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 (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. On the basis of WGS of 3,366 chickpea germplasm accessions, we report here a rich map of the genetic variation in chickpea. We provide a chickpea pan-genome and offer insights into species divergence, the migration of the cultigen (C. arietinum), rare allele burden and fitness loss in chickpea. We propose three genomic breeding approaches— haplotype-based breeding, genomic prediction and OCS—for developing tailor-made high-yielding and climate-resilient chickpea varieties. We sequenced 3,366 chickpea germplasm lines, including 3,171 cultivated and 195 wild accessions at an average coverage of around 12× (Methods, Extended Data Fig. 1, Supplementary Data 1 Tables 1, 2). Alignment of WGS data to the CDC Frontier reference genome11 identified 3.94 million and 19.57 million single-nucleotide polymorphisms (SNPs) in 3,171 cultivated and 195 wild accessions, respectively (Extended Data Table 1, Supplementary Data 1 Tables 3–7, Supplementary Notes). This SNP dataset was used to assess linkage disequilibrium (LD) decay (Supplementary Data 2 Tables 1, 2, Extended Data Fig. 2, Supplementary Notes) and identify private and population-enriched SNPs (Supplementary Data 3 Tables 1–4, Supplementary Notes). These private and population-enriched SNPs suggest rapid adaptation and can enhance the genetic foundation in the elite gene pool

Integrated transcriptome, small RNA and degradome sequencing approaches provide insights into Ascochyta blight resistance in chickpea

Ascochyta blight (AB) is one of the major biotic stresses known to limit the chickpea production worldwide. To dissect the complex mechanisms of AB resistance in chickpea, three approaches, namely, transcriptome, small RNA and degradome sequencing were used. The transcriptome sequencing of 20 samples including two resistant genotypes, two susceptible genotypes and one introgression line under control and stress conditions at two time points (3rd and 7th day post inoculation) identified a total of 6767 differentially expressed genes (DEGs). These DEGs were mainly related to pathogenesis�related proteins, disease resistance genes like NBS�LRR, cell wall biosynthesis and various secondary metabolite synthesis genes. The small RNA sequencing of the samples resulted in the identification of 651 miRNAs which included 478 known and 173 novel miRNAs. A total of 297 miRNAs were differentially expressed between different genotypes, conditions and time points. Using degradome sequencing and in silico approaches, 2131 targets were predicted for 629 miRNAs. The combined analysis of both small RNA and transcriptome datasets identified 12 miRNA�mRNA interaction pairs that exhibited contrasting expression in resistant and susceptible genotypes and also, a subset of genes that might be post�transcriptionally silenced during AB infection. The comprehensive integrated analysis in the study provides better insights into the transcriptome dynamics and regulatory network components associated with AB stress in chickpea and, also offers candidate genes for chickpea improvement

Genomic-enabled prediction models using multi-environment trials to estimate the effect of genotype × environment interaction on prediction accuracy in chickpea

Genomic selection (GS) by selecting lines prior to field phenotyping using genotyping data has the potential to enhance the rate of genetic gains. Genotype × environment (G × E) interaction inclusion in GS models can improve prediction accuracy hence aid in selection of lines across target environments. Phenotypic data on 320 chickpea breeding lines for eight traits for three seasons at two locations were recorded. These lines were genotyped using DArTseq (1.6 K SNPs) and Genotyping-by-Sequencing (GBS; 89 K SNPs). Thirteen models were fitted including main effects of environment and lines, markers, and/or naïve and informed interactions to estimate prediction accuracies. Three cross-validation schemes mimicking real scenarios that breeders might encounter in the fields were considered to assess prediction accuracy of the models (CV2: incomplete field trials or sparse testing; CV1: newly developed lines; and CV0: untested environments). Maximum prediction accuracies for different traits and different models were observed with CV2. DArTseq performed better than GBS and the combined genotyping set (DArTseq and GBS) regardless of the cross validation scheme with most of the main effect marker and interaction models. Improvement of GS models and application of various genotyping platforms are key factors for obtaining accurate and precise prediction accuracies, leading to more precise selection of candidates

Whole genome resequencing and phenotyping of MAGIC population for high resolution mapping of drought tolerance in chickpea

Terminal drought is one of the major constraints to crop production in chickpea (Cicer arietinum L.). In order to map drought tolerance related traits at high resolution, we sequenced multi-parent advanced generation intercross (MAGIC) population using whole genome resequencing approach and phenotyped it under drought stress environments for two consecutive years (2013-14 and 2014-15). A total of 52.02 billion clean reads containing 4.67 TB clean data were generated on the 1136 MAGIC lines and eight parental lines. Alignment of clean data on to the reference genome enabled identification of a total, 932,172 of SNPs, 35,973 insertions, and 35,726 deletions among the parental lines. A high-density genetic map was constructed using 57,180 SNPs spanning a map distance of 1606.69 cM. Using compressed mixed linear model, genome-wide association study (GWAS) enabled us to identify 737 markers significantly associated with days to 50% flowering, days to maturity, plant height, 100 seed weight, biomass, and harvest index. In addition to the GWAS approach, an identity-by-descent (IBD)-based mixed model approach was used to map quantitative trait loci (QTLs). The IBD-based mixed model approach detected major QTLs that were comparable to those from the GWAS analysis as well as some exclusive QTLs with smaller effects. The candidate genes like FRIGIDA and CaTIFY4b can be used for enhancing drought tolerance in chickpea. The genomic resources, genetic map, marker-trait associations, and QTLs identified in the study are valuable resources for the chickpea community for developing climate resilient chickpeas

Development and Application of High-Density Axiom Cajanus SNP Array with 56K SNPs to Understand the Genome Architecture of Released Cultivars and Founder Genotypes

As one of the major outputs of next-generation sequencing (NGS), a large number of genome-wide single-nucleotide polymorphisms (SNPs) have been developed in pigeonpea [Cajanus cajan (L.) Huth.]. However, SNPs require a genotyping platform or assay to be used in different evolutionary studies or in crop improvement programs. Therefore, we developed an Axiom Cajanus SNP array with 56K SNPs uniformly distributed across the genome and assessed its utility in a genetic diversity study. From the whole-genome resequencing (WGRS) data on 104 pigeonpea lines, ∼2 million sequence variations (SNPs and insertion–deletions [InDels]) were identified, from which a subset of 56,512 unique and informative sequence variations were selected to develop the array. The Axiom Cajanus SNP array developed was used for genotyping 103 pigeonpea lines encompassing 63 cultivars released between 1960 and 2014 and 40 breeding, germplasm, and founder lines. Genotyping data thus generated on 103 pigeonpea lines provided 51,201 polymorphic SNPs and InDels. Genetic diversity analysis provided in-depth insights into the genetic architecture and trends in temporal diversity in pigeonpea cultivars. Therefore, the continuous use of the high-density Axiom Cajanus SNP array developed will accelerate high-resolution trait mapping, marker-assisted breeding, and genomic selection efforts in pigeonpea

Genomic resources in plant breeding for sustainable agriculture

Climate change during the last 40 years has had a serious impact on agriculture and threatens global food and nutritional security. From over half a million plant species, cereals and legumes are the most important for food and nutritional security. Although systematic plant breeding has a relatively short history, conventional breeding coupled with advances in technology and crop management strategies has increased crop yields by 56 % globally between 1965-85, referred to as the Green Revolution. Nevertheless, increased demand for food, feed, fiber, and fuel necessitates the need to break existing yield barriers in many crop plants. In the first decade of the 21st century we witnessed rapid discovery, transformative technological development and declining costs of genomics technologies. In the second decade, the field turned towards making sense of the vast amount of genomic information and subsequently moved towards accurately predicting gene-to-phenotype associations and tailoring plants for climate resilience and global food security. In this review we focus on genomic resources, genome and germplasm sequencing, sequencing-based trait mapping, and genomics-assisted breeding approaches aimed at developing biotic stress resistant, abiotic stress tolerant and high nutrition varieties in six major cereals (rice, maize, wheat, barley, sorghum and pearl millet), and six major legumes (soybean, groundnut, cowpea, common bean, chickpea and pigeonpea). We further provide a perspective and way forward to use genomic breeding approaches including marker-assisted selection, marker-assisted backcrossing, haplotype based breeding and genomic prediction approaches coupled with machine learning and artificial intelligence, to speed breeding approaches. The overall goal is to accelerate genetic gains and deliver climate resilient and high nutrition crop varieties for sustainable agriculture

Pearl millet genome sequence provides a resource to improve agronomic traits in arid environments

Pearl millet [Pennisetum glaucum (L.) R. Br., syn. Cenchrus americanus (L.) Morrone], is a staple food for over 90 million poor farmers in arid and semi-arid regions of sub-Saharan Africa and South Asia. We report the ~1.79 Gb genome sequence of reference genotype Tift 23D2B1-P1-P5, which contains an estimated 38,579 genes. Resequencing analysis of 994 (963 inbreds of the highly cross-pollinated cultigen, and 31 wild accessions) provides insights into population structure, genetic diversity, evolution and domestication history. In addition we demonstrated the use of re-sequence data for establishing marker trait associations, genomic selection and prediction of hybrid performance and defining heterotic pools. The genome wide variations and abiotic stress proteome data are useful resources for pearl millet improvement through deploying modern breeding tools for accelerating genetic gains in pearl millet.publishersversionPeer reviewe

Transcriptome profiling reveals the expression and regulation of genes associated with Fusarium wilt resistance in chickpea (Cicer arietinum L.)

Abstract Fusarium wilt (FW) is one of the most significant biotic stresses limiting chickpea production worldwide. To dissect the molecular mechanism of FW resistance in chickpea, comparative transcriptome analyses of contrasting resistance sources of chickpea genotypes under control and Fusarium oxysporum f. sp. ciceris (Foc) inoculated conditions were performed. The high‐throughput transcriptome sequencing generated about 1137 million sequencing reads from 24 samples representing two resistant genotypes, two susceptible genotypes, and two near‐isogenic lines under control and stress conditions at two‐time points (7th‐ and 12th‐day post‐inoculation). The analysis identified 5182 differentially expressed genes (DEGs) between different combinations of chickpea genotypes. Functional annotation of these genes indicated their involvement in various biological processes such as defense response, cell wall biogenesis, secondary metabolism, and disease resistance. A significant number (382) of transcription factor encoding genes exhibited differential expression patterns under stress. Further, a considerable number of the identified DEGs (287) co‐localized with previously reported quantitative trait locus for FW resistance. Several resistance/susceptibility‐related genes, such as SERINE/THREONINE PROTEIN KINASE, DIRIGENT, and MLO exhibiting contrasting expression patterns in resistant and susceptible genotypes upon Foc inoculation, were identified. The results presented in the study provide valuable insights into the transcriptional dynamics associated with FW stress response in chickpea and provide candidate genes for the development of disease‐resistant chickpea cultivars