RNA-Sequencing Reveals the Involvement of Sesquiterpene Biosynthesis Genes and Transcription Factors during an Early Response to Mechanical Wounding of Aquilaria sinensis

Plants respond to wounding by reprogramming the expression of genes involved in secondary metabolism. Aquilaria trees produce many bioactive secondary metabolites in response to wounding, but the regulatory mechanism of agarwood formation in the early response to mechanical wounding has remained unclear. To gain insights into the process of transcriptome changes and to determine the regulatory networks of Aquilaria sinensis to an early response (15 days) to mechanical wounding, we collected A. sinensis samples from the untreated (Asc1) and treated (Asf1) xylem tissues and performed RNA sequencing (RNA-seq). This generated 49,102,523 (Asc1) and 45,180,981 (Asf1) clean reads, which corresponded to 18,927 (Asc1) and 19,258 (Asf1) genes, respectively. A total of 1596 differentially expressed genes (DEGs) were detected in Asf1 vs. Asc1 (|log(2) (fold change)| ≥ 1, P(adj) ≤ 0.05), of which 1088 were up-regulated and 508 genes were down-regulated. GO and KEGG enrichment analysis of DEGs showed that flavonoid biosynthesis, phenylpropanoid biosynthesis, and sesquiterpenoid and triterpenoid biosynthesis pathways might play important roles in wound-induced agarwood formation. Based on the transcription factor (TF)-gene regulatory network analysis, we inferred that the bHLH TF family could regulate all DEGs encoding for farnesyl diphosphate synthase, sesquiterpene synthase, and 1-deoxy-D-xylulose-5-phosphate synthase (DXS), which contribute to the biosynthesis and accumulation of agarwood sesquiterpenes. This study provides insight into the molecular mechanism regulating agarwood formation in A. sinensis, and will be helpful in selecting candidate genes for improving the yield and quality of agarwood.