Resequencing of 429 chickpea accessions from 45 countries provides insights into genome diversity, domestication and agronomic traits

Chromosome-level reference genome and resequencing of 322 accessions reveal evolution, genomic imprint and key agronomic traits in adzuki bean

Adzuki bean (Vigna angularis) is an important legume crop cultivated in over 30 countries worldwide. We developed a high-quality chromosome-level reference genome of adzuki bean cultivar Jingnong6 by combining PacBio Sequel long-read sequencing with short-read and Hi-C technologies. The assembled genome covers 97.8% of the adzuki bean genome with a contig N50 of approximately 16 Mb and a total of 32 738 protein-coding genes. We also generated a comprehensive genome variation map of adzuki bean by whole-genome resequencing (WGRS) of 322 diverse adzuki beans accessions including both wild and cultivated. Furthermore, we have conducted comparative genomics and a genome-wide association study (GWAS) on key agricultural traits to investigate the evolution and domestication. GWAS identified several candidate genes, including VaCycA3;1, VaHB15, VaANR1 and VaBm, that exhibited significant associations with domestication traits. Furthermore, we conducted functional analyses on the roles of VaANR1 and VaBm in regulating seed coat colour. We provided evidence for the highest genetic diversity of wild adzuki (Vigna angularis var. nipponensis) in China with the presence of the most original wild adzuki bean, and the occurrence of domestication process facilitating transition from wild to cultigen. The present study elucidates the genetic basis of adzuki bean domestication traits and provides crucial genomic resources to support future breeding efforts in adzuki bean.

Unlocking the nutritional potential of chickpea: strategies for biofortification and enhanced multinutrient quality

Chickpea (Cicer arietinum L.) is a vital grain legume, offering an excellent balance of protein, carbohydrates, fats, fiber, essential micronutrients, and vitamins that can contribute to addressing the global population's increasing food and nutritional demands. Chickpea protein offers a balanced source of amino acids with high bioavailability. Moreover, due to its balanced nutrients and affordable price, chickpea is an excellent alternative to animal protein, offering a formidable tool for combating hidden hunger and malnutrition, particularly prevalent in low-income countries. This review examines chickpea's nutritional profile, encompassing protein, amino acids, carbohydrates, fatty acids, micronutrients, vitamins, antioxidant properties, and bioactive compounds of significance in health and pharmaceutical domains. Emphasis is placed on incorporating chickpeas into diets for their myriad health benefits and nutritional richness, aimed at enhancing human protein and micronutrient nutrition. We discuss advances in plant breeding and genomics that have facilitated the discovery of diverse genotypes and key genomic variants/regions/quantitative trait loci contributing to enhanced macro- and micronutrient contents and other quality parameters. Furthermore, we explore the potential of innovative breeding tools such as CRISPR/Cas9 in enhancing chickpea's nutritional profile. Envisioning chickpea as a nutritionally smart crop, we endeavor to safeguard food security, combat hunger and malnutrition, and promote dietary diversity within sustainable agrifood systems.

Cicer super-pangenome provides insights into species evolution and agronomic trait loci for crop improvement in chickpea

Chickpea (Cicer arietinum L.)-an important legume crop cultivated in arid and semiarid regions-has limited genetic diversity. Efforts are being undertaken to broaden its diversity by utilizing its wild relatives, which remain largely unexplored. Here, we present the Cicer super-pangenome based on the de novo genome assemblies of eight annual Cicer wild species. We identified 24,827 gene families, including 14,748 core, 2,958 softcore, 6,212 dispensable and 909 species-specific gene families. The dispensable genome was enriched for genes related to key agronomic traits. Structural variations between cultivated and wild genomes were used to construct a graph-based genome, revealing variations in genes affecting traits such as flowering time, vernalization and disease resistance. These variations will facilitate the transfer of valuable traits from wild Cicer species into elite chickpea varieties through marker-assisted selection or gene-editing. This study offers valuable insights into the genetic diversity and potential avenues for crop improvement in chickpea.

Celebrating Professor Rajeev K. Varshney's transformative research odyssey from genomics to the field on his induction as Fellow of the Royal Society

Professor Rajeev K. Varshney's transformative impact on crop genomics, genetics, and agriculture is the result of his passion, dedication, and unyielding commitment to harnessing the potential of genomics to address the most pressing challenges faced by the global agricultural community. Starting from a small town in India and reaching the global stage, Professor Varshney's academic and professional trajectory has inspired many scientists active in research today. His ground-breaking work, especially his effort to list orphan tropical crops to genomic resource-rich entities, has been transformative. Beyond his scientific achievements, Professor Varshney is recognized by his colleagues as an exemplary mentor, fostering the growth of future researchers, building institutional capacity, and strengthening scientific capability. His focus on translational genomics and strengthening seed system in developing countries for the improvement of agriculture has made a tangible impact on farmers' lives. His skills have been best utilized in roles at leading research centres where he has applied his expertise to deliver a new vision for crop improvement. These efforts have now been recognized by the Royal Society with the award of the Fellowship (FRS). As we mark this significant milestone in his career, we not only celebrate Professor Varshney's accomplishments but also his wider contributions that continue to transform the agricultural landscape.

Uncovering DNA methylation landscapes to decipher evolutionary footprints of phenotypic diversity in chickpea

Genetic diversity and environmental factors are long believed to be the dominant contributors to phenotypic diversity in crop plants. However, it has been recently established that, besides genetic variation, epigenetic variation, especially variation in DNA methylation, plays a significant role in determining phenotypic diversity in crop plants. Therefore, assessing DNA methylation diversity in crop plants becomes vital, especially in the case of crops like chickpea, which has a narrow genetic base. Thus, in the present study, we employed whole-genome bisulfite sequencing to assess DNA methylation diversity in wild and cultivated (desi and kabuli) chickpea. This revealed extensive DNA methylation diversity in both wild and cultivated chickpea. Interestingly, the methylation diversity was found to be significantly higher than genetic diversity, suggesting its potential role in providing vital phenotypic diversity for the evolution and domestication of the Cicer gene pool. The phylogeny based on DNA methylation variation also indicates a potential complementary role of DNA methylation variation in addition to DNA sequence variation in shaping chickpea evolution. Besides, the study also identified diverse epi-alleles of many previously known genes of agronomic importance. The Cicer MethVarMap database developed in this study enables researchers to readily visualize methylation variation within the genes and genomic regions of their interest (http://223.31.159.7/cicer/public/). Therefore, epigenetic variation like DNA methylation variation can potentially explain the paradox of high phenotypic diversity despite the narrow genetic base in chickpea and can potentially be employed for crop improvement.

Unique genetic architecture of prolificacy in 'Sikkim Primitive' maize unraveled through whole-genome resequencing-based DNA polymorphism

KEY MESSAGE

'Sikkim Primitive' maize landrace, unique for prolificacy (7-9 ears per plant) possesses unique genomic architecture in branching and inflorescence-related gene(s), and locus Zm00001eb365210 encoding glycosyltransferases was identified as the putative candidate gene underlying QTL (qProl-SP-8.05) for prolificacy. The genotype possesses immense usage in breeding high-yielding baby-corn genotypes. 'Sikkim Primitive' is a native landrace of North Eastern Himalayas, and is characterized by having 7-9 ears per plant compared to 1-2 ears in normal maize. Though 'Sikkim Primitive' was identified in the 1960s, it has not been characterized at a whole-genome scale. Here, we sequenced the entire genome of an inbred (MGUSP101) derived from 'Sikkim Primitive' along with three non-prolific (HKI1128, UMI1200, and HKI1105) and three prolific (CM150Q, CM151Q and HKI323) inbreds. A total of 942,417 SNPs, 24,160 insertions, and 27,600 deletions were identified in 'Sikkim Primitive'. The gene-specific functional mutations in 'Sikkim Primitive' were classified as 10,847 missense (54.36%), 402 non-sense (2.015%), and 8,705 silent (43.625%) mutations. The number of transitions and transversions specific to 'Sikkim Primitive' were 666,021 and 279,950, respectively. Among all base changes, (G to A) was the most frequent (215,772), while (C to G) was the rarest (22,520). Polygalacturonate 4-α-galacturonosyltransferase enzyme involved in pectin biosynthesis, cell-wall organization, nucleotide sugar, and amino-sugar metabolism was found to have unique alleles in 'Sikkim Primitive'. The analysis further revealed the Zm00001eb365210 gene encoding glycosyltransferases as the putative candidate underlying QTL (qProl-SP-8.05) for prolificacy in 'Sikkim Primitive'. High-impact nucleotide variations were found in ramosa3 (Zm00001eb327910) and zeaxanthin epoxidase1 (Zm00001eb081460) genes having a role in branching and inflorescence development in 'Sikkim Primitive'. The information generated unraveled the genetic architecture and identified key genes/alleles unique to the 'Sikkim Primitive' genome. This is the first report of whole-genome characterization of the 'Sikkim Primitive' landrace unique for its high prolificacy.

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.

Unraveling the genetics of heat tolerance in chickpea landraces (Cicer arietinum L.) using genome-wide association studies

Chickpea, being an important grain legume crop, is often confronted with the adverse effects of high temperatures at the reproductive stage of crop growth, drastically affecting yield and overall productivity. The current study deals with an extensive evaluation of chickpea genotypes, focusing on the traits associated with yield and their response to heat stress. Notably, we observed significant variations for these traits under both normal and high-temperature conditions, forming a robust basis for genetic research and breeding initiatives. Furthermore, the study revealed that yield-related traits exhibited high heritability, suggesting their potential suitability for marker-assisted selection. We carried out single-nucleotide polymorphism (SNP) genotyping using the genotyping-by-sequencing (GBS) method for a genome-wide association study (GWAS). Overall, 27 marker-trait associations (MTAs) linked to yield-related traits, among which we identified five common MTAs displaying pleiotropic effects after applying a stringent Bonferroni-corrected p-value threshold of <0.05 [-log10(p) > 4.95] using the BLINK (Bayesian-information and linkage-disequilibrium iteratively nested keyway) model. Through an in-depth in silico analysis of these markers against the CDC Frontier v1 reference genome, we discovered that the majority of the SNPs were located at or in proximity to gene-coding regions. We further explored candidate genes situated near these MTAs, shedding light on the molecular mechanisms governing heat stress tolerance and yield enhancement in chickpeas such as indole-3-acetic acid-amido synthetase GH3.1 with GH3 auxin-responsive promoter and pentatricopeptide repeat-containing protein, etc. The harvest index (HI) trait was associated with marker Ca3:37444451 encoding aspartic proteinase ortholog sequence of Oryza sativa subsp. japonica and Medicago truncatula, which is known for contributing to heat stress tolerance. These identified MTAs and associated candidate genes may serve as valuable assets for breeding programs dedicated to tailoring chickpea varieties resilient to heat stress and climate change.

Greater ecophysiological stress tolerance in the core environment than in extreme environments of wild chickpea (Cicer reticulatum)

Global climate change and land use change underlie a need to develop new crop breeding strategies, and crop wild relatives (CWR) have become an important potential source of new genetic material to improve breeding efforts. Many recent approaches assume adaptive trait variation increases towards the relative environmental extremes of a species range, potentially missing valuable trait variation in more moderate or typical climates. Here, we leveraged distinct genotypes of wild chickpea (Cicer reticulatum) that differ in their relative climates from moderate to more extreme and perform targeted assessments of drought and heat tolerance. We found significance variation in ecophysiological function and stress tolerance between genotypes but contrary to expectations and current paradigms, it was individuals from more moderate climates that exhibited greater capacity for stress tolerance than individuals from warmer and drier climates. These results indicate that wild germplasm collection efforts to identify adaptive variation should include the full range of environmental conditions and habitats instead of only environmental extremes, and that doing so may significantly enhance the success of breeding programs broadly.

A genomic variation map provides insights into peanut diversity in China and associations with 28 agronomic traits

Peanut (Arachis hypogaea L.) is an important allotetraploid oil and food legume crop. China is one of the world's largest peanut producers and consumers. However, genomic variations underlying the migration and divergence of peanuts in China remain unclear. Here we reported a genome-wide variation map based on the resequencing of 390 peanut accessions, suggesting that peanuts might have been introduced into southern and northern China separately, forming two cultivation centers. Selective sweep analysis highlights asymmetric selection between the two subgenomes during peanut improvement. A classical pedigree from South China offers a context for the examination of the impact of artificial selection on peanut genome. Genome-wide association studies identified 22,309 significant associations with 28 agronomic traits, including candidate genes for plant architecture and oil biosynthesis. Our findings shed light on peanut migration and diversity in China and provide valuable genomic resources for peanut improvement.

Whole genome resequencing and phenotyping of MAGIC population for high resolution mapping of drought tolerance in chickpea

Terminal drought is one of the major constraints to crop production in chickpea (Cicer arietinum L.). In order to map drought tolerance related traits at high resolution, we sequenced multi-parent advanced generation intercross (MAGIC) population using whole genome resequencing approach and phenotyped it under drought stress environments for two consecutive years (2013-14 and 2014-15). A total of 52.02 billion clean reads containing 4.67 TB clean data were generated on the 1136 MAGIC lines and eight parental lines. Alignment of clean data on to the reference genome enabled identification of a total, 932,172 of SNPs, 35,973 insertions, and 35,726 deletions among the parental lines. A high-density genetic map was constructed using 57,180 SNPs spanning a map distance of 1606.69 cM. Using compressed mixed linear model, genome-wide association study (GWAS) enabled us to identify 737 markers significantly associated with days to 50% flowering, days to maturity, plant height, 100 seed weight, biomass, and harvest index. In addition to the GWAS approach, an identity-by-descent (IBD)-based mixed model approach was used to map quantitative trait loci (QTLs). The IBD-based mixed model approach detected major QTLs that were comparable to those from the GWAS analysis as well as some exclusive QTLs with smaller effects. The candidate genes like FRIGIDA and CaTIFY4b can be used for enhancing drought tolerance in chickpea. The genomic resources, genetic map, marker-trait associations, and QTLs identified in the study are valuable resources for the chickpea community for developing climate resilient chickpeas.

Genome-wide association study, population structure, and genetic diversity of the tea plant in Guizhou Plateau

BACKGROUND

Guizhou Plateau, as one of the original centers of tea plant, has a profound multi-ethnic cultural heritage and abundant tea germplasm resources. However, the impact of indigenous community factors on the genetic diversity, population structure and geographical distribution of tea plant is still unclear.

RESULTS

Using the genotyping-by-sequencing (GBS) approach, we collected 415 tea plant accessions from the study sites, estimated genetic diversity, developed a core collection, and conducted a genome-wide association study (GWAS) based on 99,363 high-quality single-nucleotide polymorphisms (SNPs). A total of 415 tea accessions were clustered into six populations (GP01, GP02, GP03, GP04, GP05 and GP06), and the results showed that GP04 and GP05 had the highest and lowest genetic diversity (Pi = 0.214 and Pi = 0.145, respectively). Moreover, 136 tea accessions (33%) were selected to construct the core set that can represent the genetic diversity of the whole collection. By analyzing seven significant SNP markers associated with the traits such as the germination period of one bud and two leaves (OTL) and the germination period of one bud and three leaves (OtL), four candidate genes possibly related to OTL and OtL were identified.

CONCLUSIONS

This study revealed the impact of indigenous communities on the population structure of 415 tea accessions, indicating the importance of cultural practices for protection and utilization of tea plant genetic resources. Four potential candidate genes associated with the OTL and OtL of tea plant were also identified, which will facilitate genetic research, germplasm conservation, and breeding.

Genome-wide association studies for earliness, MYMIV resistance, and other associated traits in mungbean (Vigna radiata L. Wilczek) using genotyping by sequencing approach

Yellow mosaic disease (YMD) remains a major constraint in mungbean (Vigna radiata (L.)) production; while short-duration genotypes offer multiple crop cycles per year and help in escaping terminal heat stress, especially during summer cultivation. A comprehensive genotyping by sequencing (GBS)-based genome-wide association studies (GWAS) analysis was conducted using 132 diverse mungbean genotypes for traits like flowering time, YMD resistance, soil plant analysis development (SPAD) value, trichome density, and leaf area. The frequency distribution revealed a wide range of values for all the traits. GBS studies identified 31,953 high-quality single nucleotide polymorphism (SNPs) across all 11 mungbean chromosomes and were used for GWAS. Structure analysis revealed the presence of two genetically distinct populations based on ΔK. The linkage disequilibrium (LD) varied throughout the chromosomes and at r2 = 0.2, the mean LD decay was estimated as 39.59 kb. Two statistical models, mixed linear model (MLM) and Bayesian-information and Linkage-disequilibrium Iteratively Nested Keyway (BLINK) identified 44 shared SNPs linked with various candidate genes. Notable candidate genes identified include FPA for flowering time (VRADI10G01470; chr. 10), TIR-NBS-LRR for mungbean yellow mosaic India virus (MYMIV) resistance (VRADI09G06940; chr. 9), E3 ubiquitin-protein ligase RIE1 for SPAD value (VRADI07G28100; chr. 11), WRKY family transcription factor for leaf area (VRADI03G06560; chr. 3), and LOB domain-containing protein 21 for trichomes (VRADI06G04290; chr. 6). In-silico validation of candidate genes was done through digital gene expression analysis using Arabidopsis orthologous (compared with Vigna radiata genome). The findings provided valuable insight for marker-assisted breeding aiming for the development of YMD-resistant and early-maturing mungbean varieties.

The effect of environmental factors on the genetic differentiation of Cucurbita ficifolia populations based on whole-genome resequencing

BACKGROUND

Cucurbita ficifolia is one of the squash species most resistant to fungal pathogens, and has especially high resistance to melon Fusarium wilt. This species is therefore an important germplasm resource for the breeding of squash and melon cultivars.

RESULTS

Whole-genome resequencing of 223 individuals from 32 populations in Yunnan Province, the main cucurbit production area in China, was performed and 3,855,120 single-nucleotide polymorphisms (SNPs) and 1,361,000 InDels were obtained. SNP analysis suggested that levels of genetic diversity in C. ficifolia were high, but that different populations showed no significant genetic differentiation or geographical structure, and that individual C. ficifolia plants with fruit rinds of a similar color did not form independent clusters. A Mantel test conducted in combination with geographical distance and environmental factors suggested that genetic distance was not correlated with geographical distance, but had a significant correlation with environmental distance. Further associations between the genetic data and five environmental factors were analyzed using whole-genome association analysis. SNPs associated with each environmental factor were investigated and genes 250 kb upstream and downstream from associated SNPs were annotated. Overall, 15 marker-trait-associated SNPs (MTAs) and 293 genes under environmental selection were identified. The identified genes were involved in cell membrane lipid metabolism, macromolecular complexes, catalytic activity and other related aspects. Ecological niche modeling was used to simulate the distribution of C. ficifolia across time, from the present and into the future. We found that the area suitable for C. ficifolia changed with the changing climate in different periods.

CONCLUSIONS

Resequencing of the C. ficifolia accessions has allowed identification of genetic markers, such as SNPs and InDels. The SNPs identified in this study suggest that environmental factors mediated the formation of the population structure of C. ficifolia in China. These SNPs and Indels might also contribute to the variation in important pathways of genes for important agronomic traits such as yield, disease resistance and stress tolerance. Moreover, the genome resequencing data and the genetic markers identified from 223 accessions provide insight into the genetic variation of the C. ficifolia germplasm and will facilitate a broad range of genetic studies.

Dissecting chickpea genomic loci associated with the root penetration responsive traits in compacted soil

MAIN CONCLUSION

Soil compaction reduces root exploration in chickpea. We found genes related to root architectural traits in chickpea that can help understand and improve root growth in compacted soils. Soil compaction is a major concern for modern agriculture, as it constrains plant root growth, leading to reduced resource acquisition. Phenotypic variation for root system architecture (RSA) traits in compacted soils is present for various crops; however, studies on genetic associations with these traits are lacking. Therefore, we investigated RSA traits in different soil compaction levels and identified significant genomic associations in chickpea. We conducted a Genome-Wide Association Study (GWAS) of 210 chickpea accessions for 13 RSA traits under three bulk densities (BD) (1.1BD, 1.6BD, and 1.8BD). Soil compaction decreases root exploration by reducing 12 RSA traits, except average diameter (AD). Further, AD is negatively correlated with lateral root traits, and this correlation increases in 1.8BD, suggesting the negative effect of AD on lateral root traits. Interestingly, we identified probable candidate genes such as GLP3 and LRX for lateral root traits and CRF1-like for total length (TL) in 1.6BD soil. In heavy soil compaction, DGK2 is associated with lateral root traits. Reduction in laterals during soil compaction is mainly due to delayed seedling establishment, thus making lateral root number a critical trait. Interestingly, we also found a higher contribution of the  GxE component of the number of root tips (Tips) to the total variation than the other lateral traits. We also identified a pectin esterase, PPE8B, associated with Tips in high soil compaction and a significantly associated SNP with the relative change in Tips depicting a trade-off between Tips and AD. Identified genes and loci would help develop soil-compaction-resistant chickpea varieties.

Quantifying the impact of dynamic plant-environment interactions on metabolic regulation

A plant's genome encodes enzymes, transporters and many other proteins which constitute metabolism. Interactions of plants with their environment shape their growth, development and resilience towards adverse conditions. Although genome sequencing technologies and applications have experienced triumphantly rapid development during the last decades, enabling nowadays a fast and cheap sequencing of full genomes, prediction of metabolic phenotypes from genotype × environment interactions remains, at best, very incomplete. The main reasons are a lack of understanding of how different levels of molecular organisation depend on each other, and how they are constituted and expressed within a setup of growth conditions. Phenotypic plasticity, e.g., of the genetic model plant Arabidopsis thaliana, has provided important insights into plant-environment interactions and the resulting genotype x phenotype relationships. Here, we summarize previous and current findings about plant development in a changing environment and how this might be shaped and reflected in metabolism and its regulation. We identify current challenges in the study of plant development and metabolic regulation and provide an outlook of how methodological workflows might support the application of findings made in model systems to crops and their cultivation.

A comprehensive metabolomics and lipidomics atlas for the legumes common bean, chickpea, lentil and lupin

Legumes represent an important component of human and livestock diets; they are rich in macro- and micronutrients such as proteins, dietary fibers and polyunsaturated fatty acids. Whilst several health-promoting and anti-nutritional properties have been associated with grain content, in-depth metabolomics characterization of major legume species remains elusive. In this article, we used both gas chromatography-mass spectrometry (GC-MS) and liquid chromatography-mass spectrometry (LC-MS) to assess the metabolic diversity in the five legume species commonly grown in Europe, including common bean (Phaseolus vulgaris), chickpea (Cicer arietinum), lentil (Lens culinaris), white lupin (Lupinus albus) and pearl lupin (Lupinus mutabilis), at the tissue level. We were able to detect and quantify over 3400 metabolites covering major nutritional and anti-nutritional compounds. Specifically, the metabolomics atlas includes 224 derivatized metabolites, 2283 specialized metabolites and 923 lipids. The data generated here will serve the community as a basis for future integration to metabolomics-assisted crop breeding and facilitate metabolite-based genome-wide association studies to dissect the genetic and biochemical bases of metabolism in legume species.

Impact of Terminal Heat and Combined Heat-Drought Stress on Plant Growth, Yield, Grain Size, and Nutritional Quality in Chickpea (Cicer arietinum L.)

Chickpea is the third most consumed pulse and provides a kit of essential nutrients for an exponential population. High temperatures and drought stress are two major abiotic stresses that cause serious effects on chickpea growth and development. The comprehension of abiotic stresses' impact on chickpea productivity and nutritional quality will permit the selection of promising genotypes. The current study aimed to assess the impact of heat and drought stresses on plant growth, grain yield and its components, grain size, and nutritional quality in chickpea. For this purpose, 43 international chickpea genotypes were evaluated under normal, heat, and combined heat-drought stress conditions. The findings revealed a significant decrease of over 50% in plant height, biological yield, and seed yield under both stress conditions. Grain size and hundred-seed weight were the most heritable traits under normal, heat, and combined heat-drought stress. Proteins were accumulated under both stresses, evolving from 20.26% for normal conditions to 22.19% for heat stress and to 21.94% for combined heat-drought stress. For minerals, significant variation between treatments was observed for Mn, Mg, and Na. Our results also showed a significant impact of genotype and genotype-environment interaction factors only on K content. Using selection indices, 22 genotypes were identified as highly tolerant to the combined heat-drought stress, while eleven genotypes were heat-tolerant. Mineral profile analysis according to the contrasting tolerance clusters revealed decreased potassium content in susceptible genotypes, indicating genetic potential in the studied chickpea collection, ensuring tolerance to both stresses while maintaining good grain quality.

Genome-wide association study of fiber yield-related traits uncovers the novel genomic regions and candidate genes in Indian upland cotton (Gossypium hirsutum L.)

Upland cotton (Gossypium hirsutum L.) is a major fiber crop that is cultivated worldwide and has significant economic importance. India harbors the largest area for cotton cultivation, but its fiber yield is still compromised and ranks 22nd in terms of productivity. Genetic improvement of cotton fiber yield traits is one of the major goals of cotton breeding, but the understanding of the genetic architecture underlying cotton fiber yield traits remains limited and unclear. To better decipher the genetic variation associated with fiber yield traits, we conducted a comprehensive genome-wide association mapping study using 117 Indian cotton germplasm for six yield-related traits. To accomplish this, we generated 2,41,086 high-quality single nucleotide polymorphism (SNP) markers using genotyping-by-sequencing (GBS) methods. Population structure, PCA, kinship, and phylogenetic analyses divided the germplasm into two sub-populations, showing weak relatedness among the germplasms. Through association analysis, 205 SNPs and 134 QTLs were identified to be significantly associated with the six fiber yield traits. In total, 39 novel QTLs were identified in the current study, whereas 95 QTLs overlapped with existing public domain data in a comparative analysis. Eight QTLs, qGhBN_SCY_D6-1, qGhBN_SCY_D6-2, qGhBN_SCY_D6-3, qGhSI_LI_A5, qGhLI_SI_A13, qGhLI_SI_D9, qGhBW_SCY_A10, and qGhLP_BN_A8 were identified. Gene annotation of these fiber yield QTLs revealed 2,509 unique genes. These genes were predominantly enriched for different biological processes, such as plant cell wall synthesis, nutrient metabolism, and vegetative growth development in the gene ontology (GO) enrichment study. Furthermore, gene expression analysis using RNAseq data from 12 diverse cotton tissues identified 40 candidate genes (23 stable and 17 novel genes) to be transcriptionally active in different stages of fiber, ovule, and seed development. These findings have revealed a rich tapestry of genetic elements, including SNPs, QTLs, and candidate genes, and may have a high potential for improving fiber yield in future breeding programs for Indian cotton.

The Potential Genetic Effect for Yield and Foliar Disease Resistance in Faba Bean (Vicia faba L.) Assessed via Morphological and SCoT Markers

Faba bean is considered one of the most prominent grain legumes, with high protein content for human food consumption and livestock feed. The present study evaluated the nature of gene action and determined the genetic diversity among different populations of three crosses for resistance to foliar diseases at the molecular level. Analysis of variance exposed significant differences among the generations for all measured traits. Both dominance and additive gene effects were essential, but dominance genes, for the most part, exhibited greater effects than additive ones. This indicates an essential role for dominant genes alongside the additives one in inheriting such traits. The third cross (Marina × Giza 40) gave desired significant and positive (additive × additive) values for the number of pods/plant, seeds/plant, and seed yield/plant, in addition to desirable negative values for chocolate spot and rust characteristics. Furthermore, assessing the lines under study using seven SCoT primers disclosed three bands with recorded molecular weights of 260, 207, and 178 bp, generated by SCoT-1, SCoT-4, and SCoT-7 primers, respectively. These bands exist in the resistant parent (Marina), which could be attributed to the high-disease-resistance phenotypes, and they are absent in the sensitive parent (Giza 40) and other putative sensitive lines. Based on the molecular profiles and the genetic similarity between parents and the selected lines, the highest similarity value (0.91) was detected between Marina genotype and BC1, revealing a high foliar disease resistance. Meanwhile, Giza 40 (susceptible to foliar diseases) exhibited the maximum value (0.93) with F2. Additionally, cluster analysis based on genetic relationships was performed, and a high level of correlation between the results of PCR-based SCoT analysis and the foliar disease reactions was observed in the field. Consequently, this study concluded that SCoT markers created reliable banding profiles for evaluating genetic polymorphism among faba bean lines, which could be a foundation for developing an efficient breeding program.

Resequencing of global Lotus corniculatus accessions reveals population distribution and genetic loci, associated with cyanogenic glycosides accumulation and growth traits

BACKGROUND

Lotus corniculatus is a widely distributed perennial legume whose great adaptability to different environments and resistance to barrenness make it an excellent forage and ecological restoration plant. However, its molecular genetics and genomic relationships among populations are yet to be uncovered.

RESULT

Here we report on a genomic variation map from worldwide 272 L. corniculatus accessions by genome resequencing. Our analysis suggests that L. corniculatus accessions have high genetic diversity and could be further divided into three subgroups, with the genetic diversity centers were located in Transcaucasia. Several candidate genes and SNP site associated with CNglcs content and growth traits were identified by genome-wide associated study (GWAS). A non-synonymous in LjMTR was responsible for the decreased expression of CNglcs synthesis genes and LjZCD was verified to positively regulate CNglcs synthesis gene CYP79D3. The LjZCB and an SNP in LjZCA promoter were confirmed to be involved in plant growth.

CONCLUSION

This study provided a large number of genomic resources and described genetic relationship and population structure among different accessions. Moreover, we attempt to provide insights into the molecular studies and breeding of CNglcs and growth traits in L. corniculatus.

Plant genome resequencing and population genomics: Current status and future prospects

Advances in DNA sequencing technology have sparked a genomics revolution, driving breakthroughs in plant genetics and crop breeding. Recently, the focus has shifted from cataloging genetic diversity in plants to exploring their functional significance and delivering beneficial alleles for crop improvement. This transformation has been facilitated by the increasing adoption of whole-genome resequencing. In this review, we summarize the current progress of population-based genome resequencing studies and how these studies affect crop breeding. A total of 187 land plants from 163 countries have been resequenced, comprising 54 413 accessions. As part of resequencing efforts 367 traits have been surveyed and 86 genome-wide association studies have been conducted. Economically important crops, particularly cereals, vegetables, and legumes, have dominated the resequencing efforts, leaving a gap in 49 orders, including Lycopodiales, Liliales, Acorales, Austrobaileyales, and Commelinales. The resequenced germplasm is distributed across diverse geographic locations, providing a global perspective on plant genomics. We highlight genes that have been selected during domestication, or associated with agronomic traits, and form a repository of candidate genes for future research and application. Despite the opportunities for cross-species comparative genomics, many population genomic datasets are not accessible, impeding secondary analyses. We call for a more open and collaborative approach to population genomics that promotes data sharing and encourages contribution-based credit policy. The number of plant genome resequencing studies will continue to rise with the decreasing DNA sequencing costs, coupled with advances in analysis and computational technologies. This expansion, in terms of both scale and quality, holds promise for deeper insights into plant trait genetics and breeding design.

Insights into the Transcriptomics of Crop Wild Relatives to Unravel the Salinity Stress Adaptive Mechanisms

The narrow genomic diversity of modern cultivars is a major bottleneck for enhancing the crop's salinity stress tolerance. The close relatives of modern cultivated plants, crop wild relatives (CWRs), can be a promising and sustainable resource to broaden the diversity of crops. Advances in transcriptomic technologies have revealed the untapped genetic diversity of CWRs that represents a practical gene pool for improving the plant's adaptability to salt stress. Thus, the present study emphasizes the transcriptomics of CWRs for salinity stress tolerance. In this review, the impacts of salt stress on the plant's physiological processes and development are overviewed, and the transcription factors (TFs) regulation of salinity stress tolerance is investigated. In addition to the molecular regulation, a brief discussion on the phytomorphological adaptation of plants under saline environments is provided. The study further highlights the availability and use of transcriptomic resources of CWR and their contribution to pangenome construction. Moreover, the utilization of CWRs' genetic resources in the molecular breeding of crops for salinity stress tolerance is explored. Several studies have shown that cytoplasmic components such as calcium and kinases, and ion transporter genes such as Salt Overly Sensitive 1 (SOS1) and High-affinity Potassium Transporters (HKTs) are involved in the signaling of salt stress, and in mediating the distribution of excess Na⁺ ions within the plant cells. Recent comparative analyses of transcriptomic profiling through RNA sequencing (RNA-Seq) between the crops and their wild relatives have unraveled several TFs, stress-responsive genes, and regulatory proteins for generating salinity stress tolerance. This review specifies that the use of CWRs transcriptomics in combination with modern breeding experimental approaches such as genomic editing, de novo domestication, and speed breeding can accelerate the CWRs utilization in the breeding programs for enhancing the crop's adaptability to saline conditions. The transcriptomic approaches optimize the crop genomes with the accumulation of favorable alleles that will be indispensable for designing salt-resilient crops.

Genetic diversity, population structure, and genome-wide association study for the flowering trait in a diverse panel of 428 moth bean (Vigna aconitifolia) accessions using genotyping by sequencing

BACKGROUND

Moth bean (Vigna aconitifolia) is an underutilized, protein-rich legume that is grown in arid and semi-arid areas of south Asia and is highly resistant to abiotic stresses such as heat and drought. Despite its economic importance, the crop remains unexplored at the genomic level for genetic diversity and trait mapping studies. To date, there is no report of SNP marker discovery and association mapping of any trait in this crop. Therefore, this study aimed to dissect the genetic diversity, population structure and marker-trait association for the flowering trait in a diversity panel of 428 moth bean accessions using genotyping by sequencing (GBS) approach.

RESULTS

A total of 9078 high-quality single nucleotide polymorphisms (SNPs) were discovered by genotyping of 428 moth bean accessions. Model-based structure analysis and PCA grouped the moth bean accessions into two subpopulations. Cluster analysis revealed accessions belonging to the Northwestern region of India had higher variability than accessions from the other regions suggesting that this region represents its center of diversity. AMOVA revealed more variations within individuals (74%) and among the individuals (24%) than among the populations (2%). Marker-trait association analysis using seven multi-locus models including mrMLM, FASTmrEMMA FASTmrEMMA, ISIS EM-BLASSO, MLMM, BLINK and FarmCPU revealed 29 potential genomic regions for the trait days to 50% flowering, which were consistently detected in three or more models. Analysis of the allelic effect of the major genomic regions explaining phenotypic variance of more than 10% and those detected in at least 2 environments showed 4 genomic regions with significant phenotypic effect on this trait. Further, we also analyzed genetic relationships among the Vigna species using SNP markers. The genomic localization of moth bean SNPs on genomes of closely related Vigna species demonstrated that maximum numbers of SNPs were getting localized on Vigna mungo. This suggested that the moth bean is most closely related to V. mungo.

CONCLUSION

Our study shows that the north-western regions of India represent the center of diversity of the moth bean. Further, the study revealed flowering-related genomic regions/candidate genes which can be potentially exploited in breeding programs to develop early-maturity moth bean varieties.

Genetic analysis of global faba bean diversity, agronomic traits and selection signatures

KEY MESSAGE

We identified marker-trait associations for key faba bean agronomic traits and genomic signatures of selection within a global germplasm collection. Faba bean (Vicia faba L.) is a high-protein grain legume crop with great potential for sustainable protein production. However, little is known about the genetics underlying trait diversity. In this study, we used 21,345 high-quality SNP markers to genetically characterize 2678 faba bean genotypes. We performed genome-wide association studies of key agronomic traits using a seven-parent-MAGIC population and detected 238 significant marker-trait associations linked to 12 traits of agronomic importance. Sixty-five of these were stable across multiple environments. Using a non-redundant diversity panel of 685 accessions from 52 countries, we identified three subpopulations differentiated by geographical origin and 33 genomic regions subjected to strong diversifying selection between subpopulations. We found that SNP markers associated with the differentiation of northern and southern accessions explained a significant proportion of agronomic trait variance in the seven-parent-MAGIC population, suggesting that some of these traits were targets of selection during breeding. Our findings point to genomic regions associated with important agronomic traits and selection, facilitating faba bean genomics-based breeding.

The developmental dynamics in cool season legumes with focus on chickpea

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

Multi-locus genome-wide association study of chickpea reference set identifies genetic determinants of Pratylenchus thornei resistance

Pratylenchus thornei is an economically important species of root-lesion nematode adversely affecting chickpea (Cicer arietinum) yields globally. Integration of resistant crops in farming systems is recognised as the most effective and sustainable management strategy for plant-parasitic nematodes. However, breeding for P. thornei resistance in chickpea is limited by the lack of genetic diversity. We deployed a genome-wide association approach to identify genomic regions and candidate genes associated with P. thornei resistance in 285 genetically diverse chickpea accessions. Chickpea accessions were phenotyped for P. thornei resistance in replicated glasshouse experiments performed for two years (2018 and 2020). Whole genome sequencing data comprising 492,849 SNPs were used to implement six multi-locus GWAS models. Fourteen chickpea genotypes were found to be resistant to P. thornei. Of the six multi-locus GWAS methods deployed, FASTmrMLM was found to be the best performing model. In all, 24 significant quantitative trait nucleotides (QTNs) were identified, of which 13 QTNs were associated with lower nematode population density and 11 QTNs with higher nematode population density. These QTNs were distributed across all of the chickpea chromosomes, except chromosome 8. We identified, receptor-linked kinases (RLKs) on chromosomes 1, 4 and 6, GDSL-like Lipase/Acylhydrolase on chromosome 3, Aspartic proteinase-like and Thaumatin-like protein on chromosome 4, AT-hook DNA-binding and HSPRO2 on chromosome 6 as candidate genes for P. thornei resistance in the chickpea reference set. New sources of P. thornei resistant genotypes were identified that can be harnessed into breeding programs and putative candidate P. thornei resistant genes were identified that can be explored further to develop molecular markers and accelerate the incorporation of improved P. thornei resistance into elite chickpea cultivars.

Genome-wise association study identified genomic regions associated with drought tolerance in mungbean (Vigna radiata (L.) R. Wilczek)

A total of 282 mungbean accessions were resequenced to identify genome-wide variants and construct a highly precise variant map, and drought tolerance-related loci and superior alleles were identified by GWAS. Mungbean (Vigna radiata (L.) R. Wilczek) is an important food legume crop that is highly adapted to drought environments, but severe drought significantly curtails mungbean production. Here, we resequenced 282 mungbean accessions to identify genome-wide variants and constructed a highly precise map of mungbean variants. A genome-wide association study was performed to identify genomic regions for 14 drought tolerance-related traits in plants grown under stress and well-watered conditions over three years. One hundred forty-six SNPs associated with drought tolerance were detected, and 26 candidate loci associated with more than two traits were subsequently selected. Two hundred fifteen candidate genes were identified at these loci, including eleven transcription factor genes, seven protein kinase genes and other protein coding genes that may respond to drought stress. Furthermore, we identified superior alleles that were associated with drought tolerance and positively selected during the breeding process. These results provide valuable genomic resources for molecular breeding and will accelerate future efforts aimed at mungbean improvement.

Transcriptome-based variations effectively untangling the intraspecific relationships and selection signals in Xinyang Maojian tea population

As one of the world's top three popular non-alcoholic beverages, tea is economically and culturally valuable. Xinyang Maojian, this elegant green tea, is one of the top ten famous tea in China and has gained prominence for thousands of years. However, the cultivation history of Xinyang Maojian tea population and selection signals of differentiation from the other major variety Camellia sinensis var. assamica (CSA) remain unclear. We newly generated 94 Camellia sinensis (C. sinensis) transcriptomes including 59 samples in the Xinyang area and 35 samples collected from 13 other major tea planting provinces in China. Comparing the very low resolution of phylogeny inferred from 1785 low-copy nuclear genes with 94 C. sinensis samples, we successfully resolved the phylogeny of C. sinensis samples by 99,115 high-quality SNPs from the coding region. The sources of tea planted in the Xinyang area were extensive and complex. Specifically, Shihe District and Gushi County were the two earliest tea planting areas in Xinyang, reflecting a long history of tea planting. Furthermore, we identified numerous selection sweeps during the differentiation of CSA and CSS and these positive selection genes are involved in many aspects such as regulation of secondary metabolites synthesis, amino acid metabolism, photosynthesis, etc. Numerous specific selective sweeps of modern cultivars were annotated with functions in various different aspects, indicating the CSS and CSA populations possibly underwent independent specific domestication processes. Our study indicated that transcriptome-based SNP-calling is an efficient and cost-effective method in untangling intraspecific phylogenetic relationships. This study provides a significant understanding of the cultivation history of the famous Chinese tea Xinyang Maojian and unravels the genetic basis of physiological and ecological differences between the two major tea subspecies.

Improving phosphorus acquisition efficiency through modification of root growth responses to phosphate starvation in legumes

Phosphorus (P) is one of the essential macronutrients for plant growth and development, and it is an integral part of the major organic components, including nucleic acids, proteins and phospholipids. Although total P is abundant in most soils, a large amount of P is not easily absorbed by plants. Inorganic phosphate (Pi) is the plant-available P, which is generally immobile and of low availability in soils. Hence, Pi starvation is a major constraint limiting plant growth and productivity. Enhancing plant P efficiency can be achieved by improving P acquisition efficiency (PAE) through modification of morpho-physiological and biochemical alteration in root traits that enable greater acquisition of external Pi from soils. Major advances have been made to dissect the mechanisms underlying plant adaptation to P deficiency, especially for legumes, which are considered important dietary sources for humans and livestock. This review aims to describe how legume root growth responds to Pi starvation, such as changes in the growth of primary root, lateral roots, root hairs and cluster roots. In particular, it summarizes the various strategies of legumes to confront P deficiency by regulating root traits that contribute towards improving PAE. Within these complex responses, a large number of Pi starvation-induced (PSI) genes and regulators involved in the developmental and biochemical alteration of root traits are highlighted. The involvement of key functional genes and regulators in remodeling root traits provides new opportunities for developing legume varieties with maximum PAE needed for regenerative agriculture.

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.

Resequencing of 410 Sesame Accessions Identifies SINST1 as the Major Underlying Gene for Lignans Variation

Sesame is a promising oilseed crop that produces specific lignans of clinical importance. Hence, a molecular description of the regulatory mechanisms of lignan biosynthesis is essential for crop improvement. Here, we resequence 410 sesame accessions and identify 5.38 and 1.16 million SNPs (single nucleotide polymorphisms) and InDels, respectively. Population genomic analyses reveal that sesame has evolved a geographic pattern categorized into northern (NC), middle (MC), and southern (SC) groups, with potential origin in the southern region and subsequent introduction to the other regions. Selective sweeps analysis uncovers 120 and 75 significant selected genomic regions in MC and NC groups, respectively. By screening these genomic regions, we unveiled 184 common genes positively selected in these subpopulations for exploitation in sesame improvement. Genome-wide association study identifies 17 and 72 SNP loci for sesamin and sesamolin variation, respectively, and 11 candidate causative genes. The major pleiotropic SNP^C/A locus for lignans variation is located in the exon of the gene SiNST1. Further analyses revealed that this locus was positively selected in higher lignan content sesame accessions, and the "C" allele is favorable for a higher accumulation of lignans. Overexpression of SiNST1^C in sesame hairy roots significantly up-regulated the expression of SiMYB58, SiMYB209, SiMYB134, SiMYB276, and most of the monolignol biosynthetic genes. Consequently, the lignans content was significantly increased, and the lignin content was slightly increased. Our findings provide insights into lignans and lignin regulation in sesame and will facilitate molecular breeding of elite varieties and marker-traits association studies.

Exchanges of economic plants along the land silk road

BACKGROUD

The greatest contribution of the Silk Road is to communicate among different countries and nationalities, and promote two-way cultural exchanges between the East and the West. We now have clearer understanding about how material civilization and religious culture of Central Asia and West Asia spread eastward along the Land Silk Road. However, there is controversial about how crops migrate along the Land Silk Road.

RESULTS

We summarize archaeology, genetics, and genomics data to explore crop migration patterns. Of the 207 crops that were domesticated along the Land Silk Road, 19 for which genomic evidence was available were selected for discussion.

CONCLUSIONS

There were conflicting lines of evidence for the domestication of Tibetan barley, mustard, lettuce, buckwheat, and chickpea. The main reasons for the conflicting results may include incomplete early knowledge, record differences in different period, sample sizes, and data analysis techniques. There was strong evidence that Tibetan barley, barley, wheat, and jujube were introduced into China before the existence of the Land Silk Road; and mustard, lettuce, buckwheat, chickpea, alfalfa, walnut, cauliflower, grape, spinach, apple, cucumber, mulberry, and pea spread to China via trade and human migration along the Land Silk Road.

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.

The Key to the Future Lies in the Past: Insights from Grain Legume Domestication and Improvement Should Inform Future Breeding Strategies

Crop domestication is a co-evolutionary process that has rendered plants and animals significantly dependent on human interventions for survival and propagation. Grain legumes have played an important role in the development of Neolithic agriculture some 12,000 years ago. Despite being early companions of cereals in the origin and evolution of agriculture, the understanding of grain legume domestication has lagged behind that of cereals. Adapting plants for human use has resulted in distinct morpho-physiological changes between the wild ancestors and domesticates, and this distinction has been the focus of several studies aimed at understanding the domestication process and the genetic diversity bottlenecks created. Growing evidence from research on archeological remains, combined with genetic analysis and the geographical distribution of wild forms, has improved the resolution of the process of domestication, diversification and crop improvement. In this review, we summarize the significance of legume wild relatives as reservoirs of novel genetic variation for crop breeding programs. We describe key legume features, which evolved in response to anthropogenic activities. Here, we highlight how whole genome sequencing and incorporation of omics-level data have expanded our capacity to monitor the genetic changes accompanying these processes. Finally, we present our perspective on alternative routes centered on de novo domestication and re-domestication to impart significant agronomic advances of novel crops over existing commodities. A finely resolved domestication history of grain legumes will uncover future breeding targets to develop modern cultivars enriched with alleles that improve yield, quality and stress tolerance.

High-quality genome assembly and pan-genome studies facilitate genetic discovery in mung bean and its improvement

Mung bean is an economically important legume crop species that is used as a food, consumed as a vegetable, and used as an ingredient and even as a medicine. To explore the genomic diversity of mung bean, we assembled a high-quality reference genome (Vrad_JL7) that was ∼479.35 Mb in size, with a contig N50 length of 10.34 Mb. A total of 40,125 protein-coding genes were annotated, representing ∼96.9% of the genetic region. We also sequenced 217 accessions, mainly landraces and cultivars from China, and identified 2,229,343 high-quality single-nucleotide polymorphisms (SNPs). Population structure revealed that the Chinese accessions diverged into two groups and were distinct from non-Chinese lines. Genetic diversity analysis based on genomic data from 750 accessions in 23 countries supported the hypothesis that mung bean was first domesticated in south Asia and introduced to east Asia probably through the Silk Road. We constructed the first pan-genome of mung bean germplasm and assembled 287.73 Mb of non-reference sequences. Among the genes, 83.1% were core genes and 16.9% were variable. Presence/absence variation (PAV) events of nine genes involved in the regulation of the photoperiodic flowering pathway were identified as being under selection during the adaptation process to promote early flowering in the spring. Genome-wide association studies (GWASs) revealed 2,912 SNPs and 259 gene PAV events associated with 33 agronomic traits, including a SNP in the coding region of the SWEET10 homolog (jg24043) involved in crude starch content and a PAV event in a large fragment containing 11 genes for color-related traits. This high-quality reference genome and pan-genome will provide insights into mung bean breeding.

Unlocking the hidden variation from wild repository for accelerating genetic gain in legumes

The fluctuating climates, rising human population, and deteriorating arable lands necessitate sustainable crops to fulfil global food requirements. In the countryside, legumes with intriguing but enigmatic nitrogen-fixing abilities and thriving in harsh climatic conditions promise future food security. However, breaking the yield plateau and achieving higher genetic gain are the unsolved problems of legume improvement. Present study gives emphasis on 15 important legume crops, i.e., chickpea, pigeonpea, soybean, groundnut, lentil, common bean, faba bean, cowpea, lupin, pea, green gram, back gram, horse gram, moth bean, rice bean, and some forage legumes. We have given an overview of the world and India's area, production, and productivity trends for all legume crops from 1961 to 2020. Our review article investigates the importance of gene pools and wild relatives in broadening the genetic base of legumes through pre-breeding and alien gene introgression. We have also discussed the importance of integrating genomics, phenomics, speed breeding, genetic engineering and genome editing tools in legume improvement programmes. Overall, legume breeding may undergo a paradigm shift once genomics and conventional breeding are integrated in the near future.

The interaction between drought stress and nodule formation under multiple environments in chickpea

Environmental stresses, particularly drought, limit symbiotic nitrogen fixation in legumes, resulting in decreased yielding capacity. Drought is one of the most important constraints limiting yield potential in crops and it is the major abiotic stress that can cause more than 70% yield loss in chickpea. In this study, a total of two hundred four chickpea (Cicer arietinum L.) genotypes were selected to study the interaction between drought stress and nodule formation. This interaction was assessed by using morphological, yield and yield components. The field experiments were laid out in two locations (Terbol and Kfardan stations, Bekaa valley, Lebanon) using Alpha lattice design with two replications and two watering treatments (irrigation and rainfed) during 2016 and 2017 seasons. Parameters that were measured include days to 50% flowering (DFL), day to maturity (DM), plant height (PLH), nodule biomass (NB), nodule fresh weight (NFW), nodule dry weight (NDW), grain yield (GY), Biological yield (BY), 100 seed weight (100SW) and drought tolerance stress (DTS). The results indicated a significant variation between genotypes, environments and other morphological, yield and yield components traits. Drought stress reduced significantly the yield and the nodule's characteristics, biological and grain yield. The genotypes with the highest levels of drought tolerance, such as IG70399, IG8256, IG71832, IG70270, and IG70272, showed a minimal decrease in yield and nodule biomass. Nodule observations significantly and positively correlated with GY (0.36-0.38) under drought stress treatment. The correlation values for nodule characteristics with DFL and DM were higher under drought stress compared to irrigated conditions. This is a comparative study between drought stress and nodule formation traits associated with morphological, yield and yield components traits.

Genome-wide association mapping of seed oligosaccharides in chickpea

Chickpea (Cicer arietinum L.) is one of the major pulse crops, rich in protein, and widely consumed all over the world. Most legumes, including chickpeas, possess noticeable amounts of raffinose family oligosaccharides (RFOs) in their seeds. RFOs are seed oligosaccharides abundant in nature, which are non-digestible by humans and animals and cause flatulence and severe abdominal discomforts. So, this study aims to identify genetic factors associated with seed oligosaccharides in chickpea using the mini-core panel. We have quantified the RFOs (raffinose and stachyose), ciceritol, and sucrose contents in chickpea using high-performance liquid chromatography. A wide range of variations for the seed oligosaccharides was observed between the accessions: 0.16 to 15.13 mg g^-1 raffinose, 2.77 to 59.43 mg g^-1 stachyose, 4.36 to 90.65 mg g^-1 ciceritol, and 3.57 to 54.12 mg g^-1 for sucrose. Kabuli types showed desirable sugar profiles with high sucrose, whereas desi types had high concentrations RFOs. In total, 48 single nucleotide polymorphisms (SNPs) were identified for all the targeted sugar types, and nine genes (Ca_06204, Ca_04353, and Ca_20828: Phosphatidylinositol N-acetylglucosaminyltransferase; Ca_17399 and Ca_22050: Remorin proteins; Ca_11152: Protein-serine/threonine phosphatase; Ca_10185, Ca_14209, and Ca_27229: UDP-glucose dehydrogenase) were identified as potential candidate genes for sugar metabolism and transport in chickpea. The accessions with low RFOs and high sucrose contents may be utilized in breeding specialty chickpeas. The identified candidate genes could be exploited in marker-assisted breeding, genomic selection, and genetic engineering to improve the sugar profiles in legumes and other crop species.

Resequencing of 558 Chinese mungbean landraces identifies genetic loci associated with key agronomic traits

Mungbean is a warm-season annual food legume and plays important role in supplying food and nutritional security in many tropical countries. However, the genetic basis of its agronomic traits remains poorly understood. Therefore, we resequenced 558 Chinese mungbean landraces and produced a comprehensive map of mungbean genomic variation. We phenotyped all landraces in six different environments. Genome-wide association studies (GWAS) produced 110 signals significantly associated with nine agronomic traits, for which several candidate genes were identified. Overall, this study provides new insight into the genetic architecture of mungbean agronomic traits. Moreover, the genome-wide variations identified here should be valuable resources for future breeding studies of this important food legume.

Hybridization, missing wild ancestors and the domestication of cultivated diploid bananas

Hybridization and introgressions are important evolutionary forces in plants. They contribute to the domestication of many species, including understudied clonal crops. Here, we examine their role in the domestication of a clonal crop of outmost importance, banana (Musa ssp.). We used genome-wide SNPs generated for 154 diploid banana cultivars and 68 samples of the wild M. acuminata to estimate and geo-localize the contribution of the different subspecies of M. acuminata to cultivated banana. We further investigated the wild to domesticate transition in New Guinea, an important domestication center. We found high levels of admixture in many cultivars and confirmed the existence of unknown wild ancestors with unequal contributions to cultivated diploid. In New Guinea, cultivated accessions exhibited higher diversity than their direct wild ancestor, the latter recovering from a bottleneck. Introgressions, balancing selection and positive selection were identified as important mechanisms for banana domestication. Our results shed new lights on the radiation of M. acuminata subspecies and on how they shaped banana domestication. They point candidate regions of origin for two unknown ancestors and suggest another contributor in New Guinea. This work feed research on the evolution of clonal crops and has direct implications for conservation, collection, and 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.

Genome-wide association analysis of stress tolerance indices in an interspecific population of chickpea

Chickpea is a cool season crop that is highly vulnerable to abiotic stresses such as heat and drought. High temperature during early flowering and pod development stages significantly reduces the crop yield. The wild relatives of chickpeas can be potential donors for the introgression of heat and drought tolerance into cultivated chickpeas for crop improvement. Initially, 600 interspecific lines were derived from crosses between two elite cultivars, CDC Leader (kabuli chickpea) and CDC Consul (desi chickpea), and 20 accessions of Cicer reticulatum. The F₅ interspecific lines were tested for agronomic and seed quality traits including reaction to ascochyta blight disease under field conditions at two locations in 2018. A subset of 195 lines were selected based on resistance to ascochyta blight and acceptable seed quality. These lines were evaluated for their performance under suboptimal conditions at Lucky Lake (2019 and 2020) and Moose Jaw (2019), Saskatchewan, Canada, and Yuma, Arizona, United States (2019-2020). The lines were grown and evaluated at two seeding dates, normal (SD1) and late (SD2) seeding dates, at each location and year. The same lines were genotyped using Cicer60K Axiom® SNP chip. The population structure was determined based on 35,431 informative SNPs using fastStructure, and the interspecific lines were clustered at a k-value of 15. Significant marker-trait associations were identified for seed yield from SD1 and SD2 seeding dates, and stress tolerance indices (ATI, K₁STI, MP, SSPI, and TOL) using phenotypic values both from individual locations and combined analyses based on BLUP values. SNP marker Ca2_34600347 was significantly associated with yield from both the seeding dates. This and other SNP markers identified in this study may be useful for marker-assisted introgression of abiotic stress tolerance in chickpea.

Genome-Wide Association Analysis Reveals Trait-Linked Markers for Grain Nutrient and Agronomic Traits in Diverse Set of Chickpea Germplasm

Chickpea is an inexpensive source of protein, minerals, and vitamins to the poor people living in arid and semi-arid regions of Southern Asia and Sub-Saharan Africa. New chickpea cultivars with enhanced levels of protein, Fe and Zn content are a medium-term strategy for supplying essential nutrients for human health and reducing malnutrition. In the current study, a chickpea reference set of 280 accessions, including landraces, breeding lines, and advanced cultivars, was evaluated for grain protein, Fe, Zn content and agronomic traits over two seasons. Using a mid-density 5k SNP array, 4603 highly informative SNPs distributed across the chickpea genome were used for GWAS analysis. Population structure analysis revealed three subpopulations (K = 3). Linkage disequilibrium (LD) was extensive, and LD decay was relatively low. A total of 20 and 46 marker-trait associations (MTAs) were identified for grain nutrient and agronomic traits, respectively, using FarmCPU and BLINK models. Of which seven SNPs for grain protein, twelve for Fe, and one for Zn content were distributed on chromosomes 1, 4, 6, and 7. The marker S4_4477846 on chr4 was found to be co-associated with grain protein over seasons. The markers S1_11613376 and S1_2772537 co-associated with grain Fe content under NSII and pooled seasons and S7_9379786 marker under NSI and pooled seasons. The markers S4_31996956 co-associated with grain Fe and days to maturity. SNP annotation of associated markers were found to be related to gene functions of metal ion binding, transporters, protein kinases, transcription factors, and many more functions involved in plant metabolism along with Fe and protein homeostasis. The identified significant MTAs has potential use in marker-assisted selection for developing nutrient-rich chickpea cultivars after validation in the breeding populations.

Development of High Yielding Fusarium Wilt Resistant Cultivar by Pyramiding of “Genes” Through Marker-Assisted Backcrossing in Chickpea (Cicer arietinum L.)

Pusa 391, a mega desi chickpea variety with medium maturity duration is extensively cultivated in the Central Zone of India. Of late, this variety has become susceptible to Fusarium wilt (FW), which has drastic impact on its yield. Presence of variability in the wilt causing pathogen, Fusarium oxysporum f.sp. ciceri (foc) across geographical locations necessitates the role of pyramiding for FW resistance for different races (foc 1,2,3,4 and 5). Subsequently, the introgression lines developed in Pusa 391 genetic background were subjected to foreground selection using three SSR markers (GA16, TA 27 and TA 96) while 48 SSR markers uniformly distributed on all chromosomes, were used for background selection to observe the recovery of recurrent parent genome (RPG). BC1F1 lines with 75–85% RPG recovery were used to generate BC2F1. The plants that showed more than 90% RPG recovery in BC2F1 were used for generating BC3F1. The plants that showed more than 96% RPG recovery were selected and selfed to generate BC3F3. Multi-location evaluation of advanced introgression lines (BC2F3) in six locations for grain yield (kg/ha), days to fifty percent flowering, days to maturity, 100 seed weight and disease incidence was done. In case of disease incidence, the genotype IL1 (BGM 20211) was highly resistant to FW in Junagarh, Indore, New Delhi, Badnapur and moderately resistant at Sehore and Nandyal. GGE biplot analysis revealed that IL1(BGM20211) was the most stable genotype at Junagadh, Sehore and Nandyal. GGE biplot analysis revealed that IL1(BGM 20211) and IL4(BGM 20212) were the top performers in yield and highly stable across six environments and were nominated for Advanced Varietal Trials (AVT) of AICRP (All India Coordinated Research Project on Chickpea) in 2018–19. BGM20211 and BGM 20212 recorded 29 and 28.5% average yield gain over the recurrent parent Pusa 391, in the AVT-1 and AVT-2 over five environments. Thus, BGM20211 was identified for release and notified as Pusa Manav/Pusa Chickpea 20211 for Madhya Pradesh, Gujarat and Maharashtra, Southern Rajasthan, Bundhelkhand region of Uttar Pradesh states by the Central Sub-Committees on Crop Standards, Notification and Release of Varieties of Agricultural Crops, Ministry of Agriculture and Farmers Welfare, Government of India, for commercial cultivation in India (Gazette notification number S.O.500 (E) dt. 29-1-2021).Such pyramided lines give resilience to multiple races of fusarium wilt with added yield advantage.

Epigenomics as Potential Tools for Enhancing Magnitude of Breeding Approaches for Developing Climate Resilient Chickpea

Epigenomics has become a significant research interest at a time when rapid environmental changes are occurring. Epigenetic mechanisms mainly result from systems like DNA methylation, histone modification, and RNA interference. Epigenetic mechanisms are gaining importance in classical genetics, developmental biology, molecular biology, cancer biology, epidemiology, and evolution. Epigenetic mechanisms play important role in the action and interaction of plant genes during development, and also have an impact on classical plant breeding programs, inclusive of novel variation, single plant heritability, hybrid vigor, plant-environment interactions, stress tolerance, and performance stability. The epigenetics and epigenomics may be significant for crop adaptability and pliability to ambient alterations, directing to the creation of stout climate-resilient elegant crop cultivars. In this review, we have summarized recent progress made in understanding the epigenetic mechanisms in plant responses to biotic and abiotic stresses and have also tried to provide the ways for the efficient utilization of epigenomic mechanisms in developing climate-resilient crop cultivars, especially in chickpea, and other legume crops.

Unraveling Origin, History, Genetics, and Strategies for Accelerated Domestication and Diversification of Food Legumes

Domestication is a dynamic and ongoing process of transforming wild species into cultivated species by selecting desirable agricultural plant features to meet human needs such as taste, yield, storage, and cultivation practices. Human plant domestication began in the Fertile Crescent around 12,000 years ago and spread throughout the world, including China, Mesoamerica, the Andes and Near Oceania, Sub-Saharan Africa, and eastern North America. Indus valley civilizations have played a great role in the domestication of grain legumes. Crops, such as pigeon pea, black gram, green gram, lablab bean, moth bean, and horse gram, originated in the Indian subcontinent, and Neolithic archaeological records indicate that these crops were first domesticated by early civilizations in the region. The domestication and evolution of wild ancestors into today’s elite cultivars are important contributors to global food supply and agricultural crop improvement. In addition, food legumes contribute to food security by protecting human health and minimize climate change impacts. During the domestication process, legume crop species have undergone a severe genetic diversity loss, and only a very narrow range of variability is retained in the cultivars. Further reduction in genetic diversity occurred during seed dispersal and movement across the continents. In general, only a few traits, such as shattering resistance, seed dormancy loss, stem growth behavior, flowering–maturity period, and yield traits, have prominence in the domestication process across the species. Thus, identification and knowledge of domestication responsive loci were often useful in accelerating new species’ domestication. The genes and metabolic pathways responsible for the significant alterations that occurred as an outcome of domestication might aid in the quick domestication of novel crops. Further, recent advances in “omics” sciences, gene-editing technologies, and functional analysis will accelerate the domestication and crop improvement of new crop species without losing much genetic diversity. In this review, we have discussed about the origin, center of diversity, and seed movement of major food legumes, which will be useful in the exploration and utilization of genetic diversity in crop improvement. Further, we have discussed about the major genes/QTLs associated with the domestication syndrome in pulse crops and the future strategies to improve the food legume crops.

Genome-wide association study in quinoa reveals selection pattern typical for crops with a short breeding history

Quinoa germplasm preserves useful and substantial genetic variation, yet it remains untapped due to a lack of implementation of modern breeding tools. We have integrated field and sequence data to characterize a large diversity panel of quinoa. Whole-genome sequencing of 310 accessions revealed 2.9 million polymorphic high confidence SNP loci. Highland and Lowland quinoa were clustered into two main groups, with F_ST divergence of 0.36 and LD decay of 6.5 and 49.8 Kb, respectively. A genome-wide association study using multi-year phenotyping trials uncovered 600 SNPs stably associated with 17 traits. Two candidate genes are associated with thousand seed weight, and a resistance gene analog is associated with downy mildew resistance. We also identified pleiotropically acting loci for four agronomic traits important for adaptation. This work demonstrates the use of re-sequencing data of an orphan crop, which is partially domesticated to rapidly identify marker-trait association and provides the underpinning elements for genomics-enabled quinoa breeding.