SLiM: simulating evolution with selection and linkage.
Ontology highlight
ABSTRACT:
SUBMITTER: Messer PW
PROVIDER: S-EPMC3730910 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description: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, insertion-deletions, and selection under various fitness models and population size dynamics. The software is designed to be intuitive to work with for a wide range of users and executable in a cross-platform manner.
| S-EPMC6407609 | biostudies-literature
Project description:The influenza virus is an important human pathogen, with a rapid rate of evolution in the human population. The rate of homologous recombination within genes of influenza is essentially zero. As such, where two alleles within the same gene are in linkage disequilibrium, interference between alleles will occur, whereby selection acting upon one allele has an influence upon the frequency of the other. We here measured the relative importance of selection and interference effects upon the evolution of influenza. We considered time-resolved allele frequency data from the global evolutionary history of the haemagglutinin gene of human influenza A/H3N2, conducting an in-depth analysis of sequences collected since 1996. Using a model that accounts for selection-caused interference between alleles in linkage disequilibrium, we estimated the inherent selective benefit of individual polymorphisms in the viral population. These inherent selection coefficients were in turn used to calculate the total selective effect of interference acting upon each polymorphism, considering the effect of the initial background upon which a mutation arose, and the subsequent effect of interference from other alleles that were under selection. Viewing events in retrospect, we estimated the influence of each of these components in determining whether a mutant allele eventually fixed or died in the global viral population. Our inherent selection coefficients, when combined across different regions of the protein, were consistent with previous measurements of dN/dS for the same system. Alleles going on to fix in the global population tended to be under more positive selection, to arise on more beneficial backgrounds, and to avoid strong negative interference from other alleles under selection. However, on average, the fate of a polymorphism was determined more by the combined influence of interference effects than by its inherent selection coefficient.
| S-EPMC3531508 | 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
Project description:Androgen receptor (AR) mutation is closely associated with prostate cancer (PCa) and is one of the mechanisms of resistance to PCa therapies such as AR antagonists. Although sequencing technologies like next-generation sequencing (NGS) contributes to the high-throughput and precise detection of AR mutations carried by PCa patients, the lack of interpretations of these clinical genetic variants has still been a roadblock for PCa-targeted precision medicine. Here, we established a designer yeast reporter assay to simulate natural androgen receptor (AR) selection using AR antagonists. Yeast HIS3 gene transactivation was associated with the ligand-induced recruitment of steroid receptor coactivator-1 (SRC-1) by AR mutants, where yeast growth in histidine-free medium was determined as the outcome. This assay is applicable to determine a wide range of clinical AR mutants including those with loss of function relating to androgen insensitivity syndrome (AIS), and those associated with PCa conferring resistance to AR antagonists such as enzalutamide (ENZ), bicalutamide (BIC), and cyproterone acetate (CPA). One clinical AR mutant previously reported to confer ENZ-resistance, F877L, was found to confer partial resistance to CPA as well using designer yeast. Our simple and efficient assay can enable precise one-pot screening of AR mutants, providing a reference for tailored medicine.
| S-EPMC9386396 | biostudies-literature
Project description:Existing simulation methods usually simulate linkage disequilibrium (LD) structures starting with an initial population that is randomly generated according to specified allele frequencies. These at random based methods might be unstable because the LD level of the initial population is generally extremely low. This study presents a new algorithm, SIMLD, to simulate genome populations with real LD structures. SIMLD begins from an initial population with possibly the highest LD level, and then the LD decays to fit the desired level through processes of mating and recombination over generations. SIMLD can produce case-control samples according to various disease models. Using empirical SNP marker information from three populations of HapMap data, we implement the proposed algorithm and demonstrate a set of experimental results.
| S-EPMC4116680 | biostudies-literature
Project description:BACKGROUND:To evaluate statistical methods for genome-wide genetic analyses, one needs to be able to simulate realistic genotypes. We here describe a method, applicable to a broad range of association study designs, that can simulate autosome-wide single-nucleotide polymorphism data with realistic linkage disequilibrium and with spiked in, user-specified, single or multi-SNP causal effects. RESULTS:Our construction uses existing genome-wide association data from unrelated case-parent triads, augmented by including a hypothetical complement triad for each triad (same parents but with a hypothetical offspring who carries the non-transmitted parental alleles). We assign offspring qualitative or quantitative traits probabilistically through a specified risk model and show that our approach destroys the risk signals from the original data. Our method can simulate genetically homogeneous or stratified populations and can simulate case-parents studies, case-control studies, case-only studies, or studies of quantitative traits. We show that allele frequencies and linkage disequilibrium structure in the original genome-wide association sample are preserved in the simulated data. We have implemented our method in an R package (TriadSim) which is freely available at the comprehensive R archive network. CONCLUSION:We have proposed a method for simulating genome-wide SNP data with realistic linkage disequilibrium. Our method will be useful for developing statistical methods for studying genetic associations, including higher order effects like epistasis and gene by environment interactions.
| S-EPMC5749028 | biostudies-literature
Project description:Sex-biased genes--genes that are differentially expressed within males and females--are nonrandomly distributed across animal genomes, with sex chromosomes and autosomes often carrying markedly different concentrations of male- and female-biased genes. These linkage patterns are often gene- and lineage-dependent, differing between functional genetic categories and between species. Although sex-specific selection is often hypothesized to shape the evolution of sex-linked and autosomal gene content, population genetics theory has yet to account for many of the gene- and lineage-specific idiosyncrasies emerging from the empirical literature. With the goal of improving the connection between evolutionary theory and a rapidly growing body of genome-wide empirical studies, we extend previous population genetics theory of sex-specific selection by developing and analyzing a biologically informed model that incorporates sex linkage, pleiotropy, recombination, and epistasis, factors that are likely to vary between genes and between species. Our results demonstrate that sex-specific selection and sex-specific recombination rates can generate, and are compatible with, the gene- and species-specific linkage patterns reported in the genomics literature. The theory suggests that sexual selection may strongly influence the architectures of animal genomes, as well as the chromosomal distribution of fixed substitutions underlying sexually dimorphic traits.
| S-EPMC2998557 | biostudies-literature
Project description:The intricate adjustment of organisms to their environment demonstrates the effectiveness of natural selection. But Darwin himself recognized that certain biological features could limit this effectiveness, features that generally reduce the efficiency of natural selection or yield suboptimal adaptation. Genetic linkage is known to be one such feature, and here we show theoretically that it can introduce a more sinister flaw: when there is complete linkage between loci affecting fitness and loci affecting mutation rate, positive natural selection and recurrent mutation can drive mutation rates in an adapting population to intolerable levels. We discuss potential implications of this finding for the early establishment of recombination, the evolutionary fate of asexual populations, and immunological clearance of clonal pathogens.
| S-EPMC1851075 | biostudies-literature
Project description:Physiological and molecular evidences have shown earlier that low CO2 might have been a major driving force during the evolution of C4 photosynthesis, not only being a selection pressure but also as a signaling agent. However, a mechanistic linkage between low CO2 and C4 emergence is missing. In this study, using transcriptomics study in model plant Arabidopsis thaliana, we demonstrated that under long-term low CO2 treatments, the up-regulation of C4 related genes were linked to the up-regulation of genes involved in photorespiration, nitrogen assimilation and glycolysis. Plants under low CO2 also showed altered d13C. The carbonic anhydrase (CA), phosphoenolpyruvate carboxylase (PEPC), malate dehydrogenase (MDH), glutamine oxoglutarate aminotransferase (GS-GOGAT), which are required to reassimilate ammonia, were up-regulated under low CO2. Furthermore, under low CO2, genes involved in PEP regeneration from glycolysis were up-regulated while pyruvate orthophosphate dikinase (PPDK) was down-regulated, suggesting a route of PEP regeneration from glycolysis. All these results suggested that under long-term low CO2 treatment, the selection pressure to recapture the released ammonia from the increased photorespiration might have promoted the evolution of mechanisms for generation PEP from glycolysis and enhancement of enzymes catalyzing formation of oxaloacetate, one intermediate which can accept ammonia residue. These adjustments in metabolism provide a mechanistic linkage between low CO2 and evolution of C4 photosynthesis.
2016-12-31 | GSE64547 | GEO
Project description:This work describes the lightning nitric oxide (LNO) production schemes in the Community Multiscale Air Quality (CMAQ) model. We first document the existing LNO production scheme and vertical distribution algorithm. We then describe updates that were made to the scheme originally based on monthly National Lightning Detection Network (mNLDN) observations. The updated scheme uses hourly NLDN (hNLDN) observations. These NLDN-based schemes are good for retrospective model applications when historical lightning data are available. For applications when observed data are not available (i.e., air quality forecasts and climate studies that assume similar climate conditions), we have developed a scheme that is based on linear and log-linear parameters derived from regression of multiyear historical NLDN (pNLDN) observations and meteorological model simulations. Preliminary assessment for total column LNO production reveals that the mNLDN scheme overestimates LNO by over 40% during summer months compared with the updated hNLDN scheme that reflects the observed lightning activity more faithfully in time and space. The pNLDN performance varies with year, but it generally produced LNO columns that are comparable to hNLDN and mNLDN, and in most cases it outperformed mNLDN. Thus, when no observed lightning data are available, pNLDN can provide reasonable estimates of LNO emissions over time and space for this important natural NO source that influences air quality regulations.
| S-EPMC7087390 | biostudies-literature