Multi-omics analysis of small RNA, transcriptome, and degradome to identify putative miRNAs linked to MeJA regulated and oridonin biosynthesis in Isodon rubescens.

Isodon rubescens has garnered much attention due to its anti-tumor or anti-cancer properties. However, little is known about the molecular mechanism of oridonin biosynthesis leveraging the regulatory network between small RNAs and mRNAs. In this study, the regulatory networks of miRNAs and targets were examined by combining mRNA, miRNA, and degradome. A total of 348 miRNAs, including 287 known miRNAs and 61 novel miRNAs, were identified. Among them, 51 miRNAs were significantly expressed, and 36 miRNAs responded to MeJA. A total of 3066 target genes were associated with 228 miRNAs via degradome sequencing. Multi-omics analysis demonstrated that 27 miRNA-mRNA pairs were speculated to be involved in MeJA regulation, and 36 miRNA-mRNA pairs were hypothesized to be involved in the genotype-dependence of I. rubescens. Furthermore, 151 and 7 miRNA-mRNA modules were likely engaged in oridonin biosynthesis as identified by psRNATarget and degradome sequencing, respectively. Some miRNA-mRNA modules were confirmed via RT-qPCR. Moreover, miRNAs targeting plant hormone signal transduction pathway genes were identified, such as miR156, miR167, miR393, and PC-3p-19822_242. Collectively, our results demonstrate for the first time that miRNAs are identified in I. rubescens, and laid a solid foundation for further research on the molecular mechanism of oridonin biosynthesis mediated by miRNA.

Identification and functional characterization of three cytochrome P450 genes for the abietane diterpenoid biosynthesis in Isodon lophanthoides.

MAIN CONCLUSION: We identify two ferruginol synthases and a 11-hydroxyferruginol synthase from a traditional Chinese medicinal herb Isodon lophanthoides and propose their involvement in two independent abietane diterpenoids biosynthetic pathways. Isodon lophanthoides is a traditional Chinese medicinal herb rich in highly oxidized abietane-type diterpenoids. These compounds exhibit a wide range of pharmaceutical activities, yet the biosynthesis is barely known. Here, we describe the screening and functional characterization of P450s that oxidize the abietane skeleton abietatriene. We mainly focused on CYP76 family and identified 12 CYP76AHs by mining the RNA-seq data of I. lophanthoides. Among the 12 CYP76AHs, 6 exhibited similar transcriptional expression features as upstream diterpene synthases, including root or leaf-preferential expression pattern and highly MeJA inducibility. These six P450s were considered as first-tier candidates and functionally characterized in yeast and plant cells. In yeast assays showed that both CYP76AH42 and CYP76AH43 were ferruginol synthases hydroxylating the C12 position of abietatriene, whereas CYP76AH46 was characterized as a 11-hydroxyferruginol synthase which catalyzes two successive oxidations at C12 and C11 of abietatriene. Heterologous expression of three CYP76AHs in Nicotiana benthamiana resulted in the formation of ferruginol. qPCR analysis showed CYP76AH42 and CYP76AH43 were mainly expressed in the root, which was consistent with the distribution of ferruginol in the root periderms. CYP76AH46 was primarily expressed in the leaves where barely ferruginol or 11-hydroxyferruginol was detected. In addition to distinct organ-specific expression pattern, three CYP76AHs exhibited different genomic structures (w or w/o introns), low protein sequence identities (51-63%) and were placed in separate subclades in the phylogenetic tree. These results suggest that the identified CYP76AHs may be involved in at least two independent abietane biosynthetic pathways in the aerial and underground parts of I. lophanthoides.

A chromosome-level genome assembly reveals that tandem-duplicated CYP706V oxidase genes control oridonin biosynthesis in the shoot apex of Isodon rubescens.

The ent-kaurenoids (e.g., oridonin and enmein) from the Isodon genus (Lamiaceae) are one class of diterpenoids with rich structural diversity and intriguing pharmaceutical activity. In contrast to the well-established gibberellin pathway, oxidative modifications diversifying the ent-kaurene skeleton in Isodon have remained undetermined for half a century. Here we report a chromosome-level genome assembly of I. rubescens, a well-recognized oridonin producer long favored by Asian people as a traditional herb with antitumor effects. The shoot apex was confirmed to be the actual region actively producing ent-kaurene diterpenoids. Through comparative genomics and phylogenetic analyses, we discovered a cluster of tandem-duplicated CYP706V oxygenase-encoding genes located on an ancient genomic block widely distributed in eudicots, whereas almost exclusively emerged in Isodon plants. In the shoot apex, IrCYP706V2 and IrCYP706V7 oxidized the ent-kaurene core in the initial stage of oridonin biosynthesis. Loss of CYP706Vs in other Lamiaceae plants offered an explanation for the specific kaurenoid production in Isodon plants. Moreover, we found that the Isodon genomes encode multiple diterpenoid synthases that are potentially involved in generating diterpenoid diversity. These findings provided new insights into the evolution of the lineage-specific diterpenoid pathway and laid a foundation for improving production of bioactive ent-kaurene-type diterpenoids by molecular breeding and synthetic biology approaches.

Comparative analysis of chloroplast genomes reveals phylogenetic relationships and intraspecific variation in the medicinal plant Isodon rubescens.

Isodon rubescens (Hemsley) H. Hara (Lamiaceae) is a traditional Chinese medicine plant that has been used to treat various human diseases and conditions such as inflammation, respiratory and gastrointestinal bacterial infections, and malignant tumors. However, the contents of the main active components of I. rubescens from different origins differ significantly, which greatly affected its quality. Therefore, a molecular method to identify and classify I. rubescens is needed. Here, we report the DNA sequence of the chloroplast genome of I. rubescens collected from Lushan, Henan province. The genome is 152,642 bp in length and has a conserved structure that includes a pair of IR regions (25,726 bp), a LSC region (83,527 bp) and a SSC region (17,663 bp). The chloroplast genome contains 113 unique genes, four rRNA genes, 30 tRNA genes, and 79 protein-coding genes, 23 of which contain introns. The protein-coding genes account for a total of 24,412 codons, and most of them are A/T biased usage. We identified 32 simple sequence repeats (SSRs) and 48 long repeats. Furthermore, we developed valuable chloroplast molecular resources by comparing chloroplast genomes from three Isodon species, and both mVISTA and DnaSP analyses showed that rps16-trnQ, trnS-trnG, and ndhC-trnM are candidate regions that will allow the identification of intraspecific differences within I. rubescens. Also 14 candidate fragments can be used to identify interspecific differences between species in Isodon. A phylogenetic analysis of the complete chloroplast genomes of 24 species in subfamily Nepetoideae was performed using the maximum likelihood method, and shows that I. rubescens clustered closer to I. serra than I. lophanthoides. Interestingly, our analysis showed that I. rubescens (MW018469.1) from Xianyang, Shaanxi Province (IR-X), is closer to I. serra than to the other two I. rubescens accessions. These results strongly indicate that intraspecific diversity is present in I. rubescens. Therefore, our results provide further insight into the phylogenetic relationships and interspecific diversity of species in the genus Isodon.

Discovery and Functional Characterization of a Diverse Diterpene Synthase Family in the Medicinal Herb Isodon lophanthoides Var. gerardiana.

Isodon lophanthoides var. gerardiana (Lamiaceae), also named xihuangcao, is a traditional Chinese medicinal herb that exhibits a broad range of pharmacological activities. Abietane-type diterpenoids are the characteristic constituents of I. lophanthoides, yet their biosynthesis has not been elucidated. Although the aerial parts are the most commonly used organs of I. lophanthoides, metabolite profiling by gas chromatography-mass spectrometry showed the underground parts also contain large amounts of labdane diterpenoids including abietatriene, miltiradiene and ferruginol, which are distinct from the 13-hydroxy-8(14)-abietene detected in the aerial parts. Comparative transcriptome analysis of root and leaf samples identified a diverse diterpene synthase family including 6 copalyl diphosphate synthase (IlCPS1-6) and 5 kaurene synthase-like (IlKSL1-5). Here we report the functional characterization of six of these enzymes using yeast heterologous expression system. Both IlCPS1 and IlCPS3 synthesized (+)-copalyl diphosphate (CPP), in combination with IlKSL1 resulted in miltiradiene, precursor of abietane-type diterpenoids, while coupling with IlKSL5 led to the formation of hydroxylated diterpene scaffold nezukol. Expression profiling and phylogenetic analysis further support the distinct evolutionary relationship and spatial distribution of IlCPS1 and IlCPS3. IlCPS2 converted GGPP into labda-7,13E-dien-15-ol diphosphate. IlCPS6 was identified as ent-CPS, indicating a role in gibberellin metabolism. We further identified a single residue that determined the water addition of nezukol synthase IlKSL5. Substitution of alanine 513 with isoleucine completely altered the product outcome from hydroxylated nezukol to isopimara-7,15-diene. Together, these findings elucidated the early steps of bioactive abietane-type diterpenoid biosynthesis in I. lophanthoides and the catalytic mechanism of nezukol synthase.

An alternative splicing alters the product outcome of a class I terpene synthase in Isodon rubescens.

The labdane-related diterpenoids are an important superfamily of natural products. Their structural diversity mainly depends on diterpene synthases, which generate the hydrocarbon skeletal structures. Isodon rubescens contains an expanded family of class I terpene synthases with different functions. Here we report a novel class I terpene synthase cDNA (IrKSL3a) with loss of 18 nucleotides compared with the reported cDNA sequence (IrKSL3). Inspection of IrKSL3 genomic sequence indicated that IrKSL3a and IrKSL3 transcripts may be generated by an alternative splicing event that utilizes different 3' splice site. In vitro assays showed that IrKSL3a produced isopimaradiene and miltiradiene, while IrKSL3 only produced miltiradiene. Protein sequence alignment found the six residues encoded by the alternative exon was unique to IrKSL3, which are 17 residues away from the conserved DDXXD motif. A deletion mutant of IrKSL3 showed that maintaining two residues within the six-amino acid is sufficient for miltiradiene production, while the other mutants lost nearly all enzymatic function. Our results illustrated how product outcomes can be changed by alternative splicing, and further gave an interesting example for studying the loop conformation in tuning product outcome in class I terpene synthase.

Diterpene synthases facilitating production of the kaurane skeleton of eriocalyxin B in the medicinal plant Isodon eriocalyx.

The Isodon plants (Lamiaceae) have been used in traditional Chinese medicine to alleviate sufferings from inflammations and cancers. This feature has been attributed to the presence of pharmacologically active ent-kaurane diterpenoids such as eriocalyxin B and oridonin. The Isodon eriocalyx (Dunn) Kudô species native to southwest China can accumulate a particularly high content of ent-kaurane diterpenoids (∼1.5% w/w of dried leaves). We previously identified diterpene synthases IeCPS1 and IeCPS2 as ent-copalyl diphosphate synthases (ent-CPS) potentially involved in Isodon ent-kaurane diterpenoids biosynthesis. In this study, analysis of RNA-seq transcriptome of the I. eriocalyx plant revealed three other diterpene synthase genes (IeCPS3, IeKS1, and IeKSL1). Their functional characterization through coupled in vitro enzyme assays has confirmed that IeCPS3 is an ent-CPS specifically producing ent-copalyl diphosphate (ent-CPP). IeKS1 accepted ent-CPP to produce exclusively ent-kaurene and may thus be defined as an ent-kaurene synthase (ent-KS). When IeKSL1 was combined with IeCPS2 or IeCPS3, no product was detected. Based on tissue-specific expression and metabolic localization studies, the IeCPS3 and IeKS1 transcripts were significantly accumulated in leaves where the ent-kaurane diterpenoid eriocalyxin B dominates, whereas weak expression of both were observed in germinating seeds in which gibberellin biosynthetic pathway is normally active. Our findings suggest that both IeCPS3 and IeKS1 possess dual roles in general (gibberellins) and specialized diterpenoid metabolism, such as that of the Isodon ent-kaurane diterpenoids.

Comparative transcriptome analysis of roots, stems and leaves of Isodon amethystoides reveals candidate genes involved in Wangzaozins biosynthesis.

BACKGROUND: Isodon amethystoides (Ben-th) Cy Wu et Hsuan is an important traditional medicinal plant endowed with pharmacological properties effective in the treatment of various diseases, including pulmonary tuberculosis. The tetracyclic diterpenoids, Wangzaozins (Wangzaozin A, glaucocalyxin A, glaucocalyxin B), are the major bioactive compounds of I. amethystoides. However, the molecular information about the biosynthesis of these compounds still remains unclear. RESULTS: An examination of the accumulated levels of Wangzaozins in I. amethystoides revealed considerable variations in the root, stem, and leaf tissues of this plant, indicating possible differences in metabolite biosynthesis and accumulation among various tissues. To better elucidate the tetracyclic diterpenoid biosynthesis pathway, we generated transcriptome sequences from the root, stem, and leaf tissues, and performed de novo sequence assembly, yielding 230,974 transcripts and 114,488 unigenes, with average N50 lengths of 1914 and 1241 bp, respectively. Putative functions could be assigned to 73,693 transcripts (31.9%) based on BLAST searches against annotation databases, including GO, KEGG, Swiss-Prot, NR, and Pfam. Moreover, the candidate genes involving in the diterpenoid biosynthesis, such as CPS, KSL, were also analyzed. The expression profiles of eight transcripts, involving the tetracyclic diterpenoid biosynthesis, were validated in different I. amethystoides tissues by qRT-PCR, unraveling the gene expression profile of the pathway. The differential expressions of ISPD, ISPF and ISPH (MEP pathway), and IaCPS and IaKSL (diterpenoid pathway) candidate genes in leaves and roots, may contribute to the high accumulation of Wangzaozins in I. amethystoides leaves. CONCLUSION: The genomic dataset and analyses reported here lay the foundations for further research on this important medicinal plant.

Biosynthesis of the oxygenated diterpene nezukol in the medicinal plant Isodon rubescens is catalyzed by a pair of diterpene synthases.

Plants produce an immense diversity of natural products (i.e. secondary or specialized metabolites) that offer a rich source of known and potentially new pharmaceuticals and other desirable bioproducts. The Traditional Chinese Medicinal plant Isodon rubescens (Lamiaceae) contains an array of bioactive labdane-related diterpenoid natural products. Of these, the ent-kauranoid oridonin is the most prominent specialized metabolite that has been extensively studied for its potent antimicrobial and anticancer efficacy. Mining of a previously established transcriptome of I. rubescens leaf tissue identified seven diterpene synthase (diTPSs) candidates. Here we report the functional characterization of four I. rubescens diTPSs. IrTPS5 and IrTPS3 were identified as an ent-copalyl diphosphate (CPP) synthase and a (+)-CPP synthase, respectively. Distinct transcript abundance of IrTPS5 and the predicted ent-CPP synthase IrTPS1 suggested a role of IrTPS5 in specialized ent-kaurene metabolism possibly en route to oridonin. Nicotiana benthamiana co-expression assays demonstrated that IrTPS4 functions sequentially with IrTPS3 to form miltiradiene. In addition, IrTPS2 converted the IrTPS3 product (+)-CPP into the hydroxylated tricyclic diterpene nezukol not previously identified in I. rubescens. Metabolite profiling verified the presence of nezukol in I. rubescens leaf tissue. The proposed IrTPS2-catalyzed reaction mechanism proceeds via the common ionization of the diphosphate group of (+)-CPP, followed by formation of an intermediary pimar-15-en-8-yl+ carbocation and neutralization of the carbocation by water capture at C-8 to yield nezukol, as confirmed by nuclear magnetic resonance (NMR) analysis. Oxygenation activity is rare for the family of class I diTPSs and offers new catalysts for developing metabolic engineering platforms to produce a broader spectrum of bioactive diterpenoid natural products.

Functional Diversification of Kaurene Synthase-Like Genes in Isodon rubescens.

Ent-kaurene diterpenoids are the largest group of known Isodon diterpenoids. Among them, oridonin is accumulated in the leaves, and is the most frequently studied compound because of its antitumor and antibacterial activities. We have identified five copalyl diphosphate synthase (CPS) and six kaurene synthase-like (KSL) genes by transcriptome profiling of Isodon rubescens leaves. An in vitro assay assigns ten of them to five different diterpene biosynthesis pathways, except IrCPS3 that has a mutation in the catalytic motif. The Lamiaceae-specific clade genes (IrCPS1 and IrCPS2) synthesize the intermediate copalyl diphosphate (normal-CPP), while IrCPS4 and IrCPS5 synthesize the intermediate ent-copalyl diphosphate (ent-CPP). IrKSL2, IrKSL4, and IrKSL5 react with ent-CPP to produce an ent-isopimaradiene-like compound, ent-atiserene and ent-kaurene, respectively. Correspondingly, the Lamiaceae-specific clade genes IrKSL1 or IrKSL3 combined with normal-CPP led to the formation of miltiradiene. The compound then underwent aromatization and oxidization with a cytochrome P450 forming two related compounds, abietatriene and ferruginol, which were detected in the root bark. IrKSL6 reacts with normal-CPP to produce isopimaradiene. IrKSL3 and IrKSL6 have the γßα tridomain structure, as these proteins tend to possess the bidomain structure of IrKSL1, highlighting the evolutionary history of KSL gene domain loss and further elucidating chemical diversity evolution from a macroevolutionary stance in Lamiaceae.

Phylogeny and historical biogeography of Isodon (Lamiaceae): rapid radiation in south-west China and Miocene overland dispersal into Africa.

Rapid organismal radiations occurring on the Qinghai-Tibetan Plateau (QTP) and the mechanisms underlying Asia-Africa intercontinental disjunctions have both attracted much attention from evolutionary biologists. Here we use the genus Isodon (Lamiaceae), a primarily East Asian lineage with disjunct species in central and southern Africa, as a case study to shed light upon these processes. The molecular phylogeny and biogeographic history of Isodon were reconstructed using sequences of three plastid markers, the nuclear ribosomal internal transcribed spacer (nrITS), and a low-copy nuclear gene (LEAFY intron II). The evolution of chromosome numbers in this genus was also investigated using probabilistic models. Our results support a monophyletic Isodon that includes the two disjunct African species, both of which likely formed through allopolyploidy. An overland migration from Asia to Africa through Arabia during the early Miocene is proposed as the most likely explanation for the present disjunct distribution of Isodon. The opening of the Red Sea in the middle Miocene may appear to have had a major role in disrupting floristic exchange between Asia and Africa. In addition, a rapid radiation of Isodon was suggested to occur in the late Miocene. It corresponds with one of the major uplifts of the QTP and subsequent aridification events. Our results support the hypothesis that geological and climatic events play important roles in driving biological diversification of organisms distributed in the QTP area.

Selection for high oridonin yield in the Chinese medicinal plant Isodon (Lamiaceae) using a combined phylogenetics and population genetics approach.

Oridonin is a diterpenoid with anti-cancer activity that occurs in the Chinese medicinal plant Isodon rubescens and some related species. While the bioactivity of oridonin has been well studied, the extent of natural variation in the production of this compound is poorly known. This study characterizes natural variation in oridonin production in order to guide selection of populations of Isodon with highest oridonin yield. Different populations of I. rubescens and related species were collected in China, and their offspring were grown in a greenhouse. Samples were examined for oridonin content, genotyped using 11 microsatellites, and representatives were sequenced for three phylogenetic markers (ITS, rps16, trnL-trnF). Oridonin production was mapped on a molecular phylogeny of the genus Isodon using samples from each population as well as previously published Genbank sequences. Oridonin has been reported in 12 out of 74 species of Isodon examined for diterpenoids, and the phylogeny indicates that oridonin production has arisen at least three times in the genus. Oridonin production was surprisingly consistent between wild-collected parents and greenhouse-grown offspring, despite evidence of gene flow between oridonin-producing and non-producing populations of Isodon. Additionally, microsatellite genetic distance between individuals was significantly correlated with chemical distance in both parents and offspring. Neither heritability nor correlation with genetic distance were significant when the comparison was restricted to only populations of I. rubescens, but this result should be corroborated using additional samples. Based on these results, future screening of Isodon populations for oridonin yield should initially prioritize a broad survey of all species known to produce oridonin, rather than focusing on multiple populations of one species, such as I. rubescens. Of the samples examined here, I. rubescens or I. japonicus from Henan province would provide the best source of oridonin.

IeCPS2 is potentially involved in the biosynthesis of pharmacologically active Isodon diterpenoids rather than gibberellin.

The traditional Chinese medicinal plant, Isodon L., is remarkably rich in pharmacologically active ent-kaurane diterpenoids of diverse carbon skeletons. In an effort to create a resource for gene discovery and elucidate the biosynthesis of Isodonent-kaurane diterpenoids, three cDNAs (named IeCPS1, IeCPS2 and IeCPS2a) were isolated putatively encoding copalyl diphosphate synthases from Isodoneriocalyx leaves. Recombinant proteins of IeCPS1 and IeCPS2 were expressed, respectively, in Escherichia coli, and were shown to specifically convert geranylgeranyl diphosphate to copalyl diphosphate as demonstrated by GC-MS analyses. Based on tissue-specific expression and metabolic localization studies, the IeCPS2 transcripts were detected in young and mature leaves where the dominant ent-kaurane diterpenoid maoecrystal B accumulates, whereas no detectable expression of IeCPS2 was observed in germinating seeds where the gibberellin biosynthetic pathway is usually active. In addition, no evidence for maoecrystal B was found in germinating seeds. On the other hand, IeCPS1 transcripts significantly accumulated in germinating seeds as well as in leaves. The biochemical and molecular genetic evidence thus indicated that IeCPS2 is a copalyl diphosphate synthase potentially involved in the biosynthesis of Isodon diterpenoids in leaves, while IeCPS1 is more probably relevant to gibberellin formation and may, in addition, participate in Isodonent-kaurane diterpenoid production.

[Identification of herba rabdosiae serrae from different plant resources by RAPD method].

OBJECTIVE: To identify the different plant resources of Herba Rabdosiae Serrae by using Random Amplified Polymorphic DNA ( RAPD ) Analysis. METHODS: The mini spin columns were used to extract the genomic DNA from five different plants of Herba Rabdosiae Serrae. With the DNA extracted from these plants as template,the 85 oligo nucleic acids (10 bp)as random primers,the polymer chain reaction (PCR) was done and the results were analysed by electro-pharoses. RESULTS: 12 primers were selected with polymorphism and 7 of them showed good polymorphism in RAPD map. CONCLUSION: RAPD method can be used to identify the plant resources of Herba Rabdosiae Serrae.

Development of microsatellite markers for the medicinal plant Isodon rubescens (Lamiaceae) and related species.

PREMISE OF THE STUDY: Microsatellite markers were developed for the medicinal plant Isodon rubescens to investigate genetic variability of the species and, in future studies, to assess its relation to the content of pharmacologically active chemicals produced by the plant. METHODS AND RESULTS: Eleven primer pairs were identified and tested in multiple populations of I. rubescens and related species (I. henryi, I. enanderianus, I. lophanthoides) from the People's Republic of China. The primers amplified dinucleotide repeats and had between two and 12 alleles per locus in a given population. CONCLUSIONS: These microsatellites will be useful for understanding patterns of phytochemical variation in I. rubescens and have the potential to be applied to research on evolutionary processes in other species of the genus Isodon.