Microsatellite and trinucleotide-repeat evolution: evidence for mutational bias and different rates of evolution in different lineages

Genomic Variation, Evolvability, and the Paradox of Mental Illness

Twentieth-century genetics was hard put to explain the irregular behavior of neuropsychiatric disorders. Autism and schizophrenia defy a principle of natural selection; they are highly heritable but associated with low reproductive success. Nevertheless, they persist. The genetic origins of such conditions are confounded by the problem of variable expression, that is, when a given genetic aberration can lead to any one of several distinct disorders. Also, autism and schizophrenia occur on a spectrum of severity, from mild and subclinical cases to the overt and disabling. Such irregularities reflect the problem of missing heritability; although hundreds of genes may be associated with autism or schizophrenia, together they account for only a small proportion of cases. Techniques for higher resolution, genomewide analysis have begun to illuminate the irregular and unpredictable behavior of the human genome. Thus, the origins of neuropsychiatric disorders in particular and complex disease in general have been illuminated. The human genome is characterized by a high degree of structural and behavioral variability: DNA content variation, epistasis, stochasticity in gene expression, and epigenetic changes. These elements have grown more complex as evolution scaled the phylogenetic tree. They are especially pertinent to brain development and function. Genomic variability is a window on the origins of complex disease, neuropsychiatric disorders, and neurodevelopmental disorders in particular. Genomic variability, as it happens, is also the fuel of evolvability. The genomic events that presided over the evolution of the primate and hominid lineages are over-represented in patients with autism and schizophrenia, as well as intellectual disability and epilepsy. That the special qualities of the human genome that drove evolution might, in some way, contribute to neuropsychiatric disorders is a matter of no little interest.

Tandem repeats mediating genetic plasticity in health and disease

Accumulating evidence suggests that many classes of DNA repeats exhibit attributes that distinguish them from other genetic variants, including the fact that they are more liable to mutation; this enables them to mediate genetic plasticity. The expansion of tandem repeats, particularly of short tandem repeats, can cause a range of disorders (including Huntington disease, various ataxias, motor neuron disease, frontotemporal dementia, fragile X syndrome and other neurological disorders), and emerging data suggest that tandem repeat polymorphisms (TRPs) can also regulate gene expression in healthy individuals. TRPs in human genomes may also contribute to the missing heritability of polygenic disorders. A better understanding of tandem repeats and their associated repeatome, as well as their capacity for genetic plasticity via both germline and somatic mutations, is needed to transform our understanding of the role of TRPs in health and disease.

Microsatellite Mutation Rate in Atlantic Sturgeon (Acipenser oxyrinchus)

Understanding mutation rates can greatly extend the utility of population and conservation genetic analyses. Herein, we present an estimate of genome-wide microsatellite mutation rate in Atlantic sturgeon (Acipenser oxyrinchus) based on parent-offspring transmission patterns. We screened 307 individuals for parentage and mutation-rate analysis applying 43 variable markers. Out of 13228 allele transfers, 11 mutations were detected, producing a mutation rate of 8.3 × 10-4 per locus per generation (95% confidence interval: 1.48 × 10-3, 4.15 × 10-4). Single-step mutations predominated and there were trends toward mutations in loci with greater polymorphism and allele length. Two of the detected mutations were most probably cluster mutations, being identified in 12 and 28 sibs, respectively. Finally, we observed evidences of polyploidy based on the sporadic presence of 3 or 4 alleles per locus in the genotyped individuals, supporting previous reports of incomplete diploidization in Atlantic sturgeon.

Genome-wide comparative analysis of simple sequence coding repeats among 25 insect species

We present a detailed genome-scale comparative analysis of simple sequence repeats within protein coding regions among 25 insect genomes. The repetitive sequences in the coding regions primarily represented single codon repeats and codon pair repeats. The CAG triplet is highly repetitive in the coding regions of insect genomes. It is frequently paired with the synonymous codon CAA to code for polyglutamine repeats. The codon pairs that are least repetitive code for polyalanine repeats. The frequency of hexanucleotide and dinucleotide motifs of codon pair repeats is significantly (p<0.001) different in the Drosophila species compared to the non-Drosophila species. However, the frequency of synonymous and non-synonymous codon pair repeats varies in a correlated manner (r(2)=0.79) among all the species. Results further show that perfect and imperfect repeats have significant association with the trinucleotide and hexanucleotide coding repeats in most of these insects. However, only select species show significant association between the numbers of perfect/imperfect hexamers and repeat coding for single amino acid/amino acid pair runs. Our data further suggests that genes containing simple sequence coding repeats may be under negative selection as they tend to be poorly conserved across species. The sequences of coding repeats of orthologous genes vary according to the known phylogeny among the species. In conclusion, the study shows that simple sequence coding repeats are important features of genome diversity among insects.

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.

Population genetic structure of Aldabra giant tortoises

Evolution of population structure on islands is the result of physical processes linked to volcanism, orogenic events, changes in sea level, as well as habitat variation. We assessed patterns of genetic structure in the giant tortoise of the Aldabra atoll, where previous ecological studies suggested population subdivisions as a result of landscape discontinuity due to unsuitable habitat and island separation. Analysis of mitochondrial DNA (mtDNA) control region sequences and allelic variation at 8 microsatellite loci were conducted on tortoises sampled in 3 locations on the 2 major islands of Aldabra. We found no variation in mtDNA sequences. This pattern corroborated earlier work supporting the occurrence of a founding event during the last interglacial period and a further reduction in genetic variability during historical time. On the other hand, significant population structure recorded at nuclear loci suggested allopatric divergence possibly due to geographical barriers among islands and ecological partitions hindering tortoise movements within islands. This is the first attempt to study the population genetics of Aldabra tortoises, which are now at carrying capacity in an isolated terrestrial ecosystem where ecological factors appear to have a strong influence on population dynamics.

Polyglutamine gene function and dysfunction in the ageing brain

The coordinated regulation of gene expression and protein interactions determines how mammalian nervous systems develop and retain function and plasticity over extended periods of time such as a human life span. By studying mutations that occur in a group of genes associated with chronic neurodegeneration, the polyglutamine (polyQ) disorders, it has emerged that CAG/glutamine stretches play important roles in transcriptional regulation and protein-protein interactions. However, it is still unclear what the many structural and functional roles of CAG and other low-complexity sequences in eukaryotic genomes are, despite being the most commonly shared peptide fragments in such proteomes. In this review we examine the function of genes responsible for at least 10 polyglutamine disorders in relation to the nervous system and how expansion mutations lead to neuronal dysfunction, by particularly focusing on Huntington's disease (HD). We argue that the molecular and cellular pathways that turn out to be dysfunctional during such diseases, as a consequence of a CAG expansion, are also involved in the ageing of the central nervous system. These are pathways that control protein degradation systems (including molecular chaperones), axonal transport, redox-homeostasis and bioenergetics. CAG expansion mutations confer novel properties on proteins that lead to a slow-progressing neuronal pathology and cell death similar to that found in other age-related conditions such as Alzheimer's and Parkinson's diseases.

Dynamic mutations as digital genetic modulators of brain development, function and dysfunction

A substantial portion of the human genome has been found to consist of simple sequence repeats, including microsatellites and minisatellites. Microsatellites, tandem repeats of 1-6 nucleotides, form the template for dynamic mutations, which involve heritable changes in the lengths of repeat sequences. In recent years, a large number of human disorders have been found to be caused by dynamic mutations, the most common of which are trinucleotide repeat expansion diseases. Dynamic mutations are common to numerous nervous system disorders, including Huntington's disease, various spinocerebellar ataxias, fragile X syndrome, fragile X tremor/ataxia syndrome, Friedreich ataxia and other neurodegenerative disorders. The involvement of dynamic mutations in brain disorders will be reviewed, with a focus on the large group caused by CAG/glutamine repeat expansions. We will also outline a proposed role of tandem repeat polymorphisms (TRPs), with unique 'digital' genetic distributions, in modulating brain development and normal function, so as to generate additional mutational diversity upon which natural selection may act.

Contrasting patterns of mitochondrial and microsatellite genetic structure among Western European populations of tawny owls (Strix aluco)

A recent study of mitochondrial phylogeography of tawny owls (Strix aluco) in western Europe suggested that this species survived the Pleistocene glaciations in three allopatric refugia located in Iberia, Italy, and the Balkans, and the latter was likely the predominant source of postglacial colonization of northern Europe. New data from seven microsatellite loci from 184 individual owls distributed among 14 populations were used to assess the genetic congruence between nuclear and mitochondrial DNA (mtDNA) markers. Microsatellites corroborated the major phylogeographical conclusions reached on the basis of the mtDNA sequences, but also showed important differences leading to novel inferences. Microsatellites corroborated the three major refugia and supported the Balkan origin of northern populations. When corrected for differences in effective population size, microsatellites and mtDNA yielded generally congruent overall estimates of population structure (N*ST=0.12 vs. RST=0.16); however, there was substantial heterogeneity in the RST among the seven nuclear loci that was not correlated with heterozygosity. Populations representing the Balkans postglacial expansion interact with populations from the other two refugia forming two clines near the Alps and the Pyrenees. In both cases, the apparent position of the contact zones differed substantially between markers due to the genetic composition of populations sampled in northern Italy and Madrid. Microsatellite data did not corroborate the lower genetic diversity of northern, recently populated regions as was found with mtDNA; this discrepancy was taken as evidence for a recent bottleneck recovery. Finally, this study suggests that congruence among genetic markers should be more likely in cases of range expansion into new areas than when populations interact across contact zones.

Repeat expansion in spinocerebellar ataxia type 17 alleles of the TATA-box binding protein gene: an evolutionary approach

The variability and mutational changes of the CAG microsatellite in the TATA-box binding protein gene (TBP) were studied. We sequenced the microsatellite of the TBP gene of 25 unrelated individuals from northern Germany (10 SCA17 patients and 15 unaffected control individuals). In addition, the microsatellites were sequenced from individuals of 10 northern German families with at least one family member affected by SCA17. To study also the evolutionary history of this CAG/CAA microsatellite in nonhuman primates, the homologous regions were analysed from Pan troglodytes, Gorilla gorilla, Pongo pygmaeus, P. abellii, Hylobates lar, Nomascus leucogenys, Symphalangus syndactylus, Macaca mulatta, Papio hamadryas, Colobus polykomos and Callithrix jacchus. Three major conclusions were drawn: (i) Patterns of synonymous CAA interruptions in the microsatellite are characteristic and likely to result from selection for stabilizing the repetitive region; (ii) Interspecific comparisons indicate that SCA17 is likely to be a human trait. The most common allele in humans (37 repeats) is close to the threshold value upon which neurodegenerative changes can occur and may act as a repository for expanded, pathogenic alleles; (iii) The cassette-like structure of five out of 17 expanded alleles can be attributed to unequal crossing over. This can explain the rare and sporadic de novo generation of SCA17 alleles.

Microsatellite length differences between humans and chimpanzees at autosomal Loci are not found at equivalent haploid Y chromosomal Loci

When homologous microsatellites are compared between species, significant differences in mean length are often noted. A dominant cause of these length differences is ascertainment bias due to selection for maximum repeat number and repeat purity when the markers are being developed. However, even after ascertainment bias has been allowed for through reciprocal comparisons, significant length differences remain, suggesting that the average microsatellite mutation rate differs between species. Two classes of mechanism have been proposed: rapid evolution of enzymes involved in the generation and repair of slippage products (enzyme evolution model) and heterozygote instability, whereby interchromosomal events at heterozygous sites offer extra opportunities for mutations to occur (heterozygote instability model). To examine which of these hypotheses is most likely, we compared ascertainment bias and species length differences between humans and chimpanzees in autosomal and Y chromosomal microsatellites. We find that levels of ascertainment bias are indistinguishable, but that interspecies length differences are significantly greater for autosomal loci compared with haploid Y chromosomal loci. Such a pattern is consistent with predictions from the heterozygote instability model and is not expected under models of microsatellite evolution that do not include interchromosomal events such as the enzyme evolution model.

Neurological proteins are not enriched for repetitive sequences

Proteins associated with disease and development of the nervous system are thought to contain repetitive, simple sequences. However, genome-wide surveys for simple sequences within proteins have revealed that repetitive peptide sequences are the most frequent shared peptide segments among eukaryotic proteins, including those of Saccharomyces cerevisiae, which has few to no specialized developmental and neurological proteins. It is therefore of interest to determine if these specialized proteins have an excess of simple sequences when compared to other sets of compositionally similar proteins. We have determined the relative abundance of simple sequences within neurological proteins and find no excess of repetitive simple sequence within this class. In fact, polyglutamine repeats that are associated with many neurodegenerative diseases are no more abundant within neurological specialized proteins than within nonneurological collections of proteins. We also examined the codon composition of serine homopolymers to determine what forces may play a role in the evolution of extended homopolymers. Codon type homogeneity tends to be favored, suggesting replicative slippage instead of selection as the main force responsible for producing these homopolymers.

Mitochondrial genome size variation in New World and Old World populations of Drosophila melanogaster

Drosophila melanogaster originated in Africa, spread to Europe and Asia, and is believed to have colonized the New World in the past few hundred years. Levels of genetic variation are typically reduced in New World populations, consistent with a founder event following range expansion out of Africa and the Old World. We describe the patterns of mtDNA length variation within and among several populations of Drosophila melanogaster from the Old and New World. MtDNA length variation is due to insertion and deletion of tandem repeats in the control region (D-loop) of D. melanogaster mitochondrial genome. The distinct mutational dynamics of this system provide an opportunity to compare the patterns of variation in this marker to those of other markers with different mutational pressures and linkage relationships. The data show significantly more length variation in African and Asian samples than in New World samples. New World samples also show more pronounced skew of the length distribution. Our results are distinct from an earlier study that showed significantly higher levels of length variation and heteroplasmy. The level of heteroplasmy is highly correlated with the number of years that samples have been maintained in laboratory culture, suggesting that relaxed selection in small populations permits the accumulation of mtDNA length variation and heteroplasmy. Together, the data indicate that mtDNA length variants retain a signature of founder events and selection, and suggest that further investigation into the mutation-selection dynamics of the D-loop region of mtDNA would provide a distinct and informative marker for analysis of the recent history of populations.

Y-chromosomal microsatellite mutation rates: differences in mutation rate between and within loci

Precise estimates of mutation rates at Y-chromosomal microsatellite STR (short tandem repeat) loci make an important basis for paternity diagnostics and dating of Y chromosome lineage origins. There are indications of considerable locus mutation rate variability between (inter-) and within (intra-) loci. We have studied nine Y-STR loci-DYS19, DYS389I/II, DYS390, DYS391, DYS392, DYS393, DYS385, and DYS388-in 1,766 father-son pairs of confirmed paternity (a total of 15,894 meioses). Five biallelic markers were also analyzed in the fathers-Tat, YAP, 12f2, SRY1532, and 92R7-defining haplogroups 1, 2, 3, 4, 9, and 16, respectively. A total of 36 fragment length mutations were observed: 24 gains (22 single-step, two double-step) and 12 single-step losses. Thus, there was a significant surplus of gains (p=0.045). Overall, the mutation rate was positively correlated to STR repeat length and there was a significant relative excess of losses in long alleles and gains in short alleles (p=0.043). In contrast to the situation in autosomal STR loci and in MSY-1, no noteworthy correlation between mutation rate and the father's age at the child's birth was observed. We observed significant interlocus differences in Y-STR mutation rates (p<0.01). The number of observed mutations ranged from zero in DYS392 to eight in DYS391 and DYS390. We have also demonstrated obvious differences in mutation rates between the haplogroups studied (p=0.024), a phenomenon that is a reflection of the dependence of mutation rate on allele size. Our study has thus demonstrated the necessity of not only locus-specific, but even allele-specific, mutation rate estimates for forensic and population genetic purposes, and provides a considerable basis for such estimates.

An Analysis of Microsatellite Loci in Arabidopsis thaliana: Mutational Dynamics and Application

Microsatellite loci are among the most commonly used molecular markers. These loci typically exhibit variation for allele frequency distribution within a species. However, the factors contributing to this variation are not well understood. To expand on the current knowledge of microsatellite evolution, 20 microsatellite loci were examined for 126 accessions of the flowering plant, Arabidopsis thaliana. Substantial variability in mutation pattern among loci was found, most of which cannot be explained by the assumptions of the traditional stepwise mutation model or infinite alleles model. Here it is shown that the degree of locus diversity is strongly correlated with the number of contiguous repeats, more so than with the total number of repeats. These findings support a strong role for repeat disruptions in stabilizing microsatellite loci by reducing the substrate for polymerase slippage and recombination. Results of cluster analyses are also presented, demonstrating the potential of microsatellite loci for resolving relationships among accessions of A. thaliana.

The mutational spectrum of human autosomal tetranucleotide microsatellites

We studied by multiplex amplification and single-run electrophoretic analysis 10 microsatellite loci, composed of nine tetranucleotide-repeats (D1S1612, D3S2387, D4S2431, D5S2501, D10S1237, D15S657, D16S2622, D18S1270, and IFNAR-ALU) and one trinucleotide repeat (D2S1353). After elimination of proven null allele events involving D1S1612 and D5S2501 and of all data of D3S2387, in which we suspected but could not prove the occurrence of null alleles, we were left with nine loci, encompassing 24,224 meioses and 23 mutations. Twenty-two of the mutations (96%) were single-step events. Moreover, 18 of the mutations were paternal, four were maternal, and one was indeterminate. There was no significant difference between the number of additions and deletions in the mutants. Our findings are compatible with a simple model in which tetranucleotide microsatellites mutate primarily in paternal germinative cells by DNA slippage, such that the vast majority of mutations are equiprobable additions or deletions of a single-repeat unit. By combining the data from our tetranucleotide loci with literature information of highly and lowly mutable microsatellites, we observed a very highly significant correlation between mutation rate and the geometric mean of the length of the longest perfect repeat region (LRPR), compatible with a power or exponential relationship. The variation of the length of the LRPR explained as much as 80% of the variance of the mutation rate of autosomal tetranucleotide microsatellites.

Domain-level differences in microsatellite distribution and content result from different relative rates of insertion and deletion mutations

Microsatellites (short tandem polynucleotide repeats) are found throughout eukaryotic genomes at frequencies many orders of magnitude higher than the frequencies predicted to occur by chance. Most of these microsatellites appear to have evolved in a generally neutral manner. In contrast, microsatellites are generally absent from bacterial genomes except in locations where they provide adaptive functional variability, and these appear to have evolved under selection. We demonstrate a mutational bias towards deletion (repeat contraction) in a native chromosomal microsatellite of the bacterium Mycoplasma gallisepticum, through the collection and analysis of independent mutations in the absence of natural selection. Using this and similar existing data from two other bacterial species and four eukaryotic species, we find strong evidence that deletion biases resulting in repeat contraction are common in bacteria, while eukaryotic microsatellites generally experience unbiased mutation or a bias towards insertion (repeat expansion). This difference in mutational bias suggests that eukaryotic microsatellites should generally expand wherever selection does not exclude them, whereas bacterial microsatellites should be driven to extinction by mutational pressure wherever they are not maintained by selection. This is consistent with observed bacterial and eukaryotic microsatellite distributions. Hence, mutational biases that differ between eukaryotes and bacteria can account for many of the observed differences in microsatellite DNA content and distribution found in these two groups of organisms.

Mutation patterns at dinucleotide microsatellite loci in humans

Microsatellites are a major type of molecular markers in genetics studies. Their mutational dynamics are not clear. We investigated the patterns and characteristics of 97 mutation events unambiguously identified, from 53 multigenerational pedigrees with 630 subjects, at 362 autosomal dinucleotide microsatellite loci. A size-dependent mutation bias (in which long alleles are biased toward contraction, whereas short alleles are biased toward expansion) is observed. There is a statistically significant negative relationship between the magnitude (repeat numbers changed during mutation) and direction (contraction or expansion) of mutations and standardized allele size. Contrasting with earlier findings in humans, most mutation events (63%) in our study are multistep events that involve changes of more than one repeat unit. There was no correlation between mutation rate and recombination rate. Our data indicate that mutational dynamics at microsatellite loci are more complicated than the generalized stepwise mutation models.

Relative rates of insertion and deletion mutations in dinucleotide repeats of various lengths in mismatch repair proficient mouse and mismatch repair deficient human cells

Microsatellites are DNA elements composed of short tandem repeats of 1-5bp. These sequences are particularly prone to frameshift mutation by insertion-deletion loop formation during replication. The mismatch repair system is responsible for correcting these replication errors, and microsatellite mutation rates are significantly elevated in the absence of mismatch repair. We have investigated the effect of varying the number of repeats in a (CA)n microsatellite on mutation rates in cultured mammalian cells proficient or deficient in mismatch repair. We have also compared the relative rates of single-repeat insertions and deletions in these cells. Two plasmid vectors were constructed for each repeat unit number (n=8, 17, and 30), such that the microsatellites, placed upstream of a bacterial neomycin resistance gene (neo), disrupted the reading frame of the gene in the (-1) or (+1) direction. Plasmids were introduced separately into the cells, where they integrated into the cellular genome. Mutation rates were determined by selection of clones with frameshift mutations in the microsatellite that restored the reading frame of the neo gene. We found that mutation rates were significantly higher for (CA)17 and (CA)30 tracts than for (CA)8 tracts in both mismatch repair proficient (mouse) and deficient (human) cells. A mutational bias favoring insertions was generally observed. In both (CA)17 and (CA)30 tracts, single-repeat insertion rates were higher than single-repeat deletion rates with or without mismatch repair; deletions of multiple repeat units (> or =8bp) were observed in these tracts, where as deletions this large were not found in the (CA)8 tract. Single-repeat mutations of both types were made at similar rates in (CA)8 tracts in human mismatch repair deficient (MMR-) cells, but single-repeat insertion rates were higher than single-repeat deletion rates in mouse mismatch repair proficient (MMR+) cells. Results of these direct studies on microsatellite mutations in cultured cells should be useful for refinement of mathematical models for microsatellite evolution.

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

Recently Kruglyak, Durrett, Schug, and Aquadro showed that microsatellite equilibrium distributions can result from a balance between polymerase slippage and point mutations. Here, we introduce an elaboration of their model that keeps track of all parts of a perfect repeat and a simplification that ignores point mutations. We develop a detailed mathematical theory for these models that exhibits properties of microsatellite distributions, such as positive skewness of allele lengths, that are consistent with data but are inconsistent with the predictions of the stepwise mutation model. We use our theoretical results to analyze the successes and failures of the genetic distances (delta(mu))(2) and D(SW) when used to date four divergences: African vs. non-African human populations, humans vs. chimpanzees, Drosophila melanogaster vs. D. simulans, and sheep vs. cattle. The influence of point mutations explains some of the problems with the last two examples, as does the fact that these genetic distances have large stochastic variance. However, we find that these two features are not enough to explain the problems of dating the human-chimpanzee split. One possible explanation of this phenomenon is that long microsatellites have a mutational bias that favors contractions over expansions.

A Role for Selection in Regulating the Evolutionary Emergence of Disease-Causing and Other Coding CAG Repeats in Humans and Mice

The evolutionary expansion of CAG repeats in human triplet expansion disease genes is intriguing because of their deleterious phenotype. In the past, this expansion has been suggested to reflect a broad genomewide expansion of repeats, which would imply that mutational and evolutionary processes acting on repeats differ between species. Here, we tested this hypothesis by analyzing repeat- and flanking-sequence evolution in 28 repeat-containing genes that had been sequenced in humans and mice and by considering overall lengths and distributions of CAG repeats in the two species. We found no evidence that these repeats were longer in humans than in mice. We also found no evidence for preferential accumulation of CAG repeats in the human genome relative to mice from an analysis of the lengths of repeats identified in sequence databases. We then investigated whether sequence properties, such as base and amino acid composition and base substitution rates, showed any relationship to repeat evolution. We found that repeat-containing genes were enriched in certain amino acids, presumably as the result of selection, but that this did not reflect underlying biases in base composition. We also found that regions near repeats showed higher nonsynonymous substitution rates than the remainder of the gene and lower nonsynonymous rates in genes that contained a repeat in both the human and the mouse. Higher rates of nonsynonymous mutation in the neighborhood of repeats presumably reflect weaker purifying selection acting in these regions of the proteins, while the very low rate of nonsynonymous mutation in proteins containing a CAG repeat in both species presumably reflects a high level of purifying selection. Based on these observations, we propose that the mutational processes giving rise to polyglutamine repeats in human and murine proteins do not differ. Instead, we propose that the evolution of polyglutamine repeats in proteins results from an interplay between mutational processes and selection.