A chickpea genetic variation map based on the sequencing of 3,366 genomes.

Zero hunger and good health could be realized by 2030 through effective conservation, characterization and utilization of germplasm resources1. So far, few chickpea (Cicer arietinum) 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.

Genetic Analysis of Partially Resistant and Susceptible Chickpea Cultivars in Response to Ascochyta rabiei Infection.

The molecular mechanism involved in chickpea (Cicer arietinum L.) resistance to the necrotrophic fungal pathogen Ascochyta rabiei is not well documented. A. rabiei infection can cause severe damage in chickpea, resulting in significant economic losses. Understanding the resistance mechanism against ascochyta blight can help to define strategies to develop resistant cultivars. In this study, differentially expressed genes from two partially resistant cultivars (CDC Corinne and CDC Luna) and a susceptible cultivar (ICCV 96029) to ascochyta blight were identified in the early stages (24, 48 and 72 h) of A. rabiei infection using RNA-seq. Altogether, 3073 genes were differentially expressed in response to A. rabiei infection across different time points and cultivars. A larger number of differentially expressed genes (DEGs) were found in CDC Corinne and CDC Luna than in ICCV 96029. Various transcription factors including ERF, WRKY, bHLH and MYB were differentially expressed in response to A. rabiei infection. Genes involved in pathogen detection and immune signalings such as receptor-like kinases (RLKs), Leucine-Rich Repeat (LRR)-RLKs, and genes associated with the post-infection defence response were differentially expressed among the cultivars. GO functional enrichment and pathway analysis of the DEGs suggested that the biological processes such as metabolic process, response to stimulus and catalytic activity were overrepresented in both resistant and susceptible chickpea cultivars. The expression patterns of eight randomly selected genes revealed by RNA-seq were confirmed by quantitative PCR (qPCR) analysis. The results provide insights into the complex molecular mechanism of the chickpea defence in response to the A. rabiei infection.

Novel Alleles from Cicer reticulatum L. for Genetic Improvement of Cultivated Chickpeas Identified through Genome Wide Association Analysis.

The availability of wild chickpea (Cicer reticulatum L.) accessions has the potential to be used for the improvement of important traits in cultivated chickpeas. The main objectives of this study were to evaluate the phenotypic and genetic variations of chickpea progeny derived from interspecific crosses between C. arietinum and C. reticulatum, and to establish the association between single nucleotide polymorphism (SNP) markers and a series of important agronomic traits in chickpea. A total of 486 lines derived from interspecific crosses between C. arietinum (CDC Leader) and 20 accessions of C. reticulatum were evaluated at different locations in Saskatchewan, Canada in 2017 and 2018. Significant variations were observed for seed weight per plant, number of seeds per plant, thousand seed weight, and plant biomass. Path coefficient analysis showed significant positive direct effects of the number of seeds per plant, thousand seed weight, and biomass on the total seed weight. Cluster analysis based on the agronomic traits generated six groups that allowed the identification of potential heterotic groups within the interspecific lines for yield improvement and resistance to ascochyta blight disease. Genotyping of the 381 interspecific lines using a modified genotyping by sequencing (tGBS) generated a total of 14,591 SNPs. Neighbour-joining cluster analysis using the SNP data grouped the lines into 20 clusters. The genome wide association analysis identified 51 SNPs that had significant associations with different traits. Several candidate genes associated with early flowering and yield components were identified. The candidate genes and the significant SNP markers associated with different traits have a potential to aid the trait introgression in the breeding program.

Integrating genomics for chickpea improvement: achievements and opportunities.

KEY MESSAGE: Integration of genomic technologies with breeding efforts have been used in recent years for chickpea improvement. Modern breeding along with low cost genotyping platforms have potential to further accelerate chickpea improvement efforts. The implementation of novel breeding technologies is expected to contribute substantial improvements in crop productivity. While conventional breeding methods have led to development of more than 200 improved chickpea varieties in the past, still there is ample scope to increase productivity. It is predicted that integration of modern genomic resources with conventional breeding efforts will help in the delivery of climate-resilient chickpea varieties in comparatively less time. Recent advances in genomics tools and technologies have facilitated the generation of large-scale sequencing and genotyping data sets in chickpea. Combined analysis of high-resolution phenotypic and genetic data is paving the way for identifying genes and biological pathways associated with breeding-related traits. Genomics technologies have been used to develop diagnostic markers for use in marker-assisted backcrossing programmes, which have yielded several molecular breeding products in chickpea. We anticipate that a sequence-based holistic breeding approach, including the integration of functional omics, parental selection, forward breeding and genome-wide selection, will bring a paradigm shift in development of superior chickpea varieties. There is a need to integrate the knowledge generated by modern genomics technologies with molecular breeding efforts to bridge the genome-to-phenome gap. Here, we review recent advances that have led to new possibilities for developing and screening breeding populations, and provide strategies for enhancing the selection efficiency and accelerating the rate of genetic gain in chickpea.

QTL sequencing strategy to map genomic regions associated with resistance to ascochyta blight in chickpea.

Whole-genome sequencing-based bulked segregant analysis (BSA) for mapping quantitative trait loci (QTL) provides an efficient alternative approach to conventional QTL analysis as it significantly reduces the scale and cost of analysis with comparable power to QTL detection using full mapping population. We tested the application of next-generation sequencing (NGS)-based BSA approach for mapping QTLs for ascochyta blight resistance in chickpea using two recombinant inbred line populations CPR-01 and CPR-02. Eleven QTLs in CPR-01 and six QTLs in CPR-02 populations were mapped on chromosomes Ca1, Ca2, Ca4, Ca6 and Ca7. The QTLs identified in CPR-01 using conventional biparental mapping approach were used to compare the efficiency of NGS-based BSA in detecting QTLs for ascochyta blight resistance. The QTLs on chromosomes Ca1, Ca4, Ca6 and Ca7 overlapped with the QTLs previously detected in CPR-01 using conventional QTL mapping method. The QTLs on chromosome Ca4 were detected in both populations and overlapped with the previously reported QTLs indicating conserved region for ascochyta blight resistance across different chickpea genotypes. Six candidate genes in the QTL regions identified using NGS-based BSA on chromosomes Ca2 and Ca4 were validated for their association with ascochyta blight resistance in the CPR-02 population. This study demonstrated the efficiency of NGS-based BSA as a rapid and cost-effective method to identify QTLs associated with ascochyta blight in chickpea.

Genome-wide SNP discovery for development of high-density genetic map and QTL mapping of ascochyta blight resistance in chickpea (Cicer arietinum L.).

KEY MESSAGE: A high-density linkage map of chickpea using 3430 SNPs was constructed and used to identify QTLs and candidate genes for ascochyta blight resistance in chickpea. Chickpea cultivation in temperate conditions is highly vulnerable to ascochyta blight infection. Cultivation of resistant cultivars in combination with fungicide application within an informed disease management package is the most effective method to control ascochyta blight in chickpeas. Identifying new sources of resistance is critical for continued improvement in ascochyta blight resistance in chickpea. The objective of this study was to identify genetic loci and candidate genes controlling the resistance to ascochyta blight in recombinant inbred lines derived from crossing cultivars Amit and ICCV 96029. The RILs were genotyped using the genotyping-by-sequencing procedure and Illumina® GoldenGate array. The RILs were evaluated in the field over three site-years and in three independent greenhouse experiments. A genetic map with eight linkage groups was constructed using 3430 SNPs. Eight QTLs for resistance were identified on chromosomes 2, 3, 4, 5 and 6. The QTLs individually explained 7-40% of the phenotypic variations. The QTLs on chromosomes 2 and 6 were associated with the resistance at vegetative stage only. The QTLs on chromosomes 2 and 4 that were previously reported to be conserved across diverse genetic backgrounds and against different isolates of Ascochyta rabiei were confirmed in this study. Candidate genes were identified within the QTL regions. Their co-localization with the underlying QTLs was confirmed by genetic mapping. The candidate gene-based SNP markers would lead to more efficient marker-assisted selection for ascochyta blight resistance and would provide a framework for fine mapping and subsequent cloning of the genes associated with the resistance.

Correction to: Genome-wide SNP discovery for development of high-density genetic map and QTL mapping of ascochyta blight resistance in chickpea (Cicer arietinum L.).

Unfortunately there was an error in the name of the gene "Ca-AKL18" in the discussion section.

Construction of high-density linkage maps for mapping quantitative trait loci for multiple traits in field pea (Pisum sativum L.).

BACKGROUND: The objective of this research was to map quantitative trait loci (QTLs) of multiple traits of breeding importance in pea (Pisum sativum L.). Three recombinant inbred line (RIL) populations, PR-02 (Orb x CDC Striker), PR-07 (Carerra x CDC Striker) and PR-15 (1-2347-144 x CDC Meadow) were phenotyped for agronomic and seed quality traits under field conditions over multiple environments in Saskatchewan, Canada. The mapping populations were genotyped using genotyping-by-sequencing (GBS) method for simultaneous single nucleotide polymorphism (SNP) discovery and construction of high-density linkage maps. RESULTS: After filtering for read depth, segregation distortion, and missing values, 2234, 3389 and 3541 single nucleotide polymorphism (SNP) markers identified by GBS in PR-02, PR-07 and PR-15, respectively, were used for construction of genetic linkage maps. Genetic linkage groups were assigned by anchoring to SNP markers previously positioned on these linkage maps. PR-02, PR-07 and PR-15 genetic maps represented 527, 675 and 609 non-redundant loci, and cover map distances of 951.9, 1008.8 and 914.2 cM, respectively. Based on phenotyping of the three mapping populations in multiple environments, 375 QTLs were identified for important traits including days to flowering, days to maturity, lodging resistance, Mycosphaerella blight resistance, seed weight, grain yield, acid and neutral detergent fiber concentration, seed starch concentration, seed shape, seed dimpling, and concentration of seed iron, selenium and zinc. Of all the QTLs identified, the most significant in terms of explained percentage of maximum phenotypic variance (PVmax) and occurrence in multiple environments were the QTLs for days to flowering (PVmax = 47.9%), plant height (PVmax = 65.1%), lodging resistance (PVmax = 35.3%), grain yield (PVmax = 54.2%), seed iron concentration (PVmax = 27.4%), and seed zinc concentration (PVmax = 43.2%). CONCLUSION: We have identified highly significant and reproducible QTLs for several agronomic and seed quality traits of breeding importance in pea. The QTLs identified will be the basis for fine mapping candidate genes, while some of the markers linked to the highly significant QTLs are useful for immediate breeding applications.

The Chickpea Early Flowering 1 (Efl1) Locus Is an Ortholog of Arabidopsis ELF3.

In climates that experience short growing seasons due to drought, heat, or end-of-season frost, early flowering is a highly desirable trait for chickpea (Cicer arietinum). In this study, we mapped, sequenced, and characterized Early flowering3 (Efl3), an ortholog of Arabidopsis (Arabidopsis thaliana) EARLY FLOWERING3 (ELF3) that confers early flowering in chickpea. In a recombinant inbred line population derived from a cross between CDC Frontier and ICCV 96029, this gene was mapped to the site of a quantitative trait locus on Ca5 that explained 59% of flowering time variation under short days. Sequencing of ELF3 in ICCV 96029 revealed an 11-bp deletion in the first exon that was predicted to result in a premature stop codon. The effect of this mutation was tested by transgenic complementation in the Arabidopsis elf3-1 mutant, with the CDC Frontier form of CaELF3a partially complementing elf3-1 while the ICCV 96029 form had no effect on flowering time. While induction of FLOWERING LOCUS T homologs was very early in ICCV 96029, an analysis of circadian clock function failed to show any clear loss of rhythm in the expression of clock genes in ICCV 96029 grown under continuous light, suggesting redundancy with other ELF3 homologs or possibly an alternative mode of action for this gene in chickpea. The 11-bp deletion was associated with early flowering in global chickpea germplasm but was not widely distributed, indicating that this mutation arose relatively recently.

Genotypic variation in the response of chickpea to arbuscular mycorrhizal fungi and non-mycorrhizal fungal endophytes.

Plant roots host symbiotic arbuscular mycorrhizal (AM) fungi and other fungal endophytes that can impact plant growth and health. The impact of microbial interactions in roots may depend on the genetic properties of the host plant and its interactions with root-associated fungi. We conducted a controlled condition experiment to investigate the effect of several chickpea (Cicer arietinum L.) genotypes on the efficiency of the symbiosis with AM fungi and non-AM fungal endophytes. Whereas the AM symbiosis increased the biomass of most of the chickpea cultivars, inoculation with non-AM fungal endophytes had a neutral effect. The chickpea cultivars responded differently to co-inoculation with AM fungi and non-AM fungal endophytes. Co-inoculation had additive effects on the biomass of some cultivars (CDC Corrine, CDC Anna, and CDC Cory), but non-AM fungal endophytes reduced the positive effect of AM fungi on Amit and CDC Vanguard. This study demonstrated that the response of plant genotypes to an AM symbiosis can be modified by the simultaneous colonization of the roots by non-AM fungal endophytes. Intraspecific variations in the response of chickpea to AM fungi and non-AM fungal endophytes indicate that the selection of suitable genotypes may improve the ability of crop plants to take advantage of soil ecosystem services.

Gene-based SNP discovery in tepary bean (Phaseolus acutifolius) and common bean (P. vulgaris) for diversity analysis and comparative mapping.

BACKGROUND: Common bean (Phaseolus vulgaris) is an important grain legume and there has been a recent resurgence in interest in its relative, tepary bean (P. acutifolius), owing to this species' ability to better withstand abiotic stresses. Genomic resources are scarce for this minor crop species and a better knowledge of the genome-level relationship between these two species would facilitate improvement in both. High-throughput genotyping has facilitated large-scale single nucleotide polymorphism (SNP) identification leading to the development of molecular markers with associated sequence information that can be used to place them in the context of a full genome assembly. RESULTS: Transcript-based SNPs were identified from six common bean and two tepary bean accessions and a subset were used to generate a 768-SNP Illumina GoldenGate assay for each species. The tepary bean assay was used to assess diversity in wild and cultivated tepary bean and to generate the first gene-based map of the tepary bean genome. Genotypic analyses of the diversity panel showed a clear separation between domesticated and cultivated tepary beans, two distinct groups within the domesticated types, and P. parvifolius was confirmed to be distinct. The genetic map of tepary bean was compared to the common bean genome assembly to demonstrate high levels of collinearity between the two species with differences limited to a few intra-chromosomal rearrangements. CONCLUSIONS: The development of the first set of genomic resources specifically for tepary bean has allowed for greater insight into the structure of this species and its relationship to its agriculturally more prominent relative, common bean. These resources will be helpful in the development of efficient breeding strategies for both species and will facilitate the introgression of agriculturally important traits from one crop into the other.

QTL mapping of early flowering and resistance to ascochyta blight in chickpea.

In western Canada, chickpea (Cicer arietinum L.) production is challenged by short growing seasons and infestations with ascochyta blight. Research was conducted to determine the genetic basis of the association between flowering time and reaction to ascochyta blight in chickpea. Ninety-two chickpea recombinant inbred lines (RILs) developed from a cross between ICCV 96029 and CDC Frontier were evaluated for flowering responses and ascochyta blight reactions in growth chambers and fields at multiple locations and during several years. A wide range of variation was exhibited by the RILs for days to flower, days to maturity, node of first flowering, plant height, and ascochyta blight resistance. Moderate to high broad sense heritability was estimated for ascochyta blight reaction (H(2) = 0.14-0.34) and for days to flowering (H(2) = 0.45-0.87) depending on the environments. Negative correlations were observed among the RILs for days to flowering and ascochyta blight resistance, ranging from r = -0.21 (P < 0.05) to -0.58 (P < 0.0001). A genetic linkage map consisting of eight linkage groups was developed using 349 SNP markers. Seven QTLs for days to flowering were identified that individually explained 9%-44% of the phenotypic variation. Eight QTLs were identified for ascochyta blight resistance that explained phenotypic variation ranging from 10% to 19%. Clusters of QTLs for days to flowering and ascochyta blight resistances were found on chromosome 3 at the interval of 8.6-23.11 cM and on chromosome 8 at the interval of 53.88-62.33 cM.

Genotype-specific variation in the structure of root fungal communities is related to chickpea plant productivity.

Increasing evidence supports the existence of variations in the association of plant roots with symbiotic fungi that can improve plant growth and inhibit pathogens. However, it is unclear whether intraspecific variations in the symbiosis exist among plant cultivars and if they can be used to improve crop productivity. In this study, we determined genotype-specific variations in the association of chickpea roots with soil fungal communities and evaluated the effect of root mycota on crop productivity. A 2-year field experiment was conducted in southwestern Saskatchewan, the central zone of the chickpea growing region of the Canadian prairie. The effects of 13 cultivars of chickpea, comprising a wide range of phenotypes and genotypes, were tested on the structure of root-associated fungal communities based on internal transcribed spacer (ITS) and 18S rRNA gene markers using 454 amplicon pyrosequencing. Chickpea cultivar significantly influenced the structure of the root fungal community. The magnitude of the effect varied with the genotypes evaluated, and effects were consistent across years. For example, the roots of CDC Corrine, CDC Cory, and CDC Anna hosted the highest fungal diversity and CDC Alma and CDC Xena the lowest. Fusarium sp. was dominant in chickpea roots but was less abundant in CDC Corrine than the other cultivars. A bioassay showed that certain of these fungal taxa, including Fusarium species, can reduce the productivity of chickpea, whereas Trichoderma harzianum can increase chickpea productivity. The large variation in the profile of chickpea root mycota, which included growth-promoting and -inhibiting species, supports the possibility of improving the productivity of chickpea by improving its root mycota in chickpea genetic improvement programs using traditional breeding techniques.

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.

Genetic characterization of the acetohydroxyacid synthase (AHAS) gene responsible for resistance to imidazolinone in chickpea (Cicer arietinum L.).

KEY MESSAGE: A point mutation in the AHAS1 gene leading to resistance to imidazolinone in chickpea was identified. The resistance is inherited as a single gene. A KASP marker targeting the mutation was developed. Weed control in chickpea (Cicer arietinum L.) is challenging due to poor crop competition ability and limited herbicide options. A chickpea genotype with resistance to imidazolinone (IMI) herbicides has been identified, but the genetic inheritance and the mechanism were unknown. In many plant species, resistance to IMI is caused by point mutation(s) in the acetohydroxyacid synthase (AHAS) gene resulting in an amino acid substitution preventing herbicide attachment to the molecule. The main objective of this research was to characterize the resistance to IMI herbicides in chickpea. Two homologous AHAS genes namely AHAS1 and AHAS2 sharing 80 % amino acid sequence similarity were identified in the chickpea genome. Cluster analysis indicated independent grouping of AHAS1 and AHAS2 across legume species. A point mutation in the AHAS1 gene at C675 to T675 resulting in an amino acid substitution from Ala205 to Val205 confers the resistance to IMI in chickpea. A KASP marker targeting the point mutation was developed and effectively predicted the response to IMI herbicides in a recombinant inbred (RI) population of chickpea. The RI population was used in molecular mapping where the major locus for the reaction to IMI herbicide was mapped to chromosome 5. Segregation analysis across an F2 population and RI population demonstrated that the resistance is inherited as a single gene in a semi-dominant fashion. The simple genetic inheritance and the availability of KASP marker generated in this study would speed up development of chickpea varieties with resistance to IMI herbicides.

Gene-based SNP discovery and genetic mapping in pea.

KEY MESSAGE: Gene-based SNPs were identified and mapped in pea using five recombinant inbred line populations segregating for traits of agronomic importance. Pea (Pisum sativum L.) is one of the world's oldest domesticated crops and has been a model system in plant biology and genetics since the work of Gregor Mendel. Pea is the second most widely grown pulse crop in the world following common bean. The importance of pea as a food crop is growing due to its combination of moderate protein concentration, slowly digestible starch, high dietary fiber concentration, and its richness in micronutrients; however, pea has lagged behind other major crops in harnessing recent advances in molecular biology, genomics and bioinformatics, partly due to its large genome size with a large proportion of repetitive sequence, and to the relatively limited investment in research in this crop globally. The objective of this research was the development of a genome-wide transcriptome-based pea single-nucleotide polymorphism (SNP) marker platform using next-generation sequencing technology. A total of 1,536 polymorphic SNP loci selected from over 20,000 non-redundant SNPs identified using deep transcriptome sequencing of eight diverse Pisum accessions were used for genotyping in five RIL populations using an Illumina GoldenGate assay. The first high-density pea SNP map defining all seven linkage groups was generated by integrating with previously published anchor markers. Syntenic relationships of this map with the model legume Medicago truncatula and lentil (Lens culinaris Medik.) maps were established. The genic SNP map establishes a foundation for future molecular breeding efforts by enabling both the identification and tracking of introgression of genomic regions harbouring QTLs related to agronomic and seed quality traits.

Analysis of acetohydroxyacid synthase1 gene in chickpea conferring resistance to imazamox herbicide.

Chickpea (Cicer arietinum L.) production in the Canadian prairies is challenging due to a lack of effective weed management mainly because of poor competition ability of the crop and limited registered herbicide options. Chickpea genotype with resistance to imidazolinone (IMI) herbicides has been identified. A point mutation in the acetohydroxyacid synthase1 (AHAS1) gene at C581 to T581, resulting in an amino acid substitution from Ala194 to Val194 (position 205, standardized to arabidopsis), confers the resistance to imazamox in chickpea. However, the molecular mechanism leading to the resistance is not fully understood. In many plant species, contrasting transcription levels of AHAS gene has been implicated in the resistant and susceptible genotypes in response to IMI. The objectives of this research were to compare the AHAS gene expression and AHAS enzyme activity in resistant and susceptible chickpea cultivars in response to imazamox herbicide treatment. Results from RT-qPCR indicated that there is no significant change in the transcript levels of AHAS1 between the susceptible and the resistant genotypes in response to imazamox treatment. Protein hydrophobic cluster analysis, protein-ligand docking analysis, and AHAS enzyme activity assay all indicated that the resistance to imazamox in chickpea is due to the alteration of interaction of the AHAS1 enzyme with the imazamox herbicide.

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.

Freezing stress response of wild and cultivated chickpeas.

Chickpea is an economically and nutritionally important grain legume globally, however, cold stress has adverse effects on its growth. In cold countries, like Canada where the growing season is short, having cold stress-tolerant varieties is crucial. Crop wild relatives of chickpea, especially Cicer reticulatum, can survive in suboptimal environments and are an important resource for crop improvement. In this study, we explored the performance of eleven C. reticulatum wild accessions and two chickpea cultivars, CDC Leader and CDC Consul, together with a cold sensitive check ILC533 under freezing stress. Freezing tolerance was scored based on a 1-9 scale. The wild relatives, particularly Kesen_075 and CudiA_152, had higher frost tolerance compared to the cultivars, which all died after frost treatment. We completed transcriptome analysis via mRNA sequencing to assess changes in gene expression in response to freezing stress and identified 6,184 differentially expressed genes (DEGs) in CDC Consul, and 7,842 DEGs in Kesen_075. GO (gene ontology) analysis of the DEGs revealed that those related to stress responses, endogenous and external stimuli responses, secondary metabolite processes, and photosynthesis were significantly over-represented in CDC Consul, while genes related to endogenous stimulus responses and photosynthesis were significantly over-represented in Kesen_075. These results are consistent with Kesen_075 being more tolerant to freezing stress than CDC Consul. Moreover, our data revealed that the expression of CBF pathway-related genes was impacted during freezing conditions in Kesen_075, and expression of these genes is believed to alleviate the damage caused by freezing stress. We identified genomic regions associated with tolerance to freezing stress in an F2 population derived from a cross between CDC Consul and Kesen_075 using QTL-seq analysis. Eight QTLs (P<0.05) on chromosomes Ca3, Ca4, Ca6, Ca7, Ca8, and two QTLs (P<0.01) on chromosomes Ca4 and Ca8, were associated with tolerance to freezing stress. Interestingly, 58 DEGs co-located within these QTLs. To our knowledge, this is the first study to explore the transcriptome and QTLs associated with freezing tolerance in wild relatives of chickpea under controlled conditions. Altogether, these findings provide comprehensive information that aids in understanding the molecular mechanism of chickpea adaptation to freezing stress and further provides functional candidate genes that can assist in breeding of freezing-stress tolerant varieties.