Analysis of the transcriptome of Isodon rubescens and key enzymes involved in terpenoid biosynthesis

5-epi-ent-Kaurane diterpenoids from the aerial parts of Isodon eriocalyx and their anti-atherosclerotic potential

The phytochemical investigation of the EtOAc extract from the aerial parts of Isodon eriocalyx afforded seventeen diterpenoids, including eight undescribed compounds. Eriocalyxins H-L have unique structural characteristics featuring a 5-epi-ent-kaurane diterpenoid scaffold with eriocalyxins H-K also possess an unusual 6,11-epoxyspiro-lactone ring while eriocalyxin L, a 1,7:3,20-diepoxy-ent kaurene, features an 1,7-oxygen linkage. The structures of these compounds were elucidated by spectroscopic data interpretation, and the absolute configurations of eriocalyxins H, I, L, and M were confirmed by single-crystal X-ray diffraction. The isolates were screened for their inhibitory activities against VCAM-1 and ICAM-1 at 5 μM. While eriocalyxin O, coetsoidin A and laxiflorin P were found to significantly inhibit both VCAM-1 and ICAM-1, 8 (17),13-ent-labdadien-15 → 16-lactone-19-oic acid displayed evidently inhibitory effect against ICAM-1.

Oridonin Induces Apoptosis in Esophageal Squamous Cell Carcinoma by Inhibiting Cytoskeletal Protein LASP1 and PDLIM1

Esophageal squamous cell carcinoma is a severe malignancy for its high mortality and poor prognosis. Mainstay chemotherapies cause serious side effects for their ways of inducing cell death. Oridonin is the main bioactive constituent from natural plants that has anticancer ability and weak side effects. The proteomics method is efficient to understand the anticancer mechanism. However, proteins identified by proteomics aimed at understanding oridonin's anticancer mechanism is seldom overlapped by different groups. This study used proteomics based on two-dimensional electrophoresis sodium dodecyl sulfate-polyacrylamide gel electrophoresis (2-DE SDS-PAGE) integrated with mass spectrometry and Gene Set Enrichment Analysis (GSEA) to understand the anticancer mechanism of oridonin on esophageal squamous cell carcinoma (ESCC). The results showed that oridonin induced ESCC cell death via apoptosis by decreasing the protein expression of LASP1 and PDLIM1.

The complete chloroplast genome of Isodon rubescens, a traditional Chinese herb

Isodon rubescens (Hemsley) H. Hara is a kind of traditional medical herb which can be used for cancer treatment. In this study, the complete chloroplast genome sequence of I. rubescens was assembled. Its complete circular chloroplast DNA length was 152,761 bp. The genome was made up of a large single-copy region of 83,655 bp, a small single-copy region of 17,660 bp, and a pair of inverted repeat regions of 25,723 bp. The genome totally encoded 129 genes, containing 85 protein-coding genes, 36 tRNA genes, and eight rRNA genes. Phylogenetic analysis indicated that I. rubescens had a close relationship with basil (Ocimum basilicum).

Cross-Species Annotation of Expressed Genes and Detection of Different Functional Gene Modules Between 10 Cold- and 10 Hot-Propertied Chinese Herbal Medicines

According to the traditional Chinese medicine (TCM) system, Chinese herbal medicines (HMs) can be divided into four categories: hot, warm, cold, and cool. A cool nature usually is categorized as a cold nature, and a warm nature is classified as a hot nature. However, the detectable characteristics of the gene expression profile associated with the cold and hot properties have not been studied. To address this question, a strategy for the cross-species annotation of conserved genes was established in the present study by using transcriptome data of 20 HMs with cold and hot properties. Functional enrichment analysis was performed on group-specific expressed genes inferred from the functional genome of the reference species (i.e., Arabidopsis). Results showed that metabolic pathways relevant to chrysoeriol, luteolin, paniculatin, and wogonin were enriched for cold-specific genes, and pathways of inositol, heptadecane, lauric acid, octanoic acid, hexadecanoic acid, and pentadecanoic acid were enriched for hot-specific genes. Six functional modules were identified in the HMs with the cold property: nucleotide biosynthetic process, peptidy-L-cysteine S-palmitoylation, lipid modification, base-excision repair, dipeptide transport, and response to endoplasmic reticulum stress. For the hot HMs, another six functional modules were identified: embryonic meristem development, embryonic pattern specification, axis specification, regulation of RNA polymerase II transcriptional preinitiation complex assembly, mitochondrial RNA modification, and cell redox homeostasis. The research provided a new insight into HMs’ cold and hot properties from the perspective of the gene expression profile of plants.

Development of EST-based SSR and SNP markers in Gastrodia elata (herbal medicine) by sequencing, de novo assembly and annotation of the transcriptome

Tianma (Gastrodia elata Blume) has unique biological characteristics and high medicinal value. The wild resource of G. elata is being overutilized and should be conserved as it is already included in the list of endangered species in China. The population size of cultivated G. elata is small because of domestication bottleneck. Therefore, it is of utmost importance to evolve high-quality varieties and conserve wild resources of G. elata. In this study, we sequenced tuber transcriptomes of three major cultivated sub-species of Gastrodia elata, namely G. elata BI. f. elata, G. elata Bl. f. glauca S. Chow, and G. elata Bl. f. Viridis, and obtained about 7.8G clean data. The assembled high-quality reads of three sub-species were clustered into 56,884 unigenes. Of these, 31,224 (54.89%), 25,733 (45.24%), 22,629 (39.78%), and 11,856 (20.84%) unigenes were annotated by Nr, Swiss-Port, Eukaryotic Ortholog Groups (KOG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases, respectively. Here, a total of 3766 EST-SSRs and 128,921 SNPs were identified from the unigenes. The results not only offer huge number of genes that were responsible for the growth, development, and metabolism of bioactive components, but also a large number of molecular markers were detected for future studies on the conservation genetics and molecular breeding of G. elata.

Comparative transcriptome analyses on terpenoids metabolism in field- and mountain-cultivated ginseng roots

BACKGROUND

There exist differences in morphological traits and phytochemical compositions between field- and mountain-cultivated Panax ginseng (FCG and MCG), which might be attributed to variations of terpenoids metabolism adapting to different growth conditions. The present work aims to uncover these variations.

RESULTS

Among 26,648 differentially expressed genes, 496 genes distributed in seven dominant terpenoids pathways were identified. Diterpenoids and triterpenoids biosynthesis genes were significantly higher-expressed in FCG root. Conversely, biosynthesis of carotenoids was significantly more active in MCG root. Additionally, terpenoids backbones, monoterpenoids, sesquiterpenoids, and terpenoid-quinones biosyntheses were neither obviously inclined. Our determination also revealed that there were more gibberellins and steroids accumulated in FCG root which might be responsible for its quick vegetative growth, and enriched abscisic acid and germacrenes as well as protopanaxatriol-type ginsenosides might be major causes of enhanced stress-resistance in MCG root.

CONCLUSIONS

The study firstly provided an overview of terpenoids metabolism in roots of FCG and MCG in elucidating the underlying mechanisms for their different morphological appearances and phytochemical compositions.

Identification, expression, and phylogenetic analyses of terpenoid biosynthesis-related genes in secondary xylem of loblolly pine (Pinus taeda L.) based on transcriptome analyses

Loblolly pine (Pinus taeda L.) is one of the most important species for oleoresin (a mixture of terpenoids) in South China. The high oleoresin content of loblolly pine is associated with resistance to bark beetles and other economic benefits. In this study, we conducted transcriptome analyses of loblolly pine secondary xylem to gain insight into the genes involved in terpenoid biosynthesis. A total of 372 unigenes were identified as being critical for oleoresin production, including genes for ATP-binding cassette (ABC) transporters, the cytochrome P450 (CYP) protein family, and terpenoid backbone biosynthesis enzymes. Six key genes involved in terpenoid biosynthetic pathways were selected for multiple sequence alignment, conserved motif prediction, and phylogenetic and expression profile analyses. The protein sequences of all six genes exhibited a higher degree of sequence conservation, and upstream genes were relatively more conserved than downstream genes in terpenoid biosynthetic pathways. The N-terminal regions of these sequences were less conserved than the C-terminal ends, as the N-terminals were quite diverse in both length and composition. The phylogenetic analyses revealed that most genes originated from gene duplication after species divergence, and partial genes exhibited incomplete lineage sorting. In addition, the expression profile analyses showed that all six genes exhibited high expression levels during the high-oleoresin-yielding phase.

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.

Molecular mechanism of seed dormancy release induced by fluridone compared with cod stratification in Notopterygium incisum

BACKGROUND

Notopterygium incisum is an important Chinese medicinal plant. Its mature seeds have underdeveloped embryos and are physiological dormant. We found the seeds with full developed embryos can germinate after treated by fluridone (FL), an inhibitor of abscisic acid (ABA). In order to understand the molecular mechanisms underlying seed dormancy release by FL, we compared the transcriptomic changes in dormancy release induced by two different methods, FL and cold stratification (CS) in N. incisum. We further analyzed the gene expression patterns involved in seed germination and dormancy using quantitative reverse-transcription PCR.

RESULTS

RNA-sequence analysis revealed more dramatic changes in the transcriptomes of FL than those in CS, particularly for genes involved in the biosynthesis and regulation of gibberellins (GAs) and ABA. The down-regulation of ABA biosynthesis genes and the dramatic up-regulation of NiCYP707As, an ABA catabolic gene, contributed to the reduced ABA levels in FL. The increased GA3 levels in CS-treated seeds were due to the up-regulation of NiGA3OX. Both NiABI5 (a positive ABA regulator) and NiGAI (a negative regulator of GA) were down-regulated in FL and CS. The upregulation of strigolactones (SLs; the metabolites with the same precursor as ABA) biosynthesis and regulatory genes in both FL- and CS-treated seeds indicates that SLs contribute positively to seed dormancy release in N. incisum.

CONCLUSIONS

Our results indicated that FL- and CS-seed dormancy release possibly depends on two totally different mechanisms: alleviation of the effects of ABA and potentiation of the effects of GA, respectively. However, NiABI5 and NiGAI probably function as common factors integrating the effects of ABA and GA on seed dormancy release.

Physiological and de novo transcriptome analysis of the fermentation mechanism of Cerasus sachalinensis roots in response to short-term waterlogging

Background

Cerasus sachalinensis is widely used in cool regions as a sweet cherry rootstock and is known for its sensitivity to soil waterlogging and waterlogging stress. However, the limited availability of Cerasus genomic resources has considerably restricted the exploration of its waterlogging response mechanism. To understand its reaction to short-term waterlogging, we analyzed the physiology and transcriptomes of C. sachalinensis roots in response to different waterlogging durations.

Results

In this study, 12,487 differentially expressed genes (DEGs) were identified from Cerasus sachalinensis roots under different waterlogging durations. Carbon metabolism and energy maintenance formed the first coping mechanism stage of C. sachalinensis in response to low oxygen conditions. Root energy processes, including root respiration and activities of the fermentation enzymes alcohol dehydrogenase, pyruvate decarboxylase, and lactate dehydrogenase, showed unique changes after 0 h, 3 h, 6 h, and 24 h of waterlogging exposure. Ribonucleic acid sequencing was used to analyze transcriptome changes in C. sachalinensis roots treated with 3 h, 6 h, and 24 h of waterlogging stress. After de novo assembly, 597,474 unigenes were recognized, of which 355,350 (59.47%) were annotated. To identify the most important pathways represented by DEGs, Gene Ontology and Kyoto Encyclopedia of Genes and Genomes databases were used to compare these genes. The first stage of root reaction to waterlogging stress was activation of carbohydrate metabolism to produce more glucose and maintain energy levels. At 3 h, the glycolytic and fermentation pathways were activated to maintain adenosine triphosphate production. At 24 h, pathways involved in the translation of proteins were activated to further assist the plant in tolerating waterlogging stress. These findings will facilitate a further understanding of the potential mechanisms of plant responses to waterlogging at physiological and transcriptome levels.

Conclusions

Carbon metabolism and energy maintenance formed the first coping mechanism C. sachalinensis in response to low oxygen conditions, and they may be responsible for its short-term waterlogging response. Our study not only provides the assessment of genomic resources of Cerasus but also paves the way for probing the metabolic and molecular mechanisms underlying the short-term waterlogging response in C. sachalinensis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-4055-1) contains supplementary material, which is available to authorized users.

De Novo RNA Sequencing and Transcriptome Analysis of Monascus purpureus and Analysis of Key Genes Involved in Monacolin K Biosynthesis

Monascus purpureus is an important medicinal and edible microbial resource. To facilitate biological, biochemical, and molecular research on medicinal components of M. purpureus, we investigated the M. purpureus transcriptome by RNA sequencing (RNA-seq). An RNA-seq library was created using RNA extracted from a mixed sample of M. purpureus expressing different levels of monacolin K output. In total 29,713 unigenes were assembled from more than 60 million high-quality short reads. A BLAST search revealed hits for 21,331 unigenes in at least one of the protein or nucleotide databases used in this study. The 22,365 unigenes were categorized into 48 functional groups based on Gene Ontology classification. Owing to the economic and medicinal importance of M. purpureus, most studies on this organism have focused on the pharmacological activity of chemical components and the molecular function of genes involved in their biogenesis. In this study, we performed quantitative real-time PCR to detect the expression of genes related to monacolin K (mokA-mokI) at different phases (2, 5, 8, and 12 days) of M. purpureus M1 and M1-36. Our study found that mokF modulates monacolin K biogenesis in M. purpureus. Nine genes were suggested to be associated with the monacolin K biosynthesis. Studies on these genes could provide useful information on secondary metabolic processes in M. purpureus. These results indicate a detailed resource through genetic engineering of monacolin K biosynthesis in M. purpureus and related species.