Dynamics of microsatellite divergence under stepwise mutation and proportional slippage/point mutation models

Ancient and Modern Genomes Reveal Microsatellites Maintain a Dynamic Equilibrium Through Deep Time

Microsatellites are widely used in population genetics, but their evolutionary dynamics remain poorly understood. It is unclear whether microsatellite loci drift in length over time. This is important because the mutation processes that underlie these important genetic markers are central to the evolutionary models that employ microsatellites. We identify more than 27 million microsatellites using a novel and unique dataset of modern and ancient Adélie penguin genomes along with data from 63 published chordate genomes. We investigate microsatellite evolutionary dynamics over 2 timescales: one based on Adélie penguin samples dating to ∼46.5 ka and the other dating to the diversification of chordates aged more than 500 Ma. We show that the process of microsatellite allele length evolution is at dynamic equilibrium; while there is length polymorphism among individuals, the length distribution for a given locus remains stable. Many microsatellites persist over very long timescales, particularly in exons and regulatory sequences. These often retain length variability, suggesting that they may play a role in maintaining phenotypic variation within populations.

Considerations on the application of a mutation model for Y-STR interpretation

If Y-STR profiling is to be more effective in criminal casework, the methods used to evaluate evidential weight require improvement. Many forensic scientists assign an evidential weight by estimating the number of times a Y-STR profile obtained from a questioned sample has been observed in YHRD datasets. More sophisticated models have been suggested but not yet implemented into routine casework, e.g. Andersen & Balding [1]. Mutation is inherent to STR meiosis (or inheritance) and is encountered in practice. We evaluated a mutation model that can be incorporated into a method for assigning evidential weight to Y-STR profiles, an essential part of bringing any method into practice. Since an important part of implementation to casework is communication, the article is written in an accessible format for practitioners as well as statisticians. The mutation component within the MUTEA model by Willems et al. [2] incorporates the potential for multistep mutations and a tendency for alleles to revert towards a central length, reflecting observed mutation data, e.g. [3]. We have estimated the parameters in this model and in a simplified symmetric version of this model, using sequence data from father/son pairs [4] and deep-rooted pedigrees [5]. Both datasets contain multistep mutations, which may have an effect on models based on simulations [1]. We introduce Beta-Binomial and Beta-Geometric conjugate analyses for estimating rate and step parameters for the mutation models presented here, which require only summations and multiplications. We proved mathematically that the parameters can be estimated independently. We show the importance of reporting the variability of the parameters and not only a point estimate. The parameters can be easily incorporated into statistical models, and updated sequentially as more data becomes available. We recommend fuller publication of data to enable the development and evaluation of a wider range of mutation models.

The complex Y-chromosomal history of gorillas

Studies of the evolutionary relationships among gorilla populations using autosomal and mitochondrial sequences suggest that male-mediated gene flow may have been important in the past, but data on the Y-chromosomal relationships among the gorilla subspecies are limited. Here, we genotyped blood and noninvasively collected fecal samples from 12 captives and 257 wild male gorillas of known origin representing all four subspecies (Gorilla gorilla gorilla, G. g. diehli, G. beringei beringei, and G. b. graueri) at 10 Y-linked microsatellite loci resulting in 102 unique Y-haplotypes for 224 individuals. We found that western lowland gorilla (G. g. gorilla) haplotypes were consistently more diverse than any other subspecies for all measures of diversity and comprised several genetically distinct groups. However, these did not correspond to geographical proximity and some closely related haplotypes were found several hundred kilometers apart. Similarly, our broad sampling of eastern gorillas revealed that mountain (G. b. beringei) and Grauer's (G. b. graueri) gorilla Y-chromosomal haplotypes did not form distinct clusters. These observations suggest structure in the ancestral population with subsequent mixing of differentiated haplotypes by male dispersal for western lowland gorillas, and postisolation migration or incomplete lineage sorting due to short divergence times for eastern gorillas.

Origin and dispersion pathways of guava in the Galapagos Islands inferred through genetics and historical records

Guava (Psidium guajava) is an aggressive invasive plant in the Galapagos Islands. Determining its provenance and genetic diversity could explain its adaptability and spread, and how this relates to past human activities. With this purpose, we analyzed 11 SSR markers in guava individuals from Isabela, Santa Cruz, San Cristobal, and Floreana islands in the Galapagos, as well as from mainland Ecuador. The mainland guava population appeared genetically differentiated from the Galapagos populations, with higher genetic diversity levels found in the former. We consistently found that the Central Highlands region of mainland Ecuador is one of the most likely origins of the Galapagos populations. Moreover, the guavas from Isabela and Floreana show a potential genetic input from southern mainland Ecuador, while the population from San Cristobal would be linked to the coastal mainland regions. Interestingly, the proposed origins for the Galapagos guava coincide with the first human settlings of the archipelago. Through approximate Bayesian computation, we propose a model where San Cristobal was the first island to be colonized by guava from the mainland, and then, it would have spread to Floreana and finally to Santa Cruz; Isabela would have been seeded from Floreana. An independent trajectory could also have contributed to the invasion of Floreana and Isabela. The pathway shown in our model agrees with the human colonization history of the different islands in the Galapagos. Our model, in conjunction with the clustering patterns of the individuals (based on genetic distances), suggests that guava introduction history in the Galapagos archipelago was driven by either a single event or a series of introduction events in rapid succession. We thus show that genetic analyses supported by historical sources can be used to track the arrival and spread of invasive species in novel habitats and the potential role of human activities in such processes.

Distribution patterns and variation analysis of simple sequence repeats in different genomic regions of bovid genomes

As the first examination of distribution, guanine-cytosine (GC) pattern, and variation analysis of microsatellites (SSRs) in different genomic regions of six bovid species, SSRs displayed nonrandomly distribution in different regions. SSR abundances are much higher in the introns, transposable elements (TEs), and intergenic regions compared to the 3′-untranslated regions (3′UTRs), 5′UTRs and coding regions. Trinucleotide perfect SSRs (P-SSRs) were the most frequent in the coding regions, whereas, mononucleotide P-SSRs were the most in the introns, 3′UTRs, TEs, and intergenic regions. Trifold P-SSRs had more GC-contents in the 5′UTRs and coding regions than that in the introns, 3′UTRs, TEs, and intergenic regions, whereas mononucleotide P-SSRs had the least GC-contents in all genomic regions. The repeat copy numbers (RCN) of the same mono- to hexanucleotide P-SSRs showed significantly different distributions in different regions (P < 0.01). Except for the coding regions, mononucleotide P-SSRs had the most RCNs, followed by the pattern: di- > tri- > tetra- > penta- > hexanucleotide P-SSRs in the same regions. The analysis of coefficient of variability (CV) of SSRs showed that the CV variations of RCN of the same mono- to hexanucleotide SSRs were relative higher in the intronic and intergenic regions, followed by the CV variation of RCN in the TEs, and the relative lower was in the 5′UTRs, 3′UTRs, and coding regions. Wide SSR analysis of different genomic regions has helped to reveal biological significances of their distributions.

Challenges in analysis and interpretation of microsatellite data for population genetic studies

Advancing technologies have facilitated the ever-widening application of genetic markers such as microsatellites into new systems and research questions in biology. In light of the data and experience accumulated from several years of using microsatellites, we present here a literature review that synthesizes the limitations of microsatellites in population genetic studies. With a focus on population structure, we review the widely used fixation (F ST) statistics and Bayesian clustering algorithms and find that the former can be confusing and problematic for microsatellites and that the latter may be confounded by complex population models and lack power in certain cases. Clustering, multivariate analyses, and diversity-based statistics are increasingly being applied to infer population structure, but in some instances these methods lack formalization with microsatellites. Migration-specific methods perform well only under narrow constraints. We also examine the use of microsatellites for inferring effective population size, changes in population size, and deeper demographic history, and find that these methods are untested and/or highly context-dependent. Overall, each method possesses important weaknesses for use with microsatellites, and there are significant constraints on inferences commonly made using microsatellite markers in the areas of population structure, admixture, and effective population size. To ameliorate and better understand these constraints, researchers are encouraged to analyze simulated datasets both prior to and following data collection and analysis, the latter of which is formalized within the approximate Bayesian computation framework. We also examine trends in the literature and show that microsatellites continue to be widely used, especially in non-human subject areas. This review assists with study design and molecular marker selection, facilitates sound interpretation of microsatellite data while fostering respect for their practical limitations, and identifies lessons that could be applied toward emerging markers and high-throughput technologies in population genetics.

How old are chimpanzee communities? Time to the most recent common ancestor of the Y-chromosome in highly patrilocal societies

Many human societies are patrilineal, with males passing on their name or descent group affiliation to their offspring. Y-chromosomes are also passed on from father to son, leading to the simple expectation that males sharing the same surname or descent group membership should have similar Y-chromosome haplotypes. Although several studies in patrilineal human societies have examined the correspondence between Y-chromosome variation and surname or descent group membership, similar studies in non-human animals are lacking. Chimpanzees represent an excellent species for examining the relationship between descent group membership and Y-chromosome variation because they live in strongly male philopatric communities that arise by a group-fissioning process. Here we take advantage of recent analytical advances in the calculation of the time to the most recent common male ancestor and a large sample size of 273 Y-chromosome short tandem repeat haplotypes to inform our understanding of the potential ages of eight communities of chimpanzees. We find that the times to the most recent common male ancestor of chimpanzee communities are several hundred to as much as over two thousand years. These genetic estimates of the great time depths of chimpanzee communities accord well with behavioral observations suggesting that community fissions are a very rare event and are similar to genetic estimates of the time depth of patrilineal human groups.

Revisiting an important component of plant genomes: microsatellites

Microsatellites are some of the most highly variable repetitive DNA tracts in genomes. Few studies focus on whether the characteristic instability of microsatellites is linked to phenotypic effects in plants. We summarise recent data to investigate how microsatellite variations affect gene expression and hence phenotype. We discuss how the basic characteristics of microsatellites may contribute to phenotypic effects. In summary, microsatellites in plants are universal and highly mutable, they coexist and coevolve with transposable elements, and are under selective pressure. The number of motif nucleotides, the type of motif and transposon activity all contribute to the nonrandom generation and decay of microsatellites, and to conservation and distribution biases. Although microsatellites are generated by accident, they mature through responses to environmental change before final decay. This process is mediated by organism adjustment mechanisms, which maintain a balance between birth versus death and growth versus decay in microsatellites. Close relationships also exist between the physical structure, variation and functionality of microsatellites: in most plant species, sequences containing microsatellites are associated with catalytic activity and binding functions, are expressed in the membrane and organelles, and participate in the developmental and metabolic processes. Microsatellites contribute to genome structure and functional plasticity, and may be considered to promote species evolution in plants in response to environmental changes. In conclusion, the generation, loss, functionality and evolution of microsatellites can be related to plant gene expression and functional alterations. The effect of microsatellites on phenotypic variation may be as significant in plants as it is in animals.

Towards improvements in the estimation of the coalescent: implications for the most effective use of Y chromosome short tandem repeat mutation rates

Over the past two decades, many short tandem repeat (STR) microsatellite loci on the human Y chromosome have been identified together with mutation rate estimates for the individual loci. These have been used to estimate the coalescent age, or the time to the most recent common ancestor (TMRCA) expressed in generations, in conjunction with the average square difference measure (ASD), an unbiased point estimator of TMRCA based upon the average within-locus allele variance between haplotypes. The ASD estimator, in turn, depends on accurate mutation rate estimates to be able to produce good approximations of the coalescent age of a sample. Here, a comparison is made between three published sets of per locus mutation rate estimates as they are applied to the calculation of the coalescent age for real and simulated population samples. A novel evaluation method is developed for estimating the degree of conformity of any Y chromosome STR locus of interest to the strict stepwise mutation model and specific recommendations are made regarding the suitability of thirty-two commonly used Y-STR loci for the purpose of estimating the coalescent. The use of the geometric mean for averaging ASD and across loci is shown to improve the consistency of the resulting estimates, with decreased sensitivity to outliers and to the number of STR loci compared or the particular set of mutation rates selected.

Microsatellites as targets of natural selection

The ability to survey polymorphism on a genomic scale has enabled genome-wide scans for the targets of natural selection. Theory that connects patterns of genetic variation to evidence of natural selection most often assumes a diallelic locus and no recurrent mutation. Although these assumptions are suitable to selection that targets single nucleotide variants, fundamentally different types of mutation generate abundant polymorphism in genomes. Moreover, recent empirical results suggest that mutationally complex, multiallelic loci including microsatellites and copy number variants are sometimes targeted by natural selection. Given their abundance, the lack of inference methods tailored to the mutational peculiarities of these types of loci represents a notable gap in our ability to interrogate genomes for signatures of natural selection. Previous theoretical investigations of mutation-selection balance at multiallelic loci include assumptions that limit their application to inference from empirical data. Focusing on microsatellites, we assess the dynamics and population-level consequences of selection targeting mutationally complex variants. We develop general models of a multiallelic fitness surface, a realistic model of microsatellite mutation, and an efficient simulation algorithm. Using these tools, we explore mutation-selection-drift equilibrium at microsatellites and investigate the mutational history and selective regime of the microsatellite that causes Friedreich's ataxia. We characterize microsatellite selective events by their duration and cost, note similarities to sweeps from standing point variation, and conclude that it is premature to label microsatellites as ubiquitous agents of efficient adaptive change. Together, our models and simulation algorithm provide a powerful framework for statistical inference, which can be used to test the neutrality of microsatellites and other multiallelic variants.

Measuring microsatellite conservation in mammalian evolution with a phylogenetic birth-death model

Microsatellites make up ∼3% of the human genome, and there is increasing evidence that some microsatellites can have important functions and can be conserved by selection. To investigate this conservation, we performed a genome-wide analysis of human microsatellites and measured their conservation using a binary character birth--death model on a mammalian phylogeny. Using a maximum likelihood method to estimate birth and death rates for different types of microsatellites, we show that the rates at which microsatellites are gained and lost in mammals depend on their sequence composition, length, and position in the genome. Additionally, we use a mixture model to account for unequal death rates among microsatellites across the human genome. We use this model to assign a probability-based conservation score to each microsatellite. We found that microsatellites near the transcription start sites of genes are often highly conserved, and that distance from a microsatellite to the nearest transcription start site is a good predictor of the microsatellite conservation score. An analysis of gene ontology terms for genes that contain microsatellites near their transcription start site reveals that regulatory genes involved in growth and development are highly enriched with conserved microsatellites.

Microevolution of pandemic Vibrio parahaemolyticus assessed by the number of repeat units in short sequence tandem repeat regions

The emergence of the pandemic strain Vibrio parahaemolyticus O3:K6 in 1996 caused a large increase of diarrhea outbreaks related to seafood consumption in Southeast Asia, and later worldwide. Isolates of this strain constitutes a clonal complex, and their effectual differentiation is possible by comparison of their variable number tandem repeats (VNTRs). The differentiation of the isolates by the differences in VNTRs will allow inferring the population dynamics and microevolution of this strain but this requires knowing the rate and mechanism of VNTRs' variation. Our study of mutants obtained after serial cultivation of clones showed that mutation rates of the six VNTRs examined are on the order of 10(-4) mutant per generation and that difference increases by stepwise addition of single mutations. The single stepwise mutation (SSM) was deduced because mutants with 1, 2, 3, or more repeat unit deletions or insertions follow a geometric distribution. Plausible phylogenetic trees are obtained when, according to SSM, the genetic distance between clusters with different number of repeats is assessed by the absolute differences in repeats. Using this approach, mutants originated from different isolates of pandemic V. parahaemolyticus after serial cultivation are clustered with their parental isolates. Additionally, isolates of pandemic V. parahaemolyticus from Southeast Asia, Tokyo, and northern and southern Chile are clustered according their geographical origin. The deepest split in these four populations is observed between the Tokyo and southern Chile populations. We conclude that proper phylogenetic relations and successful tracing of pandemic V. parahaemolyticus requires measuring the differences between isolates by the absolute number of repeats in the VNTRs considered.

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.

Genetic diversity of microsatellite loci in hierarchically structured populations

Microsatellite loci are widely used for investigating patterns of genetic variation within and among populations. Those patterns are in turn determined by population sizes, migration rates, and mutation rates. We provide exact expressions for the first two moments of the allele frequency distribution in a stochastic model appropriate for studying microsatellite evolution with migration, mutation, and drift under the assumption that the range of allele sizes is bounded. Using these results, we study the behavior of several measures related to Wright's F(ST), including Slatkin's R(ST). Our analytical approximations for F(ST) and R(ST) show that familiar relationships between N(e)m and F(ST) or R(ST) hold when the migration and mutation rates are small. Using the exact expressions for F(ST) and R(ST), our numerical results show that, when the migration and mutation rates are large, these relationships no longer hold. Our numerical results also show that the diversity measures most closely related to F(ST) depend on mutation rates, mutational models (stepwise versus two-phase), migration rates, and population sizes. Surprisingly, R(ST) is relatively insensitive to the mutation rates and mutational models. The differing behaviors of R(ST) and F(ST) suggest that properties of the among-population distribution of allele frequencies may allow the roles of mutation and migration in producing patterns of diversity to be distinguished, a topic of continuing investigation.

Evidence for an invasive aphid "superclone": extremely low genetic diversity in Oleander aphid (Aphis nerii) populations in the southern United States

Background

The importance of genetic diversity in successful biological invasions is unclear. In animals, but not necessarily plants, increased genetic diversity is generally associated with successful colonization and establishment of novel habitats. The Oleander aphid, Aphis nerii, though native to the Mediterranean region, is an invasive pest species throughout much of the world. Feeding primarily on Oleander (Nerium oleander) and Milkweed (Asclepias spp.) under natural conditions, these plants are unlikely to support aphid populations year round in the southern US. The objective of this study was to describe the genetic variation within and among US populations of A. nerii, during extinction/recolonization events, to better understand the population ecology of this invasive species.

Methodology/Principal Findings

We used five microsatellite markers to assess genetic diversity over a two year period within and among three aphid populations separated by small (100 km) and large (3,700 km) geographic distances on two host plant species. Here we provide evidence for A. nerii “superclones”. Genotypic variation was absent in all populations (i.e., each population consisted of a single multilocus genotype (MLG) or “clone”) and the genetic composition of only one population completely changed across years. There was no evidence of sexual reproduction or host races on different plant species.

Conclusions/Significance

Aphis nerii is a well established invasive species despite having extremely low genetic diversity. As this aphid appears to be obligatorily asexual, it may share more similarities with clonally reproducing invasive plants, than with other animals. Patterns of temporal and geographic genetic variation, viewed in the context of its population dynamics, have important implications for the management of invasive pests and the evolutionary biology of asexual species.

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.

Phylogeography, genetic structure and population divergence time of cheetahs in Africa and Asia: evidence for long-term geographic isolates

The cheetah (Acinonyx jubatus) has been described as a species with low levels of genetic variation. This has been suggested to be the consequence of a demographic bottleneck 10 000–12 000 years ago (ya) and also led to the assumption that only small genetic differences exist between the described subspecies. However, analysing mitochondrial DNA and microsatellites in cheetah samples from most of the historic range of the species we found relatively deep phylogeographic breaks between some of the investigated populations, and most of the methods assessed divergence time estimates predating the postulated bottleneck. Mitochondrial DNA monophyly and overall levels of genetic differentiation support the distinctiveness of Northern-East African cheetahs (Acinonyx jubatus soemmeringii). Moreover, combining archaeozoological and contemporary samples, we show that Asiatic cheetahs (Acinonyx jubatus venaticus) are unambiguously separated from African subspecies. Divergence time estimates from mitochondrial and nuclear data place the split between Asiatic and Southern African cheetahs (Acinonyx jubatus jubatus) at 32 000–67 000 ya using an average mammalian microsatellite mutation rate and at 4700–44 000 ya employing human microsatellite mutation rates. Cheetahs are vulnerable to extinction globally and critically endangered in their Asiatic range, where the last 70–110 individuals survive only in Iran. We demonstrate that these extant Iranian cheetahs are an autochthonous monophyletic population and the last representatives of the Asiatic subspecies A. j. venaticus. We advocate that conservation strategies should consider the uncovered independent evolutionary histories of Asiatic and African cheetahs, as well as among some African subspecies. This would facilitate the dual conservation priorities of maintaining locally adapted ecotypes and genetic diversity.

Mutation biases and mutation rate variation around very short human microsatellites revealed by human-chimpanzee-orangutan genomic sequence alignments

I have studied mutation patterns around very short microsatellites, focusing mainly on sequences carrying only two repeat units. By using human-chimpanzee-orangutan alignments, inferences can be made about both the relative rates of mutations and which bases have mutated. I find remarkable non-randomness, with mutation rate depending on a base's position relative to the microsatellite, the identity of the base itself and the motif in the microsatellite. Comparing the patterns around AC2 with those around other four-base combinations reveals that AC2 does not stand out as being special in the sense that non-repetitive tetramers also generate strong mutation biases. However, comparing AC2 and AC3 with AC4 reveals a step change in both the rate and nature of mutations occurring, suggesting a transition state, AC4 exhibiting an alternating high-low mutation rate pattern consistent with the sequence patterning seen around longer microsatellites. Surprisingly, most changes in repeat number occur through base substitutions rather than slippage, and the relative probability of gaining versus losing a repeat in this way varies greatly with repeat number. Slippage mutations reveal rather similar patterns of mutability compared with point mutations, being rare at two repeats where most cause the loss of a repeat, with both mutation rate and the proportion of expansion mutations increasing up to 6-8 repeats. Inferences about longer repeat tracts are hampered by uncertainties about the proportion of multi-species alignments that fail due to multi-repeat mutations and other rearrangements.

Comparative genomics and molecular dynamics of DNA repeats in eukaryotes

Repeated elements can be widely abundant in eukaryotic genomes, composing more than 50% of the human genome, for example. It is possible to classify repeated sequences into two large families, "tandem repeats" and "dispersed repeats." Each of these two families can be itself divided into subfamilies. Dispersed repeats contain transposons, tRNA genes, and gene paralogues, whereas tandem repeats contain gene tandems, ribosomal DNA repeat arrays, and satellite DNA, itself subdivided into satellites, minisatellites, and microsatellites. Remarkably, the molecular mechanisms that create and propagate dispersed and tandem repeats are specific to each class and usually do not overlap. In the present review, we have chosen in the first section to describe the nature and distribution of dispersed and tandem repeats in eukaryotic genomes in the light of complete (or nearly complete) available genome sequences. In the second part, we focus on the molecular mechanisms responsible for the fast evolution of two specific classes of tandem repeats: minisatellites and microsatellites. Given that a growing number of human neurological disorders involve the expansion of a particular class of microsatellites, called trinucleotide repeats, a large part of the recent experimental work on microsatellites has focused on these particular repeats, and thus we also review the current knowledge in this area. Finally, we propose a unified definition for mini- and microsatellites that takes into account their biological properties and try to point out new directions that should be explored in a near future on our road to understanding the genetics of repeated sequences.

The SMM model as a boundary value problem using the discrete diffusion equation

A generalized single-step stepwise mutation model (SMM) is developed that takes into account an arbitrary initial state to a certain partial difference equation. This is solved in both the approximate continuum limit and the more exact discrete form. A time evolution model is developed for Y DNA or mtDNA that takes into account the reflective boundary modeling minimum microsatellite length and the original difference equation. A comparison is made between the more widely known continuum Gaussian model and a discrete model, which is based on modified Bessel functions of the first kind. A correction is made to the SMM model for the probability that two individuals are related that takes into account a reflecting boundary modeling minimum microsatellite length. This method is generalized to take into account the general n-step model and exact solutions are found. A new model is proposed for the step distribution.

Mining microsatellites in eukaryotic genomes

During recent decades, microsatellites have become the most popular source of genetic markers. More recently, the availability of enormous sequence data for a large number of eukaryotic genomes has accelerated research aimed at understanding the origin and functions of microsatellites and searching for new applications. This review presents recent developments of in silico mining of microsatellites to reveal various facets of the distribution and dynamics of microsatellites in eukaryotic genomes. Two aspects of microsatellite search strategies--using a suitable search tool and accessing a relevant microsatellite database--have been explored. Judicious microsatellite mining not only helps in addressing biological questions but also facilitates better exploitation of microsatellites for diverse applications.

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.

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.

Microsatellite evolution: Markov transition functions for a suite of models

This paper takes from the collection of models considered by Whittaker et al. [2003. Likelihood-based estimation of microsatellite mutation rates. Genetics 164, 781-787] derived from direct observation of microsatellite mutation in parent-child pairs and provides analytical expressions for the probability distributions for the change in number of repeats over any given number of generations. The mathematical framework for this analysis is the theory of Markov processes. We find these expressions using two approaches, approximating by circulant matrices and solving a partial differential equation satisfied by the generating function. The impact of the differing choice of models is examined using likelihood estimates for time to most recent common ancestor. The analysis presented here may play a role in elucidating the connections between these two approaches and shows promise in reconciling differences between estimates for mutation rates based on Whittaker's approach and methods based on phylogenetic analyses.

Some probabilistic results on the nonrandomness of simple sequence repeats in DNA sequences

Some probabilistic results on simple sequence repeats (SSRs) in DNA sequences are derived and used to quantify the nonrandomness of SSRs as an index of nonrandomness. The applicability of the index of nonrandomness is illustrated using several examples from the literature on selected human diseased genes.

Long repeats in a huge genome: microsatellite loci in the grasshopper Chorthippus biguttulus

It is commonly believed that both the average length and the frequency of microsatellites correlate with genome size. We have estimated the frequency and the average length for 69 perfect dinucleotide microsatellites in an insect with an exceptionally large genome: Chorthippus biguttulus (Orthoptera, Acrididae). Dinucleotide microsatellites are not more frequent in C. biguttulus, but repeat arrays are 1.4 to 2 times longer than in other insect species. The average repeat number in C. biguttulus lies in the range of higher vertebrates. Natural populations are highly variable. At least 30 alleles per locus were found and the expected heterozygosity is above 0.95 at all three loci studied. In contrast, the observed heterozygosity is much lower (< or = 0.51), which could be caused by long null alleles.

Inference on microsatellite mutation processes in the invasive mite, Varroa destructor, using reversible jump Markov chain Monte Carlo

Varroa destructor is a parasitic mite of the Eastern honeybee Apis cerana. Fifty years ago, two distinct evolutionary lineages (Korean and Japanese) invaded the Western honeybee Apis mellifera. This haplo-diploid parasite species reproduces mainly through brother-sister matings, a system which largely favors the fixation of new mutations. In a worldwide sample of 225 individuals from 21 locations collected on Western honeybees and analyzed at 19 microsatellite loci, a series of de novo mutations was observed. Using historical data concerning the invasion, this original biological system has been exploited to compare three mutation models with allele size constraints for microsatellite markers: stepwise (SMM) and generalized (GSM) mutation models, and a model with mutation rate increasing exponentially with microsatellite length (ESM). Posterior probabilities of the three models have been estimated for each locus individually using reversible jump Markov Chain Monte Carlo. The relative support of each model varies widely among loci, but the GSM is the only model that always receives at least 9% support, whatever the locus. The analysis also provides robust estimates of mutation parameters for each locus and of the divergence time of the two invasive lineages (67,000 generations with a 90% credibility interval of 35,000-174,000). With an average of 10 generations per year, this divergence time fits with the last post-glacial Korea-Japan land separation.

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.

Multiple origins of Ashkenazi Levites: Y chromosome evidence for both Near Eastern and European ancestries

Previous Y chromosome studies have shown that the Cohanim, a paternally inherited Jewish priestly caste, predominantly share a recent common ancestry irrespective of the geographically defined post-Diaspora community to which they belong, a finding consistent with common Jewish origins in the Near East. In contrast, the Levites, another paternally inherited Jewish caste, display evidence for multiple recent origins, with Ashkenazi Levites having a high frequency of a distinctive, non-Near Eastern haplogroup. Here, we show that the Ashkenazi Levite microsatellite haplotypes within this haplogroup are extremely tightly clustered, with an inferred common ancestor within the past 2,000 years. Comparisons with other Jewish and non-Jewish groups suggest that a founding event, probably involving one or very few European men occurring at a time close to the initial formation and settlement of the Ashkenazi community, is the most likely explanation for the presence of this distinctive haplogroup found today in >50% of Ashkenazi Levites.

Two distinct modes of microsatellite mutation processes: evidence from the complete genomic sequences of nine species

We surveyed microsatellite distribution in 10 completely sequenced genomes. Using a permutation-based statistic, we assessed for all 10 genomes whether the microsatellite distribution significantly differed from expectations. Consistent with previous reports, we observed a highly significant excess of long microsatellites. Focusing on short microsatellites containing only a few repeat units, we demonstrate that this repeat class is significantly underrepresented in most genomes. This pattern was observed across different repeat types. Computer simulations indicated that neither base substitutions nor a combination of length-dependent slippage and base substitutions could explain the observed pattern of microsatellite distribution. When we introduced one additional mutation process, a length-independent slippage (indel slippage) operating at repeats with few repetitions, our computer simulations captured the observed pattern of microsatellite distribution.

Active Alu element "A-tails": size does matter

Long and short interspersed elements (LINEs and SINEs) are retroelements that make up almost half of the human genome. L1 and Alu represent the most prolific human LINE and SINE families, respectively. Only a few Alu elements are able to retropose, and the factors determining their retroposition capacity are poorly understood. The data presented in this paper indicate that the length of Alu "A-tails" is one of the principal factors in determining the retropositional capability of an Alu element. The A stretches of the Alu subfamilies analyzed, both old (Alu S and J) and young (Ya5), had a Poisson distribution of A-tail lengths with a mean size of 21 and 26, respectively. In contrast, the A-tails of very recent Alu insertions (disease causing) were all between 40 and 97 bp in length. The L1 elements analyzed displayed a similar tendency, in which the "disease"-associated elements have much longer A-tails (mean of 77) than do the elements even from the young Ta subfamily (mean of 41). Analysis of the draft sequence of the human genome showed that only about 1000 of the over one million Alu elements have tails of 40 or more adenosine residues in length. The presence of these long A stretches shows a strong bias toward the actively amplifying subfamilies, consistent with their playing a major role in the amplification process. Evaluation of the 19 Alu elements retrieved from the draft sequence of the human genome that are identical to the Alu Ya5a2 insert in the NF1 gene showed that only five have tails with 40 or more adenosine residues. Sequence analysis of the loci with the Alu elements containing the longest A-tails (7 of the 19) from the genomes of the NF1 patient and the father revealed that there are at least two loci with A-tails long enough to serve as source elements within our model. Analysis of the A-tail lengths of 12 Ya5a2 elements in diverse human population groups showed substantial variability in both the Alu A-tail length and sequence homogeneity. On the basis of these observations, a model is presented for the role of A-tail length in determining which Alu elements are active.

Mismatch repair-driven mutational bias in D. melanogaster

We have used microsatellite sequences to evaluate the influence of the mismatch repair system on mutation bias in D. melanogaster. While mismatch-proficient cells have the highest mutation rate at (GT)(n) repeats, (AT)(n) repeats were the least stable ones in spel1(-/-) flies lacking functional mismatch repair. Furthermore, the mutation spectrum of long microsatellite alleles in spel1(-/-) was slightly upward biased, resulting in a gain of repeats, whereas wild-type flies have a strong downward bias. Interestingly, this mismatch repair-mediated downward mutation bias is reflected in the genome composition of D. melanogaster. When compared to other species, D. melanogaster has significantly shorter microsatellites. Our results suggest that the mismatch repair system may have an important role in shaping genome composition.

Microsatellite evolution inferred from human-chimpanzee genomic sequence alignments

Most studies of microsatellite evolution utilize long, highly mutable loci, which are unrepresentative of the majority of simple repeats in the human genome. Here we use an unbiased sample of 2,467 microsatellite loci derived from alignments of 5.1 Mb of genomic sequence from human and chimpanzee to investigate the mutation process of tandemly repetitive DNA. The results indicate that the process of microsatellite evolution is highly heterogeneous, exhibiting differences between loci of different lengths and motif sizes and between species. We find a highly significant tendency for human dinucleotide repeats to be longer than their orthologues in chimpanzees, whereas the opposite trend is observed in mononucleotide repeat arrays. Furthermore, the rate of divergence between orthologues is significantly higher at longer loci, which also show significantly greater mutability per repeat number. These observations have important consequences for understanding the molecular mechanisms of microsatellite mutation and for the development of improved measures of genetic distance.

Heterogeneity in the rate and pattern of germline mutation at individual microsatellite loci

There is a lack of information on how individual microsatellite loci differ with respect to their mutation properties. Such variation will have an important bearing on our understanding of the ubiquitous occurrence of simple repeat sequences in eukaryotic genomes and on deriving proper mutation models that can be incorporated into genetic distance estimates. We genotyped approximately 100 families of the bird barn swallow (Hirundo rustica) for two hypervariable (heterozygosity >95%) microsatellite markers: HrU6, an (AAAG)n tetranucleotide repeat, and HrU10, an (AAGAG)n pentanucleotide repeat. A total of 27 germline mutation events were documented, corresponding to mutation rates of 0.57% (HrU6) and 1.56% (HrU10). The mutation rate increased with allele size, at approximately 0.1% per repeat unit over the observed range of allele sizes (approximately 10-100 repeat units). Single repeat unit changes dominated, with 21/27 mutations representing the gain or loss of one repeat unit. There was no clear difference in the number of gains versus losses nor was there an effect of allele size on the magnitude or direction of mutation. Unexpectedly, the mutation rate of females (maternally transmitted mutations) was 2.5-5 times higher than that of males. Contrasting these observations with mutation data from other microsatellite loci reveals differences not only in the mutation rate, but also in the magnitude, direction and effect of sex on mutation. Thus, microsatellite mutation and evolution may be viewed as a dynamic and variable process.