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

CaLAP1 and CaLAP2 orchestrate anthocyanin biosynthesis in the seed coat of Cicer arietinum

MAIN CONCLUSION

Our findings shed light on the regulation of anthocyanin and proanthocyanidin biosynthesis in chickpea seed coats. Expression of R2R3-MYB transcription factors CaLAP1 and CaLAP2 enhanced the anthocyanins and proanthocyanidins content in chickpea. The seed coat color is a major economic trait in leguminous crop chickpea (Cicer arietinum). Anthocyanins and proanthocyanidins (PAs) are two classes of flavonoids that mainly contribute to the flower, seed coat and color of Desi chickpea cultivars. Throughout the land plant lineage, the accumulation of anthocyanins and PAs is regulated by MYB and bHLH transcription factors (TFs), which form an MBW (MYB, bHLH, and WD40) complex. Here, we report two R2R3-MYB TFs in chickpea belonging to the anthocyanin-specific subgroup-6, CaLAP1 (Legume Anthocyanin Production 1), and CaLAP2 (Legume Anthocyanin Production 2), which are mainly expressed in the flowers and developmental stages of the seeds. CaLAP1 and CaLAP2 interact with TT8-like CabHLH1 and WD40, forming the MBW complex, and bind to the promoter sequences of anthocyanin- and PA biosynthetic genes CaCHS6, CaDFR2, CaANS, and CaANR, leading to anthocyanins and PA accumulation in the seed coat of chickpea. Moreover, these CaLAPs partially complement the anthocyanin-deficient phenotype in the Arabidopsis thaliana sextuple mutant seedlings. Overexpression of CaLAPs in chickpea resulted in significantly higher expression of anthocyanin and PA biosynthetic genes leading to a darker seed coat color with higher accumulation of anthocyanin and PA. Our findings show that CaLAPs positively modulate anthocyanin and PA content in seed coats, which might influence plant development and resistance to various biotic and abiotic stresses.

Genetic dissection of domestication traits in interspecific chickpea populations

Chickpea (Cicer arietinum) is a pulse crop that provides an integral source of nutrition for human consumption. The close wild relatives Cicer reticulatum and Cicer echinospermum harbor untapped genetic diversity that can be exploited by chickpea breeders to improve domestic varieties. Knowledge of genomic loci that control important chickpea domestication traits will expedite the development of improved chickpea varieties derived from interspecific crosses. Therefore, we set out to identify genomic loci underlying key chickpea domestication traits by both association and quantitative trait locus (QTL) mapping using interspecific F2 populations. Diverse phenotypes were recorded for various agronomic traits. A total of 11 high-confidence markers were detected on chromosomes 1, 3, and 7 by both association and QTL mapping; these were associated with growth habit, flowering time, and seed traits. Furthermore, we identified candidate genes linked to these markers, which advanced our understanding of the genetic basis of domestication traits and validated known genes such as the FLOWERING LOCUS gene cluster that regulates flowering time. Collectively, this study has elucidated the genetic basis of chickpea domestication traits, which can facilitate the development of superior chickpea varieties.

Functional allele of a MATE gene selected during domestication modulates seed color in chickpea

Seed color is one of the key target traits of domestication and artificial selection in chickpeas due to its implications on consumer preference and market value. The complex seed color trait has been well dissected in several crop species; however, the genetic mechanism underlying seed color variation in chickpea remains poorly understood. Here, we employed an integrated genomics strategy involving QTL mapping, high-density mapping, map-based cloning, association analysis, and molecular haplotyping in an inter-specific RIL mapping population, association panel, wild accessions, and introgression lines (ILs) of Cicer gene pool. This delineated a MATE gene, CaMATE23, encoding a Transparent Testa (TT) and its natural allele (8-bp insertion) and haplotype underlying a major QTL governing seed color on chickpea chromosome 4. Signatures of selective sweep and a strong purifying selection reflected that CaMATE23, especially its 8-bp insertion natural allelic variant, underwent selection during chickpea domestication. Functional investigations revealed that the 8-bp insertion containing the third cis-regulatory RY-motif element in the CaMATE23 promoter is critical for enhanced binding of CaFUSCA3 transcription factor, a key regulator of seed development and flavonoid biosynthesis, thereby affecting CaMATE23 expression and proanthocyanidin (PA) accumulation in the seed coat to impart varied seed color in chickpea. Consequently, overexpression of CaMATE23 in Arabidopsis tt12 mutant partially restored the seed color phenotype to brown pigmentation, ascertaining its functional role in PA accumulation in the seed coat. These findings shed new light on the seed color regulation and evolutionary history, and highlight the transcriptional regulation of CaMATE23 by CaFUSCA3 in modulating seed color in chickpea. The functionally relevant InDel variation, natural allele, and haplotype from CaMATE23 are vital for translational genomic research, including marker-assisted breeding, for developing chickpea cultivars with desirable seed color that appeal to consumers and meet global market demand.

Genome-wide association study as a powerful tool for dissecting competitive traits in legumes

Legumes are extremely valuable because of their high protein content and several other nutritional components. The major challenge lies in maintaining the quantity and quality of protein and other nutritional compounds in view of climate change conditions. The global need for plant-based proteins has increased the demand for seeds with a high protein content that includes essential amino acids. Genome-wide association studies (GWAS) have evolved as a standard approach in agricultural genetics for examining such intricate characters. Recent development in machine learning methods shows promising applications for dimensionality reduction, which is a major challenge in GWAS. With the advancement in biotechnology, sequencing, and bioinformatics tools, estimation of linkage disequilibrium (LD) based associations between a genome-wide collection of single-nucleotide polymorphisms (SNPs) and desired phenotypic traits has become accessible. The markers from GWAS could be utilized for genomic selection (GS) to predict superior lines by calculating genomic estimated breeding values (GEBVs). For prediction accuracy, an assortment of statistical models could be utilized, such as ridge regression best linear unbiased prediction (rrBLUP), genomic best linear unbiased predictor (gBLUP), Bayesian, and random forest (RF). Both naturally diverse germplasm panels and family-based breeding populations can be used for association mapping based on the nature of the breeding system (inbred or outbred) in the plant species. MAGIC, MCILs, RIAILs, NAM, and ROAM are being used for association mapping in several crops. Several modifications of NAM, such as doubled haploid NAM (DH-NAM), backcross NAM (BC-NAM), and advanced backcross NAM (AB-NAM), have also been used in crops like rice, wheat, maize, barley mustard, etc. for reliable marker-trait associations (MTAs), phenotyping accuracy is equally important as genotyping. Highthroughput genotyping, phenomics, and computational techniques have advanced during the past few years, making it possible to explore such enormous datasets. Each population has unique virtues and flaws at the genomics and phenomics levels, which will be covered in more detail in this review study. The current investigation includes utilizing elite breeding lines as association mapping population, optimizing the choice of GWAS selection, population size, and hurdles in phenotyping, and statistical methods which will analyze competitive traits in legume breeding.

The developmental dynamics in cool season legumes with focus on chickpea

Chickpea is one of the most widely consumed grain legume world-wide. Advances in next-generation sequencing and genomics tools have led to genetic dissection and identification of potential candidate genes regulating agronomic traits in chickpea. However, the developmental particularities and its potential in reforming the yield and nutritional value remain largely unexplored. Studies in crops such as rice, maize, tomato and pea have highlighted the contribution of key regulator of developmental events in yield related traits. A comprehensive knowledge on the development aspects of a crop can pave way for new vistas to explore. Pea and Medicago are the close relatives of genus Cicer and the basic developmental events in these legumes are similar. However, there are some distinct developmental features in chickpea which hold potential for future crop improvement endeavours. The global chickpea germplasm encompasses wide range of diversities in terms of morphology at both vegetative and reproductive stages. There is an immediate need for understanding the genetic and molecular basis of this diversity and utilizing them for the yield contributing trait improvement. The review discusses some of the key developmental events which have potential in yield enhancement and the lessons which can be learnt from model legumes in this regard.

The R2R3-MYB gene family in Cicer arietinum: genome-wide identification and expression analysis leads to functional characterization of proanthocyanidin biosynthesis regulators in the seed coat

We identified 119 typical CaMYB encoding genes and reveal the major components of the proanthocyanidin regulatory network. CaPARs emerged as promising targets for genetic engineering toward improved agronomic traits in C. arietinum. Chickpea (Cicer arietinum) is among the eight oldest crops and has two main types, i.e., desi and kabuli, whose most obvious difference is the color of their seeds. We show that this color difference is due to differences in proanthocyanidin content of seed coats. Using a targeted approach, we performed in silico analysis, metabolite profiling, molecular, genetic, and biochemical studies to decipher the transcriptional regulatory network involved in proanthocyanidin biosynthesis in the seed coat of C. arietinum. Based on the annotated C. arietinum reference genome sequence, we identified 119 typical CaMYB encoding genes, grouped in 32 distinct clades. Two CaR2R3-MYB transcription factors, named CaPAR1 and CaPAR2, clustering with known proanthocyanidin regulators (PARs) were identified and further analyzed. The expression of CaPAR genes correlated well with the expression of the key structural proanthocyanidin biosynthesis genes CaANR and CaLAR and with proanthocyanidin levels. Protein-protein interaction studies suggest the in vivo interaction of CaPAR1 and CaPAR2 with the bHLH-type transcription factor CaTT8. Co-transfection analyses using Arabidopsis thaliana protoplasts showed that the CaPAR proteins form a MBW complex with CaTT8 and CaTTG1, able to activate the promoters of CaANR and CaLAR in planta. Finally, transgenic expression of CaPARs in the proanthocyanidin-deficient A. thaliana mutant tt2-1 leads to complementation of the transparent testa phenotype. Taken together, our results reveal main components of the proanthocyanidin regulatory network in C. arietinum and suggest that CaPARs are relevant targets of genetic engineering toward improved agronomic traits.

Genetic Augmentation of Legume Crops Using Genomic Resources and Genotyping Platforms for Nutritional Food Security

Recent advances in next generation sequencing (NGS) technologies have led the surge of genomic resources for the improvement legume crops. Advances in high throughput genotyping (HTG) and high throughput phenotyping (HTP) enable legume breeders to improve legume crops more precisely and efficiently. Now, the legume breeder can reshuffle the natural gene combinations of their choice to enhance the genetic potential of crops. These genomic resources are efficiently deployed through molecular breeding approaches for genetic augmentation of important legume crops, such as chickpea, cowpea, pigeonpea, groundnut, common bean, lentil, pea, as well as other underutilized legume crops. In the future, advances in NGS, HTG, and HTP technologies will help in the identification and assembly of superior haplotypes to tailor the legume crop varieties through haplotype-based breeding. This review article focuses on the recent development of genomic resource databases and their deployment in legume molecular breeding programmes to secure global food security.

Harnessing the hidden allelic diversity of wild Cicer to accelerate genomics-assisted chickpea crop improvement

Chickpea, commonly called Bengal gram or Garbanzo bean, faces a productivity crisis around the globe due to numerous biotic and abiotic stresses. The eroded genetic base of the cultivated Cicer gene pool is becoming a significant bottleneck in developing stress-resilient chickpea cultivars. In this scenario, the crop wild relatives (CWR) of chickpea, with the useful genomic wealth of their wild adaptation, give a ray of hope to improve the genetic background of the cultivated Cicer gene pool. To extrapolate these unearthed genomic diversities of wild, we require a thorough understanding of the pre-historic domestication episodes that are changing their shape with the expansion of the available scientific evidence. Keeping aforesaid in view, the current review article provides a glimpsed overview on several efforts done so far to reveal the mysterious origin and evolution of the Cicer gene pool, along with the constraints in their utilization for chickpea crop improvement. It encapsulates various stress-resilient CWR of chickpea and their use in several pre-breeding programs to develop numerous breeding populations for crop genetic enhancement. Further, this review will recapitulate the significant contributions of structural, functional and comparative genomics, pan-genomics and diverse genomics-assisted breeding strategy in dissecting the untapped trait-specific allelic/gene diversity and domestication pattern behind the CWR of chickpea, along with their potential and promises. We expect the newly explored genetic variations may be used in the breeding programs for re-wilding the cultigens' genomic background to open a new avenue for genetic gain and crop improvement capacity of chickpea.

Genome resequencing reveals DNA polymorphisms associated with seed size/weight determination in chickpea

Crop productivity in legumes is determined by number and size/weight of seeds. To understand the genetic basis of seed size/weight in chickpea, we performed genome resequencing of 13 small- and 5 large-seeded genotypes using Illumina platform. Single nucleotide polymorphisms (SNPs) and insertions/deletions (InDels) differentiating small- and large-seeded genotypes were identified. A total of 17,902 SNPs and 2594 InDels located in promoter and/or coding regions that may contribute to seed size/weight were detected. Of these, 266 SNPs showed significant association with seed size/weight trait. Twenty-three genes including those involved in cell growth/division, encoding transcription factors and located within QTLs associated with seed size/weight harbored SNPs within transcription factor binding motif(s) and/or coding region. The non-synonymous SNPs were found to affect the mutational sensitivity and stability of the encoded proteins. Overall, we provided a high-quality SNP map for large-scale genotyping applications and identified candidate genes that determine seed size/weight in chickpea.

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.

Is Chickpea a Potential Substitute for Soybean? Phenolic Bioactives and Potential Health Benefits

Legume seeds are rich sources of protein, fiber, and minerals. In addition, their phenolic compounds as secondary metabolites render health benefits beyond basic nutrition. Lowering apolipoprotein B secretion from HepG2 cells and decreasing the level of low-density lipoprotein (LDL)-cholesterol oxidation are mechanisms related to the prevention of cardiovascular diseases (CVD). Likewise, low-level chronic inflammation and related disorders of the immune system are clinical predictors of cardiovascular pathology. Furthermore, DNA-damage signaling and repair are crucial pathways to the etiology of human cancers. Along CVD and cancer, the prevalence of obesity and diabetes is constantly increasing. Screening the ability of polyphenols in inactivating digestive enzymes is a good option in pre-clinical studies. In addition, in vivo studies support the role of polyphenols in the prevention and/or management of diabetes and obesity. Soybean, a well-recognized source of phenolic isoflavones, exerts health benefits by decreasing oxidative stress and inflammation related to the above-mentioned chronic ailments. Similar to soybeans, chickpeas are good sources of nutrients and phenolic compounds, especially isoflavones. This review summarizes the potential of chickpea as a substitute for soybean in terms of health beneficial outcomes. Therefore, this contribution may guide the industry in manufacturing functional foods and/or ingredients by using an undervalued feedstock.

Co-localization of genomic regions associated with seed morphology and composition in a desi chickpea (Cicer arietinum L.) population varying in seed protein concentration

Major QTL on LG 1 and 3 control seed filling and seed coat development, thereby affecting seed shape, size, color, composition and weight, key determinants of crop yield and quality. A chickpea (Cicer arietinum L.) population consisting of 189 recombinant inbred lines (RILs) derived from a cross between medium-protein ICC 995 and high-protein ICC 5912 genotypes of the desi market class was analyzed for seed properties. Seed from the parental lines and RILs was produced in four different environments for determination of seed shape (SS), 100-seed weight (100-SW), protein (PRO) and starch (STA) concentration. Polymorphic genetic markers for the population were identified by Genotyping by Sequencing and assembled into a 522.5 cM genetic map. Phenotype data from the different growth environments were analyzed by QTL mapping done by single and multi-environment analyses and in addition, single marker association mapping. The analyses identified in total 11 QTL, of which the most significant (P < 0.05) loci were located on LG 1 (q-1.1), LG 2 (q-2.1), LG 3 (q-3.2, q-3.3), LG 4 (q-4.2), and LG 5 (q-5.1). STA was mostly affected by q-1.1, which explained 19.0% of the phenotypic variance for the trait. The largest QTL effects were demonstrated by q-3.2 that explained 52.5% of the phenotypic variances for 100-SW, 44.3% for PRO, and 14.6% for SS. This locus was also highly associated with flower color (COL; 95.2% explained) and showed q-3.2 alleles from the ICC 5912 parent conferred the blue flower color and production of small, round seeds with relatively high protein concentration. Genes affecting seed filling at q-1.1 and seed coat development at q-3.2, respectively, were considered to underlie differences in seed composition and morphology in the RIL population.

Influence of the preparation process on the chemical composition and nutritional value of canned purée of kabuli and Apulian black chickpeas

Chickpeas are classified into two main commercial seed types: desi, and kabuli. Furthermore, the Apulian black is another chickpea type, less common, which has peculiar phenotypic and genetic features and is the object of an increasing attention by geneticists to avoid the risk of genetic erosion. A strategy to increase the consumption of pulses consists in proposing ready-to-eat gastronomic preparations which, however, must keep their natural features and nutritional value as intact as possible. In this paper the influence of the preparation process on the chemical composition and nutritional value of ready-to-eat canned purée of Apulian black chickpeas has been evaluated, in comparison with purée of kabuli chickpeas. Total dietary fiber content was high enough to consider the kabuli chickpea purée as "source of fiber", and the black chickpea purée as "high fiber", in accordance with the current European Regulation on nutrition claims. Along the preparation process, an increase in lipid content was observed. Protein content, instead, showed a different behaviour, i.e. increased in Apulian black chickpea purée and remained constant in kabuli chickpea purée.

The preparation process strongly influenced fatty acid composition. In particular, unsaturated fatty acids decreased in both Apulian black and kabuli chickpea purées, whereas saturated fatty acids significantly increased during processing.

Apulian black chickpeas are particularly rich of bioactive compounds, but the preparation process of purée caused a strong decrease of total carotenoids, anthocyanins and phenolic compounds. However, even after processing, this purée could still be a good source of bioactive compounds.

All these features make canned purée of chickpeas a healthy ready-to-eat food, which is at the same time rich in fiber and bioactive compounds, able to fulfill the time-saving needs of modern lifestyle. These findings could promote a greater use of Apulian black chickpeas and contribute to reduce the risk of genetic erosion.

Population structure and association analysis of heat stress relevant traits in chickpea (Cicer arietinum L.)

Understanding genetic diversity and population structure is prerequisite to broaden the cultivated base of any crop. In the current investigation, we report discovery of a total of 319 alleles by assaying 81 SSRs on 71 chickpea genotypes. The cluster analysis based on Jaccard coefficient and unweighted neighbor joining algorithm categorized all genotypes into two major clusters. Cultivars grown within the same agro-climatic zones were clustered together, whereas the remaining genotypes particularly advanced breeding lines and accessions assigned to another cluster. Population structure analysis separated the entire collection into two subpopulations (K = 2) and the clustering pattern remained in close agreement with those of distance-based methods. Importantly, we also discovered marker trait association for membrane stability index (MSI) and leaf chlorophyll content measured as SPAD chlorophyll meter reading (SCMR), the two important physiological parameters indicative of heat stress (HS) tolerance in chickpea. Association analysis using both general linear and mixed linear models of the mean phenotypic data of traits recorded in 2016 and 2017 uncovered significant association of NCPGR206 and H2L102 with the MSI trait. Likewise, SSR markers GA9, TR31 and TA113 exhibited significant association with SCMR trait. The genomic regions putatively linked with two traits may be investigated in greater detail to further improve knowledge about the genetic architecture of HS tolerance in chickpea.

Population structure and association analysis of heat stress relevant traits in chickpea (Cicer arietinum L.)

Understanding genetic diversity and population structure is prerequisite to broaden the cultivated base of any crop. In the current investigation, we report discovery of a total of 319 alleles by assaying 81 SSRs on 71 chickpea genotypes. The cluster analysis based on Jaccard coefficient and unweighted neighbor joining algorithm categorized all genotypes into two major clusters. Cultivars grown within the same agro-climatic zones were clustered together, whereas the remaining genotypes particularly advanced breeding lines and accessions assigned to another cluster. Population structure analysis separated the entire collection into two subpopulations (K = 2) and the clustering pattern remained in close agreement with those of distance-based methods. Importantly, we also discovered marker trait association for membrane stability index (MSI) and leaf chlorophyll content measured as SPAD chlorophyll meter reading (SCMR), the two important physiological parameters indicative of heat stress (HS) tolerance in chickpea. Association analysis using both general linear and mixed linear models of the mean phenotypic data of traits recorded in 2016 and 2017 uncovered significant association of NCPGR206 and H2L102 with the MSI trait. Likewise, SSR markers GA9, TR31 and TA113 exhibited significant association with SCMR trait. The genomic regions putatively linked with two traits may be investigated in greater detail to further improve knowledge about the genetic architecture of HS tolerance in chickpea.

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.

A Distinct Genetic Cluster in Cultivated Chickpea as Revealed by Genome-wide Marker Discovery and Genotyping

The accurate description of plant biodiversity is of utmost importance to efficiently address efforts in conservation genetics and breeding. Herein, we report the successful application of a genotyping-by-sequencing (GBS) approach in chickpea ( L.), resulting in the characterization of a cultivated germplasm collection with 3187 high-quality single nucleotide polymorphism (SNP) markers. Genetic structure inference, principal component analysis, and hierarchical clustering all indicated the identification of a genetic cluster corresponding to black-seeded genotypes traditionally cultivated in Southern Italy. Remarkably, this cluster was clearly distinct at both genetic and phenotypic levels from germplasm groups reflecting commercial chickpea classification into and seed types. Fixation index estimates for individual polymorphisms pointed out loci and genomic regions that might be of significance for the diversification of agronomic and commercial traits. Overall, our findings provide information on genetic relationships within cultivated chickpea and highlight a gene pool of great interest for the scientific community and chickpea breeding, which is limited by the low genetic diversity available in the primary gene pool.

Genetic dissection of plant growth habit in chickpea

A combinatorial genomics-assisted breeding strategy encompassing association analysis, genetic mapping and expression profiling is found most promising for quantitative dissection of complex traits in crop plants. The present study employed GWAS (genome-wide association study) using 24,405 SNPs (single nucleotide polymorphisms) obtained with genotyping-by-sequencing (GBS) of 92 sequenced desi and kabuli accessions of chickpea. This identified eight significant genomic loci associated with erect (E)/semi-erect (SE) vs. spreading (S)/semi-spreading (SS)/prostrate (P) plant growth habit (PGH) trait differentiation regardless of diverse desi and kabuli genetic backgrounds of chickpea. These associated SNPs in combination explained 23.8% phenotypic variation for PGH in chickpea. Five PGH-associated genes were validated successfully in E/SE and SS/S/P PGH-bearing parental accessions and homozygous individuals of three intra- and interspecific RIL (recombinant inbred line) mapping populations as well as 12 contrasting desi and kabuli chickpea germplasm accessions by selective genotyping through Sequenom MassARRAY. The shoot apical, inflorescence and floral meristems-specific expression, including upregulation (seven-fold) of five PGH-associated genes especially in germplasm accessions and homozygous RIL mapping individuals contrasting with E/SE PGH traits was apparent. Collectively, this integrated genomic strategy delineated diverse non-synonymous SNPs from five candidate genes with strong allelic effects on PGH trait variation in chickpea. Of these, two vernalization-responsive non-synonymous SNP alleles carrying SNF2 protein-coding gene and B3 transcription factor associated with PGH traits were found to be the most promising in chickpea. The SNP allelic variants associated with E/SE/SS/S PGH trait differentiation were exclusively present in all cultivated desi and kabuli chickpea accessions while wild species/accessions belonging to primary, secondary and tertiary gene pools mostly contained prostrate PGH-associated SNP alleles. This indicates strong adaptive natural/artificial selection pressure (Tajima's D 3.15 to 4.57) on PGH-associated target genomic loci during chickpea domestication. These vital leads thus have potential to decipher complex transcriptional regulatory gene function of PGH trait differentiation and for understanding the selective sweep-based PGH trait evolution and domestication pattern in cultivated and wild chickpea accessions adapted to diverse agroclimatic conditions. Collectively, the essential inputs generated will be of profound use in marker-assisted genetic enhancement to develop cultivars with desirable plant architecture of erect growth habit types in chickpea.

Potential Uses of Wild Germplasms of Grain Legumes for Crop Improvement

Challenged by population increase, climatic change, and soil deterioration, crop improvement is always a priority in securing food supplies. Although the production of grain legumes is in general lower than that of cereals, the nutritional value of grain legumes make them important components of food security. Nevertheless, limited by severe genetic bottlenecks during domestication and human selection, grain legumes, like other crops, have suffered from a loss of genetic diversity which is essential for providing genetic materials for crop improvement programs. Illustrated by whole-genome-sequencing, wild relatives of crops adapted to various environments were shown to maintain high genetic diversity. In this review, we focused on nine important grain legumes (soybean, peanut, pea, chickpea, common bean, lentil, cowpea, lupin, and pigeonpea) to discuss the potential uses of their wild relatives as genetic resources for crop breeding and improvement, and summarized the various genetic/genomic approaches adopted for these purposes.

A Combined Comparative Transcriptomic, Metabolomic, and Anatomical Analyses of Two Key Domestication Traits: Pod Dehiscence and Seed Dormancy in Pea (Pisum sp.)

The origin of the agriculture was one of the turning points in human history, and a central part of this was the evolution of new plant forms, domesticated crops. Seed dispersal and germination are two key traits which have been selected to facilitate cultivation and harvesting of crops. The objective of this study was to analyze anatomical structure of seed coat and pod, identify metabolic compounds associated with water-impermeable seed coat and differentially expressed genes involved in pea seed dormancy and pod dehiscence. Comparative anatomical, metabolomics, and transcriptomic analyses were carried out on wild dormant, dehiscent Pisum elatius (JI64, VIR320) and cultivated, indehiscent Pisum sativum non-dormant (JI92, Cameor) and recombinant inbred lines (RILs). Considerable differences were found in texture of testa surface, length of macrosclereids, and seed coat thickness. Histochemical and biochemical analyses indicated genotype related variation in composition and heterogeneity of seed coat cell walls within macrosclereids. Liquid chromatography-electrospray ionization/mass spectrometry and Laser desorption/ionization-mass spectrometry of separated seed coats revealed significantly higher contents of proanthocyanidins (dimer and trimer of gallocatechin), quercetin, and myricetin rhamnosides and hydroxylated fatty acids in dormant compared to non-dormant genotypes. Bulk Segregant Analysis coupled to high throughput RNA sequencing resulted in identification of 770 and 148 differentially expressed genes between dormant and non-dormant seeds or dehiscent and indehiscent pods, respectively. The expression of 14 selected dormancy-related genes was studied by qRT-PCR. Of these, expression pattern of four genes: porin (MACE-S082), peroxisomal membrane PEX14-like protein (MACE-S108), 4-coumarate CoA ligase (MACE-S131), and UDP-glucosyl transferase (MACE-S139) was in agreement in all four genotypes with Massive analysis of cDNA Ends (MACE) data. In case of pod dehiscence, the analysis of two candidate genes (SHATTERING and SHATTERPROOF) and three out of 20 MACE identified genes (MACE-P004, MACE-P013, MACE-P015) showed down-expression in dorsal and ventral pod suture of indehiscent genotypes. Moreover, MACE-P015, the homolog of peptidoglycan-binding domain or proline-rich extensin-like protein mapped correctly to predicted Dpo1 locus on PsLGIII. This integrated analysis of the seed coat in wild and cultivated pea provides new insight as well as raises new questions associated with domestication and seed dormancy and pod dehiscence.

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.