De Novo Transcriptome Analysis of Medicinally Important Plantago ovata Using RNA-Seq

A chromosome-level genome assembly of Plantago ovata

Plantago ovata is cultivated for production of its seed husk (psyllium). When wet, the husk transforms into a mucilage with properties suitable for pharmaceutical industries, utilised in supplements for controlling blood cholesterol levels, and food industries for making gluten-free products. There has been limited success in improving husk quantity and quality through breeding approaches, partly due to the lack of a reference genome. Here we constructed the first chromosome-scale reference assembly of P. ovata using a combination of 5.98 million PacBio and 636.5 million Hi-C reads. We also used corrected PacBio reads to estimate genome size and transcripts to generate gene models. The final assembly covers ~ 500 Mb with 99.3% gene set completeness. A total of 97% of the sequences are anchored to four chromosomes with an N50 of ~ 128.87 Mb. The P. ovata genome contains 61.90% repeats, where 40.04% are long terminal repeats. We identified 41,820 protein-coding genes, 411 non-coding RNAs, 108 ribosomal RNAs, and 1295 transfer RNAs. This genome will provide a resource for plant breeding programs to, for example, reduce agronomic constraints such as seed shattering, increase psyllium yield and quality, and overcome crop disease susceptibility.

Integrative Application of Transcriptomics and Metabolomics Provides Insights into Unsynchronized Growth in Sea Cucumber (Stichopus monotuberculatus)

Ever-increasing consumer demand for sea cucumbers mainly leads to huge damage to wild sea cucumber resources, including Stichopus monotuberculatus, which in turn exerts negative impacts on marine environments due to the lack of ecological functions performed by sea cucumbers. Aquaculture of sea cucumbers is an effective way to meet consumer demand and restore their resources. Unsynchronous growth is a prominent problem in the aquaculture of sea cucumbers which has concealed unelucidated molecular mechanisms until now. In this study, we carried out an integrative analysis of transcriptomics and metabolomics on fast-growing (SMF) and slow-growing (SMS) groups of S. monotuberculatus cultured in the same environmental conditions. The results revealed that a total of 2054 significantly differentially expressed genes (DEGs) were identified, which are mainly involved in fat digestion and absorption, histidine metabolism, arachidonic acid metabolism, and glutathione metabolism. 368 differential metabolites (DMs) were screened out between the SMF group and the SMS group; these metabolites are mainly involved in glycerophospholipid metabolism, purine metabolism, biosynthesis of unsaturated fatty acids, pyrimidine metabolism, arachidonic acid metabolism, and other metabolic pathways. The integrative analysis of transcriptomics and metabolomics of S. monotuberculatus suggested that the SMF group had a higher capacity for lipid metabolism and protein synthesis, and had a more frequent occurrence of apoptosis events, which are likely to be related to coping with environmental stresses. The results of this study provide potential values for the aquaculture of sea cucumbers which may promote their resource enhancement.

The comparisons of expression pattern reveal molecular regulation of fruit metabolites in S. nigrum and S. lycopersicum

Solanum nigrum, known as black nightshade, is a medicinal plant that contains many beneficial metabolites in its fruit. The molecular mechanisms underlying the synthesis of these metabolites remain uninvestigated due to limited genetic information. Here, we identified 47,470 unigenes of S. nigrum from three different tissues by de novo transcriptome assembly, and 78.4% of these genes were functionally annotated. Moreover, gene ontology (GO) analysis using 18,860 differentially expressed genes (DEGs) revealed tissue-specific gene expression regulation. We compared gene expression patterns between S. nigrum and tomato (S. lycopersicum) in three tissue types. The expression patterns of carotenoid biosynthetic genes were different between the two species. Comparison of the expression patterns of flavonoid biosynthetic genes showed that 9 out of 14 enzyme-coding genes were highly upregulated in the fruit of S. nigrum. Using CRISPR-Cas9-mediated gene editing, we knocked out the R2R3-MYB transcription factor SnAN2 gene, an ortholog of S. lycopersicum ANTHOCYANIN 2. The mutants showed yellow/green fruits, suggesting that SnAN2 plays a major role in anthocyanin synthesis in S. nigrum. This study revealed the connection between gene expression regulation and corresponding phenotypic differences through comparative analysis between two closely related species and provided genetic resources for S. nigrum.

Comparative transcriptome mining for terpenoid biosynthetic pathway genes in wild and cultivated species of Plantago

Plantagos are important economical and medicinal plants that possess several bioactive secondary metabolites, such as phenolics, iridoids, triterpenes, and alkaloids. Triterpenoids are the ubiquitous and dynamic secondary metabolites that are deployed by plants for chemical interactions and protection under biotic/abiotic stress. Plantago ovata, a cultivated species, is the source of psyllium, while Plantago major, a wild species, has significant therapeutic potential. Wild species are considered more tolerant to stressful conditions in comparison to their cultivated allies. In view of this, the present study aimed to decipher the terpenoid biosynthetic pathway operative in P. ovata and P. major using a comparative transcriptomics approach. Majority of terpenoid biosynthetic genes were observed as upregulated in P. major including rate limiting genes of MVA (HMGR) and MEP (DXR) pathways and genes (α-AS, BAS, SM, and CYP716) involved in ursolic acid biosynthesis, an important triterpenoid prevalent in Plantago species. The HPLC output further confirmed the higher concentration of ursolic acid in P. major as compared to P. ovata leaf samples, respectively. In addition to terpenoid biosynthesis, KEGG annotation revealed the involvement of differentially expressed unigenes in several metabolic pathways, aminoacyl-tRNA biosynthesis, biosynthesis of antibiotics, and biosynthesis of secondary metabolites. MYB was found as the most abundant transcription factor family in Plantago transcriptome. We have been able to generate valuable information which can help in improving terpenoid production in Plantago. Additionally, the present study has laid a strong foundation for deciphering other important metabolic pathways in Plantago.

De novo transcriptome and tissue specific expression analysis of genes associated with biosynthesis of secondary metabolites in Operculina turpethum (L.)

This study reported the first-ever de novo transcriptome analysis of Operculina turpethum, a high valued endangered medicinal plant, using the Illumina HiSeq 2500 platform. The de novo assembly generated a total of 64,259 unigenes and 20,870 CDS (coding sequence) with a mean length of 449 bp and 571 bp respectively. Further, 20,218 and 16,458 unigenes showed significant similarity with identified proteins of NR (non-redundant) and UniProt database respectively. The homology search carried out against publicly available database found the best match with Ipomoea nil sequences (82.6%). The KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analysis identified 6538 unigenes functionally assigned to 378 modules with phenylpropanoid biosynthesis pathway as the most enriched among the secondary metabolite biosynthesis pathway followed by terpenoid biosynthesis. A total of 17,444 DEGs were identified among which majority of the DEGs (Differentially Expressed Gene) involved in secondary metabolite biosynthesis were found to be significantly upregulated in stem as compared to root tissues. The qRT-PCR validation of 9 unigenes involved in phenylpropanoid and terpenoid biosynthesis also showed a similar expression pattern. This finding suggests that stem tissues, rather than root tissues, could be used to prevent uprooting of O. turpethum in the wild, paving the way for the plant's effective conservation. Moreover, the study formed a valuable repository of genetic information which will provide a baseline for further molecular research.

Extensive Metabolite Profiling in the Unexploited Organs of Black Tiger for Their Potential Valorization in the Pharmaceutical Industry

Black tiger (Kadsura coccinea (Lem.)) has been reported to hold enormous pharmaceutical potential. The fruit and rhizome of black tiger are highly exploited in the pharmaceutical and other industries. However, the most important organs from the plant such as the leaf and stem are considered biowastes mainly because a comprehensive metabolite profile has not been reported in these organs. Knowledge of the metabolic landscape of the unexploited black tiger organs could help identify and isolate important compounds with pharmaceutical and nutritional values for a better valorization of the species. In this study, we used a widely targeted metabolomics approach to profile the metabolomes of the K. coccinea leaf (KL) and stem (KS) and compared them with the root (KR). We identified 642, 650 and 619 diverse metabolites in KL, KS and KR, respectively. A total of 555 metabolites were mutually detected among the three organs, indicating that the leaf and stem organs may also hold potential for medicinal, nutritional and industrial applications. Most of the differentially accumulated metabolites between organs were enriched in flavone and flavonol biosynthesis, phenylpropanoid biosynthesis, arginine and proline metabolism, arginine biosynthesis, tyrosine metabolism and 2-oxocarboxylic acid metabolism pathways. In addition, several important organ-specific metabolites were detected in K. coccinea. In conclusion, we provide extensive metabolic information to stimulate black tiger leaf and stem valorization in human healthcare and food.

Full-Length Transcriptome Sequencing Provides Insights into Flavonoid Biosynthesis in Fritillaria hupehensis

One of the most commonly utilized medicinal plants in China is Fritillaria hupehensis (Hsiao et K.C. Hsia). However, due to a lack of genomic resources, little is known about the biosynthesis of relevant compounds, particularly the flavonoid biosynthesis pathway. A PacBio RS II sequencing generated a total of 342,044 reads from the bulb, leaf, root, and stem, of which 316,438 were full-length (FL) non-redundant reads with an average length of 1365 bp and a N50 of 1888 bp. There were also 38,607 long non-coding RNAs and 7914 simple sequence repeats detected. To improve our understanding of processes implicated in regulating secondary metabolite biosynthesis in F. hupehensis tissues, we evaluated potential metabolic pathways. Overall, this study provides a repertoire of FL transcripts in F. hupehensis for the first time, and it will be a valuable resource for marker-assisted breeding and research into bioactive compounds for medicinal and pharmacological applications.

De Novo assembly, characterization and development of EST-SSRs from Bletilla striata transcriptomes profiled throughout the whole growing period

Bletilla striata is an endangered orchid that has been used for millennia as a medicinal herb, in cosmetics and as a horticultural plant. To construct the first nucleotide database for this species and to develop abundant EST-SSR markers for facilitating further studies, various tissues and organs of plants in the main developmental stages were harvested for mRNA isolation and subsequent RNA sequencing. A total of 106,054,784 clean reads were generated by using Illumina paired-end sequencing technology. The reads were assembled into 127,261 unigenes by the Trinity package; the unigenes had an average length of 612 bp and an N50 of 957 bp. Of these unigenes, 67,494 (51.86%) were annotated in a series of databases. Of these annotated unigenes, 41,818 and 24,615 were assigned to gene ontology categories and clusters of orthologous groups, respectively. Additionally, 20,764 (15.96%) unigenes were mapped onto 275 pathways using the KEGG database. In addition, 25,935 high-quality EST-SSR primer pairs were developed from the 15,433 unigenes by MISA mining. To validate the accuracy of the newly designed markers, 87 of 100 randomly selected primers were effectively amplified; 63 of those yielded PCR products of the expected size, and 25 yielded products with significant amounts of polymorphism among the 4 landraces. Furthermore, the transferability test of the 25 polymorphic markers was performed in 6 individuals of two closely related genus Phalaenopsis and dendrobium. Which results showed a total of 5 markers can successfully amplified among these populations. This research provides a comprehensive nucleotide database and lays a solid foundation for functional gene mining and genomic research in B. striata. The developed EST-SSR primers could facilitate phylogenetic studies and breeding.

Identification and characterization of some putative genes involved in arabinoxylan biosynthesis in Plantago ovata

Plantago ovata is an important source of Psyllium (Isabgol), which swells upon contact with water forming mucilaginous mass, largely composed of arabinoxylans. In this study, we analyzed the expression pattern of arabinoxylan biosynthetic pathway genes at different stages of seed development in P. ovata. Besides, arabinoxylans were quantified at different stages of seed development in water extractable and water unextractable fractions. The expression analysis revealed 5-8 fold increase in the levels of expression of some genes involved in arabinoxylan biosynthetic pathway such as UDP-arabinopyranose mutase, UDP-xylosyltransferase 2 and xylan glucuronosyltransferase at 15 days after pollination stage in seed. The xylose and arabinose units were analyzed at different stages of seed development and also in water-soluble (cold water and hot water), alkali and ethanolic fractions. The concentration of xylose and arabinose units increased steadily after pollination. Overall, alkali extract had high concentration of xylose (0.70 ± 0.022 mg/g) and arabinose units (0.10 ± 0.01 mg/g) at 15 days after pollination stage.

Enrichment of genomic resources and identification of simple sequence repeats from medicinally important Clausena excavata

To broaden and delve into the genomic information of Clausena excavata, an important medicinal plant in many Asian countries, RNA sequencing (RNA-seq) analysis was performed and a total of 16,638 non-redundant unigenes (≥ 300 bp) with an average length of 755 bp were generated by de novo assembly from 17,580,456 trimmed clear reads. The functional categorization of the identified unigenes by a gene ontology (GO) term resulted in 2305 genes in the cellular component, 5577 in the biological processes, and 8056 in the molecular functions, respectively. The top sub-category in biological processes was the metabolic process with 4374 genes. Among annotated genes, 3006 were mapped to 123 metabolic pathways by the Kyoto Encyclopedia of Genes and Genomes (KEGG) metabolic pathway analysis tool. The search for simple sequence repeats (SSRs) resulted in 845 SSRs from 749 SSR-containing unigenes and the most abundant SSR motifs was AAG/CTT with 179 occurrences. Twelve SSR markers were tested for cross transferability among five Clausena species; eight of them exhibited polymorphism. Taken together, these data provide valuable resources for genomic or genetic studies of Clausena species and other relative studies. The transcriptome shotgun assembly data have been deposited at DDBJ/EMBL/GenBank under the accession GGEM00000000.

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.

Transcriptomic resources and marker validation for diploid and polyploid Veronica (Plantaginaceae) from New Zealand and Europe

PREMISE OF THE STUDY

Polyploidy may generate novel variation, leading to adaptation and species diversification. An excellent natural system to study polyploid evolution in a comparative framework is Veronica (Plantaginaceae), which comprises several parallel, recently evolved polyploid series.

METHODS

Over 105 million Illumina paired-end sequence reads were generated from cDNA libraries of leaf tissue from eight individuals representing three European and four New Zealand species. Forty-eight simple sequence repeat (SSR) and 48 low-copy nuclear (LCN) markers were developed and validated with Fluidigm microfluidic PCR and Illumina MiSeq amplicon sequencing on 48 different individuals each.

RESULTS

Individual Trinity assemblies were similar regarding annotated transcripts (13,009-14,271), mean contig length (635-742 bp), N50 value (916-1133 bp), E90N50 value (1099-1308 bp), contigs with positive BLAST hits (42-63%), and gene ontology terms. Analyses of 29,738 single-nucleotide polymorphisms (8746 phylogenetically informative) mined from these transcriptomes plus two outgroups (Picrorhiza kurrooa and Plantago ovata) showed moderate to high bootstrap support for all branches and reticulation among sampled European Veronica.

DISCUSSION

The transcriptome sequences themselves, as well as the validated SSR (40/48) and LCN (11/48) markers derived from them, show inter- and intraspecific genetic variation. These resources will be invaluable for future population genetic, phylogenetic, and functional genetic investigations in polyploid Veronica.