Correction: Altered Heart Rate Variability During Rest in Schizophrenia: A State Marker.

[This corrects the article DOI: 10.7759/cureus.44145.].

Complete biosynthesis of QS-21 in engineered yeast.

QS-21 is a potent vaccine adjuvant and remains the only saponin-based adjuvant that has been clinically approved for use in humans1,2. However, owing to the complex structure of QS-21, its availability is limited. Today, the supply depends on laborious extraction from the Chilean soapbark tree or on low-yielding total chemical synthesis3,4. Here we demonstrate the complete biosynthesis of QS-21 and its precursors, as well as structural derivatives, in engineered yeast strains. The successful biosynthesis in yeast requires fine-tuning of the host's native pathway fluxes, as well as the functional and balanced expression of 38 heterologous enzymes. The required biosynthetic pathway spans seven enzyme families-a terpene synthase, P450s, nucleotide sugar synthases, glycosyltransferases, a coenzyme A ligase, acyl transferases and polyketide synthases-from six organisms, and mimics in yeast the subcellular compartmentalization of plants from the endoplasmic reticulum membrane to the cytosol. Finally, by taking advantage of the promiscuity of certain pathway enzymes, we produced structural analogues of QS-21 using this biosynthetic platform. This microbial production scheme will allow for the future establishment of a structure-activity relationship, and will thus enable the rational design of potent vaccine adjuvants.

Complete biosynthesis of the potent vaccine adjuvant QS-21.

QS-21 is a potent vaccine adjuvant currently sourced by extraction from the Chilean soapbark tree. It is a key component of human vaccines for shingles, malaria, coronavirus disease 2019 and others under development. The structure of QS-21 consists of a glycosylated triterpene scaffold coupled to a complex glycosylated 18-carbon acyl chain that is critical for immunostimulant activity. We previously identified the early pathway steps needed to make the triterpene glycoside scaffold; however, the biosynthetic route to the acyl chain, which is needed for stimulation of T cell proliferation, was unknown. Here, we report the biogenic origin of the acyl chain, characterize the series of enzymes required for its synthesis and addition and reconstitute the entire 20-step pathway in tobacco, thereby demonstrating the production of QS-21 in a heterologous expression system. This advance opens up unprecedented opportunities for bioengineering of vaccine adjuvants, investigating structure-activity relationships and understanding the mechanisms by which these compounds promote the human immune response.

Discovery of isoflavone phytoalexins in wheat reveals an alternative route to isoflavonoid biosynthesis.

Isoflavones are a group of phenolic compounds mostly restricted to plants of the legume family, where they mediate important interactions with plant-associated microbes, including in defense from pathogens and in nodulation. Their well-studied health promoting attributes have made them a prime target for metabolic engineering, both for bioproduction of isoflavones as high-value molecules, and in biofortification of food crops. A key gene in their biosynthesis, isoflavone synthase, was identified in legumes over two decades ago, but little is known about formation of isoflavones outside of this family. Here we identify a specialized wheat-specific isoflavone synthase, TaCYP71F53, which catalyzes a different reaction from the leguminous isoflavone synthases, thus revealing an alternative path to isoflavonoid biosynthesis and providing a non-transgenic route for engineering isoflavone production in wheat. TaCYP71F53 forms part of a biosynthetic gene cluster that produces a naringenin-derived O-methylated isoflavone, 5-hydroxy-2',4',7-trimethoxyisoflavone, triticein. Pathogen-induced production and in vitro antimicrobial activity of triticein suggest a defense-related role for this molecule in wheat. Genomic and metabolic analyses of wheat ancestral grasses further show that the triticein gene cluster was introduced into domesticated emmer wheat through natural hybridization ~9000 years ago, and encodes a pathogen-responsive metabolic pathway that is conserved in modern bread wheat varieties.

Altered Heart Rate Variability During Rest in Schizophrenia: A State Marker.

BACKGROUND: Autonomic nervous system (ANS) imbalance has been reported in a number of psychiatric disorders such as depression, schizophrenia, panic disorder, etc. Autonomic dysfunction in schizophrenia has been associated with the symptoms and manifestation of psychosis. Heart rate variability (HRV) as a tool has been widely used to assess ANS activity and the effect of disease on the sympathovagal balance. Therefore, in the present study, HRV derived from electrocardiogram (ECG) lead II at rest was investigated in order to understand the changes in frequency domain measures in patients with schizophrenia and their first-degree relatives compared to healthy controls. METHODS: Twenty-five patients with schizophrenia, 24 first-degree relatives of patients, and 24 healthy controls (Diagnostic and Statistical Manual of Mental Disorders (DSM)-5; 18-45 years) were included in the study. HRV of the subjects was measured after five minutes of rest. ECG lead II was recorded for five minutes and HRV was analysed in the frequency domain: low frequency (LF), high frequency (HF), total power, and LF/HF ratio. HRV parameters and heart rate were statistically analysed for group comparisons using general linear model multivariate analysis. RESULTS: Patients had significantly higher minimum heart rate and lower HF (normalized units (nu)) compared to their first-degree relatives. A trend was observed in HF (nu) with the lowest in patients followed by healthy controls and first-degree relatives and LF/HF ratio was the highest in patients followed by healthy controls and first-degree relatives, although not statistically significant. No significant difference was found between first-degree relatives and healthy controls. CONCLUSION: The alteration of HRV in schizophrenia could be attributed to reduction in vagal tone and sympathetic dominance, which in turn could serve as state markers of schizophrenia.

Elucidation of the pathway for biosynthesis of saponin adjuvants from the soapbark tree.

The Chilean soapbark tree (Quillaja saponaria) produces soap-like molecules called QS saponins that are important vaccine adjuvants. These highly valuable compounds are sourced by extraction from the bark, and their biosynthetic pathway is unknown. Here, we sequenced the Q. saponaria genome. Through genome mining and combinatorial expression in tobacco, we identified 16 pathway enzymes that together enable the production of advanced QS pathway intermediates that represent a bridgehead for adjuvant bioengineering. We further identified the enzymes needed to make QS-7, a saponin with excellent therapeutic properties and low toxicity that is present in low abundance in Q. saponaria bark extract. Our results enable the production of Q. saponaria vaccine adjuvants in tobacco and open the way for new routes to access and engineer natural and new-to-nature immunostimulants.

Biosynthesis of saponin defensive compounds in sea cucumbers.

Soft-bodied slow-moving sea creatures such as sea stars and sea cucumbers lack an adaptive immune system and have instead evolved the ability to make specialized protective chemicals (glycosylated steroids and triterpenes) as part of their innate immune system. This raises the intriguing question of how these biosynthetic pathways have evolved. Sea star saponins are steroidal, while those of the sea cucumber are triterpenoid. Sterol biosynthesis in animals involves cyclization of 2,3-oxidosqualene to lanosterol by the oxidosqualene cyclase (OSC) enzyme lanosterol synthase (LSS). Here we show that sea cucumbers lack LSS and instead have two divergent OSCs that produce triterpene saponins and that are likely to have evolved from an ancestral LSS by gene duplication and neofunctionalization. We further show that sea cucumbers make alternate sterols that confer protection against self-poisoning by their own saponins. Collectively, these events have enabled sea cucumbers to evolve the ability to produce saponins and saponin-resistant sterols concomitantly.

Pathogen-induced biosynthetic pathways encode defense-related molecules in bread wheat.

Wheat is a widely grown food crop that suffers major yield losses due to attack by pests and pathogens. A better understanding of biotic stress responses in wheat is thus of major importance. The recently assembled bread wheat genome coupled with extensive transcriptomic resources provides unprecedented new opportunities to investigate responses to pathogen challenge. Here, we analyze gene coexpression networks to identify modules showing consistent induction in response to pathogen exposure. Within the top pathogen-induced modules, we identify multiple clusters of physically adjacent genes that correspond to six pathogen-induced biosynthetic pathways that share a common regulatory network. Functional analysis reveals that these pathways, all of which are encoded by biosynthetic gene clusters, produce various different classes of compounds­namely, flavonoids, diterpenes, and triterpenes, including the defense-related compound ellarinacin. Through comparative genomics, we also identify associations with the known rice phytoalexins momilactones, as well as with a defense-related gene cluster in the grass model plant Brachypodium distachyon. Our results significantly advance the understanding of chemical defenses in wheat and open up avenues for enhancing disease resistance in this agriculturally important crop. They also exemplify the power of transcriptional networks to discover the biosynthesis of chemical defenses in plants with large, complex genomes.

UGT86C11 is a novel plant UDP-glycosyltransferase involved in labdane diterpene biosynthesis.

Glycosyltransferases constitute a large family of enzymes across all domains of life, but knowledge of their biochemical function remains largely incomplete, particularly in the context of plant specialized metabolism. The labdane diterpenes represent a large class of phytochemicals with many pharmacological benefits, such as anti-inflammatory, hepatoprotective, and anticarcinogenic. The medicinal plant kalmegh (Andrographis paniculata) produces bioactive labdane diterpenes; notably, the C19-hydroxyl diterpene (andrograpanin) is predominantly found as C19-O-glucoside (neoandrographolide), whereas diterpenes having additional hydroxylation(s) at C3 (14-deoxy-11,12-didehydroandrographolide) or C3 and C14 (andrographolide) are primarily detected as aglycones, signifying scaffold-selective C19-O-glucosylation of diterpenes in planta. Here, we analyzed UDP-glycosyltransferase (UGT) activity and diterpene levels across various developmental stages and tissues and found an apparent correlation of UGT activity with the spatiotemporal accumulation of neoandrographolide, the major diterpene C19-O-glucoside. The biochemical analysis of recombinant UGTs preferentially expressed in neoandrographolide-accumulating tissues identified a previously uncharacterized UGT86 member (ApUGT12/UGT86C11) that catalyzes C19-O-glucosylation of diterpenes with strict scaffold selectivity. ApUGT12 localized to the cytoplasm and catalyzed diterpene C19-O-glucosylation in planta. The substrate selectivity demonstrated by the recombinant ApUGT12 expressed in plant and bacterium hosts was comparable to native UGT activity. Recombinant ApUGT12 showed significantly higher catalytic efficiency using andrograpanin compared with 14-deoxy-11,12-didehydroandrographolide and trivial activity using andrographolide. Moreover, ApUGT12 silencing in plants led to a drastic reduction in neoandrographolide content and increased levels of andrograpanin. These data suggest the involvement of ApUGT12 in scaffold-selective C19-O-glucosylation of labdane diterpenes in plants. This knowledge of UGT86 function might help in developing plant chemotypes and synthesis of pharmacologically relevant diterpenes.

Prospective Randomized Control Trial Comparing Effect of Dexamethasone Versus Fentanyl as Adjuvants to Intrathecal Bupivacaine for Orthopedic Surgery.

Introduction Spinal anesthesia is the most consistent block for lower limb orthopedic surgeries. We conducted this randomized prospective study to evaluate comparative efficacy of intrathecal dexamethasone with fentanyl and normal saline as adjuvants to hyperbaric bupivacaine in spinal anesthesia administered to patients scheduled for lower limb orthopedic surgery. Materials and methods 105 patients scheduled for lower limb orthopedic surgeries under spinal anesthesia were included in this clinical trial. After randomization, patients received an intrathecal injection of hyperbaric bupivacaine (12.5 mg) with 4 mg of dexamethasone in group I, hyperbaric bupivacaine (12.5 mg) with 25 ug fentanyl with 0.5 ml of normal saline in group II and hyperbaric bupivacaine (12.5 mg) with normal saline (1 ml) in group III, so as to make volume of drug equal in all three groups. The observer evaluated the sensory and motor blocks and other parameters like time to self-void, stay in post-anesthesia care unit (PACU) and complications.  Results The total duration of sensory blockade was found to be 311.43, 197.86 and 115.29 minutes and motor blockade of 223.43, 163.86 and 83.0 minutes in groups I, II and III respectively. The PACU stay was 233.14, 173.86 and 93.00 minutes in groups I, II and III, respectively. The average time to self-void was 400.00, 315.29 and 203.00 in three groups, respectively. Conclusion Intrathecal dexamethasone seemed to be an effective adjuvant to spinal bupivacaine as it prolongs the duration of analgesia, stable hemodynamic profile with minimal side effects. Further studies are required to evaluate the optimum dose and long-term safety of intrathecal dexamethasone.

Virus-Induced Gene Silencing in Sweet Basil (Ocimum basilicum).

Virus-induced gene silencing (VIGS) is a powerful reverse genetic tool for rapid functional analysis of plant genes. Over the last decade, VIGS has been widely used for conducting rapid gene knockdown experiment in plants and played a crucial role in advancing applied and basic research in plant science. VIGS was studied extensively in model plants Arabidopsis and tobacco. Moreover, several non-model plants such as Papaver (Hileman et al., Plant J 44:334-341, 2005), Aquilegia (Gould and Kramer, Plant Methods 3:6, 2007), Catharanthus (Liscombe and O'Connor, Phytochemistry 72:1969-1977, 2011), Withania (Singh et al., Plant Biol J 13:1287-1299, 2015), and Ocimum (Misra et al., New Phytol 214:706-720, 2017) were also successfully explored. We have recently developed a robust protocol for VIGS in sweet basil (Ocimum basilicum). Sweet basil, a popular medicinal/aromatic herb, is being studied for the diversity of specialized metabolites produced in it.

Transcriptome analysis and functional characterization of oxidosqualene cyclases of the arjuna triterpene saponin pathway.

Arjuna (Terminalia arjuna) tree has been popular in Indian traditional medicine to treat cardiovascular ailments. The tree accumulates bioactive triterpene glycosides (saponins) and aglycones (sapogenins), in a tissue-preferential manner. Oleanane triterpenes/saponins (derived from ß-amyrin) with potential cardioprotective function predominantly accumulate in the bark. However, arjuna triterpene saponin pathway enzymes remain to be identified and biochemically characterized. Here, we employed a combined transcriptomics, metabolomics and biochemical approach to functionally define a suite of oxidosqualene cyclases (OSCs) that catalyzed key reactions towards triterpene scaffold diversification. De novo assembly of 131 millions Illumina NextSeq500 sequencing reads obtained from leaf and stem bark samples led to a total of 156,650 reference transcripts. Four distinct OSCs (TaOSC1-4) with 54-71 % sequence identities were identified and functionally characterized. TaOSC1, TaOSC3 and TaOSC4 were biochemically characterized as ß-amyrin synthase, cycloartenol synthase and lupeol synthase, respectively. However, TaOSC2 was found to be a multifunctional OSC producing both α-amyrin and ß-amyrin, but showed a preference for α-amyrin product. Both TaOSC1 and TaOSC2 produced ß-amyrin, the direct precursor for oleanane triterpene/saponin biosynthesis; but, TaOSC1 transcript expressed preferentially in bark, suggesting a major role of TaOSC1 in the biosynthesis of oleanane triterpenes/saponins in bark.

Oxidosqualene cyclase and CYP716 enzymes contribute to triterpene structural diversity in the medicinal tree banaba.

Pentacyclic triterpenes (PCTs) represent a major class of bioactive metabolites in banaba (Lagerstroemia speciosa) leaves; however, biosynthetic enzymes and their involvement in the temporal accumulation of PCTs remain to be studied. We use an integrated approach involving transcriptomics, metabolomics and gene function analysis to identify oxidosqualene cyclases (OSCs) and cytochrome P450 monooxygenases (P450s) that catalyzed sequential cyclization and oxidative reactions towards PCT scaffold diversification. Four monofunctional OSCs (LsOSC1,3-5) converted the triterpene precursor 2,3-oxidosqualene to either lupeol, ß-amyrin or cycloartenol, and a multifunctional LsOSC2 formed α-amyrin as a major product along with ß-amyrin. Two CYP716 family P450s (CYP716A265, CYP716A266) catalyzed C-28 oxidation of α-amyrin, ß-amyrin and lupeol to form ursolic acid, oleanolic acid and betulinic acid, respectively. However, CYP716C55 catalyzed C-2α hydroxylation of ursolic acid and oleanolic acid to produce corosolic acid and maslinic acid, respectively. Besides, combined transcript and metabolite analysis suggested major roles for the LsOSC2, CYP716A265 and CYP716C55 in determining leaf ursane and oleanane profiles. Combinatorial expression of OSCs and CYP716s in Saccharomyces cerevisiae and Nicotiana benthamiana led to PCT pathway reconstruction, signifying the utility of banaba enzymes for bioactive PCT production in alternate plant/microbial hosts that are more easily tractable than the tree species.

Two CYP716A subfamily cytochrome P450 monooxygenases of sweet basil play similar but nonredundant roles in ursane- and oleanane-type pentacyclic triterpene biosynthesis.

The medicinal plant sweet basil (Ocimum basilicum) accumulates bioactive ursane- and oleanane-type pentacyclic triterpenes (PCTs), ursolic acid and oleanolic acid, respectively, in a spatio-temporal manner; however, the biosynthetic enzymes and their contributions towards PCT biosynthesis remain to be elucidated. Two CYP716A subfamily cytochrome P450 monooxygenases (CYP716A252 and CYP716A253) are identified from a methyl jasmonate-responsive expression sequence tag collection and functionally characterized, employing yeast (Saccharomyces cerevisiae) expression platform and adapting virus-induced gene silencing (VIGS) in sweet basil. CYP716A252 and CYP716A253 catalyzed sequential three-step oxidation at the C-28 position of α-amyrin and ß-amyrin to produce ursolic acid and oleanolic acid, respectively. Although CYP716A253 was more efficient than CYP716A252 for amyrin C-28 oxidation in yeast, VIGS revealed essential roles for both of these CYP716As in constitutive biosynthesis of ursolic acid and oleanolic acid in sweet basil leaves. However, CYP716A253 played a major role in elicitor-induced biosynthesis of ursolic acid and oleanolic acid. Overall, the results suggest similar as well as distinct roles of CYP716A252 and CYP716A253 for the spatio-temporal biosynthesis of PCTs. CYP716A252 and CYP716A253 might be useful for the alternative and sustainable production of PCTs in microbial host, besides increasing plant metabolite content through genetic modification.

Use of Dexmedetomidine and Esmolol for Hypotension in Lumbar Spine Surgery.

BACKGROUND: The importance of decreasing bleeding in spine surgery is not only important to maintain the patient's hemodynamic balance but also allow a better view of the surgical field. OBJECTIVES: The current study aimed to compare dexmedetomidine and Esmolol™ as agents to induce hypotension in lumbar spine surgeries. PATIENTS AND METHODS: A total of 50 patients aged 20 to 65 years belonging to the American society of anaesthesiologist (ASA) class I - II scheduled for decompression and fixation of the lumbar spine were included and divided into two groups namely, Group I, who received Esmolol and group II, who received dexmedetomidine, intravenously. The patients were compared for intraoperative hemodynamic parameters, estimated blood loss, operation time, intraoperative analgesic (fentanyl) consumption, and total fall in haemoglobin (Hb) during the perioperative period. RESULTS: The study results showed that dexmedetomidine had lower (100.8 µg) fentanyl and sevoflurane consumption (1.2%), and less blood loss (278 mL) in comparison to the Esmolol group. CONCLUSIONS: Both dexmedetomidine and Esmolol can be used as agents to control hypotension in patients undergoing lumbar spine decompression and fixation surgery; the dexmedetomidine group, however, was associated with better intraoperative hemodynamic stability and reduced intraoperative analgesic and volatile anaesthetic requirement.

A thaumatin-like protein of Ocimum basilicum confers tolerance to fungal pathogen and abiotic stress in transgenic Arabidopsis.

Plant often responds to fungal pathogens by expressing a group of proteins known as pathogenesis-related proteins (PRs). The expression of PR is mediated through pathogen-induced signal-transduction pathways that are fine-tuned by phytohormones such as methyl jasmonate (MeJA). Here, we report functional characterization of an Ocimum basilicum PR5 family member (ObTLP1) that was identified from a MeJA-responsive expression sequence tag collection. ObTLP1 encodes a 226 amino acid polypeptide that showed sequence and structural similarities with a sweet-tasting protein thaumatin of Thaumatococcus danielli and also with a stress-responsive protein osmotin of Nicotiana tabacum. The expression of ObTLP1 in O. basilicum was found to be organ-preferential under unstressed condition, and responsive to biotic and abiotic stresses, and multiple phytohormone elicitations. Bacterially-expressed recombinant ObTLP1 inhibited mycelial growth of the phytopathogenic fungi, Scleretonia sclerotiorum and Botrytis cinerea; thereby, suggesting its antifungal activity. Ectopic expression of ObTLP1 in Arabidopsis led to enhanced tolerance to S. sclerotiorum and B. cinerea infections, and also to dehydration and salt stress. Moreover, induced expression of the defense marker genes suggested up-regulation of the defense-response pathways in ObTLP1-expressing Arabidopsis upon fungal challenge. Thus, ObTLP1 might be useful for providing tolerance to the fungal pathogens and abiotic stresses in crops.

Effect of Moderate Aerobic Exercise Training on Autonomic Functions and its Correlation with the Antioxidant Status.

Aerobic exercise is a proven measure to enhance the cardio-respiratory efficiency. This improvement isprimarily due to improvement in autonomic function where in there is increase in the parasympathetic function and decrease in sympathetic function. Exercise also affects many metabolic activities in our body and may be one of the factors which reduce the oxidative stress. The aim of the present study was to see the effect of moderate aerobic exercise training on autonomic function and its correlation with antioxidant status. 30 healthy volunteers in the age group of 18-22 years were screened. Autonomic function tests included activity (tone) and reactivity parameters. Antioxidant status was assessed by the level of malondialdehyde in plasma. We observed a significant change in SBP, DBP, LF nu, HF nu, delta value of DBP in CPT and MDA. Our findings are consistent with earlier findings that short duration physical training is known to reduce blood pressure and that there is a relationship between HF (in HRV) and training response Physical exercise also provides a favorable change in the biochemical parameters such as MDA.

Involvement of an ent-copalyl diphosphate synthase in tissue-specific accumulation of specialized diterpenes in Andrographis paniculata.

Ent-labdane-related diterpene (ent-LRD) specialized (i.e. secondary) metabolites of the medicinal plant kalmegh (Andrographis paniculata) have long been known for several pharmacological activities. However, our understanding of the ent-LRD biosynthetic pathway has remained largely incomplete. Since ent-LRDs accumulate in leaves, we carried out a comparative transcriptional analysis using leaf and root tissues, and identified 389 differentially expressed transcripts, including 223 transcripts that were preferentially expressed in leaf tissue. Analysis of the transcripts revealed various specialized metabolic pathways, including transcripts of the ent-LRD biosynthetic pathway. Two class II diterpene synthases (ApCPS1 and ApCPS2) along with one (ApCPS1') and two (ApCPS2' and ApCPS2″) transcriptional variants that were the outcomes of alternative splicing of the precursor mRNA and alternative transcriptional termination, respectively, were identified. ApCPS1 and ApCPS2 encode for 832- and 817-amino acids proteins, respectively, and are phylogenetically related to the dicotyledons ent-copalyl diphosphate synthases (ent-CPSs). The spatio-temporal patterns of ent-LRD metabolites accumulation and gene expression suggested a likely role for ApCPS1 in general (i.e. primary) metabolism, perhaps by providing precursor for the biosynthesis of phytohormone gibberellin (GA). However, ApCPS2 is potentially involved in tissue-specific accumulation of ent-LRD specialized metabolites. Bacterially expressed recombinant ApCPS2 catalyzed the conversion of (E,E,E)-geranylgeranyl diphosphate (GGPP), the general precursor of diterpenes to ent-copalyl diphosphate (ent-CPP), the precursor of ent-LRDs. Taken together, these results advance our understanding of the tissue-specific accumulation of specialized ent-LRDs of medicinal importance.

Andrographis paniculata transcriptome provides molecular insights into tissue-specific accumulation of medicinal diterpenes.

BACKGROUND: Kalmegh (Andrographis paniculata) has been widely exploited in traditional medicine for the treatment of infectious diseases and health disorders. Ent-labdane-related diterpene (ent-LRD) specialized (i.e., secondary) metabolites of kalmegh such as andrographolide, neoandrographolide and 14-deoxy-11,12-didehydroandrographolide, are known for variety of pharmacological activities. However, due to the lack of genomic and transcriptomic information, underlying molecular basis of ent-LRDs biosynthesis has remained largely unknown. To identify candidate genes of the ent-LRD biosynthetic pathway, we performed comparative transcriptome analysis using leaf and root tissues that differentially accumulate ent-LRDs. RESULTS: De novo assembly of Illumina HiSeq2000 platform-generated paired-end sequencing reads resulted into 69,011 leaf and 64,244 root transcripts which were assembled into a total of 84,628 unique transcripts. Annotation of these transcripts to the Uniprot, Kyoto Encyclopedia of Genes and Genomes (KEGG) and Carbohydrate-Active Enzymes (CAZy) databases identified candidate transcripts of the ent-LRD biosynthetic pathway. These included transcripts that encode enzymes of the plastidial 2C-methyl-D-erythritol-4-phosphate pathway which provides C5 isoprenoid precursors for the ent-LRDs biosynthesis, geranylgeranyl diphosphate synthase, class II diterpene synthase (diTPS), cytochrome P450 monooxygenase and glycosyltransferase. Three class II diTPSs (ApCPS1, ApCPS2 and ApCPS3) that showed distinct tissue-specific expression profiles and are phylogenetically related to the dicotyledon ent-copalyl diphosphate synthases, are identified. ApCPS1, ApCPS2 and ApCPS3 encode for 832-, 817- and 797- amino acids proteins of 55-63 % identity, respectively. Spatio-temporal patterns of transcripts and ent-LRDs accumulation are consistent with the involvement of ApCPS1 in general (i.e., primary) metabolism for the biosynthesis of phytohormone gibberellin, ApCPS2 in leaf specialized ent-LRDs biosynthesis and ApCPS3 in root diterpene biosynthesis. Moreover, simple sequence repeats (SSRs) that might assist in genotyping and developing specific chemotypes were identified in transcripts of the specialized metabolic pathways, including ent-LRDs. CONCLUSIONS: Comparative analysis of root and leaf transcriptomes disclosed novel genes of the ent-LRD biosynthetic pathway, including three class II diTPSs that showed discrete spatio-temporal expression patterns; thus, suggesting their participation into distinct diterpene metabolic pathways of kalmegh. Overall, these results will be useful in understanding molecular basis of the medicinal ent-LRDs biosynthesis and developing breeding strategies for improving their yields.