De novo sequencing transcriptome of endemic Gentiana straminea (Gentianaceae) to identify genes involved in the biosynthesis of active ingredients

Cloning and Expression Analysis of Key Enzyme Gene CoGPPS Involved in Iridoid Glycoside Synthesis in Cornus officinalis

Cornus iridoid glycosides (CIGs), including loganin and morroniside, are the main active components of Cornus officinalis. As one of the key enzymes in the biosynthesis of CIGs, geranyl pyrophosphate synthase (GPPS) catalyzes the formation of geranyl pyrophosphate, which is the direct precursor of CIGs. In this study, the C. officinalis geranyl pyrophosphate synthase (CoGPPS) sequence was cloned from C. officinalis and analyzed. The cDNA sequence of the CoGPPS gene was 915 bp (GenBank No. OR725699). Phylogenetic analysis showed that CoGPPS was closely related to the GPPS sequence of Actinidia chinensis and Camellia sinensis, but relatively distantly related to Paeonia lactiflora and Tripterygium wilfordii. Results from the quantitative real-time PCR showed the spatiotemporal expression pattern of CoGPPS; that is, CoGPPS was specifically expressed in the fruits. Subcellular localization assay proved that CoGPPS was specifically found in chloroplasts. Loganin and morroniside contents in the tissues were detected by high-performance liquid chromatography, and both compounds were found to be at higher levels in the fruits than in leaves. Thus, this study laid the foundation for further studies on the synthetic pathway of CIGs.

Transcriptomic divergence of the Rheum palmatum complex derived from top-geoherb and non-geoherb areas provides the insights into geoherbalism properties of rhubarb

Geoherb usually represents high-quality medicinal herbs with better clinical therapeutic effects, and elucidating the geoherbalism is essential for the quality improvement of traditional Chinese Medicine. However, few researches were conducted to clarify the geoherbalism based on a large scale of transcriptomics. In the present study, we compared the transcriptomes of Rheum palmatum complex derived from top-geoherb and non-geoherb areas to show the geoherbalism properties of rhubarb. A total of 412.32 Gb clean reads were obtained with unigene numbers of 100,615 after assembly. Based on the obtained transcriptome datasets, key enzyme-encoding genes involved in the anthraquinones biosynthesis were also obtained. We also found that 21 anthraquinone-related unigenes were differentially expressed between two different groups, and some of these DEGs were correlated to the content accumulation of five free anthraquinones, indicating that the gene expression profiles may promote the geoherbalism formation of rhubarb. In addition, the selective pressure analyses indicated that most paired orthologous genes between these two groups were subject to negative selection, and only a low proportion of orthologs under positive selection were detected. Functional annotation analyses indicated that these positive-selected genes related to the functions such as gene expression, substance transport, stress response and metabolism, indicating that discrepant environment also enhanced the formation of geoherbalism. Our study not only provided insights for the genetic mechanism of geoherbalism of rhubarb, but also laid more genetic cues for the future rhubarb germplasms improvement and utilization.

Comparative transcriptome analyses of three Gentiana species provides signals for the molecular footprints of selection effects and the phylogenetic relationships

Species in Gentiana section Cruciata are important alpine plants with a center of diversity and speciation in Qinghai-Tibet Plateau (QTP), and some of these species are sympatrically distributed in northeastern QTP. Studies on genome features and natural selection signatures of sympatric species in section Crucata have been impeded by a lack of genomic resources. Here, we showed transcript characterizations and molecular footprints of selection effects on G. straminea, G. dahurica and G. officinalis based on the comparative transcriptome. A total of 62.97 Gb clean reads were obtained with unigene numbers per species ranging from 141,819 to 236,408 after assembly. We found that these three species had similar distribution of functional categories in different databases, and key enzyme-encoding genes involved in the iridoids biosynthesis were also obtained. The selective pressure analyses indicated that most paired orthologs between these three species were subject to negative selection, and only a low proportion of the orthologs that underwent positive selection were detected. We found that some positive selected genes were involved in "catalytic activity", "metabolic process", "response to stimulus" and "response to stress". Besides, large numbers of SSR primer pairs with transferabilities were successfully designed based on the available transcriptome datasets of three Gentiana species. The phylogenetic relationships reconstructed based on 352 single-copy nuclear genes provided a rough phylogenetic framework for this genus and confirmed the monophyly of section Cruciata. Our study not only provides insights for the natural selection effects on sympatric Gentiana species, but also enhances future genetic breeding or evolutionary studies on Qinjiao species.

De Novo Transcriptome Analysis Reveals Flowering-Related Genes That Potentially Contribute to Flowering-Time Control in the Japanese Cultivated Gentian Gentiana triflora

Japanese cultivated gentians are perennial plants that flower in early summer to late autumn in Japan, depending on the cultivar. Several flowering-related genes, including GtFT1 and GtTFL1, are known to be involved in regulating flowering time, but many such genes remain unidentified. In this study, we obtained transcriptome profiling data using the Gentiana triflora cultivar ‘Maciry’, which typically flowers in late July. We conducted deep RNA sequencing analysis using gentian plants grown under natural field conditions for three months before flowering. To investigate diurnal changes, the plants were sampled at 4 h intervals over 24 h. Using these transcriptome data, we determined the expression profiles of leaves based on homology searches against the Flowering-Interactive Database of Arabidopsis. In particular, we focused on transcription factor genes, belonging to the BBX and MADS-box families, and analyzed their developmental and diurnal variation. The expression levels of representative BBX genes were also analyzed under long- and short-day conditions using in-vitro-grown seedlings, and the expression patterns of some BBX genes differed. Clustering analysis revealed that the transcription factor genes were coexpressed with GtFT1. Overall, these expression profiles will facilitate further analysis of the molecular mechanisms underlying the control of flowering time in gentians.

De Novo Transcriptome Analysis Reveals Putative Genes Involved in Anthraquinone Biosynthesis in Rubia yunnanensis

Rubia yunnanensis Diels (R. yunnanensis), a Chinese perennial plant, is well-known for its medicinal values such as rheumatism, contusion, and anemia. It is rich in bioactive anthraquinones, but the biosynthetic pathways of anthraquinones in R. yunnanensis remain unknown. To investigate genes involved in anthraquinone biosynthesis in R. yunnanensis, we generated a de novo transcriptome of R. yunnanensis using the Illumina HiSeq 2500 sequencing platform. A total of 636,198 transcripts were obtained, in which 140,078 transcripts were successfully annotated. A differential gene expression analysis identified 15 putative genes involved in anthraquinone biosynthesis. Additionally, the hairy roots of R. yunnanensis were treated with 200 µM Methyl Jasmonate (MeJA). The contents of six bioactive anthraquinones and gene expression levels of 15 putative genes were measured using ultra performance liquid chromatography coupled with mass spectrometry (UPLC-MS/MS) and real-time quantitative polymerase chain reaction (RT-qPCR), respectively. The results showed that the expressions levels for 11 of the 15 genes and the contents of two of six anthraquinones significantly increased by MeJA treatment. Pearson’s correlation analyses indicated that the expressions of 4 of the 15 putative genes were positively correlated with the contents of rubiquinone (Q3) and rubiquinone-3-O-β-d-xylopranosyl-(1→6)-β-d-glucopyranoside (Q20). This study reported the first de novo transcriptome of R. yunnanensis and shed light on the anthraquinone biosynthesis and genetic information for R. yunnanensis.

The chromosome-level genome assembly of Gentiana dahurica (Gentianaceae) provides insights into gentiopicroside biosynthesis

Gentiana dahurica Fisch. is a perennial herb of the family Gentianaceae. This species is used as a traditional Tibetan medicine because of its rich gentiopicroside constituents. Here, we generate a high-quality, chromosome-level genome of G. dahurica with a total length of 1,416.54 Mb. Comparative genomic analyses showed that G. dahurica shared one whole-genome duplication (WGD) event with Gelsemium sempervirens of the family Gelsemiaceaei and had one additional species-specific WGD after the ancient whole-genome triplication with other eudicots. Further transcriptome analyses identified numerous enzyme coding genes and the transcription factors related to gentiopicroside biosynthesis. A set of candidate cytochrome P450 genes were identified for being involved in biosynthetic shifts from swertiamarin to gentiopicroside. Both gene expressions and the contents measured by high-performance liquid chromatography indicated that the gentiopicrosides were mainly synthesized in the rhizomes with the highest contents. In addition, we found that two above-mentioned WGDs, contributed greatly to the identified candidate genes involving in gentiopicroside biosynthesis. The first reference genome of Gentianaceae we generated here will definitely accelerate evolutionary, ecological, and pharmaceutical studies of this family.

Tissue-specific transcriptome for Rheum tanguticum reveals candidate genes related to the anthraquinones biosynthesis

Rheum tanguticum (Maxim. ex Regel) Maxim. ex Balf. is a herbaceous perennial plant indigenous to China, and its root and rhizomes were usually used as an important traditional Chinese medicine. However, the genomic resources are still scarce for R. tanguticum and even for Rheum genus. Transcriptome datasets from different tissues of R. tanguticum were obtained to screen the genes related to anthraquinones biosynthesis, and five free anthraquinones were also determined. Nine cDNA libraries of roots, stems and leaves were generated, and a total of 272 million high-quality reads were assembled into 257,942 unigenes. Based on the functional annotation, A total of 227 candidate enzyme genes involved in the MVA, MEP, shikimate and polyketide pathways were identified, and several differentially expressed genes found functionally associated with anthraquinones biosynthesis showed distinct tissue-specific expression patterns. Especially, we found that the expression levels of PKS III genes might result in the content differences of free anthraquinones in different tissues of R. tanguticum. Besides, 137,400 SSR loci were identified, and 64,081 SSR primer pairs were successfully designed based on these loci. Our results not only provide cues for the genetic mechanism of anthraquinone content differences in different tissues of R. tanguticum, but also lay genomic foundation for the subsequent genetic engineering and breeding for Rheum species.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01099-8.

Transcriptome Analysis and Microsatellite Markers Development of a Traditional Chinese Medicinal Herb Halenia elliptica D. Don (Gentianaceae)

Halenia elliptica is a popular Chinese medicinal herb that is used to treat jaundice disease and virus hepatitis, and its wild populations have been reduced significantly due to overharvesting recently. However, effective conservation could not be implemented because of the lack of genomic information and genetic markers. In this study, a de novo transcriptome of H elliptica was sequenced using the NGS Illumina, and 132 695 unigenes with the length >200 bp (base pairs) were obtained. Among them, a total of 32 109 unigenes were scanned to develop simple sequence repeats (SSRs). Based on NCBI (National Center for Biotechnology Information) nonredundant database (Nr), these SSR sequences were annotated and assigned into gene ontology categories. In addition, we designed 126 pairs of SSR primers for polymerase chain reaction amplification, of which 12 pairs were identified to be polymorphic among 40 individuals from 8 populations. We then used the 12 polymorphic SSRs to construct a UPGMA dendrogram of the 40 individuals. In addition, a significant correlation between the genetic relationship and the geographic distance was found, suggesting a phylogeographic structure in H elliptica. Moreover, 2 of these SSRs were also successfully amplified in a related species Veratrilla baillonii, suggesting their cross-species transferability. Generally, the SSR markers with high polymorphisms identified in this study provide valuable genetic resources and represent an initial step for exploring the genetic diversity and population histories of H elliptica and its related species.

Immunomodulatory Activity of Ganoderma atrum Polysaccharide on Purified T Lymphocytes through Ca2+/CaN and Mitogen-Activated Protein Kinase Pathway Based on RNA Sequencing

Our previous study has demonstrated that Ganoderma atrum polysaccharide (PSG-1) has immunomodulatory activity on spleen lymphocytes. However, how PSG-1 exerts its effect on purified lymphocytes is still obscure. Thus, this study aimed to investigate the immunomodulatory activity of PSG-1 on purified T lymphocytes and further elucidate the underlying mechanism based on RNA sequencing (RNA-seq). Our results showed that PSG-1 promoted T lymphocytes proliferation and increased the production of IL-2, IFN-γ, and IL-12. Meanwhile, RNA-seq analysis found 394 differentially expressed genes. KEGG pathway analysis identified 20 significant canonical pathways and seven biological functions. Furthermore, PSG-1 elevated intracellular Ca²⁺ concentration and calcineurin (CaN) activity and raised the p-ERK, p-JNK, and p-p38 expression levels. T lymphocytes proliferation and the production of IL-2, IFN-γ, and IL-12 were decreased by the inhibitors of calcium channel and mitogen-activated protein kinases (MAPKs). These results indicated that PSG-1 possesses immunomodulatory activity on purified T lymphocytes, in which Ca²⁺/CaN and MAPK pathways play essential roles.

Using Small RNA-seq Data to Detect siRNA Duplexes Induced by Plant Viruses

Small interfering RNA (siRNA) duplexes are short (usually 21 to 24 bp) double-stranded RNAs (dsRNAs) with several overhanging nucleotides at both 5'- and 3'-ends. It has been found that siRNA duplexes bind the RNA-induced silencing complex (RISC) and cleave the sense strands with endonucleases. In this study, for the first time, we detected siRNA duplexes induced by plant viruses on a large scale using next-generation sequencing (NGS) data. In addition, we used the detected 21 nucleotide (nt) siRNA duplexes with 2 nt overhangs to construct a dataset for future data mining. The analytical results of the features in the detected siRNA duplexes were consistent with those from previous studies. The investigation of siRNA duplexes is useful for a better understanding of the RNA interference (RNAi) mechanism. It can also help to improve the virus detection based on the small RNA sequencing (sRNA-seq) technologies and to rationally design siRNAs for RNAi experiments.

Deep sequencing and transcriptome analyses to identify genes involved in secoiridoid biosynthesis in the Tibetan medicinal plant Swertia mussotii

Swertia mussotii Franch. is an important traditional Tibetan medicinal plant with pharmacological properties effective in the treatment of various ailments including hepatitis. Secoiridoids are the major bioactive compounds in S. mussotii. To better understand the secoiridoid biosynthesis pathway, we generated transcriptome sequences from the root, leaf, stem, and flower tissues, and performed de novo sequence assembly, yielding 98,613 unique transcripts with an N50 of 1,085 bp. Putative functions could be assigned to 35,029 transcripts (35.52%) based on BLAST searches against annotation databases including GO and KEGG. The expression profiles of 39 candidate transcripts encoding the key enzymes for secoiridoid biosynthesis were examined in different S. mussotii tissues, validated by qRT-PCR, and compared with the homologous genes from S. japonica, a species in the same family, unveiling the gene expression, regulation, and conservation of the pathway. The examination of the accumulated levels of three bioactive compounds, sweroside, swertiamarin, and gentiopicroside, revealed their considerable variations in different tissues, with no significant correlation with the expression profiles of key genes in the pathway, suggesting complex biological behaviours in the coordination of metabolite biosynthesis and accumulation. The genomic dataset and analyses presented here lay the foundation for further research on this important medicinal plant.

De Novo transcriptome characterization of Dracaena cambodiana and analysis of genes involved in flavonoid accumulation during formation of dragon's blood

Dragon’s blood is a red resin mainly extracted from Dracaena plants, and has been widely used as a traditional medicine in East and Southeast Asia. The major components of dragon’s blood are flavonoids. Owing to a lack of Dracaena plants genomic information, the flavonoids biosynthesis and regulation in Dracaena plants remain unknown. In this study, three cDNA libraries were constructed from the stems of D. cambodiana after injecting the inducer. Approximately 266.57 million raw sequencing reads were de novo assembled into 198,204 unigenes, of which 34,873 unique sequences were annotated in public protein databases. Many candidate genes involved in flavonoid accumulation were identified. Differential expression analysis identified 20 genes involved in flavonoid biosynthesis, 27 unigenes involved in flavonoid modification and 68 genes involved in flavonoid transport that were up-regulated in the stems of D. cambodiana after injecting the inducer, consistent with the accumulation of flavonoids. Furthermore, we have revealed the differential expression of transcripts encoding for transcription factors (MYB, bHLH and WD40) involved in flavonoid metabolism. These de novo transcriptome data sets provide insights on pathways and molecular regulation of flavonoid biosynthesis and transport, and improve our understanding of molecular mechanisms of dragon’s blood formation in D. cambodiana.

Genome-wide analysis of Dongxiang wild rice (Oryza rufipogon Griff.) to investigate lost/acquired genes during rice domestication

BACKGROUND

It is widely accepted that cultivated rice (Oryza sativa L.) was domesticated from common wild rice (Oryza rufipogon Griff.). Compared to other studies which concentrate on rice origin, this study is to genetically elucidate the substantially phenotypic and physiological changes from wild rice to cultivated rice at the whole genome level.

RESULTS

Instead of comparing two assembled genomes, this study directly compared the Dongxiang wild rice (DXWR) Illumina sequencing reads with the Nipponbare (O. sativa) complete genome without assembly of the DXWR genome. Based on the results from the comparative genomics analysis, structural variations (SVs) between DXWR and Nipponbare were determined to locate deleted genes which could have been acquired by Nipponbare during rice domestication. To overcome the limit of the SV detection, the DXWR transcriptome was also sequenced and compared with the Nipponbare transcriptome to discover the genes which could have been lost in DXWR during domestication. Both 1591 Nipponbare-acquired genes and 206 DXWR-lost transcripts were further analyzed using annotations from multiple sources. The NGS data are available in the NCBI SRA database with ID SRP070627.

CONCLUSIONS

These results help better understanding the domestication from wild rice to cultivated rice at the whole genome level and provide a genomic data resource for rice genetic research or breeding. One finding confirmed transposable elements contribute greatly to the genome evolution from wild rice to cultivated rice. Another finding suggested the photophosphorylation and oxidative phosphorylation system in cultivated rice could have adapted to environmental changes simultaneously during domestication.