skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Fixation probabilities in graph-structured populations under weak selection
A population’s spatial structure affects the rate of genetic change and the outcome of natural selection. These effects can be modeled mathematically using the Birth-death process on graphs. Individuals occupy the vertices of a weighted graph, and reproduce into neighboring vertices based on fitness. A key quantity is the probability that a mutant type will sweep to fixation, as a function of the mutant’s fitness. Graphs that increase the fixation probability of beneficial mutations, and decrease that of deleterious mutations, are said to amplify selection. However, fixation probabilities are difficult to compute for an arbitrary graph. Here we derive an expression for the fixation probability, of a weakly-selected mutation, in terms of the time for two lineages to coalesce. This expression enables weak-selection fixation probabilities to be computed, for an arbitrary weighted graph, in polynomial time. Applying this method, we explore the range of possible effects of graph structure on natural selection, genetic drift, and the balance between the two. Using exhaustive analysis of small graphs and a genetic search algorithm, we identify families of graphs with striking effects on fixation probability, and we analyze these families mathematically. Our work reveals the nuanced effects of graph structure on natural selection and neutral drift. In particular, we show how these notions depend critically on the process by which mutations arise.  more » « less
Award ID(s):
1715315
PAR ID:
10297328
Author(s) / Creator(s):
; ; ; ; ; ;
Editor(s):
Traulsen, Arne
Date Published:
Journal Name:
PLOS Computational Biology
Volume:
17
Issue:
2
ISSN:
1553-7358
Page Range / eLocation ID:
e1008695
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Barton, N (Ed.)
    Abstract In rapidly evolving populations, numerous beneficial and deleterious mutations can arise and segregate within a population at the same time. In this regime, evolutionary dynamics cannot be analyzed using traditional population genetic approaches that assume that sites evolve independently. Instead, the dynamics of many loci must be analyzed simultaneously. Recent work has made progress by first analyzing the fitness variation within a population, and then studying how individual lineages interact with this traveling fitness wave. However, these “traveling wave” models have previously been restricted to extreme cases where selection on individual mutations is either much faster or much slower than the typical coalescent timescale Tc. In this work, we show how the traveling wave framework can be extended to intermediate regimes in which the scaled fitness effects of mutations (Tcs) are neither large nor small compared to one. This enables us to describe the dynamics of populations subject to a wide range of fitness effects, and in particular, in cases where it is not immediately clear which mutations are most important in shaping the dynamics and statistics of genetic diversity. We use this approach to derive new expressions for the fixation probabilities and site frequency spectra of mutations as a function of their scaled fitness effects, along with related results for the coalescent timescale Tc and the rate of adaptation or Muller’s ratchet. We find that competition between linked mutations can have a dramatic impact on the proportions of neutral and selected polymorphisms, which is not simply summarized by the scaled selection coefficient Tcs. We conclude by discussing the implications of these results for population genetic inferences. 
    more » « less
  2. Abstract Population structure affects the outcome of natural selection. These effects can be modeled using evolutionary games on graphs. Recently, conditions were derived for a trait to be favored under weak selection, on any weighted graph, in terms of coalescence times of random walks. Here we consider isothermal graphs, which have the same total edge weight at each node. The conditions for success on isothermal graphs take a simple form, in which the effects of graph structure are captured in the ‘effective degree’—a measure of the effective number of neighbors per individual. For two update rules (death-Birth and birth-Death), cooperative behavior is favored on a large isothermal graph if the benefit-to-cost ratio exceeds the effective degree. For two other update rules (Birth-death and Death-birth), cooperation is never favored. We relate the effective degree of a graph to its spectral gap, thereby linking evolutionary dynamics to the theory of expander graphs. Surprisingly, we find graphs of infinite average degree that nonetheless provide strong support for cooperation. 
    more » « less
  3. Abstract Cancer is a genetic disease that results from accumulation of unfavorable mutations. As soon as genetic and epigenetic modifications associated with these mutations become strong enough, the uncontrolled tumor cell growth is initiated, eventually spreading through healthy tissues. Clarifying the dynamics of cancer initiation is thus critically important for understanding the molecular mechanisms of tumorigenesis. Here we present a new theoretical method to evaluate the dynamic processes associated with the cancer initiation. It is based on a discrete-state stochastic description of the formation of tumors as a fixation of cancerous mutations in tissues. Using a first-passage analysis the probabilities for the cancer to appear and the times before it happens, which are viewed as fixation probabilities and fixation times, respectively, are explicitly calculated. It is predicted that the slowest cancer initiation dynamics is observed for neutral mutations, while it is fast for both advantageous and, surprisingly, disadvantageous mutations. The method is applied for estimating the cancer initiation times from experimentally available lifetime cancer risks for different types of cancer. It is found that the higher probability of the cancer to occur does not necessary lead to the faster times of starting the cancer. Our theoretical analysis helps to clarify microscopic aspects of cancer initiation processes. 
    more » « less
  4. Heritable variation in gene expression is common within and among species and contributes to phenotypic diversity. Mutations affecting eithercis- ortrans-regulatory sequences controlling gene expression give rise to variation in gene expression, and natural selection acting on this variation causes some regulatory variants to persist in a population for longer than others. To understand how mutation and selection interact to produce the patterns of regulatory variation we see within and among species, my colleagues and I have been systematically determining the effects of new mutations on expression of theTDH3gene inSaccharomyces cerevisiaeand comparing them to the effects of polymorphisms segregating within this species. We have also investigated the molecular mechanisms by which regulatory variants act. Over the past decade, this work has revealed properties ofcis- andtrans-regulatory mutations including their relative frequency, effects, dominance, pleiotropy and fitness consequences. Comparing these mutational effects to the effects of polymorphisms in natural populations, we have inferred selection acting on expression level, expression noise and phenotypic plasticity. Here, I summarize this body of work and synthesize its findings to make inferences not readily discernible from the individual studies alone. This article is part of the theme issue ‘Interdisciplinary approaches to predicting evolutionary biology’. 
    more » « less
  5. Abstract Adaptation dynamics on fitness landscapes is often studied theoretically in the strong-selection, weak-mutation regime. However, in a large population, multiple beneficial mutants can emerge before any of them fixes in the population. Competition between mutants is known as clonal interference, and while it is known to slow down the rate of adaptation (when compared to the strong-selection, weak-mutation model with the same parameters), how it affects the shape of long-term fitness trajectories in the presence of epistasis is an open question. Here, by considering how changes in fixation probabilities arising from weak clonal interference affect the dynamics of adaptation on fitness-parameterized landscapes, we find that the change in the shape of fitness trajectory arises only through changes in the supply of beneficial mutations (or equivalently, the beneficial mutation rate). Furthermore, a depletion of beneficial mutations as a population climbs up the fitness landscape can speed up the rescaled fitness trajectory (where adaptation speed is measured relative to its value at the start of the experiment), while an enhancement of the beneficial mutation rate does the opposite of slowing it down. Our findings suggest that by carrying out evolution experiments in both regimes (with and without clonal interference), one could potentially distinguish the different sources of macroscopic epistasis (fitness effect of mutations vs change in fraction of beneficial mutations). 
    more » « less