Genetic structure and diversity analysis of the primary gene pool of chickpea using SSR markers

Genome-wide association study as a powerful tool for dissecting competitive traits in legumes

Legumes are extremely valuable because of their high protein content and several other nutritional components. The major challenge lies in maintaining the quantity and quality of protein and other nutritional compounds in view of climate change conditions. The global need for plant-based proteins has increased the demand for seeds with a high protein content that includes essential amino acids. Genome-wide association studies (GWAS) have evolved as a standard approach in agricultural genetics for examining such intricate characters. Recent development in machine learning methods shows promising applications for dimensionality reduction, which is a major challenge in GWAS. With the advancement in biotechnology, sequencing, and bioinformatics tools, estimation of linkage disequilibrium (LD) based associations between a genome-wide collection of single-nucleotide polymorphisms (SNPs) and desired phenotypic traits has become accessible. The markers from GWAS could be utilized for genomic selection (GS) to predict superior lines by calculating genomic estimated breeding values (GEBVs). For prediction accuracy, an assortment of statistical models could be utilized, such as ridge regression best linear unbiased prediction (rrBLUP), genomic best linear unbiased predictor (gBLUP), Bayesian, and random forest (RF). Both naturally diverse germplasm panels and family-based breeding populations can be used for association mapping based on the nature of the breeding system (inbred or outbred) in the plant species. MAGIC, MCILs, RIAILs, NAM, and ROAM are being used for association mapping in several crops. Several modifications of NAM, such as doubled haploid NAM (DH-NAM), backcross NAM (BC-NAM), and advanced backcross NAM (AB-NAM), have also been used in crops like rice, wheat, maize, barley mustard, etc. for reliable marker-trait associations (MTAs), phenotyping accuracy is equally important as genotyping. Highthroughput genotyping, phenomics, and computational techniques have advanced during the past few years, making it possible to explore such enormous datasets. Each population has unique virtues and flaws at the genomics and phenomics levels, which will be covered in more detail in this review study. The current investigation includes utilizing elite breeding lines as association mapping population, optimizing the choice of GWAS selection, population size, and hurdles in phenotyping, and statistical methods which will analyze competitive traits in legume breeding.

Genome-wide discovery of di-nucleotide SSR markers based on whole genome re-sequencing data of Cicer arietinum L. and Cicer reticulatum Ladiz

Simple sequence repeats (SSRs) are valuable genetic markers due to their co-dominant inheritance, multi-allelic and reproducible nature. They have been largely used for exploiting genetic architecture of plant germplasms, phylogenetic analysis, and mapping studies. Among the SSRs, di-nucleotide repeats are the most frequent of the simple repeats distributed throughout the plant genomes. In present study, we aimed to discover and develop di-nucleotide SSR markers by using the whole genome re-sequencing (WGRS) data from Cicer arietinum L. and C. reticulatum Ladiz. A total of 35,329 InDels were obtained in C. arietinum, whereas 44,331 InDels in C. reticulatum. 3387 InDels with 2 bp length were detected in C. arietinum, there were 4704 in C. reticulatum. Among 8091 InDels, 58 di-nucleotide regions that were polymorphic between two species were selected and used for validation. We tested primers for evaluation of genetic diversity in 30 chickpea genotypes including C. arietinum, C. reticulatum, C. echinospermum P.H. Davis, C. anatolicum Alef., C. canariense A. Santos & G.P. Lewis, C. microphyllum Benth., C. multijugum Maesen, C. oxyodon Boiss. & Hohen. and C. songaricum Steph ex DC. A total of 244 alleles were obtained for 58 SSR markers giving an average of 2.36 alleles per locus. The observed heterozygosity was 0.08 while the expected heterozygosity was 0.345. Polymorphism information content was found to be 0.73 across all loci. Phylogenetic tree and principal coordinate analysis clearly divided the accessions into four groups. The SSR markers were also evaluated in 30 genotypes of a RIL population obtained from an interspecific cross between C. arietinum and C. reticulatum. Chi-square (χ²) test revealed an expected 1:1 segregation ratio in the population. These results demonstrated the success of SSR identification and marker development for chickpea with the use of WGRS data. The newly developed 58 SSR markers are expected to be useful for chickpea breeders.

SSR markers in revealing extent of genetic diversity and phylogenetic relationships among chickpea core collection accessions for Western Himalayas

BACKGROUND

The exploration of genetic diversity is the key source of germplasm conservation and potential to broaden its genetic base. The globally growing demand for chickpea suggests superior/climate-resilient varieties, which in turn necessitates the germplasm characterization to unravel underlying genetic variation.

METHODOLOGY AND RESULTS

A chickpea core collection comprising of diverse 192 accessions which include cultivated Cicer arietinum, and wild C. reticulatum, C. echinospermum, and C. microphyllum species were investigated to analyze their genetic diversity and relationship, by assaying 33 unlinked simple sequence repeat (SSR) markers. The results amplified a total of 323 alleles (Na), ranging from 2 to 8 with an average of 4.25 alleles per locus. Expected heterozygosity (He) differed from 0.46 to 0.86 with an average of 0.68. Polymorphic information content (PIC) ranged from 0.73 to 0.98 with an average of 0.89. Analysis of molecular variance (AMOVA) showed that most of the variation was among individuals (87%). Cluster analysis resulted in the formation of four distinct clusters. Cluster I represented all cultivated and clusters II, III, and IV comprised a heterogeneous group of cultivated and wild chickpea accessions.

CONCLUSION

We report considerable diversity and greater resolving power of SSR markers for assessing variability and interrelationship among the chickpea accessions. The chickpea core is expected to be an efficient resource for breeders for broadening the chickpea genetic base and could be useful for selective breeding of desirable traits and in the identification of target genes for genomics-assisted breeding.

A Comprehensive Review on Chickpea (Cicer arietinum L.) Breeding for Abiotic Stress Tolerance and Climate Change Resilience

Chickpea is one of the most important pulse crops worldwide, being an excellent source of protein. It is grown under rain-fed conditions averaging yields of 1 t/ha, far from its potential of 6 t/ha under optimum conditions. The combined effects of heat, cold, drought, and salinity affect species productivity. In this regard, several physiological, biochemical, and molecular mechanisms are reviewed to confer tolerance to abiotic stress. A large collection of nearly 100,000 chickpea accessions is the basis of breeding programs, and important advances have been achieved through conventional breeding, such as germplasm introduction, gene/allele introgression, and mutagenesis. In parallel, advances in molecular biology and high-throughput sequencing have allowed the development of specific molecular markers for the genus Cicer, facilitating marker-assisted selection for yield components and abiotic tolerance. Further, transcriptomics, proteomics, and metabolomics have permitted the identification of specific genes, proteins, and metabolites associated with tolerance to abiotic stress of chickpea. Furthermore, some promising results have been obtained in studies with transgenic plants and with the use of gene editing to obtain drought-tolerant chickpea. Finally, we propose some future lines of research that may be useful to obtain chickpea genotypes tolerant to abiotic stress in a scenario of climate change.

Genetic variability and population structure of Ethiopian chickpea (Cicer arietinum L.) germplasm

Evaluation of the genetic diversity and an understanding of the genetic structure and relationships of chickpea genotypes are valuable to design efficient germplasm conservation strategies and crop breeding programs. Information is limited, in these regards, for Ethiopian chickpea germplasms. Therefore, the present study was carried out to estimate the genetic diversity, population structure, and relationships of 152 chickpea genotypes using simple sequence repeats (SSR) markers. Twenty three SSR markers exhibited polymorphism producing a total of 133 alleles, with a mean of 5.8 alleles per locus. Analyses utilizing various genetic-based statistics included pairwise population Nei’s genetic distance, heterozygosity, Shannon’s information index, polymorphic information content, and percent polymorphism. These analyses exemplified the existence of high genetic variation within and among chickpea genotypes. The 152 genotypes were divided into two major clusters based on Nei’s genetic distances. The exotic genotypes were grouped in one cluster exclusively showing that these genotypes are distinct to Ethiopian genotypes, while the patterns of clustering of Ethiopian chickpea genotypes based on their geographic region were not consistent because of the seed exchange across regions. Model-based population structure clustering identified two discrete populations. These finding provides useful insight for chickpea collections and ex-situ conservation and national breeding programs for widening the genetic base of chickpea.

Molecular Genetic Diversity and Population Structure in Ethiopian Chickpea Germplasm Accessions

Chickpea (Cicer arietinum L.) is a cheap source of protein and rich in minerals for people living in developing countries. In order to assess the existing molecular genetic diversity and determine population structures in selected Ethiopian chickpea germplasm accessions (118), a set of 46 simple sequence repeat (SSR) markers equally distributed on the chickpea genome were genotyped. A total of 572 alleles were detected from 46 SSR markers, and the number of alleles per locus varied from 2 (ICCM0289) to 28 (TA22). The average number of alleles per locus, polymorphism information content, and expected heterozygosity were 12, 0.684, and 0.699, respectively. Phylogenetic analysis grouped the 118 chickpea genotypes from diverse sources into three evolutionary and/or biological groups (improved desi, improved kabuli, and landraces). The population structure analysis revealed six sub-populations from 118 chickpea genotypes studied. AMOVA revealed that 57%, 29%, and 14% of the total genetic variations were observed among individuals, within populations, and among populations. The insights into the genetic diversity at molecular levels in the Ethiopian germplasm lines can be used for designing conservation strategies as well as the diverse germplasm lines identified in this study can be used for trait dissection and trait improvement.

Population structure and association analysis of heat stress relevant traits in chickpea (Cicer arietinum L.)

Understanding genetic diversity and population structure is prerequisite to broaden the cultivated base of any crop. In the current investigation, we report discovery of a total of 319 alleles by assaying 81 SSRs on 71 chickpea genotypes. The cluster analysis based on Jaccard coefficient and unweighted neighbor joining algorithm categorized all genotypes into two major clusters. Cultivars grown within the same agro-climatic zones were clustered together, whereas the remaining genotypes particularly advanced breeding lines and accessions assigned to another cluster. Population structure analysis separated the entire collection into two subpopulations (K = 2) and the clustering pattern remained in close agreement with those of distance-based methods. Importantly, we also discovered marker trait association for membrane stability index (MSI) and leaf chlorophyll content measured as SPAD chlorophyll meter reading (SCMR), the two important physiological parameters indicative of heat stress (HS) tolerance in chickpea. Association analysis using both general linear and mixed linear models of the mean phenotypic data of traits recorded in 2016 and 2017 uncovered significant association of NCPGR206 and H2L102 with the MSI trait. Likewise, SSR markers GA9, TR31 and TA113 exhibited significant association with SCMR trait. The genomic regions putatively linked with two traits may be investigated in greater detail to further improve knowledge about the genetic architecture of HS tolerance in chickpea.

Population structure and association analysis of heat stress relevant traits in chickpea (Cicer arietinum L.)

Understanding genetic diversity and population structure is prerequisite to broaden the cultivated base of any crop. In the current investigation, we report discovery of a total of 319 alleles by assaying 81 SSRs on 71 chickpea genotypes. The cluster analysis based on Jaccard coefficient and unweighted neighbor joining algorithm categorized all genotypes into two major clusters. Cultivars grown within the same agro-climatic zones were clustered together, whereas the remaining genotypes particularly advanced breeding lines and accessions assigned to another cluster. Population structure analysis separated the entire collection into two subpopulations (K = 2) and the clustering pattern remained in close agreement with those of distance-based methods. Importantly, we also discovered marker trait association for membrane stability index (MSI) and leaf chlorophyll content measured as SPAD chlorophyll meter reading (SCMR), the two important physiological parameters indicative of heat stress (HS) tolerance in chickpea. Association analysis using both general linear and mixed linear models of the mean phenotypic data of traits recorded in 2016 and 2017 uncovered significant association of NCPGR206 and H2L102 with the MSI trait. Likewise, SSR markers GA9, TR31 and TA113 exhibited significant association with SCMR trait. The genomic regions putatively linked with two traits may be investigated in greater detail to further improve knowledge about the genetic architecture of HS tolerance in chickpea.

Genetic Diversity in Lens Species Revealed by EST and Genomic Simple Sequence Repeat Analysis

Low productivity of pilosae type lentils grown in South Asia is attributed to narrow genetic base of the released cultivars which results in susceptibility to biotic and abiotic stresses. For enhancement of productivity and production, broadening of genetic base is essentially required. The genetic base of released cultivars can be broadened by using diverse types including bold seeded and early maturing lentils from Mediterranean region and related wild species. Genetic diversity in eighty six accessions of three species of genus Lens was assessed based on twelve genomic and thirty one EST-SSR markers. The evaluated set of genotypes included diverse lentil varieties and advanced breeding lines from Indian programme, two early maturing ICARDA lines and five related wild subspecies/species endemic to the Mediterranean region. Genomic SSRs exhibited higher polymorphism in comparison to EST SSRs. GLLC 598 produced 5 alleles with highest gene diversity value of 0.80. Among the studied subspecies/species 43 SSRs detected maximum number of alleles in L. orientalis. Based on Nei's genetic distance cultivated lentil L. culinaris subsp. culinaris was found to be close to its wild progenitor L. culinaris subsp. orientalis. The Prichard's structure of 86 genotypes distinguished different subspecies/species. Higher variability was recorded among individuals within population than among populations.

Repeat length variation in the 5'UTR of myo-inositol monophosphatase gene is related to phytic acid content and contributes to drought tolerance in chickpea (Cicer arietinum L.)

Myo-inositol metabolism plays a significant role in plant growth and development, and is also used as a precursor for many important metabolites, such as ascorbate, pinitol, and phytate. Phytate (inositol hexakisphosphate) is the major storage pool for phosphate in the seeds. It is utilized during seed germination and growth of the developing embryo. In addition, it is implicated in protection against oxidative stress. In the present study, a panel of chickpea accessions was used for an association analysis. Association analysis accounting for population structure and relative kinship identified alleles of a simple sequence repeat marker, NCPGR90, that are associated with both phytic acid content and drought tolerance. These alleles varied with respect to the dinucleotide CT repeats present within the marker. NCPGR90 located to the 5'UTR of chickpea myo-inositol monophosphatase gene (CaIMP) and showed transcript length variation in drought-tolerant and drought-susceptible accessions. CaIMP from a drought-tolerant accession with a smaller repeat was almost 2-fold upregulated as compared to a susceptible accession having a longer repeat, even under control non-stressed conditions. This study suggests an evolution of simple sequence repeat length variation in CaIMP, which might be regulating phytic acid levels to confer drought tolerance in natural populations of chickpea.

CTDB: An Integrated Chickpea Transcriptome Database for Functional and Applied Genomics

Chickpea is an important grain legume used as a rich source of protein in human diet. The narrow genetic diversity and limited availability of genomic resources are the major constraints in implementing breeding strategies and biotechnological interventions for genetic enhancement of chickpea. We developed an integrated Chickpea Transcriptome Database (CTDB), which provides the comprehensive web interface for visualization and easy retrieval of transcriptome data in chickpea. The database features many tools for similarity search, functional annotation (putative function, PFAM domain and gene ontology) search and comparative gene expression analysis. The current release of CTDB (v2.0) hosts transcriptome datasets with high quality functional annotation from cultivated (desi and kabuli types) and wild chickpea. A catalog of transcription factor families and their expression profiles in chickpea are available in the database. The gene expression data have been integrated to study the expression profiles of chickpea transcripts in major tissues/organs and various stages of flower development. The utilities, such as similarity search, ortholog identification and comparative gene expression have also been implemented in the database to facilitate comparative genomic studies among different legumes and Arabidopsis. Furthermore, the CTDB represents a resource for the discovery of functional molecular markers (microsatellites and single nucleotide polymorphisms) between different chickpea types. We anticipate that integrated information content of this database will accelerate the functional and applied genomic research for improvement of chickpea. The CTDB web service is freely available at http://nipgr.res.in/ctdb.html.

Genome-wide high-throughput SNP discovery and genotyping for understanding natural (functional) allelic diversity and domestication patterns in wild chickpea

We identified 82489 high-quality genome-wide SNPs from 93 wild and cultivated Cicer accessions through integrated reference genome- and de novo-based GBS assays. High intra- and inter-specific polymorphic potential (66-85%) and broader natural allelic diversity (6-64%) detected by genome-wide SNPs among accessions signify their efficacy for monitoring introgression and transferring target trait-regulating genomic (gene) regions/allelic variants from wild to cultivated Cicer gene pools for genetic improvement. The population-specific assignment of wild Cicer accessions pertaining to the primary gene pool are more influenced by geographical origin/phenotypic characteristics than species/gene-pools of origination. The functional significance of allelic variants (non-synonymous and regulatory SNPs) scanned from transcription factors and stress-responsive genes in differentiating wild accessions (with potential known sources of yield-contributing and stress tolerance traits) from cultivated desi and kabuli accessions, fine-mapping/map-based cloning of QTLs and determination of LD patterns across wild and cultivated gene-pools are suitably elucidated. The correlation between phenotypic (agromorphological traits) and molecular diversity-based admixed domestication patterns within six structured populations of wild and cultivated accessions via genome-wide SNPs was apparent. This suggests utility of whole genome SNPs as a potential resource for identifying naturally selected trait-regulating genomic targets/functional allelic variants adaptive to diverse agroclimatic regions for genetic enhancement of cultivated gene-pools.

Characterization of genetic diversity in chickpea using SSR markers, Start Codon Targeted Polymorphism (SCoT) and Conserved DNA-Derived Polymorphism (CDDP)

To evaluate the genetic diversity among 48 genotypes of chickpea comprising cultivars, landraces and internationally developed improved lines genetic distances were evaluated using three different molecular marker techniques: Simple Sequence Repeat (SSR); Start Codon Targeted (SCoT) and Conserved DNA-derived Polymorphism (CDDP). Average polymorphism information content (PIC) for SSR, SCoT and CDDP markers was 0.47, 0.45 and 0.45, respectively, and this revealed that three different marker types were equal for the assessment of diversity amongst genotypes. Cluster analysis for SSR and SCoT divided the genotypes in to three distinct clusters and using CDDP markers data, genotypes grouped in to five clusters. There were positive significant correlation (r = 0.43, P < 0.01) between similarity matrix obtained by SCoT and CDDP. Three different marker techniques showed relatively same pattern of diversity across genotypes and using each marker technique it's obvious that diversity pattern and polymorphism for varieties were higher than that of genotypes, and CDDP had superiority over SCoT and SSR markers. These results suggest that efficiency of SSR, SCOT and CDDP markers was relatively the same in fingerprinting of chickpea genotypes. To our knowledge, this is the first detailed report of using targeted DNA region molecular marker (CDDP) for genetic diversity analysis in chickpea in comparison with SCoT and SSR markers. Overall, our results are able to prove the suitability of SCoT and CDDP markers for genetic diversity analysis in chickpea for their high rates of polymorphism and their potential for genome diversity and germplasm conservation.

Deploying QTL-seq for rapid delineation of a potential candidate gene underlying major trait-associated QTL in chickpea

A rapid high-resolution genome-wide strategy for molecular mapping of major QTL(s)/gene(s) regulating important agronomic traits is vital for in-depth dissection of complex quantitative traits and genetic enhancement in chickpea. The present study for the first time employed a NGS-based whole-genome QTL-seq strategy to identify one major genomic region harbouring a robust 100-seed weight QTL using an intra-specific 221 chickpea mapping population (desi cv. ICC 7184 × desi cv. ICC 15061). The QTL-seq-derived major SW QTL (CaqSW1.1) was further validated by single-nucleotide polymorphism (SNP) and simple sequence repeat (SSR) marker-based traditional QTL mapping (47.6% R² at higher LOD >19). This reflects the reliability and efficacy of QTL-seq as a strategy for rapid genome-wide scanning and fine mapping of major trait regulatory QTLs in chickpea. The use of QTL-seq and classical QTL mapping in combination narrowed down the 1.37 Mb (comprising 177 genes) major SW QTL (CaqSW1.1) region into a 35 kb genomic interval on desi chickpea chromosome 1 containing six genes. One coding SNP (G/A)-carrying constitutive photomorphogenic9 (COP9) signalosome complex subunit 8 (CSN8) gene of these exhibited seed-specific expression, including pronounced differential up-/down-regulation in low and high seed weight mapping parents and homozygous individuals during seed development. The coding SNP mined in this potential seed weight-governing candidate CSN8 gene was found to be present exclusively in all cultivated species/genotypes, but not in any wild species/genotypes of primary, secondary and tertiary gene pools. This indicates the effect of strong artificial and/or natural selection pressure on target SW locus during chickpea domestication. The proposed QTL-seq-driven integrated genome-wide strategy has potential to delineate major candidate gene(s) harbouring a robust trait regulatory QTL rapidly with optimal use of resources. This will further assist us to extrapolate the molecular mechanism underlying complex quantitative traits at a genome-wide scale leading to fast-paced marker-assisted genetic improvement in diverse crop plants, including chickpea.

Evaluation of microsatellite-based genetic diversity, protein and mineral content in chickpea accessions grown in Kyrgyzstan

The genetic diversity of 23 chickpea accessions representing Kyrgyz landraces and cultivars, ICARDA breeding lines, Spanish and Turkish cultivars was characterized using nine microsatellite (SSR) markers which generated a total of 122 alleles. The number of alleles (Na) per locus varied from 9 to 20. The observed heterozygosity (Ho) ranged between 0.05 and 0.43 (average 0.13) whereas both the expected heterozygosity (He) and polymorphic information content (PIC) ranged from 0.71 to 0.90 (average 0.83). Analysis of molecular variance (AMOVA) showed that 62% of the total genetic variation was found within accessions while the remaining 38% was found among accessions. Principal coordinate analysis (PCoA) indicated the presence of two groups. The two Kyrgyz cultivars were found apart from these groups. Cluster analysis generally confirmed the results of PCoA and also separated the Kyrgyz cultivars from the subcluster formed by Kyrgyz landraces and the subclusters formed by breeding lines from ICARDA along with landraces from Turkey and Spain. In addition, protein content and mineral concentration were determined. Protein content and mineral concentrations for Ca, S, Mg, P, K, Fe, Mn, Cu and Zn varied significantly among accessions. The results show that Kyrgyz germplasm provides a source of diversity for improvement of chickpea.

An integrated genomic approach for rapid delineation of candidate genes regulating agro-morphological traits in chickpea

The identification and fine mapping of robust quantitative trait loci (QTLs)/genes governing important agro-morphological traits in chickpea still lacks systematic efforts at a genome-wide scale involving wild Cicer accessions. In this context, an 834 simple sequence repeat and single-nucleotide polymorphism marker-based high-density genetic linkage map between cultivated and wild parental accessions (Cicer arietinum desi cv. ICC 4958 and Cicer reticulatum wild cv. ICC 17160) was constructed. This inter-specific genetic map comprising eight linkage groups spanned a map length of 949.4 cM with an average inter-marker distance of 1.14 cM. Eleven novel major genomic regions harbouring 15 robust QTLs (15.6–39.8% R² at 4.2–15.7 logarithm of odds) associated with four agro-morphological traits (100-seed weight, pod and branch number/plant and plant hairiness) were identified and mapped on chickpea chromosomes. Most of these QTLs showed positive additive gene effects with effective allelic contribution from ICC 4958, particularly for increasing seed weight (SW) and pod and branch number. One robust SW-influencing major QTL region (qSW4.2) has been narrowed down by combining QTL mapping with high-resolution QTL region-specific association analysis, differential expression profiling and gene haplotype-based association/LD mapping. This enabled to delineate a strong SW-regulating ABI3VP1 transcription factor (TF) gene at trait-specific QTL interval and consequently identified favourable natural allelic variants and superior high seed weight-specific haplotypes in the upstream regulatory region of this gene showing increased transcript expression during seed development. The genes (TFs) harbouring diverse trait-regulating QTLs, once validated and fine-mapped by our developed rapid integrated genomic approach and through gene/QTL map-based cloning, can be utilized as potential candidates for marker-assisted genetic enhancement of chickpea.

Genomics-assisted breeding for drought tolerance in chickpea

Terminal drought is one of the major constraints in chickpea (Cicer arietinum L.), causing more than 50% production losses. With the objective of accelerating genetic understanding and crop improvement through genomics-assisted breeding, a draft genome sequence has been assembled for the CDC Frontier variety. In this context, 544.73Mb of sequence data were assembled, capturing of 73.8% of the genome in scaffolds. In addition, large-scale genomic resources including several thousand simple sequence repeats and several million single nucleotide polymorphisms, high-density diversity array technology (15360 clones) and Illumina GoldenGate assay genotyping platforms, high-density genetic maps and transcriptome assemblies have been developed. In parallel, by using linkage mapping approach, one genomic region harbouring quantitative trait loci for several drought tolerance traits has been identified and successfully introgressed in three leading chickpea varieties (e.g. JG 11, Chefe, KAK 2) by using a marker-assisted backcrossing approach. A multilocation evaluation of these marker-assisted backcrossing lines provided several lines with 10-24% higher yield than the respective recurrent parents.Modern breeding approaches like marker-assisted recurrent selection and genomic selection are being deployed for enhancing drought tolerance in chickpea. Some novel mapping populations such as multiparent advanced generation intercross and nested association mapping populations are also being developed for trait mapping at higher resolution, as well as for enhancing the genetic base of chickpea. Such advances in genomics and genomics-assisted breeding will accelerate precision and efficiency in breeding for stress tolerance in chickpea.

Natural allelic diversity, genetic structure and linkage disequilibrium pattern in wild chickpea

Characterization of natural allelic diversity and understanding the genetic structure and linkage disequilibrium (LD) pattern in wild germplasm accessions by large-scale genotyping of informative microsatellite and single nucleotide polymorphism (SNP) markers is requisite to facilitate chickpea genetic improvement. Large-scale validation and high-throughput genotyping of genome-wide physically mapped 478 genic and genomic microsatellite markers and 380 transcription factor gene-derived SNP markers using gel-based assay, fluorescent dye-labelled automated fragment analyser and matrix-assisted laser desorption ionization-time of flight (MALDI-TOF) mass array have been performed. Outcome revealed their high genotyping success rate (97.5%) and existence of a high level of natural allelic diversity among 94 wild and cultivated Cicer accessions. High intra- and inter-specific polymorphic potential and wider molecular diversity (11-94%) along with a broader genetic base (13-78%) specifically in the functional genic regions of wild accessions was assayed by mapped markers. It suggested their utility in monitoring introgression and transferring target trait-specific genomic (gene) regions from wild to cultivated gene pool for the genetic enhancement. Distinct species/gene pool-wise differentiation, admixed domestication pattern, and differential genome-wide recombination and LD estimates/decay observed in a six structured population of wild and cultivated accessions using mapped markers further signifies their usefulness in chickpea genetics, genomics and breeding.

Genome wide SNP identification in chickpea for use in development of a high density genetic map and improvement of chickpea reference genome assembly

BACKGROUND

In the whole genome sequencing, genetic map provides an essential framework for accurate and efficient genome assembly and validation. The main objectives of this study were to develop a high-density genetic map using RAD-Seq (Restriction-site Associated DNA Sequencing) genotyping-by-sequencing (RAD-Seq GBS) and Illumina GoldenGate assays, and to examine the alignment of the current map with the kabuli chickpea genome assembly.

RESULTS

Genic single nucleotide polymorphisms (SNPs) totaling 51,632 SNPs were identified by 454 transcriptome sequencing of Cicer arietinum and Cicer reticulatum genotypes. Subsequently, an Illumina GoldenGate assay for 1,536 SNPs was developed. A total of 1,519 SNPs were successfully assayed across 92 recombinant inbred lines (RILs), of which 761 SNPs were polymorphic between the two parents. In addition, the next generation sequencing (NGS)-based GBS was applied to the same population generating 29,464 high quality SNPs. These SNPs were clustered into 626 recombination bins based on common segregation patterns. Data from the two approaches were used for the construction of a genetic map using a population derived from an intraspecific cross. The map consisted of 1,336 SNPs including 604 RAD recombination bins and 732 SNPs from Illumina GoldenGate assay. The map covered 653 cM of the chickpea genome with an average distance between adjacent markers of 0.5 cM. To date, this is the most extensive genetic map of chickpea using an intraspecific population. The alignment of the map with the CDC Frontier genome assembly revealed an overall conserved marker order; however, a few local inconsistencies within the Cicer arietinum pseudochromosome 1 (Ca1), Ca5 and Ca8 were detected. The map enabled the alignment of 215 unplaced scaffolds from the CDC Frontier draft genome assembly. The alignment also revealed varying degrees of recombination rates and hotspots across the chickpea genome.

CONCLUSIONS

A high-density genetic map using RAD-Seq GBS and Illumina GoldenGate assay was developed and aligned with the existing kabuli chickpea draft genome sequence. The analysis revealed an overall conserved marker order, although some localized inversions between draft genome assembly and the genetic map were detected. The current analysis provides an insight of the recombination rates and hotspots across the chickpea genome.

Glycolytic enzyme activities and gene expression in Cicer arietinum exposed to water-deficit stress

The specific activities and transcript levels of glycolytic enzymes were examined in shoots of chickpea (Cicer arietinum L.) cultivars, Pusa362 (drought tolerant) and SBD377 (drought sensitive), subjected to water-deficit stress 30 days after sowing. Water-deficit stress resulted in decrease in relative water content, chlorophyll content, plant dry weight, and NADP/NADPH ratio and increase in NAD/NADH ratio in both the cultivars. A successive decline in the specific activities of fructose-1,6-bisphosphate aldolase (aldolase), 3-phosphoglycerate kinase (PGK), and NADP-glyceraldehyde-3-phosphate dehydrogenase (NADP-GAPDH) and elevation in the specific activities of phosphoglycerate mutase (PGM) and triosephosphate isomerase (TPI) was observed in both the cultivars under stress as compared to their respective control plants. The specific activities of hexokinase, fructose-6-phosphate kinase (PFK), and NAD-GAPDH were least affected. The transcript levels of PGK and NADP-GAPDH decreased and that of glucose-6-phosphate isomerase (GPI), PGM, and PFK increased in response to water-deficit stress while water-deficit stress had no effect on the steady-state transcript levels of hexokinase, aldolase, TPI, and NAD-GAPDH. The results suggest that under water-deficit stress, the activities and transcript levels of most of the glycolytic enzymes are not significantly affected, except the increased activity and transcript level of PGM and decreased activities and transcript levels of PGK and NADP-GAPDH. Further, the glycolytic enzymes do not show much variation between the tolerant and sensitive cultivars under water deficit.

Genetic diversity and geographical differentiation of Iranian landrace, cultivars, and exotic chickpea lines as revealed by morphological and microsatellite markers

Assessment of the extent of genetic variability within chickpea is fundamental for chickpea breeding and conservation of genetic resources and is particularly useful as a general guide in the choice of parents for breeding hybrids. To establish genetic diversity among 60 accessions of chickpea comprising landraces, internationally developed improved lines, and cultivars, genetic distances were evaluated using 14 simple sequence repeat markers. These markers showed a high level of polymorphism; a total of 59 different alleles were detected, with a mean of 4.2 alleles per locus. The polymorphic information content (PIC) value ranged from 0.31 to 0.89. All the markers, with the exception of TAA170, TA110, GA34, and Ts35, were considered to be informative (PIC > 0.5), indicating their potential usefulness for cultivar identification. Based on the UNJ clustering method, all accessions were clustered in five groups, which indicated the probable origin and region similarity of Iranian landraces over the other cultivars. It also represents a wide diversity among available germplasm. The result has firmly established that introduction of genetic materials from exotic sources has broadened the genetic base of the national chickpea breeding program. As further implications of the findings, this study can be useful for selective breeding for specific traits and in enhancing the genetic base of breeding programs.

Molecular characterization of primary gene pool of chickpea based on ISSR markers

Genetic diversity and relationships within and among members of the primary gene pool of chickpea, including 38 accessions of Cicer arietinum, six of C. reticulatum,, and four of C. echinospermum, were investigated using 31 ISSR markers. The study revealed moderate diversity, detecting 141 fragments, of which 79 (56%) were polymorphic. Averages were 0.125 for polymorphic information content, 0.350 for marker index, and 0.715 for resolving power. The UPGMA dendrogram and the principal coordinate analysis revealed a clear differentiation between wild and cultivated accessions. The clustering pattern did not strictly follow the grouping of accessions by geographic origin but was in good agreement with the pedigree data and the seed type. The study demonstrates that ISSRs provide promising marker tools in revealing genetic diversity and relationships in chickpea and can contribute to efficient identification, conservation, and utilization of germplasm for plant improvement through conventional as well as molecular breeding approaches.

Achievements and prospects of genomics-assisted breeding in three legume crops of the semi-arid tropics

Advances in next-generation sequencing and genotyping technologies have enabled generation of large-scale genomic resources such as molecular markers, transcript reads and BAC-end sequences (BESs) in chickpea, pigeonpea and groundnut, three major legume crops of the semi-arid tropics. Comprehensive transcriptome assemblies and genome sequences have either been developed or underway in these crops. Based on these resources, dense genetic maps, QTL maps as well as physical maps for these legume species have also been developed. As a result, these crops have graduated from 'orphan' or 'less-studied' crops to 'genomic resources rich' crops. This article summarizes the above-mentioned advances in genomics and genomics-assisted breeding applications in the form of marker-assisted selection (MAS) for hybrid purity assessment in pigeonpea; marker-assisted backcrossing (MABC) for introgressing QTL region for drought-tolerance related traits, Fusarium wilt (FW) resistance and Ascochyta blight (AB) resistance in chickpea; late leaf spot (LLS), leaf rust and nematode resistance in groundnut. We critically present the case of use of other modern breeding approaches like marker-assisted recurrent selection (MARS) and genomic selection (GS) to utilize the full potential of genomics-assisted breeding for developing superior cultivars with enhanced tolerance to various environmental stresses. In addition, this article recommends the use of advanced-backcross (AB-backcross) breeding and development of specialized populations such as multi-parents advanced generation intercross (MAGIC) for creating new variations that will help in developing superior lines with broadened genetic base. In summary, we propose the use of integrated genomics and breeding approach in these legume crops to enhance crop productivity in marginal environments ensuring food security in developing countries.