{'fraction_non_alphanumeric': 0.04502763367578747, 'fraction_numerical': 0.028311802800664563, 'mean_word_length': 4.601899335232669, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 2, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 2, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
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)

Vigna reflexo-pilosa, which includes a neglected crop, is the only one tetraploid species in genus Vigna. The ancestral species that make up this allotetraploid species have not conclusively been identified, although previous studies suggested that a donor genome of V. reflexo-pilosa is V. trinervia. In this study, 1,429 azuki bean EST-SSR markers were developed of which 38 EST-SSR primer pairs that amplified one product in diploid species and two discrete products in tetraploid species were selected to analyze 268 accessions from eight taxa of seven Asian Vigna species including V. reflexo-pilosa var. glabra, V. reflexo-pilosa var. reflexo-pilosa, V. exilis, V. hirtella, V. minima, V. radiata var. sublobata, V. tenuicaulis and V. trinervia to identify genome donor of V. reflexo-pilosa. Since both diploid and tetraploid species were analyzed and each SSR primer pair detected two loci in the tetraploid species, we separated genomes of the tetraploid species into two different diploid types, viz. A and B. In total, 445 alleles were detected by 38 EST-SSR markers. The highest gene diversity was observed in V. hirtella. By assigning the discrete PCR products of V. reflexo-pilosa into two distinguished genomes, we were able to identify the two genome donor parents of créole bean. Phylogenetic and principal coordinate analyses suggested that V. hirtella is a species complex and may be composed of at least three distinct taxa. Both analyses also clearly demonstrated that V. trinervia and one taxon of V. hirtella are the genome donors of V. reflexo-pilosa. Gene diversity indicates that the evolution rate of EST-SSRs on genome B of créole bean might be faster than that on genome A. Species relationship among the Vigna species in relation to genetic data, morphology and geographical distribution are presented.

Cultivated and wild forms of créole bean are classified as V. reflexo-pilosa var. glabra and V. reflexo-pilosa var. reflexo-pilosa, respectively [6]. Wild créole bean is widely distributed in East, Southeast and South Asia, and across the islands from the west to the north Pacific islands. It is also found in Papua New Guinea and northern Australia [7] [8] [9]. The cultivated créole bean was formerly recognized as a glabrous variety of mungbean, V. radiata var. glabra [1]. Then, it was treated as a distinct species, V. glabrescens [3]. It is differentiated from its wild progenitor principally by thick glabrous stem and erect growth habit and reported to be cultivated as pulse in Vietnam and the Philippines or as forage in India, Mauritius and Tanzania [2]. This crop shows resistance to several insect pests and diseases such as bruchids, bean fly, powdery mildew, and cucumber mosaic virus [10], and is partially cross-compatible with mungbean [11]. Thus, créole bean has potential to be a gene source for breeding other Vigna crops. In addition, cultivated créole bean could be considered as a novel crop for the future and wild créole bean as a wider genepool to improve cultivated créole bean.
The origin and genome donors of créole bean have been the subject of debate. Based on the study on isozyme banding patterns and interspecific hybridization among the Ceratotropis species, Egawa et al [10] [12] proposed that V. reflexo-pilosa var. reflexopilosa is the descendant of natural interspecific hybridization between V. trinervia and V. minima followed by spontaneous chromosome doubling. An accession used as V. minima in Egawa et al. [10] [12] was correctly identified as V. hirtella by Konarev et al. [13], and was used in the present study (No. 21, JP108851, from Malaysia). Proteinase inhibitors polymorphism study in Asian Vigna by Konarev et al. [13] supported that one genome donor of créole bean is V. trinervia and suggested that the other genome donor is most likely V. hirtella or its closely related species. The phylogenetic studies based on plastid sequence data also supported that V. trinervia is a genome donor of créole bean [14] [15]. In contrast, rDNA-ITS sequence variation showed high similarity between créole bean and V. exilis, V. hirtella and V. umbellata [16], suggesting that one of these species is the genome donor of V. reflexo-pilosa.
The objective of this study was to determine the putative genome donor species of tetraploid créole bean using EST-SSR markers. To do so, we developed EST-SSR from azuki bean and used them to analyze accessions of créole bean and other Vigna species that are candidate genome donor species.

Plant Materials
A Japanese azuki bean cultivar 'Erimo-shouzu' (V. angularis var. angularis, accession no. JP37752) obtained from the Genebank, National Institute of Agrobiological Sciences (NIAS), Tsukuba, Japan was used for development of EST-SSR markers. Eight accessions consisting of four major Vigna crop species; azuki bean, rice bean, black gram and mungbean were used for assessing transferability and polymorphisms of the EST-SSR markers (Table 1). Two hundred and sixty-eight accessions from eight taxa of seven Asian Vigna species (Figure 1, Table S1) including 7 of V. reflexo-pilosa var. glabra, 51 of V. reflexo-pilosa var. reflexopilosa, 13 of V. exilis, 47 of V. hirtella, 49 of V. minima, 13 of V. radiata var. sublobata, 42 of V. tenuicaulis, and 46 of V. trinervia were used to analyze the genome origin of V. reflexo-pilosa.

Development of Azuki Bean EST-SSR Markers
Total RNA was extracted from seedlings and young pods of Erimo-shouzu using Plant RNA Purification Reagent (Invitrogen, CA, USA). Purification of polyadenylated RNA and conversion to cDNA were performed as described by [17]. Synthesized cDNA was resolved by 1% agarose gel electrophoresis, and fragments ranging from 1 to 3 kb were recovered. The recovered fragments were cloned into the Eco RI-Xho I site of the pBluescript II SKplasmid vector (Stratagene, CA, USA) and introduced into the E. coli ElectroTen-Blue strain (Stratagene, CA, USA) by electroporation. To generate ESTs, plasmid DNAs were amplified from the colonies using TempliPhi (GE Healthcare UK Ltd, Buckinghamshire, England) and subjected to sequencing using a BigDye Terminator Cycle Sequencing Ready Reaction Kit (Applied Biosystems, CA, USA). The reaction mixtures were run on an automated ABI PRISM 3730 DNA Analyzer (Applied Biosystems).
Sequencing chromatograms were converted into nucleotide bases with Phred [18] [19], and the sequences derived from the vector and linkers were removed with CROSSMATCH [19]. The EST reads were quality-trimmed with TRIM2 [20] using the Phred quality score $20, and ambiguous regions including more than ten X or N bases were trimmed. Contiguous, high-quality reads $100 bp were submitted to the DDBJ/EMBL/GenBank databases under the accession numbers HX939204 to HX950377 (11,174 entries). The PHRAP program with default parameters was used to cluster and identify non-redundant azuki bean ESTs [19].
Simple sequence repeats (SSRs) $15 nucleotides in length, which contained all possible combinations of di-nucleotide (NN), tri-nucleotide (NNN), and tetra-nucleotide (NNNN) repeats, were identified from the non-redundant azuki bean ESTs using fuzznuc program in EMBOSS [21] for SSRs within two mismatches. Primer pairs for amplification of SSR-containing regions were designed based on the flanking sequences of each SSR with the aid of the Primer3 [22] so that the amplified fragment sizes were between 90 bp and 300 bp in length. The newly developed markers were designated as VES (Vigna EST-derived SSR) markers.

EST-SSR Marker Analysis
Total genomic DNA was extracted from young leaves of each Vigna accession using the method described by [23] with a slight modification. The DNA was quantified against a lambda DNA on 1.0% agarose gel stained with ethidium bromide and diluted to 5 ng/ml for PCR amplification.
Amplification, transferability and polymorphisms of all the EST-SSR markers were confirmed using eight accessions, two each from azuki bean complex, blackgram, rice bean and mungbean (Table 1). Five microliters of PCR reaction mixture including 0.5 ng of total genomic DNA, 0.02 U BIOTAQ DNA polymerase (BIOLINE, UK), 1x PCR buffer (BIOLINE, UK), 3 mM MgCl 2 , 0.2 mM dNTPs and 2 pmol of the forward and reverse EST-SSR primers. The PCR thermal cycling was performed as follows: 94uC for 1 min; 3 cycles of 94uC for 30 s and 70uC for 30 s, followed by 3 rounds of the same program in which the annealing and extension temperatures were decreased by 2uC every 3 cycles; 3 cycles of 94uC for 30 s, 62uC for 30 s and 72uC for 30 s, followed by 2 rounds of the same program in which the annealing temperatures were decreased by 2uC every 2 cycles; 30 cycles of 94uC for 30 s, 55uC for 30 s and 72uC for 30 s, and a final cycle at 72uC for 10 min. The PCR products were separated by 10% polyacrylamide gel electrophoresis in tris-borate-ethylene diamine tetraacetic acid (TBE) buffer according to the standard protocol, and banding patterns in the gel were stored as pictures. The characteristics of PCR products such as intensity, banding pattern and size range were recorded for each accession and EST-SSR markers. Polymorphism between wild and cultivated accessions in each species was recorded. Out of 1,429 azuki bean EST-SSR markers, 175 markers with good transferability in all the four major Vigna crop species were labeled with fluorescent dyes and further examined for polymorphism in a panel of eight taxa of Asian Vigna consisting of 16 accessions, two accessions for each species randomly chosen from the 268 accessions (Table S1). The 59-end of the reverse primers were labeled with one of the following fluorescent dyes; 6-FAM (blue), HEX (green), and NED (yellow) (Applied Biosystems). Five microliters of PCR reaction mixture contained 5 ng of genomic DNA, 16 QIAGEN Multiplex PCR Master Mix, and 5 pmol of the forward and reverse primers. PCR reactions were performed in a GeneAmp PCR System 9700 (Applied Biosystems). The PCR thermal cycling was programmed as follows: 95uC for 15 min followed by 40 cycles of 94uC for 30 s, 55uC for 90 s, 72uC for 60 s, and a final cycle at 72uC for 10 min. One microliter of 10 times dilution PCR product was mixed with 8.5 mL of Hi-Di formamide and 0.125 mL of ROX size standard (Applied Biosystems). The mixture was denatured at 95uC for 5 min and then run on an ABI PRISM 3130xl DNA Analyzer (Applied Biosystems). Allele size for the highest stutter peak with the height ranging between 500 and 10,000 RFU were recorded. The genotyping was performed using GeneMapper 4.0 (Applied Biosystems).
Thirty-eight EST-SSR primer pairs amplified one product in diploid species and two discrete products in tetraploid species with height ranging RFU (Table S2) were selected for further analysis in the 268 Vigna accessions. Three primers with different labels and size products were mixed into a single PCR reaction mixture and amplified as multiplex PCR. After genotyping, a single allele size was scored for each marker in each accession corresponding to the strongest peak.

Data Analysis
Since both diploid (1 genome set) and tetraploid (2 genome set) species were analyzed together in this study and each SSR primer pair detected 2 loci in the tetraploid species, we separated the genome of V. reflexo-pilosa, into 2 different diploid types, viz. A and B. Therefore 10 genomes were recognized from the 8 taxa of 7 Vigna species.
Genetic distance (D A ) [24] for all possible pairs of accessions was calculated using software POPULATIONS 1.2.28 (available at www.cnrs-gif.fr/pge/bioinfo/populations). D A among the 10 genomes was also calculated using the same software. A phylogenic tree was constructed to reveal relationships among accessions based on D A using neighbor-joining clustering method by software MEGA 5.05 [25] with bootstrap support (1,000 replicates) obtained by re-sampling the allelic frequency data. In order to confirm the results from neighbor-joining clustering, principal coordinate analysis (PCoA) was also performed using software PAST [26]   the relationship of individuals based on PCoA analysis enabled 12 sub-genome groups to be determined. Using average D A genetic distance between 12 sub-genomes, an unrooted dendrogram showing relationships between these genomes was constructed by neighbor-joining method using software POPULATIONS 1.2.32. (available at www.cnrs-gif.fr/pge/bioinfo/populations) with bootstrap support (400 replicates) obtained by re-sampling the allelic frequency data. The trees were visualized with MEGA ver. 5.05 [25].

Azuki Bean EST-SSR Markers
A total of 11,167 cDNA clones were sequenced consisting of 7,534 clones from a seedling library and 3,633 clones from a young pod library. After clustering, 4,896 potential non-redundant EST sequences, including 2,350 contigs and 2,546 singletons, were generated with a total of 4,284,693 qualified bases (Table S3). By using the fuzznuc program in EMBOSS, a total of 1,188 SSRs were identified in the 4,896 non-redundant EST sequences. Of the 1,188 SSRs, di-, tri-, and tetra-nucleotide SSRs accounted for 71.6%, 26.7%, and 1.7%, respectively (Table S4). Assuming that total length of the non-redundant azuki bean EST sequences is 4.3 Mbp, the frequency of occurrence of the SSRs in transcribed regions was estimated to be one in every 3.6 kb. Altogether105 primer pairs were initially designed on the flanking regions of the 132 perfect SSR motifs (Table S4). To increase the number of candidate EST-SSR markers, additional 1,324 primer pairs were designed on the flanking regions of 1,545 imperfect SSR motifs allowing one or two base mismatching. As a result, a total of 1,429 EST-SSR markers were designed (Table S4), of which 149, 28, 5 and 5 markers identified 2, 3, 4 and #5 SSRs in the regions between the primer pairs, respectively. Thus the total number of identified SSRs by the 1,429 EST-SSR markers was 1,677 that consisted of 137 (8.2%) di-nucleotide repeats, 1,400 (83.5%) trinucleotide repeats and 140 (8.3%) tetra-nucleotide repeats (Table  S4). Among the di-nucleotide repeats, poly(AG) n (n = 88, 5.2% of total) were most frequently observed, followed by poly(AT) n (n = 35, 2.1%) and poly (AC) n (n = 14, 0.8%). Among the ten types of tri-nucleotide repeats observed, poly(AAG) n (n = 322, 19.2%) were the most abundant, followed by poly(GGA) n (n = 195, 11.6%) and poly(ATC) n (n = 180, 10.7%). Among the thirteen tetranucleotide repeats, poly(AAAG) n (n = 45, 2.7%) were the most frequently observed, followed by poly(AAAT) n (n = 30, 1.8%), and poly(AAAC) n (n = 20, 1.2%). The details of the designed azuki bean EST-SSR primers, along with the corresponding SSR motif, product size, and primer sequence, are available at http://marker. kazusa.or.jp/Azuki and in Table S5.

Amplification and Transferability of Azuki Bean EST-SSR Markers
Transferability and polymorphism of the 1,429 azuki bean EST-SSR markers were initially examined by polyacrylamide gel electrophoresis using eight accessions consisting of four major Vigna crop species; azuki bean complex, black gram, rice bean and mungbean (Table 1, Table S1). These wild and cultivated parental accessions had been used for a linkage map construction in each species ([27] [28] [29] [30]. Amplification of azuki bean EST-SSR markers in the same species, V. angularis, and the closely related species, V. umbellata, were 91.3 and 92.9% which were higher than those in V. radiata and V. mungo (Table 1). However, the size and pattern of PCR products of EST-SSR markers did not always reveal simple banding pattern with expected size, even in the same species. Among the amplified markers, 81.8 to 83.6% of them revealed simple banding pattern suitable for further application. When polymorphism between cultivated and wild accessions was examined, V. radiata possessed the highest polymorphism in which 277 markers (23.1% of amplified markers) were polymorphic. Further, 175 markers with good transferability in the four Vigna species were tested for polymorphism in a panel of 16 accessions from 8 taxa of Asian Vigna (Table S1). The test revealed that 22 markers (12.6%) failed to amplify, 36 (20.6%) amplified some accessions, whereas 117 (66.8%) successfully amplified all the 16 accessions (Table 2). Among the 117 amplifiable markers, 2 (1.1%) gave product size larger than 500 bp, 4 (2.3%) amplified multiple products, 45 (25.7%) amplified one product in tetraploid species, 10 (5.7%) were monomorphic, and 56 (32%) were polymorphic, amplifying one product in diploid species and two discrete products in tetraploid species and were considered to be suitable for analyzing origin of genome donor of the tetraploid species.

EST-SSR Polymorphism and Genetic Diversity
Thirty-eight primer pairs of azuki bean EST-SSR were analyzed in 268 accessions of 10 genomes from eight taxa of Asian Vigna. Based on discrimination of the diploid species we distinguished two discrete products detected by each marker of V. reflexo-pilosa into two types, viz. A and B. The product size that is the same or very similar to V. trinervia was scored as type A, while the other product size was scored as type B. In total, 445 alleles were detected in 10 genomes by the 38 EST-SSR loci (Tables 3,  4). The number of alleles detected per locus ranged between 3 (VES0777 and VES1271) and 34 (VES1172) with a mean of 11.7 alleles per locus ( Table 4) The allelic richness was between 2.0 for marker VES0777 and 9.9 for marker VES1172 with a mean of 4.7. Fifteen primers had allelic richness higher than 5.0 (Table 4). V. hirtella possessed the highest gene diversity (0.478), while V. reflexo-pilosa showed very low gene diversity (0.053 to 0.094) ( Table 3). All genomes of V. reflexo-pilosa showed no observed heterozygosity, while the other genomes showed very low observed heterozygosity from 0.008 to 0.028 (Table 3).

Genetic Relationship among Genomes
Genetic distance (D A ) among the ten genomes is shown in Table 5a. In most cases D A among genomes was high (.0.6). Nonetheless, D A between V. trinervia and V. reflexo-pilosa var. glabra (A) or V. reflexo-pilosa var. reflexo-pilosa (A) was low being 0.333 and 0.313, respectively. While, both V. reflexo-pilosa var. glabra (B) and V. reflexo-pilosa var. reflexo-pilosa (B) showed lowest D A with V. hirtella being 0.585 and 0.589 (Table 5a) or with V. hirtella (a1), being 0.429 and 0.422 (Table 6), respectively. Interestingly, V. trinervia, which belongs to section Angulares showed lower D A with V. radiata var. sublobata (section Ceratotropis) than with those species in the section Angulares.
The principle coordinate analysis (PCoA) was also conducted to confirm species relationship based on the neighbor-joining tree using the same genetic distance estimates. The first three PCs together accounted up to 56.6% of the total variation. The first, second and third PCs accounted for 38.1%, 10.6% and 7.9%, respectively ( Figure 4). A three-dimensional PC plot (PC1, PC2 and PC3) of the twelve genome groups is shown in Figure 4      Difference between the largest and smallest fragments amplified by each primer is shown in parentheses.
See Table 3 for the abbreviations of genome number. doi:10.1371/journal.pone.0104990.t004 Table 5. Genetic distance(D A ) within and among 10 genomes with (Vigna hirtella is not divided).

Transferability of Azuki Bean EST-SSRs
One of the advantages of SSRs derived from EST (EST-SSR markers) to other DNA marker systems is the high transferability across species and genera because of primer pairs are designed from conserved transcribed regions, although they are generally less polymorphic than genomic SSR markers (gSSR). In this study, more than 83% of the azuki bean EST-SSR markers were able to amplify DNA from 8 Asian Vigna species including four major crop species, viz. V. angularis, V. umbellata, V. radiata and V. mungo of the subgenus Ceratotropis. This result is in agreement with that of [31] who reported that more than 80% of mungbean EST-SSR markers were transferable to other species in the subgenus Ceratotropis. Transferability of the EST-SSR markers from azuki bean (in this study) and mungbean [31] is greater than that of gSSR markers from azuki bean (.67% transferable; [28]) and mungbean (.59% transferable; [32]). The high transferability of azuki bean EST-SSRs and of mungbean EST-SSRs to the related Vigna species indicates high genome conservation among the species in the subgenus Ceratotropis. Therefore, the azuki bean EST-SSR markers developed in this study will be useful for comparative genomic and genetic diversity studies in Vigna species. Nonetheless, as expected, polymorphism rate of the azuki bean EST-SSR markers detected between cultivated and wild accessions of azuki bean, rice bean, mungbean and black gram is low (14% in black gram and 23% in mungbean). In the same plant materials, Chaitieng et al. [28] reported that polymorphism rate of the genomic azuki bean SSRs was as high as 50% in black gram to 63% in mungbean. Thus low polymorphism rate of the azuki bean EST-SSRs can be an undesirable character for the use of these markers in comparative genome mapping in the genus Vigna. However, EST-SSRs represent true genetic diversity which may be more directly associated with traits of interest in breeding as compared to gSSRs [33], and they may reflect better relationships among the related species in genetic diversity study.
Evolution and Domestication of Tetraploid V. reflexopilosa V. reflexo-pilosa is the only tetraploid species of the genus Vigna. Although previous studies on Vigna species using DNA markers and plastid sequences clearly demonstrated that V. trinervia is a genome donor of V. reflexo-pilosa, those studies were unable to identify the other genome donor. In the present study, five putative genome donor species of V. reflexo-pilosa belong to the section Angulares and one species representing the section Ceratotropis were analyzed with azuki bean EST-SSRs. By assigning the discrete PCR products of V. reflexo-pilosa into two distinguished genomes (designated as genome A and B), we were able to identify the two genome donor parents of créole bean. Phylogenetic tree ( Figure 2) and PCoA plot (Figure 4) clearly demonstrated that genome A of V. reflexo-pilosa is received from V. trinervia. This result confirms the previous results reported by [7] [10] [13] [14] [15] that V. trinervia is a diploid genome donor. V. trinervia and V. reflexo-pilosa shares several similar morphological characters such as seed shape (rectangular), non-protruded hilum and large golden flower. The phylogenetic tree and the PCoA plot also unambiguously demonstrated that the other genome (genome B) of V. reflexo-pilosa is originated from V. hirtella. This result is in agreement with the result reported by Tateishi [7] [10] [13] [16] that V. hirtella is one of candidate genome donor of créole bean. Results from phylogenetic analysis of plastid DNA sequences of [14] [15] suggested that V. trinervia is the maternal genome donor of V. reflexo-pilosa. Since the Table 6. Genetic distance(D A ) within and among 10 genomes (Vigna hirtella is divided into three genomes).  Table 3 for the abbreviations of genome number.
Genetic distances within and among genomes are shown on the diagonal and below the diagonal, respectively.
The accessions that belong to genome a1, a2 and b are shown in Figure 3. previous and our results demonstrated close genetic relationship between tetraploid V. reflexo-pilosa and diploid V. trinervia/V. hirtella, we propose that créole bean evolved from interspecific hybridization between V. trinervia as female parent and V. hirtella as male parent, followed by genome duplication. Phylogenetic tree ( Figure 2) and PCoA plot ( Figure 4) showed lower divergence of genome A of V. reflexo-pilosa compared to genome B. Very low gene diversity within both genomes A and B suggested that this tetraploid species has a monophyletic (single) origin and only evolved recently. Gene diversity within the genome A in both V. reflexo-pilosa var. reflexo-pilosa (wild form) and V. reflexo-pilosa var. glabra (cultivated form) is very similar (0.058 vs. 0.053), while gene diversity of the genome B in wild créole bean is about twice of that in the cultivated one (0.94 vs. 0.52) ( Table 3). This indicates that the evolution rate of EST-SSRs on genome B of créole bean might be faster than that on genome A. The marked morphological differences between the cultivated and wild créole beans are thicker and erect glabrous stem in the former. Moreover, cultivated V. reflexo-pilosa still possesses a relatively high degree of pod shattering (Somta and Chankaew, personal observation), an important domestication trait for legume crops [34]. Therefore, the cultivated V. reflexo-pilosa can be treated as a semi-domesticated form. Based on cultivation, utilization and distribution of diploid genome donor species, créole bean appears to have been domesticated in Southeast Asia. The crop is now very rarely cultivated, although it is produced in northern mountainous villages in Vietnam under the same name and consumed in the same way as mungbean [2].
V. trinervia is widely distributed across Asia and also found in Papua New Guinea, Madagascar and East Africa. It is the second most widely distributed species in the subgenus Ceratotropis after V. radiata var. sublobata [2]. Nevertheless, gene diversity of V. trinervia was not high as compared to the other related species (Table 3). In Southeast Asia, especially in Thailand and the peninsular Malaysia, V. trinervia is often found on rural roadside habitats and in or near to rubber or oil palm plantations. It seems that some of the populations in those areas were introduced recently as a cover plant in the plantations [2]. This may account for the low diversity of V. trinervia.
V. tenuicaulis distributes in open wet habitats in northern Southeast Asia, and appeared to be closely related to V. hirtella (b) and V. angularis [13] [35] [6] [36]. Tomooka et al. suggested that there are no major barriers to hybridization between V. tenuicaulis and V. hirtella (b) [2]. Our results confirmed the close relationship between these two species (Figure 2, 4). Previous studies based on molecular, biochemical and morphological variations all showed high level of distinctness and variation  Table 6. doi:10.1371/journal.pone.0104990.g003 within V. tenuicaulis. EST-SSR marker variation in our study revealed contrasting results. Although we used as many as 30 accessions of V. tenuicaulis from three countries, viz. Laos, Myanmar and Thailand. (Figure 1, Table S1), the species showed relatively low gene diversity (0.309).
V. exilis has been reported only in Thailand and Myanmar [2] [14]. In Thailand, this species is restricted to rocky limestone mountains. Its habitats suggest that it may be useful as gene sources for resistant to alkaline soil and drought conditions. SSR analysis revealed high level of intra-specific diversity in V. exilis [37]. Our results also supported this despite only 13 accessions from narrow geographical origin (west of Thailand) of this species were used (Figure 1, Table S1), the species showed similar level of gene diversity to V. tenuicaulis and V. trinervia (Table 3).
V. minima has broad environments adaptation and grows well in shaded deciduous forest floors and open-wet habitats in East and Southeast Asia [2]. It is the only species in section Angulares that is found on the forest floor. Among the Asian Vigna species analyzed in this study, V. minima is the second most diverse species after V. hirtella (Table 3). The high differentiation of V. minima is due to wide geographical distribution of the analyzed accessions. In Southeast Asia, isolation of V. minima in forests in different mountainous regions across Thailand and Myanmar and its sporadic occurrence in patches of forests in those regions may account for high level of population divergence [35]. Accessions of V. minima in the present study were most closely related with V. hirtella, especially V. hirtella (b), followed by V. exilis (Figure 2, 4). Similar finding was reported by [35]. This agrees with their morphological appearance [2] that V. minima is sometimes confused with V. hirtella. V. minima can be distinguished from V. hirtella by smaller bracteole and more protruding hilum with well-developed rim-aril [2].
In summary, one of the advantages of EST-SSR markers to other DNA marker systems is their co-dominant nature and high transferability to closely related species. In this study, we used EST-SSR primer pairs, which amplified one PCR product in diploid species and two discrete PCR products in tetraploid species to separate genomes of V. reflexo-pilosa into 2 different diploid types, viz. A and B, that could successfully determine different sets of the donor genomes. Phylogenetic analyses revealed that V. trinervia and V. hirtella are the genome donors of V. reflexopilosa. Both genomes of cultivated and wild créole bean accessions showed low divergence suggesting that domestication of V. reflexopilosa is a relatively recent event.

Supporting Information
Table S1 Origin of 268 accessions used to analyze origin of genome of V. reflexo-pilosa. (XLSX)