Estimates Of Mutation Rates Research Articles

Modeling age-dependent radiation-induced second cancer risks and estimation of mutation rate: an evolutionary approach.

Although the survival rate of cancer patients has significantly increased due to advances in anti-cancer therapeutics, one of the major side effects of these therapies, particularly radiotherapy, is the potential manifestation of radiation-induced secondary malignancies. In this work, a novel evolutionary stochastic model is introduced that couples short-term formalism (during radiotherapy) and long-term formalism (post-treatment). This framework is used to estimate the risks of second cancer as a function of spontaneous background and radiation-induced mutation rates of normal and pre-malignant cells. By fitting the model to available clinical data for spontaneous background risk together with data of Hodgkin's lymphoma survivors (for various organs), the second cancer mutation rate is estimated. The model predicts a significant increase in mutation rate for some cancer types, which may be a sign of genomic instability. Finally, it is shown that the model results are in agreement with the measured results for excess relative risk (ERR) as a function of exposure age and that the model predicts a negative correlation of ERR with increase in attained age. This novel approach can be used to analyze several radiotherapy protocols in current clinical practice and to forecast the second cancer risks over time for individual patients.

Estimation of the spontaneous mutation rate in Heliconius melpomene.

We estimated the spontaneous mutation rate in Heliconius melpomene by genome sequencing of a pair of parents and 30 of their offspring, based on the ratio of number of de novo heterozygotes to the number of callable site-individuals. We detected nine new mutations, each one affecting a single site in a single offspring. This yields an estimated mutation rate of 2.9 × 10−9 (95% confidence interval, 1.3 × 10−9–5.5 × 10−9), which is similar to recent estimates in Drosophila melanogaster, the only other insect species in which the mutation rate has been directly estimated. We infer that recent effective population size of H. melpomene is about 2 million, a substantially lower value than its census size, suggesting a role for natural selection reducing diversity. We estimate that H. melpomene diverged from its Müllerian comimic H. erato about 6 Ma, a somewhat later date than estimates based on a local molecular clock.

Haplotype and mutation analysis for newly suggested Y-STRs in Korean father-son pairs.

In this study, 363 Korean father-son haplotype transfers in 351 families were analyzed using an in-house multiplex PCR system for 14 Y-STRs (DYS385a/b, DYF387S1, DYS391, DYS449, DYS460, DYS481, DYS518, DYS533, DYS549, DYS570, DYS576, DYS627 and DYS643), that included 11 loci newly added to the PowerPlex Y23 system or the Yfiler Plus system. The Y-STRs showed gene diversity values ranging from 0.2499 to 0.9612; the multicopy Y-STR loci DYS385 and DYF387S1 had high gene diversity of 0.9612 and 0.9457, respectively. In addition, DYF387S1, which has two copies, showed three alleles in seven individuals, and micro-variant alleles were observed in 14 individuals at four loci (DYS448, DYS518, DYS570 and DYS627). Among 351 haplotypes for the 11 newly added Y-STRs, 350 different haplotypes were observed, with an overall haplotype diversity of 0.9999 and discrimination capacity of 99.72%. In 363 haplotype transfers from 351 pedigrees, 29 single-step mutations were observed at 11 Y-STRs. Locus-specific mutation rate estimates varied from 0.0 to 1.93×10(-2), with an average estimated mutation rate of 6.66×10(-3). Two father-son pairs had mutations at two different loci in 11 Y-STRs. The number of pairs with mutations at multiple loci increased to five when the mutation event was investigated for haplotype transfer at 28 Y-STRs including 17 Yfiler loci and 11 Y-STRs examined in this study: four father-son pairs had mutations at two loci, and one pair had mutations at three loci. Overall, mutations were frequently observed at DYS449, DYS576 and DYS627 loci, which are known to be rapidly mutating Y-STRs. Mutation rate estimates at most loci were not significantly different from rates in other populations, but estimates for DYF387S1, DYS518 and DYS570 were considerably lower in the Korean population than in other populations.

Evolutionary ecotoxicology of perfluoralkyl substances (PFASs) inferred from multigenerational exposure: A case study with Chironomus riparius (Diptera, Chironomidae)

A multigeneration toxicity test on Chironomus riparius was performed with the aim of investigating the evolutionary consequences of exposure to perfluoralkyl substances (perfluorooctane sulfonic acid, PFOS; perfluorooctanoic acid, PFOA; perfluorobutane sulfonate, PFBS). Six-hundred larvae were bred per treatment and per generation until emergence and egg deposition under a nominal concentration of 10μg/L of contaminants. Newborn larvae were used to start the next generation. Evolution of genetic variability was evaluated along a total of 10 consecutive generations based on 5 microsatellite loci. Analysis of life-history traits (survival, sex ratio and reproduction) was also carried out. Rapid genetic variability reduction was observed in all treatments, including controls, across generations due to the test conditions. Nevertheless, an increased mutation rate determined a stronger conservation of genetic variability in PFOS and, at minor extent, in PFBS exposed populations compared to controls. No significant effects were induced by exposure to PFOA. Direct mutagenicity or induced stress conditions may be at the base of increased mutation rate, indicating the potential risk of mutational load caused by exposure to PFOS and PFBS.The test provided the opportunity to evaluate the use of approximate Bayesian computation (ABC) and coalescent approaches in evolutionary ecotoxicology. A weak performance was evidenced for ABC, either in terms of bias or dispersion of effective population sizes and of estimates of mutation rate. On the contrary, coalescent simulations proved the sensitivity of traditional genetic endpoints (i.e. heterozygosity and number of alleles) to the alteration of mutation rate, but not to erosion of genetic effective size.

Estimation of the RNU2 macrosatellite mutation rate by BRCA1 mutation tracing.

Large tandem repeat sequences have been poorly investigated as severe technical limitations and their frequent absence from the genome reference hinder their analysis. Extensive allelotyping of this class of variation has not been possible until now and their mutational dynamics are still poorly known. In order to estimate the mutation rate of a macrosatellite, we analysed in detail the RNU2 locus, which displays at least 50 different alleles containing 5-82 copies of a 6.1 kb repeat unit. Mining data from the 1000 Genomes Project allowed us to precisely estimate copy numbers of the RNU2 repeat unit using read depth of coverage. This further revealed significantly different mean values in various recent modern human populations, favoring a scenario of fast evolution of this locus. Its proximity to a disease gene with numerous founder mutations, BRCA1, within the same linkage disequilibrium block, offered the unique opportunity to trace RNU2 arrays over a large timescale. Analysis of the transmission of RNU2 arrays associated with one ‘private’ mutation in an extended kindred and four founder mutations in multiple kindreds gave an estimation by maximum likelihood of 5 × 10−3 mutations per generation, which is close to that of microsatellites.

Unbiased estimation of mutation rates under fluctuating final counts.

Estimation methods for mutation rates (or probabilities) in Luria-Delbrück fluctuation analysis usually assume that the final number of cells remains constant from one culture to another. We show that this leads to systematically underestimate the mutation rate. Two levels of information on final numbers are considered: either the coefficient of variation has been independently estimated, or the final number of cells in each culture is known. In both cases, unbiased estimation methods are proposed. Their statistical properties are assessed both theoretically and through Monte-Carlo simulation. As an application, the data from two well known fluctuation analysis studies on Mycobacterium tuberculosis are reexamined.

A mutation burst during the acute phase of Helicobacter pylori infection in humans and rhesus macaques.

The evolution rate and genetic changes that occur during chronic infection with Helicobacter pylori have been analysed, but little is known about the genomic changes during the initial, acute bacterial infection phase. Here we analyse the rate and pattern of genome evolution in H. pylori from the genomes of two input strains isolated from human volunteers with asymptomatic infection, and the genomes of two output strains collected 20 and 44 days after re-infection. Similarly, we analyse genome evolution in bacteria from the genome sequences of input and output strains sequentially taken after experimental infection of a rhesus macaque. The estimated mutation rate reveals a mutation burst during the acute infection phase that is over 10 times faster than the mutation rate during chronic infection, and orders of magnitude faster than mutation rates in any other bacteria. The elevated frequency of mutations in outer membrane protein genes suggests that the mutation burst facilitates rapid host adaptation of the bacteria.

Effects of new mutations on fitness: insights from models and data.

The rates and properties of new mutations affecting fitness have implications for a number of outstanding questions in evolutionary biology. Obtaining estimates of mutation rates and effects has historically been challenging, and little theory has been available for predicting the distribution of fitness effects (DFE); however, there have been recent advances on both fronts. Extreme-value theory predicts the DFE of beneficial mutations in well-adapted populations, while phenotypic fitness landscape models make predictions for the DFE of all mutations as a function of the initial level of adaptation and the strength of stabilizing selection on traits underlying fitness. Direct experimental evidence confirms predictions on the DFE of beneficial mutations and favors distributions that are roughly exponential but bounded on the right. A growing number of studies infer the DFE using genomic patterns of polymorphism and divergence, recovering a wide range of DFE. Future work should be aimed at identifying factors driving the observed variation in the DFE. We emphasize the need for further theory explicitly incorporating the effects of partial pleiotropy and heterogeneity in the environment on the expected DFE.

Precise estimates of mutation rate and spectrum in yeast.

Mutation is the ultimate source of genetic variation. The most direct and unbiased method of studying spontaneous mutations is via mutation accumulation (MA) lines. Until recently, MA experiments were limited by the cost of sequencing and thus provided us with small numbers of mutational events and therefore imprecise estimates of rates and patterns of mutation. We used whole-genome sequencing to identify nearly 1,000 spontaneous mutation events accumulated over ∼311,000 generations in 145 diploid MA lines of the budding yeast Saccharomyces cerevisiae. MA experiments are usually assumed to have negligible levels of selection, but even mild selection will remove strongly deleterious events. We take advantage of such patterns of selection and show that mutation classes such as indels and aneuploidies (especially monosomies) are proportionately much more likely to contribute mutations of large effect. We also provide conservative estimates of indel, aneuploidy, environment-dependent dominant lethal, and recessive lethal mutation rates. To our knowledge, for the first time in yeast MA data, we identified a sufficiently large number of single-nucleotide mutations to measure context-dependent mutation rates and were able to (i) confirm strong AT bias of mutation in yeast driven by high rate of mutations from C/G to T/A and (ii) detect a higher rate of mutation at C/G nucleotides in two specific contexts consistent with cytosine methylation in S. cerevisiae.

Estimation of the in vivo recombination rate for a plant RNA virus

Phylogenomic evidence suggested that recombination is an important evolutionary force for potyviruses, one of the larger families of plant RNA viruses. However, mixed-genotype potyvirus infections are marked by low levels of cellular coinfection, precluding template switching and recombination events between virus genotypes during genomic RNA replication. To reconcile these conflicting observations, we evaluated the in vivo recombination rate (rg) of Tobacco etch virus (TEV; genus Potyvirus, family Potyviridae) by coinfecting plants with pairs of genotypes marked with engineered restriction sites as neutral markers. The recombination rate was then estimated using two different approaches: (i) a classical approach that assumed recombination between marked genotypes can occur in the whole virus population, rendering an estimate of rg = 7.762 × 10(-8) recombination events per nucleotide site per generation, and (ii) an alternative method that assumed recombination between marked genotypes can occur only in coinfected cells, rendering a much higher estimate of rg = 3.427 × 10(-5) recombination events per nucleotide site per generation. This last estimate is similar to the TEV mutation rate, suggesting that recombination should be at least as important as point mutation in creating variability. Finally, we compared our mutation and recombination rate estimates to those reported for animal RNA viruses. Our analysis suggested that high recombination rates may be an unavoidable consequence of selection for fast replication at the cost of low fidelity.

Formae speciales of cereal powdery mildew: close or distant relatives?

Powdery mildew is an important disease of cereals, affecting both grain yield and end-use quality. The causal agent of powdery mildew on cereals, Blumeria graminis, has been classified into eight formae speciales (ff.spp.), infecting crops and wild grasses. Advances in research on host specificity and resistance, and on pathogen phylogeny and origins, have brought aspects of the subspecific classification system of B. graminis into ff.spp. into question, because it is based on adaptation to certain hosts rather than strict host specialization. Cereals therefore cannot be considered as typical non-hosts to non-adapted ff.spp. We introduce the term 'non-adapted resistance' of cereals to inappropriate ff.spp. of B. graminis, which involves both pathogen-associated molecular pattern-triggered immunity (PTI) and effector-triggered immunity (ETI). There is no clear distinction between the mechanisms of resistance to adapted and non-adapted ff.spp. Molecular evolutionary data suggest that the taxonomic grouping of B. graminis into different ff.spp. is not consistent with the phylogeny of the fungus. Imprecise estimates of mutation rates and the lack of genetic variation in introduced populations may explain the uncertainty with regard to divergence times, in the Miocene or Holocene epochs, of ff.spp. of B. graminis which infect cereal crop species. We propose that most evidence favours divergence in the Holocene, during the course of early agriculture. We also propose that the forma specialis concept should be retained for B. graminis pathogenic on cultivated cereals to include clades of the fungus which are strongly specialized to these hosts, i.e. ff.spp. hordei, secalis and tritici, as well as avenae from cultivated A. sativa, and that the forma specialis concept should no longer be applied to B. graminis from most wild grasses.

Estimation of the spontaneous mutation rate per nucleotide site in a Drosophila melanogaster full-sib family.

We employed deep genome sequencing of two parents and 12 of their offspring to estimate the mutation rate per site per generation in a full-sib family of Drosophila melanogaster recently sampled from a natural population. Sites that were homozygous for the same allele in the parents and heterozygous in one or more offspring were categorized as candidate mutations and subjected to detailed analysis. In 1.23 × 10(9) callable sites from 12 individuals, we confirmed six single nucleotide mutations. We estimated the false negative rate in the experiment by generating synthetic mutations using the empirical distributions of numbers of nonreference bases at heterozygous sites in the offspring. The proportion of synthetic mutations at callable sites that we failed to detect was <1%, implying that the false negative rate was extremely low. Our estimate of the point mutation rate is 2.8 × 10(-9) (95% confidence interval = 1.0 × 10(-9) - 6.1 × 10(-9)) per site per generation, which is at the low end of the range of previous estimates, and suggests an effective population size for the species of ∼1.4 × 10(6). At one site, point mutations were present in two individuals, indicating that there had been a premeiotic mutation cluster, although surprisingly one individual had a G→A transition and the other a G→T transversion, possibly associated with error-prone mismatch repair. We also detected three short deletion mutations and no insertions, giving a deletion mutation rate of 1.2 × 10(-9) (95% confidence interval = 0.7 × 10(-9) - 11 × 10(-9)).

Molecular phylogenetic analysis of the genus Abies (Pinaceae) based on the nucleotide sequence of chloroplast DNA

A phylogenetic study of firs (Abies Mill.) was conducted using nucleotide sequences of several chloroplast DNA regions with a total length of 5580 bp. The analysis included 37 taxa, which represented the main evolutionary lineages of the genus, and Keteleeria daviana. According to phylogenetic reconstruction the Abies species were subdivided into six main groups, generally corresponding to their geographic distribution. The phylogenetic tree had three basal clades. All of these clades contained American species, and only one of them contained Eurasian species. The divergence time calibrations, based on paleobotanical data and the chloroplast DNA mutation rate estimates in Pinaceae, produced similar results..The age of diversification among the clades of the present-day Abies was estimated as the end of the Oligocene-beginning of Miocene. The age of the separation of Mediterranean firs from the Asian-North American branch corresponds to the Miocene. The age of diversification within the young groups of Mediterranean, Asian, and boreal American firs (A. lasiocarpa, A. balsamea, A. fraseri) was estimated as the Pliocene-Pleistocene. Based on the phylogenetic reconstruction obtained, the most plausible biogeographic scenarios were suggested. It is noted that the existing systematic classification of the genus Abies strongly contradicts with phylogenetic reconstruction and requires revision.

Fluctuation analysis: can estimates be trusted?

The estimation of mutation rates and relative fitnesses in fluctuation analysis is based on the unrealistic hypothesis that the single-cell times to division are exponentially distributed. Using the classical Luria-Delbrück distribution outside its modelling hypotheses induces an important bias on the estimation of the relative fitness. The model is extended here to any division time distribution. Mutant counts follow a generalization of the Luria-Delbrück distribution, which depends on the mean number of mutations, the relative fitness of normal cells compared to mutants, and the division time distribution of mutant cells. Empirical probability generating function techniques yield precise estimates both of the mean number of mutations and the relative fitness of normal cells compared to mutants. In the case where no information is available on the division time distribution, it is shown that the estimation procedure using constant division times yields more reliable results. Numerical results both on observed and simulated data are reported.

Modelling and inference reveal nonlinear length-dependent suppression of somatic instability for small disease associated alleles in myotonic dystrophy type 1 and Huntington disease

More than 20 human genetic diseases are associated with inheriting an unstable expanded DNA simple sequence tandem repeat, for example, CTG (cytosine-thymine-guanine) repeats in myotonic dystrophy type 1 (DM1) and CAG (cytosine-adenine-guanine) repeats in Huntington disease (HD). These sequences mutate by changing the number of repeats not just between generations, but also during the lifetime of affected individuals. Levels of somatic instability contribute to disease onset and progression but as changes are tissue-specific, age- and repeat length-dependent, interpretation of the level of somatic instability in an individual is confounded by these considerations. Mathematical models, fitted to CTG repeat length distributions derived from blood DNA, from a large cohort of DM1-affected or at risk individuals, have recently been used to quantify inherited repeat lengths and mutation rates. Taking into account age, the estimated mutation rates are lower than predicted among individuals with small alleles (inherited repeat lengths less than 100 CTGs), suggesting that these rates may be suppressed at the lower end of the disease-causing range. In this study, we propose that a length-specific effect operates within this range and tested this hypothesis using a model comparison approach. To calibrate the extended model, we used data derived from blood DNA from DM1 individuals and, for the first time, buccal DNA from HD individuals. In a novel application of this extended model, we identified individuals whose effective repeat length, with regards to somatic instability, is less than their actual repeat length. A plausible explanation for this distinction is that the expanded repeat tract is compromised by interruptions or other unusual features. We quantified effective length for a large cohort of DM1 individuals and showed that effective length better predicts age of onset than inherited repeat length, thus improving the genotype-phenotype correlation. Under the extended model, we removed some of the bias in mutation rates making them less length-dependent. Consequently, rates adjusted in this way will be better suited as quantitative traits to investigate cis- or trans-acting modifiers of somatic mosaicism, disease onset and progression.

Mutation and mutation rates at Y chromosome specific Short Tandem Repeat Polymorphisms (STRs): A reappraisal

Mutation is a topic of intense research and raises important problems in forensics. Since the markers of choice in current forensic genetics analyses are microsatellites or Short Tandem Repeat Polymorphisms (STRs), mutation is sufficiently common to cause difficulties in evaluating DNA evidence in a significant proportion of cases but at the same time rare enough to turn the estimation of the corresponding probability of occurrence into a hard task. We address these issues using the simplest model of transmission: the Y chromosome specific STRs. Within this model, and under an explicit set of definitions and involved assumptions, we developed the theoretical framework required for the study of allelic transitions in gametogenesis, identifying the required parameters and associated probabilities and finally we discuss the estimation of these parameters and their application in forensics. We conclude that (i) for forensic casework the relevant parameter for incorporation in a likelihood ratio is biallelic specific (i.e. the mutation rate estimate corresponds to the probability of the specific allelic transition observed) and (ii) for these estimates as well as in order to provide data for testing mutation models the absolute frequency of mutated and non-mutated transmissions per allele, along with the description of the observed mutations should be reported.

Using HSV-1 Genome Phylogenetics to Track Past Human Migrations

We compared 31 complete and nearly complete globally derived HSV-1 genomic sequences using HSV-2 HG52 as an outgroup to investigate their phylogenetic relationships and look for evidence of recombination. The sequences were retrieved from NCBI and were then aligned using Clustal W. The generation of a maximum likelihood tree resulted in a six clade structure that corresponded with the timing and routes of past human migration. The East African derived viruses contained the greatest amount of genetic diversity and formed four of the six clades. The East Asian and European/North American derived viruses formed separate clades. HSV-1 strains E07, E22 and E03 were highly divergent and may each represent an individual clade. Possible recombination was analyzed by partitioning the alignment into 5 kb segments, performing individual phylogenetic analysis on each partition and generating a.phylogenetic network from the results. However most evidence for recombination spread at the base of the tree suggesting that recombination did not significantly disrupt the clade structure. Examination of previous estimates of HSV-1 mutation rates in conjunction with the phylogenetic data presented here, suggests that the substitution rate for HSV-1 is approximately 1.38 × 10(-7) subs/site/year. In conclusion, this study expands the previously described HSV-1 three clade phylogenetic structures to a minimum of six and shows that the clade structure also mirrors global human migrations. Given that HSV-1 has co-evolved with its host, sequencing HSV-1 isolated from various populations could serve as a surrogate biomarker to study human population structure and migration patterns.

Factors Influencing Ascertainment Bias of Microsatellite Allele Sizes: Impact on Estimates of Mutation Rates

Microsatellite loci play an important role as markers for identification, disease gene mapping, and evolutionary studies. Mutation rate, which is of fundamental importance, can be obtained from interspecies comparisons, which, however, are subject to ascertainment bias. This bias arises, for example, when a locus is selected on the basis of its large allele size in one species (cognate species 1), in which it is first discovered. This bias is reflected in average allele length in any noncognate species 2 being smaller than that in species 1. This phenomenon was observed in various pairs of species, including comparisons of allele sizes in human and chimpanzee. Various mechanisms were proposed to explain observed differences in mean allele lengths between two species. Here, we examine the framework of a single-step asymmetric and unrestricted stepwise mutation model with genetic drift. Analysis is based on coalescent theory. Analytical results are confirmed by simulations using the simuPOP software. The mechanism of ascertainment bias in this model is a tighter correlation of allele sizes within a cognate species 1 than of allele sizes in two different species 1 and 2. We present computations of the expected average allele size difference, given the mutation rate, population sizes of species 1 and 2, time of separation of species 1 and 2, and the age of the allele. We show that when the past demographic histories of the cognate and noncognate taxa are different, the rate and directionality of mutations affect the allele sizes in the two taxa differently from the simple effect of ascertainment bias. This effect may exaggerate or reverse the effect of difference in mutation rates. We reanalyze literature data, which indicate that despite the bias, the microsatellite mutation rate estimate in the ancestral population is consistently greater than that in either human or chimpanzee and the mutation rate estimate in human exceeds or equals that in chimpanzee with the rate of allele length expansion in human being greater than that in chimpanzee. We also demonstrate that population bottlenecks and expansions in the recent human history have little impact on our conclusions.

Accumulation of spontaneous mutations in the ciliate Tetrahymena thermophila.

Knowledge of the rate and fitness effects of mutations is essential for understanding the process of evolution. Mutations are inherently difficult to study because they are rare and are frequently eliminated by natural selection. In the ciliate Tetrahymena thermophila, mutations can accumulate in the germline genome without being exposed to selection. We have conducted a mutation accumulation (MA) experiment in this species. Assuming that all mutations are deleterious and have the same effect, we estimate that the deleterious mutation rate per haploid germline genome per generation is U=0.0047 (95% credible interval: 0.0015, 0.0125), and that germline mutations decrease fitness by s=11% when expressed in a homozygous state (95% CI: 4.4%, 27%). We also estimate that deleterious mutations are partially recessive on average (h=0.26; 95% CI: -0.022, 0.62) and that the rate of lethal mutations is <10% of the deleterious mutation rate. Comparisons between the observed evolutionary responses in the germline and somatic genomes and the results from individual-based simulations of MA suggest that the two genomes have similar mutational parameters. These are the first estimates of the deleterious mutation rate and fitness effects from the eukaryotic supergroup Chromalveolata and are within the range of those of other eukaryotes.

Why to account for finite sites in population genetic studies and how to do this with Jaatha 2.0

With the advent of next-generation sequencing technologies, large data sets of several thousand loci from multiple conspecific individuals are available. Such data sets should make it possible to obtain accurate estimates of population genetic parameters, even for complex models of population history. In the analyses of large data sets, it is difficult to consider finite-sites mutation models (FSMs). Here, we use extensive simulations to demonstrate that the inclusion of FSMs is necessary to avoid severe biases in the estimation of the population mutation rate θ, population divergence times, and migration rates. We present a new version of Jaatha, an efficient composite-likelihood method for estimating demographic parameters from population genetic data and evaluate the usefulness of Jaatha in two biological examples. For the first application, we infer the speciation process of two wild tomato species, Solanum chilense and Solanum peruvianum. In our second application example, we demonstrate that Jaatha is readily applicable to NGS data by analyzing genome-wide data from two southern European populations of Arabidopsis thaliana. Jaatha is now freely available as an R package from the Comprehensive R Archive Network (CRAN).