Evolutionary dynamics of microsatellite DNA

Within the past decade microsatellites have developed into one of the most popular genetic markers. Despite the widespread use of microsatellite analysis, an integral picture of the mutational dynamics of microsatellite DNA is just beginning to emerge. Here, I review both generally agreed and controversial results about the mutational dynamics of microsatellite DNA. Microsatellites are short DNA sequence stretches in which a motif of one to six bases is tandemly repeated. It has been known for some time that these sequences can differ in repeat number among individuals. With the advent of polymerase chain reaction (PCR) technology this property of microsatellite DNA was converted into a highly versatile genetic marker (Litt and Luty 1989; Tautz 1989; Weber and May 1989). Polymerase chain reaction products of different length can be amplified with primers flanking the variable microsatellite region. Due to the availability of high-throughput capillary sequencers or mass spectrography the sizing of alleles is no longer a bottleneck in microsatellite analysis. The almost random distribution of microsatellites and their high level of polymorphism greatly facilitated the construction of genetic maps (Dietrich et al. 1994; Dib et al. 1996) and enabled subsequent positional cloning of several genes. Almost at the same time, microsatellites were established as the marker of choice for the identification of individuals and paternity testing. The high sensitivity of PCR-based microsatellite analysis was not only of great benefit for forensics, but opened completely new research areas, such as the analysis of samples with limited DNA amounts (e.g., many social insects) or degraded DNA (e.g., feces, museum material) (Schl?tterer and Pemberton 1998). More recently, microsatellite analysis has also been employed in population genetics (Goldstein and Schl?tterer 1999). Compared with allozymes, microsatellites offer the advantage that, in principle, several thousand potentially polymorphic markers are available. Nevertheless, the application of microsatellites to population genetic questions requires a more detailed understanding of the mutation processes of microsatellite DNA as the evolutionary time frames covered in population genetics are often too long to allow novel microsatellite mutations to be ignored. Additional interest in the evolution of microsatellite DNA comes from the discovery that trinucleotide repeats, a special class of microsatellites, are involved in human neurodegenerative diseases (e.g., fragile X and Huntington's disease). A detailed understanding of the processes underlying microsatellite instability is therefore an important contribution toward a better understanding of these human neurodegenerative diseases.     

Mutational Dynamics of Microsatellites

Microsatellites are a ubiquitous class of simple repetitive DNA sequences, which are widespread in both eukaryotic and prokaryotic genomes. The use of microsatellites as polymorphic DNA markers has considerably increased both in the number of studies and in the number of organisms, primarily for genetic mapping, studying genomic instability in cancer, population genetics, forensics, conservation biology, molecular anthropology and in the studies of human evolutionary history. Although simple sequence repeats have been extensively used in studies encompassing varied areas of genetics, the mutation dynamics of these genome regions is still not well understood. The present review focuses on the mutational dynamics of microsatellite DNA with special reference to mutational mechanisms and their role in microsatellite evolution.     

Microsatellite repeats are not randomly distributed within Norway spruce (Picea abies K.) expressed sequences.

A Norway spruce (Picea abies K.) cDNA library obtained from vegetative bud tissue was screened for the presence of (AG)n and (AC)n microsatellite repeats. Ten (AG)n and six (AC)n microsatellites were found, with an average length of 25.5 repeat units. Most of the microsatellites are simple perfect repeats. The microsatellite distribution within the clones is clearly non-random, with different classes of repeats lying in different positions relative to the coding region and in a highly conserved orientation. An estimate of the frequency of dinucleotide microsatellites in expressed regions was obtained, showing that SSRs (simple sequence repeats) are found in genes about 20 times less frequently than in random genomic clones, with (AG)n repeats more frequent than (AC)n repeats. Potential applications of these sequences as expressed region-based molecular markers are shown by developing six SSR markers for the detection of natural variation in Norway spruce populations and testing two of them for the identification of illegitimate progenies from a mapping population.     

Evidence for widespread convergent evolution around human microsatellites

Microsatellites are a major component of the human genome, and their evolution has been much studied. However, the evolution of microsatellite flanking sequences has received less attention, with reports of both high and low mutation rates and of a tendency for microsatellites to cluster. From the human genome we generated a database of many thousands of (AC)_n flanking sequences within which we searched for common characteristics. Sequences flanking microsatellites of similar length show remarkable levels of convergent evolution, indicating shared mutational biases. These biases extend 25–50 bases either side of the microsatellite and may therefore affect more than 30% of the entire genome. To explore the extent and absolute strength of these effects, we quantified the observed convergence. We also compared homologous human and chimpanzee loci to look for evidence of changes in mutation rate around microsatellites. Most models of DNA sequence evolution assume that mutations are independent and occur randomly. Allowances may be made for sites mutating at different rates and for general mutation biases such as the faster rate of transitions over transversions. Our analysis suggests that these models may be inadequate, in that proximity to even very short microsatellites may alter the rate and distribution of mutations that occur. The elevated local mutation rate combined with sequence convergence, both of which we find evidence for, also provide a possible resolution for the apparently contradictory inferences of mutation rates in microsatellite flanking sequences.     

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 diversity of loquat germplasm (Eriobotrya japonica (Thunb) Lindl) assessed by SSR markers.

Genetic relationships among 40 loquat (Eriobotrya japonica (Thunb) Lindl) accessions that originated from different countries and that are part of the germplasm collection of the Instituto Valenciano de Investigaciones Agrarias (IVIA) (Valencia, Spain) were evaluated using microsatellites. Thirty primer pairs flanking microsatellites previously identified in Malus x domestica (Borkh.) were assayed. Thirteen of them amplified polymorphic products and unambiguously distinguished 34 genotypes from the 40 accessions analyzed. Six accessions showing identical marker patterns were Spanish local varieties thought to have been derived from 'Algerie' by a mutational process very common in loquat species. A total of 39 alleles were detected in the population studied, with a mean value of 2.4 alleles per locus. The expected and observed heterozygosities were 0.46 and 51% on average, respectively, leading to a negative value of the Wright's fixation index (-0.20). The values of these parameters indicate a smaller degree of genetic diversity in the set of loquat accessions analyzed than in other members of the Rosaceae family. Unweighted pair-group method (UPGMA) cluster analysis, based on Nei's genetic distance, generally grouped genotypes according to their geographic origins and pedigrees. The high number of alleles and the high expected heterozygosity detected with SSR markers developed in Malus x domestica (Borkh.) make them a suitable tool for loquat cultivar identification, confirming microsatellite marker transportability among genera in the Rosaceae family.     

Linkage disequilibrium on the bovine X chromosome: characterization and use in quantitative trait locus mapping

We herein demonstrate that in the Holstein-Friesian dairy cattle population, microsatellites are as polymorphic on the X chromosome as on the autosomes but that the level of linkage disequilibrium between these markers is higher on the X chromosome than on the autosomes. The latter observation is not compatible with the small male-to-female ratio that prevails in this population and results in a higher gonosomal than autosomal effective population size. It suggests that the X chromosome undergoes distinct selective or mutational forces. We describe and characterize a novel Markovian approach to exploit this linkage disequilibrium to compute the probability that two chromosomes are identical-by-descent conditional on flanking marker data. We use the ensuing probabilities in a restricted maximum-likelihood approach to search for quantitative trait loci (QTL) affecting 48 traits of importance to the dairy industry and provide evidence for the presence of QTL affecting 5 of these traits on the bovine X chromosome.     

Heterosis, marker mutational processes and population inbreeding history.

Genotype-fitness correlations (GFC) have previously been studied using allozyme markers and have often focused on short-term processes such as recent inbreeding. Thus, models of GFC usually neglect marker mutation and only use heterozygosity as a genotypic index. Recently, GFC have also been reported (i) with DNA markers such as microsatellites, characterized by high mutation rates and specific mutational processes and (ii) using new individual genotypic indices assumed to be more precise than heterozygosity. The aim of this article is to evaluate the theoretical impact of marker mutation on GFC. We model GFC due to short-term processes generated by the current breeding system (partial selfing) and to long-term processes generated by past population history (hybridization). Various mutation rates and mutation models corresponding to different kinds of molecular markers are considered. Heterozygosity is compared to other genotypic indices designed for specific marker types. Highly mutable markers (such as microsatellites) are particularly suitable for the detection of GFC that evolve in relation to short-term processes, whereas GFC due to long-term processes are best observed with intermediate mutation rates. Irrespective of the marker type and population scenario, heterozygosity usually provides higher correlations than other genotypic indices under most biologically plausible conditions.     

IRAP and REMAP for retrotransposon-based genotyping and fingerprinting

Retrotransposons can be used as markers because their integration creates new joints between genomic DNA and their conserved ends. To detect polymorphisms for retrotransposon insertion, marker systems generally rely on PCR amplification between these ends and some component of flanking genomic DNA. We have developed two methods, retrotransposon-microsatellite amplified polymorphism (REMAP) analysis and inter-retrotransposon amplified polymorphism (IRAP) analysis, that require neither restriction enzyme digestion nor ligation to generate the marker bands. The IRAP products are generated from two nearby retrotransposons using outward-facing primers. In REMAP, amplification between retrotransposons proximal to simple sequence repeats (microsatellites) produces the marker bands. Here, we describe protocols for the IRAP and REMAP techniques, including methods for PCR amplification with a single primer or with two primers and for agarose gel electrophoresis of the product using optimal electrophoresis buffers and conditions. This protocol can be completed in 1-2 d.     

Conserved flanking microsatellite sequences (ReFS) differentiate between Lepidoptera species,and provide insight into microsatellite evolution

High rates of mutation and homoplasy mean that microsatellites generally are not considered to be useful molecular markers for inferring systematic relationships between species. However, an earlier pilot study suggested that conserved flanking microsatellite sequences, also known as repetitive flanking sequences (ReFS), may form a basis for a dominant marker that can differentiate between species of Lepidoptera. We present data that demonstrate that ReFS are quick and easy to use, and generate highly repeatable banding patterns from a range of Lepidoptera species. Sequence data from a subset of ReFS‐amplified bands revealed microsatellite families with flanking sequences that are more conserved within than among species: this is probably attributable to recombination‐mediated events, transposition of mobile elements or a combination of the two. Our data support the use of ReFS as dominant interspecific molecular markers, and add to the growing literature on the evolution of microsatellites in Lepidoptera.     

Heterogeneous Nature and Distribution of Interruptions in Dinucleotides May Indicate the Existence of Biased Substitutions Underlying Microsatellite Evolution

Some aspects of microsatellite evolution, such as the role of base substitutions, are far from being fully understood. To examine the significance of base substitutions underlying the evolution of microsatellites we explored the nature and the distribution of interruptions in dinucleotide repeats from the human genome. The frequencies that we inferred in the repetitive sequences were statistically different from the frequencies observed in other noncoding sequences. Additionally, we detected that the interruptions tended to be towards the ends of the microsatellites and 5'-3' asymmetry. In all the estimates nucleotides forming the same repetitive motif seem to be affected by different base substitution rates in AC and AG. This tendency itself could generate patterning and similarity in flanking sequences and reconcile these phenomena with the high mutation rate found in flanking sequences without invoking convergent evolution. Nevertheless, our data suggest that there is a regional bias in the substitution pattern of microsatellites. The accumulation of random substitutions alone cannot explain the heterogeneity and the asymmetry of interruptions found in this study or the relative frequency of different compound microsatellites in the human genome. Therefore, we cannot rule out the possibility of a mutational bias leading to convergent or parallel evolution in flanking sequences.     

Linkage disequilibrium and founder effect analysis of the NF1 gene in French Canadians from the Quebec population

We genotyped 19 neurofibromatosis type 1 (NF1) families from French Canadians of the Quebec population with four intragenic microsatellites (IVS26-2.3, IVS27AC28.4, IVS27AC33.1, and IVS38GT53.0). Linkage analysis of the four microsatellite markers among the 19 NF1 families indicates that the four microsatellites are strongly linked with NF1 disease (LOD = 2.76-3.64). The four markers are associated (P = 0-0.077) except marker pair IVS26-2.3/IVS27AC33.1 (P = 0.18 or 0.17). However, perhaps due to the high mutation rate of the NF1 gene, no founder effect for NF1 was detected in the Quebec French Canadians.     

Comparative genomic analysis reveals species-dependent complexities that explain difficulties with microsatellite marker development in molluscs

Reliable population DNA molecular markers are difficult to develop for molluscs, the reasons for which are largely unknown. Identical protocols for microsatellite marker development were implemented in three gastropods. Success rates were lower for Gibbula cineraria compared to Littorina littorea and L. saxatilis. Comparative genomic analysis of 47.2 kb of microsatellite containing sequences (MCS) revealed a high incidence of cryptic repetitive DNA in their flanking regions. The majority of these were novel, and could be grouped into DNA families based upon sequence similarities. Significant inter-specific variation in abundance of cryptic repetitive DNA and DNA families was observed. Repbase scans show that a large proportion of cryptic repetitive DNA was identified as transposable elements (TEs). We argue that a large number of TEs and their transpositional activity may be linked to differential rates of DNA multiplication and recombination. This is likely to be an important factor explaining inter-specific variation in genome stability and hence microsatellite marker development success rates. Gastropods also differed significantly in the type of TEs classes (autonomous vs non-autonomous) observed. We propose that dissimilar transpositional mechanisms differentiate the TE classes in terms of their propensity for transposition, fixation and/or silencing. Consequently, the phylogenetic conservation of non-autonomous TEs, such as CvA, suggests that dispersal of these elements may have behaved as microsatellite-inducing elements. Results seem to indicate that, compared to autonomous, non-autonomous TEs maybe have a more active role in genome rearrangement processes. The implications of the findings for genomic rearrangement, stability and marker development are discussed.     

GinMicrosatDb: a genome-wide microsatellite markers database for sesame (<Emphasis Type="Italic">Sesamum indicum</Emphasis> L.)

Molecular breeding in sesame is still at infancy due to limited number of microsatellite markers available and the low level of polymorphism exhibited by them. Therefore, whole genome sequencing was used for development of microsatellite markers so as to ensure availability of substantial number of polymorphic markers for use in marker assisted breeding programs. Whole genome sequencing of sesame variety ‘Swetha’ was done using Illumina paired-end sequencing and Roche 454 shotgun sequencing technologies (GCA_000975565.1 in GenBank). ‘GinMicrosatDb’, a genome-wide microsatellite marker database has been developed using the whole genome sequence data of sesame variety ‘Swetha’. The database consists of microsatellites localized on both linkage groups and scaffolds with their genomic co-ordinates. It provides five sets of forward and reverse primers for each of the microsatellite loci along with the flanking sequences, primer GC content, product size and melting temperature etc. The distribution of microsatellites can be viewed and selected through a genome browser as well as through a physical map. The newly identified microsatellite markers are expected to help sesame breeders in developing marker tags for traits of economic importance thereby bringing about greater efficiency in marker-assisted selection programs.     

Eleven polymorphic microsatellites isolated from red swamp crayfish, Procambarus clarkii

Eleven microsatellites were isolated from a genomic DNA library enriched for CA- and GA-repeats, and were characterized in 48 individuals of red swamp crayfish Procambarus clarkii. All 11 microsatellites were polymorphic with an average of allele number of 6.6 per locus. Genotypic distributions of these 11 markers were found to be in conformance with the expectations of Hardy-Weinberg equilibrium. All 11 markers were unlinked. These markers are being used to study the invasion routine, genetic diversity and population structure of the species P. clarkii.     

Quantifying ascertainment bias and species-specific length differences in human and chimpanzee microsatellites using genome sequences

Surveys of variability of homologous microsatellite loci among species reveal an ascertainment bias for microsatellite length where microsatellite loci isolated in one species tend to be longer than homologous loci in related species. Here, we take advantage of the availability of aligned human and chimpanzee genome sequences to compare length difference of homologous microsatellites for loci identified in humans to length difference for loci identified in chimpanzees. We are able to quantify ascertainment bias for a range of motifs and microsatellite lengths. Because ascertainment bias should not exist if a microsatellite selected in one species is as likely to be longer as it is to be shorter than its homologue, we propose that the nature of ascertainment bias can provide evidence for understanding how microsatellites evolve. We show that bias is greater for longer microsatellites but also that many long microsatellites have short homologues. These results are consistent with the notion that growth of long microsatellites is constrained by an upper length boundary that, when reached, sometimes results in large deletions. By evaluating ascertainment bias separately for interrupted and uninterrupted repeats we also show that long microsatellites tend to become interrupted, thereby contributing a second component of ascertainment bias. Having accounted for ascertainment bias, in agreement with results published elsewhere, we find that microsatellites in humans are longer on average than those in chimpanzees. This length difference is similar among repeat motifs but surprisingly comprises two roughly equal components, one associated with the repeats themselves and one with the flanking sequences. The differences we find can only be explained if microsatellites are both evolving directionally under a biased mutation process and are doing so at different rates in different closely related species.     

Development of a genome-wide anchored microsatellite map for common bean (<Emphasis Type="Italic">Phaseolus vulgaris</Emphasis> L.)

A total of 150 microsatellite markers developed for common bean (Phaseolus vulgaris L.) were tested for parental polymorphism and used to determine the positions of 100 genetic loci on an integrated genetic map of the species. The value of these single-copy markers was evident in their ability to link two existing RFLP-based genetic maps with a base map developed for the Mesoamerican × Andean population, DOR364 × G19833. Two types of microsatellites were mapped, based respectively on gene-coding and anonymous genomic-sequences. Gene-based microsatellites proved to be less polymorphic (46.3%) than anonymous genomic microsatellites (64.3%) between the parents of two inter-genepool crosses. The majority of the microsatellites produced single bands and detected single loci, however four of the gene-based and three of the genomic microsatellites produced consistent double or multiple banding patterns and detected more than one locus. Microsatellite loci were found on each of the 11 chromosomes of common bean, the number per chromosome ranging from 5 to 17 with an average of ten microsatellites each. Total map length for the base map was 1,720 cM and the average chromosome length was 156.4 cM, with an average distance between microsatellite loci of 19.5 cM. The development of new microsatellites from sequences in the Genbank database and the implication of these results for genetic mapping, quantitative trait locus analysis and marker-assisted selection in common bean are described.Communicated by H.F. Linskens     

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.