Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Plant Sci., 25 September 2025

Sec. Plant Breeding

Volume 16 - 2025 | https://doi.org/10.3389/fpls.2025.1661440

This article is part of the Research TopicAdvances in Breeding for Quantitative Disease Resistance - Volume IIView all 4 articles

Discovery and development of single-nucleotide polymorphism markers for resistance to Striga gesnerioides in cowpea (Vigna unguiculata)

  • 1International Institute of Tropical Agriculture (IITA), Kano, Nigeria
  • 2International Institute of Tropical Agriculture (IITA), Ibadan, Nigeria

Introduction: The parasitic weed [Striga gesnerioides (Willd.) Vatke] is a principal biotic constraint to cowpea [Vigna unguiculata (L.) Walp.] production in West and Central Africa, causing severe yield reductions. Multiple races of S. gesnerioides exist across the cowpea-growing areas of the sub-region. Past efforts identified some resistant sources and race-specific genes underpinning Striga resistance, but deployment of associated markers in breeding is limited. Here, we utilized a 51K cowpea iSelect single-nucleotide polymorphisms (SNPs) to decipher genomic regions underlying Striga resistance and explore marker conversion and validation for easy deployment.

Method: The study used two-year phenotypic data on a minicore panel of 368 cowpea genotypes screened at two sites in Northern Nigeria. SNPs performances were verified and validated using two independent sets of 60 and 20 diverse genotypes respectively.

Results: The minicore displayed apparent differences in response to the S. gesnerioides attack. A genome-wide scan uncovered a primary gene effect signal on chromosome Vu11 and minor regions on chromosomes Vu02, Vu03, Vu07, Vu09 and Vu10. The major effect region on Vu11 harbored a coil-coil nucleotide-binding site leucine-rich repeat (CC-NBS-LRR) protein, encoded by the RSG3–301 gene, previously implicated in race-specific resistance to S. gesnerioides in cowpea. The associated SNPs were successfully converted into Kompetitive Allele-Specific PCR (KASP) assays and validated using 20 independent diverse cowpea genotypes. Five KASP markers, snpVU00075, snpVU00076, snpVU00077, snpVU00078, and snpVU00079, depicted consistent and significant associations with the phenotype in the validation set.

Discussion: The markers provide valuable tools for efficient marker-assisted selection (MAS) in breeding programs focused on developing Striga-resistant cowpea varieties.

1 Introduction

Cowpea (Vigna unguiculata L. Walp.) is a vital legume crop widely grown in sub-Saharan Africa, Asia, and Latin America (Abate et al., 2012; Singh et al., 2014). It is a self-pollinated diploid with 2n=22 chromosomes and a genome size of 640.6 Mbp (Lonardi et al., 2019). Cowpea has synonymous names worldwide, including black-eyed pea, lubia, Kathir pea, China pea, crowder pea, niebe, and southern pea (Heuzé et al., 2015). Recent evidence traces the cowpea center of domestication to Nigeria in West Africa (WA) (Vaillancourt and Weeden, 1992; Panzeri et al., 2022). The crop is grown in over 88 countries worldwide, with Nigeria being the largest producer, accounting for 46% of world production. Cowpea is an essential food and income source for more than 200 million smallholder farmers, particularly in semi-arid regions where other crops may not thrive due to limited water availability (Abate et al., 2012; Singh, 2014; da Silva et al., 2018). Cowpea is a highly nutritious crop rich in protein, dietary fiber, vitamins, and minerals (Singh et al., 2014; Jayathilake et al., 2018). The crop is also an essential source of nitrogen for the soil, making it a valuable component of sustainable cropping systems (Blade et al., 1997; Mortimore et al., 1997); it serves as feed for livestock and as a cover crop to protect and enrich the soil (Abate et al., 2012).

Despite its importance, cowpea production is often constrained by a range of biotic and abiotic stresses, including parasitic weeds (Striga and Alectra), insect pests and diseases, drought, and poor soil fertility (Atokple et al., 1995; Fatokun et al., 2002; Li et al., 2009; Boukar et al., 2013; Nkomo et al., 2021). Among these constraints, the parasitic weed [Striga gesnerioides (Willd.) Vatke], is quite devastating on cowpea, especially in the Sudano-Sahelian belt of West Africa, where the crop is most widely cultivated (Parker and Riches, 1993; Mohamed et al., 2001; Li et al., 2009; CABI, 2020). Cowpea yield loss to S. gesnerioides ranges from moderate to total crop loss in some parts of Nigeria, Niger, and Burkina Faso (CABI, 2020).

Striga gesnerioides belongs to the Orobanchaceae family (CABI, 2020). It is an obligate parasite with tiny seeds, unable to establish itself without the help of a host plant CABI, 2020). Germination depends on a period of moist conditioning and exposure to germination stimulants in the host plant’s root exudates, the most important of which is alectrol, a stimulant for both the Striga and related parasite Alectra vogelii (Müller et al., 1992; CABI, 2020). Over 28 Striga species and six subspecies have been characterized, among which purple witchweed [Striga hermonthica (Delile) Benth.], Asiatic witchweed [Striga asiatica (L.) Kuntze], and S. gesnerioides are the most economically important (Mohamed et al., 2001). Striga hermonthica and S. asiatica primarily infect cereals in the Poaceae family, while the primary hosts of S. gesnerioides are cowpea and wild legume species (Ohlson and Timko, 2020). However, S. gesnerioides can parasitize hosts across genera, including Ipomea, Jaquemontia, Merremia, Euphorbia, and Nicotiana (Mohamed et al., 2001).

Several control strategies have been developed for parasitic weeds, including improved cultural practices, breeding using wild and cultivated germplasm as sources of resistance, and the use of chemical control, but the use of resistant cultivars is still considered the most effective (Parker and Riches, 1993; Singh et al., 2006; Li et al., 2009). Given its significance in cowpea production, S. gesnerioides resistance is required for varieties released for WA’s Sahelian and Sudan Savanna zones. Hence, Striga resistance is a “must-have trait” in the cowpea product profiles. Consequently, significant efforts have been made to study the Striga race structure and identify resistance sources in cowpea. Molecular profiling and host differential response studies using S. gesnerioides isolates across the WA region revealed seven distinct parasite races (Botanga and Timko, 2006). An additional report by (Li et al., 2009; Ohlson and Timko, 2020) confirmed the seven races of S. gesnerioides. The authors observed that some cowpea cultivars were differentially resistant to various geographic isolates of the parasite. The races have been designated as SG1 (Burkina Faso), SG2 (Mali), SG3 (Nigeria and Niger), SG4 (Benin), SG4z (localized to the Zakpota region of Benin), SG5 (Cameroon), and SG6 (Senegal) (Li et al., 2009). However, Ohlson and Timko (2020) reported that SG6 from Senegal was related to the SG1 race, and the authors revealed a novel race (designated as SG6) in Nigeria that was able to overcome more cowpea resistance genes than any previously reported race.

Screening efforts have identified sources of resistance to S. gesnerioides. These include B301 (Botswana landrace), IT82D-849, IT81D-994, and Wango-1 (from Burkina Faso), all having monogenic dominant inheritance mode and are resistant to races SG1, SG2, and SG3 (Singh and Emechebe, 1990; Atokple et al., 1993; Li et al., 2009). Genotype HTR (from Niger) was reported to carry one or two dominant genes and is resistant to race SG1, while Suvita-2 has a single recessive gene against SG3 and a single dominant gene against SG1 and SG2 (Li et al., 2009). In addition, Omoigui et al. (2010) screened some cowpea genotypes and found B301, IT03K-338-1, and IT99K-573-2–1 to be free of emerged Striga and Alectra shoots, while IT98K-1092–1 and IT97K-205–8 were resistant to Striga but supported the emergence of some Alectra shoots.

Past efforts also identified some molecular markers associated with resistance to S. gesnerioides in cowpea. Notably, Ouédraogo et al. (2001), (2002); Li et al. (2009) reported amplified fragment length polymorphism (AFLP) markers on linkage group 1 (LG1) linked to race-specific genes: Rsg2–1 in the cowpea line IT82D-849, Rsg1–1 in B301, and Rsg4–3 in Tvu14676. Additionally, race-specific genes Rsg3–1 and Rsg994-1, present in Suvita-2 and IT91D-994 respectively, were mapped on LG6 (Ouédraogo et al., 2001; Li et al., 2009). Two sequence-characterized amplified regions (SCARs) markers 61R (E-ACT/M-CAA) and SEACTMCAC83/8, linked to S. gesnerioides resistance, were developed and deployed for marker-assisted selection (Ouédraogo et al., 2001; Li et al., 2009; Ouédraogo et al., 2012). Four AFLP markers, E-ACT/M-CTC115, E-ACT/M-CAC115, E-ACA/M-CAG108, and E-AAG/E-CTA190, were found to be associated with the Rsg1 gene in a resistant line IT93K- 693–2 that confers resistance to SG3 race (Boukar et al., 2004). Some Simple Sequence Repeat (SSR) markers associated with resistance to S. gesnerioides race 3 (SG3) were identified and deployed in breeding for resistance (Li et al., 2009; Asare et al., 2013; Omoigui et al., 2017; Essem et al., 2019). However, these old marker technologies have known limitations for wide-scale routine applications in breeding, including difficulty in handling and scoring, and a lack of automation to allow high-throughput genotyping. Consequently, there has been a significant shift toward using SNPs, given their genomic abundance and rapid emergence of novel, faster, and cheaper methods of genotyping (Tsuchihashi and Dracopoli, 2002).

Recent progress in developing next-generation genomics and genetic resources for cowpea has generated more robust molecular marker platforms than those described earlier. Notably, Muchero et al. (2009) developed a cowpea GoldenGate assay consisting of 1536 SNP markers, being utilized for linkage mapping and QTL analyses (Lucas et al., 2011; Muchero et al., 2013; Pottorff et al., 2014) and assessment of genetic diversity (Huynh et al., 2013). Illumina Cowpea iSelect Consortium Array having 51,128 SNPs, was also developed (Muñoz-Amatriaín et al., 2017) and has been deployed in genome-wide mapping of several traits in cowpea (Herniter et al., 2018, 2019; Lo et al., 2018, 2019; Miesho et al., 2019; Paudel et al., 2021; Ongom et al., 2022). Recently, cowpea researchers have developed low-density Kompetitive Allele Specific PCR (KASP) assays (Ongom et al., 2021; Wu et al., 2021) and medium-density DArTag genotyping panel (Ongom et al., 2024, 2024). Several cowpea genetic resources have also been developed; among them are the UCR minicore, a diverse set of genotypes that have been genotyped with the iSelect SNPs panel (Muñoz-Amatriaín et al., 2021), and the IITA minicore, a set of genotypes representing the cowpea germplasm maintained at the IITA Genetic Resources Center that have been genotyped based on genotyping by sequencing (GBS) (Fatokun et al., 2018). The cowpea genomic and genetic resources so far developed have opened doors for QTL/gene discovery and development of robust molecular markers to enhance breeding for essential traits, including Striga resistance.

Resistance to Striga gesnerioides is notoriously difficult to evaluate in the field due to a complex interplay of confounding factors: seasonal shifts, soil heterogeneity, subjective scoring criteria, and significant parasite race diversity (Li et al., 2009). These challenges highlight a critical need for a reliable marker-assisted selection (MAS) platform, especially given the prevalence of multiple S. gesnerioides races that may circumvent resistance derived from a single source. In response, this study aimed to: (i) identify SNPs robustly associated with Striga resistance; and (ii) develop Kompetitive Allele Specific PCR (KASP) assays suitable for routine use in cowpea breeding. The study harnessed the diversity in the UCR minicore panel, combining two years of Striga resistance phenotypic data with the high-density iSelect SNPs to pinpoint genomic regions associated with S. gesnerioides in cowpea. Ultimately, the development and validation of KASP markers linked to Striga resistance will advance the integration of MAS into breeding pipelines. These markers offer breeders a precise, cost-effective tool for selecting Striga resistance, expediting the development of resilient cowpea varieties and accelerating their deployment in farmers’ fields.

2 Materials and methods

2.1 Genetic materials

The study used 368 UCR minicore genotypes described by (Muñoz-Amatriaín et al., 2021). The UCR minicore contains diverse landraces and breeding materials from 50 countries covering Africa, Asia, North and South America, Europe, and Australia. The minicore population, previously genotyped with high-density SNPs, represents the existing genetic and phenotypic diversity of cultivated cowpea while maintaining a sample size that can be managed by most researchers and breeders for evaluating traits of interest (Muñoz-Amatriaín et al., 2021). For marker development, the present study used 60 diverse cowpea genotypes along with 46 F1 progenies for SNPs technical verification and another independent set of 20 genotypes consisting of some known resistant and susceptible sources commonly used as standard checks in breeding programs for marker validation.

2.2 Study sites

The study was conducted at three sites in Northern Nigeria where there is a widespread of Striga gesnerioides. The first site was a Striga hot spot at IITA Minjibir Research Farm, Kano, Latitude 12° 14′ 35.30″ N and Longitude 8° 66′ 62.10″ E. The second Striga hot spot site was at Malam Madori, Jigawa, located at Latitude 12° 33’ 36.32” N and Longitude 9° 59’ 9.56” E. Previously, the samples of Striga races collected from these two sites have been characterized as belonging to race SG3, which is a dominant S. gesnerioides race in Nigeria among other less frequent races labeled as SG1, SG5 and SG6 (Ohlson and Timko, 2020). The third site designated for validation screening is situated at the IITA station in Kano, located at latitude 11°58’50.0”N and longitude 8°33’26.8”E. This facility includes a nursery that has been deliberately infested with Striga gesnerioides race SG3 to facilitate controlled research on cowpea resistance to this parasitic weed.

2.3 Striga resistance phenotyping

The minicore genotypes were evaluated in the two locations in Northern Nigeria for two years in Minjibir and one year in Malam Madori. The genotypes were planted in a two-row plot of 4 meters long and at a spacing of 0.75m between the rows and 0.2m within the row. The trials were established as alpha lattice designs with two replications. Variation in Striga infestation was assessed using two criteria: First, the Striga score was evaluated on a scale of 0-3, where 0 = no Striga emergence; 1 = few Striga emerged; 2 = moderate Striga emergence; 3 = heavy Striga infestation The second method was the presence-absence rating, where ‘0’ indicates that Striga is absent in the plot and ‘1’ indicates Striga is present. In addition, a validation set of 20 cowpea genotypes was screened in an artificially infested nursery using a randomized complete block design experiment. For the past several years, this nursery has been dedicated to Striga screening following heavy inoculation of the soil with seeds of Striga collected from infested fields at Malam Madori, Jigawa State, Nigeria. Prior to soil inoculation, Striga seeds were preconditioned through surface-sterilization using 10% sodium hypochlorite for ~10 minutes followed by incubation in the dark at 29°C for 10 days. The soil was then inoculated with Striga seeds at a rate of 2 g/m2, evenly distributing the seeds across the soil surface. The 20 cowpea genotypes were each planted in a 2-row plot of 4m and spaced at 0.75m by 0.2m between and within rows, respectively, and the experiment had two replications. The presence or absence of Striga was assessed for each plot, including the number of Striga plants that emerged, and the data was analyzed with Analysis of Variance (ANOVA) and means compared using LSD test.

2.4 SNP genotype data

The SNP data was obtained from (Muñoz-Amatriaín et al., 2017). The genotyping was done at the University of Southern California Molecular Genomics Core lab (Los Angeles, California, USA) following the procedure described by (Muñoz-Amatriaín et al., 2021). Briefly, genomic DNA was extracted from the young leaves of individual plants using DNeasy Plant Kit (Qiagen, Valencia, California, USA). The Cowpea iSelect Consortium Array including 51,128 SNPs (Muñoz-Amatriaín et al., 2017) was used for genotyping each DNA sample. SNPs were called in GenomeStudio (Illumina Inc., San Diego, California, USA) and manually curated to remove those with high levels (>20%) of missing data and/or heterozygous calls. The genotype data contained 47129 SNPs after removing the contigs with no chromosomal assignment. Further filtering was conducted, setting the minimum minor allele frequency at 0.05 and maximum heterozygous proportion at 0.1, resulting in 41510 SNPs for downstream analyses.

2.5 Data analysis

The phenotypic data were analyzed using a general linear model tailored for an alpha lattice design, implemented through the agricolae package in R (de Mendiburu, 2020). The model applied is represented as:

Y_ijk=μ+τ_i+γ_j+ρ_k(j) +ϵ_ijk

Where: Y_ijk denotes the observed phenotypic value for the ith treatment in the kth block within the jth replication., μ is the general mean, τ_i represents the fixed effect of the ith treatment, γ_j signifies the effect of the jth replication, and ρ_k(j)  accounts for the effect of the kth incomplete block nested within the jth replication, ϵ_ijk is the random error term associated with the observation. The mean values for each minicore accession were extracted from the model using the LSD.test () function and later utilized in GWAS analysis.

GWAS was conducted using the 41510 filtered SNPs in TASSEL v 5.2.20 (Bradbury et al., 2007) and rMVP package v 1.03 (Yin et al., 2021), utilizing the Mixed Linear Model (MLM) (Zhang et al., 2010) and the Fixed and random model Circulating Probability Unification (FarmCPU) (Liu et al., 2016) on the following traits: Striga score and Striga presence-absence rating. The FarmCPU model uses a multiple loci linear mixed model (MLMM) and incorporates multiple markers simultaneously as covariates in a stepwise MLM to partially remove the confounding between testing markers and kinship (Liu et al., 2016). A genomic PCA matrix (P) and kinship matrix (K) were used to capture the population structure and relatedness among individuals in the panel (Kang et al., 2008). In TASSEL, kinship, and principal components were computed and fitted together with the SNPs in the MLM model to account for relatedness and population structure, respectively. The MLM statistics were also utilized separately in CMplot package 3.1.3 (Yin, 2016) to generate customized GWAS Manhattan and QQ plots. Decisions on significant GWAS signals were based on both the conservative Bonferroni (Henry, 2015) and a less conservative false discovery rate (FDR) (Glickman et al., 2014; Ongom et al., 2022, 2024) corrections of multiple statistical tests to reduce the risk of a type I error.

2.6 Candidate gene search

Candidate genes linked to Striga resistance were pinpointed by aligning the positions of significant SNPs from the GWAS with the cowpea reference genome (version 1.1) available on Phytozome 13 Genome Browser (https://phytozome-next.jgi.doe.gov/info/Vunguiculata_v1_1, accessed on 26 January 2025). A length of 200 kb was added or removed from either end of the significant marker to locate potential regions for comparison based on the LD rate of the cowpea minicore population (Ongom et al., 2022). In selecting candidate genes, the following criteria were used: (i) genes of known function in cowpea related to the trait under study, (ii) genes with function-known orthologs in Arabidopsis related to the trait under study, and (iii) genes pinpointed by the peak SNPs. Putative candidate genes were subsequently researched in the literature for verification. We then used chromoMap v4.1.1 package (Anand and Rodriguez Lopez, 2022) to visualize the distribution of selected genes on target chromosomes where GWAS signals were discovered.

2.7 KASP marker development

Technical verification of KASP assay: Seventeen SNP markers spanning Striga resistance gene regions were selected to develop KASP assays. The design sequences of the selected SNPs are presented in Table 1. The validity of these KASP assays was verified using 60 cowpea genotypes and 46 heterozygous F1 progenies. The samples for this SNP technical verification were collected in three replicates per genotype. The leaf samples and SNP design sequences were sent to Intertek Lab Sweden for assay development. Leaf sampling followed the standard procedure previously described (Ongom et al., 2021). Genotyping data were visualized using SNPviewer version 5.1.1.27582 (LGC Genomics, 2025) and SNP cluster plots were generated.

Table 1
www.frontiersin.org

Table 1. Candidate SNP markers for Striga resistance used to develop the KASP assay.

Testing the marker-trait association: To validate the performance of these markers for tracking Striga resistance, we deployed the candidate markers in a panel of 20 cowpea validation genotypes consisting of some known resistant and susceptible sources. These cowpea genotypes were independently phenotyped in a an artificially Striga-infested nursery to confirm their reaction status. The number of Striga plants that emerged in each plot was counted, and the data was used to categorize the 20 genotypes into either resistant (R) or susceptible (S) with zero Striga count = R and any count greater than zero = S. The 20 cowpea genotypes were then genotyped with the 17 candidate KASP markers. A single marker analysis was then conducted to verify the consistency of these candidate markers in tracking Striga resistance in cowpea. A chi-square test of independence and a t-test were performed in R for each marker with the null hypothesis that each marker genotype is independent of the observed variation in the phenotypic status of the cowpea genotypes. The results of these tests were visualized using plots generated by ggplot2 and ggstatsplot packages.

3 Results

3.1 Phenotypic assessment

Phenotypic assessment

The phenotypic data from both Minjibir and Malam Madori sites depicted skewed distributions for Striga ratings, with a greater proportion of the minicore genotypes showing susceptibility to Striga. This is portrayed by the stacked bar chart in Figure 1A, where only 2.7 to 6.5% of the minicore showed resistance across the two testing sites. In contrast, a mean Striga score of 2.35 (scale 1-3) was registered for the combined data across the two locations (Figure 1B). The susceptible genotypes triggered extensive growth of Striga plants, causing severe yellowing and death of cowpea plants. At the same time, no emergence of Striga was observed in plots planted with resistant genotypes (Figure 1C).

Figure 1
Panel A shows a bar chart comparing the presence and absence of Striga across three environments, with high proportions present (red). Panel B illustrates a histogram of density distribution with a mean of 2.35 and standard deviation of 0.35. Panel C consists of four images: (i) displays resistant and susceptible cowpea plants; (ii) shows Striga-infested plants; (iii) depicts yellow and wilted cowpea plants; (iv) features a dead cowpea plant with a label.

Figure 1. Density curve displaying the phenotypic distribution and reaction of minicore genotypes to Striga gesnerioides. (A) the distribution for presence-absence rating at Minjibir with the x-axis indicating a Striga rating of 0 (Striga absent) and 1 (Striga present), (B) the distribution for presence-absence rating at Malam Madori with the x-axis indicating Striga rating of 0 (Striga absent) and 1 (Striga present) (C) the distribution for Striga score (on the x-axis) for combined data across the two locations, (D) the response of minicore genotypes to Striga infestation in the field, with panel (i) depicting resistant and susceptible genotypes and panels (ii) to (iv) displaying the phases of plant growth under Striga infestation from green healthy-looking but infested plants through severe yellowing and wilting up to completely dead cowpea plants.

3.2 Association analysis

GWAS discovered one major association signal for Striga resistance on chromosome Vu11 that was consistently significant in both locations and five other minor signals on chromosomes Vu02, Vu03, Vu07, Vu09, and Vu10 (Figure 2A). The major association signal on chromosome Vu11 and two minor regions on chromosomes Vu07 and Vu10 were detectable in both Minjibir and Malam Madori locations, while the regions on Vu03, Vu02, and Vu09 were significant only in one location (Figure 2A). The QQ plot shows that the observed p-values largely align with the expected distribution under the null hypothesis, indicating no evidence of inflation. However, the upward deflection at the tail suggests the presence of a small number of SNPs exhibiting stronger association signals than expected by chance.” (Figure 2B). It was also evident from the QQ plot that the GWAS signals were stronger in Malam Madori than in Minjibir and the combined location data set for Striga scores (Figure 2B).

Figure 2
Manhattan plot and Q-Q plot for genetic association study data. The Manhattan plot (A) displays −log₁₀(p-value) across chromosomes Vu01 to Vu11, showing multiple colored points denoting different datasets. Horizontal lines indicate p-value thresholds. The Q-Q plot (B) compares observed versus expected −log₁₀(p-values), with colors corresponding to datasets, demonstrating deviations from the null hypothesis.

Figure 2. Genome-wide association signals based on the mixed linear model (MLM). (A) Manhattan plot depicting a major GWAS signal on chromosome Vu11 and five minor ones on chromosomes Vu02, Vu03, Vu07, Vu09 and Vu10. The data is presented for Striga presence-absence rating in Minjibir (stgrpaMJB_Y1.MLM and stgrpaMJB_Y2.MLM) and Malam Madori (stgrpaMD_Y1.MLM), Striga score in Malam Madori (stgrscoreMD_Y1), Striga score in Minjibir (stgrscoreMJB_Y1.MLM) and Striga score for data combined across the two locations (atgrscore_COMB.MLM). Beneath each chromosome in the Manhattan plot is the SNP density and distribution defined by the colored key on the top right, with gray color representing low SNP density while red depicts high SNP density regions. (B) the QQ plot for the six experimental data sets.

A total of 309 significant SNPs with the Boneforrini threshold of -log10(p) > 5.92 were found to be associated with Striga resistance (Supplementary Table 1). Out of these, seventeen (17) representative SNPs that had R2 values above 7.5% were selected for KASP assay development (Table 2). The strongest GWAS signal on chromosome Vu11 was represented by five significant SNPs, accounting for 15.8% to 19.8% of phenotypic variance (Table 2). The other minor signals on chromosomes Vu02, Vu03, Vu07, Vu09, and Vu10 explained 6.8% to 14.1% of phenotypic variance (Supplementary Table 1).

Table 2
www.frontiersin.org

Table 2. Representative SNPs that were significantly associated with Striga resistance in the cowpea minicore population.

3.3 Candidate genes

Searching candidate genes within 200kb of peak GWAS signals identified 178 genes on four target chromosomes, with 28 on Vu11, 38 on Vu10, 56 on Vu07, and 57 on Vu02 (Supplementary Table 2). Sixty-four of the 178 genes were clustered closely to the peak SNPs on the target chromosomes (Figure 3). Further examinations revealed 20 unique annotated proteins associated with the 64 genes (Figure 3), given that most of these gene loci encoded similar functional proteins. Of the 64 candidate genes identified, 11 were positioned within approximately -167461 bp to 191285 bp of the lead SNP on Chromosome Vu11, five genes were within ~ -139927 bp to 83997 bp of the association peak on Vu10, 36 genes were within ~ -190061 bp to 202188 bp of the signal on Vu07 and 12 were within ~ -145716 bp to 125678 bp of the association locus on Vu02 (Table 3). Several of these clusters contain genes with related biological functions. For instance, up to five genes on Chromosome Vu02 belonged to the Pentatricopeptide repeat (PPR) superfamily protein. On chromosome Vu07, up to 28 genes had functions related to Cysteine-rich receptor-kinase-like protein, while 10 genes on chromosome Vu11 had functions related to LEUCINE-RICH REPEAT-CONTAINING PROTEIN (Table 3). The rest of the genes on each chromosome had unique functions, including, among others, C2H2-like zinc finger protein, myb transcription factor, MADS-box transcription factor family protein, ETHYLENE-RESPONSIVE TRANSCRIPTION FACTOR ERF003, and WRKY DNA -binding domain (WRKY) (Table 3, Figure 3).

Table 3
www.frontiersin.org

Table 3. Selected genes found within 200kb of the representative GWAS SNPs with implicated functions in plant immunity and defense signaling.

Figure 3
Four horizontal black bars labeled Vu02, Vu07, Vu10, and Vu11, each representing a chromosome segment. Small colored blocks within the bars indicate gene locations. A color key on the right identifies gene types, including leucine-rich repeat and cytochrome P450 protein. Horizontal axis shows ruler markers from zero to forty megabase pairs.

Figure 3. Heatmap of the position of key genes on chromosomes. It depicts the selected proximal genes to the peak SNPs with likely functional involvement in plant defense and immune signaling on chromosomes Vu02, Vu07, Vu10, and Vu11. The chromosomal positions of the genes are marked in green, while a cluster plot of genes is displayed above the marked regions. The color-coded legend on the right identifies 20 unique proteins associated with genes found within the QTL signal region.

3.4 KASP marker assay development

KASP assays were successfully developed for the 17 representative SNPs having strong associations with Striga resistance in cowpea. The SNP calls were verified using a unique technical validation set of cowpea genotypes, including homozygous genotypes and highly heterozygous F1. The raw genotype calls displaying the KASP assay performance of the 17 candidate SNPs in the cowpea technical validation set are provided in Supplementary Table 3. The KASP assay technical verification results revealed 12 SNPs as having good quality assays, allowing easy scoring of the alleles (Figure 4A). Two of the SNPs were of medium quality, with some ambiguity in differentiating between homozygotes and heterozygotes (Figure 4B). Two other SNPs were rated as having inconclusive quality due to their sensitivity to DNA concentration and the resultant difficulty in scoring the marker genotypes (Figure 4C), while one SNP did not form scorable clusters and was regarded as having a bad quality (Figure 4D). A detailed report on the quality assessment of each SNP is presented in Table 4. SNPs that were easy to score formed three genotype clusters representing the two Mendelian homozygotes and one heterozygote. Those with inconclusive results formed either one or two genotype clusters, while the bad-quality SNPs did not create any meaningful clusters and were considered failed SNPs (Table 4). Overall, the performance of 14 out of the 17 SNPs was deemed acceptable and selected as candidates for further validation.

Figure 4
Four scatterplots of SNP genotyping data from Biosearch Technologies show different SNP markers. Each plot (A-D) displays clusters of colored dots representing genotypes with axes labeled. The grids below detail additional sample information.

Figure 4. Visualizations of KASP assays technical verification of the candidate SNPs. (A) represents the 12 SNPs that had good quality KASP assays, (B) represents the 2 SNPs that had medium KASP assay quality and were still scorable, (C) represents the 2 SNPs that had inclusive results (only one or two genotype clusters), (D) the SNP that could not be scored since no clear genotype clusters were formed. In this figure, the blue and red clusters are the two Mendelian homozygous genotypes, while the green cluster is the heterozygous genotype. The black data points represent the no template controls (NTC). The pink data points marked with “?” represent samples that did not generate consistent signals or failed to amplify. The DNA sample plate layout is shown below the cluster plot.

Table 4
www.frontiersin.org

Table 4. KASP assay technical verification report for 17 candidate SNPs for Striga resistance.

3.5 KASP marker validation

The validation exercise involved screening 20 cowpea genotypes in an artificially Striga-infested nursery and comparing the Striga phenotypic data with the marker genotypes on these cowpea genotypes. The phenotypic distribution of the 20 cowpea genotypes in response to Striga infestation is presented in Figures 5A, B. The resistant cowpea genotypes had zero Striga emergence, as portrayed by zero mean and median. In contrast, the susceptible cowpea genotypes showed dispersed distributions with the mean and median Striga count clearly above zero (Figures 5A, B). Consequently, the counts of Striga emergence as influenced by the 20 cowpea genotypes were significantly different, with each genotype recording varying mean Striga count (Table 5). Five cowpea genotypes were entirely immune to Striga (no Striga emerged in these plots), while the other genotypes differed in the average number of Striga ranging from 3-10.5, indicating different response levels to Striga (Table 5). As expected, an old variety named Achishiru registered the highest number of Striga, followed by the landraces TVu-867 and Vita7, among the susceptible genotypes compared to low Striga counts registered in some of the released varieties like IT07K-318-33, IT98K-1111-1, IT08K-150-12, IT89KD-288 and IT98K-589-2 (Table 4), suggesting that these are Striga tolerant genotypes.

Figure 5
Panel A shows a bar chart comparing phenotype R and S by Striga count, with R having a mean of zero from five samples, and S having a mean of four point eight three from fifteen samples. Panel B displays a violin plot comparing the same phenotypes, with R and S showing different distributions in Striga count.

Figure 5. Distribution of Striga counts for the 20 cowpea genotypes evaluated in an artificially Striga-infested nursery at IITA Kano, Nigeria, grouped by their phenotypic status. (A) histogram showing the distribution of the resistant genotypes (R) in red and susceptible genotypes (S) in blue. The red and blue vertical dash lines depict the mean of the resistant and susceptible groups, respectively. (B) boxplot portraying the resistant genotypes (R) dispersion in red and susceptible genotypes (S) in blue. The black line at the center of the boxplot indicates the median Striga count.

Table 5
www.frontiersin.org

Table 5. Reactions of 20 selected cowpea genotypes to Striga infestation, evaluated in an artificially infested Striga nursery at IITA Kano station, Nigeria.

Using a chi-square test of independence, the study compared the association between each of the 14 candidate SNPs with the phenotypic status of the 20 cowpea genotypes (Table 6). Five SNPs in close proximity to the major association signal on chromosome Vu11 displayed highly significant deviations (P ≤ 0.01) from the null expectation of no marker–phenotype correlation, supporting their strong association with genomic regions influencing Striga resistance in cowpea (Table 6). Five other SNPs positioned close to each other on chromosome Vu02 showed moderate statistical significance (P ≤ 0.05). The rest of the SNPs on Chromosomes Vu07 and Vu10 revealed no significant chi-square values (Table 6).

Table 6
www.frontiersin.org

Table 6. Chi-square test of independence between the candidate SNPs and the Striga resistance in cowpea.

Further efforts to group the phenotypic reactions by marker alleles revealed a strong marker-phenotype association for the candidate markers on chromosome Vu11 and a weak relation for the other SNPs (Figure 6, Supplementary Figures 1A–D). For instance, the GG allele of the marker snpVU00077 on Vu11 was responsible for 100% of the cowpea genotypes that had a Striga-resistant phenotype, while the alternative allele CC was carried by 94% of phenotypically susceptible genotypes, with only 6% phenotypic misclassification (Figure 6). This observation was supported by a highly significant chi-square test and high Cramer’s correlation (V^Cramer =0.86). Similar results were observed for all five proximal SNPs on Vu11 (Supplementary Figure 1D), suggesting that these markers are strongly associated with genes underlying Striga resistance in cowpea. However, a relatively high percentage of misclassifications and V^Cramer less than 0.5 were observed for the remaining candidate markers on chromosomes Vu02, Vu07, and Vu10 (Figure 6, Supplementary Figures 1A–C), suggesting weak marker-trait association.

Figure 6
Four bar charts titled “Chromosome Vu02,” “Chromosome Vu07,” “Chromosome Vu10,” and “Chromosome Vu11” compare phenotypes S and R across different SNPs. Each chart presents two bars showing varying percentages of phenotypes S (green) and R (orange) with sample sizes below each bar. Statistical data including chi-square and p-values is noted above each chart.

Figure 6. Stacked bar chart depicting the grouping of phenotypic reactions to Striga by marker alleles. Four markers: snpVU00064, snpVU00070, snpVU00074, and snpVU0007, represent the candidate regions on chromosomes Vu02, Vu07, Vu10, and Vu11, respectively. Results for all 14 markers are presented in Supplementary Figure 1. The 20 cowpea genotypes are grouped based on their phenotypic reaction to Striga into resistant (R) and susceptible (S) classes, and the stacked bar charts are plotted to reflect the percentage of genotypes in each phenotypic category that carry the two alleles of each SNP marker. Summary statistics are presented at the top of each stacked bar chart that includes the Pearson chi-square test (χ2Pearson), probability (p) presented for each marker allele class, and for the overall chi-square test of independence, effect size measured by Cramer’s correlation (V^Cramer), 5% confidence interval (CI95%) and number of observations (nobs).

Similarly, a t-test was conducted to compare the differences in the mean Striga count between cowpea genotypes carrying the two alleles of each marker. Results for four markers representing the candidates on chromosomes Vu02, Vu07, Vu10, and Vu11 are presented in Figure 7. Highly significant differences (p 0.001) between the mean of the allelic groups were found only for the markers on chromosome Vu11, while the other markers were not significant, strongly supporting the chi-square test results. The marker alleles conferring resistance to Striga had significantly lower mean Striga counts than the alternative alleles. The results of the t-test for all 14 markers specifying the alleles underlying resistance and susceptibility and how they differ have been presented in Supplementary Figures 2A–D.

Figure 7
Violin plots display Striga count data across genetic variants for chromosomes Vu02, Vu07, Vu10, and Vu11. Each plot includes statistical information: t Welch test, p-value, g Hedges, confidence intervals, and sample size. Significant differences and mean values are marked.

Figure 7. The violin chart depicting the difference in mean Striga counts between cowpea genotypes carrying the two alleles of each candidate marker. Results of four markers, snpVU00064, snpVU00070, snpVU00074, and snpVU0007, are presented to represent the candidate regions on chromosomes Vu02, Vu07, Vu10, and Vu11, respectively. Results for all 14 markers are presented in Supplementary Figure 2. The mean (μ^mean) of each allelic group is indicated by a red dot. Summary statistics are presented at the top of the violin chart, including t-statistical test (twelch), probability (p) for the t-test, effect size measured by Hedges’ g (ĝHedges), a 5% confidence interval (CI95%) and number of observations (nobs).

4 Discussion

Striga gesnerioides, commonly known as cowpea witchweed, is a parasitic plant that poses a significant threat to cowpea (Vigna unguiculata) production, particularly in West and Central Africa. This parasitic weed attaches to cowpea roots, extracting essential nutrients and water, leading to substantial yield reductions and, in severe infestations, complete crop failure (Li et al., 2009; Sadda et al., 2021).​ The multiplication of S. gesnerioides is facilitated by its prolific seed production, with a single capsule containing hundreds of seeds that can remain viable in the soil for several years (Mohamed et al., 2001). Addressing the challenges posed by S. gesnerioides requires integrated management strategies, among which developing and cultivating Striga-resistant cowpea varieties is pivotal (Abdullahi et al., 2022). The application of molecular markers has been emphasized as the best approach to facilitate breeding for resistance to this parasitic weed, given its race complexity (Ouédraogo et al., 2001; Boukar et al., 2004; Li et al., 2009; Ouédraogo et al., 2012; Essem et al., 2019). Successfully deploying markers in breeding requires discovering and validating markers that closely tag the resistance loci. In the case of S. gesnerioides, which has multiple races (Botanga and Timko, 2006; Ohlson and Timko, 2020), it is vital to develop molecular markers linked to resistance against each race. These markers would facilitate the pyramiding of resistance genes, hence the development of cowpea varieties suitable for cultivation across the Striga-affected area of WA. The present study significantly contributes to these efforts by uncovering key loci and developing molecular markers underlying resistance to Nigeria’s dominant Striga race SG3.

4.1 Genomic regions associated with Striga resistance

The present study deciphered genomic regions controlling S. gesnerioides resistance in cowpea by utilizing a high-density single nucleotide polymorphism marker and diverse cowpea minicore genotypes. The phenotypic data from the two test sites revealed significant genetic variation in resistance to Striga gesnerioides with skewed distribution, indicating that major-effect genes likely govern resistance. The finding aligns with prior studies, which suggested that resistance to S. gesnerioides in cowpea is often controlled by specific resistance (R) genes exhibiting major effects. Notably, the cowpea genotypes B301 and IT82D-849 exhibited resistance controlled by a single dominant gene effective against Striga races SG1, SG2, and SG3 found in Nigeria (Singh and Emechebe, 1990; Atokple et al., 1993; Li et al., 2009). In a cross of cowpea genotypes HTR (from Niger) and Wango-1 (from Burkina Faso), it was reported that resistance to S. gesnerioides race SG1 in HTR was controlled by one or two dominant genes that are non-allelic to the genes in B301 and IT82D-849 (Timko et al., 2007; Li et al., 2009). In addition, monogenic but recessive inheritance of resistance to S. gesnerioides race SG3 from Niger was also reported (Touré et al., 1997). However, the interaction between cowpea and S. gesnerioides is complex due to the existence of multiple parasite races. Seven distinct races (SG1 to SG7) were identified, each capable of overcoming specific resistance genes in cowpea (Li et al., 2009; Ohlson and Timko, 2020). This variability necessitates discovering key loci involved, followed by stacking multiple resistance genes to achieve broad-spectrum and durable resistance.

A genome-wide association analysis revealed a major association signal on chromosome Vu11, alongside f significant minor signals on chromosomes Vu02, Vu07, and Vu10 linked to Striga resistance. The prominent Vu11 signal aligns closely with the a coiled-coil nucleotide-binding site leucine-rich repeat (CC-NBS-LRR) protein, encoded by the RSG3–301 gene, was implicated in race-specific resistance to S. gesnerioides in cowpea (Li and Timko, 2009). This proximity suggests that similar R-genes may underlie the Vu11 association signal.​ Additionally, a BLAST analysis using primers from earlier marker technologies (e.g., AFLP, SSR, and SCAR) returned best matches within the same Vu11 region (Table 7), further supporting the colocalization of historical markers with our GWAS-identified signals. The previously reported genes spanning the same area on chromosome Vu11 included Rsg3 (conferring resistance to Striga race SG3 and SG5), Rsg2–1 (effective against Striga race SG1), and Rsg4–3 (effective against Striga race SG3) (Ouédraogo et al., 2001; Essem et al., 2019). Earlier linkage mapping studies identified resistance loci to Striga in cowpea via bi-parental populations. Ouédraogo et al. (2002) mapped a race-specific resistance gene (designated Rsg3 and Rsg994) to linkage group 1 (equivalent to chromosome Vu10 in the reference genome), conferring resistance to Striga race SG1. Boukar et al. (2004) mapped another resistance locus, Rsg1, to linkage group 1 (chromosome Vu08), conferring race-specific resistance to SG3 The studies above used traditional linkage mapping in F2 bi-parental populations to identify the Striga resistance loci. In comparing these classical mapping results with our GWAS findings, we observed a strong correspondence between previously reported loci and our newly identified association signals (see Table 7). This concordance supports the involvement of these genomic regions, now validated by independent GWAS in mediating Striga resistance.

Table 7
www.frontiersin.org

Table 7. Previously mapped Striga resistance QTL in cowpea with some of the loci occupying the same genomic region discovered in the present GWAS.

4.2 Candidate genes for Striga resistance

Through our candidate gene analysis, we identified 20 unique proteins that have been implicated in plant defense and immune response signaling, collectively annotated to 64 cowpea genes spanning the GWAS-identified association intervals. Notably, among these is a leucine-rich repeat (LRR)-containing protein mapped near the major association peak on chromosome Vu11.LRR domains are well-established as pivotal in pathogen recognition and activation of plant defense responses, as seen in both extracellular receptor kinases and intracellular NB-LRR immune receptors (Jones and Jones, 1997; McHale et al., 2006). An earlier study in cowpea pinned leucine-rich repeat (CC-NBS-LRR) protein to a gene-for-gene resistance mechanism in the interactions between Striga and cowpea, with a corresponding gene in cowpea named RSG3-301 (Li and Timko, 2009). The authors also unlocked the hypothesis of race specificity resistance by silencing the RSG3–301 gene in cultivar B301, which rendered it susceptible to race RG3 but remained resistant to races SG2 and SG5. Another gene PTHR22952:SF183 - TRANSCRIPTION FACTOR TGA5 on Vu11, has been found to work along with TGA2 and TGA6, playing a crucial role in plant systemic acquired resistance (SAR), a broad-spectrum defense mechanism induced after a local infection by avirulent pathogens (Zhang et al., 2003). On chromosome Vu10, the likely genes involved in plant defense were ETHYLENE-RESPONSIVE TRANSCRIPTION FACTOR ERF003, WRKY DNA-binding domain, and PF03763 - Remorin, C-terminal region (Remorin_C). Ethylene-responsive transcription factors (ERFs) are integral components of the APETALA2/ERF superfamily, playing pivotal roles in regulating plant responses to biotic and abiotic stresses (Müller and Munné-Bosch, 2015). ERFs regulate molecular response to pathogen attack by binding to specific cis-acting elements in the promoters of stress-responsive genes, thereby modulating their expression (Müller and Munné-Bosch, 2015). WRKY transcription factors, on the other hand, have been associated with responses to various pathogens in cowpea, particularly exhibiting differential expressions when the plants were challenged with Fusarium oxysporum, a pathogen responsible for Fusarium wilt (Hao et al., 2024). This finding underscores the significant role of WRKY genes in cowpea’s defense mechanisms against biotic stresses. Given the conserved nature of WRKY transcription factors across plant species, it is plausible that manipulating specific WRKY genes in cowpea could enhance resistance to Striga. However, targeted studies are required to identify which WRKY genes are involved and to elucidate their mechanisms of action in the context of Striga resistance. In the context of plant defense, remorins have been observed to interact with receptor-like kinases (RLKs) and pathogen effectors, suggesting a role in the early stages of immune signaling (Yu, 2020). Chromosome Vu07 harbored genes encoding MADS-box transcription factors, cysteine-rich receptor-like kinases, MYB-like proteins, and GATA-type zinc finger transcription factors, all of which play roles in plant development and stress responses (Abdullah-Zawawi et al., 2021; Zhang et al., 2023a, b). The genes on chromosome Vu02, included, among others, pentatricopeptide repeat (PPR) proteins, C2H2-like zinc finger proteins, glucuronosyltransferase PGSIP8-like, receptor-like kinases, and lipid transfer proteins that contribute to the complex network of plant defense mechanisms, each playing distinct roles in responding to biotic and abiotic stresses (Goff and Ramonell, 2007; Finkina et al., 2016; Zhang et al., 2020; Huang et al., 2022; Meng et al., 2024).

4.3 Marker development

Seventeen markers tagging GWAS-identified association intervals for Striga resistance were converted into KASP assays and independently validated across 20 cowpea genotypes with known resistance profiles. Five of these markers, all situated near the major association signal on chromosome Vu11, were confirmed to be associated with resistance. KASP is a globally recognized technology for SNP genotyping, known for being user-friendly, relatively cheap, and can be automated, ensuring high throughput genotyping (Ongom et al., 2021; Dipta et al., 2024). The five Striga-associated SNP markers designated by Intertek KASP assay IDs; snpVU00075, snpVU00076, snpVU00077, snpVU00078, and snpVU00079 are all positioned proximally to each other within the major QTL region on Vu11. Previously reported marker systems for Striga resistance were mostly AFLPs, SCARs, and SSRs (Ouédraogo et al., 2001, 2002; Boukar et al., 2004; Essem et al., 2019). Although some of these old marker technologies were reportedly effective in marker-assisted selection (MAS) for Striga resistance (Larweh et al., 2017; Omoigui et al., 2017), they have limitations, including labor-intensive gel electrophoresis or capillary electrophoresis, lack of automation and cumbersomeness for high-throughput applications, high cost per data point, and low reproducibility and accuracy (Mut et al., 2008; Dipta et al., 2024). The present study underscores the potential of new KASP-based SNP markers for efficient MAS in cowpea breeding programs aimed at developing Striga-resistant varieties.

5 Conclusion

This study deepens our understanding of the genetic architecture of Striga resistance in cowpea by pinpointing key association signals and candidate genes. A dominant signal on chromosome Vu11 aligns with a known leucine-rich repeat (LRR) resistance gene, reinforcing our confidence in targeting this region for breeding applications. SNPs within these associated intervals were converted to breeder-friendly KASP markers for routine breeding applications. These markers, particularly those on chromosome Vu11, provide valuable tools for breeding programs focused on developing Striga-resistant cowpea varieties, thereby contributing to improved cowpea productivity in Striga-infested regions. Ongoing efforts aim to validate these markers across diverse genetic backgrounds using segregating bi-parental populations to broaden their applicability. Since multiple Striga races affect cowpea across West Africa, pyramiding resistance alleles will be essential for achieving durable, broad-spectrum resistance. Continued research is also needed to identify and validate resistance loci against other predominant Striga races.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Author contributions

PO: Supervision, Formal Analysis, Writing – review & editing, Methodology, Writing – original draft, Data curation, Investigation, Visualization, Conceptualization, Software, Validation. CF: Investigation, Writing – review & editing, Validation, Supervision. OB: Resources, Funding acquisition, Project administration, Supervision, Validation, Writing – review & editing.

Funding

The author (s) declare financial support was received for the research and/or publication of this article. This research was funded in part by Bill & Melinda Gates foundation through the Accelerated Varietal Improvement and Seed Delivery of Legumes and Cereals in Africa (AVISA project, Grant#OPP1198373, and the Foundation for Food and Agriculture Research (FFAR), Grant ID ICRC20-0000000048.

Acknowledgments

We wish to express our sincere gratitude to the International Institute of Tropical Agriculture (IITA) for graciously hosting this study and granting access to their phenotyping facilities. This support was instrumental in enabling the thorough evaluation of cowpea-Striga gesnerioides interactions using the minicore collection, which IITA maintains as part of its global germplasm repository. We also wish to thank the team at the University of California, Riverside (UCR) for genotyping the cowpea minicore population using the high-density iSelect SNP array and making this valuable data publicly available at https://doi.org/10.1002/leg3.95. This open-access resource, covering over 51,000 SNP loci, greatly facilitated our genome-wide association analyses. Finally, we extend our appreciation to all field technicians and support personnel involved in data generation and analysis. Their collective commitment has been central to the discovery and development of KASP markers, which will underpin future marker-assisted breeding efforts for Striga-resistant cowpea varieties.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2025.1661440/full#supplementary-material

References

Abate, T., Alene, A. D., Bergvinson, D., Shiferaw, B., Silim, S., Orr, A., et al. (2012). Tropical grain legumes in Africa and South Asia: Knowledge and opportunities (Nairobi, Kenya: International Crops Research Institute for the Semi-Arid Tropics). Available online at: https://core.ac.uk/download/pdf/211013778.pdf.

Google Scholar

Abdullahi, W. M., Dianda, M., Boukar, O., Dieng, I., Mohammed, G. S., Belko, N., et al. (2022). Integrated management of Striga gesnerioides in cowpea using resistant varieties, improved crop nutrition and rhizobium inoculants. Plant Soil 473, 197–213. doi: 10.1007/S11104-022-05295-7/FIGURES/7

Crossref Full Text | Google Scholar

Abdullah-Zawawi, M. R., Ahmad-Nizammuddin, N. F., Govender, N., Harun, S., Mohd-Assaad, N., and Mohamed-Hussein, Z. A. (2021). Comparative genome-wide analysis of WRKY, MADS-box and MYB transcription factor families in Arabidopsis and rice. Sci. Rep. 11, 1–18. doi: 10.1038/s41598-021-99206-y

PubMed Abstract | Crossref Full Text | Google Scholar

Anand, L. and Rodriguez Lopez, C. M. (2022). ChromoMap: an R package for interactive visualization of multi-omics data and annotation of chromosomes. BMC Bioinf. 23, 1–9. doi: 10.1186/S12859-021-04556-Z/FIGURES/5

PubMed Abstract | Crossref Full Text | Google Scholar

Asare, A. T., Gowda, B. S., Galyuon, I. K. A., Aboagye, L. M., Takrama, J. F., Padi, F. K., et al. (2013). Identification of potential sources of striga resistance in cowpea [Vigna unguiculata (L.) Walp.] accessions from Ghana. J. Microbiol. Biotechnol. Res. Scholars Res. Library J. Microbiol. Biotech. Res. 3, 14–22.

Google Scholar

Atokple, I. D. K., Singh, B. B., and Emechebe, A. M. (1993). Independent inheritance of Striga and Alectra resistance in cowpea genotype B301. Crop Sci. 33, 714–715. doi: 10.2135/CROPSCI1993.0011183X003300040015X

Crossref Full Text | Google Scholar

Atokple, I. D. K., Singh, B. B., and Emechebe, A. M. (1995). Genetics of resistance to Striga and Alectra in Cowpea. J. Heredity. 33 (4), 714–715. doi: 10.1093/oxfordjournals.jhered.a111524

Crossref Full Text | Google Scholar

Blade, S. F., Shetty, S. V. R., Terao, T., and Singh4, B. B. (1997). “Recent developments in cowpea cropping systems research,” in Advances in cowpea research. Eds. Singh, B. B., Raj, M., Dashiell, K. E., and Jackai, L. E. N. (Co-publication of International Institute of Tropical Agriculture (IITA) and Japan International Research Center for Agricultural Sciences (JIRCAS, Ibadan), 114–128. Available online at: https://www.cabi.org/gara/FullTextPDF/2022/20220318355.pdf.

Google Scholar

Botanga, C. J. and Timko, M. P. (2006). Phenetic relationships among different races of Striga gesnerioides (Willd.) Vatke from West Africa. Genome 49, 1351–1365. doi: 10.1139/G06-086

PubMed Abstract | Crossref Full Text | Google Scholar

Boukar, O., Bhattacharjee, R., Fatokun, C., Kumar, P. L., and Gueye, B. (2013). “Cowpea,” in Genetic and Genomic Resources of Grain Legume Improvement. Eds. Singh, M., Bisht, H. D. U., and Singh, I. (Elsevier Inc, London, Uk), 137–156. doi: 10.1016/B978-0-12-397935-3.00006-2

Crossref Full Text | Google Scholar

Boukar, O., Kong, L., Singh, B. B., Murdock, L., and Ohm, H. W. (2004). AFLP and AFLP-derived SCAR markers associated with Striga gesnerioides resistance in Cowpea. Crop Sci. 44, 1259–1264. doi: 10.2135/CROPSCI2004.1259

Crossref Full Text | Google Scholar

Bradbury, P. J., Zhang, Z., Kroon, D. E., Casstevens, T. M., Ramdoss, Y., and Buckler, E. S. (2007). TASSEL: Software for association mapping of complex traits in diverse samples. Bioinformatics 23, 2633–2635. doi: 10.1093/bioinformatics/btm308

PubMed Abstract | Crossref Full Text | Google Scholar

CABI (2020). Striga gesnerioides (cowpea witchweed) (Wallingford, UK: CABI Compendium). doi: 10.1079/cabicompendium.51785

Crossref Full Text | Google Scholar

da Silva, A. C., da C. Santos, D., Junior, D. L. T., da Silva, P. B., dos S. Siviero, R. C., da Silva, A. C., et al. (2018). “Cowpea: A strategic legume species for food security and health,” in Legume Seed Nutraceutical Research (London, United Kingdom: IntechOpen). doi: 10.5772/INTECHOPEN.79006

Crossref Full Text | Google Scholar

de Mendiburu, F. (2020). agricolae: Statistical Procedures for Agricultural Research. R Package. Available online at: https://cran.r-project.org/web/packages/agricolae/index.html (Accessed May 15, 2025).

Google Scholar

Dipta, B., Sood, S., Mangal, V., Bhardwaj, V., Thakur, A. K., Kumar, V., et al. (2024). KASP: a high-throughput genotyping system and its applications in major crop plants for biotic and abiotic stress tolerance. Mol. Biol. Rep. 51, 508. doi: 10.1007/S11033-024-09455-Z

PubMed Abstract | Crossref Full Text | Google Scholar

Essem, F., Ohlson, E. W., Asare, A. T., and Timko, M. P. (2019). Genetic markers linked to Striga gesnerioides resistance for the improvement of Ghanaian cowpea (Vigna unguiculata) cultivars. Plant Breed. 138, 599–604. doi: 10.1111/PBR.12702

Crossref Full Text | Google Scholar

Fatokun, C., Girma, G., Abberton, M., Gedil, M., Unachukwu, N., Oyatomi, O., et al. (2018). Genetic diversity and population structure of a minicore subset from the world cowpea (Vigna unguiculata (L.) Walp.) germplasm collection. Sci. Rep. 8, 16035. doi: 10.1038/s41598-018-34555-9

PubMed Abstract | Crossref Full Text | Google Scholar

Fatokun, C. A., Tarawali, S. A., Singh, B. B., Kormawa, P. M., Tamo, M., and Tamò, M. (2002). “Challenges and opportunities for enhancing sustainable cowpea production,” in Proceedings of the World Cowpea Conference III held at the International Institute of Tropical Agriculture (IITA) (IITA, Ibadan), 7–396. Available online at: https://hdl.handle.net/10568/99857.

Google Scholar

Finkina, E. I., Melnikova, D. N., Bogdanov, I. V., and Ovchinnikova, T. V. (2016). Lipid transfer proteins as components of the plant innate immune system: structure, functions, and applications. Acta Naturae 8, 47. doi: 10.32607/20758251-2016-8-2-47-61

PubMed Abstract | Crossref Full Text | Google Scholar

Glickman, M. E., Rao, S. R., and Schultz, M. R. (2014). False discovery rate control is a recommended alternative to Bonferroni-type adjustments in health studies. J. Clin. Epidemiol. 67, 850–857. doi: 10.1016/j.jclinepi.2014.03.012

PubMed Abstract | Crossref Full Text | Google Scholar

Goff, K. E. and Ramonell, K. M. (2007). The role and regulation of receptor-like kinases in plant defense. Gene Regul. Syst. Bio 1, 167. doi: 10.1177/117762500700100015

PubMed Abstract | Crossref Full Text | Google Scholar

Hao, Y., Liu, R., Mao, Z., Yang, Q., Zheng, S., Lu, X., et al. (2024). Identification and analysis of WRKY transcription factors in response to cowpea fusarium wilt in cowpea. Plants 13, 2273. doi: 10.3390/PLANTS13162273/S1

PubMed Abstract | Crossref Full Text | Google Scholar

Henry, R. (2015). Etymologia: bonferroni correction. Emerg. Infect. Dis. 21, 289. doi: 10.3201/EID2102.ET2102

PubMed Abstract | Crossref Full Text | Google Scholar

Herniter, I. A., Lo, R., Muñoz-Amatriaín, M., Lo, S., Guo, Y.-N., Huynh, B.-L., et al. (2019). Seed coat pattern QTL and development in cowpea (Vigna unguiculata [L.] Walp.). Front. Plant Sci. 0. doi: 10.3389/FPLS.2019.01346

PubMed Abstract | Crossref Full Text | Google Scholar

Herniter, I. A., Muñoz-Amatriaín, M., Lo, S., Guo, Y. N., and Close, T. J. (2018). Identification of candidate genes controlling black seed coat and pod tip color in cowpea (Vigna unguiculata [L.] Walp). G3 Genes|Genomes|Genetics 8, 3347–3355. doi: 10.1534/G3.118.200521

PubMed Abstract | Crossref Full Text | Google Scholar

Heuzé, V., Tran, G., Nozière, P., Bastianelli, D., and Lebas, F. (2015). Cowpea (Vigna unguiculata) forage (Rome, Italy: Feedipedia, a programme by INRAE, CIRAD, AFZ and FAO). Available online at: https://www.feedipedia.org/node/233.

Google Scholar

Huang, Y., Du, L., Wang, M., Ren, M., Yu, S., and Yang, Q. (2022). Multifaceted roles of zinc finger proteins in regulating various agronomic traits in rice. Front. Plant Sci. 13. doi: 10.3389/FPLS.2022.974396/BIBTEX

PubMed Abstract | Crossref Full Text | Google Scholar

Huynh, B.-L., Close, T. J., Roberts, P. A., Hu, Z., Wanamaker, S., Lucas, M. R., et al. (2013). Gene pools and the genetic architecture of domesticated cowpea. Plant Genome 6, 1–8. doi: 10.3835/plantgenome2013.03.0005

Crossref Full Text | Google Scholar

Jayathilake, C., Visvanathan, R., Deen, A., Bangamuwage, R., Jayawardana, B. C., Nammi, S., et al. (2018). Cowpea: an overview on its nutritional facts and health benefits. J. Sci. Food Agric. 98, 4793–4806. doi: 10.1002/jsfa.9074

PubMed Abstract | Crossref Full Text | Google Scholar

Jones, D. A. and Jones, J. D. G. (1997). The role of leucine-rich repeat proteins in plant defences. Adv. Bot. Res. 24, 89–167. doi: 10.1016/S0065-2296(08)60072-5

Crossref Full Text | Google Scholar

Kang, H., Zaitlen, N., Wade, C., Kirby, A., Heckerman, D., Daly, M., et al. (2008). Efficient control of population structure in model organism association mapping. Genetics 178, 1709–1723. doi: 10.1534/genetics.107.080101

PubMed Abstract | Crossref Full Text | Google Scholar

Larweh, V., Akromah, R., Amoah, S., Asibuo, J. Y., Prempeh, R., and Kusi, F. (2017). Marker assisted selection for resistance to Striga gesnerioides in Cowpea (Vigna unguiculata L. Walp). J. Biosci. Biotechnol. Discov. 2, 97–103. doi: 10.31248/JBBD2017.047

Crossref Full Text | Google Scholar

LGC Genomics (2025). SNPviewer | LGC Biosearch Technologies. Available online at: https://www.biosearchtech.com/products/pcr-reagents-kits-and-instruments/pcr-instruments-and-software/genotyping-and-lims-software/snpviewer (Accessed January 24, 2025).

Google Scholar

Li, J., Lis, K. E., and Timko, M. P. (2009). Molecular genetics of race-specific resistance of cowpea to Striga gesnerioides (Willd.). Pest Manag Sci. 65 (5), 520–7. doi: 10.1002/ps.1722

PubMed Abstract | Crossref Full Text | Google Scholar

Li, J. and Timko, M. P. (2009). Gene-for-gene resistance in Striga-cowpea associations. Science 325, 1094. doi: 10.1126/SCIENCE.1174754

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, X., Huang, M., Fan, B., Buckler, E. S., and Zhang, Z. (2016). Iterative Usage of fixed and random effect models for powerful and efficient genome-wide association studies. PloS Genet. 12, e1005767. doi: 10.1371/JOURNAL.PGEN.1005767

PubMed Abstract | Crossref Full Text | Google Scholar

Lo, S., Muñoz-Amatriaín, M., Boukar, O., Herniter, I., Cisse, N., Guo, Y. N., et al. (2018). Identification of QTL controlling domestication-related traits in cowpea (Vigna unguiculata L. Walp). Sci. Rep. 8, 1–9. doi: 10.1038/s41598-018-24349-4

PubMed Abstract | Crossref Full Text | Google Scholar

Lo, S., Muñoz-Amatriaín, M., Hokin, S. A., Cisse, N., Roberts, P. A., Farmer, A. D., et al. (2019). A genome-wide association and meta-analysis reveal regions associated with seed size in cowpea [Vigna unguiculata (L.) Walp. Theor. Appl. Genet. 132, 3079–3087. doi: 10.1007/s00122-019-03407-z

PubMed Abstract | Crossref Full Text | Google Scholar

Lonardi, S., Muñoz-Amatriaín, M., Liang, Q., Shu, S., Wanamaker, S. I., Lo, S., et al. (2019). The genome of cowpea (Vigna unguiculata [L.] Walp.). Plant J. 98, 767–782. doi: 10.1111/tpj.14349

PubMed Abstract | Crossref Full Text | Google Scholar

Lucas, M. R., Diop, N.-N., Wanamaker, S., Ehlers, J. D., Roberts, P. A., and Close, T. J. (2011). Cowpea–soybean synteny clarified through an improved genetic map. Plant Genome 4, 218–225. doi: 10.3835/PLANTGENOME2011.06.0019

Crossref Full Text | Google Scholar

McHale, L., Tan, X., Koehl, P., and Michelmore, R. W. (2006). Plant NBS-LRR proteins: Adaptable guards. Genome Biol. 7, 1–11. doi: 10.1186/GB-2006-7-4-212/FIGURES/5

PubMed Abstract | Crossref Full Text | Google Scholar

Meng, L., Du, M., Zhu, T., Li, G., Ding, Y., and Zhang, Q. (2024). PPR proteins in plants: roles, mechanisms, and prospects for rice research. Front. Plant Sci. 15. doi: 10.3389/FPLS.2024.1416742/BIBTEX

PubMed Abstract | Crossref Full Text | Google Scholar

Miesho, B., Hailay, M., Msiska, U., Bruno, A., Malinga, G. M., Obia Ongom, P., et al. (2019). Identification of candidate genes associated with resistance to bruchid (Callosobruchus maculatus) in cowpea. Plant Breed. 138, 605–613. doi: 10.1111/pbr.12705

Crossref Full Text | Google Scholar

Mohamed, K. I., Musselman, J., and Riches, C. R. (2001). The genus Striga (Scrophulariaceae) in Africa. Ann. Missouri Botanical Garden 88, 60–103. doi: 10.2307/2666132

Crossref Full Text | Google Scholar

Mortimore, M., Singh, B., Harris, F., and Blade, S. (1997). “Cowpea in traditional cropping systems,” in Advances in cowpea research. Eds. Singh, B., Mohan, R. D., Dashiel, K., and Jackai, L. (Co-publication of International Institute of Tropical Agriculture (IITA) and Japan International Research Center for Agricultural Sciences (JIRCAS, Ibadan), 99–113.

Google Scholar

Muchero, W., Diop, N. N., Bhat, P. R., Fenton, R. D., Wanamaker, S., Pottorff, M., et al. (2009). A consensus genetic map of cowpea [Vigna unguiculata (L) Walp.] and synteny based on EST-derived SNPs. Proc. Natl. Acad. Sci. U.S.A. 106, 18159–18164. doi: 10.1073/pnas.0905886106

PubMed Abstract | Crossref Full Text | Google Scholar

Muchero, W., Roberts, P. A., Diop, N. N., Drabo, I., Cisse, N., Close, T. J., et al. (2013). Genetic architecture of delayed senescence, biomass, and grain yield under drought stress in cowpea. PloS One 8, e70041. doi: 10.1371/JOURNAL.PONE.0070041

PubMed Abstract | Crossref Full Text | Google Scholar

Müller, S., Hauck, C., and Schildknecht, H. (1992). Germination stimulants produced by Vigna unguiculata Walp cv Saunders Upright. J. Plant Growth Regul. 11, 77–84. doi: 10.1007/BF00198018/METRICS

Crossref Full Text | Google Scholar

Müller, M. and Munné-Bosch, S. (2015). Ethylene response factors: A key regulatory hub in hormone and stress signaling. Plant Physiol. 169, 32. doi: 10.1104/PP.15.00677

PubMed Abstract | Crossref Full Text | Google Scholar

Muñoz-Amatriaín, M., Lo, S., Herniter, I. A., Boukar, O., Fatokun, C., Carvalho, M., et al. (2021). The UCR Minicore: a resource for cowpea research and breeding. Legume Sci. 3, e95. doi: 10.1002/LEG3.95

Crossref Full Text | Google Scholar

Muñoz-Amatriaín, M., Mirebrahim, H., Xu, P., Wanamaker, S. I., Luo, M. C., Alhakami, H., et al. (2017). Genome resources for climate-resilient cowpea, an essential crop for food security. Plant J. 89, 1042–1054. doi: 10.1111/tpj.13404

PubMed Abstract | Crossref Full Text | Google Scholar

Mut, E., Tiedemann, A., Koopmann, B., Muiru, W., and Kimenju, J. (2008). A review of the Amplified Fragment Length Polymorphism (AFLP) technique in genotyping and DNA fingerprinting studies. J. Appl. Biosci. 9, 389–395.

Google Scholar

Nkomo, G. V., Sedibe, M. M., and Mofokeng, M. A. (2021). Production constraints and improvement strategies of Cowpea (Vigna unguiculata L. Walp.) genotypes for drought tolerance. Int. J. Agron. 2021, 1–9. doi: 10.1155/2021/5536417

Crossref Full Text | Google Scholar

Ohlson, E. W. and Timko, M. P. (2020). Race structure of cowpea witchweed (Striga gesnerioides) in West Africa and its implications for Striga resistance breeding of cowpea. Weed Sci. 68, 125–133. doi: 10.1017/WSC.2020.3

Crossref Full Text | Google Scholar

Omoigui, L. O., Kamara, A. Y., Ishiyaku, M. F., and Boukar, O. (2010). Comparative responses of cowpea breeding lines to Striga and Alectra in the dry savanna of northeast Nigeria. Afr J. Agric. Res. 7, 747–754. doi: 10.5897/AJAR11.1341

Crossref Full Text | Google Scholar

Omoigui, L. O., Kamara, A. Y., Moukoumbi, Y. D., Ogunkanmi, L. A., and Timko, M. P. (2017). Breeding cowpea for resistance to Striga gesnerioides in the Nigerian dry savannas using marker-assisted selection. Plant Breed. 136, 393–399. doi: 10.1111/PBR.12475

Crossref Full Text | Google Scholar

Ongom, P., Fatokun, C., Togola, A., Garcia-Oliveira, A., Kilian, A., Lonardi, S., et al. (2024). A Mid-density single-nucleotide polymorphism panel for molecular applications in Cowpea (Vigna unguiculata (L.) Walp). Int. J. Genomics 2024, 9912987. doi: 10.1155/2024/9912987

PubMed Abstract | Crossref Full Text | Google Scholar

Ongom, P. O., Fatokun, C., Togola, A., Salvo, S., Oyebode, O. G., Ahmad, M. S., et al. (2021). Molecular fingerprinting and hybridity authentication in cowpea using single nucleotide polymorphism based Kompetitive Allele-Specific PCR Assay. Front. Plant Sci. 12. doi: 10.3389/FPLS.2021.734117/BIBTEX

PubMed Abstract | Crossref Full Text | Google Scholar

Ongom, P. O., Togola, A., Fatokun, C., and Boukar, O. (2022). A genome-wide scan divulges key loci involved in resistance to aphids (Aphis craccivora) in cowpea (Vigna unguiculata). Genes (Basel) 13, 2002. doi: 10.3390/GENES13112002/S1

PubMed Abstract | Crossref Full Text | Google Scholar

Ouédraogo, J. T., Maheshwari, V., Berner, D. K., St-Pierre, C. A., Belzile, F., and Timko, M. P. (2001). Identification of AFLP markers linked to resistance of cowpea (Vigna unguiculata L.) to parasitism by Striga gesnerioides. Theor. Appl. Genet. 102, 1029–1036. doi: 10.1007/S001220000499/METRICS

Crossref Full Text | Google Scholar

Ouédraogo, J. T., Ouedraogo, M., Gowda, B. S., and Timko, M. P. (2012). Development of sequence characterized amplified region (SCAR) markers linked to race-specific resistance to Striga gesnerioides in cowpea (Vigna unguiculata L.). Afr J. Biotechnol. 11, 12555–12562.

Google Scholar

Ouédraogo, J. T., Tignegre, J. B., Timko, M. P., and Belzile, F. J. (2002). AFLP markers linked to resistance against Striga gesnerioides race 1 in cowpea (Vigna unguiculata). Genome 45, 787–793. doi: 10.1139/G02-043

PubMed Abstract | Crossref Full Text | Google Scholar

Panzeri, D., Nissim, W. G., Labra, M., and Grassi, F. (2022). Revisiting the domestication process of African vigna species (Fabaceae): background, perspectives and challenges. Plants 11, 532. doi: 10.3390/PLANTS11040532

PubMed Abstract | Crossref Full Text | Google Scholar

Parker, C. and Riches, C. R. (1993). Parasitic Weeds of the World: Biology and Control (Wallingford, UK: CAB International). Available online at: https://books.google.com/books/about/Parasitic_Weeds_of_the_World.html?id=KIJFAQAAIAAJ.

Google Scholar

Paudel, D., Dareus, R., Rosenwald, J., Muñoz-Amatriaín, M., and Rios, E. F. (2021). Genome-wide association study reveals candidate genes for flowering time in cowpea (Vigna unguiculata [L.] Walp.). Front. Genet. 12. doi: 10.3389/FGENE.2021.667038/FULL

PubMed Abstract | Crossref Full Text | Google Scholar

Pottorff, M. O., Li, G., Ehlers, J. D., Close, T. J., and Roberts, P. A. (2014). Genetic mapping, synteny, and physical location of two loci for Fusarium oxysporum f. sp. tracheiphilum race 4 resistance in cowpea [Vignaunguiculata (L.) Walp. Mol. Breed 33, 779–791. doi: 10.1007/S11032-013-9991-0

PubMed Abstract | Crossref Full Text | Google Scholar

Sadda, A. S., d’Eeckenbrugge, G. C., Saidou, A. A., Diouf, A., Jangorzo, N. S., Issoufou, H. B. A., et al. (2021). The witchweed Striga gesnerioides and the cultivated cowpea: A geographical and historical analysis of their West African distribution points to the prevalence of agro-ecological factors and the parasite’s multilocal evolution potential. PLoS One 16, e0254803. doi: 10.1371/JOURNAL.PONE.0254803

PubMed Abstract | Crossref Full Text | Google Scholar

Singh, B. B. (2014). Cowpea: The Food Legume of the 21st Century (Madison, WI, USA: Crop Science Society of America). doi: 10.2135/2014.COWPEA

Crossref Full Text | Google Scholar

Singh, B. B. and Emechebe, A. M. (1990). Inheritance of striga resistance in cowpea genotype B301. Crop Sci. 30, 879–881. doi: 10.2135/CROPSCI1990.0011183X003000040023X

Crossref Full Text | Google Scholar

Singh, B. B., Olufajo, O. O., Ishiyaku, M. F., Adeleke, R. A., Ajeigbe, H. A., and Mohammed, S. G. (2006). Registration of Six Improved Germplasm Lines of Cowpea with Combined Resistance to Striga gesnerioides and Alectra vogelii. Crop Sci. 46, 2332–2333. doi: 10.2135/CROPSCI2006.03.0148

Crossref Full Text | Google Scholar

Singh, B. B., Timko, M. P., and Aragao, F. J. L. (2014). “Advances in cowpea improvement and genomics,” in Legumes in the Omic Era. Eds. Gupta, S., Nadarajan, N., and Gupta, D. S. (Springer Science, New York), 131–153. doi: 10.1007/978-1-4614-8370-0_7/COVER

Crossref Full Text | Google Scholar

Timko, M. P., Gowda, B. S., Ouedraogo, J., and Ousmane, B. (2007). “Molecular markers for analysis of resistance to Striga Gesnerioides in cowpea,” in Integrating New Technologies for Striga Control: Towards Ending the Witch-Hunt, (London, UK: World Scientific Publishing Co. Pte. Ltd) 115–128. doi: 10.1142/9789812771506_0009

Crossref Full Text | Google Scholar

Touré, M., Olivier, A., Ntare, B. R., Lane, J. A., and St-Pierre, C. A. (1997). Inheritance of resistance to Striga gesnerioides biotypes from Mali and Niger in cowpea (Vigna unguiculata (L.) Walp.). Euphytica 94, 273–278. doi: 10.1023/A:1002927705470/METRICS

Crossref Full Text | Google Scholar

Tsuchihashi, Z. and Dracopoli, N. C. (2002). Progress in high throughput SNP genotyping methods. Pharmacogenomics J. 2, 103–110. doi: 10.1038/sj.tpj.6500094

PubMed Abstract | Crossref Full Text | Google Scholar

Vaillancourt, R. E. and Weeden, N. F. (1992). Chloroplast DNA polymorphism suggests Nigerian center of domestication for the Cowpea, Vigna unguiculata (Leguminosae). Am. J. Bot. 79, 1194–1199. doi: 10.1002/J.1537-2197.1992.TB13716.X

PubMed Abstract | Crossref Full Text | Google Scholar

Wu, X., Wang, B., Wu, S., Li, S., Zhang, Y., Wang, Y., et al. (2021). Development of a core set of single nucleotide polymorphism markers for genetic diversity analysis and cultivar fingerprinting in cowpea. Legume Sci. 3, e93. doi: 10.1002/LEG3.93

Crossref Full Text | Google Scholar

Yin, L. (2016). CMplot: Circular and Rectangular Manhattan Plot (Wuhan, Hubei, China: Github). Available online at: https://github.com/YinLiLin/CMplot?tab=readme-ov-file.

Google Scholar

Yin, L., Zhang, H., Tang, Z., Xu, J., Yin, D., Zhang, Z., et al. (2021). rMVP: A memory-efficient, visualization-enhanced, and parallel-accelerated tool for genome-wide association study. Genomics Proteomics Bioinf. 19, 619–628. doi: 10.1016/J.GPB.2020.10.007

PubMed Abstract | Crossref Full Text | Google Scholar

Yu, Y. (2020). Remorins: essential regulators in plant-microbe interaction and cell death induction. Plant Physiol. 183, 435. doi: 10.1104/PP.20.00490

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, Z., Ersoz, E., Lai, C.-Q., Todhunter, R. J., Tiwari, H. K., Gore, M. A., et al. (2010). Mixed linear model approach adapted for genome-wide association studies. Nat. Genet. 42, 355–360. doi: 10.1038/ng.546

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, Y., Held, M. A., and Showalter, A. M. (2020). Elucidating the roles of three β-glucuronosyltransferases (GLCATs) acting on arabinogalactan-proteins using a CRISPR-Cas9 multiplexing approach in Arabidopsis. BMC Plant Biol. 20, 1–20. doi: 10.1186/S12870-020-02420-5/FIGURES/15

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, X., Ma, J., Yang, S., Yao, W., Zhang, N., Hao, X., et al. (2023a). Analysis of GATA transcription factors and their expression patterns under abiotic stress in grapevine (Vitis vinifera L.). BMC Plant Biol. 23, 1–17. doi: 10.1186/S12870-023-04604-1/FIGURES/8

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, Y., Tessaro, M. J., Lassner, M., and Li, X. (2003). Knockout analysis of arabidopsis transcription factors TGA2, TGA5, and TGA6 reveals their redundant and essential roles in systemic acquired resistance. Plant Cell 15, 2647. doi: 10.1105/TPC.014894

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, Y., Tian, H., Chen, D., Zhang, H., Sun, M., Chen, S., et al. (2023b). Cysteine-rich receptor-like protein kinases: emerging regulators of plant stress responses. Trends Plant Sci. 28, 776–794. doi: 10.1016/J.TPLANTS.2023.03.028

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: cowpea (Vigna unguiculata), Striga gesnerioides, genome-wide association study, candidate genes, SNP-based KASP markers, marker validation

Citation: Ongom PO, Fatokun CA and Boukar O (2025) Discovery and development of single-nucleotide polymorphism markers for resistance to Striga gesnerioides in cowpea (Vigna unguiculata). Front. Plant Sci. 16:1661440. doi: 10.3389/fpls.2025.1661440

Received: 21 July 2025; Accepted: 15 September 2025;
Published: 25 September 2025.

Edited by:

Anna Maria Mastrangelo, Council for Agricultural and Economics Research (CREA), Italy

Reviewed by:

Ofere Francis Emeriewen, Julius Kühn Institute (JKI), Germany
Goitsemang Dikane, Central University of Technology, South Africa

Copyright © 2025 Ongom, Fatokun and Boukar. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Patrick Obia Ongom, UC5PbmdvbUBjZ2lhci5vcmc=

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.