Next-generation sequencing (NGS) transcriptomes reveal association of multiple genes and pathways contributing to secondary metabolites accumulation in tuberous roots of Aconitum heterophyllum Wall

The known, unknown, and the intriguing about members of a critically endangered traditional medicinal plant genus Aconitum

Humanity will always be indebted to plants. In the ongoing scientific era, the 'Herbal Revolution' has helped discover several valuable medicinal plants and associated novel secondary metabolites from the diverse unexplored ecosystems, treating several diseases via phytotherapy. The Aconitum genus comprises several economically-important poisonous mountainous medicinal plant species whose unique biodiversity is on the verge of extinction due to illegal human intervention triggered habitat loss, over-harvesting, and unrestricted trading. Owing to its vast diversity of diterpene alkaloids, most species are extensively used to treat several ailments in rural parts of the world. Irrespective of this, many unexplored and intriguing prospects exist to understand and utilize this critical plant for human benefit. This systematic review tries to fill this gap by compiling information from the sporadically available literature known for ~300 Aconitum spp. regarding its nomenclature and classification, endangerment, plant morphology, ploidy, secondary metabolites, drug pharmacokinetics, conservation, and omics-based computational studies. We also depicted the disparity in the studied model organisms for this diverse genus. The absence of genomic/metagenomic data is becoming a limiting factor in understanding its plant physiology, metabolic pathways, and plant-microbes interactions, and therefore must be promoted. Additionally, government support and public participation are crucial in establishing conservation protocols to save this plant from endangerment.

A review on efforts for improvement in medicinally important chemical constituents inAconitum through biotechnological interventions

The genus Aconitum belongs to the family Ranunculaceae, is endowed with more than 350 species on the earth. Medicinally important aconitine type of diterpenoid alkaloids are the characteristic compounds in most of the Aconitum species. The present review endeavored the major research carried out in the field of genetic resource characterization, pharmacological properties, phytochemistry, major factors influencing quantity, biosynthetic pathways and processing methods for recovery of active ingredients, variety improvement, propagation methods, and important metabolite production through cell/organ culture of various Aconitum species. More than 450 derivatives of aconitine-type C19 and C20-diterpenoid alkaloids along with a few other non-alkaloidal compounds, such as phenylpropanoids, flavonoids, terpenoids, and fatty acids, have been identified in the genus. A few Aconitum species and their common diterpenoid alkaloid compounds are also well characterized for analgesic, inflammatory and cytotoxic properties. However, the different isolated compound needs to be validated for supporting other traditional therapeutical uses of the plant species. Aconitine alkaloids shared common biosynthesis pathway, but their diversification mechanism remains unexplored in the genus. Furthermore, the process needs to be developed on secondary metabolite recovery, mass-scale propagation methods, and agro-technologies for maintaining the quality of products. Many species are losing their existence in nature due to over-exploitation or anthropogenic factors; thus, temporal monitoring of the population status in its habitat, and suitable management programs for ascertaining conservation needs to be developed.

Neuropharmacological Potential of Diterpenoid Alkaloids

This study provides a narrative review of diterpenoid alkaloids (DAs), a family of extremely important natural products found predominantly in some species of Aconitum and Delphinium (Ranunculaceae). DAs have long been a focus of research attention due to their numerous intricate structures and diverse biological activities, especially in the central nervous system (CNS). These alkaloids originate through the amination reaction of tetra or pentacyclic diterpenoids, which are classified into three categories and 46 types based on the number of carbon atoms in the backbone structure and structural differences. The main chemical characteristics of DAs are their heterocyclic systems containing β-aminoethanol, methylamine, or ethylamine functionality. Although the role of tertiary nitrogen in ring A and the polycyclic complex structure are of great importance in drug-receptor affinity, in silico studies have emphasized the role of certain sidechains in C13, C14, and C8. DAs showed antiepileptic effects in preclinical studies mostly through Na⁺ channels. Aconitine (1) and 3-acetyl aconitine (2) can desensitize Na⁺ channels after persistent activation. Lappaconitine (3), N-deacetyllapaconitine (4), 6-benzoylheteratisine (5), and 1-benzoylnapelline (6) deactivate these channels. Methyllycaconitine (16), mainly found in Delphinium species, possesses an extreme affinity for the binding sites of α7 nicotinic acetylcholine receptors (nAChR) and contributes to a wide range of neurologic functions and the release of neurotransmitters. Several DAs such as bulleyaconitine A (17), (3), and mesaconitine (8) from Aconitum species have a drastic analgesic effect. Among them, compound 17 has been used in China for decades. Their effect is explained by increasing the release of dynorphin A, activating the inhibitory noradrenergic neurons in the β-adrenergic system, and preventing the transmission of pain messages by inactivating the Na⁺ channels that have been stressed. Acetylcholinesterase inhibitory, neuroprotective, antidepressant, and anxiolytic activities are other CNS effects that have been investigated for certain DAs. However, despite various CNS effects, recent advances in developing new drugs from DAs were insignificant due to their neurotoxicity.

TcbZIP60 positively regulates pyrethrins biosynthesis in Tanacetum cinerariifolium

Pyrethrins, synthesized in the perennial plant Tanacetum cinerariifolium, are a class of terpene mixtures with high insecticidal activity and low human toxicity, which are widely used in plant-derived pesticides. Numerous studies have identified multiple pyrethrins biosynthesis enzymes, which can be enhanced by exogenous hormones such as methyl jasmonate (MeJA). However, the mechanism by which hormone signaling regulates pyrethrins biosynthesis and the potential involvement of certain transcription factors (TFs) remain unclear. In this study, we found that the expression level of a TF in T. cinerariifolium was significantly increased after treatment with plant hormones (MeJA, abscisic acid). Subsequent analysis identified this TF as a member of the basic region/leucine zipper (bZIP) family and was thus named TcbZIP60. TcbZIP60 was localized in the nucleus, suggesting that it is involved in the transcription process. The expression profiles of TcbZIP60 were similar to those of pyrethrins synthesis genes in different flower organs and at different flowering stages. Furthermore, TcbZIP60 could directly bind to the E-box/G-box motifs in the promoters of the pyrethrins synthesis genes TcCHS and TcAOC to activate their expression. Transient overexpression of TcbZIP60 increased the expression levels of pyrethrins biosynthesis genes, leading to the significant accumulation of pyrethrins. Silencing of TcbZIP60 significantly downregulated pyrethrins accumulation and the expression of related genes. Overall, our results reveal a novel TF, TcbZIP60, that regulates both the terpenoid and jasmonic acid pathways of pyrethrins biosynthesis in T. cinerariifolium.

The clonal growth in Aconitum carmichaelii Debx

Aconitum carmichaelii Debx. is used as traditional herbal medicine in China, Japan, and other Asian countries. A. carmichaelii has two modes for reproduction: sexual reproduction with seed and vegetative reproduction with vegetative propagules. The vegetative propagules are belowground and invisible. To date, only a handful of studies for the clonal growth are available. In this study, we investigated the clonal growth by anatomical and morphological changes. Results revealed that the axillary bud appeared on the rhizome. Furthermore, the axillary meristem in the axillary bud differentiated a bud upwards and an adventitious root (AR) downwards. The AR expanded to a tuberous root in order to provide the bud nutrients for the new plant. The AR branched LRs. In addition, some lateral roots (LRs) on the AR also swelled. Both the AR and LR were found to follow a similar pattern of development. However, high lignification in the stele region of LRs inhibited further expansion. AR development was attributed to activities of the cambium and meristem cell, starch accumulation, stele lignification, and a polyarch stele. Our study not only provides a better understanding of clonal growth but also provides clues to explore the regulatory mechanisms underlying AR development in A. carmichaelii.

De novo and comparative transcriptomic analysis explain morphological differences in Panax notoginseng taproots

Background

Panax notoginseng (Burk.) F. H. Chen (PN) belonging to the genus Panax of family Araliaceae is widely used in traditional Chinese medicine to treat various diseases. PN taproot, as the most vital organ for the accumulation of bioactive components, presents a variable morphology (oval or long), even within the same environment. However, no related studies have yet explained the molecular mechanism of phenotypic differences. To investigate the cause of differences in the taproot phenotype, de novo and comparative transcriptomic analysis on PN taproot was performed.

Results

A total of 133,730,886 and 114,761,595 paired-end clean reads were obtained based on high-throughput sequencing from oval and long taproot samples, respectively. 121,955 unigenes with contig N50 = 1,774 bp were generated by using the de novo assembly transcriptome, 63,133 annotations were obtained with the BLAST. And then, 42 genes belong to class III peroxidase (PRX) gene family, 8 genes belong to L-Ascorbate peroxidase (APX) gene family, and 55 genes belong to a series of mitogen-activated protein kinase (MAPK) gene family were identified based on integrated annotation results. Differentially expressed genes analysis indicated substantial up-regulation of PnAPX3 and PnPRX45, which are related to reactive oxygen species metabolism, and the PnMPK3 gene, which is related to cell proliferation and plant root development, in long taproots compared with that in oval taproots. Furthermore, the determination results of real-time quantitative PCR, enzyme activity, and H₂O₂ content verified transcriptomic analysis results.

Conclusion

These results collectively demonstrate that reactive oxygen species (ROS) metabolism and the PnMPK3 gene may play vital roles in regulating the taproot phenotype of PN. This study provides further insights into the genetic mechanisms of phenotypic differences in other species of the genus Panax.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08283-w.

Functional identification of the terpene synthase family involved in diterpenoid alkaloids biosynthesis in Aconitum carmichaelii

Aconitum carmichaelii is a high-value medicinal herb widely used across China, Japan, and other Asian countries. Aconitine-type diterpene alkaloids (DAs) are the characteristic compounds in Aconitum. Although six transcriptomes, based on short-read next generation sequencing technology, have been reported from the Aconitum species, the terpene synthase (TPS) corresponding to DAs biosynthesis remains unidentified. We apply a combination of Pacbio isoform sequencing and RNA sequencing to provide a comprehensive view of the A. carmichaelii transcriptome. Nineteen TPSs and five alternative splicing isoforms belonging to TPS-b, TPS-c, and TPS-e/f subfamilies were identified. In vitro enzyme reaction analysis functional identified two sesqui-TPSs and twelve diTPSs. Seven of the TPS-c subfamily genes reacted with GGPP to produce the intermediate ent-copalyl diphosphate. Five AcKSLs separately reacted with ent-CPP to produce ent-kaurene, ent-atiserene, and ent-13-epi-sandaracopimaradie: a new diterpene found in Aconitum. AcTPSs gene expression in conjunction DAs content analysis in different tissues validated that ent-CPP is the sole precursor to all DAs biosynthesis, with AcKSL1, AcKSL2s and AcKSL3-1 responsible for C20 atisine and napelline type DAs biosynthesis, respectively. These data clarified the molecular basis for the C20-DAs biosynthetic pathway in A. carmichaelii and pave the way for further exploration of C19-DAs biosynthesis in the Aconitum species.

Transcriptome repository of North-Western Himalayan endangered medicinal herbs: a paramount approach illuminating molecular perspective of phytoactive molecules and secondary metabolism

Medicinal plants of the North-Western Himalayan region are known for their unprecedented biodiversity and valuable secondary metabolites that are unique to this dynamic geo-climatic region. From ancient times these medicinal herbs have been used traditionally for their therapeutic potentials. But from the last 2 decades increasing pharmaceutical demand, illegal and unorganized trade of these medicinal plants have accelerated the rate of over-exploitation in a non-scientific manner. In addition, climate change and anthropogenic activities also affected their natural habitat and driving most of these endemic plant species to critically endangered that foresee peril of mass extinction from this eco-region. Hence there is an urgent need for developing alternative sustainable approaches and policies to utilize this natural bioresource ensuring simultaneous conservation. Hither, arise the advent of sequencing-based transcriptomic studies significantly contributes to better understand the background of important metabolic pathways and related genes/enzymes of high-value medicinal herbs, in the absence of genomic information. The use of comparative transcriptomics in conjunction with biochemical techniques in North-Western Himalayan medicinal plants has resulted in significant advances in the identification of the molecular players involved in the production of secondary metabolic pathways over the last decade. This information could be used to further engineer metabolic pathways and breeding programs, ultimately leading to the development of in vitro systems dedicated to the production of pharmaceutically important secondary metabolites at the industrial level. Collectively, successful adoption of these approaches can certainly ensure the sustainable utilization of Himalayan bioresource by reducing the pressure on the wild population of these critically endangered medicinal herbs. This review provides novel insight as a transcriptome-based bioresource repository for the understanding of important secondary metabolic pathways genes/enzymes and metabolism of endangered high-value North-Western Himalayan medicinal herbs, so that researchers across the globe can effectively utilize this information for devising effective strategies for the production of pharmaceutically important compounds and their scale-up for sustainable usage and take a step forward in omics-based conservation genetics.

Structural diversity, bioactivities, and biosynthesis of natural diterpenoid alkaloids

Covering: 2009 to 2018. Diterpenoid alkaloids, originating from the amination of natural tetracyclic diterpenes, are a diverse class of compounds having complex structural features with many stereocenters. The important pharmacological activities and structural complexity of the diterpenoid alkaloids have long interested scientists due to their medicinal uses, infamous toxicity, and unique biosynthesis. Since 2009, 373 diterpenoid alkaloids, assigned to 46 skeletons, have been isolated and identified from plants mostly in the Ranunculaceae family. The names, classes, molecular weight, molecular formula, NMR data, and plant sources of these diterpene alkaloids are collated here. This review will be a detailed update of the naturally occurring diterpene alkaloids reported from the plant kingdom from 2009-2018, providing an in-depth discussion of their diversity, biological activities, pharmacokinetics, toxicity, application, evolution, and biosynthesis.

Application of DNA barcoding and metabarcoding for species identification in salmon products

Mislabelling is a significant manifestation of food fraud. Traditional Sanger sequencing technology is the gold standard for seafood species identification. However, this method is not suitable for analysing processed samples that may contain more than one species. This study tested the feasibility of next-generation sequencing in identifying mixed salmon products. Salmon samples containing up to eight species were amplified using 16S rRNA mini-barcode primers, and sequenced on an Illumina HiSeq2500 platform. All species were accurately identified, and mixtures as low as 1% (w/w) could be detected. Furthermore, this study conducted a market survey of 32 products labelled as salmon. For pure and mixed fish products, Sanger and next-generation sequencing techniques were respectively used for species identification, and for NGS results, we also used real-time PCR method to cross-validate the mixed products to further verify the accuracy of the DNA metabarcoding technology established in this study. DNA barcoding and metabarcoding of commercial salmon food products revealed the presence of mislabelling in 16 of 32 (50%) samples. The developed DNA barcoding and metabarcoding methods are useful for the identification of salmon species in food and can be used for quality control of various types of salmon products.

Comparative transcriptome analysis infers bulb derived in vitro cultures as a promising source for sipeimine biosynthesis in Fritillaria cirrhosa D. Don (Liliaceae, syn. Fritillaria roylei Hook.) - High value Himalayan medicinal herb

Fritillaria cirrhosa D. Don (Liliaceae, syn. Fritillaria roylei Hook.) is a critically endangered medicinal herb of immense importance due to its pharmaceutical bioactive compound, especially sipeimine, used for the treatment of chronic respiratory disorders. However, the industrial demand for sipeimine solely depends on its endangered natural habitat. Therefore; there is an utmost need for its biodiversity conservation as well as for the sustainable utilization of phytochemicals. Plant cell culture and transcriptomics-based molecular bioprospection of key regulatory genes involved in sipeimine biosynthesis as such will play a crucial role in exploring the unexplored traits, that are in supply crisis or nearly in extinction stage. De novo comparative transcriptome sequencing of the bulb (in vivo), callus, and regenerated plantlets (in vitro) resulted in more than 150 million high-quality paired-end clean reads that assembled into final 31,428 transcripts. Functional annotation and unigenes classification with multiple public databases such as KEGG, Refseq, Uniprot, TAIR, GO, and COG, etc. along with chemical structures and functional biocatalytic activity analysis of different steroidal alkaloids facilitated the identification of 30 unigenes specific to sipeimine biosynthesis. Additionally, ABC transporters and TFs like bHLH, MYC, MYB, and WRKY suggests their possible role in metabolite translocation and regulation in vivo as well as in vitro tissues. Differential gene expression and quantitative analysis revealed that the MVA pathway probably the predominant route for 5C intermediate (IPP & DMAPP) biosynthesis. Further, the genes involved in the downstream biosynthesis pathway viz. SQLE, CAS1, SMT1, SMO1, SMO2, SC5DL, DHCR7, DHCR24, CYP710A, 3β-HSD, CYP90D2, and CYP374A6 shown similar expression pattern with RNA-Seq and qRT-PCR findings. The positive correlation between higher expression of proposed biosynthetic pathway genes and relatively higher accumulation of sipeimine in differentiated naturally grown bulb tissues (in vivo), undifferentiated cells (callus), and de-differentiated tissues i.e. regenerated plantlets (in vitro) has been evident from the present study. Comprehensive genomic resources created in F. cirrhosa will provide strong evidence of bulb derived in vitro culture as an alternative promising source for steroidal alkaloids biosynthesis and metabolite upscaling through genetic and metabolic engineering.

De novo RNA sequencing and analysis reveal the putative genes involved in diterpenoid biosynthesis in Aconitum vilmorinianum roots

In this study, the putative genes involved in diterpenoid alkaloids biosynthesis in A. vilmorinianum roots were revealed by transcriptome sequencing. 59.39 GB of clean bases and 119,660 unigenes were assembled, of which 69,978 unigenes (58.48%) were annotated. We identified 27 classes of genes (139 candidate genes) involved in the synthesis of diterpenoid alkaloids, including the mevalonate (MVA) pathway, the methylerythritol 4-phosphate (MEP) pathway, the farnesyl diphosphate regulatory pathway, and the diterpenoid scaffold synthetic pathway. 12 CYP450 genes were identified. We found that hydroxymethylglutaryl-CoA reductase was the key enzyme in MVA metabolism, which was regulated by miR6300. Transcription factors, such as bHLH, AP2/EREBP, and MYB, used to synthesize the diterpenes were analyzed.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02646-6.

Identification and characterization of SlbHLH, SlDof and SlWRKY transcription factors interacting with SlDPD gene involved in costunolide biosynthesis in Saussurea lappa

The genes involved in costunolide biosynthesis in Saussurea lappa have been identified recently by our lab. However, the study of transcriptional regulators of these genes was lacking for better opportunities for engineering the pharmacologically important biosynthetic pathway. Therefore, we cloned the promoter region of diphosphomevalonate decarboxylase gene (DPD) and analyzed its cis-acting regulatory elements to reveal the potential transcription factor (TF) binding sites for Dof, bHLH and WRKY family proteins in the gene promoter. The transcriptome study approach followed by the hidden Markov model based search, digital gene expression, co-expression network analysis, conserved domain properties and evolutionary analyses were carried out to screen out seven putative TFs for the DPD-TF interaction studies. Yeast one-hybrid assays were performed and three TFs were reported, namely, SlDOF2, SlbHLH3 and SlWRKY2 from Dof, bHLH and WRKY families, respectively that interacted positively with the DPD gene of the costunolide biosynthetic pathway. The tissue specific relative gene expression studies also supported the linked co-expression of the gene and its interacting TFs The present report will improve the understanding of transcriptional regulation pattern of costunolide biosynthetic pathway.

De novo transcriptome analysis of Lantana camara L. revealed candidate genes involved in phenylpropanoid biosynthesis pathway

Lantana camara L. is an economically important essential oil producing plant belonging to family Verbenaceae. It is used in medication for treating various diseases like cancer, ulcers, tumor, asthma and fever. The plant is a useful source of essential bioactive compounds such as steroids, flavonoids and phenylpropanoid glycosides etc. Nonetheless, very little is known about the genomic or transcriptomic resources of L. camara, and this might be the reason of hindering molecular studies leading to identification of improved lines. Here we used Illumina sequencing platform and performed the L. camara leaf (LCL) and root (LCR) de novo transcriptome analyses. A total of 70,155,594 and 84,263,224 clean reads were obtained and de novo assembly generated 72,877 and 513,985 unigenes from leaf (LCL) and root (LCR) respectively. Furthermore, the pathway analysis revealed the presence of 229 and 943 genes involved in the phenylpropanoid biosynthesis in leaf and root tissues respectively. Similarity search was performed against publically available genome databases and best matches were found with Sesamum indicum (67.5%) that were much higher than that of Arabidopsis thaliana (3.9%). To the best of our knowledge, this is the first comprehensive transcriptomic analysis of leaf and root tissues of this non-model plant from family Verbenaceae and may serve as a baseline for further molecular studies.

De novo leaf transcriptome assembly of Bougainvillea spectabilis for the identification of genes involves in the secondary metabolite pathways

Bougainvillea spectabilis is known as a vital medicinal, ornamental as well as an essential oil producing plant. It is also a rich source of important secondary metabolites with several therapeutic properties. Various studies on its pharmacological and toxicological aspects have been published but there is no genomic or transcriptomic resource available in the public databases. To address this important issue, the de-novo transcriptome assembly of B. spectabilis leaf tissue has been done for the identification of genes involved in various important secondary metabolites, Single nucleotide polymorphism (SNPs) and Simple sequence repeats (SSRs). The transcriptome sequencing of B. spectabilis leaf tissue generated 79,811,024 raw reads with GC value 42.77%. The transcriptomic assembly was performed by Trinity software which generated 100,374 transcripts and 99,793 unigenes with minimum and maximum length of 201 bp and 13,237 bp and N50 value of 1470 and 1472 respectively. Annotation of these unigenes was performed using seven databases including NR, PFAM, GO and KEGG. Approximately, 44,302 unigenes were annotated in GO database. The KEGG pathway analysis revealed 23,102 unigenes in which 19,054 genes were assigned to five groups in KEGG and 130 biochemical pathways. The highest group among the five groups was Metabolism with 9230 unigenes. Moreover, about 63,226 SNPs and 30,333 SSRs in the leaf transcriptome of B. spectabilis were identified. To the best of my understanding it will be the first comprehensive transcriptome analysis of B. spectabilis from family Nyctaginaceae which will help as a reference line for further genomic and transcriptomic studies.

Biotic elicitors enhance diosgenin production in Helicteres isora L. suspension cultures via up-regulation of CAS and HMGR genes

In an attempt to find an alternative and potent source of diosgenin, a steroidal saponin in great demand for its pharmaceutical importance, Helicteres isora suspension cultures were explored for diosgenin extraction. The effect of biotic elicitors on the biosynthesis of diosgenin, in suspension cultures of H. isora was studied. Bacterial as well as fungal elicitors such as Escherichia coli, Bacillus subtilis, Saccharomyces cerevisiae and Aspergillus niger were applied at varying concentrations to investigate their effects on diosgenin content. The HPLC based quantification of the treated samples proved that amongst the biotic elicitors, E. coli (1.5%) proved best with a 9.1-fold increase in diosgenin content over respective control cultures. Further, the scaling-up of the suspension culture to shake-flask and ultimately to bioreactor level were carried out for production of diosgenin. During all the scaling-up stages, diosgenin yield obtained was in the range between 7.91 and 8.64 mg l^-1, where diosgenin content was increased with volume of the medium. The quantitative real-time PCR (qRT-PCR) analysis showed biotic elicitors induced the expression levels of regulatory genes in diosgenin biosynthetic pathway, the 3-hydroxy-3-methylglutaryl-CoA reductase (HMGR) and cycloartenol synthase (CAS), which can be positively correlated with elicited diosgenin contents in those cultures. The study holds significance as H. isora represents a cleaner and easy source of diosgenin where unlike other traditional sources, it is not admixed with other steroidal saponins, and the scaled-up levels of diosgenin achieved herein have the potential to be explored commercially.

Comparative transcriptome and metabolome analyses provide new insights into the molecular mechanisms underlying taproot thickening in Panax notoginseng

BACKGROUND

Taproot thickening is a complex biological process that is dependent on the coordinated expression of genes controlled by both environmental and developmental factors. Panax notoginseng is an important Chinese medicinal herb that is characterized by an enlarged taproot as the main organ of saponin accumulation. However, the molecular mechanisms of taproot enlargement are poorly understood.

RESULTS

A total of 29,957 differentially expressed genes (DEGs) were identified during the thickening process in the taproots of P. notoginseng. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment revealed that DEGs associated with "plant hormone signal transduction," "starch and sucrose metabolism," and "phenylpropanoid biosynthesis" were predominantly enriched. Further analysis identified some critical genes (e.g., RNase-like major storage protein, DA1-related protein, and Starch branching enzyme I) and metabolites (e.g., sucrose, glucose, fructose, malate, and arginine) that potentially control taproot thickening. Several aspects including hormone crosstalk, transcriptional regulation, homeostatic regulation between sugar and starch, and cell wall metabolism, were identified as important for the thickening process in the taproot of P. notoginseng.

CONCLUSION

The results provide a molecular regulatory network of taproot thickening in P. notoginseng and facilitate the further characterization of the genes responsible for taproot formation in root medicinal plants or crops.

Transcriptome analysis reveals in vitro-cultured regeneration bulbs as a promising source for targeted Fritillaria cirrhosa steroidal alkaloid biosynthesis

The bulbs of Fritillaria cirrhosa is wildly used in traditional Chinese medicine to treat lung-related disease, which has recently been found to have antitussive, anti-inflammatory, antihypertensive and anti-tumor activity. Steroidal alkaloids are the major effective ingredients of F. cirrhosa. In the current study, we demonstrated an efficient strategy for F. cirrhosa bulb regeneration in vitro by cytokinin/auxin induction. Our data showed that the regenerated bulbs accumulated higher alkaloid content that the wild ones. We further performed RNA-seq and bioinformatics analysis to study the gene expression profile, especially those related to alkaloids biosynthesis. KEGG pathway annotation identified genes related to "Metabolic pathways" were the most abundant (2644, 26.0%), followed by those for "Biosynthesis of secondary metabolites" (1319, 13.0%) among the 113,865 unigenes identified. Further analysis suggested MEP pathway, other than MVA pathway, might be the major route for steroidal alkaloid biosynthesis of F. cirrhosa, as all the key genes in this pathway were found to be unregulated in our study. We also showed that accumulation of different phytochemicals was linked to plant hormone addition. Our current study demonstrated that in vitro cultivation is a promising strategy for mass production of F. cirrhosa steroidal alkaloids for pharmacological industry.

Comparative transcriptomics uncovers differences in photoautotrophic versus photoheterotrophic modes of nutrition in relation to secondary metabolites biosynthesis in Swertia chirayita

Swertia chirayita is a high-value medicinal herb exhibiting antidiabetic, hepatoprotective, anticancer, antiediematogenic and antipyretic properties. Scarcity of its plant material has necessitated in vitro production of therapeutic metabolites; however, their yields were low compared to field grown plants. Possible reasons for this could be differences in physiological and biochemical processes between plants grown in photoautotrophic versus photoheterotrophic modes of nutrition. Comparative transcriptomes of S. chirayita were generated to decipher the crucial molecular components associated with the secondary metabolites biosynthesis. Illumina HiSeq sequencing yielded 57,460 and 43,702 transcripts for green house grown (SCFG) and tissue cultured (SCTC) plants, respectively. Biological role analysis (GO and COG assignments) revealed major differences in SCFG and SCTC transcriptomes. KEGG orthology mapped 351 and 341 transcripts onto secondary metabolites biosynthesis pathways for SCFG and SCTC transcriptomes, respectively. Nineteen out of 30 genes from primary metabolism showed higher in silico expression (FPKM) in SCFG versus SCTC, possibly indicating their involvement in regulating the central carbon pool. In silico data were validated by RT-qPCR using a set of 16 genes, wherein 10 genes showed similar expression pattern across both the methods. Comparative transcriptomes identified differentially expressed transcription factors and ABC-type transporters putatively associated with secondary metabolism in S. chirayita. Additionally, functional classification was performed using NCBI Biosystems database. This study identified the molecular components implicated in differential modes of nutrition (photoautotrophic vs. photoheterotrophic) in relation to secondary metabolites production in S. chirayita.

De Novo RNA Sequencing and Expression Analysis of Aconitum carmichaelii to Analyze Key Genes Involved in the Biosynthesis of Diterpene Alkaloids

Aconitum carmichaelii is an important medicinal herb used widely in China, Japan, India, Korea, and other Asian countries. While extensive research on the characterization of metabolic extracts of A. carmichaelii has shown accumulation of numerous bioactive metabolites including aconitine and aconitine-type diterpene alkaloids, its biosynthetic pathway remains largely unknown. Biosynthesis of these secondary metabolites is tightly controlled and mostly occurs in a tissue-specific manner; therefore, transcriptome analysis across multiple tissues is an attractive method to identify the molecular components involved for further functional characterization. In order to understand the biosynthesis of secondary metabolites, Illumina-based deep transcriptome profiling and analysis was performed for four tissues (flower, bud, leaf, and root) of A. carmichaelii, resulting in 5.5 Gbps clean RNA-seq reads assembled into 128,183 unigenes. Unigenes annotated as possible rate-determining steps of an aconitine-type biosynthetic pathway were highly expressed in the root, in accordance with previous reports describing the root as the accumulation site for these metabolites. We also identified 21 unigenes annotated as cytochrome P450s and highly expressed in roots, which represent candidate unigenes involved in the diversification of secondary metabolites. Comparative transcriptome analysis of A. carmichaelii with A. heterophyllum identified 20,232 orthogroups, representing 30,633 unigenes of A. carmichaelii, gene ontology enrichment analysis of which revealed essential biological process together with a secondary metabolic process to be highly enriched. Unigenes identified in this study are strong candidates for aconitine-type diterpene alkaloid biosynthesis, and will serve as useful resources for further validation studies.

Modular Design of Picroside-II Biosynthesis Deciphered through NGS Transcriptomes and Metabolic Intermediates Analysis in Naturally Variant Chemotypes of a Medicinal Herb, Picrorhiza kurroa

Picroside-II (P-II), an iridoid glycoside, is used as an active ingredient of various commercial herbal formulations available for the treatment of liver ailments. Despite this, the knowledge of P-II biosynthesis remains scarce owing to its negligence in Picrorhiza kurroa shoots which sets constant barrier for function validation experiments. In this study, we utilized natural variation for P-II content in stolon tissues of different P. kurroa accessions and deciphered its metabolic route by integrating metabolomics of intermediates with differential NGS transcriptomes. Upon navigating through high vs. low P-II content accessions (1.3-2.6%), we have established that P-II is biosynthesized via degradation of ferulic acid (FA) to produce vanillic acid (VA) which acts as its immediate biosynthetic precursor. Moreover, the FA treatment in vitro at 150 μM concentration provided further confirmation with 2-fold rise in VA content. Interestingly, the cross-talk between different compartments of P. kurroa, i.e., shoots and stolons, resolved spatial complexity of P-II biosynthesis and consequently speculated the burgeoning necessity to bridge gap between VA and P-II production in P. kurroa shoots. This work thus, offers a forward looking strategy to produce both P-I and P-II in shoot cultures, a step toward providing a sustainable production platform for these medicinal compounds via-à-vis relieving pressure from natural habitat of P. kurroa.

Abscisic-acid-dependent basic leucine zipper (bZIP) transcription factors in plant abiotic stress

One of the major causes of significant crop loss throughout the world is the myriad of environmental stresses including drought, salinity, cold, heavy metal toxicity, and ultraviolet-B (UV-B) rays. Plants as sessile organisms have evolved various effective mechanism which enable them to withstand this plethora of stresses. Most of such regulatory mechanisms usually follow the abscisic-acid (ABA)-dependent pathway. In this review, we have primarily focussed on the basic leucine zipper (bZIP) transcription factors (TFs) activated by the ABA-mediated signalosome. Upon perception of ABA by specialized receptors, the signal is transduced via various groups of Ser/Thr kinases, which phosphorylate the bZIP TFs. Following such post-translational modification of TFs, they are activated so that they bind to specific cis-acting sequences called abscisic-acid-responsive elements (ABREs) or GC-rich coupling elements (CE), thereby influencing the expression of their target downstream genes. Several in silico techniques have been adopted so far to predict the structural features, recognize the regulatory modification sites, undergo phylogenetic analyses, and facilitate genome-wide survey of TF under multiple stresses. Current investigations on the epigenetic regulation that controls greater accessibility of the inducible regions of DNA of the target gene to the bZIP TFs exclusively under stress situations, along with the evolved stress memory responses via genomic imprinting mechanism, have been highlighted. The potentiality of overexpression of bZIP TFs, either in a homologous or in a heterologous background, in generating transgenic plants tolerant to various abiotic stressors have also been addressed by various groups. The present review will provide a coherent documentation on the functional characterization and regulation of bZIP TFs under multiple environmental stresses, with the major goal of generating multiple-stress-tolerant plant cultivars in near future.

Transcriptome-wide mining suggests conglomerate of genes associated with tuberous root growth and development in Aconitum heterophyllum Wall

Tuberous roots of Aconitum heterophyllum constitute storage organ for secondary metabolites, however, molecular components contributing to their formation are not known. The transcriptomes of A. heterophyllum were analyzed to identify possible genes associated with tuberous root development by taking clues from genes implicated in other plant species. Out of 18 genes, eight genes encoding GDP-mannose pyrophosphorylase (GMPase), SHAGGY, Expansin, RING-box protein 1 (RBX1), SRF receptor kinase (SRF), β-amylase, ADP-glucose pyrophosphorylase (AGPase) and Auxin responsive factor 2 (ARF2) showed higher transcript abundance in roots (13-171 folds) compared to shoots. Comparative expression analysis of those genes between tuberous root developmental stages showed 11-97 folds increase in transcripts in fully developed roots compared to young rootlets, thereby implying their association in biosynthesis, accumulation and storage of primary metabolites towards root biomass. Cluster analysis revealed a positive correlation with the gene expression data for different stages of tuberous root formation in A. heterophyllum. The outcome of this study can be useful in genetic improvement of A. heterophyllum for root biomass yield.

Renaissance in phytomedicines: promising implications of NGS technologies

Medicinal plant research is growing significantly in faith to discover new and more biologically compatible phytomedicines. Deposition of huge genome/trancriptome sequence data assisted by NGS technologies has revealed the new possibilities for producing upgraded bioactive molecules in medicinal plants. Growing interest of investors and consumers in the herbal drugs raises the need for extensive research to open the facts and details of every inch of life canvas of medicinal plants to produce improved quality of phytomedicines. As in agriculture crops, knowledge emergence from medicinal plant's genome/transcriptome, can be used to assure their amended quality and these improved varieties are then transported to the fields for cultivation. Genome studies generate huge sequence data which can be exploited further for obtaining information regarding genes/gene clusters involved in biosynthesis as well as regulation. This can be achieved rapidly at a very large scale with NGS platforms. Identification of new RNA molecules has become possible, which can lead to the discovery of novel compounds. Sequence information can be combined with advanced phytochemical and bioinformatics tools to discover functional herbal drugs. Qualitative and quantitative analysis of small RNA species put a light on the regulatory aspect of biosynthetic pathways for phytomedicines. Inter or intra genomic as well as transcriptomic interactive processes for biosynthetic pathways can be elucidated in depth. Quality management of herbal material will also become rapid and high throughput. Enrichment of sequence information will be used to engineer the plants to get more efficient phytopharmaceuticals. The present review comprises of role of NGS technologies to boost genomic studies of pharmaceutically important plants and further, applications of sequence information aiming to produce enriched phytomedicines. Emerging knowledge from the medicinal plants genome/transcriptome can give birth to deep understanding of the processes responsible for biosynthesis of medicinally important compounds.

Molecular dissection of pathway components unravel atisine biosynthesis in a non-toxic Aconitum species, A. heterophyllum Wall

Aconitum heterophyllum is an important component for various herbal drug formulations due to the occurrence of non-toxic aconites including marker compound, atisine. Despite huge pharmacological potential, the reprogramming of aconites production is limited due to lack of understanding on their biosynthesis. To address this problem, we have proposed here the complete atisine biosynthetic pathway for the first time connecting glycolysis, MVA/MEP, serine biosynthesis and diterpene biosynthetic pathways. The transcript profiling revealed phosphorylated pathway as a major contributor towards serine production in addition to repertoire of genes in glycolysis (G6PI, PFK, ALD and ENO), serine biosynthesis (PGDH and PSAT) and diterpene biosynthesis (KO and KH) sharing a similar pattern of expression (2-4-folds) in roots compared to shoots vis-à-vis atisine content (0-0.37 %). Quantification of steviol and comparative analysis of shortlisted genes between roots of high (0.37 %) vs low (0.14 %) atisine content accessions further confirmed the route of atisine biosynthesis. The results showed 6-fold increase in steviol content and 3-62-fold up-regulation of all the selected genes in roots of high content accession ascertaining their association towards atisine production. Moreover, significant positive correlations were observed between selected genes suggesting their co-expression and crucial role in atisine biosynthesis. This study, thus, offers unprecedented opportunities to explore the selected candidate genes for enhanced production of atisine in cultivated plant cells.