124 research outputs found
Multidimensional Epistasis and the Advantage of Sex
Kondrashov and Kondrashov (2001) suggest that there is usually a disadvantage for sex in systems with multidimensional epistasis. They define systems of 'unidimensional epistasis' to be those where the fitness of a genotype is a function of the number of mutations it carries, and in contrast describe a system where the fitness of a genotype is a function of the numbers of mutations in two (or more) disjoint subsets of loci creating 'multidimensional epistasis'. In an example landscape an asexual population evolves fit genotypes about twice as fast as a sexual one. Here we examine other landscapes with multidimensional epistasis and find cases where an asexual population evolves fit genotypes 20 and 180 times slower than a sexual population
Evaluating the structure coefficient theorem of evolutionary game theory
In order to accommodate the empirical fact that population structures are
rarely simple, modern studies of evolutionary dynamics allow for complicated
and highly-heterogeneous spatial structures. As a result, one of the most
difficult obstacles lies in making analytical deductions, either qualitative or
quantitative, about the long-term outcomes of evolution. The "structure
coefficient theorem" is a well-known approach to this problem for
mutation-selection processes under weak selection, but a general method of
evaluating the terms it comprises is lacking. Here, we provide such a method
for populations of fixed (but arbitrary) size and structure, using easily
interpretable demographic measures. This method encompasses a large family of
evolutionary update mechanisms and extends the theorem to allow for asymmetric
contests to provide a better understanding of the mutation-selection balance
under more realistic circumstances. We apply the method to study social goods
produced and distributed among individuals in spatially-heterogeneous
populations, where asymmetric interactions emerge naturally and the outcome of
selection varies dramatically depending on the nature of the social good, the
spatial topology, and frequency with which mutations arise.Comment: 49 page
Recommended from our members
Extending Coalescent Theory to Autotetraploids
We develop coalescent models for autotetraploid species with tetrasomic inheritance. We show that the ancestral genetic process in a large population without recombination may be approximated using Kingman’s standard coalescent, with a coalescent effective population size 4N. Numerical results suggest that this approximation is accurate for population sizes on the order of hundreds of individuals. Therefore, existing coalescent simulation programs can be adapted to study population history in autotetraploids simply by interpreting the timescale in units of 4N generations. We also consider the possibility of double reduction, a phenomenon unique to polysomic inheritance, and show that its effects on gene genealogies are similar to partial self-fertilization.Organismic and Evolutionary Biolog
Loss and Recovery of Genetic Diversity in Adapting Populations of HIV
The evolution of drug resistance in HIV occurs by the fixation of specific, well-known, drug-resistance mutations, but the underlying population genetic processes are not well understood. By analyzing within-patient longitudinal sequence data, we make four observations that shed a light on the underlying processes and allow us to infer the short-term effective population size of the viral population in a patient. Our first observation is that the evolution of drug resistance usually occurs by the fixation of one drug-resistance mutation at a time, as opposed to several changes simultaneously. Second, we find that these fixation events are accompanied by a reduction in genetic diversity in the region surrounding the fixed drug-resistance mutation, due to the hitchhiking effect. Third, we observe that the fixation of drug-resistance mutations involves both hard and soft selective sweeps. In a hard sweep, a resistance mutation arises in a single viral particle and drives all linked mutations with it when it spreads in the viral population, which dramatically reduces genetic diversity. On the other hand, in a soft sweep, a resistance mutation occurs multiple times on different genetic backgrounds, and the reduction of diversity is weak. Using the frequency of occurrence of hard and soft sweeps we estimate the effective population size of HIV to be ( confidence interval ). This number is much lower than the actual number of infected cells, but much larger than previous population size estimates based on synonymous diversity. We propose several explanations for the observed discrepancies. Finally, our fourth observation is that genetic diversity at non-synonymous sites recovers to its pre-fixation value within 18 months, whereas diversity at synonymous sites remains depressed after this time period. These results improve our understanding of HIV evolution and have potential implications for treatment strategies
Pacific Salmon and the Coalescent Effective Population Size
Pacific salmon include several species that are both commercially important and endangered. Understanding the causes of loss in genetic variation is essential for designing better conservation strategies. Here we use a coalescent approach to analyze a model of the complex life history of salmon, and derive the coalescent effective population (CES). With the aid of Kronecker products and a convergence theorem for Markov chains with two time scales, we derive a simple formula for the CES and thereby establish its existence. Our results may be used to address important questions regarding salmon biology, in particular about the loss of genetic variation. To illustrate the utility of our approach, we consider the effects of fluctuations in population size over time. Our analysis enables the application of several tools of coalescent theory to the case of salmon
Discordance of Species Trees with Their Most Likely Gene Trees
Because of the stochastic way in which lineages sort during speciation, gene trees may differ in topology from each other and from species trees. Surprisingly, assuming that genetic lineages follow a coalescent model of within-species evolution, we find that for any species tree topology with five or more species, there exist branch lengths for which gene tree discordance is so common that the most likely gene tree topology to evolve along the branches of a species tree differs from the species phylogeny. This counterintuitive result implies that in combining data on multiple loci, the straightforward procedure of using the most frequently observed gene tree topology as an estimate of the species tree topology can be asymptotically guaranteed to produce an incorrect estimate. We conclude with suggestions that can aid in overcoming this new obstacle to accurate genomic inference of species phylogenies
Correlation of coalescence times in a diploid Wright-Fisher model with recombination and selfing
The correlation among the gene genealogies at different loci is crucial in
biology, yet challenging to understand because such correlation depends on many
factors including genetic linkage, recombination, natural selection and
population structure. Based on a diploid Wright-Fisher model with a single
mating type and partial selfing for a constant large population with size ,
we quantify the combined effect of genetic drift and two competing factors,
recombination and selfing, on the correlation of coalescence times at two
linked loci for samples of size two. Recombination decouples the genealogies at
different loci and decreases the correlation while selfing increases the
correlation. We obtain explicit asymptotic formulas for the correlation for
four scaling scenarios that depend on whether the selfing probability and the
recombination probability are of order or as tends to
infinity. Our analytical results confirm that the asymptotic lower bound in
[King, Wakeley, Carmi (Theor. Popul. Biol. 2018)] is sharp when the loci are
unlinked and when there is no selfing, and provide a number of new formulas for
other scaling scenarios that have not been considered before. We present
asymptotic results for the variance of Tajima's estimator of the population
mutation rate for infinitely many loci as tends to infinity. When the
selfing probability is of order and is equal to a positive constant
for all and if the samples at both loci are in the same individual, then
the variance of the Tajima's estimator tends to (hence remains positive)
even when the recombination rate, the number of loci and the population size
all tend to infinity.Comment: 39 pages, 6 figure
Recommended from our members
Bayesian Nonparametric Inference of Population Size Changes from Sequential Genealogies
Sophisticated inferential tools coupled with the coalescent model have recently emerged for estimating past population sizes from genomic data. Recent methods that model recombination require small sample sizes, make constraining assumptions about population size changes, and do not report measures of uncertainty for estimates. Here, we develop a Gaussian process-based Bayesian nonparametric method coupled with a sequentially Markov coalescent model that allows accurate inference of population sizes over time from a set of genealogies. In contrast to current methods, our approach considers a broad class of recombination events, including those that do not change local genealogies. We show that our method outperforms recent likelihood-based methods that rely on discretization of the parameter space. We illustrate the application of our method to multiple demographic histories, including population bottlenecks and exponential growth. In simulation, our Bayesian approach produces point estimates four times more accurate than maximum-likelihood estimation (based on the sum of absolute differences between the truth and the estimated values). Further, our method’s credible intervals for population size as a function of time cover 90% of true values across multiple demographic scenarios, enabling formal hypothesis testing about population size differences over time. Using genealogies estimated with ARGweaver, we apply our method to European and Yoruban samples from the 1000 Genomes Project and confirm key known aspects of population size history over the past 150,000 years
- …