Genome-wide single nucleotide polymorphism discovery and validation in adzuki bean

Adzuki bean, also known as red bean (Vigna angularis), with 2n = 22 chromosomes, is an important legume crop in East Asian countries, including China, Japan, and Korea. For single nucleotide polymorphism (SNP) discovery, we used Vigna accessions, V. angularis IT213134 and its wild relative V. nakashimae IT178530, because of the lack of DNA sequence polymorphism in the cultivated species. Short read sequences of IT213134 and IT178530 of approximately 37 billion and 35 billion bp were produced using the Illumina HiSeq 2000 system to a sequencing depth of 61.5× and 57.7×, respectively. After de novo assembly was carried out with trimmed HiSeq reads from IT213134, 98,441 contigs of various sizes were produced with N50 of 13,755 bp. Using Burrows–Wheeler Aligner software, trimmed short reads of V. nakashimae IT178530 were successfully mapped to IT213134 contigs. All sequence variations at the whole-genome level were examined between the two Vigna species. Of the 1,565,699 SNPs, 59.4 % were transitions and 40.6 % were transversions. A total of 213,758 SNPs, consisting of 122,327 non-synonymous and 91,431 synonymous SNPs, were identified in coding sequences. For SNP validation, 96 SNPs in the genic region were chosen from among IT213134 contigs longer than 10 kb. Of these 96 SNPs, 88 were confirmed by Sanger sequencing of 10 adzuki bean genotypes from various geographic origins as well as IT213134 and its wild relative IT178530. These genome-wide SNP markers will enrich the existing Vigna resources and, specifically, could be of value for constructing a genetic map and evaluating the genetic diversity of adzuki bean.

[1]  J. Batley,et al.  Single nucleotide polymorphism discovery from wheat next-generation sequence data. , 2012, Plant biotechnology journal.

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

[3]  T. Lumpkin,et al.  Azuki bean: botany, production and uses. , 1994 .

[4]  Nigel Maxted,et al.  The Asian Vigna:: Genus Vigna subgenus Ceratotropis Genetic Resources , 2003 .

[5]  Liisa Holm,et al.  The Pfam protein families database , 2009, Nucleic acids research.

[6]  L. Holm,et al.  The Pfam protein families database , 2005, Nucleic Acids Res..

[7]  Enrique Blanco,et al.  Computational gene annotation in new genome assemblies using GeneID. , 2009, Methods in molecular biology.

[8]  Robert D. Finn,et al.  The Pfam protein families database , 2004, Nucleic Acids Res..

[9]  Yasukazu Nakamura,et al.  SNP Discovery and Linkage Map Construction in Cultivated Tomato , 2010, DNA research : an international journal for rapid publication of reports on genes and genomes.

[10]  Rolf Apweiler,et al.  InterProScan: protein domains identifier , 2005, Nucleic Acids Res..

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

[12]  Steven J. M. Jones,et al.  Abyss: a Parallel Assembler for Short Read Sequence Data Material Supplemental Open Access , 2022 .

[13]  T. Ishii,et al.  A genetic linkage map of azuki bean constructed with molecular and morphological markers using an interspecific population (Vigna angularis x V. nakashimae) , 1996, Theoretical and Applied Genetics.

[14]  R. Gaur,et al.  EST-derived genic molecular markers: development and utilization for generating an advanced transcript map of chickpea , 2012, Theoretical and Applied Genetics.

[15]  Kyujung Van,et al.  Genome-wide mapping of NBS-LRR genes and their association with disease resistance in soybean , 2012, BMC Plant Biology.