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.

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.

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.

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.

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.

Employing genome-wide SNP discovery and genotyping strategy to extrapolate the natural allelic diversity and domestication patterns in chickpea

The genome-wide discovery and high-throughput genotyping of SNPs in chickpea natural germplasm lines is indispensable to extrapolate their natural allelic diversity, domestication, and linkage disequilibrium (LD) patterns leading to the genetic enhancement of this vital legume crop. We discovered 44,844 high-quality SNPs by sequencing of 93 diverse cultivated desi, kabuli, and wild chickpea accessions using reference genome- and de novo-based GBS (genotyping-by-sequencing) assays that were physically mapped across eight chromosomes of desi and kabuli. Of these, 22,542 SNPs were structurally annotated in different coding and non-coding sequence components of genes. Genes with 3296 non-synonymous and 269 regulatory SNPs could functionally differentiate accessions based on their contrasting agronomic traits. A high experimental validation success rate (92%) and reproducibility (100%) along with strong sensitivity (93-96%) and specificity (99%) of GBS-based SNPs was observed. This infers the robustness of GBS as a high-throughput assay for rapid large-scale mining and genotyping of genome-wide SNPs in chickpea with sub-optimal use of resources. With 23,798 genome-wide SNPs, a relatively high intra-specific polymorphic potential (49.5%) and broader molecular diversity (13-89%)/functional allelic diversity (18-77%) was apparent among 93 chickpea accessions, suggesting their tremendous applicability in rapid selection of desirable diverse accessions/inter-specific hybrids in chickpea crossbred varietal improvement program. The genome-wide SNPs revealed complex admixed domestication pattern, extensive LD estimates (0.54-0.68) and extended LD decay (400-500 kb) in a structured population inclusive of 93 accessions. These findings reflect the utility of our identified SNPs for subsequent genome-wide association study (GWAS) and selective sweep-based domestication trait dissection analysis to identify potential genomic loci (gene-associated targets) specifically regulating important complex quantitative agronomic traits in chickpea. The numerous informative genome-wide SNPs, natural allelic diversity-led domestication pattern, and LD-based information generated in our study have got multidimensional applicability with respect to chickpea genomics-assisted breeding.

A combinatorial approach of comprehensive QTL-based comparative genome mapping and transcript profiling identified a seed weight-regulating candidate gene in chickpea

High experimental validation/genotyping success rate (94-96%) and intra-specific polymorphic potential (82-96%) of 1536 SNP and 472 SSR markers showing in silico polymorphism between desi ICC 4958 and kabuli ICC 12968 chickpea was obtained in a 190 mapping population (ICC 4958 × ICC 12968) and 92 diverse desi and kabuli genotypes. A high-density 2001 marker-based intra-specific genetic linkage map comprising of eight LGs constructed is comparatively much saturated (mean map-density: 0.94 cM) in contrast to existing intra-specific genetic maps in chickpea. Fifteen robust QTLs (PVE: 8.8-25.8% with LOD: 7.0-13.8) associated with pod and seed number/plant (PN and SN) and 100 seed weight (SW) were identified and mapped on 10 major genomic regions of eight LGs. One of 126.8 kb major genomic region harbouring a strong SW-associated robust QTL (Caq'SW1.1: 169.1-171.3 cM) has been delineated by integrating high-resolution QTL mapping with comprehensive marker-based comparative genome mapping and differential expression profiling. This identified one potential regulatory SNP (G/A) in the cis-acting element of candidate ERF (ethylene responsive factor) TF (transcription factor) gene governing seed weight in chickpea. The functionally relevant molecular tags identified have potential to be utilized for marker-assisted genetic improvement of chickpea.

Genome-wide conserved non-coding microsatellite (CNMS) marker-based integrative genetical genomics for quantitative dissection of seed weight in chickpea

Phylogenetic footprinting identified 666 genome-wide paralogous and orthologous CNMS (conserved non-coding microsatellite) markers from 5'-untranslated and regulatory regions (URRs) of 603 protein-coding chickpea genes. The (CT)n and (GA)n CNMS carrying CTRMCAMV35S and GAGA8BKN3 regulatory elements, respectively, are abundant in the chickpea genome. The mapped genic CNMS markers with robust amplification efficiencies (94.7%) detected higher intraspecific polymorphic potential (37.6%) among genotypes, implying their immense utility in chickpea breeding and genetic analyses. Seventeen differentially expressed CNMS marker-associated genes showing strong preferential and seed tissue/developmental stage-specific expression in contrasting genotypes were selected to narrow down the gene targets underlying seed weight quantitative trait loci (QTLs)/eQTLs (expression QTLs) through integrative genetical genomics. The integration of transcript profiling with seed weight QTL/eQTL mapping, molecular haplotyping, and association analyses identified potential molecular tags (GAGA8BKN3 and RAV1AAT regulatory elements and alleles/haplotypes) in the LOB-domain-containing protein- and KANADI protein-encoding transcription factor genes controlling the cis-regulated expression for seed weight in the chickpea. This emphasizes the potential of CNMS marker-based integrative genetical genomics for the quantitative genetic dissection of complex seed weight in chickpea.

Pre-breeding for diversification of primary gene pool and genetic enhancement of grain legumes

The narrow genetic base of cultivars coupled with low utilization of genetic resources are the major factors limiting grain legume production and productivity globally. Exploitation of new and diverse sources of variation is needed for the genetic enhancement of grain legumes. Wild relatives with enhanced levels of resistance/tolerance to multiple stresses provide important sources of genetic diversity for crop improvement. However, their exploitation for cultivar improvement is limited by cross-incompatibility barriers and linkage drags. Pre-breeding provides a unique opportunity, through the introgression of desirable genes from wild germplasm into genetic backgrounds readily used by the breeders with minimum linkage drag, to overcome this. Pre-breeding activities using promising landraces, wild relatives, and popular cultivars have been initiated at International Crops Research Institute for the Semi-Arid Tropics (ICRISAT) to develop new gene pools in chickpea, pigeonpea, and groundnut with a high frequency of useful genes, wider adaptability, and a broad genetic base. The availability of molecular markers will greatly assist in reducing linkage drags and increasing the efficiency of introgression in pre-breeding programs.

Association mapping of maturity and plant height using SNP markers with the sorghum mini core collection

Plant height and maturity are two critical traits in sorghum breeding. To develop molecular tools and to identify genes underlying the traits for molecular breeding, we developed 14,739 SNP markers used to genotype the complete sorghum [Sorghum bicolor (L.) Moench] mini core collection. The collection was evaluated in four rainy and three post-rainy season environments for plant height and maturity. Association analysis identified six marker loci linked to height and ten to maturity in at least two environments with at least two SNPs in each locus. Of these, 14 were in close proximity to previously mapped height/maturity QTL in sorghum. Candidate genes for maturity or plant height close to the marker loci include a sugar transporter (SbSUC9), an auxin response factor (SbARF3), an FLC and FT regulator (SbMED12), and a photoperiod response gene (SbPPR1) for maturity and peroxidase 53, and an auxin transporter (SbLAX4) for plant height. Linkage disequilibrium analysis showed that SbPPR1 and SbARF3 were in regions with reduced sequence variation among early-maturing accessions, suggestive of past purifying selection. We also found a linkage disequilibrium block that existed only among the accessions with short plant height in rainy season environments. The block contains a gene homologous to the Arabidopsis flowering time gene, LUMINIDEPENDENS (LD). Functional LD promotes early maturity while mutation delays maturity, affecting plant height. Previous studies also found reduced sequence variations within this gene. These newly-mapped SNP markers will facilitate further efforts to identify plant height or maturity genes in sorghum.

Achievements and prospects of genomics-assisted breeding in three legume crops of the semi-arid tropics

Advances in next-generation sequencing and genotyping technologies have enabled generation of large-scale genomic resources such as molecular markers, transcript reads and BAC-end sequences (BESs) in chickpea, pigeonpea and groundnut, three major legume crops of the semi-arid tropics. Comprehensive transcriptome assemblies and genome sequences have either been developed or underway in these crops. Based on these resources, dense genetic maps, QTL maps as well as physical maps for these legume species have also been developed. As a result, these crops have graduated from 'orphan' or 'less-studied' crops to 'genomic resources rich' crops. This article summarizes the above-mentioned advances in genomics and genomics-assisted breeding applications in the form of marker-assisted selection (MAS) for hybrid purity assessment in pigeonpea; marker-assisted backcrossing (MABC) for introgressing QTL region for drought-tolerance related traits, Fusarium wilt (FW) resistance and Ascochyta blight (AB) resistance in chickpea; late leaf spot (LLS), leaf rust and nematode resistance in groundnut. We critically present the case of use of other modern breeding approaches like marker-assisted recurrent selection (MARS) and genomic selection (GS) to utilize the full potential of genomics-assisted breeding for developing superior cultivars with enhanced tolerance to various environmental stresses. In addition, this article recommends the use of advanced-backcross (AB-backcross) breeding and development of specialized populations such as multi-parents advanced generation intercross (MAGIC) for creating new variations that will help in developing superior lines with broadened genetic base. In summary, we propose the use of integrated genomics and breeding approach in these legume crops to enhance crop productivity in marginal environments ensuring food security in developing countries.

Phenotyping chickpeas and pigeonpeas for adaptation to drought

The chickpea and pigeonpea are protein-rich grain legumes used for human consumption in many countries. Grain yield of these crops is low to moderate in the semi-arid tropics with large variation due to high GxE interaction. In the Indian subcontinent chickpea is grown in the post-rainy winter season on receding soil moisture, and in other countries during the cool and dry post winter or spring seasons. The pigeonpea is sown during rainy season which flowers and matures in post-rainy season. The rainy months are hot and humid with diurnal temperature varying between 25 and 35°C (maximum) and 20 and 25°C (minimum) with an erratic rainfall. The available soil water during post-rainy season is about 200-250 mm which is bare minimum to meet the normal evapotranspiration. Thus occurrence of drought is frequent and at varying degrees. To enhance productivity of these crops cultivars tolerant to drought need to be developed. ICRISAT conserves a large number of accessions of chickpea (>20,000) and pigeonpea (>15,000). However only a small proportion (<1%) has been used in crop improvement programs mainly due to non-availability of reliable information on traits of economic importance. To overcome this, core and mini core collections (10% of core, 1% of entire collection) have been developed. Using the mini core approach, trait-specific donor lines were identified for agronomic, quality, and stress related traits in both crops. Composite collections were developed both in chickpea (3000 accessions) and pigeonpea (1000 accessions), genotyped using SSR markers and genotype based reference sets of 300 accessions selected for each crop. Screening methods for different drought-tolerant traits such as early maturity (drought escape), large and deep root system, high water-use efficiency, smaller leaflets, reduced canopy temperature, carbon isotope discrimination, high leaf chlorophyll content (drought avoidance), and breeding strategies for improving drought tolerance have been discussed.

Ascochyta blight of chickpea (Cicer arietinum L.): a review of biology, pathogenicity, and disease management

Ascochyta blight (AB), caused by Ascochyta rabiei is a major disease of chickpea (Cicer arietinum L.), especially in areas where cool, cloudy, and humid weather persists during the crop season. Several epidemics of AB causing complete yield loss have been reported. The fungus mainly survives between seasons through infected seed and in infected crop debris. Despite extensive pathological and molecular studies, the nature and extent of pathogenic variability in A. rabiei have not been clearly established. Accumulation of phenols, phytoalexins (medicarpin and maackiain), and hydrolytic enzymes has been associated with host-plant resistance (HPR). Seed treatment and foliar application of fungicides are commonly recommended for AB management, but further information on biology and survival of A. rabiei is needed to devise more effective management strategies. Recent studies on inheritance of AB resistance indicate that several quantitative trait loci (QTLs) control resistance. In this paper we review the biology of A. rabiei, HPR, and management options, with an emphasis on future research priorities.