Distribution and abundance of microsatellites in the yeast genome can Be explained by a balance between slippage events and point mutations

We fit a Markov chain model of microsatellite evolution introduced by Kruglyak et al. to data on all di-, tri-, and tetranucleotide repeats in the yeast genome. Our results suggest that many features of the distribution of abundance and length of microsatellites can be explained by this simple model, which incorporates a competition between slippage events and base pair substitutions, with no need to invoke selection or constraints on the lengths. Our results provide some new information on slippage rates for individual repeat motifs, which suggest that AT-rich trinucleotide repeats have higher slippage rates. As our model predicts, we found that many repeats were adjacent to shorter repeats of the same motif. However, we also found a significant tendency of microsatellites of different motifs to cluster.     

The structure of interrupted human AC microsatellites

Microsatellite lengths change over evolutionary time through a process of replication slippage. A recently proposed model of this process holds that the expansionary tendencies of slippage mutation are balanced by point mutations breaking longer microsatellites into smaller units and that this process gives rise to the observed frequency distributions of uninterrupted microsatellite lengths. We refer to this as the slippage/point-mutation theory. Here we derive the theory's predictions for interrupted microsatellites comprising regions of perfect repeats, labeled segments, separated by dinucleotide interruptions containing point mutations. These predictions are tested by reference to the frequency distributions of segments of AC microsatellite in the human genome, and several predictions are shown not to be supported by the data, as follows. The estimated slippage rates are relatively low for the first four repeats, and then rise initially linearly with length, in accordance with previous work. However, contrary to expectation and the experimental evidence, the inferred slippage rates decline in segments above 10 repeats. Point mutation rates are also found to be higher within microsatellites than elsewhere. The theory provides an excellent fit to the frequency distribution of peripheral segment lengths but fails to explain why internal segments are shorter. Furthermore, there are fewer microsatellites with many segments than predicted. The frequencies of interrupted microsatellites decline geometrically with microsatellite size measured in number of segments, so that for each additional segment, the number of microsatellites is 33.6% less. Overall we conclude that the detailed structure of interrupted microsatellites cannot be reconciled with the existing slippage/point-mutation theory of microsatellite evolution, and we suggest that microsatellites are stabilized by processes acting on interior rather than on peripheral segments.     

The relationship between microsatellite slippage mutation rate and the number of repeat units

Microsatellite markers are widely used for genetic studies, but the relationship between microsatellite slippage mutation rate and the number of repeat units remains unclear. In this study, microsatellite distributions in the human genome are collected from public sequence databases. We observe that there is a threshold size for slippage mutations. We consider a model of microsatellite mutation consisting of point mutations and single stepwise slippage mutations. From two sets of equations based on two stochastic processes and equilibrium assumptions, we estimate microsatellite slippage mutation rates without assuming any relationship between microsatellite slippage mutation rate and the number of repeat units. We use the least squares method with constraints to estimate expansion and contraction mutation rates. The estimated slippage mutation rate increases exponentially as the number of repeat units increases. When slippage mutations happen, expansion occurs more frequently for short microsatellites and contraction occurs more frequently for long microsatellites. Our results agree with the length-dependent mutation pattern observed from experimental data, and they explain the scarcity of long microsatellites.     

Dinucleotide repeats in the Drosophila and human genomes have complex,length-dependent mutation processes

We use methods of maximum likelihood estimation to fit several microsatellite mutation models to the observed length distribution of dinucletoide repeats in the Drosophila and human genomes. All simple models are rejected by this procedure. Two new models, one with quadratic and another with piecewise linear slippage rates, have the best fits and agree with recent experimental studies by predicting that long microsatellites have a bias toward contractions.     

Evolution of chloroplast mononucleotide microsatellites in Arabidopsis thaliana

The level of variation and the mutation rate were investigated in an empirical study of 244 chloroplast microsatellites in 15 accessions of Arabidopsis thaliana. In contrast to SNP variation, microsatellite variation in the chloroplast was found to be common, although less common than microsatellite variation in the nucleus. No microsatellite variation was found in coding regions of the chloroplast. To evaluate different models of microsatellite evolution as possible explanations for the observed pattern of variation, the length distribution of microsatellites in the published DNA sequence of the A. thaliana chloroplast was subsequently used. By combining information from these two analyses we found that the mode of evolution of the chloroplast mononucleotide microsatellites was best described by a linear relation between repeat length and mutation rate, when the repeat lengths exceeded about 7 bp. This model can readily predict the variation observed in non-coding chloroplast DNA. It was found that the number of uninterrupted repeat units had a large impact on the level of chloroplast microsatellite variation. No other factors investigated—such as the position of a locus within the chromosome, or imperfect repeats—appeared to affect the variability of chloroplast microsatellites. By fitting the slippage models to the Genbank sequence of chromosome 1, we show that the difference between microsatellite variation in the nucleus and the chloroplast is largely due to differences in slippage rate. Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

Creation of a chloroplast microsatellite reporter for detection of replication slippage in Chlamydomonas reinhardtii

Microsatellites are composed of short tandem direct repeats; deletions or duplications of those repeats through the process of replication slippage result in microsatellite instability relative to other genomic loci. Variation in repeat number occurs so frequently that microsatellites can be used for genotyping and forensic analysis. However, an accurate assessment of the rates of change can be difficult because the presence of many repeats makes it difficult to determine whether changes have occurred through single or multiple events. The current study was undertaken to experimentally assess the rates of replication slippage that occur in vivo in the chloroplast DNA of Chlamydomonas reinhardtii. A reporter construct was created in which a stretch of AAAG repeats was inserted into a functional gene to allow changes to be observed when they occurred at the synthetic microsatellite. Restoration of the reading frame occurred through replication slippage in 15 of every million viable cells. Since only one-third of the potential insertion/deletion events would restore the reading frame, the frequency of change could be deduced to be 4.5 x 10(-5). Analysis of the slippage events showed that template slippage was the primary event, resulting in deletions rather than duplications. These findings contrasted with events observed in Escherichia coli during maintenance of the plasmid, where duplications were the rule.     

Long microsatellite alleles in Drosophila melanogaster have a downward mutation bias and short persistence times, which cause their genome-wide underrepresentation

Microsatellites are short tandemly repeated DNA sequence motifs that are highly variable in most organisms. In contrast to mammals, long microsatellites (>15 repeats) are extremely rare in the Drosophila melanogaster genome. To investigate this paucity of long microsatellites in Drosophila, we studied 19 loci with exceptionally long microsatellite alleles. Inter- and intraspecific analysis showed that long microsatellite alleles arose in D. melanogaster only very recently. This lack of old alleles with many repeats indicated that long microsatellite alleles have short persistence times. The size distribution of microsatellite mutations in mutation-accumulation lines suggests that long alleles have a mutation bias toward a reduction in the number of repeat units. This bias causes the short persistence times of long microsatellite alleles. We propose that species-specific, size-dependent mutation spectra of microsatellite alleles may provide a general mechanism to account for the observed differences in microsatellite length between species.     

Microsatellite Interruptions Stabilize Primate Genomes and Exist as Population-Specific Single Nucleotide Polymorphisms within Individual Human Genomes

Interruptions of microsatellite sequences impact genome evolution and can alter disease manifestation. However, human polymorphism levels at interrupted microsatellites (iMSs) are not known at a genome-wide scale, and the pathways for gaining interruptions are poorly understood. Using the 1000 Genomes Phase-1 variant call set, we interrogated mono-, di-, tri-, and tetranucleotide repeats up to 10 units in length. We detected ∼26,000–40,000 iMSs within each of four human population groups (African, European, East Asian, and American). We identified population-specific iMSs within exonic regions, and discovered that known disease-associated iMSs contain alleles present at differing frequencies among the populations. By analyzing longer microsatellites in primate genomes, we demonstrate that single interruptions result in a genome-wide average two- to six-fold reduction in microsatellite mutability, as compared with perfect microsatellites. Centrally located interruptions lowered mutability dramatically, by two to three orders of magnitude. Using a biochemical approach, we tested directly whether the mutability of a specific iMS is lower because of decreased DNA polymerase strand slippage errors. Modeling the adenomatous polyposis coli tumor suppressor gene sequence, we observed that a single base substitution interruption reduced strand slippage error rates five- to 50-fold, relative to a perfect repeat, during synthesis by DNA polymerases α, β, or η. Computationally, we demonstrate that iMSs arise primarily by base substitution mutations within individual human genomes. Our biochemical survey of human DNA polymerase α, β, δ, κ, and η error rates within certain microsatellites suggests that interruptions are created most frequently by low fidelity polymerases. Our combined computational and biochemical results demonstrate that iMSs are abundant in human genomes and are sources of population-specific genetic variation that may affect genome stability. The genome-wide identification of iMSs in human populations presented here has important implications for current models describing the impact of microsatellite polymorphisms on gene expression.     

Length, orientation, and plant host influence the mutation frequency in microsatellites.

Microsatellites are simple, tandem DNA repeats that represent unstable regions of the genome. They undergo frequent changes in tract length by base additions or deletions due to DNA polymerase slippage during replication. To characterize factors affecting the frequency of spontaneous mutations occurring in microsatellites in plants, a reporter system was used in Arabidopsis thaliana and tomato (Lycopersicon esculentum). The beta-glucuronidase (GUS) reporter system was used to measure the mutation frequency in various microsatellites (G(7), G(10), G(13), G(16), and C(16)) in somatic tissues. Our results indicate that this frequency increases with the number of repeats: a G(16) tract was almost 80-fold more mutable than a G(7) tract. Furthermore, the frequency of mutations depends on repeat orientation, as G(16) was 3-fold more mutable than C(16). The mutation rate was also found to differ markedly in Arabidopsis and tomato for an identical microsatellite. Indeed, Arabidopsis showed a 5-fold higher mutation frequency than tomato with the same G(7) reporter construct. Finally, mutation in a G(16) tract was frequent enough that mutations transmitted germinally to the next generation could be detected at a relatively high frequency.     

Microsatellite length differences between humans and chimpanzees at autosomal Loci are not found at equivalent haploid Y chromosomal Loci

When homologous microsatellites are compared between species, significant differences in mean length are often noted. A dominant cause of these length differences is ascertainment bias due to selection for maximum repeat number and repeat purity when the markers are being developed. However, even after ascertainment bias has been allowed for through reciprocal comparisons, significant length differences remain, suggesting that the average microsatellite mutation rate differs between species. Two classes of mechanism have been proposed: rapid evolution of enzymes involved in the generation and repair of slippage products (enzyme evolution model) and heterozygote instability, whereby interchromosomal events at heterozygous sites offer extra opportunities for mutations to occur (heterozygote instability model). To examine which of these hypotheses is most likely, we compared ascertainment bias and species length differences between humans and chimpanzees in autosomal and Y chromosomal microsatellites. We find that levels of ascertainment bias are indistinguishable, but that interspecies length differences are significantly greater for autosomal loci compared with haploid Y chromosomal loci. Such a pattern is consistent with predictions from the heterozygote instability model and is not expected under models of microsatellite evolution that do not include interchromosomal events such as the enzyme evolution model.     

The mutation rates of di-, tri- and tetranucleotide repeats in Drosophila melanogaster

In a recent study, we reported that the combined average mutation rate of 10 di-, 6 tri-, and 8 tetranucleotide repeats in Drosophila melanogaster was 6.3 x 10(-6) mutations per locus per generation, a rate substantially below that of microsatellite repeat units in mammals studied to date (range = 10(-2)-10(-5) per locus per generation). To obtain a more precise estimate of mutation rate for dinucleotide repeat motifs alone, we assayed 39 new dinucleotide repeat microsatellite loci in the mutation accumulation lines from our earlier study. Our estimate of mutation rate for a total of 49 dinucleotide repeats is 9.3 x 10(-6) per locus per generation, only slightly higher than the estimate from our earlier study. We also estimated the relative difference in microsatellite mutation rate among di-, tri-, and tetranucleotide repeats in the genome of D. melanogaster using a method based on population variation, and we found that tri- and tetranucleotide repeats mutate at rates 6.4 and 8.4 times slower than that of dinucleotide repeats, respectively. The slower mutation rates of tri- and tetranucleotide repeats appear to be associated with a relatively short repeat unit length of these repeat motifs in the genome of D. melanogaster. A positive correlation between repeat unit length and allelic variation suggests that mutation rate increases as the repeat unit lengths of microsatellites increase.      

Plasmodium vivax genetic diversity: microsatellite length matters

The Plasmodium vivax genome is very diverse but has a relatively low abundance of microsatellites. Leclerc et al. had shown that these di-nucleotide repeats have a low level of polymorphism, suggesting a recent bottleneck event in the evolutionary history of P. vivax. By contrast, in a recent paper, Imwong et al. show that there is a very high level of microsatellite diversity. The difference in these results is probably due to the set array lengths chosen by each group. Longer arrays are more diverse than are shorter ones because slippage mutations become exponentially more common with an increase in array length. These studies highlight the need to consider carefully the application and design of studies involving microsatellites.     

Silene tatarica microsatellites are frequently located in repetitive DNA

The genomic distribution of microsatellites can be explained by DNA slippage, slippage like processes and base substitutions. Nevertheless, microsatellites are also frequently associated with repetitive DNA, raising the question of the relative contributions of these processes to microsatellite genesis. We show that in Silene tatarica about 50% of the microsatellites isolated by an enrichment cloning protocol are associated with repetitive DNA. Based on the flanking sequences, we distinguished seven different classes of repetitive DNA. PCR primers designed for the flanking sequences of an individual clone amplified a heterogeneous family of repetitive DNA. Despite considerable variation in the flanking sequence (pi = 0.108), the microsatellite repeats did not show any evidence for decay. Rather, we observed the emergence of a new repeat type that probably arose by mutation and was spread by replication slippage. In fact, a complete repeat type switch could be observed among the analysed clones. We propose that the analysis of microsatellite sequences embedded in repetitive DNA provides a hitherto largely unexplored tool to study microsatellite evolution.     

Genome-wide analysis of conservation and divergence of microsatellites in rice

Studies on microsatellite distribution and divergence in related genomes contribute towards understanding of genome evolution in eukaryotes. Despite the availability of whole genome sequences of four rice genomes, occurrence and significance of microsatellites in the rice genome has remained a relatively unexplored area of research. We have aligned genomes of two rice subspecies i.e. indica and japonica to understand the trends of microsatellite conservation and divergence in the rice genome. Nearly 62% of the indica microsatellites were also found in the japonica genome. Occurrence of microsatellites showed a negative association with that of retrotransposons. Microsatellites repeat unit length and sequence showed direct influence on the microsatellite locus length. Further, microsatellite allele length was also influenced by the sequence characteristics of the neighbouring regions. CCG repeats were most conserved microsatellite sequences across the different syntenic regions in the two rice genomes and often showed association with CpG islands. Our study suggested that microsatellite distribution is not only governed by a balance between replication slippage and point mutations as proposed earlier, but also by the microsatellite motif sequence and characteristics of microsatellite neighbouring regions in the genome. Thus, this study is likely to prove an important reference for understanding the process of microsatellite evolution and dynamics in the two rice subspecies.     

A Phylogenetic Perspective on Sequence Evolution in Microsatellite Loci

We examined the evolution of the repeat regions of three noncoding microsatellite loci in 58 species of the Polistinae, a subfamily of wasps that diverged over 140 million years ago. A phylogenetic approach allows two new kinds of approaches to studying microsatellite evolution: character mapping and comparative analysis. The basic repeat structure of the loci was highly conserved, but was often punctuated with imperfections that appear to be phylogenetically informative. Repeat numbers evolved more rapidly than other changes in the repeat region. Changes in number of repeats among species seem consistent with the stepwise mutation model, which is based on slippage during replication as the main source of mutations. Changes in repeat numbers can occur even when there are very few tandem repeats but longer repeats, especially perfect repeats led to greater rates of evolutionary change. Species phylogenetically closer to the one from which we identified the loci had longer stretches of uninterrupted repeats and more different motifs, but not longer total repeat regions. The number of perfect repeats increased more often than it decreased. However, there was no evidence that some species have consistently greater numbers of repeats across loci than other species have, once ascertainment bias is eliminated. We also found no evidence for a population size effect posited by one form of the directionality hypothesis. Overall, phylogenetic variation in repeat regions can be explained by adding neutral evolution to what is already known about the mutation process. The life cycle of microsatellites appears to reflect a balance between growth by slippage and degradation by an essentially irreversible accumulation of imperfections. Received: 13 April 1999 / Accepted: 8 September 1999     

Instability of interstitial telomeric sequences in the human genome

The length variability of four human interstitial telomeric sequences (ITs) is described. Three of the ITs contain short telomeric stretches ranging between 53 and 84 bp and are localized in 21q22, 2q31, and 7q36; the fourth IT derives from the subtelomeric domain of chromosome 6p and contains a tract of a few hundred basepairs of exact and degenerate repeats. Using primers flanking the repeats, we amplified the genomic DNA from unrelated individuals and from family members, and we found that all the loci are polymorphic. At the 21q22 IT locus, two equally frequent alleles were found, while the number of alleles at the 2q31, 7q36, and 6pter IT loci was 8, 6, and 4, respectively. Sequence analysis revealed that in the three loci containing short ITs the alleles differ from one another for multiples of the hexanucleotide; it is likely that the mechanism leading to the polymorphism is DNA polymerase slippage. These loci were also unstable in gastric tumor cells characterized by microsatellite instability. At the 6pter IT locus, the four alleles range in length from about 500 to about 700 bp; this variability is probably due to unequal exchange or gene conversion. Our data indicate that stretches of exact internal telomeric repeats can be highly unstable, like microsatellites with shorter units, and that they can be useful polymorphic markers for linkage analysis, for forensic applications, and for the detection of genetic instability in tumors.     

Microsatellite mutation models: insights from a comparison of humans and chimpanzees

Using genomic data from homologous microsatellite loci of pure AC repeats in humans and chimpanzees, several models of microsatellite evolution are tested and compared using likelihood-ratio tests and the Akaike information criterion. A proportional-rate, linear-biased, one-phase model emerges as the best model. A focal length toward which the mutational and/or substitutional process is linearly biased is a crucial feature of microsatellite evolution. We find that two-phase models do not lead to a significantly better fit than their one-phase counterparts. The performance of models based on the fit of their stationary distributions to the empirical distribution of microsatellite lengths in the human genome is consistent with that based on the human-chimp comparison. Microsatellites interrupted by even a single point mutation exhibit a twofold decrease in their mutation rate when compared to pure AC repeats. In general, models that allow chimps to have a larger per-repeat unit slippage rate and/or a shorter focal length compared to humans give a better fit to the human-chimp data as well as the human genomic data.     

Characterization and isolation of novel microsatellites from the Drosophila dunni subgroup

We have isolated and characterized 77 novel microsatellites from two species, Drosophila dunni and Drosophila nigrodunni, which are closely related Caribbean-island endemics from the Drosophila cardini species group. These species are very distantly related to all other Drosophila from which microsatellites have previously been characterized. We find that the average length of microsatellites isolated in these species is quite small, with an overall mean length of 9.8 repeat units for dinucleotide microsatellites in the two study species. The nucleotide composition of dinucleotides differs between the two species: D. nigrodunni has a predominance of (AC/GT)n repeats, whereas D. dunni has equal numbers of (AC/GT)n and (AG/CT)n repeats. Tri- and tetranucleotide repeats are not abundant in either species. We assayed the variability of eight microsatellites in a closely related third species, Drosophila arawakana, using wild-caught individuals from the island of Guadeloupe. We found the microsatellites to be extremely variable in this population, with observed heterozygosities ranging from 0.541 to 0.889. DNA amplification trials suggest that these eight microsatellites are widely conserved across the D. cardini group, with five of the eight producing amplification products in every species tested. However, the loci are very poorly conserved over greater phylogenetic distances. DNA amplification of the microsatellite loci was unreliable in members of the closely related Drosophila quinaria, Drosophila calloptera, Drosophila guarani and Drosophila tripunctata species groups. Furthermore, these microsatellites could not be detected in the genome of Drosophila melanogaster, despite the conservation of microsatellite flanking regions at some loci. These data indicate that Drosophila microsatellite loci are quite short lived over evolutionary timescales relative to many other taxa.     

Multiple levels of single-strand slippage at cetacean tri- and tetranucleotide repeat microsatellite loci.

Between three and six tri- and tetranucleotide repeat microsatellite loci were analyzed in 3720 samples collected from four different species of baleen whales. Ten of the 18 species/locus combinations had imperfect allele arrays, i.e., some alleles differed in length by other than simple integer multiples of the basic repeat length. The estimate of the average number of alleles and heterozygosity was higher at loci with imperfect allele arrays relative to those with perfect allele arrays. Nucleotide sequences of 23 different alleles at one tetranucleotide repeat microsatellite locus in fin whales, Balaenoptera physalus, and humpback whales, Megaptera novaeangliae, revealed sequence changes including perfect repeats only, multiple repeats, and partial repeats. The relative rate of the latter two categories of mutation was estimated at 0.024 of the mutation rate involving perfect repeats only. It is hypothesized that single-strand slippage of partial repeats may provide a mechanism for counteracting the continuous expansion of microsatellite loci, which is the logical consequence of recent reports demonstrating directional mutations. Partial-repeat mutations introduce imperfections in the repeat array, which subsequently could reduce the rate of single-strand slippage. Limited computer simulations confirmed this predicted effect of partial-repeat mutations.     

Low abundance of Escherichia coli microsatellites is associated with an extremely low mutation rate

It is widely assumed that microsatellites are generated by replication slippage, a mutation process specific to repetitive DNA. Consistent with their high mutation rate, microsatellites are highly abundant in most eukaryotic genomes. In Escherichia coli, however, microsatellites are rare and short despite the fact that a high microsatellite mutation rate was described. We show that this high microsatellite instability depends on the presence of the F-plasmid. E. coli cells lacking the F-plasmid have extremely low microsatellite mutation rates. This result provides a possible explanation for the genome-wide low density of microsatellites in E. coli. Furthermore, we show that the F-plasmid induced microsatellite instability is independent of the mismatch repair pathway.