Development of a dense genetic map and QTL analysis for pod borer Helicoverpa armigera (Hübner) resistance component traits in chickpea (Cicer arietinum L.)

Genetic resources and genes/QTLs for gram pod borer (Helicoverpa armigera Hübner) resistance in chickpea from the Western Himalayas

Helicoverpa armigera (also known as gram pod borer) is a serious threat to chickpea production in the world. A set of 173 chickpea genotypes were evaluated for H. armigera resistance, including mean larval population (MLP), percentage pod damage (PPD), and pest resistance (PR) for 2 consecutive years (year 2020 and 2021). The same core set was also genotyped with 50K Axiom CicerSNP Array. The trait data and 50,000 single nucleotide polymorphism genotypic data were used together to work out marker-trait associations (MTAs) using different genome-wide association studies models. For MLP, a total of 53 MTAs were identified, including 25 MTAs in year 2020 and 28 MTAs in year 2021. A set of three MTAs was found common in both environments. For PPD, two MTAs in year 2020 and five MTAs in year 2021 were identified. A set of two MTAs were common in both environments. Similarly, for PR, only two MTAs common in both environments were identified. Interestingly, a common MTA (Affx_123255526) on chromosome 2 (Ca2) was found to be associated with all the three component traits (MLP, PPD, and PR) of pod borer resistance in chickpea. Further, we report key genes that encode SCAMPs (that facilitates the secretion of defense-related molecules), quinone oxidoreductase (enables the production of reactive oxygen species that promotes diapause of gram pod borer), and NB-LRR proteins that have been implicated in plant defense against H. armigera. The resistant chickpea genotypes, MTAs, and key genes reported in the present study may prove useful in the future for developing pod borer-resistant chickpea varieties.

Field screening and identification of biochemical indices of pod borer (Helicoverpa armigera) resistance in chickpea mutants

Chickpea pod borer (CPB) (Helicoverpa armigera) is one of the major pests, causing significant yield losses. The objectives were to screen chickpea mutants for pod borer resistance/tolerance under field conditions and identification of biochemical markers of tolerance. Chickpea mutant CM216-A/15 had highest leaf (25 trichomes/mm2) and stem trichome density (17 trichomes/mm2) with least pod damage at Kallur Kot and highest pod weight per plant (22.8 ± 2.6g) at AZRI. Higher total phenolic contents (TPCs) and antioxidant capacity were detected in tolerant mutants, i.e., CM216-A/15 and CM664/15. TPC was positively associated with pod yield and had negative correlation with pod damage. Mutants CM216-A/15, CM664/15, and CM766/15 depicted the highest resilience to CPB, owing to higher hairiness, better antioxidant defense response, and lower levels of hydrolytic enzymes and sugars. Identified biochemical markers like TPC, total oxidant status, superoxide dismutase, and pigments can be used for screening of CPB-tolerant/resistant mutants.

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.

A comprehensive investigation of lipid-transfer proteins from Cicer arietinum disentangles their role in plant defense against Helicoverpa armigera-infestation

Lipid Transfer Proteins (LTPs) play a crucial role in synthesizing lipid barrier polymers and are involved in defense signaling during pest and pathogen attacks. Although LTPs are conserved with multifaceted roles in plants, these are not yet identified and characterized in Cicer arietinum. In this study, a genome-wide analysis of LTPs was executed and their physiochemical properties, biochemical function, gene structure analysis, chromosomal localization, promoter analysis, gene duplication, and evolutionary analysis were performed using in silico tools. Furthermore, tissue-specific expression analysis and gene expression analysis during pest attack was also conducted for the LTPs. A total of 48 LTPs were identified and named as CaLTPs. They were predicted to be small unstable proteins with "Glycolipid transfer protein" and "Alpha-Amylase Inhibitors, Lipid Transfer and Seed Storage" domains, that are translocated to the extracellular region. CaLTPs were predicted to possess 3-4 introns and were located on all the eight chromosomes of chickpea with half of the CaLTPs being localized on chromosomes 4, 5, and 6, and found to be closely related to LTPs of Arabidopsis thaliana and Medicago trancatula. Gene duplication and synteny analysis revealed that most of the CaLTPs have evolved due to tandem or segmental gene duplication and were subjected to purifying selection during evolution. The promoters of CaLTPs had development-related, phytohormone-responsive, and abiotic and biotic stress-related cis-acting elements. A few CaLTP transcripts exhibited differential expression in diverse tissue types, while others showed no/very low expression. Out of 20 jasmonate-regulated CaLTPs, 14 exhibited differential expression patterns during Helicoverpa armigera-infestation, indicating their role in plant defense response. This study identified and characterized CaLTPs from an important legume, C. arietinum, and indicated their involvement in plant defense against H. armigera-infestation, which can be further utilized to explore lipid signaling during plant-pest interaction and pest management.

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.

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.

Karyotype Differentiation in Cultivated Chickpea Revealed by Oligopainting Fluorescence in situ Hybridization

Chickpea (Cicer arietinum L.) is one of the main sources of plant proteins in the Indian subcontinent and West Asia, where two different morphotypes, desi and kabuli, are grown. Despite the progress in genome mapping and sequencing, the knowledge of the chickpea genome at the chromosomal level, including the long-range molecular chromosome organization, is limited. Earlier cytogenetic studies in chickpea suffered from a limited number of cytogenetic landmarks and did not permit to identify individual chromosomes in the metaphase spreads or to anchor pseudomolecules to chromosomes in situ. In this study, we developed a system for fast molecular karyotyping for both morphotypes of cultivated chickpea. We demonstrate that even draft genome sequences are adequate to develop oligo-fluorescence in situ hybridization (FISH) barcodes for the identification of chromosomes and comparative analysis among closely related chickpea genotypes. Our results show the potential of oligo-FISH barcoding for the identification of structural changes in chromosomes, which accompanied genome diversification among chickpea cultivars. Moreover, oligo-FISH barcoding in chickpea pointed out some problematic, most probably wrongly assembled regions of the pseudomolecules of both kabuli and desi reference genomes. Thus, oligo-FISH appears as a powerful tool not only for comparative karyotyping but also for the validation of genome assemblies.