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

Novel gurmarin-like peptides from Gymnema sylvestre and their interactions with the sweet taste receptor T1R2/T1R3

Gymnema sylvestre (GS) is a traditional medicinal plant known for its hypoglycemic and hypolipidemic effects. Gurmarin (hereafter Gur-1) is the only known active peptide in GS. Gur-1 has a suppressive sweet taste effect in rodents but no or only a very weak effect in humans. Here, 8 gurmarin-like peptides (Gur-2 to Gur-9) and their isoforms are reported in the GS transcriptome. The molecular mechanism of sweet taste suppression by Gur-1 is still largely unknown. Therefore, the complete architecture of human and mouse sweet taste receptors T1R2/T1R3 and their interaction with Gur-1 to Gur-9 were predicted by AlphaFold-Multimer (AF-M) and validated. Only Gur-1 and Gur-2 interact with the T1R2/T1R3 receptor. Indeed, Gur-1 and Gur-2 bind to the region of the cysteine-rich domain (CRD) and the transmembrane domain (TMD) of the mouse T1R2 subunit. In contrast, only Gur-2 binds to the TMD of the human T1R2 subunit. This result suggests that Gur-2 may have a suppressive sweet taste effect in humans. Furthermore, AF-M predicted that Gα-gustducin, a protein involved in sweet taste transduction, interacts with the intracellular domain of the T1R2 subunit. These results highlight an unexpected diversity of gurmarin-like peptides in GS and provide the complete predicted architecture of the human and mouse sweet taste receptor with the putative binding sites of Gur-1, Gur-2, and Gα-gustducin. In addition, gurmarin-like peptides may serve as promising drug scaffolds for the development of antidiabetic molecules.

Identification of a putative novel cholivirus in the transcriptome of Gymnema sylvestre R. Br

Gymnema sylvestre is a tropical climber species that is widely used in traditional medicine since ages. In the present study, the transcriptome datasets of G. sylvestre available in public domain were screened for the presence of novel plant viral sequences and a putative novel virus tentatively named as Gymnema sylvestre virus 1 (GysV1) was identified. Coding-complete genome segments of GysV1 that are 6.35 kb (RNA1) and 3.98 kb (RNA2) long possessed a single large open reading frame coding for a polyprotein. BLASTp, sequence identity and phylogenetic analyses revealed the relatedness of GysV1 to the members of the subgenus Cholivirus (genus Sadwavirus; family Secoviridae; order Picornavirales). Based on the species demarcation criteria of the family Secoviridae, GysV1 can be regarded as a new cholivirus member.

Characterization of squalene synthase gene from Gymnema sylvestre R. Br

Background

Squalene synthase (SQS) is a rate-limiting enzyme necessary to produce pentacyclic triterpenes in plants. It is an important enzyme producing squalene molecules required to run steroidal and triterpenoid biosynthesis pathways working in competitive inhibition mode. Reports are available on information pertaining to SQS gene in several plants, but detailed information on SQS gene in Gymnema sylvestre R. Br. is not available. G. sylvestre is a priceless rare vine of central eco-region known for its medicinally important triterpenoids. Our work aims to characterize the GS-SQS gene in this high-value medicinal plant.

Results

Coding DNA sequences (CDS) with 1245 bp length representing GS-SQS gene predicted from transcriptome data in G. sylvestre was used for further characterization. The SWISS protein structure modeled for the GS-SQS amino acid sequence data had MolProbity Score of 1.44 and the Clash Score 3.86. The quality estimates and statistical score of Ramachandran plots analysis indicated that the homology model was reliable. For full-length amplification of the gene, primers designed from flanking regions of CDS encoding GS-SQS were used to get amplification against genomic DNA as template which resulted in approximately 6.2-kb sized single-band product. The sequencing of this product through NGS was carried out generating 2.32 Gb data and 3347 number of scaffolds with N50 value of 457 bp. These scaffolds were compared to identify similarity with other SQS genes as well as the GS-SQSs of the transcriptome. Scaffold_3347 representing the GS-SQS gene harbored two introns of 101 and 164 bp size. Both these intronic regions were validated by primers designed from adjoining outside regions of the introns on the scaffold representing GS-SQS gene. The amplification took place when the template was genomic DNA and failed when the template was cDNA confirmed the presence of two introns in GS-SQS gene in Gymnema sylvestre R. Br.

Conclusion

This study shows GS-SQS gene was very closely related to Coffea arabica and Gardenia jasminoides and this gene harbored two introns of 101 and 164 bp size.

Transcriptome analysis for molecular landscaping of genes controlling diterpene andrographolide biosynthesis in Andrographis paniculata (Burm . f.) Nees

Kalmegh [Andrographis paniculata (Burm. f.) Nees.] is one of the essential medicinal plants due to an important terpenoid, i.e. andrographolide possesses immense therapeutic and pharmacological uses. The experiment was performed to elucidate the expression of candidate genes associated with andrographolide biosynthesis. Based on results obtained in chromatography for andrographolide content analysis of six genotypes, two contrast genotypes, i.e. IC-520361 (maximum andrographolide content-2.33%) and Anand Local (lowest andrographolide content-1.01%) were selected for the transcriptome analysis. A total of 1.04 Gb of raw data were produced using MiSeq Illumina platform, in which IC 520361 generated 645 million base pairs sequence along with 4,524,251 raw reads and Anand Local produced 419 million base pairs sequence along with 3,021,316 raw reads. The combined assembly of high quality reads generated for both the samples had 33,247,454 bp of total assembled bases and 38,292 of transcripts. The GC percent of assembled transcripts was 44.79%, an average read length was 800 bp and N50 value was 1186 bp. Species-specific distribution using BLAST X (Nr), showed the highest Blast hits with Sesamum indicum. Out of 23,346 transcripts, 87% of transcripts annotated in UniProt KB (Universal Protein Resource KnowledgeBase) database and only 0.21% of transcripts were annotated in TAIR (The Arabidopsis Information Resources). Biological processes gene ontology classified based on Blast2GO showed, out of 6853 transcripts, 1370 of transcripts were represented by terpenoid biosynthetic pathway, which involved in secondary metabolite andrographolide biosynthesis. The heat map showed 1016 transcripts were differentially expressed between two kalmegh genotypes, in which nine important differentially expressed transcripts related to MEP (2C methyl-d-erythritol 4-phosphate) and MVA (Mevalonic acid) andrographolide biosynthesis pathways such as, geranyl diphosphate synthase small subunit, Isopentenyl-diphosphate delta-isomerase i-like, 4, 13-hydroxy-3-methylglutaryl-coenzyme a reductase etc. were upregulated in IC 520361 as compared to Anand Local, which were validated through RT-qPCR. The highest expression of gene 13-hydroxy-3-methylglutaryl-coenzyme a reductase (HMGR) was reported, which is responsible for accumulation of andrographolide in leaf. This comparative transcriptome analysis confirmed the expression level of genes were higher in accession IC 520361 as compare to Anand Local related to andrographolide biosynthesis pathways i.e. MEP and MVA. These up-regulated genes could be over-expressed to enhance the andrographolide content using genetic engineering of these metabolic pathways. It will also give an idea to the breeder for development of molecular markers for direct screening of the genotypes.

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.