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.

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.

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.

Mapping Quantitative Trait Loci for Carotenoid Concentration in Three F2 Populations of Chickpea.

CORE IDEAS: Quantitative trait locus (QTL) analyses for carotenoids in chickpea were completed for three F2 populations. A moderate number of QTLs and candidate genes associated with carotenoid concentration in chickpea seeds were identified. Green cotyledon color is positively associated with provitamin A carotenoids. Three F2 populations derived from crosses between cultivars with green and yellow cotyledon colors were used to identify quantitative trait loci (QTLs) associated with carotenoid components in chickpea (Cicer arietinum L.) seeds developed by the Crop Development Centre (CDC). Carotenoids including violaxanthin, lutein, zeaxanthin, ß-cryptoxanthin, and ß-carotene were assessed in the F2:3 seeds via high-performance liquid chromatography (HPLC). In the 'CDC Jade' × 'CDC Frontier' population, 1068 bin markers derived from the 50K Axiom CicerSNP array were mapped onto eight linkage groups (LGs). Eight QTLs, including two each for ß-carotene and zeaxanthin and one each for total carotenoids, ß-cryptoxanthin, ß-carotene, and violaxanthin were identified in this population. In the 'CDC Cory' × 'CDC Jade' population, 694 bin markers were mapped onto eight LGs and one partial LG. Quantitative trait loci for ß-cryptoxanthin, ß-carotene, violaxanthin, lutein, and total carotenoids were identified on LG8. A map with eight LGs was developed from 581 bin markers in the third population derived from the 'ICC4475' × 'CDC Jade' cross. One QTL for ß-carotene and four QTLs, one each for ß-cryptoxanthin, ß-carotene, lutein, and total carotenoids, were identified in this population. The highest phenotypic variation explained by the QTLs was for ß-carotene, which ranged from 58 to 70% in all three populations. A major gene for cotyledon color was mapped on LG8 in each population. A significant positive correlation between cotyledon color and carotenoid concentration was observed. Potential candidate genes associated with carotenoid components were obtained and their locations on the chickpea genome are presented.

Genetic Analysis of NBS-LRR Gene Family in Chickpea and Their Expression Profiles in Response to Ascochyta Blight Infection.

Ascochyta blight is one of the major diseases of chickpea worldwide. The genetic resistance to ascochyta blight in chickpea is complex and governed by multiple QTLs. However, the molecular mechanism of quantitative disease resistance to ascochyta blight and the genes underlying these QTLs are still unknown. Most often disease resistance is determined by resistance (R) genes. The most predominant R-genes contain nucleotide binding site and leucine rich repeat (NBS-LRR) domains. A total of 121 NBS-LRR genes were identified in the chickpea genome. Ninety-eight of these genes contained all essential conserved domains while 23 genes were truncated. The NBS-LRR genes were grouped into eight distinct classes based on their domain architecture. Phylogenetic analysis grouped these genes into two major clusters based on their structural variation, the first cluster with toll or interleukin-1 like receptor (TIR) domain and the second cluster either with or without a coiled-coil domain. The NBS-LRR genes are distributed unevenly across the eight chickpea chromosomes and nearly 50% of the genes are present in clusters. Thirty of the NBS-LRR genes were co-localized with nine of the previously reported ascochyta blight QTLs and were tested as potential candidate genes for ascochyta blight resistance. Expression pattern of these genes was studied in two resistant (CDC Corinne and CDC Luna) and one susceptible (ICCV 96029) genotypes at different time points after ascochyta blight infection using real-time quantitative PCR. Twenty-seven NBS-LRR genes showed differential expression in response to ascochyta blight infection in at least one genotype at one time point. Among these 27 genes, the majority of the NBS-LRR genes showed differential expression after inoculation in both resistant and susceptible genotypes which indicates the involvement of these genes in response to ascochyta blight infection. Five NBS-LRR genes showed genotype specific expression. Our study provides a new insight of NBS-LRR gene family in chickpea and the potential involvement of NBS-LRR genes in response to ascochyta blight infection.

Genome-Wide Analysis of the Aquaporin Gene Family in Chickpea (Cicer arietinum L.).

Aquaporins (AQPs) are essential membrane proteins that play critical role in the transport of water and many other solutes across cell membranes. In this study, a comprehensive genome-wide analysis identified 40 AQP genes in chickpea (Cicer arietinum L.). A complete overview of the chickpea AQP (CaAQP) gene family is presented, including their chromosomal locations, gene structure, phylogeny, gene duplication, conserved functional motifs, gene expression, and conserved promoter motifs. To understand AQP's evolution, a comparative analysis of chickpea AQPs with AQP orthologs from soybean, Medicago, common bean, and Arabidopsis was performed. The chickpea AQP genes were found on all of the chickpea chromosomes, except chromosome 7, with a maximum of six genes on chromosome 6, and a minimum of one gene on chromosome 5. Gene duplication analysis indicated that the expansion of chickpea AQP gene family might have been due to segmental and tandem duplications. CaAQPs were grouped into four subfamilies including 15 NOD26-like intrinsic proteins (NIPs), 13 tonoplast intrinsic proteins (TIPs), eight plasma membrane intrinsic proteins (PIPs), and four small basic intrinsic proteins (SIPs) based on sequence similarities and phylogenetic position. Gene structure analysis revealed a highly conserved exon-intron pattern within CaAQP subfamilies supporting the CaAQP family classification. Functional prediction based on conserved Ar/R selectivity filters, Froger's residues, and specificity-determining positions suggested wide differences in substrate specificity among the subfamilies of CaAQPs. Expression analysis of the AQP genes indicated that some of the genes are tissue-specific, whereas few other AQP genes showed differential expression in response to biotic and abiotic stresses. Promoter profiling of CaAQP genes for conserved cis-acting regulatory elements revealed enrichment of cis-elements involved in circadian control, light response, defense and stress responsiveness reflecting their varying pattern of gene expression and potential involvement in biotic and abiotic stress responses. The current study presents the first detailed genome-wide analysis of the AQP gene family in chickpea and provides valuable information for further functional analysis to infer the role of AQP in the adaptation of chickpea in diverse environmental conditions.

The CarERF genes in chickpea (Cicer arietinum L.) and the identification of CarERF116 as abiotic stress responsive transcription factor.

The AP2/ERF family is one of the largest transcription factor gene families that are involved in various plant processes, especially in response to biotic and abiotic stresses. Complete genome sequences of one of the world's most important pulse crops chickpea (Cicer arietinum L.), has provided an important opportunity to identify and characterize genome-wide ERF genes. In this study, we identified 120 putative ERF genes from chickpea. The genomic organization of the chickpea ERF genes suggested that the gene family might have been expanded through the segmental duplications. The 120 member ERF family was classified into eleven distinct groups (I-X and VI-L). Transcriptional factor CarERF116, which is differentially expressed between drought tolerant and susceptible chickpea cultivar under terminal drought stress has been identified and functionally characterized. The CarERF116 encodes a putative protein of 241 amino acids and classified into group IX of ERF family. An in vitro CarERF116 protein-DNA binding assay demonstrated that CarERF116 protein specifically interacts with GCC box. We demonstrate that CarERF116 is capable of transactivation activity of and show that the functional transcriptional domain lies at the C-terminal region of the CarERF116. In transgenic Arabidopsis plants overexpressing CarERF116, significant up-regulation of several stress related genes were observed. These plants also exhibit resistance to osmotic stress and reduced sensitivity to ABA during seed germination. Based on these findings, we conclude that CarERF116 is an abiotic stress responsive gene, which plays an important role in stress tolerance. In addition, the present study leads to genome-wide identification and evolutionary analyses of chickpea ERF gene family, which will facilitate further research on this important group of genes and provides valuable resources for comparative genomics among the grain legumes.

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.

Construction and analysis of an SSH cDNA library of early heat-induced genes of Vigna aconitifolia variety RMO-40.

Moth bean ( Vigna aconitifolia (Jacq.) Marechal) is an important grain legume crop grown in rain fed areas of hot desert regions of Thar, India, under scorching sun rays with very little supplementation of water. An SSH cDNA library was generated from leaf tissues of V. aconitifolia var. RMO-40 exposed to an elevated temperature of 42 °C for 5 min to identify early-induced genes. A total of 488 unigenes (114 contigs and 374 singletons) were derived by cluster assembly and sequence alignment of 738 ESTs; out of 206 ESTs (28%) of unknown proteins, 160 ESTs (14%) were found to be novel to moth bean. Only 578 ESTs (78%) showed significant BLASTX similarity (<1 × 10(-6)) in the NCBI non-redundant database. Gene ontology functional classification terms were retrieved for 479 (65%) sequences, and 339 sequences were annotated with 165 EC codes and mapped to 68 different KEGG pathways. Four hundred and fifty-two ESTs were further annotated with InterProScan (IPS), and no IPS was assigned to 153 ESTs. In addition, the expression level of 27 ESTs in response to heat stress was evaluated through semiquantitative RT-PCR assay. Approximately 20 different signaling genes and 16 different transcription factors have been shown to be associated with heat stress in moth bean for the first time.

Comparative analysis of expressed sequence tags (ESTs) between drought-tolerant and -susceptible genotypes of chickpea under terminal drought stress.

BACKGROUND: Chickpea (Cicer arietinum L.) is an important grain-legume crop that is mainly grown in rainfed areas, where terminal drought is a major constraint to its productivity. We generated expressed sequence tags (ESTs) by suppression subtraction hybridization (SSH) to identify differentially expressed genes in drought-tolerant and -susceptible genotypes in chickpea. RESULTS: EST libraries were generated by SSH from root and shoot tissues of IC4958 (drought tolerant) and ICC 1882 (drought resistant) exposed to terminal drought conditions by the dry down method. SSH libraries were also constructed by using 2 sets of bulks prepared from the RNA of root tissues from selected recombinant inbred lines (RILs) (10 each) for the extreme high and low root biomass phenotype. A total of 3062 unigenes (638 contigs and 2424 singletons), 51.4% of which were novel in chickpea, were derived by cluster assembly and sequence alignment of 5949 ESTs. Only 2185 (71%) unigenes showed significant BLASTX similarity (<1E-06) in the NCBI non-redundant (nr) database. Gene ontology functional classification terms (BLASTX results and GO term), were retrieved for 2006 (92.0%) sequences, and 656 sequences were further annotated with 812 Enzyme Commission (EC) codes and were mapped to 108 different KEGG pathways. In addition, expression status of 830 unigenes in response to terminal drought stress was evaluated using macro-array (dot blots). The expression of few selected genes was validated by northern blotting and quantitative real-time PCR assay. CONCLUSION: Our study compares not only genes that are up- and down-regulated in a drought-tolerant genotype under terminal drought stress and a drought susceptible genotype but also between the bulks of the selected RILs exhibiting extreme phenotypes. More than 50% of the genes identified have been shown to be associated with drought stress in chickpea for the first time. This study not only serves as resource for marker discovery, but can provide a better insight into the selection of candidate genes (both up- and downregulated) associated with drought tolerance. These results can be used to identify suitable targets for manipulating the drought-tolerance trait in chickpea.

A comprehensive resource of drought- and salinity- responsive ESTs for gene discovery and marker development in chickpea (Cicer arietinum L.).

BACKGROUND: Chickpea (Cicer arietinum L.), an important grain legume crop of the world is seriously challenged by terminal drought and salinity stresses. However, very limited number of molecular markers and candidate genes are available for undertaking molecular breeding in chickpea to tackle these stresses. This study reports generation and analysis of comprehensive resource of drought- and salinity-responsive expressed sequence tags (ESTs) and gene-based markers. RESULTS: A total of 20,162 (18,435 high quality) drought- and salinity- responsive ESTs were generated from ten different root tissue cDNA libraries of chickpea. Sequence editing, clustering and assembly analysis resulted in 6,404 unigenes (1,590 contigs and 4,814 singletons). Functional annotation of unigenes based on BLASTX analysis showed that 46.3% (2,965) had significant similarity (< or =1E-05) to sequences in the non-redundant UniProt database. BLASTN analysis of unique sequences with ESTs of four legume species (Medicago, Lotus, soybean and groundnut) and three model plant species (rice, Arabidopsis and poplar) provided insights on conserved genes across legumes as well as novel transcripts for chickpea. Of 2,965 (46.3%) significant unigenes, only 2,071 (32.3%) unigenes could be functionally categorised according to Gene Ontology (GO) descriptions. A total of 2,029 sequences containing 3,728 simple sequence repeats (SSRs) were identified and 177 new EST-SSR markers were developed. Experimental validation of a set of 77 SSR markers on 24 genotypes revealed 230 alleles with an average of 4.6 alleles per marker and average polymorphism information content (PIC) value of 0.43. Besides SSR markers, 21,405 high confidence single nucleotide polymorphisms (SNPs) in 742 contigs (with > or = 5 ESTs) were also identified. Recognition sites for restriction enzymes were identified for 7,884 SNPs in 240 contigs. Hierarchical clustering of 105 selected contigs provided clues about stress- responsive candidate genes and their expression profile showed predominance in specific stress-challenged libraries. CONCLUSION: Generated set of chickpea ESTs serves as a resource of high quality transcripts for gene discovery and development of functional markers associated with abiotic stress tolerance that will be helpful to facilitate chickpea breeding. Mapping of gene-based markers in chickpea will also add more anchoring points to align genomes of chickpea and other legume species.