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.

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.

Cowpea (Vigna unguiculata L. Walp., 2n = 2x = 22) is a protein-rich crop that complements staple cereals for humans and serves as fodder for livestock. It is widely grown in Africa and other developing countries as the primary source of protein in the diet; therefore, it is necessary to identify the protein-related loci to improve cowpea breeding. In the current study, we conducted a genome-wide association study (GWAS) on 161 cowpea accessions (151 USDA germplasm plus 10 Arkansas breeding lines) with a wide range of seed protein contents (21.8~28.9%) with 110,155 high-quality whole-genome single-nucleotide polymorphisms (SNPs) to identify markers associated with protein content, then performed genomic prediction (GP) for future breeding. A total of seven significant SNP markers were identified using five GWAS models (single-marker regression (SMR), the general linear model (GLM), Mixed Linear Model (MLM), Fixed and Random Model Circulating Probability Unification (FarmCPU), and Bayesian-information and Linkage-disequilibrium Iteratively Nested Keyway (BLINK), which are located at the same locus on chromosome 8 for seed protein content. This locus was associated with the gene Vigun08g039200, which was annotated as the protein of the thioredoxin superfamily, playing a critical function for protein content increase and nutritional quality improvement. In this study, a genomic prediction (GP) approach was employed to assess the accuracy of predicting seed protein content in cowpea. The GP was conducted using cross-prediction with five models, namely ridge regression best linear unbiased prediction (rrBLUP), Bayesian ridge regression (BRR), Bayesian A (BA), Bayesian B (BB), and Bayesian least absolute shrinkage and selection operator (BL), applied to seven random whole genome marker sets with different densities (10 k, 5 k, 2 k, 1 k, 500, 200, and 7), as well as significant markers identified through GWAS. The accuracies of the GP varied between 42.9% and 52.1% across the seven SNPs considered, depending on the model used. These findings not only have the potential to expedite the breeding cycle through early prediction of individual performance prior to phenotyping, but also offer practical implications for cowpea breeding programs striving to enhance seed protein content and nutritional quality.

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.

Guar [Cyamopsis tetragonoloba (L.) Taub] was domesticated in India and Pakistan. It is mainly self-pollinated, bushy, and deeply tap-rooted. Guar seed endosperm contains galactomannan gum, which is used in many food products, pharmaceuticals, cosmetics, explosives, meat products, and pet foods, and in the textile industry, yet its genetic diversity remains largely underexplored. Using 7000 high-quality single nucleotide polymorphism (SNP) markers acquired from genotyping by sequencing (GBS), we analyzed the genetic diversity and population structure in 225 guar accessions from India, Pakistan, and the United States. Structure Harvester revealed that K = 3 had the best delta K, whereas K = 2 had the second-highest delta K. Three major genetic clusters (K = 3) were identified using population structure analysis, utilizing an admixture model: 156 accessions (69.3%) were classified into Q1, 23 accessions (10.2%) in Q2, and 16 accessions (7.1%) in Q3. The remaining 30 accessions (13.3%) were included in the admixture. In all three of the subpopulations at K = 3, most of the guar accessions came from India. We also found that these clusters mostly correlated with geographic origins. Results showed that the Q2 and Q3 subpopulations included only guar accessions from India. Genetic resources from Q2 and Q3 may represent an untapped reservoir for introducing beneficial variety into the U.S. breeding populations. This genetic diversity and population structure analysis of the guar gene pool will be of interest to conduct allele mining and donor parent selection for the development of new and better guar germplasm for desired traits.