Comparisons of mutation rate variation at genome-wide microsatellites: evolutionary insights from two cultivated rice and their wild relatives

Background  

Mutation rate (μ) per generation per locus is an important parameter in the models of population genetics. Studies on mutation rate and its variation are of significance to elucidate the extent and distribution of genetic variation, further infer evolutionary relationships among closely related species, and deeply understand genetic variation of genomes. However, patterns of rate variation of microsatellite loci are still poorly understood in plant species. Furthermore, how their mutation rates vary in di-, tri-, and tetra-nucleotide repeats within the species is largely uninvestigated across related plant genomes.     

Rate and pattern of mutation at microsatellite loci in maize

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

High mutation rate of a long microsatellite allele in Drosophila melanogaster provides evidence for allele-specific mutation rates

Within recent years, microsatellite have become one of the most powerful genetic markers in biology. For several mammalian species, microsatellite mutation rates have been estimated on the order of 10(- 3)-10(-5). A recent study, however, demonstrated mutation rates in Drosophila melanogaster of at least one order of magnitude lower than those in mammals. To further test this result, we examined mutation rates of different microsatellite loci using a larger sample size. We screened 24 microsatellite loci in 119 D. melanogaster lines maintained for approximately 250 generations and detected 9 microsatellite mutations. The average mutation rate of 6.3 x 10(-6) is identical to the mutation rate from a previous study. Most interestingly, all nine mutations occurred at the same allele of one locus (DROYANETSB). This hypermutable allele has 28 dinucleotide repeats and is among the longest microsatellite reported in D. melanogaster. The allele-specific mutation rate of 3.0 x 10(-4) per generation is within the range of mammalian mutation rates. Future microsatellite analyses will have to account for the dramatic differences in allele-specific mutation rates.      

The mutation rates of di-, tri- and tetranucleotide repeats in Drosophila melanogaster

In a recent study, we reported that the combined average mutation rate of 10 di-, 6 tri-, and 8 tetranucleotide repeats in Drosophila melanogaster was 6.3 x 10(-6) mutations per locus per generation, a rate substantially below that of microsatellite repeat units in mammals studied to date (range = 10(-2)-10(-5) per locus per generation). To obtain a more precise estimate of mutation rate for dinucleotide repeat motifs alone, we assayed 39 new dinucleotide repeat microsatellite loci in the mutation accumulation lines from our earlier study. Our estimate of mutation rate for a total of 49 dinucleotide repeats is 9.3 x 10(-6) per locus per generation, only slightly higher than the estimate from our earlier study. We also estimated the relative difference in microsatellite mutation rate among di-, tri-, and tetranucleotide repeats in the genome of D. melanogaster using a method based on population variation, and we found that tri- and tetranucleotide repeats mutate at rates 6.4 and 8.4 times slower than that of dinucleotide repeats, respectively. The slower mutation rates of tri- and tetranucleotide repeats appear to be associated with a relatively short repeat unit length of these repeat motifs in the genome of D. melanogaster. A positive correlation between repeat unit length and allelic variation suggests that mutation rate increases as the repeat unit lengths of microsatellites increase.      

Human Microsatellites: Mutation and Evolution

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

Simultaneous estimation of all the parameters of a stepwise mutation model.

Minisatellite and microsatellite are short tandemly repetitive sequences dispersed in eukaryotic genomes, many of which are highly polymorphic due to copy number variation of the repeats. Because mutation changes copy numbers of the repeat sequences in a generalized stepwise fashion, stepwise mutation models are widely used for studying the dynamics of these loci. We propose a minimum chi-square (MCS) method for simultaneous estimation of all the parameters in a stepwise mutation model and the ancestral allelic type of a sample. The MCS estimator requires knowing the mean number of alleles of a certain size in a sample, which can be estimated using Monte Carlo samples generated by a coalescent algorithm. The method is applied to samples of seven (CA)n repeat loci from eight human populations and one chimpanzee population. The estimated values of parameters suggest that there is a general tendency for microsatellite alleles to expand in size, because (1) each mutation has a slight tendency to cause size increase and (2) the mean size increase is larger than the mean size decrease for a mutation. Our estimates also suggest that most of these CA-repeat loci evolve according to multistep mutation models rather than single-step mutation models. We also introduced several quantities for measuring the quality of the estimation of ancestral allelic type, and it appears that the majority of the estimated ancestral allelic types are reasonably accurate. Implications of our analysis and potential extensions of the method are discussed.SINCE the discovery that a large number of loci with tandemly repeated sequences in human and many eukaryote species are highly polymorphic because of copy number variation of the repeats in different individuals (Jeffreys 1985; Litt and Luty 1989; Weber and May 1989), allele size data from such loci are rapidly becoming the dominant source of genetic markers for genome mapping, forensic testing, and population studies. Loci with repeat sequences longer than 5 bp are generally referred to as minisatellite or variable number tandem repeat loci, and those with repeat sequences between 2 to 5 bp are referred to as microsatellite or short tandem repeat loci (Tautz 1993). Because mutations change the copy number of such loci in a stepwise fashion, rapid accumulation of population samples from minisatellite and microsatellite loci has resurrected the interest of the stepwise mutation model (SMM), which was popular in the 1970s.     

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

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

Temporal changes in allele frequencies in two reciprocally selected maize populations

The effects of breeding on allele frequency changes at 82 restriction fragment length polymorphism (RFLP) loci were examined in two maize (Zea mays L.) populations undergoing reciprocal recurrent selection, Iowa Stiff Stalk Synthetic and Iowa Corn Borer Synthetic #1. After 12 cycles of selection, approximately 30% of the alleles were extinct and 10% near fixation in each population. A test of selective neutrality identified several loci in each population whose allele frequency changes cannot be explained by genetic drift; interpopulation mean expected heterozygosity increased for that subset of 28 loci but not for the remaining 54 loci. Mean expected heterozygosity within the two subpopulations decreased 39%, while the between-population component of genetic variation increased from 0.5% to 33.4% of the total. Effective population size is a key parameter for discerning allele frequency changes due to genetic drift versus those resulting from selection and genetic hitchhiking. Empirical estimates of effective population size for each population were within the range predicted by the breeding method. Received: 10 June 1998 / Accepted: 29 April 1999     

Human microsatellites: mutation and evolution

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

High density of long dinucleotide microsatellites in Drosophila subobscura

We isolated 96 dinucleotide repeats with five or more tandemly repeated units from a subgenomic Drosophila subobscura library. The mean repeat unit length of microsatellite clones in D. subobscura is 15, higher than that observed in other Drosophila species. Population variation was assayed in 32-40 chromosomes from Barcelona, Spain, using 18 randomly chosen microsatellite loci. Positive correlation between measures of variation and perfect repeat length measures (mean size, most common, and longest allele) is consistent with a higher mutation rate in loci with longer repeat units. Levels of microsatellite variation measured as variance in repeat number and heterozygosity in D. subobscura were similar to those of Drosophila pseudoobscura and higher than those of Drosophila melanogaster and Drosophila simulans. Our data suggest that higher levels of microsatellite variation, and possibly density, in D. subobscura compared with D. melanogaster are due to both a higher average effective population and a higher intrinsic slippage rate in the former species.     

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 five trimeric and tetrameric tandem repeat loci in four human population groups.

Trimeric and tetrameric short tandem repeats (STRs) represent a rich source of highly polymorphic markers in the human genome that may be studied with the polymerase chain reaction (PCR). We report the analysis of a multilocus genotype survey of 97-380 chromosomes in U.S. Black, White, Mexican-American, and Asian populations at five STR loci located on chromosomes 1, 4, 11, and X. The heterozygote frequencies of the loci ranged from 0.36 to 0.91 and the number of alleles from 6 to 20 for the 20 population and locus combinations. Relative allele frequencies exhibited differences between populations and unimodal, bimodal, and complex distributions. Although deviations were noted at some locus-population test combinations, genotype data from the loci were consistent overall with Hardy-Weinberg equilibrium by three tests. Population subheterogeneity within each ethnic group was not detected by two additional tests. No mutations were detected in a total of 860 meioses for two loci studied in the CEPH kindreds and five loci studied in other families. An indirect estimate of the mutation rates gave values from 2.3 x 10(-5) to 15.9 x 10(-5) for the five loci. Higher mutation rates appear to be associated with greater numbers of tandem repeats in the core motif. The most frequent genotype for all five loci combined appears to have a frequency of 7.59 x 10(-4). Together, these results suggest that trimeric and tetrameric STR loci are useful markers for the study of new mutations and genetic linkage analysis and for application to personal identification in the medical and forensic sciences.     

Variance of Neutral Genetic Variances within and between Populations for a Quantitative Character

The variances of genetic variances within and between finite populations were systematically studied using a general multiple allele model with mutation in terms of identity by descent measures. We partitioned the genetic variances into components corresponding to genetic variances and covariances within and between loci. We also analyzed the sampling variance. Both transient and equilibrium results were derived exactly and the results can be used in diverse applications. For the genetic variance within populations, sigma 2 omega, the coefficient of variation can be very well approximated as [formula: see text] for a normal distribution of allelic effects, ignoring recurrent mutation in the absence of linkage, where m is the number of loci, N is the effective population size, theta 1(0) is the initial identity by descent measure of two genes within populations and t is the generation number. The first term is due to genic variance, the second due to linkage disequilibrium, and third due to sampling. In the short term, the variation is predominantly due to linkage disequilibrium and sampling; but in the long term it can be largely due to genic variance. At equilibrium with mutation [formula: see text] where u is the mutation rate. The genetic variance between populations is a parameter. Variance arises only among sample estimates due to finite sampling of populations and individuals. The coefficient of variation for sample gentic variance between populations, sigma 2b, can be generally approximated as [formula: see text] when the number of loci is large where S is the number of sampling populations.     

Toward fully automated genotyping: allele assignment, pedigree construction, phase determination, and recombination detection in Duchenne muscular dystrophy.

Human genetic maps have made quantum leaps in the past few years, because of the characterization of > 2,000 CA dinucleotide repeat loci: these PCR-based markers offer extraordinarily high PIC, and within the next year their density is expected to reach intervals of a few centimorgans per marker. These new genetic maps open new avenues for disease gene research, including large-scale genotyping for both simple and complex disease loci. However, the allele patterns of many dinucleotide repeat loci can be complex and difficult to interpret, with genotyping errors a recognized problem. Furthermore, the possibility of genotyping individuals at hundreds or thousands of polymorphic loci requires improvements in data handling and analysis. The automation of genotyping and analysis of computer-derived haplotypes would remove many of the barriers preventing optimal use of dense and informative dinucleotide genetic maps. Toward this end, we have automated the allele identification, genotyping, phase determinations, and inheritance consistency checks generated by four CA repeats within the 2.5-Mbp, 10-cM X-linked dystrophin gene, using fluorescein-labeled multiplexed PCR products analyzed on automated sequencers. The described algorithms can deconvolute and resolve closely spaced alleles, despite interfering stutter noise; set phase in females; propagate the phase through the family; and identify recombination events. We show the implementation of these algorithms for the completely automated interpretation of allele data and risk assessment for five Duchenne/Becker muscular dystrophy families. The described approach can be scaled up to perform genome-based analyses with hundreds or thousands of CA-repeat loci, using multiple fluorophors on automated sequencers.     

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.     

Differentiation of strains of Xylella fastidiosa by a variable number of tandem repeat analysis

Short sequence repeats (SSRs) with a potential variable number of tandem repeat (VNTR) loci were identified in the genome of the citrus pathogen Xylella fastidiosa and used for typing studies. Although mono- and dinucleotide repeats were absent, we found several intermediate-length 7-, 8-, and 9-nucleotide repeats, which we examined for allelic polymorphisms using PCR. Five genuine VNTR loci were highly polymorphic within a set of 27 X. fastidiosa strains from different hosts. The highest average Nei's measure of genetic diversity (H) estimated for VNTR loci was 0.51, compared to 0.17 derived from randomly amplified polymorphic DNA (RAPD) analysis. For citrus X. fastidiosa strains, some specific VNTR loci had a H value of 0.83, while the maximum value given by specific RAPD loci was 0.12. Our approach using VNTR markers provides a high-resolution tool for epidemiological, genetic, and ecological analysis of citrus-specific X. fastidiosa strains.     

Microsatellite evolution at two hypervariable loci revealed by extensive avian pedigrees

Genealogies generated through a long-term study of superb fairy-wrens (Malurus cyaneus) were used to investigate mutation within two hypervariable microsatellite loci. Of 3,230 meioses examined at the tetranucleotide locus (Mcy micro 8), 45 mutations were identified, giving a mutation rate of 1.4%. At the dinucleotide locus (Mcy micro 4) 30 mutations were recorded from 2,750 meioses giving a mutation rate of 1.1%. Mutations at both loci primarily (80%; 60/75) involved the loss or gain of a single repeat unit. Unlike previous studies, there was no significant bias toward additions over deletions. The mutation rate at Mcy micro 8 increased with allele size, and very long alleles (>70 repeats) mutated at a rate of almost 20%. The length of the mutating allele and allele span, however, were strongly correlated so it was not possible to isolate the causative factor. Allele size did not appear to affect mutation rate at Mcy micro 4, but the repeat number was considerably lower at this locus. The gender of the mutating parent was significant only at Mcy micro 8, where mutations occurred more frequently in maternal alleles. However, at both loci we found that alleles inherited from the mother were on average larger than those from the father, and this in part drove the higher mutation rate among maternally inherited alleles at Mcy micro 8.     

Detecting population growth, selection and inherited fertility from haplotypic data in humans

The frequency of a rare mutant allele and the level of allelic association between this allele and one or several closely linked markers are frequently measured in genetic epidemiology. Both quantities are related to the time elapsed since the appearance of the mutation in the population and the intrinsic growth rate of the mutation (which may be different from the average population growth rate). Here, we develop a method that uses these two kinds of genetic data to perform a joint estimation of the age of the mutation and the minimum growth rate that is compatible with its present frequency. In absence of demographic data, it provides a useful estimate of population growth rate. When such data are available, contrasts among estimates from several loci allow demographic processes, affecting all loci similarly, to be distinguished from selection, affecting loci differently. Testing these estimates on populations for which data are available for several disorders shows good congruence with demographic data in some cases whereas in others higher growth rates are obtained, which may be the result of selection or hidden demographic processes.     

Biased mutations and microsatellite variation

Mutation bias is one of the forces that may constrain the variation at microsatellite loci. Here, we study the dynamics of population statistics and the genetic distance between two populations under multiple stepwise mutations with linear bias and random drift. Expressions are derived for these statistics as functions of time, as well as at mutation-drift equilibrium. Applying these expressions to published data on humans and chimpanzees, the regression coefficient of mutation bias on allele size was estimated to be at least between - 0.0064 and -0.013. The assumption of mutational bias produces larger estimates of divergence times than are obtained in its absence; in particular, the time of split between African and non-African human populations is estimated to be between 183,000 and 222,000 years, assuming one-step mutations and no selection. With multistep mutations, the divergence time is estimated to be lower.      

Demographic Genetics of Brown Trout (Salmo Trutta) and Estimation of Effective Population Size from Temporal Change of Allele Frequencies

We studied temporal allele frequency shifts over 15 years and estimated the genetically effective size of four natural populations of brown trout (Salmo trutta L.) on the basis of the variation at 14 polymorphic allozyme loci. The allele frequency differences between consecutive cohorts were significant in all four populations. There were no indications of natural selection, and we conclude that random genetic drift is the most likely cause of temporal allele frequency shifts at the loci examined. Effective population sizes were estimated from observed allele frequency shifts among cohorts, taking into consideration the demographic characteristics of each population. The estimated effective sizes of the four populations range from 52 to 480 individuals, and we conclude that the effective size of natural brown trout populations may differ considerably among lakes that are similar in size and other apparent characteristics. In spite of their different effective sizes all four populations have similar levels of genetic variation (average heterozygosity) indicating that excessive loss of genetic variability has been retarded, most likely because of gene flow among neighboring populations.