This study compared the transcriptome of Atractylodes lancea rhizome at different development stages and explored genes encoding the key enzymes of the sesquiterpenoid biosynthesis pathway. Specifically, Illumina NovaSeq 6000 was employed for sequencing the cDNA libraries of A. lancea rhizome samples at the growth stage(SZ), flowering stage(KH), and harvesting stage(CS), respectively. Finally, a total of 388 201 748 clean reads were obtained, and 16 925, 8 616, and 13 702 differentially expressed genes(DEGs) were identified between SZ and KH, KH and CS, and SZ and CS, separately. Among them, 53 genes were involved in the sesquiterpenoid biosynthesis pathways: 9 encoding 6 enzymes of the mevalonic acid(MVA) pathway, 15 encoding 7 enzymes of the 2-C-methyl-D-erythritol-4-phosphate(MEP) pathway, and 29 of sesquiterpenoid and triterpenoid biosynthesis pathway. Weighted gene co-expression network analysis(WGCNA) yielded 12 genes related to sesquiterpenoid biosynthesis for the SZ, 1 gene for the KH, and 1 gene for CS, and several candidate genes for sesquiterpenoid biosynthesis were discovered based on the co-expression network. This study laid a solid foundation for further research on the sesquiterpenoid biosynthesis pathway, analysis of the regulation mechanism, and mechanism for the accumulation of sesquiterpenoids in A. lancea.
Keywords: Atractylodes lancea rhizome; sesquiterpenoid biosynthesis; transcriptome analysis.