Complexity of gene paralogues resolved in biosynthetic pathway of hepatoprotective iridoid glycosides in a medicinal herb, Picrorhiza kurroa through differential NGS transcriptomes

A Transcriptomic and Metabolomic Study on the Biosynthesis of Iridoids in Phlomoides rotata from the Qinghai-Tibet Plateau

Phlomoides rotata is a traditional Chinese herbal medicine that grows in the Qinghai-Tibet Plateau region at a 3100-5000 m altitude. Iridoid compounds are the main active compounds of the P. rotata used as medical ingredients and display anti-inflammatory, analgesic, and hepatoprotective properties. To better understand the biological mechanisms of iridoid compounds in this species, we performed a comprehensive analysis of the transcriptome and metabolome of P. rotata leaves from four different regions (3540-4270 m). Global metabolome profiling detected 575 metabolites, and 455 differentially accumulated metabolites (DAMs) were detected in P. rotata leaves from the four regions. Eight major DAMs related to iridoid metabolism in P. rotata leaves were investigated: shanzhiside methyl ester, 8-epideoxyloganic acid, barlerin, shanzhiside, geniposide, agnuside, feretoside, and catalpin. In addition, five soil physical and chemical indicators in P. rotata rhizosphere soils were analyzed. Four significant positive correlations were observed between alkaline nitrogen and geniposide, exchangeable calcium and geniposide, available potassium and shanzhiside, and available phosphorus and shanzhiside methyl ester. The transcriptome data showed 12 P. rotata cDNA libraries with 74.46 Gb of clean data, which formed 29,833 unigenes. Moreover, 78.91% of the unigenes were annotated using the eight public databases. Forty-one candidate genes representing 23 enzymes involved in the biosynthesis of iridoid compounds were identified in P. rotata leaves. Moreover, the DXS1, IDI1, 8-HGO1, and G10H2 genes associated with iridoid biosynthesis were specifically expressed in P. rotata. The integration of transcriptome and metabolome analyses highlights the crucial role of soil physical and chemical indicators and major gene expression related to iridoid metabolism pathways in P. rotata from different areas. Our findings provide a theoretical foundation for exploring the molecular mechanisms underlying iridoid compound accumulation in P. rotata.

ABC transporters mined through comparative transcriptomics associate with organ-specific accumulation of picrosides in a medicinal herb, Picrorhiza kurroa

Picrorhiza kurroa Royle ex Benth is a valuable medicinal herb of North-Western Himalayas due to presence of two major bioactive compounds, picroside-I and picroside-II used in the preparation of several hepatoprotective herbal drugs. These compounds accumulate in stolons/rhizomes; however, biosynthesized in different organs, viz., picroside-I in shoots and picroside-II in roots. As of today, no information exists on what transporters are transporting these metabolites from shoots and roots to the final storage organ, stolon, which ultimately transforms into rhizome. The ATP-binding cassette (ABC) transporters are reported to transport majority of secondary metabolites, including terpenoids in plants, therefore, we mined P. kurroa transcriptomes to identify and shortlist potential candidates. A total of 99 ABC transporter-encoding transcripts were identified in 3 differential transcriptomes, PKSS (shoots), PKSTS (stolons), and PKSR (roots) of P. kurroa, based on in silico comparative analysis and transcript abundance. 15 of these transcripts were further validated for their association using qRT-PCR in shoots, roots and stolon tissues in P. kurroa accessions varying for picroside-I and picroside-II contents. Organ-specific expression analysis revealed that PkABCA1, PkABCG1, and PkABCB5 had comparatively elevated expression in shoots; PkABCB2 and PkABCC2 in roots; PkABCB3 and PkABCC1 in stolon tissues of P. kurroa. Co-expression network analysis using ABC genes as hubs further unravelled important interactions with additional components of biosynthetic machinery. Our study has provided leads, first to our knowledge as of today, on putative ABC transporters possibly involved in long distance and local transport of picrosides in P. kurroa organs, thus opening avenues for designing a suitable genetic intervention strategy.

Picrorhiza kurroa, Royle ex Benth:Traditional uses, phytopharmacology, and translational potential in therapy of fatty liver disease

Picrorhiza kurroa Royle ex Benth, Kutki (P.kurroa) is an important medicinal plant, traditionally recommended and used in Ayurveda for millennia, with certain cautions. There has been a significant revival of keen interest in its pharmacology, pharmacognosy, and phytochemistry for the last few decades. The evidence of its hepatoprotective activity, in experimental and clinical studies, accelerated the correlation of the specific phytochemical constituents of P.kurroa with precise pharmacological activities. Iridoid glycosides, particularly picrosides, emerged as the active molecules. For effective translation of traditional remedies into modern therapy, value addition by mechanistic understanding of molecular actions, drug targets, the degrees of efficacy and safety as well as convenient dosage forms is needed. Reverse pharmacology approach and phytopharmaceutical drug category facilitate such a translation. The present review illustrates how a potential translation of traditional practices of using P.kurroa into a phytochemically standardized, clinically targeted natural product for global unmet medical needs viz. Fatty liver disease can be attained.

Capturing acyltransferase(s) transforming final step in the biosynthesis of a major Iridoid Glycoside, (Picroside-II) in a Himalayan Medicinal Herb, Picrorhiza kurroa

BACKGROUND

Picrorhiza kurroa has been reported as an age-old ayurvedic hepato-protection to treat hepatic disorders due to the presence of iridoids such as picroside-II (P-II), picroside-I, and kutkoside. The acylation of catalpol and vanilloyl coenzyme A by acyltransferases (ATs) is critical step in P-II biosynthesis. Since accumulation of P-II occurs only in roots, rhizomes and stolons in comparison to leaves uprooting of this critically endangered herb has been the only source of this compound. Recently, we reported that P-II acylation likely happen in roots, while stolons serve as the vital P-II storage compartment. Therefore, developing an alternate engineered platform for P-II biosynthesis require identification of P-II specific AT/s.

METHODS AND RESULTS

In that direction, egg-NOG function annotated 815 ATs from de novo RNA sequencing of tissue culture based 'shoots-only' system and nursery grown shoots, roots, and stolons varying in P-II content, were cross-compared in silico to arrive at ATs sequences unique and/or common to stolons and roots. Verification for organ and accession-wise upregulation in gene expression of these ATs by qRT-PCR has shortlisted six putative 'P-II-forming' ATs. Further, six-frame translation, ab initio protein structure modelling and protein-ligand molecular docking of these ATs signified one MBOAT domain containing AT with preferential binding to the vanillic acid CoA thiol ester as well as with P-II, implying that this could be potential AT decorating final structure of P-II.

CONCLUSIONS

Organ-wise comparative transcriptome mining coupled with reverse transcription real time qRT-PCR and protein-ligand docking led to the identification of an acyltransferases, contributing to the final structure of P-II.