- Open Access
Comparative transcriptome analysis reveals sesquiterpenoid biosynthesis among 1-, 2- and 3-year old Atractylodes chinensis
BMC Plant Biology volume 21, Article number: 354 (2021)
Atractylodes chinensis (DC.) Koidz is a well-known medicinal plant containing the major bioactive compound, atractylodin, a sesquiterpenoid. High-performance liquid chromatography (HPLC) analysis demonstrated that atractylodin was most abundant in 3-year old A. chinensis rhizome, compared with those from 1- and 2-year old rhizomes, however, the molecular mechanisms underlying accumulation of atractylodin in rhizomes are poorly understood.
In this study, we characterized the transcriptomes from rhizomes of 1-, 2- and 3-year old (Y1, Y2 and Y3, respectively) A. chinensis, to identify differentially expressed genes (DEGs). We identified 240, 169 and 131 unigenes encoding the enzyme genes in the mevalonate (MVA), methylerythritol phosphate (MEP), sesquiterpenoid and triterpenoid biosynthetic pathways, respectively. To confirm the reliability of the RNA sequencing analysis, eleven key gene encoding factors involved in the sesquiterpenoid and triterpenoid biosynthetic pathway, as well as in pigment, amino acid, hormone and transcription factor functions, were selected for quantitative real time PCR (qRT-PCR) analysis. The results demonstrated similar expression patterns to those determined by RNA sequencing, with a Pearson’s correlation coefficient of 0.9 between qRT-PCR and RNA-seq data. Differential gene expression analysis of rhizomes from different ages revealed 52 genes related to sesquiterpenoid and triterpenoid biosynthesis. Among these, seven DEGs were identified in Y1 vs Y2, Y1 vs Y3 and Y2 vs Y3, of which five encoded four key enzymes, squalene/phytoene synthase (SS), squalene-hopene cyclase (SHC), squalene epoxidase (SE) and dammarenediol II synthase (DS). These four enzymes directly related to squalene biosynthesis and subsequent catalytic action. To validate the result of these seven DEGs, qRT-PCR was performed and indicated most of them displayed lower relative expression in 3-year old rhizome, similar to transcriptomic analysis.
The enzymes SS, SHC, SE and DS down-regulated expression in 3-year old rhizome. This data corresponded to the higher content of sesquiterpenoid in 3-year old rhizome, and confirmed by qRT-PCR. The results of comparative transcriptome analysis and identified key enzyme genes laid a solid foundation for investigation of production sesquiterpenoid in A. chinensis.
Atractylodes lancea and Atractylodes chinensis (typically referred to as “Mao Cang Zhu” and “Bei Cang Zhu” in Chinese), together constitute the rhizome atractylodes, and belong to the Asteraceae family. The rhizome atractylodes are widely used in East Asia, and have great economic and medicinal value. A. lancea is currently on the verge of extinction, therefore, A. chinensis is the main source of the rhizome atractylodes that widely distributed in most areas of Northern China. The main bioactive compounds in A. chinensis rhizome are used to treat digestive disorders, rheumatic diseases and night blindness . Modern pharmacological studies have demonstrated that A. chinensis also has anti-inflammatory, anti-bacterial [11, 20] and antitumor  properties. Although the sesquiterpenoid components of A. chinensis rhizome have important pharmacological activities, the molecular mechanisms underlying accumulation of bioactive sesquiterpenoids are poorly understood. In plants, sesquiterpenoids are generally synthesized via MVA and MEP biosynthetic pathways.
Natural populations of A. chinensis currently being rapidly depleted, due to heavy use and the weak reproductive capacity of perennial herbs. Thus, artificial cultivation is urgently needed to protect the natural populations and ensure sustainable utilization. A crucial question is how to ensure, or even improve, rhizome atractylodes quality, in terms of sesquiterpenoid content. Although the phytochemistry [11, 33], pharmacology [5, 14, 20, 24] and cultivation [27, 32, 37, 38] of A. chinensis have been studied, the molecular mechanisms underlying their accumulation of bioactive compounds remains unclear, largely due to a lack of genomic and transcriptomic data.
Transcriptome analysis is an effective approach for analysis of secondary metabolites biosynthesis, and has been used to determine the functions of genes in medicinal plants, including Danshen (Salvia miltiorrhiza) , Renshen (Panax ginseng) , Sanqi (Panax notoginseng)  and Yunnan chonglou (Paris polyphylla var. yunnanensis) , among others. Recently, understanding of the molecular processes involved in sesquiterpenoid biosynthesis has improved, with various genes involved in this biosynthetic pathway investigated by transcriptome analysis in the genus, Atractylodes [2, 12, 36]. Sesquiterpenoids are the main bioactive components in the rhizomes of A. lancea and A. chinensis; however, there are differences between them in the composition and content of sesquiterpenoids. In addition, the content of bioactive components in perennial medicinal herbs is influenced by the year of cultivation [1, 16, 29]. To date, one study has reported the transcriptome profiles of 3-year old A. chinensis rhizome ; however, there are no data regarding the molecular mechanism involved in the relationship between sesquiterpenoid accumulation and year of cultivation. Elucidating factors involved in the biosynthesis and accumulation of bioactive components and identifying key enzyme genes in the biosynthetic pathway will be important steps toward improvements in sesquiterpenoid production.
Here, rhizomes from 1-, 2- and 3-year old A. chinensis were subjected to high throughput transcriptome sequencing, enabling us to characterize the transcriptomes and differential expression profiles of A. chinensis rhizomes cultivated for different ages, to profile differentially expressed genes (DEGs) among rhizomes from different years of cultivation, and to identify DEGs related to biosynthesis and accumulation of sesquiterpenoids. Discovering the key enzyme genes in the sesquiterpenoid biosynthetic pathway is necessary to improve atractylodin production. This study could provide insights into the relationship between changes in atrctylodin content and year of cultivation, and contribute to uncovering the underlying molecular mechanisms in A. chinensis.
Materials and methods
A. chinensis seeds were collected from cultivation base of Qinhuangdao Tongsheng Pharmaceutical Co., Ltd, Qinhuangdao City, Hebei Province, China. To ensure a similar physical environment, seeds were sown separately in 2016, 2017 and 2018, in the same open field at the Hebei Normal University of Science & Technology. No permission is required to collect wild resources of A. chinensis. A. chinensis rhizome, which is the part of the plant used in medicinal preparations, serves as a store for photosynthetic products and bioactive compounds. For use as a medicine, A. chinensis is optimally harvested the 3rd to 4th rhizomes during late October to early November in Hebei province, therefore, 1-, 2- and 3-year old rhizomes with 3 biological replicates (15 rhizomes for each biological replicates) were collected in early November 2018 (as seen in Fig. 7). After collection, rhizomes were cleaned in running water, and divided into two groups. One group was dried at 60 °C for atractylodin extraction and HPLC analysis. The other group was immediately frozen in liquid nitrogen and stored at -80 °C prior to RNA extraction and sequencing.
Atractylodin extraction and HPLC analysis
The dried rhizomes were and then ground into a powder, 0.2 g of which was immersed in 50 mL methanol (purity ≥ 99.9%, Grade/Application information: ACS. Reag. Ph Eur, CSA-No. 67–56-1) and ultrasonically extracted (Power 250 W, Frequency 40 kHz) for 1 h. Next, 1 mL of the supernatant was collected and passed through a 0.22 μm microporous filter membrane (JTSF0311, Tianjin Jinteng Experiment Equipment Co. Ltd.).
Determination of atractylodin content in 1-, 2- and 3-year old rhizomes was conducted using a Thermo Fisher UltiMate 3000 UPLC system, equipped with a Uv–vis detector, on C18 Column (4.6 × 250 mm, 5 μm, Thermo Fisher). The mobile phase was methanol:water (79:21), with a flow rate of 1.0 mL∙min−1. The HPLC chromatogram was monitored at 340 nm, and the column temperature was set at 30 °C.
The standard atractylodin (number 111924–201,806, gbw114.com, China) solution of 20 μg∙ml−1 was prepared with methanol, then diluted to 1 µg∙ml−1, 2 µg∙ml−1, 5 µg∙ml−1and 10 µg∙ml−1, separately. According to the relationship between standard atractylodin content (y) and its cover area (x), the standard curve was obtained y = 4.2254 x (R2 = 0.9999). Atractylodin content was determined by standard curve. The content of atractylodin (%) = y * 50 mL (methanol) / 0.2 g (rhizome powder). The mean values of three biological replicates were calculated. Statistical significance of atractylodin contents among 1-, 2- and 3-year old rhizomes was analyzed by DPS (version 14.5).
RNA sequencing and functional annotation of unigenes
To extract total RNA, three biological replicates of rhizomes from 1-, 2- and 3-year old A. chinensis were extracted using TRIzol Reagent (Invitrogen) separately, then treated with DNase I (TaKaRa). RNA quality was tested by 1% agarose gel electrophoresis and the concentration determined using Nanodrop spectrophotometer (Thermo). RNA pools were prepared by mixing equal amounts of the three biological replicates for each age rhizome. Transcriptome data of 1-, 2- and 3-year old rhizomes were acquired using based on the Illumina hiseq Xten PE150 platform, by Novogene Co. (Beijing, China). The raw paired end reads were trimmed and quality controlled by SeqPrep (https://github.com/jstjohn/SeqPrep) and Sickle (https://github.com/najoshi/sickle) with default parameters. The sequencing data with high-quality reads were gathered using fastp (version 0.19.5, https://github.com/OpenGene/fastp). In short, adapter’s contamination, bases with low quality, reads having (≥ 10%) ambiguous bases, and reads having ˂20 bp were removed. Then clean data from all three age rhizomes were used to do de novo assembly with Trinity (version 2.8.5, http://trinityrnaseq.sourceforge.net/). The datasets generated and analyzed during the current study are available in the Sequence Read Archive (SRA) repository (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/sra/PRJNA698794). The transcriptome data have been uploaded to SRA (BioProject ID PRJNA698794, https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/sra/PRJNA698794).
Sequence annotation and classification were referenced to the method of Ramya et al. . For functional annotation, all the assembled transcripts were searched against the Nr (NCBI protein non-redundant), COG (Clusters of Orthologous Groups of proteins) and KEGG (Kyoto Encyclopedia of Genes and Genomes) databases using BLASTx to identify the proteins that had the highest sequence similarity with the given transcripts to retrieve their function annotations and a typical cut-off E-values less than 1.0 × 10–5 was set. BLAST2GO (http://www.blast2go.com/b2ghome) program was used to get GO (Gene Ontology) annotations of unique assembled transcripts. Metabolic pathway analysis was performed using the KEGG (http://www.genome.jp/kegg/). For other sequences not involved in the BLAST searches, we used the TransDecoder program (https://github.com/TransDecoder/TransDecoder) to predict coding sequence (CDS) and orientation.
Analysis of differentially expressed genes
Differential expression of unigenes among three age rhizomes of A. chinensis were determined using DESeq2 (Version 1.24.0) software. DESeq2 provides statistical routines for determining differential expression in digital gene expression data using a model based on the negative binomial distribution. For functional-enrichment analysis including GO and KEGG were performed to identify which DEGs were significantly enriched in GO terms and KEGG metabolic pathways at Bonferroni-corrected P-adjust < 0.05 compared with the whole-transcriptome background. GO functional enrichment and KEGG pathway analysis were carried out by Goatools (version 0.6.5, https://github.com/tanghaibao/Goatools) and KOBAS (http://kobas.cbi.pku.edu.cn/kobas3/?t=1).
Differences in gene expression were evaluated using the chi-square test and the false discovery rate (FDR) was controlled. The FPKM (fragments per kilobase of exons model per million mapped and reads) were assigned to normalize reads expression values. This study took FDR value of (≤ 0.05) and (log2FC ≥ 1) as a criterion for screening DEGs. Corrected P-adjust < 0.05 were used as thresholds for “enriched” DEGs. Heat maps were generated to display genes with significantly altered expression at the three ages.
Quantitative real-time PCR
To confirm the reliability of the RNA sequencing analysis, qRT-PCR analyses were performed using samples from the same 1-, 2- and 3-year old rhizomes as used for RNA-seq. Eleven genes (cluster-15114.3, cluster-8388.71372, cluster-8388.203329, cluster-388.168445, cluster-8388.64828, cluster-8388.299573, cluster-8388.162261, cluster-8388.157231, cluster-8388.172353, cluster-8388.295361 and cluster-8388.295722), with key functions in sesquiterpenoid biosynthetic pathway, as well as in pigment, amino acid, hormone, and transcription factor functions, were randomly selected for qRT-PCR analysis. Primers for qRT-PCR were designed using Primer v5.0 and synthesized by Sangon Biotech (Shanghai, China) Co., Ltd. cDNAs were reverse-transcribed from total RNA using the PrimeScript RT reagent Kit (TaKaRa), and qRT-PCR analyses were performed on an BIO-RAD CFX Connect™ Real-Time PCR detect System (Singapore Biosystems). Relative expression data were normalized to those of the UBQ2 gene, which was used as an internal control . Each qRT-PCR experiment was repeated three times. The relative expression of each gene was calculated using the 2−∆∆Ct method  and the SD was calculated among three biological replicates. All primers used are listed in Supplementary Table S1.
Validation of the seven DEGs related to sesquiterpenoid and triterpenoid biosynthesis using qRT-PCR according to method of reliability confirmation for the RNA sequencing. The primers are listed in Supplementary Table S2.
The atrctylodin content in 1-, 2- and 3-year old A. chinensis rhizomes
As atrctylodin is the main and index bioactive constituent in A. chinensis, its levels in rhizome from 1-, 2- and 3-year old A. chinensis plants were measured by HPLC analysis, with atrctylodine contents (%) recorded as 0.2252, 0.2378 and 0.2939, respectively (Table 1, Supplementary Fig. S1). The atractylodin content in 3-year old rhizome was significantly higher than 1- and 2-year old rhizomes (Table 1). These data showed that cultivation year had marked effect on atrctylodin content of A. chinensis rhizome, however, the molecular mechanisms underlying the higher atrctylodin content in 3-year old rhizome is unclear.
Sequencing analysis and de novo assembly
To study the molecular mechanisms involved in the relationship between increased atrctylodin content and A. chinensis cultivation year, transcriptome sequencing was conducted. After filtering out adapter sequences and reads ≤ 50 bp, 58,394,019, 51,583,471, and 59,107,505 high-quality (HQ) reads were obtained from 1-, 2- and 3-year old A. chinensis rhizomes, respectively. Reads from the three samples were also pooled and the above steps repeated, resulting in identification of 143,616 unigenes (mean length 825 bp, N50 length 1121 bp). The GC content of the reads and unigenes was in the range 45.18% ~ 45.56% (Table 2). Analysis of length distribution demonstrated that 16% unigenes were > 1 kb.
Functional annotation and classification
A total of 56,759 unigenes (39.52% of the total assembled unigenes) had matches in the Nr database, with 37,475 (26.09%), 31,272 (21.77%), 6,540 (4.55%), 39,372 (27.41%) and 31,424 (14.92%) unigenes showing significant similarity to sequences in the Swiss-Prot, Pfam, COG, GO and KEGG databases, respectively (Table 3).
The counts of the unigenes in the three GO categories included biological process (57,008 unigenes), cellular component (55,906 unigenes) and molecular function (47,730 unigenes) with 47 sub categories (Fig. 1A). In total, 6,540 unigenes were annotated and grouped into 24 COG classifications (Fig. 1B), among which, the cluster for “translation, ribosomal structure and biogenesis” (464 unigenes, 7.10% of the total COG annotated unigenes) accounted for the largest proportion, followed by “posttranslational modification, protein turnover, chaperones” (377 unigenes, 5.80%) and “general function prediction only” (373 unigenes, 5.70%).
KEGG pathway analysis was performed to functionally classify biochemical pathways active in 1-, 2- and 3-year old A. chinensis rhizomes. A total of 31,424 unigenes were assigned to 6 KEGG categories with 130 sub categories: “metabolism”, “genetic information processing”, “environmental information processing”, “cellular processes”, “organismal systems” and “Human diseases” (Fig. 2).
Of KEGG secondary metabolic pathways, most unigenes were assigned to “phenylpropanoid biosynthesis” (723 unigenes, ko00940), “terpenoid backbone biosynthesis” (517 unigenes, ko00900), “carotenoid biosynthesis” (394 unigenes, ko00906), and “sesquiterpenoid and triterpenoid biosynthesis” (204 unigenes, ko00909) (Table 4). A total of 240 unigenes were identified as key enzyme genes in the MVA pathway, with 169 unigenes in the MEP pathway, and 131 unigenes in the sesquiterpenoid and triterpenoid biosynthetic pathway (Table 5). The discovery of these genes related to sesquiterpenoid and triterpenoid biosynthetic pathway may help us to elucidate the molecular mechanisms underlying the higher atrctylodin content in 3-year old rhizomes.
Differential expression of transcripts in A. chinensis rhizomes from different cultivation year
To compare the unigenes from different age A. chinensis rhizomes, a Venn diagram was constructed (Fig. 3). The results showed that 31,895 unigenes were shared by all three age rhizomes. A total of 7,027, 8,879 and 6,109 unigenes were specific to 1-, 2- and 3-year old A. chinensis rhizomes, respectively, with the 2-year-old A. chinensis rhizome having the highest number of unique unigenes.
To identify DEGs among the three age rhizomes, the tag frequencies of 1- vs 2-year-old (Y1 vs Y2) rhizome, 2- vs 3-year old (Y2 vs Y3) rhizome and 1- vs 3-year old (Y1 vs Y3) rhizomes were assessed, with 6,699, 2,840, and 3,633 DEGs detected between the three pair comparisons, respectively (Fig. 4). Y1 vs Y2, Y2 vs Y3 and Y1 vs Y3, revealed 3,880, 2,214 and 2,292 up-regulated genes. There were more up-regulated than down-regulated genes in Y1 vs Y2, with the opposite detected in the Y2 vs Y3 and Y1 vs Y3 comparisons.
KEGG pathway enrichment analysis of all DEGs was performed to characterize the complex biological behaviors. The enriched pathways are presented in Fig. 5, and reflected the preferential biological functions of samples from different age rhizomes. Hierarchical clustering of all DEGs indicated that overall unigenes enrichment characteristics were similar between the Y1 vs Y2 and Y1 vs Y3 comparisons (Fig. 5A and 5B), with genes involved in “Carbohydrate metabolism”, “Signal transduction”, “Amino acid metabolism”, “Lipid metabolism” and “Biosynthesis of secondary metabolites” over-expressed. In Y2 vs Y3, genes involved in “Lipid metabolism”, “Amino acid metabolism”, “Biosynthesis of secondary metabolites” and “Replication and repair” were overexpressed (Fig. 5C).
Pathways involved in bioactive compounds metabolism are of particular interest in medicinal plants. DEGs involved in “Biosynthesis of secondary metabolites” were overexpressed in all three pair comparisons; 52 genes related to sesquiterpenoid and triterpenoid biosynthesis were detected, of which seven were differentially expressed in Y1 vs Y2, Y1 vs Y3 and Y2 vs Y3 (Fig. 6).
Heatmap trees were constructed based on gene expression levels, to further investigate the seven differentially expressed genes involved in sesquiterpenoid and triterpenoid biosynthetic pathway, including NAD-dependent epimerase/dehydratase (NDE), squalene/phytoene synthase (SS), squalene-hopene cyclases (SHC), squalene epoxidase (SE), dammarenediol-II synthase (DS), and serine/threonine-protein kinase SRK2E (SPK). All of these seven DEGs down-regulated expression in 3-year old rhizome, comparing with 1- and 2-year old rhizomes (Fig. 7). Notably, five of these DEGs encoded four key enzymes: SS, SHC, SE and DS. These four enzymes directly related to squalene biosynthesis and subsequent catalytic action. According to putative pathway, squalene is the first precursor in triterpenoid biosynthetic pathway (Fig. 8).
Validation of RNA-seq analysis by qRT-PCR
To confirm the reliability of the RNA sequencing analysis, eleven genes representing key genes in sesquiterpenoid and triterpenoid biosynthetic pathway, as well as in pigment, amino acid, hormone and transcription factor functions, were selected for qRT-PCR analysis. The result demonstrated similar expression patterns to those determined by RNA sequencing, with a Pearson’s correlation co-efficient between qRT-PCR and RNA-seq data of 0.9 (Fig. 9).
Further validation of seven DEGs related to sesquiterpenoid and triterpenoid biosynthetic pathway was performed by qRT-PCR. The relative expression levels of these seven DEGs noted in 3-year old rhizome were significantly lower than those in 1- and 2-year old rhizomes (Fig. 10). These results are consistent with the data of transcriptomic sequencing analysis.
As genome data for the Atractylodes genus are not yet available, Illumina-based RNA sequencing was performed to characterize the A. chinensis transcriptome. We obtained 143, 616 unigenes, of which 40.71% could be functionally annotated based on public databases. In addition, the qRT-PCR results demonstrated similar expression patterns of eleven randomly selected genes to those determined by RNA-seq analysis, demonstrating the reliability of our A. chinensis transcriptome data.
Transcriptomic analysis to investigate sesquiterpenoid accumulation patterns in different tissues (leaf and rhizome) of A. lancea discovered 69 unigenes in the MVA pathway, including nine key enzymes, and 28 unigenes in the MEP pathway, involving seven key enzymes . In this study, we investigated the sesquiterpenoid accumulation patterns in A. chinensis of different cultivation year and discovered 240 unigenes in the MVA pathway, involving twelve key enzymes, and 169 unigenes in the MEP pathway, involving six crucial enzymes. It inferred more differences between different age rhizomes than different tissues from the same individual. These data will facilitate further study of the molecular mechanisms underlying sesquiterpenoid accumulation.
In this study, we found that atractylodin content in A. chinensis rhizome increased with the increase of cultivation year. Li et al.  confirmed that the age of cultivation medicinal plants was important in increasing saponins production in Panax notoginseng rhizomes. Based on this natural phenomena, we performed differential expression analysis using transcriptome data from 1-, 2- and 3-year old rhizomes, to identify candidate DEGs encoding key enzymes in sesquiterpenoid and triterpenoid biosynthetic pathway. Differential gene expression patterns were further investigated to profile global gene expression differences between Y1 vs Y2, Y2 vs Y3 and Y1 vs Y3. Most DEGs between Y1 vs Y2 and Y1 vs Y3 were assigned to 19 metabolic pathways, including signal transduction, primary metabolic pathways (carbohydrate metabolism, amino acid metabolism and lipid metabolism), and biosynthesis of other secondary metabolites. In Y2 vs Y3, DEGs were assigned to 10 metabolic pathways, of which lipid metabolism, amino acid metabolism, replication and repair, and biosynthesis of other secondary metabolites comprised a higher percentage. These data indicate that the metabolic characteristics of 2-year old rhizome are more similar to those of 3-year old rhizome, relative to 1-year old rhizome. Further, the metabolic characteristics of DEGs were consistent with the rhizome’s physiological function as a storage organ for photosynthetic products and bioactive compounds. These data demonstrated that vitality of medicinal plants and the production of secondary metabolic became increased over the cultivation year, likely because they are crucial for defense against stress in older plants.
Further analysis of DEGs provided information crucial for investigation of the molecular mechanisms involved in sesquiterpenoid biosynthesis and accumulation in A. chinensis. Seven key genes related to sesquiterpenoid and triterpenoid biosynthesis were discovered by analysis DEGs between Y1 vs Y2, Y1 vs Y3, and Y2 vs Y3. Of the seven key genes, five encoding four enzymes: SE, SHC, SS and DS. The biological production of sesquiterpenoid and triterpenoid is an extremely complicated process, with synthesis occurring via MEP and MVA pathway. Many enzymes are involved in the process of isoprenenyl diphosphate (IPP) biosynthesis catalysis, which was then catalyzed toward two biosynthesis branch, sesquiterpenoids biosynthesis branch and squalene biosynthesis branch.
The identified four enzymes, SE, SHC, SS and DS, play important role in squalene biosynthesis and the subsequent catalytic reactions in this biosynthesis branch. The enzyme SS as a key enzyme in the terpenoid biosynthetic pathway catalyzes the synthesis of the first precursor of terpenoid compounds, squalene [7, 15, 23, 35, 39]. The SHC enzyme can catalyze the formation of hopene from its precursor squalene [21, 25], toward triterpenoid or steroid biosynthesis. The enzyme, SE, which catalyzes the oxidation of squalene to 2, 3-oxysqualene, is a rate-limiting enzyme in the sterol biosynthesis . DS enzyme is the first dedicated enzyme for ginsenoside biosynthesis, one of triterpenoid compounds .
The enzymes SS and SE has been found the rate-limiting enzymes involved in sterol and cholesterol biosynthesis [7, 26]. The substrate for SS is farnesyl diphosphate (FPP), which was the important substrate of sesquiterpenoid and triterpenoid biosynthesis. In the case of suppression of enzyme SS activity was observed induction of sesquiterpenoid cyclase, toward the sesquiterpenoid biosynthetic pathway . Evidence suggested that inhibition of either SS or SE was found to trigger a severalfold increase in enzyme activity of HMGR . HMGR is the first rate-controlling enzyme for the synthesis of variety of isoprenoids via MVA pathway [8, 10, 30]. This study revealed that the four enzymes SS, SHC, SE and DS down-regulated expression in 3-year old rhizome. This data corresponded to the higher content of sesquiterpenoid in 3-year old rhizome, and confirmed by qRT-PCR. We would like to infer that the sesquiterpenoid biosynthesis branch is the main biosynthetic pathway in 3-year old rhizome of A. chinensis. This study reported the results of comparative transcriptome analysis and identified key enzyme genes, laid a solid foundation for investigation of production sesquiterpenoid in A. chinensis.
Availability of data and materials
The datasets analyzed during the current study are available in the SRA (BioProject ID PRJNA698794, https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/sra/PRJNA698794) repository. The public access to the databases is open.
- A. chinensis :
Atractylodes chinensis (DC.) Koids.
Differentially expressed genes
Quantitative real-time PCR
3-Hydroxy-3-methylglutaryl coenzyme A synthase
3-Hydroxy-3-methylglutaryl coenzyme A reductase
Geranyl diphosphate synthase
Farnesyl diphosphate synthase
Mevalonate pyrophosphate decarboxylase
Isopentenyl diphosphate isomerase
Farnesyl diphosphate synthase
2-C-Methyl-D-erythritol 4-phosphate cytidylyl transferase
4-Hydroxy-3-methyl but-2-(E)-enyl-diphosphate synthase
4-Hydroxy-3-methyl but-2-(E)-enyl-diphosphate reductase
Germacrene D synthase
Germacrene A synthase
Agnieszka G, Mariola D, Anna P, Piotr K, Natalia W, Aneta S , et al. Qualitative and quantitative analyses of bioactive compounds from ex vitro Chamaenerion angustifolium (L.) (Epilobium augustifolium) herb in different harvest times. Ind. Crop. Prod. 2018;123:208–220.
Ahmed S, Zhan CS, Yang YY, Wang XK, Yang TW, Zhao ZY, et al. The transcript profile of a traditional chinese medicine, Atractylodes lancea, Revealing its sesquiterpenoid biosynthesis of the major active components. Plos One. 2016;11(3):e0151975.
Chen F, Wei YX, Zhang JM, Sang XM, Dai CC. Transcriptomics analysis investigates sesquiterpenoids accumulation pattern in different tissues of Atractylodes lancea (Thunb.) DC. plantlet. Plant Cell Tiss. Organ. Cult. 2017;130:73–90.
Chen S, Luo H, Li Y, Sun Y, Wu Q, Niu Y, et al. 454 EST analysis detects genes putatively involved in ginsenside biosynthesis in Panax ginseng. Plant Cell Rep. 2011;30(9):1593–601.
Cheng Y, Chen TY, Yang XL, Xue JH, Chen JJ. Atractylodin induces apoptosis and suppresses metastasis in hepatic cancer cells and inhibits growth in vivo. Cancer Manage Res. 2019;11:5883–94.
Committee SP. Pharmacopoeia of the People’s Republic of China. Beijing: People’s Medical Publishing House; 2020. p. 168–9.
Dan J, Rong QX, Chen YJ, Yuan QJ, Ye S, Guo J, et al. Molecular cloning and functional analysis of squalene synthase (SS) in Panax notoginseng. Int J Biol Macromol. 2017;95:658–66.
Devi K, Patar L, Modi M, Sen P. An insight into structure, function, and expression analysis of 3-Hydroxy-3-Methylglutaryl-CoA Reductase of Cymbopogon winterianus. Bioinform Biol Insights. 2017;11:1–11.
Gao XY, Zhang X, Chen W, Li J, Yang WJ, Zhang XW, et al. Transcriptome analysis of Paris polyphylla var. yunnanenesis illuminates the biosynthesis and accumulation of steroidal saponins in rhizomes and leaves. Phytochemistry. 2020;178:112460.
Gu W, Geng C, Xue W, Wu Q, Chao J, Xu F, et al. Characterization and function of the 3-hydroxy-3-methylglutaryl-CoA reductase gene in Alisma orientale (Sam.) Juz. and its relationship with protostane triterpene production. Plant Physiol. Biochem. 2015;97:378–389.
Hossen MJ, Chou JY, Li SM, Fu XQ, Yin C, Guo H, et al. An ethanol extract of the rhizome of Atractylodes chinensis exerts anti gastritis activities and inhibits Akt/NF-kappa B signaling. J Ethnopharmacol. 2019;228:18–25.
Huang QQ, Huang X, Deng J, Liu HG, Liu YW, Yu K, et al. Differential gene expression between leaf and rhizome in Atractylodes lancea: a comparative transcriptome analysis. Front Plant Sci. 2016;7:348.
Ishii T, Okuyama T, Noguchi N, Nishidono Y, Okumura T, Kaibori M, et al. Antiinflammatory constituents of Atractylodes chinensis rhizome improve glomerular lesions in immunoglobulin A nephropathy model mice. J Nat Med. 2020;74(3):616.
Kim JK, Doh EJ, Lee G. Chemical differentiation of genetically identified Atractylodes japonica, A. macrocephala, and A. chinensis rhizomes using high-performance liquid chromatography with chemometric analysis. J. Evidence-based complementary Altern. Med. 2018;2018:4860371.
Kim OT, Seong NS, Kim MY, Hwang B. Isolation and characterization of squalene synthase cDNA from Centella asiatica (L) urban. J Plant Biol. 2005;48:263–9.
Kong D, Li Y, Bai M, Deng Y, Liang G, Wu H. A comparative study of the dynamic accumulation of polyphenol components and the changes in their antioxidant activities in diploid and tetraploid Lonicera japonica. Plant Physiol Biochem. 2017;112:87–96.
Li J, Ma L, Zhang ST, Zuo CL, Song N, Zhi SS, Wu JS. Transcriptome analysis of 1- and 3-year-old Panax notoginseng rhizomes and functional characterization of saponin biosynthetic genes DS and CYP716A47 like. Planta. 2019;249:1229–37.
Liu MH, Yang BR, Cheung WF, Yang KY, Zhou HF, Kwok JSL, et al. Transcriptome analysis of leaves, roots and flowers of Panax notoginseng identifies genes involved in ginsenoside and alkaloid biosynthesis. BMC Genomics. 2015;16:265.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2-∆∆C(T) method. Methods. 2001;25:402–8.
Lyu Z, Ji XF, Chen G, An BY. Atractylodin ameliorates lipopolysaccharide and D-galactosamine-induced acute liver failure via the suppression of inflammation and oxidative stress. Inter Immunopharmacol. 2019;72:348–57.
Nakano C, Watanabe T, Minamino M, Hoshino T. Enzymatic syntheses of novel carbocyclic scaffolds with a 6,5 + 5,5 ring system by squalene hopene cyclase. Org Biomol Chem. 2019;17:9375–89.
Ramya M, Park PH, Chuang YC, Kwon OK, An HR, Park PM, et al. RNA sequencing analysis of Cymbidium goeringii identifies floral scent biosynthesis related genes. BMC Plant Biol. 2019;19:337.
Shao CM, Wang CK, Zhang SX, Shi YY, MA KL, Yang QS, et al. Transcriptome analysis of Clinopodium gracile (Benth.) Matsum and identification of genes related to Triterpenoid Saponin biosynthesis. BMC Genomics. 2020;21:49.
Shimato Y, Ota M, Asai K, Atsumi T, Tabuchi Y, Makino T. Comparison of byakujutsu (Atractylodes rhizome) and sojutsu (Atractylodes lancea rhizome) on anti-inflammatory and immunostimulative effects in vitro. J Nat Med. 2018;72:192–201.
Siedenburg G, Jendrossek D. Squalene-hopene cyclases. Appl Environ Microbiol. 2011;77:3905–15.
Sui ZH, Zhou JH, Cheng ZJ, Lu PH. Squalene epoxidase (SQLE), one of the rate-limiting enzymes in the cholesterol biosynthesis, recently has been found to be involved in the tumorigenesis. Tumor Biol. 2015;36:6173–9.
Sun WM, Wen XL, Qi HX, Feng LN, Cao J, Han ZL, et al. First Report of Anthracnose of Atractylodes chinensis (DC.) Koidz. caused by Colletotrichum chlorophyti in China. Plant disease. 2019;103:764–764.
Tansakul P, Shibuya M, Kushiro T, Ebizuka Y. Dammarenediol-II synthase, the first dedicated enzyme for ginsenoside biosynthesis, in Panax ginseng. FEBS Lett. 2006;580:5143–9.
Wang YZ, Li P. Effect of cultivation years on saponins in Paris Polyphylla var. yunnanensis using ultra-high liquid chromatography–tandem mass spectrometry and Fourier transform infrared spectroscopy. Plant Growth Regul. 2018;84:373–381.
Wei H, Xu C, Movahedi A, Sun WB, L DW, Zhuge Q. Characterization and function of 3-Hydroxy-3-Methylglutaryl-CoA Reductase in Populus trichocarpa: overexpression of PtHMGR enhances terpenoids in transgenic poplar. Front. Plant Sci. 2019;10:1476.
Wentzinger LF, Bach TJ, Hartmann MA. Inhibition of squalene synthase and squalene epoxidase in tobacco cells triggers an up-regulation of 3-Hydroxy-3-Methylglutaryl coenzyme A reductase. Plant Physiol. 2002;130:334–46.
Xu HJ, Zhou RJ, Fu JF, Yuan Y, Ge XX, Damm U. Colletotrichum atractylodicola sp nov.: the anthracnose pathogen of Atractylodes chinensis in China. Mycological progress. 2018;17:393–402.
Xu J, Chen D, Liu C, Wu XZ, Dong CX, Zhou J. Structural characterization and anti-tumor effects of an inulin-type fructan from Atractylodes chinensis. Int J Boil Macromol. 2016;82:765–71.
Yang L, Ding G, Lin H, Cheng HN, Kong Y, Wei YK, et al. Transcriptome analysis of medicinal plant Salvia miltiorrhiza and identification of genes related to tanshinone biosynthesis. Plos One. 2013;8:e80464.
Ye Y, Wang RF, Jin L, Shen JH, Li XT, Yang T, et al. Molecular cloning and differential expression analysis of a squalene synthase gene from Dioscorea zingiberensis an important pharmaceutical plant. Mol Biol Rep. 2014;41:6097–104.
Zhao JH, Zhao CY, Sun CZ, Shi FY, Chen LN, Zheng JS. Transcriptomic analysis of Atractylodes chinensis and elucidation of genes in sesquiterpenes biosynthesis. Plant Physiol. J. 2020;56:1458–1466. (in Chinses with English abstract, doi: https://0-doi-org.brum.beds.ac.uk/10.13592/j.cnki.ppj.2019.0560)
Zheng JS, Wang WP, Li YS. Effects of temperature and substrate water content on seed germination and seedling morphogenesis of Atractylodes chinensis. J. Chinses Medi. Materi. 2018;41:1282–1284. (in Chinses with English abstract, doi: https://0-doi-org.brum.beds.ac.uk/10.13863/j.issn1001-4454.2018.06.007)
Zheng JS., Wang WP, Wu YG. Effects of different sowing depth and seedling substrate on seed emergence of Atractylodes chinensis. J. Chinses Medi. Materi. 2019; 41:2501–2504. (in Chinses with English abstract, doi: https://0-doi-org.brum.beds.ac.uk/10.13863/j.issn1001-4454.2018.12.005 )
Zheng ZJ, Cao XY, Li CG, Chen YQ, Yuan B, Xu Y, et al. Molecular cloning and expression analysis of a squalene synthase gene from a medicinal plant Euphorbia pekinensis Rupr. Acta Physiol Plant. 2013;35:3007–14.
Zook MN, Kuć JA. Induction of sesquiterpene cyclase and suppression of squalene synthase activity in elicitor-treated or fungal-infected potato tuber tissue. Physiol Mol Plant Pathol. 1991;39:377–90.
Thanks for the help from laboratory of college of horticulture science and technology, Hebei Normal University of Science & Technology in atractylodin extraction and HPLC analysis.
This work was supported financially by National Fund of Hebei Province, China (Project No. H2019407120), which provide support for the test of transcriptomic sequencing, Key Research and Development Project of Hebei Province, China (Project No.19226354D) and Science and Technology Project Hebei Education Department, China (Project No. QN2019162), which provide support for the test of atrctylodin. There is no role of the funding body in the design of the study.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Primers of qRT-PCR for validation of the reliability of RNA-seq analysis. Table S2. Primers of qRT-PCR for validation of the seven DEGs involved in sesquiterpenoid and triterpenoid biosynthetic pathway.
The content of atractylodin in A. chinensis rhizomes. A. standard (20 mg∙ml−1); B. 1-year old rhizome; C. 2-year old rhizome; D. 3-year old rhizome.
About this article
Cite this article
Zhao, J., Sun, C., Shi, F. et al. Comparative transcriptome analysis reveals sesquiterpenoid biosynthesis among 1-, 2- and 3-year old Atractylodes chinensis. BMC Plant Biol 21, 354 (2021). https://0-doi-org.brum.beds.ac.uk/10.1186/s12870-021-03131-1
- Atractylodes chinensis (DC.) Koids.
- Differentially expressed genes