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.     

Microsatellite mutations during the polymerase chain reaction: mean field approximations and their applications

We develop a novel mathematical model for microsatellite mutations during polymerase chain reaction (PCR). Based on the model, we study the first- and second-order moments of the number of repeat units in a randomly chosen molecule after n PCR cycles and their corresponding mean field approximations. We give upper bounds for the approximation errors and show that the approximation errors are small when the mutation rate is low. Based on the theoretical results, we develop a moment estimation method to estimate the mutation rate per-repeat-unit per PCR cycle and the probability of expansion when mutations occur. Simulation studies show that the moment estimation method can accurately recover the true mutation rate and probability of expansion. Finally, the method is applied to experimental data from single-molecule PCR experiments.     

Microsatellite evolutionary rate and pattern in Schistocerca gregaria inferred from direct observation of germline mutations

Unravelling variation among taxonomic orders regarding the rate of evolution in microsatellites is crucial for evolutionary biology and population genetics research. The mean mutation rate of microsatellites tends to be lower in arthropods than in vertebrates, but data are scarce and mostly concern accumulation of mutations in model species. Based on parent–offspring segregations and a hierarchical Bayesian model, the mean rate of mutation in the orthopteran insect Schistocerca gregaria was estimated at 2.1e^?4 per generation per untranscribed dinucleotide locus. This is close to vertebrate estimates and one order of magnitude higher than estimates from species of other arthropod orders, such as Drosophila melanogaster and Daphnia pulex. We also found evidence of a directional bias towards expansions even for long alleles and exceptionally large ranges of allele sizes. Finally, at transcribed microsatellites, the mean rate of mutation was half the rate found at untranscribed loci and the mutational model deviated from that usually considered, with most mutations involving multistep changes that avoid disrupting the reading frame. Our direct estimates of mutation rate were discussed in the light of peculiar biological and genomic features of S. gregaria, including specificities in mismatch repair and the dependence of its activity to allele length. Shedding new light on the mutational dynamics of grasshopper microsatellites is of critical importance for a number of research fields. As an illustration, we showed how our findings improve microsatellite application in population genetics, by obtaining a more precise estimation of S. gregaria effective population size from a published data set based on the same microsatellites.     

Determining microsatellite genotyping reliability and mutation detection ability: an approach using small-pool PCR from sperm DNA

Microsatellite genotyping from trace DNA is now common in fields as diverse as medicine, forensics and wildlife genetics. Conversely, small-pool PCR (SP-PCR) has been used to investigate microsatellite mutation mechanisms in human DNA, but has had only limited application to non-human species. Trace DNA and SP-PCR studies share many challenges, including problems associated with allelic drop-out, false alleles and other PCR artefacts, and the need to reliably identify genuine alleles and/or mutations. We provide a framework for the validation of such studies without a multiple tube approach and demonstrate the utility of that approach with an analysis of microsatellite mutations in the tammar wallaby (Macropus eugenii). Specifically, we amplified three autosomal microsatellites from somatic DNA to characterise efficiency and reliability of PCR from low-template DNA. Reconstruction experiments determined our ability to discriminate mutations from parental alleles. We then developed rules to guide data interpretation. We estimated mutation rates in sperm DNA to range from 1.5 × 10⁻² to 2.2 × 10⁻³ mutations per locus per generation. Large multi-step mutations were observed, providing evidence for complex mutation processes at microsatellites and potentially violating key assumptions in the stepwise mutation model. Our data demonstrate the necessity of actively searching for large mutation events when investigating microsatellite evolution and highlight the need for a thorough understanding of microsatellite amplification characteristics before embarking on SP-PCR or trace DNA studies.     

Characterization of Demographic Expansions From Pairwise Comparisons of Linked Microsatellite Haplotypes

This work extends the methods of demographic inference based on the distribution of pairwise genetic differences between individuals (mismatch distribution) to the case of linked microsatellite data. Population genetics theory describes the distribution of mutations among a sample of genes under different demographic scenarios. However, the actual number of mutations can rarely be deduced from DNA polymorphisms. The inclusion of mutation models in theoretical predictions can improve the performance of statistical methods. We have developed a maximum-pseudolikelihood estimator for the parameters that characterize a demographic expansion for a series of linked loci evolving under a stepwise mutation model. Those loci would correspond to DNA polymorphisms of linked microsatellites (such as those found on the Y chromosome or the chloroplast genome). The proposed method was evaluated with simulated data sets and with a data set of chloroplast microsatellites that showed signal for demographic expansion in a previous study. The results show that inclusion of a mutational model in the analysis improves the estimates of the age of expansion in the case of older expansions.     

Empirical evaluation reveals best fit of a logistic mutation model for human Y-chromosomal microsatellites

The rate of microsatellite mutation is dependent upon both the allele length and the repeat motif, but the exact nature of this relationship is still unknown. We analyzed data on the inheritance of human Y-chromosomal microsatellites in father-son duos, taken from 24 published reports and comprising 15,285 directly observable meioses. At the six microsatellites analyzed (DYS19, DYS389I, DYS390, DYS391, DYS392, and DYS393), a total of 162 mutations were observed. For each locus, we employed a maximum-likelihood approach to evaluate one of several single-step mutation models on the basis of the data. For five of the six loci considered, a novel logistic mutation model was found to provide the best fit according to Akaike's information criterion. This implies that the mutation probability at the loci increases (nonlinearly) with allele length at a rate that differs between upward and downward mutations. For DYS392, the best fit was provided by a linear model in which upward and downward mutation probabilities increase equally with allele length. This is the first study to empirically compare different microsatellite mutation models in a locus-specific fashion.     

Characteristics and frequency of germline mutations at microsatellite loci from the human Y chromosome, as revealed by direct observation in father/son pairs

A number of applications of analysis of human Y-chromosome microsatellite loci to human evolution and forensic science require reliable estimates of the mutation rate and knowledge of the mutational mechanism. We therefore screened a total of 4,999 meioses from father/son pairs with confirmed paternity (probability >/=99. 9%) at 15 Y-chromosomal microsatellite loci and identified 14 mutations. The locus-specific mutation-rate estimates were 0-8. 58x10-3, and the average mutation rate estimates were 3.17x10-3 (95% confidence interval [CI] 1.89-4.94x10-3) across 8 tetranucleotide microsatellites and 2.80x10-3 (95% CI 1.72-4.27x10-3) across all 15 Y-chromosomal microsatellites studied. Our data show a mutational bias toward length increase, on the basis of observation of more repeat gains than losses (10:4). The data are in almost complete agreement with the stepwise-mutation model, with 13 single-repeat changes and 1 double-repeat change. Sequence analysis revealed that all mutations occurred in uninterrupted homogenous arrays of >/=11 repeats. We conclude that mutation rates and characteristics of human Y-chromosomal microsatellites are consistent with those of autosomal microsatellites. This indicates that the general mutational mechanism of microsatellites is independent of recombination.     

鲤鱼微卫星突变速率和模式(英文)

利用150个微卫星分子标记在F1代家系的基因型分析过程中,共有27600个等位基因从亲本向子代传递,其中在5个微卫星座位上检测到6个突变的等位基因。对突变的等位基因数目进行统计分析后得出:鲤鱼平均每个世代每个微卫星座位的突变速率为2.53×10-4。在发现突变的5个位点中,经测序发现,突变序列中插入1个以上的重复单元就导致了突变的发生。这些突变表明,鲤鱼的微卫星突变没有遵循严格的渐变突变模型(stepwise mutation model,SMM)。该文关于鲤鱼微卫星突变速率和模式的研究将会对统计鲤鱼有效群体的统计提供有效参数。     

Microsatellites can be misleading: an empirical and simulation study

Abstract. It has been long recognized that highly polymorphic genetic markers can lead to underestimation of divergence between populations when migration is low. Microsatellite loci, which are characterized by extremely high mutation rates, are particularly likely to be affected. Here, we report genetic differentiation estimates in a contact zone between two chromosome races of the common shrew ( Sorex araneus ), based on 10 autosomal microsatellites, a newly developed Y-chromosome microsatellite, and mitochondrial DNA. These results are compared to previous data on proteins and karyotypes. Estimates of genetic differentiation based on F - and R -statistics are much lower for autosomal microsatellites than for all other genetic markers. We show by simulations that this discrepancy stems mainly from the high mutation rate of microsatellite markers for F -statististics and from deviations from a single-step mutation model for R -statistics. The sex-linked genetic markers show that all gene exchange between races is mediated by females. The absence of male-mediated gene flow most likely results from male hybrid sterility.     

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.     

DNA dinucleotide evolution in humans: fitting theory to facts

We examine length distributions of approximately 6000 human dinucleotide microsatellite loci, representing chromosomes 1-22, from the GDB database. Under the stepwise mutation model, results from theory and simulation are compared with the empirical data. In both constant and expanding population scenarios, a simple single-step model with parameters chosen to account for the observed variance of microsatellite lengths produces results inconsistent with the observed heterozygosity and the dispersion of length skewness. Complicating the model by allowing a variable mutation rate accounts for the homozygosity, and introducing a small probability of a large mutation step accounts for the dispersion in skewnesses. We discuss these results in light of the long-term evolution of microsatellites.     

Moderately and Highly Polymorphic Microsatellites Provide Discordant Estimates of Population Divergence in Sockeye Salmon, Oncorhynchus nerka

Mutation rate can vary widely among microsatellite loci. This variation may cause discordant single-locus and multi-locus estimates of F_ST, the commonly used measure of population divergence. We use 16 microsatellite and five allozyme loci from 14 sockeye salmon populations to address two questions about the affect of mutation rate on estimates of F_ST: (1) does mutation rate influence F_ST estimates from all microsatellites to a similar degree relative to allozymes?; (2) does the influence of mutation rate on F_ST estimates from microsatellites vary with geographic scale in spatially structured populations? For question one we find that discordant estimates of F_ST among microsatellites as well as between the two marker classes are correlated with mean within-population heterozygosity (H_S) and thus are likely due to differences in mutation rate. Highly polymorphic microsatellites (H_S > 0.84) provide significantly lower estimates of F_ST than moderately polymorphic microsatellites and allozymes (H_S < 0.60). Estimates of F_ST from binned allele frequency data and R_ST provide more accurate measures of population divergence for highly polymorphic but not for moderately polymorphic microsatellites. We conclude it is more important to pool loci of like H_S rather than marker class when estimating F_ST. For question two we find the F_ST values for moderately and highly polymorphic loci, while significantly different, are positively correlated for geographically proximate but not geographically distant population pairs. These results are consistent with expectations from the equilibrium approximation of Wright's infinite island model and confirm that the influence of mutation on estimates of F_ST can vary in spatially structured populations presumably because the rate of migration varies inversely with geographic scale.     

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.     

Microsatellite variation in North American populations of Drosophila melanogaster.

Computer database searching for microsatellites can be particularly effective for organisms like Drosophila melanogaster for which there are extensive sequence data. Here we demonstrate that 17 out of 18 such microsatellites are also highly polymorphic in natural populations of Drosophila, and that this variation is easily scorable with PCR followed by electrophoresis on high-resolution agarose. This form of variation is likely to be of great value in studies of the genomic distribution of polymorphism, population structure, the relation between intraspecific polymorphism and interspecific divergence and the mutation rate and pattern of mutations of microsatellites. In this preliminary survey of 15 lines, we find that the variance in repeat count is most strongly correlated with the maximum count, that perfect repeats are significantly more variable than imperfect repeats and that repeats which are split by an imperfection have unexpectedly low variance given the size of the perfectly repeated portion.     

Complex mutations in a high proportion of microsatellite loci from the protozoan parasite Plasmodium falciparum

Microsatellite loci are generally assumed to evolve via a stepwise mutational process and a battery of statistical techniques has been developed in recent years based on this or related mutation models. It is therefore important to investigate the appropriateness of these models in a wide variety of taxa. We used two approaches to examine mutation patterns in the malaria parasite Plasmodium falciparum: (i) we examined sequence variation at 12 tri-nucleotide repeat loci; and (ii) we analysed patterns of repeat structure and heterozygosity at 114 loci using data from 12 laboratory parasite lines. The sequencing study revealed complex patterns of mutation in five of the 12 loci studied. Alleles at two loci contain indels of 24 bp and 57 bp in flanking regions, while in the other three loci, blocks of imperfect microsatellites appear to be duplicated or inserted; these loci essentially consist of minisatellite repeats, with each repeat unit containing four to eight microsatellites. The survey of heterozygosity revealed a positive relationship between repeat number and microsatellite variability for both di- and trinucleotides, indicating a higher mutation rate in loci with longer repeat arrays. Comparisons of levels of variation in different repeat types indicate that the mutation rate of dinucleotide-bearing loci is 1.6-2.1 times faster than trinucleotides, consistent with the lower mean number of repeats in trinucleotide-bearing loci. However, despite the evidence that microsatellite arrays themselves are evolving in a manner consistent with stepwise mutation model in P. falciparum, the high frequency of complex mutations precludes the use of analytical tools based on this mutation model for many microsatellite-bearing loci in this protozoan. The results call into question the generality of models based on stepwise mutation for analysing microsatellite data, but also demonstrate the ease with which loci that violate model assumptions can be detected using minimal sequencing effort.     

Size homoplasy and mutational processes of interrupted microsatellites in two bee species, Apis mellifera and Bombus terrestris (Apidae)

Similar microsatellite electromorphs (PCR products of the same size) can arise from independent mutational events. Such alleles are not identical by descent. This phenomenon, termed size homoplasy, was studied by sequencing electromorphs of two microsatellite loci in which the stretch of basic repeats is interrupted by different short (1-2 bp) DNA motifs. The number and position of these interruptions were established for electromorphs from closely and distantly related populations of honeybees and bumblebees. No sequence difference was found when electromorphs came from the same subspecies or from closely related subspecies, suggesting that they were probably identical by descent. In contrast, sequence differences were often detected in distantly related subspecies, showing that size homoplasy frequently occurs at this level of population differentiation. Size homoplasy is increased by limits to free length variation of alleles, a phenomenon that seems to act on interrupted microsatellites when comparing distantly related taxa, that is, honeybee subspecies from different evolutionary lineages. Electromorph sequences suggest that, within the scope of these limits, large mutation events have occurred frequently at both interrupted loci studied. In good agreement with the molecular data, computations based on the observed heterozygosity and number of electromorphs and simulation studies showed that neither locus fits the one-step stepwise mutant model (SMM). We speculate that interrupted microsatellites in general could be characterized by a higher variance in repeat number and consequently a lower homoplasy rate than pure ones. Hence, interrupted microsatellites should be most appropriate for investigating population differentiation and evolutionary relationship between relatively distant populations.      

A maximum-likelihood approach to fitting equilibrium models of microsatellite evolution

Here, we develop a new approach to Markov chain modeling of microsatellite evolution through polymerase slippage and introduce new models: a "constant-slippage-rate" model, in which there is no dependence of slippage rate on microsatellite length, as envisaged by Moran; and a "linear-with-constant" model, in which slippage rate increases linearly with microsatellite length, but the line of best fit is not constrained to go through the origin. We show how these and a linear no-constant model can be fitted to data hierarchically using maximum likelihood. This has advantages over previous methods in allowing statistical comparisons between models. When applied to a previously analyzed data set, the method allowed us to statistically establish that slippage rate increases with microsatellite length for dinucleotide microsatellites in humans, mice, and fruit flies, and suggested that no slippage occurs in very short microsatellites of one to four repeats. The suggestion that slippage rates are zero or close to zero for very short microsatellites of one to four repeats has important implications for understanding the mechanism of polymerase slippage.     

Estimation of the mutation rate during error-prone polymerase chain reaction.

Error-prone polymerase chain reaction (PCR) is widely used to introduce point mutations during in vitro evolution experiments. Accurate estimation of the mutation rate during error-prone PCR is important in studying the diversity of error-prone PCR product. Although many methods for estimating the mutation rate during PCR are available, all the existing methods depend on the assumption that the mutation rate is low and mutations occur at different places whenever they occur. The available methods may not be applicable to estimate the mutation rate during error-prone PCR. We develop a mathematical model for error-prone PCR and present methods to estimate the mutation rate during error-prone PCR without assuming low mutation rate. We also develop a computer program to simulate error-prone PCR. Using the program, we compare the newly developed methods with two other methods. We show that when the mutation rate is relatively low (< 10(-3) per base per PCR cycle), the newly developed methods give roughly the same results as previous methods. When the mutation rate is relatively high (> 5 x 10(-3) per base per PCR cycle, the mutation rate for most error-prone PCR experiments), the previous methods underestimate the mutation rate and the newly developed methods approximate the true mutation rate.     

Evidence for Nonindependent Evolution of Adjacent Microsatellites in the Human Genome

Microsatellites are short tandem repeats that evolve predominantly through a stepwise mutation model. Despite intensive study, many aspects of their evolution remain unresolved, particularly the question of how compound microsatellites containing two different motifs evolve. Previous work described profound asymmetries in the probability that any given second motif lies either 3′ or 5′ of an AC repeat tract. Here we confirm and extend this analysis to examine the length dependence of these asymmetries. We then use the differences in length between homologous human and chimpanzee microsatellites as a surrogate measure of the slippage-based mutation rate to explore factors that influence this process. We find that the dominant predictor of mutation rate is the length of the tract being considered, which is a stronger predictor than the length of the two tracts combined, but other factors also have a significant impact, including the length of the second tract and which of the two tracts lies upstream. We conclude that compound microsatellites rarely arise through random point mutations generating a second motif within a previously pure tract. Instead, our analyses point toward a model in which poorly understood mutation biases, probably affecting both slippage and point mutations and often showing 3′-5′ polarity, promote the formation of compound microsatellites. The result is convergent evolution. We suggest that, although their exact nature remains unclear, these biases are likely attributable to structural features, such as the propensity of AC tracts to form Z-DNA.     

A coalescent approach to the polymerase chain reaction.

A versatile algorithm is developed to model PCR on a computer. The method is based on a modification of the coalescent process and provides a general framework to analyse data from PCR. It allows for incorporation of the dynamics of the replication process as described in terms of the number of starting template molecules and cycle-dependent PCR efficiency. The simulation method generates, as a first step, the genealogy of a set of sequences sampled from a final PCR product. In a second step a mutation process is superimposed and the resulting data set is analysed. The efficiency of our algorithm enables us to get reliable approximations of various sample distributions. We demonstrate the relevance of our method with two applications: maximum likelihood estimation of the error rate in PCR and a test of homogeneity of the template.