Printer Friendly

Assessment of inter simple sequence repeat markers to differentiate sympatric wild and domesticated populations of common bean.

EFFICIENT MANAGEMENT of genetic resources, whether for conservation or for utilization in plant breeding programs, requires accurate and fast assessment of levels of genetic diversity and degrees of genetic relatedness. In the last decade, increasing use has been made of molecular markers based on DNA analyses (Bretting and Widrlechner, 1995; Gepts, 1995). An enduring concern with these markers is that they are fairly laborious and costly. Hence, there has been a trend to markers that, while highly polymorphic, are relatively easy and fast to analyze. In common bean, these include randomly amplified polymorphic DNA (RAPD) (e.g., Freyre et al., 1996) and amplified fragment length polymorphisms (AFLPs) (Tohme et al., 1996; Papa and Gepts, 2003). The ISSR marker technique involves polymerase chain reaction (PCR) amplification of DNA using a single primer composed of a microsatellite sequence, such as [(GACA).sub.4], anchored at the 3' or 5' end by two to four arbitrary, often degenerate, nucleotides' (Zietkiewicz et al., 1994). The ISSRs have proven to be a rapid, simple, and inexpensive way to assess genetic diversity (Kantety et al., 1995; Tsumura et al., 1996; Esselman et al., 1999; Gilbert et al., 1999; Sarla et al., 2003; Tanyolac, 2003; Jin et al., 2003), to identify closely related cultivars (Fang and Roose, 1997; Fang et al., 1997, 1998; Eiadthong et al., 1999; Prevost and Wilkinson, 1999; Martins et al., 2003), and to study evolutionary processes, such as reproductive systems (Liston et al., 2003) and gene flow (Wolfe et al., 1998). A recent paper by Galvan et al. (2003) describes differences between 10 Argentinean and three French bean cultivars, in particular broad differences between the Andean and Mesoamerican gene pools.

In the current research, the level of genetic diversity and differentiation was determined among one wild and four domesticated sympatric populations of common bean (Phaseolus vulgaris; 2n = 2x = 22) from the Sierra Norte de Puebla. Most studies of genetic diversity in common bean have been conducted in samples representing broad geographical ranges to identify major subdivisions in its germplasm, such as gene pools and races (Singh et al., 1991; Gepts, 1998). Fewer studies have been conducted to determine genetic relationships on smaller geographic levels (Cattan-Toupance et al., 1998; Papa and Gepts, 2003). These studies are useful because they provide information on evolutionary factors that can shape genetic diversity in common bean populations, such as gene flow. For example, Cattan-Toupance et al. (1998) showed the importance of a disease pressure in determining local population differentiation, in addition to genetic drift. They also considered partial outbreeding to be a potential factor affecting association of traits. Papa and Gepts (2003) showed that, in sympatry, gene flow from domesticated to wild types is three times more important than in the opposite direction. Furthermore, spatial autocorrelation studies showed that wild beans exhibited a stronger local population structure than domesticated beans.

In addition to studying the level of genetic diversity among sympatric populations of common bean, we sought to determine the usefulness of ISSR markers in common bean with respect to their reproducibility, variation in DNA amplification among primers, relative linkage map location, and levels of polymorphism detected. The ISSRs were able to distinguish sympatric bean populations but also identify within-population polymorphism.


Plant Materials

Seeds of a wild population of P. vulgaris were collected in early winter of 1995 by Francisco Basurto (FB1934) in the municipalities of Xochitlan de Romero Rubio and Huapalejcan in the Sierra Notre de Puebla (state of Puebla, Mexico). This wild population grows within bean fields, but surprisingly not outside fields. Seeds of four commonly grown landraces--labeled according to their seed color: Beige (B), Negro Brillante (NB, shiny black), Negro Opaco (NO, dull black), and Rojo (R, red)--were collected from local, indigenous farmers from the same area. As far as we can tell, these are not bred materials. All landraces were indeterminate strong climbers, classified as Type IV (Singh, 1982).

DNA Extraction and ISSR Analysis

A total of 615 seeds (Table 1) were planted in a greenhouse at Davis. These seeds included 79, 156, 138, 34, and 208 individuals for the Beige, Negro Brillante, Negro Opaco, Rojo, and wild populations, respectively. DNA was extracted using the CTAB method as described by Doyle and Doyle (1987). Inter simple sequence repeats were analyzed using a total of 32 primers synthesized by the Protein Sequencing Laboratory of the University of California, Davis. From the 30 primers tested originally, four primers were selected that produced a high level of polymorphism among both landraces and wild populations (Table 1). Two of these primers, [(GACA).sub.3] RT and [(GACAC).sub.2], were run on the entire sample, whereas two additional primers, YR [(GACA).sub.3] and [(ACTG).sub.3] RG, were only run on a randomly chosen subset of 417 individuals, which included 61, 100, 117, 32, and 107 individuals of the Beige, Negro Brillante, Negro Opaco, Rojo, and wild populations, respectively. Each 20-[micro]L amplification reaction consisted of 10 mM Tris-HCl (pH 9.0), 50 mM KCl, 0.1% Triton X-100, 2 mM Mg[Cl.sub.2], 200 [micro]M of each deoxyribonucleotide phosphate, 1 [micro]M of primer, 1 unit of Taq polymerase (Promega, Madison, WI), and 20 ng of template DNA. Each reaction mixture was overlaid with mineral oil. Amplification was performed in a 96-well Ericomp thermal cycler (Ericomp, Inc., San Diego, CA) under the following conditions: 4 min at 94[degrees]C for 1 cycle, followed by 2 min at 94[degrees]C, 1 min at 42[degrees]C, and 2 min at 70[degrees]C for 35 cycles, and 5 min at 72[degrees]C for final extension. Amplification products were separated on 320 by 380 by 0.4 mm of 6% nondenaturing 30:1 bis-acrylamide gels containing 3 M urea and 1 x tris-borate buffer (Zietkiewicz et al., 1994). A 123-bp molecular marker standard was included in each gel. Gels were run for 12 h at 300 V, and PCR products were detected by silver staining (Bassam et al., 1991) (Fig. 1). Each gel was loaded with PCR products of at least eight individuals of each accession to be able to compare the amplification pattern among accessions within the same gel. The PCR fragments generated by the ISSR primers were labeled with a code consisting of a description of the primer used for the reaction followed by the molecular size of the fragment in base pairs [i.e., [(ACTG).sub.3] RG.434: fragment size of 434 bp].


Linkage Mapping

The same four ISSRs primers used to generate fingerprinting in the five populations (Table 1) were used to generate PCR-amplified DNA fragments in 75 lines of the BAT93 x Jalo EEP556 recombinant inbred population, which had reached the [F.sub.8] generation. This population has been used previously to develop a RFLP genetic linkage map (Nodari et al., 1993) and to develop a core linkage map for the common bean genome (Freyre et al., 1998). Segregation was scored for presence or absence of the ISSRs bands. Segregation of polymorphic fragments was analyzed by a chi-square test for goodness-of-fit to a 1:1 ratio. Linkage analysis of the polymorphic set of markers was done using Mapmaker 3.0b (Lander et al., 1987) as described in Menendez et al. (1997) (Fig. 2). Briefly, markers already previously established in the common bean map (Freyre et al., 1998) were used to anchor linkage groups. The Assign command was then used to allocate ISSRs markers into linkage groups. To establish the most likely order of markers within each linkage group, the command Order was applied under the following conditions: (i) LOD score above or equal to 3.00, and (ii) a maximum distance of 30 cM. Markers not placed under this limit were placed using lower LOD scores (between 3 and 2).


Analyses of Genetic Diversity

The PCR products amplified using the ISSR primers were scored as present (1) or absent (0) for the five populations. Analyses were initially conducted on individuals with results from four or two primers. Because results were similar between the two groups, only the results based on individuals with complete data sets (four primers) are presented here. Of a total of 20 850 data points (417 individuals x 50 markers), 1079 (5%) were considered as missing data. Individuals with such missing data were excluded without effect on the results (data not shown). Thus, our final sample contained 329 individuals, including 51, 85, 92, 24, and 77 individuals from the Beige, Negro Brillante, Negro Opaco, Rojo, and wild populations, respectively. When considering all 50 fragments together, 170 multilocus combinations (haplotypes) could be identified. Unless otherwise mentioned, standard indices of genetic diversity were calculated using the Arlequin (Schneider et al., 2000) or the POPGENE (Yeh et al., 1997) softwares. Indices of genetic similarity were calculated using Dice's coefficient (Dice, 1945) to estimate the relationship between accessions and to observe the placement of putative outcrossed individuals. From the similarity matrix, a sequential, agglomerative, hierarchical, and nested (SAHN) cluster analysis was performed using the unweighted pair group method with arithmetic means (UPGMA) algorithm computed using NTSYS-pc (Rohlf, 1997). A dendrogram was generated using NTSYS-pc (Rohlf, 1997) to show the genetic relationships and distances of each accession (Fig. 3).



Development and Mapping of ISSR Markers in Common Bean

The pattern of DNA amplification using ISSR primers was very reproducible based on results from at least 20 different DNA extractions and PCR reactions from each landrace, which produced reproducible banding patterns (Fig. 1). All primers tested amplified some bands in common bean except primers [(AT).sub.8]YG and [(AT).sub.8]RG, with R = G, A and Y = C, T (data not shown). The four ISSR primers selected amplified 50 bands, which were strongly amplified and well-stained across gels. The bands ranged in size from 250 to 1300 bp. To examine the distribution of these ISSR markers in the bean genome, a mapping study was conducted in the BAT93 x Jalo EEP558 recombinant inbred population, the core mapping population in common bean (Freyre et al., 1998). The core map established in this population includes some 560 markers, including 120 RFLPs and 430 RAPDs. The four ISSR primers used in the gene flow studies amplified 54 polymorphic fragments in this population. The segregation data of these fragments were added to those of existing framework markers to determine their location on the map. The Chi-square test (1:1 expectation, P = 0.99) indicated that nine markers (16.6%) showed a segregation distortion that was significantly different from the expected 1:1 ratio at the 1% level. Fifty of the amplified bands mapped onto nine of the 11 linkage groups of the common bean linkage map, indicating a broad coverage of the genome by the ISSRs (Fig. 2). None of the ISSRs markers scored mapped on linkage groups 4 or 9 (Fig. 2) and four bands were unlinked. The mapped ISSRs markers were distributed mainly on linkage group 7 (11 markers), 10 (9 markers), and 11 (9 markers). Further inspection of Fig. 2 suggests some degree of clustering of these markers on individual linkage groups.

Polymorphisms in Phaseolus vulgaris of the Sierra Norte de Puebla

Ninety-two (46/50) percent of the analyzed markers displayed variation in this study (Tables 1 and 2). Furthermore, the five populations analyzed in this study showed contrasting levels of polymorphism. The wild population was the most diverse, whether measured by the number of polymorphic fragments, Nei's (1973) gene diversity, or the number of haplotypes (Table 1). As a whole, the domesticated group (n = 252) had higher or similar levels of polymorphism compared with the sample of the wild population (n = 77) (Table 1). None of the regressions between population size (independent variable) and measures of genetic diversity as dependent variables (number of polymorphic bands, Nei's gene diversity, and number of haplotypes) were statistically significant at the P = 0.05. Thus, differences in genetic diversity among populations were not due to differences in population sizes.

Genetic Distances and Differentiation

The ISSR markers clearly distinguished the wild from the domesticated populations included in this study (Fig. 3). Among the 50 markers used for the final analysis, 22% showed highly contrasting frequencies in the wild and domesticated gene pools: A12, A13, A40, B23, B26, B27, C6, C18, D9, D28, and D33. Frequencies of these markers in one gene pool were below 20% and above 80% in the other gene pool (Table 2). Fourteen percent were monomorphic (frequencies [greater than or equal to] 0.95) among these two groups (i.e., markers A24 [[(GACA).sub.3]RG.703], A29 [[(GACA).sub.3]RG.616], B17 [YR[(GACA).sub.3]], C8 [[(GACAC).sub.2], 887], C10 [[(GACAC).sub.2].821], C12 [[(GACAC).sub.2].727], and D5 [[(ACTG).sub.3]RG.1013] in Table 2). The remaining percentage consisted of bands with various frequency patterns in the wild and domesticated components (Table 2). In addition, ISSR markers were also able to clearly distinguish the four domesticated populations, including the Negro Opaco and Negro Brillante populations, in spite of the phenotypic similarities of the latter (Fig. 3). Although the main difference between the two black-seeded populations--the shininess of the seed--is controlled by a single gene (Asp; Beninger et al., 2000), there were four markers that were able to differentiate these two populations: D12 [[(ACTG).sub.3]RG.652], D25 [[(ACTG).sub.3]RG.485], D27 [[(ACTG).sub.3]RG.434], and D29 [[(ACTG).sub.3]RG.415]. A more heterogeneous fingerprinting pattern was observed in the wild accessions, but nevertheless, some bands were found that had a high frequency in one group of accessions (either wild or domesticated) and were present at a much lower frequency in the other group (Table 2).

All ISSRs markers were present at some frequency in the wild and the two black-seeded landraces (Table 2). The beige and red landraces did not have unique markers associated with them, but were missing some of the bands observed in the other three accessions. Within each population, markers were generally almost completely fixed. Their frequency was either close to one or zero, especially in the domesticated populations, reflecting the codominant nature of ISSRs (Fig. 4). Among domesticated populations, most of the variation was distributed among and not within populations as shown by a [G.sub.ST] value (Nei, 1973) of 0.70 ([H.sub.T] = 0.23; [H.sub.S] = 0.07). Adding the wild population increased total genetic diversity but did not change the distribution of genetic diversity markedly ([G.sub.ST] = 0.67; [H.sub.T] = 0.30; [H.sub.S] = 0.10). These results are consistent with the predominantly selling nature of the species, as further confirmed by the low values for Nm, the average number of migrants per generation (around 0.24). The Nm values of 1 are generally considered to be the threshold above which migration will counteract population differentiation (Neigel, 1997). The wild population, on one hand, and all the domesticated populations, considered as a single group, showed low differentiation ([F.sub.ST] = 0.06) and, concurrently, a high Nm value of 8.47. In contrast, among the different domesticated populations, a much higher differentiation was observed, ranging from [F.sub.ST] = 0.49 between the two black-seeded populations to [F.sub.ST] [approximately equal to] 0.85 among the Beige, Negro Opaco, and Rojo populations (Table 3).


A UPGMA dendrogram resulting from a SAHN cluster analysis based on the Dice similarity index is depicted in Fig. 3. It includes 329 individuals that were evaluated with the four ISSR primers and had a complete data set. The phenogram contains six clusters, while only five populations were included in this study. All the beige individuals (n = 51) grouped in a single cluster, which also included two red-seeded individuals (R23 and R24). With the exception of the latter two individuals, all the red-seeded individuals grouped in a single cluster, which did not contain any individuals from the other populations, Most of the individuals of the wild population grouped in a single cluster that included also two individuals from the Negro Brillante population (NB84 and NB85). The Negro Brillante cluster contained two Negro Opaco individuals (NO51 and NO56) and the Negro Opaco cluster included one Negro Brillante individual (NB29). The mixed cluster contained seven Negro Brillante (NB24, 49, 63, 64, 73, 81, 83), nine Negro Opaco (NO16, 58, 64, 66, 67, 82, 83, 85, 92), and two wild individuals (W62, 66). Thus, there were 25/339 (8%) individuals for which the seed phenotypes did not match the cluster membership based on a UPGMA analysis of molecular data.


Inter simple sequence repeats have been used to explore genome structure, to assess genetic diversity within germplasm collections as well as to determine the frequency of simple sequence repeats in several species (Blair et al., 1999; Gilbert et al., 1999; Prevost and Wilkinson, 1999). The ISSRs have proven to be a reliable, easy to generate, and versatile set of markers that do not require previous knowledge of the genome sequence to generate DNA markers, unlike SSRs (Zietkiewicz et al., 1994; Gupta et al., 1994). The flexibility to design primers containing a di-, tri-, or tetra-nucleotide repetitive motif anchored by one or more nucleotides at the 3' or 5' end make them ideal to explore the genome of any species, including those without previous knowledge of DNA sequence. Several reports have compared the level of polymorphism detected using RFLPs, RAPDs, and ISSRs (Nagaoka and Ogihara, 1997), and more recently ISSRs and AFLPs. The consensus is that ISSRs are very powerful to detect polymorphism, scan the whole genome, are inexpensive, and easy to generate. They can detect even more polymorphism than RFLPs (Kantety et al., 1995) and more than AFLPs in rice (Blair et al., 1999). However, it is now becoming clear that the polymorphism detected by ISSRs is dependent on the plant species being investigated as well as the type of simple sequence repeat incorporated in the ISSRs primer used to amplify PCR products. Our results confirm these earlier observations. Several characteristics make ISSRs also a very useful marker in common bean. They are able to distinguish among closely related common bean populations. The ISSR primers amplify fragments that are dispersed throughout the genome as demonstrated by mapping ISSR markers in the [F.sub.8] generation of a recombinant inbred population. The ISSR markers were able to detect individual resulting from potential introgression among wild and domesticated populations (although independent evidence is necessary for confirmation, see below).

Microsatellite fingerprinting in the genus Phaseolus was investigated by Hamann et al. (1995), who concluded that the repetitive element [(GATA).sub.4] and, to a lesser extent, [(GACA).sub.4] were well represented in P. vulgaris and P. lunatus. They also concluded that that the dinucleotide motif [(CA).sub.8] seemed to be less represented in the of P. vulgaris. On this basis, several of our primers contained [(GACA).sub.3] and [(GATA).sub.3] as the repeat sequence. Among the five accessions tested in this study, we found that the primers containing [(GACA).sub.3] as the repetitive motif produced better amplification and higher polymorphism than those containing [(GATA).sub.3]. Although a nonexhaustive effort was attempted to optimize and characterize all the primers designed for amplifying ISSRs, all the primers tested amplified fragments in beans, indicating the potential of this technique to measure genetic variability in bean germplasm collections. The combination of the four ISSR primers was sufficient to distinguish all the populations included in this study. Generation of ISSR bands was based on four primers, two of which included [(GACA).sub.4] as the repeat sequence anchored either at the 3' or 5' end (Table 1). The double pentanucleotide motif [(GACAC).sub.2] was not anchored but produced a very reliable and useful fingerprinting that differentiated all the populations included in this study. Coverage of the whole genome by ISSRs has been claimed previously by others (Blair et al., 1999; Fang and Roose, 1999). Nevertheless, ISSRs were difficult to map in einkorn wheat (Kojima et al., 1998). In our mapping, we found that a sample of some 50 markers mapped to 9 of the 11 linkage groups of the existing core genetic map of bean (Freyre et al., 1998). Although a tendency for clustering was observed among some of the markers (Fig. 2), this has been also observed in some attempts to map AFLPs markers in beans (R. Papa and P. Gepts, 2004, unpublished results).

There was a strong association between molecular haplotypes and seed color as 304/329 (92%) of the individuals clustered with individuals showing the same seed color (Fig. 3). Conversely, 25/329 (8%) of the individuals grouped with individuals with seeds of a different color. A strong association is expected in predominantly self-pollinated species (Hedrick et al., 1978). This reproductive system leads to multilocus associations, whether the loci involved are linked or not. Our results show that ISSR markers in common bean are distributed over the entire bean genome, although individual markers can be linked (Fig. 2). Seed color and color pattern genes are also distributed throughout the bean genome (McClean et al., 2002). Thus, in the absence of outcrossing and the ensuing recombination, multilocus associations between seed color and ISSR markers can be expected.

Given this expectation, what is the origin of the individuals for which the association broke down? We suggest that there are two possible reasons. First, these individuals could represent incomplete lineage sorting, the latter being defined as the progressive loss of polymorphisms in related lineages of populations eventually leading to the absence of shared polymorphisms. In the transitional period between the separation of the populations from their common ancestor and complete lineage sorting, certain polymorphisms will remain shared unless they arose after the separation such as the seed color alleles. The landraces included in this study are probably not derived directly from the sympatric wild population as the Mesoamerican bean domestication has been traced genetically to the west-central part of Mexico (Jalisco and Guanajuato; Gepts, 1988). The archaeological record of common bean in Mesoamerica is currently not older than 2300 yr (Kaplan and Lynch, 1999). Thus, on an evolutionary time scale, the separation between wild and domesticated beans and among the different landraces may have been too short to achieve full lineage sorting.

An alternative explanation is that the unusual individuals actually result from hybridization one or more generations before the observations reported in this study. On the basis of the dendrogram, the extent of introgression would reach about 8%. However, dendrograms may underestimate the frequency of hybrids because separate clusters only appear when a sufficient number of markers show significant differences in frequency. For example, the two most closely related populations (Negro Opaco and Negro Brillante; Table 3) could be distinguished by six markers (B25, B27, D12, D26, D27, and D29) with a frequency difference of at least 0.48 (Table 2). Introgression, however, may affect only a minor portion of the genome, which would only be detected with a small number ([greater than or equal to] 1) of markers. These would not necessarily be detected in a dendrogram. A different way of analyzing this issue is to consider individual markers, especially those occurring at low frequency (<0.10 or <0.05) in the different populations, but occur at a higher frequency in the other populations analyzed (for examples, see Fig. 4). Individuals exhibiting one or more of these markers could have resulted from introgression (although incomplete lineage sorting cannot be excluded, as pointed out earlier). At the more stringent level (5%), some 20 to 25% of the individuals carried one or more loci with low-frequency markers. For the less stringent criterion of 10%, the frequency ranged between 23% for the Beige population to 58% for the wild population (data not shown). Overall, these estimates of potentially outcrossed or introgressed individuals based on individual loci were therefore higher that the estimates based on the entire set of loci, obtained either statistically or graphically.

Common bean is considered a self-pollinated species, but several lines of evidence suggest that outcrossing may be occasionally important (Wells et al., 1988). Ibarra-Perez et al. (1997) recently summarized the literature reporting on experiments conducted to directly measure the level of outcrossing. Although the outcrossing rates were generally below 5%, some studies reported higher levels, reaching 60 to 80% in some instances (Table 1 in Ibarra-Perez et al., 1997). Differences in outcrossing rates appeared to vary according to environmental conditions prevailing during the study, the bean genotypes involved, and the presence of pollinating insects. Debouck et al. (1993) and Freyre et al. (1996) identified several cases of gene flow from domesticated to wild beans in Colombia, Ecuador, and Bolivia. Beebe et al. (1997) mentioned that extensive wild-weedy-domesticated complexes of P. vulgaris were observed in regions of Peru and Colombia where wild and domesticated beans are sympatric. In the Mesoamerican gene pool, Vanderborght (1983) described the existence of weedy forms in Mexico and introgression rates of up to 50% in wild populations. Crosses between wild and domesticated common beans yield viable and fertile progenies (Burkart and Briicher, 1953; Miranda Colin, 1979; Evans, 1980; Koenig and Gepts, 1989; Kornegay et al., 1993; Singh et al., 1995; Pereira et al., 1996; Mumba and Galwey, 1998, 1999). With the exception of a preliminary report by Triana et al. (1993), all estimates of outcrossing in beans involved cultivars. Although the Triana et al. (1993) study involved wild as well as domesticated beans, their study was conducted on a research station outside the natural distribution area of wild beans. Thus, further investigations into levels of outcrossing in common bean in different locations and years is of great interest. Additional research is needed to determine whether the individuals in which the correlation between seed color and ISSR markers has apparently been broken or which carry unusual markers result from incomplete lineage sorting or hybridization, or a combination of both.
Table 1. Levels of polymorphism in wild and domesticated populations
of common bean from the Sierra Norte de Puebla, Mexico.

                                       Number of polymorphic
                                          ISSR ([dagger])

                           No. of         [(GACA).sub.3]RT
Population               individuals   (12) ([double dagger])

Beige                         51                  6
Negro Brillante               85                  8
Negro Opaco                   92                  8
Rojo                          24                  9
Domesticated (#)             252                 12
Wild                          77                 10
All ([dagger][dagger])       329                 12

                             Number of polymorphic
                           ISSR ([dagger]) fragments

                                B                  C
                         YR[(GACA).sub.3]   [(GACAC).sub.2]
Population                     (12)               (9)

Beige                            4                 1
Negro Brillante                  7                 6
Negro Opaco                      8                 7
Rojo                             5                 5
Domesticated (#)                 9                 7
Wild                             8                 6
All ([dagger][dagger])          10                 8

                         Number of polymorphic
                            ISSR ([dagger])

                         [(ACTG).sub.3]RG   Total
Population                     (17)         (50)

Beige                            4           15
Negro Brillante                 15           36
Negro Opaco                     11           34
Rojo                            10           29
Domesticated (#)                15           43
Wild                            16           40
All ([dagger][dagger])          16           46

                         Nei's gene        No. of
Population               diversity       haplotypes

Beige                       0.02       14 ([section])
Negro Brillante             0.10       46 ([paragraph])
Negro Opaco                 0.06       28 ([paragraph])
Rojo                        0.10       10 ([section])
Domesticated (#)            0.19       95
Wild                        0.20       75
All ([dagger][dagger])      0.29      170

([dagger]) ISSR, inter simple sequence repeat.

([double dagger]) ISSR primer (see Materials and Methods) and
total number of fragments in parenthesis.

([section]) 1 haplotype shared between Beige and Rojo.

([paragraph]) 2 haplotypes shared between Negro Brillante and
Negro Opaco.

(#) Beige + Negro Brillante + Negro Opaco + Rojo - 3 shared

([dagger][dagger]) Domesticated + wild.

Table 2. Frequency of inter simple sequence repeat markers
within the different populations. ([dagger])

Fragment no.                    Marker

A11            [(GACA).sub.3] RG.1000 ([double dagger])
A12            [(GACA).sub.3] RG.957
A13            [(GACA).sub.3] RG.936
A16            [(GACA).sub.3] RG.896
A18            [(GACA).sub.3] RG.820
A19            [(GACA).sub.3] RG.802
A21            [(GACA).sub.3] RG.785
A24            [(GACA).sub.3] RG.703
A26            [(GACA).sub.3] RG.658
A29            [(GACA).sub.3] RG.616
A37            [(GACA).sub.3] RG.500
A40            [(GACA).sub.3] RG.429
A43            [(GACA).sub.3] RG.359
A44            [(GACA).sub.3] RG.355
B15            YR [(GACA).sub.3].873
B16            YR [(GACA).sub.3].856
B17            YR [(GACA).sub.3].759
B19            YR [(GACA).sub.3].715
B20            YR [(GACA).sub.3].680
B23            YR [(GACA).sub.3].634
B25            YR [(GACA).sub.3].591
B26            YR [(GACA).sub.3].568
B27            YR [(GACA).sub.3].551
B31            YR [(GACA).sub.3].438
B33            YR [(GACA).sub.3].253
C6             [(GACAC).sub.2].963
C7             [(GACAC).sub.2].947
C8             [(GACAC).sub.2].887
C10            [(GACAC).sub.2].821
C11            [(GACAC).sub.2].743
C12            [(GACAC).sub.2].727
C14            [(GACAC).sub.2].623
C18            [(GACAC).sub.2].479
C20            [(GACAC).sub.2].424
D1             [(ACTG).sub.3] RG.1208
D2             [(ACTG).sub.3] RG.1181
D5             [(ACTG).sub.3] RG.1013
D9             [(ACTG).sub.3] RG.795
D12            [(ACTG).sub.3] RG.652
D14            [(ACTG).sub.3] RG.638
D16            [(ACTG).sub.3] RG.635
D18            [(ACTG).sub.3] RG.621
D19            [(ACTG).sub.3] RG.604
D25            [(ACTG).sub.3] RG.485
D26            [(ACTG).sub.3] RG.479
D27            [(ACTG).sub.3] RG.434
D28            [(ACTG).sub.3] RG.420
D29            [(ACTG).sub.3] RG.415
D32            [(ACTG).sub.3] RG.345
D33            [(ACTG).sub.3] RG.333

Fragment no.    B      NB     NO     R      W

A11            0.98   1.00   1.00   0.18   0.95
A12            0.98   0.96   0.98   0.00   0.13
A13            0.02   0.05   0.06   1.00   0.98
A16            1.00   0.92   0.96   1.00   0.96
A18            1.00   1.00   1.00   0.00   1.00
A19            0.03   1.00   0.96   1.00   0.15
A21            1.00   0.08   0.12   0.86   0.93
A24            1.00   1.00   1.00   1.00   1.00
A26            1.00   1.00   1.00   0.14   1.00
A29            1.00   1.00   1.00   1.00   1.00
A37            1.00   0.92   1.00   1.00   0.28
A40            0.03   0.11   0.06   0.11   0.90
A43            0.02   0.21   0.33   0.18   0.73
A44            0.00   0.05   0.09   0.04   0.34
B15            0.88   0.09   0.04   0.27   0.49
B16            1.00   1.00   1.00   0.88   0.99
B17            1.00   1.00   1.00   1.00   0.99
B19            0.00   0.06   0.10   0.00   0.62
B20            0.00   0.91   0.97   0.04   0.07
B23            0.00   0.03   0.03   0.04   0.79
B25            0.98   0.68   0.05   1.00   0.04
B26            1.00   1.00   1.00   1.00   0.18
B27            0.00   0.52   0.04   0.04   0.99
B31            0.98   1.00   0.99   0.12   0.38
B33            0.05   0.97   0.99   0.96   0.98
C6             1.00   1.00   0.95   0.00   0.09
C7             0.00   0.02   0.06   1.00   0.48
C8             1.00   1.00   1.00   1.00   0.95
C10            1.00   0.99   0.98   1.00   1.00
C11            0.00   0.08   0.01   1.00   0.04
C12            1.00   1.00   1.00   1.00   1.00
C14            0.00   0.03   0.02   1.00   0.12
C18            0.02   0.10   0.04   0.00   0.99
C20            0.00   0.98   0.97   0.00   0.30
D1             1.00   0.02   0.05   0.00   0.87
D2             0.02   0.99   0.99   1.00   0.32
D5             1.00   1.00   1.00   1.00   0.96
D9             0.15   0.08   0.10   0.00   0.92
D12            1.00   0.97   0.09   0.00   0.36
D14            0.00   0.04   0.03   0.10   0.36
D16            0.97   0.98   1.00   0.80   0.94
D18            1.00   0.99   1.00   1.00   0.60
D19            1.00   0.99   1.00   1.00   0.79
D25            0.00   0.27   0.95   1.00   0.04
D26            1.00   0.89   0.11   0.03   0.98
D27            1.00   0.90   0.05   0.00   0.03
D28            0.02   0.07   0.04   0.00   0.85
D29            0.03   0.16   0.89   1.00   0.10
D32            1.00   0.99   1.00   0.97   0.47
D33            0.00   0.04   0.04   0.07   0.85

([dagger]) B = Beige, NB = Negro Brillante, NO = Negro Opaco,
R = Rojo, W = wild.

([double dagger]) The fragment name consists of the primer sequence
followed by the amplified fragment size expressed in bp.

Table 3. Genetic differentiation (FST) and number of migrants
(Nm) among populations included in this study. ([dagger])

Populations       Beige   Negro Brillante   Negro Opaco   Rojo   Wild

Beige                          0.24            0.10       0.08   0.22
Negro Brillante   0.68                         0.52       0.19   0.25
Negro Opaco       0.83         0.49                       0.17   0.19
Rojo              0.86         0.72            0.75              0.27
Wild              0.69         0.67            0.73       0.65

([dagger]) Above the diagonal: Nm; below the diagonal: [F.sub.ST].

Abbreviations: AFLP, amplified fragment length polymorphism: ISSR, inter simple sequence repeat; PCR, polymerase chain reaction: RAPD, randomly amplified polymorphic DNA: RFLP, restriction fragment length polymorphism; SAHN, sequential, agglomerative, hierarchical, and nested.


We are grateful to the cooperation and assistance of Don Vicente (Nauzontla, Puebla), farmer-owner of the field where the experiment took place. Field work would not have been possible without the assistance of Gabriel Flores Franco, Sara Fuentes Soriano, Pedro Mercado Ruaro, Francisco Basurto, and Leticia Torres Colin. This project was funded by the McKnight Foundation, Minneapolis, MN, USA.


Bassam, B., G. Caetano-Anolles, and P. Gresshoff. 1991. Fast and sensitive silver staining of DNA in polyacrylamide gels. Anal. Biochem. 196:80-83.

Beebe, S., O. Toro, A. Gonzalez, M. Chacon, and D. Debouck. 1997. Wild-weed-crop complexes of common bean (Phaseolus vulgaris L., Fabaceae) in the Andes of Peru and Colombia, and their implications for conservation and breeding. Genet. Res. Crop Evol. 44: 73-91.

Beninger, C., G. Hosfield, M. Bassett, and S. Owens. 2000. Chemical and morphological expression of the B and Asp seedcoat genes in Phaseolus vulgaris. J. Am. Soc. Hortic. Sci. 125:52-58.

Blair, M., O. Panaud, and S. McCouch. 1999. Inter-simple sequence repeats (ISSR) amplification for analysis of microsatellite motif frequency and fingerprinting of rice (Oryza sativa L.). Theor. Appl. Genet. 98:780-792.

Bretting, P.K., and M.P. Widrlechner. 1995. Genetic markers and plant genetic resource management. Plant Breed. Rev. 13:11-86.

Burkart, A., and H. Brucher. 1953. Phaseolus aborigineus Burkart, die mutmassliche andine Stammform der Kulturbohne. Zuchter 23:65-72.

Cattan-Toupance, I., Y. Michalakis, and C. Neema. 1998. Genetic structure of wild bean populations in their South-Andean centre of origin. Theor. Appl. Genet. 96:844-851.

Debouck, D.G., O. Toro, O.M. Paredes, W.C. Johnson, and P. Gepts. 1993. Genetic diversity and ecological distribution of Phaseolus vulgaris in northwestern South America. Econ. Bot. 47:408-423.

Dice, L.R. 1945. Measures of the amount of ecologic association between species. Ecology 26:297-302.

Doyle, J., and J. Doyle. 1987. A rapid DNA isolation procedure fro small quantities of fresh leaf tissue. Phytochem. Bull. 19:11-15.

Eiadthong, W., K. Yonemori, A. Sugiura, N. Utsunomiya, and S. Subhadrabandhu. 1999. Identification of mango cultivars of Thailand and evaluation of their genetic variation using the amplified fragments by simple sequence repeat-(SSR-) anchored primers. Sci. Hortic. 82:57-66.

Esselman, E., L. Jianqiang, D. Crawford, J. Windus, and A. Wolfe. 1999. Clonal diversity in the rare Calamagrostis porteri ssp. insperata (Poaceae): Comparative results for allozymes and random amplified polymorphic DNA (RAPD) and intersimple sequence repeat (ISSR) markers. Mol. Ecol. 8:443-451.

Evans, A.M. 1980. Structure, variation, evolution, and classification in Phaseolus. p. 337-347. In R.J. Summerfield and A.H. Bunting (ed.) Advances in legume science. Royal Botanic Gardens, Kew, UK.

Fang, D., and M. Roose. 1997. Identification of closely related citrus cultivars with inter-simple sequence repeat markers. Theor. Appl. Genet. 95:408-417.

Fang, D., and M. Roose. 1999. Inheritance of inter simple sequence repeat markers in citrus. J. Hered. 90:247-249.

Fang, D., R. Krueger, and M. Roose. 1998. Phylogenetic relationships among selected Citrus germplasm accessions revealed by intersimple sequence repeat (ISSR) markers. J. Am. Soc. Hortic. Sci. 123:612-617.

Fang, D., M. Roose, R. Krueger, and C. Federici. 1997. Fingerprinting trifoliate orange germplasm accessions with isozymes, RFLPs, and inter-simple sequence repeat markers. Theor. Appl. Genet. 95: 211-219.

Freyre, R., R. Rios, L. Guzman, D. Debouck, and P. Gepts. 1996. Ecogeographic distribution of Phaseolus spp. (Fabaceae) in Bolivia. Econ. Bot. 50:195-215.

Freyre, R., P. Skroch, V. Geffroy, A.-F. Adam-Blondon, A. Shirmohamadali, W. Johnson, V. Llaca, R. Nodari, P. Pereira, S.-M. Tsai, J. Tohme, M. Dron, J. Nienhuis, C. Vallejos, and P. Gepts. 1998. Towards an integrated linkage map of common bean. 4. Development of a core map and alignment of RFLP maps. Theor. Appl. Genet. 97:847-856.

Galvan, M.Z., B. Bornet, P.A. Balatti, and M. Branchard. 2003. Inter simple sequence repeat (ISSR) markers as a tool for the assessment of both genetic diversity and gene pool origin in common bean (Phaseolus vulgaris L.). Euphytica 132:297-301.

Gepts, P. (ed.) 1988. Phaseolin as an evolutionary marker, p. 215-241. In Genetic resources of Phaseolus beans. Kluwer, Dordrecht, the Netherlands.

Gepts, P. 1995. Genetic markers and core collections, p. 127-146. In T. Hodgkin et al. (ed.) Core collections of plant genetic resources. John Wiley & Sons, New York.

Gepts, P. 1998. Origin and evolution of common bean: Past events and recent trends. HortScience 33:1124-1130.

Gilbert, J., R. Lewis, M. Wilkinson, and P. Caligari. 1999. Developing an appropriate strategy to assess genetic variability in plant germplasm collections. Theor. Appl. Genet. 98:1125-1131.

Gupta, M., Y. Chyi, J. Romero-Severson, and J. Owen. 1994. Amplification of DNA markers from evolutionarily diverse genomes using single primers of simple-sequence repeats. Theor. Appl. Genet. 89: 998-1006.

Hamann, A., D. Zink, and W. Nagl. 1995. Microsatellite fingerprinting in the genus Phaseolus. Genome 38:507-515.

Hedrick, P., S. Jain, and L. Holden. 1978. Multilocus systems in evolution. Evol. Biol. 11:101-182.

Ibarra-Perez, F., B. Ehdaie, and G. Waines. 1997. Estimation of outcrossing rate in common bean. Crop Sci. 37:60-65.

Jin, Y., T. He, and B.-R. Lu. 2003. Fine scale genetic structure in a wild soybean (Glycine soja) population and the implications for conservation. New Phytol. 159:513-519.

Kantety, R., X. Zeng, J. Bennetzen, and B. Zehr. 1995. Assessment of genetic diversity in Dent and Popcorn (Zea mays L.) inbred lines using Inter-Simple Sequence Repeat (ISSR) amplification. Mol. Breed. 1:365-373.

Kaplan, L., and T. Lynch. 1999. Phaseolus (Fabaceae) in archaeology: AMS radiocarbon dates and their significance for pre-Columbian agriculture. Econ. Bot. 53:261-272.

Koenig, R., and P. Gepts. 1989. Segregation and linkage of genes for seed proteins, isozymes, and morphological traits in common bean (Phaseolus vulgaris). J. Hered. 80:455-459.

Kojima, T., T. Nagaoka, K. Noda, and Y. Ogihara. 1998. Genetic linkage map of ISSR and RAPD markers in einkorn wheat in relation to that of RFLP markers. Theor. Appl. Genet. 96:37-45.

Kornegay, J., C. Cardona, and C.E. Posso. 1993. Inheritance of resistance to Mexican bean weevil in common bean, determined by bioassay and biochemical tests. Crop Sci. 33:589-594.

Lander, E.S., P. Green, J. Abrahamson, A. Barlow, M. Daly, S.E. Lincoln, and L. Newburg. 1987. MAPMAKER: An interactive computer package for constructing primary genetic linkage maps of experimental and natural populations. Genomics 1:174-181.

Liston, A., B.L. Wilson, W.A. Robinson, P.S. Doescher, N.R. Harris, and T. Svejcar. 2003. The relative importance of sexual reproduction versus clonal spread in an aridland bunchgrass. Oecologia 137:216-225.

Martins, M., R. Tenreiro, and M.M. Oliveira. 2003. Genetic relatedness of Portuguese almond cultivars assessed by RAPD and ISSR markers. Plant Cell Rep. 22:71-78.

McClean, P., R. Lee, C. Otto, P. Gepts, and M. Bassett. 2002. Molecular and phenotypic mapping of genes controlling seed coat pattern and color in common bean (Phaseolus vulgaris L.). J. Hered. 93:148-152.

Menendez, C., A. Hall, and P. Gepts. 1997. A genetic linkage map of cowpea (Vigna unguiculata) developed from a cross between two inbred, domesticated lines. Theor. Appl. Genet. 95:1210-1217.

Miranda Colin, S. 1979. Evolucion de Phaseolus vulgaris y P. coccineus, p. 83-99. In E.M. Engleman (ed.) Contribuciones al conocimiento del frijol (Phaseolus) en Mexico. Colegio de Postgraduados, Chapingo, Mexico.

Mumba, L., and N. Galwey. 1998. Compatibility of crosses between gene pools and evolutionary classes in the common bean (Phaseolus vulgaris L.). Genet. Res. Crop Evol. 45:69-80.

Mumba, L., and N. Galwey. 1999. Compatibility between wild and cultivated common bean (Phaseolus vulgaris L.) genotypes of the Mesoamerican and Andean gene pools: Evidence from the inheritance of quantitative characters. Euphytica 108:105-119.

Nagaoka, T., and Y. Ogihara. 1997. Applicability of inter-simple sequence repeat polymorphisms in wheat for use as DNA markers in comparison to RFLP and RAPD markers. Theor. Appl. Genet. 94:597-602.

Nei, M. 1973. Analysis of gene diversity in subdivided populations. Proc. Natl. Acad. Sci. USA 70:3321-3323.

Neigel, J.E. 1997. A comparison of alternative strategies for estimating gene flow from genetic markers. Annu. Rev. Ecol. Syst. 28:105-128.

Nodari, R.O., S.M. Tsai, R.L. Gilbertson, and P. Gepts. 1993. Towards an integrated linkage map of common bean. II. Development of an RFLP-based linkage map. Theor. Appl. Genet. 85:513-520.

Papa, R., and P. Gepts. 2003. Asymmetry of gene flow and differential geographical structure of molecular diversity in wild and domesticated common bean (Phaseolus vulgaris L.) from Mesoamerica. Theor. Appl. Genet. 106:239-250.

Pereira, P.A.A., J.E. de Souza Carneiro, H. Torres da Silva, M.J. Del Peloso, and P. Gepts. 1996. Introgressao de genes de feijoes silvestres em feijao cultivado, p. 393-396. In Anais V Reuniao Nacional de Pesquisa de Feijao, Goiania, GO. 14-18 Oct. 1996. EMBRAPA-CNPAF, Goiania, GO, Brazil.

Prevost, A, and M. Wilkinson. 1999. A new system of comparing PCR primers applied to ISSR fingerprinting of potato cultivars. Theor. Appl. Genet. 98:107-112.

Rohlf, F.J. 1997. NTSYS-pc. Numerical taxonomy and multivariate analysis system, v. 2.02. Exeter Software, Setauket, NY.

Sarla, N., S. Bobba, and E.A. Siddiq. 2003. ISSR and SSR markers based on AG and GA repeats delineate geographically diverse Oryza nivara accessions and reveal rare alleles. Curr. Sci. 84:683-690.

Schneider, S., D. Roessli, and L. Excoffier. 2000. Arlequin ver. 2.000: A software for population genetics data analysis [Online]. available at[verified 12 Nov. 2004]. Genetics and Biometry Laboratory, Univ. of Geneva, Switzerland.

Singh, S.P. 1982. A key for identification of different growth habits of Phaseolus vulgaris L. Annu. Rept. Bean Improv. Coop. 25:92-95.

Singh, S.P., P. Gepts, and D.G. Debouck. 1991. Races of common bean (Phaseolus vulgaris L., Fabaceae). Econ. Bot. 45:379-396.

Singh, S.P., A. Molina, and P. Gepts. 1995. Potential of wild common bean for seed yield improvement of cultivars in the tropics. Can. J. Plant Sci. 75:807 813.

Tanyolac, B. 2003. Inter-simple sequence repeat (ISSR) and RAPD variation among wild barley (Hordeum vulgare subsp, spontaneum) populations from west Turkey. Genet. Res. Crop Evol. 50:611-614.

Tohme, J., D.O. Gonzalez, S. Beebe, and M.C. Duque. 1996. AFLP analysis of gene pools of a wild bean core collection. Crop Sci. 36:1375-1384.

Triana, B., M. Iwanaga, H. Rubiano, and M. Andrade. 1993. A study of allogamy in wild Phaseolus vulgaris. Annu. Rep. Bean Improv. Coop. 36:20-21.

Tsumura, Y., K. Ohba, and S. Strauss. 1996. Diversity and inheritance of inter-simple sequence repeat polymorphism in douglas-fir (Pseudotsuga menziesii) and sugi (Cryptomeria japonica). Theor. Appl. Genet. 92:40-45.

Vanderborght, T. 1983. Evaluation of Phaseolus vulgaris wild types and weedy forms. Plant Genet. Res. Newsl. 54:18-25.

Wells, W.C., W.H. Isom, and J.G. Waines. 1988. Outcrossing rates of six common bean lines. Crop Sci. 28:177-178.

Wolfe, A., Q. Xiang, S. Kephart, and Q. Xiang. 1998. Assessing hybridization in natural populations of Penstemon (Scrophulariaceae) using hypervariable intersimple sequence repeat (ISSR) bands. Mol. Ecol. 7:1107-1125.

Yeh, F.C., R.-C. Yang, T.J.B. Boyle, Z.-H. Ye, and J.X.P. Mao. 1997. POPGENE, the user-friendly shareware for population genetic analysis. Molecular Biology and Biotechnology Centre, Univ. of Alberta, Edmonton, AB, Canada.

Zietkiewicz, E., A. Rafalski, and D. Labuda. 1994. Genome finger-printing by Simple Sequence Repeat (SSR)-anchored polymerase chain reaction amplification. Genomics 20:176-183.

A. Gonzalez, A. Wong, A. Delgado-Salinas, R. Papa, and P. Gepts *

A. Gonzalez and P. Gepts, Dep. of Agronomy and Range Science, Univ. of California, 1 Shields Ave., Davis, CA 95616-8515, USA (A. Gonzalez, current address: CSIRO Plant Industry, Horticulture Unit, PMB 44, Winnellie 0822, Darwin, NT, Australia); A. Wong and A. Delgado-Salinas, Inst. de Biologia, Univ. Nacional Autonoma de Mexico, Mexico, DF, Mexico; R. Papa, Dipartimento di Scienze degli Alimenti, Univ. Politecnica delle Marche, Ancona, Italy. This research was funded by the McKnight Foundation Collaborative Crop Research Program. Received 28 May 2003. Plant Genetic Resources. * Corresponding author (
COPYRIGHT 2005 Crop Science Society of America
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2005 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Title Annotation:Plant Genetic Resources
Author:Gonzalez, A.; Wong, A.; Delgado-Salinas, A.; Papa, R.; Gepts, P.
Publication:Crop Science
Date:Mar 1, 2005
Previous Article:Effects of abscisic acid on antioxidant systems of Stylosanthes guianensis (Aublet) Sw. under chilling stress.
Next Article:Contrasting toxic-endophyte contamination between endophyte-free and Nontoxic-Endophyte tall fescue pastures.

Terms of use | Privacy policy | Copyright © 2021 Farlex, Inc. | Feedback | For webmasters |