De novo transcriptomic analysis of cowpea (Vigna unguiculata L. Walp.) for genic SSR marker development

BackgroundCowpea [Vigna unguiculata (L.) Walp.] is one of the most important legumes in tropical and semi-arid regions. However, there is relatively little genomic information available for genetic research on and breeding of cowpea. The objectives of this study were to analyse the cowpea transcriptome and develop genic molecular markers for future genetic studies of this genus.ResultsApproximately 54 million high-quality cDNA sequence reads were obtained from cowpea based on Illumina paired-end sequencing technology and were de novo assembled to generate 47,899 unigenes with an N50 length of 1534 bp. Sequence similarity analysis revealed 36,289 unigenes (75.8%) with significant similarity to known proteins in the non-redundant (Nr) protein database, 23,471 unigenes (49.0%) with BLAST hits in the Swiss-Prot database, and 20,654 unigenes (43.1%) with high similarity in the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. Further analysis identified 5560 simple sequence repeats (SSRs) as potential genic molecular markers. Validating a random set of 500 SSR markers yielded 54 polymorphic markers among 32 cowpea accessions.ConclusionsThis transcriptomic analysis of cowpea provided a valuable set of genomic data for characterizing genes with important agronomic traits in Vigna unguiculata and a new set of genic SSR markers for further genetic studies and breeding in cowpea and related Vigna species.

[1]  Rajeev K. Varshney,et al.  Integration of novel SSR and gene-based SNP marker loci in the chickpea genetic map and establishment of new anchor points with Medicago truncatula genome , 2010, Theoretical and Applied Genetics.

[2]  T. Sharma,et al.  Exploring genetic variability within lentil (Lens culinaris Medik.) and across related legumes using a newly developed set of microsatellite markers , 2014, Molecular Biology Reports.

[3]  S. Gupta,et al.  Advances in genome mapping in orphan grain legumes of genus Vigna , 2013 .

[4]  W. Mao,et al.  Genetic Diversity Analysis of Faba Bean (Vicia faba L.) Based on EST-SSR Markers , 2011 .

[5]  S. Cannon,et al.  A Comprehensive Transcriptome Assembly of Pigeonpea (Cajanus cajan L.) using Sanger and Second-Generation Sequencing Platforms , 2012, Molecular plant.

[6]  R. Varshney,et al.  Comparative assessment of EST-SSR, EST-SNP and AFLP markers for evaluation of genetic diversity and conservation of genetic resources using wild, cultivated and elite barleys , 2007 .

[7]  M. Fisher,et al.  Illuminating Choices for Library Prep: A Comparison of Library Preparation Methods for Whole Genome Sequencing of Cryptococcus neoformans Using Illumina HiSeq , 2014, PloS one.

[8]  M. Blair,et al.  BAC-End Microsatellites from Intra and Inter-Genic Regions of the Common Bean Genome and Their Correlation with Cytogenetic Features , 2014, PloS one.

[9]  Rajeev K. Varshney,et al.  Genome sequence of mungbean and insights into evolution within Vigna species , 2014, Nature Communications.

[10]  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.

[11]  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.

[12]  M. Blair,et al.  Transcriptome Sequencing of Mung Bean (Vigna radiate L.) Genes and the Identification of EST-SSR Markers , 2015, PloS one.

[13]  C. Town.  Annotating the genome of Medicago truncatula. , 2006, Current opinion in plant biology.

[14]  Andreas Graner,et al.  Genic microsatellite markers in plants: features and applications. , 2005, Trends in biotechnology.

[15]  S. Jackson,et al.  Integration of physical and genetic maps of common bean through BAC-derived microsatellite markers , 2010, BMC Genomics.

[16]  S. Tabata,et al.  Genome Sequencing and Genome Resources in Model Legumes , 2007, Plant Physiology.

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

[18]  M. Timko,et al.  Genetic relationship of cowpea (Vigna unguiculata) varieties from Senegal based on SSR markers. , 2012, Genetics and molecular research : GMR.

[19]  M. Blair,et al.  Development of microsatellite markers for common bean (Phaseolus vulgaris L.) based on screening of non-enriched, small-insert genomic libraries. , 2009, Genome.

[20]  James K. Hane,et al.  Draft genome sequence of chickpea (Cicer arietinum) provides a resource for trait improvement , 2013, Nature Biotechnology.

[21]  T. Close,et al.  An improved genetic linkage map for cowpea (Vigna unguiculata L.) combining AFLP, RFLP, RAPD, biochemical markers, and biological resistance traits. , 2002, Genome.

[22]  M. Melotto,et al.  Development, characterization, and comparative analysis of polymorphism at common bean SSR loci isolated from genic and genomic sources. , 2007, Genome.

[23]  Hongwei Jiang,et al.  Identification and characterization of SSRs from soybean (Glycine max) ESTs , 2012, Molecular Biology Reports.

[24]  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.

[25]  J. Souframanien,et al.  De novo Assembly, Characterization of Immature Seed Transcriptome and Development of Genic-SSR Markers in Black Gram [Vigna mungo (L.) Hepper] , 2015, PloS one.

[26]  S. Tabata,et al.  Lotus japonicus as a platform for legume research. , 2006, Current opinion in plant biology.

[27]  S. Gupta,et al.  Development of unigene-derived SSR markers in cowpea (Vigna unguiculata) and their transferability to other Vigna species. , 2010, Genome.

[28]  Jin Wang,et al.  De novo assembly and Characterisation of the Transcriptome during seed development, and generation of genic-SSR markers in Peanut (Arachis hypogaea L.) , 2012, BMC Genomics.

[29]  Fengxiang Yin,et al.  Development of Gene-Based SSR Markers in Rice Bean (Vigna umbellata L.) Based on Transcriptome Data , 2016, PloS one.

[30]  B. Leikind,et al.  An Assessment of , 2011 .

[31]  T. Sharma,et al.  RESEARCH ARTICLE Open Access Development of genic-SSR markers by deep transcriptome sequencing in pigeonpea [Cajanus cajan (L.) Millspaugh] , 2022 .

[32]  N. Sethy,et al.  Development of chickpea EST-SSR markers and analysis of allelic variation across related species , 2009, Theoretical and Applied Genetics.

[33]  Linhai Wang,et al.  Characterization of the sesame (Sesamum indicum L.) global transcriptome using Illumina paired-end sequencing and development of EST-SSR markers , 2011, BMC Genomics.

[34]  Colin N. Dewey,et al.  De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis , 2013, Nature Protocols.

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

[36]  M. Gerstein,et al.  RNA-Seq: a revolutionary tool for transcriptomics , 2009, Nature Reviews Genetics.

[37]  S. Gupta,et al.  Development and characterization of genic SSR markers for mungbean (Vigna radiata (L.) Wilczek) , 2013, Euphytica.

[38]  B. Mou,et al.  Genetic Diversity and Population Structure of Cowpea (Vigna unguiculata L. Walp) , 2016, PloS one.

[39]  R. Wing,et al.  Gene-based SSR markers for common bean (Phaseolus vulgaris L.) derived from root and leaf tissue ESTs: an integration of the BMc series , 2011, BMC Plant Biology.

[40]  Juan Miguel García-Gómez,et al.  BIOINFORMATICS APPLICATIONS NOTE Sequence analysis Manipulation of FASTQ data with Galaxy , 2005 .

[41]  K. Hilu,et al.  Microsatellites and RAPD Markers to Study Genetic Relationships Among Cowpea Breeding Lines and Local Varieties in Senegal , 2004, Genetic Resources and Crop Evolution.

[42]  S. Isobe,et al.  Detection of Genome Donor Species of Neglected Tetraploid Crop Vigna reflexo-pilosa (Créole Bean), and Genetic Structure of Diploid Species Based on Newly Developed EST-SSR Markers from Azuki Bean (Vigna angularis) , 2014, PloS one.

[43]  M. Timko,et al.  Assessment of the genetic diversity in cowpea (Vigna unguiculata L. Walp.) germplasm from Ghana using simple sequence repeat markers , 2010, Plant Genetic Resources.

[44]  The transferability and polymorphism of mung bean SSR markers in rice bean germplasm , 2015, Molecular Breeding.

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

[46]  M. Nelson,et al.  Next generation DNA sequencing technology delivers valuable genetic markers for the genomic orphan legume species, Bituminaria bituminosa , 2011, BMC Genetics.

[47]  N. Young,et al.  Comparative genome analysis of mungbean (Vigna radiata L. Wilczek) and cowpea (V. unguiculata L. Walpers) using RFLP mapping data , 1993, Theoretical and Applied Genetics.

[48]  Akiyasu C. Yoshizawa,et al.  KAAS: an automatic genome annotation and pathway reconstruction server , 2007, Environmental health perspectives.

[49]  Liping Liu,et al.  Development and Validation of EST-SSR Markers from the Transcriptome of Adzuki Bean (Vigna angularis) , 2015, PloS one.

[50]  Ju-Kyung Yu,et al.  Nonrandom distribution and frequencies of genomic and EST-derived microsatellite markers in rice, wheat, and barley , 2005, BMC Genomics.

[51]  M. Blair,et al.  Assessment of genetic diversity and population structure of mung bean (Vigna radiata) germplasm using EST-based and genomic SSR markers. , 2015, Gene.

[52]  T. Close,et al.  Genomic Tools in Cowpea Breeding Programs: Status and Perspectives , 2016, Front. Plant Sci..

[53]  D. Diouf.  Recent advances in cowpea [ Vigna unguiculata (L.) Walp.] “omics” research for genetic improvement , 2011 .

[54]  P. Gepts,et al.  Genetic diversity in cowpea [Vigna unguiculata (L.) Walp.] as revealed by RAPD markers , 2004, Genetic Resources and Crop Evolution.