SANTA-SIM: simulating viral sequence evolution dynamics under selection and recombination.
Ontology highlight
ABSTRACT:
SUBMITTER: Jariani A
PROVIDER: S-EPMC6407609 | biostudies-literature | 2019 Jan
REPOSITORIES: biostudies-literature
Publications
SANTA-SIM: simulating viral sequence evolution dynamics under selection and recombination.
Jariani Abbas A Warth Christopher C Deforche Koen K Libin Pieter P Drummond Alexei J AJ Rambaut Andrew A Matsen Iv Frederick A FA Theys Kristof K
Virus evolution 20190101 1
Simulations are widely used to provide expectations and predictive distributions under known conditions against which to compare empirical data. Such simulations are also invaluable for testing and comparing the behaviour and power of inference methods. We describe SANTA-SIM, a software package to simulate the evolution of a population of gene sequences forwards through time. It models the underlying biological processes as discrete components: replication, recombination, point mutations, insert ...[more]
Similar Datasets
Project description:SLiM is an efficient forward population genetic simulation designed for studying the effects of linkage and selection on a chromosome-wide scale. The program can incorporate complex scenarios of demography and population substructure, various models for selection and dominance of new mutations, arbitrary gene structure, and user-defined recombination maps.
| S-EPMC3730910 | biostudies-literature
Project description:The within-host evolution of influenza is a vital component of its epidemiology. A question of particular interest is the role that selection plays in shaping the viral population over the course of a single infection. We here describe a method to measure selection acting upon the influenza virus within an individual host, based upon time-resolved genome sequence data from an infection. Analysing sequence data from a transmission study conducted in pigs, describing part of the haemagglutinin gene (HA1) of an influenza virus, we find signatures of non-neutrality in six of a total of sixteen infections. We find evidence for both positive and negative selection acting upon specific alleles, while in three cases, the data suggest the presence of time-dependent selection. In one infection we observe what is potentially a specific immune response against the virus; a non-synonymous mutation in an epitope region of the virus is found to be under initially positive, then strongly negative selection. Crucially, given the lack of homologous recombination in influenza, our method accounts for linkage disequilibrium between nucleotides at different positions in the haemagglutinin gene, allowing for the analysis of populations in which multiple mutations are present at any given time. Our approach offers a new insight into the dynamics of influenza infection, providing a detailed characterisation of the forces that underlie viral evolution.
| S-EPMC4117419 | biostudies-literature
Project description:We study the emergence of conditional cooperation in the presence of both intra-group and inter-group selection. Individuals play public goods games within their groups using conditional strategies, which are represented as piecewise linear response functions. Accordingly, groups engage in conflicts with a certain probability. In contrast to previous studies, we consider continuous contribution levels and a rich set of conditional strategies, allowing for a wide range of possible interactions between strategies. We find that the existence of conditional strategies enables the stabilization of cooperation even under strong intra-group selection. The strategy that eventually dominates in the population has two key properties: (i) It is unexploitable with strong intra-group selection; (ii) It can achieve full contribution to outperform other strategies in the inter-group selection. The success of this strategy is robust to initial conditions as well as changes to important parameters. We also investigate the influence of different factors on cooperation levels, including group conflicts, group size, and migration rate. Their effect on cooperation can be attributed to and explained by their influence on the relative strength of intra-group and inter-group selection.
| S-EPMC4786827 | biostudies-other
Project description:MotivationConventional phylogenetic analysis for characterizing the relatedness between taxa typically assumes that a single relationship exists between species at every site along the genome. This assumption fails to take into account recombination which is a fundamental process for generating diversity and can lead to spurious results. Recombination induces a localized phylogenetic structure which may vary along the genome. Here, we generalize a hidden Markov model (HMM) to infer changes in phylogeny along multiple sequence alignments while accounting for rate heterogeneity; the hidden states refer to the unobserved phylogenic topology underlying the relatedness at a genomic location. The dimensionality of the number of hidden states (topologies) and their structure are random (not known a priori) and are sampled using Markov chain Monte Carlo algorithms. The HMM structure allows us to analytically integrate out over all possible changepoints in topologies as well as all the unknown branch lengths.ResultsWe demonstrate our approach on simulated data and also to the genome of a suspected HIV recombinant strain as well as to an investigation of recombination in the sequences of 15 laboratory mouse strains sequenced by Perlegen Sciences. Our findings indicate that our method allows us to distinguish between rate heterogeneity and variation in phylogeny caused by recombination without being restricted to 4-taxa data.
| S-EPMC2639012 | biostudies-literature
Project description:Multigene families-immunity genes or sensory receptors, for instance-are often subject to diversifying selection. Allelic diversity may be favored not only through balancing or frequency-dependent selection at individual loci but also by associating different alleles in multicopy gene families. Using a combination of analytical calculations and simulations, we explored a population genetic model of epistatic selection and unequal recombination, where a trade-off exists between the benefit of allelic diversity and the cost of copy abundance. Starting from the neutral case, where we showed that gene copy number is Gamma distributed at equilibrium, we derived also the mean and shape of the limiting distribution under selection. Considering a more general model, which includes variable population size and population substructure, we explored by simulations mean fitness and some summary statistics of the copy number distribution. We determined the relative effects of selection, recombination, and demographic parameters in maintaining allelic diversity and shaping the mean fitness of a population. One way to control the variance of copy number is by lowering the rate of unequal recombination. Indeed, when encoding recombination by a rate modifier locus, we observe exactly this prediction. Finally, we analyzed the empirical copy number distribution of 3 genes in human and estimated recombination and selection parameters of our model.
| S-EPMC9252282 | biostudies-literature
Project description:Plant microRNAs (miRNA) guide cleavage of target mRNAs by DICER-like proteins, thereby reducing mRNA abundance. Native precursor miRNAs can be redesigned to target RNAs of interest, and one application of such artificial microRNA (amiRNA) technology is to generate plants resistant to pathogenic viruses. Transgenic Arabidopsis plants expressing amiRNAs designed to target the genome of two unrelated viruses were resistant, in a highly specific manner, to the appropriate virus. Here, we pursued two different goals. First, we confirmed that the 21-nt target site of viral RNAs is both necessary and sufficient for resistance. Second, we studied the evolutionary stability of amiRNA-mediated resistance against a genetically plastic RNA virus, TuMV. To dissociate selective pressures acting upon protein function from those acting at the RNA level, we constructed a chimeric TuMV harboring a 21-nt, amiRNA target site in a non-essential region. In the first set of experiments designed to assess the likelihood of resistance breakdown, we explored the effect of single nucleotide mutation within the target 21-nt on the ability of mutant viruses to successfully infect amiRNA-expressing plants. We found non-equivalency of the target nucleotides, which can be divided into three categories depending on their impact in virus pathogenicity. In the second set of experiments, we investigated the evolution of the virus mutants in amiRNA-expressing plants. The most common outcome was the deletion of the target. However, when the 21-nt target was retained, viruses accumulated additional substitutions on it, further reducing the binding/cleavage ability of the amiRNA. The pattern of substitutions within the viral target was largely dominated by G to A and C to U transitions.
| S-EPMC2642722 | biostudies-other
Project description:Infectious disease is a great enemy of humankind. The ravages of COVID-19 are leading to profound crises across the world. There is an urgent requirement for analyzing the current pandemic situation, predicting trends over time, and assessing the effectiveness of containment measures. Thus, numerous statistical models, primarily based on the susceptible-exposed-infected-recovered or removed (SEIR) model, have been established. However, these models are highly technical, which are difficult for the public and governing bodies to understand and use. To address this issue, we developed a simple operating software based on our improved K-SEIR model termed as the kernelkernel SEIR simulator (K-SEIR-Sim). This software includes natural propagation parameters, containment measure parameters, and certain characteristic parameters that can deduce the effects of natural propagation and containment measures. Further, the applicability of the proposed software was demonstrated using the example of the COVID-19 outbreak in the United States and the city of Wuhan, China. Operating results verified the potency of the proposed software in evaluating the epidemic situation and human intervention during COVID-19. Importantly, the software can perform real-time, backward-looking, and forward-looking analysis by functioning in data-driven and model-driven ways. All of them have considerable practical values in their applications according to the actual needs of personal use. Conclusively, K-SEIR-Sim is the first simple customized operating software that is highly valuable for the global fight against COVID-19 and other infectious diseases.
| S-EPMC8025586 | biostudies-literature
Project description:The evolutionary dynamics of HIV during the chronic phase of infection is driven by the host immune response and by selective pressures exerted through drug treatment. To understand and model the evolution of HIV quantitatively, the parameters governing genetic diversification and the strength of selection need to be known. While mutation rates can be measured in single replication cycles, the relevant effective recombination rate depends on the probability of coinfection of a cell with more than one virus and can only be inferred from population data. However, most population genetic estimators for recombination rates assume absence of selection and are hence of limited applicability to HIV, since positive and purifying selection are important in HIV evolution. Yet, little is known about the distribution of selection differentials between individual viruses and the impact of single polymorphisms on viral fitness. Here, we estimate the rate of recombination and the distribution of selection coefficients from time series sequence data tracking the evolution of HIV within single patients. By examining temporal changes in the genetic composition of the population, we estimate the effective recombination to be rho = 1.4+/-0.6 x 10(-5) recombinations per site and generation. Furthermore, we provide evidence that the selection coefficients of at least 15% of the observed non-synonymous polymorphisms exceed 0.8% per generation. These results provide a basis for a more detailed understanding of the evolution of HIV. A particularly interesting case is evolution in response to drug treatment, where recombination can facilitate the rapid acquisition of multiple resistance mutations. With the methods developed here, more precise and more detailed studies will be possible as soon as data with higher time resolution and greater sample sizes are available.
| S-EPMC2813257 | biostudies-literature
Project description:We use first principles of population genetics to model the evolution of proteins under persistent positive selection (PPS). PPS may occur when organisms are subjected to persistent environmental change, during adaptive radiations, or in host-pathogen interactions. Our mutation-selection model indicates protein evolution under PPS is an irreversible Markov process, and thus proteins under PPS show a strongly asymmetrical distribution of selection coefficients among amino acid substitutions. Our model shows the criteria ω>1 (where ω is the ratio of nonsynonymous over synonymous codon substitution rates) to detect positive selection is conservative and indeed arbitrary, because in real proteins many mutations are highly deleterious and are removed by selection even at positively selected sites. We use a penalized-likelihood implementation of the PPS model to successfully detect PPS in plant RuBisCO and influenza HA proteins. By directly estimating selection coefficients at protein sites, our inference procedure bypasses the need for using ω as a surrogate measure of selection and improves our ability to detect molecular adaptation in proteins.
| S-EPMC8760937 | biostudies-literature
Project description:MotivationSimulation is an essential technique for generating biomolecular data with a 'known' history for use in validating phylogenetic inference and other evolutionary methods. On longer time scales, simulation supports investigations of equilibrium behavior and provides a formal framework for testing competing evolutionary hypotheses. Twenty years of molecular evolution research have produced a rich repertoire of simulation methods. However, current models do not capture the stringent constraints acting on the domain insertions, duplications, and deletions by which multidomain architectures evolve. Although these processes have the potential to generate any combination of domains, only a tiny fraction of possible domain combinations are observed in nature. Modeling these stringent constraints on domain order and co-occurrence is a fundamental challenge in domain architecture simulation that does not arise with sequence and gene family simulation.ResultsHere, we introduce a stochastic model of domain architecture evolution to simulate evolutionary trajectories that reflect the constraints on domain order and co-occurrence observed in nature. This framework is implemented in a novel domain architecture simulator, DomArchov, using the Metropolis-Hastings algorithm with data-driven transition probabilities. The use of a data-driven event module enables quick and easy redeployment of the simulator for use in different taxonomic and protein function contexts. Using empirical evaluation with metazoan datasets, we demonstrate that domain architectures simulated by DomArchov recapitulate properties of genuine domain architectures that reflect the constraints on domain order and adjacency seen in nature. This work expands the realm of evolutionary processes that are amenable to simulation.Availability and implementationDomArchov is written in Python 3 and is available at http://www.cs.cmu.edu/~durand/DomArchov. The data underlying this article are available via the same link.Supplementary informationSupplementary data are available at Bioinformatics online.
| S-EPMC9236583 | biostudies-literature