Rate and directionality of mutations and effects of allele size constraints at anonymous, gene-associated, and disease-causing trinucleotide loci.

We studied the patterns of within- and between-population variation at 29 trinucleotide loci in a random sample of 200 healthy individuals from four diverse populations: Germans, Nigerians, Chinese, and New Guinea highlanders. The loci were grouped as disease-causing (seven loci with CAG repeats), gene-associated (seven loci with CAG/CCG repeats and eight loci with AAT repeats), or anonymous (seven loci with AAT repeats). We used heterozygosity and variance of allele size (expressed in units of repeat counts) as measures of within-population variability and GST (based on heterozygosity as well as on allele size variance) as the measure of genetic differentiation between populations. Our observations are: (1) locus type is the major significant factor for differences in within-population genetic variability; (2) the disease-causing CAG repeats (in the nondisease range of repeat counts) have the highest within-population variation, followed by the AAT-repeat anonymous loci, the AAT-repeat gene-associated loci, and the CAG/CTG-repeat gene-associated loci; (3) an imbalance index beta, the ratio of the estimates of the product of effective population size and mutation rate based on allele size variance and heterozygosity, is the largest for disease-causing loci, followed by AAT- and CAG/CCG-repeat gene-associated loci and AAT-repeat anonymous loci; (4) mean allele size correlates positively with allele size variance for AAT- and CAG/CCG-repeat gene-associated loci and negatively for anonymous loci; and (5) GST is highest for the disease-causing loci. These observations are explained by specific differences of rates and patterns of mutations in these four groups of trinucleotide loci, taking into consideration the effects of the past demographic history of the modern human population.     

Genetic variation at nine short tandem repeat loci among islanders of the eastern Adriatic coast of Croatia

We have analyzed the extent of genetic variation at nine autosomal short tandem repeat loci (D3S1358, VWA, FGA, TH01, TPOX, CSF1PO, D5S818, D13S317, D7S820) among six populations from Croatia: five distributed in the islands of the eastern Adriatic coast and one from the mainland. The purpose is to investigate the usefulness of these loci in detecting regional genetic differentiation in the studied populations. Significant heterogeneity among the island and mainland populations is revealed in the distributions of allele frequencies; however, the absolute magnitude of the coefficient of gene differentiation is small but significant. The summary measures of genetic variation, namely, heterozygosity, number of alleles, and allele size variance, do not indicate reduced genetic variation in the island populations compared to the mainland population. In contrast to the two measures of genetic variation, allele size variance and within-locus heterozygosity, the imbalance index (beta) indicates evidence of recent expansion of population sizes in all islands and in the mainland. High mutation rates of the studied loci together with local drift effects are likely explanations for interisland genetic variation and the observed lack of reduced genetic diversity among the island populations.     

Allele Frequencies at Microsatellite Loci: The Stepwise Mutation Model Revisited

We summarize available data on the frequencies of alleles at microsatellite loci in human populations and compare observed distributions of allele frequencies to those generated by a simulation of the stepwise mutation model. We show that observed frequency distributions at 108 loci are consistent with the results of the model under the assumption that mutations cause an increase or decrease in repeat number by one and under the condition that the product Nu, where N is the effective population size and u is the mutation rate, is larger than one. We show that the variance of the distribution of allele sizes is a useful estimator of Nu and performs much better than previously suggested estimators for the stepwise mutation model. In the data, there is no correlation between the mean and variance in allele size at a locus or between the number of alleles and mean allele size, which suggests that the mutation rate at these loci is independent of allele size.     

Dependence of expected heterozygosity on locus number with stabilizing selection and drift

I determine expected levels of heterozygosity in two allele multilocus models with mutation, stabilizing selection and drift. In the range 2 to 32 loci, the per locus heterozygosity can depend on the locus number. The per locus heterozygosity for ten loci can be as low as three fourths of the per locus heterozygosity in the limit, as the number of loci gets large. Simulations indicate that this dependence on locus number is not due to the population approaching equilibria at which the mean differs from the optimum, but is due to changes in the substitution rate as a function of the number of loci.     

Characteristics of loci and individuals are associated with germline microsatellite mutation rates in lesser kestrels (Falco naumanni)

Although microsatellites are one of the most popular tools in genetic studies, their mutational dynamics and evolution remain unclear. Here, we apply extensive pedigree genotyping to identify and analyze the patterns and factors associated with de novo germline mutations across nine microsatellite loci in a wild population of lesser kestrels (Falco naumanni). A total of 10 germline mutations events were unambiguously identified in four loci, yielding an average mutation rate of 2.96x10(-3). Across loci, mutation rate was positively correlated with locus variability and average allele size. Mutations were primarily compatible with a stepwise mutation model, although not exclusively involved single-step changes. Unexpectedly, we found an excess of maternally transmitted mutations (male-to-female ratio of 0.1). One of the analyzed loci (Fn2.14) resulted hypermutable (mutation rate=0.87%). This locus showed a size-dependent mutation bias, with longer alleles displaying deletions or additions of a small number of repeat than shorter alleles. Mutation probability at Fn2.14 was higher for females and increased with parental (maternal) age but was not associated with individual physical condition, multilocus heterozygosity, allele length or allele span. Overall, our results do not support the male-biased mutation rate described in other organisms and suggest that mutation dynamics at microsatellite loci are a complex process which requires further research.     

Hitchhiking and associative overdominance at a microsatellite locus

The possible effects of a selected locus on a closely linked microsatellite locus are discussed and analyzed in terms of coalescent theory and models of the mutation process. Background selection caused by recurrent deleterious mutations will reduce the variance of allele size at a microsatellite locus. The occasional substitution of advantageous alleles (genetic hitchhiking) will also reduce the variance, but a high mutation rate at a microsatellite locus can restore the variance relatively rapidly. Overdominance at the selected locus will increase the variance at the microsatellite locus and create partitioning of the variation in allele size among gametes carrying one or the other of the overdominant alleles. These results suggest that neutral microsatellite loci can provide indicators of selective processes at closely linked loci.      

Adult survival and microsatellite diversity in possums: effects of major histocompatibility complex-linked microsatellite diversity but not multilocus inbreeding estimators

Adult survival is perhaps the fitness parameter most important to population growth in long-lived species. Intrinsic and extrinsic covariates of survival are therefore likely to be important drivers of population dynamics. We used long-term mark-recapture data to identify genetic, individual and environmental covariates of local survival in a natural population of mountain brushtail possums (Trichosurus cunninghami). Rainfall and intra-individual diversity at microsatellite DNA markers were associated with increased local survival of adults and juveniles. We contrasted the performance of several microsatellite heterozygosity measures, including internal relatedness (IR), homozygosity by loci (HL) and the mean multilocus estimate of the squared difference in microsatellite allele sizes within an individual (mean d ²). However, the strongest effect on survival was not associated with multilocus microsatellite diversity (which would indicate a genome-wide inbreeding effect), but a subset of two loci. This included a major histocompatibility complex (MHC)-linked marker and a putatively neutral microsatellite locus. For both loci, diversity measures incorporating allele size information had stronger associations with survival than measures based on heterozygosity, whether or not allele frequency information was included (such as IR). Increased survival was apparent among heterozygotes at the MHC-linked locus, but the benefits of heterozygosity to survival were reduced in heterozygotes with larger differences in allele size. The effect of heterozygosity on fitness-related traits was supported by data on endoparasites in a subset of the individuals studied in this population. There was no apparent density dependence in survival, nor an effect of sex, age or immigrant status. Our findings suggest that in the apparent absence of inbreeding, variation at specific loci can generate strong associations between fitness and diversity at linked markers.     

Usefulness of molecular markers for detecting population bottlenecks via monitoring genetic change

It is important to detect population bottlenecks in threatened and managed species because bottlenecks can increase the risk of population extinction. Early detection is critical and can be facilitated by statistically powerful monitoring programs for detecting bottleneck-induced genetic change. We used Monte Carlo computer simulations to evaluate the power of the following tests for detecting genetic changes caused by a severe reduction in a population's effective size ( N _e): a test for loss of heterozygosity, two tests for loss of alleles, two tests for change in the distribution of allele frequencies, and a test for small N _e based on variance in allele frequencies (the 'variance test'). The variance test was most powerful; it provided an 85% probability of detecting a bottleneck of size N _e = 10 when monitoring five microsatellite loci and sampling 30 individuals both before and one generation after the bottleneck. The variance test was almost 10-times more powerful than a commonly used test for loss of heterozygosity, and it allowed for detection of bottlenecks before 5% of a population's heterozygosity had been lost. The second most powerful tests were generally the tests for loss of alleles. However, these tests had reduced power for detecting genetic bottlenecks caused by skewed sex ratios. We provide guidelines for the number of loci and individuals needed to achieve high-power tests when monitoring via the variance test. We also illustrate how the variance test performs when monitoring loci that have widely different allele frequency distributions as observed in five wild populations of mountain sheep ( Ovis canadensis ).     

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

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

Measures of Variation at DNA Repeat Loci under a General Stepwise Mutation Model

Polymorphisms at tandem repeat loci are caused by mutations with allele sizes occasionally altered by more than one repeat unit in both forward and backward directions. Such mutational changes may occur with asymmetric probabilities. Therefore, a one-step symmetric stepwise mutation model may not be appropriate for studying the population dynamics at all repeat loci. In this work, we evaluated the expectation and variance of the within-population variance of the allele size distribution in a finite population, and the expected homozygosity at a locus by the coalescence approach under a general stepwise mutation model, where mutational transitions of allele sizes can be arbitrary, including being asymmetric. Under the special cases of symmetric one-step, two-step, and multi-step geometric distributions of mutations, our general results reduce to the corresponding results obtained by earlier investigators. The general results indicate that in a finite population, which has reached a steady state under the (general stepwise) mutation and drift balance, the within-population variance of allele sizes has a simple expectation (i.e., proportional toNν, the product of the mutation rate,ν, and effective population size,N). However, its stochastic variance is a quadratic function of this composite parameter,Nν. Furthermore, this second-order variance does not decay with the number of alleles sampled from a population. Application of this theory to data on allele size distributions in unrelated Caucasians from the CEPH pedigree (obtained from the Genome Data Base) shows that the relationship of the variance and mean of within-population variance of allele sizes at tandem repeat loci, grouped by their chromosomal assignment, has a trend compatible with the theory. However, there is an indication that the second-order variance is generally underestimated. One reason for this departure might be that the CEPH sample may not represent a single homogeneous population that reached equilibrium at all tandem repeat loci.     

Evidence for heterozygote instability in microsatellite loci in house wrens

Microsatellite loci have high mutation rates and high levels of allelic variation, but the factors influencing their mutation rate are not well understood. The proposal that heterozygosity may increase mutation rates has profound implications for understanding the evolution of microsatellite loci, but currently has limited empirical support. We examined 20 microsatellite mutations identified in an analysis of 12 260 meiotic events across three loci in two populations of a songbird, the house wren (Troglodytes aedon). We found that for an allele of a given length, mutation was significantly more likely when there was a relatively large difference in size between the allele and its homologue (i.e. a large ‘allele span’). Our results support the proposal of heterozygote instability at microsatellite loci.     

Genetic analysis of microsatellite polymorphism in the Elliot's Pheasant (Syrmaticus ellioti) in China

In this study, we reported the population genetic analyses in the Elliot's Pheasant(Syrnaticus ellioti) using seven polymorphism microsatellite loci based on 105 individuals from 4 geographical populations. Departures from Hardy-Weinberg equilibrium were found in four geographical populations. The average number of alleles was 8.86, with a total of 62 alleles across 7 loci; observed heterozygosity (HO) was generally low and the average number was 0.504. For the seven microsatellite loci, the polymorphism information content ranged from 0.549 to 0.860, with an average number 0.712. Population bottlenecks of the four geographical populations were tested by infinite allele mutation model, step-wise mutation model and two-phase mutation model, which found that each population had experienced bottleneck effect during the recent period. Fst analysis across all geographical populations indicated that the genetic differentiaton between the Guizhou geographical population and the Hunan geographical population was highly significant (P<0.001), a finding supported by the far genetic relationship showed by the neighbor-joining tree of four geographical populations based on Nei's unbiased genetic distances. Using hierarchical analysis of molecular variance (Guizhou geographical population relative to all others pooled), we found a low level of the genetic variation among geographical populations and that between groups. However, differences among populations relative to the total sample explained most of the genetic variance (92.84%), which was significant.     

Towards Improvements in the Estimation of the Coalescent: Implications for the Most Effective Use of Y Chromosome Short Tandem Repeat Mutation Rates

Over the past two decades, many short tandem repeat (STR) microsatellite loci on the human Y chromosome have been identified together with mutation rate estimates for the individual loci. These have been used to estimate the coalescent age, or the time to the most recent common ancestor (TMRCA) expressed in generations, in conjunction with the average square difference measure (ASD), an unbiased point estimator of TMRCA based upon the average within-locus allele variance between haplotypes. The ASD estimator, in turn, depends on accurate mutation rate estimates to be able to produce good approximations of the coalescent age of a sample. Here, a comparison is made between three published sets of per locus mutation rate estimates as they are applied to the calculation of the coalescent age for real and simulated population samples. A novel evaluation method is developed for estimating the degree of conformity of any Y chromosome STR locus of interest to the strict stepwise mutation model and specific recommendations are made regarding the suitability of thirty-two commonly used Y-STR loci for the purpose of estimating the coalescent. The use of the geometric mean for averaging ASD and across loci is shown to improve the consistency of the resulting estimates, with decreased sensitivity to outliers and to the number of STR loci compared or the particular set of mutation rates selected.     

Microsatellite mutations and inferences about human demography

Microsatellites have been widely used as tools for population studies. However, inference about population processes relies on the specification of mutation parameters that are largely unknown and likely to differ across loci. Here, we use data on somatic mutations to investigate the mutation process at 14 tetranucleotide repeats and carry out an advanced multilocus analysis of different demographic scenarios on worldwide population samples. We use a method based on less restrictive assumptions about the mutation process, which is more powerful to detect departures from the null hypothesis of constant population size than other methods previously applied to similar data sets. We detect a signal of population expansion in all samples examined, except for one African sample. As part of this analysis, we identify an "anomalous" locus whose extreme pattern of variation cannot be explained by variability in mutation size. Exaggerated mutation rate is proposed as a possible cause for its unusual variation pattern. We evaluate the effect of using it to infer population histories and show that inferences about demographic histories are markedly affected by its inclusion. In fact, exclusion of the anomalous locus reduces interlocus variability of statistics summarizing population variation and strengthens the evidence in favor of demographic growth.     

Estimation of effective population size and migration rate from one- and two-locus identity measures

Standard methods for inferring demographic parameters from genetic data are based mainly on one-locus theory. However, the association of genes at different loci (e.g., two-locus identity disequilibrium) may also contain some information about demographic parameters of populations. In this article, we define one- and two-locus parameters of population structure as functions of one- and two-locus probabilities for the identity in state of genes. Since these parameters are known functions of demographic parameters in an infinite island model, we develop moment-based estimators of effective population size and immigration rate from one- and two-locus parameters. We evaluate this method through simulation. Although variance and bias may be quite large, increasing the number of loci on which the estimates are derived improves the method. We simulate an infinite allele model and a K allele model of mutation. Bias and variance are smaller with increasing numbers of alleles per locus. This is, to our knowledge, the first attempt of a joint estimation of local effective population size and immigration rate.     

Vntr Allele Frequency Distributions under the Stepwise Mutation Model: A Computer Simulation Approach

Variable numbers of tandem repeats (VNTRs) are a class of highly informative and widely dispersed genetic markers. Despite their wide application in biological science, little is known about their mutational mechanisms or population dynamics. The objective of this work was to investigate four summary measures of VNTR allele frequency distributions: number of alleles, number of modes, range in allele size and heterozygosity, using computer simulations of the one-step stepwise mutation model (SMM). We estimated these measures and their probability distributions for a wide range of mutation rates and compared the simulation results with predictions from analytical formulations of the one-step SMM. The average heterozygosity from the simulations agreed with the analytical expectation under the SMM. The average number of alleles, however, was larger in the simulations than the analytical expectation of the SMM. We then compared our simulation expectations with actual data reported in the literature. We used the sample size and observed heterozygosity to determine the expected value, 5th and 95th percentiles for the other three summary measures, allelic size range, number of modes and number of alleles. The loci analyzed were classified into three groups based on the size of the repeat unit: microsatellites (1-2 base pair (bp) repeat unit), short tandem repeats [(STR) 3-5 bp repeat unit], and minisatellites (15-70 bp repeat unit). In general, STR loci were most similar to the simulation results under the SMM for the three summary measures (number of alleles, number of modes and range in allele size), followed by the microsatellite loci and then by the minisatellite loci, which showed deviations in the direction of the infinite allele model (IAM). Based on these differences, we hypothesize that these three classes of loci are subject to different mutational forces.     

MICROSATELLITE EVOLUTION IN VERTEBRATES: INFERENCE FROM AC DINUCLEOTIDE REPEATS

Abstract We analyze published data from 592 AC microsatellite loci from 98 species in five vertebrate classes including fish, reptiles, amphibians, birds, and mammals. We use these data to address nine major questions about microsatellite evolution. First, we find that larger genomes do not have more microsatellite loci and therefore reject the hypothesis that microsatellites function primarily to package DNA into chromosomes. Second, we confirm that microsatellite loci are relatively rare in avian genomes, but reject the hypothesis that this is due to physical constraints imposed by flight. Third, we find that microsatellite variation differs among species within classes, possibly relating to population dynamics. Fourth, we reject the hypothesis that microsatellite structure (length, number of alleles, allele dispersion, range in allele sizes) differs between poikilotherms and homeotherms. The difference is found only in fish, which have longer microsatellites and more alleles than the other classes. Fifth, we find that the range in microsatellite allele size at a locus is largely due to the number of alleles and secondarily to allele dispersion. Sixth, length is a major factor influencing mutation rate. Seventh, there is a directional mutation toward an increase in microsatellite length. Eighth, at the species level, microsatellite and allozyme heterozygosity covary and therefore inferences based on large-scale studies of allozyme variation may also reflect microsatellite genetic diversity. Finally, published microsatellite loci (isolated using conventional hybridization methods) provide a biased estimate of the actual mean repeat length of microsatellites in the genome.     

Identification of selective sweeps in closely related populations of the house mouse based on microsatellite scans

Genome scans of polymorphisms promise to provide insights into the patterns and frequencies of positive selection under natural conditions. The use of microsatellites as markers has the potential to focus on very recent events, since in contrast to SNPs, their high mutation rates should remove signatures of older events. We assess this concept here in a large-scale study. We have analyzed two population pairs of the house mouse, one pair of the subspecies Mus musculus domesticus and the other of M. m. musculus. A total of 915 microsatellite loci chosen to cover the whole genome were assessed in a prescreening procedure, followed by individual typing of candidate loci. Schlötterer's ratio statistics (lnRH) were applied to detect loci with significant deviations from patterns of neutral expectation. For eight loci from each population pair we have determined the size of the potential sweep window and applied a second statistical procedure (linked locus statistics). For the two population pairs, we find five and four significant sweep loci, respectively, with an average estimated window size of 120 kb. On the basis of the analysis of individual allele frequencies, it is possible to identify the most recent sweep, for which we estimate an onset of 400–600 years ago. Given the known population history for the French–German population pair, we infer that the average frequency of selective sweeps in these populations is higher than 1 in 100 generations across the whole genome. We discuss the implications for adaptation processes in natural populations.     

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

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