Impact of crowding on the diversity of expanding populations.

Crowding effects critically impact the self-organization of densely packed cellular assemblies, such as biofilms, solid tumors, and developing tissues. When cells grow and divide, they push each other apart, remodeling the structure and extent of the population's range. Recent work has shown that crowding has a strong impact on the strength of natural selection. However, the impact of crowding on neutral processes, which controls the fate of new variants as long as they are rare, remains unclear. Here, we quantify the genetic diversity of expanding microbial colonies and uncover signatures of crowding in the site frequency spectrum. By combining Luria-Delbrück fluctuation tests, lineage tracing in a novel microfluidic incubator, cell-based simulations, and theoretical modeling, we find that the majority of mutations arise behind the expanding frontier, giving rise to clones that are mechanically "pushed out" of the growing region by the proliferating cells in front. These excluded-volume interactions result in a clone-size distribution that solely depends on where the mutation first arose relative to the front and is characterized by a simple power law for low-frequency clones. Our model predicts that the distribution depends on a single parameter-the characteristic growth layer thickness-and hence allows estimation of the mutation rate in a variety of crowded cellular populations. Combined with previous studies on high-frequency mutations, our finding provides a unified picture of the genetic diversity in expanding populations over the whole frequency range and suggests a practical method to assess growth dynamics by sequencing populations across spatial scales.

Evolutionary rescue of resistant mutants is governed by a balance between radial expansion and selection in compact populations.

Mutation-mediated treatment resistance is one of the primary challenges for modern antibiotic and anti-cancer therapy. Yet, many resistance mutations have a substantial fitness cost and are subject to purifying selection. How emerging resistant lineages may escape purifying selection via subsequent compensatory mutations is still unclear due to the difficulty of tracking such evolutionary rescue dynamics in space and time. Here, we introduce a system of fluorescence-coupled synthetic mutations to show that the probability of evolutionary rescue, and the resulting long-term persistence of drug resistant mutant lineages, is dramatically increased in dense microbial populations. By tracking the entire evolutionary trajectory of thousands of resistant lineages in expanding yeast colonies we uncover an underlying quasi-stable equilibrium between the opposing forces of radial expansion and natural selection, a phenomenon we term inflation-selection balance. Tailored computational models and agent-based simulations corroborate the fundamental nature of the observed effects and demonstrate the potential impact on drug resistance evolution in cancer. The described phenomena should be considered when predicting multi-step evolutionary dynamics in any mechanically compact cellular population, including pathogenic microbial biofilms and solid tumors. The insights gained will be especially valuable for the quantitative understanding of response to treatment, including emerging evolution-based therapy strategies.

Collective motion conceals fitness differences in crowded cellular populations.

Many cellular populations are tightly packed, such as microbial colonies and biofilms, or tissues and tumours in multicellular organisms. The movement of one cell in these crowded assemblages requires motion of others, so that cell displacements are correlated over many cell diameters. Whenever movement is important for survival or growth, these correlated rearrangements could couple the evolutionary fate of different lineages. However, little is known about the interplay between mechanical forces and evolution in dense cellular populations. Here, by tracking slower-growing clones at the expanding edge of yeast colonies, we show that the collective motion of cells prevents costly mutations from being weeded out rapidly. Joint pushing by neighbouring cells generates correlated movements that suppress the differential displacements required for selection to act. This mechanical screening of fitness differences allows slower-growing mutants to leave more descendants than expected under non-mechanical models, thereby increasing their chance for evolutionary rescue. Our work suggests that, in crowded populations, cells cooperate with surrounding neighbours through inevitable mechanical interactions. This effect has to be considered when predicting evolutionary outcomes, such as the emergence of drug resistance or cancer evolution.

Emergence of evolutionary driving forces in pattern-forming microbial populations.

Evolutionary dynamics are controlled by a number of driving forces, such as natural selection, random genetic drift and dispersal. In this perspective article, we aim to emphasize that these forces act at the population level, and that it is a challenge to understand how they emerge from the stochastic and deterministic behaviour of individual cells. Even the most basic steric interactions between neighbouring cells can couple evolutionary outcomes of otherwise unrelated individuals, thereby weakening natural selection and enhancing random genetic drift. Using microbial examples of varying degrees of complexity, we demonstrate how strongly cell-cell interactions influence evolutionary dynamics, especially in pattern-forming systems. As pattern formation itself is subject to evolution, we propose to study the feedback between pattern formation and evolutionary dynamics, which could be key to predicting and potentially steering evolutionary processes. Such an effort requires extending the systems biology approach from the cellular to the population scale.This article is part of the theme issue 'Self-organization in cell biology'.

SCWISh network is essential for survival under mechanical pressure.

Cells that proliferate within a confined environment build up mechanical compressive stress. For example, mechanical pressure emerges in the naturally space-limited tumor environment. However, little is known about how cells sense and respond to mechanical compression. We developed microfluidic bioreactors to enable the investigation of the effects of compressive stress on the growth of the genetically tractable model organism Saccharomyces cerevisiae We used this system to determine that compressive stress is partly sensed through a module consisting of the mucin Msb2 and the cell wall protein Sho1, which act together as a sensor module in one of the two major osmosensing pathways in budding yeast. This signal is transmitted via the MAPKKK kinase Ste11. Thus, we term this mechanosensitive pathway the "SMuSh" pathway, for Ste11 through Mucin/Sho1 pathway. The SMuSh pathway delays cells in the G1 phase of the cell cycle and improves cell survival in response to growth-induced pressure. We also found that the cell wall integrity (CWI) pathway contributes to the response to mechanical compressive stress. These latter results are confirmed in complimentary experiments in Mishra et al. [Mishra R, et al. (2017) Proc Natl Acad Sci USA, 10.1073/pnas.1709079114]. When both the SMuSh and the CWI pathways are deleted, cells fail to adapt to compressive stress, and all cells lyse at relatively low pressure when grown in confinement. Thus, we define a network that is essential for cell survival during growth under pressure. We term this mechanosensory system the SCWISh (survival through the CWI and SMuSh) network.

Excess of mutational jackpot events in expanding populations revealed by spatial Luria-Delbrück experiments.

The genetic diversity of growing cellular populations, such as biofilms, solid tumours or developing embryos, is thought to be dominated by rare, exceptionally large mutant clones. Yet, the emergence of these mutational jackpot events is only understood in well-mixed populations, where they stem from mutations that arise during the first few cell divisions. To study jackpot events in spatially structured populations, we track mutant clones in microbial populations using fluorescence microscopy and population sequencing. High-frequency mutations are found to be massively enriched in microbial colonies compared with well-shaken liquid cultures, as a result of late-occurring mutations surfing at the edge of range expansions. Thus, jackpot events can be generated not only when mutations arise early but also when they occur at favourable locations, which exacerbates their role in adaptation and disease. In particular, because spatial competition with the wild type keeps most mutant clones in a quiescent state, strong selection pressures that kill the wild type promote drug resistance.

The small heat shock protein Hsp27 affects assembly dynamics and structure of keratin intermediate filament networks.

The mechanical properties of living cells are essential for many processes. They are defined by the cytoskeleton, a composite network of protein fibers. Thus, the precise control of its architecture is of paramount importance. Our knowledge about the molecular and physical mechanisms defining the network structure remains scarce, especially for the intermediate filament cytoskeleton. Here, we investigate the effect of small heat shock proteins on the keratin 8/18 intermediate filament cytoskeleton using a well-controlled model system of reconstituted keratin networks. We demonstrate that Hsp27 severely alters the structure of such networks by changing their assembly dynamics. Furthermore, the C-terminal tail domain of keratin 8 is shown to be essential for this effect. Combining results from fluorescence and electron microscopy with data from analytical ultracentrifugation reveals the crucial role of kinetic trapping in keratin network formation.

Neurofilament sidearms modulate parallel and crossed-filament orientations inducing nematic to isotropic and re-entrant birefringent hydrogels.

Neurofilaments are intermediate filaments assembled from the subunits neurofilament-low, neurofilament-medium and neurofilament-high. In axons, parallel neurofilaments form a nematic liquid-crystal hydrogel with network structure arising from interactions between the neurofilaments' C-terminal sidearms. Here we report, using small-angle X-ray-scattering, polarized-microscopy and rheometry, that with decreasing ionic strength, neurofilament-low-high, neurofilament-low-medium and neurofilament-low-medium-high hydrogels transition from the nematic hydrogel to an isotropic hydrogel (with random, crossed-filament orientation) and to an unexpected new re-entrant liquid-crystal hydrogel with parallel filaments--the bluish-opaque hydrogel--with notable mechanical and water retention properties reminiscent of crosslinked hydrogels. Significantly, the isotropic gel phase stability is sidearm-dependent: neurofilament-low-high hydrogels exhibit a wide ionic strength range, neurofilament-low-medium hydrogels a narrow ionic strength range, whereas neurofilament-low hydrogels lack the isotropic gel phase. This suggests a dominant regulatory role for neurofilament-high sidearms in filament reorientation plasticity, facilitating organelle transport in axons. Neurofilament-inspired biomimetic hydrogels should therefore exhibit remarkable structure-dependent moduli and slow and fast water-release properties.

Towards constructing extracellular matrix-mimetic hydrogels: an elastic hydrogel constructed from tandem modular proteins containing tenascin FnIII domains.

Protein-based hydrogels have been developed for various biomedical applications where they provide artificial extracellular microenvironments that mimic the physical and biochemical characteristics of natural extracellular matrices (ECMs). In natural ECMs, a large number of proteins are tandem modular proteins consisting of many individually folded functional domains that confer structural and biological functionalities. Such tandem modular proteins are promising building blocks for constructing ECM-mimetic biomaterials. However, their use for such purposes has not been explored extensively. Tenascin-C (TNC) is an ECM tandem modular protein and plays an important role in mechanotransduction by regulating important cell-matrix interactions. The third FnIII domain of TNC (TNfn3) contains an RGD sequence and is known to bind integrins. Here we use the TNfn3 domain and resilin sequence-based tandem modular protein FRF4RF4R (F represents the TNfn3 domain and R represents the resilin sequence, respectively) as a building block to construct protein-based ECM-mimetic hydrogels. The tandem modular protein-based building block FRF4RF4R closely mimics the architecture of the naturally occurring tandem modular ECM protein TNC and incorporates intact RGD-containing FnIII domains. Our results demonstrate that tandem modular proteins containing TNfn3 can be readily photochemically crosslinked into elastic hydrogels, whose Young's modulus can be tuned by the concentration of the tandem modular protein solution. In vitro studies demonstrate that none of the photochemical crosslinking reaction components are cytotoxic at the level tested, and the hydrogel supports the spread of human lung fibroblast cells. Our results demonstrate that FRF4RF4R-based hydrogel is a novel ECM-mimetic hydrogel.

Stability of a surface-bound oligonucleotide duplex inferred from molecular dynamics: a study of single nucleotide defects using DNA microarrays.

Microarray technology uses the sequence dependent hybridization (binding) affinity of surface-bound oligonucleotide strands for the quantification of complex nucleic acid mixtures. In spite of its huge potential in life science and medicine, microarray oligonucleotide hybridization remains far from being understood. Taking advantage of microarray combinatorial possibilities we show that, although surface bound, the hybridization affinities of single-base mismatched oligonucleotides can be derived from first principles using parameters from bulk.

Position dependent mismatch discrimination on DNA microarrays - experiments and model.

BACKGROUND: The propensity of oligonucleotide strands to form stable duplexes with complementary sequences is fundamental to a variety of biological and biotechnological processes as various as microRNA signalling, microarray hybridization and PCR. Yet our understanding of oligonucleotide hybridization, in particular in presence of surfaces, is rather limited. Here we use oligonucleotide microarrays made in-house by optically controlled DNA synthesis to produce probe sets comprising all possible single base mismatches and base bulges for each of 20 sequence motifs under study. RESULTS: We observe that mismatch discrimination is mostly determined by the defect position (relative to the duplex ends) as well as by the sequence context. We investigate the thermodynamics of the oligonucleotide duplexes on the basis of double-ended molecular zipper. Theoretical predictions of defect positional influence as well as long range sequence influence agree well with the experimental results. CONCLUSION: Molecular zipping at thermodynamic equilibrium explains the binding affinity of mismatched DNA duplexes on microarrays well. The position dependent nearest neighbor model (PDNN) can be inferred from it. Quantitative understanding of microarray experiments from first principles is in reach.

Impact of point-mutations on the hybridization affinity of surface-bound DNA/DNA and RNA/DNA oligonucleotide-duplexes: comparison of single base mismatches and base bulges.

BACKGROUND: The high binding specificity of short 10 to 30 mer oligonucleotide probes enables single base mismatch (MM) discrimination and thus provides the basis for genotyping and resequencing microarray applications. Recent experiments indicate that the underlying principles governing DNA microarray hybridization - and in particular MM discrimination - are not completely understood. Microarrays usually address complex mixtures of DNA targets. In order to reduce the level of complexity and to study the problem of surface-based hybridization with point defects in more detail, we performed array based hybridization experiments in well controlled and simple situations. RESULTS: We performed microarray hybridization experiments with short 16 to 40 mer target and probe lengths (in situations without competitive hybridization) in order to systematically investigate the impact of point-mutations - varying defect type and position - on the oligonucleotide duplex binding affinity. The influence of single base bulges and single base MMs depends predominantly on position - it is largest in the middle of the strand. The position-dependent influence of base bulges is very similar to that of single base MMs, however certain bulges give rise to an unexpectedly high binding affinity. Besides the defect (MM or bulge) type, which is the second contribution in importance to hybridization affinity, there is also a sequence dependence, which extends beyond the defect next-neighbor and which is difficult to quantify. Direct comparison between binding affinities of DNA/DNA and RNA/DNA duplexes shows, that RNA/DNA purine-purine MMs are more discriminating than corresponding DNA/DNA MMs. In DNA/DNA MM discrimination the affected base pair (C.G vs. A.T) is the pertinent parameter. We attribute these differences to the different structures of the duplexes (A vs. B form). CONCLUSION: We have shown that DNA microarrays can resolve even subtle changes in hybridization affinity for simple target mixtures. We have further shown that the impact of point defects on oligonucleotide stability can be broken down to a hierarchy of effects. In order to explain our observations we propose DNA molecular dynamics - in form of zipping of the oligonucleotide duplex - to play an important role.