Transcriptional Controls for Early Bolting and Flowering in Angelica sinensis

Untargeted metabolomics and functional analyses reveal that the secondary metabolite quinic acid associates with Angelica sinensis flowering

Flowering is a critical step in the plant life cycle. Angelica sinensis (Oliv.) Diels is a medicinal crop whose root is a well-known herbal medicine used in Asia. Early flowering causes changes in secondary metabolic flow and results in the loss of medicinal quality. Based on untargeted metabolomics studies, quinic acid was identified as a metabolite present in significantly higher concentrations during the early-flowering stage in A. sinensis leaves. This metabolite was subsequently investigated as a potential marker for early bolting in A. sinensis under field conditions. Moreover, quinic acid was found to accelerate flowering in the model plant Arabidopsis thaliana. Importantly, the flowering time was delayed in the quinate dehydrogenase At mutant, and this delay was reversed by quinic acid. Quinic acid upregulated the expression of the GA20OX and GID1 receptors and downregulated the expression of the inhibitor DELLA, thereby affecting the levels of FT and LFY and accelerating plant flowering. Quinic acid also significantly changed the expression of genes such as LOX, JAZ1, MYC2 and MYC3 in the jasmonic acid pathway. The trends of GID1, DELLA (GAI) and LOX2 protein expression were essentially consistent with those at the transcription level. These results suggest that quinic acid may promote plant flowering primarily by regulating the expression of genes and proteins in the gibberellin and jasmonic acid pathways.

Visualizing the Spatial Distribution of Metabolites in Angelica sinensis Roots by Matrix-Assisted Laser Desorption/Ionization Mass Spectrometry Imaging

INTRODUCTION

Angelica sinensis is one of the most popular traditional Chinese medicines (TCM) and has been extensively used to treat various diseases. Hundreds of endogenous ingredients have been isolated and identified from this herb, but their spatial distribution within the plant root is largely unknown.

OBJECTIVES

In this study, we tried to investigate and map within-tissue spatial distribution of metabolites in Angelica sinensis roots.

MATERIAL AND METHODS

After optimization of experiment conditions, the 1,5-diaminonaphthalene (1,5-DAN) was chosen as the matrix and was sprayed on the surface of root sections. Then matrix-assisted laser desorption/ionization mass spectrometry imaging (MALDI-MSI) was employed to perform in situ detection and obtain detail spatial distribution information of metabolites in Angelica sinensis roots.

RESULTS

The spatial distributions of a wide range of metabolites including organic acids, amino acids, oligosaccharides, and phospholipids were characterized and visualized in Angelica sinensis roots. Majority of these metabolites were located in the phloem and xylem, while ferulic acid was mainly present in the cork layer. The results revealed a dramatic metabolic heterogeneity among different regions of the roots and distinct spatial distribution patterns of different metabolites. Additionally, the metabolic pathways involved in the biosynthesis of choline were also successfully localized and visualized.

CONCLUSION

This study comprehensively characterized the spatial distribution of metabolites in Angelica sinensis roots, which would prompt the understanding of its chemical separation, biosynthesis, and pharmacological activities.

Light-Regulated Growth, Anatomical, Metabolites Biosynthesis and Transcriptional Changes in Angelica sinensis

Angelica sinensis is an alpine medicinal plant that has been widely used as a general blood tonic and gynecological indications over 2000 years, which depend on the bioactive metabolites (e.g., volatile oils, organic acids, and flavonoids). Although the accumulation of these metabolites is significantly affected by the environmental factors (e.g., altitude, temperature, and sunshine) as found in previous studies, the regulatory mechanism of different lights has not been clearly revealed. Here, growth parameters, contents of bioactive metabolites, and expression levels of related genes were examined when A. sinensis was exposed to different white-light (WL) and UV-B radiation treatments. The results showed that the differences in growth parameters (e.g., plant height, root length, and plant biomass) and leaf tissue characteristics (e.g., leaf thickness, stomatal density and shape, and chloroplast density) were observed under different light treatments. The contents of Z-ligustilide and ferulic acid elevated with the increase of WL (50 to 150 µmol·m2/s) and maximized under the combination of WL-100 and UV-B (107 µW/m2, UV-107) radiation, while the total flavonoids and polysaccharides contents, as well as in vitro antioxidant capacity, elevated with the increasing of WL and UV-B. mRNA transcripts encoding for the biosynthesis of volatile oils, ferulic acid, flavonoids, and polysaccharides were found to be differentially regulated under the different WL and UV-B treatments. These morphological, anatomical, and transcriptional changes are consistent with the elevated bioactive metabolites in A. sinensis under the combination of WL and UV-B. These findings will provide useful references for improving bioactive metabolite production via the cultivation and bioengineering of A. sinensis.

Stabilised pathway and its anti-oxygen/glucose deprived activity about ligustilide, as a pharmacodynamic marker for Angelica sinensis (Oliv) Diles or Ligusticum chuanxiong Hort

The inherent structural instability and very low bioavailability of ligustilide (Lig) lead to the limited application of the clinical neuroprotection. We developed a stabilised method for Lig, which involved a chemical synthesis between cyclopropylamine and the active phthalides extract included Lig from Angelica sinensis (Oliv) Diles or Ligusticum chuanxiong Hort, which generated correspond phthalide derivatives by nucleophilic additive and substituted reaction. During this process, four phthalide derivatives have been obtained and two compounds (3 and 4) stands out as an unprecedented dimeric phthalide entity. This research has meticulously delineated the stabilised pathway and polymerisation mechanics of Lig or its derivatives, an oxygen/glucose deprived PC12 cells experiment model utilised to screen the anti-stroke activity of this isolated Lig derivatives and the results supported this chemical transformation could achieve the stabilised goal of Lig and improved its anti-oxygen/glucose deprived activity.

Integrated transcriptomic and metabolomic analysis reveals the regulation mechanism of early bolting and flowering in two cultivars of Angelica sinensis

The root of Angelica sinensis is utilized in Traditional Chinese medicine to enhance blood replenishment and facilitate blood circulation. The early bolting and flowering (EBF) of A. sinensis, however, compromises the quality of the roots and restricts the yield of medicinal substances. The study was conducted to compare the transcriptomic and metabolomic profiles between EBF plants and normal plants of two cultivars of A. sinensis, followed by validation of the transcriptome results using qRT-PCR. There were 3677 DEGs in EBF plants compared to normal plants of cultivar 2 (Mingui No.2), and cultivar 4 (Mingui No.4) was 3354. The main differential metabolites in the EBF and normal plants were phenolic acids, flavonoids, lignans, and coumarins. The analysis of 5 EBF-related pathways revealed 28 genes exhibiting differential expression and 5 metabolites showing differential accumulation. The expression of the Lhcb5, Lhcb2, Lhcb6, Lhcb1, Lhca4, ATPG1, EGLC, CELB, AMY, glgA, CYCD3, SnRK2, PYL, AHK2, AUX1, BSK, FabI/K, ACACA and FabV decreased and the expression of the PsbR, PsbA, LHY, FT, CO, malQ, HK, GPI and DELLA increased in EBF plants. In addition, the Abscisic acid, d-Glucose-6P, α-d-Glucose-1P, NADP+, and ADP were more significantly enriched in EBF plants. The findings offer novel perspectives on the EBF mechanisms in A. sinensis and other medicinal plants of the Apiaceae family.

ASAP: a platform for gene functional analysis in Angelica sinensis

BACKGROUND

Angelica sinensis (Danggui), a renowned medicinal orchid, has gained significant recognition for its therapeutic effects in treating a wide range of ailments. Genome information serves as a valuable resource, enabling researchers to gain a deeper understanding of gene function. In recent times, the availability of chromosome-level genomes for A. sinensis has opened up vast opportunities for exploring gene functionality. Integrating multiomics data can allow researchers to unravel the intricate mechanisms underlying gene function in A. sinensis and further enhance our knowledge of its medicinal properties.

RESULTS

In this study, we utilized genomic and transcriptomic data to construct a coexpression network for A. sinensis. To annotate genes, we aligned them with sequences from various databases, such as the NR, TAIR, trEMBL, UniProt, and SwissProt databases. For GO and KEGG annotations, we employed InterProScan and GhostKOALA software. Additionally, gene families were predicted using iTAK, HMMER, OrholoFinder, and KEGG annotation. To facilitate gene functional analysis in A. sinensis, we developed a comprehensive platform that integrates genomic and transcriptomic data with processed functional annotations. The platform includes several tools, such as BLAST, GSEA, Heatmap, JBrowse, and Sequence Extraction. This integrated resource and approach will enable researchers to explore the functional aspects of genes in A. sinensis more effectively.

CONCLUSION

We developed a platform, named ASAP, to facilitate gene functional analysis in A. sinensis. ASAP ( www.gzybioinformatics.cn/ASAP ) offers a comprehensive collection of genome data, transcriptome resources, and analysis tools. This platform serves as a valuable resource for researchers conducting gene functional research in their projects, providing them with the necessary data and tools to enhance their studies.

A chromosome-scale genome of Peucedanum praeruptorum provide insights into Apioideae evolution and medicinal ingredient biosynthesis

Peucedanum praeruptorum Dunn, a traditional Chinese medicine rich in coumarin, belongs to the Apiaceae family. A high-quality assembled genome of P. praeruptorum is lacking, which has posed obstacles to functional identification and molecular evolution studies of genes associated with coumarin production. Here, a chromosome-scale reference genome of P. praeruptorum, an important medicinal and aromatic plant, was first sequenced and assembled using Oxford Nanopore Technologies and Hi-C sequencing. The final assembled genome size was 1.83 Gb, with a contig N50 of 11.12 Mb. The entire BUSCO evaluation and second-generation read comparability rates were 96.0 % and 99.31 %, respectively. Furthermore, 99.91 % of the genome was anchored to 11 pseudochromosomes. The comparative genomic study revealed the presence of 18,593 orthogroups, which included 476 species-specific orthogroups and 1211 expanded gene families. Two whole-genome duplication (WGD) events and one whole-genome triplication (WGT) event occurred in P. praeruptorum. In addition to the γ-WGT shared by core eudicots or most eudicots, the first WGD was shared by Apiales, while the most recent WGD was unique to Apiaceae. Our study demonstrated that WGD events that occurred in Apioideae highlighted the important role of tandem duplication in the biosynthesis of coumarins and terpenes in P. praeruptorum. Additionally, the expansion of the cytochrome P450 monooxygenase, O-methyltransferase, ATP-binding cassette (ABC) transporter, and terpene synthase families may be associated with the abundance of coumarins and terpenoids. Moreover, we identified >170 UDP-glucosyltransferase members that may be involved in the glycosylation post-modification of coumarins. Significant gene expansion was observed in the ABCG, ABCB, and ABCC subgroups of the ABC transporter family, potentially facilitating the transmembrane transport of coumarins after bolting. The P. praeruptorum genome provides valuable insights into the machinery of coumarin biosynthesis and enhances our understanding of Apiaceae evolution.

Integrating genomic and multiomic data for Angelica sinensis provides insights into the evolution and biosynthesis of pharmaceutically bioactive compounds

Angelica sinensis roots (Angelica roots) are rich in many bioactive compounds, including phthalides, coumarins, lignans, and terpenoids. However, the molecular bases for their biosynthesis are still poorly understood. Here, an improved chromosome-scale genome for A. sinensis var. Qinggui1 is reported, with a size of 2.16 Gb, contig N50 of 4.96 Mb and scaffold N50 of 198.27 Mb, covering 99.8% of the estimated genome. Additionally, by integrating genome sequencing, metabolomic profiling, and transcriptome analysis of normally growing and early-flowering Angelica roots that exhibit dramatically different metabolite profiles, the pathways and critical metabolic genes for the biosynthesis of these major bioactive components in Angelica roots have been deciphered. Multiomic analyses have also revealed the evolution and regulation of key metabolic genes for the biosynthesis of pharmaceutically bioactive components; in particular, TPSs for terpenoid volatiles, ACCs for malonyl CoA, PKSs for phthalide, and PTs for coumarin biosynthesis were expanded in the A. sinensis genome. These findings provide new insights into the biosynthesis of pharmaceutically important compounds in Angelica roots for exploration of synthetic biology and genetic improvement of herbal quality.

Apiaceae Medicinal Plants in China: A Review of Traditional Uses, Phytochemistry, Bolting and Flowering (BF), and BF Control Methods

Apiaceae plants have been widely used in traditional Chinese medicine (TCM) for the removing dampness, relieving superficies, and dispelling cold, etc. In order to exploit potential applications as well as improve the yield and quality of Apiaceae medicinal plants (AMPs), the traditional use, modern pharmacological use, phytochemistry, effect of bolting and flowering (BF), and approaches for controlling BF were summarized. Currently, about 228 AMPs have been recorded as TCMs, with 6 medicinal parts, 79 traditional uses, 62 modern pharmacological uses, and 5 main kinds of metabolites. Three different degrees (i.e., significantly affected, affected to some extent, and not significantly affected) could be classed based on the yield and quality. Although the BF of some plants (e.g., Angelica sinensis) could be effectively controlled by standard cultivation techniques, the mechanism of BF has not yet been systemically revealed. This review will provide useful references for the reasonable exploration and high-quality production of AMPs.

Investigation of Enzymes in the Phthalide Biosynthetic Pathway in Angelica sinensis Using Integrative Metabolite Profiles and Transcriptome Analysis

The roots of Angelica sinensis (Oliv.) Diels are well known for their efficacy in promoting blood circulation. Although many studies have indicated that phthalides are the main chemical components responsible for the pharmacological properties of A. sinensis, the phthalide biosynthetic pathway and enzymes that transform different phthalides are still poorly understood. We identified 108 potential candidate isoforms for phthalide accumulation using transcriptome and metabolite profile analyses. Then, six enzymes, including phospho-2-dehydro-3-deoxyheptonate aldolase 2, shikimate dehydrogenase, primary amine oxidase, polyphenol oxidase, tyrosine decarboxylase, and shikimate O-hydroxycinnamoyl transferase, were identified and proven to be involved in phthalide accumulation by heterologously expressing these proteins in Escherichia coli. We proposed a possible mechanism underlying phthalide transformation and biosynthetic pathways in A. sinensis based on our findings. The results of our study can provide valuable information for understanding the mechanisms underlying phthalide accumulation and transformation and enable further development of quality control during the cultivation of A. sinensis.

Regulatory Networks of Flowering Genes in Angelica sinensis during Vernalization

Angelica sinensis is a low-temperature and long-day perennial herb that has been widely used for cardio-cerebrovascular diseases in recent years. In commercial cultivation, up to 40% of flowering decreases the officinal yield of roots and accumulation of bioactive compounds. Although the regulatory mechanism of flowering genes during the photoperiod has been revealed, the networks during vernalization have not been mapped. Here, transcriptomics profiles of A. sinensis with uncompleted (T1), completed (T2) and avoided vernalization (T3) were performed using RNA-seq, and genes expression was validated with qRT-PCR. A total of 61,241 isoforms were annotated on KEGG, KOG, Nr and Swiss-Prot databases; 4212 and 5301 differentially expressed genes (DEGs) were observed; and 151 and 155 genes involved in flowering were dug out at T2 vs. T1 and T3 vs. T1, respectively. According to functional annotation, 104 co-expressed genes were classified into six categories: FLC expression (22; e.g., VILs, FCA and FLK), sucrose metabolism (12; e.g., TPSs, SUS3 and SPSs), hormone response (18; e.g., GID1B, RAP2s and IAAs), circadian clock (2; i.e., ELF3 and COR27), downstream floral integrators and meristem identity (15; e.g., SOC1, AGL65 and SPLs) and cold response (35; e.g., PYLs, ERFs and CORs). The expression levels of candidate genes were almost consistent with FPKM values and changes in sugar and hormone contents. Based on their functions, four pathways that regulate flowering during vernalization were mapped, including the vernalization pathway, the autonomic pathway, the age pathway and the GA (hormone) pathway. This transcriptomic analysis provides new insights into the gene-regulatory networks of flowering in A. sinensis.

Comparative Transcriptomics Unveil the Crucial Genes Involved in Coumarin Biosynthesis in Peucedanum praeruptorum Dunn

Peucedanum praeruptorum Dunn is a commonly used traditional Chinese medicine that is abundant in furano- and dihydropyrano coumarins. When P. praeruptorum reaches the bolting stage, the roots gradually lignified, and the content of coumarins declines rapidly. Non-bolting has always been a decisive factor for harvesting the P. praeruptorum materials. To evaluate the amount of coumarin components in unbolted and bolted P. praeruptorum, the variations of praeruptorin A, praeruptorin B, praeruptorin E, peucedanocoumarin I, and peucedanocoumarin II were determined. Additionally, 336,505 transcripts were obtained from the comparative transcriptome data. Among them, a total of 1,573 differentially expressed genes were screened out. To identify the critical genes involved in coumarin biosynthesis, comparative transcriptomics coupled with co-expression associated analysis was conducted. Finally, coumarin biosynthesis-related eighteen candidate genes were selected for the validation of qPCR. Additionally, a phylogenetic tree and the expression profile of ATP-binding cassette (ABC) transporters were constructed. To clarify the main genes in the regulation of coumarin biosynthesis, the interaction network of the co-expression genes from thirteen modules was constructed. Current results exhibited the significant increment of praeruptorin A, praeruptorin B and praeruptorin E in the bolted P. praeruptorum. Although, peucedanocoumarin I and peucedanocoumarin II were slightly increased. Besides the content of coumarins, the essential genes involved in the coumarin biosynthesis also exhibited an overall downward trend after bolting. Three peroxidases (PRXs) involved in the production of lignin monomers had been demonstrated to be downregulated. PAL, C4H, HCT, COMT, CCoAOMT, and some ABC transporters were dramatically downregulated at the bolting stage. These results indicated that the downregulation of coumarin biosynthetic genes in the bolted P. praeruptorum ultimately reduced the formation of coumarins. However, the mechanism through which bolting indirectly affects the formation of coumarin still needs extra functional verification.

Transcriptomic Analysis Reveals LncRNAs Associated with Flowering of Angelica sinensis during Vernalization

Angelica sinensis is a “low-temperature and long-day” perennial plant that produces bioactive compounds such as phthalides, organic acids, and polysaccharides for various types of clinical agents, including those with cardio-cerebrovascular, hepatoprotective, and immunomodulatory effects. To date, the regulatory mechanism of flowering under the photoperiod has been revealed, while the regulatory network of flowering genes during vernalization, especially in the role of lncRNAs, has yet to be identified. Here, lncRNAs associated with flowering were identified based on the full-length transcriptomic analysis of A. sinensis at vernalization and freezing temperatures, and the coexpressed mRNAs of lncRNAs were validated by qRT-PCR. We obtained a total of 2327 lncRNAs after assessing the protein-coding potential of coexpressed mRNAs, with 607 lncRNAs aligned against the TAIR database of model plant Arabidopsis, 345 lncRNAs identified, and 272 lncRNAs characterized on the SwissProt database. Based on the biological functions of coexpressed mRNAs, the 272 lncRNAs were divided into six categories: (1) chromatin, DNA/RNA and protein modification; (2) flowering; (3) stress response; (4) metabolism; (5) bio-signaling; and (6) energy and transport. The differential expression levels of representatively coexpressed mRNAs were almost consistent with the flowering of A. sinensis. It can be concluded that the flowering of A. sinensis is positively or negatively regulated by lncRNAs, which provides new insights into the regulation mechanism of the flowering of A. sinensis.

Comparative Proteomics Reveals the Difference in Root Cold Resistance between Vitis. riparia × V. labrusca and Cabernet Sauvignon in Response to Freezing Temperature

Grapevines, bearing fruit containing large amounts of bioactive metabolites that offer health benefits, are widely cultivated around the world. However, the cold damage incurred when grown outside in extremely low temperatures during the overwintering stage limits the expansion of production. Although the morphological, biochemical, and molecular levels in different Vitis species exposed to different temperatures have been investigated, differential expression of proteins in roots is still limited. Here, the roots of cold-resistant (Vitis. riparia × V. labrusca, T1) and cold-sensitive varieties (Cabernet Sauvignon, T3) at -4 °C, and also at -15 °C for the former (T2), were measured by iTRAQ-based proteomic analysis. Expression levels of genes encoding candidate proteins were validated by qRT-PCR, and the root activities during different treatments were determined using a triphenyl tetrazolium chloride method. The results show that the root activity of the cold-resistant variety was greater than that of the cold-sensitive variety, and it declined with the decrease in temperature. A total of 25 proteins were differentially co-expressed in T2 vs. T1 and T1 vs. T3, and these proteins were involved in stress response, bio-signaling, metabolism, energy, and translation. The relative expression levels of the 13 selected genes were consistent with their fold-change values of proteins. The signature translation patterns for the roots during spatio-temporal treatments of different varieties at different temperatures provide insight into the differential mechanisms of cold resistance of grapevine.

Cool Temperature Enhances Growth, Ferulic Acid and Flavonoid Biosynthesis While Inhibiting Polysaccharide Biosynthesis in Angelica sinensis

Angelica sinensis, a perennial herb that produces ferulic acid and phthalides for the treatment of cardio-cerebrovascular diseases, prefers growing at an altitude of 1800-3000 m. Geographical models have predicted that high altitude, cool temperature and sunshade play determining roles in geo-authentic formation. Although the roles of altitude and light in yield and quality have been investigated, the role of temperature in regulating growth, metabolites biosynthesis and gene expression is still unclear. In this study, growth characteristics, metabolites contents and related genes expression were investigated by exposing A. sinensis to cooler (15 °C) and normal temperatures (22 °C). The results showed that plant biomass, the contents of ferulic acid and flavonoids and the expression levels of genes related to the biosynthesis of ferulic acid (PAL1, 4CLL4, 4CLL9, C3H, HCT, CCOAMT and CCR) and flavonoids (CHS and CHI) were enhanced at 15 °C compared to 22 °C. The contents of ligustilide and volatile oils exhibited slight increases, while polysaccharide contents decreased in response to cooler temperature. Based on gene expression levels, ferulic acid biosynthesis probably depends on the CCOAMT pathway and not the COMT pathway. It can be concluded that cool temperature enhances plant growth, ferulic acid and flavonoid accumulation but inhibits polysaccharide biosynthesis in A. sinensis. These findings authenticate that cool temperature plays a determining role in the formation of geo-authentic and also provide a strong foundation for regulating metabolites production of A. sinensis.

Integrated Metabolomic and Transcriptomic Analysis Reveals Differential Mechanism of Flavonoid Biosynthesis in Two Cultivars of Angelica sinensis

Angelica sinensis is a traditional Chinese medicinal plant that has been primarily used as a blood tonic. It largely relies on its bioactive metabolites, which include ferulic acid, volatile oils, polysaccharides and flavonoids. In order to improve the yield and quality of A. sinensis, the two cultivars Mingui 1 (M1), with a purple stem, and Mingui 2 (M2), with a green stem, have been selected in the field. Although a higher root yield and ferulic acid content in M1 than M2 has been observed, the differences of flavonoid biosynthesis and stem-color formation are still limited. In this study, the contents of flavonoids and anthocyanins were determined by spectrophotometer, the differences of flavonoids and transcripts in M1 and M2 were conducted by metabolomic and transcriptomic analysis, and the expression level of candidate genes was validated by qRT-PCR. The results showed that the contents of flavonoids and anthocyanins were 1.5- and 2.6-fold greater in M1 than M2, respectively. A total of 26 differentially accumulated flavonoids (DAFs) with 19 up-regulated (UR) and seven down-regulated (DR) were obtained from the 131 identified flavonoids (e.g., flavonols, flavonoid, isoflavones, and anthocyanins) in M1 vs. M2. A total 2210 differentially expressed genes (DEGs) were obtained from the 34,528 full-length isoforms in M1 vs. M2, and 29 DEGs with 24 UR and 5 DR were identified to be involved in flavonoid biosynthesis, with 25 genes (e.g., CHS1, CHI3, F3H, DFR, ANS, CYPs and UGTs) mapped on the flavonoid biosynthetic pathway and four genes (e.g., RL1, RL6, MYB90 and MYB114) belonging to transcription factors. The differential accumulation level of flavonoids is coherent with the expression level of candidate genes. Finally, the network of DAFs regulated by DEGs was proposed. These findings will provide references for flavonoid production and cultivars selection of A. sinensis.

Bolting reduces ferulic acid and flavonoid biosynthesis and induces root lignification in Angelica sinensis

Angelica sinensis is a perennial herbaceous species that produces the bioactive metabolites ferulic acid and alkylphthalides widely applied in the treatment of cardio-cerebrovascular diseases. While the effects of bolting on plant biomass and metabolites accumulation have been partly investigated, the mechanism of bolting reducing metabolites biosynthesis is still limited. In this study, the root biomass, accumulations of ferulic acid, flavonoids and lignin, antioxidant capacity, and related genes expression at four different bolting stages were investigated. The results showed that there was a 2.2-, 2.4- and 2.9-fold decrease of the root biomass, ferulic acid and flavonoids contents, while a 2.9-fold increase of lignin content on a per plant basis during the bolting stages. The antioxidant capacity also exhibited significant decrease with growth and development. The differential expression levels of the 20 genes, which are involved in biosynthesis of ferulic acid (e.g. AsPAL1, As4CLs and AsHCT), flavonoids (e.g. AsCHS, AsCHI and AsI3'H) and lignin (e.g. AsCAD1 and AsLACs), were consistent with changes in the above metabolites accumulation. The findings will provide useful references for improving the production of bioactive metabolites in A. sinensis.