Exploring the medicinally important secondary metabolites landscape through the lens of transcriptome data in fenugreek (Trigonella foenum graecum L.) | Scientific Reports – Nature.com

I am sending you 1π! Pi is a new digital currency developed by Stanford PhDs, with over 30 million members worldwide. To claim your Pi, follow this link https://minepi.com/Bigtoff and use my username (Bigtoff) as your invitation code.

Exploring the medicinally important secondary metabolites landscape through the lens of transcriptome data in fenugreek (Trigonella foenum graecum L.) | Scientific Reports – Nature.com

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Advertisement
Scientific Reports volume 12, Article number: 13534 (2022)
Metrics details
Fenugreek (Trigonella foenum-graecum L.) is a self-pollinated leguminous crop belonging to the Fabaceae family. It is a multipurpose crop used as herb, spice, vegetable and forage. It is a traditional medicinal plant in India attributed with several nutritional and medicinal properties including antidiabetic and anticancer. We have performed a combined transcriptome assembly from RNA sequencing data derived from leaf, stem and root tissues. Around 209,831 transcripts were deciphered from the assembly of 92% completeness and an N50 of 1382 bases. Whilst secondary metabolites of medicinal value, such as trigonelline, diosgenin, 4-hydroxyisoleucine and quercetin, are distributed in several tissues, we report transcripts that bear sequence signatures of enzymes involved in the biosynthesis of such metabolites and are highly expressed in leaves, stem and roots. One of the antidiabetic alkaloid, trigonelline and its biosynthesising enzyme, is highly abundant in leaves. These findings are of value to nutritional and the pharmaceutical industry.
Fenugreek (Trigonella foenum-graecum L.) is a dicotyledonous leguminous self-pollinated crop with chromosome number 2n = 16 belonging to the Fabaceae (Leguminosae) family1. It is native to India and eastern Mediterranean countries. The genus Trigonella includes 140 species distributed globally, among which 107 species occur in Asia2. In India mainly two types are commonly cultivated i.e., common fenugreek (Trigonella foenum-graecum L.) for leaves, stem and seeds, and Kasuri fenugreek (Trigonella corniculata L.) for leaves. Fenugreek is one of the oldest known medicinal plants of India and is traditionally used in ayurvedic medicines. It is considered to be a multipurpose crop for its food and non-food uses such as herb, spice, fodder, and its nutraceutical, pharmaceutical and therapeutic properties3. In 2017, Venkata et al.4 reviewed fenugreek for its broad pharmaceutical usage on preclinical and clinical research on human diseases and antipathogenic properties.
Although there has been considerable advancement in the field of drug development, medicinal plants are being increasingly utilised for treating various diseases owing to their therapeutic properties and safety. Fenugreek has a high content of dietary fibre, vitamin and mineral content as well as bioactive compounds such as modified amino acids, phenolic acids, alkaloids and saponins (Fig. 1). It has multiple health benefits, such as antidiabetic, anticancer, galactagogic, helps in digestion, hepatoprotective effects, regulatory functions, antioxidant properties, and works against anorexia, antilithogenic, hepatoprotective effect, antipathogenic properties and several other medicinal properties5,6,7. It is also Generally Recognized As Safe (GRAS) while using as a flavour by the U.S. Food and Drug Administration (FDA)8.
A snapshot of selected secondary metabolites produced by Trigonella foenum-graecum (shown in image is the AFG-1 variety). The abundance of these compounds make fenugreek a repertoire of nutraceutical and pharmaceutical properties.
A recent review reports the presence of more than 100 phytochemicals in the fenugreek seeds9. Among the several important secondary metabolites produced by fenugreek, diosgenin, trigonelline and 4-hydroxyisoleucine are three compounds which have been implicated in the anti-diabetic properties of the plant (Fig. 1). Diosgenin, a steroidal saponin, is an important starting material for the preparation of several steroidal drugs in the pharmaceutical industry and has shown high potential in the treatment of various types of disorders such as cancer, hypercholesterolemia, inflammation, and several types of infections10. The pharmacological activities of trigonelline, an alkaloid, include hypoglycaemic, hypolipidemic, neuroprotective, antimigraine, sedative, memory-improving, antibacterial, antiviral, and anti-tumour activities11. 4-hydroxyisoleucine (4-HIL), a modified amino acid is another bioactive compound abundant in fenugreek seeds. The antidiabetic property of 4-HIL is attributed to its ability to induce insulin secretion in a glucose-dependent manner12,13. There are several other important classes of metabolites such as volatiles (eugenol, linalool), flavonoids (quercetin, rutin, vitexin and isovitexin) which render fenugreek multiple medicinal properties7,14.
There have been a limited number of studies for identification of potential genes involved in biosynthesis of bioactive compounds attributing the nutritional and pharmaceutical properties of the fenugreek as well as differential expression of these genes across the tissues of the plant. The transcriptome studies in fenugreek will pave the way for identification of candidate genes involved in medicinal as well as environmental stress tolerance properties of the plant and these can be taken up further in molecular breeding programs15,16. Transcriptomics provide the details of genes associated with functions at the cellular and tissue level; analysis of transcriptomic data from different tissues of the organism is particularly useful in understanding differential gene expression17. Towards this goal, we have performed transcriptome sequencing from leaf, stem and root tissues of fenugreek and estimated the transcript abundance across these tissues. We use this data to explore the secondary metabolite landscape of the plant by deciphering genes involved in synthesis of several nutritionally and medicinally important metabolites. We have also quantified the expression of selected genes and performed quantification of the selected metabolites.
Seeds of T. foenum-graecum (Ajmer Fenugreek-1 (AFG-1) variety) were procured from College of Horticulture, GKVK Campus, Bengaluru with due permission. All the studies on T. foenum-graecum were carried out in accordance with relevant institutional, national and international guidelines and legislation. For transcriptome sequencing, RNA from three different tissues-leaf, stem and root, was isolated from one month old plants using the Plant Spectrum Plant RNA extraction kit (Sigma-Aldrich). These tissue samples, with two biological replicates, were sequenced using Illumina HiSeq 1000 platform with 100 base pairs (bp) read length. Trimmomatic was used to process the reads using default settings18. A total of 226.5 million reads were obtained after quality processing from six libraries. These reads were assembled de novo using Trinity with default parameters19. BUSCO was used to assess the completeness of the transcriptome assembly20.
TransDecoder (v3.0.0) was used to identify the candidate coding regions within transcript sequences obtained from Trinity19. A length cut-off of 100 bases was used to predict Open Reading Frames (ORFs). Function assignment for amino acid sequences was obtained using BLASTP21 against UniProt Viridiplantae database with an E-value cut-off of 10–522 and domain identification was performed using hmmscan of the HMMER suite (v.3.1b2)23 against Pfam (v.31.0) using an E-value cut-off of 10–524. For transcripts, BLASTX21 was performed against the above sequence database with the default parameters. Additional annotation was rendered using SignalP (v.4.0)25, TMHMM (v.2.0)26 and RNAmmer (v.1.2)27 to extract information on signal peptides, transmembrane regions and RNA sequences. Gene ontology (GO) terms28, orthologous relationships (eggNOG)29 and pathway information (KEGG)30 were also mined.
Orthology analysis of T. foenum-graecum was carried out with the well-studied model plants Arabidopsis thaliana and Oryza sativa and plants which belong to the same tribe as fenugreek (Trifolieae) viz. Trifolium pratense and Medicago truncatula using OrthoMCL31. An E-value cut-off of 10–5 was used for the all-against-all BLAST runs. Additionally, 34 different plant species proteomes from Phytozome (v.12)32 were used for large scale orthology analysis using Proteinortho (v5.11)33 (Supplementary Table 1) with an E-value cut-off of 10–5.
The concatenated nucleotide sequence alignments for ribulose-1,5-bisphosphate carboxylase/oxygenase large subunit (rbcL), maturase K (matK) and internal transcribed spacer 1 and 2 (ITS1 and ITS2) were used for constructing species phylogenetic tree for the species belonging to subclass Rosids. In case of partial sequences, the alignment was trimmed. The individual sequences were aligned using MAFFT34 with maximum allowance for aligning gapped regions and the separate alignments were concatenated. Phylogenetic reconstruction was performed using the Neighbour-Joining method implemented in MAFFT with 100 bootstraps. The orthologue distribution obtained from Proteinortho analysis along with their phylogenetic tree for comparison amongst these closely related species was plotted. Through orthogroup analysis, the in-paralogous groups and singletons of T. foenum-graecum common to both methods were identified. The sequences were compared with each other with respect to their expression profiles in different tissues, Pfam domains22 and GO terms28.
The reads from the three tissues were mapped onto the de novo transcriptome assembly using Bowtie2 (v2.3.0). The expression counts were obtained using eXpress (v1.5.1) and RSEM35,36,37. Transcripts per Million (TPM) values were calculated for all the transcripts. The average TPM value was considered for the biological replicates and compared across the three tissues.
The genes involved in synthesis of secondary metabolites belonging to alkaloids, saponins, volatiles and flavonoids, present in T. foenum-graecum were identified using strategy adapted from Joshi et al.38. The pathway information and related functional annotation was extracted from multiple sources such as PlantCyc MetaCyc, KEGG and literature11,30,39,40,41,42. The pathway images were adapted from PlantCyc or literature. The sequences for the enzymes involved in the synthesis of secondary metabolites serving as start points for sequence search were identified from UniProt22. This was followed with creation and curation of multiple sequence alignment profiles including functionally important residues (FIR), typically used for selecting trusted homologues21,43. Such profiles were implemented to search the transcriptome data for potential hits and were ascertained by presence of the FIRs. The phylogenetic trees were constructed using the neighbor-joining method from MEGA (v10)44 for 1000 bootstraps, to assess co-clustering of the candidate genes with the curated start points. The phylogenetic trees obtained were visualized and edited in Figtree (v1.4.2)45 and iToL46.
The genes identified to be involved in production of secondary metabolites trigonelline and diosgenin were validated using quantitative reverse transcription PCR (qRT-PCR). The details of methods used for RNA isolation, qRT-PCR quantification are described in the Supplementary Text (Sect. 1.1).
The quantification of diosgenin, trigonelline and 4-hydroxyisoleucine was carried out from leaf, stem, root, and seed tissue samples using HLPC-PDA method. The details of this method are described in the methods part of Supplementary Text (Sect. 1.2).
We carried out transcriptome sequencing of three tissues (leaf, stem and root) of fenugreek plant and constructed a combined assembly along with functional annotation. We investigated the orthology relationship with other plant species and also estimated the transcript abundance across tissues. Further, we explored the secondary metabolite landscape of fenugreek considering metabolites with reported medicinal properties from several classes of compounds (alkaloids, saponins, volatiles, flavonoids, etc.) and mined for the enzymes involved in their synthesis.
Combined transcriptome assembly was obtained from three different tissues (leaf, stem and root) with N50 of 1382 bases. The 169.5 Mb assembly contained 209,831 transcripts. This assembly was 92.2% complete (95.49% using partial core genes) as assessed through the BUSCO20.
The candidate coding regions within the fenugreek transcriptome were predicted using Transdecoder program. A set of 64,670 ORFs (complete: 37,294, partial: 27,376—terminally partial: 18,743, internally partial: 8633) were predicted from the assembly which varied from 200 to 2000 base pairs. We obtained annotations, using BLASTX against UniProt Viridiplantae database, for 26% of transcripts (Table 1). Whereas, among the predicted ORFs, we obtained annotation for 90% of ORFs, through BLASTP using E-value threshold of 10–5 and 50% query coverage cut-off. Most of these ORFs (76%) were associated with Medicago truncatula, a closely related species, as the best hit.
Among the predicted ORFs, 66% (42,524) could be associated with Pfam domains, representing 4400 unique Pfam domains. The leucine-rich repeat (LRR), pentatricopeptide repeat (PPR), tetratricopeptide repeats (TPR) and several transcription factor binding (WD40, F-box, AAA) domain families were the most abundant annotated domains observed. Similar to other angiosperms, several other Pfam domains such as, protein serine/threonine kinases, protein tyrosine kinases, calcium binding EF-hand and ankyrin repeats were found in abundance (Supplementary Fig. 1A). We could associate GO terms to more than 50% of total ORFs. The top enriched GO terms were classified into molecular function (MF), biological process (BP) and cellular component (CC) (Supplementary Fig. 1B). The enriched terms under MF included transporter activity, transcription regulator, antioxidant activity and nutrient reservoir. Whereas, some of the most enriched BP terms were metabolic process, response to stimuli and immune system process. It was observed that transmembrane helices were present in 19% of ORFs, signal peptides were predicted in 6% of ORFs, whereas eight rRNAs were identified. Of the total number of amino acid sequences obtained from the transcriptome assembly, 59% were mapped to functional annotation databases, namely eggNOG (orthogroups), GO and KEGG pathways (Table 1; complete annotation data available in Supplementary Data 1).
The orthologs were detected using OrthoMCL and Proteinortho for fenugreek and closely related plants (Fig. 2A,B and Supplementary Fig. 2). Among the five plants used for OrthoMCL analysis, 8822 orthogroups were common to all whereas, for Proteinortho, 114 orthogroups were common across 39 plants considered for the analysis. It was observed that the orthogroup distribution was highly similar across all the plants from subfamily Faboideae (G. max, P. vulgaris, T. pratense, M. truncatula and T. foenum-graecum). We constructed the phylogenetic tree using selected plants from subclass Rosids and O. sativa. (Fig. 2B).
Orthology analyses (A) OrthoMCL: Venn diagram representing the common and unique orthogroups for the five proteomes studied. (B) Proteinortho results showing the distribution of orthogroups between few selected plants from the subclass Rosids (including Arabidopsis thaliana and few plants from Faboideae and Rosaceae) and Oryza sativa. Each coloured section represents the number of species with which the orthogroups are shared.
We observed 20,230 fenugreek specific paragroups and singletons, among which only the longest ORFs (15,712, 24% of all predicted gene products) were assessed for their functional annotation (hereafter referred to as singletons). A total of 6300 of the singletons in fenugreek proteins show presence of 1771 unique Pfam families. PPR, LRR, protein kinases, protein tyrosine kinases, zinc knuckle and reverse transcriptase domain were among the most abundant Pfam domains in the singleton proteins (Supplementary Fig. 3A). Almost half of the fenugreek singletons (7235) were found to have Gene Ontology (GO) annotation. Nucleotide/nucleic acid binding, metal ion binding were abundant in molecular functions within GO terms (Supplementary Fig. 3B). Translation, regulation of transcription and defence response were the most abundant biological process (Supplementary Fig. 3C). The most common cellular component term within these GO terms was ‘integral component of membrane’ (Supplementary Fig. 3D).
The transcript abundance across all three tissues (leaf, stem and root) was assessed. For each tissue, the average TPM value between the replicates was considered for each transcript. The top 20 highly abundant transcripts were identified from each tissue and are represented as a heatmap (Fig. 3). While observing variation between highly abundant transcripts in each tissue, we could identify six transcripts (protease inhibitor, metallothionein, dehydrin B, thioredoxin H, plant invertase inhibitor and translationally controlled tumor-like protein), which were highly expressed across three tissues. Among the top most abundant transcripts in the leaf sample, transcripts mostly related to photosynthesis, for example, ribulose bisphosphate carboxylase, chlorophyll a–b binding protein, carbonic anhydrase etc., were observed. In the stem, apart from photosynthesis related transcripts, few cytochrome genes and pathogenesis related protein families were among the most abundant transcripts. In the root, the top three highly abundant annotated transcripts included metallothionein, dehydrin B and Nmr-A like family of genes. These genes are well documented to be involved in drought tolerance in crops47,48,49,50. The tissue-wide average TPM values for enzymes involved in synthesis of several secondary metabolites considered in this work are mentioned in the subsequent sections and documented in the Supplementary Table 3.
Heatmap for top 20 most abundant transcripts (based on average TPM values) from leaf, stem and root tissues. The * sign indicates transcripts with multiple annotations. The annotation derived from the best hit have been ascribed for such transcripts.
We explored the secondary metabolites landscape in fenugreek with an aim to identify the transcripts involved in their biosynthesis. The candidate metabolites and the enzymes involved in their biosynthesis were chosen following an extensive literature survey. These spanned several classes of secondary metabolites: alkaloids, saponins, volatiles and flavonoids.
Alkaloids are heterocyclic nitrogen containing compounds displaying a wide array of pharmacological properties. These compounds are not only associated with plant growth and defence but are a storage source of nitrogen51. We considered trigonelline and glycine-betaine from the alkaloid class found in fenugreek. We pursued identification of enzymes involved in their biosynthesis by adopting the strategy described in Joshi et al.38.
Trigonelline is a derivative of nicotinic acid, first isolated from fenugreek (Trigonella foenum-graecum). Trigonelline has several medicinal properties such as hypoglycaemic, neuroprotective, antibacterial, and anti-tumour activities52. It is synthesised through purine nucleotide cycling, wherein, nicotinate is methylated to methyl nicotinate (trigonelline) by nicotinate N-methyltransferase (NNMT) (EC: 2.1.1.7). We selected and curated protein sequences for this enzyme found in other plants, as a start point to carry out sequence search in fenugreek transcriptome. There are five functionally important residues (FIRs), a substrate-binding motif Asn 21, Tyr 120, His 124 a catalytic motif Thr 264 and a SAM binding motif in the NNMT enzyme42. Although 12 hits were obtained, a single hit (Tfoe_c19814_g2_i1_m7872) was identified based on co-clustering with Medicago truncatula sequence (NCBI RefSeq ID: XP_013464471.1) (Supplementary Fig. 4E) with a sequence identity of 92.6% and a 100% query coverage. Further this hit was validated through mapping of the FIRs (Supplementary Fig. 4A–D). Based on the average TPM values, this transcript was most abundant in the leaf tissue sample (Fig. 4, Supplementary Table 3). This trend was also similar to the expression level estimation through qRT-PCR (Supplementary Fig. 5A). Trigonelline, as a compound, was quantified across different tissues of fenugreek (described further in “Quantification of metabolites in different tissue samples”).
Trigonelline biosynthesis pathway and relative abundance across tissues of the transcript encoding the enzyme-Nicotinate N-methyltransferase (represented as heatmap of average TPM value in log2 scale, in green-white colour scale). The pathway has been adapted from Plant Metabolic Network (PMN), https://pmn.plantcyc.org/PLANT/NEW-IMAGE?type=PATHWAY&object=PWY-5110&detail-level=3, on www.plantcyc.org, Aug 12, 2021.
Glycine betaine is one of the major osmoprotectants in plants and also contributes to maintaining membrane and cell stability. It plays a protective role against inflammation, carcinogenesis and neurodegenerative diseases53,54,55. It was identified in the extract of seeds of fenugreek56. It is synthesized in plants by the oxidation of choline through a two-step process (Supplementary Data 2). The first step is catalyzed by choline monooxygenase (CHMO) (EC: 1.14.15.7) which is an unusual iron-sulfur containing enzyme. We identified one hit (Tfoe_c38951_g1_i2_m30577) from fenugreek which contained both the motifs (Rieske 2Fe-2S and non-heme binding region). The second step catalyzes the formation of glycine-betaine from betaine-aldehyde by betaine aldehyde dehydrogenase (BADH) (EC: 1.2.1.8)57. We identified one hit for BADH (Tfoe_c82300_g1_i1_m72031) and mapped the FIRs (proton acceptor: Glu260, Nucleophile: Cys294 and cation-pi interacting Trp285, Trp456). The hits for both the enzymes in fenugreek co-clustered with the M. truncatula proteins (Supplementary Data 2). Both these transcripts were relatively abundant in the leaf tissue sample (Supplementary Table 3).
Saponins are amphipathic glycosides grouped phenomenologically by the soap-like foaming they produce when shaken in aqueous solutions. These metabolites are found in the plant families-Sapindaceae, Aceraceae, Hippocastanaceae and Fabaceae58. We focussed on the biosynthetic pathway of diosgenin, a steroidal sapogenin attributed to treatment of metabolic disorders such as diabetes and hypercholesterolemia and also implicated in treatment of cancer, inflammation and several infections10. Tissue-wide diosgenin content varies in fenugreek plant and is reported to be present in highest quantity in seeds. Moreover, its quantity differs across several cultivars and is considered one of the criterion in crop improvement programs59.
From the proposed biosynthetic pathway of diosgenin, as per Ciura et al.60, (Supplementary Data 2), putative sequences corresponding to the last two enzymes (sterol-3β-glucosyltransferase and β-glucosidase) were identified from the set of protein sequences obtained from the transcriptome assembly. Sterol-3β-glucosyltransferase (EC: 2.4.1.173) contains two motifs-putative steroid binding domain (PSBD) motif and plant secondary product glucosyltransferase (PSPG) motif61. The hits obtained from the searches against the proteins sequences were filtered for the presence of these two motifs and of the four hits that contained this motif, from the FIR analysis, two transcripts (Tfoe_c26467_g1_i1_m12030 and Tfoe_c33087_g1_i3_m19456) expressing significantly in all three tissues of the plant were considered as true hits (referred henceforth as sterol-3β-glucosyltransferase-1 and -2). While sterol-3β-glucosyltransferase-1 (Tfoe_c26467_g1_i1_m12030) clustered with a sequence from a closely related plant, M. truncatula, the sterol-3β-glucosyltransferase-2 (Tfoe_c33087_g1_i3_m19456) clustered with the enzyme sequence from Glycine soja. Sterol-3β-glucosyltransferase-1 expressed well in all tissues of the plant (leaf, stem and root) while sterol-3β-glucosyltransferase-2 showed relatively higher average TPM value in leaves and this was further validated using qRT-PCR (Supplementary Fig. 5B,C, Supplementary Table 3).
The last enzyme in the biosynthetic pathway of diosgenin, β-glucosidase (EC: 3.2.1.21), is known to contain two catalytic glutamate residues which are part of the (I/V)TENG and TFNEP motifs and these have been identified from crystal structures of β-glucosidase from rice in apo and bound forms (2RGL and 2RGM)62. We identified the β-glucosidase in fenugreek (Tfoe_c120410_g1_i1_m78951) and confirmed the presence of these motifs. Although, we found a total of 22 such hits which had > 70% query coverage and > 30% sequence identity, with the query β-glucosidase sequences from Viridiplantae, the selected candidate protein Tfoe_c120410_g1_i1_m78951 clustered with the enzyme from M. truncatula and was thus considered a true hit. This transcript showed high average TPM value in root followed by stem and the same relative expression pattern was observed from qRT-PCR validation (Supplementary Fig. 5D, Supplementary Table 3). The metabolic pathway, sequence motif, phylogenetic tree and the expression data for all diosgenin synthesising enzymes is documented in Supplementary Data 2.
The volatiles are a class of compounds which typically include aromatic organic compounds. These compounds are used in perfumes, flavours, topical formulations, etc. The medicinal properties of essential volatiles have been known since ancient times. We have explored the enzymes involved in the biosynthesis of eugenol and linalool which are volatile oils found in fenugreek63,64. The therapeutic potential of eugenol and isoeugenol, which are phenolic compounds, is attributed to antioxidant potency and anti-inflammatory benefits65,66. Linalool, a terpene, has been reported to show anticancer, analgesic, anxiolytic and other neuroprotective properties67. Eugenol and isoeugenol are phenylpropene volatiles that have a phenyl ring with a propenyl side chain synthesized from coniferyl alcohol precursor in multiple enzymatic reactions (Supplementary Data 2). We studied three enzymes in the eugenol synthesis pathway: Coniferyl alcohol acetyltransferase, eugenol synthase and isoeugenol synthase68. Coniferyl alcohol acetyltransferase (EC: 2.3.1.224) associated with two proteins (Tfoe_c37673_g1_i1_m27189, Tfoe_c102741_g1_i1_m76275), however, Tfoe_c37673_g1_i1_m27189 was selected to be the true candidate gene (marked in Supplementary Data 2) since it co-clustered with the M. truncatula protein and was also supported by the presence of the FIR. Similarly, the candidate proteins from fenugreek for eugenol synthase (EC: 1.1.1.318) and the isoeugenol synthase (EC: 1.1.1.319) co-clustered with clades belonging to known enzymes (highlighted in Supplementary Data 2). The average TPM values across tissues for each of the enzymes involved in synthesis of eugenol are represented in Supplementary Table 3. The transcripts corresponding to both the enzymes were found to be most abundant in root compared to the other tissues.
In the linalool synthesis pathway we studied three enzymes: Farnesyl-pyrophosphate synthase, (3S)-linalool synthase and (3R)-linalool synthase69,70. (3S)-linalool and (3R)-linalool are two enantiomeric forms of naturally occurring monoterpenoids. The candidate protein from fenugreek for farnesyl-pyrophosphate synthase (EC: 2.5.1.1) (Tfoe_c5299_g1_i1_m1972) co-clustered with another legume Lupinus albus (Fabaceae family). The (3S)-linalool synthase (EC: 4.2.3.25) (Tfoe_c11093_g1_i1_m4195, Tfoe_c11093_g1_i2_m4196) and (3R)-linalool synthase (EC: 4.2.3.26) (Tfoe_c44084_g4_i2_m55446) co-clustered with other sequence homologues (as highlighted in the Supplementary Data 2). The average TPM values for each of the enzymes involved in synthesis of linalool are represented in Supplementary Table 3. While the first enzyme-farnesyl-pyrophosphate synthase, is expressed highest in root compared to the other tissues, (3S)-linalool synthase and (3R)-linalool synthase is expressed highly in leaves.
Flavonoids, an important class of secondary metabolites with a polyphenolic structure, widely found in fruits, vegetables and certain beverages. This class of compounds has been linked to a variety of nutraceutical, pharmaceutical, and medicinal properties, including antioxidative, anti-inflammatory, antimutagenic, and anticarcinogenic properties71. We studied biosynthetic pathways of four flavonoids (quercetin, vitexin, isovitexin and rutin) known to be present in fenugreek63. These metabolites are synthesised from coumaroyl-CoA through different enzymatic reactions. The flavonoids share a common flavonoid backbone structure with differences in their hydroxylation patterns.
Quercetin biosynthesis pathway involves six enzymes [4-Coumarate-CoA ligase (EC: 6.2.1.12), chalcone synthase (EC: 2.3.1.74), chalcone flavone isomerase (EC: 5.5.1.6), flavanone 3-hydroxylase (EC: 1.14.11.9)/flavonol synthase (EC: 1.14.20.6), tricin synthase (EC: 2.1.1.175) and flavonoid 3′-monooxygenase (EC: 1.14.14.82)] that convert coumaroyl-CoA to quercetin through several intermediates17. For each of the six enzymes involved in the pathway, multiple transcript hits were identified in fenugreek as mentioned in Supplementary Table 3. The sequence hits for these enzymes were validated as true hits, based on co-clustering with well annotated sequence of same sub-family observed in the phylogenetic tree (as described in Methods “Identification and analysis of genes involved in synthesis of secondary metabolites”). It was observed that the selected enzymes are highly abundant in roots of the plant followed by the stem (Supplementary Table 3, Supplementary Data 2).
We further studied biosynthetic pathways for vitexin and isovitexin, which are two enantiomeric metabolites. The enzymes involved their pathway are identical72, and they produce these metabolites in a mixture. A three-step reaction is known to be involved in the conversion of naringenin to vitexin/isovitexin. The first two reactions are enzymatic while the third step (hydration reaction)73, although, has been hypothesized to be an enzymatic one, there is no definitive enzyme that has been isolated for it. For the first two enzymes [naringenin 2-hydroxylase (EC: 1.14.14.162) and C-glucosyltransferase (EC: 2.4.1.360)] known in the reaction, four and 13 hits were obtained from the fenugreek transcriptome Supplementary Table 3. These hits were validated using phylogenetic analysis of the enzymes from the same sub-family (Supplementary Data 2). The highest abundance of the naringenin 2-hydroxylase in the pathway was observed in the stem, while for C-glucosyltransferase, although there were 13 hits identified fenugreek, most of them were highly expressed in leaf.
The third flavonoid, rutin, is synthesised in plants from taxol. The pathway involves two enzymatic reactions. While vitexin/isovitexin are obtained by C-glycosylation, rutin is obtained by O-glycosylation of the substrates. The enzymes involved in the pathway are flavonol-3-o-glucosyltransferase (EC: 2.4.1.91) and flavonol-3-o-glucoside rhamnosyltransferase (EC: 2.4.1.15)74. There were 11 hits identified for flavonol-3-o-glucosyltransferase, and most of them had higher expression value in leaf and root tissue. In case of flavonol-3-o-glucoside rhamnosyltransferase a single hit (Tfoe_c121323_g1_i1_m79650) was identified with a high abundance in leaf (Supplementary Table 3, Supplementary Data 2).
Fenugreek is a traditionally well-known plant for its medicinal properties and the presence of complex mixtures of key phytochemicals renders this plant such benefits4,63. In the literature, details of these chemical constituents in fenugreek tissue extracts have mainly been demonstrated by bio-analytical techniques. We quantified, presence of important metabolites- diosgenin, trigonelline and 4-hydroxyisoleucine (4-HIL) in different tissues (leaf, stem, root and seed) of fenugreek, by HPLC and mass spectrometry. Although high-quality RNA could not be purified from the seed tissue, we used it nonetheless for metabolite quantification. Diosgenin and trigonelline were quantified using HPLC–PDA and 4-HIL was quantified using LC–MS analysis from extracts of different tissues as described in Supplementary Methods Sect. 1.2. Supplementary Fig. 6A,B show the standard curve of diosgenin and trigonelline respectively. The concentration of each compound was quantitatively determined by comparison of the peak area of the standard with that of the samples. Supplementary Fig. 7A–D,E–H refers to the HPLC peaks which showed the presence of trigonelline and diosgenin, and Supplementary Fig. 7L–N refers to LC profile from LC–MS data of 4-HIL. As seen in Supplementary Table 4, seed gave the highest concentration of diosgenin i.e. 1.56 µg/mg of the tissue followed by the root (1.036 µg/mg). Leaf tissue showed the highest concentration of trigonelline i.e. 4.715 µg/mg of the tissue followed by the seed (2.698 µg/mg). 4-HIL was detected at a very high concentration of 84.77 µg/mg in the dry seed and 77.48 µg/mg from the leaf. The concentration of 4-HIL was observed to significantly decrease with water-soaked seeds post 12 h and 24 h. The concentrations of these medicinally important metabolites vary across several cultivars of fenugreek75,76.
We have surveyed here the landscape of secondary metabolites in fenugreek, which are medicinally important, through the lens of transcriptome data. Apart from being a food crop of semi-arid nature, fenugreek is not only a food source but has been used in traditional medicines in several parts of the world. The medicinal properties of this drought tolerant leguminous plant can largely be attributed to an array of secondary metabolites produced in it63. In this study we have identified candidate genes involved in the biosynthesis of such metabolites belonging to four major classes, alkaloids, saponins, volatiles and flavonoids. The process adopted for identification of such genes in the transcriptome has been previously implemented in Ocimum tenuiflorum, Moringa oleifera17,38,77. It involved creation and curation of a knowledgebase for candidate gene products from other (closely) related plants and applying it through a rigorous sequence search and validation process to ascertain the transcripts in the fenugreek transcriptome assembly. The raw data from leaf, stem and root tissue samples was processed and assembled into a set of transcripts and their level of expression was compared across these tissues in measure of average TPM values. In particular, several drought responsive genes such as metallothionein, dehydrin B and Nmr-A like family of genes as described in the results “Transcript abundance estimation” were found to be among most abundant transcripts, especially in the root tissue sample.
While identifying candidate genes involved in biosynthesis of secondary metabolites, it is vital to select appropriate start points for the rigorous sequence search and subsequent validation. We rely on multiple sources such as PlantCyc, MetaCyc, KEGG pathway and metabolites cited in literature11,39,40,41,42 to create such a knowledge base of start points. It further helps in identification of FIRs which resolves selection of true candidate genes along with the co-clustering of proteins from fenugreek transcripts with the selected start points. For instance, we considered the trigonelline synthesising enzyme NNMT (“Alkaloids”, Fig. 4), which methylates nicotinate to methyl nicotinate or trigonelline using S-adenosyl-l-methionine (SAM) as the methyl donor. We identified 12 hits from the transcriptome through the sequence search, however, we make use of the knowledge-driven filtering to ascertain the correct annotation. We identified a single hit (Tfoe_c19814_g2_i1_m7872) in which we could find all the functionally important residues viz. the catalytic residue Thr 264, substrate-binding motif Asn 21, Tyr 120, His 124 and a SAM binding motif (Supplementary Fig. 4A–D). Subsequently this sequence co-clustered with the NNMT sequence (XP_013464471.1_MEDTR) from Medicago truncatula, which is closely related to fenugreek from the Trifolieae tribe (Supplementary Fig. 4E). These multiple checkpoints helped in selecting the correct fenugreek hit for rendering the annotation. The TPM values for this transcript were compared across all the tissues and found to be higher in leaves. This trend was followed in the qRT-PCR based expression (Supplementary Fig. 5A) as well the metabolite quantification (“Quantification of metabolites in different tissue samples”), suggesting that the fenugreek leaves are a rich source of this anti-diabetic compound, trigonelline. This knowledge-driven protocol was applied in selecting the correct start points and ascertaining the true fenugreek hits for all the other alkaloids, saponins, flavonoids and volatiles (described in the Results “Identification and analysis of genes involved in synthesis of secondary metabolites”). The candidate genes identified for the production of volatiles- eugenol, isoeugenol and linalool were found to be highly abundant in roots except for the enzyme (3S)-linalool synthase and (3R)-linalool synthase which were abundant in leaves compared to other tissues. The transcripts for the enzymes involved in the biosynthesis of the flavonoids-vitexin, isovitexin and rutin were found to be highly abundant in leaves, whereas, for quercetin biosynthesis the transcripts were abundant in the root and stem.
In the case of the diosgenin pathway, which is well studied in fenugreek78, there are two enzymes which were explored in the current study. Sterol-3β-glucosyltransferase (the penultimate enzyme) is associated with four hits in fenugreek. However, based on identification of PSBD and PSPG motifs and co-clustering with sequences from closely related leguminous plants, ensured identification of two true candidate genes (Tfoe_c26467_g1_i1_m12030, Tfoe_c33087_g1_i3_m19456) (“Saponins”). The first hit was found to be abundant across all the three tissues- leaves, stem and root while the second hit was found to be high in leaves. This was observed through the average TPM values across the tissue. Additionally, the qRT-PCR quantification also corroborated with this trend. The last enzyme in diosgenin biosynthesis, β-glucosidase, despite presence of (I/V)TENG and TFNEP motifs in all of the associated 22 hits, only a single hit was selected as the candidate gene based on co-clustering with M. truncatula protein sequence. This enzyme was found to be abundant in root and stem from the TPM values, which was corroborated by the qRT-PCR studies (“Saponins”). Diosgenin quantification from different tissues also indicated predominantly high concentrations in the root (“Quantification of metabolites in different tissue samples”). There are several transcriptome studies available in NCBI for fenugreek albeit with different cultivars and under different study conditions (BioProject IDs: PRJNA544308, PRJNA508420, PRJNA383660). We confirmed the presence of both the diosgenin pathway enzymes in their corresponding SRA datasets. We observed more than 98% sequence identity and close to 90% query coverage for the candidate genes from our study. This suggested that the candidate genes identified in our study agree well with the publicly available data for fenugreek.
Several fenugreek cultivars are widely cultivated across India. Traditionally, the AFG-1 variety, described in this study, has been one of the cultivars commonly grown across India, with bold seeds and good yield79. Many of the fenugreek cultivars grown across India face the powdery mildew disease, including the AFG-1 variety80. Ongoing work in UHS, Bagalkot, India (unpublished), on Indian fenugreek cultivars for identification of biotic stress resistant genotype, has observed that among selected cultivars, Kasuri methi showed high immune response towards the biotic stress arising from the powdery mildew disease. This could be further utilized to understand disease resistance mechanisms in fenugreek cultivars. The current transcriptome data will assist in establishing the genetic determinants for the nutritional and medicinal properties among the cultivars. It will also help us to estimate the important secondary metabolites and develop varieties with high nutraceutical value. This will further pave the way to develop and breed biotic and abiotic tolerant cultivars of fenugreek.
We present here the transcriptome of the Ajmer Fenugreek-1 (AFG-1) variety of T. foenum-graecum. We identified candidate genes of enzymes involved in the biosynthesis of secondary metabolites such as alkaloids, saponins, flavonoids and volatile compounds. Some of these secondary metabolites are of immense value in the pharmaceutical and nutraceutical industries. Among these, diosgenin and trigonelline are known to be specific to fenugreek and are antidiabetic, antioxidant, anti-inflammatory, antiobesity along with several other medicinal values. The abundance of transcripts encoding these candidate genes was assessed and validated through qRT-PCR for tissue-level expression. Additionally, these metabolites were quantified across the tissues. Identification of genes involved in secondary metabolites biosynthesis in fenugreek will aid in bioengineering efforts for large scale production. Further, this transcriptome can be utilised for cross comparison of several cultivars and through molecular breeding programs lead to crop improvement.
NCBI BioProject: PRJNA734905, BioSample: (leaf replicates: SAMN19548760, SAMN19548761; stem replicates: SAMN19548765, SAMN19548764; root replicates: SAMN19548762, SAMN19548763), SRA accession: (leaf replicates: SRR14721911, SRR14721912; stem replicates: SRR14721915, SRR14721916; root replicates: SRR14721913, SRR14721914). The Supplementary Data is available at the following link: http://caps.ncbs.res.in/download/tfoe_data/. The details are indexed in Supplementary Text file (“Materials and methods”).
Trigonella foenum-graecum
Ajmer Fenugreek-1
Open reading frame
Transcripts per million
Functionally important residues
Quantitative reverse transcriptase polymerase chain reaction
4-Hydroxyleucine
Gene ontology
Nicotinate N-methyl transferase
Betaine aldehyde dehydrogenase
Choline monooxygenase
Putative sterol binding domain
Plant secondary product glycosyltransferase
High performance liquid chromatography–photo diode array
Liquid chromatography–mass spectrometry
Martin, E., Akan, H., Ekici, M. & Ayta, Z. Karyotype analyses of ten sections of Trigonella (Fabaceae). Comp. Cytogenet. 5, 105–121 (2011).
CAS  Article  Google Scholar 
Ranjbar, M. & Zahra, H. Chromosome numbers and biogeography of the genus Trigonella (Fabaceae). Caryologia https://doi.org/10.1080/00087114.2016.1169090 (2016).
Article  Google Scholar 
Ahmad, A., Alghamdi, S. S., Mahmood, K. & Afzal, M. Fenugreek a multipurpose crop: Potentialities and improvements. Saudi J. Biol. Sci. 23, 300–310 (2016).
Article  Google Scholar 
Nagulapalli Venkata, K. C., Swaroop, A., Bagchi, D. & Bishayee, A. A small plant with big benefits: Fenugreek (Trigonella foenum-graecum Linn.) for disease prevention and health promotion. Mol. Nutr. Food Res. 61, 1600950 (2017).
Article  Google Scholar 
Srinivasan, K. Fenugreek (Trigonella foenum-graecum): A review of health beneficial physiological effects. Food Rev. Int. https://doi.org/10.1080/87559120600586315 (2006).
Article  Google Scholar 
Yadav, U. C. S. & Baquer, N. Z. Pharmacological effects of Trigonella foenum-graecum L. in health and disease. Pharm. Biol. 52, 243–254 (2014).
Article  Google Scholar 
Syed, Q. A. et al. Nutritional and therapeutic properties of fenugreek (Trigonella foenum-graecum): A review. Int. J. Food Prop. https://doi.org/10.1080/10942912.2020.1825482 (2020).
Article  Google Scholar 
Vaughn, C. J. Drugs and lactation database: LactMed. J. Electron. Resour. Med. Libr. https://doi.org/10.1080/15424065.2012.735134 (2012).
Article  Google Scholar 
Yao, D. et al. Advances on application of fenugreek seeds as functional foods: Pharmacology, clinical application, products, patents and market. Crit. Rev. Food Sci. Nutr. 60, 2342–2352 (2020).
CAS  Article  Google Scholar 
Jesus, M., Martins, A. P. J., Gallardo, E. & Silvestre, S. Diosgenin: Recent highlights on pharmacology and analytical methodology. J. Anal. Methods Chem. 2016, 4156293 (2016).
Article  Google Scholar 
Zhou, J., Chan, L. & Zhou, S. Trigonelline: A plant alkaloid with therapeutic potential for diabetes and central nervous system disease. Curr. Med. Chem. 19, 3523–3531 (2012).
CAS  Article  Google Scholar 
Broca, C. et al. 4-Hydroxyisoleucine: Experimental evidence of its insulinotropic and antidiabetic properties. Am. J. Physiol.-Endocrinol. Metab. 277, E617–E623 (1999).
CAS  Article  Google Scholar 
Avalos-Soriano, A., de la Cruz-Cordero, R., Rosado, J. L. & Garcia-Gasca, T. 4-Hydroxyisoleucine from Fenugreek (Trigonella foenum-graecum): Effects on insulin resistance associated with obesity. Molecules (Basel, Switzerland) 21, 1596 (2016).
Article  Google Scholar 
Tewari, D. et al. Fenugreek (Trigonella foenum-graecum L.) seeds dietary supplementation regulates liver antioxidant defense systems in aging mice. Nutrients 12, 2552 (2020).
CAS  Article  Google Scholar 
Choudhary, S. et al. Transcriptome profiling of coriander: A dual purpose crop unravels stem gall resistance genes. J. Genet. https://doi.org/10.1007/s12041-019-1064-7 (2019).
Article  PubMed  Google Scholar 
Choudhary, S., Naika, M. B. N. & Meena, R. D. Development and characterization of genic SSR-FDM for stem gall disease resistance in coriander (Coriandrum sativum L.) and its cross species transferability. Mol. Biol. Rep. https://doi.org/10.1007/s11033-021-06396-9 (2021).
Article  PubMed  Google Scholar 
Pasha, S. N. et al. The transcriptome enables the identification of candidate genes behind medicinal value of Drumstick tree (Moringa oleifera). Genomics https://doi.org/10.1016/j.ygeno.2019.04.014 (2020).
Article  PubMed  Google Scholar 
Bolger, A. M., Lohse, M. & Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics (Oxford, England) 30, 2114–2120 (2014).
CAS  Article  Google Scholar 
Haas, B. J. et al. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat. Protoc. 8, 1494–1512 (2013).
CAS  Article  Google Scholar 
Nishimura, O., Hara, Y. & Kuraku, S. gVolante for standardizing completeness assessment of genome and transcriptome assemblies. Bioinformatics 33, 3635–3637 (2017).
CAS  Article  Google Scholar 
Altschul, S. F. BLAST algorithm. Encycl. Life Sci. https://doi.org/10.1038/npg.els.0005253 (2005).
Article  Google Scholar 
UniProt: The universal protein knowledgebase in 2021. Nucleic Acids Res. 49, D480–D489 (2021).
Finn, R. D., Clements, J. & Eddy, S. R. HMMER web server: Interactive sequence similarity searching. Nucleic Acids Res. 39, W29–W37 (2011).
CAS  Article  Google Scholar 
Sonnhammer, E. L. L., Eddy, S. R. & Durbin, R. Pfam: A comprehensive database of protein domain families based on seed alignments. Proteins Struct. Funct. Genet. 28, 405–420 (1997).
CAS  Article  Google Scholar 
Petersen, T. N., Brunak, S., von Heijne, G. & Nielsen, H. SignalP 4.0: Discriminating signal peptides from transmembrane regions. Nat. Methods 8, 785–786 (2011).
CAS  Article  Google Scholar 
Krogh, A., Larsson, B., von Heijne, G. & Sonnhammer, E. L. L. Predicting transmembrane protein topology with a hidden Markov model: Application to complete genomes. J. Mol. Biol. https://doi.org/10.1006/jmbi.2000.4315 (2001).
Article  PubMed  Google Scholar 
Lagesen, K. et al. RNAmmer: Consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 35, 3100–3108 (2007).
ADS  CAS  Article  Google Scholar 
Ashburner, M. et al. Gene ontology: Tool for the unification of biology. Nat. Genet. https://doi.org/10.1038/75556 (2000).
Article  PubMed  PubMed Central  Google Scholar 
Huerta-Cepas, J. et al. eggNOG 5.0: A hierarchical, functionally and phylogenetically annotated orthology resource based on 5090 organisms and 2502 viruses. Nucleic Acids Res. 47, D309–D314 (2019).
CAS  Article  Google Scholar 
Kanehisa, M. & Goto, S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 28, 27–30 (2000).
CAS  Article  Google Scholar 
Li, L., Stoeckert, C. J. & Roos, D. S. OrthoMCL: Identification of ortholog groups for eukaryotic genomes. Genome Res. 13, 2178–2189 (2003).
CAS  Article  Google Scholar 
Goodstein, D. M. et al. Phytozome: A comparative platform for green plant genomics. Nucleic Acids Res. 40, D1178–D1186 (2012).
CAS  Article  Google Scholar 
Lechner, M. et al. Proteinortho: Detection of (Co-)orthologs in large-scale analysis. BMC Bioinform. 12, 124 (2011).
Article  Google Scholar 
Katoh, K. & Standley, D. M. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol. Biol. Evol. 30, 772–780 (2013).
CAS  Article  Google Scholar 
Li, B. & Dewey, C. N. RSEM: Accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinform. 12, 323 (2011).
CAS  Article  Google Scholar 
Roberts, A., Trapnell, C., Donaghey, J., Rinn, J. L. & Pachter, L. Improving RNA-Seq expression estimates by correcting for fragment bias. Genome Biol. https://doi.org/10.1186/gb-2011-12-3-r22 (2011).
Article  PubMed  PubMed Central  Google Scholar 
Langmead, B. & Salzberg, S. L. Fast gapped-read alignment with Bowtie 2. Nat. Methods https://doi.org/10.1038/nmeth.1923 (2012).
Article  PubMed  PubMed Central  Google Scholar 
Joshi, A. G. et al. A knowledge-driven protocol for prediction of proteins of interest with an emphasis on biosynthetic pathways. MethodsX https://doi.org/10.1016/j.mex.2020.101053 (2020).
Article  PubMed  PubMed Central  Google Scholar 
Schläpfer, P. et al. Genome-wide prediction of metabolic enzymes, pathways, and gene clusters in plants. Plant Physiol. 173, 2041–2059 (2017).
Article  Google Scholar 
Caspi, R. et al. The MetaCyc database of metabolic pathways and enzymes and the BioCyc collection of Pathway/Genome Databases. Nucleic Acids Res. 42, D459–D471 (2014).
CAS  Article  Google Scholar 
Mizuno, K. et al. Conversion of nicotinic acid to trigonelline is catalyzed by N-methyltransferase belonged to motif B’ methyltransferase family in Coffea arabica. Biochem. Biophys. Res. Commun. 452, 1060–1066 (2014).
CAS  Article  Google Scholar 
Li, W. et al. A novel N-methyltransferase in arabidopsis appears to feed a conserved pathway for nicotinate detoxification among land plants and is associated with lignin biosynthesis. Plant Physiol. 174, 1492–1504 (2017).
CAS  Article  Google Scholar 
Sievers, F. et al. Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mol. Syst. Biol. 7, 539 (2011).
Article  Google Scholar 
Kumar, S., Stecher, G., Li, M., Knyaz, C. & Tamura, K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 35, 1547–1549 (2018).
CAS  Article  Google Scholar 
Rambaut, A. FigTree v1.4.2, A Graphical Viewer of Phylogenetic Trees. http://tree.bio.ed.ac.uk/software/figtree/ (2014).
Letunic, I. & Bork, P. Interactive Tree Of Life (iTOL): An online tool for phylogenetic tree display and annotation. Bioinformatics (Oxford, England) 23, 127–128 (2007).
CAS  Article  Google Scholar 
Joshi, R. & Karan, R. Physiological, biochemical and molecular mechanisms of drought tolerance in plants. Mol. Approach. Plant Abiotic Stress https://doi.org/10.1201/b15538-16 (2013).
Article  Google Scholar 
Yang, Z., Wu, Y., Li, Y., Ling, H. Q. & Chu, C. OsMT1a, a type 1 metallothionein, plays the pivotal role in zinc homeostasis and drought tolerance in rice. Plant Mol. Biol. https://doi.org/10.1007/s11103-009-9466-1 (2009).
Article  PubMed  Google Scholar 
Lopez, C. G., Banowetz, G. M., Peterson, C. J. & Kronstad, W. E. Dehydrin expression and drought tolerance in seven wheat cultivars. Crop Sci. https://doi.org/10.2135/cropsci2003.0577 (2003).
Article  Google Scholar 
Hassan, N. M., El-Bastawisy, Z. M., El-Sayed, A. K., Ebeed, H. T. & Nemat Alla, M. M. Roles of dehydrin genes in wheat tolerance to drought stress. J. Adv. Res. https://doi.org/10.1016/j.jare.2013.11.004 (2015).
Article  PubMed  PubMed Central  Google Scholar 
Waller, G. R. & Nowacki, E. K. The Role of Alkaloids in Plants BT—Alkaloid Biology and Metabolism in Plants (eds. Waller, G. R. & Nowacki, E. K.) 143–181 (Springer US, 1978). https://doi.org/10.1007/978-1-4684-0772-3_5.
Zhou, J., Chan, L. & Zhou, S. Trigonelline: A plant alkaloid with therapeutic potential for diabetes and central nervous system disease. Curr. Med. Chem. https://doi.org/10.2174/092986712801323171 (2012).
Article  PubMed  Google Scholar 
Zhao, G. et al. Betaine in inflammation: Mechanistic aspects and applications. Front. Immunol. https://doi.org/10.3389/fimmu.2018.01070 (2018).
Article  PubMed  PubMed Central  Google Scholar 
Lee, J. E. et al. Choline and betaine intake and the risk of colorectal cancer in men. Cancer Epidemiol. Biomark. Prev. https://doi.org/10.1158/1055-9965.EPI-09-1295 (2010).
Article  Google Scholar 
Sun, J., Wen, S., Zhou, J. & Ding, S. Association between malnutrition and hyperhomocysteine in Alzheimer’s disease patients and diet intervention of betaine. J. Clin. Lab. Anal. https://doi.org/10.1002/jcla.22090 (2017).
Article  PubMed  PubMed Central  Google Scholar 
Hassanein, R. A. Effect of heat shock on some biochemical and molecular criteria of fenugreek (Trigonella foenum-graceum L.). J. Med. Plants Res. https://doi.org/10.5897/jmpr11.1624 (2012).
Article  Google Scholar 
Rathinasabapathi, B. et al. Choline monooxygenase, an unusual iron-sulfur enzyme catalyzing the first step of glycine betaine synthesis in plants: Prosthetic group characterization and cDNA cloning. Proc. Natl. Acad. Sci. U.S.A. https://doi.org/10.1073/pnas.94.7.3454 (1997).
Article  PubMed  PubMed Central  Google Scholar 
Wink, M. Evolution of secondary metabolites in legumes (Fabaceae). S. Afr. J. Bot. https://doi.org/10.1016/j.sajb.2013.06.006 (2013).
Article  Google Scholar 
Giridhar, K., Surya Kumari, S., Rajani, A., Sarada, C. & Naram Naidu, L. Identification of potential genotypes of fenugreek in rainfed vertisols for yield and diosgenin content. Indian J. Agric. Res. https://doi.org/10.18805/ijare.v0iOF.8603 (2016).
Article  Google Scholar 
Ciura, J., Szeliga, M., Grzesik, M. & Tyrka, M. Next-generation sequencing of representational difference analysis products for identification of genes involved in diosgenin biosynthesis in fenugreek (Trigonella foenum-graecum). Planta https://doi.org/10.1007/s00425-017-2657-0 (2017).
Article  PubMed  PubMed Central  Google Scholar 
Ramirez-Estrada, K. et al. Tomato UDP-glucose sterol glycosyltransferases: A family of developmental and stress regulated genes that encode cytosolic and membrane-associated forms of the enzyme. Front. Plant Sci. https://doi.org/10.3389/fpls.2017.00984 (2017).
Article  PubMed  PubMed Central  Google Scholar 
Chuenchor, W. et al. Structural insights into rice BGlu1 β-glucosidase oligosaccharide hydrolysis and transglycosylation. J. Mol. Biol. https://doi.org/10.1016/j.jmb.2008.01.076 (2008).
Article  PubMed  Google Scholar 
Wani, S. A. & Kumar, P. Fenugreek: A review on its nutraceutical properties and utilization in various food products. J. Saudi Soc. Agric. Sci. 17, 97–106 (2018).
Google Scholar 
Al-Daghri, N. M. et al. Fenugreek extract as an inducer of cellular death via autophagy in human T lymphoma Jurkat cells. BMC Complement. Altern. Med. https://doi.org/10.1186/1472-6882-12-202 (2012).
Article  PubMed  PubMed Central  Google Scholar 
Yogalakshmi, B., Viswanathan, P. & Anuradha, C. V. Investigation of antioxidant, anti-inflammatory and DNA-protective properties of eugenol in thioacetamide-induced liver injury in rats. Toxicology 268, 204–212 (2010).
CAS  Article  Google Scholar 
Rauscher, F. M., Sanders, R. A. & Watkins, J. B. Effects of isoeugenol on oxidative stress pathways in normal and streptozotocin-induced diabetic rats. J. Biochem. Mol. Toxicol. https://doi.org/10.1002/jbt.13 (2001).
Article  PubMed  Google Scholar 
Pereira, I., Severino, P., Santos, A. C., Silva, A. M. & Souto, E. B. Linalool bioactive properties and potential applicability in drug delivery systems. Colloids Surf. B https://doi.org/10.1016/j.colsurfb.2018.08.001 (2018).
Article  Google Scholar 
Koeduka, T. et al. Eugenol and isoeugenol, characteristic aromatic constituents of spices, are biosynthesized via reduction of a coniferyl alcohol ester. Proc. Natl. Acad. Sci. U.S.A. https://doi.org/10.1073/pnas.0603732103 (2006).
Article  PubMed  PubMed Central  Google Scholar 
Raguso, R. A. & Pichersky, E. A day in the life of a linalool molecule: Chemical communication in a plant-pollinator system. Part 1: Linalool biosynthesis in flowering plants. Plant Species Biol. https://doi.org/10.1046/j.1442-1984.1999.00014.x (1999).
Article  Google Scholar 
Carrau, F. M. et al. De novo synthesis of monoterpenes by Saccharomyces cerevisiae wine yeasts. FEMS Microbiol. Lett. https://doi.org/10.1016/j.femsle.2004.11.050 (2005).
Article  PubMed  Google Scholar 
Panche, A. N., Diwan, A. D. & Chandra, S. R. Flavonoids: An overview. J. Nutr. Sci. 5, e47–e47 (2016).
CAS  Article  Google Scholar 
Du, Y., Chu, H., Chu, I. K. & Lo, C. CYP93G2 is a flavanone 2-hydroxylase required for C-glycosylflavone biosynthesis in rice. Plant Physiol. 154, 324–333 (2010).
CAS  Article  Google Scholar 
Kerscher, F. & Franz, G. Biosynthesis of vitexin and isovitexin: enzymatic synthesis of the C-glucosylflavones vitexin and isovitexin with an enzyme preparation from Fagopyrum esculentum M. Seedlings. Zeitschrift für Naturforschung C 42, 519–524 (1987).
CAS  Article  Google Scholar 
Lucci, N. & Mazzafera, P. Rutin synthase in fava d’anta: Purification and influence of stressors. Can. J. Plant Sci. 89, 895–902 (2009).
CAS  Article  Google Scholar 
Paramesha, M., Priyanka, N., Crassina, K. & Shetty, N. P. Evaluation of diosgenin content from eleven different Indian varieties of fenugreek and fenugreek leaf powder fortified bread. J. Food Sci. Technol. https://doi.org/10.1007/s13197-021-04967-z (2021).
Article  PubMed  PubMed Central  Google Scholar 
Hassanzadeh, E., Chaichi, M. R., Mazaheri, D., Rezazadeh, S. & Badi, H. A. N. Physical and chemical variabilities among domestic Iranian Fenugreek (Trigonella foenum-graceum) seeds. Asian J. Plant Sci. https://doi.org/10.3923/ajps.2011.323.330 (2011).
Article  Google Scholar 
Upadhyay, A. K. et al. Genome sequencing of herb Tulsi (Ocimum tenuiflorum) unravels key genes behind its strong medicinal properties. BMC Plant Biol. https://doi.org/10.1186/s12870-015-0562-x (2015).
Article  PubMed  PubMed Central  Google Scholar 
Ciura, J., Szeliga, M., Grzesik, M. & Tyrka, M. Next-generation sequencing of representational difference analysis products for identification of genes involved in diosgenin biosynthesis in fenugreek (Trigonella foenum-graecum). Planta 245, 977–991 (2017).
CAS  Article  Google Scholar 
ICAR-NRCSS. List of varieties/technologies developed at ICAR-NRCSS.
Rathi, A. S., Gupta, P. P. & Jhorar, B. S. Yield losses due to powdery mildew disease in fenugreek (Trigonella foenum-graecum L.). J. Spices Aromat. Crops 11, 143–145 (2002).
Google Scholar 
Download references
We thank the reviewers and the editorial office for all the comments and points to improve the manuscript. We acknowledge the initial discussions regarding the transcriptome analysis pipeline with A. Gandhimathi, Dhingra Surbhi, Mathew Oommen K., Malhotra Sony, Manoharan Malini, Mutt Eshita, Shingate Prashant N., Sukhwal Anshul, Sunitha Margaret S., Upadhyay Atul K. and Upadhyayula Raghavender S. We also thank Dr. Padma Ramakrishnan from Centre for cellular and molecular platforms, Bengaluru for use of mass spectrometry and HPLC facilities. We thank College of Horticulture, GKVK campus, Bengaluru for plant seed materials. We also thank Department of Crop Physiology, UAS, Bengaluru for use of experimental facilities.
This work was supported by the Department of Biotechnology, Government of India grant (BT/PR10550/BID/7/479/2013), the JC Bose fellowship, Science and Engineering Research Board, Department of Science and Technology, Government of India (SB/S2/JC-071/2015), Institute of Bioinformatics and Applied Biotechnology for the funding through her Mazumdar-Shaw Chair in Computational Biology (IBAB/MSCB/182/2022) and Bioinformatics Centre Grant funded by Department of Biotechnology, India (BT/PR40187/BTIS/137/9/2021) to RS. We also thank NCBS (TIFR) for infrastructural and financial support.
National Centre for Biological Sciences (TIFR), GKVK Campus, Bengaluru, India
Mahantesha B. N. Naika, Nitish Sathyanarayanan, Radha Sivarajan Sajeevan, Teerna Bhattacharyya, Pritha Ghosh, Meenakshi S. Iyer, Mahita Jarjapu, Adwait G. Joshi, K. Harini, K. Mohamed Shafi, Neha Kalmankar, Snehal D. Karpe, Bhavika Mam, Shaik Naseer Pasha & Ramanathan Sowdhamini
Department of Biotechnology and Crop Improvement, K.R.C. College of Horticulture, University of Horticultural Sciences-Bagalkot, Arabhavi, Karnataka, 591218, India
Mahantesha B. N. Naika
The University of Trans-Disciplinary Health Sciences & Technology (TDU), Yelahanka, Bengaluru, Karnataka, 560064, India
Nitish Sathyanarayanan, K. Mohamed Shafi, Neha Kalmankar & Bhavika Mam
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
R.S.: Conceived and designed experiments, provided resources and tools and critically reviewed the manuscript. M.B.N.N., S.N., R.S.S., I.M.S., K.H., K.N., S.N.P.: Data generation. M.B.N.N., S.N., R.S.S., B.T., G.P., I.M.S., J.M., J.A.G., K.H., K.M.S., K.N., K.S.D., M.B., S.N.P.: Data analysis and presentation. M.B.N.N., S.N., R.S.S., B.T., G.P., I.M.S., J.M., J.A.G., K.H., K.M.S., K.N., K.S.D., M.B., S.N.P. and R.S.: Writing of manuscript.
Correspondence to Ramanathan Sowdhamini.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Reprints and Permissions
Naika, M.B.N., Sathyanarayanan, N., Sajeevan, R.S. et al. Exploring the medicinally important secondary metabolites landscape through the lens of transcriptome data in fenugreek (Trigonella foenum graecum L.). Sci Rep 12, 13534 (2022). https://doi.org/10.1038/s41598-022-17779-8
Download citation
Received: 11 March 2022
Accepted: 30 July 2022
Published: 08 August 2022
DOI: https://doi.org/10.1038/s41598-022-17779-8
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Advertisement
Advanced search
Scientific Reports (Sci Rep) ISSN 2045-2322 (online)
© 2022 Springer Nature Limited
Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

source

How do I find SEO services
Dominate search engine result pages, crush your competition, get more sales, and be the #1 provider in your service area. Work with Top SEO Agency in Dublin that also provides quality Digital Marketing Optimisation services. 

Our Search Sngine Optimisation Consultants will improve your site performance in all major search engines by implementing the Best SEO Solutions, Strategies and Techniques. Let us help you increase organic traffic, get more leads, more customers, and grow your revenue with a customised affordable SEO package.

Contact Us

Globe Boss Digital Marketing Agency

Hampton Square
Dublin - Ireland
Tel: (+353)1 868 2345