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.

Exploring Chickpea Germplasm Diversity for Broadening the Genetic Base Utilizing Genomic Resourses

Legume crops provide significant nutrition to humans as a source of protein, omega-3 fatty acids as well as specific macro and micronutrients. Additionally, legumes improve the cropping environment by replenishing the soil nitrogen content. Chickpeas are the second most significant staple legume food crop worldwide behind dry bean which contains 17%–24% protein, 41%–51% carbohydrate, and other important essential minerals, vitamins, dietary fiber, folate, β-carotene, anti-oxidants, micronutrients (phosphorus, calcium, magnesium, iron, and zinc) as well as linoleic and oleic unsaturated fatty acids. Despite these advantages, legumes are far behind cereals in terms of genetic improvement mainly due to far less effort, the bottlenecks of the narrow genetic base, and several biotic and abiotic factors in the scenario of changing climatic conditions. Measures are now called for beyond conventional breeding practices to strategically broadening of narrow genetic base utilizing chickpea wild relatives and improvement of cultivars through advanced breeding approaches with a focus on high yield productivity, biotic and abiotic stresses including climate resilience, and enhanced nutritional values. Desirable donors having such multiple traits have been identified using core and mini core collections from the cultivated gene pool and wild relatives of Chickpea. Several methods have been developed to address cross-species fertilization obstacles and to aid in inter-specific hybridization and introgression of the target gene sequences from wild Cicer species. Additionally, recent advances in “Omics” sciences along with high-throughput and precise phenotyping tools have made it easier to identify genes that regulate traits of interest. Next-generation sequencing technologies, whole-genome sequencing, transcriptomics, and differential genes expression profiling along with a plethora of novel techniques like single nucleotide polymorphism exploiting high-density genotyping by sequencing assays, simple sequence repeat markers, diversity array technology platform, and whole-genome re-sequencing technique led to the identification and development of QTLs and high-density trait mapping of the global chickpea germplasm. These altogether have helped in broadening the narrow genetic base of chickpeas.

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.

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.

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.

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.

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.

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 Population Structure Analyses of Three Minor Millets: Kodo Millet, Little Millet, and Proso Millet

Developed genome-wide SNP marker data for kodo, proso, and little millet Marker data used to analyze genetic diversity Heritability results of various traits used to validate genetic data Millets are a diverse group of small-seeded grains that are rich in nutrients but have received relatively little advanced plant breeding research. Millets are important to smallholder farmers in Africa and Asia because of their short growing season, good stress tolerance, and high nutritional content. To advance the study and use of these species, we present genome-wide marker datasets and population structure analyses for three minor millets: kodo millet (Paspalum scrobiculatum L.), little millet (Panicum sumatrense Roth), and proso millet (Panicum miliaceum L.).We generated genome-wide marker data sets for 190 accessions of each species with genotyping-by-sequencing (GBS). After filtering, we retained between 161 and 165 accessions of each species, with 3461, 2245, and 1882 single-nucleotide polymorphisms (SNPs) for kodo, little, and proso millet, respectively. Population genetic analysis revealed seven putative subpopulations of kodo millet and eight each of proso millet and little millet. To confirm the accuracy of this genetic data, we used public phenotype data on a subset of these accessions to estimate the heritability of various agronomically relevant phenotypes. Heritability values largely agree with the prior expectation for each phenotype, indicating that these SNPs provide an accurate genome-wide sample of genetic variation. These data represent one of first genome-wide population genetics analyses-and the most extensive-in these species and the first genomic analyses of any sort for little millet and kodo millet. These data will be a valuable resource for researchers and breeders trying to improve these crops for smallholder farmers.

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

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.

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.

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.

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.

ABC Transporter-Mediated Transport of Glutathione Conjugates Enhances Seed Yield and Quality in Chickpea

The identification of functionally relevant molecular tags is vital for genomics-assisted crop improvement and enhancement of seed yield, quality, and productivity in chickpea (Cicer arietinum). The simultaneous improvement of yield/productivity as well as quality traits often requires pyramiding of multiple genes, which remains a major hurdle given various associated epistatic and pleotropic effects. Unfortunately, no single gene that can improve yield/productivity along with quality and other desirable agromorphological traits is known, hampering the genetic enhancement of chickpea. Using a combinatorial genomics-assisted breeding and functional genomics strategy, this study identified natural alleles and haplotypes of an ABCC3-type transporter gene that regulates seed weight, an important domestication trait, by transcriptional regulation and modulation of the transport of glutathione conjugates in seeds of desi and kabuli chickpea. The superior allele/haplotype of this gene introgressed in desi and kabuli near-isogenic lines enhances the seed weight, yield, productivity, and multiple desirable plant architecture and seed-quality traits without compromising agronomic performance. These salient findings can expedite crop improvement endeavors and the development of nutritionally enriched high-yielding cultivars in chickpea.

Genetic dissection of photosynthetic efficiency traits for enhancing seed yield in chickpea

Understanding the genetic basis of photosynthetic efficiency (PE) contributing to enhanced seed yield per plant (SYP) is vital for genomics-assisted crop improvement of chickpea. The current study employed an integrated genomic strategy involving photosynthesis pathway gene-based association mapping, genome-wide association study, quantitative trait loci (QTL) mapping, and expression profiling. This identified 16 potential single nucleotide polymorphism loci linked to major QTLs underlying 16 candidate genes significantly associated with PE and SYP traits in chickpea. The allelic variants were tightly linked to positively interacting QTLs regulating both enhanced PE and SYP traits as exemplified by a chlorophyll A-B binding protein-coding gene. The leaf tissue-specific pronounced up-regulated expression of 16 associated genes in germplasm accessions and homozygous individuals of mapping population was evident. Such combinatorial genomic strategy coupled with gene haplotype-specific association and in silico protein-protein interaction study delineated natural alleles and superior haplotypes from a chlorophyll A-B binding (CAB) protein-coding gene and its interacting gene, Timing of CAB Expression 1 (TOC1), which appear to be most promising candidates in modulating chickpea PE and SYP traits. These functionally pertinent molecular signatures identified have efficacy to drive marker-assisted selection for developing PE-enriched cultivars with high seed yield in chickpea.

Biochemical diversity evaluation in chickpea accessions employing mini-core collection

The seeds of chickpea provide an exceptional source of dietary proteins and is one of the important legumes in both developed and developing countries over the world. The available germplasm of cultivated chickpea is deficient in desired biochemical signatures. To identify new sources of variations for breeding, reduced subsets of germplasm such as mini-core collection can be explored as an effective resource. In the present investigation, mini-core collections consisting of 215 accessions of chickpea were extensively evaluated for tapping biochemical diversity. Analysis included ten biochemical parameters comprising total protein, total free amino acids, phytic acid, tannin, total phenolics, total flavonoids, lectin, DPPH radical scavenging activity, in vitro digestibility of protein and starch. The spectrum of diversity was documented for total protein (4.60-33.90%), total free amino acids (0.092-9.33 mg/g), phytic acid (0.009-4.06 mg/g), tannin (0.232-189.63 mg/g), total phenolics (0.15-0.81 mg/g), total flavonoids (0.04-1.57 mg/g), lectin (0.07-330.32 HU/mg), DPPH radical scavenging activity (26.74-49.11%), in vitro protein digestibility (59.45-76.22%) and in vitro starch digestibility (45.63-298.39 mg of maltose/g). The principal component analysis revealed association of chickpea higher protein content to the lower level of total phenolics and flavonoid contents. The dendrogram obtained by unweighted pair group method using arithmetic average cluster analysis grouped the chickpea accessions into two major clusters. This is the first comprehensive report on biochemical diversity analysed in the mini-core chickpea accessions. The ultimate purpose of conducting such studies was to deliver information on nutritional characteristics for effective breeding programmes. Depending on the objectives of the breeding aforesaid accessions could be employed as a parent.

Genome-wide generation and genotyping of informative SNPs to scan molecular signatures for seed yield in chickpea

We discovered 2150 desi and 2199 kabuli accessions-derived SNPs by cultivar-wise individual assembling of sequence-reads generated through genotyping-by-sequencing of 92 chickpea accessions. Subsequent large-scale validation and genotyping of these SNPs discovered 619 desi accessions-derived (DAD) SNPs, 531 kabuli accessions-derived (KAD) SNPs, 884 multiple accessions-derived (MAD) SNPs and 1083 two accessions (desi ICC 4958 and kabuli CDC Frontier)-derived (TAD) SNPs that were mapped on eight chromosomes. These informative SNPs were annotated in coding/non-coding regulatory sequence components of genes. The MAD-SNPs were efficient to detect high intra-specific polymorphic potential and wide natural allelic diversity level including high-resolution admixed-population genetic structure and precise phylogenetic relationship among 291 desi and kabuli accessions. This signifies their effectiveness in introgression breeding and varietal improvement studies targeting useful agronomic traits of chickpea. Six trait-associated genes with SNPs including quantitative trait nucleotides (QTNs) in combination explained 27.5% phenotypic variation for seed yield per plant (SYP). A pentatricopeptide repeat (PPR) gene with a synonymous-coding SNP/QTN significantly associated with SYP trait was found most-promising in chickpea. The essential information delineated can be of immense utility in genomics-assisted breeding applications to develop high-yielding chickpea cultivars.

Cross-genera transferability of rice and finger millet genomic SSRs to barnyard millet (Echinochloa spp.)

Barnyard millet (Echinochloa spp.) is an important crop from nutritional point of view, nevertheless, the genetic information is very scarce. In the present investigation, rice and finger millet genomic SSRs were used for assessing cross transferability, identification of polymorphic markers, syntenic regions, genetic diversity and population structure analysis of barnyard millet genotypes. We observed 100% cross transferability for finger millet SSRs, of which 91% were polymorphic, while 71% of rice markers were cross transferable with 48% polymorphic out of them. Twenty-nine and sixteen highly polymorphic finger millet and rice SSRs yielded a mean of 4.3 and 3.38 alleles per locus in barnyard millet genotypes, respectively. The PIC values varied from 0.27 to 0.73 at an average of 0.54 for finger millet SSRs, whereas it was from 0.15 to 0.67 at an average of 0.44 for rice SSRs. High synteny was observed for markers related to panicle length, yield-related traits, spikelet fertility, plant height, root traits, leaf senescence, blast and brown plant hopper resistance. Although the rice SSRs located on chromosome 10 followed by chromosome 6 and 11 were found to be more transferable to barnyard millet, the finger millet SSRs were more polymorphic and transferable to barnyard millet genotypes. These SSR data of finger millet and rice individually as well as combined together grouped the 11 barnyard millet genotypes into 2 major clusters. The results of population structure analysis were similar to cluster analysis.

Diversifying Food Systems in the Pursuit of Sustainable Food Production and Healthy Diets

Increasing demand for nutritious, safe, and healthy food because of a growing population, and the pledge to maintain biodiversity and other resources, pose a major challenge to agriculture that is already threatened by a changing climate. Diverse and healthy diets, largely based on plant-derived food, may reduce diet-related illnesses. Investments in plant sciences will be necessary to design diverse cropping systems balancing productivity, sustainability, and nutritional quality. Cultivar diversity and nutritional quality are crucial. We call for better cooperation between food and medical scientists, food sector industries, breeders, and farmers to develop diversified and nutritious cultivars that reduce soil degradation and dependence on external inputs, such as fertilizers and pesticides, and to increase adaptation to climate change and resistance to emerging pests.

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

Draft genome, 994 re-sequenced lines and GWAS for yield-traits provide a resource of genetics and genomics tools for pearl millet researchers and breeders.

Pearl millet [Cenchrus americanus (L.) Morrone] is a staple food for more than 90 million farmers in arid and semi-arid regions of sub-Saharan Africa, India and South Asia. We report the ∼1.79 Gb draft whole genome sequence of reference genotype Tift 23D₂B₁-P1-P5, which contains an estimated 38,579 genes. We highlight the substantial enrichment for wax biosynthesis genes, which may contribute to heat and drought tolerance in this crop. We resequenced and analyzed 994 pearl millet lines, enabling insights into population structure, genetic diversity and domestication. We use these resequencing data to establish marker trait associations for genomic selection, to define heterotic pools, and to predict hybrid performance. We believe that these resources should empower researchers and breeders to improve this important staple crop.

Identification of putative QTLs for seedling stage phosphorus starvation response in finger millet (Eleusine coracana L. Gaertn.) by association mapping and cross species synteny analysis

A germplasm assembly of 128 finger millet genotypes from 18 countries was evaluated for seedling-stage phosphorus (P) responses by growing them in P sufficient (Psuf) and P deficient (Pdef) treatments. Majority of the genotypes showed adaptive responses to low P condition. Based on phenotype behaviour using the best linear unbiased predictors for each trait, genotypes were classified into, P responsive, low P tolerant and P non-responsive types. Based on the overall phenotype performance under Pdef, 10 genotypes were identified as low P tolerants. The low P tolerant genotypes were characterised by increased shoot and root length and increased root hair induction with longer root hairs under Pdef, than under Psuf. Association mapping of P response traits using mixed linear models revealed four quantitative trait loci (QTLs). Two QTLs (qLRDW.1 and qLRDW.2) for low P response affecting root dry weight explained over 10% phenotypic variation. In silico synteny analysis across grass genomes for these QTLs identified putative candidate genes such as Ser-Thr kinase and transcription factors such as WRKY and basic helix-loop-helix (bHLH). The QTLs for response under Psuf were mapped for traits such as shoot dry weight (qHSDW.1) and root length (qHRL.1). Putative associations of these QTLs over the syntenous regions on the grass genomes revealed proximity to cytochrome P450, phosphate transporter and pectin methylesterase inhibitor (PMEI) genes. This is the first report of the extent of phenotypic variability for P response in finger millet genotypes during seedling-stage, along with the QTLs and putative candidate genes associated with P starvation tolerance.

A Multiple QTL-Seq Strategy Delineates Potential Genomic Loci Governing Flowering Time in Chickpea

Identification of functionally relevant potential genomic loci using an economical, simpler and user-friendly genomics-assisted breeding strategy is vital for rapid genetic dissection of complex flowering time quantitative trait in chickpea. A high-throughput multiple QTL-seq strategy was employed in two inter (Cicer arietinum desi accession ICC 4958 × C reticulatum wild accession ICC 17160)- and intra (ICC 4958 × C. arietinum kabuli accession ICC 8261)-specific RIL mapping populations to identify the major QTL genomic regions governing flowering time in chickpea. The whole genome resequencing discovered 1635117 and 592486 SNPs exhibiting differentiation between early- and late-flowering mapping parents and bulks, constituted by pooling the homozygous individuals of extreme flowering time phenotypic trait from each of two aforesaid RIL populations. The multiple QTL-seq analysis using these mined SNPs in two RIL mapping populations narrowed-down two longer (907.1 kb and 1.99 Mb) major flowering time QTL genomic regions into the high-resolution shorter (757.7 kb and 1.39 Mb) QTL intervals on chickpea chromosome 4. This essentially identified regulatory as well as coding (non-synonymous/synonymous) novel SNP allelic variants from two efl1 (early flowering 1) and GI (GIGANTEA) genes regulating flowering time in chickpea. Interestingly, strong natural allelic diversity reduction (88-91%) of two known flowering genes especially mapped at major QTL intervals as compared to that of background genomic regions (where no flowering time QTLs were mapped; 61.8%) in cultivated vis-à-vis wild Cicer gene pools was evident inferring the significant impact of evolutionary bottlenecks on these loci during chickpea domestication. Higher association potential of coding non-synonymous and regulatory SNP alleles mined from efl1 (36-49%) and GI (33-42%) flowering genes for early and late flowering time differentiation among chickpea accessions was evident. The robustness and validity of two functional allelic variants-containing genes localized at major flowering time QTLs was apparent by their identification from multiple intra-/inter-specific mapping populations of chickpea. The functionally relevant molecular tags delineated can be of immense use for deciphering the natural allelic diversity-based domestication pattern of flowering time and expediting genomics-aided crop improvement to develop early flowering cultivars of chickpea.

Genotypic variation in soil water use and root distribution and their implications for drought tolerance in chickpea

Chickpeas are often grown under receding soil moisture and suffer ~50% yield losses due to drought stress. The timing of soil water use is considered critical for the efficient use of water under drought and to reduce yield losses. Therefore the root growth and the soil water uptake of 12 chickpea genotypes known for contrasts in drought and rooting response were monitored throughout the growth period both under drought and optimal irrigation. Root distribution reduced in the surface and increased in the deep soil layers below 30cm in response to drought. Soil water uptake was the maximum at 45-60cm soil depth under drought whereas it was the maximum at shallower (15-30 and 30-45cm) soil depths when irrigated. The total water uptake under drought was 1-fold less than optimal irrigation. The amount of water left unused remained the same across watering regimes. All the drought sensitive chickpea genotypes were inferior in root distribution and soil water uptake but the timing of water uptake varied among drought tolerant genotypes. Superiority in water uptake in most stages and the total water use determined the best adaptation. The water use at 15-30cm soil depth ensured greater uptake from lower depths and the soil water use from 90-120cm soil was critical for best drought adaptation. Root length density and the soil water uptake across soil depths were closely associated except at the surface or the ultimate soil depths of root presence.

Root traits confer grain yield advantages under terminal drought in chickpea (Cicer arietinum L.)

Chickpea, the second most important legume crop, suffers major yield losses by terminal drought stress (DS). Stronger root system is known to enhance drought yields but this understanding remains controversial. To understand precisely the root traits contribution towards yield, 12 chickpea genotypes with well-known drought response were field evaluated under drought and optimal irrigation. Root traits, such as root length density (RLD), total root dry weight (RDW), deep root dry weight (deep RDW) and root:shoot ratio (RSR), were measured periodically by soil coring up to 1.2 m soil depth across drought treatments. Large variations were observed for RLD, RDW, deep RDW and RSR in both the drought treatments. DS increased RLD below 30 cm soil depth, deep RDW, RSR but decreased the root diameter. DS increased the genetic variation in RDW more at the penultimate soil depths. Genetic variation under drought was the widest for RLD ∼50 DAS, for deep RDW ∼50-75 DAS and for RSR at 35 DAS. Genotypes ICC 4958, ICC 8261, Annigeri, ICC 14799, ICC 283 and ICC 867 at vegetative stage and genotypes ICC 14778, ICCV 10, ICC 3325, ICC 14799 and ICC 1882 at the reproductive phase produced greater RLD. Path- and correlation coefficients revealed strong positive contributions of RLD after 45 DAS, deep RDW at vicinity of maturity and RSR at early podfill stages to yield under drought. Breeding for the best combination of profuse RLD at surface soil depths, and RDW at deeper soil layers, was proposed to be the best selection strategy, for an efficient water use and an enhanced terminal drought tolerance in chickpea.

Development and Evaluation of a High Density Genotyping 'Axiom_Arachis' Array with 58 K SNPs for Accelerating Genetics and Breeding in Groundnut

Single nucleotide polymorphisms (SNPs) are the most abundant DNA sequence variation in the genomes which can be used to associate genotypic variation to the phenotype. Therefore, availability of a high-density SNP array with uniform genome coverage can advance genetic studies and breeding applications. Here we report the development of a high-density SNP array 'Axiom_Arachis' with 58 K SNPs and its utility in groundnut genetic diversity study. In this context, from a total of 163,782 SNPs derived from DNA resequencing and RNA-sequencing of 41 groundnut accessions and wild diploid ancestors, a total of 58,233 unique and informative SNPs were selected for developing the array. In addition to cultivated groundnuts (Arachis hypogaea), fair representation was kept for other diploids (A. duranensis, A. stenosperma, A. cardenasii, A. magna and A. batizocoi). Genotyping of the groundnut 'Reference Set' containing 300 genotypes identified 44,424 polymorphic SNPs and genetic diversity analysis provided in-depth insights into the genetic architecture of this material. The availability of the high-density SNP array 'Axiom_Arachis' with 58 K SNPs will accelerate the process of high resolution trait genetics and molecular breeding in cultivated groundnut.

Genome-wide development and deployment of informative intron-spanning and intron-length polymorphism markers for genomics-assisted breeding applications in chickpea

The discovery and large-scale genotyping of informative gene-based markers is essential for rapid delineation of genes/QTLs governing stress tolerance and yield component traits in order to drive genetic enhancement in chickpea. A genome-wide 119169 and 110491 ISM (intron-spanning markers) from 23129 desi and 20386 kabuli protein-coding genes and 7454 in silico InDel (insertion-deletion) (1-45-bp)-based ILP (intron-length polymorphism) markers from 3283 genes were developed that were structurally and functionally annotated on eight chromosomes and unanchored scaffolds of chickpea. A much higher amplification efficiency (83%) and intra-specific polymorphic potential (86%) detected by these markers than that of other sequence-based genetic markers among desi and kabuli chickpea accessions was apparent even by a cost-effective agarose gel-based assay. The genome-wide physically mapped 1718 ILP markers assayed a wider level of functional genetic diversity (19-81%) and well-defined phylogenetics among domesticated chickpea accessions. The gene-derived 1424 ILP markers were anchored on a high-density (inter-marker distance: 0.65cM) desi intra-specific genetic linkage map/functional transcript map (ICC 4958×ICC 2263) of chickpea. This reference genetic map identified six major genomic regions harbouring six robust QTLs mapped on five chromosomes, which explained 11-23% seed weight trait variation (7.6-10.5 LOD) in chickpea. The integration of high-resolution QTL mapping with differential expression profiling detected six including one potential serine carboxypeptidase gene with ILP markers (linked tightly to the major seed weight QTLs) exhibiting seed-specific expression as well as pronounced up-regulation especially in seeds of high (ICC 4958) as compared to low (ICC 2263) seed weight mapping parental accessions. The marker information generated in the present study was made publicly accessible through a user-friendly web-resource, "Chickpea ISM-ILP Marker Database". The designing of multiple ISM and ILP markers (2-5 markers/gene) from an individual gene (transcription factor) with numerous aforementioned desirable genetic attributes can widen the user-preference to select suitable primer combination for simultaneous large-scale assaying of functional allelic variation, natural allelic diversity, molecular mapping and expression profiling of genes among chickpea accessions. This will essentially accelerate the identification of functionally relevant molecular tags regulating vital agronomic traits for genomics-assisted crop improvement by optimal resource expenses in chickpea.

The Evolution of Photoperiod-Insensitive Flowering in Sorghum, A Genomic Model for Panicoid Grasses

Of central importance in adapting plants of tropical origin to temperate cultivation has been selection of daylength-neutral genotypes that flower early in the temperate summer and take full advantage of its long days. A cross between tropical and temperate sorghums [Sorghum propinquum (Kunth) Hitchc.×S. bicolor (L.) Moench], revealed a quantitative trait locus (QTL), FlrAvgD1, accounting for 85.7% of variation in flowering time under long days. Fine-scale genetic mapping placed FlrAvgD1 on chromosome 6 within the physically largest centiMorgan in the genome. Forward genetic data from “converted” sorghums validated the QTL. Association genetic evidence from a diversity panel delineated the QTL to a 10-kb interval containing only one annotated gene, Sb06g012260, that was shown by reverse genetics to complement a recessive allele. Sb06g012260 (SbFT12) contains a phosphatidylethanolamine-binding (PEBP) protein domain characteristic of members of the “FT” family of flowering genes acting as a floral suppressor. Sb06g012260 appears to have evolved ∼40 Ma in a panicoid ancestor after divergence from oryzoid and pooid lineages. A species-specific Sb06g012260 mutation may have contributed to spread to temperate regions by S. halepense (“Johnsongrass”), one of the world’s most widespread invasives. Alternative alleles for another family member, Sb02g029725 (SbFT6), mapping near another flowering QTL, also showed highly significant association with photoperiod response index (P = 1.53×10 ^− 6). The evolution of Sb06g012260 adds to evidence that single gene duplicates play large roles in important environmental adaptations. Increased knowledge of Sb06g012260 opens new doors to improvement of sorghum and other grain and cellulosic biomass crops.

Exploiting Phenylpropanoid Derivatives to Enhance the Nutraceutical Values of Cereals and Legumes

Phenylpropanoids are a diverse chemical class with immense health benefits that are biosynthesized from the aromatic amino acid L-phenylalanine. This article reviews the progress for accessing variation in phenylpropanoids in germplasm collections, the genetic and molecular basis of phenylpropanoid biosynthesis, and the development of cultivars dense in seed-phenylpropanoids. Progress is also reviewed on high-throughput assays, factors that influence phenylpropanoids, the site of phenylpropanoids accumulation in seed, Genotype × Environment interactions, and on consumer attitudes for the acceptance of staple foods rich in phenylpropanoids. A paradigm shift was noted in barley, maize, rice, sorghum, soybean, and wheat, wherein cultivars rich in phenylpropanoids are grown in Europe and North and Central America. Studies have highlighted some biological constraints that need to be addressed for development of high-yielding cultivars that are rich in phenylpropanoids. Genomics-assisted breeding is expected to facilitate rapid introgression into improved genetic backgrounds by minimizing linkage drag. More research is needed to systematically characterize germplasm pools for assessing variation to support crop genetic enhancement, and assess consumer attitudes to foods rich in phenylpropanoids.

EcoTILLING-Based Association Mapping Efficiently Delineates Functionally Relevant Natural Allelic Variants of Candidate Genes Governing Agronomic Traits in Chickpea

The large-scale mining and high-throughput genotyping of novel gene-based allelic variants in natural mapping population are essential for association mapping to identify functionally relevant molecular tags governing useful agronomic traits in chickpea. The present study employs an alternative time-saving, non-laborious and economical pool-based EcoTILLING approach coupled with agarose gel detection assay to discover 1133 novel SNP allelic variants from diverse coding and regulatory sequence components of 1133 transcription factor (TF) genes by genotyping in 192 diverse desi and kabuli chickpea accessions constituting a seed weight association panel. Integrating these SNP genotyping data with seed weight field phenotypic information of 192 structured association panel identified eight SNP alleles in the eight TF genes regulating seed weight of chickpea. The associated individual and combination of all SNPs explained 10-15 and 31% phenotypic variation for seed weight, respectively. The EcoTILLING-based large-scale allele mining and genotyping strategy implemented for association mapping is found much effective for a diploid genome crop species like chickpea with narrow genetic base and low genetic polymorphism. This optimized approach thus can be deployed for various genomics-assisted breeding applications with optimal expense of resources in domesticated chickpea. The seed weight-associated natural allelic variants and candidate TF genes delineated have potential to accelerate marker-assisted genetic improvement of chickpea.

Genetic dissection of seed-iron and zinc concentrations in chickpea

The SNP-based high-resolution QTL mapping mapped eight major genomic regions harbouring robust QTLs governing seed-Fe and Zn concentrations (39.4% combined phenotypic variation explained/PVE) on six chromosomes of an intra-specific high-density genetic linkage map (1.56 cM map-density). 24620 SNPs discovered from genome-wide GBS (genotyping-by-sequencing) and 13 known cloned Fe and Zn contents-related chickpea gene-orthologs were genotyped in a structured population of 92 sequenced desi and kabuli accessions. The large-scale 16591 SNP genotyping- and phenotyping-based GWAS (genome-wide association study) identified 16 genomic loci/genes associated (29% combined PVE) with seed-Fe and Zn concentrations. Of these, 11 trait-associated SNPs in the genes linked tightly with eight QTLs were validated by QTL mapping. The seed-specific expression, including pronounced differential-regulation of 16 trait-associated genes particularly in accessions/mapping individuals with contrasting level of seed-Fe and Zn contents was apparent. Collectively, the aforementioned rapid integrated genomic strategy led to delineate novel functional non-synonymous and regulatory SNP allelic-variants from 16 known/candidate genes, including three strong trait-associated genes (encoding late embryogenesis abundant and yellow stripe-like 1 protein, and vacuolar protein sorting-associated protein) and eight major QTLs regulating seed-Fe and Zn concentrations in chickpea. These essential inputs thus have potential to be deployed in marker-assisted genetic enhancement for developing nutritionally-rich iron/zinc-biofortified chickpea cultivars.

Identification of candidate genes for dissecting complex branch number trait in chickpea

The present study exploited integrated genomics-assisted breeding strategy for genetic dissection of complex branch number quantitative trait in chickpea. Candidate gene-based association analysis in a branch number association panel was performed by utilizing the genotyping data of 401 SNP allelic variants mined from 27 known cloned branch number gene orthologs of chickpea. The genome-wide association study (GWAS) integrating both genome-wide GBS- (4556 SNPs) and candidate gene-based genotyping information of 4957 SNPs in a structured population of 60 sequenced desi and kabuli accessions (with 350-400 kb LD decay), detected 11 significant genomic loci (genes) associated (41% combined PVE) with branch number in chickpea. Of these, seven branch number-associated genes were further validated successfully in two inter (ICC 4958 × ICC 17160)- and intra (ICC 12299 × ICC 8261)-specific mapping populations. The axillary meristem and shoot apical meristem-specific expression, including differential up- and down-regulation (4-5 fold) of the validated seven branch number-associated genes especially in high branch number as compared to the low branch number-containing parental accessions and homozygous individuals of two aforesaid mapping populations was apparent. Collectively, this combinatorial genomic approach delineated diverse naturally occurring novel functional SNP allelic variants in seven potential known/candidate genes [PIN1 (PIN-FORMED protein 1), TB1 (teosinte branched 1), BA1/LAX1 (BARREN STALK1/LIKE AUXIN1), GRAS8 (gibberellic acid insensitive/GAI, Repressor of ga13/RGA and Scarecrow8/SCR8), ERF (ethylene-responsive element-binding factor), MAX2 (more axillary growth 2) and lipase] governing chickpea branch number. The useful information generated from this study have potential to expedite marker-assisted genetic enhancement by developing high-yielding cultivars with more number of productive (pods and seeds) branches in chickpea.

Genome-Wide Scans for Delineation of Candidate Genes Regulating Seed-Protein Content in Chickpea

Identification of potential genes/alleles governing complex seed-protein content (SPC) is essential in marker-assisted breeding for quality trait improvement of chickpea. Henceforth, the present study utilized an integrated genomics-assisted breeding strategy encompassing trait association analysis, selective genotyping in traditional bi-parental mapping population and differential expression profiling for the first-time to understand the complex genetic architecture of quantitative SPC trait in chickpea. For GWAS (genome-wide association study), high-throughput genotyping information of 16376 genome-based SNPs (single nucleotide polymorphism) discovered from a structured population of 336 sequenced desi and kabuli accessions [with 150-200 kb LD (linkage disequilibrium) decay] was utilized. This led to identification of seven most effective genomic loci (genes) associated [10-20% with 41% combined PVE (phenotypic variation explained)] with SPC trait in chickpea. Regardless of the diverse desi and kabuli genetic backgrounds, a comparable level of association potential of the identified seven genomic loci with SPC trait was observed. Five SPC-associated genes were validated successfully in parental accessions and homozygous individuals of an intra-specific desi RIL (recombinant inbred line) mapping population (ICC 12299 × ICC 4958) by selective genotyping. The seed-specific expression, including differential up-regulation (>four fold) of six SPC-associated genes particularly in accessions, parents and homozygous individuals of the aforementioned mapping population with a high level of contrasting SPC (21-22%) was evident. Collectively, the integrated genomic approach delineated diverse naturally occurring novel functional SNP allelic variants in six potential candidate genes regulating SPC trait in chickpea. Of these, a non-synonymous SNP allele-carrying zinc finger transcription factor gene exhibiting strong association with SPC trait was found to be the most promising in chickpea. The informative functionally relevant molecular tags scaled-down essentially have potential to accelerate marker-assisted genetic improvement by developing nutritionally rich chickpea cultivars with enhanced SPC.

Landrace Germplasm for Improving Yield and Abiotic Stress Adaptation

Plant landraces represent heterogeneous, local adaptations of domesticated species, and thereby provide genetic resources that meet current and new challenges for farming in stressful environments. These local ecotypes can show variable phenology and low-to-moderate edible yield, but are often highly nutritious. The main contributions of landraces to plant breeding have been traits for more efficient nutrient uptake and utilization, as well as useful genes for adaptation to stressful environments such as water stress, salinity, and high temperatures. We propose that a systematic landrace evaluation may define patterns of diversity, which will facilitate identifying alleles for enhancing yield and abiotic stress adaptation, thus raising the productivity and stability of staple crops in vulnerable environments.

Selection of Host Differentials for Elucidating Pathogenic Variation in Magnaporthe grisea Populations Adapted to Finger Millet (Eleusine coracana)

Blast, caused by Pyricularia grisea (teleomorph: Magnaporthe grisea), is the most devastating disease of finger millet affecting production, utilization, and trade in Africa and Southeast Asia. An attempt was made to select a set of putative host differentials that can be used to determine virulence diversity in finger-millet-infecting populations of M. grisea. Thus, a differential set comprising eight germplasm accessions selected from finger millet core collection (IE 2911, IE 2957, IE 3392, IE 4497, IE 5091, IE 6240, IE 6337, and IE 7079) and a resistant ('GPU 28') and a susceptible ('VR 708') variety was developed. This differential set was used to study pathogenic variation in 25 isolates of M. grisea collected from Karnataka, Telangana, and Andhra Pradesh states in India. Based on the reaction (virulent = score ≥4 and avirulent = score ≤3 on a 1-to-9 scale) on host differentials, nine pathotypes were identified among 25 M. grisea isolates. Pathotype 9, represented by isolate Pg23 from Vizianagaram, was the most virulent because it could infect all of the host differentials except GPU 28. This study will be helpful in devising strategies for monitoring virulence change in M. grisea populations, and for identification of blast resistance in finger millet for use in disease resistance breeding programs.

Association mapping of germinability and seedling vigor in sorghum under controlled low-temperature conditions

Sorghum is one of the world's most important food, feed, and fiber crops as well as a potential feedstock for lignocellulosic bioenergy. Early-season planting extends sorghum's growing season and increases yield in temperate regions. However, sorghum's sensitivity to low soil temperatures adversely impacts seed germination. In this study, we evaluated the 242 accessions of the ICRISAT sorghum mini core collection for seed germination and seedling vigor at 12 °C as a measure of cold tolerance. Genome-wide association analysis was performed with approximately 162,177 single nucleotide polymorphism markers. Only one marker locus (Locus 7-2) was significantly associated with low-temperature germination and none with vigor. The linkage of Locus 7-2 to low-temperature germination was supported by four lines of evidence: strong association in three independent experiments, co-localization with previously mapped cold tolerance quantitative trait loci (QTL) in sorghum, a candidate gene that increases cold tolerance and germination rate when its wheat homolog is overexpressed in tobacco, and its syntenic region in rice co-localized with two cold tolerance QTL in rice. This locus may be useful in developing tools for molecular breeding of sorghums with improved low-temperature germinability.

A Genome-wide Combinatorial Strategy Dissects Complex Genetic Architecture of Seed Coat Color in Chickpea

The study identified 9045 high-quality SNPs employing both genome-wide GBS- and candidate gene-based SNP genotyping assays in 172, including 93 cultivated (desi and kabuli) and 79 wild chickpea accessions. The GWAS in a structured population of 93 sequenced accessions detected 15 major genomic loci exhibiting significant association with seed coat color. Five seed color-associated major genomic loci underlying robust QTLs mapped on a high-density intra-specific genetic linkage map were validated by QTL mapping. The integration of association and QTL mapping with gene haplotype-specific LD mapping and transcript profiling identified novel allelic variants (non-synonymous SNPs) and haplotypes in a MATE secondary transporter gene regulating light/yellow brown and beige seed coat color differentiation in chickpea. The down-regulation and decreased transcript expression of beige seed coat color-associated MATE gene haplotype was correlated with reduced proanthocyanidins accumulation in the mature seed coats of beige than light/yellow brown seed colored desi and kabuli accessions for their coloration/pigmentation. This seed color-regulating MATE gene revealed strong purifying selection pressure primarily in LB/YB seed colored desi and wild Cicer reticulatum accessions compared with the BE seed colored kabuli accessions. The functionally relevant molecular tags identified have potential to decipher the complex transcriptional regulatory gene function of seed coat coloration and for understanding the selective sweep-based seed color trait evolutionary pattern in cultivated and wild accessions during chickpea domestication. The genome-wide integrated approach employed will expedite marker-assisted genetic enhancement for developing cultivars with desirable seed coat color types in chickpea.

A genome-scale integrated approach aids in genetic dissection of complex flowering time trait in chickpea

A combinatorial approach of candidate gene-based association analysis and genome-wide association study (GWAS) integrated with QTL mapping, differential gene expression profiling and molecular haplotyping was deployed in the present study for quantitative dissection of complex flowering time trait in chickpea. Candidate gene-based association mapping in a flowering time association panel (92 diverse desi and kabuli accessions) was performed by employing the genotyping information of 5724 SNPs discovered from 82 known flowering chickpea gene orthologs of Arabidopsis and legumes as well as 832 gene-encoding transcripts that are differentially expressed during flower development in chickpea. GWAS using both genome-wide GBS- and candidate gene-based genotyping data of 30,129 SNPs in a structured population of 92 sequenced accessions (with 200-250 kb LD decay) detected eight maximum effect genomic SNP loci (genes) associated (34% combined PVE) with flowering time. Six flowering time-associated major genomic loci harbouring five robust QTLs mapped on a high-resolution intra-specific genetic linkage map were validated (11.6-27.3% PVE at 5.4-11.7 LOD) further by traditional QTL mapping. The flower-specific expression, including differential up- and down-regulation (>three folds) of eight flowering time-associated genes (including six genes validated by QTL mapping) especially in early flowering than late flowering contrasting chickpea accessions/mapping individuals during flower development was evident. The gene haplotype-based LD mapping discovered diverse novel natural allelic variants and haplotypes in eight genes with high trait association potential (41% combined PVE) for flowering time differentiation in cultivated and wild chickpea. Taken together, eight potential known/candidate flowering time-regulating genes [efl1 (early flowering 1), FLD (Flowering locus D), GI (GIGANTEA), Myb (Myeloblastosis), SFH3 (SEC14-like 3), bZIP (basic-leucine zipper), bHLH (basic helix-loop-helix) and SBP (SQUAMOSA promoter binding protein)], including novel markers, QTLs, alleles and haplotypes delineated by aforesaid genome-wide integrated approach have potential for marker-assisted genetic improvement and unravelling the domestication pattern of flowering time in chickpea.

Population Genetics and Structure of a Global Foxtail Millet Germplasm Collection

Foxtail millet [Setaria italica (L.) P. Beauv.] is one among the most ancient crops of dryland agriculture. It is the second most important crop among millets grown for grains or forage. Foxtail millet germplasm resources provide reservoirs of novel alleles and genes for crop improvement that have remained mostly unexplored. We genotyped a set of 190 foxtail millet germplasm accessions (including 155 accessions of the foxtail millet core collection) using genotyping-by-sequencing (GBS) for rapid single nucleotide polymorphisms (SNP) characterization to study population genetics and structure, which enable allele mining through association mapping approaches. After filtering a total 350,000 raw SNPs identified across 190 germplasm accessions for minor allele frequency (MAF), coverage for samples and coverage for sites, we retained 181 accessions with 17,714 high-quality SNPs with ≥5% MAF. Genetic structure analyses revealed that foxtail millet germplasm accessions are structured along both on the basis of races and geographic origin, and the maximum proportion of variation was due to among individuals within populations. Accessions of race indica were less diverse and are highly differentiated from those of maxima and moharia. Genome-wide linkage disequilibrium (LD) analysis showed on an average LD extends up to ∼150 kbp and varied with individual chromosomes. The utility of the data for performing genome-wide association studies (GWASs) was tested with plant pigmentation and days to flowering and identified significant marker-trait associations. This SNP data provides a foundation for exploration of foxtail millet diversity and for mining novel alleles and mapping genes for economically important traits.

Genome-wide insertion-deletion (InDel) marker discovery and genotyping for genomics-assisted breeding applications in chickpea

We developed 21,499 genome-wide insertion-deletion (InDel) markers (2- to 54-bp in silico fragment length polymorphism) by comparing the genomic sequences of four (desi, kabuli and wild C. reticulatum) chickpea [Cicer arietinum (L.)] accessions. InDel markers showing 2- to 6-bp fragment length polymorphism among accessions were abundant (76.8%) in the chickpea genome. The physically mapped 7,643 and 13,856 markers on eight chromosomes and unanchored scaffolds, respectively, were structurally and functionally annotated. The 4,506 coding (23% large-effect frameshift mutations) and regulatory InDel markers were identified from 3,228 genes (representing 11.7% of total 27,571 desi genes), suggesting their functional relevance for trait association/genetic mapping. High amplification (97%) and intra-specific polymorphic (60-83%) potential and wider genetic diversity (15-89%) were detected by genome-wide 6,254 InDel markers among desi, kabuli and wild accessions using even a simpler cost-effective agarose gel-based assay. This signifies added advantages of this user-friendly genetic marker system for manifold large-scale genotyping applications in laboratories with limited infrastructure and resources. Utilizing 6,254 InDel markers-based high-density (inter-marker distance: 0.212 cM) inter-specific genetic linkage map (ICC 4958 × ICC 17160) of chickpea as a reference, three major genomic regions harboring six flowering and maturity time robust QTLs (16.4-27.5% phenotypic variation explained, 8.1-11.5 logarithm of odds) were identified. Integration of genetic and physical maps at these target QTL intervals mapped on three chromosomes delineated five InDel markers-containing candidate genes tightly linked to the QTLs governing flowering and maturity time in chickpea. Taken together, our study demonstrated the practical utility of developing and high-throughput genotyping of such beneficial InDel markers at a genome-wide scale to expedite genomics-assisted breeding applications in chickpea.

Haploids: Constraints and opportunities in plant breeding

The discovery of haploids in higher plants led to the use of doubled haploid (DH) technology in plant breeding. This article provides the state of the art on DH technology including the induction and identification of haploids, what factors influence haploid induction, molecular basis of microspore embryogenesis, the genetics underpinnings of haploid induction and its use in plant breeding, particularly to fix traits and unlock genetic variation. Both in vitro and in vivo methods have been used to induce haploids that are thereafter chromosome doubled to produce DH. Various heritable factors contribute to the successful induction of haploids, whose genetics is that of a quantitative trait. Genomic regions associated with in vitro and in vivo DH production were noted in various crops with the aid of DNA markers. It seems that F2 plants are the most suitable for the induction of DH lines than F1 plants. Identifying putative haploids is a key issue in haploid breeding. DH technology in Brassicas and cereals, such as barley, maize, rice, rye and wheat, has been improved and used routinely in cultivar development, while in other food staples such as pulses and root crops the technology has not reached to the stage leading to its application in plant breeding. The centromere-mediated haploid induction system has been used in Arabidopsis, but not yet in crops. Most food staples are derived from genomic resources-rich crops, including those with sequenced reference genomes. The integration of genomic resources with DH technology provides new opportunities for the improving selection methods, maximizing selection gains and accelerate cultivar development. Marker-aided breeding and DH technology have been used to improve host plant resistance in barley, rice, and wheat. Multinational seed companies are using DH technology in large-scale production of inbred lines for further development of hybrid cultivars, particularly in maize. The public sector provides support to national programs or small-medium private seed for the exploitation of DH technology in plant breeding.

Genome-environment associations in sorghum landraces predict adaptive traits

Improving environmental adaptation in crops is essential for food security under global change, but phenotyping adaptive traits remains a major bottleneck. If associations between single-nucleotide polymorphism (SNP) alleles and environment of origin in crop landraces reflect adaptation, then these could be used to predict phenotypic variation for adaptive traits. We tested this proposition in the global food crop Sorghum bicolor, characterizing 1943 georeferenced landraces at 404,627 SNPs and quantifying allelic associations with bioclimatic and soil gradients. Environment explained a substantial portion of SNP variation, independent of geographical distance, and genic SNPs were enriched for environmental associations. Further, environment-associated SNPs predicted genotype-by-environment interactions under experimental drought stress and aluminum toxicity. Our results suggest that genomic signatures of environmental adaptation may be useful for crop improvement, enhancing germplasm identification and marker-assisted selection. Together, genome-environment associations and phenotypic analyses may reveal the basis of environmental adaptation.

A genome-wide SNP scan accelerates trait-regulatory genomic loci identification in chickpea

We identified 44844 high-quality SNPs by sequencing 92 diverse chickpea accessions belonging to a seed and pod trait-specific association panel using reference genome- and de novo-based GBS (genotyping-by-sequencing) assays. A GWAS (genome-wide association study) in an association panel of 211, including the 92 sequenced accessions, identified 22 major genomic loci showing significant association (explaining 23-47% phenotypic variation) with pod and seed number/plant and 100-seed weight. Eighteen trait-regulatory major genomic loci underlying 13 robust QTLs were validated and mapped on an intra-specific genetic linkage map by QTL mapping. A combinatorial approach of GWAS, QTL mapping and gene haplotype-specific LD mapping and transcript profiling uncovered one superior haplotype and favourable natural allelic variants in the upstream regulatory region of a CesA-type cellulose synthase (Ca_Kabuli_CesA3) gene regulating high pod and seed number/plant (explaining 47% phenotypic variation) in chickpea. The up-regulation of this superior gene haplotype correlated with increased transcript expression of Ca_Kabuli_CesA3 gene in the pollen and pod of high pod/seed number accession, resulting in higher cellulose accumulation for normal pollen and pollen tube growth. A rapid combinatorial genome-wide SNP genotyping-based approach has potential to dissect complex quantitative agronomic traits and delineate trait-regulatory genomic loci (candidate genes) for genetic enhancement in crop plants, including chickpea.

Ultra-high density intra-specific genetic linkage maps accelerate identification of functionally relevant molecular tags governing important agronomic traits in chickpea

We discovered 26785 and 16573 high-quality SNPs differentiating two parental genotypes of a RIL mapping population using reference desi and kabuli genome-based GBS assay. Of these, 3625 and 2177 SNPs have been integrated into eight desi and kabuli chromosomes, respectively in order to construct ultra-high density (0.20-0.37 cM) intra-specific chickpea genetic linkage maps. One of these constructed high-resolution genetic map has potential to identify 33 major genomic regions harbouring 35 robust QTLs (PVE: 17.9-39.7%) associated with three agronomic traits, which were mapped within <1 cM mean marker intervals on desi chromosomes. The extended LD (linkage disequilibrium) decay (~15 cM) in chromosomes of genetic maps have encouraged us to use a rapid integrated approach (comparative QTL mapping, QTL-region specific haplotype/LD-based trait association analysis, expression profiling and gene haplotype-based association mapping) rather than a traditional QTL map-based cloning method to narrow-down one major seed weight (SW) robust QTL region. It delineated favourable natural allelic variants and superior haplotype-containing one seed-specific candidate embryo defective gene regulating SW in chickpea. The ultra-high-resolution genetic maps, QTLs/genes and alleles/haplotypes-related genomic information generated and integrated strategy for rapid QTL/gene identification developed have potential to expedite genomics-assisted breeding applications in crop plants, including chickpea for their genetic enhancement.