Gene-based SSR markers for common bean (Phaseolus vulgarisL.) derived from root and leaf tissue ESTs: an integration of the BMc series
© Blair et al; licensee BioMed Central Ltd. 2011
Received: 26 November 2010
Accepted: 22 March 2011
Published: 22 March 2011
Sequencing of cDNA libraries for the development of expressed sequence tags (ESTs) as well as for the discovery of simple sequence repeats (SSRs) has been a common method of developing microsatellites or SSR-based markers. In this research, our objective was to further sequence and develop common bean microsatellites from leaf and root cDNA libraries derived from the Andean gene pool accession G19833 and the Mesoamerican gene pool accession DOR364, mapping parents of a commonly used reference map. The root libraries were made from high and low phosphorus treated plants.
A total of 3,123 EST sequences from leaf and root cDNA libraries were screened and used for direct simple sequence repeat discovery. From these EST sequences we found 184 microsatellites; the majority containing tri-nucleotide motifs, many of which were GC rich (ACC, AGC and AGG in particular). Di-nucleotide motif microsatellites were about half as common as the tri-nucleotide motif microsatellites but most of these were AGn microsatellites with a moderate number of ATn microsatellites in root ESTs followed by few ACn and no GCn microsatellites. Out of the 184 new SSR loci, 120 new microsatellite markers were developed in the BMc (Bean Microsatellites from cDNAs) series and these were evaluated for their capacity to distinguish bean diversity in a germplasm panel of 18 genotypes. We developed a database with images of the microsatellites and their polymorphism information content (PIC), which averaged 0.310 for polymorphic markers.
The present study produced information about microsatellite frequency in root and leaf tissues of two important genotypes for common bean genomics: namely G19833, the Andean genotype selected for whole genome shotgun sequencing from race Peru, and DOR364 a race Mesoamerica subgroup 2 genotype that is a small-red seeded, released variety in Central America. Both race Peru and Mesoamerica subgroup 2 (small red beans) have been understudied in comparison to race Nueva Granada and Mesoamerica subgroup 1 (black beans) both with regards to gene expression and as sources of markers. However, we found few differences between SSR type and frequency between the G19833 leaf and DOR364 root tissue-derived ESTs. Overall, our work adds to the analysis of microsatellite frequency evaluation for common bean and provides a new set of 120 BMc markers which combined with the 248 previously developed BMc markers brings the total in this series to 368 markers. Once we include BMd markers, which are derived from GenBank sequences, the current total of gene-based markers from our laboratory surpasses 500 markers. These markers are basic for studies of the transcriptome of common bean and can form anchor points for genetic mapping studies in the future.
Genic microsatellites are those microsatellites based on simple sequence repeats (SSRs) found within, or closely associated with, gene sequences from a given genome . These SSRs tend to be more conserved and of different motifs than SSRs located in other non-gene containing regions of the genome, which are often referred to as genomic microsatellites simply to distinguish them from genic microsatellites ; although both gene and non-gene derived microsatellites are obviously part of the overall genome. Simple sequence repeats are defined as small stretches of repeated DNA, usually of two to six nucleotides, tandemly repeated and located in a given pattern between segments of non-repeated DNA . In practice, remnant repeats can be found on either side of a stretch of SSR and in some occasions different motifs are combined together or either motif is interrupted . This differentiates microsatellites into compound or simple microsatellites in the first case, and perfect and imperfect microsatellites in the latter case .
Common bean, Phaseolus vulgaris L., is an important food legume, basic to the diet of the poor in tropical regions of the world, and a major source of income for small farmers there. Genic microsatellites have been limited in number for this crop. This is perhaps due to two main reasons: 1) a lack of funding has precluded large scale expressed sequence tag (EST) sequencing or even the sufficient construction of many cDNA libraries for the crop and 2) those ESTs and cDNA libraries that exist have not been extensively screened for gene-based SSRs with the exception of the work of Blair et al.  and Hanai et al. [7, 8]. Yet, common bean is essential for micronutrient nutrition and is adaptable to marginal areas for small-scale farm agriculture despite problems of low phosphorus soils or other abiotic constraints [9, 10] and a range of diseases and pests . Therefore a more complete toolbox of molecular tools for this crop is needed especially in the case of gene-based markers which can be based on SSRs polymorphisms as will be discussed here.
In our efforts to accumulate a larger set of genic SSRs, we previously constructed a leaf based cDNA library from Andean genotype G19833  and used a hybridization approach to discover SSRs of various di-nucleotide and tri-nucleotide motifs and develop microsatellites from this library in the BMc (Bean microsatellites from cDNAs) series . We have also recently developed two additional root based cDNA libraries under high and low phosphorus conditions from the Mesoamerican genotype DOR364, the other parent of the mapping population of Blair et al.  and sequenced ESTs from the libraries to discover new SSRs.
The EST sequencing of these libraries is used in this research as the basis for determining the frequency of SSR sequences in root expressed genes as opposed to leaf expressed genes and for adding to the BMc series of microsatellites through an in silico approach to microsatellite discovery as described by Varshney et al.  for some species of cereals. EST-SSRs are more common in cereals than they are in legumes [14–16].
Apart from our efforts, currently there are approximately 70,000 other EST sequences from common bean including collections from Ramirez et al. , Melotto et al.  and Thibivilliers et al.  along with small groups of GenBank entries and a wish-list of further EST efforts . However most of these libraries have not been screened nor compared for SSR markers. The Melotto et al.  libraries from anthracnose infected common bean leaves which contain together approximately 4,000 unigenes has been screened for microsatellites, yielding a set of 140 EST-based SSRs for Hanai et al. [7, 8], although many of these have been used for genetic mapping rather than for germplasm characterization.
The objective of this research was to evaluate the frequency of microsatellites in sequences from different leaf and root EST libraries made in our laboratory, comparing the types of microsatellites from each source tissue. From there we developed the most promising microsatellite loci as gene-based SSR markers that we added to the BMc series of markers . To validate these BMc markers we compared their ability to detect polymorphism in a standard germplasm panel of 18 mapping parent genotypes, which included Mesoamerican, Andean, wild and cultivated accessions that were useful for determining polymorphism information content of the different groups of markers. A final objective was to determine whether any difference in the ability to uncover polymorphism existed between the newly developed BMc markers found in the random EST sequencing versus BMc markers developed by our previous hybridization-based approach.
cDNA library and EST sequencing
Three cDNA libraries were searched for microsatellite containing sequences. These libraries were based on 1) mRNA from leaf and stem tissues as described in Blair et al.  and Ramirez et al.  for the genotype G19833; 2) a library that was made in the pBS-SKII vector from mRNAs of hydroponically grown DOR364 roots which were produced under low phosphorus (LP) conditions and 3) a final library also made in the pBS-SKII vector from mRNAs of hydroponically grown DOR364 roots but which were from a high phosphorus (HP) treatment.
In total, 1308 ESTs were sequenced from the G19833 leaf and stem tissue library: 540 sequences (GenBank entry BQ481427-BQ481965) from Blair et al.  and 768 sequences (HS089176-HS089943) sequenced for this study. Meanwhile, a total of 1815 ESTs were sequenced from the DOR364 root tissue libraries: these being 862 from the HP library (GenBank entries, HS103978-HS104836) and 953 from the LP library (GenBank entries, HS103028-HS103977). Clones from all cDNA libraries were sequenced from the 5'end using BigDye chemistry (Applied Biosystems by Life Technologies; Carlsbad, CA) and di-deoxy-based Sanger sequencing reactions at the Clemson University Genomics Institute (CUGI). All EST sequences were screened for microsatellites to be assigned to the BMc series as described in Blair et al.  and with the methods given below.
SSR identification, primer design and microsatellite amplification
SSRs were identified by screening the EST collections with SSR Locator  with the default option of 1 to 6 nucleotide repeats. Primers were designed using Primer3  with the following conditions: optimum primer length of 20 nucleotides (nt, minimum 18 nt - maximum 26 nt), optimum melting temperature of 50°C (min. 45°C - max. 55°C), an optimum product size of 125 base-pairs (bp, min. 100 bp - max. 350 bp) and an optimum G/C content of 50% (min. 45%- max. 55%). New markers were submitted as STS entries to GenBank and are listed in the Additional file 1 (Table S1).
PCR reaction conditions for all newly designed BMc markers and for the 248 BMc markers from Blair et al.  are as follows: 30 ng of genomic DNA, 0.16 μM of mixed forward and reverse primers, 1X Buffer (10 mM de Tris-HCl pH 8.2, 50 mM KCl, Triton 0.1%, BSA 1mg/ml), 1.5 mM MgCl2, 0.2 mM dNTPs and 1 U Taq polymerase in 12 μL reaction volumes. Amplification conditions were based on those described in Blair et al. [6, 22] with 35 cycles and 47°C annealing temperature. PCR reaction products were run on PTC-200 thermal cyclers (MJR, Bio-Rad Laboratories; Hercules, CA) and then denatured at 94°C and run on 4% polyacrylamide gels (5M urea, 0.5X TBE) in metal backed Owl T-Rex vertical S3S gel units (Thermo Fisher Scientific Inc; Waltham, MA) at constant 120 W. Silver staining was performed as described in Blair et al. [22, 23].
The set of genotypes used for the polymorphism survey in this study was based on a germplasm panel of 18 genotypes described in Blair et al.  as panel I. Both the DOR364 genotype, a Mesoamerican gene pool advanced line from the International Center for Tropical Agriculture (CIAT), and the G19833 genotype, an Andean gene pool Peruvian landrace in the FAO collection at CIAT were obtained from the gene bank in the Genetic Resources Unit (GRU), and used in a polymorphism survey since these were the sources of the EST libraries we screened for microsatellite loci. Along with these two genotypes the germplasm survey included nine more domesticated Mesoamerican accessions and varieties (G3513, G4825, G11360, G11350, G14519, G21212, BAT477, BAT881 and DOR390), four other domesticated Andean accessions or varieties (G21078, G21657, G21242, Radical Cerinza) and three wild accessions (G19892, G24390 and G24404) representing Andean, Mesoamerican and Colombian wild sub-populations) which were also provided by the GRU.
DNA extraction consisted in a CTAB based mini-prep procedure as described in Afanador et al.  using bulk leaf tissue from four greenhouse grown plants per genotype or line. Since the accessions were from lines separated by seed color and maintained at the gene bank, or from advanced lines from the CIAT collection, we assumed homozygosity for all the germplasm but noted any double banding that could indicate a heterozygote or heterogeneous mixture from the four plants. Although beans are a highly inbreeding species (95 to 99%) some outcrossing occurs occasionally so there can be some within accession or intra-population variation and this would be observable in any lanes containing more than one band, representing more than one allele in seeds of the accession.
Allele sizes were estimated for the survey panel and mapping gels based on comparison with 10 and 25 bp molecular weight ladders that were distributed twice on each silver stained gel. A neighbor-joining (NJ) dendogram was constructed with the proportion of shared alleles coefficient and matrix of alleles and genotypes for the survey panel with the software programs Darwin . Polymorphism information content (PIC) was calculated for each marker with Powermarker .
Comparisons of EST-SSR repeat types and marker development
Microsatellites, simples sequence repeat (SSR) class and motif type found with in EST collections positive for SSR loci
GenBank entries for ESTs
Blair (this study)
HP root cDNA
Blair (this study)
LP root cDNA
Blair (this study)
Percentage of SSR types across four EST collections
SSR Type/Genotype/Tissue source
G19833 set 1leaf cDNAs
G19833 set 2 leaf cDNAs
DOR364 root HP
DOR364 root LP
Total SSR and Seq.
Di-nucleotide motifs 1
In the effort to develop additional cDNA-derived microsatellites, we added 120 new BMc (bean microsatellites from cDNAs series) markers to the 248 previously developed BMc markers . Among the microsatellites, the first seventeen (BMc1 to BMc17) were developed from leaf cDNAs in the library described in [6, 12] and as shown in the Additional file 1 (Table S1). A second set of leaf cDNA derived microsatellites from our second EST sequencing effort in this library were designated as BMc18 to BMc27. Meanwhile, 47 microsatellite markers (BMc28 to BMc74 plus BMc77 to BMc109 except BMc55 and BMc59) were developed from the HP root library and 46 other microsatellite markers (BMc55, BMc59, BMc75, BMc76 and BMc78 to BMc108 as well as BMc110 to BMc120) were developed from the LP root cDNA libraries. In summary the largest number of new cDNA derived microsatellites were found in the root libraries (93 out of the 120) compared to the leaf library (27 out of the 120).
Summary of the motif and polymorphism characteristics of microsatellites found in BMc markers
BMc marker types
Leaf EST source
Root EST source
Number of SSRs
Percentage of total
Tetra, penta or hexa-nt based
Monomorphic in survey
Polymorphic in survey
Polymorphic in cultivated
Polymorphic in wild only
Genetic diversity detected
As described above, out of the 120 new BMc markers a total of 98 microsatellites amplified well in the survey panel and these were used for polymorphism survey for the germplasm panel and diversity analysis. In this final set of 98 functional markers, 59 (60.2%) were monomorphic and 39 (39.8%) were polymorphic. The average PIC value of the new polymorphic BMc markers was 0.310 and ranged from 0.099 for the least polymorphic markers to 0.657 for the most polymorphic marker (BMc70).
Polymorphism comparison of the di-nucleotide and tri-nucleotide markers showed that they had similar average PIC values (0.131 and 0.125, respectively) when considering both monomorphic and polymorphic microsatellites together. A similar situation was observed when considering only polymorphic microsatellites, where di- and tri-nucleotide based markers again had similar PIC values (0.322 and 0.301, respectively). None of the tetra-, penta- or hexa-nucleotide repeat-based markers was polymorphic.
Polymorphic markers were in similar proportion (38% in each case) for the BMc markers from leaf ESTs (8 out of 21 functioning markers) and for the BMc markers from root ESTs (30 out of 77 functioning markers). Interestingly some polymorphic root-derived BMc markers (BMc30, BMc40, BMc58, BMc60 and BMc70) showed monomorphic background bands suggesting they were members of gene families with different degrees of diversity in different homologs.
A set of five microsatellites (BMc17, BMc36, BMc44, BMc61 and BMc68) was only polymorphic in the wild accessions but not in the cultivated accessions or varieties. These markers had relatively low PIC values of 0.099 to 0.157. From the 368 current BMc markers, including the 248 from the previous study of Blair et al.  and the 120 described here, a total of 209 (56.8%) of the BMc markers yielded monomorphic results while 159 (43.2%) produced polymorphic results in the germplasm survey.
When the results of the phylogeny analysis of the newly developed makers were combined with the previous markers from Blair et al.  an even clearer picture of the associations emerged (Figure 3c). Although all dendrograms showed very highly supported nodes for the separation of the two main gene pools and the two wild accessions; in the combined analysis, we found very high bootstrap values (ranging from 90 to 100%) based on the strength of the total set of markers evaluated.
The major achievements of this research were 1) to evaluate microsatellite frequency in three cDNA libraries from root and leaf tissues with one of the root libraries developed for the abiotic stress of low phosphorus and 2) to create additional genic microsatellite markers based on low-level sequencing of these EST libraries to use in a polymorphism survey both to understand common bean genetic diversity and to understand the differences in various microsatellite types from different sources and their ability to uncover bean diversity. The creation of new genic microsatellites is especially pressing as only about 230 [2, 7, 8, 27, 28] had been reported before we started our work on the design of BMc microsatellite markers.
In total we have now designed 368 genic microsatellites in the BMc series between the efforts of this study and the previous work of Blair et al. ; all BMc markers were designed from cDNA libraries made from different tissues of the mapping population parents used by Blair et al. [12, 23]. In addition, with this study we have created BMc markers from two different genotypes including G19833 and DOR364 and from leaf tissue and root tissues subjected to low or high phosphorus conditions. The advantage of having markers developed from sequences of both genotypes resides is the fact that the Andean G19833 is being used for whole-genome shotgun sequencing and the Mesoamerican DOR364 provides a commercially useful tropical, small red seeded counterpart to the Andean genotype and to black beans which have been better studied in terms of agronomy as well as EST development .
In addition, both marker types from both genotype sources are useful for evaluation in the reference map based on DOR364 × G19833 studied by Blair et al. [2, 23] which is linked both the UC-Davis  and Univ. of Florida  genetic maps. In terms of the practical use of the microsatellites, the PCR amplification strength was similar for SSRs of different motifs and repeat lengths, which may be typical of gene-derived microsatellites and distinct from genomic microsatellites as first suggested by Blair et al. .
In our previous study of cDNA derived microsatellites  we found that uniformly strong PCR products were obtained with the specific primer sets around the SSR loci in cDNA sequences. In comparison, amplification with non-gene based microsatellites is prone to some pitfalls as discussed by Blair et al.  for AT-rich microsatellites and Blair et al.  for hybridization derived genomic microsatellites. Differences between genic and different kinds of genomic microsatellites have been observed for other marker sets as well [7, 32, 33]. Although the SSR and EST sequencing effort from most of these projects has been small it is useful to have added their sequences to GenBank to compare in the future to larger EST collections from Ramirez et al. , Melotto et al.  and Thibivilliers et al.  as well as future genomic sequences for common bean or related species. Furthermore the possible role of microsatellites as promoters or gene expression enhancers especially in root genes where many AGn microsatellites were found could be studied.
In terms of other di-nucleotide motifs, the lack of GC microsatellites has been observed before within the bean genome [6, 31], while AT-rich microsatellites were not expected to be found in genic sequences neither as di-nucleotides nor tri-nucleotides such as those studied by Blair et al. . There were only a few ACn based microsatellites which was surprising given that enrichment for this motif has yielded about the same number of markers as enrichment with AGn or GAn based probes [7, 34]. Among the tri-nucleotide motifs it appears that AAG (23), ACC (12), AGC (12), AGG (16) and ATC (12) microsatellites are the most common and this may have to do with their frequency in triplet codon use for amino acid incorporation into polypeptides. Additionally, open reading frames are known to have a higher GC percentage than non-translated regions  which might favor tri-nucleotide motifs such as ACC, AGC and AGG. Compared to the results of Blair et al. [6, 31] the ratio of tri-nucleotide to di-nucleotide motifs was fairly high (99 versus 57 in total). Perhaps this was due to a majority being located in the open reading frame rather than in untranslated regions of the original mRNA transcripts represented by the cDNA sequences.
In the second step of this study, we analyzed the potential of two different groups of BMc markers, one from cDNA clone sequencing (120 BMc markers) and one from cDNA hybridization with SSR motifs (248 BMc markers developed from 497 positive cDNA clones) to be used in phylogeny analysis. The full group of markers, therefore, included a total of 368 BMc microsatellites all evaluated against the same germplasm survey from Blair et al. . In that evaluation, genetic diversity was reliably predicted by both types of cDNA based BMc microsatellites. Both sets of markers were useful in separating the Andean and Mesoamerican genepool and accurately placing the wild accessions within each genepool. Two wild accessions (Colombian and Mexican) were separated from the cultivated accessions. Similar results were found with the same diversity panel in Blair et al. .
In summary, cDNA derived markers seem to be very useful for diversity analysis due to the fact that they are derived from genic sequences that are conserved and are highly transferable between different accessions of beans. They were critical in recent studies of diversity in both dry and snap bean cultivars of Phaseolus vulgaris [36, 37]. Therefore, in the future we plan to analyze the frequency of gene-based microsatellites in larger collections of ESTs such as those of Ramirez et al.  or Thibivilliers et al.  which surpass the numbers of ESTs evaluated in the libraries we used here. It will be interesting to see if SSR frequency is similar or different for the multiple libraries used by the first of these authors or the larger set of ESTs from a single rust-infected leaf library evaluated by the second research group. One lesson from this microsatellite evaluation is that it is important to test new markers for consistent patterns of genetic diversity detection. We also plan to test the gene-derived markers in related Phaseolus species.
In terms of the evaluation of genetic diversity we found that genic microsatellites from both EST sequencing and hybridization based approaches performed equally well in distinguishing Andean and Mesoamerican genepools and the Argentinean, Colombian and Mexican wild beans as separate accessions. Therefore, these markers can be used for diversity analysis and for breeding especially in crosses between wild and cultivated beans or between genepools. We expect that next generation sequencing will make the discovery of new transcriptome-based SSRs even easier than the two approaches used so far. Nonetheless, the utility of cDNA derived microsatellites for diversity analysis is well established and is perhaps best explained due to their conservation and slower rate of evolution than genomic microsatellites. In summary, gene-based or 'genic' microsatellites appear to be especially useful for genetic analysis of common bean and it would be ideal to have a larger set of these markers for functional diversity analysis and perhaps association mapping once they are genetically mapped which will be the subject of a separate manuscript to define the regions of the genome that are part of the transcriptome. Finally, these gene-based markers may be the keys to selection of specific traits as they represent expressed genes some of which are likely to have multiple functional alleles with diverse phenotypes as a result. Simple sequence repeats in promoter regions have sometimes been found to be important in controlling gene expression and this may be the case for some of the genic markers discovered in this study as well.
We are grateful to Agobardo Hoyos for germplasm curation and development. We also wish to thank the staff of CUGI that made the sequencing possible including Christopher Saski, Diane Cohen, Michael Atkins and Michael Palmer. Joe Tohme in CIAT and Dorrie Main in CUGI are acknowledged for advice. The funding from USAID-SLO linkage grants is gratefully recognized.
- Varshney RK, Graner A, Sorrells ME: Genic microsatellite markers in plants: features and applications. Trends Biotech. 2005, 23: 48-55. 10.1016/j.tibtech.2004.11.005.View ArticleGoogle Scholar
- Blair MW, Pedraza F, Buendia HF, Gaitán-Solís E, Beebe SE, Gepts P, Tohme J: Development of a genome-wide anchored microsatellite map for common bean (Phaseolus vulgaris L.). Theor Appl Genet. 2003, 107: 1362-1374. 10.1007/s00122-003-1398-6.PubMedView ArticleGoogle Scholar
- Hancock JM: Microsatellites and other simple sequences: genomic context and mutation mechanisms. Edited by: Goldstein DB and Schlotterer C. Microsatellites: Evolution and Applications. Oxford Univ. Press, New York; 1999, 1-9.Google Scholar
- Amos W: A comparative approach to the study of microsatellite evolution. Edited by: Goldstein DB and Schlotterer C. Microsatellites: Evolution and Applications. Oxford Univ. Press, New York 1999, 66-79.Google Scholar
- Ellegren H: Microsatellites: Simple sequences with complex evolution. Nature. 2004, 5: 435-445.Google Scholar
- Blair MW, Muñoz-Torres M, Giraldo MC, Pedraza F: Development and diversity assessment of Andean-derived, gene-based microsatellites for common bean (Phaseolus vulgaris L.). BMC Plant Bio. 2009, 9: 100-10.1186/1471-2229-9-100.View ArticleGoogle Scholar
- Hanai LR, de Campos T, Camargo LEA, Benchimol LL, de Souza AP, Melotto M, Carbonell SAM, Chioratto AF, Consoli L, Formighieri EF, Siquiera MF, Tsai SM, Vieira MLC: Development, characterization and comparative analysis of polymorphism at common bean SSR loci isolated from genic and genomic sources. Genome. 2007, 50: 266-277. 10.1139/G07-007.PubMedView ArticleGoogle Scholar
- Hanai LR, Santini L, Aranha LEC, Pelegrinelli MHF, Gepts P, Tsai SM, Carneiro ML: Extension of the core map of common bean with EST-SSR, RGA, AFLP, and putative functional markers. Mol Breeding. 2010, 25: 25-45. 10.1007/s11032-009-9306-7.View ArticleGoogle Scholar
- Broughton WJ, Hernández G, Blair MW, Beebe SE, Gepts P, Vanderleyden J: Beans (Phaseolus spp.) - Model Food Legumes. Plant Soil. 2003, 252: 55-128. 10.1023/A:1024146710611.View ArticleGoogle Scholar
- Rao IM: Role of physiology in improving crop adaptation to abiotic stresses in the tropics: The case of common bean and tropical forages. Edited by: Handbook of plant and crop physiology (Pessarakli M,). 2002, Marcel Dekker Inc, New York, USA, 583-613.Google Scholar
- Miklas PN, Kelly JD, Beebe SE, Blair MW: Common bean breeding for resistance against biotic and abiotic stresses: from classical to MAS breeding. Euphytica. 2006, 147: 105-131. 10.1007/s10681-006-4600-5.View ArticleGoogle Scholar
- Blair MW, Munoz MC, Pedraza F, Gaitan E, Tohme J, Main D, Frisch D, Wing R: Generation of expressed sequence tags (ESTs) from vegetative tissues of a common bean (Phaseolus vulgaris) mapping parent, G19833. GenBank. 2002, BQ481427-965.Google Scholar
- Varshney RK, Thiel T, Stein N, Landrige P, Graner A: In silico analysis on frequency and distribution of microsatellites in ESTs of some cereal species. Cell Mol Biol Lett. 2002, 7: 537-546.PubMedGoogle Scholar
- Gao L, Tang J, Li H, Jia J: Analysis of microsatellites in major crops assessed by computational and experimental approaches. Molecular Breeding. 2003, 12: 245-261. 10.1023/A:1026346121217.View ArticleGoogle Scholar
- Choumane W, Winter P, Baum M, Kahl G: Conservation of microsatellite flanking sequences in different taxa of Leguminosae. Euphytica. 2004, 138: 239-245. 10.1023/B:EUPH.0000047091.69486.05.View ArticleGoogle Scholar
- Kumpatla SP, Mukhopadhyay S: Mining and survey of simple sequence repeats in expressed sequence tags of dicotyledonous species. Genome. 2005, 48: 985-998. 10.1139/g05-060.PubMedView ArticleGoogle Scholar
- Ramírez M, Graham MA, Blanco-López L, Silvente S, Medrano-Soto A, Blair MW, Hernández G, Vance CP, Lara M: Sequencing and analysis of common bean ESTs: Building a foundation for functional genomics. Plant Physiol. 2005, 137: 1211-1227.PubMedPubMed CentralView ArticleGoogle Scholar
- Melotto M, Monteiro-Vitorello CB, Bruschi AG, Camargo LEA: Comparative bioinformatic analysis of genes expressed in common bean (Phaseolus vulgaris) seedlings. Genome. 2005, 48: 562-570. 10.1139/g05-010.PubMedView ArticleGoogle Scholar
- Thibivilliers S, Joshi T, Campbell KB, Scheffler B, Xu D, Cooper B, Nguyen HT, Stacey G: Generation of Phaseolus vulgaris ESTs and investigation of their regulation upon Uromyces appendiculatus infection. BMC Plant Bio. 2009, 9: 46-10.1186/1471-2229-9-46.View ArticleGoogle Scholar
- da Maia L, Palmieri D, Queiroz V, Marini M, Félix FA, Costa A: SSRLocator: Tool for simple sequence repeat discovery integrated with primer design and PCR simulation. Int J Plant Genomics. 2008, 1-9.Google Scholar
- Rozen S, Skaletsky HJ: Primer3 on the WWW for general users and for biologist programmers. Edited by: Krawetz, S., & Misener, S. Bioinformatics Methods and Protocols: Methods in Molecular Biology, New Jersey, U.S.A.: Humana Press, Ottawa, CA; 2000.Google Scholar
- Blair MW, Giraldo MC, Buendia HF, Tovar E, Duque MC, Beebe S: Microsatellite marker diversity in common bean (Phaseolus vulgaris L.). Theor Appl Genet. 2006, 113: 100-109. 10.1007/s00122-006-0276-4.PubMedView ArticleGoogle Scholar
- Blair MW, Buendia HF, Giraldo MC, Metais I, Peltier D: Characterization of AT-rich microsatellites in common bean. (Phaseolus vulgaris L) Theor Appl Genet. 2008, 118: 91-103. 10.1007/s00122-008-0879-z.PubMedView ArticleGoogle Scholar
- Afanador L, Haley S, Kelly JD: Adoption of a "mini-prep" DNA extraction method for RAPD's marker analysis in common bean. Phaseolus vulgaris Bean Imp Coop. 1993, 36: 10-11.Google Scholar
- Perrier X, Flori A, Bonnot F: Data analysis methods. Edited by: Hamon, P., Seguin, M., Perrier, X., Glaszmann, J. C. Genetic diversity of cultivated tropical plants. Enfield, Science Publishers. Montpellier; 2003,43-76.Google Scholar
- Liu K, Muse SV: PowerMarker: an integrated analysis environment for genetic markers analysis. Bioinformatics. 2005, 21: 22128-2129.Google Scholar
- Yu K, Park SJ, Poysa V: Abundance and variation of microsatellite DNA sequences in beans (Phaseolus and Vigna). Genome. 1999, 42: 27-34. 10.1139/gen-42-1-27.View ArticleGoogle Scholar
- Yu K, Park SJ, Poysa V, Gepts P: Integration of simple sequence repeat (SSR) markers into a molecular linkage map of common bean (Phaseolus vulgaris L.). J Hered. 2000, 91: 429-434. 10.1093/jhered/91.6.429.PubMedView ArticleGoogle Scholar
- Freyre R, Skroch PW, Geffory V, Adam-Blondon AF, Shirmohamadali A, Johnson WC, Llaca V, Nodari RO, Periera PA, Tsai SM, Tohme J, Dron M, Nienhuis J, Vallejos CE, Gepts P: Towards an integrated linkage map of common bean. 4 Development of a core linkage map and alignment of RFLP maps. Theor Appl Genet. 1998, 97: 847-856. 10.1007/s001220050964.View ArticleGoogle Scholar
- Vallejos CE, Sakiyama NE, Chase CD: A molecular marker based linkage map of Phaseolus vulgaris L. Genetics. 1992, 131: 733-740.PubMedPubMed CentralGoogle Scholar
- Blair MW, Muñoz M, Pedraza F, Giraldo MC, Buendía HF, Hurtado N: Development of microsatellite markers for common bean (Phaseolus vulgaris L.) based on screening of non-enriched, small-insert genomic libraries. Genome. 2009, 52: 772-782. 10.1139/G09-053.PubMedView ArticleGoogle Scholar
- Benchimol LL, de Campos T, Carbonell SAM, Colombo CA, Chioratto AF, Formighieri EF, Gouvêa LRL, de Souza AP: 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. 2007, 54: 1747-1762. 10.1007/s10722-006-9184-3.View ArticleGoogle Scholar
- Campos T, Benchimol LL, Carbonell SAM, Chioratto AF, Formighieri EF, de Souza AP: Microsatellites for genetic studies and breeding programs in common bean. Pes Agropec Bras. 2007, 42: 589-592.View ArticleGoogle Scholar
- Gaitán-Solís E, Duque MC, Edwards KJ, Tohme J: Microsatellite Repeats in Common Bean (Phaseolus vulgaris): Isolation, Characterization, and Cross-Species Amplification in Phaseolus ssp. Crop Sci. 2002, 42: 2128-2136.View ArticleGoogle Scholar
- Li YC, Korol AB, Fahima T, Nevo E: Microsatellites within genes: structure, function, and evolution. Mol Bio Evol. 2004, 21: 991-1007. 10.1093/molbev/msh073.View ArticleGoogle Scholar
- Blair MW, Gonzales LF, Kimani P, Butare L: Inter-genepool introgression, genetic diversity and nutritional quality of common bean (Phaseolus vulgaris L.) landraces from Central Africa. Theor Appl Genet. 2010, 121: 237-248. 10.1007/s00122-010-1305-x.PubMedPubMed CentralView ArticleGoogle Scholar
- Blair MW, Chaves A, Tofiño A, Calderón JF, Palacio JD: Extensive diversity and inter-genepool exchange of phaseolin alleles found in world-wide snap bean germplasm analyzed with AFLP and microsatellite markers. Theor Appl Genet. 2010, 120: 1381-1391. 10.1007/s00122-010-1262-4.PubMedView ArticleGoogle Scholar