Mutation rate varies among alleles at a microsatellite locus: phylogenetic evidence

The major histocompatibility complex gene region and sarcoidosis susceptibility in African Americans

Investigators have intensively evaluated the major histocompatibility (MHC) complex for sarcoidosis susceptibility genes with the majority of reports implicating the human leukocyte antigen (HLA)-DRB1 gene. Because most studies have been performed in white and Asian populations, we sought to determine which MHC genes might be risk factors for sarcoidosis in African Americans. We genotyped six microsatellite markers spanning 11.6 megabases that overlapped the MHC region on chromosome 6p21-22 in 225 nuclear families ascertained by African American probands with a history of sarcoidosis. Using a family-based association methods approach, we performed multiallelic tests of association between each marker and sarcoidosis. A statistically significant association was detected between sarcoidosis and the DQCAR marker (p = 0.002) less than two kilobases telomeric from the HLA-DQB1 gene. Typing two additional markers in this region revealed that DQCAR-G51152 haplotypes, spanning a 38-kilobase region across the HLA-DQB1 gene, were associated with sarcoidosis on a global level (p = 0.022). Analysis of individual DQCAR and G51152 alleles showed that the DQCAR 178 (expected = 21.0; observed = 10; p = 0.0005) and G51152 217 (expected = 25.6; observed = 14; p = 0.0009) alleles were transmitted to affected offspring less often than expected; whereas the DQCAR 182 allele was transmitted more often than expected (expected = 52.6; observed = 66; p = 0.002). Our results indicate that HLA-DQB1 and not HLA-DRB1 plays an important role in sarcoidosis susceptibility in African Americans. Identification of the specific HLA-DQB1 alleles that influence sarcoidosis susceptibility in African Americans and the study of their antigenic-binding properties may reveal why African Americans suffer disproportionately from this disease.

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.

Hitchhiking mapping: a population-based fine-mapping strategy for adaptive mutations in Drosophilamelanogaster

The identification of genes contributing to the adaptation of local populations is of great biological interest. In an attempt to characterize functionally important differences among African and non-African Drosophila melanogaster populations, we surveyed neutral microsatellite variation in an 850-kb genomic sequence. Three genomic regions were identified that putatively bear an adaptive mutation associated with the habitat expansion of D. melanogaster. A further inspection of two regions by sequence analysis of multiple fragments confirmed the presence of a recent beneficial mutation in the non-African populations. Our study suggests that hitchhiking mapping is a universal approach for the identification of ecologically important mutations.

Enhanced selection for MHC diversity in social tuco-tucos

To explore the effects of behavior and demography on balancing selection at major histocompatibility complex (MHC) loci, we examined allelic diversity at exon 2 of the MHC class II DQbeta locus in a social and a solitary species of tuco-tuco (Rodentia: Ctenomyidae: Ctenomys), both of which occur in the same valley in southwestern Argentina. By comparing patterns of diversity at this MHC gene to the diversity evident at fifteen microsatellite loci, we demonstrate that balancing selection at the DQbeta locus is enhanced in the social species compared to its solitary congener. These findings have intriguing implications for the role of behavioral and demographic parameters in maintaining diversity at MHC loci.

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.

Detection of genetically divergent clone mates in apomictic dandelions

This study aims to identify genetically diverged clone mates in apomictic dandelions. Clone mates are defined as individuals that may have diverged as a result of mutation accumulation and that have undergone only clonal reproduction since their most recent common ancestor. Based on distinctive morphology and an aberrant and rare chloroplast haplotype, northwest European individuals of Taraxacum section Naevosa are well suited for the detection of clonal lineages in which mutation has occurred. In the case of strictly clonal reproduction, nuclear genetic variability was expected to be hierarchically organized. Nucleotide polymorphisms in internal transcribed spacer (ITS) sequences, however, were incompatible with a clonal structure of the Norwegian individuals, probably due to persistent ancestral polymorphisms that pre-date the origin of the Naevosa clone. This interpretation is supported by the presence of ITS variants in section Naevosa that were also found in distantly related dandelions. In contrast to the ITS sequence data, amplified fragment length polymorphisms (AFLPs), isozymes and microsatellites strongly supported the contention of prolonged clonal reproduction and mutation accumulation in Norwegian Naevosa. Because these markers are generally considered to be more variable and more rapidly evolving than ITS sequences, mutations in these markers probably evolved after the origin of the clone. Within the Norwegian clone, a surprising number of markers distinguished the clone mates. As a consequence, incorporation of mutation in the detection of clone mates is anticipated to have a big impact on estimates of size, geographical range and age of clones as well as on experimental designs of studies of clonal plants.

The estimation of population differentiation with microsatellite markers

Microsatellite markers are routinely used to investigate the genetic structuring of natural populations. The knowledge of how genetic variation is partitioned among populations may have important implications not only in evolutionary biology and ecology, but also in conservation biology. Hence, reliable estimates of population differentiation are crucial to understand the connectivity among populations and represent important tools to develop conservation strategies. The estimation of differentiation is c from Wright's FST and/or Slatkin's RST, an FST -analogue assuming a stepwise mutation model. Both these statistics have their drawbacks. Furthermore, there is no clear consensus over their relative accuracy. In this review, we first discuss the consequences of different temporal and spatial sampling strategies on differentiation estimation. Then, we move to statistical problems directly associated with the estimation of population structuring itself, with particular emphasis on the effects of high mutation rates and mutation patterns of microsatellite loci. Finally, we discuss the biological interpretation of population structuring estimates.

Using linked markers to infer the age of a mutation

Advances in sequencing and genotyping technologies over the last decade have enabled geneticists to easily characterize genetic variation at the nucleotide level. Hundreds of genes harboring mutations associated with genetic disease have now been identified by positional cloning. Using variation at closely linked genetic markers, it is possible to predict the times in the past at which particular mutations arose. Such studies suggest that many of the rare mutations underlying human genetic disorders are relatively young. Studies of variation at genetic markers linked to particular mutations can provide insights into human geographic history, and historical patterns of natural selection and disease, that are not available from other sources. We review two approaches for estimating allele age using variation at linked genetic markers. A phylogenetic approach aims to reconstruct the gene tree underlying a sample of chromosomes carrying a particular mutation, obtaining a "direct" estimate of allele age from the age of the root of this tree. A population genetic approach relies on models of demography, mutation, and/or recombination to estimate allele age without explicitly reconstructing the gene tree. Phylogenetic methods are best suited for studies of ancient mutations, while population genetic methods are better suited for studies of recent mutations. Methods that rely on recombination to infer the ages of alleles can be fine-tuned by choosing linked markers at optimal map distances to maximize the information available about allele age. A limitation of methods that rely on recombination is the frequent lack of a fine-scale linkage map. Maximum likelihood and Bayesian methods for estimating allele age that rely on intensive numerical computation are described, as well as "composite" likelihood and moment-based methods that lead to simple estimators. The former provide more accurate estimates (particularly for large samples of chromosomes) and should be employed if computationally practical.

Immune related genetic polymorphisms and schizophrenia among the Chinese

Genetic association studies were conducted among two independent cohorts of Chinese ethnicity. The samples consisted of cases and unrelated controls, ascertained from Guangzhou, China, and Singapore. The studies were prompted by our earlier report of an association between schizophrenia and HLA DQB1 alleles (HLA DQB1*0602 and HLA DQB1*0303) in the Singapore sample. Polymorphisms of HLA DQB1 and flanking markers on chromosome 6p21.3 were investigated in the first part of the study. A significant negative association with HLA DQB1*0402 was detected in the Guangzhou sample (Odds ratio, OR 0.26, 95% confidence intervals, CI 0.1, 0.6; p < 0.02, corrected for multiple comparisons). Additional analysis of the Guangzhou and Singapore samples revealed associations at three other anonymous markers flanking HLA DQB1. In the second part of the study, three polymorphisms at the Interleukin-1 gene cluster (IL-1, chromosome 2q13-q21) were investigated in both cohorts, since associations with schizophrenia have been reported in another sample. Persuasive evidence for an association at IL-1 was not detected in either sample. Our results suggest a susceptibility locus for schizophrenia in the HLA region among the Chinese, but further clarification is necessary.

Spectrum of nonrandom associations between microsatellite loci on human chromosome 11p15

Most evidence about nonrandom association of alleles at different loci, or gametic disequilibrium, across extensive anonymous regions of the human genome is based on the analysis of overall disequilibrium between pairs of microsatellites. However, analysis of interallelic associations is also necessary for a more complete description of disequilibrium. Here, we report a study characterizing the frequency and strength of both overall and interallelic disequilibrium between pairs of 12 microsatellite loci (CA repeats) spanning 19 cM (14 Mb) on human chromosome 11p15, in a large sample (810 haplotypes deduced from 405 individuals) drawn from a single population. Characterization of disequilibrium was carried out, taking into account the sign of the observed disequilibria. This strategy facilitates detection of associations and gives more accurate estimates of their intensities. Our results demonstrate that the incidence of disequilibrium over an extensive human chromosomal region is much greater than is commonly considered for populations that have expanded in size. In total, 44% of the pairs of microsatellite loci and 18% of the pairs of alleles showed significant nonrandom association. All the loci were involved in disequilibrium, although both the frequency and strength of interallelic disequilibrium were distributed nonuniformly along 11p15. These findings are especially relevant since significant associations were detected between loci separated by as much as 17-19 cM (7 cM on average). It was also found that the overall disequilibrium masks complicated patterns of association between pairs of alleles, dependent on their frequency and size. We suggest that the complex mutational dynamics at microsatellite loci could explain the allele-dependent disequilibrium patterns. These observations are also relevant to evaluation of the usefulness of microsatellite markers for fine-scale localization of disease genes.

Sustaining genetic variation in a small population: evidence from the Mauritius kestrel

We obtained measures of genetic diversity in 10 kestrel species at a suite of 12 microsatellite loci. We estimated the relative effective size (Ne) of the species using a Markov chain Monte Carlo (MCMC) approach, which jointly estimated the locus specific mutation rates as nuisance parameters. There was surprisingly high genetic diversity found in museum specimens of the Mauritius kestrel. Being an endemic species on a small island, it is known to have a long history of small population size. Conversely, kestrels with a continental distribution had Ne estimates that were only one order of magnitude larger and similar to each other, despite having current population sizes that were between one and three orders of magnitude larger than the Mauritius kestrel. We show how many of the theoretical results describing the effective size of a subdivided population can be captured in terms of three rates which describe the branching pattern of the gene genealogy, and that they are useful in estimating the time to migration-drift and mutation-drift equilibrium. We use this approach to argue that population subdivision has helped retain genetic diversity in the Mauritius kestrel, and that the continental species' genetic diversity has yet to reach equilibrium after the range changes following the last ice age. We draw parallels with Hewitt's observation that genetic variation seems to survive species' range compression and is rather vulnerable to range expansion.

Genetic variation at microsatellite loci in Spanish sheep

Genetic variation at 18 microsatellite loci was analysed in six indigenous Spanish sheep: Churra; Latxa; Manchega; Rasa-Aragonesa; Castellana and Merino. Merinos had frequently the highest number of alleles per locus, whereas Latxas showed the lowest one at many loci. Markers ordered decreasingly according to the number of variants differentiated in the whole population were: MAF70; TGLA13; CSSM66; BM143, BM6444; MAF36; MAF64; CSSM6; TGLA53; OarFCB11; MAF33; BM4621; MAF48; MAF65; BM1258; ILSTS002; ADCYC and OarCP34. Parameters of variability such as effective number of alleles and gene diversities corroborated the high level of variation frequently displayed by microsatellite markers. Comparison of allele distributions among populations and loci did not reveal consistent shapes. Distributions were centralised in some cases, whereas in others some kind of skewness was evident. Breed-specific alleles were detected at most loci, being frequent in Merinos and rare in Churras.

Microsatellite mutations in the germline: implications for evolutionary inference

Microsatellite DNA sequences mutate at rates several orders of magnitude higher than that of the bulk of DNA. Such high rates mean that spontaneous mutations that form new-length variants can realistically be seen in pedigree analysis. Data on observed mutation events from various organisms are now accumulating, allowing inferences on DNA sequence evolution to be made through an unusually direct approach. Here I discuss and integrate microsatellite mutation data in an evolutionary context. A striking feature of the mutation process is that it seems highly heterogeneous, with distinct differences between species, repeat types, loci and alleles. Age and sex also affect the mutation rate. Within genomes at equilibrium, the microsatellite-length distribution is a delicate balance between biased mutation processes and point mutations acting towards the decay of repetitive DNA. Indeed, simple repeats do not evolve simply.

On GATAGATA and other "junk" in the barren stretch of genomic desert

(GATA)n and additional simple repeated DNA sequences have long been known but functional analyses along traditional research strategies have not yielded definitive results. Therefore these elements are commonly regarded as genomic "junk" although they represent assets for innovative methodological approaches and especially for speculations. While extensively using microsatellites for genome mapping and genetic relationship analyses, it was also realized that certain simple repeats can cause human trinucleotide block expansion diseases. Hence these simple repeat elements comprise an exceptional range of biological meaning, i.e. from genetically absolutely neutral markers for genome analysis, behavioral and population genetic studies to dominant and invariably disease causing mutations with complete penetrance. The quest into why there is so much repetitive junk in recent genomes can be answered only in part and sequentially, but sometimes--even after substantial efforts--not for the given element.

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.

Microsatellites: consensus and controversy

Microsatellite DNA loci have recently been adopted for many biological applications. Comparative studies across a wide range of species has revealed many details of their mutational properties and evolutionary life cycles. Experience shows that a full understanding of these processes is essential to ensure the effective use of microsatellites as analytical tools. In this article, we review the controversies that have arisen as biologists have taken up this new technology and the emerging consensus that has resulted from their debates. We point to the need for comparative DNA sequencing studies to produce input data for a new generation of theoretical models of microsatellite behaviour. We conclude by presenting our own conceptual model, 'Snakes and Ladders', as an aid to theory development.

Removal of microsatellite interruptions by DNA replication slippage: phylogenetic evidence from Drosophila

Microsatellites are tandem repetitions of short (1-6 bp) motifs. It is widely assumed that microsatellites degenerate through the accumulation of base substitutions in the repeat array. Using a phylogenetic framework, we studied the evolutionary dynamics of interruptions in three Drosophila microsatellite loci. For all three loci, we show that the interruptions in a microsatellite can be lost, resulting in a longer uninterrupted microsatellite stretch. These results indicate that mutations in the microsatellite array do not necessarily lead to decay but may represent only a transition state during the evolution of a microsatellite. Most likely, this purification of interrupted microsatellites is caused by DNA replication slippage.

Characteristics and frequency of germline mutations at microsatellite loci from the human Y chromosome, as revealed by direct observation in father/son pairs

A number of applications of analysis of human Y-chromosome microsatellite loci to human evolution and forensic science require reliable estimates of the mutation rate and knowledge of the mutational mechanism. We therefore screened a total of 4,999 meioses from father/son pairs with confirmed paternity (probability >/=99. 9%) at 15 Y-chromosomal microsatellite loci and identified 14 mutations. The locus-specific mutation-rate estimates were 0-8. 58x10-3, and the average mutation rate estimates were 3.17x10-3 (95% confidence interval [CI] 1.89-4.94x10-3) across 8 tetranucleotide microsatellites and 2.80x10-3 (95% CI 1.72-4.27x10-3) across all 15 Y-chromosomal microsatellites studied. Our data show a mutational bias toward length increase, on the basis of observation of more repeat gains than losses (10:4). The data are in almost complete agreement with the stepwise-mutation model, with 13 single-repeat changes and 1 double-repeat change. Sequence analysis revealed that all mutations occurred in uninterrupted homogenous arrays of >/=11 repeats. We conclude that mutation rates and characteristics of human Y-chromosomal microsatellites are consistent with those of autosomal microsatellites. This indicates that the general mutational mechanism of microsatellites is independent of recombination.

Conservation and evolution of microsatellite loci in primate taxa

Microsatellites are promising genetic markers for the study of demographic structure and phylogenetic history in populations. However, little information exists on the molecular nature of the repeats and their flanking sequences of a same microsatellite in a large range of species. In this study, we report polymorphism and consensus sequences of eight microsatellite loci using human primers in 20 primate species. The results show size polymorphism in almost all species and microsatellites. These loci are therefore useful markers for population genetic studies between populations of the same species. Insertion/deletion events are frequent in the flanking regions, the majority concerning several contiguous bases. This is in contrast with the more usual single base pair events in non-coding regions. The ranges of allele lengths in non-human primates often show no overlap with that of human, usually due to the deletion/insertion events in the flanking sequences, producing smaller allele lengths rather than smaller numbers of repeats. The use of length of PCR product will bias the inter-species interpretation reducing the number of observable alleles and treating as the same allele very divergent molecular sequences. Caution should be used when employing microsatellites in cross-species comparisons in which the species under study are separated by significant amounts of evolutionary time: in such cases allele comparison cannot be based on lengths alone.

Trans-speciation maintenance in the MHC region of a polymorphism which includes a polymorphic dinucleotide locus, and the de novo arisal of a polymorphic tetranucleotide microsatellite

Alleles and the surrounding regions of DQCAR, a dinucleotide repeat tightly linked to HLA-DQB1, were sequenced in a range of primate species including man. Three polymorphic regions can usefully be defined in the description of these sequences: the dinucleotide GT repeat itself, the anonymous region 5' of this repeat, and a variable CTGT repeat in the 3' region. The 5' sequence displayed six alleles in the individuals studied. One of these alleles was invariably associated with substitutions in the GT repeat and absence of the CTGT repeat, the others with pure, polymorphic GT repeats and variation in the numbers of CTGT repeats. Haplotypes can be classified by the allele in the 5' region. Those carrying allele 1 were only found in man, those with allele 2 in man, chimpanzee and gorilla. The third haplotype (indicated by the presence of allele 3) was found in chimpanzee, gorilla and orang-utan, the fourth in chimpanzee and gibbon, the fifth in baboon, guenon and mangabey and the sixth in guenon and macaque. The alleles in the 5' region, but from different species, are thus often more similar than alleles from the same species, a phenomenon already shown for some HLA genes. This suggests that major histocompatibility sequences and surrounding sequences shared a correlated evolutionary history. The new polymorphic tetranucleotide microsatellite (CTGT, 3rd region) has possibly arisen de novo from the pre-existing dinucleotide GT. This study provides information not only on the molecular evolution of this particular microsatellite but also of the trans-speciation maintenance of polymorphism of its surrounding sequences.

Genetic diversity of two African and sixteen South American populations determined on the basis of six hypervariable loci

A total of 582 individuals (1,164 chromosomes) from two African, eight African-derived South American, five South American Amerindian, and three Brazilian urban populations were studied at four variable number of tandem repeat (VNTR) and two short tandem repeat (STR) hypervariable loci. These two sets of loci did not show distinct allele profiles, which might be expected if different processes promoted their molecular differentiation. The two African groups showed little difference between them, and their intrapopulational variation was similar to those obtained in the African-derived South American communities. The latter showed different degrees of interpopulation variability, despite the fact that they presented almost identical average degrees of non-African admixture. The F(ST) single locus estimates differed in the five sets of populations, probably due to genetic drift, indicating the need to consider population structure in the evaluation of their total variability. A high interpopulational diversity was found among Amerindian populations in relation to Brazilian African-derived isolated communities. This is probably a consequence of the differences in the patterns of gene flow and genetic drift that each of these semi-isolated groups experienced.

Clustered microsatellite mutations in the pipefish Syngnathus typhle

Clustered mutations are copies of a mutant allele that enter a population's gene pool together due to replication from a premeiotic germline mutation and distribution to multiple successful gametes of an individual. Although the phenomenon has been studied in Drosophila and noted in a few other species, the topic has received scant attention despite claims of being of major importance to population genetics theory. Here we capitalize upon the reproductive biology of male-pregnant pipefishes to document the occurrence of clustered microsatellite mutations and to estimate their rates and patterns from family data. Among a total of 3195 embryos genetically screened from 110 families, 40% of the 35 detected de novo mutant alleles resided in documented mutational clusters. Most of the microsatellite mutations appeared to involve small-integer changes in repeat copy number, and they arose in approximately equal frequency in paternal and maternal germlines. These findings extend observations on clustered mutations to another organismal group and motivate a broader critique of the mutation cluster phenomenon. They also carry implications for the evolution of microsatellites with respect to mutational models and homoplasy among alleles.

Tracing the origin of HLA-DRB1 alleles by microsatellite polymorphism

We analyzed the origin of allelic diversity at the class II HLA-DRB1 locus, using a complex microsatellite located in intron 2, close to the polymorphic second exon. A phylogenetic analysis of human, gorilla, and chimpanzee DRB1 sequences indicated that the structure of the microsatellite has evolved, primarily by point mutations, from a putative ancestral (GT)x(GA)y-complex-dinucleotide repeat. In all contemporary DRB1 allelic lineages, with the exception of the human *04 and the gorilla *08 lineages, the (GA)y repeat is interrupted, often by a G-->C substitution. In general, the length of the 3' (GA)y repeat correlates with the allelic lineage and thus evolves more slowly than a middle (GA)z repeat, whose length correlates with specific alleles within the lineage. Comparison of the microsatellite sequence from 30 human DRB1 alleles showed the longer 5' (GT)x to be more variable than the shorter middle (GA)z and 3' (GA)y repeats. Analysis of multiple samples with the same exon sequence, derived from different continents, showed that the 5' (GT)x repeat evolves more rapidly than the middle (GA)z and the 3' (GA)y repeats, which is consistent with findings of a higher mutation rate for longer tracts. The microsatellite-repeat-length variation was used to trace the origin of new DRB1 alleles, such as the new *08 alleles found in the Cayapa people of Ecuador and the Ticuna people of Brazil.

Evidence for complex mutations at microsatellite loci in Drosophila

Fifteen lines each of Drosophila melanogaster, D. simulans, and D. sechellia were scored for 19 microsatellite loci. One to four alleles of each locus in each species were sequenced, and microsatellite variability was compared with sequence structure. Only 7 loci had their size variation among species consistent with the occurrence of strictly stepwise mutations in the repeat array, the others showing extensive variability in the flanking region compared to that within the microsatellite itself. Polymorphisms apparently resulting from complex nonstepwise mutations involving the microsatellite were also observed, both within and between species. Maximum number of perfect repeats and variance of repeat count were found to be strongly correlated in microsatellites showing an apparently stepwise mutation pattern. These data indicate that many microsatellite mutation events are more complex than represented even by generalized stepwise mutation models. Care should therefore be taken in inferring population or phylogenetic relationships from microsatellite size data alone. The analysis also indicates, however, that evaluation of sequence structure may allow selection of microsatellites that more closely match the assumptions of stepwise models.

Microsatellite evolution: polarity of substitutions within repeats and neutrality of flanking sequences

Though extensively used in a variety of disciplines, the evolutionary pattern of microsatellite sequences is still unclear. We addressed several questions relating to microsatellite evolution by analysing historically accumulated mutation events in a large set of artiodactyl (CA)n repeats, through sequence analysis of orthologous bovine and ovine loci. The substitution rate in microsatellite flanking sequences was not different from that in intron sequences, suggesting that if intron sequences in general are selectively neutral, sequences close to microsatellites are similarly so. This observation thus does not support the idea that successful heterologous amplification of microsatellites across distantly related taxa would be due to flanking sequences generally being under some form of selection. Interestingly, the substitution rate at the first nucleotide positions flanking repeats was significantly higher than in sequences further away. Moreover, the substitution rate in repeat units in the very end of microsatellites was significantly higher than that in the middle of repeat regions. Together these observations suggest a relative instability close to the boundary between repetitive and unique sequences. We present three models that potentially could explain such a feature, all involving inefficiency of mismatch repair systems.

Mutability of microsatellites developed for the ant Camponotus consobrinus

Five highly polymorphic (GA)n microsatellite loci are reported for the formicine ant Camponotus consobrinus. The occurrence of many nests with a simple family structure enabled a search for new mutations, 11 of which were found from 3055 informative typings. These mutations were not randomly distributed across loci, 10 of them occurring at the locus Ccon70. The spectrum of mutations across alleles at Ccon70 was also nonrandom, with all of them occurring in alleles in the upper half of the allele size distribution. Six of the Ccon70 mutations decreased allele size. The mutations observed fit the stepwise mutation model well, i.e. mutations could always be assigned to an allele which differed in size from them by one repeat unit. The parental origins of the Ccon70 mutations were established and appear more female biased than vertebrate mutations, significantly so compared with human haemophilia A and primate intron mutations. This result may indicate that the lack of meiosis in males (which are haploid in ants) reduces the mutation rate in that sex relative to species in which both sexes are diploid.

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 single nucleotide polymorphisms in the HLA-DQ region

Sequencing studies were performed in three previously described microsatellite and minisatellite markers located within the HLA-DQ region, DQCAR, DQCARII and G51152. Multiple nucleotide substitutions that did not change size polymorphisms were observed in all three markers. In all loci, the number of core repeats did not correlate with neighboring DQ allele sequence motifs while single nucleotide changes within or flanking the microsatellite sequence did. This result indicates higher mutation rates for microsatellite expansions/contractions than for nucleotide substitutions in these loci. Further analysis indicated an almost complete phylogenetic correspondence between DQCAR single nucleotide polymorphisms (SNPs) and DQB1 sequences on one side (1.0-1.5 kb apart) and a complete relationship between DQCARII and DQA1 sequences on the other (4.5 kb apart). In contrast, G51152 sequences did not correspond perfectly with DQB1 allelic sequences, thus suggesting the existence of several ancestral crossovers between this marker and DQB1 (20-25 kb). Sequencing microsatellites might be useful in disease mapping studies by increasing marker informativeness and by helping in the interpretation of association study results. It is also proposed that SNPs within the flanking region of CA repeats could be used to develop biallelic markers from already available mapped microsatellite markers.

Mutation rate in human microsatellites: influence of the structure and length of the tandem repeat

In 10,844 parent/child allelic transfers at nine short-tandem-repeat (STR) loci, 23 isolated STR mismatches were observed. The parenthood in each of these cases was highly validated (probability >99.97%). The event was always repeat related, owing to either a single-step mutation (n=22) or a double-step mutation (n=1). The mutation rate was between 0 and 7 x 10(-3) per locus per gamete per generation. No mutations were observed in three of the nine loci. Mutation events in the male germ line were five to six times more frequent than in the female germ line. A positive exponential correlation between the geometric mean of the number of uninterrupted repeats and the mutation rate was observed. Our data demonstrate that mutation rates of different loci can differ by several orders of magnitude and that different alleles at one locus exhibit different mutation rates.

Signatures of population expansion in microsatellite repeat data

To examine the signature of population expansion on genetic variability at microsatellite loci, we consider a population that evolves according to the time-continuous Moran model, with growing population size and mutations that follow a general asymmetric stepwise mutation model. We present calculations of expected allele-size variance and homozygosity at a locus in such a model for several variants of growth, including stepwise, exponential, and logistic growth. These calculations in particular prove that population bottleneck followed by growth in size causes an imbalance between allele size variance and heterozygosity, characterized by the variance being transiently higher than expected under equilibrium conditions. This effect is, in a sense, analogous to that demonstrated before for the infinite allele model, where the number of alleles transiently increases after a stepwise growth of population. We analyze a set of data on tetranucleotide repeats that reveals the imbalance expected under the assumption of bottleneck followed by population growth in two out of three major racial groups. The imbalance is strongest in Asians, intermediate in Europeans, and absent in Africans. This finding is consistent with previous findings by others concerning the population expansion of modern humans, with the bottleneck event being most ancient in Africans, most recent in Asians, and intermediate in Europeans. Nevertheless, the imbalance index alone cannot reliably estimate the time of initiation of population expansion.

DQ microsatellite association studies in three ethnic groups

Polymorphism at the level of three microsatellite markers (DQCAR, DQCARII, G51152) located in the HLA-DQ region was characterized in 78 10th International Histocompatibility Workshop B-cell lines, 718 random Japanese Asians, 99 Norwegian Caucasians and 95 New Guinean Aborigines with established HLA-DRB1, -DQA1 and -DQB1 typing. DQCAR, DQCARII, and G51152 result in 13, 13, and 11 alleles respectively. All three markers were in tight linkage disequilibrium with HLA-DRB1, -DQA1 and -DQB1. DRB1, DQA1, DQCARII, DQCAR, DQB1, and G51152 haplotypes could be defined for all subjects. In fact, DQ microsatellite typing data could predict DQA1 and DQB1 genotypes with high accuracy and may be used as a simple first pass HLA-DQ typing method. The haplotype data was also used to determine recombination in the DRB1-DQA1 (about 80 kb), DQA1-DQCARII (about 4.5 kb), DQCARII-DQCAR (about 7.5 kb), DQCAR-DQB1 (about 1-1.5 kb) and DQB1-G51152 (about 20-25 kb) genomic segments and the relative rate of slippage microsatellite mutations for DQCAR, DQCARII, and G51152. This led us to conclude that recombination is more frequent in the DRB1-DQA1 and DQCAR-DQCARII segments, thus suggesting cross-overs within small genomic segments are not proportional to genetic distance. We also observed that DQCAR had a higher mutation rate than DQCARII or G51152 and that 1 or 2 CA slippage mutations were arising more frequently from large size microsatellite alleles.

Phylogenetic assessment of length variation at a microsatellite locus

Sixty-six haplotypes at a locus containing a simple dinucleotide (CA)n microsatellite repeat were isolated by PCR-single-strand conformational polymorphism from populations of the horseshoe crab Limulus polyphemus. These haplotypes were sequenced to assess nucleotide variation directly. Thirty-four distinct sequences (alleles) were identified in a region 570 bp long that included the microsatellite motif. In the repeat region itself, CA-number varied in integer values from 5 to 11 across alleles, except that a (CA)8 class was not observed. Differences among alleles were due also to polymorphisms at 22 sites in regions immediately flanking the microsatellite repeats. Nucleotide substitutions in these regions were used to estimate phylogenetic relationships among alleles, and the gene phylogeny was used to trace the evolution of length variation and CA repeat numbers. A low correlation between size variation and genealogical relationships among alleles suggests that absolute fragment size (as normally scored in microsatellite assays) is an unreliable indicator of historical affinities among alleles. This finding on the molecular fine structure of microsatellite variation suggests the need for caution in the use of repeat counts at microsatellite loci as secure indicators of allelic relationships.

The complex mutation pattern of a microsatellite

DQCAR is a (CA)n microsatellite located in the HLA class II region and tightly linked to HLA-DQB1. Previous studies showed a strikingly low level of size variation in DQCAR alleles within an extensive subfamily of HLA-DQ subtypes (DQ1). DQCAR alleles in non-DQ1 subtypes showed a higher degree of size polymorphism. In this study sequence analysis demonstrates that DQ1-associated DQCAR alleles have a single C-->A nucleotide substitution interrupting the CA repeat array. Frequent CA-->GA mutations are also observed in DQ1-associated microsatellites with identical allele sizes. In contrast, DQCAR alleles associated with non-DQ1 haplotypes display a perfect CA repeat sequence and the variation in allele size is attributable only to differences in the number of CA repeats. Our results imply that several mutational mechanisms are involved in the generation of allelic diversity within the same microsatellite locus. The possibility of different mutation rates in the same locus should to be taken into account when using these markers in evolutionary and disease studies.

Microsatellite diversity and the demographic history of modern humans

We have examined differences in diversity at 60 microsatellite loci among human population samples from three major continental groups to evaluate the hypothesis of greater African diversity in this rapidly evolving class of loci. Application of a statistical test that assumes equal mutation rates at all loci fails to demonstrate differences in microsatellite diversity, while a randomization test that does not make this assumption finds that Africans have significantly greater microsatellite diversity (P < 10(-8)) than do Asians and Europeans. Greater African diversity is most apparent at loci with smaller overall variance in allele size, suggesting that the record of population history has been erased at repeat loci with higher mutation rates. A power analysis shows that only 35-40 microsatellites are needed to establish this difference statistically, demonstrating the considerable evolutionary information contained in these systems. On average, African populations have approximately 20% greater microsatellite diversity than do Asian and European populations. A comparison of continental diversity differences in microsatellites and mtDNA sequences suggests earlier demographic expansion of the ancestors of Africans.