A draft genome sequence of the pulse crop chickpea (Cicer arietinum L.)

Physical Seed Dormancy in Legumes: Molecular Advances and Perspectives

Physical dormancy of seeds is a form of dormancy due to the presence of an impermeable seed coat layer, and it represents a feature for plants to adapt to environmental changes over an extended period of phylogenetic evolution. However, in agricultural practice, physical dormancy is problematic. because it prevents timely and uniform seed germination. Therefore, physical dormancy is an important agronomical trait to target in breeding and domestication, especially for many leguminous crops. Compared to the well-characterized physiological dormancy, research progress on physical dormancy at the molecular level has been limited until recent years, due to the lack of suitable research materials. This review focuses on the structure of seed coat, factors affecting physical dormancy, genes controlling physical dormancy, and plants suitable for studying physical dormancy at the molecular level. Our goal is to provide a plethora of information for further molecular research on physical dormancy.

Convergent and/or parallel evolution of RNA-binding proteins in angiosperms after polyploidization

Increasing studies suggest that the biased retention of stress-related transcription factors (TFs) after whole-genome duplications (WGDs) could rewire gene transcriptional networks, facilitating plant adaptation to challenging environments. However, the role of posttranscriptional factors (e.g. RNA-binding proteins, RBPs) following WGDs has been largely ignored. Uncovering thousands of RBPs in 21 representative angiosperm species, we integrate genomic, transcriptomic, regulatomic, and paleotemperature datasets to unravel their evolutionary trajectories and roles in adapting to challenging environments. We reveal functional enrichments of RBP genes in stress responses and identify their convergent retention across diverse angiosperms from independent WGDs, coinciding with global cooling periods. Numerous RBP duplicates derived from WGDs are then identified as cold-induced. A significant overlap of 29 orthogroups between WGD-derived and cold-induced RBP genes across diverse angiosperms highlights a correlation between WGD and cold stress. Notably, we unveil an orthogroup (Glycine-rich RNA-binding Proteins 7/8, GRP7/8) and relevant TF duplicates (CCA1/LHY, RVE4/8, CBF2/4, etc.), co-retained in different angiosperms post-WGDs. Finally, we illustrate their roles in rewiring circadian and cold-regulatory networks at both transcriptional and posttranscriptional levels during global cooling. Altogether, we underline the adaptive evolution of RBPs in angiosperms after WGDs during global cooling, improving our understanding of plants surviving periods of environmental turmoil.

Efficient Single Nucleotide Polymorphism Marker-Assisted Selection to Fusarium Wilt in Chickpea

Fusarium wilt is one of the most destructive chickpea diseases worldwide. Race 5 (Foc5) is the most harmful in the Mediterranean basin. The primary objective of this study is to validate a block of six SNP markers previously mapped in Ca2 in a diverse panel of cultivars, advanced and inbred lines phenotyped for resistance to fusarium wilt. Additionally, we aim to assess the effectiveness of using these markers in the selection of resistant Foc5 lines in an ongoing breeding program. The results showed a 100% coincidence between phenotype and expected haplotype in plant material evaluated for Foc5. We also analyzed 67 inbred lines previously phenotyped by different authors for fusarium wilt reaction, though the specific race was not specified. In these accessions, 65.8% of the analyzed lines exhibited complete correspondence between the phenotype and haplotype. Our results suggest that in early generations it is possible to select resistant materials with reliability, leading to the removal of a significant number of lines, thereby reducing costs and facilitating the handling of materials for additional trait evaluations. Functional annotation of genes delimited by the SNP block revealed several genes in the "response to stimulus" category with potential roles in the resistance reaction.

A chromosome-level genome assembly for Onobrychis viciifolia reveals gene copy number gain underlying enhanced proanthocyanidin biosynthesis

Sainfoin (Onobrychis viciifolia), which belongs to subfamily Papilionoideae of Leguminosae, is a vital perennial forage known as "holy hay" due to its high contents of crude proteins and proanthocyanidins (PAs, also called condensed tannins) that have various pharmacological properties in animal feed, such as alleviating rumen tympanic disease in ruminants. In this study, we select an autotetraploid common sainfoin (2n = 4x = 28) and report its high-quality chromosome-level genome assembly with 28 pseudochromosomes and four haplotypes (~1950.14 Mb, contig N50 = 10.91 Mb). The copy numbers of genes involved in PA biosynthesis in sainfoin are significantly greater than those in four selected Fabales species, namely, autotetraploid Medicago sativa and three other diploid species, Lotus japonicus, Medicago truncatula, and Glycine max. Furthermore, gene expansion is confirmed to be the key contributor to the increased expression of these genes and subsequent PA enhancement in sainfoin. Transcriptomic analyses reveal that the expression of genes involved in the PA biosynthesis pathway is significantly increased in the lines with high PA content compared to the lines with medium and low PA content. The sainfoin genome assembly will improve our understanding of leguminous genome evolution and biosynthesis of secondary metabolites in sainfoin.

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.

MicroRNA397 regulates tolerance to drought and fungal infection by regulating lignin deposition in chickpea root

Plants deposit lignin in the secondary cell wall as a common response to drought and pathogen attacks. Cell wall localised multicopper oxidase family enzymes LACCASES (LACs) catalyse the formation of monolignol radicals and facilitate lignin formation. We show an upregulation of the expression of several LAC genes and a downregulation of microRNA397 (CamiR397) in response to natural drought in chickpea roots. CamiR397 was found to target LAC4 and LAC17L out of twenty annotated LACs in chickpea. CamiR397 and its target genes are expressed in the root. Overexpression of CamiR397 reduced expression of LAC4 and LAC17L and lignin deposition in chickpea root xylem causing reduction in xylem wall thickness. Downregulation of CamiR397 activity by expressing a short tandem target mimic (STTM397) construct increased root lignin deposition in chickpea. CamiR397-overexpressing and STTM397 chickpea lines showed sensitivity and tolerance, respectively, towards natural drought. Infection with a fungal pathogen Macrophomina phaseolina, responsible for dry root rot (DRR) disease in chickpea, induced local lignin deposition and LAC gene expression. CamiR397-overexpressing and STTM397 chickpea lines showed more sensitivity and tolerance, respectively, to DRR. Our results demonstrated the regulatory role of CamiR397 in root lignification during drought and DRR in an agriculturally important crop chickpea.

Genetic mapping of quantitative trait loci associated with drought tolerance in chickpea (Cicer arietinum L.)

Elucidation of the genetic basis of drought tolerance is vital for genomics-assisted breeding of drought tolerant crop varieties. Here, we used genotyping-by-sequencing (GBS) to identify single nucleotide polymorphisms (SNPs) in recombinant inbred lines (RILs) derived from a cross between a drought tolerant chickpea variety, Pusa 362 and a drought sensitive variety, SBD 377. The GBS identified a total of 35,502 SNPs and subsequent filtering of these resulted in 3237 high-quality SNPs included in the eight linkage groups. Fifty-one percent of these SNPs were located in the genic regions distributed throughout the genome. The high density linkage map has total map length of 1069 cm with an average marker interval of 0.33 cm. The linkage map was used to identify 9 robust and consistent QTLs for four drought related traits viz. membrane stability index, relative water content, seed weight and yield under drought, with percent variance explained within the range of 6.29%-90.68% and LOD scores of 2.64 to 6.38, which were located on five of the eight linkage groups. A genomic region on LG 7 harbors quantitative trait loci (QTLs) explaining > 90% phenotypic variance for membrane stability index, and > 10% PVE for yield. This study also provides the first report of major QTLs for physiological traits such as membrane stability index and relative water content for drought stress in chickpea. A total of 369 putative candidate genes were identified in the 6.6 Mb genomic region spanning these QTLs. In-silico expression profiling based on the available transcriptome data revealed that 326 of these genes were differentially expressed under drought stress. KEGG analysis resulted in reduction of candidate genes from 369 to 99, revealing enrichment in various signaling pathways. Haplotype analysis confirmed 5 QTLs among the initially identified 9 QTLs. Two QTLs, qRWC1.1 and qYLD7.1, were chosen based on high SNP density. Candidate gene-based analysis revealed distinct haplotypes in qYLD7.1 associated with significant phenotypic differences, potentially linked to pathways for secondary metabolite biosynthesis. These identified candidate genes bolster defenses through flavonoids and phenylalanine-derived compounds, aiding UV protection, pathogen resistance, and plant structure.The study provides novel genomic regions and candidate genes which can be utilized in genomics-assisted breeding of superior drought tolerant chickpea cultivars.

Association mapping in bambara groundnut [Vigna subterranea (L.) Verdc.] reveals loci associated with agro-morphological traits

BACKGROUND

Genome-wide association studies (GWAS) are important for the acceleration of crop improvement through knowledge of marker-trait association (MTA). This report used DArT SNP markers to successfully perform GWAS on agro-morphological traits using 270 bambara groundnut [Vigna subterranea (L.) Verdc.] landraces sourced from diverse origins. The study aimed to identify marker traits association for nine agronomic traits using GWAS and their candidate genes. The experiment was conducted at two different locations laid out in alpha lattice design. The cowpea [Vigna unguiculata (L.) Walp.] reference genome (i.e. legume genome most closely related to bambara groundnut) assisted in the identification of candidate genes.

RESULTS

The analyses showed that linkage disequilibrium was found to decay rapidly with an average genetic distance of 148 kb. The broadsense heritability was relatively high and ranged from 48.39% (terminal leaf length) to 79.39% (number of pods per plant). The GWAS identified a total of 27 significant marker-trait associations (MTAs) for the nine studied traits explaining 5.27% to 24.86% of phenotypic variations. Among studied traits, the highest number of MTAs was obtained from seed coat colour (6) followed by days to flowering (5), while the least is days to maturity (1), explaining 5.76% to 11.03%, 14.5% to 19.49%, and 11.66% phenotypic variations, respectively. Also, a total of 17 candidate genes were identified, varying in number for different traits; seed coat colour (6), days to flowering (3), terminal leaf length (2), terminal leaf width (2), number of seed per pod (2), pod width (1) and days to maturity (1).

CONCLUSION

These results revealed the prospect of GWAS in identification of SNP variations associated with agronomic traits in bambara groundnut. Also, its present new opportunity to explore GWAS and marker assisted strategies in breeding of bambara groundnut for acceleration of the crop improvement.

Development of an Agrobacterium-delivered codon-optimized CRISPR/Cas9 system for chickpea genome editing

Chickpea is considered recalcitrant to in vitro tissue culture amongst all edible legumes. The clustered, regularly interspaced short palindromic repeats/CRISPR-associated protein 9 (CRISPR/Cas9)-based genome editing in chickpea can remove the bottleneck of limited genetic variation in this cash crop, which is rich in nutrients and protein. However, generating stable mutant lines using CRISPR/Cas9 requires efficient and highly reproducible transformation protocols. As an attempt to solve this problem, we developed a modified and optimized protocol for chickpea transformation. This study transformed the single cotyledon half-embryo explants using CaMV35S promoter to drive two marker genes (β-glucuronidase gene; GUS and green fluorescent protein; GFP) through binary vectors pBI101.2 and modified pGWB2, respectively. These vectors were delivered in the explants through three different strains of Agrobacterium tumefaciens, viz., GV3101, EHA105, and LBA4404. We found better efficiency with the strain GV3101 (17.56%) compared with two other strains, i.e., 8.54 and 5.43%, respectively. We recorded better regeneration frequencies in plant tissue culture for the constructs GUS and GFP, i.e., 20.54% and 18.09%, respectively. The GV3101 was further used for the transformation of the genome editing construct. For the development of genome-edited plants, we used this modified protocol. We also used a modified binary vector pPZP200 by introducing a CaMV35S-driven chickpea codon-optimized SpCas9 gene. The promoter of the Medicago truncatula U6.1 snRNA gene was used to drive the guide RNA cassettes. This cassette targeted and edited the chickpea phytoene desaturase (CaPDS) gene. A single gRNA was found sufficient to achieve high efficiency (42%) editing with the generation of PDS mutants with albino phenotypes. A simple, rapid, highly reproducible, stable transformation and CRISPR/Cas9-based genome editing system for chickpea was established. This study aimed to demonstrate this system's applicability by performing a gene knockout of the chickpea PDS gene using an improved chickpea transformation protocol for the first time.

The Progression in Developing Genomic Resources for Crop Improvement

Sequencing technologies have rapidly evolved over the past two decades, and new technologies are being continually developed and commercialized. The emerging sequencing technologies target generating more data with fewer inputs and at lower costs. This has also translated to an increase in the number and type of corresponding applications in genomics besides enhanced computational capacities (both hardware and software). Alongside the evolving DNA sequencing landscape, bioinformatics research teams have also evolved to accommodate the increasingly demanding techniques used to combine and interpret data, leading to many researchers moving from the lab to the computer. The rich history of DNA sequencing has paved the way for new insights and the development of new analysis methods. Understanding and learning from past technologies can help with the progress of future applications. This review focuses on the evolution of sequencing technologies, their significant enabling role in generating plant genome assemblies and downstream applications, and the parallel development of bioinformatics tools and skills, filling the gap in data analysis techniques.

Comparative transcriptome analysis reveals molecular regulation of salt tolerance in two contrasting chickpea genotypes

Salinity is a major abiotic stress that causes substantial agricultural losses worldwide. Chickpea (Cicer arietinum L.) is an important legume crop but is salt-sensitive. Previous physiological and genetic studies revealed the contrasting response of two desi chickpea varieties, salt-sensitive Rupali and salt-tolerant Genesis836, to salt stress. To understand the complex molecular regulation of salt tolerance mechanisms in these two chickpea genotypes, we examined the leaf transcriptome repertoire of Rupali and Genesis836 in control and salt-stressed conditions. Using linear models, we identified categories of differentially expressed genes (DEGs) describing the genotypic differences: salt-responsive DEGs in Rupali (1,604) and Genesis836 (1,751) with 907 and 1,054 DEGs unique to Rupali and Genesis836, respectively, salt responsive DEGs (3,376), genotype-dependent DEGs (4,170), and genotype-dependent salt-responsive DEGs (122). Functional DEG annotation revealed that the salt treatment affected genes involved in ion transport, osmotic adjustment, photosynthesis, energy generation, stress and hormone signalling, and regulatory pathways. Our results showed that while Genesis836 and Rupali have similar primary salt response mechanisms (common salt-responsive DEGs), their contrasting salt response is attributed to the differential expression of genes primarily involved in ion transport and photosynthesis. Interestingly, variant calling between the two genotypes identified SNPs/InDels in 768 Genesis836 and 701 Rupali salt-responsive DEGs with 1,741 variants identified in Genesis836 and 1,449 variants identified in Rupali. In addition, the presence of premature stop codons was detected in 35 genes in Rupali. This study provides valuable insights into the molecular regulation underpinning the physiological basis of salt tolerance in two chickpea genotypes and offers potential candidate genes for the improvement of salt tolerance in chickpeas.

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.

Medicago truncatula resources to study legume biology and symbiotic nitrogen fixation

Medicago truncatula is a chosen model for legumes towards deciphering fundamental legume biology, especially symbiotic nitrogen fixation. Current genomic resources for M. truncatula include a completed whole genome sequence information for R108 and Jemalong A17 accessions along with the sparse draft genome sequences for other 226 M. truncatula accessions. These genomic resources are complemented by the availability of mutant resources such as retrotransposon (Tnt1) insertion mutants in R108 and fast neutron bombardment (FNB) mutants in A17. In addition, several M. truncatula databases such as small secreted peptides (SSPs) database, transporter protein database, gene expression atlas, proteomic atlas, and metabolite atlas are available to the research community. This review describes these resources and provide information regarding how to access these resources.

Inheritance of Early and Late Ascochyta Blight Resistance in Wide Crosses of Chickpea

Chickpea (Cicer arietinum) is a globally important food legume but its yield is negatively impacted by the fungal pathogen Ascochyta blight (Ascochyta rabiei) causing necrotic lesions leading to plant death. Past studies have found that Ascochyta resistance is polygenic. It is important to find new resistance genes from the wider genepool of chickpeas. This study reports the inheritance of Ascochyta blight resistance of two wide crosses between the cultivar Gokce and wild chickpea accessions of C. reticulatum and C. echinospermum under field conditions in Southern Turkey. Following inoculation, infection damage was scored weekly for six weeks. The families were genotyped for 60 SNPs mapped to the reference genome for quantitative locus (QTL) mapping of resistance. Family lines showed broad resistance score distributions. A late responding QTL on chromosome 7 was identified in the C. reticulatum family and three early responding QTLs on chromosomes 2, 3, and 6 in the C. echinospermum family. Wild alleles mostly showed reduced disease severity, while heterozygous genotypes were most diseased. Interrogation of 200k bp genomic regions of the reference CDC Frontier genome surrounding QTLs identified nine gene candidates involved in disease resistance and cell wall remodeling. This study identifies new candidate chickpea Ascochyta blight resistance QTLs of breeding potential.

Mining legume germplasm for genetic gains: An Indian perspective

Legumes play a significant role in food and nutritional security and contribute to environmental sustainability. Although legumes are highly beneficial crops, it has not yet been possible to enhance their yield and production to a satisfactory level. Amid a rising population and low yield levels, per capita average legume consumption in India has fallen by 71% over the last 50 years, and this has led to protein-related malnutrition in a large segment of the Indian population, especially women and children. Several factors have hindered attempts to achieve yield enhancement in grain legumes, including biotic and abiotic pressures, a lack of good ideotypes, less amenability to mechanization, poorer responsiveness to fertilizer input, and a poor genetic base. Therefore, there is a need to mine the approximately 0.4 million ex situ collections of legumes that are being conserved in gene banks globally for identification of ideal donors for various traits. The Indian National Gene Bank conserves over 63,000 accessions of legumes belonging to 61 species. Recent initiatives have been undertaken in consortia mode with the aim of unlocking the genetic potential of ex situ collections and conducting large-scale germplasm characterization and evaluation analyses. We assume that large-scale phenotyping integrated with omics-based science will aid the identification of target traits and their use to enhance genetic gains. Additionally, in cases where the genetic base of major legumes is narrow, wild relatives have been evaluated, and these are being exploited through pre-breeding. Thus far, >200 accessions of various legumes have been registered as unique donors for various traits of interest.

Unravelling Differential DNA Methylation Patterns in Genotype Dependent Manner under Salinity Stress Response in Chickpea

DNA methylation is one of the epigenetic mechanisms that govern gene regulation in response to abiotic stress in plants. Here, we analyzed the role of epigenetic variations by exploring global DNA methylation and integrating it with differential gene expression in response to salinity stress in tolerant and sensitive chickpea genotypes. Genome-wide DNA methylation profiles showed higher CG methylation in the gene body regions and higher CHH methylation in the TE body regions. The analysis of differentially methylated regions (DMRs) suggested more hyper-methylation in response to stress in the tolerant genotype compared to the sensitive genotype. We observed higher enrichment of CG DMRs in genes and CHH DMRs in transposable elements (TEs). A positive correlation of gene expression with CG gene body methylation was observed. The enrichment analysis of DMR-associated differentially expressed genes revealed they are involved in biological processes, such as lateral root development, transmembrane transporter activity, GTPase activity, and regulation of gene expression. Further, a high correlation of CG methylation with CHG and CHH methylation under salinity stress was revealed, suggesting crosstalk among the methylation contexts. Further, we observed small RNA-mediated CHH hypermethylation in TEs. Overall, the interplay between DNA methylation, small RNAs, and gene expression provides new insights into the regulatory mechanism underlying salinity stress response in chickpeas.

Comprehensive Transcriptome Profiling Uncovers Molecular Mechanisms and Potential Candidate Genes Associated with Heat Stress Response in Chickpea

Chickpea (Cicer arietinum L.) production is highly susceptible to heat stress (day/night temperatures above 32/20 °C). Identifying the molecular mechanisms and potential candidate genes underlying heat stress response is important for increasing chickpea productivity. Here, we used an RNA-seq approach to investigate the transcriptome dynamics of 48 samples which include the leaf and root tissues of six contrasting heat stress responsive chickpea genotypes at the vegetative and reproductive stages of plant development. A total of 14,544 unique, differentially expressed genes (DEGs) were identified across different combinations studied. These DEGs were mainly involved in metabolic processes, cell wall remodeling, calcium signaling, and photosynthesis. Pathway analysis revealed the enrichment of metabolic pathways, biosynthesis of secondary metabolites, and plant hormone signal transduction, under heat stress conditions. Furthermore, heat-responsive genes encoding bHLH, ERF, WRKY, and MYB transcription factors were differentially regulated in response to heat stress, and candidate genes underlying the quantitative trait loci (QTLs) for heat tolerance component traits, which showed differential gene expression across tolerant and sensitive genotypes, were identified. Our study provides an important resource for dissecting the role of candidate genes associated with heat stress response and also paves the way for developing climate-resilient chickpea varieties for the future.

CicerSpTEdb2.0: An Upgrade of Cicer Species Transposable Elements Database

To meet the critical demand of LTR-RTs data-driven research, we updated the CicerSpTEdb database to version 2.0, which includes more accurate intact LTR-RT elements with annotation of internal domains. We also added the ability to BLAST against TEs of Cicer species. As a result, 3701 intact LTR-RTs were detected in the studied genomes, including 2840 Copia and 861 Gypsy elements. Of the 3701 intact LTR-RTs, 588 were in C. arietinum, including 475 Copia and 113 Gypsy. While 1373 were detected in C. reticulatum, including 1041 Copia and 332 Gypsy. Furthermore, 1740 were found in C. echinospermum, including 1324 Copia and 416 Gypsy. Based on LTR-RT clades, the analysis classified the 3701 identified intact LTR-RTs in the studied genomes as Ale (850), SIRE (740), unknown (455), Ikeros (323), Reina (290), Tork (290), Ivana (282), Tekay (197), Athila (128), TAR (99), CRM (31), and Ogre (16) elements. The newly updated CicerSpTEdb2.0 will be a valuable resource for TEs of Cicer species and their comparative genomics.Database URL: http://cicersptedb.easyomics.org/index.php.

Advancing Grain Legumes Domestication and Evolution Studies with Genomics

Grain legumes were domesticated in parallel with cereals in several regions of the world and formed the economic basis of early farming cultures. Since then, legumes have played a vital role in human and animal diets and in fostering agrobiodiversity. Increasing grain legume cultivation will be crucial to safeguard nutritional security and the resilience of agricultural ecosystems across the globe. A better understanding of the molecular underpinnings of domestication and crop evolution of grain legumes may be translated into practical approaches in modern breeding programs to stabilize yield, which is threatened by evolving pathogens and changing climates. During recent decades, domestication research in all crops has greatly benefited from the fast progress in genomic technologies. Yet still, many questions surrounding the domestication and diversification of legumes remain unanswered. In this review, we assess the potential of genomic approaches in grain legume research. We describe the centers of origin and the crucial domestication traits of grain legumes. In addition, we survey the effect of domestication on both above-ground and below-ground traits that have economic importance. Finally, we discuss open questions in grain legume domestication and diversification and outline how to bridge the gap between the preservation of historic crop diversity and their utilization in modern plant breeding.

Legume Pangenome: Status and Scope for Crop Improvement

In the last decade, legume genomics research has seen a paradigm shift due to advances in genome sequencing technologies, assembly algorithms, and computational genomics that enabled the construction of high-quality reference genome assemblies of major legume crops. These advances have certainly facilitated the identification of novel genetic variants underlying the traits of agronomic importance in many legume crops. Furthermore, these robust sequencing technologies have allowed us to study structural variations across the whole genome in multiple individuals and at the species level using 'pangenome analysis.' This review updates the progress of constructing pangenome assemblies for various legume crops and discusses the prospects for these pangenomes and how to harness the information to improve various traits of economic importance through molecular breeding to increase genetic gain in legumes and tackle the increasing global food crisis.

Chickpea (Cicer arietinum L.) Biology and Biotechnology: From Domestication to Biofortification and Biopharming

Chickpea (Cicer arietinum L.), the world's second most consumed legume crop, is cultivated in more than 50 countries around the world. It is a boon for diabetics and is an excellent source of important nutrients such as vitamins A, C, E, K, B1-B3, B5, B6, B9 and minerals (Fe, Zn, Mg and Ca) which all have beneficial effects on human health. By 2050, the world population can cross 9 billion, and in order to feed the teaming millions, chickpea production should also be increased, as it is a healthy alternative to wheat flour and a boon for diabetics. Moreover, it is an important legume that is crucial for food, nutrition, and health security and the livelihood of the small-scale farmers with poor resources, in developing countries. Although marvelous improvement has been made in the development of biotic and abiotic stress-resistant varieties, still there are many lacunae, and to fulfill that, the incorporation of genomic technologies in chickpea breeding (genomics-assisted breeding, high-throughput and precise-phenotyping and implementation of novel breeding strategies) will facilitate the researchers in developing high yielding, climate resilient, water use efficient, salt-tolerant, insect/pathogen resistant varieties, acceptable to farmers, consumers, and industries. This review focuses on the origin and distribution, nutritional profile, genomic studies, and recent updates on crop improvement strategies for combating abiotic and biotic stresses in chickpea.

An integrated transcriptome mapping the regulatory network of coding and long non-coding RNAs provides a genomics resource in chickpea

Large-scale transcriptome analysis can provide a systems-level understanding of biological processes. To accelerate functional genomic studies in chickpea, we perform a comprehensive transcriptome analysis to generate full-length transcriptome and expression atlas of protein-coding genes (PCGs) and long non-coding RNAs (lncRNAs) from 32 different tissues/organs via deep sequencing. The high-depth RNA-seq dataset reveal expression dynamics and tissue-specificity along with associated biological functions of PCGs and lncRNAs during development. The coexpression network analysis reveal modules associated with a particular tissue or a set of related tissues. The components of transcriptional regulatory networks (TRNs), including transcription factors, their cognate cis-regulatory motifs, and target PCGs/lncRNAs that determine developmental programs of different tissues/organs, are identified. Several candidate tissue-specific and abiotic stress-responsive transcripts associated with quantitative trait loci that determine important agronomic traits are also identified. These results provide an important resource to advance functional/translational genomic and genetic studies during chickpea development and environmental conditions.

A full-length transcriptome and expression atlas of protein-coding genes and long non-coding RNAs is generated in chickpea. Components of transcriptional regulatory networks and candidate tissue-specific transcripts associated with quantitative trait loci are identified.

SSR markers in revealing extent of genetic diversity and phylogenetic relationships among chickpea core collection accessions for Western Himalayas

BACKGROUND

The exploration of genetic diversity is the key source of germplasm conservation and potential to broaden its genetic base. The globally growing demand for chickpea suggests superior/climate-resilient varieties, which in turn necessitates the germplasm characterization to unravel underlying genetic variation.

METHODOLOGY AND RESULTS

A chickpea core collection comprising of diverse 192 accessions which include cultivated Cicer arietinum, and wild C. reticulatum, C. echinospermum, and C. microphyllum species were investigated to analyze their genetic diversity and relationship, by assaying 33 unlinked simple sequence repeat (SSR) markers. The results amplified a total of 323 alleles (Na), ranging from 2 to 8 with an average of 4.25 alleles per locus. Expected heterozygosity (He) differed from 0.46 to 0.86 with an average of 0.68. Polymorphic information content (PIC) ranged from 0.73 to 0.98 with an average of 0.89. Analysis of molecular variance (AMOVA) showed that most of the variation was among individuals (87%). Cluster analysis resulted in the formation of four distinct clusters. Cluster I represented all cultivated and clusters II, III, and IV comprised a heterogeneous group of cultivated and wild chickpea accessions.

CONCLUSION

We report considerable diversity and greater resolving power of SSR markers for assessing variability and interrelationship among the chickpea accessions. The chickpea core is expected to be an efficient resource for breeders for broadening the chickpea genetic base and could be useful for selective breeding of desirable traits and in the identification of target genes for genomics-assisted breeding.

Genome-wide association mapping of nutritional traits for designing superior chickpea varieties

Micronutrient malnutrition is a serious concern in many parts of the world; therefore, enhancing crop nutrient content is an important challenge. Chickpea (Cicer arietinum L.), a major food legume crop worldwide, is a vital source of protein and minerals in the vegetarian diet. This study evaluated a diverse set of 258 chickpea germplasm accessions for 12 key nutritional traits. A significant variation was observed for several nutritional traits, including crude protein (16.56-24.64/100 g), β-Carotene (0.003-0.104 mg/100 g), calcium (60.69-176.55 mg/100 g), and folate (0.413-6.537 mg/kg). These data, combined with the available whole-genome sequencing data for 318,644 SNPs, were used in genome-wide association studies comprising single-locus and multi-locus models. We also explored the effect of varying the minor allele frequency (MAF) levels and heterozygosity. We identified 62 significant marker-trait associations (MTAs) explaining up to 28.63% of the phenotypic variance (PV), of which nine were localized within genes regulating G protein-coupled receptor signaling pathway, proteasome assembly, intracellular signal transduction, and oxidation-reduction process, among others. The significant effect MTAs were located primarily on Ca1, Ca3, Ca4, and Ca6. Importantly, varying the level of heterozygosity was found to significantly affect the detection of associations contributing to traits of interest. We further identified seven promising accessions (ICC10399, ICC1392, ICC1710, ICC2263, ICC1431, ICC4182, and ICC16915) with superior agronomic performance and high nutritional content as potential donors for developing nutrient-rich, high-yielding chickpea varieties. Validation of the significant MTAs with higher PV could identify factors controlling the nutrient acquisition and facilitate the design of biofortified chickpeas for the future.

Genetic mapping of QTLs for drought tolerance in chickpea (Cicer arietinum L.)

Chickpea yield is severely affected by drought stress, which is a complex quantitative trait regulated by multiple small-effect genes. Identifying genomic regions associated with drought tolerance component traits may increase our understanding of drought tolerance mechanisms and assist in the development of drought-tolerant varieties. Here, a total of 187 F8 recombinant inbred lines (RILs) developed from an interspecific cross between drought-tolerant genotype GPF 2 (Cicer arietinum) and drought-sensitive accession ILWC 292 (C. reticulatum) were evaluated to identify quantitative trait loci (QTLs) associated with drought tolerance component traits. A total of 21 traits, including 12 morpho-physiological traits and nine root-related traits, were studied under rainfed and irrigated conditions. Composite interval mapping identified 31 QTLs at Ludhiana and 23 QTLs at Faridkot locations for morphological and physiological traits, and seven QTLs were identified for root-related traits. QTL analysis identified eight consensus QTLs for six traits and five QTL clusters containing QTLs for multiple traits on linkage groups CaLG04 and CaLG06. The identified major QTLs and genomic regions associated with drought tolerance component traits can be introgressed into elite cultivars using genomics-assisted breeding to enhance drought tolerance in chickpea.

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

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

The Genome of the Mimosoid Legume Prosopis cineraria, a Desert Tree

The mimosoid legumes are a clade of ~40 genera in the Caesalpinioideae subfamily of the Fabaceae that grow in tropical and subtropical regions. Unlike the better studied Papilionoideae, there are few genomic resources within this legume group. The tree Prosopis cineraria is native to the Near East and Indian subcontinent, where it thrives in very hot desert environments. To develop a tool to better understand desert plant adaptation mechanisms, we sequenced the P. cineraria genome to near-chromosomal assembly, with a total sequence length of ~691 Mb. We predicted 77,579 gene models (76,554 CDS, 361 rRNAs and 664 tRNAs) from the assembled genome, among them 55,325 (~72%) protein-coding genes that were functionally annotated. This genome was found to consist of over 58% repeat sequences, primarily long terminal repeats (LTR-)-retrotransposons. We find an expansion of terpenoid metabolism genes in P. cineraria and its relative Prosopis alba, but not in other legumes. We also observed an amplification of NBS-LRR disease-resistance genes correlated with LTR-associated retrotransposition, and identified 410 retrogenes with an active burst of chimeric retrogene creation that approximately occurred at the same time of divergence of P. cineraria from a common lineage with P. alba~23 Mya. These retrogenes include many biotic defense responses and abiotic stress stimulus responses, as well as the early Nodulin 93 gene. Nodulin 93 gene amplification is consistent with an adaptive response of the species to the low nitrogen in arid desert soil. Consistent with these results, our differentially expressed genes show a tissue specific expression of isoprenoid pathways in shoots, but not in roots, as well as important genes involved in abiotic salt stress in both tissues. Overall, the genome sequence of P. cineraria enriches our understanding of the genomic mechanisms of its disease resistance and abiotic stress tolerance. Thus, it is a very important step in crop and legume improvement.

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.

Genome-wide analyses of the mung bean NAC gene family reveals orthologs, co-expression networking and expression profiling under abiotic and biotic stresses

BACKGROUND

Mung bean is a short-duration and essential food crop owing to its cash prominence in Asia. Mung bean seeds are rich in protein, fiber, antioxidants, and phytonutrients. The NAC transcription factors (TFs) family is a large plant-specific family, participating in tissue development regulation and abiotic and biotic stresses.

RESULTS

In this study, we perform genome-wide comparisons of VrNAC with their homologs from Arabidopsis. We identified 81 NAC transcription factors (TFs) in mung bean genome and named as per their chromosome location. A phylogenetic analysis revealed that VrNACs are broadly distributed in nine groups. Moreover, we identified 20 conserved motifs across the VrNACs highlighting their roles in different biological process. Based on the gene structure of the putative VrNAC and segmental duplication events might be playing a vital role in the expansion of mung bean genome. A comparative phylogenetic analysis of mung bean NAC together with homologs from Arabidopsis allowed us to classify NAC genes into 13 groups, each containing several orthologs and paralogs. Gene ontology (GO) analysis categorized the VrNACs into biological process, cellular components and molecular functions, explaining the functions in different plant physiology processes. A gene co-expression network analysis identified 173 genes involved in the transcriptional network of putative VrNAC genes. We also investigated how miRNAs potentially target VrNACs and shape their interactions with proteins. VrNAC1.4 (Vradi01g03390.1) was targeted by the Vra-miR165 family, including 9 miRNAs. Vra-miR165 contributes to leaf development and drought tolerance. We also performed qRT-PCR on 22 randomly selected VrNAC genes to assess their expression patterns in the NM-98 genotype, widely known for being tolerant to drought and bacterial leaf spot disease.

CONCLUSIONS

This genome-wide investigation of VrNACs provides a unique resource for further detailed investigations aimed at predicting orthologs functions and what role the play under abiotic and biotic stress, with the ultimate aim to improve mung bean production under diverse environmental conditions.

Prospects and potentials of underutilized leafy Amaranths as vegetable use for health-promotion

Climate change causes environmental variation worldwide, which is one of the most serious threats to global food security. In addition, more than 2 billion people in the world are reported to suffer from serious malnutrition, referred to as 'hidden hunger.' Dependence on only a few crops could lead to the loss of genetic diversity and high fragility of crop breeding in systems adapting to global scale climate change. The exploitation of underutilized species and genetic resources, referred to as orphan crops, could be a useful approach for resolving the issue of adaptability to environmental alteration, biodiversity preservation, and improvement of nutrient quality and quantity to ensure food security. Moreover, the use of these alternative crops will help to increase the human health benefits and the income of farmers in developing countries. In this review, we highlight the potential of orphan crops, especially amaranths, for use as vegetables and health-promoting nutritional components. This review highlights promising diversified sources of amaranth germplasms, their tolerance to abiotic stresses, and their nutritional, phytochemical, and antioxidant values for vegetable purposes. Betalains (betacyanins and betaxanthins), unique antioxidant components in amaranth vegetables, are also highlighted regarding their chemodiversity across amaranth germplasms and their stability and degradation. In addition, we discuss the physiological functions, antioxidant, antilipidemic, anticancer, and antimicrobial activities, as well as the biosynthesis pathway, molecular, biochemical, genetics, and genomic mechanisms of betalains in detail.

A Comprehensive Review on Chickpea (Cicer arietinum L.) Breeding for Abiotic Stress Tolerance and Climate Change Resilience

Chickpea is one of the most important pulse crops worldwide, being an excellent source of protein. It is grown under rain-fed conditions averaging yields of 1 t/ha, far from its potential of 6 t/ha under optimum conditions. The combined effects of heat, cold, drought, and salinity affect species productivity. In this regard, several physiological, biochemical, and molecular mechanisms are reviewed to confer tolerance to abiotic stress. A large collection of nearly 100,000 chickpea accessions is the basis of breeding programs, and important advances have been achieved through conventional breeding, such as germplasm introduction, gene/allele introgression, and mutagenesis. In parallel, advances in molecular biology and high-throughput sequencing have allowed the development of specific molecular markers for the genus Cicer, facilitating marker-assisted selection for yield components and abiotic tolerance. Further, transcriptomics, proteomics, and metabolomics have permitted the identification of specific genes, proteins, and metabolites associated with tolerance to abiotic stress of chickpea. Furthermore, some promising results have been obtained in studies with transgenic plants and with the use of gene editing to obtain drought-tolerant chickpea. Finally, we propose some future lines of research that may be useful to obtain chickpea genotypes tolerant to abiotic stress in a scenario of climate change.

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.

Beyond a reference genome: pangenomes and population genomics of underutilized and orphan crops for future food and nutrition security

Underutilized crops are, by definition, under-researched compared to staple crops yet come with traits that may be especially important given climate change and the need to feed a globally increasing population. These crops are often stress-tolerant, and this combined with unique and beneficial nutritional profiles. Whilst progress is being made by generating reference genome sequences, in this Tansley Review, we show how this is only the very first step. We advocate that going 'beyond a reference genome' should be a priority, as it is only at this stage one can identify the specific genes and the adaptive alleles that underpin the valuable traits. We sum up how population genomic and pangenomic approaches have led to the identification of stress- and disease-tolerant alleles in staple crops and compare this to the small number of examples from underutilized crops. We also demonstrate how previously underutilized crops have benefitted from genomic advances and that many breeding targets in underutilized crops are often well studied in staple crops. This cross-crop population-level resequencing could lead to an understanding of the genetic basis of adaptive traits in underutilized crops. This level of investment may be crucial for fully understanding the value of these crops before they are lost.

Insights into the gene and protein structures of the CaSWEET family members in chickpea (Cicer arietinum), and their gene expression patterns in different organs under various stress and abscisic acid treatments

'Sugars Will Eventually be Exported Transporters' (SWEETs) are a group of sugar transporters that play crucial roles in various biological processes, particularly plant stress responses. However, no information is available yet for the CaSWEET family in chickpea. Here, we identified all putative CaSWEET members in chickpea, and obtained their major characteristics, including physicochemical patterns, chromosomal distribution, subcellular localization, gene organization, conserved motifs and three-dimensional protein structures. Subsequently, we explored available transcriptome data to compare spatiotemporal transcript abundance of CaSWEET genes in various major organs. Finally, we studied the changes in their transcript levels in leaves and/or roots following dehydration and exogenous abscisic acid treatments using RT-qPCR to obtain valuable information underlying their potential roles in chickpea responses to water-stress conditions. Our results provide the first insights into the characteristics of the CaSWEET family members and a foundation for further functional characterizations of selected candidate genes for genetic engineering of chickpea.

Defining the mutation sites in chickpea nodulation mutants PM233 and PM405

BACKGROUND

Like most legumes, chickpeas form specialized organs called root nodules. These nodules allow for a symbiotic relationship with rhizobium bacteria. The rhizobia provide fixed atmospheric nitrogen to the plant in a usable form. It is of both basic and practical interest to understand the host plant genetics of legume root nodulation. Chickpea lines PM233 and PM405, which harbor the mutationally identified nodulation genes rn1 and rn4, respectively, both display nodulation-deficient phenotypes. Previous investigators identified the rn1 mutation with the chickpea homolog of Medicago truncatula nodulation gene NSP2, but were unable to define the mutant rn1 allele. We used Illumina and Nanopore sequencing reads to attempt to identify and characterize candidate mutation sites responsible for the PM233 and PM405 phenotypes.

RESULTS

We aligned Illumina reads to the available desi chickpea reference genome, and did a de novo contig assembly of Nanopore reads. In mutant PM233, the Nanopore contigs allowed us to identify the breakpoints of a ~ 35 kb deleted region containing the CaNSP2 gene, the Medicago truncatula homolog of which is involved in nodulation. In mutant PM405, we performed variant calling in read alignments and identified 10 candidate mutations. Genotyping of a segregating progeny population narrowed that pool down to a single candidate gene which displayed homology to M. truncatula nodulation gene NIN.

CONCLUSIONS

We have characterized the nodulation mutation sites in chickpea mutants PM233 and PM405. In mutant PM233, the rn1 mutation was shown to be due to deletion of the entire CaNSP2 nodulation gene, while in mutant PM405 the rn4 mutation was due to a single base deletion resulting in a frameshift mutation between the predicted RWP-RK and PB1 domains of the NIN nodulation gene. Critical to characterization of the rn1 allele was the generation of Nanopore contigs for mutant PM233 and its wild type parent ICC 640, without which the deletional boundaries could not be defined. Our results suggest that efforts of prior investigators were hampered by genomic misassemblies in the CaNSP2 region of both the desi and kabuli reference genomes.

Chickpea Roots Undergoing Colonisation by Phytophthora medicaginis Exhibit Opposing Jasmonic Acid and Salicylic Acid Accumulation and Signalling Profiles to Leaf Hemibiotrophic Models

Hemibiotrophic pathogens cause significant losses within agriculture, threatening the sustainability of food systems globally. These microbes colonise plant tissues in three phases: a biotrophic phase followed by a biotrophic-to-necrotrophic switch phase and ending with necrotrophy. Each of these phases is characterized by both common and discrete host transcriptional responses. Plant hormones play an important role in these phases, with foliar models showing that salicylic acid accumulates during the biotrophic phase and jasmonic acid/ethylene responses occur during the necrotrophic phase. The appropriateness of this model to plant roots has been challenged in recent years. The need to understand root responses to hemibiotrophic pathogens of agronomic importance necessitates further research. In this study, using the root hemibiotroph Phytophthora medicaginis, we define the duration of each phase of pathogenesis in Cicer arietinum (chickpea) roots. Using transcriptional profiling, we demonstrate that susceptible chickpea roots display some similarities in response to disease progression as previously documented in leaf plant–pathogen hemibiotrophic interactions. However, our transcriptomic results also show that chickpea roots do not conform to the phytohormone responses typically found in leaf colonisation by hemibiotrophs. We found that quantified levels of salicylic acid concentrations in root tissues decreased significantly during biotrophy while jasmonic acid concentrations were significantly induced. This study demonstrated that a wider spectrum of plant species should be investigated in the future to understand the physiological changes in plants during colonisation by soil-borne hemibiotrophic pathogens before we can better manage these economically important microbes.

Genomics Associated Interventions for Heat Stress Tolerance in Cool Season Adapted Grain Legumes

Cool season grain legumes occupy an important place among the agricultural crops and essentially provide multiple benefits including food supply, nutrition security, soil fertility improvement and revenue for farmers all over the world. However, owing to climate change, the average temperature is steadily rising, which negatively affects crop performance and limits their yield. Terminal heat stress that mainly occurred during grain development phases severely harms grain quality and weight in legumes adapted to the cool season, such as lentils, faba beans, chickpeas, field peas, etc. Although, traditional breeding approaches with advanced screening procedures have been employed to identify heat tolerant legume cultivars. Unfortunately, traditional breeding pipelines alone are no longer enough to meet global demands. Genomics-assisted interventions including new-generation sequencing technologies and genotyping platforms have facilitated the development of high-resolution molecular maps, QTL/gene discovery and marker-assisted introgression, thereby improving the efficiency in legumes breeding to develop stress-resilient varieties. Based on the current scenario, we attempted to review the intervention of genomics to decipher different components of tolerance to heat stress and future possibilities of using newly developed genomics-based interventions in cool season adapted grain legumes.

Discovery of DNA polymorphisms via whole genome resequencing and their functional relevance in salinity stress response in chickpea

Salinity stress is one of the major constraints for plant growth and yield. The salinity stress response of different genotypes of crop plants may largely be governed by DNA polymorphisms. To determine the molecular genetic factors involved in salinity stress tolerance in chickpea, we performed a whole genome resequencing data analysis of three each of salinity-sensitive and salinity-tolerant genotypes. A total of 6173 single nucleotide polymorphisms and 920 insertions and deletions differentiating the chickpea genotypes with contrasting salinity stress responses were identified. Gene ontology analysis revealed the enrichment of functional terms related to stress response and development among the genes harboring DNA polymorphisms in their promoter and/or coding regions. DNA polymorphisms located within the cis-regulatory motifs of the quantitative trait loci (QTL)-associated and abiotic stress related genes were identified, which may influence salinity stress response via modulating binding affinity of the transcription factors. Several genes including QTL-associated and abiotic stress response related genes harboring DNA polymorphisms exhibited differential expression in response to salinity stress especially at the reproductive stage of development in the salinity-tolerant genotype. Furthermore, effects of non-synonymous DNA polymorphisms on mutational sensitivity and structural integrity of the encoded proteins by the candidate QTL-associated and abiotic stress response related genes were revealed. The results suggest that DNA polymorphisms may determine salinity stress response via influencing differential gene expression in genotype and/or stage-dependent manner. Altogether, we provide a high-quality set of DNA polymorphisms and candidate genes that may govern salinity stress tolerance in chickpea.

Evolutionary and functional analysis of two-component system in chickpea reveals CaRR13, a TypeB RR, as positive regulator of symbiosis

The critical role of cytokinin in early nodulation in legumes is well known. In our study, exogenous cytokinin application to roots of the important crop legume, chickpea (Cicer arietinum L.), led to the formation of pseudo-nodules even in the absence of rhizobia. Hence, a genome-wide analysis of the cytokinin signalling, two-component system (TCS) genes, was conducted in chickpea, Medicago and Cajanus cajan. The integrated phylogenetic, evolutionary and expression analysis of the TCS genes was carried out, which revealed that histidine kinases (HKs) were highly conserved, whereas there was diversification leading to neofunctionalization at the level of response regulators (RRs) especially the TypeB RRs. Further, the functional role of the CaHKs in nodulation was established by complementation of the sln1Δ mutant of yeast and cre1 mutants of (Medicago) which led to restoration of the nodule-deficient phenotype. Additionally, the highest expressing TypeB RR of chickpea, CaRR13, was functionally characterized. Its localization in the nucleus and its Y1H assay-based interaction with the promoter of the early nodulation gene CaNSP2 indicated its role as a transcription factor regulating early nodulation. Overexpression, RNAi lines and complementation of cre1 mutants with CaRR13 revealed its critical involvement as an important signalling molecule regulating early events of nodule organogenesis in chickpea.

Chickpea tolerance to temperature stress: Status and opportunity for improvement

Chickpea is a globally important commercial crop and a key source of protein for vegetarian populations. Though chickpea was domesticated at least 3000 years ago, research into abiotic stress tolerance has been limited compared to cereal crops such as wheat. This review investigates the impacts of heat stress on chickpea, focusing on reproductive development. The fertilisation process is particularly sensitive to environmental stress, such as drought and heat that can reduce yields by up to 70%. Current research has largely focused on breeding cultivars that reach maturity faster to avoid stress rather than true thermotolerance and little is known of the impact of heat on cellular processes. This review suggests that there is ample variation within the chickpea gene pool for selective breeding to achieve improved abiotic stress tolerance. Rates of genetic progress will improve once key QTL are identified and the link between thermotolerance and pollen viability confirmed. Other benefits may arise from better understanding of heat shock proteins and molecular chaperones and their role in the protection of reproductive processes.

CicerSpTEdb: A web-based database for high-resolution genome-wide identification of transposable elements in Cicer species

Recently, Cicer species have experienced increased research interest due to their economic importance, especially in genetics, genomics, and crop improvement. The Cicer arietinum, Cicer reticulatum, and Cicer echinospermum genomes have been sequenced and provide valuable resources for trait improvement. Since the publication of the chickpea draft genome, progress has been made in genome assembly, functional annotation, and identification of polymorphic markers. However, work is still needed to identify transposable elements (TEs) and make them available for researchers. In this paper, we present CicerSpTEdb, a comprehensive TE database for Cicer species that aims to improve our understanding of the organization and structural variations of the chickpea genome. Using structure and homology-based methods, 3942 C. echinospermum, 3579 C. reticulatum, and 2240 C. arietinum TEs were identified. Comparisons between Cicer species indicate that C. echinospermum has the highest number of LTR-RT and hAT TEs. C. reticulatum has more Mutator, PIF Harbinger, Tc1 Mariner, and CACTA TEs, while C. arietinum has the highest number of Helitron. CicerSpTEdb enables users to search and visualize TEs by location and download their results. The database will provide a powerful resource that can assist in developing TE target markers for molecular breeding and answer related biological questions. Database URL: http://cicersptedb.easyomics.org/index.php.

Aldehyde Dehydrogenase 3 Is an Expanded Gene Family with Potential Adaptive Roles in Chickpea

Legumes play an important role in ensuring food security, improving nutrition and enhancing ecosystem resilience. Chickpea is a globally important grain legume adapted to semi-arid regions under rain-fed conditions. A growing body of research shows that aldehyde dehydrogenases (ALDHs) represent a gene class with promising potential for plant adaptation improvement. Aldehyde dehydrogenases constitute a superfamily of proteins with important functions as ‘aldehyde scavengers’ by detoxifying aldehydes molecules, and thus play important roles in stress responses. We performed a comprehensive study of the ALDH superfamily in the chickpea genome and identified 27 unique ALDH loci. Most chickpea ALDHs originated from duplication events and the ALDH3 gene family was noticeably expanded. Based on the physical locations of genes and sequence similarities, our results suggest that segmental duplication is a major driving force in the expansion of the ALDH family. Supported by expression data, the findings of this study offer new potential target genes for improving stress tolerance in chickpea that will be useful for breeding programs.

Genome-wide identification and functional prediction of salt- stress related long non-coding RNAs (lncRNAs) in chickpea (Cicer arietinum L.)

LncRNAs (long noncoding RNAs) are 200 bp length crucial RNA molecules, lacking coding potential and having important roles in regulating gene expression, particularly in response to abiotic stresses. In this study, we identified salt stress-induced lncRNAs in chickpea roots and predicted their intricate regulatory roles. A total of 3452 novel lncRNAs were identified to be distributed across all 08 chickpea chromosomes. On comparing salt-tolerant (ICCV 10, JG 11) and salt-sensitive cultivars (DCP 92-3, Pusa 256), 4446 differentially expressed lncRNAs were detected under various salt  treatments. We predicted 3373 lncRNAs to be regulating their target genes in cis regulating manner and 80 unique lncRNAs were observed as interacting with 136 different miRNAs, as eTMs (endogenous target mimic) targets of miRNAs and implicated them in the regulatory network of salt stress response. Functional analysis of these lncRNA revealed their association in targeting salt stress response-related genes like potassium transporter, transporter family genes, serine/threonine-protein kinase, aquaporins like TIP1-2, PIP2-5 and transcription factors like, AP2, NAC, bZIP, ERF, MYB and WRKY. Furthermore, about 614 lncRNA-SSRs (simple sequence repeats) were identified as a new generation of molecular markers with higher efficiency and specificity in chickpea. Overall, these findings will pave the understanding of comprehensive functional role of potential lncRNAs, which can help in providing insight into the molecular mechanism of salt tolerance in chickpea.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01093-0.

Identification of stable heat tolerance QTLs using inter-specific recombinant inbred line population derived from GPF 2 and ILWC 292

Heat stress during reproductive stages has been leading to significant yield losses in chickpea (Cicer arietinum L.). With an aim of identifying the genomic regions or QTLs responsible for heat tolerance, 187 F8 recombinant inbred lines (RILs) derived from the cross GPF 2 (heat tolerant) × ILWC 292 (heat sensitive) were evaluated under late-sown irrigated (January-May) and timely-sown irrigated environments (November-April) at Ludhiana and Faridkot in Punjab, India for 13 heat tolerance related traits. The pooled ANOVA for both locations for the traits namely days to germination (DG), days to flowering initiation (DFI), days to 50% flowering (DFF), days to 100% flowering (DHF), plant height (PH), pods per plant (NPP), biomass (BIO), grain yield (YLD), 100-seed weight (HSW), harvest index (HI), membrane permeability index (MPI), relative leaf water content (RLWC) and pollen viability (PV)) showed a highly significant difference in RILs. The phenotyping data coupled with the genetic map comprising of 1365 ddRAD-Seq based SNP markers were used for identifying the QTLs for heat tolerance. Composite interval mapping provided a total of 28 and 23 QTLs, respectively at Ludhiana and Faridkot locations. Of these, 13 consensus QTLs for DG, DFI, DFF, DHF, PH, YLD, and MPI have been identified at both locations. Four QTL clusters containing QTLs for multiple traits were identified on the same genomic region at both locations. Stable QTLs for days to flowering can be one of the major factors for providing heat tolerance as early flowering has an advantage of more seed setting due to a comparatively longer reproductive period. Identified QTLs can be used in genomics-assisted breeding to develop heat stress-tolerant high yielding chickpea cultivars.

Genetics and breeding for climate change in Orphan crops

Climate change is rapidly changing how we live, what we eat and produce, the crops we breed and the target traits. Previously underutilized orphan crops that are climate resilient are receiving much attention from the crops research community, as they are often the only crops left in the field after periods of extreme weather conditions. There are several orphan crops with incredible resilience to biotic and abiotic stresses. Some are nutritious, while others provide good sources of biofuel, medicine and other industrial raw materials. Despite these benefits, orphan crops are still lacking in important genetic and genomic resources that could be used to fast track their improvement and make their production profitable. Progress has been made in generating draft genomes of at least 28 orphan crops over the last decade, thanks to the reducing cost of sequencing. The implementation of a structured breeding program that takes advantage of additional modern crop improvement tools such as genomic selection, speed breeding, genome editing, high throughput phenotyping and breeding digitization would make rapid improvement of these orphan crops possible, but would require coordinated research investment. Other production challenges such as lack of adequate germplasm conservation, poor/non-existent seed systems and agricultural extension services, as well as poor marketing channels will also need to be improved if orphan crops were to be profitable. We review the importance of breeding orphan crops under the increasing effects of climate change, highlight existing gaps that need to be addressed and share some lessons to be learned from major crops.

Comprehensive analysis and identification of drought-responsive candidate NAC genes in three semi-arid tropics (SAT) legume crops

Background

Chickpea, pigeonpea, and groundnut are the primary legume crops of semi-arid tropics (SAT) and their global productivity is severely affected by drought stress. The plant-specific NAC (NAM - no apical meristem, ATAF - Arabidopsis transcription activation factor, and CUC - cup-shaped cotyledon) transcription factor family is known to be involved in majority of abiotic stresses, especially in the drought stress tolerance mechanism. Despite the knowledge available regarding NAC function, not much information is available on NAC genes in SAT legume crops.

Results

In this study, genome-wide NAC proteins – 72, 96, and 166 have been identified from the genomes of chickpea, pigeonpea, and groundnut, respectively, and later grouped into 10 clusters in chickpea and pigeonpea, while 12 clusters in groundnut. Phylogeny with well-known stress-responsive NACs in Arabidopsis thaliana, Oryza sativa (rice), Medicago truncatula, and Glycine max (soybean) enabled prediction of putative stress-responsive NACs in chickpea (22), pigeonpea (31), and groundnut (33). Transcriptome data revealed putative stress-responsive NACs at various developmental stages that showed differential expression patterns in the different tissues studied. Quantitative real-time PCR (qRT-PCR) was performed to validate the expression patterns of selected stress-responsive, Ca_NAC (Cicer arietinum - 14), Cc_NAC (Cajanus cajan - 15), and Ah_NAC (Arachis hypogaea - 14) genes using drought-stressed and well-watered root tissues from two contrasting drought-responsive genotypes of each of the three legumes. Based on expression analysis, Ca_06899, Ca_18090, Ca_22941, Ca_04337, Ca_04069, Ca_04233, Ca_12660, Ca_16379, Ca_16946, and Ca_21186; Cc_26125, Cc_43030, Cc_43785, Cc_43786, Cc_22429, and Cc_22430; Ah_ann1.G1V3KR.2, Ah_ann1.MI72XM.2, Ah_ann1.V0X4SV.1, Ah_ann1.FU1JML.2, and Ah_ann1.8AKD3R.1 were identified as potential drought stress-responsive candidate genes.

Conclusion

As NAC genes are known to play role in several physiological and biological activities, a more comprehensive study on genome-wide identification and expression analyses of the NAC proteins have been carried out in chickpea, pigeonpea and groundnut. We have identified a total of 21 potential drought-responsive NAC genes in these legumes. These genes displayed correlation between gene expression, transcriptional regulation, and better tolerance against drought. The identified candidate genes, after validation, may serve as a useful resource for molecular breeding for drought tolerance in the SAT legume crops.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07602-5.

Molecular mapping of quantitative trait loci for ascochyta blight and botrytis grey mould resistance in an inter-specific cross in chickpea (Cicer arietinum L.) using genotyping by sequencing

Ascochyta blight (AB) and botrytis grey mould (BGM) are the most devastating fungal diseases of chickpea worldwide. The wild relative of chickpea, C. reticulatum acc. ILWC 292 was found resistant to BGM whereas, GPF2 (Cicer arietinum L.) is resistant to AB. A total of 187 F₈ Recombinant Inbred Lines (RILs) developed from an inter-specific cross of GPF2 × C. reticulatum acc. ILWC 292 were used to identify quantitative trait loci (QTLs) responsible for resistance to AB and BGM. RILs along with parents were evaluated under artificial epiphytotic field/laboratory conditions for two years. Highly significant differences (P < 0.001) were observed for reaction to both pathogens in both years. Parents and RILs were genotyped-by-sequencing to identify genome wide single nucleotide polymorphism (SNPs). A total of 1365 filtered and parental polymorphic SNPs were used for linkage map construction, of which, 673 SNPs were arranged on eight linkage groups. Composite interval mapping revealed three QTLs for AB and four QTLs for BGM resistance. Out of which, two QTLs for AB and three QTLs for BGM were consistent in both years. These QTLs can be targeted for further fine mapping for deployment of resistance to AB and BGM in elite chickpea cultivars using marker-assisted-selection.

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.