Stirring does not make populations well mixed

Frequent asymmetric migrations suppress natural selection in spatially structured populations

Natural microbial populations often have complex spatial structures. This can impact their evolution, in particular the ability of mutants to take over. While mutant fixation probabilities are known to be unaffected by sufficiently symmetric structures, evolutionary graph theory has shown that some graphs can amplify or suppress natural selection, in a way that depends on microscopic update rules. We propose a model of spatially structured populations on graphs directly inspired by batch culture experiments, alternating within-deme growth on nodes and migration-dilution steps, and yielding successive bottlenecks. This setting bridges models from evolutionary graph theory with Wright-Fisher models. Using a branching process approach, we show that spatial structure with frequent migrations can only yield suppression of natural selection. More precisely, in this regime, circulation graphs, where the total incoming migration flow equals the total outgoing one in each deme, do not impact fixation probability, while all other graphs strictly suppress selection. Suppression becomes stronger as the asymmetry between incoming and outgoing migrations grows. Amplification of natural selection can nevertheless exist in a restricted regime of rare migrations and very small fitness advantages, where we recover the predictions of evolutionary graph theory for the star graph.

Coalescent dynamics of planktonic communities

Planktonic communities are extremely diverse and include a vast number of rare species. The dynamics of these rare species is best described by individual-based models. However, individual-based approaches to planktonic diversity face substantial difficulties, due to the large number of individuals required to make realistic predictions. In this paper, we study the diversity of planktonic communities by means of a spatial coalescence model that incorporates transport by oceanic currents. As a main advantage, our approach requires simulating a number of individuals equal to the size of the sample one is interested in, rather than the size of the entire community. By theoretical analysis and simulations, we explore the conditions upon which our coalescence model is equivalent to individual-based dynamics. As an application, we use our model to predict the impact of chaotic advection by oceanic currents on biodiversity. We conclude that the coalescent approach permits one to simulate marine microbial communities much more efficiently than with individual-based models.

A Prospective Method for Generating COVID-19 Dynamics

Generating dynamic operators are constructed here from the cumulative case function to recover all state dynamics of a Susceptible–Exposed–Infectious–Recovered (SEIR) model for COVID-19 transmission. In this study, recorded and unrecorded EIRs and a time-dependent infection rate are taken into account to accommodate immeasurable control and intervention processes. Generating dynamic operators are built and implemented on the cumulative cases. All infection processes, which are hidden in this cumulative function, can be recovered entirely by implementing the generating operators. Direct implementation of the operators on the cumulative function gives all recorded state dynamics. Further, the unrecorded daily infection rate is estimated from the ratio between IFR and CFR. The remaining dynamics of unrecorded states are directly obtained from the generating operators. The simulations are conducted using infection data provided by Worldometers from ten selected countries. It is shown that the higher number of daily PCR tests contributed directly to reducing the effective reproduction ratio. The simulations of all state dynamics, infection rates, and effective reproduction ratios for several countries in the first and second waves of transmissions are presented. This method directly measures daily transmission indicators, which can be effectively used for the day-to-day control of the epidemic.

Hydrodynamic flow and concentration gradients in the gut enhance neutral bacterial diversity

The gut microbiota features important genetic diversity, and the specific spatial features of the gut may shape evolution within this environment. We investigate the fixation probability of neutral bacterial mutants within a minimal model of the gut that includes hydrodynamic flow and resulting gradients of food and bacterial concentrations. We find that this fixation probability is substantially increased, compared with an equivalent well-mixed system, in the regime where the profiles of food and bacterial concentration are strongly spatially dependent. Fixation probability then becomes independent of total population size. We show that our results can be rationalized by introducing an active population, which consists of those bacteria that are actively consuming food and dividing. The active population size yields an effective population size for neutral mutant fixation probability in the gut.

Toward a Universal Model for Spatially Structured Populations

A key question in evolution is how likely a mutant is to take over. This depends on natural selection and on stochastic fluctuations. Population spatial structure can impact mutant fixation probabilities. We introduce a model for structured populations on graphs that generalizes previous ones by making migrations independent of birth and death. We demonstrate that by tuning migration asymmetry, the star graph transitions from amplifying to suppressing natural selection. The results from our model are universal in the sense that they do not hinge on a modeling choice of microscopic dynamics or update rules. Instead, they depend on migration asymmetry, which can be experimentally tuned and measured.

The Phylogeography of Potato Virus X Shows the Fingerprints of Its Human Vector

Potato virus X (PVX) occurs worldwide and causes an important potato disease. Complete PVX genomes were obtained from 326 new isolates from Peru, which is within the potato crop's main domestication center, 10 from historical PVX isolates from the Andes (Bolivia, Peru) or Europe (UK), and three from Africa (Burundi). Concatenated open reading frames (ORFs) from these genomes plus 49 published genomic sequences were analyzed. Only 18 of them were recombinants, 17 of them Peruvian. A phylogeny of the non-recombinant sequences found two major (I, II) and five minor (I-1, I-2, II-1, II-2, II-3) phylogroups, which included 12 statistically supported clusters. Analysis of 488 coat protein (CP) gene sequences, including 128 published previously, gave a completely congruent phylogeny. Among the minor phylogroups, I-2 and II-3 only contained Andean isolates, I-1 and II-2 were of both Andean and other isolates, but all of the three II-1 isolates were European. I-1, I-2, II-1 and II-2 all contained biologically typed isolates. Population genetic and dating analyses indicated that PVX emerged after potato's domestication 9000 years ago and was transported to Europe after the 15th century. Major clusters A-D probably resulted from expansions that occurred soon after the potato late-blight pandemic of the mid-19th century. Genetic comparisons of the PVX populations of different Peruvian Departments found similarities between those linked by local transport of seed potato tubers for summer rain-watered highland crops, and those linked to winter-irrigated crops in nearby coastal Departments. Comparisons also showed that, although the Andean PVX population was diverse and evolving neutrally, its spread to Europe and then elsewhere involved population expansion. PVX forms a basal Potexvirus genus lineage but its immediate progenitor is unknown. Establishing whether PVX's entirely Andean phylogroups I-2 and II-3 and its Andean recombinants threaten potato production elsewhere requires future biological studies.

The Moran process on 2-chromatic graphs

Resources are rarely distributed uniformly within a population. Heterogeneity in the concentration of a drug, the quality of breeding sites, or wealth can all affect evolutionary dynamics. In this study, we represent a collection of properties affecting the fitness at a given location using a color. A green node is rich in resources while a red node is poorer. More colors can represent a broader spectrum of resource qualities. For a population evolving according to the birth-death Moran model, the first question we address is which structures, identified by graph connectivity and graph coloring, are evolutionarily equivalent. We prove that all properly two-colored, undirected, regular graphs are evolutionarily equivalent (where “properly colored” means that no two neighbors have the same color). We then compare the effects of background heterogeneity on properly two-colored graphs to those with alternative schemes in which the colors are permuted. Finally, we discuss dynamic coloring as a model for spatiotemporal resource fluctuations, and we illustrate that random dynamic colorings often diminish the effects of background heterogeneity relative to a proper two-coloring.

Heterogeneity in environmental conditions can have profound effects on long-term evolutionary outcomes in structured populations. We consider a population evolving on a colored graph, wherein the color of a node represents the resources at that location. Using a combination of analytical and numerical methods, we quantify the effects of background heterogeneity on a population’s dynamics. In addition to considering the notion of an “optimal” coloring with respect to mutant invasion, we also study the effects of dynamic spatial redistribution of resources as the population evolves. Although the effects of static background heterogeneity can be quite striking, these effects are often attenuated by the movement (or “flow”) of the underlying resources.

Ocean currents promote rare species diversity in protists

Oceans host communities of plankton composed of relatively few abundant species and many rare species. The number of rare protist species in these communities, as estimated in metagenomic studies, decays as a steep power law of their abundance. The ecological factors at the origin of this pattern remain elusive. We propose that chaotic advection by oceanic currents affects biodiversity patterns of rare species. To test this hypothesis, we introduce a spatially explicit coalescence model that reconstructs the species diversity of a sample of water. Our model predicts, in the presence of chaotic advection, a steeper power law decay of the species abundance distribution and a steeper increase of the number of observed species with sample size. A comparison of metagenomic studies of planktonic protist communities in oceans and in lakes quantitatively confirms our prediction. Our results support that oceanic currents positively affect the diversity of rare aquatic microbes.

Evolution of cooperation on temporal networks

Population structure is a key determinant in fostering cooperation among naturally self-interested individuals in microbial populations, social insect groups, and human societies. Traditional research has focused on static structures, and yet most real interactions are finite in duration and changing in time, forming a temporal network. This raises the question of whether cooperation can emerge and persist despite an intrinsically fragmented population structure. Here we develop a framework to study the evolution of cooperation on temporal networks. Surprisingly, we find that network temporality actually enhances the evolution of cooperation relative to comparable static networks, despite the fact that bursty interaction patterns generally impede cooperation. We resolve this tension by proposing a measure to quantify the amount of temporality in a network, revealing an intermediate level that maximally boosts cooperation. Our results open a new avenue for investigating the evolution of cooperation and other emergent behaviours in more realistic structured populations.

Population structure enables emergence of cooperation among individuals, but the impact of the dynamic nature of real interaction networks is not understood. Here, the authors study the evolution of cooperation on temporal networks and find that temporality enhances the evolution of cooperation.

Turbulent coherent structures and early life below the Kolmogorov scale

Major evolutionary transitions, including the emergence of life, likely occurred in aqueous environments. While the role of water’s chemistry in early life is well studied, the effects of water’s ability to manipulate population structure are less clear. Population structure is known to be critical, as effective replicators must be insulated from parasites. Here, we propose that turbulent coherent structures, long-lasting flow patterns which trap particles, may serve many of the properties associated with compartments — collocalization, division, and merging — which are commonly thought to play a key role in the origins of life and other evolutionary transitions. We substantiate this idea by simulating multiple proposed metabolisms for early life in a simple model of a turbulent flow, and find that balancing the turnover times of biological particles and coherent structures can indeed enhance the likelihood of these metabolisms overcoming extinction either via parasitism or via a lack of metabolic support. Our results suggest that group selection models may be applicable with fewer physical and chemical constraints than previously thought, and apply much more widely in aqueous environments.

Models of the origin of life generally require a mechanism to structure emerging populations. Here, Krieger et al. develop spatial models showing that coherent structures arising in turbulent flows in aquatic environments could have provided compartments that facilitated the origin of life.

Motion, fixation probability and the choice of an evolutionary process

Seemingly minor details of mathematical and computational models of evolution are known to change the effect of population structure on the outcome of evolutionary processes. For example, birth-death dynamics often result in amplification of selection, while death-birth processes have been associated with suppression. In many biological populations the interaction structure is not static. Instead, members of the population are in motion and can interact with different individuals at different times. In this work we study populations embedded in a flowing medium; the interaction network is then time dependent. We use computer simulations to investigate how this dynamic structure affects the success of invading mutants, and compare these effects for different coupled birth and death processes. Specifically, we show how the speed of the motion impacts the fixation probability of an invading mutant. Flows of different speeds interpolate between evolutionary dynamics on fixed heterogeneous graphs and well-stirred populations; this allows us to systematically compare against known results for static structured populations. We find that motion has an active role in amplifying or suppressing selection by fragmenting and reconnecting the interaction graph. While increasing flow speeds suppress selection for most evolutionary models, we identify characteristic responses to flow for the different update rules we test. In particular we find that selection can be maximally enhanced or suppressed at intermediate flow speeds.

Whether a mutation spreads in a population or not is one of the most important questions in biology. The evolution of cancer and antibiotic resistance, for example, are mediated by invading mutants. Recent work has shown that population structure can have important consequences for the outcome of evolution. For instance, a mutant can have a higher or a lower chance of invasion than in unstructured populations. These effects can depend on seemingly minor details of the evolutionary model, such as the order of birth and death events. Many biological populations are in motion, for example due to external stirring. Experimentally this is known to be important; the performance of mutants in E. coli populations, for example, depends on the rate of mixing. Here, we focus on simulations of populations in a flowing medium, and compare the success of a mutant for different flow speeds. We contrast different evolutionary models, and identify what features of the evolutionary model affect mutant success for different speeds of the flow. We find that the chance of mutant invasion can be at its highest (or lowest) at intermediate flow speeds, depending on the order in which birth and death events occur in the evolutionary process.

Fixation probabilities in weakly compressible fluid flows

Competition between biological species in marine environments is affected by the motion of the surrounding fluid. An effective 2D compressibility can arise, for example, from the convergence and divergence of water masses at the depth at which passively traveling photosynthetic organisms are restricted to live. In this report, we seek to quantitatively study genetics under flow. To this end, we couple an off-lattice agent-based simulation of two populations in 1D to a weakly compressible velocity field-first a sine wave and then a shell model of turbulence. We find for both cases that even in a regime where the overall population structure is approximately unaltered, the flow can significantly diminish the effect of a selective advantage on fixation probabilities. We understand this effect in terms of the enhanced survival of organisms born at sources in the flow and the influence of Fisher genetic waves.