Most cancers in humans are large, measuring centimetres in diameter, and composed of many billions of cells1. An equivalent mass of normal cells would be highly heterogeneous as a result of the mutations that occur during each cell division. What is remarkable about cancers is that virtually every neoplastic cell within a large tumour often contains the same core set of genetic alterations, with heterogeneity confined to mutations that emerge late during tumour growth2–5. How such alterations expand within the spatially constrained three-dimensional architecture of a tumour, and come to dominate a large, pre-existing lesion, has been unclear. Here we describe a model for tumour evolution that shows how short-range dispersal and cell turnover can account for rapid cell mixing inside the tumour. We show that even a small selective advantage of a single cell within a large tumour allows the descendants of that cell to replace the precursor mass in a clinically relevant time frame. We also demonstrate that the same mechanisms can be responsible for the rapid onset of resistance to chemotherapy. Our model not only provides insights into spatial and temporal aspects of tumour growth, but also suggests that targeting short-range cellular migratory activity could have marked effects on tumour growth rates.

We study colonies of nonmotile, rod-shaped bacteria growing on solid substrates. In our model, bacteria interact purely mechanically, by pushing each other away as they grow, and consume a diffusing nutrient. We show that mechanical interactions control the velocity and shape of the advancing front, which leads to features that cannot be captured by established Fisher-Kolmogorov models. In particular, we find that the velocity depends on the elastic modulus of bacteria or their stickiness to the surface. Interestingly, we predict that the radius of an incompressible, strictly two-dimensional colony cannot grow linearly in time, unless it develops branches. Importantly, mechanical interactions can also account for the nonequilibrium transition between circular and branching colonies, often observed in the lab.

We define a new class of random walk processes which maximize entropy. This maximal entropy random walk is equivalent to generic random walk if it takes place on a regular lattice, but it is not if the underlying lattice is irregular. In particular, we consider a lattice with weak dilution. We show that the stationary probability of finding a particle performing maximal entropy random walk localizes in the largest nearly spherical region of the lattice which is free of defects. This localization phenomenon, which is purely classical in nature, is explained in terms of the Lifshitz states of a certain random operator.

The coupling of ecology and evolution during range expansions enables mutations to establish at expanding range margins and reach high frequencies. This phenomenon, called allele surfing, is thought to have caused revolutions in the gene pool of many species, most evidently in microbial communities. It has remained unclear, however, under which conditions allele surfing promotes or hinders adaptation. Here, using microbial experiments and simulations, we show that, starting with standing adaptive variation, range expansions generate a larger increase in mean fitness than spatially uniform population expansions. The adaptation gain results from ‘soft’ selective sweeps emerging from surfing beneficial mutations. The rate of these surfing events is shown to sensitively depend on the strength of genetic drift, which varies among strains and environmental conditions. More generally, allele surfing promotes the rate of adaptation per biomass produced, which could help developing biofilms and other resource-limited populations to cope with environmental challenges.

We show that the eigenvalue density of a product X = X1X2 · · · XM of M independent N × N Gaussian random matrices in the limit N → ∞ is rotationally symmetric in the complex plane and is given by a simple expression ρ(z,z) =M for |z| ≤ σ, and is zero for |z| > σ. The parameter σ corresponds to the radius of the circular support and is related to the amplitude of the Gaussian fluctuations. This form of the eigenvalue density is highly universal. It is identical for products of Gaussian Hermitian, non-Hermitian, real or complex random matrices. It does not change even if the matrices in the product are taken from different Gaussian ensembles. We present a self-contained derivation of this result using a planar diagrammatic technique. Additionally, we conjecture that this distribution also holds for any matrices whose elements are independent, centered random variables with a finite variance or even more generally for matrices which fulfill Pastur-Lindeberg's condition. We provide a numerical evidence supporting this conjecture.

Drug gradients are believed to play an important role in the evolution of bacteria resistant to antibiotics and tumors resistant to anticancer drugs. We use a statistical physics model to study the evolution of a population of malignant cells exposed to drug gradients, where drug resistance emerges via a mutational pathway involving multiple mutations. We show that a nonuniform drug distribution has the potential to accelerate the emergence of resistance when the mutational pathway involves a long sequence of mutants with increasing resistance, but if the pathway is short or crosses a fitness valley, the evolution of resistance may actually be slowed down by drug gradients. These predictions can be verified experimentally, and may help to improve strategies for combating the emergence of resistance.

Mechanical forces are obviously important in the assembly of three-dimensional multicellular structures, but their detailed role is often unclear. We have used growing microcolonies of the bacterium Escherichia coli to investigate the role of mechanical forces in the transition from two-dimensional growth (on the interface between a hard surface and a soft agarose pad) to three-dimensional growth (invasion of the agarose). We measure the position within the colony where the invasion transition happens, the cell density within the colony and the colony size at the transition as functions of the concentration of the agarose. We use a phenomenological theory, combined with individual-based computer simulations, to show how mechanical forces acting between the bacterial cells, and between the bacteria and the surrounding matrix, lead to the complex phenomena observed in our experiments—in particular the observation that agarose concentration non-trivially affects the colony size at transition. Matching these approaches leads to a prediction for how the friction between the bacteria and the agarose should vary with agarose concentration. Our experimental conditions mimic numerous clinical and environmental scenarios in which bacteria invade soft matrices, as well as shedding more general light on the transition between two- and three-dimensional growth in multicellular assemblies.

Fitness effects of mutations depend on environmental parameters. For example, mutations that increase fitness of bacteria at high antibiotic concentration often decrease fitness in the absence of antibiotic, exemplifying a tradeoff between adaptation to environmental extremes. We develop a mathematical model for fitness landscapes generated by such tradeoffs, based on experiments that determine the antibiotic dose-response curves of Escherichia coli strains, and previous observations on antibiotic resistance mutations. Our model generates a succession of landscapes with predictable properties as antibiotic concentration is varied. The landscape is nearly smooth at low and high concentrations, but the tradeoff induces a high ruggedness at intermediate antibiotic concentrations. Despite this high ruggedness, however, all the fitness maxima in the landscapes are evolutionarily accessible from the wild type. This implies that selection for antibiotic resistance in multiple mutational steps is relatively facile despite the complexity of the underlying landscape.