Isolation and characterization of a gene encoding cinnamoyl-CoA reductase from Isatis indigotica Fort

Progress based on a multi-omics research strategy in the biosynthesis and modernization of active ingredients of Herpetospermum pedunculosum seeds

Herpetospermum pedunculosum seeds also known as Herpetospermum caudigerum Wall. is the mature seed of the Herpetospermum pedunculosum(Ser.) C. B. Clarke,Cucurbitaceae. Modern pharmacological studies have shown that H. pedunculosum has hepatoprotective, anti-inflammatory, anti-gout and antibacterial pharmacological activities. The biologically active chemical components include lignin compounds such as Herpetin, Herpetetrone, Herpetoriol and so on. The natural product displays considerable skeletal diversity and structural complexity, offering significant opportunities for novel drug discovery. Based on the multi-omics research strategy and the 'gene-protein-metabolite' research framework, the biosynthetic pathway of terpenoids and lignans in H. pedunculosum has has been elucidated at multiple levels. These approaches provide comprehensive genetic information for cloning and identification of pertinent enzyme genes. Furthermore, the application of multi-omics integrative approaches provides a scientific means to elucidate entire secondary metabolic pathways. We investigated the biosynthetic pathways of lignin and terpene components in H. pedunculosum and conducted bioinformatics analysis of the crucial enzyme genes involved in the biosynthetic process using genomic and transcriptomic data. We identified candidate genes for six key enzymes in the biosynthetic pathway. This review reports on the current literature on pharmacological investigations of H. pedunculosum, proposing its potential as an antidiabetic agent. Moreover, we conclude, for the first time, the identification of key enzyme genes potentially involved in the biosynthesis of active compounds in H. pedunculosum. This review provides a scientific foundation for the discovery of novel therapeutic agents from natural sources.

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.

Functional analysis of a dirigent protein AtsDIR23 in Acorustatarinowii

Acorus tatarinowii (A. tatarinowii) is a medicinal plant of the Araceae family. Currently, pharmacology focuses on the study of volatile oils, but there are few reports of another important secondary metabolite, lignan. Dirigent protein is thought to play an important role in plant secondary metabolism and responds to a variety of biotic and abiotic stresses. However, the DIR gene family of A. tatarinowii has not been systematically analyzed, and it is unknown whether it affects lignan synthesis. In this study, a total of 27 AtsDIRs were identified by comprehensive analysis of the genome of the medicinal plant A. tatarinowii, and the candidate gene AtsDIR23 that may be involved in lignan synthesis was screened through bioinformatics and transcriptome analysis. It is worth noting that AtsDIR23 is significantly expressed in rhizomes and is a member of the DIR-a subfamily. Subsequently, subcellular localization revealed that AtsDIR23 was localized in chloroplasts. The functional verification of AtsDIR23 b y the transient transformation of A. tatarinowii and the stable transformation of Arabidopsis thaliana showed that the content of lignans in overexpressed plants increased. Co-expression analysis screening revealed the MYB transcription factor (AtsMYB91) that is highly correlated with AtsDIR23 expression, while yeast one-hybrid assays and double luciferase experiments showed that AtsMYB91 negatively regulated the expression of AtsDIR23 b y binding to the AtsDIR23 promoter. In conclusion, AtsDIR23 can promote the accumulation of lignans, which provides a reference for further research on the regulation of lignans by DIR genes.

Molecular cloning and characterization of Cinnamoyl-CoA reductase promoter gene from Asarum sieboldii Miq

Asarum sieboldii Miq., a perennial herb of the family Aristolochiaceae, is widely used in China to treat cold, fever, aphthous stomatitis, toothache, gingivitis, and rheumatoid arthritis. Methyleugenol is the most representative pharmacological constituent of this medicinal herb. Cinnamoyl-CoA reductase (CCR), which has been well known for occupying a critical position in the lignin biosynthesis pathway, is also shared with the biosynthesis of methyleugenol. To better understand the regulatory mechanisms of methyleugenol biosynthesis, a 1530-bp long promoter region of the AsCCR1 gene was isolated. PLACE and PlantCARE analysis affirmed the existence of the core promoter elements such as TATA and CAAT boxes, abiotic stress-responsive cis-regulation elements like abscisic acid-responsive element, G-box, and MBS in the isolated sequence. The histochemical assay suggested that it was a constitutive promoter, highly expressed in the root tissue. Moreover, the region of -200 bp to ATG (start codon) was enough to drive the expression of It GUS gene. Treatments with low temperature and high concentration of gibberellin or abscisic acid demonstrated the abiotic stress-induced expression of the AsCCR1 promoter. Overall, this study revealed the isolation and characterization of the AsCCR1 promoter. Moreover, it also provided a candidate gene for molecular breeding in A. sieboldii to enhance its pharmacological potential.

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.

SEPALLATA--like genes of Isatis indigotica can affect the architecture of the inflorescences and the development of the floral organs

BACKGROUND

The architecture of inflorescence and the development of floral organs can influence the yield of seeds and have a significant impact on plant propagation. E-class floral homeotic MADS-box genes exhibit important roles in regulation of floral transition and differentiation of floral organs. Woad (Isatis indigotica) possesses unique inflorescence, floral organs and fruit. However, very little research has been carried out to determine the function of MADS-box genes in this medicinal cruciferous plant species.

RESULTS

SEPALLATA orthologs in I. indigotica were cloned by degenerate PCR. The sequence possessing the highest identity with SEP2 and SEP4 of Arabidopsis were named as IiSEP2 and IiSEP4, respectively. Constitutive expression of IiSEP2 in Columbia (Col-0) ecotype of Arabidopsis led to early flowering, and the number of the flowers and the lateral branches was reduced, indicating an alteration in architecture of the inflorescences. Moreover, the number of the floral organs was declined, the sepals were turned into carpelloid tissues bearing stigmatic papillae and ovules, and secondary flower could be produced in apetalous terminal flowers. In 35S::IiSEP4-GFP transgenic Arabidopsis plants in Landsberg erecta (Ler) genetic background, the number of the floral organs was decreased, sepals were converted into curly carpelloid structures, accompanied by generation of ovules. Simultaneously, the size of petals, stamens and siliques was diminished. In 35S::IiSEP4-GFP transgenic plants of apetalous ap1 cal double mutant in Ler genetic background, the cauliflower phenotype was attenuated significantly, and the petal formation could be rescued. Occasionally, chimeric organs composed of petaloid and sepaloid tissues, or petaloid and stamineous tissues, were produced in IiSEP4 transgenic plants of apl cal double mutant. It suggested that overexpression of IiSEP4 could restore the capacity in petal differentiation. Silencing of IiSEP4 by Virus-Induced Gene Silencing (VIGS) can delay the flowering time, and reduce the number and size of the floral organs in woad flowers.

CONCLUSION

All the results showed that SEPALLATA-like genes could influence the architecture of the inflorescence and the determinacy of the floral meristems, and was also related to development of the floral organs.

Isatis indigotica: from (ethno) botany, biochemistry to synthetic biology

Isatis indigotica Fort. (Chinese woad) is a species with an ancient and well-documented history as an indigo dye and medicinal plant. It is often confused with Isatis tinctoria L. (European woad), a medicinal plant in Europe. Here, the differences between I. indigotica and I. tinctoria are systematically described. The usage development history, clinical applications and pharmacological activities, and chemical components of I. indigotica are also summarized. Lignans, indole alkaloids, and their corresponding derivatives have been identified as the major active ingredients of I. indigotica and are associated with anti-viral, anti-inflammatory, anti-cancer, and other health-promoting activities. Notable progress has been made in understanding the biosynthetic pathway and regulation mechanism of lignans and indole alkaloids in I. indigotica, the results from which should facilitate the process of targeted metabolic engineering or synthetic biology. Moreover, multiple biotechnology methods such as polyploid breeding and genetic engineering have been used with I. indigotica to result in, for example, greater yields, higher levels of bioactive component accumulation, and enhanced stress tolerance to salt, drought, and insects. Some issues require additional analyses, and suggestions for future research on I. indigotica are also discussed.

Overexpression of Cinnamoyl-CoA Reductase 2 in Brassica napus Increases Resistance to Sclerotinia sclerotiorum by Affecting Lignin Biosynthesis

Sclerotinia sclerotiorum causes severe yield and economic losses for many crop and vegetable species, especially Brassica napus. To date, no immune B. napus germplasm has been identified, giving rise to a major challenge in the breeding of Sclerotinia resistance. In the present study, we found that, compared with a Sclerotinia-susceptible line (J902), a Sclerotinia-resistant line (J964) exhibited better xylem development and a higher lignin content in the stems, which may limit the invasion and spread of S. sclerotiorum during the early infection period. In addition, genes involved in lignin biosynthesis were induced under S. sclerotiorum infection in both lines, indicating that lignin was deposited proactively in infected tissues. We then overexpressed BnaC.CCR2.b, which encodes the first rate-limiting enzyme (cinnamoyl-CoA reductase) that catalyzes the reaction of lignin-specific pathways, and found that overexpression of BnaC.CCR2.b increased the lignin content in the stems of B. napus by 2.28-2.76% under normal growth conditions. We further evaluated the Sclerotinia resistance of BnaC.CCR2.b overexpression lines at the flower-termination stage and found that the disease lesions on the stems of plants in the T₂ and T₃ generations decreased by 12.2-33.7% and 32.5-37.3% compared to non-transgenic control plants, respectively, at 7days post-inoculation (dpi). The above results indicate that overexpression of BnaC.CCR2.b leads to an increase in lignin content in the stems, which subsequently leads to increased resistance to S. sclerotiorum. Our findings demonstrate that increasing the lignin content in the stems of B. napus is an important strategy for controlling Sclerotinia.

Identification and expression analysis of the WRKY gene family in Isatis indigotica

The WRKY proteins, which represent one of the largest families of transcriptional regulators in plants, play pivotal roles in regulating multiple processes of growth and development, particularly in diverse stress responses. Isatis indigotica is widely used in Traditional Chinese Medicine and is famous for its use as a dye for the color indigo. However, reports of the WRKY gene family in I. indigotica are limited. In this study, 64 IiWRKY genes encoding proteins with the complete WRKY domain were identified from genome of I. indigotica. Based on their structure and phylogenetic relationships of this gene family in I. indigotica, the IiWRKY genes were classified into three groups: Group I (n = 13), Group II (n = 35) and Group III (n = 16). Sequence alignment revealed that IiWRKY proteins harbored two variants, WRKYRQK and WRKYGKK, of the highly conserved WRKYGQK motif. The number of exons in IiWRKY genes varied from two to 14, with most of IiWRKY genes containing three exons. Investigation of gene duplication demonstrated that 10 and 14 IiWRKY genes were incorporated in tandem and segmental duplication events, respectively. Finally, the expression profiles derived from transcriptome data and quantitative real-time PCR analysis showed distinct expression patterns of these IiWRKY gene in five different organs or in response to four abiotic stresses. Taken together, our results will contribute to functional analysis of IiWRKY genes, and also provide a basis for further clarification of the molecular mechanism of stress responses in this important herb.

A chromosome-scale genome assembly of Isatis indigotica, an important medicinal plant used in traditional Chinese medicine: An Isatis genome

Isatis indigotica (2n = 14) is an important medicinal plant in China. Its dried leaves and roots (called Isatidis Folium and Isatidis Radix, respectively) are broadly used in traditional Chinese medicine for curing diseases caused by bacteria and viruses such as influenza and viral pneumonia. Various classes of compounds isolated from this species have been identified as effective ingredients. Previous studies based on transcriptomes revealed only a few candidate genes for the biosynthesis of these active compounds in this medicinal plant. Here, we report a high-quality chromosome-scale genome assembly of I. indigotica with a total size of 293.88 Mb and scaffold N50 = 36.16 Mb using single-molecule real-time long reads and high-throughput chromosome conformation capture techniques. We annotated 30,323 high-confidence protein-coding genes. Based on homolog searching and functional annotations, we identified many candidate genes involved in the biosynthesis of main active components such as indoles, terpenoids, and phenylpropanoids. In addition, we found that some key enzyme-coding gene families related to the biosynthesis of these components were expanded due to tandem duplications, which likely drove the production of these major active compounds and explained why I. indigotica has excellent antibacterial and antiviral activities. Our results highlighted the importance of genome sequencing in identifying candidate genes for metabolite synthesis in medicinal plants.

Selection of reference genes for the quantitative real-time PCR normalization of gene expression in Isatis indigotica fortune

Background

Isatis indigotica, a traditional Chinese medicine, produces a variety of active ingredients. However, little is known about the key genes and corresponding expression profiling involved in the biosynthesis pathways of these ingredients. Quantitative real-time polymerase chain reaction (qRT-PCR) is a powerful, commonly-used method for gene expression analysis, but the accuracy of the quantitative data produced depends on the appropriate selection of reference genes.

Results

In this study, the systematic analysis of the reference genes was performed for quantitative real-Time PCR normalization in I. indigotica. We selected nine candidate reference genes, including six traditional housekeeping genes (ACT, α-TUB, β-TUB, UBC, CYP, and EF1-α), and three newly stable internal control genes (MUB, TIP41, and RPL) from a transcriptome dataset of I. indigotica, and evaluated their expression stabilities in different tissues (root, stem, leaf, and petiole) and leaves exposed to three abiotic treatments (low-nitrogen, ABA, and MeJA) using geNorm, NormFinder, BestKeeper, and comprehensive RefFind algorithms. The results demonstrated that MUB and EF1-α were the two most stable reference genes for all samples. TIP41 as the optimal reference gene for low-nitrogen stress and MeJA treatment, while ACT had the highest ranking for ABA treatment and CYP was the most suitable for different tissues.

Conclusions

The results revealed that the selection and validation of appropriate reference genes for normalizing data is mandatory to acquire accurate quantification results. The necessity of specific internal control for specific conditions was also emphasized. Furthermore, this work will provide valuable information to enhance further research in gene function and molecular biology on I. indigotica and other related species.

Electronic supplementary material

The online version of this article (10.1186/s12867-019-0126-y) contains supplementary material, which is available to authorized users.

Ectopic expression of IiFUL isolated from Isatis indigotica could change the reproductive growth of Arabidopsis thaliana

The coding sequence of IiFUL in Isatis indigotica was isolated and was used in transformation of Arabidopsis. IiFUL overexpressing Arabidopsis plants exhibited early flowering phenotype, accompanied with the reduction of flower number and the production of terminal flower on the top of the main stems. In development process, the flowers located on the top of the main stems generated a lot of variations in phenotype, including abnormal swelling of pistil, withering and numerical change of stamens and petals, appearance of stigmatoid tissues and naked ovules at the margin or inside of sepals. Besides, secondary flower could be formed within the flowers on the top of the main stems. These observations illustrated that IiFUL mainly affected the development of inflorescence meristems and pistils, but its ectopic expression could also disturb the normal growth of other floral organs. Moreover, the fertile siliques produced by the lateral inflorescences of IiFUL overexpressing Arabidopsis plants showed indehiscent phenotype, and the shape of the cauline leaves was changed significantly. The results of quantitative real-time PCR revealed that higher transcriptional levels of IiFUL could be detected in flowers and silicles of I. indigotica. In comprehensive consideration of the previous reports about the dehiscence phenotype of Arabidopsis siliques and the fact that the siliques of IiFUL overexpressing Arabidopsis plants were indehiscent in the present work, it can be speculated that high expression of IiFUL in pericarp is likely the reason why the silicles of I. indigotica possess an indehiscent phenotype.

Biochemical and Expression Analyses of the Rice Cinnamoyl-CoA Reductase Gene Family

Cinnamoyl-CoA reductase (CCR) is the first committed enzyme in the monolignol pathway for lignin biosynthesis and catalyzes the conversion of hydroxycinnamoyl-CoAs into hydroxycinnamaldehydes. In the rice genome, 33 genes are annotated as CCR and CCR-like genes, collectively called OsCCRs. To elucidate the functions of OsCCRs, their phylogenetic relationships, expression patterns at the transcription levels and biochemical characteristics were thoroughly analyzed. Of the 33 OsCCRs, 24 of them encoded polypeptides of lengths similar to those of previously identified plant CCRs. The other nine OsCCRs had much shorter peptide lengths. Phylogenetic tree and sequence similarities suggested OsCCR4, 5, 17, 18, 19, 20, and 21 as likely candidates for functional CCRs in rice. To elucidate biochemical functions, OsCCR1, 5, 17, 19, 20, 21, and 26 were heterologously expressed in Escherichia coli and the resulting recombinant OsCCRs were purified to apparent homogeneity. Activity assays of the recombinant OsCCRs with hydroxycinnamoyl-CoAs revealed that OsCCR17, 19, 20, and 21 were biochemically active CCRs, in which the NAD(P)-binding and NADP-specificity motifs as well as the CCR signature motif were fully conserved. The kinetic parameters of enzyme reactions revealed that feruloyl-CoA, a precursor for the guaiacyl (G)-unit of lignin, is the most preferred substrate of OsCCR20 and 21. This result is consistent with a high content (about 70%) of G-units in rice lignins. Phylogenetic analysis revealed that OsCCR19 and 20 were grouped with other plant CCRs involved in developmental lignification, whereas OsCCR17 and 21 were closely related to stress-responsible CCRs identified from other plant species. In agreement with the phylogenetic analysis, expression analysis demonstrated that OsCCR20 was constitutively expressed throughout the developmental stages of rice, showing particularly high expression levels in actively lignifying tissues, such as roots and stems. These results suggest that OsCCR20 is primarily involved in developmental deposition of lignins in secondary cell walls. As expected, the expressions of OsCCR17 and 21 were induced in response to biotic and abiotic stresses, such as Magnaporthe grisea and Xanthomonas oryzae pv. oryzae (Xoo) infections, UV-irradiation and high salinity, suggesting that these genes play a role in defense-related processes in rice.

AP2/ERF Transcription Factor, Ii049, Positively Regulates Lignan Biosynthesis in Isatis indigotica through Activating Salicylic Acid Signaling and Lignan/Lignin Pathway Genes

Lignans, such as lariciresinol and its derivatives, have been identified as effective antiviral ingredients in Isatis indigotica. Evidence suggests that the APETALA2/ethylene response factor (AP2/ERF) family might be related to the biosynthesis of lignans in I. indigotica. However, the special role played by the AP2/ERF family in the metabolism and its underlying putative mechanism still need to be elucidated. One novel AP2/ERF gene, named Ii049, was isolated and characterized from I. indigotica in this study. The quantitative real-time PCR analysis revealed that Ii049 was expressed highest in the root and responded to methyl jasmonate, salicylic acid (SA) and abscisic acid treatments to various degrees. Subcellular localization analysis indicated that Ii049 protein was localized in the nucleus. Knocking-down the expression of Ii049 caused a remarkable reduction of lignan/lignin contents and transcript levels of genes involved in the lignan/lignin biosynthetic pathway. Ii049 bound to the coupled element 1, RAV1AAT and CRTAREHVCBF2 motifs of genes IiPAL and IiCCR, the key structural genes in the lignan/lignin pathway. Furthermore, Ii049 was also essential for SA biosynthesis, and SA induced lignan accumulation in I. indigotica. Notably, the transgenic I. indigotica hairy roots overexpressing Ii049 showed high expression levels of lignan/lignin biosynthetic genes and SA content, resulting in significant accumulation of lignan/lignin. The best-engineered line (OVX049-10) produced 425.60 μg·g^-1 lariciresinol, an 8.3-fold increase compared with the wild type production. This study revealed the function of Ii049 in regulating lignan/lignin biosynthesis, which had the potential to increase the content of valuable lignan/lignin in economically significant medicinal plants.

Molecular Cloning, Characterization, and Expression Analysis of Lignin Genes from Sugarcane Genotypes Varying in Lignin Content

Sugarcane (Saccharum spp.) is one of the highest biomass-producing plant and the best lignocellulosic feedstock for ethanol production. To achieve more efficient conversion of biomass to ethanol, a better understanding of the main factors affecting biomass recalcitrance is needed. Therefore, with this objective, here, we report a systematic study on lignin content, deposition, identification, and cloning of genes involved in lignin biosynthesis and their differential expression in five sugarcane clones, EC11003, EC11010, IK 76-91, IK 76-99, and Co 86032. Lignin content among the clones varied from 26.87 to 23.19 % with the highest in the clone EC11010 and the lowest in high sugar Co86032. Lignin deposition studied through phloroglucinol staining of the cell walls implied that the sclerenchyma cells of the energy canes (EC11010 and EC11003) have more lignin deposition followed by the Erianthus (IK 76-91 and IK 76-99) clones whereas Co86032 has the minimum amount of lignin deposition. We cloned partial coding regions of important genes of lignification COMT (650 bp), CCR (332 bp), and PAL (650 bp) from Erianthus, wild relative of sugarcane followed by the expression analysis through real-time PCR. Differential expression analysis showed high level of expression for the three genes in the energy cane EC11010.

Characterization of developmental- and stress-mediated expression of cinnamoyl-CoA reductase in kenaf (Hibiscus cannabinus L.)

Cinnamoyl-CoA reductase (CCR) is an important enzyme for lignin biosynthesis as it catalyzes the first specific committed step in monolignol biosynthesis. We have cloned a full length coding sequence of CCR from kenaf (Hibiscus cannabinus L.), which contains a 1,020-bp open reading frame (ORF), encoding 339 amino acids of 37.37 kDa, with an isoelectric point (pI) of 6.27 (JX524276, HcCCR2). BLAST result found that it has high homology with other plant CCR orthologs. Multiple alignment with other plant CCR sequences showed that it contains two highly conserved motifs: NAD(P) binding domain (VTGAGGFIASWMVKLLLEKGY) at N-terminal and probable catalytic domain (NWYCYGK). According to phylogenetic analysis, it was closely related to CCR sequences of Gossypium hirsutum (ACQ59094) and Populus trichocarpa (CAC07424). HcCCR2 showed ubiquitous expression in various kenaf tissues and the highest expression was detected in mature flower. HcCCR2 was expressed differentially in response to various stresses, and the highest expression was observed by drought and NaCl treatments.

Biosynthesis of the active compounds of Isatis indigotica based on transcriptome sequencing and metabolites profiling

Backgroud

Isatis indigotica is a widely used herb for the clinical treatment of colds, fever, and influenza in Traditional Chinese Medicine (TCM). Various structural classes of compounds have been identified as effective ingredients. However, little is known at genetics level about these active metabolites. In the present study, we performed de novo transcriptome sequencing for the first time to produce a comprehensive dataset of I. indigotica.

Results

A database of 36,367 unigenes (average length = 1,115.67 bases) was generated by performing transcriptome sequencing. Based on the gene annotation of the transcriptome, 104 unigenes were identified covering most of the catalytic steps in the general biosynthetic pathways of indole, terpenoid, and phenylpropanoid. Subsequently, the organ-specific expression patterns of the genes involved in these pathways, and their responses to methyl jasmonate (MeJA) induction, were investigated. Metabolites profile of effective phenylpropanoid showed accumulation pattern of secondary metabolites were mostly correlated with the transcription of their biosynthetic genes. According to the analysis of UDP-dependent glycosyltransferases (UGT) family, several flavonoids were indicated to exist in I. indigotica and further identified by metabolic profile using UPLC/Q-TOF. Moreover, applying transcriptome co-expression analysis, nine new, putative UGTs were suggested as flavonol glycosyltransferases and lignan glycosyltransferases.

Conclusions

This database provides a pool of candidate genes involved in biosynthesis of effective metabolites in I. indigotica. Furthermore, the comprehensive analysis and characterization of the significant pathways are expected to give a better insight regarding the diversity of chemical composition, synthetic characteristics, and the regulatory mechanism which operate in this medical herb.