RNA sequencing and de novo assembly of Solanum trilobatum leaf transcriptome to identify putative transcripts for major metabolic pathways

Enhancement of production of glycoalkaloids by elicitors along with characterization of gene expression of pathways in Solanum xanthocarpum

Solanum xanthocarpum fruits are used in the treatment of cough, fever, and heart disorders. It possesses antipyretic, hypotensive, antiasthmatic, aphrodisiac and antianaphylactic properties. In the present study, 24 elicitors (both biotic and abiotic) were used to enhance the production of glycoalkaloids in cell cultures of S. xanthocarpum. Four concentrations of elicitors were added into the MS culture medium. The maximum accumulation (5.56-fold higher than control) of demissidine was induced by sodium nitroprusside at 50 mM concentration whereas the highest growth of cell biomass (4.51-fold higher than control) stimulated by systemin at 30 mM concentration. A total of 17 genes of biosynthetic pathways of glycoalkaloids were characterized from the cells of S. xanthocarpum. The greater accumulation of demissidine was confirmed with the expression analysis of 11 key biosynthetic pathway enzymes e.g., acetoacetic-CoA thiolase, 3- hydroxy 3-methyl glutaryl synthase, β-hydroxy β-methylglutaryl CoA reductase, mevalonate kinase, farnesyl diphosphate synthase, squalene synthase, squalene epoxidase, squalene-2,3- epoxide cyclase, cycloartenol synthase, UDP-glucose: solanidine glucosyltransferase and UDP-rhamnose: solanidine rhamno-galactosyl transferase. The maximum expression levels of UDP-rhamnose: solanidine rhamno-galactosyl transferase gene was recorded in this study.

Transcriptome sequencing of medical herb Salvia Rosmarinus (Rosemary) revealed the phenylpropanoid biosynthesis pathway genes and their phylogenetic relationships

BACKGROUND

The Salvia rosmarinus spenn. (rosemary) is considered an economically important ornamental and medicinal plant and is widely utilized in culinary and for treating several diseases. However, the procedure behind synthesizing secondary metabolites-based bioactive compounds at the molecular level in S. rosmarinus is not explored completely.

METHODS AND RESULTS

We performed transcriptomic sequencing of the pooled sample from leaf and stem tissues on the Illumina HiSeqTM X10 platform. The transcriptomics analysis led to the generation of 29,523,608 raw reads, followed by data pre-processing which generated 23,208,592 clean reads, and de novo assembly of S. rosmarinus obtained 166,849 unigenes. Among them, nearly 75.1% of unigenes i.e., 28,757 were interpreted against a non-redundant protein database. The gene ontology-based annotation classified them into 3 main categories and 55 sub-categories, and clusters of orthologous genes annotation categorized them into 23 functional categories. The Kyoto Encyclopedia of Genes and Genomes database-based pathway analysis confirmed the involvement of 13,402 unigenes in 183 biochemical pathways, among these unigenes, 1,186 are involved in the 17 secondary metabolite production pathways. Several key enzymes involved in producing aromatic amino acids and phenylpropanoids were identified from the transcriptome database. Among the identified 48 families of transcription factors from coding unigenes, bHLH, MYB, WRKYs, NAC, C2H2, C3H, and ERF are involved in flavonoids and other secondary metabolites biosynthesis.

CONCLUSION

The phylogenetic analysis revealed the evolutionary relationship between the phenylpropanoid pathway genes of rosemary with other members of Lamiaceae. Our work reveals a new molecular mechanism behind the biosynthesis of phenylpropanoids and their regulation in rosemary plants.

Dynamic transcriptome landscape of pulmonary tissues of rats infected with Paragonimus proliferus

Paragonimiasis (pulmonary fluke disease) is a foodborne parasitic disease caused by trematode infections. Paragonimus proliferus is a characteristic Paragonimus species that was first identified in Yunnan Province of China. No direct evidence has yet proven that P. proliferus can infect humans. However, we previously found that P. proliferus infects and damages rat lung tissues via an unclear mechanism. Here, we infected Sprague Dawley rats with P. proliferus and sequenced their lung transcriptomes at various intervals thereafter. We detected P. proliferus on the surface of rat lung tissues at 7 days post infection. It colonized by attaching and secreting dsRNA and utilized nutrients from the lung tissues for mitosis and meiosis and the dynein arm of lung tissues to develop symmetrical organs. The rats generated different types of immune responses that differed according to the stage of infection. We then analyzed P. proliferus responses to these immune strategies and the genes expressed during each stage of infection. Our findings provide a foundation for developing medical treatments for P. proliferus infection.

Green synthesis, characterization and biological activity of Solanum trilobatum-mediated silver nanoparticles

Biologically inspired synthesis of nanoparticles was found to be more attractive in metal nanoparticle synthesis. The present study reported an in-situ biogenic synthesis of silver nanoparticles (AgNPs) using Solanum trilobatum aqueous leaf extract. On this basis, the aqueous leaf extract of S. trilobatum acted as a reducing agent and stabilizing agent to synthesize highly stable AgNPs at ambient temperature. Eventually, the synthesized and stabilized AgNPs surface plasmon resonance was near 430 nm through a UV–visible (UV–vis) spectrophotometer. Here, the stability of the silver colloids monitored through zeta potential and mean particle size was evaluated through diffraction light scattering (DLF). Further, the average particle size was found to be 27.6 nm and spherical, confirmed with transmission electron microscopy (TEM). Also, colloidal AgNPs and aqueous extract are found to be rich sources of antioxidants and exhibit higher free radical scavenging ability. Thus, efficient inhibition with COX1 and COX2 enzymes and the protective effect with human red blood cell (HRBC) membrane stability showed significant results. These features are promising, suggesting the possibility of the AgNPs to be useful to disease-modifying for treating inflammatory disorders and associated complications.

De novo transcriptome analysis of Justicia adhatoda reveals candidate genes involved in major biosynthetic pathway

BACKGROUND

Justicia adhatoda is an important medicinal plant traditionally used in the Indian system of medicine and the absence of molecular-level studies in this plant hinders its wide use, hence the study was aimed to analyse the genes involved in its various pathways.

METHODS AND RESULTS

The RNA isolated was subjected to Illumina sequencing. De novo assembly was performed using TRINITY software which produced 171,064 transcripts with 55,528 genes and N50 value of 2065 bp, followed by annotation of unigenes against NCBI, KEGG and Gene ontology databases resulted in 105,572 annotated unigenes and 40,288 non-annotated unigenes. A total of 5980 unigenes were mapped to 144 biochemical pathways, including the metabolism and biosynthesis pathways. The pathway analysis revealed the major transcripts involved in the tryptophan biosynthesis with TPM values of 6.0903, 33.6854, 11.527, 1.6959, and 8.1662 for Anthranilate synthase alpha, Anthranilate synthase beta, Arogenate/Prephenate dehydratase, Chorismate synthase and Chorismate mutase, respectively. The qRT-PCR validation of the key enzymes showed up-regulation in mid mature leaf when compared to root and young leaf tissue. A total of 16,154 SSRs were identified from the leaf transcriptome of J. Adhatoda ,which could be helpful in molecular breeding.

CONCLUSIONS

The study aimed at identifying transcripts involved in the tryptophan biosynthesis pathway for its medicinal properties, as it acts as a precursor to the acridone alkaloid biosynthesis with major key enzymes and their validation. This is the first study that reports transcriptome assembly and annotation of J. adhatoda plant.

Transcriptome and de novo analysis of Rosa xanthina f. spontanea in response to cold stress

BACKGROUND

Rose is one of most popular ornamental plants worldwide and is of high economic value and great cultural importance. However, cold damage restricts its planting application in cold areas. To elucidate the metabolic response of rose under low temperature stress, we conducted transcriptome and de novo analysis of Rosa xanthina f. spontanea.

RESULTS

A total of 124,106 unigenes from 9 libraries were generated by de novo assembly, with N50 length was 1470 bp, under 4 °C and - 20 °C stress (23 °C was used as a control). Functional annotation and prediction analyses identified 55,084 unigenes, and 67.72% of these unigenes had significant similarity (BLAST, E ≤ 10^- 5) to those in the public databases. A total of 3031 genes were upregulated and 3891 were downregulated at 4 °C compared with 23 °C, and 867 genes were upregulated and 1763 were downregulated at - 20 °C compared with 23 °C. A total of 468 common DEGs were detected under cold stress, and the matched DEGs were involved in three functional categories: biological process (58.45%), cellular component (11.27%) and molecular function (30.28%). Based on KEGG functional annotations, four pathways were significantly enriched: metabolic pathway, response to plant pathogen interaction (32 genes); starch and sucrose metabolism (21 genes); circadian rhythm plant (8 genes); and photosynthesis antenna proteins (7 genes).

CONCLUSIONS

Our study is the first to report the response to cold stress at the transcriptome level in R. xanthina f. spontanea. The results can help to elucidate the molecular mechanism of cold resistance in rose and provide new insights and candidate genes for genetically enhancing cold stress tolerance.

Transcriptome analysis reveals important candidate gene families related to oligosaccharides biosynthesis in Morindaofficinalis

Morinda officinalis How (MO) is one of the best-known traditional herbs and is widely cultivated in subtropical and tropical areas for many years, especially in southern China. Oligosaccharides are the major constituents in the roots of MO, which is well known for its therapeutic effects with anti-depression, anti-osteoporosis, memory-enhancing, ect. To date, the main gene families that regulate the biosynthetic pathway of MO oligosaccharides metabolism yet have been published. In our study, six cDNA libraries generated from six plants of MO were sequenced utilizing an Illumina HiSeq 4000 platform. Corresponding totals of more than 132.60 million clean reads were obtained from the six libraries and assembled into 25,812 unigenes with an average length of 1288 bp. Moreover, 6036 unigenes were found to be allocated to 26 pathways maps using several public databases, and 2538 differential expression genes (DEGs) were screened. Among them, 25 genes from three families were selected as the mainly candidate genes related to MO oligosaccharides biosynthesis. Then, the expression patterns of six DEGs closely related to MO oligosaccharides biosynthesis were verified by quantitative real-time PCR (qRT-PCR). Besides, the MO was clustered more closely to Coffea arabica of Rubiaceae. In summary, the transcriptomic analysis was used to investigate the differences in expression genes of oligosaccharides biosynthesis, with the notable outcome that several key gene families were closely linked to oligosaccharides biosynthesis.

Genome skimming-based simple sequence repeat (SSR) marker discovery and characterization in Grevillea robusta

Proteaceae, a largely southern hemisphere family consisting of 80 genera distributed in Australia and southern Africa as its centres of greatest diversity, also extends well in northern and southern America. Under this family, Grevillea robusta is a fast-growing species got popularity in farm and avenue plantations. Despite the ecological and economic importance, the species has not yet been investigated for its genetic improvement and genome-based studies. Only a few molecular markers are available for the species or its close relatives, which hinders  genomic and population genetics studies. Genetic markers have been intensively applied for the main strategies in breeding programs, especially for the economically important traits. Hence, it is of utmost priority to develop genomic database resources and species-specific markers for studying quantitative genetics in G. robusta. Given this, the present study aimed to develop de novo genome sequencing, robust microsatellites markers, sequence annotation and their validation in different stands of G. robusta in northern India. Library preparation and sequencing were carried out using Illumina paired-end sequencing technology. Approximately, ten gigabases (Gb) sequence data with 70.87 million raw reads assembled into 425,923 contigs (read mapped to 76.48%) comprising 455 Mb genome size (23 × coverage) generated through genome skimming approach. In total, 9421 simple sequence repeat (SSR) primer pairs were successfully designed from 13,335 microsatellite repeats. Afterward, a subset of 161 primer pairs was randomly selected, synthesized and validated. All the tested primers showed successful amplification but only 13 showed polymorphisms. The polymorphic SSRs were further used to estimate the measures of genetic diversity in 12 genotypes each from the states of Punjab, Haryana, Himachal Pradesh and Uttarakhand. Importantly, the average number of alleles (Na), observed heterozygosity (Ho), expected heterozygosity (He), and the polymorphism information content (PIC) were recorded as 2.69, 0.356, 0.557 and 0.388, respectively. The availability of sequence information and newly developed SSR markers could potentially be used in various genetic analyses and improvements through molecular breeding strategies for G. robusta.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01035-w.

Comparative evaluation of full-length isoform quantification from RNA-Seq

Background

Full-length isoform quantification from RNA-Seq is a key goal in transcriptomics analyses and has been an area of active development since the beginning. The fundamental difficulty stems from the fact that RNA transcripts are long, while RNA-Seq reads are short.

Results

Here we use simulated benchmarking data that reflects many properties of real data, including polymorphisms, intron signal and non-uniform coverage, allowing for systematic comparative analyses of isoform quantification accuracy and its impact on differential expression analysis. Genome, transcriptome and pseudo alignment-based methods are included; and a simple approach is included as a baseline control.

Conclusions

Salmon, kallisto, RSEM, and Cufflinks exhibit the highest accuracy on idealized data, while on more realistic data they do not perform dramatically better than the simple approach. We determine the structural parameters with the greatest impact on quantification accuracy to be length and sequence compression complexity and not so much the number of isoforms. The effect of incomplete annotation on performance is also investigated. Overall, the tested methods show sufficient divergence from the truth to suggest that full-length isoform quantification and isoform level DE should still be employed selectively.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04198-1.

Dynamic transcriptome landscape of Paragonimus proliferus developmental stages in the rat lungs

Paragonimus proliferus, a lung fluke of the genus Paragonimus, was first reported in Yunnan province, China. P. proliferus can infect Sprague-Dawley (SD) rats and cause lung damage, but there is still no direct evidence of human infection. Until now, there has been a lack of studies on P. proliferus parasitism and development in mammalian lung tissue. The aim of this study was to perform transcriptomic profiling of P. proliferus at different developmental stages. SD rats were infected with P. proliferus metacercariae obtained from crabs; worms isolated from the lungs at different time points as well as metacercariae were subjected to whole transcriptome sequencing. Overall, 34,403 transcripts with the total length of 33,223,828 bp, average length of 965 bp, and N50 of 1833 bp were assembled. Comparative analysis indicated that P. proliferus, similar to other Paragonimus spp., expressed genes related to catabolism, whereas P. proliferus-specific transcripts were related to the maintenance of cellular redox homeostasis, sensitivity to bacteria, and immune response. Transcriptional dynamics analysis revealed that genes involved in the regulation of catabolism and apoptosis had stable expression over the P. proliferus life cycle, whereas those involved in development and immune response showed time-dependent changes. High expression of genes associated with immune response corresponded to that of genes regulating the sensitivity to bacteria and immune protection. We constructed a P. proliferus developmental model, including the development of the body, suckers, blood cells, reproductive and tracheal systems, lymph, skin, cartilage, and other tissues and organs, and an immune response model, which mainly involved T cells and macrophages. Our study provides a foundation for further research into the molecular biology and infection mechanism of P. proliferus.

Reference transcriptomes and comparative analyses of six species in the threatened rosewood genus Dalbergia

Dalbergia is a pantropical genus with more than 250 species, many of which are highly threatened due to overexploitation for their rosewood timber, along with general deforestation. Many Dalbergia species have received international attention for conservation, but the lack of genomic resources for Dalbergia hinders evolutionary studies and conservation applications, which are important for adaptive management. This study produced the first reference transcriptomes for 6 Dalbergia species with different geographical origins and predicted ~ 32 to 49 K unique genes. We showed the utility of these transcriptomes by phylogenomic analyses with other Fabaceae species, estimating the divergence time of extant Dalbergia species to ~ 14.78 MYA. We detected over-representation in 13 Pfam terms including HSP, ALDH and ubiquitin families in Dalbergia. We also compared the gene families of geographically co-occurring D. cochinchinensis and D. oliveri and observed that more genes underwent positive selection and there were more diverged disease resistance proteins in the more widely distributed D. oliveri, consistent with reports that it occupies a wider ecological niche and has higher genetic diversity. We anticipate that the reference transcriptomes will facilitate future population genomics and gene-environment association studies on Dalbergia, as well as contributing to the genomic database where plants, particularly threatened ones, are currently underrepresented.

De novo transcriptome analysis of Lantana camara L. revealed candidate genes involved in phenylpropanoid biosynthesis pathway

Lantana camara L. is an economically important essential oil producing plant belonging to family Verbenaceae. It is used in medication for treating various diseases like cancer, ulcers, tumor, asthma and fever. The plant is a useful source of essential bioactive compounds such as steroids, flavonoids and phenylpropanoid glycosides etc. Nonetheless, very little is known about the genomic or transcriptomic resources of L. camara, and this might be the reason of hindering molecular studies leading to identification of improved lines. Here we used Illumina sequencing platform and performed the L. camara leaf (LCL) and root (LCR) de novo transcriptome analyses. A total of 70,155,594 and 84,263,224 clean reads were obtained and de novo assembly generated 72,877 and 513,985 unigenes from leaf (LCL) and root (LCR) respectively. Furthermore, the pathway analysis revealed the presence of 229 and 943 genes involved in the phenylpropanoid biosynthesis in leaf and root tissues respectively. Similarity search was performed against publically available genome databases and best matches were found with Sesamum indicum (67.5%) that were much higher than that of Arabidopsis thaliana (3.9%). To the best of our knowledge, this is the first comprehensive transcriptomic analysis of leaf and root tissues of this non-model plant from family Verbenaceae and may serve as a baseline for further molecular studies.

Transcriptomic Features of Echinococcus granulosus Protoscolex during the Encystation Process

Cystic echinococcosis (CE) is a zoonotic infection caused by Echinococcus granulosus larvae. It seriously affects the development of animal husbandry and endangers human health. Due to a poor understanding of the cystic fluid formation pathway, there is currently a lack of innovative methods for the prevention and treatment of CE. In this study, the protoscoleces (PSCs) in the encystation process were analyzed by high-throughput RNA sequencing. A total of 32,401 transcripts and 14,903 cDNAs revealed numbers of new genes and transcripts, stage-specific genes, and differently expressed genes. Genes encoding proteins involved in signaling pathways, such as putative G-protein coupled receptor, tyrosine kinases, and serine/threonine protein kinase, were predominantly up-regulated during the encystation process. Antioxidant enzymes included cytochrome c oxidase, thioredoxin glutathione, and glutathione peroxidase were a high expression level. Intriguingly, KEGG enrichment suggested that differentially up-regulated genes involved in the vasopressin-regulated water reabsorption metabolic pathway may play important roles in the transport of proteins, carbohydrates, and other substances. These results provide valuable information on the mechanism of cystic fluid production during the encystation process, and provide a basis for further studies on the molecular mechanisms of growth and development of PSCs.

De novo transcriptome of Gymnema sylvestre identified putative lncRNA and genes regulating terpenoid biosynthesis pathway

Gymnema sylvestre is a highly valuable medicinal plant in traditional Indian system of medicine and used in many polyherbal formulations especially in treating diabetes. However, the lack of genomic resources has impeded its research at molecular level. The present study investigated functional gene profile of G. sylvestre via RNA sequencing technology. The de novo assembly of 88.9 million high quality reads yielded 23,126 unigenes, of which 18116 were annotated against databases such as NCBI nr database, gene ontology (GO), KEGG, Pfam, CDD, PlantTFcat, UniProt & GreeNC. Total 808 unigenes mapped to 78 different Transcription Factor families, whereas 39 unigenes assigned to CYP450 and 111 unigenes coding for enzymes involved in the biosynthesis of terpenoids including transcripts for synthesis of important compounds like Vitamin E, beta-amyrin and squalene. Among them, presence of six important enzyme coding transcripts were validated using qRT-PCR, which showed high expression of enzymes involved in methyl-erythritol phosphate (MEP) pathway. This study also revealed 1428 simple sequence repeats (SSRs), which may aid in molecular breeding studies. Besides this, 8 putative long non-coding RNAs (lncRNAs) were predicted from un-annotated sequences, which may hold key role in regulation of essential biological processes in G. sylvestre. The study provides an opportunity for future functional genomic studies and to uncover functions of the lncRNAs in G. sylvestre.