Person:

Nelson, David

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.

We study theoretically the mechanical failure of a simple model of double stranded DNA under an applied shear. Starting from a more microscopic Hamiltonian that describes a sheared DNA, we arrive at a nonlinear generalization of a ladder model of shear unzipping proposed earlier by deGennes [deGennes P. G. C. R. Acad. Sci., Ser. IV; Phys., Astrophys. 2001, 1505]. Using this model and a combination of analytical and numerical methods, we study the DNA "unzipping" transition when the shearing force exceeds a critical threshold at zero temperature. We also explore the effects of sequence heterogeneity and finite temperature and discuss possible applications to determine the strength of colloidal nanoparticle assemblies functionalized by DNA.

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.

As an example of life at high Reynolds number, 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 remarkably long persistence times), sinks in the turbulent field. An important consequence is a large reduction in the carrying capacity of the fluid medium. We analytically determine the regimes where this quasi-localized behavior occurs and test our predictions by numerical simulations.

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 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 (DP) process, with a phase diagram strongly influenced by an exceptional symmetric DP (DP2) transition.

Although the effects of kinetics on crystal growth are well understood, the role of substrate curvature is not yet established. We studied rigid, two-dimensional colloidal crystals growing on spherical droplets to understand how the elastic stress induced by Gaussian curvature affects the growth pathway. In contrast to crystals grown on flat surfaces or compliant crystals on droplets, these crystals formed branched, ribbon-like domains with large voids and no topological defects. We show that this morphology minimizes the curvature-induced elastic energy. Our results illustrate the effects of curvature on the ubiquitous process of crystallization, with practical implications for nanoscale disorder-order transitions on curved manifolds, including the assembly of viral capsids, phase separation on vesicles, and crystallization of tetrahedra in three dimensions.

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 ℓ 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 among different cell types, ranging from bacteria and yeast to human tissue. The thickness ℓ decreases with increasing ν, increasing cell volume fraction ϕ, and decreasing ambient nutrient concentration (ψ_∞). The theoretical results are compared with numerical simulations and experiments. In the latter studies, colonies of budding yeast, 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 ℓ.

We study the competition between domain coarsening in a symmetric binary mixture below critical temperature and turbulent fluctuations. We find that the coarsening process is arrested in the presence of turbulence. The physics of the process shares remarkable similarities with the behavior of diluted turbulent emulsions and the arrest length scale can be estimated with an argument similar to the one proposed by Kolmogorov and Hinze for the maximal stability diameter of droplets in turbulence. Although, in the absence of flow, the microscopic diffusion constant is negative, turbulence does effectively arrest the inverse cascade of concentration fluctuations by making the low wavelength diffusion constant positive for scales above the Hinze length.

Colloidal capsules can sustain an external osmotic pressure; however, for a sufficiently large pressure, they will ultimately buckle. This process can be strongly influenced by structural inhomogeneities in the capsule shells. We explore how the time delay before the onset of buckling decreases as the shells are made more inhomogeneous; this behavior can be quantitatively understood by coupling shell theory with Darcy’s law. In addition, we show that the shell inhomogeneity can dramatically change the folding pathway taken by a capsule after it buckles.

On microscopic scales, the crystallinity of flexible tethered or cross-linked membranes determines their mechanical response. We show that by controlling the type, number, and distribution of defects on a spherical elastic shell, it is possible to direct the morphology of these structures. Our numerical simulations show that by deflating a crystalline shell with defects, we can create elastic shell analogs of the classical platonic solids. These morphologies arise via a sharp buckling transition from the sphere which is strongly hysteretic in loading or unloading. We construct a minimal Landau theory for the transition using quadratic and cubic invariants of the spherical harmonic modes. Our approach suggests methods to engineer shape into soft spherical shells using a frozen defect topology.