Genetic diversity in cultivated carioca common beans based on molecular marker analysis

RESEARCH ARTICLE

Genetic diversity in cultivated carioca common beans based on molecular marker analysis

Juliana Morini Küpper Cardoso Perseguini^I; Alisson Fernando Chioratto^II; Maria Imaculada Zucchi^I; Carlos Augusto Colombo^I; Sérgio Augusto Moraes Carbonell^II; Jorge Mauricio Costa Mondego^I; Rodrigo Gazaffi^III; Antonio Augusto Franco Garcia^III; Tatiana de Campos^{IV, V}; Anete Pereira de Souza^{IV, V}; Luciana Benchimol Rubiano^I

^ICentro de Recursos Genéticos Vegetais, Instituto Agronômico de Campinas, Campinas, SP, Brazil

^IICentro de Grãos e Fibras, Instituto Agronômico de Campinas, Campinas, SP, Brazil

^IIIDepartamento de Genética, Escola Superior de Agricultura "Luiz de Queiroz", Universidade de São Paulo, Piracicaba, SP, Brazil

^IVCentro de Biologia Molecular e Engenharia Genética, Universidade Estadual de Campinas, Campinas, SP, Brazil

^VDepartamento de Genética e Evolução, Instituto de Biologia, Universidade Estadual de Campinas, Campinas, SP, Brazil

^{Send correspondence to}

ABSTRACT

A wide array of molecular markers has been used to investigate the genetic diversity among common bean species. However, the best combination of markers for studying such diversity among common bean cultivars has yet to be determined. Few reports have examined the genetic diversity of the carioca bean, commercially one of the most important common beans in Brazil. In this study, we examined the usefulness of two molecular marker systems (simple sequence repeats - SSRs and amplified fragment length polymorphisms - AFLPs) for assessing the genetic diversity of carioca beans. The amount of information provided by Roger's modified genetic distance was used to analyze SSR data and Jaccards similarity coefficient was used for AFLP data. Seventy SSRs were polymorphic and 20 AFLP primer combinations produced 635 polymorphic bands. Molecular analysis showed that carioca genotypes were quite diverse. AFLPs revealed greater genetic differentiation and variation within the carioca genotypes (Gst = 98% and Fst = 0.83, respectively) than SSRs and provided better resolution for clustering the carioca genotypes. SSRs and AFLPs were both suitable for assessing the genetic diversity of Brazilian carioca genotypes since the number of markers used in each system provided a low coefficient of variation. However, fingerprint profiles were generated faster with AFLPs, making them a better choice for assessing genetic diversity in the carioca germplasm.

Key words: AFLPs, genetic structure, genetic variability, Phaseolus vulgaris L., SSRs.

Introduction

The genus Phaseolus originated in the Americas and contains 55 species, five of which are widely cultivated: P. vulgaris L., P. lunatus L., P. coccineus L., P. acutifolius A., Gray var. latifolius Freeman and P. polyanthus Greenman (Debouck, 1993). Of these five species, the common bean (Phaseolus vulgaris L.) is the most important because its high protein content has made it the major staple food crop in Africa and Latin America (Yu and Bliss, 1978).

The cultivated common bean has two major gene pools and several races within these pools (Beebe et al., 2000). The carioca cream-striped grain type belongs to the Mesoamerican gene pool and shows marked genetic variability that accounts for many traits such as disease resistance. The common bean is essentially an autogamous species with breeding based on self-crossings and the identification of segregant populations from which superior lines are selected. For this reason, it is important to know the extent of parental genetic dissimilarity in order to direct future crosses and enhance the chances of recovering superior genotypes in segregating generations (Gepts and Debouk, 1991).

In 1967, the Agronomic Institute (IAC, Campinas, SP, Brazil) released a new carioca variety named the 'original carioca'. This new cream-striped grain type cultivar quickly became popular and was widely cultivated in Brazil and abroad, especially in Africa, because of its excellent grain quality and high grain yield. Recombination events moved the genetic structure of the original carioca forward and, with time, introduced new variability into the new carioca cultivars that were named after the original carioca because of their similar cream-striped grain type.

Molecular techniques have altered the way plant breeding is being done. Molecular markers have great potential to help breeders develop new improved varieties since they may be used to estimate the genetic diversity and level of heterozygosity among plants and animals (Dani et al., 2008; Kumar et al., 2008). Molecular markers have been used for genetic mapping (Grisi et al., 2007), marker-assisted selection (Ender et al., 2008), and to measure spatial and temporal gene flow within and among populations (Papa and Gepts, 2003).

A range of molecular data can be used to assess crop genetic diversity and has been applied to study the population structure of Mesoamerican and Andean bean gene pools (Dìaz and Blair, 2006; Blair et al., 2007, 2009; Kwak and Gepts, 2009). Talukder et al. (2010) analyzed the genetic diversity and seed mineral content of a set of common beans in order to enhance seed Zn and Fe content. Santalla et al. (2010) used simple sequence repeats (SSRs) to elucidate how adaptation to environmental conditions has sculpted the common bean genomes in southern Europe. In addition, molecular markers have been used to evaluate genetic diversity in snap bean varieties from Europe (Métais et al., 2002; Masi et al., 2003), wild populations from Mexico (Payró de la Cruz et al., 2005) and dry bean genotypes from Italy (Marotti et al., 2007), Bulgaria (Svetleva et al., 2006), Nicaragua (Gomez et al., 2004), Slovenia (Maras et al., 2006) and East Africa (Asfaw et al., 2009).

Among molecular markers, microsatellites or SSRs deserve special attention as tools for analyzing diversity. These codominant, multiallelic markers are widely distributed throughout genomes and can be highly polymorphic (Chin et al., 1996). SSRs have been used to distinguish among Mesoamerican and Andean accessions (Blair et al., 2006; Benchimol et al., 2007; Campos et al., 2007) and have identified greater genetic variability among Andean than among Mesoamerican genotypes (Gepts et al., 2008). SSRs have been successfully used to evaluate genetics, pedigree, phylogeny and/or identify various traits and/or germplasm accessions (McCouch et al., 2001); they have been especially important in assessing the genetic diversity and genetic maps of common beans (Yu et al., 1999; Guo et al., 2000; Métais et al., 2002; Gaitán-Solís et al., 2002; Masi et al., 2003; Blair et al., 2003, 2006; Benchimol et al., 2007; Grisi et al., 2007; Hanai et al., 2007; Zhang et al., 2008).

Amplified fragment length polymorphism (AFLP) is a polymerase chain reaction (PCR)-based molecular marker assay (Vos et al., 1995) that can detect a higher number of polymorphic loci in a single assay than restriction fragment length polymorphism (RFLP) or random amplification of polymorphic DNA (RAPD) (Powell et al., 1996). AFLP is highly polymorphic and shows considerable reproducibility within a laboratory. The effectiveness and reliability of AFLP has led to its increasing use in diversity studies, phylogeny, genomic linkage mapping and identification of varieties (Tohme et al., 1996; Papa and Gepts, 2003; Rosales-Serna et al., 2005). AFLPs can detect a large number of polymorphic bands in a single lane rather than high levels of polymorphism at each locus such as occurs with SSRs. AFLP has been used to distinguish very closely related genotypes belonging to the same commercial class, such as the yellow bean class (Pallottini et al., 2004). AFLPs are also reliable for distinguishing closely related cocoa varieties (Saunders et al., 2001).

Single nucleotide polymorphisms (SNPs), including insertion/deletions (Indels), are informative genetic markers. Although at any given site SNPs could, in principle, involve four different nucleotide variants, in practice they are generally biallelic. When compared to multiallelic markers such as SSRs, biallelic SNPs are less polymorphic (the expected heterozygosity is lower), but this disadvantage is compensated for by the relative abundance of SNPs (Oraguzie et al., 2007). SNP markers are useful in a variety of applications, including the construction of high resolution genetic maps, mapping traits, genetic diagnostics, analysis of the genetic structure of populations and phylogenetic analysis (Rafalski, 2002). Efforts to develop SNP- and Indel-based markers for the common bean (Phaseolus vulgaris L.) have been reported (Galeano et al., 2009).

A comparison of the different marker systems used to estimate crop genetic diversity is important in order to assess their usefulness in germplasm conservation and as plant breeding tools. Garcia et al. (2004) used different marker systems (AFLP, RAPD, RFLP and SSR) to examine the diversity of inbred tropical maize lines and concluded that AFLP was the best molecular assay for fingerprinting and assessing genetic relationships because of its high accuracy. Geleta et al. (2005) reported that both AFLP and SSR markers were efficient tools in assessing the genetic variability among sorghum genotypes.

In this study, we sought to determine the most suitable molecular marker system (SSRs or AFLPs) for assessing the genetic diversity in commercial carioca common beans and examined the genetic variability among various genotypes.

Materials and Methods

Plant material and DNA extraction

Sixty carioca genotypes from the Agronomic Institute Germplasm Bank (IAC, Campinas, SP, Brazil) were used in this study (Table 1). Total genomic DNA was extracted from powdered lyophilized young leaves using the CTAB method (Hoisington et al., 1994).

SSR analysis

A total of 85 SSRs (Table 2) were developed from two microsatellite-enriched libraries, one for the 'CAL-143' line and another for the 'IAC-UNA' variety (Benchimol et al., 2007; Campos et al., 2007; Cardoso et al., 2008; Oblessuc et al., 2009). Sixty-five of the 85 SSRs were genotyped using 6% polyacrylamide silver stained gels (Creste et al., 2001) whereas the remaining 20 SSRs were genotyped using a fluorescent labeling method that allowed high-throughput genotyping (Schuelke, 2000).

To explore the potential nature of the SSR loci assessed in this study, the Gene Ontology functional annotation tool Blast2GO (Conesa et al., 2005) was used to assign GO IDs, enzyme commission (EC) numbers and INTERPRO codes. Sequence similarities were researched using BLASTX (Altschul et al., 1990) against a databank of non-redundant protein sequences (NR, E-value cutoff = 1e^-05) and BLASTN against a databank of expressed sequenced tags (dbEST; E-value cutoff = 1e^-05). Further functional manual annotation was done using AmiGO and PFAM.

For polyacrylamide assays, amplifications were done in a final volume of 25 μL containing 50 ng of DNA, 1X buffer, 0.2 μM of each forward and reverse primer, 100 μM of each dNTP, 2.0 mM MgCl₂, 10 mM Tris-HCl (pH 8,0), 50 mM KCl, and 0.5 U of Taq DNA polymerase. The reactions were run using the following conditions: 1 min at 94 °C, then 30 cycles of 1 min at 94 °C, 1 min at the specific annealing temperature for each SSR and 1 min at 72 °C, with a final cycle of 5 min at 72 °C.

Amplifications using fluorescent SSRs were done with the M13 universal primer. The reactions were amplified in a final volume of 15 μL containing 30 ng of DNA, 1 U of Taq DNA polymerase, 1.5 mM MgCl₂, 0.15 mM of each dNTP, 1X buffer, 0.8 pmol/μL of reverse primer and labeled M13 and 0.2 pmol of forward primer/μL. The reactions were run using the following conditions: 1 min at 94 °C, then 30 cycles of 1 min at 94 °C, 1 min at the specific annealing temperature (56 °C or 60 °C) for each SSR and 1 min at 72 °C, followed by 8 cycles of 1 min at 94 °C, 1 min at 53 °C and 1 min at 72 °C, with a final cycle of 10 min at 72 °C. The amplicons sequenced with a 3730 DNA Analyzer (Applied Biosystems) and analyzed with GeneMapper^® v 3.7 software.

AFLP analysis

AFLP analysis was done as described by Vos et al. (1995), with some modifications. DNA (100 ng) was digested with EcoRI and MseI restriction enzymes, linked to specific adapters, and fragments were amplified. The pre-amplification reactions were done with primers containing one selective nucleotide. The pre-amplification PCR consisted of 3 μL of digested DNA and with the adapters already linked, 0.2 μM of the primer combinations, 100 μM of each dNTP, 2.0 mM MgCl₂, 10 mM Tris-HCl (pH 8,0), 50 mM KCl, and 0.5 U of Taq DNA polymerase in a final volume of 25 μL. The reactions were run using the following conditions: 1 min at 94 °C, then 26 cycles of 1 min at 94 °C, 1 min at 56 °C, 1 min at 72 °C, with a final cycle of 5 min at 72 °C. The pre-amplification products were diluted 1:9 in water.

Selective amplification was done with primers that had three selective nucleotides (for EcoRI and MseI primers). Twenty primer combinations were used in this step (Table 3). PCR for selective amplification was done as described for the pre-amplification reaction, using 3 μL of the pre-amplification dilution. The reactions were done using the following conditions: 2 min at 94 °C, then 12 cycles of 0.30 s at 94 °C, 0.30 s at 65 °C and 1 min at 72 °C, followed by 23 cycles of 0.30 s at 94 °C, 0.30 s at 56 °C and 1 min at 72 °C, with a final cycle of 5 min at 72 °C. The PCR amplification products were separated on 7% denaturing polyacrylamide silver stained gels. The AFLP bands were scored manually.

Data analysis

The AFLP data were scored for presence (1) or absence (0) of bands, whereas for SSRs the results were transformed into genotypic data in order to identify locus and allele frequencies. Pair-wise comparisons were used to estimate Jaccard's similarity coefficient (Jaccard, 1908) for AFLP data using NTSYS software, version 2.02E (Rohlf, 1993). Genetic distances (GDs) were calculated from the SSR data for all possible inbred pairs using Rogers modified genetic distance (RMD; Goodman and Stuber, 1983) version 1.3 (Miller, 1997). Cluster analyses for SSRs and AFLPs were done using UPGMA (Unweighted Pair-Group Method with Arithmetic Averages).

Polymorphism information content (PIC) values were calculated using the formula:

where fi is the frequency of the i^th allele (marker) for the i^th SSR locus (Lynch and Walsh, 1998). Discrimination power (DP) analysis values for the k^th primer were calculated using the formula:

where N is the number of individuals and p_j is the frequency of the j^th pattern (Tessier et al., 1999). PIC was used to measure the information of a given marker locus for the pool of genotypes, while DP was used to measure the efficiency of SSRs and AFLPs in order to identify varieties by taking into account the probability of two randomly chosen individuals having different patterns.

The bootstrap procedure (Efron and Tibshirani, 1993) was used to verify whether the number of polymorphic SSRs and AFLPs used to estimate genetic similarity was large enough to supply a precise estimation of molecular markers among the genotypes (Tivang et al., 1994). The polymorphic markers were submitted to sampling with the replacement of markers to create new samples from the original data. Genetic similarities for each of these subsets were calculated by obtaining 1000 bootstraps estimates of SSRs and AFLPs for each of these combinations. The coefficients of variation (CV) were used to construct box plots for each sample size. These analyses were done using the R software, which is a language and environmental asset for statistical computing. For each marker system (SSRs and AFLPs), the exponential function was adjusted to estimate the number of loci needed to obtain a CV of 10%. The median and maximum CV values were used to evaluate the accuracy of the genetic distance estimates because although the mean CV is often used in the literature, caution is needed when dealing with molecular marker data for which there is no assurance that the CV values have a symmetrical distribution.

Principal coordinate analysis (PCO; Gower, 1966) was done with the SSR MRD distance matrix and the first three principal coordinates were used to describe the dispersion of the 60 accessions according to their allele data.

Wrights F statistics for SSRs were estimated using the GDA program (Lewis and Zaykin, 2000). This set of statistics was used to test the structure of genetic diversity of the carioca genotypes evaluated. AMOVA (Analysis of Molecular Variance) was used to test the structure of the genetic diversity of the genotypes based on the AFLP data, and the analyses were done using Arlequin 2.0 software (Schneider and Excoffier, 1999). The significance of the fixation indices was tested by a permutation procedure with 10,000 permutations. Arlequin 2.0 software was also used to estimate the diversity fraction (F_ST) generated by AFLP analysis. The genetic differentiation coefficient, F_ST, measures the relative degree of gene differentiation among subpopulations such that each genotype can be considered a subpopulation. Wright's F_ST is considered to be identical to G_ST (Nei, 1978) for fully homozygous diploids in AFLP analysis and was calculated as: F_ST = G_ST = D_ST/h_T, where H_T = H_S + D_ST, and D_ST = (H_T - H_S).

Results and Discussion

Survey of molecular marker polymorphism

Seventy (82.4%) of the SSRs (Table 2) were polymorphic and produced 196 polymorphic alleles. The allele number ranged from 2 to 6 with a mean of 2.8 alleles per SSR. The highest numbers of alleles observed were found for SSR-IAC10, SSR-IAC62, SSR-IAC144 and SSR-IAC251. AFLP markers produced 725 bands of which 635 (87.6%) were polymorphic. The number of fragments amplified per primer pair varied from 20 (E-ACC/M-CGG) to 71 (E-TAA/M-GAA) (Table 3).

The highest polymorphic information content (PIC) found for SSRs was 0.70 (SSR-IAC144) and the lowest was 0.03 (SSR-IAC155), with a mean value of 0.47. The highest PIC value for AFLP was 0.37 (E-AAG/M-CGG) and the lowest was 0.03 (E-TAT/M-GTG), with a mean value of 0.29. Benchimol et al. (2007) analyzed genotypes from the Andean and Mesoamerican gene pools (P. vulgaris L.) with 123 polymorphic microsatellites and found PIC values from 0.04 to 0.83. These genotypes were expected to be more diverse than those of the carioca commercial type (Mesoamerican). Fifty-three SSRs used by Benchimol et al. (2007) were also used here with the carioca genotypes, and six of these SSRs showed similar PIC values (SSR-IAC-01, SSR-IAC32, SSR-IAC55, SSR-IAC64, SSR-IAC70 and SSR-IAC-83). Thus, most of the SSRs differed in their PIC, a situation that probably reflected the type and number of genotypes that were evaluated. Moreover, the higher PIC values found for the SSR markers were related to the size of the SSR motifs; perfect SSRs with a higher number of repetitions and compound SSRs showed higher PIC (Table 2). The PIC values were compatible with those obtained for SSRs used in genotyping carioca cultivars.

Discrimination power (DP) analyses for SSRs (Table 2) yielded values from 0.01 (SSR-IAC155) to 0.82 (SSR-IAC68) with most being around 0.40. Based on parameters such as the number of alleles, DP and PIC, the polymorphism estimated for the carioca genotypes with a subset of 85 SSRs showed that these loci could generate sufficient polymorphism to allow their use as molecular markers to establish genetic relationships among very closely related genotypes. These SSRs could be also useful for a wide range of genetic investigations such as linkage map construction and association mapping studies. The DP values for AFLPs were higher than for SSRs, ranging from 0.70 (E-TCA/M-GAC) to 0.94 (E-ACC/M-CCC), with most of the values being around 0.92 (Table 3).

To understand the genomic nature (functional or non-functional regions) of the SSR markers used, the 85 contigs were compared with non-redundant GenBank databases using BLASTX and BLASTN (dbEST). Of the sequences analyzed, 22% were found in the BLASTX - NR DB and 49% in the BLASTN dbEST database. GO analyses yielded a hierarchy of terms that varied at many levels and were divided into functional groups identified as "molecular function", "biological process" and "cellular component". Nineteen functional SSRs were identified (Table S1), and the SSRs that shared similarity with ESTs deposited in dbEST are shown in Table S2.

Genetic diversity of the carioca genotypes

Various studies have shown that the accuracy of genetic distance measurements is enhanced by using a mean CV of 10% (dos Santos et al., 1994; Halldén et al., 1994; Thormann et al., 1994; Tivang et al., 1994). However, according to Garcia et al. (2004), the use of an average value is not a good indicator of central tendency for skewed data. Hence, in this work, the minimum number of SSR and AFLP loci necessary for an accurate representation of genetic distances was calculated in order to construct an exponential function based on the average, median and maximum CV values of genetic distances obtained by bootstrap sampling data for each SSR and AFLP. The CV values used to calculate the median were used as follows: the sample size (number of loci) required for 50% of the genetic distances to have CV values < 10% (n_median), the sample size required for no genetic distance to have a CV > 10% (n_maximum), and the sample size required for all genetic distances to have an average CV of 10% (n_mean) (Figure 1). The results obtained based on these criteria are shown in Table 4.

The choice of the appropriate number of polymorphic loci required for a reliable estimate of genetic distance is influenced by the criteria used. According to Garcia et al. (2004), the median CV value is the best choice for evaluating the precision of genetic distance estimates based on different molecular systems. The extremely high (almost 100%) coefficients of determination for the adjusted equations for SSRs and AFLPs indicated that extrapolation to outlying points could be done. The results of the 10% CV analysis (n_median) showed that 44 SSRs were necessary to assess the genetic variability among genotypes, whereas 100 AFLPs were required for the same analysis (Table 4). Thus, AFLPs were more suitable for analyzing genetic diversity in the set of common bean genotypes evaluated here since they provided information more rapidly and more accurately. Similar results were reported by Pallottini et al. (2004) for AFLPs used to assess the yellow bean class.

Maluf et al. (2005) characterized the genetic diversity of Coffea arabica genotypes using various molecular markers (AFLP, RAPD and SSR). RAPD and SSR were more efficient in kinship analysis; however, despite the high discriminatory power of SSRs, these markers were not efficient in separating the genotypes into well-defined groups. The authors concluded that AFLP markers were more appropriate for assessing the genetic diversity of kin groups. More recently, Kumar et al. (2008) concluded that AFLP markers were very useful for assessing the genetic diversity among narrow-based common bean accessions.

Of the 70 SSRs used here, 37 were mapped in the 'IAC-UNA' x 'CAL-143' genetic map (UC map; Campos et al., 2010). These SSRs were anchored in eight out of the 11 linkage groups (GLs) generated in the UC map (B1 - SSR-IAC01, SSR-IAC21 and SSR-IAC226; B6 - SSR-IAC47, SSR-IAC128 and SSR-IAC183; B9 - SSR-IAC55, SSR-IAC62 and SSR-IAC242). The linkage group B2 had 10 SSRs (SSR-IAC18a, SSR-IAC24, SSR-IAC46, SSR-IAC51, SSR-IAC57, SSR-IAC70, SSR-IAC134, SSR-IAC141, SSR-IAC166 and SSR-IAC251), the B3 linkage group had one SSR (SSR-IAC77), with one SSR in the B8 linkage group (SSR-IAC22). The B4 linkage group had four SSRs (SSR-IAC25, SSR-IAC66, SSR-IAC67 and SSR-IAC179), with the same number of SSRs being observed in the B5 linkage group (SSR-IAC10, SSR-IAC88, SSR-IAC96 and SSR-IAC159). The B7 linkage group had six SSRs (SSR-IAC18b, SSR-IAC64, SSR-IAC101, SSR-IAC143, SSR-IAC272 and SSR-IAC239) and the B10 linkage group had two SSRs (SSR-IAC155 and SSR-IAC244). These results show that the number of SSRs used was not sufficient to scan the entire bean genome since not all of the 11 chromosomal pairs were assessed and some linkage groups were irregularly scanned. Indeed, the SSR map positions in the UC map were not known at the time these SSRs were chosen for this study.

The SSR dendrogram (Figure 2) showed high genetic variability, with genetic distances varying from 0.37 to 0.63. However, a low genetic structure was observed as some groups were not well-defined. This could be explained by the fact that the SSRs used to assess the genetic variability of common genotypes in 'carioca' beans did not cover the entire genome of the common bean since the SSRs used in the analysis were not distributed in all linkage groups of the common bean.

The AFLP dendrogram provided a better clustering pattern for the carioca genotypes, which formed two major groups (Figure 3); five accesses ('Campeão II', 'LP9979', 'IAC-Aruã', 'Aporé' and 'CariocaMG') could not be classified with either of the two major groups. Genetic distances varied from 0.09 to 0.63 and indicated marked genetic variability. The AFLP average genetic distance (GD_AFLP = 0.88) was higher than the average genetic distance for SSRs (GD_SSR = 0.58). Some of the genotypes clustered according to the institutional breeding program from which they were derived. This was observed for the IAPAR (The Agronomic Institute of the State of Paraná, in southern Brazil) genotypes ('IAPAR 81', 'IAPAR57', 'IAPAR 14', 'IAPAR 80', 'Juriti' and 'LP88-175'), for the IAC (The Agronomic Institute, Campinas, SP, Brazil) genotypes ('IAC-Carioca', 'IAC-Akytã' and 'IAC-Pyatã'; 'IAC-Tybatã', 'IAC-Votuporanga' and 'IAC-Ybaté') and for those from EMBRAPA ('A-449', 'Branquinho', 'BRS-Horizonte', 'BRS-Cometa', 'BRS-Pontal', 'BRS-Requinte' and 'BRSMG-Talismã'). This clustering pattern could be explained by the fact that each breeding Institution usually displays and uses a different genetic core collection and each breeding program has its own goals.

The Pearson correlation coefficient for SSR and AFLP genetic distances was negative (r = -0.08) indicating lack of correlation between these marker systems. The differences between SSRs and AFLPs reflect the extent of genome coverage and different evolutionary properties (Tautz and Schlötterer, 1994). The genomic distribution of SSRs is non-random (Li et al., 2004), whereas it is reasonable to expect that AFLPs should have complete genome coverage (Nowosielski et al., 2002). Saini et al. (2004) reported a low correlation (0.50) between these two classes of markers in rice and assumed that this reflected different genomic fractions that involved repeat and/or unique sequences; these sequences may have had different patterns of involvement or preservation during natural or human selection.

Menezes et al. (2004) and Palomino et al. (2005) studied the genetic diversity within carioca cultivars using RAPD markers and verified high genetic variability within carioca cultivars that differed from the 'original carioca'. Although RAPD markers have been extensively used for common bean fingerprinting and genetic analysis, AFLP analysis provides a higher level of polymorphism than RAPD (Pejic et al., 1998). AFLP markers are advantageous because they reveal a larger number of reproducible markers, thereby increasing the probability of identifying polymorphic markers, even among closely related genotypes such as in common beans (Tohme et al., 1996; Beebe et al., 2001).

Principal coordinate analysis (PCO) failed to provide a clear division of the genotypes in the first three axes for both markers. For SSRs, 20.95% of the total variance was explained in three axes (Figure 4). According to SSR PCO, the 'original carioca', which is the ancestor of most of the carioca genotypes that are currently being used, was detached from the other cariocas. This finding supports the hypothesis that carioca genotypes have diverged from the 'original carioca' ancestor and shows that there is still considerable genetic variability to be exploited among carioca common bean accessions. Indeed, carioca common bean accessions have derived from many recombination events and from crosses with different genotypes, with much of this variability being preserved. Other factors, such as retrotransposons, may also be involved in generating and maintaining carioca genetic diversity.

GDA analysis of the SSR data attributed 83% of the total molecular variance to F_ST, indicating a high level of genetic differentiation in the carioca samples. Based on the AFLP data, AMOVA showed 98% G_ST. Meanwhile, the observed heterozygosity was very low for all genotypes (0.16) and reflected the inbreeding index characteristic of an autogamous species (Beebe et al., 2000). Both groups of markers revealed high genetic variability among the carioca genotypes.

Common bean breeding programs usually explore only a low proportion of the available genetic diversity. Since most breeding programs use only a limited number of cultivars, the genetic base of elite cultivars tends to be narrow (Cooper et al., 2001). However, our findings suggest that genetic diversity is still preserved in the carioca collection. The degree of variability detected here supports the idea that the carioca group has a broader genetic base than expected, with a large number of unique alleles and low gene flow among them.

The success of common bean breeding programs is intimately related to the appropriate choice of divergent parents. To make this choice researchers need to know the genetic diversity contained in the available germplasm since the best parental combination can be obtained by combining two parents with a high degree of divergence; this will allow the exploitation of heterosis to improve breeding programs. The 'carioca' genotypes of common beans are currently the most sold in Brazil, largely because of their agronomic and commercial characteristics. Most of the genotypes have a high yield and resistance to various bean diseases, such as rust (examples: 'Aporé', 'Carioca MG', 'Rudá'), golden mosaic virus (examples: 'IAPAR 72', 'MD-806' and 'IPR-Aurora'), to bacterial diseases (examples: 'IAC-Tybatã', 'LP88-175', 'IAC-Carioca Pyatã', 'IAC-Carioca Akytã', 'IAC-Ybaté', 'BRS-Pontal'), anthracnose (examples: 'IAC-Carioca Tybatã', 'IAC-Carioca Pyatã', 'IAC-Carioca Akytã', 'IAC-Ybaté', 'H96A28', 'H96A102', 'IAPAR-31', 'BRS-Requinte', 'BRS-Pontal', 'IAC-Alvorada', 'H96A31', 'LP 9979') and angular leaf spot (examples: 'BRS-Pontal', 'IAC-Carioca Tybatã', 'Pérola', 'IAPAR-31', 'IAC-Carioca Aruã', 'BRS-Requinte'). Knowledge of the appropriate agronomic information and genetic distances among a set of 'carioca' genotypes can be used by bean breeding programs to assist in selecting superior parents and allow the exploitation of heterosis.

Information on possible crosses based on genetic diversity can also be used to assist breeding programs. Thus, for example, based on the AFLP genetic distances and the agronomic traits associated with each genotype, the following crosses could be suggested: 'IAC-Alvorada x H96A102' (GD_AFLP = 0.89), 'IAC-Alvorada' x 'LP 0181' (GD_AFLP = 0.87), 'IAC-Alvorada' x 'LP 9979' (GD_AFLP = 0.89), 'IAC-Alvorada' x 'H96A28' (GD_AFLP = 0.87), 'IAC-Alvorada' x 'IPR-Aurora' (GD_AFLP = 0.90), 'BRS-Requinte' x 'H96A102' (GD_AFLP = 0.90), 'BRS-Pontal' x 'H96A102' (GD_AFLP = 0.87), 'Branquinho' x 'H96A102' (GD_AFLP = 0.85), 'BRS-Talismã' x 'H96A102' (GD_AFLP = 0.83) and 'H96A31' x 'H96A102'(GD_AFLP = 0.80).

A large number of markers have been used to assess genetic diversity in plants. For wide-scale use in germplasm characterization and breeding, it is important that these marker technologies be exchanged between laboratories and be standardized to yield reproducible results that will allow the direct comparison of data among laboratories and studies, thereby reducing project expenses. Jones et al. (1997) described a network experiment involving several European laboratories in which the reproducibility of three popular molecular marker techniques (AFLP, RAPD and SSR) was examined. The authors concluded that the costs involved in developing SSR markers meant that AFLP was more useful for studying genetic diversity - the major limitation inherent with SSRs resides in the logistic difficulty of increasing the number of useful loci for assignment tests. Indeed, developing and applying large numbers of SSRs may be technically challenging, expensive and time-consuming (Goldstein and Pollock, 1997).

AFLP provides a practical alternative for studying genetic diversity since the technique generates a large number of loci, and the cost and time required are lower than for SSRs (Vos et al., 1995; Rieseberg, 1998; Mueller and Wolfenbarger, 1999). Furthermore, the high reproducibility of AFLP means that multiple bands can be generated in a single assay (Oliveira et al., 2004). AFLP has been successfully used to estimate genetic diversity in the common bean (Tohme et al., 1996; Caicedo et al., 1999; Maciel et al., 2003; Pallottini et al., 2004; Rosales-Serna et al., 2005).

Comparisons of the efficiency of SSRs and AFLPs in assessing the genetic diversity of plants have generally shown that AFLP is the preferred technique. Barbosa et al. (2003) investigated the genetic diversity of tropical maize using SSRs and AFLPs and concluded that the former was less efficient than the latter for assigning lines to heterotic groups in tropical maize, and for predicting single cross performance in this culture. Similarly, Garcia et al. (2004), in a comparison of AFLP, RAPD, RFLP and SSR markers to evaluate genetic diversity in tropical maize, found that AFLP was the best-suited molecular assay for accurate fingerprinting and assessment of genetic relationships among tropical maize inbred lines.

The advantages and disadvantages of different marker systems for analyzing diversity in breeding populations or germplasm bank accessions, such as the carioca commercial type, must be considered in genotyping data for crop improvement and for the ex-situ conservation of plant genetic resources. Other studies that have compared the efficiencies of AFLPs and SSRs have concluded that both marker systems are useful for assessing the genetic structure and diversity of common bean accessions (Maras et al., 2008; Masi et al., 2009).

In conclusion, this is the first study to use powerful molecular markers such as AFLP and SSRs to assess the genetic variability of carioca commercial beans. Other studies of cream-striped carioca genotypes have used RAPD markers, which are less powerful in revealing the extent of genetic diversity (Menezes et al., 2004; Palomino et al., 2005). AFLP and SSRs and the number of markers they generated were appropriate for assessing the genetic diversity among carioca genotypes. In practice, AFLP markers were easier to use when screening the whole genome for genetic diversity and were more suitable in distinguishing among very closely related genotypes belonging to the same commercial class, such as the carioca genotypes.

Acknowledgments

The authors thank Dr. Marcos Machado (IAC, Cordeirópolis, SP, Brazil) for making available the automatic sequencer. This research was supported by grants from Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP, grants 2006/61848-3 and 2006/59354-2).

Internet Resources

Supplementary Material

The following online material is available for this article:

Table S1 - Functional annotation off SSRs sharing similarity with proteins in GenBank.

Table S2- SSRs with similarity to ESTs deposited in the dbEST databank.

This material is available as part of the online article from http://www.scielo.br/gmb

Received: January 25, 2010; Accepted: July 29, 2010.

Associate Editor: Dario Grattapaglia

License information: This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Anexos

Table S1 -Clique para ampliar

Table S2 -Clique para ampliar

• Altschul SF, Gish W, Miller W, Myers EW and Lipman DJ (1990) Basic local alignment search tool. J Mol Biol 215:403-410.
• Asfaw A, Blair MW and Almekinders C (2009) Genetic diversity and population structure of common bean (Phaseolus vulgaris L.) landraces from the East African highlands. Theor Appl Genet 120:1-12.
• Barbosa AMM, Geraldi IO, Benchimol LL, Garcia AAF, Souza Jr CL and Souza AP (2003) Relationship of intra and interpopulation maize single crosses hybrid performance and genetic distances computed from AFLP and SSR markers. Euphytica 130:87-99.
• Beebe S, Skroch PW, Tohme J, Duque MC, Pedraza F and Nienhuis J (2000) Structure of genetic diversity among common bean landraces of Middle American origin based on correspondence analysis of RAPD. Crop Sci 40:264-273.
• Beebe S, Rengifo J, Gaitán E, Duque MC and Tohme J (2001) Diversity and origin of Andean landraces of common bean. Crop Sci 41:854-862.
• Benchimol LL, Campos T, Carbonell SAM, Colombo CA, Chioratto AF, Formighieri EF and Souza AP (2007) Structure of genetic diversity among common bean (Phaseolus vulgaris L.) varieties of Mesoamerican and Andean origins using new developed microsatellite markers. Genet Resour Crop Evol 54:1747-1762.
• Blair MW, Pedraza F, Buedia HF, Gaitán-Solís E, Beebe SE, Gepts P and Tohme J (2003) Development of a genome-wide anchored microsatellite map for common bean (Phaseolus vulgaris L.). Theor Appl Genet 107:1362-1374.
• Blair MW, Giraldo MC, Buendía HF, Tovar E, Duque MC and Beebe SE (2006) Microsatellite marker diversity in common bean (Phaseolus vulgaris L.). Theor Appl Genet 113:100-109.
• Blair MW, Dìaz JM, Hidalgo R, Dìaz LM and Duque MC (2007) Microsatellite characterization of Andean races of common bean (Phaseolus vulgaris L.). Theor Appl Genet 116:29-43.
• Blair MW, Dìaz LM, Buendia HF and Duque MC (2009) Genetic diversity, seed size associations and population structure of a core collection of common beans (Phaseolus vulgaris L.). Theor Appl Genet 119:955-972.
• Caicedo AL, Gaitán E, Duque MC, Chica OT, Debouck DG and Tohme J (1999) AFLP fingerprinting of Phaseolus lunatus L. and related wild species from South America. Crop Sci 39:1497-1507.
• Campos T, Benchimol LL, Carbonell SAM, Chioratto AF, Formighieri EF and Souza AP (2007) Microsatellites for genetic studies and breeding programs in common bean. Pesq Agropec Bras 42:401-408.
• Campos T, Oblessuc PR, Sforça DA, Cardoso JMK, Baroni RM, Sousa ACB, Carbonell SAM, Chioratto AF, Garcia AAF, Rubiano LB et al. (2010) Inheritance of growth habit detected by genetic linkage analysis using microsatellites in the common bean (Phaseolus vulgaris L.). Mol Breed, Epub.
• Cardoso JMK, Oblessuc PR, Campos T, Sforça DA, Carbonell SAM, Chioratto AF, Formighieri EF, Souza AP and Benchimol LL (2008) New microsatellite markers developed from an enriched microsatellite common bean library. Pesq Agropec Bras 7:929-936.
• Chin ECL, Senior ML and Shu H (1996) Maize simple repetitive DNA sequences: Abundance and allele variation. Genome 39:866-873.
• Conesa A, Götz S, García-Gómez JM, Terol J, Talón M and Robles M (2005) Blast2GO: A universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 18:3674-3676.
• Cooper HD, Spillan EC and Hodgkin T (2001) Broadening the genetic base of crops: An overview. In: Cooper HD, Spillane C and Hodgkin T (eds) Broadening the Genetic Base of Crop Production. CAB International, Wallingford, pp 1-23.
• Creste S, Tulmann A and Figueira A (2001) Detection of single sequence repeat polymorphism in denaturing polyacrylamide sequencing gels by silver staining. Plant Mol Bio Rep 19:299-306.
• Dani MAC, Heinneman MB and Dani SU (2008) Brazilian Nelori cattle: A melting pot unfolded by molecular genetics. Genet Mol Res 7:1127-1137.
• Debouck DG, Toro O, Paredes OM, Johnson WC and Gepts P (1993) Genetic diversity and ecological distribution of Phaseolus vulgaris in northwestern South America. Econ Bot 47:408-423.
• Díaz LM and Blair MW (2006) Race structure within the Mesoamerican gene pool of common bean (Phaseolus vulgaris L.) as determined by microsatellite markers. Theor Appl Genet 114:143-154.
• dos Santos JB, Nienhuis J, Skroch P, Tivang J and Slocum MK (1994) Comparison of RFLP genetic markers in determining genetic similarity among Brassica oleracea L. genotypes. Theor Appl Genet 87:909-915.
• Efron B and Tibshirani RJ (1993) An Introduction to the Bootstrap. v. 57 of Monographs on Statistics and Applied Probability. Chapman and Hall, New York.
• Ender M, Terpstra K and Kelly JD (2008) Marker assisted selection for white mold resistance in common bean. Mol Breed 2:149-157.
• Gaitan-Solís E, Duque MC, Edwards KJ and Tohme J (2002) Microsatellite in common bean (Phaseolus vulgaris): Isolation, characterization, and cross-species amplification in Phaseolus sp. Crop Sci 42:2128-2136.
• Galeano CH, Fernández AC, Gómez M and Blair MW (2009) Single strand conformation polymorphism based SNP and Indel markers for genetic mapping and synteny analysis of common bean, (Phaseolus vulgaris L.) BMC Genomics 10:629.
• Garcia AAF, Benchimol LL, Barbosa AMM, Geraldi IO, Júnior CLS and Souza AP (2004) Comparison of RAPD, RFLP, AFLP and SSR markers for diversity studies in tropical maize inbred lines. Genet Mol Biol 27:579-588.
• Geleta LF, Labuschagne MT and Viljoen CD (2005) Genetic variability in pepper (Capsicum annuum L.) estimated by morphological data and amplified fragment length polymorphism markers. Biodivers Conserv 14:2361-2375.
• Gepts P and Debouk D (1991) Origin, domestication, and evolution of the common bean (Phaseolus vulgaris L.). In: Van Schoonhoven A and Voyset O (eds) Common Beans: Research for Crop Improvement. CAB International, Wallingford, pp 7-53.
• Gepts P, Aragão F, de Barros E, Blair MW, Brondani R, Broughton W, Galasso I, Hernández G, Kami J, Lariguet P, et.al (2008) Genomics of Phaseolus beans, a major source of dietary protein and micronutrients in the tropics. In: Moore P and Ming R (eds) Genomics of Tropical Crop Plants, Springer, Heidelberg, pp 113-143.
• Goldstein DB and Pollock DD (1997) Launching microsatellites: A review of mutation processes and methods of phylogenetic inference. J Hered 88:335-342.
• Goodman MM and Stuber CW (1983) Races of maize VI. Isozyme variation among races of maize in Bolívia. Maydica 28:169-187.
• Gower JC (1966) Some distant properties of latent root and vector methods used in multivariate analysis. Biometrika 53:325-338.
• Gómez OJ, Blair MW, Frankow-Lindberg BE and Gullberg U (2004) Molecular and phenotypic diversity of common bean landraces from Nicaragua. Crop Sci 44:1412-1418.
• Grisi MCM, Blair MW, Gepts P, Brondani C, Pereira PAA and Brondadi RPV (2007) Genetic mapping of a new set of microsatellite markers in a reference common bean (Phaseolus vulgaris) population BAT93 x Jalo EEP558. Genet Mol Res 3:691-706.
• Guo JC, Hu XW, Yanagihara S and Yoshinobu E (2000) Isolation and characterization of microsatellites in snap bean. Acta Bot Sin 42:1179-1183.
• Halldén C, Nilsson NO, Rading IM and Säll T (1994) Evaluation of RFLP and RAPD markers in comparison of Brassica napus breeding lines. Theor Appl Genet 88:123-128.
• Hanai LL, Campos T, Camargo LEA, Benchimol LL, de Souza AP, Melotto M, Carbonell SA, Chioratto AF, Consoli L, Formighieri EF, et. al (2007) Development, characterization and comparative analysis of polymorphism at common bean-SSR loci isolated from genic and genomic sources. Genome 50:266-277.
• Hoisington D, Khairallah M and Gonzalez-de-Leon D (1994) Laboratory Protocols: CIMMYT Applied Molecular Genetics Laboratory. 2nd edition. CIMMYT, Mexico DF, pp 51.
• Jaccard P (1908) Nouvelles rechearches sur la distribuition florale. Bull Soc Vaudoice Sci Nat 44:223-270.
• Jones CJ, Edwards KJ, Castaglione S, Winfield MO, Sala F, Wiel C, Bredemeijer G, Vosman B, Matthes M, Daly A, et. al (1997) Reproducibility testing of RAPD, AFLP and SSR markers in plants by a network of European laboratories. Mol Breed 3:381-390.
• Kumar V, Sharma S, Kero S, Sharma S, Sharma K, Kumar M and Bhat KV (2008) Assessment of genetic diversity in common bean (Phaseolus vulgaris L.) germplasm using amplified fragment length polymorphism (AFLP). Scientia Hort 116:138-143.
• Kwak M and Gepts P (2009) Structure of genetic diversity in the two major gene pools of common bean (Phaseolus vulgaris L., Fabaceae). Theor Appl Genet 118:979-992.
• Li YC, Korol AB, Fahima T and Nevo E (2004) Microsatellites within genes: Structure, function, and evolution. Mol Biol Evol 21:991-1007.
• Lynch M and Walsh JB (1998) Genetics and Analysis of Quantitative Traits. Sinauer Associates, Sunderland, 980 pp.
• Maciel FL, Echeverrigaray S, Gerald LTS and Gobbi-Grazziotin F (2003) Genetic relationships and diversity among Brazilian cultivars and landraces of common beans (Phaseolus vulgaris L.) revealed by AFLP markers. Gen Res Crop Evol 50:887-893.
• Maluf MP, Silvestrine M, Ruggiero LMC, Guerreiro-Filho O and Colombo CA (2005) Genetic diversity of cultivated Coffea arabica inbred lines assessed by RAPD, AFLP and SSR marker systems. Sci Agric 62:366-373.
• Maras M, Susnik S, Sustar-Vozlic J and Meglic V (2006) Temporal changes in genetic diversity of common bean (Phaseolus vulgaris L.) accessions cultivated between 1800 and 2000. Russian J Genet 42:775-782.
• Maras M, Sustar-Vozlic J, Javornik B and Meglic V (2008) The efficiency of AFLP and SSR markers in genetic diversity estimation and gene pool classification of common bean (Phaseolus vulgaris L.). Acta Agric Slov 91:87-96.
• Marotti I, Bonetti A, Minelli M, Catizone P and Dinelli G (2007) Characterization of some Italian common bean (Phaseolus vulgaris L.) landraces by RAPD, semi-random and ISSR molecular markers. Genet Resour Crop Evol 54:175-188.
• Masi P, Spagnoletti-Zeuli PL and Donini P (2003) Development and analysis of multiplex microsatellite markers sets in common bean (Phaseolus vulgaris L.). Mol Breed 11:303-313.
• Masi P, Logozzo G, Donini P and Spagnoletti ZP (2009) Analysis of genetic structure in widely distributed common bean landraces with different plant growth habits using SSR and AFLP markers. Crop Sci 49:187-199.
• McCouch SR, Temnykh S, Lukashova A, Coburn J, Declerck G, Cartinhour S, Harrington S, Thomson M, Septiningsi E, Semon M, et. al (2001) Microsatellite markers in rice: Abundance, diversity and applications. In: Khush GS, Brar DS and Hardy B (eds) Rice Genetics IV. IRRI, Manila, pp 117-135.
• Menezes MVCO, Santos JB and Menezes M (2004) Genetic diversity in Carioca and Pérola cultivars of common bean based on RAPD markers. Crop Breed Appl Biotechnol 4:178-182.
• Métais I, Hamon B, Jalouzot R and Peltier D (2002) Structure and level of genetic diversity in various bean types evidenced with microsatellite markers isolated from a genomic enriched library. Theor Appl Genet 104:1346-1352.
• Miller M (1997) TFPGA - Tools for Population Genetic Analyses, v. 1.3. Northern Arizona University, Flagstaff, pp 33.
• Mueller UG and Wolfenbarger LL (1999) AFLP genotyping and fingerprinting. Trends Ecol Evol 14:389-394.
• Nei M (1978) Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics 89:583-590.
• Nowosielski J, Podyma W and Nowosielska D (2002) Molecular research on the genetic diversity of Polish varieties and landraces of Phaseolus coccineus L. and Phaseolus vulgaris L. using the RAPD and AFLP methods. Cell Mol Biol Lett 7:753-762.
• Oblessuc PO, Campos T, Cardoso JMK, Sforça DA, Baroni RM, Souza AP and Benchimol LL (2009) Adaptation of fluorescent technique for genotyping with new microsatellite markers in common bean. Pesq Agropec Bras 44:638-644.
• Oliveira KM, Laborda PR, Garcia AAF, Zagatto-Paterniani MEAG and Souza AP (2004) Evaluating genetic relationships between tropical maize inbred lines by means of AFLP profiling. Hereditas 140:24-33.
• Oraguzie NC, Rikkerink EHA, Gardiner SE and De Silva HN (2007) SNP Applications in Plants - Association Mapping in Plants. Springer, New York, pp 277.
• Palomino EC, Mori ES, Zimback L, Tamborussi EV and Moraes CB (2005) Genetic diversity of common bean genotypes of Carioca commercial group using RAPD markers. Crop Breed Appl Biotechnol 5:80-85.
• Pallottini L, Garcia E, Kami J, Barcaccia G and Gepts P (2004) The genetic anatomy of a patented yellow bean. Crop Sci 44:968-977.
• Papa R and Gepts P (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.
• Payró de la Cruz E, Gepts P, Colunga-GarciaMarín P and Zizumbo-Villareal D (2005) Spatial distribution of genetic diversity in wild populations of Phaseolus vulgaris L. from Guanajuato and Michoacán, Mexico. Genet Resour Crop Evol 52:589-599.
• Pejic I, Ajmone-Marsan P, Morgante M, Kozumplick V, Castiglioni P, Taramino G and Motto M (1998) Comparative analysis of genetic similarity among maize inbred lines detected by RFLPs, RAPDs, SSRs, and AFLPs. Theor Appl Genet 97:1248-1255.
• Powell W, Morgante M, Andre C, Hanafey M, Vogel J, Tingey S and Rafalski A (1996) The comparison of RFLP, RAPD, AFLP, and SSR (microsatellite) markers for germplasm analysis. Mol Breed 2:225-238.
• Rafalski AJ (2002) Novel genetic mapping tools in plants: SNPs and LD-based approaches. Plant Sci 162:329-333.
• Rieseberg LH (1998) Genetic mapping as a tool for studying speciation. In: Soltis PS and Doyle JJ (eds) Molecular Systematics of Plants. Chapman & Hall, New York, pp 459-487.
• Rohlf FJ (1993) NTSYS-PC Numerical taxonomy and multivariate analysis system, v. 2.0j. Exeter Publ, Setauket, pp 31.
• Rosales-Serna R, Hernandez-Delgado S, Gonzalez-Paz M, Acosta-Gallegos JA and Mayek-Perez N (2005) Genetic relationships and diversity revealed by AFLP markers in Mexican common bean bred cultivars. Crop Sci 45:1951-1957.
• Saini N, Jain N, Jain S and Jain RK (2004) Assessment of genetic diversity among Basmati and non-Basmati rice varieties using AFLP, ISSR and SSR markers. Euphytica 140:133-146.
• Santalla M, De Ron AM and De La Fuente M (2010) Integration of genome and phenotypic scanning gives evidence of genetic structure in Mesoamerican common bean (Phaseolus vulgaris L.) landraces from the southwest of Europe. Theor Appl Genet 120:1635-1651.
• Saunders JA, Hemeida AA and Mischke S (2001) USDA DNA fingerprinting programme for identification of Theobroma cacao accessions. In: Bekeleetal F (ed) Proceedings of the International Workshop on New Technology in Cocoa Breeding 2000. INGENIC Press, London, pp 108-144.
• Schneider S and Excoffier L (1999) Estimation of demographic parameters from the distribution of pairwise differences when the mutation rates vary among sites: Application to human mitochondrial DNA. Genetics 152:1079-1089.
• Schuelke M (2000) An economic method for the fluorescent labeling of PCR fragments. Nat Biotechnol 18:233-234.
• Svetleva D, Pereira G, Carlier J, Cabrita L, Leitão J and Genchev D (2006) Molecular characterization of Phaseolus vulgaris L. genotypes included in Bulgarian collection by ISSR and AFLP analyses. Sci Hort 109:198-206.
• Talukder ZI, Anderson E, Miklas PN, Blair MW, Osorno J, Dilawari M and Hossain KG (2010) Genetic diversity and selection of genotypes to enhance Zn and Fe content in common bean. Can J Plant Sci 90:49-60.
• Tautz D and Schlötterer C (1994) Simple sequences. Curr Opin Genet Dev 4:834-837.
• Tessier C, David J, This P, Boursiquot JM and Charrier A (1999) Optimization of the choice of molecular markers for varietal identification in Vitis vinifera L. Theor Appl Genet 98:171-177.
• Tivang JG, Nienhuis J and Smith OS (1994) Estimation of sampling variance of molecular marker data using the bootstrap procedure. Theor Appl Genet 89:259-264.
• Thormann CE, Ferreira ME, Camargo LEA, Tivang JG and Osborn TC (1994) Comparison of RFLP and RAPD markers to estimating genetic relationships within and among cruciferous species. Theor Appl Genet 88:973-980.
• Tohme J, Orlando-Gonzalez D, Beebe S and Duque MC (1996) AFLP analysis of gene pools of a wild bean core collection. Crop Sci 36:1375-1384.
• Vos P, Hogers R, Bleeker M, Reijans M, Van De Lee T, Hornes M, Frijters A, Pot J, Peleman J, Kuiper M et al. (1995) AFLP a new technique for DNA fingerprinting. Nucleic Acids Res 23:4407-4414.
• Yu MA and Bliss FA (1978) Seed proteins of common bean. Crop Sci 18:431-437.
• Yu K, Park SJ and Poysa V (1999) Abundance and variation of microsatellite DNA sequences in beans (Phaseolus and Vigna). Genome 42:27-34.
• Zhang X, Blair MW and Wang S (2008) Genetic diversity of Chinese common bean (Phaseolus vulgaris L.) landraces assessed with simple sequence repeat markers. Theor Appl Genet 117:629-640.
• AmiGo, http://amigo.geneontology.org (September 10, 2009).
• Lewis PO and Zaykin D (2000) Genetic data analysis: Computer program for the analysis of allelic data, v. 1.0 (d15). http://alleyn.eeb.uconn.edu/gda/2000 (June 15, 2009).
• NCBI, http://www.ncbi.nlm.nih.gov (September 10, 2009).
  » link
• NTSYS-pc Software, http://pritch.bsd.uchicago.edu (November 20, 2009).
  » link
• R-software: http://www.r-project.org (November 22, 2009).
  » link
• PFAM, http://pfam.sanger.ac.uk (September 10, 2009).
  » link
• TFPGA Software, http:// herb.bio.nau.edu/~miller (November 20, 2009).
  » link

Publication Dates

• Publication in this collection
  28 Jan 2011
• Date of issue
  2011

History

• Received
  25 Jan 2010
• Accepted
  29 July 2010