Next Article in Journal
Synthesis and Antimicrobial Activity of Some Novel Cross-Linked Chitosan Hydrogels
Previous Article in Journal
Commercial Bee Pollen with Different Geographical Origins: A Comprehensive Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Short Note

Isolation and Characterization of 13 New Polymorphic Microsatellite Markers in the Phaseolus vulgaris L. (Common Bean) Genome

1
Hubei Provincial Engineering Research Center of Vegetable Seeds and Seedlings, Wuhan Vegetable Research Institute, Wuhan Academy of Agricultural Science & Technology, Wuhan 430345, China
2
State Key Laboratory of Hybrid Rice, Department of Genetics, College of Life Sciences, Wuhan University, Wuhan 430072, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2012, 13(9), 11188-11193; https://doi.org/10.3390/ijms130911188
Submission received: 6 August 2012 / Revised: 3 September 2012 / Accepted: 3 September 2012 / Published: 7 September 2012
(This article belongs to the Section Molecular Pathology, Diagnostics, and Therapeutics)

Abstract

:
In this study, 13 polymorphic microsatellite markers were isolated from the Phaseolus vulgaris L. (common bean) by using the Fast Isolation by AFLP of Sequence COntaining Repeats (FIASCO) protocol. These markers revealed two to seven alleles, with an average of 3.64 alleles per locus. The polymorphic information content (PIC) values ranged from 0.055 to 0.721 over 13 loci, with a mean value of 0.492, and 7 loci having PIC greater than 0.5. The expected heterozygosity (HE) and observed heterozygosity (HO) levels ranged from 0.057 to 0.814 and from 0.026 to 0.531, respectively. Cross-species amplification of the 13 prime pairs was performed in its related specie of Vigna unguiculata L. Seven out of all these markers showed cross-species transferability. These markers will be useful for future genetic diversity and population genetics studies for this agricultural specie and its related species.

1. Introduction

Phaseolus vulgaris L. (common bean, 2X = 2n = 22) belonging to the Fabacaea family, is a self-pollinated vegetable and the most important legume for human consumption in the world, widely cultivated in North America, South America and Africa. Therefore, it is important to study the genetic diversity, population ecology, evolution and conservation of this species using molecular tools, such as microsatellite markers.
Until now, EST-SSR markers [1,2], BAC-derived microsatellite [3] and microsatellite markers based on non-enriched, small-insert genomic libraries [4] of the common bean have been developed to study the genetic diversity of the common bean population. However, EST-derived SSRs (Simple Sequence Repeats) rely on sets of expressed sequence tags (ESTs) in the public databases and can only be obtained from sequences that are expressed. In addition, BAC-derived microsatellite and microsatellite markers based on non-enriched, small-insert genomic libraries are dependent on the libraries in the public databases. In contrast, the Fast Isolation by AFLP of Sequences Containing repeats (FIASCO) method can quickly identify SSRs throughout the entire genome without knowing the genomic information of a given species. In addition, genomic SSR markers are more polymorphic than EST-SSRs. In this paper, we developed 13 microsatellites for P. vulgaris using the FIASCO method.

2. Results and Discussion

A total of 120 positive colonies with the fragment of expected size were sequenced. Within them, 55 unique microsatellites with sufficient flanking regions were chosen to design prime pairs. All of the 55 primer pairs were surveyed in 38 individuals of P. vulgaris provided by the Wuhan Vegetable Research Institute, Wuhan Academy of Agricultural Science & Technology, and collected from different areas in China. Thirteen primer pairs were identified as polymorphic, with a variable number of alleles per locus. In this study, we developed 13 effectives primers from 120 positive colonies and the effectives primer rate was 10.8% (13/120), which is a higher primer rate than the results recently reported by Blair et al. (2011) of 5.9% (184/3123) [2] and by Garcia et al. (2011) of 3.2% (302/9583) derived from EST-SSR [5], respectively. The effectives primer rate in our results is also higher than the results reported by Córdoba et al. (2010), in which the effectives primer rate is only 0.64% (623/98017) from the BAC-derived SSRs [3], and the results reported by Blair et al. (2009) based on screening of non-enriched, small-insert genomic libraries (microsatellite frequency overall appears to be less than 1% of clones) [4]. It is likely that the effective primer rate differences between various results reported by different researchers are due to the source sequences differences (coding or non-coding regions of the genome) used, or the different methods of generating the SSR library, or the different methods of identification the microsatellites, or the specific species studied. In this study, we use the FIASCO method to identify SSRs throughout the entire genome, apparently, the effectives primer rate was higher than that of other methods, this shows that our method used in this research is superior to other methods. In our results, the average number of alleles per locus was 3.64, and the allele number ranged from two to seven (Table 1). The polymorphic information content (PIC) values ranged from 0.055 to 0.721 over the 13 loci and 7 loci had a PIC greater than 0.5, with a mean value of 0.492, which was higher than that of the EST-derived SSR markers [2,5]. This is in agreement with the results of Garcia et al. [5]. The expected heterozygosity (HE) and observed heterozygosity (HO) levels ranged from 0.057 to 0.814 and from 0.026 to 0.531, with a mean values of 0.453 and 0.165 respectively (Table 1), which was similar to a previous study revealed by EST-SSRs [6].
The exact tests for Hardy–Weinberg equilibrium (HWE) revealed that 11 of the 13 microsatellites showed significant deviation from HWE (p < 0.05) [7] (Table 1). The reason for this may be that P. vulgaris is a predominantly self-pollinated species, which results in a low rate of random mating of different common bean species. Furthermore, the polymorphic estimating of these primers developed in a cultivation of small groups of 38 samples. No significant linkage disequilibrium (LD) was found between these polymorphic loci following the Bonferroni correction [8], which demonstrated that the pairs of loci developed in our study were distributed on different chromosomes. In addition, cross-species amplification of the 13 primer pairs was performed in the related species Vigna unguiculata L. The frequencies of null alleles for those loci that showed significant deviation from HWE (p < 0.05) were tested using the micro-checker software [7]. Of these primers, seven (C136, C130, C115, C132, C106, C45, and G10) could be successfully transferred to the tested species (Table 1) and the rate of transferability was 53.85%, but Garcia et al. (2011) reported a higher rate of transferability of 93.8% of the EST-SSRs [5]. The higher transferability rate of EST-SSRs can be explained by their correspondence to the transcribed component of a gene unit, which confers a high potential for inter-specific transferability [9]. After sequencing confirmation of the SSR-containing regions and the flanking sequences of the amplified bands, the microsatellite loci of V. unguiculata detected by the seven primer pairs were consistent with those detected in P. vulgaris.

3. Experimental Section

The genomic DNA of P. vulgaris was extracted from seeds using the CTAB (cetyltrimethyl ammonium bromide) method [10] with slight modifications. An (AG)n-enriched partial genomic library was then constructed following the FIASCO protocol with certain modifications. Briefly, total genomic DNA was digested with the restriction enzyme MseI (New England BioLabs). After ligating MseI adapters to the digested DNA, amplification was performed using an MseI primer. The amplified DNA fragments were denatured and hybridized to a 5′-biotinylated-(AG)n probe. The enriched fragments containing microsatellite repeats were ligated into a pGEM-T vector (Promega, Madison, USA) and transformed into Escherichia coli TOP10 competent cells. The transformants were cultured, and positive colonies were selected and PCR-tested to determine the insert size. A total of 140 positive colonies with a fragment of the expected size were sequenced using a HiSeq 2000 DNA sequencer (Illumina, University of Utah, USA). These sequences were analyzed for repeat motif regions of microsatellites using the software TANDEM REPEATS FINDER [11]. The Primer Premier 5 software [12] was used to design primer pairs for these regions.
The PCR conditions were optimized for each selected primer pair. All of the amplification reactions were performed in a single block of a PX2 ThermalCycler (Thermo Hybaid). The PCR amplifications were performed in a volume of 15 μL containing 20 ng genomic DNA, 2.5 μL 10× PCR buffer with MgCl2, 1 μM each primer, 250 μM each dNTP and 1 U Taq DNA polymerase (Biocolors). The following thermal profile was used for the PCR amplifications: initial denaturation at 95 °C for 5 min, followed by 30 cycles of 94 °C for 40 s, annealing at the optimal temperature (Table 1) for 45 s, and 72 °C for 40 s, and a final 10 min elongation step at 72 °C. The fragments were separated on 8% (w/v) denaturing polyacrylamide sequencing gels and visualized by silver nitrate staining. The size of the amplified fragments was determined by comparison with the pBR322 DNA/MspI DNA marker of 26–622 bp (Tiangen, Beijing, China).
Polymorphisms of all of the primer pairs were then surveyed in 38 individuals. The SSR band for each was treated as a locus, and a matrix for the SSR phenotypes was assembled. The matrix was used for the statistical analyses below. The data were calculated using the ARLEQUIN 3.01 software [13]. Several parameters of genetic diversity, including the HE, HO, number of alleles (Na) per locus, LD, HWE, and PIC, were calculated using Cervus version 3.0.3.

4. Conclusions

The FIASCO method, which quickly identifies SSRs throughout an entire genome, has previously been applied in the sacred lotus [14,15]. In the current study, we first confirmed that this method can also be used in the common bean and that cross-species amplification with the identified markers can be performed in the related species of cowpea. Hence, the polymorphic microsatellite markers developed in this study will be useful in the construction of genetic linkage maps, classification of germplasm and identification in the breeding of P. vulgaris, as well as in the study of the phylogeny, origin and evolution of Leguminosae crops.

Acknowledgements

The authors thank the earmarked fund for Modern Agro-industry Technology Research System for their support (Grant No. CARS-25-G-30) and Technology Program Project of Wuhan City (200820422172).

Reference

  1. Hanai, L.R.; Campos, T.; Camargo, L.E.A.; Benchimol, L.L.; de souza, A.P.; Melotto, M.; Carbonell, S.A.M. Development, characterization, and comparative analysis of polymorphism at common bean SSR loci isolated from genic and genomic sources. Genome 2007, 50, 266–277. [Google Scholar]
  2. Blair, M.W.; Hurtado, N.; Chavarro, C.M.; Munoz-Torres, M.C.; Giraldo, M.C.; Pedraza, F.; Tomkins, J.; Wing, R. Gene-based SSR markers for common bean (Phaseolus vulgaris L.) derived from root and leaf tissue ESTs: An integration of the BMc series. BMC Plant Biol 2011, 11, 50–59. [Google Scholar]
  3. Córdoba, J.M.; Chavarro, C.; Schlueter, J.A.; Jackson, S.A.; Blair, M.W. Integration of physical and genetic maps of common bean through BAC-derived microsatellite markers. BMC Genomics 2010, 11, 436–445. [Google Scholar]
  4. Blair, M.W.; Torres, M.M.; Pedraza, F.; Giraldo, M.C.; Buendía, H.F.; 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. [Google Scholar]
  5. Garcia, R.A.V.; Rangel, P.N.; Brondani, C.; Martins, W.S.; Melo, L.C.; Carneiro, M.S.; Borba, T.C.; Brondani, R.P. The characterization of a new set of EST-derived simple sequence repeat (SSR) markers as a resource for the genetic analysis of Phaseolus vulgaris. BMC Genetics 2011, 12, 41–54. [Google Scholar]
  6. Gill-Langarica, H.R.; Muruaga-Martínez, J.S.; Vargas-Vázquez, M.L.; Rosales-Serna, R.; Mayek-Pérez, N. Genetic diversity analysis of common beans based on molecular markers. Genet. Mol. Biol 2011, 34, 595–605. [Google Scholar]
  7. Van, O.C.; Hutchinson, W.F.; Wills, D.P.M.; Shipley, P. MICROCHECKER: Software for identifying and correcting genotyping errors in microsatellite data. Mol. Ecol. Notes 2004, 4, 535–538. [Google Scholar]
  8. Rice, W.R. Analyzing tables of statistical tests. Evolution 1989, 43, 223–225. [Google Scholar]
  9. Gupta, P.K.; Rustgi, S.; Sharma, S.; Singh, R.; Kumar, N.; Balyan, H.S. Transferable EST-SSR markers for the study of polymorphism and genetic diversity in bread wheat. Mol. Genet. Genomic 2003, 270, 315–323. [Google Scholar]
  10. Luo, Z.; Zhou, G.; Chen, X.; Lu, Q.; Hu, W. Isolation of high quality genomic DNA from plants (In Chinese). Bulletin of Human Medical University 2001, 26, 178–180. [Google Scholar]
  11. Benson, G. TANDEM REPEATS FINDER: A program to analyze DNA sequences. Nucleic Acids Res 1999, 27, 573–580. [Google Scholar]
  12. Rozen, S.; Skaletsky, H.J. Primer3 on the WWW for general users and for biologist programmers. In Bioinformatics methods and protocols; Krawetz, S., Misener, S., Eds.; Humana Press: Totowa, NJ, USA, 2000; pp. 365–386. [Google Scholar]
  13. Excoffier, L.; Laval, G.; Schneider, S. Arlequin (version 3.01): An integrated software package for population genetics data analysis; Computational and Molecular Population Genetics Lab (CMPG), Institute of Zoology, University of Berne, 2005; Berne, Switzerland. [Google Scholar]
  14. Lei, P.; Zhi, W.Q.; Shuang, M.L.; Hong, G.L.; Xin, F.H.; Wei, D.K.; Yi, D. Isolation and characterization of microsatellite markers in the sacred lotus (Nelumbo nucifera Gaertn.). Mol. Ecol. Res 2007, 7, 1054–1056. [Google Scholar]
  15. Kan, H.; Xing, F.H.; Wei, D.K.; Yi, D. Characterization of 11 new microsatellite loci in taro (Colocasia esculenta). Mol. Ecol. Res 2009, 9, 582–584. [Google Scholar]
Table 1. Characteristics of the 13 microsatellite primers developed in P. vulgaris.
Table 1. Characteristics of the 13 microsatellite primers developed in P. vulgaris.
LocusPrimer Sequence (5′3′)Repeat MotifSize Range (bp)HoHeNaDPICTa (°C)GenBank
C9F:ACAGAGACGAGTGCGTGAGAGTTAG(AT)15445–4700.4690.4284*0.47157JQ739888
R:AAAGACAGTTCTAGGAAGAACCGTC
C33F: CTCTTTCTGCTTCCTTTCTACGC(AG)15536–5650.5310.8147*0.72159JQ739889
R:TTCTTCACAGTCAAGGGAGTAGAAG
C42F:TGTCATAAACCTGCTGGTGAATAAC(AG)28287–3420.0260.4733NS0.44460JQ739900
R:CCTATTTCAGAATCACAGCTATGAC
C45aF: CATGAATGATTACGATTCGAGC(AG)20106–1990.1620.1512NS0.28160JQ739901
R:TGAGTGATTACTAGTGGAACCCA
C76F: AGGCGAAGCAAGAGGTTTTATCCAC(CT)20194–2330.0270.243*0.45459JQ739890
R:GAAGAGGCAGAAAGAACTTACAGCG
C106aF:TTGCAGGTAGCAGGTTGT(TCTA)3TCTG383–4310.0290.5855*0.65657JQ739891
R:CAGACAGATAGATAGAGACGG(TCTA)3TCAA(TCTA)5
C115aF: CGTAGTCTCTTTCGTCCTTTTCTGC(CT)17212–2450.1140.5333*0.57557JQ739892
R: TTACATGCCCTTTCCCTCGTTTG
C117F: GTACCTCCTTTTGAGTTTGTAAGG(CT)12150–1730.0280.4063*0.5555JQ739893
R: GTTCGTAAGCCTACTTTCTCA
C119F: CCACCATTGCTCTCAGTGTTA(CT)21251–2920.1390.453*0.45857JQ739894
R: TAGATGTGTGTTTGTGTTCCG
C130aF: CCATTTCAAAGCAAACCCCTT(CT)3CG(CT)2TG (CT)19298–3490.0280.5973*0.59060JQ739896
R: TGACCCGCGATTATTTACCAC
C132aF: CAGTGGTTATTCTGGGGATT(CT)13477–5020.1180.6745*0.68758JQ739897
R: GGTTGTTTATGGCAGTAGCA
C136aF:GTAAAAGTCTCCTTCTACTTTCCCC(CT)22272–3150.0570.3064*0.50360JQ739898
R:CTCTCAAGCATGTTTGGATTGTAGC
G10aF: TCTTCTGTCCATCCCTCCATACT(AG)27220–2730.0880.3273*0.44460JQ739899
R: GATTGGTGGAAATCGACTTGTCT
F: the forward primer sequences; R: the reverse primer sequences; HO: observed heterozygosity; He: expected heterozygosity; Na: number of alleles; D: deviation from Hardy-Weinberg equilibrium; Ta: annealing temperature; NS: not significant;
*: significant deviations from Hardy-Weinberg expectations after Bonferroni correction at α = 0.05.
The superscript a represents a primer pair that can be used for successful amplification in Vigna unguiculata (Linn.).

Share and Cite

MDPI and ACS Style

Wang, A.; Ding, Y.; Hu, Z.; Lin, C.; Wang, S.; Wang, B.; Zhang, H.; Zhou, G. Isolation and Characterization of 13 New Polymorphic Microsatellite Markers in the Phaseolus vulgaris L. (Common Bean) Genome. Int. J. Mol. Sci. 2012, 13, 11188-11193. https://doi.org/10.3390/ijms130911188

AMA Style

Wang A, Ding Y, Hu Z, Lin C, Wang S, Wang B, Zhang H, Zhou G. Isolation and Characterization of 13 New Polymorphic Microsatellite Markers in the Phaseolus vulgaris L. (Common Bean) Genome. International Journal of Molecular Sciences. 2012; 13(9):11188-11193. https://doi.org/10.3390/ijms130911188

Chicago/Turabian Style

Wang, Aihua, Yi Ding, Zhenhua Hu, Chufa Lin, Shuzhen Wang, Bingcai Wang, Hongyuan Zhang, and Guolin Zhou. 2012. "Isolation and Characterization of 13 New Polymorphic Microsatellite Markers in the Phaseolus vulgaris L. (Common Bean) Genome" International Journal of Molecular Sciences 13, no. 9: 11188-11193. https://doi.org/10.3390/ijms130911188

Article Metrics

Back to TopTop