Cicer super-pangenome provides insights into species evolution and agronomic trait loci for crop improvement in chickpea.

Chickpea (Cicer arietinum L.)-an important legume crop cultivated in arid and semiarid regions-has limited genetic diversity. Efforts are being undertaken to broaden its diversity by utilizing its wild relatives, which remain largely unexplored. Here, we present the Cicer super-pangenome based on the de novo genome assemblies of eight annual Cicer wild species. We identified 24,827 gene families, including 14,748 core, 2,958 softcore, 6,212 dispensable and 909 species-specific gene families. The dispensable genome was enriched for genes related to key agronomic traits. Structural variations between cultivated and wild genomes were used to construct a graph-based genome, revealing variations in genes affecting traits such as flowering time, vernalization and disease resistance. These variations will facilitate the transfer of valuable traits from wild Cicer species into elite chickpea varieties through marker-assisted selection or gene-editing. This study offers valuable insights into the genetic diversity and potential avenues for crop improvement in chickpea.

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.

Genome-wide association mapping identifies an SNF4 ortholog that impacts biomass and sugar yield in sorghum and sugarcane.

Sorghum is a feed/industrial crop in developed countries and a staple food elsewhere in the world. This study evaluated the sorghum mini core collection for days to 50% flowering (DF), biomass, plant height (PH), soluble solid content (SSC), and juice weight (JW), and the sorghum reference set for DF and PH, in 7-12 testing environments. We also performed genome-wide association mapping with 6 094 317 and 265 500 single nucleotide polymorphism markers in the mini core collection and the reference set, respectively. In the mini core panel we identified three quantitative trait loci for DF, two for JW, one for PH, and one for biomass. In the reference set panel we identified another quantitative trait locus for PH on chromosome 6 that was also associated with biomass, DF, JW, and SSC in the mini core panel. Transgenic studies of three genes selected from the locus revealed that Sobic.006G061100 (SbSNF4-2) increased biomass, SSC, JW, and PH when overexpressed in both sorghum and sugarcane, and delayed flowering in transgenic sorghum. SbSNF4-2 encodes a γ subunit of the evolutionarily conserved AMPK/SNF1/SnRK1 heterotrimeric complexes. SbSNF4-2 and its orthologs will be valuable in genetic enhancement of biomass and sugar yield in plants.

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.

Genome-Wide Association Mapping Identifies Novel Panicle Morphology Loci and Candidate Genes in Sorghum.

Panicle morphology is an important trait in racial classification and can determine grain yield and other agronomic traits in sorghum. In this study, we performed association mapping of panicle length, panicle width, panicle compactness, and peduncle recurving in the sorghum mini core panel measured in multiple environments with 6,094,317 single nucleotide polymorphism (SNP) markers. We mapped one locus each on chromosomes 7 and 9 to recurving peduncles and eight loci for panicle length, panicle width, and panicle compactness. Because panicle length was positively correlated with panicle width, all loci for panicle length and width were colocalized. Among the eight loci, two each were on chromosomes 1, 2, and 6, and one each on chromosomes 8 and 10. The two loci on chromosome 2, i.e., Pm 2-1 and Pm 2-2, were detected in 7 and 5 out of 11 testing environments, respectively. Pm 2-2 colocalized with panicle compactness. Candidate genes were identified from both loci. The rice Erect Panicle2 (EP2) ortholog was among the candidate genes in Pm 2-2. EP2 regulates panicle erectness and panicle length in rice and encodes a novel plant-specific protein with unknown functions. The results of this study may facilitate the molecular identification of panicle morphology-related genes and the enhancement of yield and adaptation in sorghum.

Variation for Photoperiod and Temperature Sensitivity in the Global Mini Core Collection of Sorghum.

Information on photoperiod and temperature sensitivity of sorghum germplasm is important to identify appropriate sources for developing cultivars with a broad adaptation. The sorghum mini core collection consisting of 242 accessions along with three control cultivars were evaluated for days to 50% flowering (DFL) and plant height in two long-day rainy and two short-day post-rainy seasons, and for grain yield and 100-seed weight in the two post-rainy seasons. Differences in DFL and cumulative growing degree days (CGDD) in the rainy and post-rainy seasons were used to classify the accessions for photoperiod and temperature sensitivity. Results revealed 18 mini core landraces as photoperiod and temperature insensitive (PTINS), 205 as photoperiod sensitive and temperature insensitive (PSTINS), and 19 as photoperiod and temperature-sensitive (PTS) sources. The 19 PTS sources and 80 PSTINS sources took less DFL in the long-day rainy seasons than in the short-day post-rainy season indicating their adaptation to the rainy season and a possible different mechanism than that trigger flowering in the short-day sorghums. In all three groups, several accessions with desirable combinations of agronomic traits were identified for use in the breeding programs to develop climate-resilient cultivars and for genomic studies to identify genes responsible for the photoperiod and temperature responses.

Exploring the genetic variability and diversity of pearl millet core collection germplasm for grain nutritional traits improvement.

Improving essential nutrient content in staple food crops through biofortification breeding can overcome the micronutrient malnutrition problem. Genetic improvement depends on the availability of genetic variability in the primary gene pool. This study was aimed to ascertain the magnitude of variability in a core germplasm collection of diverse origin and predict pearl millet biofortification prospects for essential micronutrients. Germplasm accessions were evaluated in field trials at ICRISAT, India. The accessions differed significantly for all micronutrients with over two-fold variation for Fe (34-90 mg kg-1), Zn (30-74 mg kg-1), and Ca (85-249 mg kg-1). High estimates of heritability (> 0.81) were observed for Fe, Zn, Ca, P, Mo, and Mg. The lower magnitude of genotype (G) × environment (E) interaction observed for most of the traits implies strong genetic control for grain nutrients. The top-10 accessions for each nutrient and 15 accessions, from five countries for multiple nutrients were identified. For Fe and Zn, 39 accessions, including 15 with multiple nutrients, exceeded the Indian cultivars and 17 of them exceeded the biofortification breeding target for Fe (72 mg kg-1). These 39 accessions were grouped into 5 clusters. Most of these nutrients were positively and significantly associated among themselves and with days to 50% flowering and 1000-grain weight (TGW) indicating the possibility of their simultaneous improvement in superior agronomic background. The identified core collection accessions rich in specific and multiple-nutrients would be useful as the key genetic resources for developing biofortified and agronomically superior cultivars.

A sorghum practical haplotype graph facilitates genome-wide imputation and cost-effective genomic prediction.

Successful management and utilization of increasingly large genomic datasets is essential for breeding programs to accelerate cultivar development. To help with this, we developed a Sorghum bicolor Practical Haplotype Graph (PHG) pangenome database that stores haplotypes and variant information. We developed two PHGs in sorghum that were used to identify genome-wide variants for 24 founders of the Chibas sorghum breeding program from 0.01x sequence coverage. The PHG called single nucleotide polymorphisms (SNPs) with 5.9% error at 0.01x coverage-only 3% higher than PHG error when calling SNPs from 8x coverage sequence. Additionally, 207 progenies from the Chibas genomic selection (GS) training population were sequenced and processed through the PHG. Missing genotypes were imputed from PHG parental haplotypes and used for genomic prediction. Mean prediction accuracies with PHG SNP calls range from .57-.73 and are similar to prediction accuracies obtained with genotyping-by-sequencing or targeted amplicon sequencing (rhAmpSeq) markers. This study demonstrates the use of a sorghum PHG to impute SNPs from low-coverage sequence data and shows that the PHG can unify genotype calls across multiple sequencing platforms. By reducing input sequence requirements, the PHG can decrease the cost of genotyping, make GS more feasible, and facilitate larger breeding populations. Our results demonstrate that the PHG is a useful research and breeding tool that maintains variant information from a diverse group of taxa, stores sequence data in a condensed but readily accessible format, unifies genotypes across genotyping platforms, and provides a cost-effective option for genomic selection.

Phenotypic responses of foxtail millet (Setaria italica) genotypes to phosphate supply under greenhouse and natural field conditions.

Food insecurity is a looming threat for the burgeoning world population. Phosphorus (P), which is absorbed from soil as inorganic phosphate (Pi), is an essential macronutrient for the growth of all agricultural crops. This study reports phenotype analysis for P responses in natural field and greenhouse conditions, using 54 genotypes of foxtail millet (Setaria italica) representing wide geographic origins. The genotype responses were assessed in natural field conditions in two different seasons (monsoon and summer) under Pi-fertilized (P+) and unfertilized (P-) soil for eight above-ground traits. Enormous variations were seen among the genotypes in phenotypic responses for all the measured parameters under low P stress conditions. Variations were significant for plant height, leaf number and length, tillering ability and seed yield traits. Genotypes ISe 1234 and ISe 1541 were P+ responders, and the genotypes ISe 1181, ISe 1655, ISe 783 and ISe 1892 showed tolerance to low P for total seed yield. Genotypes that performed well under P- conditions were almost as productive as genotypes that performed well under P+ conditions suggesting some genotypes are well adapted to nutrient-poor soils. In the greenhouse, most of the genotypes produced changes in root architecture that are characteristic of P- stress, but to differing degrees. Significant variation was seen in root hair density and root hair number and in fresh and dry weight of shoot and root under P- stress. However, there was not much difference in the shoot and root total P and Pi levels of five selected high and low responding genotypes. We noticed contrasting responses in the greenhouse and natural field experiments for most of these genotypes. The leads from the study form the basis for breeding and improvement of foxtail millet for better Pi-use efficiency.

Superior haplotypes for haplotype-based breeding for drought tolerance in pigeonpea (Cajanus cajan L.).

Haplotype-based breeding, a recent promising breeding approach to develop tailor-made crop varieties, deals with identification of superior haplotypes and their deployment in breeding programmes. In this context, whole genome re-sequencing data of 292 genotypes from pigeonpea reference set were mined to identify the superior haplotypes for 10 drought-responsive candidate genes. A total of 83, 132 and 60 haplotypes were identified in breeding lines, landraces and wild species, respectively. Candidate gene-based association analysis of these 10 genes on a subset of 137 accessions of the pigeonpea reference set revealed 23 strong marker-trait associations (MTAs) in five genes influencing seven drought-responsive component traits. Haplo-pheno analysis for the strongly associated genes resulted in the identification of most promising haplotypes for three genes regulating five component drought traits. The haplotype C. cajan_23080-H2 for plant weight (PW), fresh weight (FW) and turgid weight (TW), the haplotype C. cajan_30211-H6 for PW, FW, TW and dry weight (DW), the haplotype C. cajan_26230-H11 for FW and DW and the haplotype C. cajan_26230-H5 for relative water content (RWC) were identified as superior haplotypes under drought stress condition. Furthermore, 17 accessions containing superior haplotypes for three drought-responsive genes were identified. The identified superior haplotypes and the accessions carrying these superior haplotypes will be very useful for deploying haplotype-based breeding to develop next-generation tailor-made better drought-responsive pigeonpea cultivars.

Genomic signatures of seed mass adaptation to global precipitation gradients in sorghum.

Seed mass is a key component of adaptation in plants and a determinant of yield in crops. The climatic drivers and genomic basis of seed mass variation remain poorly understood. In the cereal crop Sorghum bicolor, globally-distributed landraces harbor abundant variation in seed mass, which is associated with precipitation in their agroclimatic zones of origin. This study aimed to test the hypothesis that diversifying selection across precipitation gradients, acting on ancestral cereal grain size regulators, underlies seed mass variation in global sorghum germplasm. We tested this hypothesis in a set of 1901 georeferenced and genotyped sorghum landraces, 100-seed mass from common gardens, and bioclimatic precipitation variables. As predicted, 100-seed mass in global germplasm varies significantly among botanical races and is correlated to proxies of the precipitation gradients. With general and mixed linear model genome-wide associations, we identified 29 and 56 of 100 a priori candidate seed size genes with polymorphisms in the top 1% of seed mass association, respectively. Eleven of these genes harbor polymorphisms associated with the precipitation gradient, including orthologs of genes that regulate seed size in other cereals. With FarmCPU, 13 significant SNPs were identified, including one at an a priori candidate gene. Finally, we identified eleven colocalized outlier SNPs associated with seed mass and precipitation that also carry signatures of selection based on FST scans and PCAdapt, which represents a significant enrichment. Our findings suggest that seed mass in sorghum was shaped by diversifying selection on drought stress, and can inform genomics-enabled breeding for climate-resilient cereals.

Hypoallergen Peanut Lines Identified Through Large-Scale Phenotyping of Global Diversity Panel: Providing Hope Toward Addressing One of the Major Global Food Safety Concerns.

Peanut allergy is one of the serious health concern and affects more than 1% of the world's population mainly in Americas, Australia, and Europe. Peanut allergy is sometimes life-threatening and adversely affect the life quality of allergic individuals and their families. Consumption of hypoallergen peanuts is the best solution, however, not much effort has been made in this direction for identifying or developing hypoallergen peanut varieties. A highly diverse peanut germplasm panel was phenotyped using a recently developed monoclonal antibody-based ELISA protocol to quantify five major allergens. Results revealed a wide phenotypic variation for all the five allergens studied i.e., Ara h 1 (4-36,833 µg/g), Ara h 2 (41-77,041 µg/g), Ara h 3 (22-106,765 µg/g), Ara h 6 (829-103,892 µg/g), and Ara h 8 (0.01-70.12 µg/g). The hypoallergen peanut genotypes with low levels of allergen proteins for Ara h 1 (4 µg/g), Ara h 2 (41 µg/g), Ara h 3 (22 µg/g), Ara h 6 (829 µg/g), and Ara h 8 (0.01 µg/g) have paved the way for their use in breeding and genomics studies. In addition, these hypoallergen peanut genotypes are available for use in cultivation and industry, thus opened up new vistas for fighting against peanut allergy problem across the world.

Genome-wide cis-regulatory signatures for modulation of agronomic traits as exemplified by drought yield index (DYI) in chickpea.

Developing functional molecular tags from the cis-regulatory sequence components of genes is vital for their deployment in efficient genetic dissection of complex quantitative traits in crop plants including chickpea. The current study identified 431,194 conserved non-coding SNP (CNSNP) from the cis-regulatory element regions of genes which were annotated on a chickpea genome. These genome-wide CNSNP marker resources are made publicly accessible through a user-friendly web-database ( http://www.cnsnpcicarbase.com ). The CNSNP-based quantitative trait loci (QTL) and expression QTL (eQTL) mapping and genome-wide association study (GWAS) were further integrated with global gene expression landscapes, molecular haplotyping, and DNA-protein interaction study in the association panel and recombinant inbred lines (RIL) mapping population to decode complex genetic architecture of one of the vital seed yield trait under drought stress, drought yield index (DYI), in chickpea. This delineated two constituted natural haplotypes and alleles from a histone H3 protein-coding gene and its transcriptional regulator NAC transcription factor (TF) harboring the major QTLs and trans-acting eQTL governing DYI in chickpea. The effect of CNSNPs in TF-binding cis-element of a histone H3 gene in altering the binding affinity and transcriptional activity of NAC TF based on chromatin immunoprecipitation-quantitative PCR (ChIP-qPCR) assay was evident. The CNSNP-led promising molecular tags scanned will essentially have functional significance to decode transcriptional gene regulatory function and thus can drive translational genomic analysis in chickpea.

Resequencing of 429 chickpea accessions from 45 countries provides insights into genome diversity, domestication and agronomic traits.

We report a map of 4.97 million single-nucleotide polymorphisms of the chickpea from whole-genome resequencing of 429 lines sampled from 45 countries. We identified 122 candidate regions with 204 genes under selection during chickpea breeding. Our data suggest the Eastern Mediterranean as the primary center of origin and migration route of chickpea from the Mediterranean/Fertile Crescent to Central Asia, and probably in parallel from Central Asia to East Africa (Ethiopia) and South Asia (India). Genome-wide association studies identified 262 markers and several candidate genes for 13 traits. Our study establishes a foundation for large-scale characterization of germplasm and population genomics, and a resource for trait dissection, accelerating genetic gains in future chickpea breeding.

Genome wide association analysis of sorghum mini core lines regarding anthracnose, downy mildew, and head smut.

In previous studies, a sorghum mini core collection was scored over several years for response to Colletotrichum sublineola, Peronosclerospora sorghi, and Sporisorium reilianum, the causal agents of the disease anthracnose, downy mildew, and head smut, respectively. The screening results were combined with over 290,000 Single nucleotide polymorphic (SNP) loci from an updated version of a publicly available genotype by sequencing (GBS) dataset available for the mini core collection. GAPIT (Genome Association and Prediction Integrated Tool) R package was used to identify chromosomal locations that differ in disease response. When the top scoring SNPs were mapped to the most recent version of the published sorghum genome, in each case, a nearby and most often the closest annotated gene has precedence for a role in host defense.

CLAVATA signaling pathway genes modulating flowering time and flower number in chickpea.

KEY MESSAGE: A combinatorial genomic strategy delineated functionally relevant natural allele of a CLAVATA gene and its marker (haplotype)-assisted introgression led to development of the early-flowering chickpea cultivars with high flower number and enhanced yield/productivity. Unraveling the genetic components involved in CLAVATA (CLV) signaling is crucial for modulating important shoot apical meristem (SAM) characteristics and ultimately regulating diverse SAM-regulated agromorphological traits in crop plants. A genome-wide scan identified 142 CLV1-, 28 CLV2- and 6 CLV3-like genes, and their comprehensive genomic constitution and phylogenetic relationships were deciphered in chickpea. The QTL/fine mapping and map-based cloning integrated with high-resolution association analysis identified SNP loci from CaCLV3_01 gene within a major CaqDTF1.1/CaqFN1.1 QTL associated with DTF (days to 50% flowering) and FN (flower number) traits in chickpea, which was further ascertained by quantitative expression profiling. Molecular haplotyping of CaCLV3_01 gene, expressed specifically in SAM, constituted two major haplotypes that differentiated the early-DTF and high-FN chickpea accessions from late-DTF and low-FN. Enhanced accumulation of transcripts of superior CaCLV3_01 gene haplotype and known flowering promoting genes was observed in the corresponding haplotype-introgressed early-DTF and high-FN near-isogenic lines (NILs) with narrow SAM width. The superior haplotype-introgressed NILs exhibited early-flowering, high-FN and enhanced seed yield/productivity without compromising agronomic performance. These delineated molecular signatures can regulate DTF and FN traits through SAM proliferation and differentiation and thereby will be useful for translational genomic study to develop early-flowering cultivars with enhanced yield/productivity.

Transcriptional signatures modulating shoot apical meristem morphometric and plant architectural traits enhance yield and productivity in chickpea.

Plant height (PH) and plant width (PW), two of the major plant architectural traits determining the yield and productivity of a crop, are defined by diverse morphometric characteristics of the shoot apical meristem (SAM). The identification of potential molecular tags from a single gene that simultaneously modulates these plant/SAM architectural traits is therefore prerequisite to achieve enhanced yield and productivity in crop plants, including chickpea. Large-scale multienvironment phenotyping of the association panel and mapping population have ascertained the efficacy of three vital SAM morphometric trait parameters, SAM width, SAM height and SAM area, as key indicators to unravel the genetic basis of the wide PW and PH trait variations observed in desi chickpea. This study integrated a genome-wide association study (GWAS); quantitative trait locus (QTL)/fine-mapping and map-based cloning with molecular haplotyping; transcript profiling; and protein-DNA interaction assays for the dissection of plant architectural traits in chickpea. These exertions delineated natural alleles and superior haplotypes from a CabHLH121 transcription factor (TF) gene within the major QTL governing PW, PH and SAM morphometric traits. A genome-wide protein-DNA interaction assay assured the direct binding of a known stem cell master regulator, CaWUS, to the WOX-homeodomain TF binding sites of a CabHLH121 gene and its constituted haplotypes. The differential expression of CaWUS and transcriptional regulation of its target CabHLH121 gene/haplotypes were apparent, suggesting their collective role in altering SAM morphometric characteristics and plant architectural traits in the contrasting near isogenic lines (NILs). The NILs introgressed with a superior haplotype of a CabHLH121 exhibited optimal PW and desirable PH as well as enhanced yield and productivity without compromising any component of agronomic performance. These molecular signatures of the CabHLH121 TF gene have the potential to regulate both PW and PH traits through the modulation of proliferation, differentiation and maintenance of the meristematic stem cell population in the SAM; therefore, these signatures will be useful in the translational genomic study of chickpea genetic enhancement. The restructured cultivars with desirable PH (semidwarf) and PW will ensure maximal planting density in a specified cultivable field area, thereby enhancing the overall yield and productivity of chickpea. This can essentially facilitate the achievement of better remunerative outputs by farmers with rational land use, therefore ensuring global food security in the present scenario of an increasing population density and shrinking per capita land area.