*Department of Genetics, University of Wisconsin, Madison;
Crop Genetics Research and Development, DuPont Agriculture and Nutrition, Pioneer Hi-Bred International, Johnston, Iowa;
Fukui Prefectural University, Matsuoka-cho, Yoshida-gun, Fukui, Japan;
National Center for Genome Resources, Santa Fe, New Mexico
![]() |
Abstract |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Introduction |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Numerous studies in animals have calculated the mutation rate of microsatellites, and these studies have shown that the mutation rate varies greatly among species, ranging from 5 x 10-6 in Drosophila (Schug, Mackay, and Aquadro 1997
; Schlötterer et al. 1998
; Schug et al. 1998
; Vazquez et al. 2000
) to 10-3 in human (Brinkmann et al. 1998
; Xu et al. 2000
). These studies have also outlined several trends of the mutation process that are important for understanding microsatellite evolution. They have shown that the mutation rate varies widely among loci within species (Di Rienzo et al. 1998
; Harr et al. 1998
) and that the mutation rate increases with the length of the microsatellite (Primmer et al. 1996
; Brinkmann et al. 1998
). It has also been shown that there is a constraint on the size of microsatellites (Garza, Slatkin, and Freimer 1995
) which may simply be the effect of the increased probability of contraction with the size of the microsatellite (Ellegren 2000
; Harr and Schlötterer 2000
; Xu et al. 2000
).
Mutation rate is a critical parameter in population genetic models because it enables one to relate the variability at microsatellite loci to the history of a population or to the history of a portion of the genome. For example, under the hypothesis of a generalized stepwise model of mutation, knowledge of the mutation rate permits one to estimate the time of divergence between species (Wehrhahn 1975
; Goldstein et al. 1995
) and the effective population size of species (Slatkin 1995
).
In this article, we report the mutation rate for maize (Zea mays subsp. mays) microsatellites as determined in two different experiments of mutation accumulation involving a total of 142 microsatellites. We also describe the nature of the mutation process, including whether there is a bias toward mutations that increase versus decrease allele size, whether the number of repeats in the progenitor allele is correlated with the mutation rate, and whether the mutation rate differs among loci with dinucleotide versus trinucleotide or higher-repeat motifs. Finally, we apply the mutation rate that we have determined to estimate the effective population sizes of maize and its wild progenitor (Z. mays subsp. parviglumis).
![]() |
Materials and Methods |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
In experiment II, 86 recombinant inbred (RI) lines were studied. These lines were derived from two different crosses, T232 x CM37 (45 lines) and Co159 x Tx303 (41 lines) (Burr et al. 1988
). The F2 progeny of a single F1 plant from each cross were selfed for 9 to 12 generations. The average number of generations of inbreeding after the F1 was 11.3 generations for the T232 x CM37 cross and 11.0 generations for the Co159 x Tx303 cross. The 86 RI lines and the four parental inbreds were all genotyped and compared for evidence of new mutations. Our genotyping of these 86 RI lines with 98 microsatellite loci revealed no plants with nonparental alleles at multiple loci, as expected if there had been pollen contamination.
Microsatellite Genotyping
For experiment I, DNA was extracted from leaf tissue of 10-day-old seedlings or from freeze-dried tissue of young seedlings (Smith et al. 1997
). DNAs of each progeny were allocated into duplicate 96-well liquid handling plates. All PCR amplifications and gel runs were made in duplicate for each progeny, and additional replicates of the parents were amplified and electrophoresed for each microsatellite. Forty-eight microsatellite loci were used and represented a variety of repeat motif types: 6 with dinucleotide repeat motifs, 21 with trinucleotide, 15 with tetranucleotide, 4 with pentanucleotide, 1 with hexanucleotide, and 1 with a di-tetra motif (table 1
). Microsatellite genotyping was performed on ABI automated sequencers, using procedures that have been described previously (Smith et al. 1997
). DNA samples showing putative mutants were amplified a third time for the microsatellite in question.
|
|
In experiment II, the mutant alleles and the two parental inbreds for its RI line were sequenced. Microsatellite loci were amplified by PCR using PCR Supermix (BRL) and 10 pmol of each primer, as described by Matsuoka et al. (2002a)
. Two PCRs were performed for each analysis, and the combined PCR products were purified on Qiagen columns. For lines in which the mutant allele was homozygous, the PCR products were directly sequenced with the primers used for the PCR. For heterozygous individuals, the PCR products were cloned into a plasmid vector, using the TOPO cloning kit (Invitrogen). Plasmid clones (from two to eight) were sequenced for each mutation. The sequencing was performed using M13 forward and reverse primers with a BigDye terminator sequencing kit (ABI) at the University of Wisconsin Biotechnology Center (Madison).
Statistics
To estimate the mutation rate, one divides the number of observed mutations by the number of independent generations that the two alleles present at the last generation have experienced (i.e., the number of allele-generations). For experiment I which examined only a single generation, the number of allele-generations is simply equal to 23 kernels x 2 alleles x 4 ears x 6 inbreds x 48 loci (=52,992) minus any missing data. For experiment II, one might consider that the number of allele-generations is simply the number of generations (g) times 2 (for a diploid). However, a mutation appearing in the heterozygous state in an early generation has a certain probability to be lost by drift during the successive generation of selfing. So, the number of allele-generations is somewhat less than 2g. Thus, the number of allele-generations must be determined using the probability of coalescence of the alleles over the 11 generations since the F1.
We calculated the number of allele-generations for experiment II as follows. From the last generation, going backward one generation, the probability that the two alleles at this last generation coalesced in the generation before the last is 1/2 because the plants were self-pollinated. In this case, the total number of allele-generations that the two alleles have experienced is g + 1. If the two alleles have not coalesced at the generation before the last generation, then the probability that they coalesced two generations before the last (knowing that they have not coalesced one generation before the last generation) is 1/4. In this case, the number of allele-generations is g + 2. This process can be extended up to the F2 generation at which point the probability that the two alleles have coalesced, knowing that they have not coalesced elsewhere, is (1/2)g-1, with the number of allele-generations being 2g - 1. Finally, the probability that the two alleles do not coalesce is (1/2)g-1, with the number of allele-generations being 2g. A general formula for the expected number of allele-generations is:
|
Two main parameters are important in the mutation process of microsatellites: the mutation rate, µ, and the variance of change in the number of repeats among mutations, m2 (Slatkin 1995
; Zhivotovsky and Feldman 1995
). To estimate
m2, one needs to know the change in the number of repeats for each mutation from its progenitor allele. For experiment II, when sequence polymorphisms in the regions flanking the microsatellite repeat enabled us to identify the progenitor allele, we inferred the change in the number of repeat units for the mutation by comparison with the progenitor allele. In cases where the progenitor allele could not be unambiguously identified, the change in the number of repeat units for the mutation was inferred by comparison with the parental allele that was most similar in size to the mutant allele.
Even when no mutations are observed, it is possible to calculate an upper limit of the mutation rate using the Poisson law. The probability of zero mutations is P(X = 0) = e-Gµ. We can solve this equation for P(X = 0) = 0.05 to obtain the upper limit of the mutation rate (Schug, Mackay, and Aquadro 1997
). The 95% confidence interval for the observed mutation rate can also be calculated using the Poisson law. The probability of observing k or fewer mutations with
= Gµ is
|
For experiment II, most mutations have been confirmed by two different assays (either two PCRs or one PCR plus sequencing). However, 15 of 194 putative mutations were not analyzed by a second PCR or sequencing. These may be real mutations or errors in the first PCR. To infer the proportion of the unconfirmed putative mutations that are predicted to be real mutations, we used the proportion of real mutations among those putative mutations that had been subjected to two assays. Before doing this, we divided these unconfirmed putative mutations into two classes: those that differ by 1 bp from a parental allele and those that differ by more than 1 bp from the parental allele.
![]() |
Results |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Of the 98 microsatellites analyzed in experiment II, several failed in one of the two RI populations. Loci bnlg1839, bnlg2086, and phi096 failed for the T232 x CM37 population, and bngl1074, bnlg1257, and phi064 failed for the Co159 x Tx303 population. There were also 528 cases where the PCR failed for one or a few plants but worked well for the population as a whole. Excluding these missing data, a total of 7,683 successful PCRs were performed (table 2 ). Of these, 194 gave a nonparental allele or potential mutation. Of these 194 potential mutations, 31 differed by 1 bp from the parental allele and 163 by 2 bp or more. A second PCR or sequencing was performed on 179 of these potential mutant alleles. Of these, 72 were confirmed and were classified as real mutations. Of these 72 real mutations, 3 differed by only 1 bp from the parental allele and 69 by more than 1 bp (table 2 , supplementary material available at MBE web site: www.molbiolevol.org).
The nc009 locus gave an unexpected result for six RI lines of the Co159 x Tx303 population, each of which possessed the same nonparental allele of 133 bp. These six nonparental alleles were sequenced, and all have the same number of repeats (AG)18. Sequence polymorphism in the 3' flanking region identified the 151-bp parental allele as the progenitor of the 133-bp nonparental allele, indicating that there was an 18-bp deletion or loss of nine repeat units. Because one F1 plant was used to form the F2 population, it seems most likely that a premeiotic somatic mutation in either the ear or tassel cell lineage of the F1 plant gave rise to the 133-bp allele which was then inherited by the six RI lines. This type of event has been reported previously (Jones et al. 1999
). For this reason, we have interpreted this result as a single mutation rather than six independent mutations.
To further analyze the nature of the mutations, 54 of the 72 mutations from experiment II were sequenced. All mutations confirmed by two PCRs were again verified by sequencing. Of the three 1-bp mutants confirmed by two PCRs, all were 1-bp changes in the length of a mononucleotide tract flanking the microsatellite.
Estimation of the Mutation Rate
In experiment I, a single mutation was observed in one of the seven dinucleotide loci assayed. The number of allele-generations in this case is 7,718 after subtracting 10 missing data points. This yields a mutation rate per generation for dinucleotide microsatellites of 1.3 x 10-4 (table 3
). The 95% confidence interval for this rate is 3.1 x 10-5 to 7.2 x 10-4. For microsatellites with repeat units of greater than 2 bp in length, no mutations were observed, and thus we can only calculate the upper bound of the mutation rate. Here, the number of allele-generations is 44,568 which gives an upper bound of the mutation rate of µ = 6.7 x 10-5.
|
In experiment II, one might also consider mutations that were hidden because one parental allele mutated to the other parental allele. This problem would only be significant when the two parental alleles differ by one repeat because mutations of one repeat are the most common. Twenty-one of the loci-by-RI combinations have a one-repeat difference between the two parental alleles, and among these, we observed three one-repeat mutations. Because 50% of the mutations will be to the other parental allele and thus hidden, we estimate three hidden mutations. This will only increase the mutation rate from 8.3 x 10-4 to 8.7 x 10-4.
If experiments I and II are grouped together by adding up the number of mutations and the number of allele-generations, the mutation rate for dinucleotide microsatellites is 7.7 x 10-4, with a confidence interval of 5.2 x 10-4 to 1.1 x 10-3 (table 3 ). For microsatellites with repeats of more than 2 bp in length, the upper bound of the mutation rate for the combined experiments is 5.1 x 10-5, considering that 59,100 allele-generations were analyzed.
We also used the data from experiment II to calculate that the variance of change in the number of repeats (m2) is 3.17 (table 4 ). This parameter is strongly influenced by outliers with a large change in the number of repeats. For example, if the mutation in nc009 that caused a decrease of nine repeats is excluded, then
m2 drops to 2.03. Accordingly, the value given here should be taken with some caution.
|
|
|
|
![]() |
Discussion |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
One clear factor governing the rate variation among loci within maize is the length of the repeat motif. We observed no mutations among loci with repeat motifs of more than 2 bp in length, as compared with 70 mutations among loci with dinucleotide repeats. This result is consistent with what has been seen in some other organisms, including Drosophila (Schug et al. 1998
) and humans (Chakraborty et al. 1997
). Nevertheless, there are counterexamples where the opposite relationship has been observed such that loci with dinucleotide repeat motifs mutate at a slower rate than those with larger repeat motifs (Weber and Wong 1993
; Eckert and Yan 2000
). These conflicting reports highlight the idiosyncratic nature of the mutational process at microsatellites and caution against applying results from one organism to another. It is also possible that the difference in the mutation rate between di-, tri-, and tetranucleotide microsatellites can be explained by a difference in the average number of repeats (Schug et al. 1998
; Harr and Schlötterer 2000
), with dinucleotide loci having a higher average number of repeats and thus a higher mutation rate. Because we have not determined the average length in size for all the studied loci, we cannot test this hypothesis in maize.
We can also ask whether the mutation rate that we calculated for each locus applies to "natural" populations of maize. If it does, the number of observed mutations at a locus should be positively correlated with the number of alleles and heterozygosity in maize. Matsuoka et al. (2002b)
have investigated genetic diversity for the 98 loci used in experiment II among a sample of 193 maize plants. Heterozygosity and the number of alleles were estimated using this data for dinucleotide loci. We found that the number of mutations observed in our experiment is correlated with both heterozygosity (nonparametric Spearman correlation, Rs = 0.42, P < 0.001) and the number of alleles (Rs = 0.45, P < 0.001).
Trends in the Mutational Process
In addition to estimating the mutation rate, our data revealed several features of the mutational process for microsatellites in maize. Similar to what has been observed in animals (Amos et al. 1996
; Primmer et al. 1996
), we found that mutations of a single repeat in length are far more common than mutations of multiple repeats. Of the mutations observed in dinucleotide microsatellites, 83% are one repeat, 13% are two repeats, and 4% are more than two repeats in length (table 4
). We also observed heterogeneity in the rate among dinucleotide loci (fig. 1
), and because our dinucleotide loci are mostly (AG)n repeats, this heterogeneity cannot be explained by diversity in the sequence composition of the repeat motif (see Bachtrog et al. 2000
). However, this heterogeneity may be partly attributed to variation in the number of repeats in the parental alleles because we observed a correlation between the number of repeats in the parental allele and the number of observed mutations (fig. 3a
). Alleles with a greater number of repeats appear to be more mutable.
We have also observed that there is a higher probability to mutate to a larger allele than to a smaller one (fig. 2
), with 56 of the 71 observed mutations from experiment II having caused an increase in the size of the allele. This same bias has been observed with other organisms (Amos et al. 1996
; Ellegren 2000
; Xu et al. 2000
). Given this bias, previous authors have asked: why don't microsatellites increase infinitely in size? One possible explanation is that there is an equilibrium between mutations that alter the size of the microsatellite and base substitutions that lead to the degradation of the microsatellite (Kruglyak et al. 1998
). In Drosophila, however, such an equilibrium process is inadequate to explain the underrepresentation of large microsatellites in the genome (Harr and Schlötterer 2000
). Another nonexclusive possible explanation is that the larger the allele, the greater the probability that a mutation will cause a contraction in size. Our observation that mutations causing a decrease in size are on average larger than those that cause an increase is consistent with this mechanism (fig. 2
), partially explaining why microsatellites do not increase infinitely in size (also see Harr and Schlötterer 2000
). However, we did not find a significant negative correlation between the standardized size of the progenitor allele and magnitude of the mutation as detected in humans (Ellegren 2000
; Xu et al. 2000
), although the largest contractions are associated with large alleles.
Short-Term Versus Long-Term Mutational Processes
In this study, we evaluated the short-term mutational pattern at microsatellite loci in maize in two different experiments of mutation accumulation. We determined the DNA sequence of the mutant and progenitor alleles for 55 of 73 new mutations. The sequence analysis revealed that all mutations were changes in the number of repeats in the microsatellite or in the length of mononucleotide tracts flanking the microsatellite. We did not observe any indels in the flanking regions, except for the aforementioned changes in mononucleotide tracts.
Contrary to our results, Matsuoka et al. (2002a)
observed that microsatellite alleles among lines of maize and teosinte typically differ by indels of 2 to 50 bp (or larger) in the regions flanking the microsatellite repeat. Because we have observed no such mutations in our short-term evolutionary study, this class of indels likely arises only over longer evolutionary periods at a rate far below our estimated rate for dinucleotide microsatellites. Using the Poisson law, the 95% upper bound for the rate for such indels is 2 x 10-5, given that we have examined a total of more than 153,000 allele-generations without observing any indels in the flanking sequences.
Microsatellites are assayed by screening for length polymorphisms in a DNA region between a pair of primers that flank the microsatellite repeat. Our results combined with those of Matsuoka et al. (2002a)
indicate that the observed length polymorphism is the result of several processes that proceed at different rates. First, there can be changes in the number of repeats in the microsatellite that can proceed at an average rate of 7.7 x 10-4 mutations/generation for dinucleotide repeat loci. Second, there can be multistep mutations with the variance of change in the number of repeats among mutations (
m2) being 3.2. Third, there can be indels of 2 to 50 bp or more in the flanking regions that accumulate at rates below 2 x 10-5 mutations/generation. This mixed mutational pattern cautions against the casual use of models based on a simple stepwise mutation process with a single mutation rate. However, with knowledge of these complexities, it is possible to identify a subset of microsatellites that more closely follow a stepwise process and have a more uniform mutation rate.
Estimation of Effective Population Sizes
For estimation of the effective population size of maize, one needs a set of markers that behave in a stepwise manner and have a known mutation rate. Matsuoka et al. (2002b)
have investigated genetic diversity for the 98 loci used in experiment II among a sample of 264 maize and teosinte plants. For that sample, 33 of the loci had allelic distributions with less than 10% nonstepwise alleles (table 1
). Using the data from our study, we calculated the mutation rate for these 33 loci to be 4.3 x 10-4 and the variance of change in the size of the allele (
m2) to be 2.08. Using the Poisson distribution, we cannot reject the null hypothesis that the rate is homogeneous among these 33 loci (
2 = 0.89, P = 0.35), and so this rate can be applied to all 33 loci.
The variance of allele size for this set of 33 microsatellite loci was calculated for 193 maize plants and 34 teosinte (Z. mays subsp. parviglumis) plants, using the data presented in Matsuoka et al. (2002b)
. For maize, the average variance of allele size was 23.5 repeats (range 0.8879.5), and for subsp. parviglumis, it was 26.8 (range 0.8085.1). With these data, we can calculate the effective population size, using a generalized stepwise model that allows for steps of more than one repeat in length (Slatkin 1995
). Under this model, the variance in allele size (
2), the effective population size (N), the variance of change in the size of the allele (
m2), and the mutation rate (µ) are related by the following formula:
![]() |
One can also estimate effective population size using the equilibrium expectation for heterozygosity (H) for microsatellite loci following a strict stepwise model (Kimura and Ohta 1975):
![]() |
This effective population size for maize was previously calculated using polymorphism at adh1 and estimated to be 660,000 (Gaut and Clegg 1993
). In a similar study of adh1, Eyre-Walker et al. (1998)
reported an estimate of the effective size of subsp. parviglumis at 940,000. Our estimates are more than an order of magnitude less than these. The adh1-based estimates assume a DNA substitution rate inferred from the amount of DNA sequence divergence accumulated over the 50 to 60-Myr history of the grass family. There is a concern about this rate because it represents a long-term rate over the history of the grasses and may not be appropriate to recent events in the maize lineage, where a lineage-specific rate acceleration can be anticipated because of a generation-time effect (White and Doebley 1999
).
In another study, Remington et al. (2001)
reported an effective population size of 200,000 for maize, based on the degree of linkage disequilibrium among maize inbred lines and the function C = 4Nc, where C is the population recombination parameter, and c is the recombination rate ([crossovers/bp] x generation). However, this report assumes a recombination rate of 10-8, whereas reported rates in maize genes are nearer to 10-7 (Patterson et al. 1995
; Xu et al. 1995
; Dooner and Martinez-Ferez 1997
; Okagaki and Weil 1997
). Using the latter rate would give a 4Nc-based estimate of 20,000, a value much closer to the estimate based on microsatellites. The value for C reported by Remington et al. (2001)
may also be inappropriate because it is based on a biased sample of the maize germ plasm pool. A value for C based on a more representative sample of maize is
0.02 (Tenaillon et al. 2001
). Using this value of C and the observed recombination rate for maize genes (
10-7), the effective population size for maize would be 50,000. The difficulty with these estimates is the uncertainty surrounding values for c. The cause of the differences in estimates of effective population size based on c, sequence polymorphism, and microsatellites will require further exploration.
This study has investigated the mutation rate and process for maize microsatellites. Because these markers are widely used in plants for a variety of purposes, such estimates of the mutation rate and knowledge of the mutation process are needed to clarify the origin and maintenance of genetic diversity at these loci. Whereas the mutation pattern of maize microsatellites is complex, a fuller understanding of these complexities will facilitate their application to a variety of questions in maize genetics and evolution.
![]() |
Acknowledgements |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Footnotes |
---|
Keywords: Zea mays subsp. mays
microsatellite
mutation rate
maize
teosinte, SSR
Address for correspondence and reprints: John Doebley, Department of Genetics, University of Wisconsin, Madison, Wisconsin 53706. jdoebley{at}facstaff.wisc.edu
![]() |
References |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Amos W., S. J. Sawcer, R. W. Feakes, D. C. Rubinstein, 1996 Microsatellites show mutational bias and heterozygote instability Nat. Genet 13:390-391[ISI][Medline]
Bachtrog D., M. Agis, M. Imhof, C. Schlötterer, 2000 Microsatellite variability differs between dinucleotide repeat motifsevidence from Drosophila melanogaster Mol. Biol. Evol 17:1277-1285
Bowcock A., J. Ruiz-Linares, E. Tomfohrde, J. Minch, J. Kidd, L. Cavalli-Sforza, 1994 High resolution human evolutionary trees with polymorphic microsatellites Nature 368:455-457[ISI][Medline]
Brinkmann B., M. Klintschar, F. Neuhuber, J. Hühne, B. Rolf, 1998 Mutation rate in human microsatellites: influence of the structure and length of the tandem repeat Am. J. Hum. Genet 62:1408-1415[ISI][Medline]
Burr B., F. A. Burr, K. H. Thompson, M. C. Alberston, C. W. Stuber, 1988 Gene mapping with recombinant inbreds in maize Genetics 118:519-526
Chakraborty R., M. Kimmel, D. Stivers, L. Davison, R. Deka, 1997 Relative mutation rates at di-, tri-, and tetranucleotide microsatellite loci Proc. Natl. Acad. Sci. USA 94:1041-1046
Dallas J. F., 1992 Estimation of the mutation rates in recombinant inbred strains of mouse Mamm. Genome 3:452-456[ISI][Medline]
Dietrich W., H. Katz, S. Lincoln, H. Shin, J. Friedman, N. Dracopoli, E. Lander, 1992 A genetic map of the mouse suitable for typing intraspecific crosses Genetics 131:423-427
Di Rienzo A., P. Donnelly, C. Toomajian, B. Sisk, A. Hill, M. L. Petzl-Erle, G. H. Haines, D. H. Barch, 1998 Heterogeneity of microsatellite mutations within and between loci, and implication for human demographic histories Genetics 148:1269-1284
Dooner H., I. Martinez-Ferez, 1997 Recombination occurs uniformly within the bronze gene, a meiotic recombination hotspot in the maize genome Plant Cell 9:1633-1646
Eckert K. A., G. Yan, 2000 Mutational analyses of dinucleotide and tetranucleotide microsatellites in Escherichia coli: influence of sequence on expansion mutagenesis Nucleic Acids Res 28:2831-2838
Ellegren H., 2000 Heterogenous mutation processes in human microsatellite DNA sequences Nat. Genet 24:400-402[ISI][Medline]
Eyre-Walker A., R. S. Gaut, H. Hilton, D. L. Feldman, B. S. Gaut, 1998 Investigating the bottleneck leading to the domestication of maize Proc. Natl. Acad. Sci. USA 95:4441-4446
Garza J. G., M. Slatkin, N. B. Freimer, 1995 Microsatellite allele frequencies in human and chimpanzees, with implications for constraints on allele size Mol. Biol. Evol 12:594-603[Abstract]
Gaut B. S., M. T. Clegg, 1993 Molecular evolution of the Adh1 locus in the genus Zea Proc. Natl. Acad. Sci. USA 90:5095-5099[Abstract]
Goldstein D. B., A. R. Linares, L. L. Cavalli-Sforza, M. W. Feldman, 1995 Genetic absolute dating based on microsatellites and the origin of modern human Proc. Natl. Acad. Sci. USA 92:11549-11552[Abstract]
Goldstein D. B., G. Roemer, D. Smith, D. Reich, A. Bergman, R. Wayne, 1999 The use of microsatellite variation to infer population structure and demographic history in a natural model system Genetics 151:797-801
Harr B., C. Schlötterer, 2000 Long microsatellite alleles in Drosophila melanogaster have a downward mutation bias and short persistence times, which cause their genome-wide under representation Genetics 155:1231-1220
Harr B., B. Zangerl, G. Brem, C. Schlötterer, 1998 Conservation of locus-specific microsatellite variability across species: a comparison of two Drosophila sibling species, D. melanogaster and D. simulans Mol. Biol. Evol 15:176-184[Abstract]
Jones A. G., G. Rosenqvist, A. Berglund, J. C. Avise, 1999 Clustered microsatellite mutations in the pipefish Syngnathus typhle Genetics 152:1057-1063
Kimura M., T. Ohta, 1975 Distribution of allelic frequencies in a finite population under stepwise production of neutral alleles Proc. Natl. Acad. Sci. USA 72:2761-2764[Abstract]
Kovalchuk O., Y. E. Dubrova, A. Arkhipov, B. Hohn, I. Kovalchuk, 2000 Wheat mutation rate after Chernobyl Nature 407:583-584[ISI][Medline]
Kruglyak S., R. T. Durrett, M. D. Schug, C. F. Aquadro, 1998 Equilibrium distributions of microsatellite repeat length resulting from a balance between slippage events and point mutations Proc. Natl. Acad. Sci. USA 95:10774-10778
Lagercrantz U., H. Ellegren, L. Andersson, 1993 The abundance of various polymorphic microsatellite motifs differs between plants and vertebrates Nucleic Acids Res 21:1111-1115[Abstract]
Matsuoka Y., S. E. Mitchell, S. Kresovich, M. Goodman, J. Doebley, 2002a. Microsatellites in Zeavariability, patterns of mutations, and use for evolutionary studies Theor. Appl. Genet 104:436-450[ISI][Medline]
Matsuoka Y., Y. Vigouroux, J. Sanchez, M. Goodman, E. S. Buckler IV, J. Doebley, 2002b. Maize domestication and diversification inferred from microsatellite DNA Proc. Natl. Acad. Sci. USA 99:60806084
Okagaki R., C. Weil, 1997 Analysis of recombination sites within the maize waxy locus Genetics 147:815-821
Patterson G., K. Kubo, T. Shroyer, V. Chandler, 1995 Sequences required for paramutation of the maize b gene map to a region containing the promoter and upstream sequences Genetics 140:1389-1406
Peters J. M., D. C. Queller, J. E. Strassmann, C. R. Solis, 1995 Maternity assignment and queen replacement in a social wasp Proc. R. Soc. Lond. Ser. B Biol. Sci 260:7-12[ISI][Medline]
Powell W., M. Morgante, C. Andre, M. Hanafey, J. Vogel, S. Tingey, A. Rafalski, 1996 The comparison of RFLP, RAPD, AFLP and SSR (microsatellite) markers for germplasm analysis Mol. Breed 2:225-238[ISI]
Primmer C. R., H. Ellegren, N. Saino, A. P. Møller, 1996 Directional evolution in germline microsatellite mutations Nat. Genet 13:391-393[ISI][Medline]
Remington D., J. Thornsberry, Y. Matsuoka, L. Wilson, S. Whitt, J. Doebley, S. Kresovich, M. Goodman, E. S. Buckler, 2001 Structure of linkage disequilibrium and phenotypic associations in the maize genome Proc. Natl. Acad. Sci. USA 98:11479-11484
Schlötterer C., R. Ritter, B. Harr, G. Brem, 1998 High mutation rate of a long microsatellite allele in Drosophila melanogaster provides evidence for allele-specific mutation rates Mol. Biol. Evol 15:1269-1274
Schug M. D., C. M. Hutter, K. A. Wetterstrand, M. S. Gaudette, T. F. C. Mackay, C. F. Aquadro, 1998 The mutation rates of di-, tri-, tetranucleotide repeats in Drosophila melanogaster Mol. Biol. Evol 15:1751-1760
Schug M. D., T. F. C. Mackay, C. F. Aquadro, 1997 Low mutation rates of microsatellite loci in Drosophila melanogaster Nat. Genet 15:99-102[ISI][Medline]
Slatkin M., 1995 A measure of population subdivision based on microsatellite allele frequencies Genetics 139:457-462
Smith J. S. C., E. C. L. Chin, H. Shu, O. S. Smith, S. J. Wall, M. L. Senior, S. E. Mitchell, S. Kresovich, J. Ziegle, 1997 An evaluation of the utility of SSR loci as molecular markers in maize (Zea mays L.): comparisons with data from RFLPs and pedigree Theor. Appl. Genet 95:163-173[ISI]
Tenaillon M., M. Sawkins, A. Long, R. Gaut, J. Doebley, B. Gaut, 2001 Patterns of DNA sequence polymorphism along chromosome 1 of maize (Zea mays ssp. mays) Proc. Natl. Acad. Sci. USA 98:9161-9166.
Udupa S. M., M. Baum, 2001 High mutation rate and mutational bias at (TAA)n microsatellite loci in chickpea (Cicer arietinum L) Mol. Genet. Genomics 265:1097-1103[ISI][Medline]
Vazquez F., T. Perez, J. Albornoz, A. Dominguez, 2000 Estimation of the mutation rates in Drosophila melanogaster Genet. Res 76:323-326[ISI][Medline]
Weber J. L., C. Wong, 1993 Mutation of human short tandem repeats Hum. Mol. Genet 2:1123-1128[Abstract]
Wehrhahn C. F., 1975 The evolution of selectively similar electrophoretically detectable alleles in finite natural populations Genetics 80:375-394
White S., J. Doebley, 1999 The molecular evolution of terminal ear1, a regulatory gene in the genus Zea Genetics 153:1455-1462
Xu X., A. Hsia, L. Zhang, B. Nikolau, P. Schnable, 1995 Meiotic recombination break points resolve at high rates at the 5' end of a maize coding sequence Plant Cell 7:2151-2161
Xu X., M. Peng, Z. Fang, X. Xu, 2000 The direction of microsatellite mutations is dependent upon the allele length Nat. Genet 24:396-399[ISI][Medline]
Zhivotovsky L. A., M. W. Feldman, 1995 Microsatellite variability and genetic distance Proc. Natl. Acad. Sci. USA 92:11549-11552[Abstract]