Experimental estimation of mutation rates in a wheat population with a gene genealogy approach

Microsatellite markers are extensively used to evaluate genetic diversity in natural or experimental evolving populations. Their high degree of polymorphism reflects their high mutation rates. Estimates of the mutation rates are therefore necessary when characterizing diversity in populations. As a complement to the classical experimental designs, we propose to use experimental populations, where the initial state is entirely known and some intermediate states have been thoroughly surveyed, thus providing a short timescale estimation together with a large number of cumulated meioses. In this article, we derived four original gene genealogy-based methods to assess mutation rates with limited bias due to relevant model assumptions incorporating the initial state, the number of new alleles, and the genetic effective population size. We studied the evolution of genetic diversity at 21 microsatellite markers, after 15 generations in an experimental wheat population. Compared to the parents, 23 new alleles were found in generation 15 at 9 of the 21 loci studied. We provide evidence that they arose by mutation. Corresponding estimates of the mutation rates ranged from 0 to 4.97 x 10(-3) per generation (i.e., year). Sequences of several alleles revealed that length polymorphism was only due to variation in the core of the microsatellite. Among different microsatellite characteristics, both the motif repeat number and an independent estimation of the Nei diversity were correlated with the novel diversity. Despite a reduced genetic effective size, global diversity at microsatellite markers increased in this population, suggesting that microsatellite diversity should be used with caution as an indicator in biodiversity conservation issues.     

Factors affecting polymorphism at microsatellite loci in bread wheat [<Emphasis Type="Italic">Triticum aestivum</Emphasis> (L.) Thell]: effects of mutation processes and physical distance from the centromere

The effects of factors known to influence the level of polymorphism at microsatellite loci were studied using 99 markers and seven lines of bread wheat. Mutational factors as well as indirect selective events shape diversity at these loci. Theory predicts that the selection of favorable alleles should reduce polymorphism at neutral neighboring loci in genomic areas with low recombination rates. In wheat, local recombination rate is positively correlated with physical distance from the centromere. Seventy four loci among the 99 used could be physically located on the chromosome. We studied how the following affected the diversity among a set of inbred lines: the length of the alleles, the motif (CA versus CT), the structure of the loci (perfect versus imperfect) and the chromosomal position of the loci. For each locus, we determined whether the polymorphism observed at a locus was compatible with the Stepwise Mutation Model (SMM) or the Two-Phase Model (TPM). Both the mutation rate and the compatibility with the SMM or the TPM were shown to be variable between loci. Wheat microsatellite loci were found to be more variable when segregating alleles were perfect and had long motifs (composed of many repetitions). Diversity observed at 19 loci was not compatible with the SMM. Loci located in distal regions, with presumably high recombination rates, had longer allele sizes and were more polymorphic than loci located in proximal regions. We conclude that both mutation factors and indirect selective events vary according to the local recombination rate and therefore jointly influence the level of polymorphism at microsatellite loci in wheat.Communicated by J. Dvorak     

Reciprocal translocation of small numbers of inbred individuals rescues immunogenetic diversity

Genetic rescue can reduce inbreeding depression and increase fitness of small populations, even when the donor populations are highly inbred. In a recent experiment involving two inbred island populations of the New Zealand South Island robin, Petroica australis, reciprocal translocations improved microsatellite diversity and individual fitness. While microsatellite loci may reflect patterns of genome‐wide diversity, they generally do not indicate the specific genetic regions responsible for increased fitness. We tested the effectiveness of this reciprocal translocation for rescuing diversity of two immunogenetic regions: Toll‐like receptor (TLR) and major histocompatibility complex (MHC) genes. We found that the relatively small number of migrants (seven and ten per island) effectively brought the characteristic TLR gene diversity of each source population into the recipient population. However, when migrants transmitted TLR alleles that were already present at high frequency in the recipient population, it was possible for offspring of mixed heritage to have decreased gene diversity compared to recipient population diversity prior to translocation. In contrast to TLRs, we did not observe substantial changes in MHC allelic diversity following translocation, with limited evidence of a decrease in differentiation, perhaps because most MHC alleles were observed at both sites prior to the translocation. Overall, we conclude that small numbers of migrants may successfully restore the diversity of immunogenetic loci with few alleles, but that translocating larger numbers of animals would provide additional opportunity for the genetic rescue of highly polymorphic immunity regions, such as the MHC, even when the source population is inbred.     

Microsatellite variation and evolution in the Mimulus guttatus species complex with contrasting mating systems

Mutational variability at microsatellite loci is shaped by both population history and the mating system. In turn, alternate mating systems in flowering plants can resolve aspects of microsatellite loci evolution. Five species of yellow monkeyflowers (Mimulus sect. Simiolis) differing for historical rates of inbreeding were surveyed for variation at six microsatellite loci. High levels of diversity at these loci were found in both outcrossing and selfing taxa. In line with allozyme studies, inbreeders showed more partitioning of diversity among populations, and diversity in selfing taxa was lower than expected from reductions in effective population size due to selfing alone, suggesting the presence of either population bottlenecks or background selection in selfers. Evaluation of the stepwise mutation model (a model of DNA replication slippage) suggests that these loci evolve in a stepwise fashion. Inferred coalescent times of microsatellite alleles indicate that past bottlenecks of population size or colonization events were important in reducing diversity in the inbreeding taxon.      

Mutation and evolution of microsatellite loci in Neurospora

The patterns of mutation and evolution at 13 microsatellite loci were studied in the filamentous fungal genus Neurospora. First, a detailed investigation was performed on five microsatellite loci by sequencing each microsatellite, together with its nonrepetitive flanking regions, from a set of 147 individuals from eight species of Neurospora. To elucidate the genealogical relationships among microsatellite alleles, repeat number was mapped onto trees constructed from flanking-sequence data. This approach allowed the potentially convergent microsatellite mutations to be placed in the evolutionary context of the less rapidly evolving flanking regions, revealing the complexities of the mutational processes that have generated the allelic diversity conventionally assessed in population genetic studies. In addition to changes in repeat number, frequent substitution mutations within the microsatellites were detected, as were substitutions and insertion/deletions within the flanking regions. By comparing microsatellite and flanking-sequence divergence, clear evidence of interspecific allele length homoplasy and microsatellite mutational saturation was observed, suggesting that these loci are not appropriate for inferring phylogenetic relationships among species. In contrast, little evidence of intraspecific mutational saturation was observed, confirming the utility of these loci for population-level analyses. Frequency distributions of alleles within species were generally consistent with the stepwise mutational model. By comparing variation within species at the microsatellites and the flanking-sequence, estimated microsatellite mutation rates were approximately 2500 times greater than mutation rates of flanking DNA and were consistent with estimates from yeast and fruit flies. A positive relationship between repeat number and variance in repeat number was significant across three genealogical depths, suggesting that longer microsatellite alleles are more mutable than shorter alleles. To test if the observed patterns of microsatellite variation and mutation could be generalized, an additional eight microsatellite loci were characterized and sequenced from a subset of the same Neurospora individuals.     

Genetic diversity and bottleneck analysis of Indian bellary sheep by microsatellite markers

Bellary sheep population variability and structure was investigated genetically utilizing FAO recommended microsatellite markers. Genetic variation at 20 microsatellite loci, population structure, and genetic bottleneck hypothesis were examined. Estimates of genetic variability such as effective number of alleles and gene diversities revealed substantial genetic variation frequently displayed by microsatellite markers. A total of 133 alleles were detected. Average polymorphism across the studied loci and expected gene diversity in the population were 1.419 +/- 0.405 and 0.684 +/- 0.140, respectively. No significant genotypic linkage disequilibrium was detected across population, suggesting no evidence of linkage between loci. The population was observed to be significantly differentiated into different groups, showed fairly high level of inbreeding (f = 0.253 +/- 0.050) and global heterozygote deficit. Population structure analysis indicated the intermixing/introduction of unique/rare alleles in these migrating flocks. A normal L-shaped distribution of mode-shift test, non-significant heterozygosity excess on the basis of different models, as revealed from Sign, Standardized differences and Wilcoxon sign rank tests suggested that there was no recent bottleneck. The study revealed that even breed with increasing population trend needs genetic management for the conservation and improvement.     

Ecological-genomic diversity of microsatellites in wild barley,Hordeum spontaneum,populations in Jordan

We analyzed the ecological-genomic diversity of microsatellites of wild barley, Hordeum spontaneum (C. Koch) Thell., at 18 loci in 306 individuals of 16 populations from Jordan across a southward transect of increasing aridity. The 18 microsatellites revealed a total of 249 alleles, with an average of 13.8 alleles per locus (range 3-29), with nonrandom distribution. The proportion of polymorphic loci per population averaged 0.91 (range 0.83-1.00); gene diversity, He, averaged 0.512 (range 0.38-0.651). We compared the number of alleles of the 18 loci to those found in Israel populations by Turpeinen et al. Out of the 280 alleles, 138 (49.3%) were unique (i.e. occurred in only one of the countries). The percentage of unique alleles in Jordan and Israel populations was 43.0% and 17.9%, respectively, suggesting that Jordan is an important center of origin and diversity of wild barley. Estimates of mean gene diversity were highest in the populations collected near the Golan Heights, such as Shuni North, Shuni South and Jarash. Sixty nine percent of the microsatellite variation was partitioned within populations and 31% between populations. Associations between ecogeographical values and gene diversity were established for eight microsatellite loci. The cluster produced by simple sequence repeat (SSR) data is mostly coincidence with the result of the dendrogram of the Spalax ehrenbergi superspecies of subterranean mole rats in Jordan based on allozyme gene loci. The major soil type in the wild barley habitat of each ecological group was different. Stepwise multiple regression analysis indicated that the variance of gene diversity was explained by altitude (R(2) = 0.362**). These observations suggest that microsatellites are at least partly adaptive and subject to natural selection.     

Population structure of the African savannah elephant inferred from mitochondrial control region sequences and nuclear microsatellite loci

Two hundred and thirty-six mitochondrial DNA nucleotide sequences were used in combination with polymorphism at four nuclear microsatellite loci to assess the amount and distribution of genetic variation within and between African savannah elephants. They were sampled from 11 localities in eastern, western and southern Africa. In the total sample, 43 haplotypes were identified and an overall nucleotide diversity of 2.0% was observed. High levels of polymorphism were also observed at the microsatellite loci both at the level of number of alleles and gene diversity. Nine to 14 alleles per locus across populations and 44 alleles in the total sample were found. The gene diversity ranged from 0.51 to 0.72 in the localities studied. An analysis of molecular variance showed significant genetic differentiation between populations within regions and also between regions. The extent of subdivision between populations at the mtDNA control region was approximately twice as high as shown by the microsatellite loci (mtDNA F(ST) = 0.59; microsatellite R(ST) = 0.31). We discuss our results in the light of Pleistocene refugia and attribute the observed pattern to population divergence in allopatry accompanied by a recent population admixture following a recent population expansion.     

Isolation and characterization of microsatellite loci in Swietenia humilis (Meliaceae): an endangered tropical hardwood species

Microsatellites are now firmly established as an informative marker system, with increasing popularity as a tool amongst molecular ecologists. We have developed a method of constructing an enriched microsatellite library for the tropical tree species Swietenia humilis Zucc. (Meliaceae). This method is based on a precloning enrichment of SSRs using synthetic oligonucleotide probes, bound to magnetic beads and hybridizing to complementary microsatellite core sequences in digested genomic DNA. Here we describe the isolation and characterization of 10 microsatellite loci that have been used to survey the genetic diversity within a natural population of S. humilis. A total of 97 alleles were identified with an average of 9.7 alleles over all loci. Very high levels of allelic polymorphism were detected at individual loci, with 23 alleles observed at the most variable. The mean observed heterozygosity was 0.415 (range 0.038-0.815) exceeding levels of diversity detected in related species which used isozymes as the marker system. Subpopulation differentiation at a microgeographical scale was low (F_ST= 0.036) and the values of Nm, calculated from the allelic frequencies, were greater than 1 thus reflecting the extent of gene flow occurring between individual trees.     

Low diversity in the major histocompatibility complex class II DRB1 gene of the Spanish ibex, Capra pyrenaica

During the last two centuries, the Spanish ibex (Capra pyrenaica) has shown a significant demographic decline as a result of the progressive destruction of its natural habitat, disease epidemics, and uncontrolled hunting. Partial sequencing of the class II MHC DRB1 gene revealed that the Spanish ibex has remarkably low levels of genetic variation at this locus, with only six different DRB1 alleles and an observed heterozygosity of 0.429-0.579. The rates of nonsynonymous vs synonymous substitutions were significantly different in the peptide-binding region (dN/dS=5.347, P=0.002), a feature that indicates that the DRB1 gene is under positive selection. A phylogenetic analysis of the Spanish ibex and a set of domestic goat DRB1 alleles revealed that the reported sequences represent four major allelic lineages. The limited allelic repertoire of the DRB1 gene in the Spanish ibex is likely the direct result of the recent history of population bottlenecks and marked demographic decline of this species. A genetic survey of 13 microsatellite loci was consistent with this idea. The Spanish ibex subspecies C. p. hispanica and C. p. victoriae consistently showed considerably lower levels of microsatellite heterozygosity (Ho=0.184-0.231) and allelic diversity (mean number of alleles per locus=2-2.4) than those reported in other wild ruminants. This study demonstrates the significance of both natural selection and the demographic history of populations in determining patterns of genetic variation at MHC loci. In addition, our results emphasize the importance of locally adapted populations for the preservation of genetic diversity.     

Evidence for convergent nucleotide evolution and high allelic turnover rates at the complementary sex determiner gene of Western and Asian honeybees

Our understanding of the impact of recombination, mutation, genetic drift, and selection on the evolution of a single gene is still limited. Here we investigate the impact of all these evolutionary forces at the complementary sex determiner (csd) gene that evolves under a balancing mode of selection. Females are heterozygous at the csd gene and males are hemizygous; diploid males are lethal and occur when csd is homozygous. Rare alleles thus have a selective advantage, are seldom lost by the effect of genetic drift, and are maintained over extended periods of time when compared with neutral polymorphisms. Here, we report on the analysis of 17, 19, and 15 csd alleles of Apis cerana, Apis dorsata, and Apis mellifera honeybees, respectively. We observed great heterogeneity of synonymous (piS) and nonsynonymous (piN) polymorphisms across the gene, with a consistent peak in exons 6 and 7. We propose that exons 6 and 7 encode the potential specifying domain (csd-PSD) that has accumulated elevated nucleotide polymorphisms over time by balancing selection. We observed no direct evidence that balancing selection favors the accumulation of nonsynonymous changes at csd-PSD (piN/piS ratios are all <1, ranging from 0.6 to 0.95). We observed an excess of shared nonsynonymous changes, which suggest that strong evolutionary constraints are operating at csd-PSD resulting in the independent accumulation of the same nonsynonymous changes in different alleles across species (convergent evolution). Analysis of csd-PSD genealogy revealed relatively short average coalescence times ( approximately 6 Myr), low average synonymous nucleotide diversity (piS < 0.09), and a lack of trans-specific alleles that substantially contrasts with previously analyzed loci under strong balancing selection. We excluded the possibility of a burst of diversification after population bottlenecking and intragenic recombination as explanatory factors, leaving high turnover rates as the explanation for this observation. By comparing observed allele richness and average coalescence times with a simplified model of csd-coalescence, we found that small long-term population sizes (i.e., N(e) < 10(4)), but not high mutation rates, can explain short maintenance times, implicating a strong historical impact of genetic drift on the molecular evolution of highly social honeybees.     

Comparison of microsatellite and blood protein diversity in sheep: inconsistencies in fragmented breeds

Finnsheep, Romanov, Oxford Down and three local breeds from Finland or northwestern Russia were assessed at 15 microsatellite and 7 protein loci. A novel albumin allele was identified. Diversity patterns were mostly concordant between marker types, but discrepancies appeared for the local Viena and Vepsia sheep, both demonstrating frequent linkage disequilibria for both marker types and excess of homozygotes for microsatellites, and in the case of Vepsia also for proteins as signs of breed fragmentation. On the basis of microsatellite data, the neighbour-joining tree and two-dimensional map constructed from DA distances suggested that difference in longitude of breed origin would relate to breed relationship, whereas on the basis of protein data latitude would have this quality. These different impressions resulted because genetic distances involving Vepsia sheep were relatively low for protein variation compared with microsatellites. Microsatellite variation correlated positively with protein variation, but for the local Viena sheep protein variation was comparatively low. Populations had significant differences in allelic richness, but not in genetic diversity. Analysis implied that at least 30 polymorphic loci were needed to detect a difference in diversity between populations using a paired t-test, if the true mean diversity difference was 0.2. In the total sample, proteins demonstrated larger theta-values, but this was reversed for Finnsheep, for which model-based clustering of microsatellite genotypes revealed a structure associated with coat colour. Imported and rare sheep exhibited lowered allelic variability and increased frequency of pairwise disequilibria between unlinked markers. Our results emphasize that more loci are required for studying fragmented breeds.     

Selection and recombination drive the evolution of MHC class II DRB diversity in ungulates

Major histocompatibility complex (MHC) antigen-presenting genes are the most variable loci in vertebrate genomes. Host-parasite co-evolution is assumed to maintain the excessive polymorphism in the MHC loci. However, the molecular mechanisms underlying the striking diversity in the MHC remain contentious. The extent to which recombination contributes to the diversity at MHC loci in natural populations is still controversial, and there have been only few comparative studies that make quantitative estimates of recombination rates. In this study, we performed a comparative analysis for 15 different ungulates species to estimate the population recombination rate, and to quantify levels of selection. As expected for all species, we observed signatures of strong positive selection, and identified individual residues experiencing selection that were congruent with those constituting the peptide-binding region of the human DRB gene. However, in addition for each species, we also observed recombination rates that were significantly different from zero on the basis of likelihood-permutation tests, and in other non-quantitative analyses. Patterns of synonymous and non-synonymous sequence diversity were consistent with differing demographic histories between species, but recent simulation studies by other authors suggest inference of selection and recombination is likely to be robust to such deviations from standard models. If high rates of recombination are common in MHC genes of other taxa, re-evaluation of many inference-based phylogenetic analyses of MHC loci, such as estimates of the divergence time of alleles and trans-specific polymorphism, may be required.     

三江源和祁连山国家公园雪豹种群的遗传结构分析

掌握遗传信息对濒危物种的保护和管理具有重要意义。本研究在我国雪豹重要分布区祁连山和三江源国家公园分别采集粪便样品,利用mtDNA的cyt b基因、微卫星多态性位点进行了雪豹的物种鉴定、个体识别和种群遗传结构评估。在采集286份疑似雪豹粪便样品中,成功的对86份雪豹样品进行了扩增鉴定,利用微卫星位点进行个体识别获得41只雪豹个体,其中祁连山国家公园26只,三江源国家公园15只。通过等位基因数、有效等位基因数、观测杂合度、期望杂合度、多态信息含量等指标进行种群遗传多样性评估,认为雪豹种群遗传多样性相对较低,但祁连山国家公园雪豹种群遗传多样性相对较高。STRUCTURE进行群体遗传结构分析表明,4个种群可以划分为3个遗传类群,祁连山国家公园的种群（YCW和QLS）与三江源国家种群(DC和SJ)的遗传差异,可能与种群间的地理隔离存在明显的相似性。     

Genetic diversity and bottleneck analysis of Indian bellary sheep by microsatellite markers

Bellary sheep population variability and structure was investigated genetically utilizing FAO recommended microsatellite markers. Genetic variation at 20 microsatellite loci, population structure, and genetic bottleneck hypothesis were examined. Estimates of genetic variability such as effective number of alleles and gene diversities revealed substantial genetic variation frequently displayed by microsatellite markers. A total of 133 alleles were detected. Average polymorphism across the studied loci and expected gene diversity in the population were 1.419 ± 0.405 and 0.684 ± 0.140, respectively. No significant genotypic linkage disequilibrium was detected across population, suggesting no evidence of linkage between loci. The population was observed to be significantly differentiated into different groups, showed fairly high level of inbreeding (f = 0.253 ± 0.050) and global heterozygote deficit. Population structure analysis indicated the intermixing/introduction of unique/rare alleles in these migrating flocks. A normal L-shaped distribution of mode-shift test, non-significant heterozygosity excess on the basis of different models, as revealed from sign, standardized differences and Wilcoxon sign rank tests suggested that there was no recent bottleneck. The study revealed that even a breed with increasing population trend needs genetic management for the conservation and improvement. The text was submitted by the authors in English.     

Footprints of divergent selection in natural populations of Castanopsis fargesii (Fagaceae)

Given predicted rapid climate change, an understanding of how environmental factors affect genetic diversity in natural populations is important. Future selection pressures are inherently unpredictable, so forest management policies should maintain both overall diversity and identify genetic markers associated with the environmental factors expected to change most rapidly, like temperature and rainfall. In this study, we genotyped 648 individuals in 28 populations of Castanopsis fargesii (Fagaceae) using 32 expressed sequence tag (EST)-derived microsatellite markers. After removing six loci that departed from Hardy–Weinberg equilibrium, we measured genetic variation, population structure and identified candidate loci putatively under selection by temperature and precipitation. We found that C. fargesii populations possessed high genetic diversity and moderate differentiation among them, indicating predominant outcrossing and few restrictions to gene flow. These patterns reduce the possible impact of stochastic effects or the influence of genetic isolation. Clear footprints of divergent selection at four loci were discovered. Frequencies of five alleles at these loci were strongly correlated with environmental factors, particularly extremes in precipitation. These alleles varied from being near fixation at one end of the gradient to being completely absent at the other. Our study species is an important forest tree in the subtropical regions of China and could have a major role in future management and reforestation plans. Our results demonstrate that the gene flow is widespread and abundant in natural populations, maintaining high diversity, while diversifying selection is acting on specific genomic regions.     

Nucleotide variation at Msn and Alas2, two genes flanking the centromere of the X chromosome in humans

The centromeric region of the X chromosome in humans experiences low rates of recombination over a considerable physical distance. In such a region, the effects of selection may extend to linked sites that are far away. To investigate the effects of this recombinational environment on patterns of nucleotide variability, we sequenced 4581 bp at Msn and 4697 bp at Alas2, two genes situated on either side of the X chromosome centromere, in a worldwide sample of 41 men, as well as in one common chimpanzee and one orangutan. To investigate patterns of linkage disequilibrium (LD) across the centromere, we also genotyped several informative sites from each gene in 120 men from sub-Saharan Africa. By studying X-linked loci in males, we were able to recover haplotypes and study long-range patterns of LD directly. Overall patterns of variability were remarkably similar at these two loci. Both loci exhibited (i) very low levels of nucleotide diversity (among the lowest seen in the human genome); (ii) a strong skew in the distribution of allele frequencies, with an excess of both very-low and very-high-frequency derived alleles in non-African populations; (iii) much less variation in the non-African than in the African samples; (iv) very high levels of population differentiation; and (v) complete LD among all sites within loci. We also observed significant LD between Msn and Alas2 in Africa, despite the fact that they are separated by approximately 10 Mb. These observations are difficult to reconcile with a simple demographic model but may be consistent with positive and/or purifying selection acting on loci within this large region of low recombination.     

On the genealogy of a duplicated microsatellite

When a microsatellite locus is duplicated in a diploid organism, a single pair of PCR primers may amplify as many as four distinct alleles. To study the evolution of a duplicated microsatellite, we consider a coalescent model with symmetric stepwise mutation. Conditional on the time of duplication and a mutation rate, both in a model of completely unlinked loci and in a model of completely linked loci, we compute the probabilities for a sampled diploid individual to amplify one, two, three, or four distinct alleles with one pair of microsatellite PCR primers. These probabilities are then studied to examine the nature of their dependence on the duplication time and the mutation rate. The mutation rate is observed to have a stronger effect than the duplication time on the four probabilities, and the unlinked and linked cases are seen to behave similarly. Our results can be useful for helping to interpret genetic variation at microsatellite loci in species with a very recent history of gene and genome duplication.     

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.