Estimating the distribution of fitness effects from DNA sequence data: implications for the molecular clock.
Ontology highlight
ABSTRACT:
SUBMITTER: Piganeau G
PROVIDER: S-EPMC193562 | biostudies-literature | 2003 Sep
REPOSITORIES: biostudies-literature
Publications
Estimating the distribution of fitness effects from DNA sequence data: implications for the molecular clock.
Piganeau Gwenaël G Eyre-Walker Adam A
Proceedings of the National Academy of Sciences of the United States of America 20030818 18
We present a method for estimating the distribution of fitness effects of new amino acid mutations when those mutations can be assumed to be slightly advantageous, slightly deleterious, or strongly deleterious. We apply the method to mitochondrial data from several different species. In the majority of the data sets, the shape of the distribution is approximately exponential. Our results provide an estimate of the distribution of fitness effects of weakly selected mutations and provide a possibl ...[more]
Similar Datasets
Project description:The evolution of HIV during acute infection is often considered a neutral process. Recent analysis of sequencing data from this stage of infection, however, showed high levels of shared mutations between independent viral populations. This suggests that selection might play a role in the early stages of HIV infection. We adapted an existing model for random evolution during acute HIV-infection to include selection. Simulations of this model were used to fit a global mutational fitness effects distribution to previously published sequencing data of the env gene of individuals with acute HIV infection. Measures of sharing between viral populations were used as summary statistics to compare the data to the simulations. We confirm that evolution during acute infection is significantly different from neutral. The distribution of mutational fitness effects is best fit by a distribution with a low, but significant fraction of beneficial mutations and a high fraction of deleterious mutations. While most mutations are neutral or deleterious in this model, about 5% of mutations are beneficial. These beneficial mutations will, on average, result in a small but significant increase in fitness. When assuming no epistasis, this indicates that, at the moment of transmission, HIV is near, but not on the fitness peak for early infection.
| S-EPMC6172364 | biostudies-literature
Project description:Accurately dating when the first bilaterally symmetrical animals arose is crucial to our understanding of early animal evolution. The earliest unequivocally bilaterian fossils are approximately 555 million years old. In contrast, molecular-clock analyses calibrated by using the fossil record of vertebrates estimate that vertebrates split from dipterans (Drosophila) approximately 900 million years ago (Ma). Nonetheless, comparative genomic analyses suggest that a significant rate difference exists between vertebrates and dipterans, because the percentage difference between the genomes of mosquito and fly is greater than between fish and mouse, even though the vertebrate divergence is almost twice that of the dipteran. Here we show that the dipteran rate of molecular evolution is similar to other invertebrate taxa (echinoderms and bivalve molluscs) but not to vertebrates, which significantly decreased their rate of molecular evolution with respect to invertebrates. Using a data set consisting of the concatenation of seven different amino acid sequences from 23 ingroup taxa (giving a total of 11 different invertebrate calibration points scattered throughout the bilaterian tree and across the Phanerozoic), we estimate that the last common ancestor of bilaterians arose somewhere between 573 and 656 Ma, depending on the value assigned to the parameter scaling molecular substitution rate heterogeneity. These results are in accord with the known fossil record and support the view that the Cambrian explosion reflects, in part, the diversification of bilaterian phyla.
| S-EPMC404080 | biostudies-literature
Project description:BACKGROUND:Phylogenetic analysis of DNA from modern and ancient samples allows the reconstruction of important demographic and evolutionary processes. A critical component of these analyses is the estimation of evolutionary rates, which can be calibrated using information about the ages of the samples. However, the reliability of these rate estimates can be negatively affected by among-lineage rate variation and non-random sampling. Using a simulation study, we compared the performance of three phylogenetic methods for inferring evolutionary rates from time-structured data sets: regression of root-to-tip distances, least-squares dating, and Bayesian inference. We also applied these three methods to time-structured mitogenomic data sets from six vertebrate species. RESULTS:Our results from 12 simulation scenarios show that the three methods produce reliable estimates when the substitution rate is high, rate variation is low, and samples of similar ages are not all grouped together in the tree (i.e., low phylo-temporal clustering). The interaction of these factors is particularly important for least-squares dating and Bayesian estimation of evolutionary rates. The three estimation methods produced consistent estimates of rates across most of the six mitogenomic data sets, with sequence data from horses being an exception. CONCLUSIONS:We recommend that phylogenetic studies of ancient DNA sequences should use multiple methods of inference and test for the presence of temporal signal, among-lineage rate variation, and phylo-temporal clustering in the data.
| S-EPMC5956955 | biostudies-literature
Project description:The distribution of fitness effects (DFE) encompasses the fraction of deleterious, neutral, and beneficial mutations. It conditions the evolutionary trajectory of populations, as well as the rate of adaptive molecular evolution (α). Inferring DFE and α from patterns of polymorphism, as given through the site frequency spectrum (SFS) and divergence data, has been a longstanding goal of evolutionary genetics. A widespread assumption shared by previous inference methods is that beneficial mutations only contribute negligibly to the polymorphism data. Hence, a DFE comprising only deleterious mutations tends to be estimated from SFS data, and α is then predicted by contrasting the SFS with divergence data from an outgroup. We develop a hierarchical probabilistic framework that extends previous methods to infer DFE and α from polymorphism data alone. We use extensive simulations to examine the performance of our method. While an outgroup is still needed to obtain an unfolded SFS, we show that both a DFE, comprising both deleterious and beneficial mutations, and α can be inferred without using divergence data. We also show that not accounting for the contribution of beneficial mutations to polymorphism data leads to substantially biased estimates of the DFE and α We compare our framework with one of the most widely used inference methods available and apply it on a recently published chimpanzee exome data set.
| S-EPMC5676230 | biostudies-literature
Project description:A common challenge arising from the observation of an evolutionary system over time is to infer the magnitude of selection acting upon a specific genetic variant, or variants, within the population. The inference of selection may be confounded by the effects of genetic drift in a system, leading to the development of inference procedures to account for these effects. However, recent work has suggested that deterministic models of evolution may be effective in capturing the effects of selection even under complex models of demography, suggesting the more general application of deterministic approaches to inference. Responding to this literature, we here note a case in which a deterministic model of evolution may give highly misleading inferences, resulting from the nondeterministic properties of mutation in a finite population. We propose an alternative approach that acts to correct for this error, and which we denote the delay-deterministic model. Applying our model to a simple evolutionary system, we demonstrate its performance in quantifying the extent of selection acting within that system. We further consider the application of our model to sequence data from an evolutionary experiment. We outline scenarios in which our model may produce improved results for the inference of selection, noting that such situations can be easily identified via the use of a regular deterministic model.
| S-EPMC5937181 | biostudies-literature
Project description:Despite the success of highly active antiretroviral therapy (HAART) in the management of human immunodeficiency virus (HIV)-1 infection, virological failure due to drug resistance development remains a major challenge. Resistant mutants display reduced drug susceptibilities, but in the absence of drug, they generally have a lower fitness than the wild type, owing to a mutation-incurred cost. The interaction between these fitness costs and drug resistance dictates the appearance of mutants and influences viral suppression and therapeutic success. Assessing in vivo viral fitness is a challenging task and yet one that has significant clinical relevance. Here, we present a new computational modelling approach for estimating viral fitness that relies on common sparse cross-sectional clinical data by combining statistical approaches to learn drug-specific mutational pathways and resistance factors with viral dynamics models to represent the host-virus interaction and actions of drug mechanistically. We estimate in vivo fitness characteristics of mutant genotypes for two antiretroviral drugs, the reverse transcriptase inhibitor zidovudine (ZDV) and the protease inhibitor indinavir (IDV). Well-known features of HIV-1 fitness landscapes are recovered, both in the absence and presence of drugs. We quantify the complex interplay between fitness costs and resistance by computing selective advantages for different mutants. Our approach extends naturally to multiple drugs and we illustrate this by simulating a dual therapy with ZDV and IDV to assess therapy failure. The combined statistical and dynamical modelling approach may help in dissecting the effects of fitness costs and resistance with the ultimate aim of assisting the choice of salvage therapies after treatment failure.
| S-EPMC4222584 | biostudies-literature
Project description:The distribution of fitness effects (DFE) of new mutations is a key parameter in determining the course of evolution. This fact has motivated extensive efforts to measure the DFE or to predict it from first principles. However, just as the DFE determines the course of evolution, the evolutionary process itself constrains the DFE. Here, we analyze a simple model of genome evolution in a constant environment in which natural selection drives the population toward a dynamic steady state where beneficial and deleterious substitutions balance. The distribution of fitness effects at this steady state is stable under further evolution and provides a natural null expectation for the DFE in a population that has evolved in a constant environment for a long time. We calculate how the shape of the evolutionarily stable DFE depends on the underlying population genetic parameters. We show that, in the absence of epistasis, the ratio of beneficial to deleterious mutations of a given fitness effect obeys a simple relationship independent of population genetic details. Finally, we analyze how the stable DFE changes in the presence of a simple form of diminishing-returns epistasis.
| S-EPMC4423373 | biostudies-literature
Project description:The molecular clock of neutral mutations, which represents linear mutation fixation over generations, is theoretically explained by genetic drift in fitness-steady evolution or hitchhiking in adaptive evolution. The present study is the first experimental demonstration for the molecular clock of neutral mutations in a fitness-increasing evolutionary process. The dynamics of genome mutation fixation in the thermal adaptive evolution of Escherichia coli were evaluated in a prolonged evolution experiment in duplicated lineages. The cells from the continuously fitness-increasing evolutionary process were subjected to genome sequencing and analyzed at both the population and single-colony levels. Although the dynamics of genome mutation fixation were complicated by the combination of the stochastic appearance of adaptive mutations and clonal interference, the mutation fixation in the population was simply linear over generations. Each genome in the population accumulated 1.6 synonymous and 3.1 non-synonymous neutral mutations, on average, by the spontaneous mutation accumulation rate, while only a single genome in the population occasionally acquired an adaptive mutation. The neutral mutations that preexisted on the single genome hitchhiked on the domination of the adaptive mutation. The successive fixation processes of the 128 mutations demonstrated that hitchhiking and not genetic drift were responsible for the coincidence of the spontaneous mutation accumulation rate in the genome with the fixation rate of neutral mutations in the population. The molecular clock of neutral mutations to the fitness-increasing evolution suggests that the numerous neutral mutations observed in molecular phylogenetic trees may not always have been fixed in fitness-steady evolution but in adaptive evolution.
| S-EPMC4503671 | biostudies-literature
Project description:Telomeres play a key role in replicative ageing and undergo age-dependent attrition in vivo. Here, we report a novel method, TelSeq, to measure average telomere length from whole genome or exome shotgun sequence data. In 260 leukocyte samples, we show that TelSeq results correlate with Southern blot measurements of the mean length of terminal restriction fragments (mTRFs) and display age-dependent attrition comparably well as mTRFs.
| S-EPMC4027178 | biostudies-literature
Project description:Observations that rates of molecular evolution vary widely within and among lineages have cast doubts on the existence of a single "molecular clock." Differences in the timing of evolutionary events estimated from genetic and fossil evidence have raised further questions about the accuracy of molecular clocks. Here, we present a model of nucleotide substitution that combines theory on metabolic rate with the now-classic neutral theory of molecular evolution. The model quantitatively predicts rate heterogeneity and may reconcile differences in molecular- and fossil-estimated dates of evolutionary events. Model predictions are supported by extensive data from mitochondrial and nuclear genomes. By accounting for the effects of body size and temperature on metabolic rate, this model explains heterogeneity in rates of nucleotide substitution in different genes, taxa, and thermal environments. This model also suggests that there is indeed a single molecular clock, as originally proposed by Zuckerkandl and Pauling [Zuckerkandl, E. & Pauling, L. (1965) in Evolving Genes and Proteins, eds. Bryson, V. & Vogel, H. J. (Academic, New York), pp. 97-166], but that it "ticks" at a constant substitution rate per unit of mass-specific metabolic energy rather than per unit of time. This model therefore links energy flux and genetic change. More generally, the model suggests that body size and temperature combine to control the overall rate of evolution through their effects on metabolism.
| S-EPMC544068 | biostudies-literature