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

Haplotypes of ATP-Binding Cassette CaABCC6 in Chickpea from Kazakhstan Are Associated with Salinity Tolerance and Leaf Necrosis via Oxidative Stress

Salinity tolerance was studied in chickpea accessions from a germplasm collection and in cultivars from Kazakhstan. After NaCl treatment, significant differences were found between genotypes, which could be arranged into three groups. Those that performed poorest were found in group 1, comprising five ICC accessions with the lowest chlorophyll content, the highest leaf necrosis (LN), Na+ accumulation, malondialdehyde (MDA) content, and a low glutathione ratio GSH/GSSG. Two cultivars, Privo-1 and Tassay, representing group 2, were moderate in these traits, while the best performance was for group 3, containing two other cultivars, Krasnokutsky-123 and Looch, which were found to have mostly green plants and an exact opposite pattern of traits. Marker-trait association (MTA) between 6K DArT markers and four traits (LN, Na+, MDA, and GSH/GSSG) revealed the presence of four possible candidate genes in the chickpea genome that may be associated with the three groups. One gene, ATP-binding cassette, CaABCC6, was selected, and three haplotypes, A, D1, and D2, were identified in plants from the three groups. Two of the most salt-tolerant cultivars from group 3 were found to have haplotype D2 with a novel identified SNP. RT-qPCR analysis confirmed that this gene was strongly expressed after NaCl treatment in the parental- and breeding-line plants of haplotype D2. Mass spectrometry of seed proteins showed a higher accumulation of glutathione reductase and S-transferase, but not peroxidase, in the D2 haplotype. In conclusion, the CaABCC6 gene was hypothesized to be associated with a better response to oxidative stress via glutathione metabolism, while other candidate genes are likely involved in the control of chlorophyll content and Na+ accumulation.

Zinc finger knuckle genes are associated with tolerance to drought and dehydration in chickpea (Cicer arietinum L.)

Chickpea (Cicer arietinum L.) is a very important food legume and needs improved drought tolerance for higher seed production in dry environments. The aim of this study was to determine diversity and genetic polymorphism in zinc finger knuckle genes with CCHC domains and their functional analysis for practical improvement of chickpea breeding. Two CaZF-CCHC genes, Ca04468 and Ca07571, were identified as potentially important candidates associated with plant responses to drought and dehydration. To study these genes, various methods were used including Sanger sequencing, DArT (Diversity array technology) and molecular markers for plant genotyping, gene expression analysis using RT-qPCR, and associations with seed-related traits in chickpea plants grown in field trials. These genes were studied for genetic polymorphism among a set of chickpea accessions, and one SNP was selected for further study from four identified SNPs between the promoter regions of each of the two genes. Molecular markers were developed for the SNP and verified using the ASQ and CAPS methods. Genotyping of parents and selected breeding lines from two hybrid populations, and SNP positions on chromosomes with haplotype identification, were confirmed using DArT microarray analysis. Differential expression profiles were identified in the parents and the hybrid populations under gradual drought and rapid dehydration. The SNP-based genotypes were differentially associated with seed weight per plant but not with 100 seed weight. The two developed and verified SNP molecular markers for both genes, Ca04468 and Ca07571, respectively, could be used for marker-assisted selection in novel chickpea cultivars with improved tolerance to drought and dehydration.

Meta QTL analysis for dissecting abiotic stress tolerance in chickpea

BACKGROUND

Chickpea is prone to many abiotic stresses such as heat, drought, salinity, etc. which cause severe loss in yield. Tolerance towards these stresses is quantitative in nature and many studies have been done to map the loci influencing these traits in different populations using different markers. This study is an attempt to meta-analyse those reported loci projected over a high-density consensus map to provide a more accurate information on the regions influencing heat, drought, cold and salinity tolerance in chickpea.

RESULTS

A meta-analysis of QTL reported to be responsible for tolerance to drought, heat, cold and salinity stress tolerance in chickpeas was done. A total of 1512 QTL responsible for the concerned abiotic stress tolerance were collected from literature, of which 1189 were projected on a chickpea consensus genetic map. The QTL meta-analysis predicted 59 MQTL spread over all 8 chromosomes, responsible for these 4 kinds of abiotic stress tolerance in chickpea. The physical locations of 23 MQTL were validated by various marker-trait associations and genome-wide association studies. Out of these reported MQTL, CaMQAST1.1, CaMQAST4.1, CaMQAST4.4, CaMQAST7.8, and CaMQAST8.2 were suggested to be useful for different breeding approaches as they were responsible for high per cent variance explained (PVE), had small intervals and encompassed a large number of originally reported QTL. Many putative candidate genes that might be responsible for directly or indirectly conferring abiotic stress tolerance were identified in the region covered by 4 major MQTL- CaMQAST1.1, CaMQAST4.4, CaMQAST7.7, and CaMQAST6.4, such as heat shock proteins, auxin and gibberellin response factors, etc. CONCLUSION: The results of this study should be useful for the breeders and researchers to develop new chickpea varieties which are tolerant to drought, heat, cold, and salinity stresses.

Integration of genome-wide association studies, metabolomics, and transcriptomics reveals phenolic acid- and flavonoid-associated genes and their regulatory elements under drought stress in rapeseed flowers

Introduction

Biochemical and metabolic processes help plants tolerate the adverse effects of drought. In plants accumulating bioactive compounds, understanding the genetic control of the biosynthesis of biochemical pathways helps the discovery of candidate gene (CG)-metabolite relationships.

Methods

The metabolic profile of flowers in 119 rapeseed (Brassica napus) accessions was assessed over two irrigation treatments, one a well-watered (WW) condition and the other a drought stress (DS) regime. We integrated information gained from 52,157 single-nucleotide polymorphism (SNP) markers, metabolites, and transcriptomes to identify linked SNPs and CGs responsible for the genetic control of flower phenolic compounds and regulatory elements.

Results

In a genome-wide association study (GWAS), of the SNPs tested, 29,310 SNPs were qualified to assess the population structure and linkage disequilibrium (LD), of which several SNPs for radical scavenging activity (RSA) and total flavanol content (TFLC) were common between the two irrigation conditions and pleiotropic SNPs were found for chlorogenic and coumaric acids content. The principal component analysis (PCA) and stepwise regression showed that chlorogenic acid and epicatechin in WW and myricetin in DS conditions were the most important components for RSA. The hierarchical cluster analysis (HCA) showed that vanillic acid, myricetin, gallic acid, and catechin were closely associated in both irrigation conditions. Analysis of GWAS showed that 60 CGs were identified, of which 18 were involved in stress-induced pathways, phenylpropanoid pathway, and flavonoid modifications. Of the CGs, PAL1, CHI, UGT89B1, FLS3, CCR1, and CYP75B137 contributed to flavonoid biosynthetic pathways. The results of RNA sequencing (RNA-seq) revealed that the transcript levels of PAL, CHI, and CYP75B137 known as early flavonoid biosynthesis-related genes and FLS3, CCR1, and UGT89B1 related to the later stages were increased during drought conditions. The transcription factors (TFs) NAC035 and ERF119 related to flavonoids and phenolic acids were upregulated under drought conditions.

Discussion

These findings expand our knowledge on the response mechanisms to DS, particularly regarding the regulation of key phenolic biosynthetic genes in rapeseed. Our data also provided specific linked SNPs for marker-assisted selection (MAS) programs and CGs as resources toward realizing metabolomics-associated breeding of rapeseed.

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.

Natural alleles of Mediator subunit genes modulate plant height in chickpea

Plant height (PH) is an important plant architectural trait targeted during Green Revolution to enhance crop yields. Identification of genes and natural alleles governing plant height without compromising agronomic performance can fill the lacuna of knowledge connecting ideal plant architecture with maximum achievable yield in chickpea. Through coherent strategy involving genome‐wide association study, QTL/fine mapping, map‐based cloning, molecular haplotyping, and downstream functional genomics, the current study identified two Mediator subunit genes namely, CaMED23 and CaMED5b and their derived natural alleles/haplotypes underlying the major QTLs and trans‐acting eQTLs regulating plant height in chickpea. Differential accumulation of haplotype‐specific transcripts of these two Mediator genes in corresponding haplotype‐introgressed near‐isogenic lines (NILs) correlates negatively with the plant height trait. Quantitative as well as qualitative estimation based on histology, scanning electron microscopy, and histochemical assay unraveled the reduced lengths and cell sizes of internodes along with compromised lignin levels in dwarf/semi‐dwarf chickpea NILs introgressed with superior CaMED23 and CaMED5b gene haplotypes. This observation, supported by global transcriptome profiling‐based diminished expression of various phenylpropanoid pathway genes upstream of lignin biosynthesis in dwarf/semi‐dwarf NILs, essentially links plant height with lignin accumulation. The identified molecular signatures in the Mediator subunit genes can be efficiently utilized to develop desirable dwarf/semi‐dwarf‐type chickpea cultivars without affecting their yield per plant via modulating lignin/phenylpropanoid biosynthesis.

Fatty acid composition and genome-wide associations of a chickpea (Cicer arietinum L.) diversity panel for biofortification efforts

Chickpea is a nutritionally dense pulse crop with high levels of protein, carbohydrates, micronutrients and low levels of fats. Chickpea fatty acids are associated with a reduced risk of obesity, blood cholesterol, and cardiovascular diseases in humans. We measured four primary chickpea fatty acids; palmitic acid (PA), linoleic acid (LA), alpha-linolenic acid (ALA), and oleic acid (OA), which are crucial for human health and plant stress responses in a chickpea diversity panel with 256 accessions (Kabuli and desi types). A wide concentration range was found for PA (450.7-912.6 mg/100 g), LA (1605.7-3459.9 mg/100 g), ALA (416.4-864.5 mg/100 g), and OA (1035.5-1907.2 mg/100 g). The percent recommended daily allowances also varied for PA (3.3-6.8%), LA (21.4-46.1%), ALA (34.7-72%), and OA (4.3-7.9%). Weak correlations were found among fatty acids. Genome-wide association studies (GWAS) were conducted using genotyping-by-sequencing data. Five significant single nucleotide polymorphisms (SNPs) were identified for PA. Admixture population structure analysis revealed seven subpopulations based on ancestral diversity in this panel. This is the first reported study to characterize fatty acid profiles across a chickpea diversity panel and perform GWAS to detect associations between genetic markers and concentrations of selected fatty acids. These findings demonstrate biofortification of chickpea fatty acids is possible using conventional and genomic breeding techniques, to develop superior cultivars with better fatty acid profiles for improved human health and plant stress responses.

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.

Genetic dissection of thousand-seed weight in linseed (Linum usitatissimum L.) using multi-locus genome-wide association study

Flaxseed/linseed is an important oilseed crop having applications in the food, nutraceutical, and paint industry. Seed weight is one of the most crucial determinants of seed yield in linseed. Here, quantitative trait nucleotides (QTNs) associated with thousand-seed weight (TSW) have been identified using multi-locus genome-wide association study (ML-GWAS). Field evaluation was carried out in five environments in multi-year-location trials. SNP genotyping information of the AM panel of 131 accessions comprising 68,925 SNPs was employed for ML-GWAS. From the six ML-GWAS methods employed, five methods helped identify a total of 84 unique significant QTNs for TSW. QTNs identified in ≥ 2 methods/environments were designated as stable QTNs. Accordingly, 30 stable QTNs have been identified for TSW accounting up to 38.65% trait variation. Alleles with positive effect on trait were analyzed for 12 strong QTNs with r ² ≥ 10.00%, which showed significant association of specific alleles with higher trait value in three or more environments. A total of 23 candidate genes have been identified for TSW, which included B3 domain-containing transcription factor, SUMO-activating enzyme, protein SCARECROW, shaggy-related protein kinase/BIN2, ANTIAUXIN-RESISTANT 3, RING-type E3 ubiquitin transferase E4, auxin response factors, WRKY transcription factor, and CBS domain-containing protein. In silico expression analysis of candidate genes was performed to validate their possible role in different stages of seed development process. The results from this study provide significant insight and elevate our understanding on genetic architecture of TSW trait in linseed.

A superior gene allele involved in abscisic acid signaling enhances drought tolerance and yield in chickpea

Identifying potential molecular tags for drought tolerance is essential for achieving higher crop productivity under drought stress. We employed an integrated genomics-assisted breeding and functional genomics strategy involving association mapping, fine mapping, map-based cloning, molecular haplotyping and transcript profiling in the introgression lines (ILs)- and near isogenic lines (NILs)-based association panel and mapping population of chickpea (Cicer arietinum). This combinatorial approach delineated a bHLH (basic helix-loop-helix) transcription factor, CabHLH10 (Cicer arietinum bHLH10) underlying a major QTL, along with its derived natural alleles/haplotypes governing yield traits under drought stress in chickpea. CabHLH10 binds to a cis-regulatory G-box promoter element to modulate the expression of RD22 (responsive to desiccation 22), a drought/abscisic acid (ABA)-responsive gene (via a trans-expression QTL), and two strong yield-enhancement photosynthetic efficiency (PE) genes. This, in turn, upregulates other downstream drought-responsive and ABA signaling genes, as well as yield-enhancing PE genes, thus increasing plant adaptation to drought with reduced yield penalty. We showed that a superior allele of CabHLH10 introgressed into the NILs improved root and shoot biomass and PE, thereby enhancing yield and productivity during drought without compromising agronomic performance. Furthermore, overexpression of CabHLH10 in chickpea and Arabidopsis (Arabidopsis thaliana) conferred enhanced drought tolerance by improving root and shoot agro-morphological traits. These findings facilitate translational genomics for crop improvement and the development of genetically tailored, climate-resilient, high-yielding chickpea cultivars.

Genome wide association analysis for grain micronutrients and anti-nutritional traits in mungbean [Vigna radiata (L.) R. Wilczek] using SNP markers

Mungbean is an important food grain legume for human nutrition and nutritional food due to its nutrient-dense seed, liked palatability, and high digestibility. However, anti-nutritional factors pose a significant risk to improving nutritional quality for bio-fortification. In the present study, genetic architecture of grain micronutrients (grain iron and zinc concentration) and anti-nutritional factors (grain phytic acid and tannin content) in association mapping panel of 145 diverse mungbean were evaluated. Based on all four parameters genotypes PUSA 1333 and IPM 02-19 were observed as desired genotypes as they had high grain iron and zinc concentration but low grain phytic acid and tannin content. The next generation sequencing (NGS)-based genotyping by sequencing (GBS) identified 14,447 genome-wide SNPs in a diverse selected panel of 127 mungbean genotypes. Population admixture analysis revealed the presence of four different ancestries among the genotypes and LD decay of ∼57.6 kb kb physical distance was noted in mungbean chromosomes. Association mapping analysis revealed that a total of 20 significant SNPs were shared by both GLM and Blink models associated with grain micronutrient and anti-nutritional factor traits, with Blink model identifying 35 putative SNPs. Further, this study identified the 185 putative candidate genes. Including potential candidate genes Vradi07g30190, Vradi01g09630, and Vradi09g05450 were found to be associated with grain iron concentration, Vradi10g04830 with grain zinc concentration, Vradi08g09870 and Vradi01g11110 with grain phytic acid content and Vradi04g11580 and Vradi06g15090 with grain tannin content. Moreover, two genes Vradi07g15310 and Vradi09g05480 showed significant variation in protein structure between native and mutated versions. The identified SNPs and candidate genes are potential powerful tools to provide the essential information for genetic studies and marker-assisted breeding program for nutritional improvement in mungbean.

Genome-wide analysis suggests the potential role of lncRNAs during seed development and seed size/weight determination in chickpea

The integrated transcriptome data analyses suggested the plausible roles of lncRNAs during seed development in chickpea. The candidate lncRNAs associated with QTLs and those involved in miRNA-mediated seed size/weight determination in chickpea have been identified. Long non-coding RNAs (lncRNAs) are important regulators of various biological processes. Here, we identified lncRNAs at seven successive stages of seed development in small-seeded and large-seeded chickpea cultivars. In total, 4751 lncRNAs implicated in diverse biological processes were identified. Most of lncRNAs were conserved between the two cultivars, whereas only a few of them were conserved in other plants, suggesting their species-specificity. A large number of lncRNAs differentially expressed between the two chickpea cultivars associated with seed development-related processes were identified. The lncRNAs acting as precursors of miRNAs and those mimicking target protein-coding genes of miRNAs involved in seed size/weight determination, including HAIKU1, BIG SEEDS1, and SHB1, were also revealed. Further, lncRNAs located within seed size/weight associated quantitative trait loci were also detected. Overall, we present a comprehensive resource and identified candidate lncRNAs that may play important roles during seed development and seed size/weight determination in chickpea.

Identification of QTNs Associated With Flowering Time, Maturity, and Plant Height Traits in Linum usitatissimum L. Using Genome-Wide Association Study

Early flowering, maturity, and plant height are important traits for linseed to fit in rice fallows, for rainfed agriculture, and for economically viable cultivation. Here, Multi-Locus Genome-Wide Association Study (ML-GWAS) was undertaken in an association mapping panel of 131 accessions, genotyped using 68,925 SNPs identified by genotyping by sequencing approach. Phenotypic evaluation data of five environments comprising 3 years and two locations were used. GWAS was performed for three flowering time traits including days to 5%, 50%, and 95% flowering, days to maturity, and plant height by employing five ML-GWAS methods: FASTmrEMMA, FASTmrMLM, ISIS EM-BLASSO, mrMLM, and pLARmEB. A total of 335 unique QTNs have been identified for five traits across five environments. 109 QTNs were stable as observed in ≥2 methods and/or environments, explaining up to 36.6% phenotypic variance. For three flowering time traits, days to maturity, and plant height, 53, 30, and 27 stable QTNs, respectively, were identified. Candidate genes having roles in flower, pollen, embryo, seed and fruit development, and xylem/phloem histogenesis have been identified. Gene expression of candidate genes for flowering and plant height were studied using transcriptome of an early maturing variety Sharda (IC0523807). The present study unravels QTNs/candidate genes underlying complex flowering, days to maturity, and plant height traits in linseed.

Genetic variability and population structure of Ethiopian chickpea (Cicer arietinum L.) germplasm

Evaluation of the genetic diversity and an understanding of the genetic structure and relationships of chickpea genotypes are valuable to design efficient germplasm conservation strategies and crop breeding programs. Information is limited, in these regards, for Ethiopian chickpea germplasms. Therefore, the present study was carried out to estimate the genetic diversity, population structure, and relationships of 152 chickpea genotypes using simple sequence repeats (SSR) markers. Twenty three SSR markers exhibited polymorphism producing a total of 133 alleles, with a mean of 5.8 alleles per locus. Analyses utilizing various genetic-based statistics included pairwise population Nei’s genetic distance, heterozygosity, Shannon’s information index, polymorphic information content, and percent polymorphism. These analyses exemplified the existence of high genetic variation within and among chickpea genotypes. The 152 genotypes were divided into two major clusters based on Nei’s genetic distances. The exotic genotypes were grouped in one cluster exclusively showing that these genotypes are distinct to Ethiopian genotypes, while the patterns of clustering of Ethiopian chickpea genotypes based on their geographic region were not consistent because of the seed exchange across regions. Model-based population structure clustering identified two discrete populations. These finding provides useful insight for chickpea collections and ex-situ conservation and national breeding programs for widening the genetic base of chickpea.

Identifying Anti-Oxidant Biosynthesis Genes in Pearl Millet [Pennisetum glaucum (L.) R. Br.] Using Genome-Wide Association Analysis

Pearl millet [Pennisetum glaucum (L.) R Br.] is an important staple food crop in the semi-arid tropics of Asia and Africa. It is a cereal grain that has the prospect to be used as a substitute for wheat flour for celiac patients. It is an important antioxidant food resource present with a wide range of phenolic compounds that are good sources of natural antioxidants. The present study aimed to identify the total antioxidant content of pearl millet flour and apply it to evaluate the antioxidant activity of its 222 genotypes drawn randomly from the pearl millet inbred germplasm association panel (PMiGAP), a world diversity panel of this crop. The total phenolic content (TPC) significantly correlated with DPPH (1,1-diphenyl-2-picrylhydrazyl) radical scavenging activity (% inhibition), which ranged from 2.32 to 112.45% and ferric-reducing antioxidant power (FRAP) activity ranging from 21.68 to 179.66 (mg ascorbic acid eq./100 g). Genome-wide association studies (GWAS) were conducted using 222 diverse accessions and 67 K SNPs distributed across all the seven pearl millet chromosomes. Approximately, 218 SNPs were found to be strongly associated with DPPH and FRAP activity at high confidence [-log (p) > 3.0-7.4]. Furthermore, flanking regions of significantly associated SNPs were explored for candidate gene harvesting. This identified 18 candidate genes related to antioxidant pathway genes (flavanone 7-O-beta-glycosyltransferase, GDSL esterase/lipase, glutathione S-transferase) residing within or near the association signal that can be selected for further functional characterization. Patterns of genetic variability and the associated genes reported in this study are useful findings, which would need further validation before their utilization in molecular breeding for high antioxidant-containing pearl millet cultivars.

Evolutionary Genomics of Plant Gametophytic Selection

It has long been recognized that natural selection during the haploid gametophytic phase of the plant life cycle may have widespread importance for rates of evolution and the maintenance of genetic variation. Recent theoretical advances have further highlighted the significance of gametophytic selection for diverse evolutionary processes. Genomic approaches offer exciting opportunities to address key questions about the extent and effects of gametophytic selection on plant evolution and adaptation. Here, we review the progress and prospects for integrating functional and evolutionary genomics to test theoretical predictions, and to examine the importance of gametophytic selection on genetic diversity and rates of evolution. There is growing evidence that selection during the gametophyte phase of the plant life cycle has important effects on both gene and genome evolution and is likely to have important pleiotropic effects on the sporophyte. We discuss the opportunities to integrate comparative population genomics, genome-wide association studies, and experimental approaches to further distinguish how differential selection in the two phases of the plant life cycle contributes to genetic diversity and adaptive evolution.

Genome-wide association study for phosphate deficiency responsive root hair elongation in chickpea

Root hairs (RHs) are single-celled elongated epidermal cells and play a vital role in nutrient absorption, particularly for immobile minerals like phosphorus (P). As an adaptive response to P deficiency, an increase in RH length enhances root-soil contact and absorptive area for P absorption. Genetic variations have been reported for RH length and its response to P deficiency in plants. However, only a few association studies have been conducted to identify genes and genetic loci associated with RH length. Here, we screened desi chickpea accessions for RH length and its plasticity under P deficiency. Further, the genome-wide association study (GWAS) was conducted to identify the genetic loci associated with RH length in P deficient and sufficient conditions. Although high variability was observed in terms of RH length in diverse genotypes, majority of the accessions showed typical response of increase in RH length in low P. Genome-wide association mapping identified many SNPs with significant associations with RH length in P-sufficient and P-deficient conditions. A few candidate genes for RH length in P deficient (SIZ1-like and HAD superfamily protein) and sufficient (RSL2-like and SMAP1-like) conditions were identified which have known roles in RH development and P deficiency response or both. Highly associated loci and candidate genes identified in this study would be useful for genomic-assisted breeding to develop P-efficient chickpea.

Integrating genomics for chickpea improvement: achievements and opportunities

Integration of genomic technologies with breeding efforts have been used in recent years for chickpea improvement. Modern breeding along with low cost genotyping platforms have potential to further accelerate chickpea improvement efforts. The implementation of novel breeding technologies is expected to contribute substantial improvements in crop productivity. While conventional breeding methods have led to development of more than 200 improved chickpea varieties in the past, still there is ample scope to increase productivity. It is predicted that integration of modern genomic resources with conventional breeding efforts will help in the delivery of climate-resilient chickpea varieties in comparatively less time. Recent advances in genomics tools and technologies have facilitated the generation of large-scale sequencing and genotyping data sets in chickpea. Combined analysis of high-resolution phenotypic and genetic data is paving the way for identifying genes and biological pathways associated with breeding-related traits. Genomics technologies have been used to develop diagnostic markers for use in marker-assisted backcrossing programmes, which have yielded several molecular breeding products in chickpea. We anticipate that a sequence-based holistic breeding approach, including the integration of functional omics, parental selection, forward breeding and genome-wide selection, will bring a paradigm shift in development of superior chickpea varieties. There is a need to integrate the knowledge generated by modern genomics technologies with molecular breeding efforts to bridge the genome-to-phenome gap. Here, we review recent advances that have led to new possibilities for developing and screening breeding populations, and provide strategies for enhancing the selection efficiency and accelerating the rate of genetic gain in chickpea.

Method for Genome-Wide Association Study: A Soybean Example

Genome-wide association studies (GWAS) are a valuable approach to identify single nucleotide polymorphisms (SNPs) associated with a phenotype of interest. There are now a variety of R-packages and command line tools available to perform GWAS. Here, we provide an example downloading and filtering SNP data, followed by GWAS analysis using the R-package rMVP.

Updates on Genomic Resources in Chickpea for Crop Improvement

In recent years, rapid advancement has been done in generation of genomic resources for the important legume crop chickpea. Here, we provide an update on important advancements made on availability of genomic resources for this crop. The availability of reference genome and transcriptome sequences, and resequencing of several accessions have enabled the discovery of gene space and molecular markers in chickpea. These resources have helped in elucidating evolutionary relationship and identification of quantitative trait loci for important agronomic traits. Gene expression in different tissues/organs during development and under abiotic/biotic stresses has been interrogated. In addition, single-base resolution DNA methylation patterns in different organs have been analyzed to understand gene regulation. Overall, we provide a consolidated overview of available genomic resources of chickpea that may help in fulfilling the promises for improvement of this important crop.

Genotyping for Species Identification and Diversity Assessment Using Double-Digest Restriction Site-Associated DNA Sequencing (ddRAD-Seq)

Genotyping-by-sequencing (GBS) is a powerful approach for studying the genetic diversity of legume species. By using restriction enzymes or other methods to generate a reduced representation of the genome for sequencing, GBS can provide genome-wide single nucleotide polymorphisms (SNP) for diversity analysis at high throughput and low cost. Here we describe a novel double-digest restriction site-associated DNA sequencing (ddRAD-seq) approach. We also describe the downstream bioinformatic analysis of the sequencing data, including alignment to a reference genome, de novo assembly, SNP calling, phylogenetic analysis, and structure analysis.

The Mediator subunit OsMED15a is a transcriptional co-regulator of seed size/weight-modulating genes in rice

Although several transcription factors (TFs) that regulate seed size/weight in plants are known, the molecular landscape regulating this important trait is unclear. Here, we report that a Mediator subunit, OsMED15a, links rice grain size/weight-regulating TFs to their target genes. Expression analysis and high-resolution quantitative trait loci (QTL) mapping suggested that OsMED15a is involved in rice seed development. OsMED15a has an N-terminal, three-helical KIX domain. Two of these helices, α1 and α3, and three amino acids, 76LRC78, within OsMED15a helix α3 were important for its interaction with several proteins, including interactions with the transactivation domains of two NAC-type TFs, OsNAC024 and OsNAC025. Moreover, OsMED15a, OsNAC024, and OsNAC025 all exhibited increased expression during seed development, and we identified several grain size/weight-associated SNPs in these genes in 509 low- and high-grain-weight rice genotypes. RNAi-mediated repression of OsMED15a expression down-regulated the expression of the grain size/weight regulating genes GW2, GW5 and DR11 and reduced grain length, weight, and yield. Of note, both OsNAC024 and OsNAC025 bound to the promoters of these three genes. We conclude that the transactivation domains of OsNAC024 and OsNAC025 target the KIX domain of OsMED15a in the regulation of grain size/weight-associated genes such as GW2, GW5, and D11. We propose that the integrated molecular-genetics approach used here could help identify networks of functional alleles of other regulator and co-regulator genes and thereby inform efforts for marker-assisted introgression of useful alleles in rice crop improvement.

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.

Towards Exploitation of Adaptive Traits for Climate-Resilient Smart Pulses

Pulses are the main source of protein and minerals in the vegetarian diet. These are primarily cultivated on marginal lands with few inputs in several resource-poor countries of the world, including several in South Asia. Their cultivation in resource-scarce conditions exposes them to various abiotic and biotic stresses, leading to significant yield losses. Furthermore, climate change due to global warming has increased their vulnerability to emerging new insect pests and abiotic stresses that can become even more serious in the coming years. The changing climate scenario has made it more challenging to breed and develop climate-resilient smart pulses. Although pulses are climate smart, as they simultaneously adapt to and mitigate the effects of climate change, their narrow genetic diversity has always been a major constraint to their improvement for adaptability. However, existing genetic diversity still provides opportunities to exploit novel attributes for developing climate-resilient cultivars. The mining and exploitation of adaptive traits imparting tolerance/resistance to climate-smart pulses can be accelerated further by using cutting-edge approaches of biotechnology such as transgenics, genome editing, and epigenetics. This review discusses various classical and molecular approaches and strategies to exploit adaptive traits for breeding climate-smart pulses.

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.

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.

Toward the sequence-based breeding in legumes in the post-genome sequencing era

Efficiency of breeding programs of legume crops such as chickpea, pigeonpea and groundnut has been considerably improved over the past decade through deployment of modern genomic tools and technologies. For instance, next-generation sequencing technologies have facilitated availability of genome sequence assemblies, re-sequencing of several hundred lines, development of HapMaps, high-density genetic maps, a range of marker genotyping platforms and identification of markers associated with a number of agronomic traits in these legume crops. Although marker-assisted backcrossing and marker-assisted selection approaches have been used to develop superior lines in several cases, it is the need of the hour for continuous population improvement after every breeding cycle to accelerate genetic gain in the breeding programs. In this context, we propose a sequence-based breeding approach which includes use of independent or combination of parental selection, enhancing genetic diversity of breeding programs, forward breeding for early generation selection, and genomic selection using sequencing/genotyping technologies. Also, adoption of speed breeding technology by generating 4-6 generations per year will be contributing to accelerate genetic gain. While we see a huge potential of the sequence-based breeding to revolutionize crop improvement programs in these legumes, we anticipate several challenges especially associated with high-quality and precise phenotyping at affordable costs, data analysis and management related to improving breeding operation efficiency. Finally, integration of improved seed systems and better agronomic packages with the development of improved varieties by using sequence-based breeding will ensure higher genetic gains in farmers' fields.

Next-Generation Sequencing Accelerates Crop Gene Discovery

The identification and isolation of genes underlying quantitative trait loci (QTLs) associated with agronomic traits in crops have been recently accelerated thanks to next-generation sequencing (NGS)-based technologies combined with plant genetics. With NGS, various revisited genetic approaches, which benefited from higher marker density, have been elaborated. These approaches improved resolution in QTL position and assisted in determining functional causative variations in genes. Examples of QTLs/genes associated with agronomic traits in crops and identified using different strategies based on whole-genome sequencing (WGS)/whole-genome resequencing (WGR) or RNA-seq are presented and discussed in this review. More specifically, we summarize and illustrate how NGS boosted bulk-segregant analysis (BSA), expression profiling, and the construction of polymorphism databases to facilitate the detection of QTLs and causative genes.

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.

Secondary Wall Regulating NACs Differentially Bind at the Promoter at a CELLULOSE SYNTHASE A4 Cis-eQTL

Arabidopsis thaliana CELLULOSE SYNTHASE A4/7/8 (CESA4/7/8) are three non-redundant subunits of the secondary cell wall cellulose synthase complex. Transcript abundance of these genes can vary among genotypes and expression quantitative trait loci (eQTL) were identified in a recombinant population of the accessions Bay-0 and Shahdara. Genetic mapping and analysis of the transcript levels of CESAs between two distinct near isogenic lines (NILs) confirmed a change in CESA4 expression that segregates within that interval. We sequenced the promoters and identified 16 polymorphisms differentiating CESA4^Sha and CESA4^Bay . In order to determine which of these SNPs could be responsible for this eQTL, we screened for transcription factor protein affinity with promoter fragments of CESA4^Bay, CESA4^Sha , and the reference genome CESA4^Col . The wall thickening activator proteins NAC SECONDARY WALL THICKENING PROMOTING FACTOR2 (NST2) and NST3 exhibited a decrease in binding with the CESA4^Sha promoter with a tracheary element-regulating cis-element (TERE) polymorphism. While NILs harboring the TERE polymorphisms exhibited significantly different CESA4 expression, cellulose crystallinity and cell wall thickness were indistinguishable. These results suggest that the TERE polymorphism resulted in differential transcription factor binding and CESA4 expression; yet A. thaliana is able to tolerate this transcriptional variability without compromising the structural elements of the plant, providing insight into the elasticity of gene regulation as it pertains to cell wall biosynthesis and regulation. We also explored available DNA affinity purification sequencing data to resolve a core binding site, C(G/T)TNNNNNNNA(A/C)G, for secondary wall NACs referred to as the VNS element.

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.

Current advances in chickpea genomics: applications and future perspectives

Chickpea genomics promises to illuminate our understanding of genome organization, structural variations, evolutionary and domestication-related insights and fundamental biology of legume crops. Unprecedented advancements of next generation sequencing (NGS) technologies have enabled in decoding of multiple chickpea genome sequences and generating huge genomic resources in chickpea both at functional and structural level. This review is aimed to update the current progress of chickpea genomics ranging from high density linkage map development, genome-wide association studies (GWAS), functional genomics resources for various traits, emerging role of abiotic stress responsive coding and non-coding RNAs after the completion of draft chickpea genome sequences. Additionally, the current efforts of whole genome re-sequencing (WGRS) approach of global chickpea germplasm to capture the global genetic diversity existing in the historically released varieties across the world and increasing the resolution of the previously identified candidate gene(s) of breeding importance have been discussed. Thus, the outcomes of these genomics resources will assist in genomics-assisted selection and facilitate breeding of climate-resilient chickpea cultivars for sustainable agriculture.

Association mapping unveils favorable alleles for grain iron and zinc concentrations in lentil (Lens culinaris subsp. culinaris)

Lentil is a major cool-season grain legume grown in South Asia, West Asia, and North Africa. Populations in developing countries of these regions have micronutrient deficiencies; therefore, breeding programs should focus more on improving the micronutrient content of food. In the present study, a set of 96 diverse germplasm lines were evaluated at three different locations in India to examine the variation in iron (Fe) and zinc (Zn) concentration and identify simple sequence repeat (SSR) markers that associate with the genetic variation. The genetic variation among genotypes of the association mapping (AM) panel was characterized using a genetic distance-based and a general model-based clustering method. The model-based analysis identified six subpopulations, which satisfactorily explained the genetic structure of the AM panel. AM analysis identified three SSRs (PBALC 13, PBALC 206, and GLLC 563) associated with grain Fe concentration explaining 9% to 11% of phenotypic variation and four SSRs (PBALC 353, SSR 317-1, PLC 62, and PBALC 217) were associated with grain Zn concentration explaining 14%, to 21% of phenotypic variation. These identified SSRs exhibited consistent performance across locations. These candidate SSRs can be used in marker-assisted genetic improvement for developing Fe and Zn fortified lentil varieties. Favorable alleles and promising genotypes identified in this study can be utilized for lentil biofortification.

Identification of Purple Acid Phosphatases in Chickpea and Potential Roles of CaPAP7 in Seed Phytate Accumulation

Purple acid phosphatases (PAPs) play important roles in phosphate (Pi) acquisition and utilization. These PAPs hydrolyze organic Phosphorus (P) containing compounds in rhizosphere as well as inside the plant cell. However, roles of PAPs in one of the most widely cultivated legumes, chickpea (Cicer arietnum L.), have not been unraveled so far. In the present study, we identified 25 putative PAPs in chickpea (CaPAPs) which possess functional PAP motifs and domains. Differential regulation of CaPAPs under different nutrient deficiencies revealed their roles under multiple nutrient stresses including Pi deficiency. Interestingly, most of the CaPAPs were prominently expressed in flowers and young pods indicating their roles in flower and seed development. Association mapping of SNPs underlying CaPAPs with seed traits revealed significant association of low Pi inducible CaPAP7 with seed weight and phytate content. Biochemical characterization of recombinant CaPAP7 established it to be a functional acid phosphatase with highest activity on most abundant organic-P substrate, phytate. Exogenous application of recombinant CaPAP7 enhanced biomass and Pi content of Arabidopsis seedlings supplemented with phytate as sole P source. Taken together, our results uncover the PAPs in chickpea and potential roles of CaPAP7 in seed phytate accumulation.

Assessing and Exploiting Functional Diversity in Germplasm Pools to Enhance Abiotic Stress Adaptation and Yield in Cereals and Food Legumes

There is a need to accelerate crop improvement by introducing alleles conferring host plant resistance, abiotic stress adaptation, and high yield potential. Elite cultivars, landraces and wild relatives harbor useful genetic variation that needs to be more easily utilized in plant breeding. We review genome-wide approaches for assessing and identifying alleles associated with desirable agronomic traits in diverse germplasm pools of cereals and legumes. Major quantitative trait loci and single nucleotide polymorphisms (SNPs) associated with desirable agronomic traits have been deployed to enhance crop productivity and resilience. These include alleles associated with variation conferring enhanced photoperiod and flowering traits. Genetic variants in the florigen pathway can provide both environmental flexibility and improved yields. SNPs associated with length of growing season and tolerance to abiotic stresses (precipitation, high temperature) are valuable resources for accelerating breeding for drought-prone environments. Both genomic selection and genome editing can also harness allelic diversity and increase productivity by improving multiple traits, including phenology, plant architecture, yield potential and adaptation to abiotic stresses. Discovering rare alleles and useful haplotypes also provides opportunities to enhance abiotic stress adaptation, while epigenetic variation has potential to enhance abiotic stress adaptation and productivity in crops. By reviewing current knowledge on specific traits and their genetic basis, we highlight recent developments in the understanding of crop functional diversity and identify potential candidate genes for future use. The storage and integration of genetic, genomic and phenotypic information will play an important role in ensuring broad and rapid application of novel genetic discoveries by the plant breeding community. Exploiting alleles for yield-related traits would allow improvement of selection efficiency and overall genetic gain of multigenic traits. An integrated approach involving multiple stakeholders specializing in management and utilization of genetic resources, crop breeding, molecular biology and genomics, agronomy, stress tolerance, and reproductive/seed biology will help to address the global challenge of ensuring food security in the face of growing resource demands and climate change induced stresses.

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.

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.

Diversity and association of phenotypic and metabolomic traits in the close model grasses Brachypodium distachyon, B. stacei and B. hybridum

BACKGROUND AND AIMS

Morphological traits in combination with metabolite fingerprinting were used to investigate inter- and intraspecies diversity within the model annual grasses Brachypodium distachyon, Brachypodium stacei and Brachypodium hybridum .

METHODS

Phenotypic variation of 15 morphological characters and 2219 nominal mass ( m / z ) signals generated using flow infusion electrospray ionization-mass spectrometry (FIE-MS) were evaluated in individuals from a total of 174 wild populations and six inbred lines, and 12 lines, of the three species, respectively. Basic statistics and multivariate principal component analysis and discriminant analysis were used to differentiate inter- and intraspecific variability of the two types of variable, and their association was assayed with the rcorr function.

KEY RESULTS

Basic statistics and analysis of variance detected eight phenotypic characters [(stomata) leaf guard cell length, pollen grain length, (plant) height, second leaf width, inflorescence length, number of spikelets per inflorescence, lemma length, awn length] and 434 tentatively annotated metabolite signals that significantly discriminated the three species. Three phenotypic traits (pollen grain length, spikelet length, number of flowers per inflorescence) might be genetically fixed. The three species showed different metabolomic profiles. Discriminant analysis significantly discriminated the three taxa with both morphometric and metabolome traits and the intraspecific phenotypic diversity within B. distachyon and B. stacei . The populations of B. hybridum were considerably less differentiated.

CONCLUSIONS

Highly explanatory metabolite signals together with morphological characters revealed concordant patterns of differentiation of the three taxa. Intraspecific phenotypic diversity was observed between northern and southern Iberian populations of B. distachyon and between eastern Mediterranean/south-western Asian and western Mediterranean populations of B. stacei . Significant association was found for pollen grain length and lemma length and ten and six metabolomic signals, respectively. These results would guide the selection of new germplasm lines of the three model grasses in ongoing genome-wide association studies.

Genome resources for climate-resilient cowpea, an essential crop for food security

Cowpea (Vigna unguiculata L. Walp.) is a legume crop that is resilient to hot and drought-prone climates, and a primary source of protein in sub-Saharan Africa and other parts of the developing world. However, genome resources for cowpea have lagged behind most other major crops. Here we describe foundational genome resources and their application to the analysis of germplasm currently in use in West African breeding programs. Resources developed from the African cultivar IT97K-499-35 include a whole-genome shotgun (WGS) assembly, a bacterial artificial chromosome (BAC) physical map, and assembled sequences from 4355 BACs. These resources and WGS sequences of an additional 36 diverse cowpea accessions supported the development of a genotyping assay for 51 128 SNPs, which was then applied to five bi-parental RIL populations to produce a consensus genetic map containing 37 372 SNPs. This genetic map enabled the anchoring of 100 Mb of WGS and 420 Mb of BAC sequences, an exploration of genetic diversity along each linkage group, and clarification of macrosynteny between cowpea and common bean. The SNP assay enabled a diversity analysis of materials from West African breeding programs. Two major subpopulations exist within those materials, one of which has significant parentage from South and East Africa and more diversity. There are genomic regions of high differentiation between subpopulations, one of which coincides with a cluster of nodulin genes. The new resources and knowledge help to define goals and accelerate the breeding of improved varieties to address food security issues related to limited-input small-holder farming and climate stress.

Functional natural allelic variants of flavonoid 3',5'-hydroxylase gene governing catechin traits in tea plant and its relatives

Functional allelic variants of the flavonoid 3',5'-hydroxylase (F3'5'H) gene provides new information of F3'5'H function of tea plant and its relatives. This insight may serve as the foundation upon which to advance molecular breeding in the tea plant. Catechins are the active components of tea that determine its quality and health attributes. This study established the first integrated genomic strategy for deciphering the genetic basis of catechin traits of tea plant. With the RNA-sequencing analysis of bulked segregants representing the tails of a F1 population segregated for total catechin content, we identified a flavonoid 3',5'-hydroxylase (F3'5'H) gene. F3'5'H had one copy in the genomic DNA of tea plant. Among 202 tea accessions, we identified 120 single nucleotide polymorphisms (SNPs) at F3'5'H locus. Seventeen significant marker-trait associations were identified by association mapping in multiple environments, which were involved in 10 SNP markers, and the traits including the ratio of di/tri-hydroxylated catechins and catechin contents. The associated individual and combination of SNPs explained 4.5-25.2 and 53.0-63.0% phenotypic variations, respectively. In the F1 population (validation population), the catechin trait variation percentages explained by F3'5'H diplotype were 6.9-74.3%. The genotype effects of ten functional SNPs in the F1 population were all consistent with the association population. Furthermore, the function of SNP-711/-655 within F3'5'H was validated by gene expression analysis. Altogether, our work indicated functional SNP allelic variants within F3'5'H governing the ratio of di/tri-hydroxylated catechins and catechin contents. The strong catechin-associated SNPs identified in this study can be used for future marker-assisted selection to improve tea quality.

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.

Genetic diversity analysis of Gossypium arboreum germplasm accessions using genotyping-by-sequencing

The diploid cotton species Gossypium arboreum possesses many favorable agronomic traits such as drought tolerance and disease resistance, which can be utilized in the development of improved upland cotton cultivars. The USDA National Plant Germplasm System maintains more than 1600 G. arboreum accessions. Little information is available on the genetic diversity of the collection thereby limiting the utilization of this cotton species. The genetic diversity and population structure of the G. arboreum germplasm collection were assessed by genotyping-by-sequencing of 375 accessions. Using genome-wide single nucleotide polymorphism sequence data, two major clusters were inferred with 302 accessions in Cluster 1, 64 accessions in Cluster 2, and nine accessions unassigned due to their nearly equal membership to each cluster. These two clusters were further evaluated independently resulting in the identification of two sub-clusters for the 302 Cluster 1 accessions and three sub-clusters for the 64 Cluster 2 accessions. Low to moderate genetic diversity between clusters and sub-clusters were observed indicating a narrow genetic base. Cluster 2 accessions were more genetically diverse and the majority of the accessions in this cluster were landraces. In contrast, Cluster 1 is composed of varieties or breeding lines more recently added to the collection. The majority of the accessions had kinship values ranging from 0.6 to 0.8. Eight pairs of accessions were identified as potential redundancies due to their high kinship relatedness. The genetic diversity and genotype data from this study are essential to enhance germplasm utilization to identify genetically diverse accessions for the detection of quantitative trait loci associated with important traits that would benefit upland cotton improvement.

Single nucleotide polymorphism in sugar pathway and disease resistance genes in sugarcane

Single nucleotide polymorphism in sugar pathway and disease resistance genes showing genetic association with sugar content and red rot resistance would be useful in marker-assisted genetic improvement of sugarcane. Validation and genotyping of potential sequence variants in candidate genes are necessary to understand their functional significance and trait association potential. We discovered, characterized, validated and genotyped SNPs and InDels in sugar pathway and disease resistance genes of Saccharum complex and sugarcane varieties using amplicon sequencing and CAPS assays. The SNPs were abundant in the non-coding 3'UTRs than 5'UTRs and coding sequences depicting a strong bias toward C to T transition substitutions than transversions. Sequencing of cloned amplicons validated 61.6 and 45.2 % SNPs detected in silico in 21 sugar pathway and 16 disease resistance genes, respectively. Sixteen SNPs in four sugar pathway genes and 10 SNPs in nine disease resistance genes were validated through cost-effective CAPS assay. Functional and adaptive significance of SNP and protein haplotypes identified in sugar pathway and disease resistance genes was assessed by correlating their allelic variation with missense amino acid substitutions in the functional domains, alteration in protein structure models and possible modulation of catalytic enzyme activity in contrasting high and low sugar and moderately red rot resistant and highly susceptible sugarcane genotypes. A strong genetic association of five SNPs in the sugar pathway and disease resistance genes, and an InDel marker in the promoter sequence of sucrose synthase-2 gene, with sugar content and red rot resistance, was evident. The functionally relevant SNPs and InDels, detected and validated in sugar pathway and disease resistance genes, and genic CAPS markers designed, would be of immense use in marker-assisted genetic improvement of sugarcane for sugar content and disease resistance.

Identification of candidate genes and natural allelic variants for QTLs governing plant height in chickpea

In the present study, molecular mapping of high-resolution plant height QTLs was performed by integrating 3625 desi genome-derived GBS (genotyping-by-sequencing)-SNPs on an ultra-high resolution intra-specific chickpea genetic linkage map (dwarf/semi-dwarf desi cv. ICC12299 x tall kabuli cv. ICC8261). The identified six major genomic regions harboring six robust QTLs (11.5-21.3 PVE), associated with plant height, were mapped within <0.5 cM average marker intervals on six chromosomes. Five SNPs-containing genes tightly linked to the five plant height QTLs, were validated based upon their high potential for target trait association (12.9-20.8 PVE) in 65 desi and kabuli chickpea accessions. The vegetative tissue-specific expression, including higher differential up-regulation (>5-fold) of five genes especially in shoot, young leaf, shoot apical meristem of tall mapping parental accession (ICC8261) as compared to that of dwarf/semi-dwarf parent (ICC12299) was apparent. Overall, combining high-resolution QTL mapping with genetic association analysis and differential expression profiling, delineated natural allelic variants in five candidate genes (encoding cytochrome-c-biosynthesis protein, malic oxidoreductase, NADH dehydrogenase iron-sulfur protein, expressed protein and bZIP transcription factor) regulating plant height in chickpea. These molecular tags have potential to dissect complex plant height trait and accelerate marker-assisted genetic enhancement for developing cultivars with desirable plant height ideotypes in chickpea.

Emerging Genomic Tools for Legume Breeding: Current Status and Future Prospects

Legumes play a vital role in ensuring global nutritional food security and improving soil quality through nitrogen fixation. Accelerated higher genetic gains is required to meet the demand of ever increasing global population. In recent years, speedy developments have been witnessed in legume genomics due to advancements in next-generation sequencing (NGS) and high-throughput genotyping technologies. Reference genome sequences for many legume crops have been reported in the last 5 years. The availability of the draft genome sequences and re-sequencing of elite genotypes for several important legume crops have made it possible to identify structural variations at large scale. Availability of large-scale genomic resources and low-cost and high-throughput genotyping technologies are enhancing the efficiency and resolution of genetic mapping and marker-trait association studies. Most importantly, deployment of molecular breeding approaches has resulted in development of improved lines in some legume crops such as chickpea and groundnut. In order to support genomics-driven crop improvement at a fast pace, the deployment of breeder-friendly genomics and decision support tools seems appear to be critical in breeding programs in developing countries. This review provides an overview of emerging genomics and informatics tools/approaches that will be the key driving force for accelerating genomics-assisted breeding and ultimately ensuring nutritional and food security in developing countries.

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.