Bacterial growth presents many beautiful phenomena that pose new theoretical challenges to statistical physicists, and are also amenable to laboratory experimentation. This review provides some of ...the essential biological background, discusses recent applications of statistical physics in this field, and highlights the potential for future research.
Bacterial conglomerates such as biofilms and microcolonies are ubiquitous in nature and play an important role in industry and medicine. In contrast to well-mixed cultures routinely used in microbial ...research, bacteria in a microcolony interact mechanically with one another and with the substrate to which they are attached. Here, we use a computer model of a microbial colony of rod-shaped cells to investigate how physical interactions between cells determine their motion in the colony and how this affects biological evolution. We show that the probability that a faster-growing mutant ‘surfs’ at the colony's frontier and creates a macroscopic sector depends on physical properties of cells (shape, elasticity and friction). Although all these factors contribute to the surfing probability in seemingly different ways, their effects can be summarized by two summary statistics that characterize the front roughness and cell alignment. Our predictions are confirmed by experiments in which we measure the surfing probability for colonies of different front roughness. Our results show that physical interactions between bacterial cells play an important role in biological evolution of new traits, and suggest that these interactions may be relevant to processes such as de novo evolution of antibiotic resistance.
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.
Recently available cancer sequencing data have revealed a complex view of the cancer genome containing a multitude of mutations, including drivers responsible for cancer progression and neutral ...passengers. Measuring selection in cancer and distinguishing drivers from passengers have important implications for development of novel treatment strategies. It has recently been argued that a third of cancers are evolving neutrally, as their mutational frequency spectrum follows a 1/f power law expected from neutral evolution in a particular intermediate frequency range. We study a stochastic model of cancer evolution and derive a formula for the probability distribution of the cancer cell frequency of a subclonal driver, demonstrating that driver frequency is biased towards 0 and 1. We show that it is difficult to capture a driver mutation at an intermediate frequency, and thus the calling of neutrality due to a lack of such driver will significantly overestimate the number of neutrally evolving tumors. Our approach provides quantification of the validity of the 1/f statistic across the entire range of relevant parameter values. We also show that our conclusions remain valid for non-exponential models: spatial 3d model and sigmoidal growth, relevant for early- and late stages of cancer growth.
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
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.
Evolutionary pathways describe trajectories of biological evolution in the space of different variants of organisms (genotypes). The probability of existence and the number of evolutionary pathways ...that lead from a given genotype to a better-adapted genotype are important measures of accessibility of local fitness optima and the reproducibility of evolution. Both quantities have been studied in simple mathematical models where genotypes are represented as binary sequences of two types of basic units, and the network of permitted mutations between the genotypes is a hypercube graph. However, it is unclear how these results translate to the biologically relevant case in which genotypes are represented by sequences of more than two units, for example four nucleotides (DNA) or 20 amino acids (proteins), and the mutational graph is not the hypercube. Here we investigate accessibility of the best-adapted genotype in the general case of K > 2 units. Using computer generated and experimental fitness landscapes we show that accessibility of the global fitness maximum increases with K and can be much higher than for binary sequences. The increase in accessibility comes from the increase in the number of indirect trajectories exploited by evolution for higher K. As one of the consequences, the fraction of genotypes that are accessible increases by three orders of magnitude when the number of units K increases from 2 to 16 for landscapes of size N ∼ 106 genotypes. This suggests that evolution can follow many different trajectories on such landscapes and the reconstruction of evolutionary pathways from experimental data might be an extremely difficult task.
The rate of adaptation in theoretical models of biological evolution generally increases with the mutation rate. However, there is a threshold beyond which mutations into lethal states lead to ...extinction. It would be logical to assume that eliminating such lethal states could be advantageous for evolution. Here, we demonstrate that lethal mutations actually accelerate adaptation on rugged fitness landscapes with multiple peaks and valleys in the presence of competition for resources. We investigate a modified stochastic version of the quasispecies model, incorporating two types of genotypes—viable and lethal—and show that higher rates of lethal mutations result in shorter evolution times towards the best-fit genotype. This phenomenon can be attributed to an increased frequency of traversing fitness valleys, facilitated by reduced selection pressure against less-fit variants.
Stochastic phenotype switching has been suggested to play a beneficial role in microbial populations by leading to the division of labour among cells, or ensuring that at least some of the population ...survives an unexpected change in environmental conditions. Here we use a computational model to investigate an alternative possible function of stochastic phenotype switching: as a way to adapt more quickly even in a static environment. We show that when a genetic mutation causes a population to become less fit, switching to an alternative phenotype with higher fitness (growth rate) may give the population enough time to develop compensatory mutations that increase the fitness again. The possibility of switching phenotypes can reduce the time to adaptation by orders of magnitude if the "fitness valley" caused by the deleterious mutation is deep enough. Our work has important implications for the emergence of antibiotic-resistant bacteria. In line with recent experimental findings, we hypothesise that switching to a slower growing - but less sensitive - phenotype helps bacteria to develop resistance by providing alternative, faster evolutionary routes to resistance.
Cancer evolution is predominantly studied by focusing on differences in the genetic characteristics of malignant cells within tumors. However, the spatiotemporal dynamics of clonal outgrowth that ...underlie evolutionary trajectories remain largely unresolved. Here, we sought to unravel the clonal dynamics of colorectal cancer (CRC) expansion in space and time by using a color-based clonal tracing method. This method involves lentiviral red-green-blue (RGB) marking of cell populations, which enabled us to track individual cells and their clonal outgrowth during tumor initiation and growth in a xenograft model. We found that clonal expansion largely depends on the location of a clone, as small clones reside in the center and large clones mostly drive tumor growth at the border. These dynamics are recapitulated in a computational model, which confirms that the clone position within a tumor rather than cell-intrinsic features, is crucial for clonal outgrowth. We also found that no significant clonal loss occurs during tumor growth and clonal dispersal is limited in most models. Our results imply that, in addition to molecular features of clones such as (epi-)genetic differences between cells, clone location and the geometry of tumor growth are crucial for clonal expansion. Our findings suggest that either microenvironmental signals on the tumor border or differences in physical properties within the tumor, are major contributors to explain heterogeneous clonal expansion. Thus, this study provides further insights into the dynamics of solid tumor growth and progression, as well as the origins of tumor cell heterogeneity in a relevant model system.
Phenotypic delay-the time delay between genetic mutation and expression of the corresponding phenotype-is generally neglected in evolutionary models, yet recent work suggests that it may be more ...common than previously assumed. Here, we use computer simulations and theory to investigate the significance of phenotypic delay for the evolution of bacterial resistance to antibiotics. We consider three mechanisms which could potentially cause phenotypic delay: effective polyploidy, dilution of antibiotic-sensitive molecules and accumulation of resistance-enhancing molecules. We find that the accumulation of resistant molecules is relevant only within a narrow parameter range, but both the dilution of sensitive molecules and effective polyploidy can cause phenotypic delay over a wide range of parameters. We further investigate whether these mechanisms could affect population survival under drug treatment and thereby explain observed discrepancies in mutation rates estimated by Luria-Delbrück fluctuation tests. While the effective polyploidy mechanism does not affect population survival, the dilution of sensitive molecules leads both to decreased probability of survival under drug treatment and underestimation of mutation rates in fluctuation tests. The dilution mechanism also changes the shape of the Luria-Delbrück distribution of mutant numbers, and we show that this modified distribution provides an improved fit to previously published experimental data.