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

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

Delineation of genes for a major QTL governing heat stress tolerance in chickpea.

Chickpea (Cicer arietinum) is a cool season grain legume experiencing severe yield loss during heat stress due to the intensifying climate changes and its associated gradual increase of mean temperature. Hence, understanding the genetic architecture regulating heat stress tolerance has emerged as an important trait to be addressed for enhancing yield and productivity of chickpea under heat stress. The present study is intended to identify the major genomic region(s) governing heat stress tolerance in chickpea. For this, an integrated genomics-assisted breeding strategy involving NGS-based high-resolution QTL-seq assay, QTL region-specific association analysis and molecular haplotyping was deployed in a population of 206 mapping individuals and a diversity panel of 217 germplasm accessions of chickpea. This combinatorial strategy delineated a major 156.8 kb QTL genomic region, which was subsequently narrowed-down to a functional candidate gene CaHSFA5 and its natural alleles associated strongly with heat stress tolerance in chickpea. Superior natural alleles and haplotypes delineated from the CaHSFA5 gene have functional significance in regulating heat stress tolerance in chickpea. Histochemical staining, interaction studies along with differential expression profiling of CaHSFA5 and ROS scavenging genes suggest a cross talk between CaHSFA5 with ROS homeostasis pertaining to heat stress tolerance in chickpea. Heterologous gene expression followed by heat stress screening further validated the functional significance of CaHSFA5 for heat stress tolerance. The salient outcomes obtained here can have potential to accelerate multiple translational genomic analysis including marker-assisted breeding and gene editing in order to develop high-yielding heat stress tolerant chickpea varieties.

Functional allele of a MATE gene selected during domestication modulates seed color in chickpea.

Seed color is one of the key target traits of domestication and artificial selection in chickpeas due to its implications on consumer preference and market value. The complex seed color trait has been well dissected in several crop species; however, the genetic mechanism underlying seed color variation in chickpea remains poorly understood. Here, we employed an integrated genomics strategy involving QTL mapping, high-density mapping, map-based cloning, association analysis, and molecular haplotyping in an inter-specific RIL mapping population, association panel, wild accessions, and introgression lines (ILs) of Cicer gene pool. This delineated a MATE gene, CaMATE23, encoding a Transparent Testa (TT) and its natural allele (8-bp insertion) and haplotype underlying a major QTL governing seed color on chickpea chromosome 4. Signatures of selective sweep and a strong purifying selection reflected that CaMATE23, especially its 8-bp insertion natural allelic variant, underwent selection during chickpea domestication. Functional investigations revealed that the 8-bp insertion containing the third cis-regulatory RY-motif element in the CaMATE23 promoter is critical for enhanced binding of CaFUSCA3 transcription factor, a key regulator of seed development and flavonoid biosynthesis, thereby affecting CaMATE23 expression and proanthocyanidin (PA) accumulation in the seed coat to impart varied seed color in chickpea. Consequently, overexpression of CaMATE23 in Arabidopsis tt12 mutant partially restored the seed color phenotype to brown pigmentation, ascertaining its functional role in PA accumulation in the seed coat. These findings shed new light on the seed color regulation and evolutionary history, and highlight the transcriptional regulation of CaMATE23 by CaFUSCA3 in modulating seed color in chickpea. The functionally relevant InDel variation, natural allele, and haplotype from CaMATE23 are vital for translational genomic research, including marker-assisted breeding, for developing chickpea cultivars with desirable seed color that appeal to consumers and meet global market demand.

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

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