A novel jamming phase diagram links tumor invasion to non-equilibrium phase separation.

It is well established that the early malignant tumor invades surrounding extracellular matrix (ECM) in a manner that depends upon material properties of constituent cells, surrounding ECM, and their interactions. Recent studies have established the capacity of the invading tumor spheroids to evolve into coexistent solid-like, fluid-like, and gas-like phases. Using breast cancer cell lines invading into engineered ECM, here we show that the spheroid interior develops spatial and temporal heterogeneities in material phase which, depending upon cell type and matrix density, ultimately result in a variety of phase separation patterns at the invasive front. Using a computational approach, we further show that these patterns are captured by a novel jamming phase diagram. We suggest that non-equilibrium phase separation based upon jamming and unjamming transitions may provide a unifying physical picture to describe cellular migratory dynamics within, and invasion from, a tumor.

Geometric constraints during epithelial jamming.

As an injury heals, an embryo develops, or a carcinoma spreads, epithelial cells systematically change their shape. In each of these processes cell shape is studied extensively whereas variability of shape from cell-to-cell is regarded most often as biological noise. But where do cell shape and its variability come from? Here we report that cell shape and shape variability are mutually constrained through a relationship that is purely geometrical. That relationship is shown to govern processes as diverse as maturation of the pseudostratified bronchial epithelial layer cultured from non-asthmatic or asthmatic donors, and formation of the ventral furrow in the Drosophila embryo. Across these and other epithelial systems, shape variability collapses to a family of distributions that is common to all. That distribution, in turn, is accounted for by a mechanistic theory of cell-cell interaction showing that cell shape becomes progressively less elongated and less variable as the layer becomes progressively more jammed. These findings suggest a connection between jamming and geometry that spans living organisms and inert jammed systems, and thus transcends system details. Although molecular events are needed for any complete theory of cell shape and cell packing, observations point to the hypothesis that jamming behavior at larger scales of organization sets overriding geometrical constraints.

Scaling Physiologic Function from Cell to Tissue in Asthma, Cancer, and Development.

The formation of an integrated tissue from individual cells depends on the properties of the individual cells as well as the interaction of many cells acting as a collective. Three fundamental physiological processes govern the collective scaling from the individual cell to a working tissue: cell sorting, tissue assembly, and collective cellular migration. Mechanistically, cell sorting is governed by differential adhesion, whereas tissue assembly is controlled by the epithelial-to-mesenchymal transition and its inverse, the mesenchymal-to-epithelial transition. The mechanism driving collective cellular migration, however, is not clear. To fill that gap, here we consider cell jamming and unjamming, and their role in collective cellular migration.

Contact guidance and collective migration in the advancing epithelial monolayer.

At the edge of a confluent cell layer, cell-free empty space is a cue that can drive directed collective cellular migration. Similarly, contact guidance is also a robust mechanical cue that can drive cell migration. However, it is unclear which of the two effects is stronger, and how each mechanism affects collective migration. To address this question, here we explore the trajectories of cells migrating collectively on a substrate containing micropatterned grooves (10-20 µm in periodicity, 2 µm in height) compared with unpatterned control substrates. Compared with unpatterned controls, the micropatterned substrates attenuated path variance by close to 70% and augmented migration coordination by more than 30%. Together, these results show that contact guidance can play an appreciable role in collective cellular migration. Also, our result can provide insights into tissue repair and regeneration with the remodeling of the connective tissue matrix.

Collective motion of mammalian cell cohorts in 3D.

Collective cell migration is ubiquitous in biology, from development to cancer; it occurs in complex systems comprised of heterogeneous cell types, signals and matrices, and requires large scale regulation in space and time. Understanding how cells achieve organized collective motility is crucial to addressing cellular and tissue function and disease progression. While current two-dimensional model systems recapitulate the dynamic properties of collective cell migration, quantitative three-dimensional equivalent model systems have proved elusive. To establish such a model system, we study cell collectives by tracking individuals within cell cohorts embedded in three dimensional collagen scaffolding. We develop a custom algorithm to quantify the temporal and spatial heterogeneity of motion in cell cohorts during motility events. In the absence of external driving agents, we show that these cohorts rotate in short bursts, <2 hours, and translate for up to 6 hours. We observe, track, and analyze three dimensional motion of cell cohorts composed of 3-31 cells, and pave a path toward understanding cell collectives in 3D as a complex emergent system.

A quantitative comparison of human HT-1080 fibrosarcoma cells and primary human dermal fibroblasts identifies a 3D migration mechanism with properties unique to the transformed phenotype.

Here, we describe an engineering approach to quantitatively compare migration, morphologies, and adhesion for tumorigenic human fibrosarcoma cells (HT-1080s) and primary human dermal fibroblasts (hDFs) with the aim of identifying distinguishing properties of the transformed phenotype. Relative adhesiveness was quantified using self-assembled monolayer (SAM) arrays and proteolytic 3-dimensional (3D) migration was investigated using matrix metalloproteinase (MMP)-degradable poly(ethylene glycol) (PEG) hydrogels ("synthetic extracellular matrix" or "synthetic ECM"). In synthetic ECM, hDFs were characterized by vinculin-containing features on the tips of protrusions, multipolar morphologies, and organized actomyosin filaments. In contrast, HT-1080s were characterized by diffuse vinculin expression, pronounced ß1-integrin on the tips of protrusions, a cortically-organized F-actin cytoskeleton, and quantitatively more rounded morphologies, decreased adhesiveness, and increased directional motility compared to hDFs. Further, HT-1080s were characterized by contractility-dependent motility, pronounced blebbing, and cortical contraction waves or constriction rings, while quantified 3D motility was similar in matrices with a wide range of biochemical and biophysical properties (including collagen) despite substantial morphological changes. While HT-1080s were distinct from hDFs for each of the 2D and 3D properties investigated, several features were similar to WM239a melanoma cells, including rounded, proteolytic migration modes, cortical F-actin organization, and prominent uropod-like structures enriched with ß1-integrin, F-actin, and melanoma cell adhesion molecule (MCAM/CD146/MUC18). Importantly, many of the features observed for HT-1080s were analogous to cellular changes induced by transformation, including cell rounding, a disorganized F-actin cytoskeleton, altered organization of focal adhesion proteins, and a weakly adherent phenotype. Based on our results, we propose that HT-1080s migrate in synthetic ECM with functional properties that are a direct consequence of their transformed phenotype.

Decoupling directed and passive motion in dynamic systems: particle tracking microrheology of sputum.

Probing the physical properties of heterogeneous materials is essential to understand the structure, function and dynamics of complex fluids including cells, mucus, and polymer solutions. Particle tracking microrheology is a useful method to passively probe viscoelastic properties on micron length scales by tracking the thermal motion of beads embedded in the sample. However, errors associated with active motion have limited the implementation to dynamic systems. We present a simple method to decouple active and Brownian motion, enabling particle tracking to be applied to fluctuating heterogeneous systems. We use the movement perpendicular to the major axis of motion in time to calculate rheological properties. Through simulated data we demonstrate that this method removes directed motion and performs equally well when there is no directed motion, with an average percent error of <1%. We use this method to measure glycerol-water mixtures to show the capability to measure a range of materials. Finally, we use this technique to characterize the compliance of human sputum. We also investigate the effect of a liquefaction agent used to prepare sputum for diagnostic purposes. Our results suggest that the addition of high concentration sodium hydroxide increases sample heterogeneity by increasing the maximum observed creep compliance.

Nanomechanics of multiple units in the erythrocyte membrane skeletal network.

Erythrocytes undergo deformations when they transport O(2) and CO(2) across the membrane, yet the 3D nanomechanics of the skeletal network remains poorly understood. Expanding from our previous single isolated unit, we now simulate networks consisting of 1-10 concentric rings of repeating units in equibiaxial deformation. The networks are organized with (1) a 3D model for a single unit, (2) a wrap-around mode between Sp and actin protofilament in the intra-unit interaction, and (3) a random inter-unit connectivity. These assumptions permit efficient five-degrees-of-freedom (5DOF) simulations when up to 30 pN of radial forces are applied to the boundary spectrin (Sp) and the center and other units are analyzed. As 6 Sp balance their tensions, hexagonal units become irregular. While actin protofilaments remain tangent to the network, their yaw (Phi) angles change drastically with addition of neighboring units or an Sp unfolding. It is anticipated that during deformation, transmembrane complexes associated with the network move laterally through the lipid bilayer and increase the diffusion of molecules across the membrane. When protofilament/Sp sweeps under the lipid bilayer, they mix up the submembrane concentration gradient. Thus, the nanomechanics of actin protofilaments and Sp may enhance the transport of molecules during erythrocyte deformation.