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 (AC)₂ with those around other four-base combinations reveals that (AC)₂ does not stand out as being special in the sense that non-repetitive tetramers also generate strong mutation biases. However, comparing (AC)₂ and (AC)₃ with (AC)₄ reveals a step change in both the rate and nature of mutations occurring, suggesting a transition state, (AC)₄ 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.     

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.     

Simple sequence repeats in mycobacterial genomes

Simple sequence repeats (SSRs) or microsatellites are the repetitive nucleotide sequences of motifs of length 1–6 bp. They are scattered throughout the genomes of all the known organisms ranging from viruses to eukaryotes. Microsatellites undergo mutations in the form of insertions and deletions (INDELS) of their repeat units with some bias towards insertions that lead to microsatellite tract expansion. Although prokaryotic genomes derive some plasticity due to microsatellite mutations they have in-built mechanisms to arrest undue expansions of microsatellites and one such mechanism is constituted by post-replicative DNA repair enzymes MutL, MutH and MutS. The mycobacterial genomes lack these enzymes and as a null hypothesis one could expect these genomes to harbour many long tracts. It is therefore interesting to analyse the mycobacterial genomes for distribution and abundance of microsatellites tracts and to look for potentially polymorphic microsatellites. Available mycobacterial genomes, Mycobacterium avium, M. leprae, M. bovis and the two strains of M. tuberculosis (CDC1551 and H37Rv) were analysed for frequencies and abundance of SSRs. Our analysis revealed that the SSRs are distributed throughout the mycobacterial genomes at an average of 220–230 SSR tracts per kb. All the mycobacterial genomes contain few regions that are conspicuously denser or poorer in microsatellites compared to their expected genome averages. The genomes distinctly show scarcity of long microsatellites despite the absence of a post-replicative DNA repair system. Such severe scarcity of long microsatellites could arise as a result of strong selection pressures operating against long and unstable sequences although influence of GC-content and role of point mutations in arresting microsatellite expansions can not be ruled out. Nonetheless, the long tracts occasionally found in coding as well as non-coding regions may account for limited genome plasticity in these genomes. Supplementary Data pertaining to this article is available on the Journal of Biosciences Website at     

Co-polymer tracts in eukaryotic, prokaryotic, and organellar DNA.

Large variations in DNA base composition and noticeable strand asymmetries are known to occur between different organisms and within different regions of the genomes of single organisms. Apparently such composition and sequence biases occur to fulfill structural rather than informational requirements. Here we report the wide occurrence of a more subtle biasing of DNA sequence that can have structural consequences: an increase or a suppression of the number of long tracts of two-base co-polymers. Strong biases were observed when the DNA sequences of the longest eukaryotic, prokaryotic, and organellar entries in the GenBank data base (totaling 773 kilobases) were analyzed for the number of occurrences of tracts of the two-base co-polymers (A,T)n, (G,C)n, and (A,C)n as a function of tract length. (The expression (A,T)n is used here to denote an uninterrupted tract, n nucleotides in length, of A and T bases in any proportion or order, terminated at each end by a G or C residue.) Characteristic differences are also observed in tract biases of eukaryotic vs. prokaryotic organisms.     

Rapid divergence of microsatellite abundance among species of Drosophila

Among major taxonomic groups, microsatellites exhibit considerable variation in composition and allele length, but they also show considerable conservation within many major groups. This variation may be explained by slow microsatellite evolution so that all species within a group have similar patterns of variation, or by taxon-specific mutational or selective constraints. Unfortunately, comparing microsatellites across species and studies can be problematic because of biases that may exist among different isolation and analysis protocols. We present microsatellite data from five Drosophila species in the Drosophila subgenus: D. arizonae, D. mojavensis, and D. pachea (three cactophilic species), and D. neotestacea and D. recens (two mycophagous species), all isolated at the same time using identical protocols. For each species, we compared the relative abundance of motifs, the distribution of repeat size, and the average number of repeats. Dimers were the most abundant microsatellites for each species. However, we found considerable variation in the relative abundance of motif size classes among species, even between sister taxa. Frequency differences among motifs within size classes for the three cactophilic species, but not the two mycophagous species, are consistent with other studied Drosophila. Frequency distributions of repeat number, as well as mean size, show significant differences among motif size classes but not across species. Sizes of microsatellites in these five species are consistent with D. virilis, another species in the subgenus Drosophila, but they have consistently higher means than in D. melanogaster, in the subgenus Sophophora. These results confirm that many aspects of microsatellite variation evolve quickly but also are subject to taxon-specific constraints. In addition, the nature of microsatellite evolution is dependent on temporal and taxonomic scales, and some variation is conserved across broad taxonomic levels despite relatively high rates of mutation for these loci.     

Motif mismatches in microsatellites: insights from genome-wide investigation among 20 insect species

We present a detailed genome-wide comparative study of motif mismatches of microsatellites among 20 insect species representing five taxonomic orders. The results show that varying proportions (∼15–46%) of microsatellites identified in these species are imperfect in motif structure, and that they also vary in chromosomal distribution within genomes. It was observed that the genomic abundance of imperfect repeats is significantly associated with the length and number of motif mismatches of microsatellites. Furthermore, microsatellites with a higher number of mismatches tend to have lower abundance in the genome, suggesting that sequence heterogeneity of repeat motifs is a key determinant of genomic abundance of microsatellites. This relationship seems to be a general feature of microsatellites even in unrelated species such as yeast, roundworm, mouse and human. We provide a mechanistic explanation of the evolutionary link between motif heterogeneity and genomic abundance of microsatellites by examining the patterns of motif mismatches and allele sequences of single-nucleotide polymorphisms identified within microsatellite loci. Using Drosophila Reference Genetic Panel data, we further show that pattern of allelic variation modulates motif heterogeneity of microsatellites, and provide estimates of allele age of specific imperfect microsatellites found within protein-coding genes.     

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.     

Occurrence and analysis of imperfect microsatellites in diverse potyvirus genomes

Simple sequence repeats (SSRs) or microsatellites are known to exhibit ubiquitous across all kingdoms of life including viruses. However, imperfections in simple sequence repeats have been analyzed in genomes of human, Escherichia coli and Human Immunodeficiency virus. The assessment of compound microsatellites in plant viral genomes is yet to be studied. Potyviruses severely affect crop plant growth and reduce economic yield in diverse cropping systems worldwide. Hence, we analyze the nature and distribution of compound microsatellites present in complete genome of 45 potyvirus species. The results indicate that compound microsatellites accounted for about 0% to 15.15% of all microsatellites and have low complexity as compared to that of prokaryotic genomes. Overall, 14% of compound microsatellites were of similar motifs and such motif duplications were observed for CA, TA and AG repeats. Among all 45 potyvirus genomes analyzed, SSR couple (AG)-x-(AC) was found to be the most abundant one. Hence it is apparent that in contrast to eukaryotes, majority of compound microsatellites in potyviruses were composed of variant motifs. We also highlight the relative frequency of different classes of compound microsatellites as well as their patterns of distribution and correlate with biology of potyviruses. Further characterization of such variation is important for elucidating the origin, mutational processes, and structure of these widely used, but incompletely understood sequences.     

Complex evolution of a highly conserved microsatellite locus in several fish species

The evolutionary dynamics of a highly conserved microsatellite locus ( Dla 11 ) were studied in several fish species. The data indicated that multiple types of compound microsatellites arose through point mutations that were sometimes followed by expansion of the derived motif. Furthermore, extensive length variation was detected among species in the regions immediately flanking the repeat region.     

A new genotyping method for detecting low abundance single nucleotide mutations based on gap ligase chain reaction and quantitative PCR assay

We tested applicability of a new genotyping technique to detect a low abundance CD17 (A → T) mutation of β-globin gene. The technique utilized a combined gap ligase chain reaction (Gap-LCR) and quantitative PCR (qPCR) methods. One pair of Gap-LCR primers was modified by adding specific sequences to the 5′ end of the upstream and the 3′ end of the downstream primer which served as a combining sequence for qPCR. First, specific mutation is detected using Gap-LCR; then, ligation products are detected by qPCR. Our results show that the amount of LCR products is directly proportional to the amount of template DNA. We further demonstrate that this technique detects a low abundance mutant DNA with a mutant/normal allele ratio as low as 1:10000. This technique was applied to detect a paternally inherited CD17 mutation from 53 maternal plasma samples. The results were consistent with those obtained by PCR/reverse dot blot of amniotic fluid cell DNA. In conclusion, by combining Gap-LCR and qPCR technology we successfully established a highly sensitive technique to detect low abundance point mutations. This technique can be applied to detect fetal DNA point mutation in maternal plasma.     

Development and application of microsatellite markers for genomic analysis of papaya

Papaya has a relatively small genome, displays high levels of phenotypic diversity, and is amenable to transformation, making it attractive as a fruit tree model system. The high level of phenotypic diversity seen among papaya cultivars in the field does not correlate with the low levels of genotypic polymorphism thus far elucidated. The highly mutable nature of microsatellites or simple sequence repeats (SSRs) make them potentially powerful markers for distinguishing deoxyribonucleic acid (DNA) polymorphisms between closely related genotypes. Genomic research for papaya has resulted in a significant quantity of sequence data. We mined 28.1 Mb of bacterial artificial chromosomes end sequences, 5.8 Mb of complementary DNA, and 1.6 Mb of random genomic sequences for SSRs. We generated 938 SSR markers and tested for polymorphism among seven varieties that had been used to produce five mapping populations. The level of polymorphism was highest for Kaek Dum × 2H94 with 210 markers, followed by UH928 × SunUp with 194, AU9 × SunUp with 189, UH918 × SunUp with 177, and Kapoho × SunUp displaying the lowest level with 97. Variation in levels of polymorphism, motif predominance, and motif length between the genomic and genic fractions indicated differential selection pressures acting on the microsatellites in these two fractions. The microsatellites developed in this study will greatly assist in the genetic and physical mapping of the papaya genome as well as enhance breeders’ ability to improve the crop. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

PKD1 intron 21: triplex DNA formation and effect on replication

Although autosomal dominant polycystic kidney disease is transmitted in an autosomal dominant fashion, there is evidence that the pathophysiology of cystogenesis involves a second hit somatic mutation superimposed upon the inherited germline mutation within the renal tubule cells. The polypurine·polypyrimidine (Pu·Py) tract of PKD1 intron 21 may play a role in promoting somatic mutations. To better characterize this tract and to evaluate its potential to participate in mutagenesis, we investigated the thermodynamics of intramolecular triplex formation by 15 Pu·Py mirror repeat tracts from PKD1 intron 21 by 2D gel electrophoresis. We demonstrate that intramolecular triplexes form with modest superhelical tensions for all the tracts examined. Primer extension studies demonstrated significant polymerase arrest within the Pu·Py tracts in one direction of replication only. We found correlation between polymerization arrest and both the potential length of the triplex and superhelical tension of intramolecular triplex formation. The presence of a Pu·Py tract also led to a replication blockade and double-strand breakage using an SV40 in vitro replication assay with HeLa cell extracts. During DNA replication, the G-rich template of the PKD1 Pu·Py tracts may form a triplex structure with the nascent strand, thereby blocking replication and potentially leading to recombination and mutation.     

Microsatellite instability in yeast: dependence on repeat unit size and DNA mismatch repair genes.

We examined the stability of microsatellites of different repeat unit lengths in Saccharomyces cerevisiae strains deficient in DNA mismatch repair. The msh2 and msh3 mutations destabilized microsatellites with repeat units of 1, 2, 4, 5, and 8 bp; a poly(G) tract of 18 bp was destabilized several thousand-fold by the msh2 mutation and about 100-fold by msh3. The msh6 mutations destabilized microsatellites with repeat units of 1 and 2 bp but had no effect on microsatellites with larger repeats. These results argue that coding sequences containing repetitive DNA tracts will be preferred target sites for mutations in human tumors with mismatch repair defects. We find that the DNA mismatch repair genes destabilize microsatellites with repeat units from 1 to 13 bp but have no effect on the stability of minisatellites with repeat units of 16 or 20 bp. Our data also suggest that displaced loops on the nascent strand, resulting from DNA polymerase slippage, are repaired differently than loops on the template strand.     

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 (≤0.51), which could be caused by long null alleles. [Reviewing Editor: Dr. Dmitri Petrov] Sequence data from this article have been deposited with the EMBL/GenBank databases under accession numbers AY532396–AY532400.     

The structure of interrupted human AC microsatellites

Microsatellite lengths change over evolutionary time through a process of replication slippage. A recently proposed model of this process holds that the expansionary tendencies of slippage mutation are balanced by point mutations breaking longer microsatellites into smaller units and that this process gives rise to the observed frequency distributions of uninterrupted microsatellite lengths. We refer to this as the slippage/point-mutation theory. Here we derive the theory's predictions for interrupted microsatellites comprising regions of perfect repeats, labeled segments, separated by dinucleotide interruptions containing point mutations. These predictions are tested by reference to the frequency distributions of segments of AC microsatellite in the human genome, and several predictions are shown not to be supported by the data, as follows. The estimated slippage rates are relatively low for the first four repeats, and then rise initially linearly with length, in accordance with previous work. However, contrary to expectation and the experimental evidence, the inferred slippage rates decline in segments above 10 repeats. Point mutation rates are also found to be higher within microsatellites than elsewhere. The theory provides an excellent fit to the frequency distribution of peripheral segment lengths but fails to explain why internal segments are shorter. Furthermore, there are fewer microsatellites with many segments than predicted. The frequencies of interrupted microsatellites decline geometrically with microsatellite size measured in number of segments, so that for each additional segment, the number of microsatellites is 33.6% less. Overall we conclude that the detailed structure of interrupted microsatellites cannot be reconciled with the existing slippage/point-mutation theory of microsatellite evolution, and we suggest that microsatellites are stabilized by processes acting on interior rather than on peripheral segments.     

Microsatellite variability differs between dinucleotide repeat motifs-evidence from Drosophila melanogaster

Recently, the use of microsatellites as genetic markers has become very popular. While their evolutionary dynamics are not yet fully understood, the emerging picture is that several factors are influencing microsatellite mutation rates. Recent experiments demonstrated a significant effect of repeat motif length on microsatellite mutation rates. Here, we studied the influence of the base composition of the microsatellite. Forty-two microsatellite loci on the second chromosome with the three most abundant dinucleotide repeat motifs (TC/AG, AT/TA, GT/CA) were characterized for six different Drosophila melanogaster populations. Applying ANOVA to the variance in repeat number, we found a significant influence of repeat motif on microsatellite variability. Calculating relative mutation rates, GT/CA appears to have the highest mutation rate, and AT/TA appears to have the lowest. Similar differences in mutation rates were obtained by an alternative method which estimates microsatellite mutation rates from their genomic length distribution.     

Use of polymorphic short and clustered coding-region microsatellites to distinguish strains of Candida albicans

Abstract We describe the identification of polymorphic microsatellite loci in the pathogenic yeast, Candida albicans . A search for all coding-region microsatellites with more than four repeats that can be found in Candida sequences in GenBank was conducted. Nine such microsatellite sequences consisting of trinucleotide motifs were found. Three of these were perfect microsatellites while the remaining six sequences were found in one imperfect microsatellite and two compound microsatellites. Because of the close proximity of some of these repeats, all could be assayed with six PCR primer pairs. All of these microsatellite sequences were found in five nuclear genes, ZNF1, CCN1, CPH1, EFG1 , and MNT2 . Except for a single (CTT)₅ serine tract, all coded for polyglutamine tracts. Another locus with seven alleles, a region of the ERK1 protein kinase gene, was also examined, and may be a representative of a new class of highly polymorphic ‘clustered’ microsatellites. Such loci, in which several non-contiguous but closely linked microsatellites are clustered together, may be a useful source of DNA polymorphisms in microorganisms in which long microsatellite sequences are unavailable. All seven regions amplified were polymorphic, having between two and seven variable length alleles in the 11 strains of Candida albicans examined. The results of this and similar searches will facilitate epidemiological and evolutionary studies of Candida and other microorganisms.     

Microsatellite Allelic Homoplasy Due to Variable Flanking Sequences

Microsatellite DNA sequences have become the dominant source of nuclear genetic markers for most applications. It is important to investigate the basis of variation between alleles and to know if current assumptions about the mechanisms of microsatellite mutation (that is to say, variations involving simple changes in the number of repeat) are correct. We have characterized, by DNA sequencing, the human alleles of a new highly informative (CA)n repeat localized approximately 20 kb centromeric to the HLA-B gene. Although 12 alleles were identified based on conventional length criteria, sequencing of the alleles demonstrated that differences between alleles were found to be more complex than previously assumed: A high degree of microsatellite variability is due to variation in the region immediately flanking the repeat. These data indicate that the mutational process which generates polymorphism in this region has involved not only simple changes in the number of dinucleotide CA repeats but also perturbations in the nonrepeated 5′ and 3′ flanking sequences. Three families of alleles (not visible from the overall length of the alleles), with presumably separate evolutionary histories, exist and can yield to homoplasy of size. Effectively, we can observe alleles of the same size with different internal structures which are separated by a significant amount of variation. Although allelic homoplasy for noninterrupted microsatellite loci has been suggested between different species, it has not been unequivocally demonstrated within species. A strong association is noted between alleles defined at the sequence level and HLA-B alleles. The observation of several families of alleles at the population level provides information about the evolutionary history and mutation processes of microsatellites and may have implications for the use of these markers in phylogenetic, linkage disequilibrium studies, and gene mapping. Received: 14 May 1996 / Accepted: 9 September 1996     

Using the recognition code to swap homeodomain target specificity in cell culture

The homeodomain (HD) is a 60 amino acid-long DNA-binding domain. A large fraction of HDs binds with high affinity sequences containing the 5′-TAAT-3′ core motif. However, NK-2 class HDs recognizes sequences containing the 5′-CAAG-3′ core motif. By using a cell transfection approach, here we show that modification of residues located in the N-terminal arm (at positions 6, 7 and 8) and in the recognition helix (at position 54) is enough to swap the “in vivo” binding specificity of TTF-1 HD (which is a member of the NK-2 class HD) from 5′-CAAG-3′ to 5′-TAAT-3′-containing targets. The role of residue at position 54 is also supported by data obtained with the HD of the Drosophila engrailed protein. These data support the notion that DNA-binding specificity “in vivo” is dictated by few critical residues.