Human microsatellites: mutation and evolution

Microsatellites (MSs) are short tandem DNA repeats with the repetitive motif of two to six nucleotides, forming tracts up to hundreds of nucleotides long. Notwithstanding the active use of MSs in genetic studies of various biological problems, the reasons for their wide occurrence in the genome, their possible functions, and mutational behavior are still unclear. The mutation rate in MS repeats is on average several orders of magnitude higher than in the remaining DNA, which allows for direct estimation of evolutionary transformation rate in nucleotide sequences of the genome. Mutation process in MSs is species-specific; furthermore, within a species it differs among loci with different repeat size, among alleles of one locus, and among individuals of different sex and age. Most MS mutations are caused by DNA slippage during replication but the probability of this event depends on the locus. In this review, a number of models of MS evolution are discussed, which account for the relationship between mutation rate and allele size, different mutation direction in alleles of different size, and the appearance of point mutations within repeat tracts restricting allele size. The MS evolution is considered mainly in the context of selective neutrality, although there is evidence showing functional significance of some variants of tandem repeats and thus their possible selective value.     

Human minisatellite alleles detectable only after PCR amplification.

We present evidence that a proportion of alleles at two human minisatellite loci is undetected by standard Southern blot hybridization. In each case the missing allele(s) can be identified after PCR amplification and correspond to tandem arrays too short to detect by hybridization. At one locus, there is only one undetected allele (population frequency 0.3), which contains just three repeat units. At the second locus, there are at least five undetected alleles (total population frequency 0.9) containing 60-120 repeats; they are not detected because these tandem repeats give very poor signals when used as a probe in standard Southern blot hybridization, and also cross-hybridize with other sequences in the genome. Under these circumstances only signals from the longest tandemly repeated alleles are detectable above the nonspecific background. The structures of these loci have been compared in human and primate DNA, and at one locus the short human allele containing three repeat units is shown to be an intermediate state in the expansion of a monomeric precursor allele in primates to high copy number in the longer human arrays. We discuss the implications of such loci for studies of human populations, minisatellite isolation by cloning, and the evolution of highly variable tandem arrays.     

Patterns of instability of expanded CAG repeats at the ERDA1 locus in general populations.

A highly polymorphic CAG repeat locus, ERDA1, was recently described on human chromosome 17q21.3, with alleles as large as 50-90 repeats and without any disease association in the general population. We have studied allelic distribution at this locus in five human populations and have characterized the mutational patterns by direct observation of 731 meioses. The data show that large alleles (>/=40 CAG repeats) are generally most common in Asian populations, less common in populations of European ancestry, and least common among Africans. We have observed a high intergenerational instability (46. 3%+/-5.1%) of the large alleles. Although the mutation rate is not dependent on parental sex, paternal transmissions have predominantly resulted in contractions, whereas maternal transmissions have yielded expansions. Within this class of large alleles, the mutation rate increases concomitantly with increasing allele size, but the magnitude of repeat size change does not depend on the size of the progenitor allele. Sequencing of specific alleles reveals that the intermediate-sized alleles (30-40 repeats) have CAT/CAC interruptions within the CAG-repeat array. These results indicate that expansion and instability of trinucleotide repeats are not exclusively disease-associated phenomena. The implications of the existence of massively expanded alleles in the general populations are not yet understood.     

Effects of temperature, Mg2+ concentration and mismatches on triplet-repeat expansion during DNA replication in vitro.

The human genome contains many simple tandem repeats that are widely dispersed and highly polymorphic. At least one group of simple tandem repeats, the DNA trinucleotide repeats, can dramaticallyexpand in size during transmission from one generation to the next to cause disease by a process known as dynamic mutation. We investigated the ability of trinucleotide repeats AAT and CAG to expand in size during DNA replication using a minimal in vitro system composed of the repeat tract, with and without unique flanking sequences, and DNA polymerase. Varying Mg2+concentration and temperature gave dramatic expansions of repeat size during DNA replication in vitro. Expansions of up to 1000-fold were observed. Mismatches partially stabilized the repeat tracts against expansion. Expansions were only detected when the primer was complementary to the repeat tract rather than the flanking sequence. The results imply that cellular environment and whether the growing strand contains a nick or gap are important factors for the expansion process in vivo.     

Tandemly repeated transgenes of the human minisatellite MS32 (D1S8), with novel mouse gamma satellite integration.

The human hypervariable minisatellite MS32 has a well characterised internal repeat unit array and high mutation rates have been observed at this locus. Analysis of MS32 mutants has shown that male germline mutations are polarised to one end of the array and frequently involve complex gene conversion-like events, suggesting that tandem repeat instability may be modulated by cis-acting sequences flanking the locus. In order to investigate the processes affecting MS32 mutation rate and mechanism, we have created transgenic mice harbouring an MS32 allele. Here we describe the organisation of eight transgenic insertions. Analysis of these transgenic loci by MVR-PCR and structural analysis of the junctions between mouse flanking DNA and the transgenic loci has shed light on mechanisms of integration and rearrangement of the tandem repeated transgenes. Sequence analysis of the mouse DNA flanking these transgenes has shown that 5 of the 8 insertions have integrated into mouse gamma satellite repeated sequence. This suggests a non-random integration of the MS32 transgene construct into the mouse genome.     

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.     

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.     

Microsatellite allele frequencies in humans and chimpanzees, with implications for constraints on allele size

The distributions of allele sizes at eight simple-sequence repeat (SSR) or microsatellite loci in chimpanzees are found and compared with the distributions previously obtained from several human populations. At several loci, the differences in average allele size between chimpanzees and humans are sufficiently small that there might be a constraint on the evolution of average allele size. Furthermore, a model that allows for a bias in the mutation process shows that for some loci a weak bias can account for the observations. Several alleles at one of the loci (Mfd 59) were sequenced. Differences between alleles of different lengths were found to be more complex than previously assumed. An 8-base-pair deletion was present in the nonvariable region of the chimpanzee locus. This locus contains a previously unrecognized repeated region, which is imperfect in humans and perfect in chimpanzees. The apparently greater opportunity for mutation conferred by the two perfect repeat regions in chimpanzees is reflected in the higher variance in repeat number at Mfd 59 in chimpanzees than in humans. These data indicate that interspecific differences in allele length are not always attributable to simple changes in the number of repeats.      

Mutation in microsatellite repeats of DNA and embryonal death in humans

In the analysis of tetranucleotide DNA repeats inheritance carried out in 55 families with a history of spontaneous miscarriages and normal karyotypes in respect to 21 loci located on seven autosomes, 8 embryos (14.5%) demonstrating 12 cases of the presence of alleles absent in both parents were described. The study of chromosome segregation using other DNA markers permitted highly probable exclusion of false paternity as well as uniparental disomy as the reasons for parent/child allele mismatches. The high probability of paternity together with the presence of a "new" allele at any offspring locus points to the mutation having occurred during game-togenesis in one of the parents. Examination of mutation in spontaneous abortuses revealed an increased number of tandem repeat units at microsatellite loci in three cases and an decreased number of these repeats in six cases. In two abortuses, a third allele absent in both parents, which resulted from a somatic mutation that occurred during embryonic development, was observed. The prevalence of the male germline mutations, revealed during investigation of the mutation origin, was probably associated with an increased number of DNA replication cycles in sperm compared to the oocytes. In spontaneous abortuses, the mean mutation rate of the tetranucleotide repeat complexes analyzed was 9.8 x 10(-3) per locus per gamete per generation. This was about five times higher than the spontaneous mutation rate of these STR loci. It can be suggested that genome instability detected at the level of repeated DNA sequences can involve not only genetically neutral loci but also active genomic regions crucial for embryonic viability. This results in cell death and termination of embryonic development. Our findings indicate that the death of embryos with normal karyotypes in most cases is associated with an increased frequency of germline and somatic microsatellite mutations. The data of the present study also provide a practical tool for the quantitative evaluation of this phenomenon and for the analysis of the reasons for miscarriages and embryonic death in certain families.     

Germline mutations of tetranucleotide DNA repeats in families with normal children and reproductive pathology

We have previously reported a high rate of tetranucleotide DNA repeat mutations, including mutations of both germline and somatic origin, in spontaneous human abortuses. To analyze in more detail mutational microsatellite (MS) variability in meiosis and its possible association with disturbed embryonic development, we have conducted a comparative study of mutation rates of a complex of 15 autosomal tetranucleotide MSs in 55 families with healthy children and in 103 families that have had spontaneous abortuses with normal karyotypes. In the families with miscarriage, the gametic MS mutation rate was higher than in the families with normal reproductive function (4.36 x 10(-3) versus 2.32 x 10(-3) per locus per gamete per generation), but this difference was statistically nonsignificant (P = 0.25). No association of MS mutations with familiar miscarriage was found. Mutations at the MS loci studied were recorded almost 3 times as often in spermatogenesis as in oogenesis, which is likely to result from a greater number of DNA replication cycles in male germline cell precursors than in female ones. Mutations increasing and reducing the MS sequence length appeared at virtually the same rate. Changes in MS DNA sequence length per one repeated element, i.e., single-step mutations (93% of cases) exceeded all other events of allele length change. The highest number of mutations (81.2%) was found in longer alleles. This distribution of mutations by size, direction, and parental origin corresponds to the multistep mutation model of their emergence via mechanism of DNA strand slippage during replication.     

Rate and pattern of mutation at microsatellite loci in maize

Microsatellites are important tools for plant breeding, genetics, and evolution, but few studies have analyzed their mutation pattern in plants. In this study, we estimated the mutation rate for 142 microsatellite loci in maize (Zea mays subsp. mays) in two different experiments of mutation accumulation. The mutation rate per generation was estimated to be 7.7 x 10(-4) for microsatellites with dinucleotide repeat motifs, with a 95% confidence interval from 5.2 x 10(-4) to 1.1 x 10(-3). For microsatellites with repeat motifs of more than 2 bp in length, no mutations were detected; so we could only estimate the upper 95% confidence limit of 5.1 x 10(-5) for the mutation rate. For dinucleotide repeat microsatellites, we also determined that the variance of change in the number of repeats (sigma(m)2) is 3.2. We sequenced 55 of the 73 observed mutations, and all mutations proved to be changes in the number of repeats in the microsatellite or in mononucleotide tracts flanking the microsatellite. There is a higher probability to mutate to an allele of larger size. There is heterogeneity in the mutation rate among dinucleotide microsatellites and a positive correlation between the number of repeats in the progenitor allele and the mutation rate. The microsatellite-based estimate of the effective population size of maize is more than an order of magnitude less than previously reported values based on nucleotide sequence variation.     

High somatic instability of a microsatellite locus in a clonal tree,<Emphasis Type="Italic"> Robinia pseudoacacia</Emphasis>

Robinia pseudoacacia L. is a clonal tree species. To investigate a mutation within eight microsatellite loci of R. pseudoacacia, we analyzed DNA samples obtained from different leaf samples within each ramet, leaves from ramets within the genet, and seeds. Of the eight loci, locus Rops15 (AG motif) displayed hypermutability. The mutation rates of Rops15 within each ramet, among ramets within the genet, and offspring were 6.27% (ranging from 0 to 31.1%), 6.11% (from 0 to 25.0%) and 3.78% (from 0 to 10.9%), respectively. The mutation rate increased with allele size (13–71 repeat units). The mutation patterns observed in Rops15 were distinctive in two ways. First, there was a significant bias toward additions over deletions, and both addition and deletion of single repeats were dominant at alleles with lengths less than 232 bp (63 repeats). Second, for the longest allele of 248 bp (71 repeats), the number of losses was higher than the number of gains. These observations suggest that the mutation patterns of microsatellites in R. pseudoacacia may follow a generalized stepwise mutation model, and that the tendency of long alleles to mutate to shorter lengths acts to prevent infinite growth. Finally, the observation of somatic hypermutability at locus Rops15 highlights the need for caution when using highly polymorphic microsatellites for population genetic structure and paternity analysis in tree species.Communicated by H.F. Linskens     

A variable tandem repeat locus mapped to chromosome band 10q26 is amplified and rearranged in leukocyte DNAs of two cancer patients.

A highly polymorphic locus associated with the variable tandem repetition of a 35 bp consensus sequence was mapped to chromosome 10, band q26. Examination of leukocyte DNA from a cancer patient revealed the twenty-fold amplification of one allelic fragment of this locus, while the other allelic fragment demonstrated a normal copy number. In another patient, Southern blotting of leukocyte DNA detected the deletion of the 3'-flanking region from one tandem repeat allele. These results indicate that variable tandem repeats may mark highly unstable regions of DNA in the human genome which can be altered by changes more extensive than simple tandem repeat variation.     

Rapid, high fidelity analysis of simple sequence repeats on an electronically active DNA microchip

We describe a method for the discrimination of short tandem repeat (STR) alleles based on active microarray hybridization. An essential factor in this method is electronic hybridization of the target DNA, at high stringency, in <5 min. High stringency is critical to avoid slippage of hybrids along repeat tracts at allele-specific test sites in the array. These conditions are attainable only with hybridization kinetics realized by electronic concentration of DNA. A sandwich hybrid is assembled, in which proper base stacking of juxtaposed terminal nucleotides results in a thermodynamically favored complex. The increased stability of this complex relative to non-stacked termini and/or base pair mismatches is used to determine the identification of STR alleles. This method is capable of simultaneous and precise identification of alleles containing different numbers of repeats, as well as mutations within these repeats. Given the throughput capabilities of microarrays our system has the potential to enhance the use of microsatellites in forensic criminology, diagnostics and genetic mapping.     

Mutation rate in the hypervariable VNTR g3 (D7S22) is affected by allele length and a flanking DNA sequence polymorphism near the repeat array.

The hypervariable human minisatellite locus D7S22 (g3) is highly polymorphic. The allelic distribution in D7S22 features a size clustering of the alleles and a comparably low allelic diversity among small alleles. This reduced diversity could reflect a situation where some alleles are less likely to mutate than others. Several factors could explain such an effect, including allele size, variation in repeat composition, and allelic differences in nearby cis-acting elements affecting the mutation rate. We have characterized 40 de novo mutations found on Southern blots in a large amount of paternity-testing material. There is a significant excess of paternal mutations, and small size changes are most frequent. Mutation rate is affected by allele length, with highest rates in larger alleles. Alleles of the family groups with D7S22 mutations and 50 small alleles were analyzed by nucleotide sequencing. Two hundred thirty-six base pairs of the immediate flanking region upstream of the repeat array were PCR amplified and screened for point mutations by DNA sequencing of the PCR products. Two base substitution polymorphisms were identified: one C/G transversion and one A/G transition, 54 bp and 173 bp upstream of the repeat array, respectively. There is a significant association between mutation and occurrence of 54C, while association is not obvious between mutation rate and the 173A/G variants. There is a marked association between different flanking haplotypes and allele size, and within the smallest allele-size group, all alleles had the 54G/173A haplotype. Both allele size and allelic state at site 54 remain associated with mutation rate when the other factor is controlled. Possible mechanisms behind the variation in mutation rate in D7S22 are discussed.     

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.     

VNTRs in review

In the last decade the study of human genetic variation has made a quantum leap from the analysis of protein and antigen intermediaries to the investigation of DNA itself.¹ The DNA double helix codes genetic information as a sequence of four different nucleotides: adenine (A), guanine (G), cytosine (C), and thymine (T). Nucleotides are nitrogenous bases that bind the complementary strands of the double helix, giving rise to the use of base pairs (bp) as a unit of DNA length. So far so good. Within the human genome there are DNA sequences that do not code for proteins and that consist of short runs of nucleotides, say GTGGACAGG, repeated in tandem hundreds, or even thousands of times. This particular sequence, known as MS1 for minisatellite 1, is found on human chromosome 1 at a locus called D1S7. Although it is old news that there are a lot of repetitive DNAs in the human genome, it is new and very interesting to find that many repetitive DNA loci have arrays of different numbers of repeats in different individuals. These loci are referred to as VNTRs, shorthand for “variable numbers of tandem repeats” or, more flippantly for “very nasty types of repeat.” The finding of hundreds of VNTR loci distributed across all chromosomes exposes a richness of genetic diversity for anthropologists studying human evolutionary history.     

RecJ, ExoI and RecG are required for genome maintenance but not for generation of genetic diversity by repeat-mediated phase variation in Haemophilus influenzae

High levels of genetic diversity are generated in Haemophilus influenzae populations through DNA repeat-mediated phase variation and recombination with DNA fragments acquired by uptake from the external milieu. Conversely, multiple pathways for maintenance of the genome sequence are encoded in H. influenzae genomes. In Escherichia coli, mutations in single-stranded-DNA exonucleases destabilise tandem DNA repeats whilst inactivation of recG can stabilise repeat tracts. These enzymes also have varying effects on recombination. Deletion mutations were constructed in H. influenzae genes encoding homologs of ExoI, RecJ and RecG whilst ExoVII was refractory to mutation. Inactivation of RecJ and RecG, but not ExoI, increased sensitivity to irradiation with ultraviolet light. An increase in spontaneous mutation rate was not observed in single mutants but only when both RecJ and ExoI were mutated. None of the single- or double-mutations increased or decreased the rates of slippage in tetranucleotide repeat tracts. Furthermore, the exonuclease mutants did not exhibit significant defects in horizontal gene transfer. We conclude that RecJ, ExoI and RecG are required for maintenance of the H. influenzae genome but none of these enzymes influence the generation of genetic diversity through mutations in the tetranucleotide repeat tracts of this species.     

Allelic and population variation of microsatellite loci in aspen (Populus tremuloides)

To develop a robust basis for inferences about population genetics and evolution, this work assayed 192 aspens (Populus tremuloides) from 11 sites in Wisconsin, USA, for allelic and population variation at 16 microsatellite loci distributed across the Populus genome. Frequency distributions of fluorochrome-labeled alleles resolved by capillary electrophoresis were analyzed for relationships to repeat size and number. Population-level statistics were compared with those of other studies, especially in Populus. All loci were polymorphic, varying widely in the number of alleles per locus (mean = 8.25, range 2-20). Expected and observed heterozygosities were high (0.45 and 0.41, respectively), with little differentiation among populations (F(ST) = 0.006-0.045) and a moderate level of inbreeding (F(IS) = 0.09), intermediate among levels reported in studies based on isozymes. Contrary to several other reports, allele frequencies clustered tightly around the modal frequency, and the genetic diversity (measured as alleles per locus or as expected heterozygosity) was not related to either the repeat unit size or to the number of repeats.