Mutation rates and intraspecific divergence of the mitochondrial genome of Pristionchus pacificus

Evolutionary reconstruction of the natural history of an organism ultimately requires knowledge about the development, population genetics, ecology, and phylogeny of the species. Such investigations would benefit from studies of mutational processes because mutations are the source of natural variation. The nematode Pristionchus pacificus has been developed as a model organism in evolutionary biology by comparing its development with Caenorhabditis elegans. Pristionchus pacificus and related species are associated with scarab beetles, and their ecology and phylogeny are well known. More than 200 P. pacificus isolates from all over the world are available for this cosmopolitan species. We generated mutation accumulation (MA) lines in P. pacificus to study spontaneous mutation rates in the mitochondrial genome and compared mutation rate estimates with natural variation between nine representative isolates of the species. The P. pacificus mitochondrial genome is 15,955 bp in length and is typical for nematodes. Pristionchus pacificus has all known mitochondrial genes and contains an unusual suppressor transfer RNA (tRNA) for the codon UAA. This has most likely influenced the spectrum of observable mutations because 6 of 12 mutations found in the 82 MA lines analyzed are nonsense mutations that can be suppressed by the suppressor tRNA. The overall mutation rate in P. pacificus is 7.6 × 10?? per site per generation and is less than one order of magnitude different from estimates in C. elegans and Drosophila. Using this mutation rate estimate in a comparison of the mitochondrial genome of nine P. pacificus isolates, we calculate the minimum time to the most recent common ancestor at 10?-10? generations. The combination of mutation rate analysis with intraspecific divergence provides a powerful tool for the reconstruction of the natural history of P. pacificus, and we discuss the ecological implication of these findings.     

The relationship between microsatellite slippage mutation rate and the number of repeat units

Microsatellite markers are widely used for genetic studies, but the relationship between microsatellite slippage mutation rate and the number of repeat units remains unclear. In this study, microsatellite distributions in the human genome are collected from public sequence databases. We observe that there is a threshold size for slippage mutations. We consider a model of microsatellite mutation consisting of point mutations and single stepwise slippage mutations. From two sets of equations based on two stochastic processes and equilibrium assumptions, we estimate microsatellite slippage mutation rates without assuming any relationship between microsatellite slippage mutation rate and the number of repeat units. We use the least squares method with constraints to estimate expansion and contraction mutation rates. The estimated slippage mutation rate increases exponentially as the number of repeat units increases. When slippage mutations happen, expansion occurs more frequently for short microsatellites and contraction occurs more frequently for long microsatellites. Our results agree with the length-dependent mutation pattern observed from experimental data, and they explain the scarcity of long microsatellites.     

Cumulative effects of spontaneous mutations for fitness in Caenorhabditis: role of genotype, environment and stress

It is often assumed that the mutation rate is an evolutionarily optimized property of a taxon. The relevant mutation rate is for mutations that affect fitness, U, but the strength of selection on the mutation rate depends on the average effect of a mutation. Determination of U is complicated by the possibility that mutational effects depend on the particular environmental context in which the organism exists. It has been suggested that the effects of deleterious mutations are typically magnified in stressful environments, but most studies confound genotype with environment, so it is unclear to what extent environmental specificity of mutations is specific to a particular starting genotype. We report a study designed to separate effects of species, genotype, and environment on the degradation of fitness resulting from new mutations. Mutations accumulated for >200 generations at 20 degrees in two strains of two species of nematodes that differ in thermal sensitivity. Caenorhabditis briggsae and C. elegans have similar demography at 20 degrees, but C. elegans suffers markedly reduced fitness at 25 degrees. We find little evidence that mutational properties differ depending on environmental conditions and mutational correlations between environments are close to those expected if effects were identical in both environments.     

Isolation and characterization of the first eight microsatellite loci in Gammarus fossarum (Crustacea, Amphipoda) and cross-amplification in Gammarus pulex and Gammarus orinos

We characterized the first microsatellite markers for Gammarus fossarum. Eight loci gave satisfactory amplification patterns in two stream populations (Southern France) with number of alleles ranging from 2 to 10 and expected heterozygosity from 0.076 to 0.857. We performed cross-amplification in two closely related gammarid species, Gammarus pulex and Gammarus orinos. Among the eight tested microsatellite loci, four correctly amplified in G. pulex and three in G. orinos.     

Allelic variation,fragment length analyses and population genetic models: a case study on Drosophila microsatellites*

The allelic variation of 16 microsatellite loci from selected species of the Drosophila melanogaster and D. obscura group was determined. Intra‐ and interspecific sequence comparisons allowed discrimination of mutations affecting the repetitive microsatellite from those affecting the flanking regions. The hypotheses that slippage needs a minimum number of repeats in order to become efficient with respect to microsatellite variability, and of an increased mutation rate with increased length of the microsatellite are supported by the results of our analyses. There is in particular at the interrupted complex microsatellite locus BICOID in the species of the D. obscura group, extensive variation in the flanking regions in addition to length and sequence variation of the repetitive microsatellite. The allelic variation at this locus can hardly be explained by slippage alone. Estimates of microsatellite variability by fragment length analyses will pick up only a minor fraction of allelic variation at such loci, and conclusions that are based on the stepwise mutation model will not hold.     

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.      

Multiple levels of single-strand slippage at cetacean tri- and tetranucleotide repeat microsatellite loci.

Between three and six tri- and tetranucleotide repeat microsatellite loci were analyzed in 3720 samples collected from four different species of baleen whales. Ten of the 18 species/locus combinations had imperfect allele arrays, i.e., some alleles differed in length by other than simple integer multiples of the basic repeat length. The estimate of the average number of alleles and heterozygosity was higher at loci with imperfect allele arrays relative to those with perfect allele arrays. Nucleotide sequences of 23 different alleles at one tetranucleotide repeat microsatellite locus in fin whales, Balaenoptera physalus, and humpback whales, Megaptera novaeangliae, revealed sequence changes including perfect repeats only, multiple repeats, and partial repeats. The relative rate of the latter two categories of mutation was estimated at 0.024 of the mutation rate involving perfect repeats only. It is hypothesized that single-strand slippage of partial repeats may provide a mechanism for counteracting the continuous expansion of microsatellite loci, which is the logical consequence of recent reports demonstrating directional mutations. Partial-repeat mutations introduce imperfections in the repeat array, which subsequently could reduce the rate of single-strand slippage. Limited computer simulations confirmed this predicted effect of partial-repeat mutations.     

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.     

Determining microsatellite genotyping reliability and mutation detection ability: an approach using small-pool PCR from sperm DNA

Microsatellite genotyping from trace DNA is now common in fields as diverse as medicine, forensics and wildlife genetics. Conversely, small-pool PCR (SP-PCR) has been used to investigate microsatellite mutation mechanisms in human DNA, but has had only limited application to non-human species. Trace DNA and SP-PCR studies share many challenges, including problems associated with allelic drop-out, false alleles and other PCR artefacts, and the need to reliably identify genuine alleles and/or mutations. We provide a framework for the validation of such studies without a multiple tube approach and demonstrate the utility of that approach with an analysis of microsatellite mutations in the tammar wallaby (Macropus eugenii). Specifically, we amplified three autosomal microsatellites from somatic DNA to characterise efficiency and reliability of PCR from low-template DNA. Reconstruction experiments determined our ability to discriminate mutations from parental alleles. We then developed rules to guide data interpretation. We estimated mutation rates in sperm DNA to range from 1.5 × 10⁻² to 2.2 × 10⁻³ mutations per locus per generation. Large multi-step mutations were observed, providing evidence for complex mutation processes at microsatellites and potentially violating key assumptions in the stepwise mutation model. Our data demonstrate the necessity of actively searching for large mutation events when investigating microsatellite evolution and highlight the need for a thorough understanding of microsatellite amplification characteristics before embarking on SP-PCR or trace DNA studies.     

Reticulate evolution of the Daphnia pulex complex as revealed by nuclear markers

The study of species complexes is of particular interest to understand how evolutionary young species maintain genomic integrity. The Daphnia pulex complex has been intensively studied as it includes species that dominate freshwater environments in the Northern hemisphere and as it is the sole North American complex that shows transitions to obligate parthenogenesis. Past studies using mitochondrial markers have revealed the presence of 10 distinct lineages in the complex. This study is the first to examine genetic relationships among seven species of the complex at nuclear markers (nine microsatellite loci and one protein-coding gene). Clones belonging to the seven species of the Daphnia pulex complex were characterized at the mitochondrial NADH dehydrogenase (ND5) gene and at the Lactate dehydrogenase (LDH) locus. K-means, principal coordinate analyses and phylogenetic network analyses on the microsatellite data all separated European D. pulicaria, D. tenebrosa, North American D. pulex, D. pulicaria and their hybrids into distinct clusters. The hybrid cluster was composed of diploid and polyploid hybrids with D. pulex mitochondria and some clones with D. pulicaria mitochondria. By contrast, the phylogeny of the D. pulex complex using Rab4 was not well resolved but still showed clusters consisting mostly of D. pulex alleles and others of D. pulicaria alleles. Incomplete lineage sorting and hybridization may obscure genetic relationships at this locus. This study shows that hybridization and introgression have played an important role in the evolution of this complex.     

Genome size evolution and polyploidy in the Daphnia pulex complex (Cladocera: Daphniidae)

Genome size was estimated in 49 clones of the Daphnia pulex complex from temperate and subarctic locations using flow cytometry and microsatellite DNA analyses. Significant genome size differences were found in diploid species belonging to the two genetically distinct groups (the pulicaria and the tenebrosa groups), with clones from the tenebrosa group having genome sizes 22% larger than those in the pulicaria group. Combined flow cytometry and microsatellite DNA analyses revealed that nearly all polyploid clones in the D. pulex complex are triploid and not tetraploid, as was previously suggested. Sequencing analyses of the ND5 gene to position clones in their respective clades within the D. pulex complex have uncovered three triploid clones of Daphnia middendorffiana with a D. pulex maternal parent. This result was unexpected because Daphnia pulicaria has always been identified as the maternal parent of these hybrid polyploid clones. Triploid clones likely owe their origins to interactions between sexual and asexual populations. Further interactions in the tenebrosa group have generated tetraploid clones but these events have been rare.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 97 , 68–79.     

Inference on microsatellite mutation processes in the invasive mite, Varroa destructor, using reversible jump Markov chain Monte Carlo

Varroa destructor is a parasitic mite of the Eastern honeybee Apis cerana. Fifty years ago, two distinct evolutionary lineages (Korean and Japanese) invaded the Western honeybee Apis mellifera. This haplo-diploid parasite species reproduces mainly through brother-sister matings, a system which largely favors the fixation of new mutations. In a worldwide sample of 225 individuals from 21 locations collected on Western honeybees and analyzed at 19 microsatellite loci, a series of de novo mutations was observed. Using historical data concerning the invasion, this original biological system has been exploited to compare three mutation models with allele size constraints for microsatellite markers: stepwise (SMM) and generalized (GSM) mutation models, and a model with mutation rate increasing exponentially with microsatellite length (ESM). Posterior probabilities of the three models have been estimated for each locus individually using reversible jump Markov Chain Monte Carlo. The relative support of each model varies widely among loci, but the GSM is the only model that always receives at least 9% support, whatever the locus. The analysis also provides robust estimates of mutation parameters for each locus and of the divergence time of the two invasive lineages (67,000 generations with a 90% credibility interval of 35,000-174,000). With an average of 10 generations per year, this divergence time fits with the last post-glacial Korea-Japan land separation.     

Factors affecting germline mutations in a hypervariable microsatellite: a comparative analysis of six species of swallows (Aves: Hirundinidae)

Microsatellites mutate frequently by replication slippage. Empirical evidence shows that the probability of such slippage mutations may increase with the length of the repeat region as well as exposure to environmental mutagens, but the mutation rate can also differ between the male and female germline. It has been hypothesized that more intense sexual selection or sperm competition can also lead to elevated mutation rates, but the empirical evidence is inconclusive. Here, we analyzed the occurrence of germline slippage mutations in the hypervariable pentanucleotide microsatellite locus HrU10 across six species of swallow (Aves: Hirundinidae). These species exhibit marked differences in the length range of the microsatellite, as well as differences in the intensity of sperm competition. We found a strong effect of microsatellite length on the probability of mutation, but no residual effect of species or their level of sperm competition when the length effect was accounted for. Neither could we detect any difference in mutation rate between tree swallows (Tachycineta bicolor) breeding in Hamilton Harbour, Ontario, an industrial site with previous documentation of elevated mutation rates for minisatellite DNA, and a rural reference population. However, our cross-species analysis revealed two significant patterns of sex differences in HrU10 germline mutations: (1) mutations in longer alleles occurred typically in the male germline, those in shorter alleles in the female germline, and (2) male germline mutations were more often expansions than contractions, whereas no directional bias was evident in the female germline. These results indicate some fundamental differences in male and female gametogenesis affecting the probability of slippage mutations. Our study also reflects the value of a comparative, multi-species approach for locus-specific mutation analyses, through which a wider range of influential factors can be assessed than in single-species studies.     

Evidence for complex mutations at microsatellite loci in Drosophila.

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

鲤鱼微卫星突变速率和模式(英文)

利用150个微卫星分子标记在F1代家系的基因型分析过程中,共有27600个等位基因从亲本向子代传递,其中在5个微卫星座位上检测到6个突变的等位基因。对突变的等位基因数目进行统计分析后得出:鲤鱼平均每个世代每个微卫星座位的突变速率为2.53×10-4。在发现突变的5个位点中,经测序发现,突变序列中插入1个以上的重复单元就导致了突变的发生。这些突变表明,鲤鱼的微卫星突变没有遵循严格的渐变突变模型(stepwise mutation model,SMM)。该文关于鲤鱼微卫星突变速率和模式的研究将会对统计鲤鱼有效群体的统计提供有效参数。     

Body temperature predicts maximum microsatellite length in mammals

A long-standing mystery in genome evolution is why short tandem repeats vary so much in length and frequency. Here, we test the hypothesis that body temperature acts to influence the rate and nature of slippage-based mutations. Using the data from both 28 species where genome sequencing is advanced and 76 species from which marker loci have been published, we show that in mammals, maximum repeat number is inversely correlated with body temperature, with warmer-blooded species having shorter 'long' microsatellites. Our results support a model of microsatellite evolution in which maximum length is limited by a temperature-dependent stability threshold.     

Microsatellite length differences between humans and chimpanzees at autosomal Loci are not found at equivalent haploid Y chromosomal Loci

When homologous microsatellites are compared between species, significant differences in mean length are often noted. A dominant cause of these length differences is ascertainment bias due to selection for maximum repeat number and repeat purity when the markers are being developed. However, even after ascertainment bias has been allowed for through reciprocal comparisons, significant length differences remain, suggesting that the average microsatellite mutation rate differs between species. Two classes of mechanism have been proposed: rapid evolution of enzymes involved in the generation and repair of slippage products (enzyme evolution model) and heterozygote instability, whereby interchromosomal events at heterozygous sites offer extra opportunities for mutations to occur (heterozygote instability model). To examine which of these hypotheses is most likely, we compared ascertainment bias and species length differences between humans and chimpanzees in autosomal and Y chromosomal microsatellites. We find that levels of ascertainment bias are indistinguishable, but that interspecies length differences are significantly greater for autosomal loci compared with haploid Y chromosomal loci. Such a pattern is consistent with predictions from the heterozygote instability model and is not expected under models of microsatellite evolution that do not include interchromosomal events such as the enzyme evolution model.     

Mutation rates, spectra and hotspots in mismatch repair-deficient Caenorhabditis elegans

Although it is clear that postreplicative DNA mismatch repair (MMR) plays a critical role in maintaining genomic stability in nearly all forms of life surveyed, much remains to be understood about the genome-wide impact of MMR on spontaneous mutation processes and the extent to which MMR-deficient mutation patterns vary among species. We analyzed spontaneous mutation processes across multiple genomic regions using two sets of mismatch repair-deficient (msh-2 and msh-6) Caenorhabditis elegans mutation-accumulation (MA) lines and compared our observations to mutation spectra in a set of wild-type (WT), repair-proficient C. elegans MA lines. Across most sequences surveyed in the MMR-deficient MA lines, mutation rates were approximately 100-fold higher than rates in the WT MA lines, although homopolymeric nucleotide-run (HP) loci composed of A:T base pairs mutated at an approximately 500-fold greater rate. In contrast to yeast and humans where mutation spectra vary substantially with respect to different specific MMR-deficient genotypes, mutation rates and patterns were overall highly similar between the msh-2 and msh-6 C. elegans MA lines. This, along with the apparent absence of a Saccharomyces cerevisiae MSH3 ortholog in the C. elegans genome, suggests that C. elegans MMR surveillance is carried out by a single Msh-2/Msh-6 heterodimer.     

Empirical evaluation reveals best fit of a logistic mutation model for human Y-chromosomal microsatellites

The rate of microsatellite mutation is dependent upon both the allele length and the repeat motif, but the exact nature of this relationship is still unknown. We analyzed data on the inheritance of human Y-chromosomal microsatellites in father-son duos, taken from 24 published reports and comprising 15,285 directly observable meioses. At the six microsatellites analyzed (DYS19, DYS389I, DYS390, DYS391, DYS392, and DYS393), a total of 162 mutations were observed. For each locus, we employed a maximum-likelihood approach to evaluate one of several single-step mutation models on the basis of the data. For five of the six loci considered, a novel logistic mutation model was found to provide the best fit according to Akaike's information criterion. This implies that the mutation probability at the loci increases (nonlinearly) with allele length at a rate that differs between upward and downward mutations. For DYS392, the best fit was provided by a linear model in which upward and downward mutation probabilities increase equally with allele length. This is the first study to empirically compare different microsatellite mutation models in a locus-specific fashion.     

Microsatellite mutations during the polymerase chain reaction: mean field approximations and their applications

We develop a novel mathematical model for microsatellite mutations during polymerase chain reaction (PCR). Based on the model, we study the first- and second-order moments of the number of repeat units in a randomly chosen molecule after n PCR cycles and their corresponding mean field approximations. We give upper bounds for the approximation errors and show that the approximation errors are small when the mutation rate is low. Based on the theoretical results, we develop a moment estimation method to estimate the mutation rate per-repeat-unit per PCR cycle and the probability of expansion when mutations occur. Simulation studies show that the moment estimation method can accurately recover the true mutation rate and probability of expansion. Finally, the method is applied to experimental data from single-molecule PCR experiments.