De novo transcriptome sequencing facilitates genomic resource generation in Tinospora cordifolia

Genome sequencing and functional analysis of a multipurpose medicinal herb Tinospora cordifolia (Giloy)

Tinospora cordifolia (Willd.) Hook.f. & Thomson, also known as Giloy, is among the most important medicinal plants that have numerous therapeutic applications in human health due to the production of a diverse array of secondary metabolites. To gain genomic insights into the medicinal properties of T. cordifolia, the genome sequencing was carried out using 10× Genomics linked read and Nanopore long-read technologies. The draft genome assembly of T. cordifolia was comprised of 1.01 Gbp, which is the genome sequenced from the plant family Menispermaceae. We also performed the genome size estimation for T. cordifolia, which was found to be 1.13 Gbp. The deep sequencing of transcriptome from the leaf tissue was also performed. The genome and transcriptome assemblies were used to construct the gene set, resulting in 17,245 coding gene sequences. Further, the phylogenetic position of T. cordifolia was also positioned as basal eudicot by constructing a genome-wide phylogenetic tree using multiple species. Further, a comprehensive comparative evolutionary analysis of gene families contraction/expansion and multiple signatures of adaptive evolution was performed. The genes involved in benzyl iso-quinoline alkaloid, terpenoid, lignin and flavonoid biosynthesis pathways were found with signatures of adaptive evolution. These evolutionary adaptations in genes provide genomic insights into the presence of diverse medicinal properties of this plant. The genes involved in the common symbiosis signalling pathway associated with endosymbiosis (Arbuscular Mycorrhiza) were found to be adaptively evolved. The genes involved in adventitious root formation, peroxisome biogenesis, biosynthesis of phytohormones, and tolerance against abiotic and biotic stresses were also found to be adaptively evolved in T. cordifolia.

De novo genome assembly and annotation of the medicinal plant Tinospora cordifolia (Willd.) Miers ex Hook. f. & Thom's

Indian natural climbing shrub Tinospora cordifolia, often known as "Guduchi" and "Amrita," is a highly esteemed medicinal plant in the Indian system of medicine (ISM). It is a member of the Menispermaceae family which consists of a rich source of protein, micronutrients, and rich source of bioactive components which are used in treating various systemic diseases. The current study was designed to know the biological characterization of the plant genome and biosynthesis of plant metabolites essential for its medicinal applications. Tinospora cordifolia's complete genome was sequenced using Illumina HiSeq2500 sequencing technology. The draft genome was assembled through a de novo method. An integrative genome annotation approach was used to perform functional gene prediction. The pathway analysis was carried out using the KEGG database. The total genome size obtained after genome assembly was 894 Mb with an N50 of 9148 bp. The integrative annotation approach resulted in 35,111 protein-coding genes. In addition, genes responsible for the synthesis of syringin, a secondary metabolite found in plants, were identified. In comparison to the standard drug (dopamine, rasagiline, and selegiline), syringin's molecular docking exhibited a greater binding affinity from the range of - 4.3 to - 6.6 kcal/mol for all the targets of Parkinson's disease and for Alzheimer's targets; it has shown the maximum potency from the range of - 6.5 to - 7.4 kcal/mol with respect to the standard drug (donepezil, galantamine, and rivastigmine). This study provides the genomic information of Tinospora cordifolia which is helpful in understanding genomic insights and metabolic pathways connected to the corresponding plant genome and predicts the possible useful effect for the molecular characterization of therapeutic drugs.

De Novo Transcriptome Profiling for the Generation and Validation of Microsatellite Markers, Transcription Factors, and Database Development for Andrographis paniculata

Andrographis paniculata belongs to the family Acanthaceae and is known for its medicinal properties owing to the presence of unique constituents belonging to the lactones, diterpenoids, diterpene glycosides, flavonoids, and flavonoid glycosides groups of chemicals. Andrographolide, a major therapeutic constituent of A. paniculata, is extracted primarily from the leaves of this plant and exhibits antimicrobial and anti-inflammatory activities. Using 454 GS-FLX pyrosequencing, we have generated a whole transcriptome profile of entire leaves of A. paniculata. A total of 22,402 high-quality transcripts were generated, with an average transcript length and N50 of 884 bp and 1007 bp, respectively. Functional annotation revealed that 19,264 (86%) of the total transcripts showed significant similarity with the NCBI-Nr database and were successfully annotated. Out of the 19,264 BLAST hits, 17,623 transcripts were assigned GO terms and distributed into three major functional categories: molecular function (44.62%), biological processes (29.19%), and cellular component (26.18%) based on BLAST2GO. Transcription factor analysis showed 6669 transcripts, belonging to 57 different transcription factor families. Fifteen TF genes that belong to the NAC, MYB, and bHLH TF categories were validated by RT PCR amplification. In silico analysis of gene families involved in the synthesis of biochemical compounds having medicinal values, such as cytochrome p450, protein kinases, heat shock proteins, and transporters, was completed and a total of 102 different transcripts encoding enzymes involved in the biosynthesis of terpenoids were predicted. Out of these, 33 transcripts belonged to terpenoid backbone biosynthesis. This study also identified 4254 EST-SSRs from 3661 transcripts, representing 16.34% of the total transcripts. Fifty-three novel EST-SSR markers generated from our EST dataset were used to assess the genetic diversity among eighteen A. paniculata accessions. The genetic diversity analysis revealed two distinct sub-clusters and all accessions based on the genetic similarity index were distinct from each other. A database based on EST transcripts, EST-SSR markers, and transcription factors has been developed using data generated from the present study combined with available transcriptomic resources from a public database using Meta transcriptome analysis to make genomic resources available in one place to the researchers working on this medicinal plant.

Molecular Characterization of Tinospora cordifolia (Willd.) Miers Using Novel g-SSR Markers and Their Comparison with EST-SSR and SCoT Markers for Genetic Diversity Study

In the present study, novel genomic-SSR (g-SSR) markers generated in our laboratory were used to characterize Tinospora cordifolia and related species. The g-SSR marker was also compared with EST-SSR and SCoT markers used earlier in our laboratory to assess the genetic diversity of T. cordifolia. A total of 26 accessions of T. cordifolia and 1 accession each of Tinospora rumphii and Tinospora sinensis were characterized using 65 novel g-SSR markers. A total of 125 alleles were detected with 49 polymorphic g-SSR markers. The number of alleles per locus varied from 1–4 with a mean value of 2.55 alleles per locus. Mean PIC, gene diversity, and heterozygosity were estimated to be 0.33, 0.41, and 0.65, respectively. The two species, namely T. rumphii and T. sinensis, showed cross-species transferability of g-SSRs developed in T. cordifolia. The success rate of cross-species transferability in T. rumphii was 95.3% and 93.8% in T. sinensis, proving the usefulness of this marker in genetic diversity studies of related species. The Tinospora accessions were also used for molecular characterization using SCoT and EST-SSR markers and compared for genetic diversity and cross-species transferability. The PIC, gene diversity, heterozygosity, and principal coordinate analysis showed that g-SSR is the better maker for a genetic diversity study of T. cordifolia. Additionally, high cross-species transferability of g-SSRs was found (95.3% and 93.8%) compared to EST-SSRs (68.8% and 67.7%) in T. rumphii and T. sinensis, respectively.

TinoTranscriptDB: A Database of Transcripts and Microsatellite Markers of Tinospora cordifolia, an Important Medicinal Plant

Tinospora cordifolia, commonly known as “Giloe” in India, is a shrub belonging to the family Menispermaceae. It is an important medicinal plant known for its antipyretic, anti-inflammatory, antispasmodic, and antidiabetic properties and is used in the treatment of jaundice, gout, and rheumatism. Despite its economic importance, the limited information related to its genomic resources prohibits its judicious exploitation through molecular breeding or biotechnological approaches. In this study, we generated a meta-transcriptome assembly of 43,090 non-redundant transcripts by merging the RNASeq data obtained from Roche 454 GS-FLX, and Illumina platforms, and report the first transcriptome-based database for simple sequence repeats and transcription factors (“TinoTranscriptDB” (Tinospora cordifolia Transcriptome Database)). We annotated 26,716 (62%) of the total transcripts successfully from National Center for Biotechnology Information non-redundant protein (NCBI-NR), gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), Swiss-Prot, and Pfam databases. This database contains information of 2620 perfect simple sequence repeats (P-SSRs) with a relative abundance of 340.12 (loci/Mb), and relative density of 6309.29 (bp/Mb). Excluding mono-nucleotides, the most abundant SSR motifs were tri-nucleotides (54.31%), followed by di-nucleotides (37.51%), tetra-nucleotides (4.54%), penta-nucleotides (3.16%) and hexa-nucleotides (0.45%). Additionally, we also identified 4,311 transcription factors (TFs) and categorized them into 55 sub-families. This database is expected to fill the gap in genomic resource availability in T. cordifolia and thus accelerate molecular breeding and related functional and other applied studies aimed towards genetic improvements of T. cordifolia and related species.

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.

Estimation of genetic diversity and population structure in Tinospora cordifolia using SSR markers

Thirty polymorphic SSRs, derived from RNA sequencing of Tinospora cordifolia (willd.), were utilized for genetic diversity and population structure evaluation among 96 accessions collected from ten different geographical regions of India. A total of 7611 SSRs were identified from 268149 transcripts. Of all SSR loci, 69.07% of them were tri-nucleotide repeat motifs, followed by di-nucleotide repeat motifs (12.82%). A total of 230 alleles were generated by 30 SSRs with an average of 7.67 alleles per locus with comparatively higher polymorphic information content (average 0.68). The expected (He) and observed (Ho) heterozygosity means were 0.71 and 0.12, respectively. All the loci showed significant deviation from Hardy-Weinberg Equilibrium (HWE). The neighbor joining clustering based on jaccard's coefficient grouped all the 96 accessions into three major cluster which was also in congruence with model-based structure plot. The result of molecular variance (AMOVA) revealed higher genetic variance within populations than among populations. The result reflects an existence of high level of genetic diversity in the collected accessions of T. cordifolia. The accessions Tc131, Tc31, Tc129, Tc38, Tc16, Tc59, Tc60, Tc17, Tc106 and Tc130 was found to be potential and diverse in nature and the SSRs TCSSR-18, TCSSR-37, TCTSSR-59, TCTSSR-92, TCTSSR-123 and TCTSSR-126 as potential markers. These accessions and newly developed SSR markers provide valuable resource and could be strategically utilized for further genetic improvement of T. cordifolia.

Tinospora cordifolia as a potential neuroregenerative candidate against glutamate induced excitotoxicity: an in vitro perspective

BACKGROUND

Glutamate, the major excitatory neurotransmitter of CNS acts as a neurotoxin at higher concentrations. Prolonged activation of glutamate receptors results in progressive neuronal damage by aggravating calcium influx, inducing mitochondrial damage and oxidative stress. Excitotoxic cell death is associated with the pathogenesis of various neurodegenerative disorders such as trauma, brain injury and neurodegenerative diseases. The current study was designed to investigate the neuroprotective and neuroregenerative potential of Tinospora cordifolia against glutamate-induced excitotoxicity using primary cerebellar neuronal cultures as a model system.

METHODS

Monosodium salt of glutamate was used to induce neurotoxic injury in primary cerebellar neurons. Four extracts including Hexane extract, Chloroform extract, Ethyl acetate, and Butanol extract were obtained from fractionation of previously reported aqueous ethanolic extract of T. cordifolia and tested for neuroprotective activity. Out of the four fractions, Butanol extract of T. cordifolia (B-TCE) exhibited neuroprotective potential by preventing degeneration of neurons induced by glutamate. Expression of different neuronal, apoptotic, inflammatory, cell cycle regulatory and plasticity markers was studied by immunostaining and Western blotting. Neurite outgrowth and migration were also studied using primary explant cultures, wound scratch and gelatin zymogram assay.

RESULTS

At molecular level, B-TCE pretreatment of glutamate-treated cultures normalized the stress-induced downregulation in the expression of neuronal markers (MAP-2, GAP-43, NF200) and anti-apoptotic marker (Bcl-xL). Further, cells exposed to glutamate showed enhanced expression of inflammatory (NF-κB, AP-1) and senescence markers (HSP70, Mortalin) as well as the extent of mitochondrial damage. However, B-TCE pretreatment prevented this increase and inhibited glutamate-induced onset of inflammation, stress and mitochondrial membrane damage. Furthermore, B-TCE was observed to promote regeneration, migration and plasticity of cerebellar neurons, which was otherwise significantly inhibited by glutamate treatment.

CONCLUSION

These results suggest that B-TCE may have neuroprotective and neuroregenerative potential against catastrophic consequences of glutamate-mediated excitotoxicity and could be a potential therapeutic candidate for neurodegenerative diseases.

Comparative transcriptome analysis to identify putative genes involved in thymol biosynthesis pathway in medicinal plant Trachyspermum ammi L

Thymol, as a dietary monoterpene, is a phenol derivative of cymene, which is the major component of the essential oil of Trachyspermum ammi (L.). It shows multiple biological activities: antifungal, antibacterial, antivirus and anti-inflammatory. T. ammi, commonly known as ajowan, belongs to Apiaceae and is an important medicinal seed spice. To identify the putative genes involved in thymol and other monoterpene biosynthesis, we provided transcriptomes of four inflorescence tissues of two ajowan ecotypes, containing different thymol yield. This study has detected the genes encoding enzymes for the go-between stages of the terpenoid biosynthesis pathways. A large number of unigenes, differentially expressed between four inflorescence tissues of two ajowan ecotypes, was revealed by a transcriptome analysis. Furthermore, differentially expressed unigenes encoding dehydrogenases, transcription factors, and cytochrome P450s, which might be associated with terpenoid diversity in T. ammi, were identified. The sequencing data obtained in this study formed a valuable repository of genetic information for an understanding of the formation of the main constituents of ajowan essential oil and functional analysis of thymol-specific genes. Comparative transcriptome analysis led to the development of new resources for a functional breeding of ajowan.

Candidate genes of flavonoid biosynthesis in Selaginella bryopteris (L.) Baker identified by RNA-Seq

In the present study, de novo transcriptome analysis of Selaginella bryopteris in frond and root was performed to understand the regulation of flavonoid (FL) biosynthesis. High-quality data of 5.84 and 5.86 Gb was generated for frond and root, respectively, that assembled into 94,713 and 81,567 transcripts. A total of 87,471 and 73,395 unigenes were obtained from frond and root, respectively. A total of 41,267 and 31,048 CDS of frond and root, respectively, were annotated by BLASTX, which showed maximum hits against S. moellendorffii. Out of 11,285 differentially expressed genes, a total of 5639 genes were found to be down-regulated and 5628 genes up-regulated in frond as compared to those in root. In silico analysis of expression of genes in frond as compared to that in root was done for those related to phenylpropanoid (PP)/FL biosynthesis along with transcription factors (TFs) after DESeq and MapMan-based information. Results showed that genes of PP/FL biosynthesis pathway namely SbCHS, SbCHI, SbF3H, SbF3'H, SbDFR, SbUF3GT, SbCCOAMT, and SbCATOMT and TFs (SbMYB1, SbMYB2, SbMYB3, SbBHLH1, and SbWD40-5) were up-regulated in frond in comparison to those in root. Further, this in silico expression data was validated by RT-PCR analysis which showed predominant expression of most of these genes in frond and indicated their importance in the biosynthesis of flavonoids in S. bryopteris. A total of 9074 simple sequence repeats (SSRs) were also identified for frond and 3811 SSRs for root; these can be used for experimental validation.

Exploring drought stress-regulated genes in senna (Cassia angustifolia Vahl.): a transcriptomic approach

De novo assembly of reads produced by next-generation sequencing (NGS) technologies offers a rapid approach to obtain expressed gene sequences for non-model organisms. Senna (Cassia angustifolia Vahl.) is a drought-tolerant annual undershrub of Caesalpiniaceae, a subfamily of Fabaceae. There are insufficient transcriptomic and genomic data in public databases for understanding the molecular mechanism underlying the drought tolerance of senna. Therefore, the main purpose of this study was to know the transcriptome profile of senna, with special reference to drought stress. RNA from two different stages of leaf development was extracted and sequenced separately using the Illumina technology. A total of 200 million reads were generated, and a de novo assembly of processed reads in the pooled transcriptome using Trinity yielded 43,413 transcripts which were further annotated using NCBI BLAST with "green plant database (txid 33090)," Swiss Prot, Kyoto Encyclopedia of Genes and Genomes (KEGG), Clusters of Orthologous Groups (COG), and Gene Ontology (GO). Out of the total transcripts, 42,280 (95.0 %) were annotated by BLASTX against the green plant database of NCBI. Senna transcriptome showed the highest similarity to Glycine max (41 %), followed by Phaseolus vulgaris (16 %), Cicer arietinum (15 %), and Medicago trancatula (5 %). The highest number of GO terms were enriched for the molecular functions category; of these "catalytic activity" (GO: 0003824) (25.10 %) and "binding activity" (GO: 0005488) (20.10 %) were most abundantly represented. We used InterProscan to see protein similarity at domain level; a total of 33,256 transcripts were annotated against the Pfam domains. The transcripts were assigned with various KEGG pathways. Coding DNA sequences (CDS) encoding various drought stress-regulated pathways such as signaling factors, protein-modifying/degrading enzymes, biosynthesis of phytohormone, phytohormone signaling, osmotically active compounds, free radical scavengers, chlorophyll metabolism, leaf cuticular wax, polyamines, and protective proteins were identified through BLASTX search. The lucine-rich repeat kinase family was the most abundantly found group of protein kinases. Orphan, bHLH, and bZIP family TFs were the most abundantly found in senna. Six genes encoding MYC2 transcription factor, 9-cis-epoxycarotenoid dioxygenase (NCED), l -ascorbate peroxidase (APX), aminocyclopropane carboxylate oxidase (ACO), abscisic acid 8'-hydroxylase (ABA), and WRKY transcription factor were confirmed through reverse transcriptase-PCR (RT-PCR) and Sanger sequencing for the first time in senna. The potential drought stress-related transcripts identified in this study provide a good start for further investigation into the drought adaptation in senna. Additionally, our transcriptome sequences are the valuable resource for accelerated genomics-assisted genetic improvement programs and facilitate manipulation of biochemical pathways for developing drought-tolerant genotypes of crop plants.