Characteristics of loci and individuals are associated with germline microsatellite mutation rates in lesser kestrels (Falco naumanni)

Although microsatellites are one of the most popular tools in genetic studies, their mutational dynamics and evolution remain unclear. Here, we apply extensive pedigree genotyping to identify and analyze the patterns and factors associated with de novo germline mutations across nine microsatellite loci in a wild population of lesser kestrels (Falco naumanni). A total of 10 germline mutations events were unambiguously identified in four loci, yielding an average mutation rate of 2.96x10(-3). Across loci, mutation rate was positively correlated with locus variability and average allele size. Mutations were primarily compatible with a stepwise mutation model, although not exclusively involved single-step changes. Unexpectedly, we found an excess of maternally transmitted mutations (male-to-female ratio of 0.1). One of the analyzed loci (Fn2.14) resulted hypermutable (mutation rate=0.87%). This locus showed a size-dependent mutation bias, with longer alleles displaying deletions or additions of a small number of repeat than shorter alleles. Mutation probability at Fn2.14 was higher for females and increased with parental (maternal) age but was not associated with individual physical condition, multilocus heterozygosity, allele length or allele span. Overall, our results do not support the male-biased mutation rate described in other organisms and suggest that mutation dynamics at microsatellite loci are a complex process which requires further research.     

MICROSATELLITE EVOLUTION IN VERTEBRATES: INFERENCE FROM AC DINUCLEOTIDE REPEATS

Abstract We analyze published data from 592 AC microsatellite loci from 98 species in five vertebrate classes including fish, reptiles, amphibians, birds, and mammals. We use these data to address nine major questions about microsatellite evolution. First, we find that larger genomes do not have more microsatellite loci and therefore reject the hypothesis that microsatellites function primarily to package DNA into chromosomes. Second, we confirm that microsatellite loci are relatively rare in avian genomes, but reject the hypothesis that this is due to physical constraints imposed by flight. Third, we find that microsatellite variation differs among species within classes, possibly relating to population dynamics. Fourth, we reject the hypothesis that microsatellite structure (length, number of alleles, allele dispersion, range in allele sizes) differs between poikilotherms and homeotherms. The difference is found only in fish, which have longer microsatellites and more alleles than the other classes. Fifth, we find that the range in microsatellite allele size at a locus is largely due to the number of alleles and secondarily to allele dispersion. Sixth, length is a major factor influencing mutation rate. Seventh, there is a directional mutation toward an increase in microsatellite length. Eighth, at the species level, microsatellite and allozyme heterozygosity covary and therefore inferences based on large-scale studies of allozyme variation may also reflect microsatellite genetic diversity. Finally, published microsatellite loci (isolated using conventional hybridization methods) provide a biased estimate of the actual mean repeat length of microsatellites in the genome.     

雌核发育银鲫微卫星突变速率和突变模式

银鲫是天然雌核发育的三倍体两性型种群, 因其遗传背景和生殖方式的特殊性, 已经成为研究单性和多倍体脊椎动物进化遗传学的理想模式鱼类。利用33个微卫星序列对雌核发育银鲫的突变速率和突变模式进行研究, 结果表明: (1) 22个子代个体中, 检测到1尾个体在15个微卫星位点具有18个突变等位基因; (2) 每个微卫星位点每代总体平均突变率是1.16×10-2, 95%置信区间是6.87×10-3~1.83×10-2, 与其他鱼类相比, 雌核发育银鲫的突变率明显偏高, 这与天然雌核发育鱼类处在单性生殖和两性生殖的过度阶段有密切关系; (3) 具有突变等位基因的13个位点的重复单元数目都在10次以上, 11个复合型微卫星位点的突变率(1.31×10-2)与21个完美型位点(1.00×10-2)的突变率没有明显差异(P = 0.67), 微卫星突变率受到重复单元数目的影响, 然而与重复结构类型和侧翼序列GC碱基含量无相关性; (4) 序列分析表明, 雌核发育银鲫的微卫星突变模式并不严格遵守逐步突变模型。     

Mutation and evolution of microsatellite loci in Neurospora

The patterns of mutation and evolution at 13 microsatellite loci were studied in the filamentous fungal genus Neurospora. First, a detailed investigation was performed on five microsatellite loci by sequencing each microsatellite, together with its nonrepetitive flanking regions, from a set of 147 individuals from eight species of Neurospora. To elucidate the genealogical relationships among microsatellite alleles, repeat number was mapped onto trees constructed from flanking-sequence data. This approach allowed the potentially convergent microsatellite mutations to be placed in the evolutionary context of the less rapidly evolving flanking regions, revealing the complexities of the mutational processes that have generated the allelic diversity conventionally assessed in population genetic studies. In addition to changes in repeat number, frequent substitution mutations within the microsatellites were detected, as were substitutions and insertion/deletions within the flanking regions. By comparing microsatellite and flanking-sequence divergence, clear evidence of interspecific allele length homoplasy and microsatellite mutational saturation was observed, suggesting that these loci are not appropriate for inferring phylogenetic relationships among species. In contrast, little evidence of intraspecific mutational saturation was observed, confirming the utility of these loci for population-level analyses. Frequency distributions of alleles within species were generally consistent with the stepwise mutational model. By comparing variation within species at the microsatellites and the flanking-sequence, estimated microsatellite mutation rates were approximately 2500 times greater than mutation rates of flanking DNA and were consistent with estimates from yeast and fruit flies. A positive relationship between repeat number and variance in repeat number was significant across three genealogical depths, suggesting that longer microsatellite alleles are more mutable than shorter alleles. To test if the observed patterns of microsatellite variation and mutation could be generalized, an additional eight microsatellite loci were characterized and sequenced from a subset of the same Neurospora individuals.     

Genetic variation and de novo mutations in the parthenogenetic Caucasian rock lizard Darevskia unisexualis

Unisexual all-female lizards of the genus Darevskia that are well adapted to various habitats are known to reproduce normally by true parthenogenesis. Although they consist of unisexual lineages and lack effective genetic recombination, they are characterized by some level of genetic polymorphism. To reveal the mutational contribution to overall genetic variability, the most straightforward and conclusive way is the direct detection of mutation events in pedigree genotyping. Earlier we selected from genomic library of D. unisexualis two polymorphic microsatellite containing loci Du281 and Du215. In this study, these two loci were analyzed to detect possible de novo mutations in 168 parthenogenetic offspring of 49 D. unisexualis mothers and in 147 offspring of 50 D. armeniaca mothers. No mutant alleles were detected in D. armeniaca offspring at both loci, and in D. unisexualis offspring at the Du215 locus. There were a total of seven mutational events in the germ lines of four of the 49 D. unisexualis mothers at the Du281 locus, yielding the mutation rate of 0.1428 events per germ line tissue. Sequencing of the mutant alleles has shown that most mutations occur via deletion or insertion of single microsatellite repeat being identical in all offspring of the family. This indicates that such mutations emerge at the early stages of embryogenesis. In this study we characterized single highly unstable (GATA)(n) containing locus in parthenogenetic lizard species D. unisexualis. Besides, we characterized various types of mutant alleles of this locus found in the D. unisexualis offspring of the first generation. Our data has shown that microsatellite mutations at highly unstable loci can make a significant contribution to population variability of parthenogenetic lizards.     

A high incidence of clustered microsatellite mutations revealed by parent-offspring analysis in the African freshwater snail, <Emphasis Type="BoldItalic">Bulinus forskalii</Emphasis> (Gastropoda,Pulmonata)

Genotyping of 11 microsatellites in 432 offspring from 28 families of the hermaphroditic, freshwater snail Bulinus forskalii detected 10 de novo mutant alleles. This gave an estimated mutation rate of 1.1 × 10^–3 per locus per gamete per generation. There was a trend towards repeat length expansion and, unlike most studies, multi-step mutations predominated, suggesting that the microsatellite mutation process does not conform to a strict stepwise mutation model. Interestingly, the ten mutant alleles appear to have arisen from only six independent germline mutation events within the microsatellite array, with seven of them residing in three mutational clusters. Our results extend observations of clustered microsatellite mutations to another taxonomic group and type of mating system, self-fertile gastropods, and provide compelling evidence of premeiotic germline mutations, a phenomenon that could greatly impact upon our understanding of mutation dynamics but which has received little attention.     

Microsatellite mutations in the germline: implications for evolutionary inference

Microsatellite DNA sequences mutate at rates several orders of magnitude higher than that of the bulk of DNA. Such high rates mean that spontaneous mutations that form new-length variants can realistically be seen in pedigree analysis. Data on observed mutation events from various organisms are now accumulating, allowing inferences on DNA sequence evolution to be made through an unusually direct approach. Here I discuss and integrate microsatellite mutation data in an evolutionary context. A striking feature of the mutation process is that it seems highly heterogeneous, with distinct differences between species, repeat types, loci and alleles. Age and sex also affect the mutation rate. Within genomes at equilibrium, the microsatellite-length distribution is a delicate balance between biased mutation processes and point mutations acting towards the decay of repetitive DNA. Indeed, simple repeats do not evolve simply.     

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.     

Vntr Allele Frequency Distributions under the Stepwise Mutation Model: A Computer Simulation Approach

Variable numbers of tandem repeats (VNTRs) are a class of highly informative and widely dispersed genetic markers. Despite their wide application in biological science, little is known about their mutational mechanisms or population dynamics. The objective of this work was to investigate four summary measures of VNTR allele frequency distributions: number of alleles, number of modes, range in allele size and heterozygosity, using computer simulations of the one-step stepwise mutation model (SMM). We estimated these measures and their probability distributions for a wide range of mutation rates and compared the simulation results with predictions from analytical formulations of the one-step SMM. The average heterozygosity from the simulations agreed with the analytical expectation under the SMM. The average number of alleles, however, was larger in the simulations than the analytical expectation of the SMM. We then compared our simulation expectations with actual data reported in the literature. We used the sample size and observed heterozygosity to determine the expected value, 5th and 95th percentiles for the other three summary measures, allelic size range, number of modes and number of alleles. The loci analyzed were classified into three groups based on the size of the repeat unit: microsatellites (1-2 base pair (bp) repeat unit), short tandem repeats [(STR) 3-5 bp repeat unit], and minisatellites (15-70 bp repeat unit). In general, STR loci were most similar to the simulation results under the SMM for the three summary measures (number of alleles, number of modes and range in allele size), followed by the microsatellite loci and then by the minisatellite loci, which showed deviations in the direction of the infinite allele model (IAM). Based on these differences, we hypothesize that these three classes of loci are subject to different mutational forces.     

Joint inference of microsatellite mutation models, population history and genealogies using transdimensional Markov Chain Monte Carlo

We provide a framework for Bayesian coalescent inference from microsatellite data that enables inference of population history parameters averaged over microsatellite mutation models. To achieve this we first implemented a rich family of microsatellite mutation models and related components in the software package BEAST. BEAST is a powerful tool that performs Bayesian MCMC analysis on molecular data to make coalescent and evolutionary inferences. Our implementation permits the application of existing nonparametric methods to microsatellite data. The implemented microsatellite models are based on the replication slippage mechanism and focus on three properties of microsatellite mutation: length dependency of mutation rate, mutational bias toward expansion or contraction, and number of repeat units changed in a single mutation event. We develop a new model that facilitates microsatellite model averaging and Bayesian model selection by transdimensional MCMC. With Bayesian model averaging, the posterior distributions of population history parameters are integrated across a set of microsatellite models and thus account for model uncertainty. Simulated data are used to evaluate our method in terms of accuracy and precision of estimation and also identification of the true mutation model. Finally we apply our method to a red colobus monkey data set as an example.     

Escherichia coli DNA polymerase IV contributes to spontaneous mutagenesis at coding sequences but not microsatellite alleles

Slipped strand mispairing during DNA synthesis is one proposed mechanism for microsatellite or short tandem repeat (STR) mutation. However, the DNA polymerase(s) responsible for STR mutagenesis have not been determined. In this study, we investigated the effect of the Escherichia colidinB gene product (Pol IV) on mononucleotide and dinucleotide repeat stability, using an HSV-tk gene episomal reporter system for microsatellite mutations. For the control vector (HSV-tk gene only) we observed a statistically significant 3.5-fold lower median mutation frequency in dinB(-) than dinB(+) cells (p<0.001, Wilcoxon Mann Whitney Test). For vectors containing an in-frame mononucleotide allele ([G/C](10)) or either of two dinucleotide alleles ([GT/CA](10) and [TC/AG](11)) we observed no statistically significant difference in the overall HSV-tk mutation frequency observed between dinB(+) and dinB(-) strains. To determine if a mutational bias exists for mutations made by Pol IV, mutational spectra were generated for each STR vector and strain. No statistically significant differences between strains were observed for either the proportion of mutational events at the STR or STR specificity among the three vectors. However, the specificity of mutational events at the STR alleles in each strain varied in a statistically significant manner as a consequence of microsatellite sequence. Our results indicate that while Pol IV contributes to spontaneous mutations within the HSV-tk coding sequence, Pol IV does not play a significant role in spontaneous mutagenesis at [G/C](10), [GT/CA](10), or [TC/AG](11) microsatellite alleles. Our data demonstrate that in a wild type genetic background, the major factor influencing microsatellite mutagenesis is the allelic sequence composition.     

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.     

Microsatellite evolution at two hypervariable loci revealed by extensive avian pedigrees

Genealogies generated through a long-term study of superb fairy-wrens (Malurus cyaneus) were used to investigate mutation within two hypervariable microsatellite loci. Of 3,230 meioses examined at the tetranucleotide locus (Mcy micro 8), 45 mutations were identified, giving a mutation rate of 1.4%. At the dinucleotide locus (Mcy micro 4) 30 mutations were recorded from 2,750 meioses giving a mutation rate of 1.1%. Mutations at both loci primarily (80%; 60/75) involved the loss or gain of a single repeat unit. Unlike previous studies, there was no significant bias toward additions over deletions. The mutation rate at Mcy micro 8 increased with allele size, and very long alleles (>70 repeats) mutated at a rate of almost 20%. The length of the mutating allele and allele span, however, were strongly correlated so it was not possible to isolate the causative factor. Allele size did not appear to affect mutation rate at Mcy micro 4, but the repeat number was considerably lower at this locus. The gender of the mutating parent was significant only at Mcy micro 8, where mutations occurred more frequently in maternal alleles. However, at both loci we found that alleles inherited from the mother were on average larger than those from the father, and this in part drove the higher mutation rate among maternally inherited alleles at Mcy micro 8.     

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 Polymorphism in Natural Populations of the Wild Plant Arabidopsis Thaliana

Variation in repeat number at 20 microsatellite loci of Arabidopsis thaliana was studied in a worldwide sample of 42 ecotypes to investigate the pattern and level of polymorphism in repetitive sequences in natural plant populations. There is a substantial amount of variation at microsatellite loci despite the selfing nature of this plant species. The average gene diversity was 0.794 and the average number of alleles per locus was 10.6. The distribution of alleles was centered around the mean of repeat number at most loci, but could not be regarded as normal. There was a significantly positive correlation between the number of repeats and the amount of variation. For most loci, the observed number of alleles was between the expected values of the infinite allele and stepwise mutation models. The two models were rejected by the sign test. Linkage disequilibrium was detected in 12.1% of the pairwise comparisons between loci. In phylogenetic tree, there was no association between ecotype and geographic origin. This result is consistent with the recent expansion of A. thaliana throughout the world.     

Simultaneous estimation of all the parameters of a stepwise mutation model.

Minisatellite and microsatellite are short tandemly repetitive sequences dispersed in eukaryotic genomes, many of which are highly polymorphic due to copy number variation of the repeats. Because mutation changes copy numbers of the repeat sequences in a generalized stepwise fashion, stepwise mutation models are widely used for studying the dynamics of these loci. We propose a minimum chi-square (MCS) method for simultaneous estimation of all the parameters in a stepwise mutation model and the ancestral allelic type of a sample. The MCS estimator requires knowing the mean number of alleles of a certain size in a sample, which can be estimated using Monte Carlo samples generated by a coalescent algorithm. The method is applied to samples of seven (CA)n repeat loci from eight human populations and one chimpanzee population. The estimated values of parameters suggest that there is a general tendency for microsatellite alleles to expand in size, because (1) each mutation has a slight tendency to cause size increase and (2) the mean size increase is larger than the mean size decrease for a mutation. Our estimates also suggest that most of these CA-repeat loci evolve according to multistep mutation models rather than single-step mutation models. We also introduced several quantities for measuring the quality of the estimation of ancestral allelic type, and it appears that the majority of the estimated ancestral allelic types are reasonably accurate. Implications of our analysis and potential extensions of the method are discussed.SINCE the discovery that a large number of loci with tandemly repeated sequences in human and many eukaryote species are highly polymorphic because of copy number variation of the repeats in different individuals (Jeffreys 1985; Litt and Luty 1989; Weber and May 1989), allele size data from such loci are rapidly becoming the dominant source of genetic markers for genome mapping, forensic testing, and population studies. Loci with repeat sequences longer than 5 bp are generally referred to as minisatellite or variable number tandem repeat loci, and those with repeat sequences between 2 to 5 bp are referred to as microsatellite or short tandem repeat loci (Tautz 1993). Because mutations change the copy number of such loci in a stepwise fashion, rapid accumulation of population samples from minisatellite and microsatellite loci has resurrected the interest of the stepwise mutation model (SMM), which was popular in the 1970s.     

Single-molecule analysis of the hypermutable tetranucleotide repeat locus D21S1245 through sperm genotyping: a heterogeneous pattern of mutation but no clear male age effect

Single molecule genotyping of the hypermutable microsatellite locus D21S1245 was used for studying how the rate and pattern of mutation varied between alleles and different age groups. In total, 203 mutation events were scored from the genotyping of DNA corresponding to an estimated 8623 sperm cells from eight different men. Allele-specific mutation rates ranged from 0.007 to 0.052, a heterogeneity related in part to variation in the mutation rate among three allelic lineages identified after allele sequencing. Alleles from these lineages differed in the overall repeat structure of this complex microsatellite locus. Also, the pattern of mutation varied between lineages in that they differed in the relative proportions of expansion and contraction mutations. Surprisingly, a group of four men aged 18-23 years showed a higher mean mutation rate than a group of four men aged 48-56 years. To some extent this age difference can probably be explained by a bias in the distribution of alleles from the three allelic lineages among the age groups. However, the absence of a clear male age effect is at odds with the idea of an increasing male mutation rate with age, which is thought to arise from the continuous replication of germline cells throughout adulthood.     

Evidence for complex mutations at microsatellite loci in Drosophila.

Fifteen lines each of Drosophila melanogaster, D. simulans, and D. sechellia were scored for 19 microsatellite loci. One to four alleles of each locus in each species were sequenced, and microsatellite variability was compared with sequence structure. Only 7 loci had their size variation among species consistent with the occurrence of strictly stepwise mutations in the repeat array, the others showing extensive variability in the flanking region compared to that within the microsatellite itself. Polymorphisms apparently resulting from complex nonstepwise mutations involving the microsatellite were also observed, both within and between species. Maximum number of perfect repeats and variance of repeat count were found to be strongly correlated in microsatellites showing an apparently stepwise mutation pattern. These data indicate that many microsatellite mutation events are more complex than represented even by generalized stepwise mutation models. Care should therefore be taken in inferring population or phylogenetic relationships from microsatellite size data alone. The analysis also indicates, however, that evaluation of sequence structure may allow selection of microsatellites that more closely match the assumptions of stepwise models.     

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.     

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.