Employing genome-wide SNP discovery and genotyping strategy to extrapolate the natural allelic diversity and domestication patterns in chickpea

Uncovering DNA methylation landscapes to decipher evolutionary footprints of phenotypic diversity in chickpea

Genetic diversity and environmental factors are long believed to be the dominant contributors to phenotypic diversity in crop plants. However, it has been recently established that, besides genetic variation, epigenetic variation, especially variation in DNA methylation, plays a significant role in determining phenotypic diversity in crop plants. Therefore, assessing DNA methylation diversity in crop plants becomes vital, especially in the case of crops like chickpea, which has a narrow genetic base. Thus, in the present study, we employed whole-genome bisulfite sequencing to assess DNA methylation diversity in wild and cultivated (desi and kabuli) chickpea. This revealed extensive DNA methylation diversity in both wild and cultivated chickpea. Interestingly, the methylation diversity was found to be significantly higher than genetic diversity, suggesting its potential role in providing vital phenotypic diversity for the evolution and domestication of the Cicer gene pool. The phylogeny based on DNA methylation variation also indicates a potential complementary role of DNA methylation variation in addition to DNA sequence variation in shaping chickpea evolution. Besides, the study also identified diverse epi-alleles of many previously known genes of agronomic importance. The Cicer MethVarMap database developed in this study enables researchers to readily visualize methylation variation within the genes and genomic regions of their interest (http://223.31.159.7/cicer/public/). Therefore, epigenetic variation like DNA methylation variation can potentially explain the paradox of high phenotypic diversity despite the narrow genetic base in chickpea and can potentially be employed for crop improvement.

Understanding genetic diversity in drought-adaptive hybrid parental lines in pearl millet

Information on genetic diversity and population structure is helpful to strategize enhancing the genetic base of hybrid parental lines in breeding programs. The present study determined the population structure and genetic diversity of 109 pearl millet hybrid parental lines, known for their better adaptation and performance in drought-prone environments, using 16,472 single nucleotide polymorphic (SNP) markers generated from GBS (genotyping-by-sequencing) platforms. The SNPs were distributed uniformly across the pearl millet genome and showed considerable genetic diversity (0.337), expected heterozygosity (0.334), and observed heterozygosity (0.031). Most of the pairs of lines (78.36%) had Identity-by-State (IBS) based genetic distances of more than 0.3, indicating a significant amount of genetic diversity among the parental lines. Bayesian model-based population stratification, neighbor-joining phylogenetic analysis, and principal coordinate analysis (PCoA) differentiated all hybrid parental lines into two clear-cut major groups, one each for seed parents (B-lines) and pollinators (R-lines). Majority of parental lines sharing common parentages were found grouped in the same cluster. Analysis of molecular variance (AMOVA) revealed 7% of the variation among subpopulations, and 93% of the variation was attributable to within sub-populations. Chromosome 3 had the highest number of LD regions. Genomic LD decay distance was 0.69 Mb and varied across the different chromosomes. Genetic diversity based on 11 agro-morphological and grain quality traits also suggested that the majority of the B- and R-lines were grouped into two major clusters with few overlaps. In addition, the combined analysis of phenotypic and genotypic data showed similarities in the population grouping patterns. The present study revealed the uniqueness of most of the inbred lines, which can be a valuable source of new alleles and help breeders to utilize these inbred lines for the development of hybrids in drought-prone environments.

Genome-wide association studies for earliness, MYMIV resistance, and other associated traits in mungbean (Vigna radiata L. Wilczek) using genotyping by sequencing approach

Yellow mosaic disease (YMD) remains a major constraint in mungbean (Vigna radiata (L.)) production; while short-duration genotypes offer multiple crop cycles per year and help in escaping terminal heat stress, especially during summer cultivation. A comprehensive genotyping by sequencing (GBS)-based genome-wide association studies (GWAS) analysis was conducted using 132 diverse mungbean genotypes for traits like flowering time, YMD resistance, soil plant analysis development (SPAD) value, trichome density, and leaf area. The frequency distribution revealed a wide range of values for all the traits. GBS studies identified 31,953 high-quality single nucleotide polymorphism (SNPs) across all 11 mungbean chromosomes and were used for GWAS. Structure analysis revealed the presence of two genetically distinct populations based on ΔK. The linkage disequilibrium (LD) varied throughout the chromosomes and at r2 = 0.2, the mean LD decay was estimated as 39.59 kb. Two statistical models, mixed linear model (MLM) and Bayesian-information and Linkage-disequilibrium Iteratively Nested Keyway (BLINK) identified 44 shared SNPs linked with various candidate genes. Notable candidate genes identified include FPA for flowering time (VRADI10G01470; chr. 10), TIR-NBS-LRR for mungbean yellow mosaic India virus (MYMIV) resistance (VRADI09G06940; chr. 9), E3 ubiquitin-protein ligase RIE1 for SPAD value (VRADI07G28100; chr. 11), WRKY family transcription factor for leaf area (VRADI03G06560; chr. 3), and LOB domain-containing protein 21 for trichomes (VRADI06G04290; chr. 6). In-silico validation of candidate genes was done through digital gene expression analysis using Arabidopsis orthologous (compared with Vigna radiata genome). The findings provided valuable insight for marker-assisted breeding aiming for the development of YMD-resistant and early-maturing mungbean varieties.

Natural alleles of Mediator subunit genes modulate plant height in chickpea

Plant height (PH) is an important plant architectural trait targeted during Green Revolution to enhance crop yields. Identification of genes and natural alleles governing plant height without compromising agronomic performance can fill the lacuna of knowledge connecting ideal plant architecture with maximum achievable yield in chickpea. Through coherent strategy involving genome‐wide association study, QTL/fine mapping, map‐based cloning, molecular haplotyping, and downstream functional genomics, the current study identified two Mediator subunit genes namely, CaMED23 and CaMED5b and their derived natural alleles/haplotypes underlying the major QTLs and trans‐acting eQTLs regulating plant height in chickpea. Differential accumulation of haplotype‐specific transcripts of these two Mediator genes in corresponding haplotype‐introgressed near‐isogenic lines (NILs) correlates negatively with the plant height trait. Quantitative as well as qualitative estimation based on histology, scanning electron microscopy, and histochemical assay unraveled the reduced lengths and cell sizes of internodes along with compromised lignin levels in dwarf/semi‐dwarf chickpea NILs introgressed with superior CaMED23 and CaMED5b gene haplotypes. This observation, supported by global transcriptome profiling‐based diminished expression of various phenylpropanoid pathway genes upstream of lignin biosynthesis in dwarf/semi‐dwarf NILs, essentially links plant height with lignin accumulation. The identified molecular signatures in the Mediator subunit genes can be efficiently utilized to develop desirable dwarf/semi‐dwarf‐type chickpea cultivars without affecting their yield per plant via modulating lignin/phenylpropanoid biosynthesis.

Narrow genetic diversity in germplasm from the Guinean and Sudano-Guinean zones in Benin indicates the need to broaden the genetic base of sweet fig banana (Musa acuminata cv Sotoumon)

Sweet fig (M. acuminata cv. Sotoumon) is an economically important dessert banana in Benin, with high nutritional, medicinal, and cultural values. Nevertheless, its productivity and yield are threatened by biotic and abiotic stresses. Relevant knowledge of the genetic diversity of this economically important crop is essential for germplasm conservation and the development of breeding programs. However, very little is known about the genetic makeup of this cultivar in Benin. To advance the understanding of genetic diversity in sweet fig banana germplasm, a Genotype-By-Sequencing (GBS) was performed on a panel of 273 accessions collected in different phytogeographical zones of Benin. GBS generated 8,457 quality SNPs, of which 1992 were used for analysis after filtering. The results revealed a low diversity in the studied germplasm (He = 0.0162). Genetic differentiation was overall very low in the collection as suggested by the negative differentiation index (Fstg = -0.003). The Analysis of Molecular Variance (AMOVA) indicated that the variation between accessions within populations accounted for 83.8% of the total variation observed (P < 0.001). The analysis of population structure and neighbor-joining tree partitioned the germplasm into three clusters out of which a predominant major one contained 98.1% of all accessions. These findings demonstrate that current sweet fig banana genotypes shared a common genetic background, which made them vulnerable to biotic and abiotic stress. Therefore, broadening the genetic base of the crop while maintaining its quality attributes and improving yield performance is of paramount importance. Moreover, the large genetic group constitutes an asset for future genomic selection studies in the crop and can guide the profiling of its conservation strategies.

Genome-wide association studies for phenological and agronomic traits in mungbean (Vigna radiata L. Wilczek)

Mungbean (Vigna radiata L. Wilczek) is one of the important warm-season food legumes, contributing substantially to nutritional security and environmental sustainability. The genetic complexity of yield-associated agronomic traits in mungbean is not well understood. To dissect the genetic basis of phenological and agronomic traits, we evaluated 153 diverse mungbean genotypes for two phenological (days to heading and days to maturity) and eight agronomic traits (leaf nitrogen status using SPAD, plant height, number of primary branches, pod length, number of pods per plant, seeds per pod, 100-seed weight, and yield per plant) under two environmental conditions. A wide array of phenotypic variability was apparent among the studied genotypes for all the studied traits. The broad sense of heritability of traits ranged from 0.31 to 0.95 and 0.21 to 0.94 at the Delhi and Ludhiana locations, respectively. A total of 55,634 genome-wide single nucleotide polymorphisms (SNPs) were obtained by the genotyping-by-sequencing method, of which 15,926 SNPs were retained for genome-wide association studies (GWAS). GWAS with Bayesian information and linkage-disequilibrium iteratively nested keyway (BLINK) model identified 50 SNPs significantly associated with phenological and agronomic traits. In total, 12 SNPs were found to be significantly associated with phenological traits across environments, explaining 7%-18.5% of phenotypic variability, and 38 SNPs were significantly associated with agronomic traits, explaining 4.7%-27.6% of the phenotypic variability. The maximum number of SNPs (15) were located on chromosome 1, followed by seven SNPs each on chromosomes 2 and 8. The BLAST search identified 19 putative candidate genes that were involved in light signaling, nitrogen responses, phosphorus (P) transport and remobilization, photosynthesis, respiration, metabolic pathways, and regulating growth and development. Digital expression analysis of 19 genes revealed significantly higher expression of 12 genes, viz. VRADI01G08170, VRADI11G09170, VRADI02G00450, VRADI01G00700, VRADI07G14240, VRADI03G06030, VRADI02G14230, VRADI08G01540, VRADI09G02590, VRADI08G00110, VRADI02G14240, and VRADI02G00430 in the roots, cotyledons, seeds, leaves, shoot apical meristems, and flowers. The identified SNPs and putative candidate genes provide valuable genetic information for fostering genomic studies and marker-assisted breeding programs that improve yield and agronomic traits in mungbean.

A superior gene allele involved in abscisic acid signaling enhances drought tolerance and yield in chickpea

Identifying potential molecular tags for drought tolerance is essential for achieving higher crop productivity under drought stress. We employed an integrated genomics-assisted breeding and functional genomics strategy involving association mapping, fine mapping, map-based cloning, molecular haplotyping and transcript profiling in the introgression lines (ILs)- and near isogenic lines (NILs)-based association panel and mapping population of chickpea (Cicer arietinum). This combinatorial approach delineated a bHLH (basic helix-loop-helix) transcription factor, CabHLH10 (Cicer arietinum bHLH10) underlying a major QTL, along with its derived natural alleles/haplotypes governing yield traits under drought stress in chickpea. CabHLH10 binds to a cis-regulatory G-box promoter element to modulate the expression of RD22 (responsive to desiccation 22), a drought/abscisic acid (ABA)-responsive gene (via a trans-expression QTL), and two strong yield-enhancement photosynthetic efficiency (PE) genes. This, in turn, upregulates other downstream drought-responsive and ABA signaling genes, as well as yield-enhancing PE genes, thus increasing plant adaptation to drought with reduced yield penalty. We showed that a superior allele of CabHLH10 introgressed into the NILs improved root and shoot biomass and PE, thereby enhancing yield and productivity during drought without compromising agronomic performance. Furthermore, overexpression of CabHLH10 in chickpea and Arabidopsis (Arabidopsis thaliana) conferred enhanced drought tolerance by improving root and shoot agro-morphological traits. These findings facilitate translational genomics for crop improvement and the development of genetically tailored, climate-resilient, high-yielding chickpea cultivars.

Mining legume germplasm for genetic gains: An Indian perspective

Legumes play a significant role in food and nutritional security and contribute to environmental sustainability. Although legumes are highly beneficial crops, it has not yet been possible to enhance their yield and production to a satisfactory level. Amid a rising population and low yield levels, per capita average legume consumption in India has fallen by 71% over the last 50 years, and this has led to protein-related malnutrition in a large segment of the Indian population, especially women and children. Several factors have hindered attempts to achieve yield enhancement in grain legumes, including biotic and abiotic pressures, a lack of good ideotypes, less amenability to mechanization, poorer responsiveness to fertilizer input, and a poor genetic base. Therefore, there is a need to mine the approximately 0.4 million ex situ collections of legumes that are being conserved in gene banks globally for identification of ideal donors for various traits. The Indian National Gene Bank conserves over 63,000 accessions of legumes belonging to 61 species. Recent initiatives have been undertaken in consortia mode with the aim of unlocking the genetic potential of ex situ collections and conducting large-scale germplasm characterization and evaluation analyses. We assume that large-scale phenotyping integrated with omics-based science will aid the identification of target traits and their use to enhance genetic gains. Additionally, in cases where the genetic base of major legumes is narrow, wild relatives have been evaluated, and these are being exploited through pre-breeding. Thus far, >200 accessions of various legumes have been registered as unique donors for various traits of interest.

Analysis of Genetic Diversity and Population Structure of Cowpea (Vigna unguiculata (L.) Walp) Genotypes Using Single Nucleotide Polymorphism Markers

Cowpea (Vigna unguiculata (L.) Walp) is an important legume crop with immense potential for nutritional and food security, income generation, and livestock feed in Sub-Saharan Africa. The crop is highly tolerant to heat and drought stresses which makes it an extremely important crop for improving resilience in crop production in the face of climate change. This study was carried out to assess the genetic diversity and population structure of 90 cowpea accessions using single nucleotide polymorphism (SNP) markers. Out of 11,940 SNPs used, 5864 SNPs were polymorphic and maintained for genome diversity analysis. Polymorphic information content (PIC) values ranged from 0.22 to 0.32 with a mean value of 0.27. The model-based Bayesian STRUCTURE analysis classified 90 cowpea accessions into four subpopulations at K = 4, while the distance-based cluster analysis grouped the accessions into three distinct clusters. The analysis of molecular variance (AMOVA) revealed that 59% and 69% of the total molecular variation was attributed to among individual variation for model-based and distance-based populations, respectively, and 18% was attributed to within individual variations. Furthermore, the low heterozygosity among cowpea accessions and the high inbreeding coefficient observed in this study suggests that the accessions reached an acceptable level of homozygosity. This study would serve as a reference for future selection and breeding programs of cowpea with desirable traits and systematic conservation of these plant genetic resources.

RP, R. L. R. Population structure and genetic diversity characterization of soybean for seed longevity

Seed longevity is an important trait in the context of germplasm conservation and economics of seed production. The identification of populations with high level of genetic variability for seed longevity and associated traits will become a valuable resource for superior alleles for seed longevity. In this study, Genotyping-by-sequencing (GBS)-single nucleotide polymorphism (SNP) approach, simple sequence repeats (SSR) markers and agro-morphological traits have been explored to investigate the diversity and population structure of assembled 96 genotypes. The GBS technique performed on 96 genotypes of soybean (Glycine max (L.) Merrill) resulted in 37,897 SNPs on sequences aligned to the reference genome sequence. The average genome coverage was 6.81X with a mapping rate of 99.56% covering the entire genome. Totally, 29,955 high quality SNPs were identified after stringent filtering and most of them were detected in non-coding regions. The 96 genotypes were phenotyped for eight quantitative and ten qualitative traits by growing in field by following augmented design. The STRUCTURE (Bayesian-model based algorithm), UPGMA (Un-weighed Pair Group Method with Arithmetic mean) and principal component analysis (PCA) approaches using SSR, SNP as well as quantitative and qualitative traits revealed population structure and diversity in assembled population. The Bayesian-model based STRUCTURE using SNP markers could effectively identify clusters with higher seed longevity associated with seed coat colour and size which were subsequently validated by UPGMA and PCA based on SSR and agro-morphological traits. The results of STRUCTURE, PCA and UPGMA cluster analysis showed high degree of similarity and provided complementary data that helped to identify genotypes with higher longevity. Six black colour genotypes, viz., Local black soybean, Kalitur, ACC Nos. 39, 109, 101 and 37 showed higher seed longevity during accelerated ageing. Higher coefficient of variability observed for plant height, number of pods per plant, seed yield per plant, 100 seed weight and seed longevity confirms the diversity in assembled population and its suitability for quantitative trait loci (QTL) mapping.

Genome-Wide Association Analysis Reveals Trait-Linked Markers for Grain Nutrient and Agronomic Traits in Diverse Set of Chickpea Germplasm

Chickpea is an inexpensive source of protein, minerals, and vitamins to the poor people living in arid and semi-arid regions of Southern Asia and Sub-Saharan Africa. New chickpea cultivars with enhanced levels of protein, Fe and Zn content are a medium-term strategy for supplying essential nutrients for human health and reducing malnutrition. In the current study, a chickpea reference set of 280 accessions, including landraces, breeding lines, and advanced cultivars, was evaluated for grain protein, Fe, Zn content and agronomic traits over two seasons. Using a mid-density 5k SNP array, 4603 highly informative SNPs distributed across the chickpea genome were used for GWAS analysis. Population structure analysis revealed three subpopulations (K = 3). Linkage disequilibrium (LD) was extensive, and LD decay was relatively low. A total of 20 and 46 marker-trait associations (MTAs) were identified for grain nutrient and agronomic traits, respectively, using FarmCPU and BLINK models. Of which seven SNPs for grain protein, twelve for Fe, and one for Zn content were distributed on chromosomes 1, 4, 6, and 7. The marker S4_4477846 on chr4 was found to be co-associated with grain protein over seasons. The markers S1_11613376 and S1_2772537 co-associated with grain Fe content under NSII and pooled seasons and S7_9379786 marker under NSI and pooled seasons. The markers S4_31996956 co-associated with grain Fe and days to maturity. SNP annotation of associated markers were found to be related to gene functions of metal ion binding, transporters, protein kinases, transcription factors, and many more functions involved in plant metabolism along with Fe and protein homeostasis. The identified significant MTAs has potential use in marker-assisted selection for developing nutrient-rich chickpea cultivars after validation in the breeding populations.

Exploring Chickpea Germplasm Diversity for Broadening the Genetic Base Utilizing Genomic Resourses

Legume crops provide significant nutrition to humans as a source of protein, omega-3 fatty acids as well as specific macro and micronutrients. Additionally, legumes improve the cropping environment by replenishing the soil nitrogen content. Chickpeas are the second most significant staple legume food crop worldwide behind dry bean which contains 17%–24% protein, 41%–51% carbohydrate, and other important essential minerals, vitamins, dietary fiber, folate, β-carotene, anti-oxidants, micronutrients (phosphorus, calcium, magnesium, iron, and zinc) as well as linoleic and oleic unsaturated fatty acids. Despite these advantages, legumes are far behind cereals in terms of genetic improvement mainly due to far less effort, the bottlenecks of the narrow genetic base, and several biotic and abiotic factors in the scenario of changing climatic conditions. Measures are now called for beyond conventional breeding practices to strategically broadening of narrow genetic base utilizing chickpea wild relatives and improvement of cultivars through advanced breeding approaches with a focus on high yield productivity, biotic and abiotic stresses including climate resilience, and enhanced nutritional values. Desirable donors having such multiple traits have been identified using core and mini core collections from the cultivated gene pool and wild relatives of Chickpea. Several methods have been developed to address cross-species fertilization obstacles and to aid in inter-specific hybridization and introgression of the target gene sequences from wild Cicer species. Additionally, recent advances in “Omics” sciences along with high-throughput and precise phenotyping tools have made it easier to identify genes that regulate traits of interest. Next-generation sequencing technologies, whole-genome sequencing, transcriptomics, and differential genes expression profiling along with a plethora of novel techniques like single nucleotide polymorphism exploiting high-density genotyping by sequencing assays, simple sequence repeat markers, diversity array technology platform, and whole-genome re-sequencing technique led to the identification and development of QTLs and high-density trait mapping of the global chickpea germplasm. These altogether have helped in broadening the narrow genetic base of chickpeas.

Discovery of DNA polymorphisms via whole genome resequencing and their functional relevance in salinity stress response in chickpea

Salinity stress is one of the major constraints for plant growth and yield. The salinity stress response of different genotypes of crop plants may largely be governed by DNA polymorphisms. To determine the molecular genetic factors involved in salinity stress tolerance in chickpea, we performed a whole genome resequencing data analysis of three each of salinity-sensitive and salinity-tolerant genotypes. A total of 6173 single nucleotide polymorphisms and 920 insertions and deletions differentiating the chickpea genotypes with contrasting salinity stress responses were identified. Gene ontology analysis revealed the enrichment of functional terms related to stress response and development among the genes harboring DNA polymorphisms in their promoter and/or coding regions. DNA polymorphisms located within the cis-regulatory motifs of the quantitative trait loci (QTL)-associated and abiotic stress related genes were identified, which may influence salinity stress response via modulating binding affinity of the transcription factors. Several genes including QTL-associated and abiotic stress response related genes harboring DNA polymorphisms exhibited differential expression in response to salinity stress especially at the reproductive stage of development in the salinity-tolerant genotype. Furthermore, effects of non-synonymous DNA polymorphisms on mutational sensitivity and structural integrity of the encoded proteins by the candidate QTL-associated and abiotic stress response related genes were revealed. The results suggest that DNA polymorphisms may determine salinity stress response via influencing differential gene expression in genotype and/or stage-dependent manner. Altogether, we provide a high-quality set of DNA polymorphisms and candidate genes that may govern salinity stress tolerance in chickpea.

Genome resequencing reveals DNA polymorphisms associated with seed size/weight determination in chickpea

Crop productivity in legumes is determined by number and size/weight of seeds. To understand the genetic basis of seed size/weight in chickpea, we performed genome resequencing of 13 small- and 5 large-seeded genotypes using Illumina platform. Single nucleotide polymorphisms (SNPs) and insertions/deletions (InDels) differentiating small- and large-seeded genotypes were identified. A total of 17,902 SNPs and 2594 InDels located in promoter and/or coding regions that may contribute to seed size/weight were detected. Of these, 266 SNPs showed significant association with seed size/weight trait. Twenty-three genes including those involved in cell growth/division, encoding transcription factors and located within QTLs associated with seed size/weight harbored SNPs within transcription factor binding motif(s) and/or coding region. The non-synonymous SNPs were found to affect the mutational sensitivity and stability of the encoded proteins. Overall, we provided a high-quality SNP map for large-scale genotyping applications and identified candidate genes that determine seed size/weight in chickpea.

Genome-Wide Association Analysis for Phosphorus Use Efficiency Traits in Mungbean (Vigna radiata L. Wilczek) Using Genotyping by Sequencing Approach

Mungbean (Vigna radiata L. Wilczek) is an annual grain legume crop affected by low availability of phosphorus. Phosphorus deficiency mainly affects the growth and development of plants along with changes in root morphology and increase in root-to-shoot ratio. Deciphering the genetic basis of phosphorus use efficiency (PUE) traits can benefit our understanding of mungbean tolerance to low-phosphorus condition. To address this issue, 144 diverse mungbean genotypes were evaluated for 12 PUE traits under hydroponics with optimum- and low-phosphorus levels. The broad sense heritability of traits ranged from 0.63 to 0.92 and 0.58 to 0.92 under optimum- and low-phosphorus conditions, respectively. This study, reports for the first time such a large number of genome wide Single nucleotide polymorphisms (SNPs) (76,160) in mungbean. Further, genome wide association study was conducted using 55,634 SNPs obtained by genotyping-by-sequencing method. The results indicated that total 136 SNPs shared by both GLM and MLM models were associated with tested PUE traits under different phosphorus regimes. We have identified SNPs with highest p value (-log10(p)) for some traits like, TLA and RDW with p value (-log10(p)) of more than 6.0 at LP/OP and OP condition. We have identified nine SNPs (three for TLA and six for RDW trait) which was found to be present in chromosomes 8, 4, and 7. One SNP present in Vradi07g06230 gene contains zinc finger CCCH domain. In total, 71 protein coding genes were identified, of which 13 genes were found to be putative candidate genes controlling PUE by regulating nutrient uptake and root architectural development pathways in mungbean. Moreover, we identified three potential candidate genes VRADI11G08340, VRADI01G05520, and VRADI04G10750 with missense SNPs in coding sequence region, which results in significant variation in protein structure at tertiary level. The identified SNPs and candidate genes provide the essential information for genetic studies and marker-assisted breeding program for improving low-phosphorus tolerance in mungbean.

Genome-wide association study for phosphate deficiency responsive root hair elongation in chickpea

Root hairs (RHs) are single-celled elongated epidermal cells and play a vital role in nutrient absorption, particularly for immobile minerals like phosphorus (P). As an adaptive response to P deficiency, an increase in RH length enhances root-soil contact and absorptive area for P absorption. Genetic variations have been reported for RH length and its response to P deficiency in plants. However, only a few association studies have been conducted to identify genes and genetic loci associated with RH length. Here, we screened desi chickpea accessions for RH length and its plasticity under P deficiency. Further, the genome-wide association study (GWAS) was conducted to identify the genetic loci associated with RH length in P deficient and sufficient conditions. Although high variability was observed in terms of RH length in diverse genotypes, majority of the accessions showed typical response of increase in RH length in low P. Genome-wide association mapping identified many SNPs with significant associations with RH length in P-sufficient and P-deficient conditions. A few candidate genes for RH length in P deficient (SIZ1-like and HAD superfamily protein) and sufficient (RSL2-like and SMAP1-like) conditions were identified which have known roles in RH development and P deficiency response or both. Highly associated loci and candidate genes identified in this study would be useful for genomic-assisted breeding to develop P-efficient chickpea.

A high-density SNP-based linkage map using genotyping-by-sequencing and its utilization for improved genome assembly of chickpea (Cicer arietinum L.)

Genotyping-by-sequencing (GBS) allows rapid identification of markers for use in development of linkage maps, which expedite efficient breeding programs. In the present study, we have utilized GBS approach to identify and genotype single-nucleotide polymorphism (SNP) markers in an inter-specific RIL population of Cicer arietinum L. X C. reticulatum. A total of 141,639 raw SNPs were identified using the TASSEL-GBS pipeline. After stringent filtering, 8208 candidate SNPs were identified of which ~ 37% were localized in the intragenic regions followed by genic regions (~ 30%) and intergenic regions (~ 27%). We then utilized 6920 stringent selected SNPs from present study and 6714 SNPs and microsatellite markers available from previous studies for construction of linkage map. The resulting high-density linkage map comprising of eight linkage groups contained 13,590 markers which spanned 1299.14 cM of map length with an average marker density of 0.095 cM. Further, the derived linkage map was used to improve the available assembly of desi chickpea genome by anchoring 443 previously unplaced scaffolds onto eight linkage groups. The present efforts have refined anchoring of the desi chickpea genome assembly to 55.57% of the ~ 520 Mb of assembled desi genome. To the best of our knowledge, the linkage map generated in the present study represents one of the most dense linkage map developed for the crop till date. It will serve as a valuable resource for fine mapping and positional cloning of important quantitative trait loci (QTLs) associated with agronomical traits and also for anchoring and ordering of future genome sequence assemblies.

DNA methylation reprogramming during seed development and its functional relevance in seed size/weight determination in chickpea

Seed development is orchestrated via complex gene regulatory networks and pathways. Epigenetic factors may also govern seed development and seed size/weight. Here, we analyzed DNA methylation in a large-seeded chickpea cultivar (JGK 3) during seed development stages. Progressive gain of CHH context DNA methylation in transposable elements (TEs) and higher frequency of small RNAs in hypermethylated TEs during seed development suggested a role of the RNA-dependent DNA methylation pathway. Frequency of intragenic TEs was higher in CHH context differentially methylated region (DMR) associated differentially expressed genes (DEGs). CG context hyper/hypomethylation within the gene body was observed for most of DMR-associated DEGs in JGK 3 as compared to small-seeded chickpea cultivar (Himchana 1). We identified candidate genes involved in seed size/weight determination exhibiting CG context hypermethylation within the gene body and higher expression in JGK 3. This study provides insights into the role of DNA methylation in seed development and seed size/weight determination in chickpea.

Rajkumar et al. report progressive gain of CHH context DNA methylation in transposable elements during seed development in chickpea, of which hypermethylation is associated with small RNAs. The candidate genes that determine seed size/weight in chickpea show CG context hypermethylation in the gene body and higher expression in large-seeded cultivar.

Population structure and genetic diversity in red clover (Trifolium pratense L.) germplasm

Red clover (Trifolium pratense L.) is a highly adaptable forage crop for temperate livestock agriculture. Genetic variation can be identified, via molecular techniques, and used to assess diversity among populations that may otherwise be indistinguishable. Here we have used genotyping by sequencing (GBS) to determine the genetic variation and population structure in red clover natural populations from Europe and Asia, and varieties or synthetic populations. Cluster analysis differentiated the collection into four large regional groups: Asia, Iberia, UK, and Central Europe. The five varieties clustered with the geographical area from which they were derived. Two methods (BayeScan and Samβada) were used to search for outlier loci indicating signatures of selection. A total of 60 loci were identified by both methods, but no specific genomic region was highlighted. The rate of decay in linkage disequilibrium was fast, and no significant evidence of any bottlenecks was found. Phenotypic analysis showed that a more prostrate and spreading growth habit was predominantly found among populations from Iberia and the UK. A genome wide association study identified a single nucleotide polymorphism (SNP) located in a homologue of the VEG2 gene from pea, associated with flowering time. The identification of genetic variation within the natural populations is likely to be useful for enhancing the breeding of red clover in the future.

Genome-wide SNP identification in Prunus rootstocks germplasm collections using Genotyping-by-Sequencing: phylogenetic analysis, distribution of SNPs and prediction of their effect on gene function

Genotyping-by-Sequencing (GBS) was applied in a set of 53 diploid Prunus rootstocks and five scion cultivars from three subgenera (Amygdalus, Prunus and Cerasus) for genome-wide SNP identification and to assess genetic diversity of both Chilean and Spanish germplasm collections. A group of 45,382 high quality SNPs (MAF >0.05; missing data <5%) were selected for analysis of this group of 58 accessions. These SNPs were distributed in genic and intergenic regions in the eight pseudomolecules of the peach genome (Peach v2.0), with an average of 53% located in exonic regions. The genetic diversity detected among the studied accessions divided them in three groups, which are in agreement with their current taxonomic classification. SNPs were classified based on their putative effect on annotated genes and KOG analysis was carried out to provide a deeper understanding of the function of 119 genes affected by high-impact SNPs. Results demonstrate the high utility for Prunus rootstocks identification and studies of diversity in Prunus species. Also, given the high number of SNPs identified in exonic regions, this strategy represents an important tool for finding candidate genes underlying traits of interest and potential functional markers for use in marker-assisted selection.

Genome-Wide SNP discovery and genomic characterization in avocado (Persea americana Mill.)

Modern crop breeding is based on the use of genetically and phenotypically diverse plant material and, consequently, a proper understanding of population structure and genetic diversity is essential for the effective development of breeding programs. An example is avocado, a woody perennial fruit crop native to Mesoamerica with an increasing popularity worldwide. Despite its commercial success, there are important gaps in the molecular tools available to support on-going avocado breeding programs. In order to fill this gap, in this study, an avocado 'Hass' draft assembly was developed and used as reference to study 71 avocado accessions which represent the three traditionally recognized avocado horticultural races or subspecies (Mexican, Guatemalan and West Indian). An average of 5.72 M reads per individual and a total of 7,108 single nucleotide polymorphism (SNP) markers were produced for the 71 accessions analyzed. These molecular markers were used in a study of genetic diversity and population structure. The results broadly separate the accessions studied according to their botanical race in four main groups: Mexican, Guatemalan, West Indian and an additional group of Guatemalan × Mexican hybrids. The high number of SNP markers developed in this study will be a useful genomic resource for the avocado community.

Whole Genome Diversity, Population Structure, and Linkage Disequilibrium Analysis of Chickpea (Cicer arietinum L.) Genotypes Using Genome-Wide DArTseq-Based SNP Markers

Characterization of genetic diversity, population structure, and linkage disequilibrium is a prerequisite for proper management of breeding programs and conservation of genetic resources. In this study, 186 chickpea genotypes, including advanced "Kabuli" breeding lines and Iranian landrace "Desi" chickpea genotypes, were genotyped using DArTseq-Based single nucleotide polymorphism (SNP) markers. Out of 3339 SNPs, 1152 markers with known chromosomal position were selected for genome diversity analysis. The number of mapped SNP markers varied from 52 (LG8) to 378 (LG4), with an average of 144 SNPs per linkage group. The chromosome size that was covered by SNPs varied from 16,236.36 kbp (LG8) to 67,923.99 kbp (LG5), while LG4 showed a higher number of SNPs, with an average of 6.56 SNPs per Mbp. Polymorphism information content (PIC) value of SNP markers ranged from 0.05 to 0.50, with an average of 0.32, while the markers on LG4, LG6, and LG8 showed higher mean PIC value than average. Unweighted neighbor joining cluster analysis and Bayesian-based model population structure grouped chickpea genotypes into four distinct clusters. Principal component analysis (PCoA) and discriminant analysis of principal component (DAPC) results were consistent with that of the cluster and population structure analysis. Linkage disequilibrium (LD) was extensive and LD decay in chickpea germplasm was relatively low. A few markers showed r2 ≥ 0.8, while 2961 pairs of markers showed complete LD (r2 = 1), and a huge LD block was observed on LG4. High genetic diversity and low kinship value between pairs of genotypes suggest the presence of a high genetic diversity among the studied chickpea genotypes. This study also demonstrates the efficiency of DArTseq-based SNP genotyping for large-scale genome analysis in chickpea. The genotypic markers provided in this study are useful for various association mapping studies when combined with phenotypic data of different traits, such as seed yield, abiotic, and biotic stresses, and therefore can be efficiently used in breeding programs to improve chickpea.

Multi-trait analysis of domestication genes in Cicer arietinum - Cicer reticulatum hybrids with a multidimensional approach: Modeling wide crosses for crop improvement

Domestication and subsequent breeding have eroded genetic diversity in the modern chickpea crop by ˜100-fold. Corresponding reductions to trait variation create the need, and an opportunity, to identify and harness the genetic capacity of wild species for crop improvement. Here we analyze trait segregation in a series of wild x cultivated hybrid populations to delineate the genetic underpinnings of domestication traits. Two species of wild chickpea, C. reticulatum and C. echinospermum, were crossed with the elite, early flowering C. arietinum cultivar ICCV96029. KASP genotyping of F2 parents with an FT-linked molecular marker enabled selection of 284 F3 families with reduced phenological variation: 255 F3 families of C. arietinum x reticulatum (AR) derived from 17 diverse wild parents and 29 F3 families of C. arietinum x echinospermum (AE) from 3 wild parents. The combined 284 lineages were genotyped using a genotyping-by-sequencing strategy and phenotyped for agronomic traits. 50 QTLs in 11 traits were detected from AR and 35 QTLs in 10 traits from the combined data. Using hierarchical clustering to assign traits to six correlated groups and mixed model based multi-trait mapping, four pleiotropic loci were identified. Bayesian analysis further identified four inter-trait relationships controlling the duration of vegetative growth and seed maturation, for which the underlying pleiotropic genes were mapped. A random forest approach was used to explore the most extreme trait differences between AR and AE progenies, identifying traits most characteristic of wild species origin. Knowledge of the genomic basis of traits that segregate in wild-cultivated hybrid populations will facilitate chickpea improvement by linking genetic and phenotypic variation in a quantitative genetic framework.

Genetic diversity, linkage disequilibrium, and population structure analysis of the tea plant (Camellia sinensis) from an origin center, Guizhou plateau, using genome-wide SNPs developed by genotyping-by-sequencing

BACKGROUND

To efficiently protect and exploit germplasm resources for marker development and breeding purposes, we must accurately depict the features of the tea populations. This study focuses on the Camellia sinensis (C. sinensis) population and aims to (i) identify single nucleotide polymorphisms (SNPs) on the genome level, (ii) investigate the genetic diversity and population structure, and (iii) characterize the linkage disequilibrium (LD) pattern to facilitate next genome-wide association mapping and marker-assisted selection.

RESULTS

We collected 415 tea accessions from the Origin Center and analyzed the genetic diversity, population structure and LD pattern using the genotyping-by-sequencing (GBS) approach. A total of 79,016 high-quality SNPs were identified; the polymorphism information content (PIC) and genetic diversity (GD) based on these SNPs showed a higher level of genetic diversity in cultivated type than in wild type. The 415 accessions were clustered into three groups by STRUCTURE software and confirmed using principal component analyses (PCA)-wild type, cultivated type, and admixed wild type. However, unweighted pair group method with arithmetic mean (UPGMA) trees indicated the accessions should be grouped into more clusters. Further analyses identified four groups, the Pure Wild Type, Admixed Wild Type, ancient landraces and modern landraces using STRUCTURE, and the results were confirmed by PCA and UPGMA tree method. A higher level of genetic diversity was detected in ancient landraces and Admixed Wild Type than that in the Pure Wild Type and modern landraces. The highest differentiation was between the Pure Wild Type and modern landraces. A relatively fast LD decay with a short range (kb) was observed, and the LD decays of four inferred populations were different.

CONCLUSIONS

This study is, to our knowledge, the first population genetic analysis of tea germplasm from the Origin Center, Guizhou Plateau, using GBS. The LD pattern, population structure and genetic differentiation of the tea population revealed by our study will benefit further genetic studies, germplasm protection, and breeding.

CLAVATA signaling pathway genes modulating flowering time and flower number in chickpea

A combinatorial genomic strategy delineated functionally relevant natural allele of a CLAVATA gene and its marker (haplotype)-assisted introgression led to development of the early-flowering chickpea cultivars with high flower number and enhanced yield/productivity. Unraveling the genetic components involved in CLAVATA (CLV) signaling is crucial for modulating important shoot apical meristem (SAM) characteristics and ultimately regulating diverse SAM-regulated agromorphological traits in crop plants. A genome-wide scan identified 142 CLV1-, 28 CLV2- and 6 CLV3-like genes, and their comprehensive genomic constitution and phylogenetic relationships were deciphered in chickpea. The QTL/fine mapping and map-based cloning integrated with high-resolution association analysis identified SNP loci from CaCLV3_01 gene within a major CaqDTF1.1/CaqFN1.1 QTL associated with DTF (days to 50% flowering) and FN (flower number) traits in chickpea, which was further ascertained by quantitative expression profiling. Molecular haplotyping of CaCLV3_01 gene, expressed specifically in SAM, constituted two major haplotypes that differentiated the early-DTF and high-FN chickpea accessions from late-DTF and low-FN. Enhanced accumulation of transcripts of superior CaCLV3_01 gene haplotype and known flowering promoting genes was observed in the corresponding haplotype-introgressed early-DTF and high-FN near-isogenic lines (NILs) with narrow SAM width. The superior haplotype-introgressed NILs exhibited early-flowering, high-FN and enhanced seed yield/productivity without compromising agronomic performance. These delineated molecular signatures can regulate DTF and FN traits through SAM proliferation and differentiation and thereby will be useful for translational genomic study to develop early-flowering cultivars with enhanced yield/productivity.

Towards Exploitation of Adaptive Traits for Climate-Resilient Smart Pulses

Pulses are the main source of protein and minerals in the vegetarian diet. These are primarily cultivated on marginal lands with few inputs in several resource-poor countries of the world, including several in South Asia. Their cultivation in resource-scarce conditions exposes them to various abiotic and biotic stresses, leading to significant yield losses. Furthermore, climate change due to global warming has increased their vulnerability to emerging new insect pests and abiotic stresses that can become even more serious in the coming years. The changing climate scenario has made it more challenging to breed and develop climate-resilient smart pulses. Although pulses are climate smart, as they simultaneously adapt to and mitigate the effects of climate change, their narrow genetic diversity has always been a major constraint to their improvement for adaptability. However, existing genetic diversity still provides opportunities to exploit novel attributes for developing climate-resilient cultivars. The mining and exploitation of adaptive traits imparting tolerance/resistance to climate-smart pulses can be accelerated further by using cutting-edge approaches of biotechnology such as transgenics, genome editing, and epigenetics. This review discusses various classical and molecular approaches and strategies to exploit adaptive traits for breeding climate-smart pulses.

Genome-wide cis-regulatory signatures for modulation of agronomic traits as exemplified by drought yield index (DYI) in chickpea

Developing functional molecular tags from the cis-regulatory sequence components of genes is vital for their deployment in efficient genetic dissection of complex quantitative traits in crop plants including chickpea. The current study identified 431,194 conserved non-coding SNP (CNSNP) from the cis-regulatory element regions of genes which were annotated on a chickpea genome. These genome-wide CNSNP marker resources are made publicly accessible through a user-friendly web-database ( http://www.cnsnpcicarbase.com ). The CNSNP-based quantitative trait loci (QTL) and expression QTL (eQTL) mapping and genome-wide association study (GWAS) were further integrated with global gene expression landscapes, molecular haplotyping, and DNA-protein interaction study in the association panel and recombinant inbred lines (RIL) mapping population to decode complex genetic architecture of one of the vital seed yield trait under drought stress, drought yield index (DYI), in chickpea. This delineated two constituted natural haplotypes and alleles from a histone H3 protein-coding gene and its transcriptional regulator NAC transcription factor (TF) harboring the major QTLs and trans-acting eQTL governing DYI in chickpea. The effect of CNSNPs in TF-binding cis-element of a histone H3 gene in altering the binding affinity and transcriptional activity of NAC TF based on chromatin immunoprecipitation-quantitative PCR (ChIP-qPCR) assay was evident. The CNSNP-led promising molecular tags scanned will essentially have functional significance to decode transcriptional gene regulatory function and thus can drive translational genomic analysis in chickpea.

Transcriptional signatures modulating shoot apical meristem morphometric and plant architectural traits enhance yield and productivity in chickpea

Plant height (PH) and plant width (PW), two of the major plant architectural traits determining the yield and productivity of a crop, are defined by diverse morphometric characteristics of the shoot apical meristem (SAM). The identification of potential molecular tags from a single gene that simultaneously modulates these plant/SAM architectural traits is therefore prerequisite to achieve enhanced yield and productivity in crop plants, including chickpea. Large-scale multienvironment phenotyping of the association panel and mapping population have ascertained the efficacy of three vital SAM morphometric trait parameters, SAM width, SAM height and SAM area, as key indicators to unravel the genetic basis of the wide PW and PH trait variations observed in desi chickpea. This study integrated a genome-wide association study (GWAS); quantitative trait locus (QTL)/fine-mapping and map-based cloning with molecular haplotyping; transcript profiling; and protein-DNA interaction assays for the dissection of plant architectural traits in chickpea. These exertions delineated natural alleles and superior haplotypes from a CabHLH121 transcription factor (TF) gene within the major QTL governing PW, PH and SAM morphometric traits. A genome-wide protein-DNA interaction assay assured the direct binding of a known stem cell master regulator, CaWUS, to the WOX-homeodomain TF binding sites of a CabHLH121 gene and its constituted haplotypes. The differential expression of CaWUS and transcriptional regulation of its target CabHLH121 gene/haplotypes were apparent, suggesting their collective role in altering SAM morphometric characteristics and plant architectural traits in the contrasting near isogenic lines (NILs). The NILs introgressed with a superior haplotype of a CabHLH121 exhibited optimal PW and desirable PH as well as enhanced yield and productivity without compromising any component of agronomic performance. These molecular signatures of the CabHLH121 TF gene have the potential to regulate both PW and PH traits through the modulation of proliferation, differentiation and maintenance of the meristematic stem cell population in the SAM; therefore, these signatures will be useful in the translational genomic study of chickpea genetic enhancement. The restructured cultivars with desirable PH (semidwarf) and PW will ensure maximal planting density in a specified cultivable field area, thereby enhancing the overall yield and productivity of chickpea. This can essentially facilitate the achievement of better remunerative outputs by farmers with rational land use, therefore ensuring global food security in the present scenario of an increasing population density and shrinking per capita land area.

ABC Transporter-Mediated Transport of Glutathione Conjugates Enhances Seed Yield and Quality in Chickpea

The identification of functionally relevant molecular tags is vital for genomics-assisted crop improvement and enhancement of seed yield, quality, and productivity in chickpea (Cicer arietinum). The simultaneous improvement of yield/productivity as well as quality traits often requires pyramiding of multiple genes, which remains a major hurdle given various associated epistatic and pleotropic effects. Unfortunately, no single gene that can improve yield/productivity along with quality and other desirable agromorphological traits is known, hampering the genetic enhancement of chickpea. Using a combinatorial genomics-assisted breeding and functional genomics strategy, this study identified natural alleles and haplotypes of an ABCC3-type transporter gene that regulates seed weight, an important domestication trait, by transcriptional regulation and modulation of the transport of glutathione conjugates in seeds of desi and kabuli chickpea. The superior allele/haplotype of this gene introgressed in desi and kabuli near-isogenic lines enhances the seed weight, yield, productivity, and multiple desirable plant architecture and seed-quality traits without compromising agronomic performance. These salient findings can expedite crop improvement endeavors and the development of nutritionally enriched high-yielding cultivars in chickpea.

Genetic dissection of photosynthetic efficiency traits for enhancing seed yield in chickpea

Understanding the genetic basis of photosynthetic efficiency (PE) contributing to enhanced seed yield per plant (SYP) is vital for genomics-assisted crop improvement of chickpea. The current study employed an integrated genomic strategy involving photosynthesis pathway gene-based association mapping, genome-wide association study, quantitative trait loci (QTL) mapping, and expression profiling. This identified 16 potential single nucleotide polymorphism loci linked to major QTLs underlying 16 candidate genes significantly associated with PE and SYP traits in chickpea. The allelic variants were tightly linked to positively interacting QTLs regulating both enhanced PE and SYP traits as exemplified by a chlorophyll A-B binding protein-coding gene. The leaf tissue-specific pronounced up-regulated expression of 16 associated genes in germplasm accessions and homozygous individuals of mapping population was evident. Such combinatorial genomic strategy coupled with gene haplotype-specific association and in silico protein-protein interaction study delineated natural alleles and superior haplotypes from a chlorophyll A-B binding (CAB) protein-coding gene and its interacting gene, Timing of CAB Expression 1 (TOC1), which appear to be most promising candidates in modulating chickpea PE and SYP traits. These functionally pertinent molecular signatures identified have efficacy to drive marker-assisted selection for developing PE-enriched cultivars with high seed yield in chickpea.

Analysis of genes encoding seed storage proteins (SSPs) in chickpea (Cicer arietinum L.) reveals co-expressing transcription factors and a seed-specific promoter

Improvement of the quality and quantity of chickpea seed protein can be greatly facilitated by an understanding of the genic organization and the genetic architecture of the genes encoding seed storage proteins (SSPs). The aim of this study was to provide a comprehensive analysis of the chickpea SSP genes, putative co-expressing transcription factors (TFs), and to identify a seed-specific SSP gene promoter. A genome-wide identification of SSP genes in chickpea led to the identification of 21 non-redundant SSP encoding genes located on 6 chromosomes. Phylogenetic analysis grouped SSP genes into 3 subgroups where members within the same clade demonstrated similar motif composition and intron-exon organization. Tandem duplications were identified to be the major contributors to the expansion of the SSP gene family in chickpea. Co-expression analysis revealed 14 TFs having expression profiles similar to the SSP genes that included members of important TF families that are known to regulate seed development. Expression analysis of SSP genes and TFs revealed significantly higher expression in late stages of seed development as well as in high seed protein content (HPC) genotypes. In silico analysis of the promoter regions of the SSP encoding genes revealed several seed-specific cis-regulatory elements such as RY repeats, ACGT motifs, CAANTG, and GCN4. A candidate promoter was analyzed for seed specificity by generating stable transgenics in Arabidopsis. Overall, this study provides a useful resource to explore the regulatory networks involved in SSP synthesis and/or accumulation for utilization in developing nutritionally improved chickpea genotypes.

Genome-wide discovery of DNA polymorphisms among chickpea cultivars with contrasting seed size/weight and their functional relevance

Seed size/weight is a major agronomic trait which determine crop productivity in legumes. To understand the genetic basis of seed size determination, we sought to identify DNA polymorphisms between two small (Himchana 1 and Pusa 362) and two large-seeded (JGK 3 and PG 0515) chickpea cultivars via whole genome resequencing. We identified a total of 75535 single nucleotide polymorphisms (SNPs), 6486 insertions and deletions (InDels), 1938 multi-nucleotide polymorphisms (MNPs) and 5025 complex variants between the two small and two large-seeded chickpea cultivars. Our analysis revealed 814, 244 and 72 seed-specific genes harboring DNA polymorphisms in promoter or non-synonymous and large-effect DNA polymorphisms, respectively. Gene ontology analysis revealed enrichment of cell growth and division related terms in these genes. Among them, at least 22 genes associated with quantitative trait loci, and those involved in cell growth and division and encoding transcription factors harbored promoter and/or large-effect/non-synonymous DNA polymorphisms. These also showed higher expression at late-embryogenesis and/or mid-maturation stages of seed development in the large-seeded cultivar, suggesting their role in seed size/weight determination in chickpea. Altogether, this study provided a valuable resource for large-scale genotyping applications and a few putative candidate genes that might play crucial role in governing seed size/weight in chickpea.

Genome-wide generation and genotyping of informative SNPs to scan molecular signatures for seed yield in chickpea

We discovered 2150 desi and 2199 kabuli accessions-derived SNPs by cultivar-wise individual assembling of sequence-reads generated through genotyping-by-sequencing of 92 chickpea accessions. Subsequent large-scale validation and genotyping of these SNPs discovered 619 desi accessions-derived (DAD) SNPs, 531 kabuli accessions-derived (KAD) SNPs, 884 multiple accessions-derived (MAD) SNPs and 1083 two accessions (desi ICC 4958 and kabuli CDC Frontier)-derived (TAD) SNPs that were mapped on eight chromosomes. These informative SNPs were annotated in coding/non-coding regulatory sequence components of genes. The MAD-SNPs were efficient to detect high intra-specific polymorphic potential and wide natural allelic diversity level including high-resolution admixed-population genetic structure and precise phylogenetic relationship among 291 desi and kabuli accessions. This signifies their effectiveness in introgression breeding and varietal improvement studies targeting useful agronomic traits of chickpea. Six trait-associated genes with SNPs including quantitative trait nucleotides (QTNs) in combination explained 27.5% phenotypic variation for seed yield per plant (SYP). A pentatricopeptide repeat (PPR) gene with a synonymous-coding SNP/QTN significantly associated with SYP trait was found most-promising in chickpea. The essential information delineated can be of immense utility in genomics-assisted breeding applications to develop high-yielding chickpea cultivars.

Consensus genetic linkage map construction and QTL mapping for plant height-related traits in linseed flax (Linum usitatissimum L.)

BACKGROUND

Flax is an important field crop that can be used for either oilseed or fiber production. Plant height and technical length are important characters for flax. For linseed flax, plants usually have a short technical length and plant height than those for fiber flax. As an important agronomical character for fiber and linseed flax, plant height is usually a selection target for breeding. However, because of limited technologies and methods available, there has been little research focused on discovering the molecular mechanism controlling plant height.

RESULTS

In this study, two related recombinant inbred line (RIL) populations developed from crosses of linseed and fiber parents were developed and phenotyped for plant height and technical length in four environments. A consensus linkage map based on two RIL populations was constructed using SNP markers generated by genotyping by sequencing (GBS) technology. A total of 4497 single nucleotide polymorphism (SNP) markers were included on 15 linkage groups with an average marker density of one marker every 2.71 cM. Quantitative trait locus (QTL) mapping analysis was performed for plant height and technical length using the two populations. A total of 19 QTLs were identified for plant height and technical length. For the MH population, eight plant height QTLs and seven technical length QTLs were identified, five of which were common QTLs for both traits. For the PH population, six plant height and three technical length QTLs were identified. By comparing the QTLs and candidate gene information in the two population, two common QTLs and three candidate genes were discovered.

CONCLUSIONS

This study provides a foundation for map-based cloning of QTLs and marker-assisted selection for plant height-related traits in linseed and fiber flax.

Consensus genetic linkage map construction and QTL mapping for plant height-related traits in linseed flax (Linum usitatissimum L.)

BACKGROUND

Flax is an important field crop that can be used for either oilseed or fiber production. Plant height and technical length are important characters for flax. For linseed flax, plants usually have a short technical length and plant height than those for fiber flax. As an important agronomical character for fiber and linseed flax, plant height is usually a selection target for breeding. However, because of limited technologies and methods available, there has been little research focused on discovering the molecular mechanism controlling plant height.

RESULTS

In this study, two related recombinant inbred line (RIL) populations developed from crosses of linseed and fiber parents were developed and phenotyped for plant height and technical length in four environments. A consensus linkage map based on two RIL populations was constructed using SNP markers generated by genotyping by sequencing (GBS) technology. A total of 4497 single nucleotide polymorphism (SNP) markers were included on 15 linkage groups with an average marker density of one marker every 2.71 cM. Quantitative trait locus (QTL) mapping analysis was performed for plant height and technical length using the two populations. A total of 19 QTLs were identified for plant height and technical length. For the MH population, eight plant height QTLs and seven technical length QTLs were identified, five of which were common QTLs for both traits. For the PH population, six plant height and three technical length QTLs were identified. By comparing the QTLs and candidate gene information in the two population, two common QTLs and three candidate genes were discovered.

CONCLUSIONS

This study provides a foundation for map-based cloning of QTLs and marker-assisted selection for plant height-related traits in linseed and fiber flax.

Genome-wide identification and characterization of InDels and SNPs in Glycine max and Glycine soja for contrasting seed permeability traits

BACKGROUND

Water permeability governed by seed coat is a major facet of seed crops, especially soybean, whose seeds lack physiological dormancy and experience rapid deterioration in seed viability under prolonged storage. Moreover, the physiological and chemical characteristics of soybean seeds are known to vary with seed coat color. Thus, to underpin the genes controlling water permeability in soybean seeds, we carried out an in-depth characterization of the associated genomic variation.

RESULTS

In the present study, we have analyzed genomic variation between cultivated soybean and its wild progenitor with implications on seed permeability, a trait related to seed storability. Whole genome resequencing of G.max and G. soja, identified SNPs and InDels which were further characterized on the basis of their genomic location and impact on gene expression. Chromosomal density distribution of the variation was assessed across the genome and genes carrying SNPs and InDels were characterized into different metabolic pathways. Seed hardiness is a complex trait that is affected by the allelic constitution of a genetic locus as well as by a tricky web of plant hormone interactions. Seven genes that hold a probable role in the determination of seed permeability were selected and their expression differences at different stages of water imbibition were analyzed. Variant interaction network derived 205 downstream interacting partners of 7 genes confirmed their role in seed related traits. Interestingly, genes encoding for Type I- Inositol polyphosphate 5 phosphatase1 and E3 Ubiquitin ligase could differentiate parental genotypes, revealed protein conformational deformations and were found to segregate among RILs in coherence with their permeability scores. The 2 identified genes, thus showed a preliminary association with the desirable permeability characteristics.

CONCLUSION

In the light of above outcomes, 2 genes were identified that revealed preliminary, but a relevant association with soybean seed permeability trait and hence could serve as a primary material for understanding the molecular pathways controlling seed permeability traits in soybean.

Current advances in chickpea genomics: applications and future perspectives

Chickpea genomics promises to illuminate our understanding of genome organization, structural variations, evolutionary and domestication-related insights and fundamental biology of legume crops. Unprecedented advancements of next generation sequencing (NGS) technologies have enabled in decoding of multiple chickpea genome sequences and generating huge genomic resources in chickpea both at functional and structural level. This review is aimed to update the current progress of chickpea genomics ranging from high density linkage map development, genome-wide association studies (GWAS), functional genomics resources for various traits, emerging role of abiotic stress responsive coding and non-coding RNAs after the completion of draft chickpea genome sequences. Additionally, the current efforts of whole genome re-sequencing (WGRS) approach of global chickpea germplasm to capture the global genetic diversity existing in the historically released varieties across the world and increasing the resolution of the previously identified candidate gene(s) of breeding importance have been discussed. Thus, the outcomes of these genomics resources will assist in genomics-assisted selection and facilitate breeding of climate-resilient chickpea cultivars for sustainable agriculture.

Transcriptomic signature of drought response in pearl millet (Pennisetum glaucum (L.) and development of web-genomic resources

Pearl millet, (Pennisetum glaucum L.), an efficient (C4) crop of arid/semi-arid regions is known for hardiness. Crop is valuable for bio-fortification combating malnutrition and diabetes, higher caloric value and wider climatic resilience. Limited studies are done in pot-based experiments for drought response at gene-expression level, but field-based experiment mimicking drought by withdrawal of irrigation is still warranted. We report de novo assembly-based transcriptomic signature of drought response induced by irrigation withdrawal in pearl millet. We found 19983 differentially expressed genes, 7595 transcription factors, gene regulatory network having 45 hub genes controlling drought response. We report 34652 putative markers (4192 simple sequence repeats, 12111 SNPs and 6249 InDels). Study reveals role of purine and tryptophan metabolism in ABA accumulation mediating abiotic response in which MAPK acts as major intracellular signal sensing drought. Results were validated by qPCR of 13 randomly selected genes. We report the first web-based genomic resource ( http://webtom.cabgrid.res.in/pmdtdb/ ) which can be used for candidate genes-based SNP discovery programs and trait-based association studies. Looking at climatic change, nutritional and pharmaceutical importance of this crop, present investigation has immense value in understanding drought response in field condition. This is important in germplasm management and improvement in endeavour of pearl millet productivity.

Genomic diversity and macroecology of the crop wild relatives of domesticated pea

There is growing interest in the conservation and utilization of crop wild relatives (CWR) in international food security policy and research. Legumes play an important role in human health, sustainable food production, global food security, and the resilience of current agricultural systems. Pea belongs to the ancient set of cultivated plants of the Near East domestication center and remains an important crop today. Based on genome-wide analysis, P. fulvum was identified as a well-supported species, while the diversity of wild P. sativum subsp. elatius was structured into 5 partly geographically positioned clusters. We explored the spatial and environmental patterns of two progenitor species of domesticated pea in the Mediterranean Basin and in the Fertile Crescent in relation to the past and current climate. This study revealed that isolation by distance does not explain the genetic structure of P. sativum subsp. elatius in its westward expansion from its center of origin. The genetic diversity of wild pea may be driven by Miocene-Pliocene events, while the phylogenetic diversity centers may reflect Pleisto-Holocene climatic changes. These findings help set research and discussion priorities and provide geographical and ecological information for germplasm-collecting missions, as well as for the preservation of extant diversity in ex-situ collections.

Food Legumes and Rising Temperatures: Effects, Adaptive Functional Mechanisms Specific to Reproductive Growth Stage and Strategies to Improve Heat Tolerance

Ambient temperatures are predicted to rise in the future owing to several reasons associated with global climate changes. These temperature increases can result in heat stress- a severe threat to crop production in most countries. Legumes are well-known for their impact on agricultural sustainability as well as their nutritional and health benefits. Heat stress imposes challenges for legume crops and has deleterious effects on the morphology, physiology, and reproductive growth of plants. High-temperature stress at the time of the reproductive stage is becoming a severe limitation for production of grain legumes as their cultivation expands to warmer environments and temperature variability increases due to climate change. The reproductive period is vital in the life cycle of all plants and is susceptible to high-temperature stress as various metabolic processes are adversely impacted during this phase, which reduces crop yield. Food legumes exposed to high-temperature stress during reproduction show flower abortion, pollen and ovule infertility, impaired fertilization, and reduced seed filling, leading to smaller seeds and poor yields. Through various breeding techniques, heat tolerance in major legumes can be enhanced to improve performance in the field. Omics approaches unravel different mechanisms underlying thermotolerance, which is imperative to understand the processes of molecular responses toward high-temperature stress.

A Multiple QTL-Seq Strategy Delineates Potential Genomic Loci Governing Flowering Time in Chickpea

Identification of functionally relevant potential genomic loci using an economical, simpler and user-friendly genomics-assisted breeding strategy is vital for rapid genetic dissection of complex flowering time quantitative trait in chickpea. A high-throughput multiple QTL-seq strategy was employed in two inter (Cicer arietinum desi accession ICC 4958 × C reticulatum wild accession ICC 17160)- and intra (ICC 4958 × C. arietinum kabuli accession ICC 8261)-specific RIL mapping populations to identify the major QTL genomic regions governing flowering time in chickpea. The whole genome resequencing discovered 1635117 and 592486 SNPs exhibiting differentiation between early- and late-flowering mapping parents and bulks, constituted by pooling the homozygous individuals of extreme flowering time phenotypic trait from each of two aforesaid RIL populations. The multiple QTL-seq analysis using these mined SNPs in two RIL mapping populations narrowed-down two longer (907.1 kb and 1.99 Mb) major flowering time QTL genomic regions into the high-resolution shorter (757.7 kb and 1.39 Mb) QTL intervals on chickpea chromosome 4. This essentially identified regulatory as well as coding (non-synonymous/synonymous) novel SNP allelic variants from two efl1 (early flowering 1) and GI (GIGANTEA) genes regulating flowering time in chickpea. Interestingly, strong natural allelic diversity reduction (88-91%) of two known flowering genes especially mapped at major QTL intervals as compared to that of background genomic regions (where no flowering time QTLs were mapped; 61.8%) in cultivated vis-à-vis wild Cicer gene pools was evident inferring the significant impact of evolutionary bottlenecks on these loci during chickpea domestication. Higher association potential of coding non-synonymous and regulatory SNP alleles mined from efl1 (36-49%) and GI (33-42%) flowering genes for early and late flowering time differentiation among chickpea accessions was evident. The robustness and validity of two functional allelic variants-containing genes localized at major flowering time QTLs was apparent by their identification from multiple intra-/inter-specific mapping populations of chickpea. The functionally relevant molecular tags delineated can be of immense use for deciphering the natural allelic diversity-based domestication pattern of flowering time and expediting genomics-aided crop improvement to develop early flowering cultivars of chickpea.

A repeat length variation in myo-inositol monophosphatase gene contributes to seed size trait in chickpea

Chickpea (Cicer arietinum L.) is the third most important food legume crop. Seed size is the most economically important trait for chickpea. To understand the genetic regulation of seed size in chickpea, the present study established a three-way association of CT repeat length variation of a simple sequence repeat (SSR) in myo-inositol monophosphatase gene (CaIMP) with seed weight and phytic acid content by large scale validation and genotyping in a set of genetically diverse germplasm accessions and two reciprocal intra-specific mapping populations. Germplasms and mapping individuals with CT repeat-length expansion in the 5' untranslated region of CaIMP exhibited a pronounced increase in CaIMP protein level, enzymatic activity, seed-phytate content and seed weight. A chickpea transient expression system demonstrated this repeat-length variation influenced the translation of CaIMP mRNA, apparently by facilitating translation initiation. Our analyses proposed that the SSR marker derived from 5' UTR of a CaIMP gene is a promising candidate for selection of seed size/weight for agronomic trait improvement of chickpea.

Genetic dissection of plant growth habit in chickpea

A combinatorial genomics-assisted breeding strategy encompassing association analysis, genetic mapping and expression profiling is found most promising for quantitative dissection of complex traits in crop plants. The present study employed GWAS (genome-wide association study) using 24,405 SNPs (single nucleotide polymorphisms) obtained with genotyping-by-sequencing (GBS) of 92 sequenced desi and kabuli accessions of chickpea. This identified eight significant genomic loci associated with erect (E)/semi-erect (SE) vs. spreading (S)/semi-spreading (SS)/prostrate (P) plant growth habit (PGH) trait differentiation regardless of diverse desi and kabuli genetic backgrounds of chickpea. These associated SNPs in combination explained 23.8% phenotypic variation for PGH in chickpea. Five PGH-associated genes were validated successfully in E/SE and SS/S/P PGH-bearing parental accessions and homozygous individuals of three intra- and interspecific RIL (recombinant inbred line) mapping populations as well as 12 contrasting desi and kabuli chickpea germplasm accessions by selective genotyping through Sequenom MassARRAY. The shoot apical, inflorescence and floral meristems-specific expression, including upregulation (seven-fold) of five PGH-associated genes especially in germplasm accessions and homozygous RIL mapping individuals contrasting with E/SE PGH traits was apparent. Collectively, this integrated genomic strategy delineated diverse non-synonymous SNPs from five candidate genes with strong allelic effects on PGH trait variation in chickpea. Of these, two vernalization-responsive non-synonymous SNP alleles carrying SNF2 protein-coding gene and B3 transcription factor associated with PGH traits were found to be the most promising in chickpea. The SNP allelic variants associated with E/SE/SS/S PGH trait differentiation were exclusively present in all cultivated desi and kabuli chickpea accessions while wild species/accessions belonging to primary, secondary and tertiary gene pools mostly contained prostrate PGH-associated SNP alleles. This indicates strong adaptive natural/artificial selection pressure (Tajima's D 3.15 to 4.57) on PGH-associated target genomic loci during chickpea domestication. These vital leads thus have potential to decipher complex transcriptional regulatory gene function of PGH trait differentiation and for understanding the selective sweep-based PGH trait evolution and domestication pattern in cultivated and wild chickpea accessions adapted to diverse agroclimatic conditions. Collectively, the essential inputs generated will be of profound use in marker-assisted genetic enhancement to develop cultivars with desirable plant architecture of erect growth habit types in chickpea.

Genetic Diversity, Population Structure, and Linkage Disequilibrium of a Core Collection of Ziziphus jujuba Assessed with Genome-wide SNPs Developed by Genotyping-by-sequencing and SSR Markers

Chinese jujube (Ziziphus jujuba Mill) is an economically important fruit species native to China with high nutritious and medicinal value. Genotyping-by-sequencing was used to detect and genotype single nucleotide polymorphisms (SNPs) in a core collection of 150 Chinese jujube accessions and further to characterize their genetic diversity, population structure, and linkage disequilibrium (LD). A total of 4,680 high-quality SNPs were identified, of which 38 sets of tri-allelic SNPs were detected. The average polymorphism information content (PIC) values based on bi-allelic SNPs and tri-allelic SNPs were 0.27 and 0.38, respectively. STRUCTURE and principal coordinate analyses based on SNPs revealed that the 150 accessions could be clustered into two groups. However, neighbor-joining trees indicated the accessions should be grouped into three major clusters. Our data confirm that the resolving power for genetic diversity was similar for the SSRs and SNPs. In contrast, regarding population structure, the resolving power was higher for SSRs than for SNPs. The LD pattern in Chinese jujube was investigated for the first time. We observed a relatively rapid LD decay with a short range (∼10 kb) for all pseudo-chromosomes and for individual pseudo-chromosomes. Our findings provide important information for future genome-wide association analyses and marker-assisted selective breeding of Chinese jujube.

Genome-Wide Single-Nucleotide Polymorphisms in CMS and Restorer Lines Discovered by Genotyping Using Sequencing and Association with Marker-Combining Ability for 12 Yield-Related Traits in Oryza sativa L. subsp. Japonica

Heterosis or hybrid vigor is closely related with general combing ability (GCA) of parents and special combining ability (SCA) of combinations. The evaluation of GCA and SCA facilitate selection of parents and combinations in heterosis breeding. In order to improve combining ability (CA) by molecular marker assist selection, it is necessary to identify marker loci associated with the CA. To identify the single nucleotide polymorphisms (SNP) loci associated with CA in the parental genomes of japonica rice, genome-wide discovered SNP loci were tested for association with the CA of 18 parents for 12 yield-related traits. In this study, 81 hybrids were created and evaluated to calculate the CA of 18 parents. The parents were sequenced by genotyping by sequencing (GBS) method for identification of genome-wide SNPs. The analysis of GBS indicated that the successful mapping of 9.86 × 10⁶ short reads in the Nipponbare reference genome consists of 39,001 SNPs in parental genomes at 11,085 chromosomal positions. The discovered SNPs were non-randomly distributed within and among the 12 chromosomes of rice. Overall, 20.4% (8026) of the discovered SNPs were coding types, and 8.6% (3344) and 9.9% (3951) of the SNPs revealed synonymous and non-synonymous changes, which provide valuable knowledge about the underlying performance of the parents. Furthermore, the associations between SNPs and CA indicated that 362 SNP loci were significantly related to the CA of 12 parental traits. The identified SNP loci of CA in our study were distributed genome wide and caused a positive or negative effect on the CA of traits. For the yield-related traits, such as grain thickness, days to heading, panicle length, grain length and 1000-grain weight, a maximum number of positive SNP loci of CA were found in CMS A171 and in the restorers LC64 and LR27. On an individual basis, some of associated loci that resided on chromosomes 2, 5, 7, 9, and 11 recorded maximum positive values for the CA of traits. From our results, we suggest that heterosis in japonica rice would be improved by pyramiding the favorable SNP loci of CA and eliminating the unfavorable loci from parental genomes.

Draft genome sequence of Cicer reticulatum L., the wild progenitor of chickpea provides a resource for agronomic trait improvement

Cicer reticulatum L. is the wild progenitor of the fourth most important legume crop chickpea (C. arietinum L.). We assembled short-read sequences into 416 Mb draft genome of C. reticulatum and anchored 78% (327 Mb) of this assembly to eight linkage groups. Genome annotation predicted 25,680 protein-coding genes covering more than 90% of predicted gene space. The genome assembly shared a substantial synteny and conservation of gene orders with the genome of the model legume Medicago truncatula. Resistance gene homologs of wild and domesticated chickpeas showed high sequence homology and conserved synteny. Comparison of gene sequences and nucleotide diversity using 66 wild and domesticated chickpea accessions suggested that the desi type chickpea was genetically closer to the wild species than the kabuli type. Comparative analyses predicted gene flow between the wild and the cultivated species during domestication. Molecular diversity and population genetic structure determination using 15,096 genome-wide single nucleotide polymorphisms revealed an admixed domestication pattern among cultivated (desi and kabuli) and wild chickpea accessions belonging to three population groups reflecting significant influence of parentage or geographical origin for their cultivar-specific population classification. The assembly and the polymorphic sequence resources presented here would facilitate the study of chickpea domestication and targeted use of wild Cicer germplasms for agronomic trait improvement in chickpea.

Food Legumes and Rising Temperatures: Effects, Adaptive Functional Mechanisms Specific to Reproductive Growth Stage and Strategies to Improve Heat Tolerance

Ambient temperatures are predicted to rise in the future owing to several reasons associated with global climate changes. These temperature increases can result in heat stress- a severe threat to crop production in most countries. Legumes are well-known for their impact on agricultural sustainability as well as their nutritional and health benefits. Heat stress imposes challenges for legume crops and has deleterious effects on the morphology, physiology, and reproductive growth of plants. High-temperature stress at the time of the reproductive stage is becoming a severe limitation for production of grain legumes as their cultivation expands to warmer environments and temperature variability increases due to climate change. The reproductive period is vital in the life cycle of all plants and is susceptible to high-temperature stress as various metabolic processes are adversely impacted during this phase, which reduces crop yield. Food legumes exposed to high-temperature stress during reproduction show flower abortion, pollen and ovule infertility, impaired fertilization, and reduced seed filling, leading to smaller seeds and poor yields. Through various breeding techniques, heat tolerance in major legumes can be enhanced to improve performance in the field. Omics approaches unravel different mechanisms underlying thermotolerance, which is imperative to understand the processes of molecular responses toward high-temperature stress.

Breaking RAD: an evaluation of the utility of restriction site-associated DNA sequencing for genome scans of adaptation

Understanding how and why populations evolve is of fundamental importance to molecular ecology. Restriction site-associated DNA sequencing (RADseq), a popular reduced representation method, has ushered in a new era of genome-scale research for assessing population structure, hybridization, demographic history, phylogeography and migration. RADseq has also been widely used to conduct genome scans to detect loci involved in adaptive divergence among natural populations. Here, we examine the capacity of those RADseq-based genome scan studies to detect loci involved in local adaptation. To understand what proportion of the genome is missed by RADseq studies, we developed a simple model using different numbers of RAD-tags, genome sizes and extents of linkage disequilibrium (length of haplotype blocks). Under the best-case modelling scenario, we found that RADseq using six- or eight-base pair cutting restriction enzymes would fail to sample many regions of the genome, especially for species with short linkage disequilibrium. We then surveyed recent studies that have used RADseq for genome scans and found that the median density of markers across these studies was 4.08 RAD-tag markers per megabase (one marker per 245 kb). The length of linkage disequilibrium for many species is one to three orders of magnitude less than density of the typical recent RADseq study. Thus, we conclude that genome scans based on RADseq data alone, while useful for studies of neutral genetic variation and genetic population structure, will likely miss many loci under selection in studies of local adaptation.