Skip to main content
  • Research article
  • Open access
  • Published:

High frequency of microsatellites in S. cerevisiae meiotic recombination hotspots

Abstract

Background

Microsatellites are highly abundant in eukaryotic genomes but their function and evolution are not yet well understood. Their elevated mutation rate makes them ideal markers of genetic difference, but high levels of unexplained heterogeneity in mutation rates among microsatellites at different genomic locations need to be elucidated in order to improve the power and accuracy of the many types of study that use them as genetic markers. Recombination could contribute to this heterogeneity, since while replication errors are thought to be the predominant mechanism for microsatellite mutation, meiotic recombination is involved in some mutation events. There is also evidence suggesting that microsatellites could function as recombination signals. The yeast S. cerevisiae is a useful model organism with which to further explore the link between microsatellites and recombination, since it is very amenable to genetic study, and meiotic recombination hotspots have been mapped throughout its entire genome.

Results

We examined in detail the relationship between microsatellites and hotspots of meiotic double-strand breaks, the precursors of meiotic recombination, throughout the S. cerevisiae genome. We included all tandem repeats with motif length (repeat period) between one and six base pairs. Long, short and two-copy arrays were considered separately. We found that long, mono-, di- and trinucleotide microsatellites are around twice as frequent in hot than non-hot intergenic regions. The associations are weak or absent for repeats with less than six copies, and also for microsatellites with 4–6 base pair motifs, but high-copy arrays with motif length greater than three are relatively very rare throughout the genome. We present evidence that the association between high-copy, short-motif microsatellites and recombination hotspots is not driven by effects on microsatellite distribution of other factors previously linked to both recombination and microsatellites, including transcription, GC-content and transposable elements.

Conclusion

Our findings suggest that a mutation bias relating to recombination hotspots causing repeats to form and grow, and/or regulation of a subset of hotspots by simple sequences, may be significant processes in yeast. Some previous evidence has cast doubt on both of these possibilities, and as a result they have not been explored on a large scale, but the strength of the association we report suggests that they deserve further experimental testing.

Background

Microsatellites are direct tandem repeats of 1–6 base pair sequence motifs, often strung together in long arrays. They occur much more commonly than expected by chance in the genomes of all eukaryotes [1–3]. The reasons for this are not yet fully understood, but increasing evidence indicates that many microsatellites are functionally important in regulating gene expression [4–10] and possibly also meiotic recombination [11–14]. Microsatellites are also of interest because of their widespread use as genetic markers for applications in genome mapping [15–17], gene hunting [18–20], forensics [21], deducing kinship [22], population genetics [23–25] and the study of the evolution of species [26–28]. These applications depend on assumptions about microsatellite evolution that, at present, are overly simplistic because of unexplained heterogeneity in mutation rates between loci, and an increased understanding of microsatellite evolution and mutational mechanisms is therefore being sought (reviewed in [29, 30]).

Slipped strand mispairing during DNA replication is currently thought to cause most microsatellite mutations [31], but it has also been proposed that unequal meiotic recombination could drive microsatellite evolution [32]. Recombination has been demonstrated to cause instability of some microsatellite loci implicated in human disease (reviewed in [33]), but evidence has counted against it being considered a significant factor in microsatellite evolution. Microsatellite instability was not found to be reduced in recombination deficient strains of E. coli [34] or S. cerevisiae [35] and similar microsatellite mutation rates have been reported for the non-recombining human Y chromosome and the autosomes [36–38]. Also, no association has been found between microsatellite variation and recombination rates on scales of several hundred thousand base pairs in humans [39, 40]. Recent evidence has shown, however, that meiotic recombination events predominantly occur in narrow hotspots of 1–2.5 kilo bases (kb) separated by as much as 50–100 kb of DNA that very seldom recombines [41–44]. Data about the relationship between microsatellites and recombination hotspots at this narrow scale are sparse, and there are some signs that it merits further investigation. A poly-AC array inserted near a recombination hotspot in S. cerevisiae mutated with high frequency [12], and it has recently been found that polymorphic microsatellites are over-represented in human hotspots [45]. There is also some evidence that microsatellites could have a role in regulating hotspot recombination [11–14], increasing the relevance of studying their association with hotspots, since the basis in sequence of the control of hotspot locations is not yet well understood [42–44, 46–48].

It has been shown previously that microsatellite frequencies correlate with broad scale recombination rates in rats, mice and humans [49]. Microsatellites are also associated with intermediate scale recombination rates [50], as well as hotspots in their narrowest known sense [43] in the human genome. So far, however, these studies have reported little detail about the relationship between recombination hotspots and microsatellites. An ideal model organism in which to further examine the association is the yeast S. cerevisiae, since it is the simplest eukaryote, and recombination hotspots have been mapped throughout its entire genome [42]. Factors that could complicate an association between microsatellites and recombination are likely to be less problematic in yeast since, for example, the locations of genes and their expression levels have been well-characterized, making it possible to control for the links between microsatellites, recombination and transcription. Also, transposable or other known repetitive elements are not likely to mediate a link between recombination hotspots and microsatellites in yeast, since these elements are not enriched in yeast hotspots [42], as they are in human hotspots [43].

We investigated in detail the association between microsatellites and hotspots of meiotic double-strand breaks (DSBs), the precursors of meiotic recombination, throughout the S. cerevisiae genome [42]. As well as long microsatellite arrays, we considered low copy number repeats, which have not been studied previously in relation to recombination, including those with only two copies. This allowed us to address the question of whether recombination is involved in the origin of microsatellites, which has previously been considered to occur mainly by accumulation of random point mutations [51]. An association between low-copy microsatellites and hotspots would suggest the involvement of recombination as a mutational mechanism in microsatellite evolution, since replication slippage is expected to act with significant frequency only on arrays of at least six copies [52–54], and there is no available evidence to suggest that short microsatellites have the potential to stimulate recombination.

We found several types of microsatellite to be strongly associated with recombination hotspots in S. cerevisiae, with levels of enrichment greater than two-fold. The associations are, however, stronger for longer microsatellites, and weak or absent for repeats with less than six copies. Our findings suggest that the link between microsatellites and recombination deserves further experimental exploration.

Results

We used hotspot locations mapped by Gerton and co-workers throughout the S. cerevisiae genome using microarray analysis of meiotic DSB frequency [42]. This study identified 177 hotspots, which encompassed all previously known meiotic recombination hotspots in the species, and 40 coldspots. For the purposes of our analysis, we extended the hotspots and coldspots to include the intergenic regions (IGRs) adjacent to the open reading frames (ORFs) identified by Gerton and co-workers [42], since yeast hotspots are typically centred on IGRs, in which most DSBs occur [55]. The hotspots as we defined them have a mean length of 3466 bp. The principal statistical comparisons we made were between hot and non-hot, rather than hot and cold regions, since the cold regions are too few to provide a reliable enough picture of microsatellite density, and recombination frequencies are very low in all experimentally tested regions outside hotspots [41, 44].

In general, numbers of repeats are very much lower in ORFs than IGRs, (Table 1), despite the fact that ORFs cover 73.5% of the genome. This is not surprising, since array length change mutations in microsatellites other than tri- or hexanucleotide repeats would cause frame-shifts in ORFs, destroying gene function. Short (3–5 bp) mononucleotide runs have similar frequency in ORFs and IGRs, but this is likely to be due to coding sequence such as AAA (Lys), GTTTTA (Val Leu), GGG (Gly) or AGGGTT (Arg Val), because the vast majority of the short mononucleotide repeats genome-wide are only three bp long. When making comparisons between hot and non hot regions, we accounted for the low microsatellite abundance in ORFs by comparing ORFs exclusively with other ORFs, and IGRs only with other IGRs. We found the abundance of short-motif, AT-rich repeats to be dramatically higher than other repeat types throughout the genome, so we divided microsatellites by motif length as well as by array length in order not to lose information about longer motifs. We also separated poly-A from poly-G. Nineteen physically independent categories of motif and array length were used in total (see Methods section).

Table 1 Total number of microsatellite repeats and percentage of regions with at least one repeat in the S. cerevisiae genome. The e value denotes the number of bases in any part of a repeat within which no more than one mismatch was allowed with respect to the consensus motif. A lower e value therefore results in the detection of more imperfect repeats.

High microsatellite frequencies in meiotic recombination hotspots

Microsatellite frequencies in meiotic recombination hotspots and non-hot regions of the S. cerevisiae genome can be found in Additional file 1, Tables S1 and S2. Several types of microsatellite have significantly different frequency in hot than non-hot areas (alpha, adjusting for Bonferroni's correction = 0.0026, Table 2). Repeat frequencies in the 40 coldspots are generally lower than in other non-hot regions, but these differences are not statistically significant (Additional file 1, Tables S1 and S2). The correlation between DSB intensity level, assayed for all yeast ORFs by Gerton and co-workers [42], and microsatellite frequency, is generally weak (Additional file 1, Tables S3 and S4), but several repeat types, especially long poly-A and dinucleotide microsatellites, are markedly more abundant in hotspots than non-hot regions (Figure 1, Table 2).

Figure 1
figure 1

Frequencies of high-copy, short-motif repeats in yeast intergenic regions. Mean microsatellite frequencies in S. cerevisiae IGRs divided according to DSB intensity into 473 hot, 89 cold and 5431 other regions, which were all IGRs not categorized as either hot or cold. Poly-AT arrays comprised the majority of dinucleotide repeats and are highlighted in grey. Error bars are plus and minus one SEM.

Table 2 Microsatellite types with a significant difference in frequency either between hot and non-hot IGRs, or hot and non-hot ORFs, in the S. cerevisiae genome. Significance was inferred where p < 0.0026, with the level of alpha adjusted for 19 independent classes of repeat using Bonferroni's correction. The Mann-Whitney U Test or T Test was used, depending whether samples were normally distributed. The e value denotes the number of bases in any part of a repeat within which no more than one mismatch was allowed with respect to the consensus motif. A lower e value therefore results in the detection of more imperfect repeats.

Of the types of microsatellite we investigated, mononucleotide runs are by far the most common, and long arrays are highly over-represented in hotspots. Although poly-A (n ≥ 6) is less than 28% enriched in hot IGRs, and is more common in non-hot than hot ORFs, poly-A (n ≥ 14) is between two and two and a half fold more common in hot IGRs, and poly-G (n ≥ 14) is nearly five fold over-represented, though this figure may be misleading as numbers of poly-G arrays are very low (Table 1). We used a lower limit of 14 bp to define long mononucleotide arrays, since a 14 bp poly-A tract was previously found to influence the activity of the S. cerevisiae ARG4 meiotic recombination hotspot [11]. Short poly-G runs are somewhat enriched in hotspots, and short poly-A is under-represented, but these differences can partly be explained by elevated GC content in hotspots, which has been shown previously [42], since correlations between DSB intensity and short mononucleotide runs are up to 50% weaker for IGRs, and are almost completely absent for ORFs, when controlling for GC content using partial correlation analysis (Additional file 1, Tables S3 and S4). For long microsatellites other than poly-G, correlations with DSB intensity are generally increased when controlling for GC-content (Additional file 1, Tables S3 and S4).

Dinucleotide repeats of six copies or more, and especially those with ten copies or more, are strongly associated with both hot IGRs and hot ORFs, with poly-AT the most abundant type of repeat involved (Figure 1, Table 2). Trinucleotide repeats of more than six copies are approximately twice as frequent in hot than non hot IGRs (p = 0.0027 Mann-Whitney U Test). This association is not quite significant when using the conservative Bonferroni correction for multiple hypotheses (alpha = 0.0026, see Methods section), but trinucleotide microsatellites are much scarcer than mono- or dinucleotide repeats in the yeast genome (Table 1), so statistical power to detect effects on their distribution is lower.

More marginal associations are present for some other repeat types. Long hexanucleotide microsatellites are many fold more frequent in hot than non-hot ORFs (p < 0.0001, Mann-Whitney U Test; Table 2), but this should be considered in view of the very small numbers of hexanucleotide repeats throughout the genome (Table 1). Dinucleotide repeats with between three and five copies are also significantly over-represented in hot compared with non hot IGRs, but levels of enrichment are much lower than for longer microsatellites (Table 2). Frequency of two-copy repeats is not significantly different in hot compared with non hot regions, despite the great abundance of these repeats relative to longer microsatellites, and the consequent high statistical power. Tetra- and pentanucleotide microsatellites show no significant associations at all, but these repeat types are relatively very rare throughout the yeast genome (Table 1).

Properties of hotspot-associated microsatellites

We examined repeat array length and purity (number of mismatches with respect to the consensus repeated motif) for microsatellites of at least six copies in hotspots and other regions of the yeast genome. In addition, we compared the frequencies of insertion, substitution and deletion mismatches, with respect to the consensus repeated motifs, between hotspot-associated microsatellites and those in other regions. We found that poly-A and poly-G arrays are significantly longer in hot IGRs, and mismatched dinucleotide repeats of at least six copies are significantly longer in hot ORFs, but we saw no other significant differences in repeat length (Additional file 1, Tables S7 and S8). Microsatellites in hot and non-hot regions do not differ significantly in purity, but dinucleotide repeats in non-hot regions do show an elevated proportion of deletion mismatches (p = 0.0006, Mann-Whitney U test).

We looked the sequence motifs of all microsatellites with repeat period between three and six to see if any particular motifs were associated with hotspots. No obvious associations were seen, but we did note that poly-purine/poly-pyrimidine motifs with only one G or C are clearly over-represented among the most common motifs for low copy repeats in both hot and non-hot regions (Additional file 1, Tables S9–S12). This is likely to be related to the enrichment of poly-purine/poly-pyrimidine tracts (PPTs) in the genome as a whole [56], and, as we have reported previously, PPTs with internal tandem repeats comprise only a small proportion of total PPTs [57]. The GC-content of all repeats with at least six copies is strikingly low in IGRs throughout the genome, but there are no significant differences between hot and non-hot regions for microsatellite GC-content (Additional file 1, Tables S5 and S6).

Possible complicating factors

The influence of microsatellites on transcriptional frequency [4–10], and the mutagenic effect of transcription on microsatellites [58] suggested that factors relating to gene expression could affect microsatellite distribution. Theoretically, this could drive the association between microsatellites and recombination hotspots in yeast, since transcriptional frequency (vegetative cells [59]) correlates with DSB intensity (p < 0.0001). However, looking at the "hottest" regions for transcriptional frequency (in equivalent numbers to the numbers of recombination hot regions studied), we found that the number of these that overlap with recombination hotspots is lower than random expectation, and the correlations between DSB intensity and frequency of microsatellites change very little when controlling for transcriptional frequency in partial correlation analysis (Additional file 1, Tables S3 and S4). DSBs have been shown to be more frequent in IGRs with two promoters (divergent transcription of flanking genes) than those with one (parallel transcription of flanking genes) or none (convergent transcription of flanking genes) [42]. We found that densities of some types of microsatellite do differ between IGRs with different numbers of promoters (Table S13). Significant differences are not present for longer microsatellites, however, with the exception of dinucleotide repeats, which are more common in IGRs with no promoters, though not significantly so when testing hot IGRs only. The association between poly-A and hotspots is not due to factors relating to the poly-A adenylation signal present in 3' untranslated regions (UTRs), since the level of enrichment of poly-A in hot over non hot IGRs does not differ by more than 5% between regions with zero, one and two promoters (two, one and zero 3' UTRs respectively).

Another factor that could complicate the association between hotspots and microsatellites is complex (tightly bunched or highly degenerate) repeats. Our initial analysis left open this possibility, since our repeat-finding algorithm does not allow multiple consecutive mismatches within single microsatellites. We therefore looked at numbers of repeats within five and ten bp of other repeats, and compared levels between hot and non-hot regions (Additional file 1, Tables S14 and S15). We found that numbers of microsatellites within complex repeats in IGRs are similar in hot and non-hot, or somewhat higher in non-hot, regions. Degenerate or complex repeats do not, therefore, affect the association between microsatellites and hot IGRs. In ORFs, complex repeats are generally somewhat more frequent in hot regions, however, and this is the case for one repeat type that showed significant over-abundance in hot ORFs, namely long dinucleotide repeats of at least six and at least ten copies. This raised the question of whether the association between this type of repeat and hot ORFs is due to the presence of highly mismatched repeats counted multiple times by our repeat finder. We therefore repeated the analysis with dinuclceotide microsatellites in ORFs occurring within 5 bp of other dinucleotide microsatellites grouped together as single arrays. This did not change the results for repeats with at least 10 copies, which still showed a strong association with hot ORFs (p < 0.0001). It did, however, reduce the significance of the association between dinucleotide arrays of at least six copies and hotspots, raising the p value to 0.014, which is above our alpha level. In view of this result, we removed dinucleotide repeats with at least six copies from our list of repeat types associated with hot ORFs.

Microsatellite frequencies in hotspot flanking regions

We reported previously that PPTs are enriched in hotspot flanking regions as far as two ORFs removed from hotspots [57]. We repeated the analysis for microsatellites, but found no consistent evidence for a similar regional enrichment (Additional file 1, Tables S16 and S17). This suggests that the association with recombination hotspots is less broad in scale for microsatellites than for PPTs. It is also possible, however, that the lower relative abundance of microsatellites could obscure a more general broad scale association than we were able to detect, since several repeat types have higher mean frequencies in hotspot flanking regions but are too sparse for statistical significance. Furthermore, since microsatellites are enriched in both hot IGRs and hot ORFs as defined by the DSB map by Gerton et al., [42], and recombination breakpoints mapped on the finest possible scale are concentrated almost entirely in IGRs in yeast [55], the relationship between microsatellites and recombination probably is distal to some degree.

Discussion

The level of enrichment of microsatellites in yeast recombination hotspots we have detailed here is considerably greater than has been seen for human hotspots [43, 45]. It is not clear why this should be the case, but it is notable that the association between microsatellites and recombination in mammals is quite marked when considering broad scales of several hundred thousand kilo bases or more [49, 50]. In view of evidence that humans and chimpanzees do not share a large proportion of hotspot locations in common [60, 61], one explanation for the discrepancy could be that hotspots do not stay in one place long enough, in these species, to leave strong local imprints in the form of simple sequences generated by hotspot-associated factors, but that hotspot density is more constant on a larger scale. Lower lability of yeast hotspots in evolutionary time could therefore, in theory, have resulted in the stronger associations we have seen.

A better-characterized difference between the yeast and human genomes, which could also contribute to the difference between the two species in the level of association between hotspots and microsatellite abundance, is the vastly greater amount of non-coding DNA in humans. Yeast intergenic regions are small, averaging only just over 500 bp, and 75% of them contain promoters. Potentially, this could complicate the association between recombination hotspots and microsatellites due to the links between microsatellites, transcription, and recombination. Our findings suggest that this is not the case, however. It is also unlikely that transposable, or known repetitive, elements mediate the link between recombination hotspots and microsatellites in yeast, since they are not over-represented in the yeast hotspots we studied [42].

The two most obvious factors that could contribute to the association are a mutation bias, relating to recombination, or some other property of hotspot regions, causing microsatellites to form and grow, and regulation of hotspot locations by simple sequences. We attempted to isolate evidence for a mutagenic effect of recombination on microsatellites by investigating short arrays, as these are not likely to be significantly effected by replication slippage, and there is no available evidence to suggest that they have the potential to stimulate recombination. We did not find strong associations with hotspots for low-copy repeats, however, and previous evidence suggests that long microsatellites have the potential to stimulate recombination, as well as to be mutated by it. Some previous findings have cast doubt on the possibility that these phenomena have a widespread influence, and this has limited the amount of attention they have so far been given, but other evidence, including our results, suggests that they should be tested further.

Evidence that microsatellites could play a role in regulating recombination has been found at a chromosomal level in S.cerevisiae for poly-A [11], poly-AC [12, 14] and pentanucleotide [13] arrays, and using extra-chromosomal DNA molecules for several repeat types [62–66]. The existence of hotspots without local microsatellites does not rule out a functional role for the sequences in recombination, since it has been established that mechanisms of hotspot regulation are heterogeneous [46, 48, 67]. High frequencies of microsatellites in some regions outside hotspots are also not conclusive evidence against their functional involvement, since the control of hotspot location has been shown to be complex and multi-levelled, with local and distal sequences, transcription factor binding and chromatin structure alterations all implicated (reviewed in [46, 48, 67]). The ability of microsatellites to bind transcription factors [68], and to affect chromatin structure in vitro [69] and in vivo [70], therefore suggest two ways in which they could function to potentiate recombination at a subset of hotspots. This could happen without DSBs actually occurring in microsatellites; deletion of a 14 bp poly-A tract reduced activity of the yeast ARG4 hotspot by 75% despite the fact that DSBs avoid poly-A [71, 72].

It is also plausible that recombination is involved in some proportion of microsatellite mutations. The vast, presently unexplained, differences in mutation rates between loci (reviewed in [30, 73]) suggest the involvement of heterogeneous mutational mechanisms or regional mutation biases. In model organisms, evidence has been found both for [12, 33, 74] and against [34, 35] a role for recombination, in the mutation of different types of microsatellite. Studies have shown microsatellite mutation rates on the human Y Chromosome to be similar to autosomal levels [36–38], but concluding from these that recombination does not play a role in microsatellite evolution is problematic, since the Y chromosome undergoes intramolecular recombination [75]. It is therefore possible that meiotic recombination, or other properties of its hotspots, could contribute to the variability in microsatellite mutation rates at different chromosomal locations. Although unequal crossing over, or meiotic gene conversion (recombination without exchange of flanking markers), are the most obvious mechanisms for this, other factors could be important, such as replication pausing, which has been linked to microsatellite mutations [76, 77], and may be causally involved in a subset of recombination hotspots [67].

Conclusion

We found that high-copy, short-motif microsatellites are strongly associated with S. cerevisiae meiotic recombination hotspots. The association is weak or absent for low-copy repeats. Our results add to the weight of evidence in favour of further studying the link between microsatellites and recombination hotspots. Large-scale experimental studies in yeast could be used to quantify the level of influence hotspots have on microsatellite evolution, and to explore the possible functional role of microsatellites in regulating recombination. This work could include tracking microsatellite mutations in mono-clonal yeast populations from recombining and non recombining strains. The effect on recombination frequency of deleting microsatellites from hotspots could also be tested.

Methods

Figures for transcriptional activity were from the study by Holstege and co-workers (1998) who mapped transcription frequency in vegetative cells for each yeast ORF [59]. For IGRs, we took the mean of the two adjacent ORFs.

Detection of microsatellites

We detected microsatellites in the yeast genome using an algorithm written in C [78]. The programme initially generated databases of all non-overlapping repeats of two copies or greater for repeated motif sizes two to six bp, and three copies or greater for mononucleotide arrays. Separate databases were created for perfect repeats, arrays with a maximum of one mismatch allowed per ten bp of repeat sequence, and arrays with a maximum of one mismatch per six bp. Microsatellites overlapping two regions were excluded from the analysis. This occurred for less than one percent of arrays overall.

Categorization of microsatellites

Copy number groups were two, three to five, six or more and ten or more. For mononucleotide repeats, we used 14 or more instead of ten or more, since a 14 bp poly-A tract has been shown to be a functional component of a yeast recombination hotspot [11]. We divided mononucleotide microsatellites into the equivalent motif groups A/T and G/C. Dinculeotide microsatellites were considered as a whole for statistical comparisons, but we divided them into the motif groups AT/TA, AC/CA/TG/GT, AG/GA/TC/CT and CG/GC in order to see the relative abundance of each motif type within the class. We examined sequence motifs of microsatellites with three to six bp motifs visually. We investigated compound and highly degenerate microsatellites by looking at numbers of arrays within five or ten bp of another microsatellite of the same or larger copy number group.

Statistical analysis

Statistical comparison of means (Student's T-test and Mann-Whitney U Test, 2-tailed tests in call cases) and correlation analyses (Spearman's Rho) were done using SPSS or SAS. We initially tested the distribution of each sample for normality (Kolmogorov-Smirnov Test) and subjected significantly non-normal samples only to non-parametric tests. Because repeats were divided into 19 physically independent categories for statistical testing, Bonferroni's correction was used to set the alpha level at 0.05/19 = 0.0026. For the purpose of this calculation, the number of categories did not include different mismatch types, because, within motif and size classes, these overlap substantially so are not independent from each other. For the same reason, the size class six copies and longer was not considered to be independent of the class 10 copies and longer for the purpose of calculating the number of independent categories. Bonferroni's correction is clearly very conservative for this study, because we lose statistical power with increasing numbers of categories due to the fact that there are proportionally fewer microsatellites in each category. We would therefore gain a large amount of power by limiting the categorization to a 4-way division of microsatellites into short and long mononucleotide repeats, and short and long 2–6 bp motif repeats. This would not change the main conclusions of the paper, because, for all motif lengths, long microsatellites are either more frequent in hotspots or are extremely rare (Additional file 1, Tables S1 and S2). Some interesting information would be lost with this scheme, since poly-A and dinucleotide repeats are highly predominant among long microsatellites, and two-copy repeats are vastly more frequent than 3–5 copy repeats (Additional file 1, Tables S1 and S2), so we favoured the 19-way division.

References

  1. Charlesworth B, Sniegowski P, Stephan W: The evolutionary dynamics of repetitive DNA in eukaryotes. Nature. 1994, 371: 215-220. 10.1038/371215a0.

    Article  PubMed  Google Scholar 

  2. Lander E, Linton LM, Birren B, al. : Initial sequencing and analysis of the human genome. Nature. 2001, 409: 860-921. 10.1038/35057062.

    Article  PubMed  Google Scholar 

  3. Li B, Xia Q, Lu C, Zhou Z, Xiang Z: Analysis on frequency and density of microsatellites in coding sequences of several eukaryotic genomes. Genomics Proteomics Bioinformatics. 2004, 2 (1): 24-31.

    PubMed  Google Scholar 

  4. Struhl K: Naturally occurring poly(dA-dT) sequences are upstream promoter elements for constitutive transcription in yeast. Proc Natl Acad Sci U S A. 1985, 82 (24): 8419-8423. 10.1073/pnas.82.24.8419.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Uhlemann AC, Szlezak NA, Vonthein R, Tomiuk J, Emmer SA, Lell B, Kremsner PG, Kun JF: DNA phasing by TA dinucleotide microsatellite length determines in vitro and in vivo expression of the gp91phox subunit of NADPH oxidase and mediates protection against severe malaria. J Infect Dis. 2004, 189 (12): 2227-2234. 10.1086/421242.

    Article  PubMed  Google Scholar 

  6. Curi RA, Oliveira HN, Silveira AC, Lopes CR: Effects of polymorphic microsatellites in the regulatory region of IGF1 and GHR on growth and carcass traits in beef cattle. Anim Genet. 2005, 36 (1): 58-62. 10.1111/j.1365-2052.2004.01226.x.

    Article  PubMed  Google Scholar 

  7. Contente A, Dittmer A, Koch MC, Roth J, Dobbelstein M: A polymorphic microsatellite that mediates induction of PIG3 by p53. Nat Genet. 2002, 30 (3): 315-320. 10.1038/ng836.

    Article  PubMed  Google Scholar 

  8. Borrmann L, Seebeck B, Rogalla P, Bullerdiek J: Human HMGA2 promoter is coregulated by a polymorphic dinucleotide (TC)-repeat. Oncogene. 2003, 22 (5): 756-760. 10.1038/sj.onc.1206073.

    Article  PubMed  Google Scholar 

  9. Hammock EA, Young LJ: Microsatellite instability generates diversity in brain and sociobehavioral traits. Science. 2005, 308 (5728): 1630-1634. 10.1126/science.1111427.

    Article  PubMed  Google Scholar 

  10. Hammock EA, Young LJ: Functional microsatellite polymorphism associated with divergent social structure in vole species. Mol Biol Evol. 2004, 21 (6): 1057-1063. 10.1093/molbev/msh104.

    Article  PubMed  Google Scholar 

  11. Schultes NP, Szostak JW: A poly(dA.dT) tract is a component of the recombination initiation site at the ARG4 locus in Saccharomyces cerevisiae. Mol Cell Biol. 1991, 11 (1): 322-328.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Gendrel CG, Boulet A, Dutreix M: (CA/GT)(n) microsatellites affect homologous recombination during yeast meiosis. Genes Dev. 2000, 14 (10): 1261-1268.

    PubMed  PubMed Central  Google Scholar 

  13. Kirkpatrick DT, Wang YH, Dominska M, Griffith JD, Petes TD: Control of meiotic recombination and gene expression in yeast by a simple repetitive DNA sequence that excludes nucleosomes. Mol Cell Biol. 1999, 19 (11): 7661-7671.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Treco D, Arnheim N: The evolutionarily conserved repetitive sequence d(TG.AC)n promotes reciprocal exchange and generates unusual recombinant tetrads during yeast meiosis. Mol Cell Biol. 1986, 6 (11): 3934-3947.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Dib C, Faure S, Fizames C, Samson D, Drouot N, Vignal A, Millasseau P, Marc S, Hazan J, Seboun E, Lathrop M, Gyapay G, Morissette J, Weissenbach J: A comprehensive genetic map of the human genome based on 5,264 microsatellites. Nature. 1996, 380 (6570): 152-154. 10.1038/380152a0.

    Article  PubMed  Google Scholar 

  16. Ihara N, Takasuga A, Mizoshita K, Takeda H, Sugimoto M, Mizoguchi Y, Hirano T, Itoh T, Watanabe T, Reed KM, Snelling WM, Kappes SM, Beattie CW, Bennett GL, Sugimoto Y: A comprehensive genetic map of the cattle genome based on 3802 microsatellites. Genome Res. 2004, 14 (10A): 1987-1998. 10.1101/gr.2741704.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Dietrich WF, Miller J, Steen R, Merchant MA, Damron-Boles D, Husain Z, Dredge R, Daly MJ, Ingalls KA, O'Connor TJ: A comprehensive genetic map of the mouse genome. Nature. 1996, 380 (6570): 149-152. 10.1038/380149a0.

    Article  PubMed  Google Scholar 

  18. Sibov ST, de Souza CL, Garcia AA, Silva AR, Garcia AF, Mangolin CA, Benchimol LL, de Souza AP: Molecular mapping in tropical maize (Zea mays L.) using microsatellite markers. 2. Quantitative trait loci (QTL) for grain yield, plant height, ear height and grain moisture. Hereditas. 2003, 139 (2): 107-115. 10.1111/j.1601-5223.2003.01667.x.

    Article  PubMed  Google Scholar 

  19. Goris A, Sawcer S, Vandenbroeck K, Carton H, Billiau A, Setakis E, Compston A, Dubois B: New candidate loci for multiple sclerosis susceptibility revealed by a whole genome association screen in a Belgian population. J Neuroimmunol. 2003, 143 (1-2): 65-69. 10.1016/j.jneuroim.2003.08.013.

    Article  PubMed  Google Scholar 

  20. Dirlewanger E, Cosson P, Howad W, Capdeville G, Bosselut N, Claverie M, Voisin R, Poizat C, Lafargue B, Baron O, Laigret F, Kleinhentz M, Arus P, Esmenjaud D: Microsatellite genetic linkage maps of myrobalan plum and an almond-peach hybrid--location of root-knot nematode resistance genes. Theor Appl Genet. 2004, 109 (4): 827-838. 10.1007/s00122-004-1694-9.

    Article  PubMed  Google Scholar 

  21. Tamaki K, Jeffreys AJ: Human tandem repeat sequences in forensic DNA typing. Leg Med (Tokyo). 2005, 7 (4): 244-250.

    Article  Google Scholar 

  22. Webster MS, Reichart L: Use of microsatellites for parentage and kinship analyses in animals. Methods Enzymol. 2005, 395: 222-238.

    Article  PubMed  Google Scholar 

  23. Schlotterer C, Pemberton J: The use of microsatellites for genetic analysis of natural populations. Exs. 1994, 69: 203-214.

    PubMed  Google Scholar 

  24. Rosenberg NA, Pritchard JK, Weber JL, Cann HM, Kidd KK, Zhivotovsky LA, Feldman MW: Genetic structure of human populations. Science. 2002, 298 (5602): 2381-2385. 10.1126/science.1078311.

    Article  PubMed  Google Scholar 

  25. Hayano A, Yoshioka M, Tanaka M, Amano M: Population differentiation in the Pacific white-sided dolphin Lagenorhynchus obliquidens inferred from mitochondrial DNA and microsatellite analyses. Zoolog Sci. 2004, 21 (9): 989-999. 10.2108/zsj.21.989.

    Article  PubMed  Google Scholar 

  26. Bowcock AM, Ruiz-Linares A, Tomfohrde J, Minch E, Kidd JR, Cavalli-Sforza LL: High resolution of human evolutionary trees with polymorphic microsatellites. Nature. 1994, 368 (6470): 455-457. 10.1038/368455a0.

    Article  PubMed  Google Scholar 

  27. Meyer E, Wiegand P, Rand SP, Kuhlmann D, Brack M, Brinkmann B: Microsatellite polymorphisms reveal phylogenetic relationships in primates. J Mol Evol. 1995, 41 (1): 10-14. 10.1007/BF00174036.

    Article  PubMed  Google Scholar 

  28. Schlotterer C: Genealogical inference of closely related species based on microsatellites. Genet Res. 2001, 78 (3): 209-212.

    Article  PubMed  Google Scholar 

  29. Buschiazzo E, Gemmell NJ: The rise, fall and renaissance of microsatellites in eukaryotic genomes. Bioessays. 2006, 28 (10): 1040-1050. 10.1002/bies.20470.

    Article  PubMed  Google Scholar 

  30. Ellegren H: Microsatellites: simple sequences with complex evolution. Nat Rev Genet. 2004, 5 (6): 435-445. 10.1038/nrg1348.

    Article  PubMed  Google Scholar 

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

    PubMed  Google Scholar 

  32. Richard GF, Paques F: Mini- and microsatellite expansions: the recombination connection. EMBO Rep. 2000, 1 (2): 122-126. 10.1093/embo-reports/kvd031.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Jakupciak JP, Wells RD: Genetic instabilities of triplet repeat sequences by recombination. IUBMB Life. 2000, 50 (6): 355-359.

    Article  PubMed  Google Scholar 

  34. Levinson G, Gutman GA: High frequencies of short frameshifts in poly-CA/TG tandem repeats borne by bacteriophage M13 in Escherichia coli K-12. Nucleic Acids Res. 1987, 15 (13): 5323-5338. 10.1093/nar/15.13.5323.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Henderson ST, Petes TD: Instability of simple sequence DNA in Saccharomyces cerevisiae. Mol Cell Biol. 1992, 12 (6): 2749-2757.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Gusmao L, Sanchez-Diz P, Calafell F, Martin P, Alonso CA, Alvarez-Fernandez F, Alves C, Borjas-Fajardo L, Bozzo WR, Bravo ML, Builes JJ, Capilla J, Carvalho M, Castillo C, Catanesi CI, Corach D, Di Lonardo AM, Espinheira R, Fagundes de Carvalho E, Farfan MJ, Figueiredo HP, Gomes I, Lojo MM, Marino M, Pinheiro MF, Pontes ML, Prieto V, Ramos-Luis E, Riancho JA, Souza Goes AC, Santapa OA, Sumita DR, Vallejo G, Vidal Rioja L, Vide MC, Vieira da Silva CI, Whittle MR, Zabala W, Zarrabeitia MT, Alonso A, Carracedo A, Amorim A: Mutation rates at Y chromosome specific microsatellites. Hum Mutat. 2005, 26 (6): 520-528. 10.1002/humu.20254.

    Article  PubMed  Google Scholar 

  37. Kayser M, Roewer L, Hedman M, Henke L, Henke J, Brauer S, Kruger C, Krawczak M, Nagy M, Dobosz T, Szibor R, de Knijff P, Stoneking M, Sajantila A: Characteristics and frequency of germline mutations at microsatellite loci from the human Y chromosome, as revealed by direct observation in father/son pairs. Am J Hum Genet. 2000, 66 (5): 1580-1588. 10.1086/302905.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Nebel A, Filon D, Hohoff C, Faerman M, Brinkmann B, Oppenheim A: Haplogroup-specific deviation from the stepwise mutation model at the microsatellite loci DYS388 and DYS392. Eur J Hum Genet. 2001, 9 (1): 22-26. 10.1038/sj.ejhg.5200577.

    Article  PubMed  Google Scholar 

  39. Huang QY, Xu FH, Shen H, Deng HY, Liu YJ, Liu YZ, Li JL, Recker RR, Deng HW: Mutation patterns at dinucleotide microsatellite loci in humans. Am J Hum Genet. 2002, 70 (3): 625-634. 10.1086/338997.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Payseur BA, Nachman MW: Microsatellite variation and recombination rate in the human genome. Genetics. 2000, 156 (3): 1285-1298.

    PubMed  PubMed Central  Google Scholar 

  41. Jeffreys AJ, Kauppi L, Neumann R: Intensely punctate meiotic recombination in the class II region of the major histocompatibility complex. Nat Genet. 2001, 29 (2): 217-222. 10.1038/ng1001-217.

    Article  PubMed  Google Scholar 

  42. Gerton JL, DeRisi J, Shroff R, Lichten M, Brown PO, Petes TD: Inaugural article: global mapping of meiotic recombination hotspots and coldspots in the yeast Saccharomyces cerevisiae. Proc Natl Acad Sci U S A. 2000, 97 (21): 11383-11390. 10.1073/pnas.97.21.11383.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Myers S, Bottolo L, Freeman C, McVean G, Donnelly P: A fine-scale map of recombination rates and hotspots across the human genome. Science. 2005, 310 (5746): 321-324. 10.1126/science.1117196.

    Article  PubMed  Google Scholar 

  44. Jeffreys AJ, Neumann R, Panayi M, Myers S, Donnelly P: Human recombination hot spots hidden in regions of strong marker association. Nat Genet. 2005, 37 (6): 601-606. 10.1038/ng1565.

    Article  PubMed  Google Scholar 

  45. Brandström M, Bagshaw ATM, Gemmell NJ, Ellegren H: In preparation.

  46. Nishant KT, Rao MR: Molecular features of meiotic recombination hot spots. Bioessays. 2006, 28 (1): 45-56. 10.1002/bies.20349.

    Article  PubMed  Google Scholar 

  47. Jeffreys AJ, Neumann R: Factors influencing recombination frequency and distribution in a human meiotic crossover hotspot. Hum Mol Genet. 2005, 14 (15): 2277-2287. 10.1093/hmg/ddi232.

    Article  PubMed  Google Scholar 

  48. Kauppi L, Jeffreys AJ, Keeney S: Where the crossovers are: recombination distributions in mammals. Nat Rev Genet. 2004, 5 (6): 413-424. 10.1038/nrg1346.

    Article  PubMed  Google Scholar 

  49. Jensen-Seaman MI, Furey TS, Payseur BA, Lu Y, Roskin KM, Chen CF, Thomas MA, Haussler D, Jacob HJ: Comparative recombination rates in the rat, mouse, and human genomes. Genome Res. 2004, 14 (4): 528-538. 10.1101/gr.1970304.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Kong A, Gudbjartsson DF, Sainz J, Jonsdottir GM, Gudjonsson SA, Richardsson B, Sigurdardottir S, Barnard J, Hallbeck B, Masson G, Shlien A, Palsson ST, Frigge ML, Thorgeirsson TE, Gulcher JR, Stefansson K: A high-resolution recombination map of the human genome. Nat Genet. 2002, 31 (3): 241-247.

    PubMed  Google Scholar 

  51. Messier W, Li SH, Stewart CB: The birth of microsatellites. Nature. 1996, 381 (6582): 483-10.1038/381483a0.

    Article  PubMed  Google Scholar 

  52. Brinkmann B, Klintschar M, Neuhuber F, Huhne J, Rolf B: Mutation rate in human microsatellites: influence of the structure and length of the tandem repeat. Am J Hum Genet. 1998, 62 (6): 1408-1415. 10.1086/301869.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Weber JL: Informativeness of human (dC-dA)n.(dG-dT)n polymorphisms. Genomics. 1990, 7 (4): 524-530. 10.1016/0888-7543(90)90195-Z.

    Article  PubMed  Google Scholar 

  54. Zhu Y, Queller DC, Strassmann JE: A phylogenetic perspective on sequence evolution in microsatellite loci. J Mol Evol. 2000, 50 (4): 324-338.

    PubMed  Google Scholar 

  55. Baudat F, Nicolas A: Clustering of meiotic double-strand breaks on yeast chromosome III. Proc Natl Acad Sci U S A. 1997, 94 (10): 5213-5218. 10.1073/pnas.94.10.5213.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Raghavan S, Burma PK, Brahmachari SK: Positional preferences of polypurine/polypyrimidine tracts in Saccharomyces cerevisiae genome: implications for cis regulation of gene expression. J Mol Evol. 1997, 45 (5): 485-498. 10.1007/PL00006253.

    Article  PubMed  Google Scholar 

  57. Bagshaw AT, Pitt JP, Gemmell NJ: Association of poly-purine/poly-pyrimidine sequences with meiotic recombination hot spots. BMC Genomics. 2006, 7: 179-10.1186/1471-2164-7-179.

    Article  PubMed  PubMed Central  Google Scholar 

  58. Wierdl M, Greene CN, Datta A, Jinks-Robertson S, Petes TD: Destabilization of simple repetitive DNA sequences by transcription in yeast. Genetics. 1996, 143 (2): 713-721.

    PubMed  PubMed Central  Google Scholar 

  59. Holstege FC, Jennings EG, Wyrick JJ, Lee TI, Hengartner CJ, Green MR, Golub TR, Lander ES, Young RA: Dissecting the regulatory circuitry of a eukaryotic genome. Cell. 1988, 95 (5): 717-728. 10.1016/S0092-8674(00)81641-4.

    Article  Google Scholar 

  60. Ptak SE, Hinds DA, Koehler K, Nickel B, Patil N, Ballinger DG, Przeworski M, Frazer KA, Paabo S: Fine-scale recombination patterns differ between chimpanzees and humans. Nat Genet. 2005, 37 (4): 429-434. 10.1038/ng1529.

    Article  PubMed  Google Scholar 

  61. Winckler W, Myers SR, Richter DJ, Onofrio RC, McDonald GJ, Bontrop RE, McVean GA, Gabriel SB, Reich D, Donnelly P, Altshuler D: Comparison of fine-scale recombination rates in humans and chimpanzees. Science. 2005, 308 (5718): 107-111. 10.1126/science.1105322.

    Article  PubMed  Google Scholar 

  62. Murphy KE, Stringer JR: RecA independent recombination of poly[d(GT)-d(CA)] in pBR322. Nucleic Acids Res. 1986, 14 (18): 7325-7340. 10.1093/nar/14.18.7325.

    Article  PubMed  PubMed Central  Google Scholar 

  63. Bullock P, Miller J, Botchan M: Effects of poly[d(pGpT).d(pApC)] and poly[d(pCpG).d(pCpG)] repeats on homologous recombination in somatic cells. Mol Cell Biol. 1986, 6 (11): 3948-3953.

    Article  PubMed  PubMed Central  Google Scholar 

  64. Napierala M, Dere R, Vetcher A, Wells RD: Structure-dependent recombination hot spot activity of GAA.TTC sequences from intron 1 of the Friedreich's ataxia gene. J Biol Chem. 2004, 279 (8): 6444-6454. 10.1074/jbc.M309596200.

    Article  PubMed  Google Scholar 

  65. Wahls WP, Wallace LJ, Moore PD: The Z-DNA motif d(TG)30 promotes reception of information during gene conversion events while stimulating homologous recombination in human cells in culture. Mol Cell Biol. 1990, 10 (2): 785-793.

    Article  PubMed  PubMed Central  Google Scholar 

  66. Napierala M, Parniewski P, Pluciennik A, Wells RD: Long CTG.CAG repeat sequences markedly stimulate intramolecular recombination. J Biol Chem. 2002, 277 (37): 34087-34100. 10.1074/jbc.M202128200.

    Article  PubMed  Google Scholar 

  67. Petes TD: Meiotic recombination hot spots and cold spots. Nat Rev Genet. 2001, 2 (5): 360-369. 10.1038/35072078.

    Article  PubMed  Google Scholar 

  68. Lu Q, Teare JM, Granok H, Swede MJ, Xu J, Elgin SC: The capacity to form H-DNA cannot substitute for GAGA factor binding to a (CT)n*(GA)n regulatory site. Nucleic Acids Res. 2003, 31 (10): 2483-2494. 10.1093/nar/gkg369.

    Article  PubMed  PubMed Central  Google Scholar 

  69. Wang YH, Amirhaeri S, Kang S, Wells RD, Griffith JD: Preferential nucleosome assembly at DNA triplet repeats from the myotonic dystrophy gene. Science. 1994, 265 (5172): 669-671. 10.1126/science.8036515.

    Article  PubMed  Google Scholar 

  70. Otten AD, Tapscott SJ: Triplet repeat expansion in myotonic dystrophy alters the adjacent chromatin structure. Proc Natl Acad Sci U S A. 1995, 92 (12): 5465-5469. 10.1073/pnas.92.12.5465.

    Article  PubMed  PubMed Central  Google Scholar 

  71. Liu J, Wu TC, Lichten M: The location and structure of double-strand DNA breaks induced during yeast meiosis: evidence for a covalently linked DNA-protein intermediate. Embo J. 1995, 14 (18): 4599-4608.

    PubMed  PubMed Central  Google Scholar 

  72. de Massy B, Rocco V, Nicolas A: The nucleotide mapping of DNA double-strand breaks at the CYS3 initiation site of meiotic recombination in Saccharomyces cerevisiae. Embo J. 1995, 14 (18): 4589-4598.

    PubMed  PubMed Central  Google Scholar 

  73. Ellegren H: Microsatellite mutations in the germline: implications for evolutionary inference. Trends Genet. 2000, 16 (12): 551-558. 10.1016/S0168-9525(00)02139-9.

    Article  PubMed  Google Scholar 

  74. Hashem VI, Rosche WA, Sinden RR: Genetic recombination destabilizes (CTG)n.(CAG)n repeats in E. coli. Mutat Res. 2004, 554 (1-2): 95-109.

    Article  PubMed  Google Scholar 

  75. Rozen S, Skaletsky H, Marszalek JD, Minx PJ, Cordum HS, Waterston RH, Wilson RK, Page DC: Abundant gene conversion between arms of palindromes in human and ape Y chromosomes. Nature. 2003, 423 (6942): 873-876. 10.1038/nature01723.

    Article  PubMed  Google Scholar 

  76. Fouche N, Ozgur S, Roy D, Griffith JD: Replication fork regression in repetitive DNAs. Nucleic Acids Res. 2006, 34 (20): 6044-6050. 10.1093/nar/gkl757.

    Article  PubMed  PubMed Central  Google Scholar 

  77. Hile SE, Eckert KA: Positive correlation between DNA polymerase pausing and mutagenesis within polypyrimidine/polypurine microsatellite sequences. J Mol Biol. 2004, 335 (3): 745-759. 10.1016/j.jmb.2003.10.075.

    Article  PubMed  Google Scholar 

  78. [http://repeatfinder.sourceforge.net/]

Download references

Acknowledgements

Financial support for this work came from a Royal Society of New Zealand Marsden grant (UOC 202).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Andrew TM Bagshaw.

Additional information

Authors' contributions

ATMB conceived and designed the experiments, analyzed the data and wrote the paper. JPWP wrote the computer programme. NJG contributed to the interpretation of the data and the writing of the paper. All authors read and approved the final manuscript.

Electronic supplementary material

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

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

Bagshaw, A.T., Pitt, J.P. & Gemmell, N.J. High frequency of microsatellites in S. cerevisiae meiotic recombination hotspots. BMC Genomics 9, 49 (2008). https://0-doi-org.brum.beds.ac.uk/10.1186/1471-2164-9-49

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/1471-2164-9-49

Keywords