Cell jamming in a collagen-based interface assay is tuned by collagen density and proteolysis.

Tumor cell invasion into heterogenous interstitial tissues consisting of network-, channel- or rift-like architectures involves both matrix metalloproteinase (MMP)-mediated tissue remodeling and cell shape adaptation to tissue geometry. Three-dimensional (3D) models composed of either porous or linearly aligned architectures have added to the understanding of how physical spacing principles affect migration efficacy; however, the relative contribution of each architecture to decision making in the presence of varying MMP availability is not known. Here, we developed an interface assay containing a cleft between two high-density collagen lattices, and we used this assay to probe tumor cell invasion efficacy, invasion mode and MMP dependence in concert. In silico modeling predicted facilitated cell migration into confining clefts independently of MMP activity, whereas migration into dense porous matrix was predicted to require matrix degradation. This prediction was verified experimentally, where inhibition of collagen degradation was found to strongly compromise migration into 3D collagen in a density-dependent manner, but interface-guided migration remained effective, occurring by cell jamming. The 3D interface assay reported here may serve as a suitable model to better understand the impact of in vivo-relevant interstitial tissue topologies on tumor invasion patterning and responses to molecular interventions.

Evolution of phenotypic plasticity leads to tumor heterogeneity with implications for therapy.

Cancer is a significant global health issue, with treatment challenges arising from intratumor heterogeneity. This heterogeneity stems mainly from somatic evolution, causing genetic diversity within the tumor, and phenotypic plasticity of tumor cells leading to reversible phenotypic changes. However, the interplay of both factors has not been rigorously investigated. Here, we examine the complex relationship between somatic evolution and phenotypic plasticity, explicitly focusing on the interplay between cell migration and proliferation. This type of phenotypic plasticity is essential in glioblastoma, the most aggressive form of brain tumor. We propose that somatic evolution alters the regulation of phenotypic plasticity in tumor cells, specifically the reaction to changes in the microenvironment. We study this hypothesis using a novel, spatially explicit model that tracks individual cells' phenotypic and genetic states. We assume cells change between migratory and proliferative states controlled by inherited and mutation-driven genotypes and the cells' microenvironment. We observe that cells at the tumor edge evolve to favor migration over proliferation and vice versa in the tumor bulk. Notably, different genetic configurations can result in this pattern of phenotypic heterogeneity. We analytically predict the outcome of the evolutionary process, showing that it depends on the tumor microenvironment. Synthetic tumors display varying levels of genetic and phenotypic heterogeneity, which we show are predictors of tumor recurrence time after treatment. Interestingly, higher phenotypic heterogeneity predicts poor treatment outcomes, unlike genetic heterogeneity. Our research offers a novel explanation for heterogeneous patterns of tumor recurrence in glioblastoma patients.

BIO-LGCA: A cellular automaton modelling class for analysing collective cell migration.

Collective dynamics in multicellular systems such as biological organs and tissues plays a key role in biological development, regeneration, and pathological conditions. Collective tissue dynamics-understood as population behaviour arising from the interplay of the constituting discrete cells-can be studied with on- and off-lattice agent-based models. However, classical on-lattice agent-based models, also known as cellular automata, fail to replicate key aspects of collective migration, which is a central instance of collective behaviour in multicellular systems. To overcome drawbacks of classical on-lattice models, we introduce an on-lattice, agent-based modelling class for collective cell migration, which we call biological lattice-gas cellular automaton (BIO-LGCA). The BIO-LGCA is characterised by synchronous time updates, and the explicit consideration of individual cell velocities. While rules in classical cellular automata are typically chosen ad hoc, rules for cell-cell and cell-environment interactions in the BIO-LGCA can also be derived from experimental cell migration data or biophysical laws for individual cell migration. We introduce elementary BIO-LGCA models of fundamental cell interactions, which may be combined in a modular fashion to model complex multicellular phenomena. We exemplify the mathematical mean-field analysis of specific BIO-LGCA models, which allows to explain collective behaviour. The first example predicts the formation of clusters in adhesively interacting cells. The second example is based on a novel BIO-LGCA combining adhesive interactions and alignment. For this model, our analysis clarifies the nature of the recently discovered invasion plasticity of breast cancer cells in heterogeneous environments.

Correction: Quantification of nematic cell polarity in three-dimensional tissues.

[This corrects the article DOI: 10.1371/journal.pcbi.1008412.].

Quantification of nematic cell polarity in three-dimensional tissues.

How epithelial cells coordinate their polarity to form functional tissues is an open question in cell biology. Here, we characterize a unique type of polarity found in liver tissue, nematic cell polarity, which is different from vectorial cell polarity in simple, sheet-like epithelia. We propose a conceptual and algorithmic framework to characterize complex patterns of polarity proteins on the surface of a cell in terms of a multipole expansion. To rigorously quantify previously observed tissue-level patterns of nematic cell polarity (Morales-Navarrete et al., eLife 2019), we introduce the concept of co-orientational order parameters, which generalize the known biaxial order parameters of the theory of liquid crystals. Applying these concepts to three-dimensional reconstructions of single cells from high-resolution imaging data of mouse liver tissue, we show that the axes of nematic cell polarity of hepatocytes exhibit local coordination and are aligned with the biaxially anisotropic sinusoidal network for blood transport. Our study characterizes liver tissue as a biological example of a biaxial liquid crystal. The general methodology developed here could be applied to other tissues and in-vitro organoids.

Inferring the effect of interventions on COVID-19 transmission networks.

Countries around the world implement nonpharmaceutical interventions (NPIs) to mitigate the spread of COVID-19. Design of efficient NPIs requires identification of the structure of the disease transmission network. We here identify the key parameters of the COVID-19 transmission network for time periods before, during, and after the application of strict NPIs for the first wave of COVID-19 infections in Germany combining Bayesian parameter inference with an agent-based epidemiological model. We assume a Watts-Strogatz small-world network which allows to distinguish contacts within clustered cliques and unclustered, random contacts in the population, which have been shown to be crucial in sustaining the epidemic. In contrast to other works, which use coarse-grained network structures from anonymized data, like cell phone data, we consider the contacts of individual agents explicitly. We show that NPIs drastically reduced random contacts in the transmission network, increased network clustering, and resulted in a previously unappreciated transition from an exponential to a constant regime of new cases. In this regime, the disease spreads like a wave with a finite wave speed that depends on the number of contacts in a nonlinear fashion, which we can predict by mean field theory.

Cell-cell adhesion and 3D matrix confinement determine jamming transitions in breast cancer invasion.

Plasticity of cancer invasion and metastasis depends on the ability of cancer cells to switch between collective and single-cell dissemination, controlled by cadherin-mediated cell-cell junctions. In clinical samples, E-cadherin-expressing and -deficient tumours both invade collectively and metastasize equally, implicating additional mechanisms controlling cell-cell cooperation and individualization. Here, using spatially defined organotypic culture, intravital microscopy of mammary tumours in mice and in silico modelling, we identify cell density regulation by three-dimensional tissue boundaries to physically control collective movement irrespective of the composition and stability of cell-cell junctions. Deregulation of adherens junctions by downregulation of E-cadherin and p120-catenin resulted in a transition from coordinated to uncoordinated collective movement along extracellular boundaries, whereas single-cell escape depended on locally free tissue space. These results indicate that cadherins and extracellular matrix confinement cooperate to determine unjamming transitions and stepwise epithelial fluidization towards, ultimately, cell individualization.