Open Journal of Genetics
Vol. 2  No. 4 (2012) , Article ID: 26308 , 6 pages DOI:10.4236/ojgen.2012.24022

Efficiency of microsatellite isolation from orchids via next generation sequencing

Madhav Pandey, Jyotsna Sharma

Department of Plant and Soil Science, Texas Tech University, Lubbock, USA

Email: madhav.mpandey@gmail.com, jyotsna.sharma@ttu.edu

Received 29 June 2012; revised 30 July 2012; accepted 1 September 2012

Keywords: 454 GS-FLX Sequencing; Cypripedium kentuckiense; Microsatellites; Orchids; Pogonia ophioglossoides; Sanger Sequencing

ABSTRACT

Microsatellites, or simple sequence repeats (SSRs), are highly polymorphic, co-dominant genetic markers commonly used for population genetics analyses although de novo development of species specific microsatellites is costand time-intensive. Orchidaceae is one of the most species-rich families of angiosperms with more than 30,000 species estimated. Despite its high species-diversity, microsatellites are available only for a few species and all were developed by only using Sanger sequencing methods. For the first time in orchids, we used 454 GS-FLX sequencing to isolate microsatellites in two species (Cypripedium kentuckiense and Pogonia ophioglossoides), and report preliminary results of the study. From 1/16th plate that was subjected to sequencing, 32,665 reads were generated, from which 15,473 fragments contained at least one SSR. We selected 20,697 SSRs representing di-, tri-, and tetra-nucleotides. While 3,674 microsatellites had flanking regions on both sides, useable primer pairs could be designed for 255 SSRs. The mean numbers of reads, SSRs, and SSR-containing reads useful for primer design estimated for other 15 orchid species using Sanger sequencing method were 166, 78 and 31, respectively. Results demonstrate that the efficiency of microsatellite isolation in orchids is substantially higher with 454 GS-FLX sequencing technique in comparison to the Sanger sequencing methods.

1. INTRODUCTION

Microsatellites, or simple sequence repeats (SSRs) are regions of DNA that contain short tandem repeats (STRs) of 1 to 6 nucleotides. Microsatellites are known to occur in the genomes of living organisms including bacteria, plants, and animals [1]. Although microsatellites are less abundant in coding regions than in non-coding regions, they are randomly distributed throughout the genome [2]. In general, microsatellites are considered neutral markers, but evidence exists for their role in gene expression [3,4]. Variation in SSRs is thought to arise primarily when strand slippage, accelerated by the presence of repeats, leads to stepwise mutation during replication resulting in higher mutation rates in SSR regions than in other regions of the DNA [2,5].

Over the past two decades, microsatellites have become one of the most popular markers for population genetics [6-10]. High intraspecific polymorphisms, codominant nature, and ease of use are some of the advantages [11] although high costand time-investment limit their widespread use, especially for species whose genomes are not yet well sequenced to allow for rapid searches for SSRs within the existing genomic sequence database. In such cases, de novo isolation of SSRs becomes a necessary first step. While several methods are available for microsatellite isolation [12], the most commonly used conventional method involves screening of microsatellites in genomic libraries after hybridizing with artificial probes, cloning of the hybridized DNA fragments, and sequencing of the colony DNA using the Sanger method [13-17]. Despite high costand timeinvestment, this method yields relatively few (100 to a few hundred, at the most) DNA fragments and the percentage of SSR-containing reads is typically quite low (0.04% to 12%; reviewed in Zane et al., [12]).

Recent advances in next generation sequencing have expedited the identification of SSRs in non-model organisms [18,19]. Various next generation sequencing techniques are available such as 454 GS-FLX (Roche Applied Science, Penzburg, Germany), SOLEXA (Illumina, San Diego, CA), SOLiD (Applied Bioystems, Foster City, CA), tSMS (Helicos Biosciences, Cambridge, MA) [18]. Among these, 454 GS-FLX sequencing is most useful for microsatellite isolation because it generates relatively longer reads than other techniques [20,21]. So far, only few studies have reported isolation of microsatellites using next generation sequencing in plants (reviewed in Zalapa et al. [22]).

Orchidaceae is one of the most species-rich families of angiosperms with more than 30,000 species estimated. Despite this, microsatellites isolated through Sanger sequencing technique are available only for few species, (e.g., [23-37]). To our knowledge, no published reports are available on orchid microsatellites isolated by using high-throughput sequencing.

In this study, we report preliminary data of microsatellite isolation from two orchid species, Cypripedium kentuckiense and Pogonia ophioglossoides, using microatellite enrichment followed by 454 GS-FLX sequencing. DNA of the two species was pooled together before the enrichment and sequencing to reduce per species cost of microsatellite isolation. Both are insect-pollinated and each is presumed to be outcrossing. Microsatellites are not available for either species. We also compared the efficiency of microsatellite isolation between this study, and other orchid studies that used Sanger sequencing method for microsatellite isolation.

2. MATERIALS AND METHODS

Total DNA from approximately 100 mg of leaf tissue from single plants of Cypripedium kentuckiense and Pogonia ophioglossoides was extracted by using Qiagen DNEasy Plant Minikit (Qiagen, Valencia, CA), following the manufacturer’s protocol. Quality and quantity of the DNA was determined in NanoDrop ND1000 (Thermo Scientific, Woburn, MA). Five micrograms of contaminant-free DNA (100 ng/µl concentration and 260/280 ratio of ~1.8) from each of the species was combined in a single tube and was digested with the restriction enzyme RsaI (New England Biolabs, Beverley, MA). After digestion, fragments were ligated to double stranded linkers, and they were denatured and hybridized to three artificial biotinylated microsatellite oligonucleotide mixes (mix one = (AG)12, (TG)12, (AAC)6, (AAG)8, (AAT)12, (ACT)12, (ATC)8; mix two = (AAAC)6, (AAAG)6, (AATC)6, (AATG)6, (ACAG)6, (ACCT)6, (ACTC)6, (ACTG)6; mix three = (AAAT)8, (AACT)8, (AAGT)8, (ACAT)8, (AGAT)8. Hybridized fragments were then captured on magnetic streptavidin beads (Dynal, Lake Success, NY). After washing out the unhybridized fragments, the remaining DNA was eluted from the beads and was amplified with the primer SimpleX-6. Sequencing of 1/16th of a plate of the SSR-enriched libraries was carried out on a 454 GSFLX high throughput sequencer by using titanium chemistry at Savannah River Ecology Lab (SREL), University of Georgia. Sequences were subjected to a 3’-end quality trim if any of the terminal 25 bases of the read contained a quality score (phred score) less than 20, or if it contained an ambiguous base. Reads shorter than 150 bp were discarded. Contig sequences were assembled in CAP3 [38] based on the criteria of 98% sequence identity and a minimal overlap of 75 bp. All singlet and contig DNA reads were screened in Msatcommander version 0.8.1 [39] to search for di-nucleotide SSRs with ≥8 repeats and triand tetra-nucleotide SSRs with ≥7 repeats each. Primer pairs were then designed for the SSRs that contained suitable flanking nucleotide sequences of appropriate length on either side by using PRIMER3 [40] based on the default criteria of the program.

To compare the efficiency of microsatellite isolation between this study and other orchid studies, we have selected 15 studies that have used Sanger sequencing methods for microsatellite isolation and have reported at least two of the three basic data, i.e., number of reads, total number of microsatellites, and total number of SSR-containing reads that were useful for primer design.

3. RESULTS

From the 1/16th plate that was subjected to sequencing, 32,665 unique DNA fragments (reads), including 2,982 contigs, were generated and yielded 9.24 Mb of sequence data. Mean G/C content and read length among the 32,665 reads were 39% and 283 bp, respectively. Approximately 66% of all reads were under 300 bp, and only about 1% of the reads had lengths over 500 bp. The longest read length across all DNA fragments was 702 bp (in two individual reads). Of the 32,665 DNA reads, 15,623 contained at least one SSR consisting of ≥7 or 8 repeats. We observed multiple SSRs in 5,305 individual DNA reads (33.95% of 15,623 reads). Overall, 20,697 SSRs with ≥7 or 8 repeats were identified in the combined, partial genomic library of both species. Di-nucleotide SSRs were the most abundant, followed by tetra- (28.62%), and then tri-nucleotide SSRs (12.58%) (Table 1). Approximately 18% (3,674 of 20,697 SSRs) of the SSRs contained useable flanking regions on both sides of each SSR and hence were selected as candidate SSRs for primer development. Of these, 255 SSRs and their flanking regions (6.94% of 3,674 SSRs) met the default primer design criteria in PRIMER3.

Four, six, and fifteen different types of SSR motifs were observed for di-, tri-, and tetra-nucleotides, respectively (Table 1). More than 87% (10,625 of 12,171) of the di-nucleotide SSRs belonged to type AG/CT, followed by AC/GT (1,192), and then TA (353). Three SSR motifs, i.e., ATC/GAT, AGT/ACT and CTT/AAG constituted approximately 97% (2514 of 2603) of the trinucleotide SSR repeats. More than 71% (4,220 of 5,923) of tetra-nucleotide SSR repeats were comprised of three types of repeat motifs, i.e., AGAT/ATCT, ATGT/ ACAT, and CATT/AATG.

Longest repeat length of di-, tri-, and tetra-nucleotide SSRs was 376 bp (188 repeats), 561 bp (187 repeats),

Table 1. Type and number of different SSR repeat motifs of 20,697 SSRs found in 15,623 reads out of 32,665 reads obtained from combined DNA of Cypripedium kentuckiense and Pogonia ophioglossoides.

and 620 bp (155 repeats), respectively. Mean read length of SSRs was 86 bp across the 20,697 SSRs. Among the three repeat types, the average read length was higher for tri-nucleotides (102 bp) than for the di- (86 bp) and tetranucleotides (68 bp). Majority of the SSRs were shorter than 50 repeats (di-nucleotide, 100 bp; tri-nucleotide, 150 bp; and tetra-nucleotide, 200 bp) (Figure 1). More than 13% (1,636 of 12,171 SSRs) of di-nucleotide and 7% (189 of 2603 SSRs) of the tri-nucleotide SSRs had 100 - 200 repeats (di-nucleotide, 200 - 400 bp; tri-nucleotide, 300 - 600 bp), while only one tetra-nucleotide microsatellite belonged to this size class 100 - 200 repeats (400 - 800 bp).

4. DISCUSSION

Frequency of microsatellites in combined genomes of C. kentuckiense and P. ophioglossoides can be considered relatively high because approximately 50% of the DNA reads contained at least one microsatellite. The number of reads generated in this study was about 200 times higher (32,665 reads) than the mean number of reads in the 15 orchid species (166 reads) using the Sanger sequencing method (see Table 2 for data and references). The number of SSRs found in our study was more than 260 folds higher (20,697 SSRs) than the average number of SSRs (78 SSRs) estimated for those 15 orchid species. The maximum and minimum numbers of SSRs found among the 15 studies were 270 and 12, respectively. The number of SSR-containing reads suitable for primer

Figure 1. Number of di-, tri-, and tetra nucleotide SSRs of different repeat lengths obtained from a combined enriched genomic library of two orchid species, Cypripedium kentuckiense and Pogonia ophioglossoides. Only one tetranucleotide SSR was between 151 and 200 repeats (the color is not visible in the pie-chart because of this low number). Minimum number of repeats for dinucleotide SSR motifs was 8, and it was 7 for tri-nucleotides and tetra-nucleotides.

design in those orchid studies ranged from 7 to 63 with the mean number 31, which is four folds lower than in our study. This clearly indicates that the efficiency of microsatellite isolation in orchids using the next generation sequencing is several folds higher than using the Sanger method.

When the results of this study were compared with the results of microsatellite isolation in other plant species, the number of SSRs (20,697 SSRs) found in this study was more than 12% higher than the mean number of

Table 2. Comparison of efficiency of microsatellite isolation using Sanger and next generation sequencing methods in orchids.

SSRs (18,437 SSRs) estimated for 20 studies that used the 454 GSFLX sequencing [22]. Similarly, the number of microsatellite-containing reads that are useful for primer designing in our study is more than 40% higher than in those 20 studies. The number of SSRs found in the combined genomes of C. kentuckiense and P. ophioglossoides was more than 84 folds higher than the mean number of SSRs (245 SSRs) reported for 71 different plant species using Sanger sequencing method [22]. The mean number of SSRs-containing reads (50 reads) in those 71 studies was about five times lower than in our study. Collectively, these observations further confirm that the combined genomes of C. kentuckiense and P. ophioglossoides have relatively high number of microsatellites, and the efficiency of microsatellite isolation using next generation sequencing is substantially higher than using the Sanger sequencing.

Although the number of oligonucleotide probes used for microsatellite enrichment was lower for di-nucleotide SSRs (2 probes) than for the tri-nucleotides (5 probes) and tetra-nucleotides (13 probes), the number of di-nucleotide SSRs in C. kentuckiense and P. ophioglossoides was higher than tri-and tetra-nucleotides (Table 1). These results indicate that the genomes of these species are rich in di-nucleotide repeats. Other studies have also reported higher percentage of di-nucleotide microsatellite repeats than tri-and/or tetra-nucleotide repeats in plants (e.g., [41] (55% di-nucleotides, and 45% tri-nucleotides); [42] (46% di-nucleotides, 34% tri-nucleotides, 8% tetra-nucleotides, and 12% others); [20] (89% di-nucleotides, 9% trinucleotides and 2% tetra-nucleotides). Despite the fact that we did not use TA as a probe for the microsatellite enrichment procedure, we still detected 353 SSRs comprised of this motif. This is not surprising as the most common di-nucleotide repeat motif found in plants is TA [43].

Generally, as the length of an SSR increases polymorphism also increases [44-46]. Although more than 21% (4,436 of 20,697) of SSRs had more than 50 repeats, primers could not be designed for any of them because in those reads either the repeats were prematurely terminated or did not have suitable flanking regions for primer design. This outcome is due to the limitation of 454 sequencing whereby a majority of the SSR-containing DNA fragments are shorter than 400 bp [21]. Percentage of SSR-containing reads that are useful for primer design are typically lower when using next generation sequenceing technique than with Sanger sequencing because the latter technique can generate fragments up to 1,000 bp long. However, shorter read length in high throughput sequencing is easily outweighed by the larger number of SSR-containing sequences.

REFERENCES

  1. Chambers, G.K. and Macavoy, E.S. (2000) Microsatellites: Consensus and controversy. Comparative Biochemistry and Physiology, 126, 455-476. doi:10.1016/S0305-0491(00)00233-9
  2. Hancock, J.M. (1995) The contribution of slippage-like processes to genome evolution. Journal of Molecular Evolution, 41, 1038-1047. doi:10.1007/BF00173185
  3. Gerber, H., Seipel, K., Georgiev, O., Höfferer, M., Hug, M., Rusconi, S. and Schaffner, W. (1994) Transcriptional activation modulated by homopolymeric glutamine and proline stretches. Science, 263, 808-811. doi:10.1126/science.8303297
  4. Kashi, Y., King, D. and Soller, M. (1997) Simple sequence repeats as a source of quantitative genetic variation. Trends in Genetics, 13, 74-78. doi:10.1016/S0168-9525(97)01008-1
  5. Schlötterer, C. and Tautz, D. (1992) Slippage synthesis of simple sequence DNA. Nucleic Acids Research, 20, 211-215. doi:10.1093/nar/20.2.211
  6. Buschiazzo, E. and Gemmell, N.J. (2006) The rise, fall and renaissance of microsatellites in eukaryotic genomes. Bioessays, 28, 1040-1050. doi:10.1002/bies.20470
  7. Selkoe, K.A. and Toonen, R.J. (2006) Microsatellites for ecologists: A practical guide to using and evaluating microsatellite markers. Ecology Letters, 9, 615-629. doi:10.1111/j.1461-0248.2006.00889.x
  8. Mittal, N. and Dubey, A. (2009) Microsatellite markers— A new practice of DNA based markers in molecular genetics. Pharmacognosy Reviews, 3, 235-246.
  9. Jones, A.G., Small, C.M., Paczolt, K.A. and Ratterman, N.L. (2010) A practical guide to methods of parentage analysis. Molecular Ecology Resources, 10, 6-30. doi:10.1111/j.1755-0998.2009.02778.x
  10. Guichoux, E., Lagache, L., Wagner, S., et al. (2011) Current trends in microsatellite genotyping. Molecular Ecology Resources, 11, 591-611. doi:10.1111/j.1755-0998.2011.03014.x
  11. Jarne, P. and Lagoda, P. (1996) Microsatellites, from molecules to populations and back. Trends in Ecology and Evolution, 11, 424-429. doi:10.1016/0169-5347(96)10049-5
  12. Zane, L., Bargelloni, L. and Patarnello, T. (2011) Strategies for microsatellite isolation: A review. Molecular Ecology, 11, 1-16. doi:10.1046/j.0962-1083.2001.01418.x
  13. Rassmann, K., Schlötterer, C. and Tautz, D. (1991) Isolation of simple-sequence loci for use in polymerase chain reaction-based DNA fingerprinting. Electrophoresis, 12, 113-118. doi:10.1002/elps.1150120205
  14. Kijas, J.M.H., Fowler, J.C.S., Garbett, C.A. and Thomas, M.R. (1194) Enrichment of microsatellites from the Citrus genome using biotinylated oligonucleotide sequences bound to streptavidin-coated magnetic particles. Biotechniques, 16, 656-662.
  15. Kandpal, R.P., Kandpal, G. and Weissman, S.M. (1994) Construction of libraries enriched for sequence repeats and jumping clones, and hybridization selection for regions specific markers. Proceedings of National Academy of Science USA, 91, 88-92. doi:10.1073/pnas.91.1.88
  16. Fischer, D. and Bachmann, K. (1998) Microsatellite enrichment in organisms with large genomes (Allium cepa L.). Biotechniques, 24, 796-802.
  17. Pandey, M., Gailing, O., Fischer, D., Hattemer, H.H. and Finkeldey, R. (2004) Characterization of microsatellite markers in sycamore (Acer pseudoplatanus L.). Molecular Ecology Notes, 4, 253-255. doi:10.1111/j.1471-8286.2004.00633.x
  18. Hudson, M.E. (2008) Sequencing breakthroughs for genomic ecology and evolutionary biology. Molecular Ecology Notes, 8, 3-17.
  19. Morozova, O. and Marra, M.A. (2008) Applications of nextgeneration sequencing technologies in functional genomics. Genomics, 92, 255-264. doi:10.1016/j.ygeno.2008.07.001
  20. Buehler, D., Graf, R., Holderegger, R. and Gugerli, F. (2011) Using the 454 pyrosequencing-based technique in the development of nuclear microsatellite loci in the alpine plant Arabis alpine (Brassicaceae). American Journal of Botany, 98, e103-e105. doi:10.3732/ajb.1000488
  21. Jennings, T.N., Knaus, B.J., Mullins, T.D., Haig, S.M. and Cronn, R.C. (2011) Multiplexed microsatellite recovery using massively parallel sequencing. Molecular Ecology Resources, 11, 1060-1067. doi:10.1111/j.1755-0998.2011.03033.x
  22. Zalapa, J.E., Cuevas, H., Zhu, H., Steffan, S., Senalik, D., Zeldin, E., McCown, B., Harbut, R. and Simon, P. (2012) Using next-generation sequencing approaches to isolate simple sequence repeat (SSR) loci in the plant sciences. American Journal of Botany, 99, 193-208. doi:10.3732/ajb.1100394
  23. Swarts, N.D., Sinclair, E.A. and Dixon, K.W. (2007) Characterization of microsatellite loci in the endangered grand spider orchid Caladenia huegelii (Orchidaceae). Molecular Ecology Notes, 7, 1141-1143. doi:10.1111/j.1471-8286.2007.01810.x
  24. Xia, K., Ye, X. and Zhang, M. (2008) Isolation and characterization of nine microsatellite markers for Cymbidium sinense. Horticulture Science, 43, 1925-1926.
  25. Tang, M., Zeng, C.-X., Bi, Y.-F. and Yang, J.-B. (2012) Microsatellite markers for the Chinese endangered and endemic orchid Cymbidium tortisepalum (Orchidaceae). American Journal of Botany, 99, e11-e13. doi:10.3732/ajb.1100307
  26. Pinheiro, F., Santos, M.O., Barros, F., et al. (2008) Isolation and characterization of microsatellite loci in the Brazilian orchid Epidendrum fulgens. Conservation Genetics, 9, 1661-1663. doi:10.1007/s10592-008-9514-4
  27. Gustafson, S. and Thorén, P.A. (2001) Microsatellite loci in Gymnadenia conopsea, the fragrant orchid. Molecular Ecology Notes, 1, 81-82. doi:10.1046/j.1471-8278.2001.00033.x
  28. Campbell, V.V., Rowe, G., Beebee, T.J.C. and Hutchings, M.J. (2002) Isolation and characterization of microsatellite primers for the fragrant orchid Gymnadenia conopsea (L.) R. Brown (Orchidaceae). Conservation Genetics, 3, 209-210. doi:10.1023/A:1015224531795
  29. Lombardo, V.T., Hopkins, S.E., Selosse, M.-A., Cozzolino, S. and Taylor, D.L. (2008) Isolation and characterization of new polymorphic microsatellite loci in the mixotrophic orchid Limodorum abortivum L. Swartz (Orchidaceae). Molecular Ecology Resources, 8, 1117-1120. doi:10.1111/j.1755-0998.2008.02176.x
  30. López-Roberts, M.C., Almeida, P.R.M., Oliveira, E.J.F. and van de Berg, C. (2012) Microsatellite marker development for the threatened orchid Masdevallia solomonii (Orchidaceae). American Journal of Botany, 99, e66-e68. doi:10.3732/ajb.1100364
  31. Soliva, M., Gautschi, B., Salzmann, C., Tenzer, I. and Widmer, A. (2000) Isolation and characterization of microsatellite loci in the orchid Ophrys araneola (Orchidaceae) and a test of cross-species amplification. Molecular Ecology, 9, 2178-2179. doi:10.1046/j.1365-294X.2000.105313.x
  32. Cotrim, H.C., Monteiro, F.A., Sousa, E.S., Fay, M.F., Chase, M.W. and Pais, M.S. (2009) Isolation and characterization of novel polymorphic nuclear microsatellite markers from Ophrys fusca (Orchidaceae) and crossspecies amplification. Conservation Genetics, 9, 739-742. doi:10.1007/s10592-008-9634-x
  33. Li, L.-N., Zeng, S., Zheng, F., Chen, Z.-L., Wu, K.-L., Zhang, J.-X. and Dual, J. (2010) Isolation and characterization of 10 polymorphic microsatellite loci in Paphiopedilum concolor (Batem.) Pfitzer (Orchidaceae) and cross-species amplification. Horticulture Science, 45, 1286-1287.
  34. Rodrigues, K.F. and Kumar S.V. (2009) Isolation and characterization of 24 microsatellite loci in Paphiopedilum rothschildianum, and endangered slipper orchid. Conservation Genetics, 10, 127-130. doi:10.1007/s10592-008-9533-1
  35. Pellegrino, G.D., Cafasso, A., Widmer, M., Soliva, A., Musacchio, A. and Cozzolino, S. (2001) Isolation of microsatellite loci from the orchid Serapias vomeracea (Orchidaceae) and cross-priming to the other Serapias species. Molecular Ecology Notes, 1, 279-280. doi:10.1046/j.1471-8278.2001.00106.x
  36. Phuekvilai, P., Pradit, P. and Surin, P. (2009) Development of microsatellite markers for Vanda orchid. Kasetsart Journal of Natural Science, 43, 497-506.
  37. Bory, S., Da Silva, D., Risterucci, A.M., Grisoni, M., Besse, P. and Duval, M.F. (2008) Development of microsatellite markers in cultivated Vanilla: Polymorphism and transferability to other Vanilla species. Scientia Horticulturae, 115, 420-425. doi:10.1016/j.scienta.2007.10.020
  38. Huang, X. and Madan, A. (1999) CAP3: A DNA sequence assembly program. Genome Research, 9, 868-877. doi:10.1101/gr.9.9.868
  39. Faircloth, B. (2008) Msatcommander: Detection of microsatellite repeat arrays and automated, locus-specific primer design. Molecular Ecology Resources, 8, 92-94. doi:10.1111/j.1471-8286.2007.01884.x
  40. Rozen, S. and Skaletsky, H.J. (2000) Primer3 on the WWW for general users and for biologist programmers. Methods in Molecular Biology, 132, 365-386.
  41. Toth, G., Gaspari, Z. and Jurka, J. (2000) Microsatellites in different eukaryotic genomes: Survey and analysis. Genome Research, 10, 967-981. doi:10.1101/gr.10.7.967
  42. Pan, L., Xia, Q., Quan, Z., Liu, H., Ke, W. and Ding, Y. (2010) Development of novel EST-SSRs from sacred lotus (Nelumbo nucifera Gaertn) and their utilization for the genetic diversity analysis of N. nucifera. Journal of Heredity, 101, 71-82. doi:10.1093/jhered/esp070
  43. Lagercrantz, U., Ellegren, H. and Andersson, L. (1993) The abundance of various polymorphic microsatellite motifs differs between plants and vertebrates. Nucleic Acid Research, 21, 1111-1115. doi:10.1093/nar/21.5.1111
  44. Goldstein, D.B. and Clark, A.G. (1995) Microsatellite variation in North American populations of Drosophila melanogaster. Nucleic Acid Research, 23, 3882-3886. doi:10.1093/nar/23.19.3882
  45. Bachtrog, D., Agis, M., Imhof, M. and Schloetterer, C. (2000) Microsatellite variability differs between dinucleotide repeat motifs-evidence from Drosophila melanogaster. Molecular Biology and Evolution, 170, 1277-1285. doi:10.1093/oxfordjournals.molbev.a026411
  46. Xu, Z., Gutierrez, L., Hitchens, M., Scherer, S., Sater, A.K., et al. (2008) Distribution of polymorphic and nonpolymorphic microsatellite repeats in Xenopus tropicalis. Bioinformatics and Biology Insights, 2, 157-169.