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.

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.

We study the localization of vibrational modes of frictionless granular media. We introduce a new method, motivated by earlier work on non-Hermitian quantum problems, which works well both in the localized regime where the localization length (\xi) is much less than the linear size (L) and in the regime (\xi \gsim L) when modes are extended throughout our finite system. Our very lowest frequency modes show "quasi-localized" resonances away from the jamming point; the spatial extent of these regions increases as the jamming point is approached, as expected theoretically. Throughout the remaining frequency range, our data show no signature of the nearness of the jamming point and collapse well when properly rescaled with the system size. Using Random Matrix Theory, we derive the scaling relation (\xi \sim L^{d/2}) for the regime (\xi \gg L) in (d) dimensions.

Across all kingdoms of biological life, protein-coding genes exhibit unequal usage of synonymous 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 use 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 as 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.

Hydrostatically pressurized circular rings confined to two dimensions (or cylinders constrained to have only z-independent deformations) undergo Euler-type buckling when the outside pressure exceeds a critical value. We perform a stability analysis of rings with arclength-dependent bending moduli and determine how weakened bending modulus segments affect the buckling critical pressure. Rings with a fourfold symmetric modulation are particularly susceptible to collapse. In addition we study the initial postbuckling stages of the pressurized rings to determine possible ring folding patterns.

We consider two-dimensional dispersions of droplets of isotropic phase in a liquid with an XY-like order parameter, tilt, nematic, and hexatic symmetries being included. Strong anchoring boundary conditions are assumed. Textures for a single droplet and a pair of droplets are calculated and a universal droplet-droplet pair potential is obtained. The interaction of dispersed droplets via the ordered phase is attractive at large distances and repulsive at short distances, which results in a well defined preferred separation for two droplets and topological stabilization of the emulsion. This interaction also drives self-assembly into chains. Preferred separations and energy barriers to coalescence are calculated, and the effects of thermal fluctuations and film thickness are discussed.