Showing 1–48 of 48 results for author: Nelson, R

Turbulent mixing controls fixation of growing antagonistic populations

Authors: Jonathan Bauermann, Roberto Benzi, David R. Nelson, Suraj Shankar, Federico Toschi

Abstract: Unlike coffee and cream that homogenize when stirred, growing micro-organisms (e.g., bacteria, baker's yeast) can actively kill each other and avoid mixing. How do such antagonistic interactions impact the growth and survival of competing strains, while being spatially advected by turbulent flows? By using numerical simulations of a continuum model, we study the dynamics of two antagonistic strain… ▽ More Unlike coffee and cream that homogenize when stirred, growing micro-organisms (e.g., bacteria, baker's yeast) can actively kill each other and avoid mixing. How do such antagonistic interactions impact the growth and survival of competing strains, while being spatially advected by turbulent flows? By using numerical simulations of a continuum model, we study the dynamics of two antagonistic strains that are dispersed by incompressible turbulent flows in two spatial dimensions. A key parameter is the ratio of the fluid transport time to that of biological reproduction, which determines the winning strain that ultimately takes over the whole population from an initial heterogeneous state. By quantifying the probability and mean time for fixation along with the spatial structure of concentration fluctuations, we demonstrate how turbulence raises the threshold for biological nucleation and antagonism suppresses flow-induced mixing by depleting the population at interfaces. Our work highlights the unusual biological consequences of the interplay of turbulent fluid flows with antagonistic population dynamics, with potential implications for marine microbial ecology and origins of biological chirality. △ Less

Submitted 21 August, 2024; originally announced August 2024.

Comments: 7 pages, 4 figures

MiNAA: Microbiome Network Alignment Algorithm

Authors: Reed Nelson, Rosa Aghdam, Claudia Solis-Lemus

Abstract: Our Microbiome Network Alignment Algorithm (MiNAA) aligns two microbial networks using a combination of the GRAph ALigner (GRAAL) algorithm and the Hungarian algorithm. Network alignment algorithms find pairs of nodes (one node from the first network and the other node from the second network) that have the highest similarity. Traditionally, similarity has been defined as topological similarity su… ▽ More Our Microbiome Network Alignment Algorithm (MiNAA) aligns two microbial networks using a combination of the GRAph ALigner (GRAAL) algorithm and the Hungarian algorithm. Network alignment algorithms find pairs of nodes (one node from the first network and the other node from the second network) that have the highest similarity. Traditionally, similarity has been defined as topological similarity such that the neighborhoods around the two nodes are similar. Recent implementations of network alignment methods such as NETAL and L-GRAAL also include measures of biological similarity, yet these methods are restricted to one specific type of biological similarity (e.g. sequence similarity in L-GRAAL). Our work extends existing network alignment implementations by allowing any type of biological similarity to be input by the user. This flexibility allows the user to choose whatever measure of biological similarity is suitable for the study at hand. In addition, unlike most existing network alignment methods that are tailored for protein or gene interaction networks, our work is the first one suited for microbiome networks. △ Less

Submitted 12 December, 2022; originally announced December 2022.

Recommendations for extending the GFF3 specification for improved interoperability of genomic data

Authors: Surya Saha, Scott Cain, Ethalinda K. S. Cannon, Nathan Dunn, Andrew Farmer, Zhi-Liang Hu, Gareth Maslen, Sierra Moxon, Christopher J Mungall, Rex Nelson, Monica F. Poelchau

Abstract: The GFF3 format is a common, flexible tab-delimited format representing the structure and function of genes or other mapped features (https://github.com/The-Sequence-Ontology/Specifications/blob/master/gff3.md). However, with increasing re-use of annotation data, this flexibility has become an obstacle for standardized downstream processing. Common software packages that export annotations in GFF3… ▽ More The GFF3 format is a common, flexible tab-delimited format representing the structure and function of genes or other mapped features (https://github.com/The-Sequence-Ontology/Specifications/blob/master/gff3.md). However, with increasing re-use of annotation data, this flexibility has become an obstacle for standardized downstream processing. Common software packages that export annotations in GFF3 format model the same data and metadata in different notations, which puts the burden on end-users to interpret the data model. The AgBioData consortium is a group of genomics, genetics and breeding databases and partners working towards shared practices and standards. Providing concrete guidelines for generating GFF3, and creating a standard representation of the most common biological data types would provide a major increase in efficiency for AgBioData databases and the genomics research community that use the GFF3 format in their daily operations. The AgBioData GFF3 working group has developed recommendations to solve common problems in the GFF3 format. We suggest improvements for each of the GFF3 fields, as well as the special cases of modeling functional annotations, and standard protein-coding genes. We welcome further discussion of these recommendations. We request the genomics and bioinformatics community to utilize the github repository (https://github.com/NAL-i5K/AgBioData_GFF3_recommendation) to provide feedback via issues or pull requests. △ Less

Submitted 15 February, 2022; originally announced February 2022.

Comments: 23 pages

Spatial population genetics with fluid flow

Authors: Roberto Benzi, David R. Nelson, Suraj Shankar, Federico Toschi, Xiaojue Zhu

Abstract: The growth and evolution of microbial populations is often subjected to advection by fluid flows in spatially extended environments, with immediate consequences for questions of spatial population genetics in marine ecology, planktonic diversity and origin of life scenarios. Here, we review recent progress made in understanding this rich problem in the simplified setting of two competing genetic m… ▽ More The growth and evolution of microbial populations is often subjected to advection by fluid flows in spatially extended environments, with immediate consequences for questions of spatial population genetics in marine ecology, planktonic diversity and origin of life scenarios. Here, we review recent progress made in understanding this rich problem in the simplified setting of two competing genetic microbial strains subjected to fluid flows. As a pedagogical example we focus on antagonsim, i.e., two killer microorganism strains, each secreting toxins that impede the growth of their competitors (competitive exclusion), in the presence of stationary fluid flows. By solving two coupled reaction-diffusion equations that include advection by simple steady cellular flows composed of characteristic flow motifs in two dimensions (2d), we show how local flow shear and compressibility effects can interact with selective advantage to have a dramatic influence on genetic competition and fixation in spatially distributed populations. We analyze several 1d and 2d flow geometries including sources, sinks, vortices and saddles, and show how simple analytical models of the dynamics of the genetic interface can be used to shed light on the nucleation, coexistence and flow-driven instabilities of genetic drops. By exploiting an analogy with phase separation with nonconserved order parameters, we uncover how these genetic drops harness fluid flows for novel evolutionary strategies, even in the presence of number fluctuations, as confirmed by agent-based simulations as well. △ Less

Submitted 30 June, 2022; v1 submitted 16 December, 2021; originally announced December 2021.

Comments: 29 pages, 22 figures

Journal ref: Rep. Prog. Phys. 85 096601, 2022

The collective effect of finite-sized inhomogeneities on the spatial spread of populations in two dimensions

Authors: Wolfram Möbius, Francesca Tesser, Kim M. J. Alards, Roberto Benzi, David R. Nelson, Federico Toschi

Abstract: The dynamics of a population expanding into unoccupied habitat has been primarily studied for situations in which growth and dispersal parameters are uniform in space or vary in one dimension. Here we study the influence of finite-sized individual inhomogeneities and their collective effect on front speed if randomly placed in a two-dimensional habitat. We use an individual-based model to investig… ▽ More The dynamics of a population expanding into unoccupied habitat has been primarily studied for situations in which growth and dispersal parameters are uniform in space or vary in one dimension. Here we study the influence of finite-sized individual inhomogeneities and their collective effect on front speed if randomly placed in a two-dimensional habitat. We use an individual-based model to investigate the front dynamics for a region in which dispersal or growth of individuals is reduced to zero (obstacles) or increased above the background (hotspots), respectively. In a regime where front dynamics is determined by a local front speed only, a principle of least time can be employed to predict front speed and shape. The resulting analytical solutions motivate an event-based algorithm illustrating the effects of several obstacles or hotspots. We finally apply the principle of least time to large heterogeneous environments by solving the Eikonal equation numerically. Obstacles lead to a slow-down that is dominated by the number density and width of obstacles, but not by their precise shape. Hotspots result in a speedup, which we characterise as function of hotspot strength and density. Our findings emphasise the importance of taking the dimensionality of the environment into account. △ Less

Submitted 20 September, 2021; v1 submitted 11 October, 2019; originally announced October 2019.

Comments: Authors Wolfram Möbius and Francesca Tesser contributed equally

Nucleation of antagonistic organisms and cellular competitions on curved, inflating substrates

Authors: Maxim O. Lavrentovich, David R. Nelson

Abstract: We consider the dynamics of spatially-distributed, diffusing populations of organisms with antagonistic interactions. These interactions are found on many length scales, ranging from kilometer-scale animal range dynamics with selection against hybrids to micron-scale interactions between poison-secreting microbial populations. We find that the dynamical line tension at the interface between antago… ▽ More We consider the dynamics of spatially-distributed, diffusing populations of organisms with antagonistic interactions. These interactions are found on many length scales, ranging from kilometer-scale animal range dynamics with selection against hybrids to micron-scale interactions between poison-secreting microbial populations. We find that the dynamical line tension at the interface between antagonistic organisms suppresses survival probabilities of small clonal clusters: the line tension introduces a critical cluster size that an organism with a selective advantage must achieve before deterministically spreading through the population. We calculate the survival probability as a function of selective advantage $δ$ and antagonistic interaction strength $σ$. Unlike a simple Darwinian selective advantage, the survival probability depends strongly on the spatial diffusion constant $D_s$ of the strains when $σ>0$, with suppressed survival when both species are more motile. Finally, we study the survival probability of a single mutant cell at the frontier of a growing spherical cluster of cells, such as the surface of an avascular spherical tumor. Both the inflation and curvature of the frontier significantly enhance the survival probability by changing the critical size of the nucleating cell cluster. △ Less

Submitted 23 October, 2019; v1 submitted 17 July, 2019; originally announced July 2019.

Comments: 15 pages, 10 figures, final version

Journal ref: Phys. Rev. E 100, 042406 (2019)

Non-Hermitian Quasi-Localization and Ring Attractor Neural Networks

Authors: Hidenori Tanaka, David R. Nelson

Abstract: Eigenmodes of a broad class of "sparse" random matrices, with interactions concentrated near the diagonal, exponentially localize in space, as initially discovered in 1957 by Anderson for quantum systems. Anderson localization plays ubiquitous roles in varieties of problems from electrons in solids to mechanical and optical systems. However, its implications in neuroscience (where the connections… ▽ More Eigenmodes of a broad class of "sparse" random matrices, with interactions concentrated near the diagonal, exponentially localize in space, as initially discovered in 1957 by Anderson for quantum systems. Anderson localization plays ubiquitous roles in varieties of problems from electrons in solids to mechanical and optical systems. However, its implications in neuroscience (where the connections can be strongly asymmetric) have been largely unexplored, mainly because synaptic connectivity matrices of neural systems are often "dense", which makes the eigenmodes spatially extended. Here, we explore roles that Anderson localization could be playing in neural networks by focusing on "spatially structured" disorder in synaptic connectivity matrices. Recently, neuroscientists have experimentally confirmed that the local excitation and global inhibition (LEGI) ring attractor model can functionally represent head direction cells in Drosophila melanogaster central brain. We first study a non-Hermitian (i.e. asymmetric) tight-binding model with disorder and then establish a connection to the LEGI ring attractor model. We discover that (i) Principal eigenvectors of the LEGI ring attractor networks with structured nearest neighbor disorder are "quasi-localized", even with fully dense inhibitory connections. (ii) The quasi-localized eigenvectors play dominant roles in the early time neural dynamics, and the location of the principal quasi-localized eigenvectors predict an initial location of the "bump of activity" representing, say, a head direction of an insect. Our investigations open up a new venue for explorations at the intersection between the theory of Anderson localization and neural networks with spatially structured disorder. △ Less

Submitted 18 November, 2018; originally announced November 2018.

Journal ref: Phys. Rev. E 99, 062406 (2019)

Fixation Probabilities in Weakly Compressible Fluid Flows

Authors: Abigail Plummer, Roberto Benzi, David R. Nelson, Federico Toschi

Abstract: 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 co… ▽ More 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. △ Less

Submitted 23 January, 2019; v1 submitted 21 August, 2018; originally announced August 2018.

Comments: 8 pages, 6 figures; SI: 9 pages, 3 figures

Journal ref: PNAS, vol. 116 no. 2, 373-378 (2019)

Spatial gene drives and pushed genetic waves

Authors: Hidenori Tanaka, Howard A. Stone, David R. Nelson

Abstract: Gene drives have the potential to rapidly replace a harmful wild-type allele with a gene drive allele engineered to have desired functionalities. However, an accidental or premature release of a gene drive construct to the natural environment could damage an ecosystem irreversibly. Thus, it is important to understand the spatiotemporal consequences of the super-Mendelian population genetics prior… ▽ More Gene drives have the potential to rapidly replace a harmful wild-type allele with a gene drive allele engineered to have desired functionalities. However, an accidental or premature release of a gene drive construct to the natural environment could damage an ecosystem irreversibly. Thus, it is important to understand the spatiotemporal consequences of the super-Mendelian population genetics prior to potential applications. Here, we employ a reaction-diffusion model for sexually reproducing diploid organisms to study how a locally introduced gene drive allele spreads to replace the wild-type allele, even though it possesses a selective disadvantage $s>0$. Using methods developed by N. Barton and collaborators, we show that socially responsible gene drives require $0.5<s<0.697$, a rather narrow range. In this "pushed wave" regime, the spatial spreading of gene drives will be initiated only when the initial frequency distribution is above a threshold profile called "critical propagule", which acts as a safeguard against accidental release. We also study how the spatial spread of the pushed wave can be stopped by making gene drives uniquely vulnerable ("sensitizing drive") in a way that is harmless for a wild-type allele. Finally, we show that appropriately sensitized drives in two dimensions can be stopped even by imperfect barriers perforated by a series of gaps. △ Less

Submitted 21 August, 2017; v1 submitted 11 April, 2017; originally announced April 2017.

Journal ref: PNAS, vol. 114 no. 32, 8452-8457 (2017)

Population Genetics with Fluctuating Population Sizes

Authors: Thiparat Chotibut, David R. Nelson

Abstract: Standard neutral population genetics theory with a strictly fixed population size has important limitations. An alternative model that allows independently fluctuating population sizes and reproduces the standard neutral evolution is reviewed. We then study a situation such that the competing species are neutral at the equilibrium population size but population size fluctuations nevertheless favor… ▽ More Standard neutral population genetics theory with a strictly fixed population size has important limitations. An alternative model that allows independently fluctuating population sizes and reproduces the standard neutral evolution is reviewed. We then study a situation such that the competing species are neutral at the equilibrium population size but population size fluctuations nevertheless favor fixation of one species over the other. In this case, a separation of timescales emerges naturally and allows adiabatic elimination of a fast population size variable to deduce the fluctuations-induced selection dynamics near the equilibrium population size. The results highlight the incompleteness of the standard population genetics with a strictly fixed population size. △ Less

Submitted 29 August, 2016; originally announced August 2016.

Comments: Submitted to Journal of Statistical Physics, Special Issue: Dedicated to the Memory of Leo Kadanoff. arXiv admin note: text overlap with arXiv:1412.6688

Striated Populations in Disordered Environments with Advection

Authors: Thiparat Chotibut, David R. Nelson, Sauro Succi

Abstract: Growth in static and controlled environments such as a Petri dish can be used to study the spatial population dynamics of microorganisms. However, natural populations such as marine microbes experience fluid advection and often grow up in heterogeneous environments. We investigate a generalized Fisher-Kolmogorov-Petrovsky-Piscounov (FKPP) equation describing single species population subject to a… ▽ More Growth in static and controlled environments such as a Petri dish can be used to study the spatial population dynamics of microorganisms. However, natural populations such as marine microbes experience fluid advection and often grow up in heterogeneous environments. We investigate a generalized Fisher-Kolmogorov-Petrovsky-Piscounov (FKPP) equation describing single species population subject to a constant flow field and quenched random spatially inhomogeneous growth rates with a fertile overall growth condition. We analytically and numerically demonstrate that the non-equilibrium steady-state population density develops a flow-driven striation pattern. The striations are highly asymmetric with a longitudinal correlation length that diverges linearly with the flow speed and a transverse correlation length that approaches a finite velocity-independent value. Linear response theory is developed to study the statistics of the steady states. Theoretical predictions show excellent agreement with the numerical steady states of the generalized FKPP equation obtained from Lattice Boltzmann simulations. These findings suggest that, although the growth disorder can be spatially uncorrelated, correlated population structures with striations emerge naturally at sufficiently strong advection. △ Less

Submitted 29 March, 2016; originally announced March 2016.

Comments: 30 pages, 8 figures

Non-Hermitian Localization in Biological Networks

Authors: Ariel Amir, Naomichi Hatano, David R. Nelson

Abstract: We explore the spectra and localization properties of the N-site banded one-dimensional non-Hermitian random matrices that arise naturally in sparse neural networks. Approximately equal numbers of random excitatory and inhibitory connections lead to spatially localized eigenfunctions, and an intricate eigenvalue spectrum in the complex plane that controls the spontaneous activity and induced respo… ▽ More We explore the spectra and localization properties of the N-site banded one-dimensional non-Hermitian random matrices that arise naturally in sparse neural networks. Approximately equal numbers of random excitatory and inhibitory connections lead to spatially localized eigenfunctions, and an intricate eigenvalue spectrum in the complex plane that controls the spontaneous activity and induced response. A finite fraction of the eigenvalues condense onto the real or imaginary axes. For large N, the spectrum has remarkable symmetries not only with respect to reflections across the real and imaginary axes, but also with respect to 90 degree rotations, with an unusual anisotropic divergence in the localization length near the origin. When chains with periodic boundary conditions become directed, with a systematic directional bias superimposed on the randomness, a hole centered on the origin opens up in the density-of-states in the complex plane. All states are extended on the rim of this hole, while the localized eigenvalues outside the hole are unchanged. The bias dependent shape of this hole tracks the bias independent contours of constant localization length. We treat the large-N limit by a combination of direct numerical diagonalization and using transfer matrices, an approach that allows us to exploit an electrostatic analogy connecting the "charges" embodied in the eigenvalue distribution with the contours of constant localization length. We show that similar results are obtained for more realistic neural networks that obey "Dale's Law" (each site is purely excitatory or inhibitory), and conclude with perturbation theory results that describe the limit of large bias g, when all states are extended. Related problems arise in random ecological networks and in chains of artificial cells with randomly coupled gene expression patterns. △ Less

Submitted 4 January, 2016; v1 submitted 17 December, 2015; originally announced December 2015.

Journal ref: Phys. Rev. E 93, 042310 (2016)

Evolutionary Dynamics with Fluctuating Population Sizes and Strong Mutualism

Authors: Thiparat Chotibut, David R. Nelson

Abstract: Game theory ideas provide a useful framework for studying evolutionary dynamics in a well-mixed environment. This approach, however, typically enforces a strictly fixed overall population size, deemphasizing natural growth processes. We study a competitive Lotka-Volterra model, with number fluctuations, that accounts for natural population growth and encompasses interaction scenarios typical of ev… ▽ More Game theory ideas provide a useful framework for studying evolutionary dynamics in a well-mixed environment. This approach, however, typically enforces a strictly fixed overall population size, deemphasizing natural growth processes. We study a competitive Lotka-Volterra model, with number fluctuations, that accounts for natural population growth and encompasses interaction scenarios typical of evolutionary games. We show that, in an appropriate limit, the model describes standard evolutionary games with both genetic drift and overall population size fluctuations. However, there are also regimes where a varying population size can strongly influence the evolutionary dynamics. We focus on the strong mutualism scenario and demonstrate that standard evolutionary game theory fails to describe our simulation results. We then analytically and numerically determine fixation probabilities as well as mean fixation times using matched asymptotic expansions, taking into account the population size degree of freedom. These results elucidate the interplay between population dynamics and evolutionary dynamics in well-mixed systems. △ Less

Submitted 29 August, 2016; v1 submitted 20 December, 2014; originally announced December 2014.

Comments: Updated Version, Published in Phys. Rev. E 92, 022718 on 20 August 2015

Journal ref: Phys. Rev. E 92, 022718 (2015)

Survival Probabilities at Spherical Frontiers

Authors: Maxim O. Lavrentovich, David R. Nelson

Abstract: Motivated by tumor growth and spatial population genetics, we study the interplay between evolutionary and spatial dynamics at the surfaces of three-dimensional, spherical range expansions. We consider range expansion radii that grow with an arbitrary power-law in time: $R(t)=R_0(1+t/t^*)^Θ$, where $Θ$ is a growth exponent, $R_0$ is the initial radius, and $t^*$ is a characteristic time for the gr… ▽ More Motivated by tumor growth and spatial population genetics, we study the interplay between evolutionary and spatial dynamics at the surfaces of three-dimensional, spherical range expansions. We consider range expansion radii that grow with an arbitrary power-law in time: $R(t)=R_0(1+t/t^*)^Θ$, where $Θ$ is a growth exponent, $R_0$ is the initial radius, and $t^*$ is a characteristic time for the growth, to be affected by the inflating geometry. We vary the parameters $t^*$ and $Θ$ to capture a variety of possible growth regimes. Guided by recent results for two-dimensional inflating range expansions, we identify key dimensionless parameters that describe the survival probability of a mutant cell with a small selective advantage arising at the population frontier. Using analytical techniques, we calculate this probability for arbitrary $Θ$. We compare our results to simulations of linearly inflating expansions ($Θ=1$ spherical Fisher-Kolmogorov-Petrovsky-Piscunov waves) and treadmilling populations ($Θ=0$, with cells in the interior removed by apoptosis or a similar process). We find that mutations at linearly inflating fronts have survival probabilities enhanced by factors of 100 or more relative to mutations at treadmilling population frontiers. We also discuss the special properties of "marginally inflating" $(Θ=1/2)$ expansions. △ Less

Submitted 1 June, 2015; v1 submitted 26 August, 2014; originally announced August 2014.

Comments: 35 pages, 11 figures, revised version

Journal ref: Theor. Popul. Biol. 102 (2015) 26-39

Genetic Influences on Brain Gene Expression in Rats Selected for Tameness and Aggression

Authors: Henrike O. Heyne, Susann Lautenschläger, Ronald Nelson, François Besnier, Maxime Rotival, Alexander Cagan, Rimma Kozhemyakina, Irina Z. Plyusnina, Lyudmila Trut, Örjan Carlborg, Enrico Petretto, Leonid Kruglyak, Svante Pääbo, Torsten Schöneberg, Frank W. Albert

Abstract: Inter-individual differences in many behaviors are partly due to genetic differences, but the identification of the genes and variants that influence behavior remains challenging. Here, we studied an F2 intercross of two outbred lines of rats selected for tame and aggressive behavior towards humans for more than 64 generations. By using a mapping approach that is able to identify genetic loci segr… ▽ More Inter-individual differences in many behaviors are partly due to genetic differences, but the identification of the genes and variants that influence behavior remains challenging. Here, we studied an F2 intercross of two outbred lines of rats selected for tame and aggressive behavior towards humans for more than 64 generations. By using a mapping approach that is able to identify genetic loci segregating within the lines, we identified four times more loci influencing tameness and aggression than by an approach that assumes fixation of causative alleles, suggesting that many causative loci were not driven to fixation by the selection. We used RNA sequencing in 150 F2 animals to identify hundreds of loci that influence brain gene expression. Several of these loci colocalize with tameness loci and may reflect the same genetic variants. Through analyses of correlations between allele effects on behavior and gene expression, differential expression between the tame and aggressive rat selection lines, and correlations between gene expression and tameness in F2 animals, we identify the genes Gltscr2, Lgi4, Zfp40 and Slc17a7 as candidate contributors to the strikingly different behavior of the tame and aggressive animals. △ Less

Submitted 14 April, 2014; originally announced April 2014.

Genetic drift opposes mutualism during spatial population expansion

Authors: Melanie JI Muller, Beverly I Neugeboren, David R Nelson, Andrew W Murray

Abstract: Mutualistic interactions benefit both partners, promoting coexistence and genetic diversity. Spatial structure can promote cooperation, but spatial expansions may also make it hard for mutualistic partners to stay together, since genetic drift at the expansion front creates regions of low genetic and species diversity. To explore the antagonism between mutualism and genetic drift, we grew cross-fe… ▽ More Mutualistic interactions benefit both partners, promoting coexistence and genetic diversity. Spatial structure can promote cooperation, but spatial expansions may also make it hard for mutualistic partners to stay together, since genetic drift at the expansion front creates regions of low genetic and species diversity. To explore the antagonism between mutualism and genetic drift, we grew cross-feeding strains of the budding yeast S. cerevisiae on agar surfaces as a model for mutualists undergoing spatial expansions. By supplying varying amounts of the exchanged nutrients, we tuned strength and symmetry of the mutualistic interaction. Strong mutualism suppresses genetic demixing during spatial expansions and thereby maintains diversity, but weak or asymmetric mutualism is overwhelmed by genetic drift even when mutualism is still beneficial, slowing growth and reducing diversity. Theoretical modeling using experimentally measured parameters predicts the size of demixed regions and how strong mutualism must be to survive a spatial expansion. △ Less

Submitted 24 February, 2014; originally announced February 2014.

Comments: 12 pages main text, 20 pages SI

Journal ref: Proc. Natl. Acad. Sci. USA 111, 1037-1042 (2014)

Asymmetric Mutualism in Two- and Three-Dimensional Range Expansions

Authors: Maxim O. Lavrentovich, David R. Nelson

Abstract: Genetic drift at the frontiers of two-dimensional range expansions of microorganisms can frustrate local cooperation between different genetic variants, demixing the population into distinct sectors. In a biological context, mutualistic or antagonistic interactions will typically be asymmetric between variants. By taking into account both the asymmetry and the interaction strength, we show that th… ▽ More Genetic drift at the frontiers of two-dimensional range expansions of microorganisms can frustrate local cooperation between different genetic variants, demixing the population into distinct sectors. In a biological context, mutualistic or antagonistic interactions will typically be asymmetric between variants. By taking into account both the asymmetry and the interaction strength, we show that the much weaker demixing in three dimensions allows for a mutualistic phase over a much wider range of asymmetric cooperative benefits, with mutualism prevailing for any positive, symmetric benefit. We also demonstrate that expansions with undulating fronts roughen dramatically at the boundaries of the mutualistic phase, with severe consequences for the population genetics along the transition lines. △ Less

Submitted 5 April, 2014; v1 submitted 1 September, 2013; originally announced September 2013.

Comments: 9 pages, 13 figures, this version published in Phys. Rev. Lett

Journal ref: Phys. Rev. Lett. 112, 138102 (2014)

Bending forces plastically deform growing bacterial cell walls

Authors: Ariel Amir, Farinaz Babaeipour, Dustin B. McIntosh, David R. Nelson, Suckjoon Jun

Abstract: Cell walls define a cell shape in bacteria. They are rigid to resist large internal pressures, but remarkably plastic to adapt to a wide range of external forces and geometric constraints. Currently, it is unknown how bacteria maintain their shape. In this work, we develop experimental and theoretical approaches and show that mechanical stresses regulate bacterial cell-wall growth. By applying a p… ▽ More Cell walls define a cell shape in bacteria. They are rigid to resist large internal pressures, but remarkably plastic to adapt to a wide range of external forces and geometric constraints. Currently, it is unknown how bacteria maintain their shape. In this work, we develop experimental and theoretical approaches and show that mechanical stresses regulate bacterial cell-wall growth. By applying a precisely controllable hydrodynamic force to growing rod-shaped Escherichia coli and Bacillus subtilis cells, we demonstrate that the cells can exhibit two fundamentally different modes of deformation. The cells behave like elastic rods when subjected to transient forces, but deform plastically when significant cell wall synthesis occurs while the force is applied. The deformed cells always recover their shape. The experimental results are in quantitative agreement with the predictions of the theory of dislocation-mediated growth. In particular, we find that a single dimensionless parameter, which depends on a combination of independently measured physical properties of the cell, can describe the cell's responses under various experimental conditions. These findings provide insight into how living cells robustly maintain their shape under varying physical environments. △ Less

Submitted 30 April, 2014; v1 submitted 24 May, 2013; originally announced May 2013.

Journal ref: PNAS 111 16 5778-5783 (2014)

Nutrient Shielding in Clusters of Cells

Authors: Maxim O. Lavrentovich, John H. Koschwanez, David R. Nelson

Abstract: Cellular nutrient consumption is influenced by both the nutrient uptake kinetics of an individual cell and the cells' spatial arrangement. Large cell clusters or colonies have inhibited growth at the cluster's center due to the shielding of nutrients by the cells closer to the surface. We develop an effective medium theory that predicts a thickness $\ell$ of the outer shell of cells in the cluster… ▽ More Cellular nutrient consumption is influenced by both the nutrient uptake kinetics of an individual cell and the cells' spatial arrangement. Large cell clusters or colonies have inhibited growth at the cluster's center due to the shielding of nutrients by the cells closer to the surface. We develop an effective medium theory that predicts a thickness $\ell$ of the outer shell of cells in the cluster that receives enough nutrient to grow. The cells are treated as partially absorbing identical spherical nutrient sinks, and we identify a dimensionless parameter $ν$ that characterizes the absorption strength of each cell. The parameter $ν$ can vary over many orders of magnitude between different cell types, ranging from bacteria and yeast to human tissue. The thickness $\ell$ decreases with increasing $ν$, increasing cell volume fraction $φ$, and decreasing ambient nutrient concentration $ψ_{\infty}$. The theoretical results are compared with numerical simulations and experiments. In the latter studies, colonies of budding yeast, \textit{Saccharomyces cerevisiae}, are grown on glucose media and imaged under a confocal microscope. We measure the growth inside the colonies via a fluorescent protein reporter and compare the experimental and theoretical results for the thickness $\ell$. △ Less

Submitted 13 June, 2013; v1 submitted 23 April, 2013; originally announced April 2013.

Comments: 19 pages, 15 figures

Journal ref: Phys. Rev. E 87, 062703 (2013)

Radial Domany-Kinzel Models with Mutation and Selection

Authors: Maxim O. Lavrentovich, Kirill S. Korolev, David R. Nelson

Abstract: We study the effect of spatial structure, genetic drift, mutation, and selective pressure on the evolutionary dynamics in a simplified model of asexual organisms colonizing a new territory. Under an appropriate coarse-graining, the evolutionary dynamics is related to the directed percolation processes that arise in voter models, the Domany-Kinzel (DK) model, contact process, etc. We explore the di… ▽ More We study the effect of spatial structure, genetic drift, mutation, and selective pressure on the evolutionary dynamics in a simplified model of asexual organisms colonizing a new territory. Under an appropriate coarse-graining, the evolutionary dynamics is related to the directed percolation processes that arise in voter models, the Domany-Kinzel (DK) model, contact process, etc. We explore the differences between linear (flat front) expansions and the much less familiar radial (curved front) range expansions. For the radial expansion, we develop a generalized, off-lattice DK model that minimizes otherwise persistent lattice artifacts. With both simulations and analytical techniques, we study the survival probability of advantageous mutants, the spatial correlations between domains of neutral strains, and the dynamics of populations with deleterious mutations. "Inflation" at the frontier leads to striking differences between radial and linear expansions. For a colony with initial radius $R_0$ expanding at velocity $v$, significant genetic demixing, caused by local genetic drift, only occurs up to a finite time $t^* = R_0/v$, after which portions of the colony become causally disconnected due to the inflating perimeter of the expanding front. As a result, the effect of a selective advantage is amplified relative to genetic drift, increasing the survival probability of advantageous mutants. Inflation also modifies the underlying directed percolation transition, introducing novel scaling functions and modifications similar to a finite size effect. Finally, we consider radial range expansions with deflating perimeters, as might arise from colonization initiated along the shores of an island. △ Less

Submitted 14 January, 2013; v1 submitted 18 September, 2012; originally announced September 2012.

Comments: 24 pages, 21 figures

Journal ref: Phys. Rev. E 87, 012103 (2013)

Growth, competition and cooperation in spatial population genetics

Authors: Simone Pigolotti, Roberto Benzi, Prasad Perlekar Mogens H. Jensen, Federico Toschi, David R. Nelson

Abstract: We study an individual based model describing competition in space between two different alleles. Although the model is similar in spirit to classic models of spatial population genetics such as the stepping stone model, here however space is continuous and the total density of competing individuals fluctuates due to demographic stochasticity. By means of analytics and numerical simulations, we st… ▽ More We study an individual based model describing competition in space between two different alleles. Although the model is similar in spirit to classic models of spatial population genetics such as the stepping stone model, here however space is continuous and the total density of competing individuals fluctuates due to demographic stochasticity. By means of analytics and numerical simulations, we study the behavior of fixation probabilities, fixation times, and heterozygosity, in a neutral setting and in cases where the two species can compete or cooperate. By concluding with examples in which individuals are transported by fluid flows, we argue that this model is a natural choice to describe competition in marine environments. △ Less

Submitted 18 December, 2012; v1 submitted 24 August, 2012; originally announced August 2012.

Comments: 29 pages, 14 figures; revised version including a section with results in the presence of fluid flows

Journal ref: Theoretical Population Biology 84, 72-86 (2013)

Dislocation-mediated growth of bacterial cell walls

Authors: Ariel Amir, David R. Nelson

Abstract: Recent experiments have illuminated a remarkable growth mechanism of rod-shaped bacteria: proteins associated with cell wall extension move at constant velocity in circles oriented approximately along the cell circumference (Garner et al., Science (2011), Dominguez-Escobar et al. Science (2011), van Teeffelen et al. PNAS (2011). We view these as dislocations in the partially ordered peptidoglycan… ▽ More Recent experiments have illuminated a remarkable growth mechanism of rod-shaped bacteria: proteins associated with cell wall extension move at constant velocity in circles oriented approximately along the cell circumference (Garner et al., Science (2011), Dominguez-Escobar et al. Science (2011), van Teeffelen et al. PNAS (2011). We view these as dislocations in the partially ordered peptidoglycan structure, activated by glycan strand extension machinery, and study theoretically the dynamics of these interacting defects on the surface of a cylinder. Generation and motion of these interacting defects lead to surprising effects arising from the cylindrical geometry, with important implications for growth. We also discuss how long range elastic interactions and turgor pressure affect the dynamics of the fraction of actively moving dislocations in the bacterial cell wall. △ Less

Submitted 7 May, 2012; originally announced May 2012.

Comments: to appear in PNAS

Journal ref: PNAS 109, 25, 9833-9838 (2012)

Supplementary information for "Selective sweeps in growing microbial colonies"

Authors: Kirill S Korolev, Melanie J I Müller, Nilay Karahan, Andrew W Murray, Oskar Hallatschek, David R Nelson

Abstract: This is supplementary information for arXiv:1204.4896, which also appeared in Physical Biology 9, 026008 (2012). Main text abstract: Evolutionary experiments with microbes are a powerful tool to study mutations and natural selection. These experiments, however, are often limited to the well-mixed environments of a test tube or a chemostat. Since spatial organization can significantly affect evolut… ▽ More This is supplementary information for arXiv:1204.4896, which also appeared in Physical Biology 9, 026008 (2012). Main text abstract: Evolutionary experiments with microbes are a powerful tool to study mutations and natural selection. These experiments, however, are often limited to the well-mixed environments of a test tube or a chemostat. Since spatial organization can significantly affect evolutionary dynamics, the need is growing for evolutionary experiments in spatially structured environments. The surface of a Petri dish provides such an environment, but a more detailed understanding of microbial growth on Petri dishes is necessary to interpret such experiments. We formulate a simple deterministic reaction-diffusion model, which successfully predicts the spatial patterns created by two competing species during colony expansion. We also derive the shape of these patterns analytically without relying on microscopic details of the model. In particular, we find that the relative fitness of two microbial strains can be estimated from the logarithmic spirals created by selective sweeps. The theory is tested with strains of the budding yeast Saccharomyces cerevisiae, for spatial competitions with different initial conditions and for a range of relative fitnesses. The reaction-diffusion model also connects the microscopic parameters like growth rates and diffusion constants with macroscopic spatial patterns and predicts the relationship between fitness in liquid cultures and on Petri dishes, which we confirmed experimentally. Spatial sector patterns therefore provide an alternative fitness assay to the commonly used liquid culture fitness assays. △ Less

Submitted 27 April, 2012; originally announced April 2012.

Comments: The main text is available as arXiv:1204.4896

Journal ref: Supplementary information for Physical Biology 9, 026008 (2012)

Selective sweeps in growing microbial colonies

Authors: Kirill S Korolev, Melanie J I Müller, Nilay Karahan, Andrew W Murray, Oskar Hallatschek, David R Nelson

Abstract: Evolutionary experiments with microbes are a powerful tool to study mutations and natural selection. These experiments, however, are often limited to the well-mixed environments of a test tube or a chemostat. Since spatial organization can significantly affect evolutionary dynamics, the need is growing for evolutionary experiments in spatially structured environments. The surface of a Petri dish p… ▽ More Evolutionary experiments with microbes are a powerful tool to study mutations and natural selection. These experiments, however, are often limited to the well-mixed environments of a test tube or a chemostat. Since spatial organization can significantly affect evolutionary dynamics, the need is growing for evolutionary experiments in spatially structured environments. The surface of a Petri dish provides such an environment, but a more detailed understanding of microbial growth on Petri dishes is necessary to interpret such experiments. We formulate a simple deterministic reaction-diffusion model, which successfully predicts the spatial patterns created by two competing species during colony expansion. We also derive the shape of these patterns analytically without relying on microscopic details of the model. In particular, we find that the relative fitness of two microbial strains can be estimated from the logarithmic spirals created by selective sweeps. The theory is tested with strains of the budding yeast Saccharomyces cerevisiae, for spatial competitions with different initial conditions and for a range of relative fitnesses. The reaction-diffusion model also connects the microscopic parameters like growth rates and diffusion constants with macroscopic spatial patterns and predicts the relationship between fitness in liquid cultures and on Petri dishes, which we confirmed experimentally. Spatial sector patterns therefore provide an alternative fitness assay to the commonly used liquid culture fitness assays. △ Less

Submitted 1 May, 2012; v1 submitted 22 April, 2012; originally announced April 2012.

Comments: Supplementary information available at arXiv:1204.6328

Journal ref: Physical Biology 9, 026008 (2012)

Population dynamics in compressible flows

Authors: Roberto Benzi, Mogens H. Jensen, David R. Nelson, Prasad Perlekar, Simone Pigolotti, Federico Toschi

Abstract: Organisms often grow, migrate and compete in liquid environments, as well as on solid surfaces. However, relatively little is known about what happens when competing species are mixed and compressed by fluid turbulence. In these lectures we review our recent work on population dynamics and population genetics in compressible velocity fields of one and two dimensions. We discuss why compressible tu… ▽ More Organisms often grow, migrate and compete in liquid environments, as well as on solid surfaces. However, relatively little is known about what happens when competing species are mixed and compressed by fluid turbulence. In these lectures we review our recent work on population dynamics and population genetics in compressible velocity fields of one and two dimensions. We discuss why compressible turbulence is relevant for population dynamics in the ocean and we consider cases both where the velocity field is turbulent and when it is static. Furthermore, we investigate populations in terms of a continuos density field and when the populations are treated via discrete particles. In the last case we focus on the competition and fixation of one species compared to another △ Less

Submitted 28 March, 2012; originally announced March 2012.

Comments: 16 pages, talk delivered at the Geilo Winter School 2011

A Quantitative Test of Population Genetics Using Spatio-Genetic Patterns in Bacterial Colonies

Authors: Kirill S. Korolev, Joao B. Xavier, David R. Nelson, Kevin R. Foster

Abstract: It is widely accepted that population genetics theory is the cornerstone of evolutionary analyses. Empirical tests of the theory, however, are challenging because of the complex relationships between space, dispersal, and evolution. Critically, we lack quantitative validation of the spatial models of population genetics. Here we combine analytics, on and off-lattice simulations, and experiments wi… ▽ More It is widely accepted that population genetics theory is the cornerstone of evolutionary analyses. Empirical tests of the theory, however, are challenging because of the complex relationships between space, dispersal, and evolution. Critically, we lack quantitative validation of the spatial models of population genetics. Here we combine analytics, on and off-lattice simulations, and experiments with bacteria to perform quantitative tests of the theory. We study two bacterial species, the gut microbe Escherichia coli and the opportunistic pathogen Pseudomonas aeruginosa, and show that spatio-genetic patterns in colony biofilms of both species are accurately described by an extension of the one-dimensional stepping-stone model. We use one empirical measure, genetic diversity at the colony periphery, to parameterize our models and show that we can then accurately predict another key variable: the degree of short-range cell migration along an edge. Moreover, the model allows us to estimate other key parameters including effective population size (density) at the expansion frontier. While our experimental system is a simplification of natural microbial community, we argue it is a proof of principle that the spatial models of population genetics can quantitatively capture organismal evolution. △ Less

Submitted 24 October, 2011; originally announced October 2011.

Journal ref: The American Naturalist 178, 538 (2011)

Population genetics in compressible flows

Authors: Simone Pigolotti, Roberto Benzi, Mogens H. Jensen, David R. Nelson

Abstract: We study competition between two biological species advected by a compressible velocity field. Individuals are treated as discrete Lagrangian particles that reproduce or die in a density-dependent fashion. In the absence of a velocity field and fitness advantage, number fluctuations lead to a coarsening dynamics typical of the stochastic Fisher equation. We then study three examples of compressibl… ▽ More We study competition between two biological species advected by a compressible velocity field. Individuals are treated as discrete Lagrangian particles that reproduce or die in a density-dependent fashion. In the absence of a velocity field and fitness advantage, number fluctuations lead to a coarsening dynamics typical of the stochastic Fisher equation. We then study three examples of compressible advecting fields: a shell model of turbulence, a sinusoidal velocity field and a linear velocity sink. In all cases, advection leads to a striking drop in the fixation time, as well as a large reduction in the global carrying capacity. Despite localization on convergence zones, one species goes extinct much more rapidly than in well-mixed populations. For a weak harmonic potential, one finds a bimodal distribution of fixation times. The long-lived states in this case are demixed configurations with a single boundary, whose location depends on the fitness advantage. △ Less

Submitted 17 June, 2011; originally announced June 2011.

Comments: 10 pages, 5 figures, submitted

Journal ref: Phys. Rev. Lett. 108, 128102, 2012

Competition and cooperation in one-dimensional stepping stone models

Authors: K. S. Korolev, David R. Nelson

Abstract: Cooperative mutualism is a major force driving evolution and sustaining ecosystems. Although the importance of spatial degrees of freedom and number fluctuations is well-known, their effects on mutualism are not fully understood. With range expansions of microbes in mind, we show that, even when mutualism confers a distinct selective advantage, it persists only in populations with high density and… ▽ More Cooperative mutualism is a major force driving evolution and sustaining ecosystems. Although the importance of spatial degrees of freedom and number fluctuations is well-known, their effects on mutualism are not fully understood. With range expansions of microbes in mind, we show that, even when mutualism confers a distinct selective advantage, it persists only in populations with high density and frequent migrations. When these parameters are reduced, mutualism is generically lost via a directed percolation process, with a phase diagram strongly influenced by an exceptional DP2 transition. △ Less

Submitted 17 April, 2011; v1 submitted 13 April, 2011; originally announced April 2011.

Comments: 8 pages, 4 figures

Genetic Demixing and Evolutionary Forces in the One-Dimensional Stepping Stone Model

Authors: K. S. Korolev, Mikkel Avlund, Oskar Hallatschek, David R. Nelson

Abstract: We review and extend results for mutation, selection, genetic drift, and migration in a one-dimensional continuous population. The population is described by a continuous limit of the stepping stone model, which leads to the stochastic Fisher-Kolmogorov-Petrovsky-Piscounov equation with additional terms describing mutations. Although the stepping stone model was first proposed for population genet… ▽ More We review and extend results for mutation, selection, genetic drift, and migration in a one-dimensional continuous population. The population is described by a continuous limit of the stepping stone model, which leads to the stochastic Fisher-Kolmogorov-Petrovsky-Piscounov equation with additional terms describing mutations. Although the stepping stone model was first proposed for population genetics, it is closely related to "voter models" of interest in nonequilibrium statistical mechanics. The stepping stone model can also be regarded as an approximation to the dynamics of a thin layer of actively growing pioneers at the frontier of a colony of microorganisms undergoing a range expansion on a Petri dish. We find that the population tends to segregate into monoallelic domains. This segregation slows down genetic drift and selection because these two evolutionary forces can only act at the boundaries between the domains; the effects of mutation, however, are not significantly affected by the segregation. Although fixation in the neutral well-mixed (or "zero dimensional") model occurs exponentially in time, it occurs only algebraically fast in the one-dimensional model. If selection is weak, selective sweeps occur exponentially fast in both well-mixed and one-dimensional populations, but the time constants are different. We also find an unusual sublinear increase in the variance of the spatially averaged allele frequency with time. Although we focus on two alleles or variants, q-allele Potts-like models of gene segregation are considered as well. We also investigate the effects of geometry at the frontier by considering growth of circular colonies. Our analytical results are checked with simulations, and could be tested against recent spatial experiments on range expansions off linear inoculations of Escherichia coli and Saccharomyces cerevisiae. △ Less

Submitted 13 April, 2011; v1 submitted 29 April, 2009; originally announced April 2009.

Comments: 29 pages, 20 figures; Reviews of Modern Physics, Volume 82, April-June 2010

Fisher equation with turbulence in one dimension

Authors: Roberto Benzi, David R. Nelson

Abstract: We investigate the dynamics of the Fisher equation for the spreading of micro-organisms in one dimension subject to both turbulent convection and diffusion. We show that for strong enough turbulence, bacteria, for example, track in a quasilocalized fashion (with remakably long persistance times) sinks in the turbulent field. An important consequence is a large reduction in the carrying capacity… ▽ More We investigate the dynamics of the Fisher equation for the spreading of micro-organisms in one dimension subject to both turbulent convection and diffusion. We show that for strong enough turbulence, bacteria, for example, track in a quasilocalized fashion (with remakably long persistance times) sinks in the turbulent field. An important consequence is a large reduction in the carrying capacity of the fluid medium. We determine analytically the regimes where this quasi-localized behavior occurs and test our predictions by numerical simulations. △ Less

Submitted 2 February, 2009; originally announced February 2009.

Comments: 15 pages

Genetic drift at expanding frontiers promotes gene segregation

Authors: Oskar Hallatschek, Pascal Hersen, Sharad Ramanathan, David R. Nelson

Abstract: Competition between random genetic drift and natural selection plays a central role in evolution: Whereas non-beneficial mutations often prevail in small populations by chance, mutations that sweep through large populations typically confer a selective advantage. Here, however, we observe chance effects during range expansions that dramatically alter the gene pool even in large microbial populat… ▽ More Competition between random genetic drift and natural selection plays a central role in evolution: Whereas non-beneficial mutations often prevail in small populations by chance, mutations that sweep through large populations typically confer a selective advantage. Here, however, we observe chance effects during range expansions that dramatically alter the gene pool even in large microbial populations. Initially well-mixed populations of two fluorescently labeled strains of Escherichia coli develop well-defined, sector-like regions with fractal boundaries in expanding colonies. The formation of these regions is driven by random fluctuations that originate in a thin band of pioneers at the expanding frontier. A comparison of bacterial and yeast colonies (Saccharomyces cerevisiae) suggests that this large-scale genetic sectoring is a generic phenomenon that may provide a detectable footprint of past range expansions. △ Less

Submitted 12 December, 2008; originally announced December 2008.

Comments: Please visit http://www.pnas.org/content/104/50/19926.abstract for published article

Journal ref: PNAS 2007 104:19926-19930

Life at the front of an expanding population

Authors: Oskar Hallatschek, David R. Nelson

Abstract: Recent microbial experiments suggest that enhanced genetic drift at the frontier of a two-dimensional range expansion can cause genetic sectoring patterns with fractal domain boundaries. Here, we propose and analyze a simple model of asexual biological evolution at expanding frontiers to explain these neutral patterns and predict the effect of natural selection. Our model attributes the observed… ▽ More Recent microbial experiments suggest that enhanced genetic drift at the frontier of a two-dimensional range expansion can cause genetic sectoring patterns with fractal domain boundaries. Here, we propose and analyze a simple model of asexual biological evolution at expanding frontiers to explain these neutral patterns and predict the effect of natural selection. Our model attributes the observed gradual decrease in the number of sectors at the leading edge to an unbiased random walk of sector boundaries. Natural selection introduces a deterministic bias in the wandering of domain boundaries that renders beneficial mutations more likely to escape genetic drift and become established in a sector. We find that the opening angle of those sectors and the rate at which they become established depend sensitively on the selective advantage of the mutants. Deleterious mutations, on the other hand, are not able to establish a sector permanently. They can, however, temporarily "surf" on the population front, and thereby reach unusual high frequencies. As a consequence, expanding frontiers are susceptible to deleterious mutations as revealed by the high fraction of mutants at mutation-selection balance. Numerically, we also determine the condition at which the wild type is lost in favor of deleterious mutants (genetic meltdown) at a growing front. Our prediction for this error threshold differs qualitatively from existing well-mixed theories, and sets tight constraints on sustainable mutation rates for populations that undergo frequent range expansions. △ Less

Submitted 12 December, 2008; v1 submitted 1 October, 2008; originally announced October 2008.

Comments: Update

Genome landscapes and bacteriophage codon usage

Authors: Julius B. Lucks, David R. Nelson, Grzegorz Kudla, Joshua B. Plotkin

Abstract: Across all kingdoms of biological life, protein-coding genes exhibit unequal usage of synonmous codons. Although alternative theories abound, translational selection has been accepted as an important mechanism that shapes the patterns of codon usage in prokaryotes and simple eukaryotes. Here we analyze patterns of codon usage across 74 diverse bacteriophages that infect E. coli, P. aeruginosa an… ▽ More Across all kingdoms of biological life, protein-coding genes exhibit unequal usage of synonmous codons. Although alternative theories abound, translational selection has been accepted as an important mechanism that shapes the patterns of codon usage in prokaryotes and simple eukaryotes. Here we analyze patterns of codon usage across 74 diverse bacteriophages that infect E. coli, P. aeruginosa and L. lactis as their primary host. We introduce the concept of a `genome landscape,' which helps reveal non-trivial, long-range patterns in codon usage across a genome. We develop a series of randomization tests that allow us to interrogate the significance of one aspect of codon usage, such a GC content, while controlling for another aspect, such as adaptation to host-preferred codons. We find that 33 phage genomes exhibit highly non-random patterns in their GC3-content, use of host-preferred codons, or both. We show that the head and tail proteins of these phages exhibit significant bias towards host-preferred codons, relative to the non-structural phage proteins. Our results support the hypothesis of translational selection on viral genes for host-preferred codons, over a broad range of bacteriophages. △ Less

Submitted 14 August, 2007; originally announced August 2007.

Comments: 9 Color Figures, 5 Tables, 53 References

Journal ref: Lucks JB, Nelson DR, Kudla GR, Plotkin JB (2008) Genome Landscapes and Bacteriophage Codon Usage. PLoS Computational Biology 4(2): e1000001

Gene surfing

Authors: Oskar Hallatschek, David R. Nelson

Abstract: Spatially resolved genetic data is increasingly used to reconstruct the migrational history of species. To assist such inference, we study, by means of simulations and analytical methods, the dynamics of neutral gene frequencies in a population undergoing a continual range expansion in one dimension. During such a colonization period, lineages can fix at the wave front by means of a ``surfing''… ▽ More Spatially resolved genetic data is increasingly used to reconstruct the migrational history of species. To assist such inference, we study, by means of simulations and analytical methods, the dynamics of neutral gene frequencies in a population undergoing a continual range expansion in one dimension. During such a colonization period, lineages can fix at the wave front by means of a ``surfing'' mechanism [Edmonds C.A., Lillie A.S. & Cavalli-Sforza L.L. (2004) Proc Natl Acad Sci USA 101: 975-979]. We quantify this phenomenon in terms of (i) the spatial distribution of lineages that reach fixation and, closely related, (ii) the continual loss of genetic diversity (heterozygosity) at the wave front, characterizing the approach to fixation. Our simulations show that an effective population size can be assigned to the wave that controls the (observable) gradient in heterozygosity left behind the colonization process. This effective population size is markedly higher in pushed waves than in pulled waves, and increases only sub-linearly with deme size. To explain these and other findings, we develop a versatile analytical approach, based on the physics of reaction-diffusion systems, that yields simple predictions for any deterministic population dynamics. △ Less

Submitted 17 March, 2007; originally announced March 2007.

Journal ref: Theoretical Population Biology, 73 (1), p. 158, 2008.

Sequence heterogeneity and the dynamics of molecular motors

Authors: Yariv Kafri, David R. Nelson

Abstract: The effect of sequence heterogeneity on the dynamics of molecular motors is reviewed and analyzed using a set of recently introduced lattice models. First, we review results for the influence of heterogenous tracks such as a single-strand of DNA or RNA on the dynamics of the motors. We stress how the predicted behavior might be observed experimentally in anomalous drift and diffusion of motors o… ▽ More The effect of sequence heterogeneity on the dynamics of molecular motors is reviewed and analyzed using a set of recently introduced lattice models. First, we review results for the influence of heterogenous tracks such as a single-strand of DNA or RNA on the dynamics of the motors. We stress how the predicted behavior might be observed experimentally in anomalous drift and diffusion of motors over a wide range of parameters near the stall force and discuss the extreme limit of strongly biased motors with one-way hopping. We then consider the dynamics in an environment containing a variety of different fuels which supply chemical energy for the motor motion, either on a heterogeneous or on a periodic track. The results for motion along a periodic track are relevant to kinesin motors in a solution with a mixture of different nucleotide triphosphate fuel sources. △ Less

Submitted 30 June, 2005; originally announced June 2005.

Comments: To appear in a JPhys special issue on molecular motors

Dynamics of Molecular Motors with Finite Processivity on Heterogeneous Tracks

Authors: Yariv Kafri, David K. Lubensky, David R. Nelson

Abstract: The dynamics of molecular motors which occasionally detach from a heterogeneous track like DNA or RNA is considered. Motivated by recent single molecule experiments, we study a simple model for a motor moving along a disordered track using chemical energy while an external force opposes its motion. The motors also have finite processivity, i.e., they can leave the track with a position dependent… ▽ More The dynamics of molecular motors which occasionally detach from a heterogeneous track like DNA or RNA is considered. Motivated by recent single molecule experiments, we study a simple model for a motor moving along a disordered track using chemical energy while an external force opposes its motion. The motors also have finite processivity, i.e., they can leave the track with a position dependent rate. We show that the response of the system to disorder in the hopping off rate depends on the value of the external force. For most values of the external force, strong disorder causes the motors which survive for long times on the track to be localized at preferred positions. However, near the stall force, localization occurs for any amount of disorder. Existence of localized states near the top of the band implies a stretched exponential contribution to the decay of the survival probability. To obtain these results, we study the complex eigenvalue spectrum of the time evolution operator. A similar spectral analysis also provides a very efficient method for studying the dynamics of motors with infinite processivity. △ Less

Submitted 8 September, 2004; originally announced September 2004.

Comments: 19 Pages, 21 Figures

Pause Point Spectra in DNA Constant-Force Unzipping

Authors: J. D. Weeks, J. B. Lucks, Y. Kafri, C. Danilowicz, D. R. Nelson, M. Prentiss

Abstract: Under constant applied force, the separation of double-stranded DNA into two single strands is known to proceed through a series of pauses and jumps. Given experimental traces of constant-force unzipping, we present a method whereby the locations of pause points can be extracted in the form of a pause point spectrum. A simple theoretical model of DNA constant-force unzipping is demonstrated to p… ▽ More Under constant applied force, the separation of double-stranded DNA into two single strands is known to proceed through a series of pauses and jumps. Given experimental traces of constant-force unzipping, we present a method whereby the locations of pause points can be extracted in the form of a pause point spectrum. A simple theoretical model of DNA constant-force unzipping is demonstrated to produce good agreement with the experimental pause point spectrum of lambda phage DNA. The locations of peaks in the experimental and theoretical pause point spectra are found to be nearly coincident below 6000 bp. The model only requires the sequence, temperature and a set of empirical base pair binding and stacking energy parameters, and the good agreement with experiment suggests that pause points are primarily determined by the DNA sequence. The model is also used to predict pause point spectra for the BacterioPhage PhiX174 genome. The algorithm for extracting the pause point spectrum might also be useful for studying related systems which exhibit pausing behavior such as molecular motors. △ Less

Submitted 10 June, 2004; originally announced June 2004.

Comments: 15 pages, 12 figures

Journal ref: Biophysical Journal, 88, 2752-2765, (2005).

Dynamics of Molecular Motors and Polymer Translocation with Sequence Heterogeneity

Authors: Yariv Kafri, David K. Lubensky, David R. Nelson

Abstract: The effect of sequence heterogeneity on polynucleotide translocation across a pore and on simple models of molecular motors such as helicases, DNA polymerase/exonuclease and RNA polymerase is studied in detail. Pore translocation of RNA or DNA is biased due to the different chemical environments on the two sides of the membrane, while the molecular motor motion is biased through a coupling to ch… ▽ More The effect of sequence heterogeneity on polynucleotide translocation across a pore and on simple models of molecular motors such as helicases, DNA polymerase/exonuclease and RNA polymerase is studied in detail. Pore translocation of RNA or DNA is biased due to the different chemical environments on the two sides of the membrane, while the molecular motor motion is biased through a coupling to chemical energy. An externally applied force can oppose these biases. For both systems we solve lattice models exactly both with and without disorder. The models incorporate explicitly the coupling to the different chemical environments for polymer translocation and the coupling to the chemical energy (as well as nucleotide pairing energies) for molecular motors. Using the exact solutions and general arguments we show that the heterogeneity leads to anomalous dynamics. Most notably, over a range of forces around the stall force (or stall tension for DNA polymerase/exonuclease systems) the displacement grows sublinearly as t^μwith μ<1. The range over which this behavior can be observed experimentally is estimated for several systems and argued to be detectable for appropriate forces and buffers. Similar sequence heterogeneity effects may arise in the packing of viral DNA. △ Less

Submitted 9 June, 2004; v1 submitted 20 October, 2003; originally announced October 2003.

Comments: 42 pages, 12 figures

Journal ref: Biophys. J. 86, 3373 (2004)

A Quasispecies on a Moving Oasis

Authors: Michael M. Desai, David R. Nelson

Abstract: A population evolving in an inhomogeneous environment will adapt differently to different regions. We study the conditions under which such a population can maintain adaptations to a particular region when that region is not stationary, but can move. In particular, we study a quasispecies living near a favorable patch ("oasis") in the middle of a large "desert." The population has two genetic st… ▽ More A population evolving in an inhomogeneous environment will adapt differently to different regions. We study the conditions under which such a population can maintain adaptations to a particular region when that region is not stationary, but can move. In particular, we study a quasispecies living near a favorable patch ("oasis") in the middle of a large "desert." The population has two genetic states, one of which which conveys a relative advantage while in the oasis at the cost of a disadvantage in the desert. We consider the population dynamics when the oasis is moving, or equivalently some form of "wind" is blowing the population away from the oasis. We find that the ratio of the two types of individuals exhibits sharp transitions at particular oasis velocities. We calculate an extinction velocity, and a switching velocity above which the dominance switches from the oasis-adapted genotype to the desert-adapted one. This switching velocity is analagous to the quasispecies mutational error threshold. Above this velocity, the population cannot maintain adaptations to the properties of the oasis. △ Less

Submitted 26 September, 2003; originally announced September 2003.

Comments: 12 pages, 5 figures

Statistical Physics of Unzipping DNA

Authors: David R. Nelson

Abstract: The denaturation of double-stranded DNA as function of force and temperature is discussed. At room temperature, sequence heterogeneity dominates the physics of single molecule force-extension curves starting about 7 piconewtons of below a ~15 pN unzipping transition. The dynamics of the unzipping fork exhibits anomalous drift and diffusion in a similar range above this transition. Energy barrier… ▽ More The denaturation of double-stranded DNA as function of force and temperature is discussed. At room temperature, sequence heterogeneity dominates the physics of single molecule force-extension curves starting about 7 piconewtons of below a ~15 pN unzipping transition. The dynamics of the unzipping fork exhibits anomalous drift and diffusion in a similar range above this transition. Energy barriers near the transition scale as the square root of the genome size. Recent observations of jumps and plateaus in the unzipping of lambda phage DNA at constant force are consistent with these predictions. △ Less

Submitted 24 September, 2003; originally announced September 2003.

Virus shapes and buckling transitions in spherical shells

Authors: Jack Lidmar, Leonid Mirny, David R. Nelson

Abstract: We show that the icosahedral packings of protein capsomeres proposed by Caspar and Klug for spherical viruses become unstable to faceting for sufficiently large virus size, in analogy with the buckling instability of disclinations in two-dimensional crystals. Our model, based on the nonlinear physics of thin elastic shells, produces excellent one parameter fits in real space to the full three-di… ▽ More We show that the icosahedral packings of protein capsomeres proposed by Caspar and Klug for spherical viruses become unstable to faceting for sufficiently large virus size, in analogy with the buckling instability of disclinations in two-dimensional crystals. Our model, based on the nonlinear physics of thin elastic shells, produces excellent one parameter fits in real space to the full three-dimensional shape of large spherical viruses. The faceted shape depends only on the dimensionless Foppl-von Karman number γ=YR^2/κ, where Y is the two-dimensional Young's modulus of the protein shell, κis its bending rigidity and R is the mean virus radius. The shape can be parameterized more quantitatively in terms of a spherical harmonic expansion. We also investigate elastic shell theory for extremely large γ, 10^3 < γ< 10^8, and find results applicable to icosahedral shapes of large vesicles studied with freeze fracture and electron microscopy. △ Less

Submitted 27 November, 2003; v1 submitted 30 June, 2003; originally announced June 2003.

Comments: 11 pages, 12 figures

Journal ref: Phys. Rev. E 68, 051910 (2003)

DNA unzipped under a constant force exhibits multiple metastable intermediates

Authors: Claudia Danilowicz, Vincent W. Coljee, Cedric Bouzigues, David K. Lubensky, David R. Nelson, Mara Prentiss

Abstract: Single molecule studies, at constant force, of the separation of double-stranded DNA into two separated single strands may provide information relevant to the dynamics of DNA replication. At constant applied force, theory predicts that the unzipped length as a function of time is characterized by jumps during which the strands separate rapidly, followed by long pauses where the number of separat… ▽ More Single molecule studies, at constant force, of the separation of double-stranded DNA into two separated single strands may provide information relevant to the dynamics of DNA replication. At constant applied force, theory predicts that the unzipped length as a function of time is characterized by jumps during which the strands separate rapidly, followed by long pauses where the number of separated base pairs remains constant. Here, we report previously uncharacterized observations of this striking behavior carried out on a number of identical single molecules simultaneously. When several single lphage molecules are subject to the same applied force, the pause positions are reproducible in each. This reproducibility shows that the positions and durations of the pauses in unzipping provide a sequence-dependent molecular fingerprint. For small forces, the DNA remains in a partially unzipped state for at least several hours. For larger forces, the separation is still characterized by jumps and pauses, but the double-stranded DNA will completely unzip in less than 30 min. △ Less

Submitted 7 February, 2003; v1 submitted 7 February, 2003; originally announced February 2003.

Journal ref: PNAS Volume 100 1694-1699 (2003)

Single Molecule Statistics and the Polynucleotide Unzipping Transition

Authors: David K. Lubensky, David R. Nelson

Abstract: We present an extensive theoretical investigation of the mechanical unzipping of double-stranded DNA under the influence of an applied force. In the limit of long polymers, there is a thermodynamic unzipping transition at a critical force value of order 10 pN, with different critical behavior for homopolymers and for random heteropolymers. We extend results on the disorder-averaged behavior of D… ▽ More We present an extensive theoretical investigation of the mechanical unzipping of double-stranded DNA under the influence of an applied force. In the limit of long polymers, there is a thermodynamic unzipping transition at a critical force value of order 10 pN, with different critical behavior for homopolymers and for random heteropolymers. We extend results on the disorder-averaged behavior of DNA's with random sequences to the more experimentally accessible problem of unzipping a single DNA molecule. As the applied force approaches the critical value, the double-stranded DNA unravels in a series of discrete, sequence-dependent steps that allow it to reach successively deeper energy minima. Plots of extension versus force thus take the striking form of a series of plateaus separated by sharp jumps. Similar qualitative features should reappear in micromanipulation experiments on proteins and on folded RNA molecules. Despite their unusual form, the extension versus force curves for single molecules still reveal remnants of the disorder-averaged critical behavior. Above the transition, the dynamics of the unzipping fork is related to that of a particle diffusing in a random force field; anomalous, disorder-dominated behavior is expected until the applied force exceeds the critical value for unzipping by roughly 5 pN. △ Less

Submitted 19 July, 2001; originally announced July 2001.

Comments: 40 pages, 18 figures

Pulling Pinned Polymers and Unzipping DNA

Authors: David K. Lubensky, David R. Nelson

Abstract: We study a class of micromanipulation experiments, exemplified by the pulling apart of the two strands of double-stranded DNA (dsDNA). When the pulling force is increased to a critical value, an ``unzipping'' transition occurs. For random DNA sequences with short-ranged correlations, we obtain exact results for the number of monomers liberated and the specific heat, including the critical behavi… ▽ More We study a class of micromanipulation experiments, exemplified by the pulling apart of the two strands of double-stranded DNA (dsDNA). When the pulling force is increased to a critical value, an ``unzipping'' transition occurs. For random DNA sequences with short-ranged correlations, we obtain exact results for the number of monomers liberated and the specific heat, including the critical behavior at the transition. Related systems include a random heteropolymer pulled away from an adsorbing surface and a vortex line in a type II superconductor tilted away from a fragmented columnar defect. △ Less

Submitted 2 August, 2000; v1 submitted 26 April, 2000; originally announced April 2000.

Comments: 4 pages, 1 EPS figure; revised references and very brief discussion of order of the transition added; to appear in Phys. Rev. Lett

Journal ref: Phys. Rev. Lett. _85_, 1572 (2000)

Population Dynamics and Non-Hermitian Localization

Authors: Karin A. Dahmen, David R. Nelson, Nadav M. Shnerb

Abstract: We review localization with non-Hermitian time evolution as applied to simple models of population biology with spatially varying growth profiles and convection. Convection leads to a constant imaginary vector potential in the Schroedinger-like operator which appears in linearized growth models. We illustrate the basic ideas by reviewing how convection affects the evolution of a population influ… ▽ More We review localization with non-Hermitian time evolution as applied to simple models of population biology with spatially varying growth profiles and convection. Convection leads to a constant imaginary vector potential in the Schroedinger-like operator which appears in linearized growth models. We illustrate the basic ideas by reviewing how convection affects the evolution of a population influenced by a simple square well growth profile. Results from discrete lattice growth models in both one and two dimensions are presented. A set of similarity transformations which lead to exact results for the spectrum and winding numbers of eigenfunctions for random growth rates in one dimension is described in detail. We discuss the influence of boundary conditions, and argue that periodic boundary conditions lead to results which are in fact typical of a broad class of growth problems with convection. △ Less

Submitted 18 March, 1999; originally announced March 1999.

Comments: 19 pages, 11 figures

Driven Polymer Translocation Through a Narrow Pore

Authors: David K. Lubensky, David R. Nelson

Abstract: Motivated by experiments in which a polynucleotide is driven through a proteinaceous pore by an electric field, we study the diffusive motion of a polymer threaded through a narrow channel with which it may have strong interactions. We show that there is a range of polymer lengths in which the system is approximately translationally invariant, and we develop a coarse-grained description of this… ▽ More Motivated by experiments in which a polynucleotide is driven through a proteinaceous pore by an electric field, we study the diffusive motion of a polymer threaded through a narrow channel with which it may have strong interactions. We show that there is a range of polymer lengths in which the system is approximately translationally invariant, and we develop a coarse-grained description of this regime. From this description, general features of the distribution of times for the polymer to pass through the pore may be deduced. We also introduce a more microscopic model. This model provides a physically reasonable scenario in which, as in experiments, the polymer's speed depends sensitively on its chemical composition, and even on its orientation in the channel. Finally, we point out that the experimental distribution of times for the polymer to pass through the pore is much broader than expected from simple estimates, and speculate on why this might be. △ Less

Submitted 16 February, 1999; originally announced February 1999.

Comments: 16 pages, 8 figures, RevTex and harvard citation style, submitted to Biophysical Journal

Life and Death near a Windy Oasis

Authors: Karin A. Dahmen, David R. Nelson, Nadav M. Shnerb

Abstract: We propose a simple experiment to study delocalization and extinction in inhomogeneous biological systems. The nonlinear steady state for, say, a bacteria colony living on and near a patch of nutrient or favorable illumination (``oasis'') in the presence of a drift term (``wind'') is computed. The bacteria, described by a simple generalization of the Fisher equation, diffuse, divide A -> A+A, di… ▽ More We propose a simple experiment to study delocalization and extinction in inhomogeneous biological systems. The nonlinear steady state for, say, a bacteria colony living on and near a patch of nutrient or favorable illumination (``oasis'') in the presence of a drift term (``wind'') is computed. The bacteria, described by a simple generalization of the Fisher equation, diffuse, divide A -> A+A, die A -> 0, and annihilate A+A -> 0. At high wind velocities all bacteria are blown into an unfavorable region (``desert''), and the colony dies out. At low velocity a steady state concentration survives near the oasis. In between these two regimes there is a critical velocity at which bacteria first survive. If the ``desert'' supports a small nonzero population, this extinction transition is replaced by a delocalization transition with increasing velocity. Predictions for the behavior as a function of wind velocity are made for one and two dimensions. △ Less

Submitted 30 July, 1998; originally announced July 1998.

Comments: 18 pages, 9 figures

Non-Hermitian Localization and Population Biology

Authors: David R. Nelson, Nadav M. Shnerb

Abstract: The time evolution of spatial fluctuations in inhomogeneous d-dimensional biological systems is analyzed. A single species continuous growth model, in which the population disperses via diffusion and convection is considered. Time-independent environmental heterogeneities, such as a random distribution of nutrients or sunlight are modeled by quenched disorder in the growth rate. Linearization of… ▽ More The time evolution of spatial fluctuations in inhomogeneous d-dimensional biological systems is analyzed. A single species continuous growth model, in which the population disperses via diffusion and convection is considered. Time-independent environmental heterogeneities, such as a random distribution of nutrients or sunlight are modeled by quenched disorder in the growth rate. Linearization of this model of population dynamics shows that the fastest growing localized state dominates in a time proportional to a power of the logarithm of the system size. Using an analogy with a Schrodinger equation subject to a constant imaginary vector potential, we propose a delocalization transition for the steady state of the nonlinear problem at a critical convection threshold separating localized and extended states. In the limit of high convection velocity, the linearized growth problem in $d$ dimensions exhibits singular scaling behavior described by a (d-1)-dimensional generalization of the noisy Burgers' equation, with universal singularities in the density of states associated with disorder averaged eigenvalues near the band edge in the complex plane. The Burgers mapping leads to unusual transverse spreading of convecting delocalized populations. △ Less

Submitted 26 January, 1998; v1 submitted 10 August, 1997; originally announced August 1997.

Comments: 22 pages, 11 figures