Comprehensive transcriptome analysis reveals novel genes involved in cardiac glycoside biosynthesis and mlncRNAs associated with secondary metabolism and stress response in Digitalis purpurea

The Ginkgo biloba microRNA160-ERF4 module participates in terpene trilactone biosynthesis

Terpene trilactones (TTLs) are important secondary metabolites in ginkgo (Ginkgo biloba); however, their biosynthesis gene regulatory network remains unclear. Here, we isolated a G. biloba ethylene response factor 4 (GbERF4) involved in TTL synthesis. Overexpression of GbERF4 in tobacco (Nicotiana tabacum) significantly increased terpenoid content and upregulated the expression of key enzyme genes (3-hydroxy-3-methylglutaryl-CoA reductase [HMGR], 3-hydroxy-3-methylglutaryl-CoA synthase [HMGS], 1-deoxy-D-xylulose-5-phosphate reductoisomerase [DXR], 1-deoxy-D-xylulose-5-phosphate synthase [DXS], acetyl-CoA C-acetyltransferase [AACT], and geranylgeranyl diphosphate synthase [GGPPS]) in the terpenoid pathway in tobacco, suggesting that GbERF4 functions in regulating the synthesis of terpenoids. The expression pattern analysis and previous microRNA (miRNA) sequencing showed that gb-miR160 negatively regulates the biosynthesis of TTLs. Transgenic experiments showed that overexpression of gb-miR160 could significantly inhibit the accumulation of terpenoids in tobacco. Targeted inhibition and dual-luciferase reporter assays confirmed that gb-miR160 targets and negatively regulates GbERF4. Transient overexpression of GbERF4 increased TTL content in G. biloba, and further transcriptome analysis revealed that DXS, HMGS, CYPs, and transcription factor genes were upregulated. In addition, yeast 1-hybrid and dual-luciferase reporter assays showed that GbERF4 could bind to the promoters of the HMGS1, AACT1, DXS1, levopimaradiene synthase (LPS2), and GGPPS2 genes in the TTL biosynthesis pathway and activate their expression. In summary, this study investigated the molecular mechanism of the gb-miR160-GbERF4 regulatory module in regulating the biosynthesis of TTLs. It provides information for enriching the understanding of the regulatory network of TTL biosynthesis and offers important gene resources for the genetic improvement of G. biloba with high contents of TTLs.

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.

Physiological and transcriptomic analyses of leaves from Gardenia jasminoides Ellis under waterlogging stress

Abstract Gardenia jasminoides Ellis is a Chinese herbal medicine with medicinal and economic value, but its mechanism of response to waterlogging stress remains unclear. In this study, the “double pots method” was used to simulate the waterlogging stress of Gardenia jasminoides Ellis to explore its physiological and transcriptomic response mechanism. We found no significant damage to Gardenia jasminoides Ellis membrane lipid during stress. POD played a vital antioxidant role, KEGG enrichment showed that secondary metabolites such as flavonoids might also play an antioxidant role, and PRO played a significant osmotic adjustment. Endogenous hormones regulate the Gardenia jasminoides Ellis's growth and development and play a role in signal transduction. Among them, light waterlogging stress is delayed. At the same time, there were 19631, 23693, and 15045 differentially expressed genes on the 5th, 10d, and 15d of Gardenia jasminoides Ellis under waterlogging stress. These genes were closely associated with the proteasome, endopeptidase, ribosome, MAPK signal transduction, and endogenous hormone signal transduction, plant-pathogen interaction and phenylpropanoid biosynthesis and other physiological and metabolic pathways, which regulate the turnover and transportation of protein, the reinforcement and adhesion of cell walls, the induction of stomatal closure, allergic reactions, defense reactions, leaf movements and others. It also can absorb ultraviolet rays to reduce the generation of oxygen free radicals, change the way of energy utilization and adjust the osmotic pressure of plant cells.

Multilayered regulation of secondary metabolism in medicinal plants

Medicinal plants represent a huge reservoir of secondary metabolites (SMs), substances with significant pharmaceutical and industrial potential. However, obtaining secondary metabolites remains a challenge due to their low-yield accumulation in medicinal plants; moreover, these secondary metabolites are produced through tightly coordinated pathways involving many spatiotemporally and environmentally regulated steps. The first regulatory layer involves a complex network of transcription factors; a second, more recently discovered layer of complexity in the regulation of SMs is epigenetic modification, such as DNA methylation, histone modification and small RNA-based mechanisms, which can jointly or separately influence secondary metabolites by regulating gene expression. Here, we summarize the findings in the fields of genetic and epigenetic regulation with a special emphasis on SMs in medicinal plants, providing a new perspective on the multiple layers of regulation of gene expression.

Identification of key genes involved in secondary metabolite biosynthesis in Digitalis purpurea

The medicinal plant Digitalis purpurea produces cardiac glycosides that are useful in the pharmaceutical industry. These bioactive compounds are in high demand due to ethnobotany's application to therapeutic procedures. Recent studies have investigated the role of integrative analysis of multi-omics data in understanding cellular metabolic status through systems metabolic engineering approach, as well as its application to genetically engineering metabolic pathways. In spite of numerous omics experiments, most molecular mechanisms involved in metabolic pathways biosynthesis in D. purpurea remain unclear. Using R Package Weighted Gene Co-expression Network Analysis, co-expression analysis was performed on the transcriptome and metabolome data. As a result of our study, we identified transcription factors, transcriptional regulators, protein kinases, transporters, non-coding RNAs, and hub genes that are involved in the production of secondary metabolites. Since jasmonates are involved in the biosynthesis of cardiac glycosides, the candidate genes for Scarecrow-Like Protein 14 (SCL14), Delta24-sterol reductase (DWF1), HYDRA1 (HYD1), and Jasmonate-ZIM domain3 (JAZ3) were validated under methyl jasmonate treatment (MeJA, 100 μM). Despite early induction of JAZ3, which affected downstream genes, it was dramatically suppressed after 48 hours. SCL14, which targets DWF1, and HYD1, which induces cholesterol and cardiac glycoside biosynthesis, were both promoted. The correlation between key genes and main metabolites and validation of expression patterns provide a unique insight into the biosynthesis mechanisms of cardiac glycosides in D. purpurea.

Integrated analysis of lncRNAs, mRNAs, and TFs to identify network modules underlying diterpenoid biosynthesis in Salvia miltiorrhiza

Long non-coding RNAs (lncRNAs) are transcripts of more than 200 nucleotides (nt) in length, with minimal or no protein-coding capacity. Increasing evidence indicates that lncRNAs play important roles in the regulation of gene expression including in the biosynthesis of secondary metabolites. Salvia miltiorrhiza Bunge is an important medicinal plant in China. Diterpenoid tanshinones are one of the main active components of S. miltiorrhiza. To better understand the role of lncRNAs in regulating diterpenoid biosynthesis in S. miltiorrhiza, we integrated analysis of lncRNAs, mRNAs, and transcription factors (TFs) to identify network modules underlying diterpenoid biosynthesis based on transcriptomic data. In transcriptomic data, we obtained 6,651 candidate lncRNAs, 46 diterpenoid biosynthetic pathway genes, and 11 TFs involved in diterpenoid biosynthesis. Combining the co-expression and genomic location analysis, we obtained 23 candidate lncRNA-mRNA/TF pairs that were both co-expressed and co-located. To further observe the expression patterns of these 23 candidate gene pairs, we analyzed the time-series expression of S. miltiorrhiza induced by methyl jasmonate (MeJA). The results showed that 19 genes were differentially expressed at least a time-point, and four lncRNAs, two mRNAs, and two TFs formed three lncRNA-mRNA and/or TF network modules. This study revealed the relationship among lncRNAs, mRNAs, and TFs and provided new insight into the regulation of the biosynthetic pathway of S. miltiorrhiza diterpenoids.

Ploidy Status, Nuclear DNA Content and Start Codon Targeted (SCoT) Genetic Homogeneity Assessment in Digitalis purpurea L., Regenerated In Vitro

Digitalis purpurea L. is a therapeutically important plant that synthesizes important cardiotonics such as digitoxin and digoxin. The present work reports a detailed and efficient propagation protocol for D. purpurea by optimizing various PGR concentrations in Murashige and Skoog (MS) medium. The genetic homogeneity of in vitro regenerants was assessed by the flow cytometric method (FCM) and Start Codon Targeted (SCoT) marker technique. Firstly, the seeds inoculated in full MS medium added with 0.5 mg/L GA₃ produced seedlings. Different parts such as hypocotyl, nodes, leaves and apical shoots were used as explants. The compact calli were obtained on BAP alone or in combinations with 2, 4-D/NAA. The hypocotyl-derived callus induced somatic embryos which proliferated and germinated best in 0.75 mg/L BAP-fortified MS medium. Scanning electron microscopic (SEM) images confirmed the presence of various developmental stages of somatic embryos. Shoot regeneration was obtained in which BAP at 1.0 mg/L and 2.0 mg/L BAP + 0.5 mg/L 2,4-D proved to be the best treatments of PGRs in inducing direct and indirect shoot buds. The regenerated shoots showed the highest rooting percentage (87.5%) with 24.7 ± 1.9 numbers of roots/shoot in 1.0 mg/L IBA augmented medium. The rooted plantlets were acclimatized in a greenhouse at a survival rate of 85-90%. The genome size and the 2C nuclear DNA content of field-grown, somatic embryo-regenerated and organogenic-derived plants were estimated and noted to be 3.1, 3.2 and 3.0 picogram (pg), respectively; there is no alteration in ploidy status and the DNA content, validating genetic uniformity. Six SCoT primers unveiled 94.3%-95.13% monomorphic bands across all the plant samples analyzed, further indicating genetic stability among in vitro clones and mother plants. This study describes for the first time successful induction of somatic embryos from hypocotyl callus; and flow cytometry and SCoT marker confirmed the genetic homogeneity of regenerated plants.

MicroRNAs in Medicinal Plants

Medicinal plant microRNAs (miRNAs) are an endogenous class of small RNA central to the posttranscriptional regulation of gene expression. Biosynthetic research has shown that the mature miRNAs in medicinal plants can be produced from either the standard messenger RNA splicing mechanism or the pre-ribosomal RNA splicing process. The medicinal plant miRNA function is separated into two levels: (1) the cross-kingdom level, which is the regulation of disease-related genes in animal cells by oral intake, and (2) the intra-kingdom level, which is the participation of metabolism, development, and stress adaptation in homologous or heterologous plants. Increasing research continues to enrich the biosynthesis and function of medicinal plant miRNAs. In this review, peer-reviewed papers on medicinal plant miRNAs published on the Web of Science were discussed, covering a total of 78 species. The feasibility of the emerging role of medicinal plant miRNAs in regulating animal gene function was critically evaluated. Staged progress in intra-kingdom miRNA research has only been found in a few medicinal plants, which may be mainly inhibited by their long growth cycle, high demand for growth environment, immature genetic transformation, and difficult RNA extraction. The present review clarifies the research significance, opportunities, and challenges of medicinal plant miRNAs in drug development and agricultural production. The discussion of the latest results furthers the understanding of medicinal plant miRNAs and helps the rational design of the corresponding miRNA/target genes functional modules.

Isolation of the 3β-HSD promoter from Digitalis ferruginea subsp. ferruginea and its functional characterization in Arabidopsis thaliana

BACKGROUND

Although members of the SDR gene family (short chain dehydrogenase) are distributed in kingdom of life, they have diverse roles in stress tolerance mechanism or secondary metabolite biosynthesis. Nevertheless, their precise roles in gene expression or regulation under stress are yet to be understood.

METHODS

As a case study, we isolated, sequenced and functionally characterized the 3β-HSD promoter from Digitalis ferruginea subsp. ferruginea in Arabidopsis thaliana.

RESULTS

The promoter fragment contained light and stress response elements such as Box-4, G-Box, TCT-motif, LAMP element, ABRE, ARE, WUN-motif, MYB, MYC, W box, STRE and Box S. The functional analysis of the 3β-HSD promoter in transgenic Arabidopsis seedlings showed that the promoter was expressed in cotyledon and root elongation zone in 2 days' seedlings. However, this expression was extended to hypocotyl and complete root in 6 days' seedlings. In 20 days-old seedlings, promoter expression was distributed to the whole seedling including hydathodes aperture, vascular bundle, shoot apical meristem, trichomes, midrib, leaf primordia, hypocotyl and xylem tissues. Further, expression of the promoter was enhanced or remained stable under the different abiotic stress conditions like osmotic, heat, cold, cadmium or low pH. In addition, the promoter also showed response to methyl jasmonate (MeJA) application. The expression could not be induced in wounded cotyledon most likely due to lack of interacting elements in the promoter fragment.

CONCLUSIONS

Taken together, the 3β-HSD promoter could be a candidate for the development of transgenic plants especially under changing environmental conditions.

Identification of Long Non-coding RNA Transcripts in Glycyrrhiza uralensis

Background:

Chinese liquorice (Glycyrrhiza uralensis), an important medicinal plant, contains various valuable secondary metabolites. Secondary metabolites biosynthesis is very tightly regulated; therefore, elucidation and manipulation of the biosynthetic pathways are of great interest. Recent studies have shown that lncRNAs play important regulatory roles in many biological processes, thus identification and modification of their expression is essential to metabolic pathways for biosynthesis of secondary metabolites.

Objectives:

In this study we attempted to identify non-coding RNA transcripts (lncRNAs) that may act as important regulators of diverse biological processes, including stress responses and developmental programs in Glycyrrhiza uralensis.

Materials and Methods:

Identification of potential lncRNAs in Chinese liquorice was performed using a bioinformatics pipeline from the available EST dataset of G. uralensis.

Results:

Bioinformatics analysis revealed that 1365 identical sequences in the range of 200 to 1286 base pair are putative lncRNAs. Only less than one percent of the predicted lncRNAs display sequence conservation with lncRNAs from other species. Moreover, 13 lncRNAs were detected as the potential precursors of 16 miRNAs. From this analysis, we also detected possible target genes of 16 known miRNA genes. The majority of the predicted miRNA target genes have important role in response to plant disease and a couple of them contribute to signalling and metabolic pathways.

Conclusion:

This study demonstrates the existence of lncRNAs in G. uralensis which has not been found before and provides valuable resources for further understanding and characterizing of lncRNAs and also a basis for additional investigation to reveal specific roles of lncRNAs in various biological processes and particularly in response to plant diseases.

Noncoding RNAs in Medicinal Plants and their Regulatory Roles in Bioactive Compound Production

BACKGROUND

Noncoding RNAs (ncRNAs), such as microRNAs (miRNAs), small interfering RNAs (siRNAs) and long noncoding RNAs (lncRNAs), play significant regulatory roles in plant development and secondary metabolism and are involved in plant response to biotic and abiotic stresses. They have been intensively studied in model systems and crops for approximately two decades and massive amount of information have been obtained. However, for medicinal plants, ncRNAs, particularly their regulatory roles in bioactive compound biosynthesis, are just emerging as a hot research field.

OBJECTIVE

This review aims to summarize current knowledge on herbal ncRNAs and their regulatory roles in bioactive compound production.

RESULTS

So far, scientists have identified thousands of miRNA candidates from over 50 medicinal plant species and 11794 lncRNAs from Salvia miltiorrhiza, Panax ginseng, and Digitalis purpurea. Among them, more than 30 miRNAs and five lncRNAs have been predicted to regulate bioactive compound production.

CONCLUSION

The regulation may achieve through various regulatory modules and pathways, such as the miR397-LAC module, the miR12112-PPO module, the miR156-SPL module, the miR828-MYB module, the miR858-MYB module, and other siRNA and lncRNA regulatory pathways. Further functional analysis of herbal ncRNAs will provide useful information for quality and quantity improvement of medicinal plants.

Quo vadis Cardiac Glycoside Research?

Cardiac glycosides (CGs), toxins well-known for numerous human and cattle poisoning, are natural compounds, the biosynthesis of which occurs in various plants and animals as a self-protective mechanism to prevent grazing and predation. Interestingly, some insect species can take advantage of the CG’s toxicity and by absorbing them, they are also protected from predation. The mechanism of action of CG’s toxicity is inhibition of Na⁺/K⁺-ATPase (the sodium-potassium pump, NKA), which disrupts the ionic homeostasis leading to elevated Ca²⁺ concentration resulting in cell death. Thus, NKA serves as a molecular target for CGs (although it is not the only one) and even though CGs are toxic for humans and some animals, they can also be used as remedies for various diseases, such as cardiovascular ones, and possibly cancer. Although the anticancer mechanism of CGs has not been fully elucidated, yet, it is thought to be connected with the second role of NKA being a receptor that can induce several cell signaling cascades and even serve as a growth factor and, thus, inhibit cancer cell proliferation at low nontoxic concentrations. These growth inhibitory effects are often observed only in cancer cells, thereby, offering a possibility for CGs to be repositioned for cancer treatment serving not only as chemotherapeutic agents but also as immunogenic cell death triggers. Therefore, here, we report on CG’s chemical structures, production optimization, and biological activity with possible use in cancer therapy, as well as, discuss their antiviral potential which was discovered quite recently. Special attention has been devoted to digitoxin, digoxin, and ouabain.

Sequencing, de novo assembly and annotation of Digitalis ferruginea subsp. schischkinii transcriptome

There is an increasing demand for elucidating the biosynthetic pathway of medicinal plants, which are capable of producing several metabolites with great potentials for industrial drug production. Digitalis species are important medicinal plants for the production of cardenolide compounds. Advancement on culture techniques is strictly related to our understanding of the genomic background of species. There are a limited number of genomic studies on Digitalis species. The goal of this study is to contribute to the genomic data of Digitalis ferruginea subsp. schischkinii by presenting transcriptome annotation. Digitalis ferruginea subsp. schischkinii has a limited distribution in Turkey and Transcaucasia, and has a high level of lanatoside C, an important cardenolide. In the study, we sequenced the cDNA library prepared from RNA pools of D. ferruginea subsp. schischkinii tissues treated with various stress conditions. Comprehensive bioinformatics approaches were used for de novo assembly and functional annotation of D. ferruginea subsp. schischkinii transcriptome sequence data along with TF families predictions and phylogenetic analysis. In the study, 58,369 unigenes were predicted and unigenes were annotated by analyzing the sequence data in the non-redundant (NR) protein database, the non-redundant nucleotide (NT) database, Gene Orthology (GO), EuKaryotic Orthologous Groups (KOG), Kyoto Encyclopedia of Genes and Genomes (KEGG), SwissProt, and InterPro databases. This study is the first transcriptome data for D. ferruginea subsp. schischkinii.

Knock-down of δ-aminolevulinic acid dehydratase via virus-induced gene silencing alters the microRNA biogenesis and causes stress-related reactions in citrus plants

The δ-aminolevulinic acid (δ-ALA) is an intermediate in the biosynthetic pathway of tetrapyrroles. Tetrapyrroles play vital roles in many biological processes such as photosynthesis, respiration, and light-sensing. ALA-dehydratase (ALAD) combines two molecules of δ-ALA to form porphobilinogen. In citrus, the silencing of ALAD caused discrete yellow spots and necrosis in leaves and stems. Additionally, it caused rapid death in developing new shoots. Herein, we hypothesize that the accumulation of δ-ALA results in severe stress and reduced meristem development. For that reason, we investigated the dynamic changes in the expression profiles of 23 microRNA (miRNA) identified through small RNA sequencing, from CTV-tALAD plants in comparison with healthy C. macrophylla and C. macrophylla infiltrated with CTV-wt. Furthermore, we reported the effect of ALAD silencing on the total phenolics, H₂O_2, and reactive oxygen species (ROS) levels, to examine the possibilities of miRNAs involving the regulation of these pathways. Our results showed that the total phenolics content, H₂O₂, and O₂^- levels were increased in CTV-tALAD plants. Moreover, 63 conserved miRNA members belonging to 23 different miRNA families were differentially expressed in CTV-tALAD plants compared to controls. The identified miRNAs are implicated in auxin biosynthesis and signaling, axillary shoot meristem formation and leaf morphology, starch metabolism, and oxidative stress. Collectively, our findings suggested that ALAD silencing initiates stress on citrus plants. As a result, CTV-tALAD plants exhibit reduced metabolic rate, growth, and development in order to cope with the stress that resulted from the accumulation of δ-ALA. This cascade of events led to leaf, stem, and meristem necrosis and failure of new shoot development.

Network Approach to Source Attribution of Salmonella enterica Serovar Typhimurium and Its Monophasic Variant

Salmonella enterica subspecies enterica serovar Typhimurium and its monophasic variant are among the most common Salmonella serovars associated with human salmonellosis each year. Related infections are often due to consumption of contaminated meat of pig, cattle, and poultry origin. In order to evaluate novel microbial subtyping methods for source attribution, an approach based on weighted networks was applied on 141 human and 210 food and animal isolates of pigs, broilers, layers, ducks, and cattle collected in Denmark from 2013 to 2014. A whole-genome SNP calling was performed along with cgMLST and wgMLST. Based on these genomic input data, pairwise distance matrices were built and used as input for construction of a weighted network where nodes represent genomes and links to distances. Analyzing food and animal Typhimurium genomes, the coherence of source clustering ranged from 89 to 90% for animal source, from 84 to 85% for country, and from 63 to 65% for year of isolation and was equal to 82% for serotype, suggesting animal source as the first driver of clustering formation. Adding human isolate genomes to the network, a percentage between 93.6 and 97.2% clustered with the existing component and only a percentage between 2.8 and 6.4% appeared as not attributed to any animal sources. The majority of human genomes were attributed to pigs with probabilities ranging from 83.9 to 84.5%, followed by broilers, ducks, cattle, and layers in descending order. In conclusion, a weighted network approach based on pairwise SNPs, cgMLST, and wgMLST matrices showed promising results for source attribution studies.

A novel source of the cardiac glycoside digoxin from the endophytic fungus Epicoccum nigrum: isolation, characterization, production enhancement by gamma irradiation mutagenesis and anticancer activity evaluation

AIMS

Different endophytic fungi were isolated and screened for their digoxin-producing ability. Strain improvement and different culture conditions were studied for more effective production of digoxin.

METHODS AND RESULTS

Among the isolated fungi, an isolate produced digoxin in a concentration of 2·07 mg l^-1 . The digoxin-producing fungal isolate was identified as Epicoccum nigrum Link according to the morphological features and phylogenetic analyses. The potentiality of the fungal strain for production enhancement of digoxin was performed by gamma radiation mutagenesis. Gamma irradiation dose of 1000 Gy intensified the digoxin yield by five-fold. Using this dose, a stable mutant strain with improved digoxin productivity was isolated and the stability for digoxin production was followed up across four successive generations. In the effort to increase digoxin magnitude, selection of the proper cultivation medium, addition of some elicitors to the most proper medium and several physical fermentation conditions were tested. Fermentation process carried out in malt extract autolysate medium (pH 6·5) supplemented by methyl jasmonate and inoculated with 2 ml of 6-day-old culture and incubated at 25°C for 10 days stimulated the highest production of digoxin to attain 50·14 mg l^-1 . Moreover, cytotoxicity of digoxin separated from the fungal culture was tested against five different cancer cell lines. Based on the MTT assay, digoxin inhibited the proliferation of the five different cancer cell lines and the recorded 50% inhibitory concentration ranged from 10·76 to 35·14 μg ml^-1 .

CONCLUSIONS

This is the first report on the production and enhancement of digoxin using fungal fermentation as a new and alternate source with high productivity.

SIGNIFICANCE AND IMPACT OF THE STUDY

These findings offer new and alternate sources with excellent biotechnological potential for digoxin production by fungal fermentation. Moreover, digoxin proved to be a promising anticancer agent whose anticancer potential should be assessed in prospective cancer therapy.

De novo transcriptome analysis deciphered polyoxypregnane glycoside biosynthesis pathway in Gymnema sylvestre

Gymnema sylvestre is an important medicinal plant containing antidiabetic activity. Through de novo transcriptomic study, the pathways of polyoxypregnane glycosides were explored and candidate genes of these pathways were identified in G. sylvestre. High-quality raw reads were assembled into transcripts which resulted in 193,615 unigenes. These unigenes further decoded 58,274 coding DNA sequences (CDSs). Functional annotation of predicted CDSs was carried out using the protein databases, i.e., NCBI's non-redundant, Uniprot and Pfam. Eukaryotic orthologous group (KOG) classification and transcription factor analysis has revealed most CDS-enriched categories as "Signal transduction mechanism" and "Basic Helix loop helix" (bHLH) transcription factor family, respectively. A total of 16,569 CDSs were assigned minimum one Gene Ontology (GO) term. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis disclosed 235 CDSs which represented total 27 genes of pregnane glycoside pathways and 19 CDSs represented 10 important enzymes of polyoxypregnane glycoside biosynthesis, i.e., sterol 24-C-methyltransferase, cycloeucalenol cycloisomerase, Δ¹⁴-sterol reductase, C-8,7 sterol isomerase, sterol methyltransferase 2, C-5 sterol desaturase, sterol Δ⁷ reductase, Δ²⁴ sterol reductase, 3β-hydroxysteroid dehydrogenase and progesterone 5β reductase (5βPOR). This transcriptome analysis provided an important resource for future functional genomic studies in G. sylvestre.

Integrated analysis of transcriptomic and metabolomic data reveals critical metabolic pathways involved in rotenoid biosynthesis in the medicinal plant Mirabilis himalaica

Mirabilis himalaica (Edgew.) Heimerl is among the most important genuine medicinal plants in Tibet. However, the biosynthesis mechanisms of the active compounds in this species are unclear, severely limiting its application. To clarify the molecular biosynthesis mechanism of the key representative active compounds, specifically rotenoid, which is of special medicinal value for M. himalaica, RNA sequencing and TOF-MS technologies were used to construct transcriptomic and metabolomic libraries from the roots, stems, and leaves of M. himalaica plants collected from their natural habitat. As a result, each of the transcriptomic libraries from the different tissues was sequenced, generating more than 10 Gb of clean data ultimately assembled into 147,142 unigenes. In the three tissues, metabolomic analysis identified 522 candidate compounds, of which 170 metabolites involved in 114 metabolic pathways were mapped to the KEGG. Of these genes, 61 encoding enzymes were identified to function at key steps of the pathways related to rotenoid biosynthesis, where 14 intermediate metabolites were also located. An integrated analysis of metabolic and transcriptomic data revealed that most of the intermediate metabolites and enzymes related to rotenoid biosynthesis were synthesized in the roots, stems and leaves of M. himalaica, which suggested that the use of non-medicinal tissues to extract compounds was feasible. In addition, the CHS and CHI genes were found to play important roles in rotenoid biosynthesis, especially, since CHS might be an important rate-limiting enzyme. This study provides a hypothetical basis for the screening of new active metabolites and the metabolic engineering of rotenoid in M. himalaica.

Researches on Transcriptome Sequencing in the Study of Traditional Chinese Medicine

Due to its incomparable advantages, the application of transcriptome sequencing in the study of traditional Chinese medicine attracts more and more attention of researchers, which greatly promote the development of traditional Chinese medicine. In this paper, the applications of transcriptome sequencing in traditional Chinese medicine were summarized by reviewing recent related papers.

Transcriptome and Metabolite analysis reveal candidate genes of the cardiac glycoside biosynthetic pathway from Calotropis procera

Calotropis procera is a medicinal plant of immense importance due to its pharmaceutical active components, especially cardiac glycosides (CG). As genomic resources for this plant are limited, the genes involved in CG biosynthetic pathway remain largely unknown till date. Our study on stage and tissue specific metabolite accumulation showed that CG's were maximally accumulated in stems of 3 month old seedlings. De novo transcriptome sequencing of same was done using high throughput Illumina HiSeq platform generating 44074 unigenes with average mean length of 1785 base pair. Around 66.6% of unigenes were annotated by using various public databases and 5324 unigenes showed significant match in the KEGG database involved in 133 different pathways of plant metabolism. Further KEGG analysis resulted in identification of 336 unigenes involved in cardenolide biosynthesis. Tissue specific expression analysis of 30 putative transcripts involved in terpenoid, steroid and cardenolide pathways showed a positive correlation between metabolite and transcript accumulation. Wound stress elevated CG levels as well the levels of the putative transcripts involved in its biosynthetic pathways. This result further validated the involvement of identified transcripts in CGs biosynthesis. The identified transcripts will lay a substantial foundation for further research on metabolic engineering and regulation of cardiac glycosides biosynthesis pathway genes.

Validation of Suitable Reference Genes for Assessing Gene Expression of MicroRNAs in Lonicera japonica

MicroRNAs (miRNAs), which play crucial regulatory roles in plant secondary metabolism and responses to the environment, could be developed as promising biomarkers for different varieties and production areas of herbal medicines. However, limited information is available for miRNAs from Lonicera japonica, which is widely used in East Asian countries owing to various pharmaceutically active secondary metabolites. Selection of suitable reference genes for quantification of target miRNA expression through quantitative real-time (qRT)-PCR is important for elucidating the molecular mechanisms of secondary metabolic regulation in different tissues and varieties of L. japonica. For precise normalization of gene expression data in L. japonica, 16 candidate miRNAs were examined in three tissues, as well as 21 cultivated varieties collected from 16 production areas, using GeNorm, NormFinder, and RefFinder algorithms. Our results revealed combination of u534122 and u3868172 as the best reference genes across all samples. Their specificity was confirmed by detecting the cycling threshold (C t) value ranges in different varieties of L. japonica collected from diverse production areas, suggesting the use of these two reference miRNAs is sufficient for accurate transcript normalization with different tissues, varieties, and production areas. To our knowledge, this is the first report on validation of reference miRNAs in honeysuckle (Lonicera spp.). Restuls from this study can further facilitate discovery of functional regulatory miRNAs in different varieties of L. japonica.

Renaissance in phytomedicines: promising implications of NGS technologies

Medicinal plant research is growing significantly in faith to discover new and more biologically compatible phytomedicines. Deposition of huge genome/trancriptome sequence data assisted by NGS technologies has revealed the new possibilities for producing upgraded bioactive molecules in medicinal plants. Growing interest of investors and consumers in the herbal drugs raises the need for extensive research to open the facts and details of every inch of life canvas of medicinal plants to produce improved quality of phytomedicines. As in agriculture crops, knowledge emergence from medicinal plant's genome/transcriptome, can be used to assure their amended quality and these improved varieties are then transported to the fields for cultivation. Genome studies generate huge sequence data which can be exploited further for obtaining information regarding genes/gene clusters involved in biosynthesis as well as regulation. This can be achieved rapidly at a very large scale with NGS platforms. Identification of new RNA molecules has become possible, which can lead to the discovery of novel compounds. Sequence information can be combined with advanced phytochemical and bioinformatics tools to discover functional herbal drugs. Qualitative and quantitative analysis of small RNA species put a light on the regulatory aspect of biosynthetic pathways for phytomedicines. Inter or intra genomic as well as transcriptomic interactive processes for biosynthetic pathways can be elucidated in depth. Quality management of herbal material will also become rapid and high throughput. Enrichment of sequence information will be used to engineer the plants to get more efficient phytopharmaceuticals. The present review comprises of role of NGS technologies to boost genomic studies of pharmaceutically important plants and further, applications of sequence information aiming to produce enriched phytomedicines. Emerging knowledge from the medicinal plants genome/transcriptome can give birth to deep understanding of the processes responsible for biosynthesis of medicinally important compounds.

Genome-wide identification and characterization of novel lncRNAs in Populus under nitrogen deficiency

Long non-coding RNAs (lncRNAs) have been identified as important regulatory factors of gene expression in eukaryotic species, such as Homo sapiens, Arabidopsis thaliana, and Oryza sativa. However, the systematic identification of potential lncRNAs in trees is comparatively rare. In particular, the characteristics, expression, and regulatory roles of lncRNAs in trees under nutrient stress remain largely unknown. A genome-wide strategy was used in this investigation to identify and characterize novel and low-nitrogen (N)-responsive lncRNAs in Populus tomentosa; 388 unique lncRNA candidates belonging to 380 gene loci were detected and only seven lncRNAs were found to belong to seven conserved non-coding RNA families indicating the majority of P. tomentosa lncRNAs are species-specific. In total, 126 lncRNAs were significantly altered under low-N stress; 8 were repressed, and 118 were induced. Furthermore, 9 and 5 lncRNAs were detected as precursors of 11 known and 14 novel Populus miRNAs, respectively, whereas 4 lncRNAs were targeted by 29 miRNAs belonging to 5 families, including 22 conserved and 7 non-conserved miRNAs. In addition, 15 antisense lncRNAs were identified to be generated from opposite strands of 14 corresponding protein-coding genes. In total, 111 protein-coding genes with regions complementary to 38 lncRNAs were also predicted with some lncRNAs corresponding to multiple genes and vice versa, and their functions were annotated, which further demonstrated the complex regulatory relationship between lncRNAs and protein-coding genes in plants. Moreover, an interaction network among lncRNAs, miRNAs, and mRNAs was investigated. These findings enrich our understanding of lncRNAs in Populus, expand the methods of miRNA identification. Our results present the first global characterization of lncRNAs and their potential target genes in response to nitrogen stress in trees, which provides more information on low-nutrition adaptation mechanisms in woody plants.

Plastoquinone and Ubiquinone in Plants: Biosynthesis, Physiological Function and Metabolic Engineering

Plastoquinone (PQ) and ubiquinone (UQ) are two important prenylquinones, functioning as electron transporters in the electron transport chain of oxygenic photosynthesis and the aerobic respiratory chain, respectively, and play indispensable roles in plant growth and development through participating in the biosynthesis and metabolism of important chemical compounds, acting as antioxidants, being involved in plant response to stress, and regulating gene expression and cell signal transduction. UQ, particularly UQ₁₀, has also been widely used in people's life. It is effective in treating cardiovascular diseases, chronic gingivitis and periodontitis, and shows favorable impact on cancer treatment and human reproductive health. PQ and UQ are made up of an active benzoquinone ring attached to a polyisoprenoid side chain. Biosynthesis of PQ and UQ is very complicated with more than thirty five enzymes involved. Their synthetic pathways can be generally divided into two stages. The first stage leads to the biosynthesis of precursors of benzene quinone ring and prenyl side chain. The benzene quinone ring for UQ is synthesized from tyrosine or phenylalanine, whereas the ring for PQ is derived from tyrosine. The prenyl side chains of PQ and UQ are derived from glyceraldehyde 3-phosphate and pyruvate through the 2-C-methyl-D-erythritol 4-phosphate pathway and/or acetyl-CoA and acetoacetyl-CoA through the mevalonate pathway. The second stage includes the condensation of ring and side chain and subsequent modification. Homogentisate solanesyltransferase, 4-hydroxybenzoate polyprenyl diphosphate transferase and a series of benzene quinone ring modification enzymes are involved in this stage. PQ exists in plants, while UQ widely presents in plants, animals and microbes. Many enzymes and their encoding genes involved in PQ and UQ biosynthesis have been intensively studied recently. Metabolic engineering of UQ₁₀ in plants, such as rice and tobacco, has also been tested. In this review, we summarize and discuss recent research progresses in the biosynthetic pathways of PQ and UQ and enzymes and their encoding genes involved in side chain elongation and in the second stage of PQ and UQ biosynthesis. Physiological functions of PQ and UQ played in plants as well as the practical application and metabolic engineering of PQ and UQ are also included.

Novel and conserved miRNAs in the halophyte Suaeda maritima identified by deep sequencing and computational predictions using the ESTs of two mangrove plants

BACKGROUND

Although miRNAs are reportedly involved in the salt stress tolerance of plants, miRNA profiling in plants has largely remained restricted to glycophytes, including certain crop species that do not exhibit any tolerance to salinity. Hence, this manuscript describes the results from the miRNA profiling of the halophyte Suaeda maritima, which is used worldwide to study salt tolerance in plants.

RESULTS

A total of 134 conserved miRNAs were identified from unique sRNA reads, with 126 identified using miRBase 21.0 and an additional eight identified using the Plant Non-coding RNA Database. The presence of the precursors of seven conserved miRNAs was validated in S. maritima. In addition, 13 novel miRNAs were predicted using the ESTs of two mangrove plants, Rhizophora mangle and Heritiera littoralis, and the precursors of seven miRNAs were found in S. maritima. Most of the miRNAs considered for characterization were responsive to NaCl application, indicating their importance in the regulation of metabolic activities in plants exposed to salinity. An expression study of the novel miRNAs in plants of diverse ecological and taxonomic groups revealed that two of the miRNAs, sma-miR6 and sma-miR7, were also expressed in Oryza sativa, whereas another two, sma-miR2 and sma-miR5, were only expressed in plants growing under the influence of seawater, similar to S. maritima.

CONCLUSION

The distribution of conserved miRNAs among only 25 families indicated the possibility of identifying a greater number of miRNAs with increase in knowledge of the genomes of more halophytes. The expression of two novel miRNAs, sma-miR2 and sma-miR5, only in plants growing under the influence of seawater suggested their metabolic regulatory roles specific to saline environments, and such behavior might be mediated by alterations in the expression of certain genes, modifications of proteins leading to changes in their activity and production of secondary metabolites as revealed by the miRNA target predictions. Moreover, the auxin responsive factor targeted by sma-miR7 could also be involved in salt tolerance because the target is conserved between species. This study also indicated that the transcriptome of one species can be successfully used to computationally predict the miRNAs in other species, especially those that have similar metabolism, even if they are taxonomically separated.

New insights of medicinal plant therapeutic activity-The miRNA transfer

MicroRNA (miRNA) has become the spotlight of the biomedical research around the world and is considered to be a major post-transcriptional gene regulator. This small, endogenous RNA (21-25 nucleotides long) plays an important role by targeting specific mRNAs in plants, animals and humans. Herbal medicine has been used for thousands of years, however little is known about its molecular mechanism of action. Since the discovery of plant miRNA in human tissue and sera after ingestion, the connection between the two kingdoms is presented under a new perspective. Forward pharmacology, such as miRNA therapeutics could be the next best step toward identifying novel therapeutic options involving medicinal plants. Besides reporting the latest findings regarding the cross-kingdom transfer of miRNA and its therapeutic application, this review can inform further investigations that could lead to a modern definition of herbal medicine.

Maternal and zygotic transcriptomes in the appendicularian, Oikopleura dioica: novel protein-encoding genes, intra-species sequence variations, and trans-spliced RNA leader

RNA sequencing analysis was carried out to characterize egg and larval transcriptomes in the appendicularian, Oikopleura dioica, a planktonic chordate, which is characterized by rapid development and short life cycle of 5 days, using a Japanese population of the organism. De novo transcriptome assembly matched with 16,423 proteins corresponding to 95.4% of the protein-encoding genes deposited in the OikoBase, the genome database of the Norwegian population. Nucleotide and amino acid sequence identities between the Japanese and Norwegian O. dioica were estimated to be around 91.0 and 94.8%, respectively. We discovered 175 novel protein-encoding genes: 144 unigenes were common to both the Japanese and Norwegian populations, whereas 31 unigenes were not found in the OikoBase genome reference. Among the total 12,311 unigenes, approximately 63% were detected in egg-stage RNAs, whereas 99% were detected in larval stage RNAs; 3772 genes were up-regulated, and 1336 genes were down-regulated more than four-fold in the larvae. Gene ontology analyses characterized gene activities in these two developmental stages. We found a messenger RNA (mRNA) 5' trans-spliced leader, which was observed in 40.8% of the total unique transcripts. It showed preferential linkage to adenine at the 5' ends of the downstream exons. Trans-splicing was observed more frequently in egg mRNAs compared with larva-specific mRNAs.

Identification and characterization of mRNA-like noncoding RNAs in Salvia miltiorrhiza

Identification and characterization of 5,446 mlncRNAs from Salvia miltiorrhiza showed that the majority of identified mlncRNAs were stress responsive, providing a framework for elucidating mlncRNA functions in S. miltiorrhiza. mRNA-like noncoding RNAs (mlncRNAs) are transcribed by RNA polymerase II and are polyadenylated, capped and spliced. They play important roles in plant development and defense responses. However, there is no information available for mlncRNAs in Salvia miltiorrhiza Bunge, the first Chinese medicinal material entering the international market. To perform a transcriptome-wide identification of S. miltiorrhiza mlncRNAs, we assembled over 8 million RNA-seq reads from GenBank database and 5,624 ESTs from PlantGDB into 44422 unigenes. Using a computational identification pipeline, we identified 5446 S. miltiorrhiza mlncRNA candidates from the assembled unigenes. Of the 5446 mlncRNAs, 2 are primary transcripts of conserved miRNAs, and 2030 can be grouped into 470 families with at least two members in a family. Quantitative real-time PCR analysis of mlncRNAs with at least 900 nt showed that the majority were differentially expressed in roots, stems, leaves and flowers and responsive to methyl jasmonate (MeJA) treatment in S. miltiorrhiza. Analysis of published RNA-seq data showed that a total of 3,044 mlncRNAs were expressed in hairy roots of S. miltiorrhiza and the expression of 1,904 of the 3,044 mlncRNAs was altered by yeast extract and Ag(+) treatment. The results indicate that the majority of mlncRNAs are involved in plant response to stress. This study provides a framework for understanding the roles of mlncRNAs in S. miltiorrhiza.

Identification of mRNA-like non-coding RNAs and validation of a mighty one named MAR in Panax ginseng

Increasing evidence suggests that long non-coding RNAs (lncRNAs) play significant roles in plants. However, little is known about lncRNAs in Panax ginseng C. A. Meyer, an economically significant medicinal plant species. A total of 3,688 mRNA-like non-coding RNAs (mlncRNAs), a class of lncRNAs, were identified in P. ginseng. Approximately 40% of the identified mlncRNAs were processed into small RNAs, implying their regulatory roles via small RNA-mediated mechanisms. Eleven miRNA-generating mlncRNAs also produced siRNAs, suggesting the coordinated production of miRNAs and siRNAs in P. ginseng. The mlncRNA-derived small RNAs might be 21-, 22-, or 24-nt phased and could be generated from both or only one strand of mlncRNAs, or from super long hairpin structures. A full-length mlncRNA, termed MAR (multiple-function-associated mlncRNA), was cloned. It generated the most abundant siRNAs. The MAR siRNAs were predominantly 24-nt and some of them were distributed in a phased pattern. A total of 228 targets were predicted for 71 MAR siRNAs. Degradome sequencing validated 68 predicted targets involved in diverse metabolic pathways, suggesting the significance of MAR in P. ginseng. Consistently, MAR was detected in all tissues analyzed and responded to methyl jasmonate (MeJA) treatment. It sheds light on the function of mlncRNAs in plants.

Comparative transcriptomic analysis of the response to cold acclimation in Eucalyptus dunnii

Eucalyptus dunnii is an important macrophanerophyte with high economic value. However, low temperature stress limits its productivity and distribution. To study the cold response mechanisms of E. dunnii, 5 cDNA libraries were constructed from mRNA extracted from leaves exposed to cold stress for varying lengths of time and were evaluated by RNA-Seq analysis. The assembly of the Illumina datasets was optimized using various assembly programs and parameters. The final optimized assembly generated 205,325 transcripts with an average length of 1,701 bp and N50 of 2,627 bp, representing 349.38 Mb of the E. dunnii transcriptome. Among these transcripts, 134,358 transcripts (65.4%) were annotated in the Nr database. According to the differential analysis results, most transcripts were up-regulated as the cold stress prolonging, suggesting that these transcripts may be involved in the response to cold stress. In addition, the cold-relevant GO categories, such as ‘response to stress’ and ‘translational initiation’, were the markedly enriched GO terms. The assembly of the E. dunnii gene index and the GO classification performed in this study will serve as useful genomic resources for the genetic improvement of E. dunnii and also provide insights into the molecular mechanisms of cold acclimation in E. dunnii.

Comprehensive analysis of alternative splicing in Digitalis purpurea by strand-specific RNA-Seq

Digitalis purpurea (D. purpurea) is one of the most important medicinal plants and is well known in the treatment of heart failure because of the cardiac glycosides that are its main active compounds. However, in the absence of strand specific sequencing information, the post-transcriptional mechanism of gene regulation in D. purpurea thus far remains unknown. In this study, a strand-specific RNA-Seq library was constructed and sequenced using Illumina HiSeq platforms to characterize the transcriptome of D. purpurea with a focus on alternative splicing (AS) events and the effect of AS on protein domains. De novo RNA-Seq assembly resulted in 48,475 genes. Based on the assembled transcripts, we reported a list of 3,265 AS genes, including 5,408 AS events in D. purpurea. Interestingly, both glycosyltransferases and monooxygenase, which were involved in the biosynthesis of cardiac glycosides, are regulated by AS. A total of 2,422 AS events occurred in coding regions, and 959 AS events were located in the regions of 882 unique protein domains, which could affect protein function. This D. purpurea transcriptome study substantially increased the expressed sequence resource and presented a better understanding of post-transcriptional regulation to further facilitate the medicinal applications of D. purpurea for human health.

Genome-wide identification and characterization of long intergenic non-coding RNAs in Ganoderma lucidum

Ganoderma lucidum is a white-rot fungus best-known for its medicinal activities. We have previously sequenced its genome and annotated the protein coding genes. However, long non-coding RNAs in G. lucidum genome have not been analyzed. In this study, we have identified and characterized long intergenic non-coding RNAs (lincRNA) in G. lucidum systematically. We developed a computational pipeline, which was used to analyze RNA-Seq data derived from G. lucidum samples collected from three developmental stages. A total of 402 lincRNA candidates were identified, with an average length of 609 bp. Analysis of their adjacent protein-coding genes (apcGenes) revealed that 46 apcGenes belong to the pathways of triterpenoid biosynthesis and lignin degradation, or families of cytochrome P450, mating type B genes, and carbohydrate-active enzymes. To determine if lincRNAs and these apcGenes have any interactions, the corresponding pairs of lincRNAs and apcGenes were analyzed in detail. We developed a modified 3' RACE method to analyze the transcriptional direction of a transcript. Among the 46 lincRNAs, 37 were found unidirectionally transcribed, and 9 were found bidirectionally transcribed. The expression profiles of 16 of these 37 lincRNAs were found to be highly correlated with those of the apcGenes across the three developmental stages. Among them, 11 are positively correlated (r>0.8) and 5 are negatively correlated (r<-0.8). The co-localization and co-expression of lincRNAs and those apcGenes playing important functions is consistent with the notion that lincRNAs might be important regulators for cellular processes. In summary, this represents the very first study to identify and characterize lincRNAs in the genomes of basidiomycetes. The results obtained here have laid the foundation for study of potential lincRNA-mediated expression regulation of genes in G. lucidum.

Transcriptome analysis of differentially expressed genes relevant to variegation in peach flowers

BACKGROUND

Variegation in flower color is commonly observed in many plant species and also occurs on ornamental peaches (Prunus persica f. versicolor [Sieb.] Voss). Variegated plants are highly valuable in the floricultural market. To gain a global perspective on genes differentially expressed in variegated peach flowers, we performed large-scale transcriptome sequencing of white and red petals separately collected from a variegated peach tree.

RESULTS

A total of 1,556,597 high-quality reads were obtained, with an average read length of 445 bp. The ESTs were assembled into 16,530 contigs and 42,050 singletons. The resulting unigenes covered about 60% of total predicted genes in the peach genome. These unigenes were further subjected to functional annotation and biochemical pathway analysis. Digital expression analysis identified a total of 514 genes differentially expressed between red and white flower petals. Since peach flower coloration is determined by the expression and regulation of structural genes relevant to flavonoid biosynthesis, a detailed examination detected four key structural genes, including C4H, CHS, CHI and F3H, expressed at a significantly higher level in red than in white petal. Except for the structural genes, we also detected 11 differentially expressed regulatory genes relating to flavonoid biosynthesis. Using the differentially expressed structural genes as the test objects, we validated the digital expression results by using quantitative real-time PCR, and the differential expression of C4H, CHS and F3H were confirmed.

CONCLUSION

In this study, we generated a large EST collection from flower petals of a variegated peach. By digital expression analysis, we identified an informative list of candidate genes associated with variegation in peach flowers, which offered a unique opportunity to uncover the genetic mechanisms underlying flower color variegation.

A stable and efficient Agrobacterium tumefaciens-mediated genetic transformation of the medicinal plant Digitalis purpurea L

In this study, we developed a rapid and efficient method for in vitro propagation and Agrobacterium tumefaciens-mediated transformation of Digitalis purpurea L. (syn. foxglove), an important medicinal plant. Mature leaf explants of D. purpurea were used for 100 % adventitious shoot regeneration on Murashige and Skoog (MS) medium supplemented with 1 mg L(-1) thidiazuron (TDZ) (a cytokine) and 0.1 mg L(-1) 1-naphthaleneacetic acid (NAA) (an auxin). Transformation was achieved by inoculating leaf explants with the A. tumefaciens strains GV2260/pBI121 or GV3101/pBI121. The binary vector pBI121 contained the reporter β-glucuronidase gene (GUS) and kanamycin selection marker nptII. Kanamycin-resistant shoots were regenerated directly on the selection medium 4-6 weeks after co-cultivation. Approximately, 52.2 and 60 % of kanamycin-resistant shoots transformed with Agrobacterium strains GV2260 and GV3101, respectively, showed strong GUS staining by histochemical assay. Furthermore, PCR and Southern blot analysis confirmed the presence of nptII and GUS on the chromosome of the transformed D. purpurea plants, and stable GUS expression was detected in the transformants by RT-PCR analysis. This efficient method of shoot regeneration and genetic transformation of D. purpurea will provide a powerful tool to increase and produce valuable components such as digitoxin, digoxin, and digoxigenin in D. purpurea through improved secondary metabolic pathways via a biotechnological approach.

Predicting gene regulatory networks of soybean nodulation from RNA-Seq transcriptome data

BACKGROUND

High-throughput RNA sequencing (RNA-Seq) is a revolutionary technique to study the transcriptome of a cell under various conditions at a systems level. Despite the wide application of RNA-Seq techniques to generate experimental data in the last few years, few computational methods are available to analyze this huge amount of transcription data. The computational methods for constructing gene regulatory networks from RNA-Seq expression data of hundreds or even thousands of genes are particularly lacking and urgently needed.

RESULTS

We developed an automated bioinformatics method to predict gene regulatory networks from the quantitative expression values of differentially expressed genes based on RNA-Seq transcriptome data of a cell in different stages and conditions, integrating transcriptional, genomic and gene function data. We applied the method to the RNA-Seq transcriptome data generated for soybean root hair cells in three different development stages of nodulation after rhizobium infection. The method predicted a soybean nodulation-related gene regulatory network consisting of 10 regulatory modules common for all three stages, and 24, 49 and 70 modules separately for the first, second and third stage, each containing both a group of co-expressed genes and several transcription factors collaboratively controlling their expression under different conditions. 8 of 10 common regulatory modules were validated by at least two kinds of validations, such as independent DNA binding motif analysis, gene function enrichment test, and previous experimental data in the literature.

CONCLUSIONS

We developed a computational method to reliably reconstruct gene regulatory networks from RNA-Seq transcriptome data. The method can generate valuable hypotheses for interpreting biological data and designing biological experiments such as ChIP-Seq, RNA interference, and yeast two hybrid experiments.

The phenylalanine ammonia-lyase gene family in Salvia miltiorrhiza: genome-wide characterization, molecular cloning and expression analysis

Salvia miltiorrhiza Bunge is a well-known material of traditional Chinese medicine. Hydrophilic phenolic acids, such as rosmarinic acid and salvianolic acid B, are a group of pharmaceutically important compounds in S. miltiorrhiza. The biosynthesis of rosmarinic acid requires the coordination of the phenylpropanoid pathway and the tyrosine-derived pathway. Phenylalanine ammonia-lyase (PAL) is the first key enzyme of the phenylpropanoid pathway. Systematic analysis of the SmPAL gene family has not been carried out. We report here the identification of three SmPALs through searching the recently obtained working draft of the S. miltiorrhiza genome and full-length cDNA cloning. Bioinformatic and phylogenetic analyses showed that SmPAL1 and SmPAL3 clustered in a sub-clade of dicot PALs, whereas SmPAL2 fell into the other one. Some important cis-elements were conserved in three SmPAL promoters, whereas the others were not. SmPAL1 and SmPAL3 were highly expressed in roots and leaves of S. miltiorrhiza, but SmPAL2 were predominately expressed in stems and flowers. It indicates that SmPAL1 and SmPAL3 function redundantly in rosmarinic acid biosynthesis. All SmPALs were induced in roots treated with PEG and MeJA, but the time and degree of responses were different, suggesting the complexity of SmPAL-associated metabolic network in S. miltiorrhiza. This is the first comprehensive study dedicated to SmPAL gene family characterization. The results provide a basis for elucidating the role of SmPAL genes in the biosynthesis of bioactive compounds.

Bioinformatics challenges in de novo transcriptome assembly using short read sequences in the absence of a reference genome sequence

Plant natural product research can be facilitated through genome and transcriptome sequencing approaches that generate informative sequence and expression datasets that enable characterization of biochemical pathways of interest. As the overwhelming majority of plant-derived natural products are derived from species with little, if any, sequence and/or genomic resources, the ability to perform whole genome shotgun sequencing and assembly has been and will continue to be transformative as access to a genome sequence provides molecular resources and a context for discovery and characterization of biosynthetic pathways. Due to the reduced size and complexity of the transcriptome relative to the genome, transcriptome sequencing provides a rapid, inexpensive approach to access gene sequences, gene expression abundances, and gene expression patterns in any species, including those that lack a reference genome sequence. To date, successful applications of RNA sequencing in conjunction with de novo transcriptome assembly has enabled identification of new genes in an array of biochemical pathways in plants. While sequencing technologies are well developed, challenges remain in the handling and analysis of transcriptome sequences. In this Highlight article, we provide an overview of the bioinformatics challenges associated with transcriptome analyses using short read sequences and how to address these issues in plant species that lack a reference genome.

Transcriptome analysis of buds and leaves using 454 pyrosequencing to discover genes associated with the biosynthesis of active ingredients in Lonicera japonica Thunb

BACKGROUND

Lonicera japonica Thunb. is a plant used in traditional Chinese medicine known for its anti-inflammatory, anti-oxidative, anti-carcinogenic, and antiviral pharmacological properties. The major active secondary metabolites of this plant are chlorogenic acid (CGA) and luteoloside. While the biosynthetic pathways of these metabolites are relatively well known, the genetic information available for this species, especially the biosynthetic pathways of its active ingredients, is limited.

METHODOLOGY/PRINCIPAL FINDINGS

We obtained one million reads (average length of 400 bp) in a whole sequence run using a Roche/454 GS FLX titanium platform. Altogether, 85.69% of the unigenes covering the entire life cycle of the plant were annotated and 325 unigenes were assigned to secondary metabolic pathways. Moreover, 2039 unigenes were predicted as transcription factors. Nearly all of the possible enzymes involved in the biosynthesis of CGA and luteoloside were discovered in L. japonica. Three hydroxycinnamoyl transferase genes, including two hydroxycinnamoyl-CoA quinate hydroxycinnamoyl transferase genes and one hydroxycinnamoyl-CoA shikimate/quinate hydroxycinnamoyl transferase (HCT) gene featuring high similarity to known genes from other species, were cloned. The HCT gene was discovered for the first time in L. japonica. In addition, 188 candidate cytochrome P450 unigenes and 245 glycosyltransferase unigenes were found in the expressed sequence tag (EST) dataset.

CONCLUSION

This study provides a high quality EST database for L. japonica by 454 pyrosequencing. Based on the EST annotation, a set of putative genes involved in CGA and luteoloside biosynthetic pathways were discovered. The database serves as an important source of public information on genetic markers, gene expression, genomics, and functional genomics in L. japonica.

Transcriptome analysis of the differentially expressed genes in the male and female shrub willows (Salix suchowensis)

Background

The dioecious system is relatively rare in plants. Shrub willow is an annual flowering dioecious woody plant, and possesses many characteristics that lend it as a great model for tracking the missing pieces of sex determination evolution. To gain a global view of the genes differentially expressed in the male and female shrub willows and to develop a database for further studies, we performed a large-scale transcriptome sequencing of flower buds which were separately collected from two types of sexes.

Results

Totally, 1,201,931 high quality reads were obtained, with an average length of 389 bp and a total length of 467.96 Mb. The ESTs were assembled into 29,048 contigs, and 132,709 singletons. These unigenes were further functionally annotated by comparing their sequences to different proteins and functional domain databases and assigned with Gene Ontology (GO) terms. A biochemical pathway database containing 291 predicted pathways was also created based on the annotations of the unigenes. Digital expression analysis identified 806 differentially expressed genes between the male and female flower buds. And 33 of them located on the incipient sex chromosome of Salicaceae, among which, 12 genes might involve in plant sex determination empirically. These genes were worthy of special notification in future studies.

Conclusions

In this study, a large number of EST sequences were generated from the flower buds of a male and a female shrub willow. We also reported the differentially expressed genes between the two sex-type flowers. This work provides valuable information and sequence resources for uncovering the sex determining genes and for future functional genomics analysis of Salicaceae spp.

In vitro propagation and production of cardiotonic glycosides in shoot cultures of Digitalis purpurea L. by elicitation and precursor feeding

Digitalis purpurea L. (Scrophulariaceae; Foxglove) is a source of cardiotonic glycosides such as digitoxin and digoxin which are commercially applied in the treatment to strengthen cardiac diffusion and to regulate heart rhythm. This investigation deals with in vitro propagation and elicited production of cardiotonic glycosides digitoxin and digoxin in shoot cultures of D. purpurea L. In vitro germinated seedlings were used as a primary source of explants. Multiple shoot formation was achieved for three explant types (nodal, internodal, and leaf) cultured on Murashige and Skoog (MS) medium with several treatments of cytokinins (6-benzyladenine-BA; kinetin-Kin; and thidiazuron-TDZ) and auxins (indole-3-acetic acid-IAA; α-naphthaleneacetic acid-NAA; and 2,4-dichlorophenoxy acetic acid-2,4-D). Maximum multiple shoots (12.7 ± 0.6) were produced from nodal explants on MS + 7.5 μM BA. Shoots were rooted in vitro on MS containing 15 μM IAA. Rooted plantlets were successfully acclimatized. To further maintain the multiple shoot induction, mother tissue was cut into four equal parts and repeatedly sub-cultured on fresh shoot induction liquid medium after each harvest. On adaptation of this strategy, an average of 18 shoots per explant could be produced. This strategy was applied for the production of biomass and glycosides digitoxin and digoxin in shoot cultures on MS medium supplemented with 7.5 μM BA and several treatments with plant growth regulators, incubation period, abiotic (salicylic acid, mannitol, sorbitol, PEG-6000, NaCl, and KCl), biotic (Aspergillus niger, Helminthosporium sp., Alternaria sp., chitin, and yeast extract) elicitors, and precursors (progesterone, cholesterol, and squalene). The treatment of KCl, mycelial mass of Helminthosporium sp., and progesterone were highly effective for the production of cardenolides. In the presence of progesterone (200 to 300 mg/l), digitoxin and digoxin accumulation was enhanced by 9.1- and 11.9-folds respectively.

High-throughput sequencing and characterization of the small RNA transcriptome reveal features of novel and conserved microRNAs in Panax ginseng

microRNAs (miRNAs) play vital regulatory roles in many organisms through direct cleavage of transcripts, translational repression, or chromatin modification. Identification of miRNAs has been carried out in various plant species. However, no information is available for miRNAs from Panax ginseng, an economically significant medicinal plant species. Using the next generation high-throughput sequencing technology, we obtained 13,326,328 small RNA reads from the roots, stems, leaves and flowers of P. ginseng. Analysis of these small RNAs revealed the existence of a large, diverse and highly complicated small RNA population in P. ginseng. We identified 73 conserved miRNAs, which could be grouped into 33 families, and 28 non-conserved ones belonging to 9 families. Characterization of P. ginseng miRNA precursors revealed many features, such as production of two miRNAs from distinct regions of a precursor, clusters of two precursors in a transcript, and generation of miRNAs from both sense and antisense transcripts. It suggests the complexity of miRNA production in P. gingseng. Using a computational approach, we predicted for the conserved and non-conserved miRNA families 99 and 31 target genes, respectively, of which eight were experimentally validated. Among all predicted targets, only about 20% are conserved among various plant species, whereas the others appear to be non-conserved, indicating the diversity of miRNA functions. Consistently, many miRNAs exhibited tissue-specific expression patterns. Moreover, we identified five dehydration- and ten heat-responsive miRNAs and found the existence of a crosstalk among some of the stress-responsive miRNAs. Our results provide the first clue to the elucidation of miRNA functions in P. ginseng.