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.

[1]  A. Korte,et al.  The advantages and limitations of trait analysis with GWAS: a review , 2013, Plant Methods.

[2]  S. Salzberg,et al.  Versatile and open software for comparing large genomes , 2004, Genome Biology.

[3]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

[4]  D. Botstein,et al.  Construction of a genetic linkage map in man using restriction fragment length polymorphisms. , 1980, American journal of human genetics.

[5]  Stefano Lonardi,et al.  Accurate Construction of Consensus Genetic Maps via Integer Linear Programming , 2011, IEEE/ACM Transactions on Computational Biology and Bioinformatics.

[6]  Richard Durbin,et al.  Sequence analysis Fast and accurate short read alignment with Burrows – Wheeler transform , 2009 .

[7]  J. Thomson.  The role of biotechnology for agricultural sustainability in Africa , 2008, Philosophical Transactions of the Royal Society B: Biological Sciences.

[8]  Marta T. Bokowiec,et al.  Sequencing and analysis of the gene-rich space of cowpea , 2008, BMC Genomics.

[9]  S. Jackson,et al.  Highly distinct chromosomal structures in cowpea (Vigna unguiculata), as revealed by molecular cytogenetic analysis , 2016, Chromosome Research.

[10]  Vinod Kumar,et al.  A genome-wide SNP scan accelerates trait-regulatory genomic loci identification in chickpea , 2015, Scientific Reports.

[11]  Arun Prabhu Dhanapal,et al.  Genome-Wide Association Study of Ureide Concentration in Diverse Maturity Group IV Soybean [Glycine max (L.) Merr.] Accessions , 2015, G3: Genes, Genomes, Genetics.

[12]  Jian Wang,et al.  SOAPdenovo2: an empirically improved memory-efficient short-read de novo assembler , 2012, GigaScience.

[13]  D. Verma,et al.  Identification of “nodule-specific” host proteins (nodulins) involved in the development of Rhizobium-Legume symbiosis , 1980, Cell.

[14]  Huanming Yang,et al.  Draft genome sequence of pigeonpea (Cajanus cajan), an orphan legume crop of resource-poor farmers , 2011, Nature Biotechnology.

[15]  G. Evanno,et al.  Detecting the number of clusters of individuals using the software structure: a simulation study , 2005, Molecular ecology.

[16]  Carolyn Thomas,et al.  High-throughput fingerprinting of bacterial artificial chromosomes using the snapshot labeling kit and sizing of restriction fragments by capillary electrophoresis. , 2003, Genomics.

[17]  P. Green,et al.  Base-calling of automated sequencer traces using phred. I. Accuracy assessment. , 1998, Genome research.

[18]  Timothy J. Close,et al.  A consensus genetic map of cowpea [Vigna unguiculata (L) Walp.] and synteny based on EST-derived SNPs , 2009, Proceedings of the National Academy of Sciences.

[19]  P. Salimath,et al.  Genetic Diversity in Cowpea , 2006 .

[20]  Madeline Fisher,et al.  B.B. Singh and His Quest to Make Cowpea the Food Legume of the 21st Century , 2014 .

[21]  J. Batley,et al.  Discovery of Single Nucleotide Polymorphisms in Complex Genomes Using SGSautoSNP , 2012, Biology.

[22]  A. Challinor,et al.  Climate variability and vulnerability to climate change: a review , 2014, Global change biology.

[23]  Gabor T. Marth,et al.  Haplotype-based variant detection from short-read sequencing , 2012, 1207.3907.

[24]  Mitchell R. Lucas,et al.  Gene Pools and the Genetic Architecture of Domesticated Cowpea , 2013 .

[25]  N. Talbot,et al.  The glyoxylate cycle is required for temporal regulation of virulence by the plant pathogenic fungus Magnaporthe grisea , 2003, Molecular microbiology.

[26]  Markers for breeding heat-tolerant cowpea , 2013, Molecular Breeding.

[27]  Rod A Wing,et al.  A reference genome for common bean and genome-wide analysis of dual domestications , 2014, Nature Genetics.

[28]  P. Donnelly,et al.  Inference of population structure using multilocus genotype data. , 2000, Genetics.

[29]  Sergey I. Nikolenko,et al.  SPAdes: A New Genome Assembly Algorithm and Its Applications to Single-Cell Sequencing , 2012, J. Comput. Biol..

[30]  T. Close,et al.  Genetic Architecture of Delayed Senescence, Biomass, and Grain Yield under Drought Stress in Cowpea , 2013, PloS one.

[31]  Rajeev K. Varshney,et al.  Genome sequencing of adzuki bean (Vigna angularis) provides insight into high starch and low fat accumulation and domestication , 2015, Proceedings of the National Academy of Sciences.

[32]  T. Close,et al.  Genomic regions, cellular components and gene regulatory basis underlying pod length variations in cowpea (V. unguiculata L. Walp) , 2016, Plant biotechnology journal.

[33]  M. Nei,et al.  F‐statistics and analysis of gene diversity in subdivided populations , 1977, Annals of human genetics.

[34]  Timothy J. Close,et al.  Cowpea–Soybean Synteny Clarified through an Improved Genetic Map , 2011 .

[35]  Gianfranco Ciardo,et al.  Combinatorial Pooling Enables Selective Sequencing of the Barley Gene Space , 2013, PLoS Comput. Biol..

[36]  Stefano Lonardi,et al.  A Graph-Theoretical Approach to the Selection of the Minimum Tiling Path from a Physical Map , 2013, IEEE/ACM Transactions on Computational Biology and Bioinformatics.

[37]  Stefano Lonardi,et al.  Efficient and Accurate Construction of Genetic Linkage Maps from the Minimum Spanning Tree of a Graph , 2008, PLoS genetics.

[38]  T. Close,et al.  Genetic mapping, synteny, and physical location of two loci for Fusarium oxysporum f. sp. tracheiphilum race 4 resistance in cowpea [Vignaunguiculata (L.) Walp] , 2013, Molecular Breeding.

[39]  T. Itoh,et al.  The Vigna Genome Server, 'VigGS': A Genomic Knowledge Base of the Genus Vigna Based on High-Quality, Annotated Genome Sequence of the Azuki Bean, Vigna angularis (Willd.) Ohwi & Ohashi. , 2016, Plant & cell physiology.

[40]  P Green,et al.  Base-calling of automated sequencer traces using phred. II. Error probabilities. , 1998, Genome research.

[41]  M. DePristo,et al.  The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. , 2010, Genome research.

[42]  J. Dvorak,et al.  A BAC-based physical map of Brachypodium distachyon and its comparative analysis with rice and wheat , 2009, BMC Genomics.

[43]  A. C. Brasileiro-Vidal,et al.  Intra- and interchromosomal rearrangements between cowpea [Vigna unguiculata (L.) Walp.] and common bean (Phaseolus vulgaris L.) revealed by BAC-FISH , 2015, Chromosome Research.

[44]  Edward S. Buckler,et al.  TASSEL: software for association mapping of complex traits in diverse samples , 2007, Bioinform..

[45]  Steven J. M. Jones,et al.  Circos: an information aesthetic for comparative genomics. , 2009, Genome research.

[46]  Carol Soderlund,et al.  FPC: a system for building contigs from restriction fingerprinted clones , 1997, Comput. Appl. Biosci..

[47]  M. Stephens,et al.  Inference of population structure using multilocus genotype data: dominant markers and null alleles , 2007, Molecular ecology notes.

[48]  Paul J. Rushton,et al.  CGKB: an annotation knowledge base for cowpea (Vigna unguiculata L.) methylation filtered genomic genespace sequences , 2007, BMC Bioinformatics.

[49]  Mitchell R. Lucas,et al.  Association Studies and Legume Synteny Reveal Haplotypes Determining Seed Size in Vigna unguiculata , 2013, Front. Plant Sci..

[50]  Gianfranco Ciardo,et al.  When Less is More: “Slicing” Sequencing Data Improves Read Decoding Accuracy and De Novo Assembly Quality , 2015, bioRxiv.