Molecular and morphophysiological characterization of superior cluster bean (Cymopsis tetragonoloba) varieties

Differential responses of two local and commercial guar cultivars for nutrient uptake and yield components under drought and biochar application

Drought is one of the abiotic stresses that can reduce crop yields. It has a major impact on crop yield reduction. For crops under stress, organic modifiers such as biochar can be useful. Guar (Cyamopsis tetragonoloba L.), an annual legume from the Fabaceae Family, is highly adaptable to arid and semi-arid regions, with many applications in various industries. Field experiments were carried out in a randomized complete block design with three replications using a split-split plots arrangement. The aim was to evaluate the influence of irrigation levels (Ir1 = 10, Ir2 = 14, and Ir3 = 17 days irrigation cycle) and biochar (B1 = 0, B2 = 5, and B3 = 10 tons ha-1) application on physiological traits [(chlorophyll a and b, chlorophyll index (SPAD), relative leaf water content (RWC), electrolyte leakage (EL), canopy temperature, leaf area, water use efficiency (WUE)], morphological parameters (length and diameter of the stem, pod length, fresh weight of root and plant, root length), yield components (seed yield, number of branch plant-1, number of clusters plant-1, pod plant-1, seed pod-1, seed plant-1, 1000-seed weight, and gum contents), and leaf nutrient uptake (Ca, Mg, P, Na, and K) of two commercial and local cultivars (cv1 = RGC-936 and cv2 = Saravan) of the guar plant. It was observed that the Ir3 irrigation treatment produced the highest seed yield (1921.8 kg ha-1) in terms of water stress. However, the maximum pod plant-1 (75.5), seed plant-1 (454.2), seed yield (1871.1 kg ha-1), leaf area (861.8 mm2), SPAD (92.2), Mg (49.8 mg g-1), Na (43.3 mg g-1) and P (0.49 mg g-1) were observed in RGC-936. The results also revealed that biochar was more effective than cultivars in terms of morphological traits. While yield and yield components were affected by cultivar, irrigation at different levels also had a significant effect on functional traits, physiology, and morphology. The addition of biochar appeared to have a positive effect on water stress alleviation and guar growth and leaf nutrient uptake. According to Pearson's correlation analysis, plant weight and length, root weight and length, stem diameter, seed pod-1, branches plant-1, and 1000-seed weight are moderately correlated with seed yield, while pod plant-1 and seed plant-1 are strongly associated with seed yield. On the other hand, the pod length, branches plant-1, and gum content showed a positive but not significant relationship.

Improvement in growth and yield attributes of cluster bean through optimization of sowing time and plant spacing under climate change scenario

Cluster bean (Cyamopsis tetragonoloba L.) yield has plateaued due to reduction in rainfall and rise in temperature. Therefore, its production cycle could not get appropriate water and temperature. It becomes important to standardize the sowing time and plant spacing of cluster beans in changing climate scenarios to get higher productivity. Therefore, a field study was conducted in 2019 at the Research area of MNS-University of Agriculture, Multan, Pakistan to evaluate the effect of four sowing times (15th May, 1st June, 15th June, and 1st July) and three plant spacings (10, 12 and 15 cm) on crop growth, yield, and physiological functions of cluster bean genotype BR-2017 under split plot arrangement under randomized complete block design (RCBD) with three replications. The sowing times (15th May, 1st June, 15th June, and 1st July) were placed in the main plot, while plant spacing (10, 12 and 15 cm) was maintained in subplots. The significant effect of sowing time and plant spacing was observed on pod plant⁻¹, pod length, grain yield, and 1000-grain weight. Results showed that 1st June sowing performed better over 15th May, 15th June, and 1st July, while plant spacing 15 cm about in all sowing times showed higher results on growth and yield parameters of cluster bean over plant spacing 10, 12, and 15 cm. The 1st June sowing time at 15 cm plant spacing showed 8.0, 22.7, and 28.5% higher grains pod^-1 than 15th May, 15th June, and 1st July sowing, respectively. Maximum grain yield was observed on 1st June in all three spacings (10, 12, and 15 cm). The chord diagram indicates that the crop has received optimum environmental conditions when sown 1st June over other sowing times. In conclusion, 1st June sowing with 15 cm plant spacing could be a good option to achieve maximum productivity of cluster bean under changing climate scenario.

De Novo Assembly and Species-Specific Marker Development as a Useful Tool for the Identification of Scutellaria L. Species

Scutellaria L. (family Lamiaceae) includes approximately 470 species found in most parts of the world and is commonly known as skullcaps. Scutellaria L. is a medicinal herb used as a folk remedy in Korea and East Asia, but it is difficult to identify and classify various subspecies by morphological methods. Since Scutellaria L. has not been studied genetically, to expand the knowledge of species in the genus Scutellaria L., de novo whole-genome assembly was performed in Scutellaria indica var. tsusimensis (H. Hara) Ohwi using the Illumina sequencing platform. We aimed to develop a molecular method that could be used to classify S.indica var. tsusimensis (H. Hara) Ohwi, S. indica L. and three other Scutellaria L. species. The assembly results for S.indica var. tsusimensis (H. Hara) Ohwi revealed a genome size of 318,741,328 bp and a scaffold N50 of 78,430. The assembly contained 92.08% of the conserved BUSCO core gene set and was estimated to cover 94.65% of the genome. The obtained genes were compared with previously registered Scutellaria nucleotide sequences and similar regions using the NCBI BLAST service, and a total of 279 similar nucleotide sequences were detected. By selecting the 279 similar nucleotide sequences and nine chloroplast DNA barcode genes, primers were prepared so that the size of the PCR product was 100 to 1000 bp. As a result, a species-specific primer set capable of distinguishing five species of Scutellaria L. was developed.

Recent Progress in Chemical Modification of the Natural Polysaccharide Guar Gum

Guar gum (GG) is a natural heteropolysaccharide. Due to its non-toxic, eco-friendly, and biodegradable nature, GG has found wide applications in many areas, in particular food, paper, textile, petroleum, and pharmaceutical industries. Therefore, GG is often called "Black Gold" as well. Due to the presence of hydroxyl groups, GG can be modified by various methods. The physical and biological properties of GG can be modulated by chemical modifications. In this manuscript, various methods for the chemical modifications of GG have been discussed according to the type of modifications. Mechanistic insights have also been provided whenever possible. In addition, potential applications of new GG derivatives have also been briefly mentioned.

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.

High temperature induced changes in quality and yield parameters of tomato (Solanum lycopersicum L.) and similarity coefficients among genotypes using SSR markers

High temperature induced by climatic fluctuations are an important threat for plant growth, development and quality of agricultural produces. Adaptableness to environmental changes generally derives from a large set of genetic traits affecting physio-morphological, biochemical and agronomic parameters. Therefore, the identification of genotypes with higher yield and good quality parameters at high temperatures is becoming increasingly necessary for future breeding programs. Here, we analyzed the performance of different tomato genotypes grown under elevated temperatures in terms of yield and nutritional quality of the fruit. High temperature stress was induced from flower initiation to maturity stage by keeping the pots in a temperature controlled green house facility for 45 days. The quality and yield parameters were taken at the harvesting stage. Starch and soluble sugar concentration in the leaves of tomato genotypes showed significant reduction in its amount under heat stress. Titrable acidity (TA), total soluble solids (TSS) and ascorbic acid content of tomato fruits were highest under high temperature conditions compared to ambient condition but lycopene content decreased with rise in temperature. The yield attributes viz., number of fruits/plant, fruit set %, average fruit weight (g), yield per plant (g/plant) were significantly lower for Arka Saurabh, Arka Rakshak and Pusa Rohini when compared to other genotypes under study. Molecular characterization of selected 22 tomato genotypes were assessed using 25 simple sequence repeat (SSR) markers. Phylogenetic tree was constructed by the unweighted neighbour-joining method (UPGMA) using NTSYSpc cluster analysis software. The Jaccard's similarity matrix was constructed using the SIMQUAL method using UPGMA algorithm in NTSYSpc. Jaccard's similarity matrix among these tomato genotypes ranged from a minimum of 0.22 to a maximum of 1 with an average genetic similarity of 0.67. Hence this study has importance in identifying genotypes that could maintain good quality and higher yield under high temperature condition.

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.

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.

Mining and Development of Novel SSR Markers Using Next Generation Sequencing (NGS) Data in Plants

Microsatellites, or simple sequence repeats (SSRs), are one of the most informative and multi-purpose genetic markers exploited in plant functional genomics. However, the discovery of SSRs and development using traditional methods are laborious, time-consuming, and costly. Recently, the availability of high-throughput sequencing technologies has enabled researchers to identify a substantial number of microsatellites at less cost and effort than traditional approaches. Illumina is a noteworthy transcriptome sequencing technology that is currently used in SSR marker development. Although 454 pyrosequencing datasets can be used for SSR development, this type of sequencing is no longer supported. This review aims to present an overview of the next generation sequencing, with a focus on the efficient use of de novo transcriptome sequencing (RNA-Seq) and related tools for mining and development of microsatellites in plants.

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

Clusterbean (Cyamopsis tetragonoloba L. Taub), is an important industrial, vegetable and forage crop. This crop owes its commercial importance to the presence of guar gum (galactomannans) in its endosperm which is used as a lubricant in a range of industries. Despite its relevance to agriculture and industry, genomic resources available in this crop are limited. Therefore, the present study was undertaken to generate RNA-Seq based transcriptome from leaf, shoot, and flower tissues. A total of 145 million high quality Illumina reads were assembled using Trinity into 127,706 transcripts and 48,007 non-redundant high quality (HQ) unigenes. We annotated 79% unigenes against Plant Genes from the National Center for Biotechnology Information (NCBI), Swiss-Prot, Pfam, gene ontology (GO) and KEGG databases. Among the annotated unigenes, 30,020 were assigned with 116,964 GO terms, 9984 with EC and 6111 with 137 KEGG pathways. At different fragments per kilobase of transcript per millions fragments sequenced (FPKM) levels, genes were found expressed higher in flower tissue followed by shoot and leaf. Additionally, we identified 8687 potential simple sequence repeats (SSRs) with an average frequency of one SSR per 8.75 kb. A total of 28 amplified SSRs in 21 clusterbean genotypes resulted in polymorphism in 13 markers with average polymorphic information content (PIC) of 0.21. We also constructed a database named ‘ClustergeneDB’ for easy retrieval of unigenes and the microsatellite markers. The tissue specific genes identified and the molecular marker resources developed in this study is expected to aid in genetic improvement of clusterbean for its end use.

RNA-Seq of Guar (Cyamopsis tetragonoloba, L. Taub.) Leaves: De novo Transcriptome Assembly, Functional Annotation and Development of Genomic Resources

Genetic improvement in industrially important guar (Cyamopsis tetragonoloba, L. Taub.) crop has been hindered due to the lack of sufficient genomic or transcriptomic resources. In this study, RNA-Seq technology was employed to characterize the transcriptome of leaf tissues from two guar varieties, namely, M-83 and RGC-1066. Approximately 30 million high-quality pair-end reads of each variety generated by Illumina HiSeq platform were used for de novo assembly by Trinity program. A total of 62,146 non-redundant unigenes with an average length of 679 bp were obtained. The quality assessment of assembled unigenes revealed 87.50% of complete and 97.18% partial core eukaryotic genes (CEGs). Sequence similarity analyses and annotation of the unigenes against non-redundant protein (Nr) and Gene Ontology (GO) databases identified 175,882 GO annotations. A total of 11,308 guar unigenes were annotated with various enzyme codes (EC) and categorized in six categories with 55 subclasses. The annotation of biochemical pathways resulted in a total of 11,971 unigenes assigned with 145 KEGG maps and 1759 enzyme codes. The species distribution analysis of the unigenes showed highest similarity with Glycine max genes. A total of 5773 potential simple sequence repeats (SSRs) and 3594 high-quality single nucleotide polymorphisms (SNPs) were identified. Out of 20 randomly selected SSRs for wet laboratory validation, 13 showed consistent PCR amplification in both guar varieties. In silico studies identified 145 polymorphic SSR markers in two varieties. To the best of our knowledge, this is the first report on transcriptome analysis and SNPs identification in guar till date.