SummaryCowpea (Vigna unguiculata [L.] Walp.) is a major crop for worldwide food and nutritional security, especially in sub‐Saharan Africa, that is resilient to hot and drought‐prone environments. An assembly of the single‐haplotype inbred genome of cowpea IT97K‐499‐35 was developed by exploiting the synergies between single‐molecule real‐time sequencing, optical and genetic mapping, and an assembly reconciliation algorithm. A total of 519 Mb is included in the assembled sequences. Nearly half of the assembled sequence is composed of repetitive elements, which are enriched within recombination‐poor pericentromeric regions. A comparative analysis of these elements suggests that genome size differences between Vigna species are mainly attributable to changes in the amount of Gypsy retrotransposons. Conversely, genes are more abundant in more distal, high‐recombination regions of the chromosomes; there appears to be more duplication of genes within the NBS‐LRR and the SAUR‐like auxin superfamilies compared with other warm‐season legumes that have been sequenced. A surprising outcome is the identification of an inversion of 4.2 Mb among landraces and cultivars, which includes a gene that has been associated in other plants with interactions with the parasitic weed Striga gesnerioides. The genome sequence facilitated the identification of a putative syntelog for multiple organ gigantism in legumes. A revised numbering system has been adopted for cowpea chromosomes based on synteny with common bean (Phaseolus vulgaris). An estimate of nuclear genome size of 640.6 Mbp based on cytometry is presented.

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.

Cowpea (Vigna unguiculata [L.] Walp.) is a major crop for worldwide food and nutritional security, especially in sub-Saharan Africa, that is resilient to hot and drought-prone environments. An assembly of the singlehaplotype inbred genome of cowpea IT97K-499-35 was developed by exploiting the synergies between single-molecule real-time sequencing, optical and genetic mapping, and an assembly reconciliation algorithm. A total of 519 Mb is included in the assembled sequences. Nearly half of the assembled sequence is composed of repetitive elements, which are enriched within recombination-poor pericentromeric regions. A comparative analysis of these elements suggests that genome size differences between Vigna species are mainly attributable to changes in the amount of Gypsy retrotransposons. Conversely, genes are more abundant in more distal, high-recombination regions of the chromosomes; there appears to be more duplication of genes within the NBS-LRR and the SAUR-like auxin superfamilies compared with other warm-season legumes that have been sequenced. A surprising outcome is the identification of an inversion of 4.2 Mb among landraces and cultivars, which includes a gene that has been associated in other plants with interactions with the parasitic weed Striga gesnerioides. The genome sequence facilitated the identification of a putative syntelog for multiple organ gigantism in legumes. A revised numbering system has been adopted for cowpea chromosomes based on synteny with common bean (Phaseolus vulgaris). An estimate of nuclear genome size of 640.6 Mbp based on cytometry is presented.

Key messageGenome resolution of a major QTL associated with theRklocus in cowpea for resistance to root-knot nematodes has significance for plant breeding programs and R gene characterization.Abstract Cowpea (Vigna unguiculata L. Walp.) is a susceptible host of root-knot nematodes (Meloidogyne spp.) (RKN), major plant-parasitic pests in global agriculture. To date, breeding for host resistance in cowpea has relied on phenotypic selection which requires time-consuming and expensive controlled infection assays. To facilitate marker-based selection, we aimed to identify and map quantitative trait loci (QTL) conferring the resistance trait. One recombinant inbred line (RIL) and two F2:3 populations, each derived from a cross between a susceptible and a resistant parent, were genotyped with genome-wide single nucleotide polymorphism (SNP) markers. The populations were screened in the field for root-galling symptoms and/or under growth-chamber conditions for nematode reproduction levels using M. incognita and M. javanica biotypes. One major QTL was mapped consistently on linkage group VuLG11 of each population. By genotyping additional cowpea lines and near-isogenic lines derived from conventional backcrossing, we confirmed that the detected QTL co-localized with the genome region associated with the Rk locus for RKN resistance that has been used in conventional breeding for many decades. This chromosomal location defined with flanking markers will be a valuable target in marker-assisted breeding and for positional cloning of genes controlling RKN resistance.Electronic supplementary materialThe online version of this article (doi:10.1007/s00122-015-2611-0) contains supplementary material, which is available to authorized users.

Cowpea is one of the most important food and forage legumes in drier regions of the tropics and subtropics. However, cowpea yield worldwide is markedly below the known potential due to abiotic and biotic stresses, including parasitism by root-knot nematodes (Meloidogyne spp., RKN). Two resistance genes with dominant effect, Rk and Rk 2 , have been reported to provide resistance against RKN in cowpea. Despite their description and use in breeding for resistance to RKN and particularly genetic mapping of the Rk locus, the exact genes conferring resistance to RKN remain unknown. In the present work, QTL mapping using recombinant inbred line (RIL) population 524B x IT84S-2049 segregating for a newly mapped locus and analysis of the transcriptome changes in two cowpea near-isogenic lines (NIL) were used to identify candidate genes for Rk and the newly mapped locus. A major QTL, designated QRk-vu9.1, associated with resistance to Meloidogyne javanica reproduction, was detected and mapped on linkage group LG9 at position 13.37 cM using egg production data. Transcriptome analysis on resistant and susceptible NILs 3 and 9 days after inoculation revealed up-regulation of 109 and 98 genes and down-regulation of 110 and 89 genes, respectively, out of 19,922 unique genes mapped to the common bean reference genome. Among the differentially expressed genes, four and nine genes were found within the QRk-vu9.1 and QRk-vu11.1 QTL intervals, respectively. Six of these genes belong to the TIR-NBS-LRR family of resistance genes and three were upregulated at one or more time-points. Quantitative RT-PCR validated gene expression to be positively correlated with RNA-seq expression pattern for eight genes. Future functional analysis of these cowpea genes will enhance our understanding of Rk-mediated resistance and identify the specific gene responsible for the resistance.

Background and Aims Endoparasitic root-knot nematodes (RKNs) (Meloidogyne spp.) cause considerable losses in banana (Musa spp.), with Meloidogyne incognita a predominant species in Cavendish sub-group bananas. This study investigates the root transcriptome in Musa acuminata genotypes 4297-06 (AA) and Cavendish Grande Naine (CAV; AAA) during early compatible interactions with M. incognita.Methods Roots were analysed by brightfield light microscopy over a 35 d period to examine nematode penetration and morphological cell transformation. RNA samples were extracted 3, 7 and 10 days after inoculation (DAI) with nematode J2 juveniles, and cDNA libraries were sequenced using lllumina HiSeq technology. Sequences were mapped to the M. acuminata ssp. malaccensis var. Pahang genome sequence, differentially expressed genes (DEGs) identified and transcript representation determined by gene set enrichment and pathway mapping.Key Results Microscopic analysis revealed a life cycle of M. incognita completing in 24 d in CAV and 27 d in 4279-06. Comparable numbers of DEGs were up-and downregulated in each genotype, with potential involvement of many in early host defence responses involving reactive oxygen species and jasmonate/ethylene signalling. DEGs revealed concomitant auxin metabolism and cell wall modification processes likely to be involved in giant cell formation. Notable transcripts related to host defence included those coding for leucine-rich repeat receptorlike serine/threonine-protein kinases, peroxidases, thaumatin-like pathogenesis-related proteins, and DREB, ERF, MYB, NAC and WRKY transcription factors. Transcripts related to giant cell development included indole acetic acid-amido synthetase GH3.8 genes, involved in auxin metabolism, as well as genes encoding expansins and hydrolases, involved in cell wall modification.Conclusions Expression analysis in M. acuminata during compatible interactions with RKNs provides insights into genes modulated during infection and giant cell formation. Increased understanding of both defence responses to limit parasitism during compatible interactions and effector-targeted host genes in this complex interaction will facilitate the development of genetic improvement measures for RKNs.

Cowpea (Vigna unguiculata L. Walp) is an affordable source of protein and strategic legume crop for food security in Africa and other developing regions; however, damage from infection by root-knot nematodes (RKN) suppresses cowpea yield. The deployment through breeding of resistance gene Rk in cowpea cultivars has provided protection to cowpea growers worldwide for many years. However, occurrence of more aggressive nematode isolates threatens the effectiveness of this monogenic resistance. A cowpea germplasm collection of 48 genotypes representing the cowpea gene-pool from Eastern and Southern Africa (cowpea has two major pools of genetic resources-Western Africa and Eastern/Southern Africa) was screened in replicated experiments under field, greenhouse and controlled-growth conditions to identify resistance to RKN, to determine the spectrum of resistance to RKN, the relative virulence (VI) among RKN species and isolates, and the relationship between root-galling (RG) and egg-mass production (EM). Analysis of variance of data for RG and EM per root system identified seven genotypes with broad-based resistance to Meloidogyne javanica (Mj), avirulent (Avr-Mi), and virulent (Mi) M. incognita isolates. Two of the 48 genotypes exhibited specific resistance to both Mi isolates. Most of the genotypes were resistant to Avr-Mi indicating predominance of Rk gene in the collection. Based on RG data, both Mj (VI = 50%) and Mi (VI = 42%) were fourfold more virulent than Avr-Mi (VI = 12%). Resistant genotypes had more effective resistance than the Rk-based resistance in cowpea genotype CB46 against Mj and Mi. Root-galling was correlated across isolates (Avr-Mi/Mj: r = 0.72; Mi/Mj: r = 0.98), and RG was correlated with EM (r = 0.60), indicating resistance to RG and EM is under control by the same genetic factors. These new sources of resistance identified in cowpea gene-pool two provide valuable target traits for breeders to improve cowpea production on RKN-infested fields.

The utility of combining simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) marker genotyping was determined for genetically mapping a novel aphid (Aphis craccivora) resistance locus in cowpea breeding line SARC 1‐57‐2 and for introgressing the resistance into elite cultivars by marker‐assisted backcrossing (MABC). The locus was tagged with codominant SSR marker CP 171F/172R with a recombination fraction of 5.91% in an F2 population from ‘Apagbaala’ x SARC 1‐57‐2. A SNP‐genotyped biparental recombinant inbred line population was genotyped for CP 171F/172R, which was mapped to position 11.5 cM on linkage group (LG) 10 (physical position 30.514 Mb on chromosome Vu10). Using CP 171F/172R for foreground selection and a KASP‐SNP‐based marker panel for background selection in MABC, the resistance from SARC 1‐57‐2 was introduced into elite susceptible cultivar ‘Zaayura’. Five BC4F3 lines of improved ‘Zaayura’ that were isogenic except for the resistance locus region had phenotypes similar to SARC 1‐57‐2. This study identified a novel aphid resistance locus and demonstrated the effectiveness of integrating SSR and SNP markers for trait mapping and marker‐assisted breeding.