Wolfram Möbius

Quantitative test of the barrier nucleosome model for statistical positioning of nucleosomes up- and downstream of transcription start sites.

The positions of nucleosomes in eukaryotic genomes determine which parts of the DNA sequence are readily accessible for regulatory proteins and which are not. Genome-wide maps of nucleosome positions have revealed a salient pattern around transcription start sites, involving a nucleosome-free region (NFR) flanked by a pronounced periodic pattern in the average nucleosome density. While the periodic pattern clearly reflects well-positioned nucleosomes, the positioning mechanism is less clear. A recent experimental study by Mavrich et al. argued that the pattern observed in Saccharomyces cerevisiae is qualitatively consistent with a “barrier nucleosome model,” in which the oscillatory pattern is created by the statistical positioning mechanism of Kornberg and Stryer. On the other hand, there is clear evidence for intrinsic sequence preferences of nucleosomes, and it is unclear to what extent these sequence preferences affect the observed pattern. To test the barrier nucleosome model, we quantitatively analyze yeast nucleosome positioning data both up- and downstream from NFRs. Our analysis is based on the Tonks model of statistical physics which quantifies the interplay between the excluded-volume interaction of nucleosomes and their positional entropy. We find that although the typical patterns on the two sides of the NFR are different, they are both quantitatively described by the same physical model with the same parameters, but different boundary conditions. The inferred boundary conditions suggest that the first nucleosome downstream from the NFR (the +1 nucleosome) is typically directly positioned while the first nucleosome upstream is statistically positioned via a nucleosome-repelling DNA region. These boundary conditions, which can be locally encoded into the genome sequence, significantly shape the statistical distribution of nucleosomes over a range of up to ∼1,000 bp to each side.

Toward a unified physical model of nucleosome patterns flanking transcription start sites

Recent genome-wide maps of nucleosome positions in different eukaryotes revealed patterns around transcription start sites featuring a nucleosome-free region flanked by a periodic modulation of the nucleosome density. For Saccharomyces cerevisiae, the average in vivo pattern was previously shown to be quantitatively described by a “nucleosome gas” model based on the statistical positioning mechanism. However, this simple physical description is challenged by the fact that the pattern differs quantitatively between species and by recent experiments that appear incompatible with statistical positioning, indicating important roles for chromatin remodelers. We undertake a data-driven search for a unified physical model to describe the nucleosome patterns of 12 yeast species and also consider an extension of the model to capture remodeling effects. We are led to a nucleosome gas that takes into account nucleosome breathing, i.e., transient unwrapping of nucleosomal DNA segments. This known biophysical property of nucleosomes rationalizes a “pressure”-induced dependence of the effective nucleosome size that is suggested by the data. By fitting this model to the data, we find an average energy cost for DNA unwrapping consistent with previous biophysical experiments. Although the available data are not sufficient to reconstruct chromatin remodeling mechanisms, a minimal model extension by one mechanism yields an “active nucleosome gas” that can rationalize the behavior of systems with reduced histone–DNA ratio and remodeler knockouts. We therefore establish a basis for a physical description of nucleosome patterns that can serve as a null model for sequence-specific effects at individual genes and in models of transcription regulation.

Physical and Mathematical Modeling in Experimental Papers.

An increasing number of publications include modeling. Often, such studies help us to gain a deeper insight into the phenomena studied and break down barriers between experimental and theoretical communities. However, combining experimental and theoretical work is challenging for authors, reviewers, and readers. To help maximize the usefulness and impact of combined theoretical and experimental research, this Primer describes the purpose, usefulness, and different types of models and addresses the practical aspect of integrated publications by outlining characteristics of good modeling, presentation, and fruitful collaborations.

How Obstacles Perturb Population Fronts and Alter Their Genetic Structure.

As populations spread into new territory, environmental heterogeneities can shape the population front and genetic composition. We focus here on the effects of an important building block of heterogeneous environments, isolated obstacles. With a combination of experiments, theory, and simulation, we show how isolated obstacles both create long-lived distortions of the front shape and amplify the effect of genetic drift. A system of bacteriophage T7 spreading on a spatially heterogeneous Escherichia coli lawn serves as an experimental model system to study population expansions. Using an inkjet printer, we create well-defined replicates of the lawn and quantitatively study the population expansion of phage T7. The transient perturbations of the population front found in the experiments are well described by a model in which the front moves with constant speed. Independent of the precise details of the expansion, we show that obstacles create a kink in the front that persists over large distances and is insensitive to the details of the obstacle’s shape. The small deviations between experimental findings and the predictions of the constant speed model can be understood with a more general reaction-diffusion model, which reduces to the constant speed model when the obstacle size is large compared to the front width. Using this framework, we demonstrate that frontier genotypes just grazing the side of an isolated obstacle increase in abundance, a phenomenon we call ‘geometry-enhanced genetic drift’, complementary to the founder effect associated with spatial bottlenecks. Bacterial range expansions around nutrient-poor barriers and stochastic simulations confirm this prediction. The effect of the obstacle on the genealogy of individuals at the front is characterized by simulations and rationalized using the constant speed model. Lastly, we consider the effect of two obstacles on front shape and genetic composition of the population illuminating the effects expected from complex environments with many obstacles.

The spatial and metabolic basis of colony size variation

Spatial structure impacts microbial growth and interactions, with ecological and evolutionary consequences. It is therefore important to quantitatively understand how spatial proximity affects interactions in different environments. We tested how proximity influences colony size when either Escherichia coli or Salmonella enterica are grown on various carbon sources. The importance of colony location changed with species and carbon source. Spatially explicit, genome-scale metabolic modeling recapitulated observed colony size variation. Competitors that determine territory size, according to Voronoi diagrams, were the most important drivers of variation in colony size. However, the relative importance of different competitors changed through time. Further, the effect of location increased when colonies took up resources quickly relative to the diffusion of limiting resources. These analyses made it apparent that the importance of location was smaller than expected for experiments with S. enterica growing on glucose. The accumulation of toxic byproducts appeared to limit the growth of large colonies and reduced variation in colony size. Our work provides an experimentally and theoretically grounded understanding of how location interacts with metabolism and diffusion to influence microbial interactions.

Genetic drift and selection in many-allele range expansions

We experimentally and numerically investigate the evolutionary dynamics of four competing strains of E. coli with differing expansion velocities in radially expanding colonies. We compare experimental measurements of the average fraction, correlation functions between strains, and the relative rates of genetic domain wall annihilations and coalescences to simulations modeling the population as a one-dimensional ring of annihilating and coalescing random walkers with deterministic biases due to selection. The simulations reveal that the evolutionary dynamics can be collapsed onto master curves governed by three essential parameters: (1) an expansion length beyond which selection dominates over genetic drift; (2) a characteristic angular correlation describing the size of genetic domains; and (3) a dimensionless constant quantifying the interplay between a colony’s curvature at the frontier and its selection length scale. We measure these parameters with a new technique that precisely measures small selective differences between spatially competing strains and show that our simulations accurately predict the dynamics without additional fitting. Our results suggest that the random walk model can act as a useful predictive tool for describing the evolutionary dynamics of range expansions composed of an arbitrary number of genotypes with different fitnesses.

Emergence of robust nucleosome patterns from an interplay of positioning mechanisms

Proper positioning of nucleosomes in eukaryotic cells is determined by a complex interplay of factors, including nucleosome-nucleosome interactions, DNA sequence, and active chromatin remodeling. Yet, characteristic features of nucleosome positioning, such as gene-averaged nucleosome patterns, are surprisingly robust across perturbations, conditions, and species. Here, we explore how this robustness arises despite the underlying complexity. We leverage mathematical models to show that a large class of positioning mechanisms merely affects the quantitative characteristics of qualitatively robust positioning patterns. We demonstrate how statistical positioning emerges as an effective description from the complex interplay of different positioning mechanisms, which ultimately only renormalize the model parameter quantifying the effective softness of nucleosomes. This renormalization can be species-specific, rationalizing a puzzling discrepancy between the effective nucleosome softness of S. pombe and S. cerevisiae. More generally, we establish a quantitative framework for dissecting the interplay of different nucleosome positioning determinants.

Evolution of populations expanding on curved surfaces(a)

The expansion of a population into new habitat is a transient process that leaves its footprints in the genetic composition of the expanding population. How the structure of the environment shapes the population front and the evolutionary dynamics during such a range expansion is little understood. Here, we investigate the evolutionary dynamics of populations consisting of many selectively neutral genotypes expanding on curved surfaces. Using a combination of individual-based off-lattice simulations, geometrical arguments, and lattice-based stepping-stone simulations, we characterise the effect of individual bumps on an otherwise flat surface. Compared to the case of a range expansion on a flat surface, we observe a transient relative increase, followed by a decrease, in neutral genetic diversity at the population front. In addition, we find that individuals at the sides of the bump have a dramatically increased expected number of descendants, while their neighbours closer to the bumps centre are far less lucky. Both observations can be explained using an analytical description of straight paths (geodesics) on the curved surface. Complementing previous studies of heterogeneous flat environments, the findings here build our understanding of how complex environments shape the evolutionary dynamics of expanding populations.

TurBase: A Software Platform for Research in Experimental and Numerical Fluid Dynamics

We present a software infrastructure for the research community working on turbulence and complex flows (TurBase), an easily accessible web platform for high quality data. Its main goal is to host, standardize and manage a large collections of heterogeneous experimental and numerical data sets from high-end European fluid dynamics experimental facilities and from High Performance Computational centres. TurBase offers scalable performances when accessing/uploading/searching data, providing at the same time maximum flexibility and power (through Jupyter notebooks) when doing online computation directly on big datasets