Dataset Information

Accelerated simulation of evolutionary trajectories in origin-fixation models.

ABSTRACT: We present an accelerated algorithm to forward-simulate origin-fixation models. Our algorithm requires, on average, only about two fitness evaluations per fixed mutation, whereas traditional algorithms require, per one fixed mutation, a number of fitness evaluations of the order of the effective population size, Ne Our accelerated algorithm yields the exact same steady state as the original algorithm but produces a different order of fixed mutations. By comparing several relevant evolutionary metrics, such as the distribution of fixed selection coefficients and the probability of reversion, we find that the two algorithms behave equivalently in many respects. However, the accelerated algorithm yields less variance in fixed selection coefficients. Notably, we are able to recover the expected amount of variance by rescaling population size, and we find a linear relationship between the rescaled population size and the population size used by the original algorithm. Considering the widespread usage of origin-fixation simulations across many areas of evolutionary biology, we introduce our accelerated algorithm as a useful tool for increasing the computational complexity of fitness functions without sacrificing much in terms of accuracy of the evolutionary simulation.

SUBMITTER: Teufel AI

PROVIDER: S-EPMC5332577 | biostudies-literature | 2017 Feb

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Accelerated simulation of evolutionary trajectories in origin-fixation models.

Teufel Ashley I AI Wilke Claus O CO

Journal of the Royal Society, Interface 20170201 127

We present an accelerated algorithm to forward-simulate origin-fixation models. Our algorithm requires, on average, only about two fitness evaluations per fixed mutation, whereas traditional algorithms require, per one fixed mutation, a number of fitness evaluations of the order of the effective population size, <i>N</i><sub>e</sub> Our accelerated algorithm yields the exact same steady state as the original algorithm but produces a different order of fixed mutations. By comparing several releva ...[more]

PMID: 28228542

Similar Datasets

Project description:Background The COVID-19 epidemic has shown that efficient prediction models are required, and the well-known SI, SIR, and SEIR models are not always capable of capturing the real dynamics. Modified models with novel structures could help identify unknown mechanisms of COVID-19 spread. Objective Our objective is to provide additional insights into the COVID-19 spread mechanisms based on different models’ parameterization which was performed using evolutionary algorithms and the first-wave data. Methods Data from the Our World in Data COVID-19 database was analysed, and several models—SI, SIR, SEIR, SEIUR, and Bass diffusion—and their variations were considered for the first wave of the COVID-19 pandemic. The models’ parameters were tuned with differential evolution optimization method L-SHADE to find the best fit. The algorithm for the automatic identification of the first wave was developed, and the differential evolution was applied to model parameterization. The reproduction rates (R0) for the first wave were calculated for 61 countries based on the best fits. Results The performed experiments showed that the Bass diffusion model-based modification could be superior compared to SI, SIR, SEIR and SEIUR due to the component responsible for spread from an external factor, which is not directly dependent on contact with infected individuals. The developed modified models containing this component were shown to perform better when fitting to the first-wave cumulative infections curve. In particular, the modified SEIR model was better fitted to the real-world data than the classical SEIR in 43 cases out of 61, based on Mann–Whitney U tests; the Bass diffusion model was better than SI for 57 countries. This showed the limitation of the classical models and indicated ways to improve them. Conclusions By using the modified models, the mechanism of infection spread, which is not directly dependent on contacts, was identified, which significantly influences the dynamics of the spread of COVID-19.

Project description:Understanding how novel functions evolve (genetic adaptation) is a critical goal of evolutionary biology. Among asexual organisms, genetic adaptation involves multiple mutations that frequently interact in a non-linear fashion (epistasis). Non-linear interactions pose a formidable challenge for the computational prediction of mutation effects. Here we use the recent evolution of β-lactamase under antibiotic selection as a model for genetic adaptation. We build a network of coevolving residues (possible functional interactions), in which nodes are mutant residue positions and links represent two positions found mutated together in the same sequence. Most often these pairs occur in the setting of more complex mutants. Focusing on extended-spectrum resistant sequences, we use network-theoretical tools to identify triple mutant trajectories of likely special significance for adaptation. We extrapolate evolutionary paths (n = 3) that increase resistance and that are longer than the units used to build the network (n = 2). These paths consist of a limited number of residue positions and are enriched for known triple mutant combinations that increase cefotaxime resistance. We find that the pairs of residues used to build the network frequently decrease resistance compared to their corresponding singlets. This is a surprising result, given that their coevolution suggests a selective advantage. Thus, β-lactamase adaptation is highly epistatic. Our method can identify triplets that increase resistance despite the underlying rugged fitness landscape and has the unique ability to make predictions by placing each mutant residue position in its functional context. Our approach requires only sequence information, sufficient genetic diversity, and discrete selective pressures. Thus, it can be used to analyze recent evolutionary events, where coevolution analysis methods that use phylogeny or statistical coupling are not possible. Improving our ability to assess evolutionary trajectories will help predict the evolution of clinically relevant genes and aid in protein design.

Dataset Information

Accelerated simulation of evolutionary trajectories in origin-fixation models.

Publications

Accelerated simulation of evolutionary trajectories in origin-fixation models.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets