Identification of Key Metabolic Pathways and Biomarkers Underlying Flowering Time of Guar (Cyamopsis tetragonoloba (L.) Taub.) via Integrated Transcriptome-Metabolome Analysis.

Guar (Cyamopsis tetragonoloba (L.) Taub.) is an annual legume crop native to India and Pakistan. Seeds of the plant serve as a source of galactomannan polysaccharide (guar gum) used in the food industry as a stabilizer (E412) and as a gelling agent in oil and gas fracturing fluids. There were several attempts to introduce this crop to countries of more northern latitudes. However, guar is a plant of a short photoperiod, therefore, its introduction, for example, to Russia is complicated by a long day length during the growing season. Breeding of new guar varieties insensitive to photoperiod slowed down due to the lack of information on functional molecular markers, which, in turn, requires information on guar genome. Modern breeding strategies, e.g., genomic predictions, benefit from integration of multi-omics approaches such as transcriptome, proteome and metabolome assays. Here we present an attempt to use transcriptome-metabolome integration to understand the genetic determination of flowering time variation among guar plants that differ in their photoperiod sensitivity. This study was performed on nine early- and six delayed-flowering guar varieties with the goal to find a connection between 63 metabolites and 1,067 differentially expressed transcripts using Shiny GAM approach. For the key biomarker of flowering in guar myo-inositol we also evaluated the KEGG biochemical pathway maps available for Arabidopsis thaliana. We found that the phosphatidylinositol signaling pathway is initiated in guar plants that are ready for flowering through the activation of the phospholipase C (PLC) gene, resulting in an exponential increase in the amount of myo-inositol in its free form observed on GC-MS chromatograms. The signaling pathway is performed by suppression of myo-inositol phosphate kinases (phosphorylation) and alternative overexpression of phosphatases (dephosphorylation). Our study suggests that metabolome and transcriptome information taken together, provide valuable information about biomarkers that can be used as a tool for marker-assisted breeding, metabolomics and functional genomics of this important legume crop.

Comparative transcriptome and metabolome profiling in the maturing seeds of contrasting cluster bean (Cyamopsis tetragonoloba L. Taub) cultivars identified key molecular variations leading to increased gum accumulation.

Cluster bean (Guar) is the major source of industrial gum. Knowledge on the molecular events regulating galactomannan gum accumulation in guar will pave way for accelerated development of gummy guar genotypes. RNA Seq analysis in the immature seeds of contrasting cluster bean genotypes HGS 563 (gum type) and Pusa Navbahar (vegetable type) resulted in the generation of 19,855,490 and 21,488,472 quality reads. Data analysis identified 4938 differentially expressed genes between the gummy vs vegetable genotypes. A set of 2241 genes were up-regulated and 2587 genes were down-regulated in gummy guar. Significant up-regulation of genes involved in the biosynthesis of galactomannan and cell wall storage polysaccharides was observed in the gummy HGS 563. Genes involved in carotenoids, flavonoids, non mevalonic acid, terpenoids, and wax metabolism were also up-regulated in HGS 563. Mannose and galactose were the major nucleotide sugars in Pusa Navbahar and HGS 563 immature seeds. Immature seeds of HGS 563 showed high concentration of mannose and galactose accumulation compared to Pusa Navbahar. qRT-PCR analysis of selected genes confirmed the findings of transcriptome data.

Elucidation of Galactomannan Biosynthesis Pathway Genes through Transcriptome Sequencing of Seeds Collected at Different Developmental Stages of Commercially Important Indian Varieties of Cluster Bean (Cyamopsis tetragonoloba L.).

Cyamopsis tetragonoloba (L) endosperm predominantly contains guar gum a polysaccharide, which has tremendous industrial applications in food, textile, paper, oil drilling and water treatment. In order to understand the genes controlling galactomannan biosynthesis, mRNA was isolated from seeds collected at different developmental stages; young pods, mature pods and young leaf from two guar varieties, HG365 and HG870 and subjected to Illumina sequencing. De novo assembly of fourteen individual read files from two varieties of guar representing seven developmental stages gave a total of 1,13,607 contigs with an N50 of 1,244 bases. Annotation of assemblies with GO mapping revealed three levels of distribution, namely, Biological Processes, Molecular Functions and Cellular Components. GO studies identified major genes involved in galactomannan biosynthesis: Cellulose synthase D1 (CS D1) and GAUT-like gene families. Among the polysaccharide biosynthetic process (GO:0000271) genes the transcript abundance for CS was found to be predominantly more in leaf samples, whereas, the transcript abundance for GAUT-like steadily increased from 65% to 90% and above from stage1 to stage5 indicating accumulation of galactomannan in developing seeds; and validated by qRT-PCR analysis. Galactomannan quantification by HPLC showed HG365 (12.98-20.66%) and HG870 (7.035-41.2%) gradually increasing from stage1 to stage 5 (10-50 DAA) and highest accumulation occurred in mature and dry seeds with 3.8 to 7.1 fold increase, respectively. This is the first report of transcriptome sequencing and complete profiling of guar seeds at different developmental stages, young pods, mature pods and young leaf material from two commercially important Indian varieties and elucidation of galactomannan biosynthesis pathway. It is envisaged that the data presented herein will be very useful for improvement of guar through biotechnological interventions in future.

Characterization of genes in guar gum biosynthesis based on quantitative RNA-sequencing in guar bean (Cyamopsis tetragonoloba).

Guar gum is an important raw material in the food, textile and oil industries, but the biosynthesis of guar gum remains unclear. To illuminate the genes involved in guar gum biosynthesis, guar beans from 30 and 40 days after flowering (DAF) were used for RNA sequencing in this study. A total of 2,535 and 2,724 preferentially expressed genes were found in 30 and 40 DAF endosperm, and 3,720 and 2,530 preferentially expressed genes were found in 30 and 40 DAF embryos, respectively. Of these, mannan synthase genes, α-galactosyltransferase genes and cellulose synthase genes were preferentially expressed in the endosperm from 30 and 40 DAF. The high expression level of these glycometabolism genes in endosperm is consistent with the expectation that the main component of guar gum is galactomannan. We believe that genes related to guar gum biosynthesis found in this study will be useful for both new variety development via genetic engineering and synthetic biology research on guar gum biosynthesis in the future.

Transcriptome and metabolome analysis of Ferula gummosa Boiss. to reveal major biosynthetic pathways of galbanum compounds.

Ferula gummosa Boiss. is an industrial and pharmaceutical plant that has been highly recognized for its valuable oleo-gum-resin, namely galbanum. Despite the fabulous value of galbanum, very little information on the genetic and biochemical mechanisms of its production existed. In the present study, the oleo-gum-resin and four organs (root, flower, stem, and leaf) of F. gummosa were assessed in terms of metabolic compositions and the expression of genes involved in their biosynthetic pathways. Results showed that the most accumulation of resin and essential oils were occurred in the roots (13.99 mg/g) and flowers (6.01 mg/g), respectively. While the most dominant compound of the resin was ß-amyrin from triterpenes, the most abundant compounds of the essential oils were α-pinene and ß-pinene from monoterpenes and α-eudesmol and germacrene-D from sesquiterpenes. Transcriptome analysis was performed by RNA sequencing (RNA-seq) for the plant roots and flowers. Differential gene expression analysis showed that 1172 unigenes were differential between two organs that 934 (79.6%) of them were up-regulated in the flowers and 238 (20.4%) unigenes were up-regulated in the roots (FDR ≤0.001). The most important up-regulated unigenes in the roots were involved in the biosynthesis of the major components of galbanum, including myrcene, germacrene-D, α-terpineol, and ß-amyrin. The results obtained by RNA-Seq were confirmed by qPCR. These analyses showed that different organs of F. gummosa are involved in the production of oleo-gum-resin, but the roots are more active than other organs in terms of the biosynthesis of triterpenes and some mono- and sesquiterpenes. This study provides rich molecular and biochemical resources for further studies on molecular genetics and functional genomics of oleo-gum-resin production in F. gummosa.