Whole genome resequencing and phenotyping of MAGIC population for high resolution mapping of drought tolerance in chickpea.

Terminal drought is one of the major constraints to crop production in chickpea (Cicer arietinum L.). In order to map drought tolerance related traits at high resolution, we sequenced multi-parent advanced generation intercross (MAGIC) population using whole genome resequencing approach and phenotyped it under drought stress environments for two consecutive years (2013-14 and 2014-15). A total of 52.02 billion clean reads containing 4.67 TB clean data were generated on the 1136 MAGIC lines and eight parental lines. Alignment of clean data on to the reference genome enabled identification of a total, 932,172 of SNPs, 35,973 insertions, and 35,726 deletions among the parental lines. A high-density genetic map was constructed using 57,180 SNPs spanning a map distance of 1606.69 cM. Using compressed mixed linear model, genome-wide association study (GWAS) enabled us to identify 737 markers significantly associated with days to 50% flowering, days to maturity, plant height, 100 seed weight, biomass, and harvest index. In addition to the GWAS approach, an identity-by-descent (IBD)-based mixed model approach was used to map quantitative trait loci (QTLs). The IBD-based mixed model approach detected major QTLs that were comparable to those from the GWAS analysis as well as some exclusive QTLs with smaller effects. The candidate genes like FRIGIDA and CaTIFY4b can be used for enhancing drought tolerance in chickpea. The genomic resources, genetic map, marker-trait associations, and QTLs identified in the study are valuable resources for the chickpea community for developing climate resilient chickpeas.

Transcriptome profiling reveals the expression and regulation of genes associated with Fusarium wilt resistance in chickpea (Cicer arietinum L.).

Fusarium wilt (FW) is one of the most significant biotic stresses limiting chickpea production worldwide. To dissect the molecular mechanism of FW resistance in chickpea, comparative transcriptome analyses of contrasting resistance sources of chickpea genotypes under control and Fusarium oxysporum f. sp. ciceris (Foc) inoculated conditions were performed. The high-throughput transcriptome sequencing generated about 1137 million sequencing reads from 24 samples representing two resistant genotypes, two susceptible genotypes, and two near-isogenic lines under control and stress conditions at two-time points (7th- and 12th-day post-inoculation). The analysis identified 5182 differentially expressed genes (DEGs) between different combinations of chickpea genotypes. Functional annotation of these genes indicated their involvement in various biological processes such as defense response, cell wall biogenesis, secondary metabolism, and disease resistance. A significant number (382) of transcription factor encoding genes exhibited differential expression patterns under stress. Further, a considerable number of the identified DEGs (287) co-localized with previously reported quantitative trait locus for FW resistance. Several resistance/susceptibility-related genes, such as SERINE/THREONINE PROTEIN KINASE, DIRIGENT, and MLO exhibiting contrasting expression patterns in resistant and susceptible genotypes upon Foc inoculation, were identified. The results presented in the study provide valuable insights into the transcriptional dynamics associated with FW stress response in chickpea and provide candidate genes for the development of disease-resistant chickpea cultivars.

Comprehensive Transcriptome Profiling Uncovers Molecular Mechanisms and Potential Candidate Genes Associated with Heat Stress Response in Chickpea.

Chickpea (Cicer arietinum L.) production is highly susceptible to heat stress (day/night temperatures above 32/20 °C). Identifying the molecular mechanisms and potential candidate genes underlying heat stress response is important for increasing chickpea productivity. Here, we used an RNA-seq approach to investigate the transcriptome dynamics of 48 samples which include the leaf and root tissues of six contrasting heat stress responsive chickpea genotypes at the vegetative and reproductive stages of plant development. A total of 14,544 unique, differentially expressed genes (DEGs) were identified across different combinations studied. These DEGs were mainly involved in metabolic processes, cell wall remodeling, calcium signaling, and photosynthesis. Pathway analysis revealed the enrichment of metabolic pathways, biosynthesis of secondary metabolites, and plant hormone signal transduction, under heat stress conditions. Furthermore, heat-responsive genes encoding bHLH, ERF, WRKY, and MYB transcription factors were differentially regulated in response to heat stress, and candidate genes underlying the quantitative trait loci (QTLs) for heat tolerance component traits, which showed differential gene expression across tolerant and sensitive genotypes, were identified. Our study provides an important resource for dissecting the role of candidate genes associated with heat stress response and also paves the way for developing climate-resilient chickpea varieties for the future.

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.

Genetic variation in CaTIFY4b contributes to drought adaptation in chickpea.

Chickpea production is vulnerable to drought stress. Identifying the genetic components underlying drought adaptation is crucial for enhancing chickpea productivity. Here, we present the fine mapping and characterization of 'QTL-hotspot', a genomic region controlling chickpea growth with positive consequences on crop production under drought. We report that a non-synonymous substitution in the transcription factor CaTIFY4b regulates seed weight and organ size in chickpea. Ectopic expression of CaTIFY4b in Medicago truncatula enhances root growth under water deficit. Our results suggest that allelic variation in 'QTL-hotspot' improves pre-anthesis water use, transpiration efficiency, root architecture and canopy development, enabling high-yield performance under terminal drought conditions. Gene expression analysis indicated that CaTIFY4b may regulate organ size under water deficit by modulating the expression of GRF-INTERACTING FACTOR1 (GIF1), a transcriptional co-activator of Growth-Regulating Factors. Taken together, our study offers new insights into the role of CaTIFY4b and on diverse physiological and molecular mechanisms underpinning chickpea growth and production under specific drought scenarios.

Translational Chickpea Genomics Consortium to Accelerate Genetic Gains in Chickpea (Cicer arietinum L.).

The Translational Chickpea Genomics Consortium (TCGC) was set up to increase the production and productivity of chickpea (Cicer arietinum L.). It represents research institutes from six major chickpea growing states (Madhya Pradesh, Maharashtra, Andhra Pradesh, Telangana, Karnataka and Uttar Pradesh) of India. The TCGC team has been engaged in deploying modern genomics approaches in breeding and popularizing improved varieties in farmers' fields across the states. Using marker-assisted backcrossing, introgression lines with enhanced drought tolerance and fusarium wilt resistance have been developed in the genetic background of 10 elite varieties of chickpea. Multi-location evaluation of 100 improved lines (70 desi and 30 kabuli) during 2016-2017 and 2018-2019 enabled the identification of top performing desi and kabuli lines. In total, 909 Farmer Participatory Varietal Selection trials were conducted in 158 villages in 16 districts of the five states, during 2017-2018, 2018-2019, and 2019-2020, involving 16 improved varieties. New molecular breeding lines developed in different genetic backgrounds are potential candidates for national trials under the ICAR-All India Coordinated Research Project on Chickpea. The comprehensive efforts of TCGC resulted in the development and adoption of high-yielding varieties that will increase chickpea productivity and the profitability of chickpea growing farmers.

Molecular Mapping of QTLs for Heat Tolerance in Chickpea.

Chickpea (Cicer arietinum L.), a cool-season legume, is increasingly affected by heat-stress at reproductive stage due to changes in global climatic conditions and cropping systems. Identifying quantitative trait loci (QTLs) for heat tolerance may facilitate breeding for heat tolerant varieties. The present study was aimed at identifying QTLs associated with heat tolerance in chickpea using 292 F8-9 recombinant inbred lines (RILs) developed from the cross ICC 4567 (heat sensitive) × ICC 15614 (heat tolerant). Phenotyping of RILs was undertaken for two heat-stress (late sown) and one non-stress (normal sown) environments. A genetic map spanning 529.11 cM and comprising 271 genotyping by sequencing (GBS) based single nucleotide polymorphism (SNP) markers was constructed. Composite interval mapping (CIM) analysis revealed two consistent genomic regions harbouring four QTLs each on CaLG05 and CaLG06. Four major QTLs for number of filled pods per plot (FPod), total number of seeds per plot (TS), grain yield per plot (GY) and % pod setting (%PodSet), located in the CaLG05 genomic region, were found to have cumulative phenotypic variation of above 50%. Nineteen pairs of epistatic QTLs showed significant epistatic effect, and non-significant QTL × environment interaction effect, except for harvest index (HI) and biomass (BM). A total of 25 putative candidate genes for heat-stress were identified in the two major genomic regions. This is the first report on QTLs for heat-stress response in chickpea. The markers linked to the above mentioned four major QTLs can facilitate marker-assisted breeding for heat tolerance in chickpea.

Molecular Mapping of Flowering Time Major Genes and QTLs in Chickpea (Cicer arietinum L.).

Flowering time is an important trait for adaptation and productivity of chickpea in the arid and the semi-arid environments. This study was conducted for molecular mapping of genes/quantitative trait loci (QTLs) controlling flowering time in chickpea using F2 populations derived from four crosses (ICCV 96029 × CDC Frontier, ICC 5810 × CDC Frontier, BGD 132 × CDC Frontier and ICC 16641 × CDC Frontier). Genetic studies revealed monogenic control of flowering time in the crosses ICCV 96029 × CDC Frontier, BGD 132 × CDC Frontier and ICC 16641 × CDC Frontier, while digenic control with complementary gene action in ICC 5810 × CDC Frontier. The intraspecific genetic maps developed from these crosses consisted 75, 75, 68 and 67 markers spanning 248.8 cM, 331.4 cM, 311.1 cM and 385.1 cM, respectively. A consensus map spanning 363.8 cM with 109 loci was constructed by integrating four genetic maps. Major QTLs corresponding to flowering time genes efl-1 from ICCV 96029, efl-3 from BGD 132 and efl-4 from ICC 16641 were mapped on CaLG04, CaLG08 and CaLG06, respectively. The QTLs and linked markers identified in this study can be used in marker-assisted breeding for developing early maturing chickpea.

Recent breeding programs enhanced genetic diversity in both desi and kabuli varieties of chickpea (Cicer arietinum L.).

In order to understand the impact of breeding on genetic diversity and gain insights into temporal trends in diversity in chickpea, a set of 100 chickpea varieties released in 14 countries between 1948 and 2012 were re-sequenced. For analysis, the re-sequencing data for 29 varieties available from an earlier study was also included. Copy number variations and presence absence variations identified in the present study have potential to drive phenotypic variations for trait improvement. Re-sequencing of a large number of varieties has provided opportunities to inspect the genetic and genomic changes reflecting the history of breeding, which we consider as breeding signatures and the selected loci may provide targets for crop improvement. Our study also reports enhanced diversity in both desi and kabuli varieties as a result of recent chickpea breeding efforts. The current study will aid the explicit efforts to breed for local adaptation in the context of anticipated climate changes.

Plant growth-promotion and biofortification of chickpea and pigeonpea through inoculation of biocontrol potential bacteria, isolated from organic soils.

Seven strains of bacteria [Pseudomonas plecoglossicida SRI-156, Brevibacterium antiquum SRI-158, Bacillus altitudinis SRI-178, Enterobacter ludwigii SRI-211, E. ludwigii SRI-229, Acinetobacter tandoii SRI-305 and Pseudomonas monteilii SRI-360; demonstrated previously for control of charcoal rot disease in sorghum and plant growth-promotion (PGP) in rice] were evaluated for their PGP and biofortification traits in chickpea and pigeonpea under field conditions. When treated on seed, the seven selected bacteria significantly enhanced the shoot height and root length of both chickpea and pigeonpea over the un-inoculated control. Under field conditions, in both chickpea and pigeonpea, the plots inoculated with test bacteria enhanced the nodule number, nodule weight, root and shoot weights, pod number, pod weight, leaf weight, leaf area and grain yield over the un-inoculated control plots. Among the seven bacteria, SRI-229 was found to significantly and consistently enhance all the studied PGP and yield traits including nodule number (24 and 36%), nodule weight (11 and 44%), shoot weight (22 and 20%), root weight (23 and 16%) and grain yield (19 and 26%) for both chickpea and pigeonpea, respectively. When the harvested grains were evaluated for their mineral contents, iron (up to 18 and 12%), zinc (up to 23 and 5%), copper (up to 19 and 8%), manganese (up to 2 and 39%) and calcium (up to 22 and 11%) contents in chickpea and pigeonpea, respectively, were found enhanced in test bacteria inoculated plots over the un-inoculated control plots. This study further confirms that the selected bacterial isolates not only have the potential for PGP in cereals and legumes but also have the potential for biofortification of mineral nutrients.

Whole genome re-sequencing reveals genome-wide variations among parental lines of 16 mapping populations in chickpea (Cicer arietinum L.).

BACKGROUND: Chickpea (Cicer arietinum L.) is the second most important grain legume cultivated by resource poor farmers in South Asia and Sub-Saharan Africa. In order to harness the untapped genetic potential available for chickpea improvement, we re-sequenced 35 chickpea genotypes representing parental lines of 16 mapping populations segregating for abiotic (drought, heat, salinity), biotic stresses (Fusarium wilt, Ascochyta blight, Botrytis grey mould, Helicoverpa armigera) and nutritionally important (protein content) traits using whole genome re-sequencing approach. RESULTS: A total of 192.19 Gb data, generated on 35 genotypes of chickpea, comprising 973.13 million reads, with an average sequencing depth of ~10 X for each line. On an average 92.18 % reads from each genotype were aligned to the chickpea reference genome with 82.17 % coverage. A total of 2,058,566 unique single nucleotide polymorphisms (SNPs) and 292,588 Indels were detected while comparing with the reference chickpea genome. Highest number of SNPs were identified on the Ca4 pseudomolecule. In addition, copy number variations (CNVs) such as gene deletions and duplications were identified across the chickpea parental genotypes, which were minimum in PI 489777 (1 gene deletion) and maximum in JG 74 (1,497). A total of 164,856 line specific variations (144,888 SNPs and 19,968 Indels) with the highest percentage were identified in coding regions in ICC 1496 (21 %) followed by ICCV 97105 (12 %). Of 539 miscellaneous variations, 339, 138 and 62 were inter-chromosomal variations (CTX), intra-chromosomal variations (ITX) and inversions (INV) respectively. CONCLUSION: Genome-wide SNPs, Indels, CNVs, PAVs, and miscellaneous variations identified in different mapping populations are a valuable resource in genetic research and helpful in locating genes/genomic segments responsible for economically important traits. Further, the genome-wide variations identified in the present study can be used for developing high density SNP arrays for genetics and breeding applications.

Salt sensitivity in chickpea.

The growth of chickpea (Cicer arietinum L.) is very sensitive to salinity, with the most susceptible genotypes dying in just 25 mm NaCl and resistant genotypes unlikely to survive 100 mm NaCl in hydroponics; germination is more tolerant with some genotypes tolerating 320 mm NaCl. When growing in a saline medium, Cl(-), which is secreted from glandular hairs on leaves, stems and pods, is present in higher concentrations in shoots than Na(+). Salinity reduces the amount of water extractable from soil by a chickpea crop and induces osmotic adjustment, which is greater in nodules than in leaves or roots. Chickpea rhizobia show a higher 'free-living' salt resistance than chickpea plants, and salinity can cause large reductions in nodulation, nodule size and N(2)-fixation capacity. Recent screenings of diverse germplasm suggest significant variation of seed yield under saline conditions. Both dominance and additive gene effects have been identified in the effects of salinity on chickpea and there appears to be sufficient genetic variation to enable improvement in yield under saline conditions via breeding. Selections are required across the entire life cycle with a range of rhizobial strains under salt-affected, preferably field, conditions.