Multi-environment QTL analyses for drought-related traits in a recombinant inbred population of chickpea (Cicer arietinum L.)

Genetic mapping of quantitative trait loci associated with drought tolerance in chickpea (Cicer arietinum L.)

Elucidation of the genetic basis of drought tolerance is vital for genomics-assisted breeding of drought tolerant crop varieties. Here, we used genotyping-by-sequencing (GBS) to identify single nucleotide polymorphisms (SNPs) in recombinant inbred lines (RILs) derived from a cross between a drought tolerant chickpea variety, Pusa 362 and a drought sensitive variety, SBD 377. The GBS identified a total of 35,502 SNPs and subsequent filtering of these resulted in 3237 high-quality SNPs included in the eight linkage groups. Fifty-one percent of these SNPs were located in the genic regions distributed throughout the genome. The high density linkage map has total map length of 1069 cm with an average marker interval of 0.33 cm. The linkage map was used to identify 9 robust and consistent QTLs for four drought related traits viz. membrane stability index, relative water content, seed weight and yield under drought, with percent variance explained within the range of 6.29%-90.68% and LOD scores of 2.64 to 6.38, which were located on five of the eight linkage groups. A genomic region on LG 7 harbors quantitative trait loci (QTLs) explaining > 90% phenotypic variance for membrane stability index, and > 10% PVE for yield. This study also provides the first report of major QTLs for physiological traits such as membrane stability index and relative water content for drought stress in chickpea. A total of 369 putative candidate genes were identified in the 6.6 Mb genomic region spanning these QTLs. In-silico expression profiling based on the available transcriptome data revealed that 326 of these genes were differentially expressed under drought stress. KEGG analysis resulted in reduction of candidate genes from 369 to 99, revealing enrichment in various signaling pathways. Haplotype analysis confirmed 5 QTLs among the initially identified 9 QTLs. Two QTLs, qRWC1.1 and qYLD7.1, were chosen based on high SNP density. Candidate gene-based analysis revealed distinct haplotypes in qYLD7.1 associated with significant phenotypic differences, potentially linked to pathways for secondary metabolite biosynthesis. These identified candidate genes bolster defenses through flavonoids and phenylalanine-derived compounds, aiding UV protection, pathogen resistance, and plant structure.The study provides novel genomic regions and candidate genes which can be utilized in genomics-assisted breeding of superior drought tolerant chickpea cultivars.

Identification of salt stress-tolerant candidate genes in the BC2F2 population at the seedling stages of G. hirsutum and G. darwinii using NGS-based bulked segregant analysis

Salinity is a major threat to the yield and productivity of cotton seedlings. In the present study, we developed a BC₂F₂ population of cotton plants from Gossypium darwinii (5-7) and Gossypium hirsutum (CCRI 12-4) salt-susceptible parents to identify salt-resistant candidate genes. The Illumina HiSeq™ strategy was used with bulked segregant analysis. Salt-resistant and salt-susceptible DNA bulks were pooled by using 30 plants from a BC₂F₂ population. Next-generation sequencing (NGS) technology was used for the sequencing of parents and both bulks. Four significant genomic regions were identified: the first genomic region was located on chromosome 18 (1.86 Mb), the second and third genomic regions were on chromosome 25 (1.06 Mb and 1.94 Mb, respectively), and the fourth was on chromosome 8 (1.41 Mb). The reads of bulk1 and bulk2 were aligned to the G. darwinii and G. hirsutum genomes, respectively, leading to the identification of 20,664,007 single-nucleotide polymorphisms (SNPs) and insertions/deletions (indels). After the screening, 6,573 polymorphic markers were obtained after filtration of the candidate regions. The SNP indices in resistant and susceptible bulks and Δ(SNP-index) values of resistant and susceptible bulks were measured. Based on the higher Δ(SNP-index) value, six effective polymorphic SNPs were selected in a different chromosome. Six effective SNPs were linked to five candidate genes in four genomic regions. Further validation of these five candidate genes was carried out using reverse transcription-quantitative polymerase chain reaction (RT-qPCR), resulting in an expression profile that showed two highly upregulated genes in the salt-tolerant species G. darwinii, i.e., Gohir.D05G367800 and Gohir.D12G239100; however, the opposite was shown in G. hirsutum, for which all genes, except one, showed partial expression. The results indicated that Gohir.D05G367800 and Gohir.D12G239100 may be salt-tolerant genes. We are confident that this study could be helpful for the cloning, transformation, and development of salt-resistant cotton varieties.

Molecular Breeding and Drought Tolerance in Chickpea

Cicer arietinum L. is the third greatest widely planted imperative pulse crop worldwide, and it belongs to the Leguminosae family. Drought is the utmost common abiotic factor on plants, distressing their water status and limiting their growth and development. Chickpea genotypes have the natural ability to fight drought stress using certain strategies viz., escape, avoidance and tolerance. Assorted breeding methods, including hybridization, mutation, and marker-aided breeding, genome sequencing along with omics approaches, could be used to improve the chickpea germplasm lines(s) against drought stress. Root features, for instance depth and root biomass, have been recognized as the greatest beneficial morphological factors for managing terminal drought tolerance in the chickpea. Marker-aided selection, for example, is a genomics-assisted breeding (GAB) strategy that can considerably increase crop breeding accuracy and competence. These breeding technologies, notably marker-assisted breeding, omics, and plant physiology knowledge, underlined the importance of chickpea breeding and can be used in future crop improvement programmes to generate drought-tolerant cultivars(s).

Genetic mapping of QTLs for drought tolerance in chickpea (Cicer arietinum L.)

Chickpea yield is severely affected by drought stress, which is a complex quantitative trait regulated by multiple small-effect genes. Identifying genomic regions associated with drought tolerance component traits may increase our understanding of drought tolerance mechanisms and assist in the development of drought-tolerant varieties. Here, a total of 187 F8 recombinant inbred lines (RILs) developed from an interspecific cross between drought-tolerant genotype GPF 2 (Cicer arietinum) and drought-sensitive accession ILWC 292 (C. reticulatum) were evaluated to identify quantitative trait loci (QTLs) associated with drought tolerance component traits. A total of 21 traits, including 12 morpho-physiological traits and nine root-related traits, were studied under rainfed and irrigated conditions. Composite interval mapping identified 31 QTLs at Ludhiana and 23 QTLs at Faridkot locations for morphological and physiological traits, and seven QTLs were identified for root-related traits. QTL analysis identified eight consensus QTLs for six traits and five QTL clusters containing QTLs for multiple traits on linkage groups CaLG04 and CaLG06. The identified major QTLs and genomic regions associated with drought tolerance component traits can be introgressed into elite cultivars using genomics-assisted breeding to enhance drought tolerance in chickpea.

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.

Construction of a high-density genetic map and QTL analysis for yield, yield components and agronomic traits in chickpea (Cicer arietinum L.)

Unravelling the genetic architecture underlying yield components and agronomic traits is important for enhancing crop productivity. Here, a recombinant inbred line (RIL) population, developed from ICC 4958 and DCP 92–3 cross, was used for constructing linkage map and QTL mapping analysis. The RIL population was genotyped using a high-throughput Axiom^®CicerSNP array, which enabled the development of a high-density genetic map consisting of 3,818 SNP markers and spanning a distance of 1064.14 cM. Analysis of phenotyping data for yield, yield components and agronomic traits measured across three years together with genetic mapping data led to the identification of 10 major-effect QTLs and six minor-effect QTLs explaining up to 59.70% phenotypic variance. The major-effect QTLs identified for 100-seed weight, and plant height possessed key genes, such as C3HC4 RING finger protein, pentatricopeptide repeat (PPR) protein, sugar transporter, leucine zipper protein and NADH dehydrogenase, amongst others. The gene ontology studies highlighted the role of these genes in regulating seed weight and plant height in crop plants. The identified genomic regions for yield, yield components, and agronomic traits, and the closely linked markers will help advance genetics research and breeding programs in chickpea.

Introgression of "QTL-hotspot" region enhances drought tolerance and grain yield in three elite chickpea cultivars

With an aim of enhancing drought tolerance using a marker-assisted backcrossing (MABC) approach, we introgressed the "QTL-hotspot" region from ICC 4958 accession that harbors quantitative trait loci (QTLs) for several drought-tolerance related traits into three elite Indian chickpea (Cicer arietinum L.) cultivars: Pusa 372, Pusa 362, and DCP 92-3. Of eight simple sequence repeat (SSR) markers in the QTL-hotspot region, two to three polymorphic markers were used for foreground selection with respective cross-combinations. A total of 47, 53, and 46 SSRs were used for background selection in case of introgression lines (ILs) developed in genetic backgrounds of Pusa 372, Pusa 362, and DCP 92-3, respectively. In total, 61 ILs (20 BC₃ F₃ in Pusa 372; 20 BC₂ F₃ in Pusa 362, and 21 BC₃ F₃ in DCP 92-3), with >90% recurrent parent genome recovery were developed. Six improved lines in different genetic backgrounds (e.g. BGM 10216 in Pusa 372; BG 3097 and BG 4005 in Pusa 362; IPC(L4-14), IPC(L4-16), and IPC(L19-1) in DCP 92-3) showed better performance than their respective recurrent parents. BGM 10216, with 16% yield gain over Pusa 372, has been released as Pusa Chickpea 10216 by the Central Sub-Committees on Crop Standards, Notification and Release of Varieties of Agricultural Crops, Ministry of Agriculture and Farmers Welfare, Government of India, for commercial cultivation in India. In summary, this study reports introgression of the QTL-hotspot for enhancing yield under rainfed conditions, development of several introgression lines, and release of Pusa Chickpea 10216 developed through molecular breeding in India.

Developing Climate-Resilient Chickpea Involving Physiological and Molecular Approaches With a Focus on Temperature and Drought Stresses

Chickpea is one of the most economically important food legumes, and a significant source of proteins. It is cultivated in more than 50 countries across Asia, Africa, Europe, Australia, North America, and South America. Chickpea production is limited by various abiotic stresses (cold, heat, drought, salt, etc.). Being a winter-season crop in northern south Asia and some parts of the Australia, chickpea faces low-temperature stress (0-15°C) during the reproductive stage that causes substantial loss of flowers, and thus pods, to inhibit its yield potential by 30-40%. The winter-sown chickpea in the Mediterranean, however, faces cold stress at vegetative stage. In late-sown environments, chickpea faces high-temperature stress during reproductive and pod filling stages, causing considerable yield losses. Both the low and the high temperatures reduce pollen viability, pollen germination on the stigma, and pollen tube growth resulting in poor pod set. Chickpea also experiences drought stress at various growth stages; terminal drought, along with heat stress at flowering and seed filling can reduce yields by 40-45%. In southern Australia and northern regions of south Asia, lack of chilling tolerance in cultivars delays flowering and pod set, and the crop is usually exposed to terminal drought. The incidences of temperature extremes (cold and heat) as well as inconsistent rainfall patterns are expected to increase in near future owing to climate change thereby necessitating the development of stress-tolerant and climate-resilient chickpea cultivars having region specific traits, which perform well under drought, heat, and/or low-temperature stress. Different approaches, such as genetic variability, genomic selection, molecular markers involving quantitative trait loci (QTLs), whole genome sequencing, and transcriptomics analysis have been exploited to improve chickpea production in extreme environments. Biotechnological tools have broadened our understanding of genetic basis as well as plants' responses to abiotic stresses in chickpea, and have opened opportunities to develop stress tolerant chickpea.

Altered Expression of an FT Cluster Underlies a Major Locus Controlling Domestication-Related Changes to Chickpea Phenology and Growth Habit

Flowering time is a key trait in breeding and crop evolution, due to its importance for adaptation to different environments and for yield. In the particular case of chickpea, selection for early phenology was essential for the successful transition of this species from a winter to a summer crop. Here, we used genetic and expression analyses in two different inbred populations to examine the genetic control of domestication-related differences in flowering time and growth habit between domesticated chickpea and its wild progenitor Cicer reticulatum. A single major quantitative trait locus for flowering time under short-day conditions [Days To Flower (DTF)3A] was mapped to a 59-gene interval on chromosome three containing a cluster of three FT genes, which collectively showed upregulated expression in domesticated relative to wild parent lines. An equally strong association with growth habit suggests a pleiotropic effect of the region on both traits. These results indicate the likely molecular explanation for the characteristic early flowering of domesticated chickpea, and the previously described growth habit locus Hg. More generally, they point to de-repression of this specific gene cluster as a conserved mechanism for achieving adaptive early phenology in temperate legumes.

Is Chickpea a Potential Substitute for Soybean? Phenolic Bioactives and Potential Health Benefits

Legume seeds are rich sources of protein, fiber, and minerals. In addition, their phenolic compounds as secondary metabolites render health benefits beyond basic nutrition. Lowering apolipoprotein B secretion from HepG2 cells and decreasing the level of low-density lipoprotein (LDL)-cholesterol oxidation are mechanisms related to the prevention of cardiovascular diseases (CVD). Likewise, low-level chronic inflammation and related disorders of the immune system are clinical predictors of cardiovascular pathology. Furthermore, DNA-damage signaling and repair are crucial pathways to the etiology of human cancers. Along CVD and cancer, the prevalence of obesity and diabetes is constantly increasing. Screening the ability of polyphenols in inactivating digestive enzymes is a good option in pre-clinical studies. In addition, in vivo studies support the role of polyphenols in the prevention and/or management of diabetes and obesity. Soybean, a well-recognized source of phenolic isoflavones, exerts health benefits by decreasing oxidative stress and inflammation related to the above-mentioned chronic ailments. Similar to soybeans, chickpeas are good sources of nutrients and phenolic compounds, especially isoflavones. This review summarizes the potential of chickpea as a substitute for soybean in terms of health beneficial outcomes. Therefore, this contribution may guide the industry in manufacturing functional foods and/or ingredients by using an undervalued feedstock.

Research Progress and Perspective on Drought Stress in Legumes: A Review

Climate change, food shortage, water scarcity, and population growth are some of the threatening challenges being faced in today's world. Drought stress (DS) poses a constant challenge for agricultural crops and has been considered a severe constraint for global agricultural productivity; its intensity and severity are predicted to increase in the near future. Legumes demonstrate high sensitivity to DS, especially at vegetative and reproductive stages. They are mostly grown in the dry areas and are moderately drought tolerant, but severe DS leads to remarkable production losses. The most prominent effects of DS are reduced germination, stunted growth, serious damage to the photosynthetic apparatus, decrease in net photosynthesis, and a reduction in nutrient uptake. To curb the catastrophic effect of DS in legumes, it is imperative to understand its effects, mechanisms, and the agronomic and genetic basis of drought for sustainable management. This review highlights the impact of DS on legumes, mechanisms, and proposes appropriate management approaches to alleviate the severity of water stress. In our discussion, we outline the influence of water stress on physiological aspects (such as germination, photosynthesis, water and nutrient uptake), growth parameters and yield. Additionally, mechanisms, various management strategies, for instance, agronomic practices (planting time and geometry, nutrient management), plant growth-promoting Rhizobacteria and arbuscular mycorrhizal fungal inoculation, quantitative trait loci (QTLs), functional genomics and advanced strategies (CRISPR-Cas9) are also critically discussed. We propose that the integration of several approaches such as agronomic and biotechnological strategies as well as advanced genome editing tools is needed to develop drought-tolerant legume cultivars.

Mapping resistance to Phytophthora root rot identifies independent loci from cultivated (Cicer arietinum L.) and wild (Cicer echinospermum P.H. Davis) chickpea

Major QTL for Phytophthora root rot resistance have been identified in three mapping populations with independent sources of resistance contributed by C. echinospermum and C. arietinum. Phytophthora root rot (PRR) caused by the oomycete Phytophthora medicaginis is a major soil-borne disease of chickpea in Australia. With no economic in-crop control of PRR, a genetic approach to improve resistance is the most practical management option. Moderate field resistance has been incorporated in the cultivated C. arietinum variety, Yorker, and a higher level of resistance has been identified in a derivative of wild chickpea (C. echinospermum, interspecific breeding line 04067-81-2-1-1). These genotypes and two other released varieties were used to develop one intra-specific and two interspecific F₆-derived recombinant inbred line mapping populations for genetic analysis of resistance. The Yorker × Genesis114 (YG), Rupali × 04067-81-2-1-1 (RB) and Yorker × 04067-81-2-1-1 (YB) populations were genotyped using genotyping-by-sequencing and phenotyped for PRR under three field environments with a mixture of 10 P. medicaginis isolates. Whole-genome QTL analysis identified major QTL QRBprrsi01, QYBprrsi01, QRBprrsi03 and QYBprrsi02 for PRR resistance on chromosomes 3 and 6, in RB and YB populations, respectively, with the resistance source derived from the wild Cicer species. QTL QYGprrsi02 and QYGprrsi03 were also identified on chromosomes 5 and 6 in YG population from C. arietinum. Aligning QTL regions to the corresponding chickpea reference genome suggested that the resistance source from C. arietinum and C. echinospermum may be different. The findings from this study provide tools for marker-assisted selection in chickpea breeding and information to assist the development of populations suitable for fine-mapping of resistance loci to determine the molecular basis for PRR resistance in chickpea.

Analysis of QTL-allele system conferring drought tolerance at seedling stage in a nested association mapping population of soybean [Glycine max (L.) Merr.] using a novel GWAS procedure

MAIN CONCLUSION

RTM-GWAS identified 111 DT QTLs, 262 alleles with high proportion of QEI and genetic variation accounting for 88.55-95.92% PV in NAM, from which QTL-allele matrices were established and candidate genes annotated. Drought tolerance (DT) is one of the major challenges for world soybean production. A nested association mapping (NAM) population with 403 lines comprising two recombinant inbred line (RIL) populations: M8206 × TongShan and ZhengYang × M8206 was tested for DT using polyethylene-glycol (PEG) treatment under spring and summer environments. The population was sequenced using restriction-site-associated DNA sequencing (RAD-seq) filtered with minor allele frequency (MAF) ≥ 0.01, 55,936 single nucleotide polymorphisms (SNPs) were obtained and organized into 6137 SNP linkage disequilibrium blocks (SNPLDBs). The restricted two-stage multi-locus genome-wide association studies (RTM-GWAS) identified 73 and 38 QTLs with 174 and 88 alleles contributed main effect 40.43 and 26.11% to phenotypic variance (PV) and QTL-environment interaction (QEI) effect 24.64 and 10.35% to PV for relative root length (RRL) and relative shoot length (RSL), respectively. The DT traits were characterized with high proportion of QEI variation (37.52-41.65%), plus genetic variation (46.90-58.40%) in a total of 88.55-95.92% PV. The identified QTLs-alleles were organized into main-effect and QEI-effect QTL-allele matrices, showing the genetic and QEI architecture of the three parents/NAM population. From the matrices, the possible best genotype was predicted to have a weighted average value over two indicators (WAV) of 1.873, while the top ten optimal crosses among RILs with 95th percentile WAV 1.098-1.132, transgressive over the parents (0.651-0.773) but much less than 1.873, implying further pyramiding potential. From the matrices, 134 candidate genes were annotated involved in nine biological processes. The present results provide a novel way for molecular breeding in QTL-allele-based genomic selection for optimal cross selection.

Investigating Drought Tolerance in Chickpea Using Genome-Wide Association Mapping and Genomic Selection Based on Whole-Genome Resequencing Data

Drought tolerance is a complex trait that involves numerous genes. Identifying key causal genes or linked molecular markers can facilitate the fast development of drought tolerant varieties. Using a whole-genome resequencing approach, we sequenced 132 chickpea varieties and advanced breeding lines and found more than 144,000 single nucleotide polymorphisms (SNPs). We measured 13 yield and yield-related traits in three drought-prone environments of Western Australia. The genotypic effects were significant for all traits, and many traits showed highly significant correlations, ranging from 0.83 between grain yield and biomass to -0.67 between seed weight and seed emergence rate. To identify candidate genes, the SNP and trait data were incorporated into the SUPER genome-wide association study (GWAS) model, a modified version of the linear mixed model. We found that several SNPs from auxin-related genes, including auxin efflux carrier protein (PIN3), p-glycoprotein, and nodulin MtN21/EamA-like transporter, were significantly associated with yield and yield-related traits under drought-prone environments. We identified four genetic regions containing SNPs significantly associated with several different traits, which was an indication of pleiotropic effects. We also investigated the possibility of incorporating the GWAS results into a genomic selection (GS) model, which is another approach to deal with complex traits. Compared to using all SNPs, application of the GS model using subsets of SNPs significantly associated with the traits under investigation increased the prediction accuracies of three yield and yield-related traits by more than twofold. This has important implication for implementing GS in plant breeding programs.

Improving Salt Tolerance of Chickpea Using Modern Genomics Tools and Molecular Breeding

INTRODUCTION

The high protein value, essential minerals, dietary fibre and notable ability to fix atmospheric nitrogen make chickpea a highly remunerative crop, particularly in low-input food production systems. Of the variety of constraints challenging chickpea productivity worldwide, salinity remains of prime concern owing to the intrinsic sensitivity of the crop. In view of the projected expansion of chickpea into arable and salt-stressed land by 2050, increasing attention is being placed on improving the salt tolerance of this crop. Considerable effort is currently underway to address salinity stress and substantial breeding progress is being made despite the seemingly highly-complex and environment-dependent nature of the tolerance trait.

CONCLUSION

This review aims to provide a holistic view of recent advances in breeding chickpea for salt tolerance. Initially, we focus on the identification of novel genetic resources for salt tolerance via extensive germplasm screening. We then expand on the use of genome-wide and cost-effective techniques to gain new insights into the genetic control of salt tolerance, including the responsive genes/QTL(s), gene(s) networks/cross talk and intricate signalling cascades.

QTL-seq for rapid identification of candidate genes for 100-seed weight and root/total plant dry weight ratio under rainfed conditions in chickpea

Terminal drought is a major constraint to chickpea productivity. Two component traits responsible for reduction in yield under drought stress include reduction in seeds size and root length/root density. QTL-seq approach, therefore, was used to identify candidate genomic regions for 100-seed weight (100SDW) and total dry root weight to total plant dry weight ratio (RTR) under rainfed conditions. Genomewide SNP profiling of extreme phenotypic bulks from the ICC 4958 × ICC 1882 population identified two significant genomic regions, one on CaLG01 (1.08 Mb) and another on CaLG04 (2.7 Mb) linkage groups for 100SDW. Similarly, one significant genomic region on CaLG04 (1.10 Mb) was identified for RTR. Comprehensive analysis revealed four and five putative candidate genes associated with 100SDW and RTR, respectively. Subsequently, two genes (Ca_04364 and Ca_04607) for 100SDW and one gene (Ca_04586) for RTR were validated using CAPS/dCAPS markers. Identified candidate genomic regions and genes may be useful for molecular breeding for chickpea improvement.

Identification of candidate genes and natural allelic variants for QTLs governing plant height in chickpea

In the present study, molecular mapping of high-resolution plant height QTLs was performed by integrating 3625 desi genome-derived GBS (genotyping-by-sequencing)-SNPs on an ultra-high resolution intra-specific chickpea genetic linkage map (dwarf/semi-dwarf desi cv. ICC12299 x tall kabuli cv. ICC8261). The identified six major genomic regions harboring six robust QTLs (11.5-21.3 PVE), associated with plant height, were mapped within <0.5 cM average marker intervals on six chromosomes. Five SNPs-containing genes tightly linked to the five plant height QTLs, were validated based upon their high potential for target trait association (12.9-20.8 PVE) in 65 desi and kabuli chickpea accessions. The vegetative tissue-specific expression, including higher differential up-regulation (>5-fold) of five genes especially in shoot, young leaf, shoot apical meristem of tall mapping parental accession (ICC8261) as compared to that of dwarf/semi-dwarf parent (ICC12299) was apparent. Overall, combining high-resolution QTL mapping with genetic association analysis and differential expression profiling, delineated natural allelic variants in five candidate genes (encoding cytochrome-c-biosynthesis protein, malic oxidoreductase, NADH dehydrogenase iron-sulfur protein, expressed protein and bZIP transcription factor) regulating plant height in chickpea. These molecular tags have potential to dissect complex plant height trait and accelerate marker-assisted genetic enhancement for developing cultivars with desirable plant height ideotypes in chickpea.

High-density linkage map construction and mapping of seed trait QTLs in chickpea (Cicer arietinum L.) using Genotyping-by-Sequencing (GBS)

This study reports the use of Genotyping-by-Sequencing (GBS) for large-scale SNP discovery and simultaneous genotyping of recombinant inbred lines (RILs) of an intra-specific mapping population of chickpea contrasting for seed traits. A total of 119,672 raw SNPs were discovered, which after stringent filtering revealed 3,977 high quality SNPs of which 39.5% were present in genic regions. Comparative analysis using physically mapped marker loci revealed a higher degree of synteny with Medicago in comparison to soybean. The SNP genotyping data was utilized to construct one of the most saturated intra-specific genetic linkage maps of chickpea having 3,363 mapped positions including 3,228 SNPs on 8 linkage groups spanning 1006.98 cM at an average inter marker distance of 0.33 cM. The map was utilized to identify 20 quantitative trait loci (QTLs) associated with seed traits accounting for phenotypic variations ranging from 9.97% to 29.71%. Analysis of the genomic sequence corresponding to five robust QTLs led to the identification of 684 putative candidate genes whose expression profiling revealed that 101 genes exhibited seed specific expression. The integrated approach utilizing the identified QTLs along with the available genome and transcriptome could serve as a platform for candidate gene identification for molecular breeding of chickpea.

Prioritization of candidate genes in "QTL-hotspot" region for drought tolerance in chickpea (Cicer arietinum L.)

A combination of two approaches, namely QTL analysis and gene enrichment analysis were used to identify candidate genes in the “QTL-hotspot” region for drought tolerance present on the Ca4 pseudomolecule in chickpea. In the first approach, a high-density bin map was developed using 53,223 single nucleotide polymorphisms (SNPs) identified in the recombinant inbred line (RIL) population of ICC 4958 (drought tolerant) and ICC 1882 (drought sensitive) cross. QTL analysis using recombination bins as markers along with the phenotyping data for 17 drought tolerance related traits obtained over 1–5 seasons and 1–5 locations split the “QTL-hotspot” region into two subregions namely “QTL-hotspot_a” (15 genes) and “QTL-hotspot_b” (11 genes). In the second approach, gene enrichment analysis using significant marker trait associations based on SNPs from the Ca4 pseudomolecule with the above mentioned phenotyping data, and the candidate genes from the refined “QTL-hotspot” region showed enrichment for 23 genes. Twelve genes were found common in both approaches. Functional validation using quantitative real-time PCR (qRT-PCR) indicated four promising candidate genes having functional implications on the effect of “QTL-hotspot” for drought tolerance in chickpea.

Development of gene-based markers for use in construction of the chickpea (Cicer arietinum L.) genetic linkage map and identification of QTLs associated with seed weight and plant height

Seed weight and plant height are important agronomic traits and contribute to seed yield. The objective of this study was to identify QTLs underlying these traits using an intra-specific mapping population of chickpea. A F11 population of 177 recombinant inbred lines derived from a cross between SBD377 (100-seed weight--48 g and plant height--53 cm) and BGD112 (100-seed weight--15 g and plant height--65 cm) was used. A total of 367 novel EST-derived functional markers were developed which included 187 EST-SSRs, 130 potential intron polymorphisms (PIPs) and 50 expressed sequence tag polymorphisms (ESTPs). Along with these, 590 previously published markers including 385 EST-based markers and 205 genomic SSRs were utilized. Of the 957 markers tested for analysis of parental polymorphism between the two parents of the mapping population, 135 (14.64%) were found to be polymorphic. Of these, 131 polymorphic markers could be mapped to the 8 linkage groups. The linkage map had a total length of 1140.54 cM with an average marker density of 8.7 cM. The map was further used for QTL identification using composite interval mapping method (CIM). Two QTLs each for seed weight, qSW-1 and qSW-2 (explaining 11.54 and 19.24% of phenotypic variance, respectively) and plant height, qPH-1 and qPH-2 (explaining 13.98 and 12.17% of phenotypic variance, respectively) were detected. The novel set of genic markers, the intra-specific linkage map and the QTLs identified in the present study will serve as valuable genomic resources in improving the chickpea seed yield using marker-assisted selection (MAS) strategies.

Genetic diversity and association mapping of iron and zinc concentrations in chickpea (Cicer arietinum L.)

Chickpea (Cicer arietinum L.) is the world's second most important pulse crop after common bean. Chickpea has historically been an important daily staple in the diet of millions of people, especially in the developing countries. Current chickpea breeding programs have mainly been directed toward high yield, biotic and abiotic stress resilience that has increased global production, but less attention has been directed toward improving micronutrient concentrations in seeds. In an effort to develop micronutrient-dense chickpea lines, a study to examine the variability and to identify SNP alleles associated with seed iron and zinc concentrations was conducted using 94 diverse accessions of chickpea. The results indicated that there is substantial variability present in chickpea germplasm for seed iron and zinc concentrations. In the current set of germplasm, zinc is negatively correlated with grain yield across all locations and years; whereas the negative correlation between iron and grain yield was only significant at the Elrose locality. Eight SNP loci associated with iron and (or) zinc concentrations in chickpea seeds were identified. One SNP located on chromosome 1 (chr1) is associated with both iron and zinc concentrations. On chr4, three SNPs associated with zinc concentration and two SNPs for iron concentration were identified. Two additional SNP loci, one on chr6 and the other on chr7, were also found to be associated with iron and zinc concentrations, respectively. The results show potential opportunity for molecular breeding for improvement of seed iron and zinc concentrations in chickpea.

Genotyping-by-sequencing based intra-specific genetic map refines a ''QTL-hotspot" region for drought tolerance in chickpea

To enhance the marker density in the “QTL-hotspot” region, harboring several QTLs for drought tolerance-related traits identified on linkage group 04 (CaLG04) in chickpea recombinant inbred line (RIL) mapping population ICC 4958 × ICC 1882, a genotyping-by-sequencing approach was adopted. In total, 6.24 Gb data from ICC 4958, 5.65 Gb data from ICC 1882 and 59.03 Gb data from RILs were generated, which identified 828 novel single-nucleotide polymorphisms (SNPs) for genetic mapping. Together with these new markers, a high-density intra-specific genetic map was developed that comprised 1,007 marker loci spanning a distance of 727.29 cM. QTL analysis using the extended genetic map along with precise phenotyping data for 20 traits collected over one to seven seasons identified 49 SNP markers in the “QTL-hotspot” region. These efforts have refined the “QTL-hotspot” region to 14 cM. In total, 164 main-effect QTLs including 24 novel QTLs were identified. In addition, 49 SNPs integrated in the “QTL-hotspot” region were converted into cleaved amplified polymorphic sequence (CAPS) and derived CAPS (dCAPS) markers which can be used in marker-assisted breeding.

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.

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.