Skip to main content
Download PDF

Mining non-model genomic libraries for microsatellites: BAC versus EST libraries and the generation of allelic richness

BMC Genomics volume 11, Article number: 428 (2010) Cite this article

Abstract

Background

Simple sequence repeats (SSRs) are tandemly repeated sequence motifs common in genomic nucleotide sequence that often harbor significant variation in repeat number. Frequently used as molecular markers, SSRs are increasingly identified via in silico approaches. Two common classes of genomic resources that can be mined are bacterial artificial chromosome (BAC) libraries and expressed sequence tag (EST) libraries.

Results

288 SSR loci were screened in the rapidly radiating Hawaiian swordtail cricket genus Laupala. SSRs were more densely distributed and contained longer repeat structures in BAC library-derived sequence than in EST library-derived sequence, although neither repeat density nor length was exceptionally elevated despite the relatively large genome size of Laupala. A non-random distribution favoring AT-rich SSRs was observed. Allelic diversity of SSRs was positively correlated with repeat length and was generally higher in AT-rich repeat motifs.

Conclusion

The first large-scale survey of Orthopteran SSR allelic diversity is presented. Selection contributes more strongly to the size and density distributions of SSR loci derived from EST library sequence than from BAC library sequence, although all SSRs likely are subject to similar physical and structural constraints, such as slippage of DNA replication machinery, that may generate increased allelic diversity in AT-rich sequence motifs. Although in silico approaches work well for SSR locus identification in both EST and BAC libraries, BAC library sequence and AT-rich repeat motifs are generally superior SSR development resources for most applications.

Background

Microsatellites, or simple sequence repeats (SSRs), are common features of eukaryotic genomes and can be characterized as generally short, repeated nucleotide sequence elements arrayed in tandem and flanked by non-repetitive regions (reviewed in [1]). SSRs often harbor high levels of polymorphism in terms of repeat number and have been developed into one of the most common classes of genetic markers due to their high degree of reproducibility, ubiquity, codominance, and variability among individuals [2–4]. The multi-allelic nature of SSR loci is thought to derive principally from errors occurring due to slipped-strand mispairing during DNA replication [5–8], however, SSRs may also be generated via alternative means, such as retrotransposition events, interhelical junctions forming during chromosome alignment, unequal crossing over, or gene conversion [9].

Frequency of SSR loci is markedly variable across genomes [10–12], although the broader principles of their genomic organization remain poorly understood [13]. Many features that shape genome evolution generally, such as nucleotide composition, may play a large role in the variability of microsatellite density across the genome (e.g. [14, 15]). Generally, length and density of microsatellites increase with genome size [16], but several exceptions to this rule have been observed (e.g. [13, 17]). Genomic regions of different functionality often maintain SSR loci with different properties. For example, SSRs located within coding regions tend to have an excess of trinucleotide repeats relative to other repeat classes and a specific excess of (CAG) n SSR loci [18]. This is generally attributed (1) to the fact that length variant trinucleotide SSRs maintain the appropriate reading frame within the coding region and (2) to the observation that glutamine (CAG) repeats have fewer detrimental effects within a protein than many other repeated amino acids [19]. Further, the striking variation of abundance in exonic, intronic, and intergenic regions suggests that selection might play a role in the genomic distribution of SSR loci (e.g. [20, 21]). The number of allelic length variants associated with an SSR locus typically increases with increasing average repeat number at that locus [22, 23] (but see [24]), however, and allelic diversity is thought to be primarily a consequence of physical parameters and structural properties of the SSR sequence motif [15].

Parallel to the rapid increase in availability of diverse DNA sequence data, highly labor-intensive methods for the generation of SSR genomic markers (e.g. [25, 26]) have been gradually replaced by in silico data mining approaches using genomic sequence databases (e.g. [4, 27, 28], but see [29]). Two such sources are sequence databases derived from expressed sequence tag (EST) and bacterial artificial chromosome (BAC) libraries. Both types of libraries contain sequenced genomic fragments that are effectively randomly distributed throughout the genome; however, they differ in being comprised of actively transcribed components only in the case of EST libraries versus random genomic fragments in the case of BAC libraries. Although transcribed, EST library derived SSR loci still maintain allelic variability comparable to that in non-transcribed genomic DNA [30] and serve as excellent molecular markers for many applications [31].

A number of genomic tools have been brought to bear in investigations of the rapidly radiating Hawaiian cricket genus Laupala, including an EST library with 10.17 Mb nucleotide sequence [32] and a BAC-end sequence library with an additional 1.71 Mb of genomic DNA sequence (unpub. data). Laupala is a unique evolutionary model system, with one of the highest documented rates of speciation known among invertebrates [33, 34]. Members of the genus are often morphologically and ecologically indistinguishable and species can differ by less than 0.1% nucleotide divergence at nuclear loci [35] but extensive divergence is observed in mate-recognition related behavioral characters [36, 37]. Additionally, Laupala appears to have reduced rates of DNA loss and maintains a relatively large genome size (approximately 11X that of Drosophila melanogaster [38]). Consequently, the Laupala genus has more developed genomic resources than most Orthopteran groups, allowing a broader investigation of the structural properties of allelic SSR variation in this important group of insects.

Here, we present a survey of 288 unique SSR loci in two species of Laupala, L. kohalensis and L. paranigra, identified in BAC and EST library sequence databases from L. kohalensis. The overall distribution of SSR loci within these genomic libraries is examined in addition to the distribution of allelic richness across all SSR loci to evaluate (1) whether SSR density or mean repeat number is elevated in a species with a relatively large genome size, (2) whether particular SSR motifs are particularly common in library sequence and (3) whether the same motifs are likely to harbor significant allelic variation. We further evaluate the structural properties of the SSR loci associated with elevated numbers of length variants and assess the efficacy of BAC and EST library nucleotide sequences for the development of informative molecular markers.

Results

Comparison of repeat structure in BAC and EST genomic libraries

SSRs were generally more abundant and comprised of lengthier repeat structures in BAC library sequence than EST library sequence (Table 1). In total, 186 SSR loci were identified in BAC library sequence (1.71 Mb) and 550 in EST library sequence (10.17 Mb); of these, we were able to design primers in flanking sequence for 135 and 435 loci, respectively. Primers could not be designed for 49 and 108 SSR loci from the BAC and EST libraries, respectively, that were located at a terminal sequence end and for an additional 7 SSR loci due to lack of suitable priming sites. SSRs were approximately two to three times more likely to be identified in a BAC versus an EST library sequence (X2 = 144, df = 1, p > 0.0001) and SSRs comprised nearly five-fold higher percentage of all sequenced bases in BAC library sequence (X2 = 6.8, df = 1, p = 0.0091). Further, SSRs identified in BAC library sequence were comprised of significantly more repeats than were those in EST library sequence (Mann-Whitney test p < 0.0001) (See Table 1).

Table 1 SSR genomic sequence distribution summary statistics
Full size table

Distribution of SSRs in genomic library sequence

To test whether SSRs were more likely to involve a particular sequence motif, we developed a posterior probability distribution of each di- and trinucleotide repeat type (e.g. AG, ACT, etc.) for both the EST and BAC genomic libraries. The analysis could not be extended to include tetra- and pentanucleotide SSR motifs because these comprised less than 2.5% of the dataset and thus were too infrequent to provide adequate statistical power (10 tetranucleotide SSRs in BAC library sequence and 16 tetra- and 2 pentanucleotide SSRs in EST library sequence). The posterior probability distribution was compared to the observed distribution of motifs within SSRs using the chi-square goodness of fit test. In each case, SSRs were found to be significantly non-randomly distributed across motif types (Table 2). This effect was magnified in SSRs identified from BAC library sequence and qualitatively different in EST and BAC libraries. In BAC library sequences, SSRs were strongly biased toward AT-rich repeat motifs for both di- and trinucleotide SSRs (Figure 1). While this was also true of trinucleotide SSRs in the EST library, dinucleotide SSRs from EST library sequence were chiefly characterized by a deficit of repeats containing solely purine or solely pyrimidine bases (Figure 1). Generally, biases in SSR repeat distribution were far greater in BAC library sequence than in EST library sequence. Data were additionally evaluated with respect to a uniform distribution of SSRs across all motif types with qualitatively similar results (i.e. significant deviation from uniform distribution; data not shown).

Table 2 Significant departure from expected distribution of SSR motifs
Full size table
Figure 1
figure 1

Nonrandom distribution of SSRs. Observed (black bars) and expected (white bars) counts of repeats of each sequence motif. Data are grouped by genomic library and by SSR motif size: (A) BAC library dinucleotide repeats and trinucleotide repeats, top and bottom, respectively; (B) EST library dinucleotide repeats and trinucleotide repeats, top and bottom, respectively. Observed values are based on posterior probability distribution of SSR motifs. BAC library SSRs show an excess of AT-rich repeat motifs, while EST library SSRs show a deficit of AT and CG motifs.

Full size image

Distribution of SSRs in genomic library sequence was additionally evaluated with respect to placement within open reading frames (ORFs), including all ORFs in any reading frame. SSR ORF placement was characterized as being either within an ORF or outside of an ORF. Comparisons using chi-square goodness of fit tests revealed that SSRs from BAC library sequence were randomly distributed with respect to ORFs, while those from EST library sequence were non-random, particularly among trinucleotide SSRs (Table 3). These deviations from random expectation among EST library SSR ORF distribution are primarily characterized by deficits of tri- and tetranucleotide SSRs and a slight, though non-significant, excess of dinucleotide SSRs within ORFs in EST library sequence (Table 3 and Figure 2).

Table 3 Placement of SSRs in genomic sequence ORFs
Full size table
Figure 2
figure 2

ORF placement of SSRs. Paired columns show observed and expected frequencies of SSRs (left and right columns, respectively) within each ORF category and are further subdivided by data partitioning showing pooled data, dinucleotide SSRs only, trinucleotide SSRs only, and tetranucleotide SSRs only presented left to right within each panel. Panel A shows distributions for SSRs from BAC library sequence; panel B shows distributions for SSRs from EST library sequence. Observed and expected frequencies are generally in agreement within BAC library sequence, but diverge in EST library sequence, primarily among tri- and tetranucleotide repeats.

Full size image

Distribution of allelic richness in genomic library sequence

Of the 288 putative SSR loci screened for allelic variation in L. kohalensis and L. paranigra, 35 failed to amplify entirely and amplification was successful for only a single species in two of the markers screened. Of the remaining 251 loci, 9 loci were scored as having at least one individual with a null allele and were dropped from subsequent analyses. In total, 242 putative SSR loci were assayed for allelic length variation (122 from the EST library sequence and 120 from the BAC library sequence, see Additional Files 1 and 2).

Allelic richness was evaluated only for loci lacking null alleles (i.e. those that were successfully amplified and scored in all eight Laupala individuals screened) and a simple count of allele number was used to estimate allelic richness. Mean allelic richness was nearly identical for total SSRs from BAC and EST library sequences (μ = 2.98 and 2.78, respectively. Mann-Whitney U test p = 0.837). However, distributions of SSR allelic richness within both BAC and EST library sequence is non-random. Allelic richness was unequally distributed across dinucleotide SSR sequence motifs in both libraries (Table 4). Trinucleotide SSRs were pooled by GC content categories (i.e. 0%, 33%, 67%, or 100%) to increase sample size and no significant effects on allelic richness were detected; however, this likely reflects small within-group sample size (< 8 in most cases) (Table 4). In both BAC and EST library sequence, allelic richness was higher among AT-rich than among GC-rich dinucleotide SSR motifs, visible particularly among AT-dinucleotide SSRs in BAC library sequence. (Figure 3). Although AT-rich motifs were generally overrepresented in BAC library sequence, this pattern was not mirrored in SSRs from EST library sequence. Consequently, this pattern of allelic richness is likely independent of sample size (SSR abundance) in the data set.

Table 4 Allelic richness across SSR sequence motifs, ORF positions, and motif sizes
Full size table
Figure 3
figure 3

Allelic richness across dinucleotide SSR sequence motifs. Medians, interquartile ranges, and maximum and minimum values of allelic richness are shown for each SSR sequence motif type, with GC content of SSR motif increasing left to right across each panel. Both BAC library SSRs (panel A) and EST library SSRs (panel B) generally have greater allelic richness in AT/TA SSR motifs.

Full size image

Allelic richness does not appear to be influenced by placement in ORF regions or size of SSR motif in either BAC or EST library sequence (e.g. di- versus trinucleotide repeat motifs) (Table 4), although the ORF screening methods used here may not have sufficient stringency to fully test this hypothesis.

Repeat number is correlated with allelic richness

Allelic richness is significantly positively correlated with the number of repeats present in genomic library sequence SSRs (Spearman correlation coefficient = 0.479 with p < 0.0001 for BAC library SSRs; Spearman correlation coefficient = 0.404 with p < 0.0001 for EST library SSRs). Reference sequence repeat number ranged widely in both genomic sequence libraries. Consequently, data were pooled by reference sequence repeat number within each library to increase sample size within each repeat class (5, 6, 7-8, 9-20, or 21-42 repeats for BAC library sequence; 5, 6, 7, 8, or 9 or more repeats for EST library sequence) and are presented in Figure 4 subdivided by the number of alleles observed at each locus. Allelic richness appears to increase more rapidly with increasing number of repeats present in EST library sequence than in BAC library sequence (Figure 4).

Figure 4
figure 4

Allelic richness is concentrated in high repeat number SSR loci. Distributions of allelic richness within repeat number classes as SSR locus counts. Left column shows EST library SSRs; right column shows BAC library SSRs. Repeat numbers are grouped to increase sample size within each repeat number class.

Full size image

Discussion

A ubiquitous feature in the genomes of diverse organisms, simple sequence repeat (SSR) loci, or microsatellites, are frequently overrepresented in eukaryotic genomes relative to total base composition [39, 40]. Further, SSR loci have been found to vary significantly in terms of repeat length across a broad range of taxonomic scales [13, 14, 20, 38]. In this study, we screened nearly 12 Mb of genomic DNA sequence from the Hawaiian swordtail cricket Laupala kohalensis, consisting of both BAC and EST library databases, for all SSR loci. Of the 736 total SSR loci identified, 288 were screened using four distinct families each from two Laupala species, L. kohalensis and L. paranigra, and characterized according to amplification success and allelic diversity at each locus. In our analysis, we found that identities of SSR loci are far more conserved across these two closely related species (i.e. that primers designed for L. kohalensis are able to amplify equally well in L. paranigra) than might have been expected [41, 42]. Despite this observation, these same SSR loci showed significant bias from neutral expectations, both in terms of overall genomic distribution and SSR characteristics corresponding to high levels of allelic diversity.

Consistent with the observation of SSR frequency being positively correlated with genome size and intergenic space [13, 16], SSR loci were generally more common in BAC-end sequence than in EST sequence, with a roughly two-fold higher density and double the repeat length. This density is within the range of variation encompassed by Drosophila [43, 44] and similar to that of the honeybee [45] although this is nearly an order of magnitude higher than that of most bivalves [13]. SSR loci do appear to be slightly less common in EST sequence, but with a haploid genome size similar to that of many bivalves [13] and an order of magnitude larger than most Drosophila [38], it seems unlikely that the overall density of SSR loci in Laupala is associated with genome size or size of intergenic space in the genome. SSR surveys in the grasshopper Chorthippus biguttulus indicate that the large genome size in this species is associated with elevated repeat length rather than density of SSR loci within the genome [46]; however, this does not appear to be the case in Laupala, where observed repeat lengths and densities were similar to other insect groups [43–45]. Interestingly, SSR distributions deviated from random with respect to ORF placement only among trinucleotide repeats derived from the EST library, in which variation would not result in a frame shift mutation, and were significantly underrepresented in this case (Table 3, Figure 2). Nevertheless, this deficit of SSR loci within ORF regions in EST library sequence when compared to BAC library sequence suggests that selection on actively expressed regions of the genome, particularly those in frame for translation, likely plays a significant role in the distribution of SSR loci in the Laupala genome.

(AC) n repeats are by far the most numerically dominant of SSR loci in humans and most eukaryotes [47–49]. Contrary to this observation, however, dinucleotide SSR motifs in Laupala BAC-end sequences are strongly biased toward an excess of (AT) n repeats and generally toward repeats with depressed GC base content. This general trend is similar to that reported among SSR loci in the silkworm, Bombyx mori [50]; however, the broader generality of the observation is unclear. No parallel compositional bias was observed for SSR loci identified from EST library sequence. This may reflect physical constraints on the structural properties of SSRs (see [15]); however, the contrast with BAC library-derived SSR loci under similar physical constraints suggests a role for selection constraining SSR composition within genic regions. Although SSR loci are well-known to occur within transcribed sequences and "EST-SSRs" are a useful source of molecular markers [30, 51], this study suggests that EST-derived SSRs may be strongly constrained by selective pressures.

Allelic diversity has been shown to be positively correlated with the length, or repeat number, of SSRs in many organisms [22, 52, 53]. Both biased mutation rates and selection acting on allele size have been suggested as mechanistic explanations for this observation [54, 55]. Here, we observed SSR loci having between 1 and 11 alleles total in the eight individuals sampled. Consistent with previous studies, allelic diversity was found to be significantly positively correlated with repeat number in library sequence, regardless of library of origin (BAC or EST). Although SSR loci with large numbers of repeats were not necessarily likely to have a high degree of allelic diversity, loci with high allelic diversity did show a tendency to be drawn from high-repeat number SSR loci. This was particularly true of BAC library-derived SSR loci. The two libraries used in this study are drawn from different genomic samples, EST and BAC, and the view of genomic processes shaping SSR repeat-number evolution consequently is different for each library. Biased mutation rates appear to be the primary factor driving allelic diversity in BAC library-derived SSRs while selection on allele size may play an additional role in EST library-derived SSRs. In this latter case, distributions of allelic richness may be constrained by selection against frame-shifting mutations or excessive length of repeated amino acid elements [18] and such markers may be of minimal use, and potentially misleading, for studies requiring the assumption of selective neutrality.

Mean allelic diversity in Laupala was nearly equal in SSRs derived from both the BAC and EST libraries, but showed strong effects of SSR base composition. SSRs with low GC content from both libraries were significantly more allele-rich than high GC content repeat motifs. This observation is independent of the numerical abundance of the same repeat motifs, although they may share a mechanistic origin. Different equilibrium microsatellite lengths across species have been attributed to species-specific rates of replication slippage [56]. If the Laupala replication machinery is intrinsically more likely to suffer slippage during replication of AT-rich motifs, this may explain both the prevalence of these repeat structures within the genome and their diversity across lineages. Similar patterns of base composition-biased allelic diversity have been reported in B. mori [50]. Previous studies have reported strikingly different levels of polymorphism among repeat motifs (di-, tri-, tetranucleotide repeats, etc.) (e.g. [50]); here, no parallel effect of motif size (i.e. di-, tri-, or tetranucleotide repeat) or ORF position was observed, although this may largely be a factor of limited sample size (number of loci) and low stringency of ORF identification used in this study.

Conclusions

Many applications of SSR loci as molecular markers are contingent on the presence of a number of selectively neutral length variant alleles to distinguish species or lineages [24, 25]. Here, we used a combination of in silico and laboratory-based analyses to evaluate 288 SSR loci in the Hawaiian swordtail cricket genus Laupala, the first such genomic-scale survey in the Orthoptera. Despite a relatively large genome size, SSR loci do not appear to be particularly dense or large in Laupala. SSR loci are significantly more common in BAC library than in EST library sequence and are heterogeneously distributed across all potential base compositions, with a deficit of GC-rich repeat motifs. While many SSR loci can be identified in EST library sequence, the likelihood that selective pressures shape the frequency and diversity of these loci may restrict their utility for certain applications. Allelic diversity of SSR loci is positively correlated with the repeat length found in library sequence and also influenced by repeat base composition. In Laupala, similar physical structural properties of SSRs and the DNA replication machinery likely contribute to the elevated abundance and allelic diversity of AT-rich repeat motifs, suggesting that future screens of Orthopteran molecular markers may benefit by focusing on such SSR motifs. Although allelic diversity profiles are similar in both BAC and EST library-derived SSR loci, the generally higher frequency of SSRs, larger number of repeats within those SSRs, and reduced likelihood of strong selective constraints relative to EST library-derived SSR loci make BAC library sequences far better in silico sources of SSR loci in Laupala and likely other developing model systems as well.

Methods

Sequence data and analysis

14363 EST sequences, comprising 10.17 Mb, from L. kohalensis were downloaded from the DFCI Cricket Gene Index http://compbio.dfci.harvard.edu/tgi/cgi-bin/tgi/gimain.pl?gudb=Cricket. An additional 1.71 Mb of unpublished BAC-end sequence, also from L. kohalensis were evaluated in parallel (library construction and sequencing by Amplicon Express, Pullman, WA). Files were compiled using the program BioEdit (Ibis Biosciences, Carlsbad, CA; available at http://www.mbio.ncsu.edu/BioEdit/BioEdit.html and redundancies were eliminated based on sequence identity. In total, 11.88 Mb of genomic library sequence from L. kohalensis, averaging 703 bp per sequencing read (674 and 708 bp per sequencing read for BAC and EST libraries, respectively), was available for analysis.

Both genomic libraries were screened for SSR motifs using the Msatfinder script implemented in PERL ([57]; available at http://www.genomics.ceh.ac.uk/msatfinder and all SSRs having repeat motifs of two or greater base pairs and greater than five repeating units were identified. Only simple, perfect repeat motifs were considered; compound and imperfect repeat structures were not included in this analysis. Primers flanking each SSR were designed using Primer3 ([58]; available at http://frodo.wi.mit.edu. Primers could not be designed for terminal repeats (i.e. those SSRs falling at the end of a sequencing read) and all such SSRs were subsequently dropped from the analysis. All primers were designed to produce products ranging in size from 200-350 bp with an optimal Tm of 60°C. Maximum self-complimentarity was set to 8.00 (3.00 at the 3' end) and maximum self-priming was set to 12.00. Program default settings were used for all other parameters.

The location of each identified SSR with respect to open reading frames (ORFs) was determined using the GenBank Open Reading Frame Finder http://www.ncbi.nlm.nih.gov/gorf/orfig.cgi and categorized as being either outside of all ORFs or within one or more ORFs. SSRs bridging an in-frame to out-of-frame boundary were discarded from subsequent analysis.

Animals

Four individuals each from two species of the Hawaiian cricket genus Laupala were used to screen SSR loci for length variants. The two species, L. kohalensis and L. paranigra, are both from the Kohalensis species group and show less than 0.1% nucleotide sequence divergence at nuclear loci [31, 33]. A single L. kohalensis male was drawn from each of four lines, collected from two different sampling localities each in the fourth laboratory generation (Pololu 2006 female #1, Pololu 2006 female #11, Kupehau 2006 female #6, Kupehau 2006 female #1). A single L. paranigra male was drawn from each of four lines, all collected from the same sampling locality (Kaiwiki 2006 #7, Kaiwiki 2006 #8, Kaiwiki 2006 female #1, Kaiwiki 2001 LP1), with the first three in the fourth laboratory generation and the final line in the eleventh laboratory generation.

All individuals were sacrificed by decapitation using a sterile razor. Thoracic and abdominal sections were archived in 100% ethanol, while the head was used for DNA extraction using the DNeasy Blood and Tissue Kit (Qiagen Inc., Valencia, CA) according to manufacturer's specifications. DNA concentration in each sample was evaluated using a NanoDrop spectrophotometer (Thermo Scientific, Wilmington, DE) and each was diluted to a working stock of 10 ng•nl-1.

Screening of SSRs

288 identified SSR loci were screened for repeat variation in the eight animals listed above. PCR reactions containing 1X PCR DyNAzyme II buffer, 0.05 U•μl-1 DyNAzyme II (both from Finnzymes Inc., Woburn, MA), 0.2 mM dNTPs (New England Biolabs, Ipswich, MA), 1 μM each forward and reverse primer (Sigma-Genosys, St Louis, MO), and 0.5 ng•μl-1 genomic DNA template were performed in 20 μl volumes. A denaturing step of 2 minutes at 95°C was followed by 35 cycles of 95°C for 30 seconds, 56°C for 30 seconds, and 72°C for one minute and a final 5 minute elongation step at 72°C. Although these conditions did not work for all 288 primer pairs, no attempt was made to optimize conditions for the minority that failed. Primers are presented in Supplementary Table 1.

All PCR products were evaluated for length variants on 4% 19:1 polyacrylamide:bis-acrylamide gels (SequaGel Sequencing System, National Diagnostics, Atlanta, GA). Each gel was pre-run for 15 minutes at 60 W prior to loading samples and 100 bp ladder (New England Biolabs, Ipswich, MA), then run for 3 hours at 55 W. Bands were visualized by silver staining and scored by hand. As precise determination of band size could not be attained, scores reflected number of unique alleles per SSR locus.

Statistical analysis

All statistics were performed using the SPSS11 statistics package (SPSS Inc., Chicago, IL). Mann-Whitney U tests were used to test for difference in repeat numbers present across libraries. Chi-square goodness-of-fit tests were used to compare observed SSR distribution within libraries with respect to (A) even distribution of SSRs across all sequence motifs (B) the posterior probability distribution of repeat motifs, (C) an even distribution of SSRs internal and external to identified ORFs, and (D) the posterior probability distribution of ORF inclusion. Posterior probability distributions were calculated separately for each genomic library, BAC and EST, respectively. The posterior probability distribution for SSR motif composition was calculated based on the total frequency of SSR motif within a complete genomic library. Correlation between SSR reference sequence repeat number and allelic richness was performed using a Spearman rank correlation. Distribution of allelic richness with regard to (A) library sequence repeat number, (B) placement within versus outside of ORF sequence, and (C) size of SSR motif (i.e. dinucleotide repeats versus trinucleotide repeats versus tetranucleotide repeats) was evaluated using Kruskal-Wallace tests. Similarly, the posterior probability distribution for SSR placement within ORF regions was based on the distribution of ORF sequence within those library sequences containing identified SSRs. To eliminate the possibility of a sampling bias influencing estimates of allelic richness, any loci with one or more individual failing to amplify were removed from analyses using allelic richness estimates. Consequently, all estimates of allelic richness were based on an equal number of individuals (4 L. kohalensis and 4 L. paranigra).

Abbreviations

SSR:

Simple Sequence Repeat

ORF:

Open Reading Frame

DNA:

Deoxyribonucleic Acid

EST:

Expressed Sequence Tag

BAC:

Bacterial Artificial Chromosome

PCR:

Polymerase Chain Reaction

bp:

base pairs

Mb:

megabases (=106 base pairs)

References

  1. Goldstein DB, Schlötterer C, (Eds): Microsatellites: Evolution and Application. 1999, Oxford University Press: Oxford

  2. Tautz D, Renz M: Simple sequences are ubiquitous repetitive components of eukaryotic genomes. Nucleic Acids Res. 1984, 12: 4127-4138. 10.1093/nar/12.10.4127.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  3. Schlötterer C: The evolution of molecular markers - just a matter of fashion?. Nat Rev Genet. 2004, 5: 63-69. 10.1038/nrg1249.

    Article  PubMed  Google Scholar 

  4. Sharma PC, Grover A, Kahl G: Mining microsatellites in eukaryotic genomes. Trends Biotech. 2007, 25: 490-498. 10.1016/j.tibtech.2007.07.013.

    Article  CAS  Google Scholar 

  5. Levinson G, Gutman GA: Slipped-strand mispairing: a major mechanism for DNA sequence evolution. Mol Biol Evol. 1987, 4: 203-221.

    CAS  PubMed  Google Scholar 

  6. Wolff RK, Plaeke R, Jeffreys AJ, White R: Unequal crossing over between homologous chromosomes is not the major mechanism involved in the generation of new alleles at VNTR loci. Genomics. 1991, 5: 382-384. 10.1016/0888-7543(89)90076-1.

    Article  Google Scholar 

  7. Innan H, Terauchi R, Miyashita NT: Microsatellite polymorphism in natural populations of wild plant Arabidopsis thaliana. Genetics. 1997, 146: 1441-1452.

    CAS  PubMed Central  PubMed  Google Scholar 

  8. Stephan WM, Kim Y: Persistence of microsatellite arrays in finite populations. Mol Biol Evol. 1998, 15: 1332-1336.

    Article  CAS  PubMed  Google Scholar 

  9. Wilder J, Hollocher H: Mobile elements and the genesis of microsatellites in dipterans. Mol Biol Evol. 2001, 18: 384-392.

    Article  CAS  PubMed  Google Scholar 

  10. Lagercrantz U, Ellegren H, Andersson L: The abundance of various polymorphic microsatellite motifs differs between plants and vertebrates. Nucleic Acids Res. 1993, 21: 1111-1115. 10.1093/nar/21.5.1111.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  11. Morgante M, Hanafey M, Powell W: Microsatellites are preferentially associated with nonrepetitive DNA in plant genomes. Nat Genet. 2002, 30: 194-200. 10.1038/ng822.

    Article  CAS  PubMed  Google Scholar 

  12. Dieringer D, Schlotterer C: Two distinct modes of microsatellite mutation processes: Evidence from the complete genomic sequence of nine species. Genome Res. 2003, 13: 2242-2251. 10.1101/gr.1416703.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  13. Cruz F, Pérez M, Presa P: Distribution and abundance of microsatellites in the genome of bivalves. Gene. 2005, 346: 241-247. 10.1016/j.gene.2004.11.013.

    Article  CAS  PubMed  Google Scholar 

  14. Bachtrog D, Weiss S, Zangerl B, Brem G, Schlötterer C: Distribution of dinucleotide microsatellites in the Drosophila melanogaster genome. Mol Biol Evol. 1999, 16: 602-210.

    Article  CAS  PubMed  Google Scholar 

  15. Bacolla A, Larson JE, Collins JR, Li J, Milosavljevic A, Stenson PD, Cooper DN, Wells RD: Abundance and length of simple repeats in vertebrate genomes are determined by their structural properties. Genome Res. 2008, 18: 1545-1553. 10.1101/gr.078303.108.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  16. Primmer CR, Raudsepp T, Chowdhary BP, Møller AP, Ellegren H: Low frequency of microsatellites in the avian genome. Genome Res. 1997, 7: 471-482.

    CAS  PubMed  Google Scholar 

  17. Elgar G, Clark MS, Meek S, Smith S, Warner S, Edwards YJK, Bouchireb N, Cottage A, Yeo GSH, Umrania Y, Williams G, Brenner S: Generation and analysis of25 Mb of genomic DNA from the pufferfish Fugu rubripes by sequence scanning. Genome Res. 1999, 9: 960-971. 10.1101/gr.9.10.960.

    Article  PubMed Central  PubMed  Google Scholar 

  18. Katti MV, Ranjekar PK, Gupta VS: Differential distribution of simple sequence repeats in eukaryotic genome sequences. Mol Biol Evol. 2001, 18: 1161-1167.

    Article  CAS  PubMed  Google Scholar 

  19. Perutz MF: Glutamine repeats and neurodegenerative diseases: molecular aspects. Trends Biochem Sci. 2001, 24: 58-63. 10.1016/S0968-0004(98)01350-4.

    Article  Google Scholar 

  20. Tóth G, Gáspári Z, Jurka J: Microsatellites in different eukaryotic genomes: survey and analysis. Genome Res. 2000, 10: 967-981. 10.1101/gr.10.7.967.

    Article  PubMed Central  PubMed  Google Scholar 

  21. Dong P, Wei YM, Chen GY, Li W, Wang JR, Nevo E, Zheng YL: EST-SSR diversity correlated with ecological and genetic factors of wild emmer wheat in Israel. Hereditas. 2009, 146: 1-10. 10.1111/j.1601-5223.2009.02098.x.

    Article  PubMed  Google Scholar 

  22. Wierdl M, Dominska M, Thomas DP: Microsatellite instability in yeast: dependence on the length of the microsatellite. Theor Pop Biol. 1997, 53: 272-283.

    Google Scholar 

  23. Gao H, Cai SL, Yan BL, Chen BY, Yu F: Discrepancy variation of dinucleotide microsatellite repeats in eukaryotic genomes. Biol Res. 2009, 42: 365-375. 10.4067/S0716-97602009000300011.

    Article  CAS  PubMed  Google Scholar 

  24. Chakraborty R, Kimmel M, Stivers DN, Davison LJ, Deka R: Relative mutation rates at di-, tri-, and tetranucleotide microsatellite loci. Proc Natl Acad Sci USA. 1997, 94: 1041-1046. 10.1073/pnas.94.3.1041.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  25. Avise J: Molecular markers, natural history, and evolution. 2004, Sinauer Associates Inc. Sunderland, MA, 2

    Google Scholar 

  26. Weising K, Nybom H, Wolff K, Kahl G: DNA fingerprinting in plants: principles, methods, and applications. 2005, CRC Press, Boca Raton, FL, 2

    Chapter  Google Scholar 

  27. Jayashree B, Punna R, Prasad P, Bantte K, Hash CT, Chandara S, Hoisington DA, Varshney RK: A database of simple sequence repeats from cereal and legume expressed sequence tags mined in silico: survey and evaluation. In Silico Biol. 2006, 6: 607-620.

    CAS  PubMed  Google Scholar 

  28. Korpelainen H, Kostamo K, Virtanen V: Microsatellite marker identification using genome screening and restriction-ligation. Biotechniques. 2007, 42: 479-486. 10.2144/000112415.

    Article  CAS  PubMed  Google Scholar 

  29. Merkel A, Gemmell NJ: Detecting microsatellites in genome data: variance in definitions and bioinformatic approaches cause systematic bias. Evol Bioinformatics. 2008, 4: 1-6.

    CAS  Google Scholar 

  30. Li Y, Korol AB, Fahima T, Beiles A, Nevo E: Microsatellites: genomic distribution, putative functions, and mutational mechanisms: a review. Mol Ecol. 2002, 11: 2453-2465. 10.1046/j.1365-294X.2002.01643.x.

    Article  CAS  PubMed  Google Scholar 

  31. Varshney RK, Graner A, Sorrells ME: Genic microsatellite markers in plants: features and applications. Trends Biotechnol. 2005, 23: 48-55. 10.1016/j.tibtech.2004.11.005.

    Article  CAS  PubMed  Google Scholar 

  32. Danley PD, Mullen SP, Liu F, Nene V, Quackenbush J, Shaw KL: A cricket gene index: a genomic resource for studying neurobiology, speciation, and molecular evolution. BMC Genomics. 2007, 8: 109-124. 10.1186/1471-2164-8-109.

    Article  PubMed Central  PubMed  Google Scholar 

  33. Mendelson TC, Shaw KL: Sexual behavior: rapid speciation in an arthropod. Nature. 2005, 433: 375-376. 10.1038/433375a.

    Article  CAS  PubMed  Google Scholar 

  34. Peccoud J, Simon JC, McLaughlin HJ, Moran NA: Post-Pleistocene radiation of the pea aphid complex revealed by rapidly evolving endosymbionts. Proc Natl Acad Sci USA. 106: 16315-16320. 10.1073/pnas.0905129106.

  35. Shaw KL: Conflict between nuclear and mitochondrial DNA phylogenies of a recent species radiation: what mtDNA reveals and conceals about modes of speciation in Hawaiian crickets. Proc Natl Acad Sci USA. 2002, 99: 16122-16127. 10.1073/pnas.242585899.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  36. Otte D: The Crickets of Hawaii: Origin, Systematics, and Evolution. 1994, Academy of Natural Sciences: Philadelphia

    Google Scholar 

  37. deCarvalho T, Shaw KL: Divergence of courtship and mating behaviors among endemic Hawaiian swordtail crickets. Behavior. 2010, 147: 479-504. 10.1163/000579509X12586249027577.

    Article  Google Scholar 

  38. Petrov DA, Sangster TA, Johnston JS, Hartl DL, Shaw KL: Evidence for DNA loss as a determinant of genome size. Science. 2000, 287: 1060-1062. 10.1126/science.287.5455.1060.

    Article  CAS  PubMed  Google Scholar 

  39. Neff BD, Gross MR: Microsatellite evolution in vertebrates: inference from AC dinucleotide repeats. Evolution. 2001, 55: 1717-1733.

    Article  CAS  PubMed  Google Scholar 

  40. Ross CL, Dyer KA, Erez T, Miller SJ, Jaenike J, Markow T: Rapid divergence of microsatellite abundance among species of Drosophila. Mol Biol Evol. 2003, 20: 1143-1157. 10.1093/molbev/msg137.

    Article  CAS  PubMed  Google Scholar 

  41. May B, Krueger CC, Kincaid HL: Genetic variation at microsatellite loci in sturgeon: primer sequence homology in Acipenser and Scaphirhynchus. Can J Fish Aquat Sci. 1997, 54: 1542-1547. 10.1139/cjfas-54-7-1542.

    Article  CAS  Google Scholar 

  42. Temperley ND, Webster LMI, Adam A, Keller LF, Johnson PCD: Cross-species utility of microsatellite markers in Trichostrongyloid nematodes. J Parasitol. 2009, 95: 487-489. 10.1645/GE-1624.1.

    Article  CAS  PubMed  Google Scholar 

  43. Hutter CM, Schug MD, Aquadro CF: Microsatellite variation in Drosophila melanogaster and Drosophila simulans: a reciprocal test of the ascertainment bias hypothesis. Mol Biol Evol. 1998, 15: 1620-1638.

    Article  CAS  PubMed  Google Scholar 

  44. Pasqual M, Schug MD, Aquadro CM: High density of long dinucleotide microsatellites in Drosophila pseudoobscura. Mol Biol Evol. 2000, 17: 1259-1267.

    Article  Google Scholar 

  45. Solignac M, Vautrin D, Baudry E, Mougel F, Loiseau A, Cornuet JM: A microsatellite-based linkage map of the honeybee, Apis mellifera L. Genetics. 2009, 167: 253-262. 10.1534/genetics.167.1.253.

    Article  Google Scholar 

  46. Ustinova J, Achmann R, Cremer S, Mayer F: Long repeats in a huge genome: microsatellite loci in the grasshopper Chorthippus biguttulus. J Mol Evol. 2006, 62: 158-167. 10.1007/s00239-005-0022-6.

    Article  CAS  PubMed  Google Scholar 

  47. International Human Genome Consortium: Initial sequencing and analysis of the human genome. Nature. 2001, 409: 860-921. 10.1038/35057062.

    Article  Google Scholar 

  48. Astolfi P, Bellizzi D, Sgaramella V: Frequency and coverage of trinucleotide repeats in eukaryotes. Gene. 2003, 317: 117-125. 10.1016/S0378-1119(03)00659-0.

    Article  CAS  PubMed  Google Scholar 

  49. Subramanian S, Mishra RK, Sing L: Genome-wide analysis of microsatellite repeats in humans: their abundance and density in specific genomic regions. Genome Biol. 2003, 4: 1311-1319. 10.1186/gb-2003-4-2-r13.

    Article  Google Scholar 

  50. Prasad MD, Muthulakshmi M, Madhu M, Archak S, Mita K, Nagaraju J: Survey and analysis of microsatellites in the silkworm, Bombyx mori: frequency, distribution, mutations, marker potential and their conservation in heterologous species. Genetics. 2005, 169: 197-214. 10.1534/genetics.104.031005.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  51. Choudhary S, Sethy NK, Shokeen B, Bhatia S: Development of chickpea EST-SSR markers and analysis of allelic variation across related species. Theor Appl Genet. 2009, 118: 591-608. 10.1007/s00122-008-0923-z.

    Article  CAS  PubMed  Google Scholar 

  52. Crozier RH, Kaufmann B, Carew ME, Crozier YC: Mutability of microsatellites developed for the ant Camponotus consobrinus. Mol Ecol. 1999, 8: 271-276. 10.1046/j.1365-294X.1999.00565.x.

    Article  CAS  PubMed  Google Scholar 

  53. Schlötterer C, Ritter R, Harr B, Brem G: High mutation rate of a long microsatellite allele in Drosophila melanogaster provides evidence for allele-specific mutation rates. Mol Biol Evol. 1998, 15: 1269-1274.

    Article  PubMed  Google Scholar 

  54. Garza JC, Slatkin M, Freimer NB: Microsatellite allele frequencies in humans and chimpanzees, with implications for constraints on allele size. Mol Biol Evol. 1995, 12: 594-630.

    CAS  PubMed  Google Scholar 

  55. Dermitzakis ET, Clark AG, Batargias C, Magoulas A, Zouros E: Negative covariance suggests mutation bias in a two-locus microsatellite system in fish Sparus aurata. Genetics. 1998, 150: 1567-1575.

    CAS  PubMed Central  PubMed  Google Scholar 

  56. Kruglyak S, Durrett RT, Schug MD, Aquadro CF: Distribution and abundance of microsatellites in the yeast genome can be explained by a balance between slippage events and point mutations. Mol Biol Evol. 2000, 17: 1210-1219.

    Article  CAS  PubMed  Google Scholar 

  57. Thurston MI, Field D: Msatfinder: detection and characterization of microsatellites. 2005, Oxford UK, [http://www.genomics.ceh.ac.uk/msatfinder/]

    Google Scholar 

  58. Rozen S, Skaletsky HJ: Primer3 on the WWW for general users and for biologist programmers. Bioinformatics methods and protocols: methods in molecular biology. Edited by: Krawetz S, Misener S. 2000, Humana Press, Totowa, NY, 365-386.

    Google Scholar 

Download references

Acknowledgements

This work was financed by National Science Foundation grant IOS-0818862 to KLS. The authors also thank D.J. Fergus and S. Lesnick for laboratory assistance and helpful conversations and three anonymous reviewers for insightful comments.

Author information

Authors and Affiliations

  1. Department of Neurobiology and Behavior, Cornell University, Ithaca, NY, 14850, USA

    Christopher K Ellison & Kerry L Shaw

Authors
  1. Christopher K Ellison
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Kerry L Shaw
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Christopher K Ellison.

Additional information

Authors' contributions

CKE designed the experiment, collected all data, performed data analysis and contributed to writing of the manuscript; KLS contributed to interpretation of results and writing of the manuscript. The final version of this manuscript was approved by both authors.

Electronic supplementary material

12864_2010_3022_MOESM1_ESM.PDF

Additional file 1: Laupala SSR primer and repeat data from BAC library sequences. Table showing all SSR loci screened in this study for allele number. Includes: SSR sequence identifier, primers, reference sequence repeat number, and observed allele number for SSRs derived from BAC library sequence. (PDF 39 KB)

12864_2010_3022_MOESM2_ESM.PDF

Additional file 2: Laupala SSR primer and repeat data from EST library sequences. Table showing all SSR loci screened in this study for allele number. Includes: SSR sequence identifier, primers, reference sequence repeat number, and observed allele number for SSRs derived from EST library sequence. (PDF 38 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Ellison, C.K., Shaw, K.L. Mining non-model genomic libraries for microsatellites: BAC versus EST libraries and the generation of allelic richness. BMC Genomics 11, 428 (2010). https://doi.org/10.1186/1471-2164-11-428

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-11-428

Keywords

BMC Genomics

ISSN: 1471-2164

Contact us