Integrative analysis of the metabolome and transcriptome provides insights into the mechanisms of lignan biosynthesis in Herpetospermum pedunculosum (Cucurbitaceae).

BACKGROUND: Herpetospermum pedunculosum (Ser.) C. B. Clarke is a traditional Chinese herbal medicine that heavily relies on the lignans found in its dried ripe seeds (Herpetospermum caudigerum), which have antioxidant and hepatoprotective functions. However, little is known regarding the lignan biosynthesis in H. pedunculosum. In this study, we used metabolomic (non-targeted UHPLC-MS/MS) and transcriptome (RNA-Seq) analyses to identify key metabolites and genes (both structural and regulatory) associated with lignan production during the green mature (GM) and yellow mature (YM) stages of H. pedunculosum. RESULTS: The contents of 26 lignan-related metabolites and the expression of 30 genes involved in the lignan pathway differed considerably between the GM and YM stages; most of them were more highly expressed in YM than in GM. UPLC-Q-TOF/MS confirmed that three Herpetospermum-specific lignans (including herpetrione, herpetotriol, and herpetin) were found in YM, but were not detected in GM. In addition, we proposed a lignan biosynthesis pathway for H. pedunculosum based on the fundamental principles of chemistry and biosynthesis. An integrated study of the transcriptome and metabolome identified several transcription factors, including HpGAF1, HpHSFB3, and HpWOX1, that were highly correlated with the metabolism of lignan compounds during seed ripening. Furthermore, functional validation assays revealed that the enzyme 4-Coumarate: CoA ligase (4CL) catalyzes the synthesis of hydroxycinnamate CoA esters. CONCLUSION: These results will deepen our understanding of seed lignan biosynthesis and establish a theoretical basis for molecular breeding of H. pedunculosum.

Transcriptomic Insights and Cytochrome P450 Gene Analysis in Kadsura coccinea for Lignan Biosynthesis.

Kadsura coccinea is a medicinal plant from the Schisandraceae family that is native to China and has great pharmacological potential due to its lignans. However, there are significant knowledge gaps regarding the genetic and molecular mechanisms of lignans. We used transcriptome sequencing technology to analyze root, stem, and leaf samples, focusing on the identification and phylogenetic analysis of Cytochrome P450 (CYP) genes. High-quality data containing 158,385 transcripts and 68,978 unigenes were obtained. In addition, 36,293 unigenes in at least one database, and 23,335 across five databases (Nr, KEGG, KOG, TrEMBL, and SwissProt) were successfully annotated. The KEGG pathway classification and annotation of these unigenes identified 10,825 categorized into major metabolic pathways, notably phenylpropanoid biosynthesis, which is essential for lignan synthesis. A key focus was the identification and phylogenetic analysis of 233 Cytochrome P450 (CYP) genes, revealing their distribution across 38 families in eight clans, with roots showing specific CYP gene expression patterns indicative of their role in lignan biosynthesis. Sequence alignment identified 22 homologous single genes of these CYPs, with 6 homologous genes of CYP719As and 1 of CYP81Qs highly expressed in roots. Our study significantly advances the understanding of the biosynthesis of dibenzocyclooctadiene lignans, offering valuable insights for future pharmacological research and development.

Tissue-specific transcriptome and metabolome analyses reveal candidate genes for lignan biosynthesis in the medicinal plant Schisandra sphenanthera.

Schisandra sphenanthera is an extremely important medicinal plant, and its main medicinal component is bioactive lignans. The S. sphenanthera fruit is preferred by the majority of consumers, and the root, stem, and leaf are not fully used. To better understand the lignan metabolic pathway, transcriptome and metabolome analyses were performed on the four major tissues of S. sphenanthera. A total of 167,972,229 transcripts and 91,215,760 unigenes with an average length of 752 bp were identified. Tissue-specific gene analysis revealed that the root had the highest abundance of unique unigenes (9703), and the leaves had the lowest (189). Transcription factor analysis showed that MYB-, bHLH- and ERF-transcription factors, which played important roles in the regulation of secondary metabolism, showed rich expression patterns and may be involved in the regulation of processes involved in lignan metabolism. In different tissues, lignans were preferentially enriched in fruit and roots by gene expression profiles related to lignan metabolism and relative lignan compound content. Furthermore, schisandrin B is an important compound in S. sphenanthera. According to weighted gene co-expression network analysis, PAL1, C4H-2, CAD1, CYB8, OMT27, OMT57, MYB18, bHLH3, and bHLH5 can be related to the accumulation of lignans in S. sphenanthera fruit, CCR5, SDH4, CYP8, CYP20, and ERF7 can be related to the accumulation of lignans in S. sphenanthera roots. In this study, transcriptome sequencing and targeted metabolic analysis of lignans will lay a foundation for the further study of their biosynthetic genes.

Savinin Triggers Programmed Cell Death of Ray Parenchyma Cells in Heartwood Formation of Taiwania cryptomerioides Hayata.

The purpose of this study was to investigate the relationship between lignan biosynthesis and programmed cell death (PCD) of ray parenchyma cells during the heartwood formation of Taiwania (Taiwania cryptomerioides Hayata). Since the PCD of ray parenchyma cells and the synthesis of lignans are the two main processes involved in the formation of heartwood, both of which need to be completed through gene regulation. Based on the results of genomics and bioinformatics analysis, that the PCD of tracheids are induced by genotoxic, and the PCD of ray parenchyma cells is induced by biological factors, such as fungi, bacteria, and viruses, which could induce oxidative stress. According to the results of time-of-flight secondary ion mass spectrometry (ToF-SIMS) analysis, lignans are produced in ray parenchyma cells, and the accumulation of savinin and its downstream lignans might be the cause of PCD in ray parenchyma cells. An in vitro experiment further confirmed that the accumulation of savinin could cause protoplasts of Taiwania's xylem to produce taiwanin A, which is the marker of heartwood formation in Taiwania. Resulting in an increase in reactive oxygen species (ROS) content, which could induce oxidative stress in ray parenchyma cells and potentially lead to PCD. Based on these findings, we conclude that accumulation of savinin could be induced PCD of ray parenchyma cells in heartwood formation in Taiwania.

Identification and functional characterization of the dirigent gene family in Phryma leptostachya and the contribution of PlDIR1 in lignan biosynthesis.

BACKGROUND: Furofuran lignans, the main insecticidal ingredient in Phryma leptostachya, exhibit excellent controlling efficacy against a variety of pests. During the biosynthesis of furofuran lignans, Dirigent proteins (DIRs) are thought to be dominant in the stereoselective coupling of coniferyl alcohol to form ( ±)-pinoresinol. There are DIR family members in almost every vascular plant, but members of DIRs in P. leptostachya are unknown. To identify the PlDIR genes and elucidate their functions in lignan biosynthesis, this study performed transcriptome-wide analysis and characterized the catalytic activity of the PlDIR1 protein. RESULTS: Fifteen full-length unique PlDIR genes were identified in P. leptostachya. A phylogenetic analysis of the PlDIRs classified them into four subfamilies (DIR-a, DIR-b/d, DIR-e, and DIR-g), and 12 conserved motifs were found among them. In tissue-specific expression analysis, except for PlDIR7, which displayed the highest transcript abundance in seeds, the other PlDIRs showed preferential expression in roots, leaves, and stems. Furthermore, the treatments with signaling molecules demonstrated that PlDIRs could be significantly induced by methyl jasmonate (MeJA), salicylic acid (SA), and ethylene (ETH), both in the roots and leaves of P. leptostachya. In examining the tertiary structure of the protein and the critical amino acids, it was found that PlDIR1, one of the DIR-a subfamily members, might be involved in the region- and stereo-selectivity of the phenoxy radical. Accordingly, LC-MS/MS analysis demonstrated the catalytic activity of recombinant PlDIR1 protein from Escherichia coli to direct coniferyl alcohol coupling into ( +)-pinoresinol. The active sites and hydrogen bonds of the interaction between PlDIR1 and bis-quinone methide (bisQM), the intermediate in ( +)-pinoresinol formation, were analyzed by molecular docking. As a result, 18 active sites and 4 hydrogen bonds (Asp-42, Ala-113, Leu-138, Arg-143) were discovered in the PlDIR1-bisQM complex. Moreover, correlation analysis indicated that the expression profile of PlDIR1 was closely connected with lignan accumulations after SA treatment. CONCLUSIONS: The results of this study will provide useful clues for uncovering P. leptostachya's lignan biosynthesis pathway as well as facilitate further studies on the DIR family.

The wild allotetraploid sesame genome provides novel insights into evolution and lignan biosynthesis.

INTRODUCTION: The wild tetraploid sesame (Sesamum schinzianum), an ancestral relative of diploid cultivated sesame, grows in the tropical desert of the African Plateau. As a valuable seed resource, wild sesame has several advantageous traits, such as strong environmental adaptability and an extremely high content of sesamolin in its seeds. High-quality genome assembly is essential for a detailed understanding of genome structure, genome evolution and crop improvement. OBJECTIVES: Here, we generated two high-quality chromosome-scale genomes from S. schinzianum and a cultivated diploid elite sesame (Sesamum indicum L.) to investigate the potential genetic basis underlying these traits of wild sesame. METHODS: The long-read data from PacBio Sequel II platform and high-throughput chromosome conformation capture (Hi-C) data were used to construct high-quality sesame genome. Then dissecting the molecular mechanisms of sesame evolution and lignan biosynthesis through comparative genomics and transcriptomics. RESULTS: We found evidence of divergent evolution that involved differences in the number, sequence and expression level of homologous genes between the two sets of subgenomes from allotetraploids in S. schinzianum, all of which might be driven by subfunctionalization after polyploidization. Furthermore, it was found that a great number of genes involved in the stress response have undergone positive selection and resulted from gene family expansion in the wild sesame genome compared with the cultivated sesame genome, which, overall, is associated with adaptative evolution to the environment. We hypothesized that the sole functional member CYP92B14 (SscC22g35272) could be associated with high content of sesamolin in wild sesame seeds. CONCLUSION: This study provides high-quality wild allotetraploid sesame and cultivated sesame genomes, reveals evolutionary features of the allotetraploid genome and provides novel insights into lignan synthesis pathways. Meanwhile, the wild sesame genome will be an important resource to conduct comparative genomic and evolutionary studies and plant improvement programmes.

Genome-wide analysis of AP2/ERF superfamily in Isatis indigotica.

OBJECTIVE: AP2/ERF (APETALA2/ethylene-responsive factor) superfamily is one of the largest gene families in plants and has been reported to participate in various biological processes, such as the regulation of biosynthesis of active lignan. However, few studies have investigated the genome-wide role of the AP2/ERF superfamily in Isatis indigotica. This study establishes a complete picture of the AP2/ERF superfamily in I. indigotica and contributes valuable information for further functional characterization of IiAP2/ERF genes and supports further metabolic engineering. METHODS: To identify the IiAP2/ERF superfamily genes, the AP2/ERF sequences from Arabidopsis thaliana and Brassica rapa were used as query sequences in the basic local alignment search tool. Bioinformatic analyses were conducted to investigate the protein structure, motif composition, chromosome location, phylogenetic relationship, and interaction network of the IiAP2/ERF superfamily genes. The accuracy of omics data was verified by quantitative polymerase chain reaction and heatmap analyses. RESULTS: One hundred and twenty-six putative IiAP2/ERF genes in total were identified from the I. indigotica genome database in this study. By sequence alignment and phylogenetic analysis, the IiAP2/ERF genes were classified into 5 groups including AP2, ERF, DREB (dehydration-responsive element-binding factor), Soloist and RAV (related to abscisic acid insensitive 3/viviparous 1) subfamilies. Among which, 122 members were unevenly distributed across seven chromosomes. Sequence alignment showed that I. indigotica and A. thaliana had 30 pairs of orthologous genes, and we constructed their interaction network. The comprehensive analysis of gene expression pattern in different tissues suggested that these genes may play a significant role in organ growth and development of I. indigotica. Members that may regulate lignan biosynthesis in roots were also preliminarily identified. Ribonucleic acid sequencing analysis revealed that the expression of 76 IiAP2/ERF genes were up- or down-regulated under salt or drought treatment, among which, 33 IiAP2/ERF genes were regulated by both stresses. CONCLUSION: This study undertook a genome-wide characterization of the AP2/ERF superfamily in I. indigotica, providing valuable information for further functional characterization of IiAP2/ERF genes and discovery of genetic targets for metabolic engineering.

Genome-wide analysis of AP2/ERF superfamily in Isatis indigotica / 中西医结合学报

OBJECTIVE@#AP2/ERF (APETALA2/ethylene-responsive factor) superfamily is one of the largest gene families in plants and has been reported to participate in various biological processes, such as the regulation of biosynthesis of active lignan. However, few studies have investigated the genome-wide role of the AP2/ERF superfamily in Isatis indigotica. This study establishes a complete picture of the AP2/ERF superfamily in I. indigotica and contributes valuable information for further functional characterization of IiAP2/ERF genes and supports further metabolic engineering.@*METHODS@#To identify the IiAP2/ERF superfamily genes, the AP2/ERF sequences from Arabidopsis thaliana and Brassica rapa were used as query sequences in the basic local alignment search tool. Bioinformatic analyses were conducted to investigate the protein structure, motif composition, chromosome location, phylogenetic relationship, and interaction network of the IiAP2/ERF superfamily genes. The accuracy of omics data was verified by quantitative polymerase chain reaction and heatmap analyses.@*RESULTS@#One hundred and twenty-six putative IiAP2/ERF genes in total were identified from the I. indigotica genome database in this study. By sequence alignment and phylogenetic analysis, the IiAP2/ERF genes were classified into 5 groups including AP2, ERF, DREB (dehydration-responsive element-binding factor), Soloist and RAV (related to abscisic acid insensitive 3/viviparous 1) subfamilies. Among which, 122 members were unevenly distributed across seven chromosomes. Sequence alignment showed that I. indigotica and A. thaliana had 30 pairs of orthologous genes, and we constructed their interaction network. The comprehensive analysis of gene expression pattern in different tissues suggested that these genes may play a significant role in organ growth and development of I. indigotica. Members that may regulate lignan biosynthesis in roots were also preliminarily identified. Ribonucleic acid sequencing analysis revealed that the expression of 76 IiAP2/ERF genes were up- or down-regulated under salt or drought treatment, among which, 33 IiAP2/ERF genes were regulated by both stresses.@*CONCLUSION@#This study undertook a genome-wide characterization of the AP2/ERF superfamily in I. indigotica, providing valuable information for further functional characterization of IiAP2/ERF genes and discovery of genetic targets for metabolic engineering.

Characterization of Peroxidase and Laccase Gene Families and In Silico Identification of Potential Genes Involved in Upstream Steps of Lignan Formation in Sesame.

Peroxidases and laccases are oxidative enzymes involved in physiological processes in plants, covering responses to biotic and abiotic stress as well as biosynthesis of health-promoting specialized metabolites. Although they are thought to be involved in the biosynthesis of (+)-pinoresinol, a comprehensive investigation of this class of enzymes has not yet been conducted in the emerging oil crop sesame and no information is available regarding the potential (+)-pinoresinol synthase genes in this crop. In the present study, we conducted a pan-genome-wide identification of peroxidase and laccase genes coupled with transcriptome profiling of diverse sesame varieties. A total of 83 and 48 genes have been identified as coding for sesame peroxidase and laccase genes, respectively. Based on their protein domain and Arabidopsis thaliana genes used as baits, the genes were classified into nine and seven groups of peroxidase and laccase genes, respectively. The expression of the genes was evaluated using dynamic transcriptome sequencing data from six sesame varieties, including one elite cultivar, white vs black seed varieties, and high vs low oil content varieties. Two peroxidase genes (SiPOD52 and SiPOD63) and two laccase genes (SiLAC1 and SiLAC39), well conserved within the sesame pan-genome and exhibiting consistent expression patterns within sesame varieties matching the kinetic of (+)-pinoresinol accumulation in seeds, were identified as potential (+)-pinoresinol synthase genes. Cis-acting elements of the candidate genes revealed their potential involvement in development, hormonal signaling, and response to light and other abiotic triggers. Transcription factor enrichment analysis of promoter regions showed the predominance of MYB binding sequences. The findings from this study pave the way for lignans-oriented engineering of sesame with wide potential applications in food, health and medicinal domains.

Functional Characterization of UDP-Glycosyltransferases Involved in Anti-viral Lignan Glycosides Biosynthesis in Isatis indigotica.

Isatis indigotica is a popular herbal medicine with its noticeable antiviral properties, which are primarily due to its lignan glycosides such as lariciresinol-4-O-ß-D-glucoside and lariciresinol-4,4'-bis-O-ß-D-glucosides (also called clemastanin B). UDP-glucose-dependent glycosyltransferases are the key enzymes involved in the biosynthesis of these antiviral metabolites. In this study, we systematically characterized the UGT72 family gene IiUGT1 and two UGT71B family genes, IiUGT4 and IiUGT71B5a, with similar enzymatic functions. Kinetic analysis showed that IiUGT4 was more efficient than IiUGT1 or IiUGT71B5a for the glycosylation of lariciresinol. Further knock-down and overexpression of these IiUGTs in I. indigotica's hairy roots indicates that they play different roles in planta: IiUGT71B5a primarily participates in the biosynthesis of coniferin not pinoresinol diglucoside, and IiUGT1 primarily participates in the biosynthesis of pinoresinol diglucoside, while IiUGT4 is responsible for the glycosylation of lariciresinol and plays a dominant role in the biosynthesis of lariciresinol glycosides in I. indigotica. Analysis of the molecular docking and site-mutagenesis of IiUGT4 have found that key residues for its catalytic activity are H373, W376, E397, and that F151 could be associated with substrate preference. This study elucidates the biosynthetic route of anti-viral lignan glycosides in I. indigotica, and provides the foundation for the production of anti-viral lignan glycosides via synthetic biology under the heterologous model.

Transcriptome analyses reveal the expression profile of genes related to lignan biosynthesis in Anthriscus sylvestris L. Hoffm. Gen.

Anthriscus sylvestris L. Hoffm. Gen (A. sylvestris) is a perennial herb widely used for antitussive and diuretic purposes in traditional Korean and Chinese medicine. Lignans are critical secondary metabolites with widely pharmacological activities in A. sylvestris. Using transcriptome data of A. sylvestris, we identified genes related to lignan biosynthesis. In all, 123,852 unigenes were obtained from the flowers, leaves, roots, and stems of A. sylvestris with the Illumina HiSeq 4000 platform. The average length of unigenes was 1,123 bp and 91,217 (73.65%) of them were annotated in public databases. Differentially expressed genes and root-specific genes were analyzed between roots and the other three tissue types by comparing gene expression profiles. Specifically, the key enzyme genes involved in lignan biosynthesis were identified and analyzed. The expression levels of some of these genes were highest in the roots, consistent with the accumulation of deoxypodophyllotoxin. These expression levels were experimentally verified via quantitative real-time polymerase chain reaction (qRT-PCR). This research provides valuable information on the transcriptome data of A. sylvestris and the identification of candidate genes associated with the biosynthesis of lignans, laying the foundation for further research on genomics in A. sylvestris and related species. Supplementary Information: The online version contains supplementary material available at 10.1007/s12298-022-01156-w.

Long-read transcriptome sequencing provides insight into lignan biosynthesis during fruit development in Schisandra chinensis.

BACKGROUND: Schisandra chinensis, an ancient member of the most basal angiosperm lineage which is known as the ANITA, is a fruit-bearing vine with the pharmacological effects of a multidrug system, such as antioxidant, anti-inflammatory, cardioprotective, neuroprotective, anti-osteoporosis effects. Its major bioactive compound is represented by lignans such as schisandrin. Molecular characterization of lignan biosynthesis in S. chinensis is of great importance for improving the production of this class of active compound. However, the biosynthetic mechanism of schisandrin remains largely unknown. RESULTS: To understand the potential key catalytic steps and their regulation of schisandrin biosynthesis, we generated genome-wide transcriptome data from three different tissues of S. chinensis cultivar Cheongsoon, including leaf, root, and fruit, via long- and short-read sequencing technologies. A total of 132,856 assembled transcripts were generated with an average length of 1.9 kb and high assembly completeness. Overall, our data presented effective, accurate gene annotation in the prediction of functional pathways. In particular, the annotation revealed the abundance of transcripts related to phenylpropanoid biosynthesis. Remarkably, transcriptome profiling during fruit development of S. chinensis cultivar Cheongsoon revealed that the phenylpropanoid biosynthetic pathway, specific to coniferyl alcohol biosynthesis, showed a tendency to be upregulated at the postfruit development stage. Further the analysis also revealed that the pathway forms a transcriptional network with fruit ripening-related genes, especially the ABA signaling-related pathway. Finally, candidate unigenes homologous to isoeugenol synthase 1 (IGS1) and dirigent-like protein (DIR), which are subsequently activated by phenylpropanoid biosynthesis and thus catalyze key upstream steps in schisandrin biosynthesis, were identified. Their expression was increased at the postfruit development stage, suggesting that they may be involved in the regulation of schisandrin biosynthesis in S. chinensis. CONCLUSIONS: Our results provide new insights into the production and accumulation of schisandrin in S. chinensis berries and will be utilized as a valuable transcriptomic resource for improving the schisandrin content.

Integrative analyses of targeted metabolome and transcriptome of Isatidis Radix autotetraploids highlighted key polyploidization-responsive regulators.

BACKGROUND: Isatidis Radix, the root of Isatis indigotica Fort. (Chinese woad) can produce a variety of efficacious compound with medicinal properties. The tetraploid I. indigotica plants exhibit superior phenotypic traits, such as greater yield, higher bioactive compounds accumulation and enhanced stress tolerance. In this study, a comparative transcriptomic and metabolomic study on Isatidis Radix autotetraploid and its progenitor was performed. RESULTS: Through the targeted metabolic profiling, 283 metabolites were identified in Isatidis Radix, and 70 polyploidization-altered metabolites were obtained. Moreover, the production of lignans was significantly increased post polyploidization, which implied that polyploidization-modulated changes in lignan biosynthesis. Regarding the transcriptomic shift, 2065 differentially expressed genes (DEGs) were identified as being polyploidy-responsive genes, and the polyploidization-altered DEGs were enriched in phenylpropanoid biosynthesis and plant hormone signal transduction. The further integrative analysis of polyploidy-responsive metabolome and transcriptome showed that 1584 DEGs were highly correlated with the 70 polyploidization-altered metabolites, and the transcriptional factors TFs-lignans network highlighted 10 polyploidy-altered TFs and 17 fluctuated phenylpropanoid pathway compounds. CONCLUSIONS: These results collectively indicated that polyploidization contributed to the high content of active compounds in autotetraploid roots, and the gene-lignan pathway network analysis highlighted polyploidy-responsive key functional genes and regulators.

Lignans of Sesame (Sesamum indicum L.): A Comprehensive Review.

Major lignans of sesame sesamin and sesamolin are benzodioxol--substituted furofurans. Sesamol, sesaminol, its epimers, and episesamin are transformation products found in processed products. Synthetic routes to all lignans are known but only sesamol is synthesized industrially. Biosynthesis of furofuran lignans begins with the dimerization of coniferyl alcohol, followed by the formation of dioxoles, oxidation, and glycosylation. Most genes of the lignan pathway in sesame have been identified but the inheritance of lignan content is poorly understood. Health-promoting properties make lignans attractive components of functional food. Lignans enhance the efficiency of insecticides and possess antifeedant activity, but their biological function in plants remains hypothetical. In this work, extensive literature including historical texts is reviewed, controversial issues are critically examined, and errors perpetuated in literature are corrected. The following aspects are covered: chemical properties and transformations of lignans; analysis, purification, and total synthesis; occurrence in Seseamum indicum and related plants; biosynthesis and genetics; biological activities; health-promoting properties; and biological functions. Finally, the improvement of lignan content in sesame seeds by breeding and biotechnology and the potential of hairy roots for manufacturing lignans in vitro are outlined.

IiWRKY34 positively regulates yield, lignan biosynthesis and stress tolerance in Isatis indigotica.

Yield potential, pharmaceutical compounds production and stress tolerance capacity are 3 classes of traits that determine the quality of medicinal plants. The autotetraploid Isatis indigotica has greater yield, higher bioactive lignan accumulation and enhanced stress tolerance compared with its diploid progenitor. Here we show that the transcription factor IiWRKY34, with higher expression levels in tetraploid than in diploid I. indigotica, has large pleiotropic effects on an array of traits, including biomass growth rates, lignan biosynthesis, as well as salt and drought stress tolerance. Integrated analysis of transcriptome and metabolome profiling demonstrated that IiWRKY34 expression had far-reaching consequences on both primary and secondary metabolism, reprograming carbon flux towards phenylpropanoids, such as lignans and flavonoids. Transcript-metabolite correlation analysis was applied to construct the regulatory network of IiWRKY34 for lignan biosynthesis. One candidate target Ii4CL3, a key rate-limiting enzyme of lignan biosynthesis as indicated in our previous study, has been demonstrated to indeed be activated by IiWRKY34. Collectively, the results indicate that the differentially expressed IiWRKY34 has contributed significantly to the polyploidy vigor of I. indigotica, and manipulation of this gene will facilitate comprehensive improvements of I. indigotica herb.

Dynamic metabolic and transcriptomic profiling of methyl jasmonate-treated hairy roots reveals synthetic characters and regulators of lignan biosynthesis in Isatis indigotica Fort.

A molecular description of lignan biosynthesis in Isatis indigotica displaying its synthetic characteristics and regulatory mechanism is of great importance for the improvement of the production of this class of active compounds. To discover the potential key catalytic steps and regulatory genes, I. indigotica hairy roots elicited by methyl jasmonate (MeJA) were used as a source of systematic variation for exploring the metabolic/transcriptional changes and candidate genes that might play key roles in lignan biosynthesis. The reprogramming modulated by MeJA was classified into three distinct phases, referred to as signal responding, transcriptional activation of metabolic pathways and accumulation of metabolites. Candidate genes were pooled according to the three phases and applied to co-expression network analysis. In total, 17 genes were identified as hub genes. 4CL3 was selected to validate its impact on lignan biosynthesis. RNAi of 4CL3 resulted in a significant reduction in lignan production. Taken together with its catalytic property, a major route of lignan biosynthesis in I. indigotica was highlighted, which was catalysed by 4CL3 via the esterization of caffeic acid. In conclusion, this study provides new insights into lignan biosynthesis as well as potential targets for metabolic engineering in I. indigotica.