De Novo Transcriptome Assembly of Cedar (Cedrela odorata L.) and Differential Gene Expression Involved in Herbivore Resistance

Abstract

Timber trees are targets of herbivorous attacks. The identification of genes associated with pest resistance can be accomplished through differential expression analysis using transcriptomes. We reported the de novo assembly of cedar (Cedrela odorata L.) transcriptome and the differential expression of genes involved in herbivore resistance. The assembly and annotation of the transcriptome were obtained using RNAseq from healthy cedar plants and those infested with Chrysobothris yucatanensis. A total of 325.6 million reads were obtained, and 127,031 (97.47%) sequences were successfully assembled. A total of 220 herbivory-related genes were detected, of which 170 genes were annotated using GO terms, and 161 genes with 245 functions were identified—165, 75, and 5 were molecular functions, biological processes, and cellular components, respectively. To protect against herbivorous infestation, trees produce toxins and volatile compounds which are modulated by signaling pathways and gene expression related to molecular functions and biological processes. The limited number of genes identified as cellular components suggests that there are minimal alterations in cellular structure in response to borer attack. The chitin recognition protein, jasmonate ZIM-domain (JAZ) motifs, and response regulator receiver domain were found to be overexpressed, whereas the terpene synthase, cytochrome P450, and protein kinase domain gene families were underexpressed. This is the first report of a cedar transcriptome focusing on genes that are overexpressed in healthy plants and underexpressed in infested plants. This method may be a viable option for identifying genes associated with herbivore resistance.