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.

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.

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 Diversity and Population Structure of Cowpea [Vigna unguiculata (L.) Walp.] Germplasm Collected from Togo Based on DArT Markers

Crop genetic diversity is a sine qua non for continuous progress in the development of improved varieties, hence the need for germplasm collection, conservation and characterization. Over the years, cowpea has contributed immensely to the nutrition and economic life of the people in Togo. However, the bulk of varieties grown by farmers are landraces due to the absence of any serious genetic improvement activity on cowpea in the country. In this study, the genetic diversity and population structure of 255 cowpea accessions collected from five administrative regions and the agricultural research institute of Togo were assessed using 4600 informative diversity array technology (DArT) markers. Among the regions, the polymorphic information content (PIC) ranged from 0.19 to 0.27 with a mean value of 0.25. The expected heterozygosity (He) varied from 0.22 to 0.34 with a mean value of 0.31, while the observed heterozygosity (Ho) varied from 0.03 to 0.07 with an average of 0.05. The average inbreeding coefficient (FIS) varied from 0.78 to 0.89 with a mean value of 0.83, suggesting that most of the accessions are inbred. Cluster analysis and population structure identified four groups with each comprising accessions from the six different sources. Weak to moderate differentiation was observed among the populations with a genetic differentiation index varying from 0.014 to 0.117. Variation was highest (78%) among accessions within populations and lowest between populations (7%). These results revealed a moderate level of diversity among the Togo cowpea germplasm. The findings of this study constitute a foundation for genetic improvement of cowpea in Togo.

Genetic worth of multiple sets of cowpea breeding lines destined for advanced yield testing

The objective of this study was to determine genetic potentials in eight sets of cowpea lines for grain yield (GY), hundred seed weight (HSDWT) and days to 50% flowering (DT50FL). A total of 614 F₆ genotypes constituting the sets, grouped by maturity, were evaluated across two locations in Northern Nigeria, in an alpha lattice design, two replications each. Data were recorded on GY, HSDWT and DT50FL.Variance components, genotypic coefficient of variation (GCV), and genetic advance (G_A) were used to decode the magnitude of genetic variance within and among sets. Genetic usefulness (Up) which depends on mean and variance to score the genetic merits in historically bi-parental populations was applied to groups of breeding lines with mixed parentage. Principal component analysis (PCA) was used to depict contribution of traits to observed variations. GY and DT50FL explained the variance within and between sets respectively. Genotypes were significantly different, although genotype-by-location and set-by-location interaction effects were also prominent. Genetic variance (δ² _G) and GCV were high for GY in Prelim2 (δ² _G = 45,897; GCV = 19.58%), HSDWT in Prelim11 (δ² _G = 7.137; GCV = 17.07%) and DT50F in Prelim5 (δ² _G = 4.54; GCV = 4.4%). Heritability varied among sets for GY (H = 0.21 to 0.57), HSDWT (H = 0.76 to 0.93) and DT50FL (H = 0.20 to 0.81). G_A and percentage G_A (G_APM) were high for GY in Prelim2 (G_APM = 24.59%; G_A = 269.05Kg/ha), HSDWT in Prelim11 (G_APM = 28.54%; G_A = 4.47 g), and DT50F in Prelim10 (G_APM = 6.49%; G_A = 3.01 days). These sets also registered high values of genetic usefulness, suggesting potential application in non-full sib populations. These approaches can be used during preliminary performance tests to reinforce decisions in extracting promising lines and choose among defined groups of lines.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at (10.1007/s10681-020-02763-y).

Identification of QTL for perenniality and floral scent in cowpea (Vigna unguiculata [L.] Walp.)

Perennial habit and floral scent are major traits that distinguish domesticated cowpeas from their wild relatives. However, the genetic basis of these two important traits remains largely unknown in cowpea. Plant longevity, a perenniality-related trait, and floral scent, an outcrossing trait, were investigated using a RIL population derived from a cross between a domesticated and a wild cowpea. QTL analysis revealed three significant loci, one on chromosome 8 associated with plant longevity and two, on chromosomes 1 and 11, for floral scent. Genes within the QTL regions were identified. Genes encoding an F-box protein (Vigun08g215300) and two kinases (Vigun08g217000, Vigun08g217800), and involved in physiological processes including regulation of flowering time and plant longevity, were identified within the perenniality QTL region. A cluster of O-methyltransferase genes (Vigun11g096800, Vigun11g096900, Vigun11g097000, Vigun11g097600, and Vigun11g097800) was identified within the floral scent QTL region. These O-methyltransferase cowpea genes are orthologs of the Arabidopsis N-acetylserotonin O-methyltransferase (ASMT) gene, which is involved in the biosynthesis of melatonin. Melatonin is an indole derivative, which is an essential molecule for plant interactions with pollinators. These findings lay the foundation for further exploration of the genetic mechanisms of perenniality and floral scent in cowpea. Knowledge from this study can help in the development of new extended-growth cycle lines with increased yield or lines with increased outcrossing for population breeding.

Identification of sources of resistance in cowpea mini core accessions to Aphis craccivora Koch (Homoptera: Aphididae) and their biochemical characterization

Cowpea (Vigna unguiculata (L. Walp) is an important grain legume for human and livestock nutrition, especially in sub-Saharan Africa. Aphid, Aphis craccivora Koch (Homoptera: Aphididae), is one of the most widespread and destructive insect pests of cowpea and host-plant resistance is an effective approach to minimize the pest damage at seedling stage. This study was aimed at identifying resistant sources to A. craccivora within the cowpea mini core collection, a set of accessions from the largest world cowpea germplasm collection maintained at the International Institute of Tropical Agriculture (IITA). A total of 375 lines including 373 from IITA mini core collection, one resistant (TVu-801) and one susceptible (TVx-3236) checks were evaluated through artificial infestation in screening cages during the seedling stage. In cages, genotypes were planted in single rows containing four plants. They were arranged in an augmented design in which the two checks were sown in individual cages. Scoring for aphid population and damage levels were carried out on individual plants at 7, 14, and 21 days after planting. Advanced bioassays and biochemical analyses were conducted to investigate the mechanism of resistance to A. craccivora. Overall, three genotypes TVu-6464, TVu-1583, and TVu-15445 showed good levels of resistance comparable to the resistant check TVu-801. The HPLC analyses proved that both low sucrose levels in the plant, as well as high levels of kaempferol and quercetin, aglycones of phenolic compounds, were related with high resistance to aphids. The above genotypes with promising levels of resistance to A. craccivora will be used in cowpea breeding programs to develop improved resistant lines against this pest.

Genetic dissection of yield associated traits in a cross between cowpea and yard-long bean (Vigna unguiculata (L.) Walp.) based on DArT markers

Both cowpea and yard-long bean belong to Vigna unguiculata ssp. unguiculata but have diverged through human induced evolution in sub-Saharan Africa and Asia, respectively. To map the quantitative trait loci (QTLs) for yield associated traits and derive new lines that may combine the attributes of both types, we developed a F_2:3 mapping population derived from a cross between cowpea line TVu2185 and yard-long bean line TVu6642. Using DArT markers, a total of 30 QTLs accounting for 1.8-13.0% phenotypic variation was detected for pod and seed traits. Some novel major QTLs for peduncle number per plant (qPeN2.2), pod length (qPoL3), seed breadth (qSB4), length (qSL7.2) and thickness (qST9) identified on chromosomes 2, 3, 4, 7 and 9, respectively, are particularly interesting and need to be validated. Moreover, we confirmed previously reported QTLs for pod length (qPoL8) and 100-seed weight (qSW8) on chromosome 8 and for seed number per pod (qSN9.2) on chromosome 9 suggesting usefulness for marker-assisted-selection purpose. Notably, some QTLs for these traits were clustered especially on chromosomes 5, 7, 8, 9 and 10 indicating the presence of the same QTL or linked loci in these regions. Moreover, the involvement of epistasis was observed for trait expressions, but compared with the main effect QTLs, the phenotypic effects of epistatic-QTLs detected were much less. The present QTL analysis may provide a useful tool for breeders to formulate efficientbreeding strategy for introgression of the desirable alleles for yield related traits in cowpea using molecular markers.

Introgression Breeding in Cowpea [Vigna unguiculata (L.) Walp.]

The narrow base of genetic diversity characteristic of cowpea can be attributed to it being self-pollinating, evolving from narrow wild germplasm and exhibiting very limited gene flow between wild and cultivated types. Backcrossing to introduce simply inherited desirable traits and utilization of improved breeding lines and varieties as parents in crossing programs further narrowed the genetic base of cowpea varieties. In most cowpea breeding programs, genes for resistance and market traits were pyramided into lines characterized by high levels of acceptance to farmers and consumers. Besides predisposing widely distributed improved varieties to genetic vulnerability, a narrow base of genetic variation may be contributing to the plateauing in cowpea grain yield, which compromises genetic gains. Cross compatible wild relatives have not been used in variety development because breeders shy away from them due to their tiny seed size, unattractive seed coat color and texture, pod shattering, and susceptibility to viruses. A number of wild cowpea relatives, both within and outside section Catiang of Vigna species, have been evaluated for their reaction to cowpea insect pests and diseases. Vigna vexillata lines were resistant to the legume pod borer (Maruca vitrata), the cowpea weevil (Callosobruchus maculatus), and Striga gesnerioides but are cross incompatible with cultivated cowpea. Some lines among the cross compatible wild relative V. unguiculata ssp. dekindtiana were found to be resistant to aphid in the seedling stage, while others showed good levels of drought and heat tolerance. Molecular markers are being generated to identify quantitative trait loci (QTL) with effects on some desirable attributes in cowpea. Modern breeding tools, including transgenics, can be applied for the improvement of cowpea, bypassing the natural barriers of traditional breeding. Transgenic cowpea with Bt gene cry1Ab showing resistance to M. vitrata has been released in Nigeria. Genome editing, a powerful emerging tool, can also be used for developing improved cowpea varieties with durable resistance to pests and diseases.

Evaluation of cowpea mini core accessions for resistance to flower bud thrips Megalurothrips sjostedti Trybom (Thysanoptera: Thripidae)

The flower bud thrips, Megalurothrips sjostedti Trybom (Thysanoptera: Thripidae), is an economically important pest of cowpea in sub-Saharan Africa. Varietal resistance is the most preferred, environmentally friendly, cost-effective and sustainable option for controlling this pest. The objective of this study was to identify sources of resistance to M. sjostedti among mini core accessions from the largest world cowpea germplasm collection maintained at the International Institute of Tropical Agriculture (IITA). The study was conducted during the 2015 and 2016 cropping seasons where 365 accessions were screened under field conditions. Each accession was rated visually for thrips damage score, flower abortion rate, number of pods per plant and number of thrips per flower. The resistance levels observed in genotypes TVu8631, TVu16368, TVu8671 and TVu7325 were similar to that of the resistant check "Sanzisabinli" (called Sanzi) during both seasons. In addition, 56 mini core genotypes showed moderate resistance to thrips damage. High heritability values were associated with thrips damage scores at 65 days after planting (0.60), percentage of effective peduncles (0.59), flower bud abortion rate (0.59), number of pods per plant (0.51) and number of peduncles with pods (0.5). The accessions identified with good levels of resistance to flower bud thrips will be used in cowpea breeding programs to develop improved resistant varieties.

Genetic diversity and population structure of a mini-core subset from the world cowpea (Vigna unguiculata (L.) Walp.) germplasm collection

The International Institute of Tropical Agriculture maintains the world's largest collection of cowpea germplasm of over 15,000 accessions. A sub-set of 298 lines from the loosely composed mini core collection of 370 landraces were genotyped based on genotyping by sequencing (GBS). Ward's minimum variance hierarchical cluster analysis, model-based ancestry analysis and discriminant analysis of principal component (DAPC) were carried out on this sub-set. Three clusters were identified by the different clustering methods. Principal component analysis further supported the three clusters especially when accessions are scattered along the axes of the first two principal components. The first two principal components explained a total of 22.30% of the variation. Cluster one comprises 115 accessions from the largest number of countries and has the highest gene diversity, heterozygosity and polymorphic information content (PIC) values. Cluster two is made up of 102 accessions, 90 percent of which are from West and Central Africa. Analysis of molecular variance shows that the most variation is among accessions and lowest among clusters. No cluster is made exclusively of accessions from a single country. Based on SNP markers, the sub set of cowpea mini core germplasm collection used in this study encompasses the diversity in the crop.

Identification of QTL controlling domestication-related traits in cowpea (Vigna unguiculata L. Walp)

Cowpea (Vigna unguiculata L. Walp) is a warm-season legume with a genetically diverse gene-pool composed of wild and cultivated forms. Cowpea domestication involved considerable phenotypic changes from the wild progenitor, including reduction of pod shattering, increased organ size, and changes in flowering time. Little is known about the genetic basis underlying these changes. In this study, 215 recombinant inbred lines derived from a cross between a cultivated and a wild cowpea accession were used to evaluate nine domestication-related traits (pod shattering, peduncle length, flower color, days to flowering, 100-seed weight, pod length, leaf length, leaf width and seed number per pod). A high-density genetic map containing 17,739 single nucleotide polymorphisms was constructed and used to identify 16 quantitative trait loci (QTL) for these nine traits. Based on annotations of the cowpea reference genome, genes within these regions are reported. Four regions with clusters of QTL were identified, including one on chromosome 8 related to increased organ size. This study provides new knowledge of the genomic regions controlling domestication-related traits in cowpea as well as candidate genes underlying those QTL. This information can help to exploit wild relatives in cowpea breeding programs.

A multi-parent advanced generation inter-cross (MAGIC) population for genetic analysis and improvement of cowpea (Vigna unguiculata L. Walp.)

Multi-parent advanced generation inter-cross (MAGIC) populations are an emerging type of resource for dissecting the genetic structure of traits and improving breeding populations. We developed a MAGIC population for cowpea (Vigna unguiculata L. Walp.) from eight founder parents. These founders were genetically diverse and carried many abiotic and biotic stress resistance, seed quality and agronomic traits relevant to cowpea improvement in the United States and sub-Saharan Africa, where cowpea is vitally important in the human diet and local economies. The eight parents were inter-crossed using structured matings to ensure that the population would have balanced representation from each parent, followed by single-seed descent, resulting in 305 F₈ recombinant inbred lines each carrying a mosaic of genome blocks contributed by all founders. This was confirmed by single nucleotide polymorphism genotyping with the Illumina Cowpea Consortium Array. These lines were on average 99.74% homozygous but also diverse in agronomic traits across environments. Quantitative trait loci (QTLs) were identified for several parental traits. Loci with major effects on photoperiod sensitivity and seed size were also verified by biparental genetic mapping. The recombination events were concentrated in telomeric regions. Due to its broad genetic base, this cowpea MAGIC population promises breakthroughs in genetic gain, QTL and gene discovery, enhancement of breeding populations and, for some lines, direct releases as new varieties.

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

Cowpea (Vigna unguiculata L. Walp.) is a legume crop that is resilient to hot and drought-prone climates, and a primary source of protein in sub-Saharan Africa and other parts of the developing world. However, genome resources for cowpea have lagged behind most other major crops. Here we describe foundational genome resources and their application to the analysis of germplasm currently in use in West African breeding programs. Resources developed from the African cultivar IT97K-499-35 include a whole-genome shotgun (WGS) assembly, a bacterial artificial chromosome (BAC) physical map, and assembled sequences from 4355 BACs. These resources and WGS sequences of an additional 36 diverse cowpea accessions supported the development of a genotyping assay for 51 128 SNPs, which was then applied to five bi-parental RIL populations to produce a consensus genetic map containing 37 372 SNPs. This genetic map enabled the anchoring of 100 Mb of WGS and 420 Mb of BAC sequences, an exploration of genetic diversity along each linkage group, and clarification of macrosynteny between cowpea and common bean. The SNP assay enabled a diversity analysis of materials from West African breeding programs. Two major subpopulations exist within those materials, one of which has significant parentage from South and East Africa and more diversity. There are genomic regions of high differentiation between subpopulations, one of which coincides with a cluster of nodulin genes. The new resources and knowledge help to define goals and accelerate the breeding of improved varieties to address food security issues related to limited-input small-holder farming and climate stress.

Leaf morphology in Cowpea [Vigna unguiculata (L.) Walp]: QTL analysis, physical mapping and identifying a candidate gene using synteny with model legume species

BACKGROUND

Cowpea [Vigna unguiculata (L.) Walp] exhibits a considerable variation in leaf shape. Although cowpea is mostly utilized as a dry grain and animal fodder crop, cowpea leaves are also used as a high-protein pot herb in many countries of Africa.

RESULTS

Leaf morphology was studied in the cowpea RIL population, Sanzi (sub-globose leaf shape) x Vita 7 (hastate leaf shape). A QTL for leaf shape, Hls (hastate leaf shape), was identified on the Sanzi x Vita 7 genetic map spanning from 56.54 cM to 67.54 cM distance on linkage group 15. SNP marker 1_0910 was the most significant over the two experiments, accounting for 74.7% phenotypic variance (LOD 33.82) in a greenhouse experiment and 71.5% phenotypic variance (LOD 30.89) in a field experiment. The corresponding Hls locus was positioned on the cowpea consensus genetic map on linkage group 4, spanning from 25.57 to 35.96 cM. A marker-trait association of the Hls region identified SNP marker 1_0349 alleles co-segregating with either the hastate or sub-globose leaf phenotype. High co-linearity was observed for the syntenic Hls region in Medicago truncatula and Glycine max. One syntenic locus for Hls was identified on Medicago chromosome 7 while syntenic regions for Hls were identified on two soybean chromosomes, 3 and 19. In all three syntenic loci, an ortholog for the EZA1/SWINGER (AT4G02020.1) gene was observed and is the candidate gene for the Hls locus. The Hls locus was identified on the cowpea physical map via SNP markers 1_0910, 1_1013 and 1_0992 which were identified in three BAC contigs; contig926, contig821 and contig25.

CONCLUSIONS

This study has demonstrated how integrated genomic resources can be utilized for a candidate gene approach. Identification of genes which control leaf morphology may be utilized to improve the quality of cowpea leaves for vegetable and or forage markets as well as contribute to more fundamental research understanding the control of leaf shape in legumes.

Regulatory considerations surrounding the deployment of Bt-expressing cowpea in Africa: report of the deliberations of an expert panel

Cowpea (Vigna unguiculata spp unguiculata) is adapted to the drier agro-ecological zones of West Africa where it is a major source of dietary protein and widely used as a fodder crop. Improving the productivity of cowpea can enhance food availability and security in West Africa. Insect predation--predominately from the legume pod borer (Maruca vitrata), flower thrips (Megalurothrips sjostedti) and a complex of pod-sucking bugs (e.g., Clavigralla spp)--is a major yield-limiting factor in West African cowpea production. Dramatic increases in yield are shown when M. vitrata is controlled with insecticides. However, availability, costs, and safety considerations limit pesticides as a viable option for boosting cowpea production. Development of Bt-cowpea through genetic modification (GM) to control the legume pod borer is a promising approach to cowpea improvement. Cowpea expressing the lepidopteran-active Cry1Ab protein from Bacillus thuringiensis is being developed as a first generation Bt-cowpea crop for West Africa. Appropriate stewardship of Bt-cowpea to assure its sustainability under West African conditions is critical to its successful development. A first step in this process is an environmental risk assessment to determine the likelihood and magnitude of adverse effects of the Cry1Ab protein on key environmental protection goals in West Africa. Here we describe the results of an expert panel convened in 2009 to develop the problem formulation phase for Bt-cowpea and to address specific issues around gene flow, non-target arthropods, and insect resistance management.

A consensus genetic map of cowpea [Vigna unguiculata (L) Walp.] and synteny based on EST-derived SNPs

Consensus genetic linkage maps provide a genomic framework for quantitative trait loci identification, map-based cloning, assessment of genetic diversity, association mapping, and applied breeding in marker-assisted selection schemes. Among "orphan crops" with limited genomic resources such as cowpea [Vigna unguiculata (L.) Walp.] (2n = 2x = 22), the use of transcript-derived SNPs in genetic maps provides opportunities for automated genotyping and estimation of genome structure based on synteny analysis. Here, we report the development and validation of a high-throughput EST-derived SNP assay for cowpea, its application in consensus map building, and determination of synteny to reference genomes. SNP mining from 183,118 ESTs sequenced from 17 cDNA libraries yielded approximately 10,000 high-confidence SNPs from which an Illumina 1,536-SNP GoldenGate genotyping array was developed and applied to 741 recombinant inbred lines from six mapping populations. Approximately 90% of the SNPs were technically successful, providing 1,375 dependable markers. Of these, 928 were incorporated into a consensus genetic map spanning 680 cM with 11 linkage groups and an average marker distance of 0.73 cM. Comparison of this cowpea genetic map to reference legumes, soybean (Glycine max) and Medicago truncatula, revealed extensive macrosynteny encompassing 85 and 82%, respectively, of the cowpea map. Regions of soybean genome duplication were evident relative to the simpler diploid cowpea. Comparison with Arabidopsis revealed extensive genomic rearrangement with some conserved microsynteny. These results support evolutionary closeness between cowpea and soybean and identify regions for synteny-based functional genomics studies in legumes.