Zero hunger and good health could be realized by 2030 through effective conservation, characterization and utilization of germplasm resources1. So far, few chickpea (Cicerarietinum) germplasm accessions have been characterized at the genome sequence level2. Here we present a detailed map of variation in 3,171 cultivated and 195 wild accessions to provide publicly available resources for chickpea genomics research and breeding. We constructed a chickpea pan-genome to describe genomic diversity across cultivated chickpea and its wild progenitor accessions. A divergence tree using genes present in around 80% of individuals in one species allowed us to estimate the divergence of Cicer over the last 21 million years. Our analysis found chromosomal segments and genes that show signatures of selection during domestication, migration and improvement. The chromosomal locations of deleterious mutations responsible for limited genetic diversity and decreased fitness were identified in elite germplasm. We identified superior haplotypes for improvement-related traits in landraces that can be introgressed into elite breeding lines through haplotype-based breeding, and found targets for purging deleterious alleles through genomics-assisted breeding and/or gene editing. Finally, we propose three crop breeding strategies based on genomic prediction to enhance crop productivity for 16 traits while avoiding the erosion of genetic diversity through optimal contribution selection (OCS)-based pre-breeding. The predicted performance for 100-seed weight, an important yield-related trait, increased by up to 23% and 12% with OCS- and haplotype-based genomic approaches, respectively.

Transmembrane (TMEM)-176A and 176B proteins belong to the MS4A family of proteins whose function in the immune system remains unclear. TMEM176A transcripts were previously shown to be elevated in liver cancer or kidney tissue with proteinuria, while marked changes in TMEM176B transcripts have been found in tolerated tissue allografts and neoplastic fibroblasts. To study the functional relationship between human TMEM176A and 176B and their putative link to cancer, we used polymerase chain reaction and biochemical assays. Here, we show that TMEM176A and 176B are widely expressed in all human tissues examined. Co-immunoprecipitation of heterologously expressed TMEM176A and 176B revealed direct physical interaction. To determine the relevance of such interaction to cancer pathology, we analyzed biopsied tissue samples from a variety of normal and cancer tissues. Our data reveal that human TMEM176A and 176B protein levels are significantly elevated in lymphoma, but not in normal tissues. The protein levels of TMEM176A are also significantly increased in lung carcinoma. Finally, analysis of the protein expression ratio of TMEM176A over 176B showed significant differences between normal and cancer tissues of the breast, lymph, skin, and liver, which indicates that both TMEM proteins could be potential useful markers for certain human cancers.

In the context of climate change, heat stress during the reproductive stages of chickpea (Cicer arietinum L.) leads to significant yield losses. In order to identify the genomic regions responsible for heat stress tolerance, a recombinant inbred line population derived from DCP 92-3 (heat sensitive) and ICCV 92944 (heat tolerant) was genotyped using the genotyping-by-sequencing approach and evaluated for two consecutive years (2017 and 2018) under normal and late sown or heat stress environments. A high-density genetic map comprising 788 single-nucleotide polymorphism markers spanning 1,125 cM was constructed. Using composite interval mapping, a total of 77 QTLs (37 major and 40 minor) were identified for 12 of 13 traits. A genomic region on CaLG07 harbors quantitative trait loci (QTLs) explaining >30% phenotypic variation for days to pod initiation, 100 seed weight, and for nitrogen balance index explaining >10% PVE. In addition, we also reported for the first time major QTLs for proxy traits (physiological traits such as chlorophyll content, nitrogen balance index, normalized difference vegetative index, and cell membrane stability). Furthermore, 32 candidate genes in the QTL regions that encode the heat shock protein genes, heat shock transcription factors, are involved in flowering time regulation as well as pollen-specific genes. The major QTLs reported in this study, after validation, may be useful in molecular breeding for developing heat-tolerant superior lines or varieties.

Biofortification through plant breeding is a cost-effective and sustainable approach towards addressing micronutrient malnutrition prevailing across the globe. Screening cultivars for micronutrient content and identification of quantitative trait loci (QTLs)/genes and markers help in the development of biofortified varieties in chickpea ( Cicer arietinum L.). With the aim of identifying the genomic regions controlling seed Fe and Zn concentrations, the F 2:3 population derived from a cross between MNK-1 and Annigeri 1 was genotyped using genotyping by sequencing approach and evaluated for Fe and Zn concentration. An intraspecific genetic linkage map comprising 839 single nucleotide polymorphisms (SNPs) spanning a total distance of 1,088.04 cM with an average marker density of 1.30 cM was constructed. By integrating the linkage map data with the phenotypic data of the F 2:3 population, a total of 11 QTLs were detected for seed Fe concentration on CaLG03, CaLG04, and CaLG05, with phenotypic variation explained ranging from 7.2% ( CaqFe3.4 ) to 13.4% ( CaqFe4.2 ). For seed Zn concentration, eight QTLs were identified on CaLG04, CaLG05, and CaLG08. The QTLs individually explained phenotypic variations ranging between 5.7% ( CaqZn8.1 ) and 13.7% ( CaqZn4.3 ). Three QTLs for seed Fe and Zn concentrations ( CaqFe4.4, CaqFe4.5 , and CaqZn4.1 ) were colocated in the “ QTL-hotspot ” region on CaLG04 that harbors several drought tolerance-related QTLs. We identified genes in the QTL regions that encode iron–sulfur metabolism and zinc-dependent alcohol dehydrogenase activity on CaLG03, iron ion binding oxidoreductase on CaLG04, and zinc-induced facilitator-like protein and ZIP zinc/iron transport family protein on CaLG05. These genomic regions and the associated markers can be used in marker-assisted selection to increase seed Fe and Zn concentrations in agronomically superior chickpea varieties.

Micronutrient malnutrition is a serious concern in many parts of the world; therefore, enhancing crop nutrient content is an important challenge. Chickpea (Cicer arietinum L.), a major food legume crop worldwide, is a vital source of protein and minerals in the vegetarian diet. This study evaluated a diverse set of 258 chickpea germplasm accessions for 12 key nutritional traits. A significant variation was observed for several nutritional traits, including crude protein (16.56–24.64/100 g), β-Carotene (0.003–0.104 mg/100 g), calcium (60.69–176.55 mg/100 g), and folate (0.413–6.537 mg/kg). These data, combined with the available whole-genome sequencing data for 318,644 SNPs, were used in genome-wide association studies comprising single-locus and multi-locus models. We also explored the effect of varying the minor allele frequency (MAF) levels and heterozygosity. We identified 62 significant marker-trait associations (MTAs) explaining up to 28.63% of the phenotypic variance (PV), of which nine were localized within genes regulating G protein-coupled receptor signaling pathway, proteasome assembly, intracellular signal transduction, and oxidation–reduction process, among others. The significant effect MTAs were located primarily on Ca1, Ca3, Ca4, and Ca6. Importantly, varying the level of heterozygosity was found to significantly affect the detection of associations contributing to traits of interest. We further identified seven promising accessions (ICC10399, ICC1392, ICC1710, ICC2263, ICC1431, ICC4182, and ICC16915) with superior agronomic performance and high nutritional content as potential donors for developing nutrient-rich, high-yielding chickpea varieties. Validation of the significant MTAs with higher PV could identify factors controlling the nutrient acquisition and facilitate the design of biofortified chickpeas for the future.