Spatially heterogeneous dynamics of cells in a growing tumor spheroid: comparison between theory and experiments

Emergence of cellular nematic order is a conserved feature of gastrulation in animal embryos

Cells undergo dramatic changes in morphology during embryogenesis, yet how these changes affect the formation of ordered tissues remains elusive. Here we find that the emergence of a nematic liquid crystal phase occurs in cells during gastrulation in the development of embryos of fish, frogs, and fruit flies. Moreover, the spatial correlations in all three organisms are long-ranged and follow a similar power-law decay( y ∼ x - α ) with α less than unity for the nematic order parameter, suggesting a common underlying physical mechanism unifies events in these distantly related species. All three species exhibit similar propagation of the nematic phase, reminiscent of nucleation and growth phenomena. Finally, we use a theoretical model along with disruptions of cell adhesion and cell specification to characterize the minimal features required for formation of the nematic phase. Our results provide a framework for understanding a potentially universal features of metazoan embryogenesis and shed light on the advent of ordered structures during animal development.

Contact area and tissue growth dynamics shape synthetic juxtacrine signaling patterns

Cell-cell communication through direct contact, or juxtacrine signaling, is important in development, disease, and many areas of physiology. Synthetic forms of juxtacrine signaling can be precisely controlled and operate orthogonally to native processes, making them a powerful reductionist tool with which to address fundamental questions in cell-cell communication in vivo. Here we investigate how cell-cell contact length and tissue growth dynamics affect juxtacrine signal responses through implementing a custom synthetic gene circuit in Drosophila wing imaginal discs alongside mathematical modeling to determine synthetic Notch (synNotch) activation patterns. We find that the area of contact between cells largely determines the extent of synNotch activation, leading to the prediction that the shape of the interface between signal-sending and signal-receiving cells will impact the magnitude of the synNotch response. Notably, synNotch outputs form a graded spatial profile that extends several cell diameters from the signal source, providing evidence that the response to juxtacrine signals can persist in cells as they proliferate away from source cells, or that cells remain able to communicate directly over several cell diameters. Our model suggests the former mechanism may be sufficient, since it predicts graded outputs without diffusion or long-range cell-cell communication. Overall, we identify that cell-cell contact area together with output synthesis and decay rates likely govern the pattern of synNotch outputs in both space and time during tissue growth, insights that may have broader implications for juxtacrine signaling in general.

3D Hydrogel Encapsulation Regulates Nephrogenesis in Kidney Organoids

Stem cell-derived kidney organoids contain nephron segments that recapitulate morphological and functional aspects of the human kidney. However, directed differentiation protocols for kidney organoids are largely conducted using biochemical signals to control differentiation. Here, the hypothesis that mechanical signals regulate nephrogenesis is investigated in 3D culture by encapsulating kidney organoids within viscoelastic alginate hydrogels with varying rates of stress relaxation. Tubular nephron segments are significantly more convoluted in kidney organoids differentiated in encapsulating hydrogels when compared with those in suspension culture. Hydrogel viscoelasticity regulates the spatial distribution of nephron segments within the differentiating kidney organoids. Consistent with these observations, a particle-based computational model predicts that the extent of deformation of the hydrogel-organoid interface regulates the morphology of nephron segments. Elevated extracellular calcium levels in the culture medium, which can be impacted by the hydrogels, decrease the glomerulus-to-tubule ratio of nephron segments. These findings reveal that hydrogel encapsulation regulates nephron patterning and morphology and suggest that the mechanical microenvironment is an important design variable for kidney regenerative medicine.

Modeling the extracellular matrix in cell migration and morphogenesis: a guide for the curious biologist

The extracellular matrix (ECM) is a highly complex structure through which biochemical and mechanical signals are transmitted. In processes of cell migration, the ECM also acts as a scaffold, providing structural support to cells as well as points of potential attachment. Although the ECM is a well-studied structure, its role in many biological processes remains difficult to investigate comprehensively due to its complexity and structural variation within an organism. In tandem with experiments, mathematical models are helpful in refining and testing hypotheses, generating predictions, and exploring conditions outside the scope of experiments. Such models can be combined and calibrated with in vivo and in vitro data to identify critical cell-ECM interactions that drive developmental and homeostatic processes, or the progression of diseases. In this review, we focus on mathematical and computational models of the ECM in processes such as cell migration including cancer metastasis, and in tissue structure and morphogenesis. By highlighting the predictive power of these models, we aim to help bridge the gap between experimental and computational approaches to studying the ECM and to provide guidance on selecting an appropriate model framework to complement corresponding experimental studies.

Free volume theory explains the unusual behavior of viscosity in a non-confluent tissue during morphogenesis

A recent experiment on zebrafish blastoderm morphogenesis showed that the viscosity (η) of a non-confluent embryonic tissue grows sharply until a critical cell packing fraction (ϕS). The increase in η up to ϕS is similar to the behavior observed in several glass-forming materials, which suggests that the cell dynamics is sluggish or glass-like. Surprisingly, η is a constant above ϕS. To determine the mechanism of this unusual dependence of η on ϕ, we performed extensive simulations using an agent-based model of a dense non-confluent two-dimensional tissue. We show that polydispersity in the cell size, and the propensity of the cells to deform, results in the saturation of the available free area per cell beyond a critical packing fraction. Saturation in the free space not only explains the viscosity plateau above ϕS but also provides a relationship between equilibrium geometrical packing to the dramatic increase in the relaxation dynamics.

Enhanced mechanical heterogeneity of cell collectives due to temporal fluctuations in cell elasticity

Cells are dynamic systems characterized by temporal variations in biophysical properties such as stiffness and contractility. Recent studies show that the recruitment and release of actin filaments into and out of the cell cortex-a network of proteins underneath the cell membrane-leads to cell stiffening prior to division and softening immediately afterward. In three-dimensional (3D) cell collectives, it is unclear whether the stiffness change during division at the single-cell scale controls the spatial structure and dynamics at the multicellular scale. This is an important question to understand because cell stiffness variations impact cell spatial organization and cancer progression. Using a minimal 3D model incorporating cell birth, death, and cell-to-cell elastic and adhesive interactions, we investigate the effect of mechanical heterogeneity-variations in individual cell stiffnesses that make up the cell collective-on tumor spatial organization and cell dynamics. We discover that spatial mechanical heterogeneity characterized by a spheroid core composed of stiffer cells and softer cells in the periphery emerges within dense 3D cell collectives, which may be a general feature of multicellular tumor growth. We show that heightened spatial mechanical heterogeneity enhances single-cell dynamics and volumetric tumor growth driven by fluctuations in cell elasticity. Our results could have important implications in understanding how spatiotemporal variations in single-cell stiffness determine tumor growth and spread.

Radiomics-based signature of breast cancer on preoperative contrast-enhanced MRI to predict axillary metastasis

Aim: This study aimed to predict axillary metastasis using radiology features in dynamic contrast-enhanced MRI. Methods: This study included 243 breast lesions confirmed as malignant based on axillary status. Most outcome-predictive features were selected using four machine-learning algorithms. Receiver operating characteristic analysis was used to reflect diagnostic performance. Results: Least absolute shrinkage and selection operator was used to dimensionally reduce 1137 radiomics features to three features. Three optimal radiomics features were used to model construction. The logistic regression model achieved an accuracy of 97% and 85% in the training and test groups. Clinical utility was evaluated using decision curve analysis. Conclusion: The novel combination of radiomics analysis and machine-learning algorithm could predict axillary metastasis and prevent invasive manipulation.

Cell clusters softening triggers collective cell migration in vivo

Embryogenesis, tissue repair and cancer metastasis rely on collective cell migration. In vitro studies propose that cells are stiffer while migrating in stiff substrates, but softer when plated in compliant surfaces which are typically considered as non-permissive for migration. Here we show that cells within clusters from embryonic tissue dynamically decrease their stiffness in response to the temporal stiffening of their native substrate to initiate collective cell migration. Molecular and mechanical perturbations of embryonic tissues reveal that this unexpected mechanical response involves a mechanosensitive pathway relying on Piezo1-mediated microtubule deacetylation. We further show that decreasing microtubule acetylation and consequently cluster stiffness is sufficient to trigger collective cell migration in soft non-permissive substrates. This suggests that reaching an optimal cluster-to-substrate stiffness ratio is essential to trigger the onset of this collective process. Overall, these in vivo findings challenge the current understanding of collective cell migration and its physiological and pathological roles.

Optimal pathways control fixation of multiple mutations during cancer initiation

Cancer starts after initially healthy tissue cells accumulate several specific mutations or other genetic alterations. The dynamics of tumor formation is a very complex phenomenon due to multiple involved biochemical and biophysical processes. It leads to a very large number of possible pathways on the road to final fixation of all mutations that marks the beginning of the cancer, complicating the understanding of microscopic mechanisms of tumor formation. We present a new theoretical framework of analyzing the cancer initiation dynamics by exploring the properties of effective free-energy landscape of the process. It is argued that although there are many possible pathways for the fixation of all mutations in the system, there are only a few dominating pathways on the road to tumor formation. The theoretical approach is explicitly tested in the system with only two mutations using analytical calculations and Monte Carlo computer simulations. Excellent agreement with theoretical predictions is found for a large range of parameters, supporting our hypothesis and allowing us to better understand the mechanisms of cancer initiation. Our theoretical approach clarifies some important aspects of microscopic processes that lead to tumor formation.

Adhesion strength between cells regulate nonmonotonic growth by a biomechanical feedback mechanism

We determine how intercellular interactions and mechanical pressure experienced by single cells regulate cell proliferation using a minimal computational model for three-dimensional multicellular spheroid (MCS) growth. We discover that emergent spatial variations in the cell division rate, depending on the location of the cells either at the core or periphery within the MCS, is regulated by intercellular adhesion strength (fad). Varying fad results in nonmonotonic proliferation of cells in the MCS. A biomechanical feedback mechanism coupling the fad and microenvironment-dependent pressure fluctuations relative to a threshold value (pc) determines the onset of a dormant phase, and explains the nonmonotonic proliferation response. Increasing fad from low values enhances cell proliferation because pressure on individual cells is smaller compared with pc. However, at high fad, cells readily become dormant and cannot rearrange effectively in spacetime, leading to arrested cell proliferation. Utilizing our theoretical predictions, we explain experimental data on the impact of adhesion strength on cell proliferation and find good agreement. Our work, which shows that proliferation is regulated by pressure-adhesion feedback mechanism, may be a general feature of multicellular growth.

Mechanical feedback controls the emergence of dynamical memory in growing tissue monolayers

The growth of a tissue, which depends on cell–cell interactions and biologically relevant processes such as cell division and apoptosis, is regulated by a mechanical feedback mechanism. We account for these effects in a minimal two-dimensional model in order to investigate the consequences of mechanical feedback, which is controlled by a critical pressure, p c. A cell can only grow and divide if its pressure, due to interaction with its neighbors, is less than p c. Because temperature is not a relevant variable, the cell dynamics is driven by self-generated active forces (SGAFs) that arise due to cell division. We show that even in the absence of intercellular interactions, cells undergo diffusive behavior. The SGAF-driven diffusion is indistinguishable from the well-known dynamics of a free Brownian particle at a fixed finite temperature. When intercellular interactions are taken into account, we find persistent temporal correlations in the force–force autocorrelation function (FAF) that extends over a timescale of several cell division times. The time-dependence of the FAF reveals memory effects, which increases as p c increases. The observed non-Markovian effects emerge due to the interplay of cell division and mechanical feedback and are inherently a non-equilibrium phenomenon.

A mathematical model for phenotypic heterogeneity in breast cancer with implications for therapeutic strategies

Inevitably, almost all cancer patients develop resistance to targeted therapy. Intratumour heterogeneity is a major cause of drug resistance. Mathematical models that explain experiments quantitatively are useful in understanding the origin of intratumour heterogeneity, which then could be used to explore scenarios for efficacious therapy. Here, we develop a mathematical model to investigate intratumour heterogeneity in breast cancer by exploiting the observation that HER2+ and HER2- cells could divide symmetrically or asymmetrically. Our predictions for the evolution of cell fractions are in quantitative agreement with single-cell experiments. Remarkably, the colony size of HER2+ cells emerging from a single HER2- cell (or vice versa), which occurs in about four cell doublings, also agrees with experimental results, without tweaking any parameter in the model. The theory explains experimental data on the responses of breast tumours under different treatment protocols. We then used the model to predict that, not only the order of two drugs, but also the treatment period for each drug and the tumour cell plasticity could be manipulated to improve the treatment efficacy. Mathematical models, when integrated with data on patients, make possible exploration of a broad range of parameters readily, which might provide insights in devising effective therapies.

Interparticle Adhesion Regulates the Surface Roughness of Growing Dense Three-Dimensional Active Particle Aggregates

Activity and self-generated motion are fundamental features observed in many living and nonliving systems. Given that interparticle adhesive forces can regulate particle dynamics, we investigate how interparticle adhesion strength controls the boundary growth and roughness of active particle aggregates. Using particle based simulations incorporating both activity (birth, death, and growth) and systematic physical interactions (elasticity and adhesion), we establish that interparticle adhesion strength (f^ad) controls the surface roughness of a densely packed three-dimensional(3D) active particle aggregate expanding into a highly viscous medium. We discover that the surface roughness of a 3D active particle aggregate increases in proportion to the interparticle adhesion strength (f^ad) and show that asymmetry in the radial and transverse active particle mean-squared displacement (MSD) suppresses 3D surface roughness at lower adhesion strengths. By analyzing the statistical properties of particle displacements at the aggregate periphery, we determine that the 3D surface roughness is driven by the movement of active particle toward the core at high interparticle adhesion strengths. Our results elucidate the physics controlling the expansion of adhesive 3D active particle collectives into a highly viscous medium, with implications into understanding stochastic interface growth in active matter systems characterized by self-generation of particles.

Mechanical heterogeneity along single cell-cell junctions is driven by lateral clustering of cadherins during vertebrate axis elongation

Morphogenesis is governed by the interplay of molecular signals and mechanical forces across multiple length scales. The last decade has seen tremendous advances in our understanding of the dynamics of protein localization and turnover at subcellular length scales, and at the other end of the spectrum, of mechanics at tissue-level length scales. Integrating the two remains a challenge, however, because we lack a detailed understanding of the subcellular patterns of mechanical properties of cells within tissues. Here, in the context of the elongating body axis of Xenopus embryos, we combine tools from cell biology and physics to demonstrate that individual cell-cell junctions display finely-patterned local mechanical heterogeneity along their length. We show that such local mechanical patterning is essential for the cell movements of convergent extension and is imparted by locally patterned clustering of a classical cadherin. Finally, the patterning of cadherins and thus local mechanics along cell-cell junctions are controlled by Planar Cell Polarity signaling, a key genetic module for CE that is mutated in diverse human birth defects.

Self-generated persistent random forces drive phase separation in growing tumors

A single solid tumor, composed of nearly identical cells, exhibits heterogeneous dynamics. Dynamics of cells in the core is glass-like, whereas those in the periphery undergoes diffusive or super-diffusive behavior. Quantification of heterogeneity using the mean square displacement or the self-intermediate scattering function, which involves averaging over the cell population, hides the complexity of the collective movement. Using the t-distributed stochastic neighbor embedding (t-SNE), a popular unsupervised machine learning dimensionality reduction technique, we show that the phase space structure of an evolving colony of cells, driven by cell division and apoptosis, partitions into nearly disjoint sets composed principally of the core and periphery cells. The non-equilibrium phase separation is driven by the differences in the persistence of self-generated active forces induced by cell division. Extensive heterogeneity revealed by t-SNE paves the way toward understanding the origins of intratumor heterogeneity using experimental imaging data.