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

The Kirkhouse Trust: Successes and Challenges in Twenty Years of Supporting Independent, Contemporary Grain Legume Breeding Projects in India and African Countries

This manuscript reviews two decades of projects funded by the Kirkhouse Trust (KT), a charity registered in the UK. KT was established to improve the productivity of legume crops important in African countries and in India. KT's requirements for support are: (1) the research must be conducted by national scientists in their home institution, either a publicly funded agricultural research institute or a university; (2) the projects need to include a molecular biology component, which to date has mostly comprised the use of molecular markers for the selection of one or more target traits in a crop improvement programme; (3) the projects funded are included in consortia, to foster the creation of scientific communities and the sharing of knowledge and breeding resources. This account relates to the key achievements and challenges, reflects on the lessons learned and outlines future research priorities.

The pattern of genetic variability in a core collection of 2,021 cowpea accessions

Cowpea is a highly drought-adapted leguminous crop with great promise for improving agricultural sustainability and food security. Here, we report analyses derived from array-based genotyping of 2,021 accessions constituting a core subset of the world's largest cowpea collection, held at the International Institute of Tropical Agriculture (IITA) in Ibadan, Nigeria. We used this dataset to examine genetic variation and population structure in worldwide cowpea. We confirm that the primary pattern of population structure is two geographically defined subpopulations originating in West and East Africa, respectively, and that population structure is associated with shifts in phenotypic distribution. Furthermore, we establish the cowpea core collection as a resource for genome-wide association studies by mapping the genetic basis of several phenotypes, with a focus on seed coat pigmentation patterning and color. We anticipate that the genotyped IITA Cowpea Core Collection will serve as a powerful tool for mapping complex traits, facilitating the acceleration of breeding programs to enhance the resilience of this crop in the face of rapid global climate change.

Genetic progress in cowpea [Vigna unguiculata (L.) Walp.] stemming from breeding modernization efforts at the International Institute of Tropical Agriculture

Genetic gain has been proposed as a quantifiable key performance indicator that can be used to monitor breeding programs' effectiveness. The cowpea breeding program at the International Institute of Tropical Agriculture (IITA) has developed and released improved varieties in 70 countries globally. To quantify the genetic changes to grain yield and related traits, we exploited IITA cowpea historical multi-environment trials (METs) advanced yield trial (AYT) data from 2010 to 2022. The genetic gain assessment targeted short duration (SD), medium duration (MD), and late duration (LD) breeding pipelines. A linear mixed model was used to calculate the best linear unbiased estimates (BLUE). Regressed BLUE of grain yield by year of genotype origin depicted realized genetic gain of 22.75 kg/ha/year (2.65%), 7.91 kg/ha/year (0.85%), and 22.82 kg/ha/year (2.51%) for SD, MD, and LD, respectively. No significant gain was realized in 100-seed weight (Hsdwt). We predicted, based on 2022 MET data, that recycling the best genotypes at AYT stage would result in grain yield gain of 37.28 kg/ha/year (SD), 28.00 kg/ha/year (MD), and 34.85 kg/ha/year (LD), and Hsdwt gain of 0.48 g/year (SD), 0.68 g/year (MD), and 0.55 g/year (LD). These results demonstrated a positive genetic gain trend for cowpea, indicating that a yield plateau has not yet been reached and that accelerated gain is expected with the recent integration of genomics in the breeding program. Advances in genomics include the development of the reference genome, genotyping platforms, quantitative trait loci mapping of key traits, and active implementation of molecular breeding.

Genetic Dissection of Diverse Seed Coat Patterns in Cowpea through a Comprehensive GWAS Approach

This study investigates the genetic determinants of seed coat color and pattern variations in cowpea (Vigna unguiculata), employing a genome-wide association approach. Analyzing a mapping panel of 296 cowpea varieties with 110,000 single nucleotide polymorphisms (SNPs), we focused on eight unique coat patterns: (1) Red and (2) Cream seed; (3) White and (4) Brown/Tan seed coat; (5) Pink, (6) Black, (7) Browneye and (8) Red/Brown Holstein. Across six GWAS models (GLM, SRM, MLM, MLMM, FarmCPU from GAPIT3, and TASSEL5), 13 significant SNP markers were identified and led to the discovery of 23 candidate genes. Among these, four specific genes may play a direct role in determining seed coat pigment. These findings lay a foundational basis for future breeding programs aimed at creating cowpea varieties aligned with consumer preferences and market requirements.

A view of the pan-genome of domesticated Cowpea (Vigna unguiculata [L.] Walp.)

Cowpea, Vigna unguiculata L. Walp., is a diploid warm-season legume of critical importance as both food and fodder in sub-Saharan Africa. This species is also grown in Northern Africa, Europe, Latin America, North America, and East to Southeast Asia. To capture the genomic diversity of domesticates of this important legume, de novo genome assemblies were produced for representatives of six subpopulations of cultivated cowpea identified previously from genotyping of several hundred diverse accessions. In the most complete assembly (IT97K-499-35), 26,026 core and 4963 noncore genes were identified, with 35,436 pan genes when considering all seven accessions. GO terms associated with response to stress and defense response were highly enriched among the noncore genes, while core genes were enriched in terms related to transcription factor activity, and transport and metabolic processes. Over 5 million single nucleotide polymorphisms (SNPs) relative to each assembly and over 40 structural variants >1 Mb in size were identified by comparing genomes. Vu10 was the chromosome with the highest frequency of SNPs, and Vu04 had the most structural variants. Noncore genes harbor a larger proportion of potentially disruptive variants than core genes, including missense, stop gain, and frameshift mutations; this suggests that noncore genes substantially contribute to diversity within domesticated cowpea.

Development and Assessment of SNP Genotyping Arrays for Citrus and Its Close Relatives

Rapid advancements in technologies provide various tools to analyze fruit crop genomes to better understand genetic diversity and relationships and aid in breeding. Genome-wide single nucleotide polymorphism (SNP) genotyping arrays offer highly multiplexed assays at a relatively low cost per data point. We report the development and validation of 1.4M SNP Axiom® Citrus HD Genotyping Array (Citrus 15AX 1 and Citrus 15AX 2) and 58K SNP Axiom® Citrus Genotyping Arrays for Citrus and close relatives. SNPs represented were chosen from a citrus variant discovery panel consisting of 41 diverse whole-genome re-sequenced accessions of Citrus and close relatives, including eight progenitor citrus species. SNPs chosen mainly target putative genic regions of the genome and are accurately called in both Citrus and its closely related genera while providing good coverage of the nuclear and chloroplast genomes. Reproducibility of the arrays was nearly 100%, with a large majority of the SNPs classified as the most stringent class of markers, "PolyHighResolution" (PHR) polymorphisms. Concordance between SNP calls in sequence data and array data average 98%. Phylogenies generated with array data were similar to those with comparable sequence data and little affected by 3 to 5% genotyping error. Both arrays are publicly available.

Quantitative trait loci and genomic prediction for grain sugar and mineral concentrations of cowpea [Vigna unguiculata (L.) Walp.]

Development of high yielding cowpea varieties coupled with good taste and rich in essential minerals can promote consumption and thus nutrition and profitability. The sweet taste of cowpea grain is determined by its sugar content, which comprises mainly sucrose and galacto-oligosaccharides (GOS) including raffinose and stachyose. However, GOS are indigestible and their fermentation in the colon can produce excess intestinal gas, causing undesirable bloating and flatulence. In this study, we aimed to examine variation in grain sugar and mineral concentrations, then map quantitative trait loci (QTLs) and estimate genomic-prediction (GP) accuracies for possible application in breeding. Grain samples were collected from a multi-parent advanced generation intercross (MAGIC) population grown in California during 2016-2017. Grain sugars were assayed using high-performance liquid chromatography. Grain minerals were determined by inductively coupled plasma-optical emission spectrometry and combustion. Considerable variation was observed for sucrose (0.6-6.9%) and stachyose (2.3-8.4%). Major QTLs for sucrose (QSuc.vu-1.1), stachyose (QSta.vu-7.1), copper (QCu.vu-1.1) and manganese (QMn.vu-5.1) were identified. Allelic effects of major sugar QTLs were validated using the MAGIC grain samples grown in West Africa in 2017. GP accuracies for minerals were moderate (0.4-0.58). These findings help guide future breeding efforts to develop mineral-rich cowpea varieties with desirable sugar content.

Understanding specific gender dynamics in the cowpea value chain for key traits to inform cowpea breeding programs in Malawi, Mozambique and Tanzania

Introduction

Cowpea is an important food and nutrition security crop in Malawi, Mozambique and Tanzania and it is mainly produced by women farmers mainly on a subsistence scale. The majority of these farmers use local varieties despite the availability of improved varieties in the region. Low acceptability and adoption of improved varieties have also hampered cowpea breeding efforts. The low adoption, especially among women farmers, has been attributed to the failure by breeding programs to involve farmers in the process of designing and developing improved varieties with a view to meeting their priorities and preferences. Despite women constituting the majority of cowpea farmers in these countries, no comprehensive gender analysis on cowpea value chain had been instituted to understand the traits that are gender and youth responsive and how to incorporate them in the product profiling so that the developed varieties benefit men, women and youth. The main objective of the gender study was, therefore, to identify preferred traits by different gender groups within the whole cowpea value chain to inform cowpea breeding programs in the three countries.

Methods

The study employed quantitative and qualitative methods to elucidate preferences, including value chain mapping, a quantitative survey of farmers, focus group discussions and key informant interviews targeting farmers/consumers, traders, policymakers and processors.

Results and discussion

Results showed that the top-ranking traits in order of importance across the countries and gender were; (1) high grain yield, (2) good grain taste, (3) early maturity time, (4) large grain size, (5) good leaf taste and (6) short cooking time. It was further noted that different gender groups preferred almost similar traits though minor variations were noted in terms of prioritization of these traits. These results have had two major influences on our cowpea breeding program: firstly, the breeding program changed the way it prioritizes traits to include ones that reflect the needs of men, women and the youth in the cowpea value chain. Secondly, our breeding objectives are closely aligned to gender differences in the target population of farmers and other users, by incorporating key priority traits that address the needs of both men and women, including the youth. That is to say, product targets and specific product profiles are more gender sensitive. Since the breeding work is ongoing, the expectation is that the development of improved varieties resulting from this gender sensitive process will translate into higher adoption levels of these varieties (compared with previous releases), which might have ripple effects on food, nutrition and income security in the region.

Urine and Dried Blood Spots From Children and Pregnant Women Reveal Phytochemicals, Amino Acids, and Carnitine Metabolites as Cowpea Consumption Biomarkers

SCOPE

Legumes consumption has been proven to promote health across the lifespan; cowpeas have demonstrated efficacy in combating childhood malnutrition and growth faltering, with an estimated malnutrition prevalence of 35.6% of children in Ghana. This cowpea feeding study aimed to identify a suite of metabolic consumption biomarkers in children and adults.

METHODS AND RESULTS

Urine and dried blood spots (DBS) from 24 children (9-21 months) and 21 pregnant women (>18 years) in Northern Ghana are collected before and after dose-escalated consumption of four cowpea varieties for 15 days. Untargeted metabolomics identified significant increases in amino acids, phytochemicals, and lipids. The carnitine metabolism pathway is represented by 137 urine and 43 DBS metabolites, with significant changes to tiglylcarnitine and acetylcarnitine. Additional noteworthy candidate biomarkers are mansouramycin C, N-acetylalliin, proline betaine, N2, N5-diacetylornithine, S-methylcysteine, S-methylcysteine sulfoxide, and cis-urocanate. S-methylcysteine and S-methylcysteine sulfoxide are targeted and quantified in urine.

CONCLUSION

This feeding study for cowpea biomarkers supports the utility of a suite of key metabolites classified as amino acids, lipids, and phytochemicals for dietary legume and cowpea-specific food exposures of global health importance.

A Mid-Density Single-Nucleotide Polymorphism Panel for Molecular Applications in Cowpea (Vigna unguiculata (L.) Walp)

Molecular markers are increasingly being deployed to accelerate genetic gain in crop plants. The objective of this study was to assess the potential of a mid-density genotyping panel for molecular applications in cowpea breeding. A core set of 2,602 targeted diversity array technology (DArTag) single-nucleotide polymorphisms (SNPs) was designed from an existing 51,128 Cowpea iSelect Consortium Array. The panel's usefulness was assessed using 376 genotypes from different populations of known genetic backgrounds. The panel was informative, with over 78% of SNPs exceeding a minor allele frequency of 0.20. The panel decoded three stratifications in the constituted population, as was expected. Linkage disequilibrium (LD) decay was correctly depicted as slower in a biparental subset than in other populations. A known flower and seed coat color gene region was located on chromosome Vu07, suggesting that the mid-density panel may be used to hypothesize genomic regions underlying target traits in cowpea. Unexpected heterozygosity was detected in some lines and highly among F1 progenies, divulging the panel's potential application in germplasm purity and hybridity verification. The study unveils the potential of an excellent genomic resource that can be tapped to enhance the development of improved cowpea cultivars.

Construction of a high-density genetic map for yardlong bean and identification of ANT1 as a regulator of anthocyanin biosynthesis

Because its long, tender pods supply essential proteins, vitamins, and fibers to humans, yardlong bean (Vigna unguiculata ssp. sesquipedalis) is a commonly consumed vegetable, especially in Southeast Asia. To provide insights into the genetic bases of key agricultural traits in yardlong bean, we here created a high-density bin-map with 2084 bin markers using 514 227 SNPs from a recombinant-inbred line (RIL) population. Quantitative trait loci (QTL) mapping was carried out to identify loci associated with anthocyanin content (ANT), vitamin E content (VE), total soluble protein content (TSP), pod length (PL), hundred-seed weight (HSW), seed length and width (SL and SW, respectively), and seed coat color (SCC). In total, 20 related QTLs were isolated, explaining 7.58-56.03% of the phenotypic variation. Of these, five major QTLs (qANT5, qTSP11, qVE7, qPL3, and qSCC9) were detected in 2020, 2021, and the combined environment, explaining 11.96-56.03% of the phenotypic variation. VuANT1 was identified as a causal gene for the QTL qANT5, which regulated anthocyanin content; VuANT1 was highly expressed in immature purple pods but barely detectable in white pods. VuANT1 overexpression in tobacco leaves and yardlong bean hairy roots resulted in purple coloration as a result of anthocyanin accumulation. These findings suggested that VuANT1 was a key regulator of anthocyanin accumulation in yardlong bean. Our results lay a firm foundation for target agricultural trait improvement and clarification of the genetic mechanisms underlying agricultural traits in yardlong bean.

Gender trait preferences among smallholder cowpea farmers in northern Ghana: lessons from a case study

Introduction

This case study reports on how a gender responsive breeding program contributes to meeting the trait preference of men and women for improved cowpea varieties in northern Ghana.

Methods

Fifty-eight early-maturing, medium-maturing and dual-purpose cowpea lines were planted at the CSIR-SARI research fields and women and men farmers invited for participatory plant breeding (PPB) in 2016. Selected lines from the PPB were further evaluated in 2017 using participatory varietal selection (PVS) in 5 districts in northern Ghana. In addition, 20 focus group discussions (FGDs) were held in 2018 in 10 randomly selected communities with 260 participants (130 women and 130 men) across the districts where the PVS had been held previously.

Results and discussion

The study finds drought tolerance, short cooking time and pest resistance to be the most preferred cowpea traits among both men and women. The study also finds that gender differences exist in trait preference, especially for traits such as seed coat color, earliness, pod above canopy and indeterminate growth habit. As breeding programs focus on improving genetic gains for tolerance of biotic and abiotic stresses, equal attention must be given to breeding for traits desired by women.

Cowpea Constraints and Breeding in Europe

Cowpea (Vigna unguiculata (L.) Walp.) is a legume with a constant rate of cultivation in Southern European countries. Consumer demand for cowpea worldwide is rising due to its nutritional content, while Europe is constantly attempting to reduce the deficit in the production of pulses and invest in new, healthy food market products. Although the climatic conditions that prevail in Europe are not so harsh in terms of heat and drought as in the tropical climates where cowpea is mainly cultivated, cowpea confronts with a plethora of abiotic and biotic stresses and yield-limiting factors in Southern European countries. In this paper, we summarize the main constraints for cowpea cultivation in Europe and the breeding methods that have been or can be used. A special mention is made of the availability plant genetic resources (PGRs) and their potential for breeding purposes, aiming to promote more sustainable cropping systems as climatic shifts become more frequent and fiercer, and environmental degradation expands worldwide.

A large-effect QTL introgressed from ricebean imparts resistance to Mungbean yellow mosaic India virus in blackgram (Vigna mungo (L.) Hepper)

Here, we report identification of a large effect QTL conferring Mungbean yellow mosaic India virus resistance introgressed from ricebean in blackgram variety Mash114. The tightly linked KASP markers would assist in marker-assisted-transfer of this region into Vigna species infected by MYMIV. Until recently, precise location of genes and marker-assisted selection was long thought in legumes such as blackgram due to lack of dense molecular maps. However, advances in next-generation sequencing based on high-throughput genotyping technologies such as QTL-seq have revolutionized trait mapping in marker-orphan crops. Using QTL-seq approach, we have identified a large-effect QTL for resistance to Mungbean yellow mosaic India virus (MYMIV) in blackgram variety Mash114. MYMIV is devastating disease responsible for huge yield losses in blackgram, greengram and other legumes. Mash114 showed consistent and high level of resistance to MYMIV since last nine years. Whole genome re-sequencing of MYMIV-resistant and susceptible bulks derived from RILs of cross KUG253 X Mash114 identified a large-effect QTL (qMYMIV6.1.1) spanning 3.4 Mb on chromosome 6 explaining 70% of total phenotypic variation. This region was further identified as an inter-specific introgression from ricebean. Linkage mapping using KASP markers developed from potent candidate genes involved in virus resistance identified the 500 kb genomic region equaling 1.9 cM on genetic map linked with MYMIV. The three KASP markers closely associated with MYMIV originated from serine threonine kinase, UBE2D2 and BAK1/BRI1-ASSOCIATED RECEPTOR KINASE genes. These KASPs can be used for marker-assisted transfer of introgressed segment into suitable backgrounds of Vigna species.

OMICS in Fodder Crops: Applications, Challenges, and Prospects

Biomass yield and quality are the primary targets in forage crop improvement programs worldwide. Low-quality fodder reduces the quality of dairy products and affects cattle's health. In multipurpose crops, such as maize, sorghum, cowpea, alfalfa, and oat, a plethora of morphological and biochemical/nutritional quality studies have been conducted. However, the overall growth in fodder quality improvement is not on par with cereals or major food crops. The use of advanced technologies, such as multi-omics, has increased crop improvement programs manyfold. Traits such as stay-green, the number of tillers per plant, total biomass, and tolerance to biotic and/or abiotic stresses can be targeted in fodder crop improvement programs. Omic technologies, namely genomics, transcriptomics, proteomics, metabolomics, and phenomics, provide an efficient way to develop better cultivars. There is an abundance of scope for fodder quality improvement by improving the forage nutrition quality, edible quality, and digestibility. The present review includes a brief description of the established omics technologies for five major fodder crops, i.e., sorghum, cowpea, maize, oats, and alfalfa. Additionally, current improvements and future perspectives have been highlighted.

A Genome-Wide Scan Divulges Key Loci Involved in Resistance to Aphids (Aphis craccivora) in Cowpea (Vigna unguiculata)

Cowpea aphids (Aphis craccivora Koch) double as a direct damaging pest and a virus vector to cowpea, threatening the economic yield of the crop. Given the multiple ecotypes, different alleles have been implicated in aphid resistance, necessitating the identification of key genes involved. The present study implemented a genome-wide scan using 365 cowpea mini-core accessions to decipher loci involved in resistance to aphid ecotype from Kano, Nigeria. Accessions were artificially inoculated with A. craccivora in insect-proof cages and damage severity assessed at 21 days after infestation. Significant phenotypic differences based on aphid damage severity were registered among the accessions. Skewed phenotypic distributions were depicted in the population, suggesting the involvement of major genes in the control of resistance. A genome-wide scan identified three major regions on chromosomes Vu10, Vu08 and Vu02, and two minor ones on chromosomes Vu01 and Vu06, that were significantly associated with aphid resistance. These regions harbored several genes, out of which, five viz Vigun01g233100.1, Vigun02g088900.1, Vigun06g224900.1, Vigun08g030200.1 and Vigun10g031100.1 were the most proximal to the peak single nucleotide polymorphisms (SNPs) positions. These genes are expressed under stress signaling, mechanical wounding and insect feeding. The uncovered loci contribute towards establishing a marker-assisted breeding platform and building durable resistance against aphids in cowpea.

Variations of Nutrient and Antinutrient Components of Bambara Groundnut (Vigna subterranea (L.) Verdc.) Seeds

Bambara groundnut (BGN) fits the bill when it comes to an acceptable level of nutrient and mineral composition. BGN is a balanced food that can help eradicate food and nutritional insecurity if it is incorporated into the major food system. However, there is a large degree of variation in nutrient composition and antinutritional factors among BGN accessions. Here, we show the degree of variability of nutrient and antinutrient components such as percentage ash, moisture, protein, fat, tryptophan, tannin, and phytate contents in seeds of 95 accessions of BGN. Data were subjected to analysis of variance (ANOVA), followed by correlation and principal component analysis. Clustering was done to show the relatedness between the accessions in response to the various traits. A high level of heterogeneity was observed among the accessions for the various traits studied. PC1 and PC2 show 41.2% of the total observed variations. Cluster analysis grouped accessions into four main clusters. This study was able to confirm the high level of diversity in the components of nutrients and antinutrients previously reported in BGN. The results of this study are expected to aid in identifying parent lines for improved breeding programs.

Development and optimization of NIRS prediction models for simultaneous multi-trait assessment in diverse cowpea germplasm

Cowpea (Vigna unguiculata (L.) Walp.) is one such legume that can facilitate achieving sustainable nutrition and climate change goals. Assessing nutritional traits conventionally can be laborious and time-consuming. NIRS is a technique used to rapidly determine biochemical parameters for large germplasm. NIRS prediction models were developed to assess protein, starch, TDF, phenols, and phytic acid based on MPLS regression. Higher RSQ_external values such as 0.903, 0.997, 0.901, 0.706, and 0.955 were obtained for protein, starch, TDF, phenols, and phytic acid respectively. Models for all the traits displayed RPD values of >2.5 except phenols and low SEP indicating the excellent prediction of models. For all the traits worked, p-value ≥ 0.05 implied the accuracy and reliability score >0.8 (except phenol) ensured the applicability of the models. These prediction models will facilitate high throughput screening of large cowpea germplasm in a non-destructive way and the selection of desirable chemotypes in any genetic background with huge application in cowpea crop improvement programs across the world.

Cowpea speed breeding using regulated growth chamber conditions and seeds of oven-dried immature pods potentially accommodates eight generations per year

BACKGROUND

Cowpea is a dryland crop with potential to improve food security in sub-Saharan Africa, where it is mostly produced and consumed. Contemporary plant improvement technologies, including genome editing, marker-assisted selection, and optimized transformation protocols, are being deployed to improve cowpea characteristics. Integrating speed breeding with these technologies would accelerate genetic gain in cowpea breeding. There are established speed breeding protocols for other important legumes, such as soybean, peanut, and chickpea, but none has been previously reported for cowpea.

RESULTS

With the aid of regulated growth conditions in two different chamber types, as well as the cultivation of new plant generations from seeds of oven-dried immature pods, we developed and validated, for the first time, an efficient speed breeding protocol that accommodates approximately seven to eight breeding generations per year for 3 cowpea genotypes. The 3 cowpea genotypes were evaluated under controlled growth conditions in light-emitting diode and metal halide lamp chambers to determine the effect of CO₂ supplementation on flowering and maturation durations, optimum conditions for plant growth, cross pollination, and pod development. Elevated CO₂ concentration had no influence on either flowering time or pod development. Adequate temperature, relative humidity and light intensity improved plant development and the rate of successful hand pollination, and  cultivating seeds of 11-day-old immature pods oven-dried at 39 °C for 2 days resulted in at least a 62% reduction in the time between pollination and sowing of the next plant generation. The plants cultivated from seeds of the oven-dried immature pods showed no defect at any stage of development.

CONCLUSIONS

Using the speed breeding protocol developed in this study, cowpea breeding cycles can be increased from the traditional one cycle per year in the field to as many as 8 generations per year in regulated growth chamber conditions. This protocol has no special technical requirements; hence, it can be implemented in any standard growth chamber. This would fast-track development, testing, validation, and utilization of improved cowpea cultivars.

Construction of a single nucleotide polymorphism linkage map and identification of quantitative trait loci controlling heat tolerance in cowpea, Vigna unguiculata (L.) Walp

Plant tolerance to heat or high temperature is crucial to crop production, especially in the situation of elevated temperature resulting from global climate change. Cowpea, Vigna unguiculata (L.) Walp., is an internationally important legume food crop and an excellent pool of genes for numerous traits resilient to environmental extremes, particularly heat and drought. Here, we report a single nucleotide polymorphism (SNP) genetic map for cowpea and identification of the loci controlling the heat tolerance in the species. The SNP map consists of 531 bins containing 4,154 SNPs grouped into 11 linkage groups, and collectively spans 1,084.7 cM, thus having a density of one SNP in 0.26 cM or 149 kb. The 11 linkage groups of the map were aligned to the 11 cowpea chromosomes. Quantitative trait locus (QTL) mapping identified nine QTLs responsible for the cowpea heat tolerance on seven of the 11 chromosomes, with each QTL explaining 6.5-21.8% of heat tolerance phenotypic variation. Moreover, we aligned these nine QTLs to the cowpea genome. Each of the QTLs was positioned in a genomic region ranging from 209,000 bp to 12,590,450 bp, and the QTL with the largest effect (21.8%) on heat tolerance, qHT4-1, was located within an interval of only 234,195 bp. These results provide SNP markers useful for marker-assisted selection for heat tolerance and lay a foundation for cloning, characterization, and applications of the genes controlling the cowpea heat tolerance for heat tolerance genetic improvement in cowpea and related crops.

Identification of Genomic Regions Associated with Fusarium Wilt Resistance in Cowpea

Fusarium wilt (FW), caused by the soil-borne fungal pathogen Fusarium oxysporum f. sp. Tracheiphilum, is a serious threat to cowpea production worldwide. Understanding the genetic architecture of FW resistance is a prerequisite to combatting this disease and developing FW resistance varieties. In the current study, a genetic diversity panel of 99 cowpea accessions was collected, and they were infected by a single strain, FW-HZ. The disease index (DI) based on the two indicators of leaf damage (LFD) and vascular discoloration (VD) varied highly across the population: most accessions were susceptible, and only seven accessions showed resistant phenotypes by both indicators. Through a genome-wide association study (GWAS), 3 and 7 single nucleotide polymorphisms (SNPs) significantly associated with LFD and VD were detected, respectively, which were distributed on chromosomes 3, 4, 5, 6 and 9, accounting for 0.68–13.92% of phenotypic variation. Based on the cowpea reference genome, 30 putative genes were identified and proposed as the likely candidates, including leucine-rich repeat protein kinase family protein, protein kinase superfamily protein and zinc finger family protein. These results provide novel insights into the genetic architecture of FW resistance and a basis for molecular breeding of FW resistant cultivars in cowpea.

Breeding of Vegetable Cowpea for Nutrition and Climate Resilience in Sub-Saharan Africa: Progress, Opportunities, and Challenges

Currently, the world population is increasing, and humanity is facing food and nutritional scarcity. Climate change and variability are a major threat to global food and nutritional security, reducing crop productivity in the tropical and subtropical regions of the globe. Cowpea has the potential to make a significant contribution to global food and nutritional security. In addition, it can be part of a sustainable food system, being a genetic resource for future crop improvement, contributing to resilience and improving agricultural sustainability under climate change conditions. In malnutrition prone regions of sub-Saharan Africa (SSA) countries, cowpea has become a strategic dryland legume crop for addressing food insecurity and malnutrition. Therefore, this review aims to assess the contribution of cowpea to SSA countries as a climate-resilient crop and the existing production challenges and perspectives. Cowpea leaves and immature pods are rich in diverse nutrients, with high levels of protein, vitamins, macro and micronutrients, minerals, fiber, and carbohydrates compared to its grain. In addition, cowpea is truly a multifunctional crop for maintaining good health and for reducing non-communicable human diseases. However, as a leafy vegetable, cowpea has not been researched and promoted sufficiently because it has not been promoted as a food security crop due to its low yield potential, susceptibility to biotic and abiotic stresses, quality assurance issues, policy regulation, and cultural beliefs (it is considered a livestock feed). The development of superior cowpea as a leafy vegetable can be approached in different ways, such as conventional breeding and gene stacking, speed breeding, mutation breeding, space breeding, demand-led breeding, a pan-omics approach, and local government policies. The successful breeding of cowpea genotypes that are high-yielding with a good nutritional value as well as having resistance to biotics and tolerant to abiotic stress could also be used to address food security and malnutrition-related challenges in sub-Saharan Africa.

Climate Change and Seed System Interventions Impact on Food Security and Incomes in East Africa

Climate change is challenging agricultural productivity, especially in Africa. Adoption of improved or diverse seed varieties is a promising strategy to achieve increased yields, support food security and reduce poverty under climate change in East Africa. However, rigorous impact evaluations linking the contributions of improved seeds to the welfare of households have been limited. This paper evaluates the impact of diversified seed systems on farm household production, sales, income, consumption and seed storage in Kenya and Uganda. It applies four-cell analysis to explore the intra-specific diversity of crops within farming systems, using primary data obtained from a random sampling of 207 treatment households and 87 control households. Propensity score matching was used to investigate the relationship between adoption of improved seeds and changes in production, sales, income, consumption, seed storage and food security. Econometric results indicate that treatment households using improved seeds saw a significant positive impact on income from bean seed sales, sorghum and millet consumption, bean livestock feed and maize and millet seed stored. We conclude that increasing seed diversity helps farmers cope with climate change and increases productivity, food availability, incomes and food security. Partnerships among seed improvement stakeholders need to be enhanced to ensure a continued supply of appropriate seeds to farmers.

The WRKY transcription factor family in cowpea: Genomic characterization and transcriptomic profiling under root dehydration

Cowpea [Vigna unguiculata (L.) Walp.] is one of the most tolerant legume crops to drought and salt stresses. WRKY transcription factor (TF) family members stand out among plant transcriptional regulators related to abiotic stress tolerance. However, little information is currently available on the expression of the cowpea WRKY gene family (VuWRKY) in response to water deficit. Thus, we analyzed genomic and transcriptomic data from cowpea to identify VuWRKY members and characterize their structure and transcriptional response under root dehydration stress. Ninety-two complete VuWRKY genes were found in the cowpea genome based on their domain characteristics. They were clustered into three groups: I (15 members), II (58), and III (16), while three genes were unclassified. Domain analysis of the encoded proteins identified four major variants of the conserved heptapeptide motif WRKYGQK. In silico analysis of VuWRKY gene promoters identified eight candidate binding motifs of cis-regulatory elements, regulated mainly by six TF families associated with abiotic stress responses. Ninety-seven VuWRKY modulated splicing variants associated with 55 VuWRKY genes were identified via RNA-Seq analysis available at the Cowpea Genomics Consortium (CpGC) database. qPCR analyses showed that 22 genes are induced under root dehydration, with VuWRKY18, 21, and 75 exhibiting the most significant induction levels. Given their central role in activating signal transduction cascades in abiotic stress response, the data provide a foundation for the targeted modification of specific VuWRKY family members to improve drought tolerance in this important climate-resilient legume in the developing world and beyond.

Genetic architecture of salt tolerance in a Multi-Parent Advanced Generation Inter-Cross (MAGIC) cowpea population

Background

Previous reports have shown that soil salinity is a growing threat to cowpea production, and thus the need for breeding salt-tolerant cowpea cultivars. A total of 234 Multi-Parent Advanced Generation Inter-Cross (MAGIC) lines along with their 8 founders were evaluated for salt tolerance under greenhouse conditions. The objectives of this study were to evaluate salt tolerance in a multi-parent advanced generation inter-cross (MAGIC) cowpea population, to identify single nucleotide polymorphism (SNP) markers associated with salt tolerance, and to assess the accuracy of genomic selection (GS) in predicting salt tolerance, and to explore possible epistatic interactions affecting salt tolerance in cowpea. Phenotyping was validated through the use of salt-tolerant and salt-susceptible controls that were previously reported. Genome-wide association study (GWAS) was conducted using a total of 32,047 filtered SNPs. The epistatic interaction analysis was conducted using the PLINK platform.

Results

Results indicated that: (1) large variation in traits evaluated for salt tolerance was identified among the MAGIC lines, (2) a total of 7, 2, 18, 18, 3, 2, 5, 1, and 23 were associated with number of dead plants, salt injury score, leaf SPAD chlorophyll under salt treatment, relative tolerance index for leaf SPAD chlorophyll, fresh leaf biomass under salt treatment, relative tolerance index for fresh leaf biomass, relative tolerance index for fresh stem biomass, relative tolerance index for the total above-ground fresh biomass, and relative tolerance index for plant height, respectively, with overlapping SNP markers between traits, (3) candidate genes encoding for proteins involved in ion transport such as Na+/Ca2+ K+ independent exchanger and H+/oligopeptide symporter were identified, and (4) epistatic interactions were identified.

Conclusions

These results will have direct applications in breeding programs aiming at improving salt tolerance in cowpea through marker-assisted selection. To the best of our knowledge, this study was one of the earliest reports using a MAGIC population to investigate the genetic architecture of salt tolerance in cowpea.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08332-y.

No disruption of rhizobial symbiosis during early stages of cowpea domestication

Modern agriculture intensely selects aboveground plant structures, while often neglecting belowground features, and evolutionary tradeoffs between these traits are predicted to disrupt host control over microbiota. Moreover, drift, inbreeding, and relaxed selection for symbiosis in crops might degrade plant mechanisms that support beneficial microbes. We studied the impact of domestication on the nitrogen fixing symbiosis between cowpea and root-nodulating Bradyrhizobium. We combined genome-wide analyses with a greenhouse inoculation study to investigate genomic diversity, heritability, and symbiosis trait variation among wild and early-domesticated cowpea genotypes. Cowpeas experienced modest decreases in genome-wide diversity during early domestication. Nonetheless, domesticated cowpeas responded efficiently to variation in symbiotic effectiveness, by forming more root nodules with nitrogen-fixing rhizobia and sanctioning non-fixing strains. Domesticated populations invested a larger proportion of host tissues into root nodules than wild cowpeas. Unlike soybean and wheat, cowpea showed no compelling evidence for degradation of symbiosis during domestication. Domesticated cowpeas experienced a less severe bottleneck than these crops and the low nutrient conditions in Africa where cowpea landraces were developed likely favored plant genotypes that gain substantial benefits from symbiosis. Breeders have largely neglected symbiosis traits, but artificial selection for improved plant responses to microbiota could increase plant performance and sustainability. This article is protected by copyright. All rights reserved.

Towards Sustainable Shifts to Healthy Diets and Food Security in Sub-Saharan Africa with Climate-Resilient Crops in Bread-Type Products: A Food System Analysis

Massive urbanization and increasing disposable incomes favor a rapid transition in diets and lifestyle in sub-Saharan Africa (SSA). As a result, the SSA population is becoming increasingly vulnerable to the double burden of malnutrition and obesity. This, combined with the increasing pressure to produce sufficient food and provide employment for this growing population together with the threat of climate change-induced declining crop yields, requires urgent sustainable solutions. Can an increase in the cultivation of climate-resilient crops (CRCs) and their utilization to produce attractive, convenient and nutritious bread products contribute to climate change adaptation and healthy and sustainable diets? A food system analysis of the bread food value chain in SSA indicates that replacement of refined, mostly imported, wheat in attractive bread products could (1) improve food and nutrition security, (2) bring about a shift to more nutritionally balanced diets, (3) increase economic inclusiveness and equitable benefits, and (4) improve sustainability and resilience of the food system. The food system analysis also provided systematic insight into the challenges and hurdles that need to be overcome to increase the availability, affordability and uptake of CRCs. Proposed interventions include improving the agronomic yield of CRCs, food product technology, raising consumer awareness and directing policies. Overall, integrated programs involving all stakeholders in the food system are needed.

Quantifying Nutrient Content in the Leaves of Cowpea Using Remote Sensing

Although hyperspectral remote sensing techniques have increasingly been used in the nutritional quantification of plants, it is important to understand whether the method shows a satisfactory response during the various phenological stages of the crop. The aim of this study was to quantify the levels of phosphorus (P), potassium (K), calcium (Ca) and zinc (Zn) in the leaves of Vigna Unguiculata (L.) Walp using spectral data obtained by a spectroradiometer. A randomised block design was used, with three treatments and twenty-five replications. The crop was evaluated at three growth stages: V4, R6 and R9. Single-band models were fitted using simple correlations. For the band ratio models, the wavelengths were selected by 2D correlation. For the models using partial least squares regression (PLSR), the stepwise method was used. The model showing the best fit was used to estimate the phosphorus content in the single-band (R² = 0.62; RMSE = 0.54 and RPD = 1.61), band ratio (R² = 0.66; RMSE = 0.65 and RPD = 1.52) and PLSR models, using data from each of the phenological stages (R² = 0.80; RMSE = 0.47 and RPD = 1.66). Accuracy in modelling leaf nutrients depends on the phenological stage, as well as the amount of data used, and is more accurate with a larger number of samples.

Gene activation via Cre/lox-mediated excision in cowpea (Vigna unguiculata)

Expression of Cre recombinase by AtRps5a_pro or AtDD45_pro enabled Cre/lox-mediated recombination at an early embryonic developmental stage upon crossing, activating transgenes in the hybrid cowpea and tobacco. Genetic engineering ideally results in precise spatiotemporal control of transgene expression. To activate transgenes exclusively in a hybrid upon fertilization, we evaluated a Cre/lox-mediated gene activation system with the Cre recombinase expressed by either AtRps5a or AtDD45 promoters that showed activity in egg cells and young embryos. In crosses between Cre recombinase lines and transgenic lines harboring a lox-excision reporter cassette with ZsGreen driven by the AtUbq3 promoter after Cre/lox-mediated recombination, we observed complete excision of the lox-flanked intervening DNA sequence between the AtUbq3_pro and the ZsGreen coding sequence in F₁ progeny upon genotyping but no ZsGreen expression in F₁ seeds or seedlings. The incapability to observe ZsGreen fluorescence was attributed to the activity of the AtUbq3_pro. Strong ZsGreen expression in F₁ seeds was observed after recombination when ZsGreen was driven by the AtUbq10 promoter. Using the AtDD45_pro to express Cre resulted in more variation in recombination frequencies between transgenic lines and crosses. Regardless of the promoter used to regulate Cre, mosaic F₁ progeny were rare, suggesting gene activation at an early embryo-developmental stage. Observation of ZsGreen-expressing tobacco embryos at the globular stage from crosses with the AtRps5a_proCre lines pollinated by the AtUbq3_prolox line supported the early activation mode.

Constraints and Prospects of Improving Cowpea Productivity to Ensure Food, Nutritional Security and Environmental Sustainability

Providing safe and secure food for an increasing number of people globally is challenging. Coping with such a human population by merely applying the conventional agricultural production system has not proved to be agro-ecologically friendly; nor is it sustainable. Cowpea (Vigna unguiculata (L) Walp) is a multi-purpose legume. It consists of high-quality protein for human consumption, and it is rich in protein for livestock fodder. It enriches the soil in that it recycles nutrients through the fixation of nitrogen in association with nodulating bacteria. However, the productivity of this multi-functional, indigenous legume that is of great value to African smallholder farmers and the rural populace, and also to urban consumers and entrepreneurs, is limited. Because cowpea is of strategic importance in Africa, there is a need to improve on its productivity. Such endeavors in Africa are wrought with challenges that include drought, salinity, the excessive demand among farmers for synthetic chemicals, the repercussions of climate change, declining soil nutrients, microbial infestations, pest issues, and so forth. Nevertheless, giant strides have already been made and there have already been improvements in adopting sustainable and smart biotechnological approaches that are favorably influencing the production costs of cowpea and its availability. As such, the prospects for a leap in cowpea productivity in Africa and in the enhancement of its genetic gain are good. Potential and viable means for overcoming some of the above-mentioned production constraints would be to focus on the key cowpea producer nations in Africa and to encourage them to embrace biotechnological techniques in an integrated approach to enhance for sustainable productivity. This review highlights the spectrum of constraints that limit the cowpea yield, but looks ahead of the constraints and seeks a way forward to improve cowpea productivity in Africa. More importantly, this review investigates applications and insights concerning mechanisms of action for implementing eco-friendly biotechnological techniques, such as the deployment of bio inoculants, applying climate-smart agricultural (CSA) practices, agricultural conservation techniques, and multi-omics smart technology in the spheres of genomics, transcriptomics, proteomics, and metabolomics, for improving cowpea yields and productivity to achieve sustainable agro-ecosystems, and ensuring their stability.

Molecular Fingerprinting and Hybridity Authentication in Cowpea Using Single Nucleotide Polymorphism Based Kompetitive Allele-Specific PCR Assay

Optimization of a breeding program for increased genetic gain requires quality assurance (QA) and quality control (QC) at key phases of the breeding process. One vital phase in a breeding program that requires QC and QA is the choice of parents and successful hybridizations to combine parental attributes and create variations. The objective of this study was to determine parental diversity and confirm hybridity of cowpea F₁ progenies using KASP (Kompetitive Allele-Specific PCR)-based single nucleotide polymorphism (SNP) markers. A total of 1,436 F₁ plants were derived from crossing 220 cowpea breeding lines and landraces to 2 elite sister lines IT99K-573-1-1 and IT99K-573-2-1 as male parents, constituting 225 cross combinations. The progenies and the parents were genotyped with 17 QC SNP markers via high-throughput KASP genotyping assay. The QC markers differentiated the parents with mean efficiency of 37.90% and a range of 3.4-82.8%, revealing unique fingerprints of the parents. Neighbor-Joining cladogram divided the 222 parents into 3 clusters. Genetic distances between parents ranged from 0 to 3.74 with a mean of 2.41. Principal component analysis (PCA) depicted a considerable overlap between parents and F₁ progenies with more scatters among parents than the F₁s. The differentiation among parents and F₁s was best contributed to by 82% of the markers. As expected, parents and F₁s showed a significant contrast in proportion of heterozygous individuals, with mean values of 0.02 and 0.32, respectively. KASP markers detected true hybridity with 100% success rate in 72% of the populations. Overall, 79% of the putative F₁ plants were true hybrids, 14% were selfed plants, and 7% were undetermined due to missing data and lack of marker polymorphism between parents. The study demonstrated an effective application of KASP-based SNP assay in fingerprinting, confirmation of hybridity, and early detection of false F₁ plants. The results further uncovered the need to deploy markers as a QC step in a breeding program.

Genetic, anatomical, and environmental patterns related to pod shattering resistance in domesticated cowpea [Vigna unguiculata (L.) Walp]

Pod shattering, which causes the explosive release of seeds from the pod, is one of the main sources of yield losses in cowpea in arid and semi-arid areas. Reduction of shattering has therefore been a primary target for selection during domestication and improvement of cowpea, among other species. Using a mini-core diversity panel of 368 cowpea accessions, four regions with a statistically significant association with pod shattering were identified. Two genes (Vigun03g321100 and Vigun11g100600), involved in cell wall biosynthesis, were identified as strong candidates for pod shattering. Microscopic analysis was conducted on a subset of accessions representing the full spectrum of shattering phenotypes. This analysis indicated that the extent of wall fiber deposition was highly correlated with shattering. The results from this study also demonstrate that pod shattering in cowpea is exacerbated by arid environmental conditions. Finally, using a subset of West African landraces, patterns of historical selection for shattering resistance related to precipitation in the environment of origin were identified. Together, these results shed light on sources of resistance to pod shattering, which will, in turn, improve climate resilience of a major global nutritional staple.

Interspecific Sample Prioritization Can Improve QTL Detection With Tree-Based Predictive Models

Due to increasing demand for new advanced crops, considerable efforts have been made to explore the improvement of stress and disease resistance cultivar traits through the study of wild crops. When both wild and interspecific hybrid materials are available, a common approach has been to study two types of materials separately and simply compare the quantitative trait locus (QTL) regions. However, combining the two types of materials can potentially create a more efficient method of finding predictive QTLs. In this simulation study, we focused on scenarios involving causal marker expression suppressed by trans-regulatory mechanisms, where the otherwise easily lost associated signals benefit the most from combining the two types of data. A probabilistic sampling approach was used to prioritize consistent genotypic phenotypic patterns across both types of data sets. We chose random forest and gradient boosting to apply the prioritization scheme and found that both facilitated the investigation of predictive causal markers in most of the biological scenarios simulated.

Disparate genetic variants associated with distinct components of cowpea resistance to the seed beetle Callosobruchus maculatus

Polygenic genome-wide association mapping identified two regions of the cowpea genome associated with different components of resistance to its major post-harvest pest, the seed beetle Callosobruchus maculatus. Cowpea (Vigna unguiculata) is an important grain and fodder crop in arid and semi-arid regions of Africa, Asia, and South America, where the cowpea seed beetle, Callosobruchus maculatus, is a serious post-harvest pest. Development of cultivars resistant to C. maculatus population growth in storage could increase grain yield and quality and reduce reliance on insecticides. Here, we use a MAGIC (multi-parent, advanced-generation intercross) population of cowpea consisting of 305 recombinant inbred lines (RILs) to identify genetic variants associated with resistance to seed beetles. Because inferences regarding the genetic basis of resistance may depend on the source of the pest or the assay protocol, we used two divergent geographic populations of C. maculatus and two complementary assays to measure several aspects of resistance. Using polygenic genome-wide association mapping models, we found that the cowpea RILs harbor substantial additive-genetic variation for most resistance measures. Variation in several components of resistance, including larval development time and survival, was largely explained by one or several linked loci on chromosome 5. A second region on chromosome 8 explained increased seed resistance via the induction of early-exiting larvae. Neither of these regions contained genes previously associated with resistance to insects that infest grain legumes. We found some evidence of gene-gene interactions affecting resistance, but epistasis did not contribute substantially to resistance variation in this mapping population. The combination of mostly high heritabilities and a relatively consistent and simple genetic architecture increases the feasibility of breeding for enhanced resistance to C. maculatus.

Limiting-Stress-Elimination Hypothesis: Using Non-hormonal Biostimulant to Reduce Stress and Increase Savanna Cowpea [Vigna unguiculata (L.) Walp.] Productivity

An alternative decision axiom to guide in determining the optimal intervention strategy to maximize cowpea production is proposed. According to the decrement from the maximum concept of Mitscherlich, the decrement from the maximum for each stressor must be minimized to produce the absolute maximum production. In crop production, this means all deficient nutrients must be supplemented to ensure maximum yield and laid the foundation in fertilizer formulation. However, its implementation is not economically feasible in many situations, particularly where multiple environmental factors impact crop productivity as in the case of low resource conditions. We propose and test the hypothesis that yield allocation will increase when the most limiting stressor among prevailing stressors is eliminated at least until the next limiting stressor impacts productivity. We selected drought limiting savanna conditions and cowpea (Vigna unguiculata), adapted to nitrogen dependence. To determine the limiting condition, we measured the response of cowpea to D-sorbitol, nitrogen, and non-hormonal biostimulant (nhB) treatments. The nhB treatment increased total biomass by 45% compared to nitrogen, 13%, and D-sorbitol, 17%, suggesting osmotic stress is more limiting in the observed savanna conditions. The effect of the biostimulant is due to antioxidants and key amino acids that stimulate metabolism and stress resistance. Where nitrogen becomes the next constraining factor, biostimulants can contribute organic nitrogen. The study supports the use of biostimulants as candidate intervention under conditions where crop productivity is limited by multiple or alternating constraints during crop growth.

Genome-Wide Association Study Reveals Candidate Genes for Flowering Time in Cowpea (Vigna unguiculata [L.] Walp.)

Cowpea (Vigna unguiculata [L.] Walp., diploid, 2n = 22) is a major crop used as a protein source for human consumption as well as a quality feed for livestock. It is drought and heat tolerant and has been bred to develop varieties that are resilient to changing climates. Plant adaptation to new climates and their yield are strongly affected by flowering time. Therefore, understanding the genetic basis of flowering time is critical to advance cowpea breeding. The aim of this study was to perform genome-wide association studies (GWAS) to identify marker trait associations for flowering time in cowpea using single nucleotide polymorphism (SNP) markers. A total of 368 accessions from a cowpea mini-core collection were evaluated in Ft. Collins, CO in 2019 and 2020, and 292 accessions were evaluated in Citra, FL in 2018. These accessions were genotyped using the Cowpea iSelect Consortium Array that contained 51,128 SNPs. GWAS revealed seven reliable SNPs for flowering time that explained 8-12% of the phenotypic variance. Candidate genes including FT, GI, CRY2, LSH3, UGT87A2, LIF2, and HTA9 that are associated with flowering time were identified for the significant SNP markers. Further efforts to validate these loci will help to understand their role in flowering time in cowpea, and it could facilitate the transfer of some of this knowledge to other closely related legume species.

Developing a rapid and highly efficient cowpea regeneration, transformation and genome editing system using embryonic axis explants

Cowpea (Vigna unguiculata (L.) Walp.) is one of the most important legume crops planted worldwide, but despite decades of effort, cowpea transformation is still challenging due to inefficient Agrobacterium-mediated transfer DNA delivery, transgenic selection and in vitro shoot regeneration. Here, we report a highly efficient transformation system using embryonic axis explants isolated from imbibed mature seeds. We found that removal of the shoot apical meristem from the explants stimulated direct multiple shoot organogenesis from the cotyledonary node tissue. The application of a previously reported ternary transformation vector system provided efficient Agrobacterium-mediated gene delivery, while the utilization of spcN as selectable marker enabled more robust transgenic selection, plant recovery and transgenic plant generation without escapes and chimera formation. Transgenic cowpea plantlets developed exclusively from the cotyledonary nodes at frequencies of 4% to 37% across a wide range of cowpea genotypes. CRISPR/Cas-mediated gene editing was successfully demonstrated. The transformation principles established here could also be applied to other legumes to increase transformation efficiencies.

Developing a rapid and highly efficient cowpea regeneration, transformation and genome editing system using embryonic axis explants

Cowpea (Vigna unguiculata (L.) Walp.) is one of the most important legume crops planted worldwide, but despite decades of effort, cowpea transformation is still challenging due to inefficient Agrobacterium-mediated transfer DNA delivery, transgenic selection and in vitro shoot regeneration. Here, we report a highly efficient transformation system using embryonic axis explants isolated from imbibed mature seeds. We found that removal of the shoot apical meristem from the explants stimulated direct multiple shoot organogenesis from the cotyledonary node tissue. The application of a previously reported ternary transformation vector system provided efficient Agrobacterium-mediated gene delivery, while the utilization of spcN as selectable marker enabled more robust transgenic selection, plant recovery and transgenic plant generation without escapes and chimera formation. Transgenic cowpea plantlets developed exclusively from the cotyledonary nodes at frequencies of 4% to 37% across a wide range of cowpea genotypes. CRISPR/Cas-mediated gene editing was successfully demonstrated. The transformation principles established here could also be applied to other legumes to increase transformation efficiencies.

Genotyping-by-Sequencing in Vigna unguiculata Landraces and Its Utility for Assessing Taxonomic Relationships

Genotyping by sequencing (GBS) was used to analyze relationships among cowpea and asparagus bean landraces from southern Italy and to assess the utility of this technology to study taxonomy in a wider panel, including V. unguiculata cultigroups, subspecies, and other Vigna species. The analysis of SNPs derived from GBS highlighted that, among the cowpea landraces, the African samples were separated from the other material, while, for the Italian landraces, a certain clustering depending on seed color/pattern was observed in the dendrogram. When examining the V. unguiculata species complex, a clear separation between the two groups of wild subspecies, i.e., the allogamous wild perennials and the perennial out/inbreds, could be observed, the former representing the more ancestral wild progenitors of V. unguiculata. The species V. vexillata appeared more closely related to V. unguiculata than to the other Vigna species analyzed.

Orthology and synteny analysis of receptor-like kinases "RLK" and receptor-like proteins "RLP" in legumes

BACKGROUND

Legume species are an important plant model because of their protein-rich physiology. The adaptability and productivity of legumes are limited by major biotic and abiotic stresses. Responses to these stresses directly involve plasma membrane receptor proteins known as receptor-like kinases and receptor-like proteins. Evaluating the homology relations among RLK and RLP for seven legume species, and exploring their presence among synteny blocks allow an increased understanding of evolutionary relations, physical position, and chromosomal distribution in related species and their shared roles in stress responses.

RESULTS

Typically, a high proportion of RLK and RLP legume proteins belong to orthologous clusters, which is confirmed in this study, where between 66 to 90% of the RLKs and RLPs per legume species were classified in orthologous clusters. One-third of the evaluated syntenic blocks had shared RLK/RLP genes among both legumes and non-legumes. Among the legumes, between 75 and 98% of the RLK/RLP were present in syntenic blocks. The distribution of chromosomal segments between Phaseolus vulgaris and Vigna unguiculata, two species that diverged ~ 8 mya, were highly similar. Among the RLK/RLP synteny clusters, seven experimentally validated resistance RLK/RLP genes were identified in syntenic blocks. The RLK resistant genes FLS2, BIR2, ERECTA, IOS1, and AtSERK1 from Arabidopsis and SLSERK1 from Solanum lycopersicum were present in different pairwise syntenic blocks among the legume species. Meanwhile, only the LYM1- RLP resistant gene from Arabidopsis shared a syntenic blocks with Glycine max.

CONCLUSIONS

The orthology analysis of the RLK and RLP suggests a dynamic evolution in the legume family, with between 66 to 85% of RLK and 83 to 88% of RLP belonging to orthologous clusters among the species evaluated. In fact, for the 10-species comparison, a lower number of singleton proteins were reported among RLP compared to RLK, suggesting that RLP positions are more physically conserved compared to RLK. The identification of RLK and RLP genes among the synteny blocks in legumes revealed multiple highly conserved syntenic blocks on multiple chromosomes. Additionally, the analysis suggests that P. vulgaris is an appropriate anchor species for comparative genomics among legumes.

Loci discovery, network-guided approach, and genomic prediction for drought tolerance index in a multi-parent advanced generation intercross (MAGIC) cowpea population

Cowpea is a nutrient-dense legume that significantly contributes to the population's diet in sub-Saharan Africa and other regions of the world. Improving cowpea cultivars to be more resilient to abiotic stress such as drought would be of great importance. The use of a multi-parent advanced generation intercross (MAGIC) population has been shown to be efficient in increasing the frequency of rare alleles that could be associated with important agricultural traits. In addition, drought tolerance index has been reported to be a reliable parameter for assessing crop tolerance to water-deficit conditions. Therefore, the objectives of this study were to evaluate the drought tolerance index for plant growth habit, plant maturity, flowering time, 100-seed weight, and grain yield in a MAGIC cowpea population, to conduct genome-wide association study (GWAS) and identify single nucleotide polymorphism (SNP) markers associated with the drought tolerance indices, to investigate the potential relationship existing between the significant loci associated with the drought tolerance indices, and to conduct genomic selection (GS). These analyses were performed using the existing phenotypic and genotypic data published for the MAGIC population which consisted of 305 F8 recombinant inbred lines (RILs) developed at University of California, Riverside. The results indicated that: (1) large variation in drought tolerance indices existed among the cowpea genotypes, (2) a total of 14, 18, 5, 5, and 35 SNPs were associated with plant growth habit change due to drought stress, and drought tolerance indices for maturity, flowering time, 100-seed weight, and grain yield, respectively, (3) the network-guided approach revealed clear interactions between the loci associated with the drought tolerance traits, and (4) the GS accuracy varied from low to moderate. These results could be applied to improve drought tolerance in cowpea through marker-assisted selection (MAS) and genomic selection (GS). To the best of our knowledge, this is the first report on marker loci associated with drought tolerance indices in cowpea.

Genomic resources in plant breeding for sustainable agriculture

Climate change during the last 40 years has had a serious impact on agriculture and threatens global food and nutritional security. From over half a million plant species, cereals and legumes are the most important for food and nutritional security. Although systematic plant breeding has a relatively short history, conventional breeding coupled with advances in technology and crop management strategies has increased crop yields by 56 % globally between 1965-85, referred to as the Green Revolution. Nevertheless, increased demand for food, feed, fiber, and fuel necessitates the need to break existing yield barriers in many crop plants. In the first decade of the 21st century we witnessed rapid discovery, transformative technological development and declining costs of genomics technologies. In the second decade, the field turned towards making sense of the vast amount of genomic information and subsequently moved towards accurately predicting gene-to-phenotype associations and tailoring plants for climate resilience and global food security. In this review we focus on genomic resources, genome and germplasm sequencing, sequencing-based trait mapping, and genomics-assisted breeding approaches aimed at developing biotic stress resistant, abiotic stress tolerant and high nutrition varieties in six major cereals (rice, maize, wheat, barley, sorghum and pearl millet), and six major legumes (soybean, groundnut, cowpea, common bean, chickpea and pigeonpea). We further provide a perspective and way forward to use genomic breeding approaches including marker-assisted selection, marker-assisted backcrossing, haplotype based breeding and genomic prediction approaches coupled with machine learning and artificial intelligence, to speed breeding approaches. The overall goal is to accelerate genetic gains and deliver climate resilient and high nutrition crop varieties for sustainable agriculture.

Developments and Prospects in Imperative Underexploited Vegetable Legumes Breeding: A Review

Vegetable legumes are an essential source of carbohydrates, vitamins, and minerals, along with health-promoting bioactive chemicals. The demand for the use of either fresh or processed vegetable legumes is continually expanding on account of the growing consumer awareness about their well-balanced diet. Therefore, sustaining optimum yields of vegetable legumes is extremely important. Here we seek to present d etails of prospects of underexploited vegetable legumes for food availability, accessibility, and improved livelihood utilization. So far research attention was mainly focused on pulse legumes' performance as compared to vegetable legumes. Wild and cultivated vegetable legumes vary morphologically across diverse habitats. This could make them less known, underutilized, and underexploited, and make them a promising potential nutritional source in developing nations where malnutrition still exists. Research efforts are required to promote underexploited vegetable legumes, for improving their use to feed the ever-increasing population in the future. In view of all the above points, here we have discussed underexploited vegetable legumes with tremendous potential; namely, vegetable pigeon pea (Cajanus cajan), cluster bean (Cyamopsis tetragonoloba), winged bean (Psophocarpus tetragonolobus), dolichos bean (Lablab purpureus), and cowpea (Vigna unguiculata), thereby covering the progress related to various aspects such as pre-breeding, molecular markers, quantitative trait locus (QTLs), genomics, and genetic engineering. Overall, this review has summarized the information related to advancements in the breeding of vegetable legumes which will ultimately help in ensuring food and nutritional security in developing nations.

A receptor-like protein mediates plant immune responses to herbivore-associated molecular patterns

Herbivory is fundamental to the regulation of both global food webs and the extent of agricultural crop losses. Induced plant responses to herbivores promote resistance and often involve the perception of specific herbivore-associated molecular patterns (HAMPs); however, precisely defined receptors and elicitors associated with herbivore recognition remain elusive. Here, we show that a receptor confers signaling and defense outputs in response to a defined HAMP common in caterpillar oral secretions (OS). Staple food crops, including cowpea (Vigna unguiculata) and common bean (Phaseolus vulgaris), specifically respond to OS via recognition of proteolytic fragments of chloroplastic ATP synthase, termed inceptins. Using forward-genetic mapping of inceptin-induced plant responses, we identified a corresponding leucine-rich repeat receptor, termed INR, specific to select legume species and sufficient to confer inceptin-induced responses and enhanced defense against armyworms (Spodoptera exigua) in tobacco. Our results support the role of plant immune receptors in the perception of chewing herbivores and defense.

Motivating Pulse-Centric Eating Patterns to Benefit Human and Environmental Well-Being

Pulses (e.g., lentil, common bean, chickpea, and dry pea) are linked to a myriad of positive human and environmental health impacts, making them an ideal food for wise and conscientious global citizens. In addition, pulses are affordable and shelf-stable. The combination of these factors, an elevated consumer interest in plant-based diets, and the COVID-19 pandemic resulted in increased purchasing of pulses and even empty grocery store shelves. Although pulses have many associated benefits, some consumers are hesitant to regularly eat pulses, claiming concerns of abdominal discomfort or a lack of knowledge on how to best prepare pulses. To capitalize on increased consumer interest and purchasing of pulses, now is the time for outreach efforts that address these concerns and the positive outcomes associated with pulses, thereby promoting public and environmental health. Consumers must actively decide to add pulses to their grocery lists and incorporate them into their regular eating patterns. Motivation to adopt new eating habits is essential because knowledge alone does not result in behavior change. Thus, to mitigate perceived barriers and drive consumption, we suggest application of the Information-Motivation-Behavioral Skills Model and emphasis of three main benefits of pulses as motivators: (1) culinary versatility, (2) sustainability, and (3) healthfulness.

Population Structure and Genetic Diversity in Korean Cowpea Germplasm Based on SNP Markers

Cowpea is one of the most essential legume crops providing inexpensive dietary protein and nutrients. The aim of this study was to understand the genetic diversity and population structure of global and Korean cowpea germplasms. A total of 384 cowpea accessions from 21 countries were genotyped with the Cowpea iSelect Consortium Array containing 51,128 single-nucleotide polymorphisms (SNPs). After SNP filtering, a genetic diversity study was carried out using 35,116 SNPs within 376 cowpea accessions, including 229 Korean accessions. Based on structure and principal component analysis, a total of 376 global accessions were divided into four major populations. Accessions in group 1 were from Asia and Europe, those in groups 2 and 4 were from Korea, and those in group 3 were from West Africa. In addition, 229 Korean accessions were divided into three major populations (Q1, Jeonra province; Q2, Gangwon province; Q3, a mixture of provinces). Additionally, the neighbor-joining tree indicated similar results. Further genetic diversity analysis within the global and Korean population groups indicated low heterozygosity, a low polymorphism information content, and a high inbreeding coefficient in the Korean cowpea accessions. The population structure analysis will provide useful knowledge to support the genetic potential of the cowpea breeding program, especially in Korea.

Developing the role of legumes in West Africa under climate change

West Africa is faced with significant challenges from climate change, including parts of the region becoming hotter with more variable rainfall. The Sahelian region in particular is already subject to severe droughts. To address this better adapted crop varieties (such as for cowpea) are clearly a central element, a complementary one is a greater use of resilient alternative crops especially underutilized legumes particularly Bambara groundnut, African yam bean, winged bean and Kersting's groundnut. Genetic diversity of these crops conserved in genebanks and farmer's field provides an opportunity to exploit climate resilient traits using cutting-edge genomic tools and to use genomics-assisted breeding to accelerate genetic gains in combination of rapid cycle breeding strategy to develop climate-resilient cultivars for sub-Saharan Africa.

Breaks of macrosynteny and collinearity among moth bean (Vigna aconitifolia), cowpea (V. unguiculata), and common bean (Phaseolus vulgaris)

Comparative cytogenetic mapping is a powerful approach to gain insights into genome organization of orphan crops, lacking a whole sequenced genome. To investigate the cytogenomic evolution of important Vigna and Phaseolus beans, we built a BAC-FISH (fluorescent in situ hybridization of bacterial artificial chromosome) map of Vigna aconitifolia (Vac, subgenus Ceratotropis), species with no sequenced genome, and compared with V. unguiculata (Vu, subgenus Vigna) and Phaseolus vulgaris (Pv) maps. Seventeen Pv BACs, eight Vu BACs, and 5S and 35S rDNA probes were hybridized in situ on the 11 Vac chromosome pairs. Five Vac chromosomes (Vac6, Vac7, Vac9, Vac10, and Vac11) showed conserved macrosynteny and collinearity between V. unguiculata and P. vulgaris. On the other hand, we observed collinearity breaks, identified by pericentric inversions involving Vac2 (Vu2), Vac4 (Vu4), and Vac3 (Pv3). We also detected macrosynteny breaks of translocation type involving chromosomes 1 and 8 of V. aconitifolia and P. vulgaris; 2 and 3 of V. aconitifolia and P. vulgaris; and 1 and 5 of V. aconitifolia and V. unguiculata. Considering our data and previous BAC-FISH studies, six chromosomes (1, 2, 3, 4, 5, and 8) are involved in major karyotype divergences between genera and five (1, 2, 3, 4, and 5) between Vigna subgenera, including mechanisms such as duplications, inversions, and translocations. Macrosynteny breaks between Vigna and Phaseolus suggest that the major chromosomal rearrangements have occurred within the Vigna clade. Our cytogenomic comparisons bring new light on the degree of shared macrosynteny and mechanisms of karyotype diversification during Vigna and Phaseolus evolution.