High Quality Unigenes and Microsatellite Markers from Tissue Specific Transcriptome and Development of a Database in Clusterbean (Cyamopsis tetragonoloba, L. Taub)

The chromosome-scale genome assembly of cluster bean provides molecular insight into edible gum (galactomannan) biosynthesis family genes

Cluster bean (Cyamopsis tetragonoloba (L.) Taub 2n = 14, is commonly known as Guar. Apart from being a vegetable crop, it is an abundant source of a natural hetero-polysaccharide called guar gum or galactomannan. Here, we are reporting a chromosome-scale reference genome assembly of a popular cluster bean cultivar RGC-936, by combining sequencing data from Illumina, 10X Genomics, Oxford Nanopore technologies. An initial assembly of 1580 scaffolds with an N50 value of 7.12 Mb was generated and these scaffolds were anchored to a high density SNP linkage map. Finally, a genome assembly of 550.31 Mb (94% of the estimated genome size of ~ 580 Mb (through flow cytometry) with 58 scaffolds was obtained, including 7 super scaffolds with a very high N50 value of 78.27 Mb. Phylogenetic analysis using single copy orthologs among 12 angiosperms showed that cluster bean shared a common ancestor with other legumes 80.6 MYA. No evidence of recent whole genome duplication event in cluster bean was found in our analysis. Further comparative transcriptomics analyses revealed pod-specific up-regulation of genes encoding enzymes involved in galactomannan biosynthesis. The high-quality chromosome-scale cluster bean genome assembly will facilitate understanding of the molecular basis of galactomannan biosynthesis and aid in genomics-assisted improvement of cluster bean.

Development of de-novo transcriptome assembly and SSRs in allohexaploid Brassica with functional annotations and identification of heat-shock proteins for thermotolerance

Crop Brassicas contain monogenomic and digenomic species, with no evidence of a trigenomic Brassica in nature. Through somatic fusion (Sinapis alba + B. juncea), a novel allohexaploid trigenomic Brassica (H1 = AABBSS; 2n = 60) was produced and used for transcriptome analysis to uncover genes for thermotolerance, annotations, and microsatellite markers for future molecular breeding. Illumina Novaseq 6000 generated a total of 76,055,546 paired-end raw reads, which were used for de-novo assembly, resulting in the development of 486,066 transcripts. A total of 133,167 coding sequences (CDSs) were predicted from transcripts with a mean length of 507.12 bp and 46.15% GC content. The BLASTX search of CDSs against public protein databases showed a maximum of 126,131 (94.72%) and a minimum of 29,810 (22.39%) positive hits. Furthermore, 953,773 gene ontology (GO) terms were found in 77,613 (58.28%) CDSs, which were divided into biological processes (49.06%), cellular components (31.67%), and molecular functions (19.27%). CDSs were assigned to 144 pathways by a pathway study using the KEGG database and 1,551 pathways by a similar analysis using the Reactome database. Further investigation led to the discovery of genes encoding over 2,000 heat shock proteins (HSPs). The discovery of a large number of HSPs in allohexaploid Brassica validated our earlier findings for heat tolerance at seed maturity. A total of 15,736 SSRs have been found in 13,595 CDSs, with an average of one SSR per 4.29 kb length and an SSR frequency of 11.82%. The first transcriptome assembly of a meiotically stable allohexaploid Brassica has been given in this article, along with functional annotations and the presence of SSRs, which could aid future genetic and genomic studies.

In silico Identification of miRNAs and Their Targets in Cluster Bean for Their Role in Development and Physiological Responses

Cluster bean popularly known as “guar” is a drought-tolerant, annual legume that has recently emerged as an economically important crop, owing to its high protein and gum content. The guar gum has wide range of applications in food, pharma, and mining industries. India is the leading exporter of various cluster bean-based products all across the globe. Non-coding RNAs (miRNAs) are involved in regulating the expression of the target genes leading to variations in the associated pathways or final protein concentrations. The understanding of miRNAs and their associated targets in cluster bean is yet to be used to its full potential. In the present study, cluster bean EST (Expressed Sequence Tags) database was exploited to identify the miRNA and their predicted targets associated with metabolic and biological processes especially response to diverse biotic and abiotic stimuli using in silico approach. Computational analysis based on cluster bean ESTs led to the identification of 57 miRNAs along with their targets. To the best of our knowledge, this is the first report on identification of miRNAs and their targets using ESTs in cluster bean. The miRNA related to gum metabolism was also identified. Most abundant miRNA families predicted in our study were miR156, miR172, and miR2606. The length of most of the mature miRNAs was found to be 21nt long and the range of minimal folding energy (MFE) was 5.8–177.3 (−kcal/mol) with an average value of 25.4 (−kcal/mol). The identification of cluster bean miRNAs and their targets is predicted to hasten the miRNA discovery, resulting in better knowledge of the role of miRNAs in cluster bean development, physiology, and stress responses.

De novo Transcriptome Analysis of Drought-Adapted Cluster Bean (Cultivar RGC-1025) Reveals the Wax Regulatory Genes Involved in Drought Resistance

Cluster bean (Cyamopsis tetragonoloba L.) is one of the multipurpose underexplored crops grown as green vegetable and for gum production in dryland areas. Cluster bean is known as relatively tolerant to drought and salinity stress. To elucidate the molecular mechanisms involved in the drought tolerance of cluster bean cultivar RGC-1025, RNA sequencing (RNA-seq) of the drought-stressed and control samples was performed. De novo assembly of the reads resulted in 66,838 transcripts involving 203 pathways. Among these transcripts, differentially expressed gene (DEG) analysis resulted in some of the drought-responsive genes expressing alpha dioxygenase 2, low temperature-induced 65 kDa protein (LDI65), putative vacuolar amino acid transporter, and late embryogenesis abundant protein (LEA 3). The analysis also reported drought-responsive transcription factors (TFs), such as NAC, WRKY, GRAS, and MYB families. The relative expression of genes by qRT-PCR revealed consistency with the DEG analysis. Key genes involved in the wax biosynthesis pathway were mapped using the DEG data analysis. These results were positively correlated with epicuticular wax content and the wax depositions on the leaf surfaces, as evidenced by scanning electron microscope (SEM) image analysis. Further, these findings support the fact that enhanced wax deposits on the leaf surface had played a crucial role in combating the drought stress in cluster beans under drought stress conditions. In addition, this study provided a set of unknown genes and TFs that could be a source of engineering tolerance against drought stress in cluster beans.

De novo transcriptome assembly and analysis of gene expression in different tissues of moth bean (Vigna aconitifolia) (Jacq.) Marechal

BACKGROUND

The underutilized species Vigna aconitifolia (Moth Bean) is an important legume crop cultivated in semi-arid conditions and is valued for its seeds for their high protein content. It is also a popular green manure cover crop that offers many agronomic benefits including nitrogen fixation and soil nutrients. Despite its economic potential, genomic resources for this crop are scarce and there is limited knowledge on the developmental process of this plant at a molecular level. In the present communication, we have studied the molecular mechanisms that regulate plant development in V. aconitifolia, with a special focus on flower and seed development. We believe that this study will greatly enrich the genomic resources for this plant in form of differentially expressed genes, transcription factors, and genic molecular markers.

RESULTS

We have performed the de novo transcriptome assembly using six types of tissues from various developmental stages of Vigna aconitifolia (var. RMO-435), namely, leaves, roots, flowers, pods, and seed tissue in the early and late stages of development, using the Illumina NextSeq platform. We assembled the transcriptome to get 150938 unigenes with an average length of 937.78 bp. About 79.9% of these unigenes were annotated in public databases and 12839 of those unigenes showed a significant match in the KEGG database. Most of the unigenes displayed significant differential expression in the late stages of seed development as compared with leaves. We annotated 74082 unigenes as transcription factors and identified 12096 simple sequence repeats (SSRs) in the genic regions of V.aconitifolia. Digital expression analysis revealed specific gene activities in different tissues which were validated using Real-time PCR analysis.

CONCLUSIONS

The Vigna aconitifolia transcriptomic resources generated in this study provide foundational resources for gene discovery with respect to various developmental stages. This study provides the first comprehensive analysis revealing the genes involved in molecular as well as metabolic pathways that regulate seed development and may be responsible for the unique nutritive values of moth bean seeds. Hence, this study would serve as a foundation for characterization of candidate genes which would not only provide novel insights into understanding seed development but also provide resources for improved moth bean and related species genetic enhancement.

Comparative transcriptomic analysis and antioxidant defense mechanisms in clusterbean (Cyamopsis tetragonoloba (L.) Taub.) genotypes with contrasting drought tolerance

To understand drought tolerance mechanism(s) in clusterbean (Cyamopsis tetragonoloba), we conducted physiological, biochemical, and de novo comparative transcriptome analysis of drought-tolerant (RGC-1002) and drought-sensitive (RGC-1066) genotypes subjected to 30 days of drought stress. Relative water content (RWC) was maintained in tolerant genotype but was reduced in sensitive genotype. Leaf pigment concentrations were higher in tolerant genotype. Net photosynthesis was significantly decreased in sensitive genotype but insignificant reduction was found in tolerant genotype. Enzymatic antioxidant (GR, APX, DHAR) activities were enhanced in tolerant genotype, while there were insignificant changes in these enzymes in sensitive genotype. The ratios of antioxidant molecules (ASC/DHA and GSH/GSSG) were higher in tolerant genotype as compared to sensitive genotype. In sensitive genotype, 6625 differentially expressed genes (DEGs) were upregulated and 5365 genes were downregulated. In tolerant genotype, 5206 genes were upregulated and 2793 genes were downregulated. In tolerant genotype, transketolase family protein, phosphoenolpyruvate carboxylase 3, temperature-induced lipocalin, and cytochrome oxidase were highly upregulated. Moreover, according to Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, the drought tolerance may be attributed to upregulated starch and sucrose metabolism-related genes in tolerant genotype. Finally, quantitative real-time PCR confirmed the reproducibility of the RNA-seq data.

Progress of Genomics-Driven Approaches for Sustaining Underutilized Legume Crops in the Post-Genomic Era

Legume crops, belonging to the Fabaceae family, are of immense importance for sustaining global food security. Many legumes are profitable crops for smallholder farmers due to their unique ability to fix atmospheric nitrogen and their intrinsic ability to thrive on marginal land with minimum inputs and low cultivation costs. Recent progress in genomics shows promise for future genetic gains in major grain legumes. Still it remains limited in minor legumes/underutilized legumes, including adzuki bean, cluster bean, horse gram, lathyrus, red clover, urd bean, and winged bean. In the last decade, unprecedented progress in completing genome assemblies of various legume crops and resequencing efforts of large germplasm collections has helped to identify the underlying gene(s) for various traits of breeding importance for enhancing genetic gain and contributing to developing climate-resilient cultivars. This review discusses the progress of genomic resource development, including genome-wide molecular markers, key breakthroughs in genome sequencing, genetic linkage maps, and trait mapping for facilitating yield improvement in underutilized legumes. We focus on 1) the progress in genomic-assisted breeding, 2) the role of whole-genome resequencing, pangenomes for underpinning the novel genomic variants underlying trait gene(s), 3) how adaptive traits of wild underutilized legumes could be harnessed to develop climate-resilient cultivars, 4) the progress and status of functional genomics resources, deciphering the underlying trait candidate genes with putative function in underutilized legumes 5) and prospects of novel breeding technologies, such as speed breeding, genomic selection, and genome editing. We conclude the review by discussing the scope for genomic resources developed in underutilized legumes to enhance their production and play a critical role in achieving the "zero hunger" sustainable development goal by 2030 set by the United Nations.

Deciphering Molecular Mechanisms Involved in Salinity Tolerance in Guar (Cyamopsis tetragonoloba (L.) Taub.) Using Transcriptome Analyses

Guar is a commercially important legume crop known for guar gum. Guar is tolerant to various abiotic stresses, but the mechanisms involved in its salinity tolerance are not well established. This study aimed to understand molecular mechanisms of salinity tolerance in guar. RNA sequencing (RNA-Seq) was employed to study the leaf and root transcriptomes of salt-tolerant (Matador) and salt-sensitive (PI 340261) guar genotypes under control and salinity. Our analyses identified a total of 296,114 unigenes assembled from 527 million clean reads. Transcriptome analysis revealed that the gene expression differences were more pronounced between salinity treatments than between genotypes. Differentially expressed genes associated with stress-signaling pathways, transporters, chromatin remodeling, microRNA biogenesis, and translational machinery play critical roles in guar salinity tolerance. Genes associated with several transporter families that were differentially expressed during salinity included ABC, MFS, GPH, and P-ATPase. Furthermore, genes encoding transcription factors/regulators belonging to several families, including SNF2, C₂H₂, bHLH, C3H, and MYB were differentially expressed in response to salinity. This study revealed the importance of various biological pathways during salinity stress and identified several candidate genes that may be used to develop salt-tolerant guar genotypes that might be suitable for cultivation in marginal soils with moderate to high salinity or using degraded water.

Development of SNP Set for the Marker-Assisted Selection of Guar (Cyamopsis tetragonoloba (L.) Taub.) Based on a Custom Reference Genome Assembly

Guar gum, a polysaccharide derived from guar seeds, is widely used in a variety of industrial applications, including oil and gas production. Although guar is mostly propagated in India, interest in guar as a new industrial legume crop is increasing worldwide, demanding the development of effective tools for marker-assisted selection. In this paper, we report a wide-ranging set of 4907 common SNPs and 327 InDels generated from RADseq genotyping data of 166 guar plants of different geographical origin. A custom guar reference genome was assembled and used for variant calling. A consensus set of variants was built using three bioinformatic pipelines for short variant discovery. The developed molecular markers were used for genome-wide association study, resulting in the discovery of six markers linked to the variation of an important agronomic trait—percentage of pods matured to the harvest date under long light day conditions. One of the associated variants was found inside the putative transcript sequence homologous to an ABC transporter in Arabidopsis, which has been shown to play an important role in D-myo-inositol phosphates metabolism. Earlier, we suggested that genes involved in myo-inositol phosphate metabolism have significant impact on the early flowering of guar plants. Hence, we believe that the developed SNP set allows for the identification of confident molecular markers of important agrobiological traits.

Morphological, physiological, and transcriptional responses to low nitrogen stress in Populus deltoides Marsh. clones with contrasting nitrogen use efficiency

Background

Nitrogen (N) is one of the main factors limiting the wood yield in poplar cultivation. Understanding the molecular mechanism of N utilization could play a guiding role in improving the nitrogen use efficiency (NUE) of poplar.

Results

In this study, three N-efficient genotypes (A1-A3) and three N-inefficient genotypes (C1-C3) of Populus deltoides were cultured under low N stress (5 μM NH₄NO₃) and normal N supply (750 μM NH₄NO₃). The dry matter mass, leaf morphology, and chlorophyll content of both genotypes decreased under N starvation. The low nitrogen adaptation coefficients of the leaves and stems biomass of group A were significantly higher than those of group C (p < 0.05). Interestingly, N starvation induced fine root growth in group A, but not in group C. Next, a detailed time-course analysis of enzyme activities and gene expression in leaves identified 2062 specifically differentially expressed genes (DEGs) in group A and 1118 in group C. Moreover, the sensitivity to N starvation of group A was weak, and DEGs related to hormone signal transduction and stimulus response played an important role in the low N response this group. Weighted gene co-expression network analysis identified genes related to membranes, catalytic activity, enzymatic activity, and response to stresses that might be critical for poplar’s adaption to N starvation and these genes participated in the negative regulation of various biological processes. Finally, ten influential hub genes and twelve transcription factors were identified in the response to N starvation. Among them, four hub genes were related to programmed cell death and the defense response, and PodelWRKY18, with high connectivity, was involved in plant signal transduction. The expression of hub genes increased gradually with the extension of low N stress time, and the expression changes in group A were more obvious than those in group C.

Conclusions

Under N starvation, group A showed stronger adaptability and better NUE than group C in terms of morphology and physiology. The discovery of hub genes and transcription factors might provide new information for the analysis of the molecular mechanism of NUE and its improvement in poplar.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07991-7.

Genome wide identification and characterization of abiotic stress responsive lncRNAs in Capsicum annuum

Long non-coding RNAs (lncRNAs) are a type of non-coding transcripts having length of more than 200 nucleotides lacking protein-coding ability. In the present study, 12807 lncRNAs were identified in Capsicum annuum tissues exposed to abiotic stress conditions viz. heat, cold, osmotic and salinity stress. Expression analysis of lncRNAs in different treatment conditions demonstrates their stress-specific expression. Thirty lncRNAs were found to act as precursors for 10 microRNAs (miRNAs) of C. annuum. Additionally, a total of 1807 lncRNAs were found to interact with 194 miRNAs which targeted 621 mRNAs of C. annuum. Among these, 344 lncRNAs were found to act as target mimics for 621 genes. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that out of those 621 gene sequences, 546 were tagged with GO terms, 105 Enzyme Code (EC) numbers were assigned to 246 genes and 223 genes are found to be involved in 63 biological pathways. In this report, we have highlighted the prospective role of lncRNAs in different abiotic stress conditions by interacting with miRNAs and regulating stress responsive transcription factors (TFs) such as DoF, WRKY, MYB, bZIP and ERF in C. annuum.

Exploring the edible gum (galactomannan) biosynthesis and its regulation during pod developmental stages in clusterbean using comparative transcriptomic approach

Galactomannan is a polymer of high economic importance and is extracted from the seed endosperm of clusterbean (C. tetragonoloba). In the present study, we worked to reveal the stage-specific galactomannan biosynthesis and its regulation in clusterbean. Combined electron microscopy and biochemical analysis revealed high protein and gum content in RGC-936, while high oil bodies and low gum content in M-83. A comparative transcriptome study was performed between RGC-936 (high gum) and M-83 (low gum) varieties at three developmental stages viz. 25, 39, and 50 days after flowering (DAF). Total 209,525, 375,595 and 255,401 unigenes were found at 25, 39 and 50 DAF respectively. Differentially expressed genes (DEGs) analysis indicated a total of 5147 shared unigenes between the two genotypes. Overall expression levels of transcripts at 39DAF were higher than 50DAF and 25DAF. Besides, 691 (RGC-936) and 188 (M-83) candidate unigenes that encode for enzymes involved in the biosynthesis of galactomannan were identified and analyzed, and 15 key enzyme genes were experimentally validated by quantitative Real-Time PCR. Transcription factor (TF) WRKY was observed to be co-expressed with key genes of galactomannan biosynthesis at 39DAF. We conclude that WRKY might be a potential biotechnological target (subject to functional validation) for developing high gum content varieties.

Developments and Prospects in Imperative Underexploited Vegetable Legumes Breeding: A Review

Vegetable legumes are an essential source of carbohydrates, vitamins, and minerals, along with health-promoting bioactive chemicals. The demand for the use of either fresh or processed vegetable legumes is continually expanding on account of the growing consumer awareness about their well-balanced diet. Therefore, sustaining optimum yields of vegetable legumes is extremely important. Here we seek to present d etails of prospects of underexploited vegetable legumes for food availability, accessibility, and improved livelihood utilization. So far research attention was mainly focused on pulse legumes' performance as compared to vegetable legumes. Wild and cultivated vegetable legumes vary morphologically across diverse habitats. This could make them less known, underutilized, and underexploited, and make them a promising potential nutritional source in developing nations where malnutrition still exists. Research efforts are required to promote underexploited vegetable legumes, for improving their use to feed the ever-increasing population in the future. In view of all the above points, here we have discussed underexploited vegetable legumes with tremendous potential; namely, vegetable pigeon pea (Cajanus cajan), cluster bean (Cyamopsis tetragonoloba), winged bean (Psophocarpus tetragonolobus), dolichos bean (Lablab purpureus), and cowpea (Vigna unguiculata), thereby covering the progress related to various aspects such as pre-breeding, molecular markers, quantitative trait locus (QTLs), genomics, and genetic engineering. Overall, this review has summarized the information related to advancements in the breeding of vegetable legumes which will ultimately help in ensuring food and nutritional security in developing nations.

Fate of nitrogen in agriculture and environment: agronomic, eco-physiological and molecular approaches to improve nitrogen use efficiency

Nitrogen is the main limiting nutrient after carbon, hydrogen and oxygen for photosynthetic process, phyto-hormonal, proteomic changes and growth-development of plants to complete its lifecycle. Excessive and inefficient use of N fertilizer results in enhanced crop production costs and atmospheric pollution. Atmospheric nitrogen (71%) in the molecular form is not available for the plants. For world's sustainable food production and atmospheric benefits, there is an urgent need to up-grade nitrogen use efficiency in agricultural farming system. The nitrogen use efficiency is the product of nitrogen uptake efficiency and nitrogen utilization efficiency, it varies from 30.2 to 53.2%. Nitrogen losses are too high, due to excess amount, low plant population, poor application methods etc., which can go up to 70% of total available nitrogen. These losses can be minimized up to 15-30% by adopting improved agronomic approaches such as optimal dosage of nitrogen, application of N by using canopy sensors, maintaining plant population, drip fertigation and legume based intercropping. A few transgenic studies have shown improvement in nitrogen uptake and even increase in biomass. Nitrate reductase, nitrite reductase, glutamine synthetase, glutamine oxoglutarate aminotransferase and asparagine synthetase enzyme have a great role in nitrogen metabolism. However, further studies on carbon-nitrogen metabolism and molecular changes at omic levels are required by using "whole genome sequencing technology" to improve nitrogen use efficiency. This review focus on nitrogen use efficiency that is the major concern of modern days to save economic resources without sacrificing farm yield as well as safety of global environment, i.e. greenhouse gas emissions, ammonium volatilization and nitrate leaching.

Comprehensive Stress-Based De Novo Transcriptome Assembly and Annotation of Guar (Cyamopsis tetragonoloba (L.) Taub.): An Important Industrial and Forage Crop

The forage crop Guar (Cyamopsis tetragonoloba (L.) Taub.) has the ability to endure heat, drought, and mild salinity. A complete image on its genic architecture will promote our understanding about gene expression networks and different tolerance mechanisms at the molecular level. Therefore, whole mRNA sequence approach on the Guar plant was conducted to provide a snapshot of the mRNA information in the cell under salinity, heat, and drought stresses to be integrated with previous transcriptomic studies. RNA-Seq technology was employed to perform a 2 × 100 paired-end sequencing using an Illumina HiSeq 2500 platform for the transcriptome of leaves of C. tetragonoloba under normal, heat, drought, and salinity conditions. Trinity was used to achieve a de novo assembly followed by gene annotation, functional classification, metabolic pathway analysis, and identification of SSR markers. A total of 218.2 million paired-end raw reads (~44 Gbp) were generated. Of those, 193.5M paired-end reads of high quality were used to reconstruct a total of 161,058 transcripts (~266 Mbp) with N50 of 2552 bp and 61,508 putative genes. There were 6463 proteins having >90% full-length coverage against the Swiss-Prot database and 94% complete orthologs against Embryophyta. Approximately, 62.87% of transcripts were blasted, 50.46% mapped, and 43.50% annotated. A total of 4715 InterProScan families, 3441 domains, 74 repeats, and 490 sites were detected. Biological processes, molecular functions, and cellular components comprised 64.12%, 25.42%, and 10.4%, respectively. The transcriptome was associated with 985 enzymes and 156 KEGG pathways. A total of 27,066 SSRs were gained with an average frequency of one SSR/9.825 kb in the assembled transcripts. This resulting data will be helpful for the advanced analysis of Guar to multi-stress tolerance.

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.

The genome size of clusterbean (Cyamopsis tetragonoloba) is significantly smaller compared to its wild relatives as estimated by flow cytometry

Clusterbean (C. tetragonoloba) is an important, leguminous vegetable and industrial crop with vast genetic diversity but meager genetic, cytological and genomic information. In the present study, an optimized procedure of flow cytometry was used to estimate the genome size of three clusterbean species, represented by C. tetragonoloba (cv. RGC-936) and two wild relatives (C. serreta and C. senegalensis). For accurate estimation of genomic content, singlet G₀/G₁ populations of multiple tissues such as leaves, hypocotyl, and matured seeds were determined and used along with three different plant species viz. Pisum sativum (as primary), Oryza sativa, and Glycine max (secondary), as external and internal reference standards. Seed tissue of the test sample and G. max provided the best estimate of nuclear DNA content in comparison to other sample tissues and reference standards. The genome size of C. tetragonoloba was detemined at 580.9±0.02Mbp (1C), while that of C. serreta and C. senegalensis was estimated at 979.6±0.02Mbp (1C) and 943.4±0.03Mbp (1C), respectively. Thus, the wild relatives harbor, nearly double the genome content of the cultivated cluster bean. Findings of this study will enrich genomic database of the legume family and can serve as the starting point for clusterbean evolutionary and genomics studies.

Identification and characterization of SSR, SNP and InDel molecular markers from RNA-Seq data of guar (Cyamopsis tetragonoloba, L. Taub.) roots

Background

Guar [Cyamopsis tetragonoloba, L. Taub.] is an important industrial crop because of the commercial applications of the galactomannan gum contained in its seeds. Plant breeding programmes based on marker-assisted selection require a rich resource of molecular markers. As limited numbers of such markers are available for guar, molecular breeding programmes have not been undertaken for the genetic improvement of this important crop. Hence, the present work was done to enrich the molecular markers resource of guar by identifying high quality SSR, SNP and InDel markers from the RNA-Seq data of the roots of two guar varieties.

Results

We carried out RNA-Seq analysis of the roots of two guar varieties, namely, RGC-1066 and M-83. A total of 102,479 unigenes with an average length of 1016 bp were assembled from about 30 million high quality pair-end reads generated by an Illumina HiSeq 2500 platform. The assembled unigenes had 86.55% complete and 97.71% partially conserved eukaryotic genes (CEGs). The functional annotation of assembled unigenes using BLASTX against six databases showed that the guar unigenes were most similar to Glycine max. We could assign GO terms to 45,200 unigenes using the UniProt database. The screening of 102,479 unigenes with MISA and SAMtools version 1.4 softwares resulted in the identification of 25,040 high-confidence molecular markers which consisted of 18,792 SSRs, 5999 SNPs and 249 InDels. These markers tagged most of the genes involved in root development, stress tolerance and other general metabolic activities. Each of the 25,040 molecular markers was characterized, particularly with respect to its position in the unigene. For 71% of the molecular markers, we could determine the names, products and functions of the unigenes. About 80% of the markers, from a random sample of molecular markers, showed PCR amplification.

Conclusions

We have identified and characterized 25,040 high confidence SSR, SNP and InDel molecular markers in guar. It is expected that these markers will be useful in molecular breeding programmes and will also be helpful in studying molecular mechanisms of root development, stress tolerance and gum synthesis in guar.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5205-9) contains supplementary material, which is available to authorized users.

V. AM, Solanke AU, Singh NK, Sharma TR, Gaikwad K. Genome-wide discovery of tissue-specific miRNAs in clusterbean (Cyamopsis tetragonoloba) indicates their association with galactomannan biosynthesis

Owing to the presence of 80% soluble dietary fibre, high protein content and high value gum, clusterbean (Cyamopsis tetragonoloba) has recently emerged as an economically important legume. The developing clusterbean seeds accumulate 90% galactomannans in the endosperm and, therefore, can be used as a model crop to understand galactomannan biosynthesis and its regulation. miRNAs are tiny master regulators of their corresponding target genes, resulting in variations in the amounts of their metabolic end products. To understand the role of these regulators in galactomannan biosynthesis regulation, small RNA libraries were prepared and sequenced from five tissues of clusterbean genotype RGC-936, and miRanalyzer and DSAP programs were used to identify conserved miRNAs and novel small RNAs. A total of 187 known and 171 novel miRNAs were found to be differentially expressed, of which 10 miRNAs were validated. A complicated network topology and 35% sharing of the target mRNAs between known and novel miRNAs suggest random evolution of novel miRNAs. The gene ontology (GO) annotation of potential target genes revealed the genes coding for signalling and carbohydrate metabolism (50.10%), kinases and other enzymes (20.75%), transcription factors (10.20%), transporters (8.35%) and other targets (10.6%). Two novel unigenes were annotated as ManS (mannosyltransferase/mannan synthase) and UGE (UDP- D-glucose 4-epimerase) and validated as targets for three novel miRNAs, that is Ct-miR3130, Ct-miR3135 and Ct-miR3157. Our findings reveal that these novel miRNAs could play an important role in the regulation of the galactomannan pathway in C. tetragonoloba and possibly other galactomannan-producing species.

Transcriptome Analysis of Two Rice Varieties Contrasting for Nitrogen Use Efficiency under Chronic N Starvation Reveals Differences in Chloroplast and Starch Metabolism-Related Genes

The nitrogen use efficiency (NUE) of crop plants is limited and enhancing it in rice, a major cereal crop, would be beneficial for farmers and the environment alike. Here we report the genome-wide transcriptome analysis of two rice genotypes, IR 64 (IR64) and Nagina 22 (N22) under optimal (N+) and chronic starvation (N−) of nitrogen (N) from 15-day-old root and shoot tissues. The two genotypes were found to be contrasting in their response to N−; IR64 root architecture and root dry weight remained almost equivalent to that under N+ conditions, while N22 showed high foraging ability but a substantial reduction in biomass under N−. Similarly, the photosynthetic pigments showed a drastic reduction in N22 under low N, while IR64 was more resilient. Nitrate reductase showed significantly low specific activity under N− in both genotypes. Glutamate synthase (GOGAT) and citrate synthase CS activity were highly reduced in N22 but not in IR64. Transcriptome analysis of these genotypes revealed nearly double the number of genes to be differentially expressed (DEGs) in roots (1016) compared to shoots (571). The response of the two genotypes to N starvation was distinctly different reflecting their morphological/biochemical response with just two and eight common DEGs in the root and shoot tissues. There were a total of 385 nitrogen-responsive DEGs (106 in shoots and 279 in roots) between the two genotypes. Fifty-two of the 89 DEGs identified as specific to N22 root tissues were also found to be differentially expressed between the two genotypes under N−. Most of these DEGs belonged to starch and chloroplast metabolism, followed by membrane and signaling proteins. Physical mapping of DEGs revealed 95 DEGs in roots and 76 in shoots to be present in quantitative trait loci (QTL) known for NUE.