Cloning and expression analysis of ten genes associated with picrosides biosynthesis in Picrorhiza kurrooa

Transcriptome analysis and exploration of genes involved in the biosynthesis of secoiridoids in Gentiana rhodantha

Gentiana rhodantha is a medicinally important perennial herb used as traditional Chinese and ethnic medicines. Secoiridoids are one of the major bioactive compounds in G. rhodantha. To better understand the secoiridoid biosynthesis pathway, we generated transcriptome sequences from four organs (root, leaf, stem and flower), followed by the de novo sequence assembly. We verified 8-HGO (8-hydroxygeraniol oxidoreductase), which may encode key enzymes of the secoiridoid biosynthesis by qRT-PCR. The mangiferin, swertiamarin and loganic acid contents in root, stem, leaf, and flower were determined by HPLC. The results showed that there were 47,871 unigenes with an average length of 1,107.38 bp. Among them, 1,422 unigenes were involved in 25 standard secondary metabolism-related pathways in the KEGG database. Furthermore, we found that 1,005 unigenes can be divided into 66 transcription factor (TF) families, with no family members exhibiting significant organ-specificity. There were 54 unigenes in G. rhodantha that encoded 17 key enzymes of the secoiridoid biosynthetic pathway. The qRT-PCR of the 8-HGO and HPLC results showed that the relative expression and the mangiferin, swertiamarin, and loganic acid contents of the aerial parts were higher than in the root. Six types of SSR were identified by SSR analysis of unigenes: mono-nucleoside repeat SSR, di-nucleoside repeat SSR, tri-nucleoside repeat SSR, tetra-nucleoside repeat SSR, penta-nucleoside repeat SSR, and hexa-nucleoside repeat SSR. This report not only enriches the Gentiana transcriptome database but helps further study the function and regulation of active component biosynthesis of G. rhodantha.

Iridoid glycosides from Picrorhiza genus endemic to the Himalayan region: phytochemistry, biosynthesis, pharmacological potential and biotechnological intercessions to boost production

Iridoid glycosides are monoterpenoids synthesized in several plant species known to exhibit a diverse range of pharmacological activities. They are used as important bioactive ingredients in many commercially available drug formulations and as lead compounds in pharmaceutical research. The genus Picrorhiza comprises two medicinally important herbs endemic to the Himalayan region viz. Picrorhiza kurrooa Royle and Picrorhiza scrophulariiflora Hong. The medicinal properties of these two species are mainly due to iridoid glycosides present in their root, rhizome, and leaves. Unregulated harvesting from the wild, habitat specificity, narrow distribution range, small population size and lack of organized cultivation led to the enrolling of these species in the endangered category by the International Union for Conservation of Nature and Natural Resources (IUCN). Therefore, there is a need for immediate biotechnological and molecular interventions. Such intercessions will open up new vistas for large-scale propagation, development of genomic/transcriptomic resources for understanding the biosynthetic pathway, the possibility of genetic/metabolic manipulations, and possible commercialization of iridoid glycosides. The current review article elucidates the phytochemistry and pharmacological importance of iridoid glycosides from the genus Picrorhiza. In addition, the role of biotechnological approaches and opportunities offered by next-generation sequencing technologies in overcoming challenges associated with the genetic engineering of these species are also discussed.

Transcriptome profiling reveal key hub genes in co-expression networks involved in Iridoid glycosides biosynthetic machinery in Picrorhiza kurroa

Picrorhiza kurroa is a medicinal herb rich in hepatoprotective iridoid glycosides, picroside-I (P-I) and picroside-II (P-II). The biosynthetic machinery of picrosides is poorly understood, therefore, 'no-direction' gene co-expression networks were used to extract linked/closed and separated interactions in terpenoid glycosides-specific sub-networks. Transcriptomes generated from different organs, varying for P-I and P-II contents such as shoots grown at 15 and 25 °C and nursery-grown shoots, stolons, and roots resulted in 47,726, 44,958, 40,117, 66,979, and 55,578 annotated transcripts, respectively. Occurrence of 2810 ± 136 nodes and 15,626 ± 696 edges in these networks indicated intense, co-expressed, closed loop interactions. Either deregulation/inhibition of abscisic acid (ABA) biosynthesis/signaling or constitutive degradation of ABA resulted in organ-specific accumulation of P-I and P-II. Biosynthesis, condensation and glucosylation of isoprene units may occur in shoots, roots or stolons; but addition of phenylpropanoid moiety and further modification/s of the iridoid backbone occurs mainly inside vacuoles in roots.

PkGPPS.SSU interacts with two PkGGPPS to form heteromeric GPPS in Picrorhiza kurrooa: Molecular insights into the picroside biosynthetic pathway

Geranyl geranyl pyrophosphate synthase (GGPPS) is known to form an integral subunit of the heteromeric GPPS (geranyl pyrophosphate synthase) complex and catalyzes the biosynthesis of monoterpene in plants. Picrorhiza kurrooa Royle ex Benth., a medicinally important high altitude plant is known for picroside biomolecules, the monoterpenoids. However, the significance of heteromeric GPPS in P. kurrooa still remains obscure. Here, transient silencing of PkGGPPS was observed to reduce picroside-I (P-I) content by more than 60% as well as picroside-II (P-II) by more than 75%. Thus, PkGGPPS was found to be involved in the biosynthesis of P-I and P-II besides other terpenoids. To unravel the mechanism, small subunit of GPPS (PkGPPS.SSU) was identified from P. kurrooa. Protein-protein interaction studies in yeast as well as bimolecular fluorescence complementation (BiFC) in planta have indicated that large subunit of GPPS PkGPPS.LSUs (PkGGPPS1 and PkGGPPS2) and PkGPPS.SSU form a heteromeric GPPS. Presence of similar conserved domains such as light responsive motifs, low temperature responsive elements (LTRE), dehydration responsive elements (DREs), W Box and MeJA responsive elements in the promoters of PkGPPS.LSU and PkGPPS.SSU documented their involvement in picroside biosynthesis. Further, the tissue specific transcript expression analysis vis-à-vis epigenetic regulation (DNA methylation) of promoters as well as coding regions of PkGPPS.LSU and PkGPPS.SSU has strongly suggested their role in picroside biosynthesis. Taken together, the newly identified PkGPPS.SSU formed the heteromeric GPPS by interacting with PkGPPS.LSUs to synthesize P-I and P-II in P. kurrooa.

Growth Kinetics, Metabolites Production and Expression Profiling of Picrosides Biosynthetic Pathway Genes in Friable Callus Culture of Picrorhiza kurroa Royle ex Benth

The rising demand for picrosides commercially and over-exploitation of Picrorhiza kurroa from natural habitat has to initiate alternative strategies for sustainable production of metabolites. In the present research, wild leaf explant of P. kurroa was used to produce friable callus under different culture condition, i.e., dark and light with two temperature variants (15 °C and 25 °C). Afterward, callus cell lines were screened based on growth biomass and metabolites content accumulation. The results revealed, maximum callus growth index along with antioxidant potential (IC₅₀-40.88 μg/mL) and total phenol content (41.35 μg/mg) were observed under dark 25 °C. However, under light 15 °C, highest accumulation of picroside II (0.58 μg/mg), cinnamic acid (0.15 μg/mg), p-hydroxy acetophenone (0.30 μg/mg), total flavonoids (77.30 μg/mg), nitrogen (7.06%), carbohydrates (18.03%), and protein (44.12%) were detected. Major reported metabolite in callus was picroside I (1.63 μg/mg) under dark 15 °C. For the first time, picroside III content (range 0.15-0.56 μg/mg) was also detected and quantified in leaf-derived calli. Expression profiling of picroside biosynthetic pathway genes showed a positive correlation with the observed metabolites. Furthermore, an optimized protocol of metabolites enriched callus biomass could be used as potential strategy for sustainable production of picrosides at commercial scale.

A Single Nucleotide Mutation of the IspE Gene Participating in the MEP Pathway for Isoprenoid Biosynthesis Causes a Green-Revertible Yellow Leaf Phenotype in Rice

Plant isoprenoids are dependent on two independent pathways, the cytosolic mevalonate (MVA) pathway and the plastidic methylerythritol phosphate (MEP) pathway. IspE is one of seven known enzymes in the MEP pathway. Currently, no IspE gene has been identified in rice. In addition, no virescent mutants have been reported to result from a gene mutation affecting the MEP pathway. In this study, we isolated a green-revertible yellow leaf mutant gry340 in rice. The mutant exhibited a reduced level of photosynthetic pigments, and an arrested development of chloroplasts and mitochondria in its yellow leaves. Map-based cloning revealed a missense mutation in OsIspE (LOC_Os01g58790) in gry340 mutant plants. OsIspE is constitutively expressed in all tissues, and its encoded protein is targeted to the chloroplast. Further, the mutant phenotype of gry340 was rescued by introduction of the wild-type gene. Therefore, we have successfully identified an IspE gene in monocotyledons via map-based cloning, and confirmed that the green-revertible yellow leaf phenotype of gry340 does result from a single nucleotide mutation in the IspE gene. In addition, the ispE ispF double mutant displayed an etiolation lethal phenotype, indicating that the isoprenoid precursors from the cytosol cannot efficiently compensate for the deficiency of the MEP pathway in rice chloroplasts. Furthermore, real-time quantitative reverse transcription-PCR suggested that this functional defect in OsIspE affected the expression of not only other MEP pathway genes but also that of MVA pathway genes, photosynthetic genes and mitochondrial genes.

Cloning and Expression Analysis of MEP Pathway Enzyme-encoding Genes in Osmanthus fragrans

The 2-C-methyl-d-erythritol 4-phosphate (MEP) pathway is responsible for the biosynthesis of many crucial secondary metabolites, such as carotenoids, monoterpenes, plastoquinone, and tocopherols. In this study, we isolated and identified 10 MEP pathway genes in the important aromatic plant sweet osmanthus (Osmanthus fragrans). Multiple sequence alignments revealed that 10 MEP pathway genes shared high identities with other reported proteins. The genes showed distinctive expression profiles in various tissues, or at different flower stages and diel time points. The qRT-PCR results demonstrated that these genes were highly expressed in inflorescences, which suggested a tissue-specific transcript pattern. Our results also showed that OfDXS1, OfDXS2, and OfHDR1 had a clear diurnal oscillation pattern. The isolation and expression analysis provides a strong foundation for further research on the MEP pathway involved in gene function and molecular evolution, and improves our understanding of the molecular mechanism underlying this pathway in plants.

Discerning picroside-I biosynthesis via molecular dissection of in vitro shoot regeneration in Picrorhiza kurroa

Expression analysis of primary and secondary metabolic pathways genes vis-à-vis shoot regeneration revealed developmental regulation of picroside-I biosynthesis in Picrorhiza kurroa. Picroside-I (P-I) is an important iridoid glycoside used in several herbal formulations for treatment of various disorders. P-I is synthesized in shoots of Picrorhiza kurroa and Picrorhiza scrophulariiflora. Current study reports on understanding P-I biosynthesis in different morphogenetic stages, viz. plant segment (PS), callus initiation (CI), callus mass (CM), shoot primordia (SP), multiple shoots (MS) and fully developed (FD) stages of P. kurroa. Expression analysis of genes involved in primary and secondary metabolism revealed that genes encoding HMGR, PMK, DXPS, ISPE, GS, G10H, DAHPS and PAL enzymes of MVA, MEP, iridoid and shikimate/phenylpropanoid pathways showed significant modulation of expression in SP, MS and FD stages in congruence with P-I content compared to CM stage. While HK, PK, ICDH, MDH and G6PDH showed high expression in MS and FD stages of P. kurroa, RBA, HisK and CytO showed high expression with progress in regeneration of shoots. Quantitative expression analysis of secondary metabolism genes at two temperatures revealed that 7 genes HMGR, PMK, DXPS, GS, G10H, DAHPS and PAL showed high transcript abundance (32-87-folds) in FD stage derived from leaf and root segments at 15 °C compared to 25 °C in P. kurroa. Further screening of these genes at species level showed high expression pattern in P. kurroa (6-19-folds) vis-à-vis P. scrophulariiflora that was in corroboration with P-I content. Therefore, current study revealed developmental regulation of P-I biosynthesis in P. kurroa which would be useful in designing a suitable genetic intervention study by targeting these genes for enhancing P-I production.

Exogenous feeding of immediate precursors reveals synergistic effect on picroside-I biosynthesis in shoot cultures of Picrorhiza kurroa Royle ex Benth

In the current study, we asked how the supply of immediate biosynthetic precursors i.e. cinnamic acid (CA) and catalpol (CAT) influences the synthesis of picroside-I (P-I) in shoot cultures of P. kurroa. Our results revealed that only CA and CA+CAT stimulated P-I production with 1.6-fold and 4.2-fold, respectively at 2.5 mg/100 mL concentration treatment. Interestingly, feeding CA+CAT not only directed flux towards p-Coumaric acid (p-CA) production but also appeared to trigger the metabolic flux through both shikimate/phenylpropanoid and iridoid pathways by utilizing more of CA and CAT for P-I biosynthesis. However, a deficiency in the supply of either the iridoid or the phenylpropanoid precursor limits flux through the respective pathways as reflected by feedback inhibition effect on PAL and decreased transcripts expressions of rate limiting enzymes (DAHPS, CM, PAL, GS and G10H). It also appears that addition of CA alone directed flux towards both p-CA and P-I production. Based on precursor feeding and metabolic fluxes, a current hypothesis is that precursors from both the iridoid and shikimate/phenylpropanoid pathways are a flux limitation for P-I production in shoot cultures of P. kurroa plants. This work thus sets a stage for future endeavour to elevate production of P-I in cultured plant cells.

NGS Transcriptomes and Enzyme Inhibitors Unravel Complexity of Picrosides Biosynthesis in Picrorhiza kurroa Royle ex. Benth

Picrorhiza kurroa is an important medicinal herb valued for iridoid glycosides, Picroside-I (P-I) and Picroside-II (P-II), which have several pharmacological activities. Genetic interventions for developing a picroside production platform would require knowledge on biosynthetic pathway and key control points, which does not exist as of today. The current study reports that geranyl pyrophosphate (GPP) moiety is mainly contributed by the non-mevalonate (MEP) route, which is further modified to P-I and P-II through phenylpropanoid and iridoid pathways, in total consisting of 41 and 35 enzymatic steps, respectively. The role of the MEP pathway was ascertained through enzyme inhibitors fosmidomycin and mevinolin along with importance of other integrating pathways using glyphosate, aminooxy acetic acid (AOA) and actinomycin D, which overall resulted in 17%-92% inhibition of P-I accumulation. Retrieval of gene sequences for enzymatic steps from NGS transcriptomes and their expression analysis vis-à-vis picrosides content in different tissues/organs showed elevated transcripts for twenty genes, which were further shortlisted to seven key genes, ISPD, DXPS, ISPE, PMK, 2HFD, EPSPS and SK, on the basis of expression analysis between high versus low picrosides content strains of P. kurroa so as to eliminate tissue type/ developmental variations in picrosides contents. The higher expression of the majority of the MEP pathway genes (ISPD, DXPS and ISPE), coupled with higher inhibition of DXPR enzyme by fosmidomycin, suggested that the MEP route contributed to the biosynthesis of P-I in P. kurroa. The outcome of the study is expected to be useful in designing a suitable genetic intervention strategy towards enhanced production of picrosides. Possible key genes contributing to picroside biosynthesis have been identified with potential implications in molecular breeding and metabolic engineering of P. kurroa.

Multiple genes of mevalonate and non-mevalonate pathways contribute to high aconites content in an endangered medicinal herb, Aconitum heterophyllum Wall

Aconitum heterophyllum Wall, popularly known as Atis or Patis, is an important medicinal herb of North-Western and Eastern Himalayas. No information exists on molecular aspects of aconites biosynthesis, including atisine- the major chemical constituent of A. heterophyllum. Atisine content ranged from 0.14% to 0.37% and total alkaloids (aconites) from 0.20% to 2.49% among 14 accessions of A. heterophyllum. Two accessions contained the highest atisine content with 0.30% and 0.37% as well as the highest alkaloids content with 2.22% and 2.49%, respectively. No atisine was detected in leaves and shoots of A. heterophyllum, thereby, suggesting that the biosynthesis and accumulation of aconite alkaloids occur mainly in roots. Quantitative expression analysis of 15 genes of MVA/MEP pathways in roots versus shoots, differing for atisine content (0-2.2 folds) showed 11-100 folds increase in transcript amounts of 4 genes of MVA pathway; HMGS, HMGR, PMK, IPPI, and 4 genes of MEP pathway; DXPS, ISPD, HDS, GDPS, respectively. The overall expression of 8 genes decreased to 5-12 folds after comparative expression analysis between roots of high (0.37%) versus low (0.14%) atisine content accessions, but their relative transcript amounts remained higher in high content accessions, thereby implying their role in atisine biosynthesis and accumulation. PCA analysis revealed a positive correlation between MVA/MEP pathways genes and alkaloids content. The current study provides first report wherein partial sequences of 15 genes of MVA/MEP pathways have been cloned and studied for their possible role in aconites biosynthesis. The outcome of study has potential applications in the genetic improvement of A. heterophyllum.

A phenylalanine ammonia-lyase ortholog (PkPAL1) from Picrorhiza kurrooa Royle ex. Benth: molecular cloning, promoter analysis and response to biotic and abiotic elicitors

Picrorhiza kurrooa Royle ex Benth. is a highly reputed medicinal herb utilised in the preparation of a number of herbal drug formulations, principally due to the presence of novel monoterpene iridoid glycosides kenned as picrosides. Phenylalanine ammonia-lyase catalyses an important rate-limiting step in phenylpropanoid pathway and supplies precursors like cinnamic acid, vanillic acid, ferulic acid, etc., to a variety of secondary metabolites including picrosides. The imperilled status of P. kurrooa coupled with lack of information regarding biogenesis of picrosides necessitates deciphering the biosynthetic pathway for picrosides. In the present study, a PAL gene, designated PkPAL1 was isolated from P. kurrooa. The cDNA is 2312 bp in length, consisting of an ORF of 2142 bp encoding for a 713 amino acid protein having a predicted molecular weight of 77.66 kDa and an isoelectric point of pH 6.82. qRT-PCR analysis of various tissues of P. kurrooa showed that PkPAL1 transcript levels were highest in the leaves, consistent with picroside accumulation pattern. Using Genome walking, a 718 bp promoter region was also isolated resulting in identification of distinct cis-regulatory elements including TGA-element, TGACG-motif, CGTCA-motif, etc. qRT-PCR indicated up-regulation of PkPAL1 by methyl jasmonate, salicylic acid, 2,4-dicholorophenoxy acetic acid and UV-B elicitations that corroborated positively with the identified cis-elements within the promoter region. Moreover, altitude was found to have a positive effect on the PkPAL1 transcript levels, driving the expression of PkPAL1 abundantly. Based on docking analysis, we identified eight residues as potentially essential for substrate binding in PkPAL1.

An inducible NADPH-cytochrome P450 reductase from Picrorhiza kurrooa - an imperative redox partner of cytochrome P450 enzymes

Picrorhiza kurrooa synthesizes a large array of pharmacologically important monoterpenoid iridoid glycosides called picrosides. Although chemical profile and pharmacological activities of P. kurrooa have been extensively studied, limited attempts have been made to decipher the biosynthetic route and to identify the key regulatory genes involved in picroside biosynthesis. In the present study, NADPH-cytochrome P450 reductase, a key enzyme involved in electron transfer to cytochrome P450s was identified from P. kurrooa. The full length cDNA (2679 bp) contained an open reading frame of 2133 bp, corresponding to 710 amino acids. PkCPR was heterologously expressed in Escherichia coli and the kinetic parameters of the recombinant enzyme were determined. Specific activity, V max and K m of PkCPR were found to be 5.8 ± 0.05 μmol min(-1) mg(-1), 8.1 ± 0.12 μmol min(-1) mg(-1) and 7.8 μM, respectively. PkCPR was found to be spatially regulated at transcript level, being maximally expressed in leaf tissues. Altitude was found to have a positive effect on the picroside concentration and the picroside content positively correlated with the PkCPR transcript levels in samples collected at varied altitudes. Further, transcript profiling under methyl jasmonate, salicylic acid, 2,4-dicholorophenoxy acetic acid and UV-B elicitations displayed differential transcriptional regulation of PkCPR that fully corroborated with the identified cis-elements within the PkCPR promoter. Expression of PkCPR was inducible by UV-B and phytohormone elicitation, indicating that the PkCPR is possibly related to defence reactions, including biosynthesis of secondary metabolites. Present study is so far the only report of identification and functional characterization of CPR ortholog from P. kurrooa.

The GATA and SORLIP motifs in the 3-hydroxy-3-methylglutaryl-CoA reductase promoter of Picrorhiza kurrooa for the control of light-mediated expression

Light upregulates the expression of 3-hydroxy-3-methylglutaryl-CoA reductase (HMGR) in Picrorhiza kurrooa, an endangered medicinal herb. Upstream sequences of HMGR of P. kurrooa (PropkHMGR) were analyzed in relation to its role in light-mediated regulation of gene expression. GATA motif in PropkHMGR exhibited stronger DNA-protein interaction with the nuclear extract of dark-exposed plants in contrast to SORLIP that exhibited stronger binding with the nuclear extract of light-exposed plants. Analysis of PropkHMGR (PropkHMGR-D1, -1,059/-1) and its deletion fragments PropkHMGR-D2 (-825/-1), PropkHMGR-D3 (-651/-1), PropkHMGR-D4 (-452/-1), and PropkHMGR-D5 (-101/-1) in Arabidopsis thaliana showed PropkHMGR to regulate gene expression [β-glucuronidase (GUS) was used as a reporter gene] at all the developmental stages but only in actively dividing tissues, excluding anthers. Whereas, PropkHMGR-D2 regulated GUS expression in relatively older seedlings but the expression was observed only in shoot apical meristem, root tips, and anthers. PropkHMGR-mediated gene expression was higher in dark as compared to that in the light in Arabidopsis across four temperatures studied. As opposed to the results in P. kurrooa, GATA motifs exhibited DNA-protein interaction with nuclear extract of light-exposed plants of Arabidopsis. SORLIP motifs in Arabidopsis also exhibited DNA-protein interaction with nuclear extract of light-exposed plants as in P. kurrooa. Data showed that (1) PropkHMGR regulated light-mediated gene expression and (2) GATA motif exhibited an inverse relationship between strength of DNA-protein interaction and the gene expression whereas the relationship was species specific for SORLIP.