Mechanical control of cell flow in multicellular spheroids

Proliferation-driven mechanical compression induces signalling centre formation during mammalian organ development

Localized sources of morphogens, called signalling centres, play a fundamental role in coordinating tissue growth and cell fate specification during organogenesis. However, how these signalling centres are established in tissues during embryonic development is still unclear. Here we show that the main signalling centre orchestrating development of rodent incisors, the enamel knot (EK), is specified by a cell proliferation-driven buildup in compressive stresses (mechanical pressure) in the tissue. Direct mechanical measurements indicate that the stresses generated by cell proliferation are resisted by the surrounding tissue, creating a circular pattern of mechanical anisotropy with a region of high compressive stress at its centre that becomes the EK. Pharmacological inhibition of proliferation reduces stresses and suppresses EK formation, and application of external pressure in proliferation-inhibited conditions rescues the formation of the EK. Mechanical information is relayed intracellularly through YAP protein localization, which is cytoplasmic in the region of compressive stress that establishes the EK and nuclear in the stretched anisotropic cells that resist the pressure buildup around the EK. Together, our data identify a new role for proliferation-driven mechanical compression in the specification of a model signalling centre during mammalian organ development.

Morphology as indicator of adaptive changes of model tissues in osmotically and chemically changing environments

We investigate the formation and maintenance of the homeostatic state in the case of 2D epithelial tissues following an induction of hyperosmotic conditions, using media enriched with 80 to 320 mOsm of mannitol, NaCl, and urea. We characterise the changes in the tissue immediately after the osmotic shock, and follow it until the new homeostatic state is formed. We characterise changes in cooperative motility and proliferation pressure in the tissue upon treatment with the help of a theoretical model based on the delayed Fisher-Kolmogorov formalism, where the delay in density evolution is induced by the the finite time of the cell division. Finally we explore the adaptation of the homeostatic tissue to highly elevated osmotic conditions by evaluating the morphology and topology of cells after 20 days in incubation. We find that hyperosmotic environments together with changes in the extracellular matrix induce different mechanical states in viable tissues, where only some remain functional. The perspective is a relation between tissue topology and function, which could be explored beyond the scope of this manuscript. Experimental investigation of morphological effect of change of osmotic conditions on long-term tissue morphology and topology Effect of osmotic changes on transient tissue growth behaviour Analysis of recovery process of tissues post-osmotic-shock Toxicity limits of osmolytes in mid- to long-term tissue evolution Tissue adaptation to physiological changes in environment Long-term tissue stabilisation under altered osmotic conditions.

Curvature induces active velocity waves in rotating spherical tissues

The multicellular organization of diverse systems, including embryos, intestines, and tumors relies on coordinated cell migration in curved environments. In these settings, cells establish supracellular patterns of motion, including collective rotation and invasion. While such collective modes have been studied extensively in flat systems, the consequences of geometrical and topological constraints on collective migration in curved systems are largely unknown. Here, we discover a collective mode of cell migration in rotating spherical tissues manifesting as a propagating single-wavelength velocity wave. This wave is accompanied by an apparently incompressible supracellular flow pattern featuring topological defects as dictated by the spherical topology. Using a minimal active particle model, we reveal that this collective mode arises from the effect of curvature on the active flocking behavior of a cell layer confined to a spherical surface. Our results thus identify curvature-induced velocity waves as a mode of collective cell migration, impacting the dynamical organization of 3D curved tissues.

Activity-driven tissue alignment in proliferating spheroids

We extend the continuum theory of active nematic fluids to study cell flows and tissue dynamics inside multicellular spheroids, spherical, self-assembled aggregates of cells that are widely used as model systems to study tumour dynamics. Cells near the surface of spheroids have better access to nutrients and therefore proliferate more rapidly than those in the resource-depleted core. Using both analytical arguments and three-dimensional simulations, we find that the proliferation gradients result in flows and in gradients of activity both of which can align the orientation axis of cells inside the aggregates. Depending on environmental conditions and the intrinsic tissue properties, we identify three distinct alignment regimes: spheroids in which all the cells align either radially or tangentially to the surface throughout the aggregate and spheroids with angular cell orientation close to the surface and radial alignment in the core. The continuum description of tissue dynamics inside spheroids not only allows us to infer dynamic cell parameters from experimentally measured cell alignment profiles, but more generally motivates novel mechanisms for controlling the alignment of cells within aggregates which has been shown to influence the mechanical properties and invasive capabilities of tumors.

A physical model for dynamic assembly of human salivary stem/progenitor microstructures

The in vitro reconstructions of human salivary glands in service of their eventual medical use represent a challenge for tissue engineering. Here, we present a theoretical approach to the dynamical formation of acinar structures from human salivary cells, focusing on observed stick-slip radial expansion as well as possible growth instabilities. Our findings demonstrate the critical importance of basement membrane remodeling in controlling the growth process.

Activity gradients in two- and three-dimensional active nematics

We numerically investigate how spatial variations of extensile or contractile active stress affect bulk active nematic systems in two and three dimensions. In the absence of defects, activity gradients drive flows which re-orient the nematic director field and thus act as an effective anchoring force. At high activity, defects are created and the system transitions into active turbulence, a chaotic flow state characterized by strong vorticity. We find that in two-dimensional (2D) systems active torques robustly align +1/2 defects parallel to activity gradients, with defect heads pointing towards contractile regions. In three-dimensional (3D) active nematics disclination lines preferentially lie in the plane perpendicular to activity gradients due to active torques acting on line segments. The average orientation of the defect structures in the plane perpendicular to the line tangent depends on the defect type, where wedge-like +1/2 defects align parallel to activity gradients, while twist defects are aligned anti-parallel. Understanding the response of active nematic fluids to activity gradients is an important step towards applying physical theories to biology, where spatial variations of active stress impact morphogenetic processes in developing embryos and affect flows and deformations in growing cell aggregates, such as tumours.

Spontaneous flow created by active topological defects

Topological defects are at the root of the large-scale organization of liquid crystals. In two-dimensional active nematics, two classes of topological defects of charges [Formula: see text] are known to play a major role due to active stresses. Despite this importance, few analytical results have been obtained on the flow-field and active-stress patterns around active topological defects. Using the generic hydrodynamic theory of active systems, we investigate the flow and stress patterns around these topological defects in unbounded, two-dimensional active nematics. Under generic assumptions, we derive analytically the spontaneous velocity and stall force of self-advected defects in the presence of both shear and rotational viscosities. Applying our formalism to the dynamics of monolayers of elongated cells at confluence, we show that the non-conservation of cell number generically increases the self-advection velocity and could provide an explanation for their observed role in cellular extrusion and multilayering. We finally investigate numerically the influence of the Ericksen stress. Our work paves the way to a generic study of the role of topological defects in active nematics, and in particular in monolayers of elongated cells.

Viscoelasticity, Like Forces, Plays a Role in Mechanotransduction

Viscoelasticity and its alteration in time and space has turned out to act as a key element in fundamental biological processes in living systems, such as morphogenesis and motility. Based on experimental and theoretical findings it can be proposed that viscoelasticity of cells, spheroids and tissues seems to be a collective characteristic that demands macromolecular, intracellular component and intercellular interactions. A major challenge is to couple the alterations in the macroscopic structural or material characteristics of cells, spheroids and tissues, such as cell and tissue phase transitions, to the microscopic interferences of their elements. Therefore, the biophysical technologies need to be improved, advanced and connected to classical biological assays. In this review, the viscoelastic nature of cytoskeletal, extracellular and cellular networks is presented and discussed. Viscoelasticity is conceptualized as a major contributor to cell migration and invasion and it is discussed whether it can serve as a biomarker for the cells' migratory capacity in several biological contexts. It can be hypothesized that the statistical mechanics of intra- and extracellular networks may be applied in the future as a powerful tool to explore quantitatively the biomechanical foundation of viscoelasticity over a broad range of time and length scales. Finally, the importance of the cellular viscoelasticity is illustrated in identifying and characterizing multiple disorders, such as cancer, tissue injuries, acute or chronic inflammations or fibrotic diseases.

Viscoelastic Properties in Cancer: From Cells to Spheroids

AFM-based rheology methods enable the investigation of the viscoelastic properties of cancer cells. Such properties are known to be essential for cell functions, especially for malignant cells. Here, the relevance of the force modulation method was investigated to characterize the viscoelasticity of bladder cancer cells of various invasiveness on soft substrates, revealing that the rheology parameters are a signature of malignancy. Furthermore, the collagen microenvironment affects the viscoelastic moduli of cancer cell spheroids; thus, collagen serves as a powerful proxy, leading to an increase of the dynamic moduli vs. frequency, as predicted by a double power law model. Taken together, these results shed new light on how cancer cells and tissues adapt their viscoelastic properties depending on their malignancy and the microenvironment. This method could be an attractive way to control their properties in the future, based on the similarity of spheroids with in vivo tumor models.

Hydraulic and electric control of cell spheroids

We use a theoretical approach to examine the effect of a radial fluid flow or electric current on the growth and homeostasis of a cell spheroid. Such conditions may be generated by a drain of micrometric diameter. To perform this analysis, we describe the tissue as a continuum. We include active mechanical, electric, and hydraulic components in the tissue material properties. We consider a spherical geometry and study the effect of the drain on the dynamics of the cell aggregate. We show that a steady fluid flow or electric current imposed by the drain could be able to significantly change the spheroid long-time state. In particular, our work suggests that a growing spheroid can systematically be driven to a shrinking state if an appropriate external field is applied. Order-of-magnitude estimates suggest that such fields are of the order of the indigenous ones. Similarities and differences with the case of tumors and embryo development are briefly discussed.

Extracellular matrix in multicellular aggregates acts as a pressure sensor controlling cell proliferation and motility

Imposed deformations play an important role in morphogenesis and tissue homeostasis, both in normal and pathological conditions. To perceive mechanical perturbations of different types and magnitudes, tissues need appropriate detectors, with a compliance that matches the perturbation amplitude. By comparing results of selective osmotic compressions of CT26 mouse cells within multicellular aggregates and global aggregate compressions, we show that global compressions have a strong impact on the aggregates growth and internal cell motility, while selective compressions of same magnitude have almost no effect. Both compressions alter the volume of individual cells in the same way over a shor-timescale, but, by draining the water out of the extracellular matrix, the global one imposes a residual compressive mechanical stress on the cells over a long-timescale, while the selective one does not. We conclude that the extracellular matrix is as a sensor that mechanically regulates cell proliferation and migration in a 3D environment.

Quantifying Material Properties of Cell Monolayers by Analyzing Integer Topological Defects

In developing organisms, internal cellular processes generate mechanical stresses at the tissue scale. The resulting deformations depend on the material properties of the tissue, which can exhibit long-ranged orientational order and topological defects. It remains a challenge to determine these properties on the time scales relevant for developmental processes. Here, we build on the physics of liquid crystals to determine material parameters of cell monolayers. Specifically, we use a hydrodynamic description to characterize the stationary states of compressible active polar fluids around defects. We illustrate our approach by analyzing monolayers of C2C12 cells in small circular confinements, where they form a single topological defect with integer charge. We find that such monolayers exert compressive stresses at the defect centers, where localized cell differentiation and formation of three-dimensional shapes is observed.

Mechanical Control of Cell Proliferation Increases Resistance to Chemotherapeutic Agents

While many cellular mechanisms leading to chemotherapeutic resistance have been identified, there is an increasing realization that tumor-stroma interactions also play an important role. In particular, mechanical alterations are inherent to solid cancer progression and profoundly impact cell physiology. Here, we explore the influence of compressive stress on the efficacy of chemotherapeutics in pancreatic cancer spheroids. We find that increased compressive stress leads to decreased drug efficacy. Theoretical modeling and experiments suggest that mechanical stress decreases cell proliferation which in turn reduces the efficacy of chemotherapeutics that target proliferating cells. Our work highlights a mechanical form of drug resistance and suggests new strategies for therapy.

Mathematical modelling reveals cellular dynamics within tumour spheroids

Tumour spheroids are widely used as an in vitro assay for characterising the dynamics and response to treatment of different cancer cell lines. Their popularity is largely due to the reproducible manner in which spheroids grow: the diffusion of nutrients and oxygen from the surrounding culture medium, and their consumption by tumour cells, causes proliferation to be localised at the spheroid boundary. As the spheroid grows, cells at the spheroid centre may become hypoxic and die, forming a necrotic core. The pressure created by the localisation of tumour cell proliferation and death generates an cellular flow of tumour cells from the spheroid rim towards its core. Experiments by Dorie et al. showed that this flow causes inert microspheres to infiltrate into tumour spheroids via advection from the spheroid surface, by adding microbeads to the surface of tumour spheroids and observing the distribution over time. We use an off-lattice hybrid agent-based model to re-assess these experiments and establish the extent to which the spatio-temporal data generated by microspheres can be used to infer kinetic parameters associated with the tumour spheroids that they infiltrate. Variation in these parameters, such as the rate of tumour cell proliferation or sensitivity to hypoxia, can produce spheroids with similar bulk growth dynamics but differing internal compositions (the proportion of the tumour which is proliferating, hypoxic/quiescent and necrotic/nutrient-deficient). We use this model to show that the types of experiment conducted by Dorie et al. could be used to infer spheroid composition and parameters associated with tumour cell lines such as their sensitivity to hypoxia or average rate of proliferation, and note that these observations cannot be conducted within previous continuum models of microbead infiltration into tumour spheroids as they rely on resolving the trajectories of individual microbeads.

Tumour spheroids are an experimental assay used to characterise the dynamics and response to treatment of different cancer cell lines. Previous experiments have demonstrated that the localisation of tumour cell proliferation to the spheroid edge (due to the gradient formed by nutrient diffusing from the surrounding medium) causes cells to be pushed from the proliferative rim towards the nutrient-deficient necrotic core. This movement allows inert particles to infiltrate tumour spheroids. We use a hybrid agent-based model to reproduce this data. We show further how data from individual microbead trajectories can be used to infer the composition of simulated tumour spheroids, and to estimate model parameters pertaining to tumour cell proliferation rates and their responses to hypoxia. Since these measurements are possible using modern imaging techniques, this could motivate new experiments in which spheroid composition could be inferred by observing passive infiltration of inert particles.

A minimal modeling framework of radiation and immune system synergy to assist radiotherapy planning

Recent evidence indicates the ability of radiotherapy to induce local and systemic tumor-specific immune responses as a result of immunogenic cell death. However, fractionation regimes routinely used in clinical practice typically ignore the synergy between radiation and the immune system, and instead attempt to completely eradicate tumors by the direct lethal effect of radiation on cancer cells. This paradigm is expected to change in the near future due to the potential benefits of considering radiation-induced antitumor immunity during treatment planning. Towards this goal, we propose a minimal modeling framework based on key aspects of the tumor-immune system interplay to simulate the effects of radiation on tumors and the immunological consequences of radiotherapy. The impacts of tumor-associated vasculature and intratumoral oxygen-mediated heterogeneity on treatment outcomes are ininvestigated. The model provides estimates of the minimum radiation doses required for tumor eradication given a certain number of treatment fractions. Moreover, estimates of treatment duration for disease control given predetermined fractional radiation doses can be also obtained. Although theoretical in nature, this study motivates the development and establishment of immune-based decision-support tools in radiotherapy planning.

Fluid pumping and active flexoelectricity can promote lumen nucleation in cell assemblies

We discuss the physical mechanisms that promote or suppress the nucleation of a fluid-filled lumen inside a cell assembly or a tissue. We discuss lumen formation in a continuum theory of tissue material properties in which the tissue is described as a 2-fluid system to account for its permeation by the interstitial fluid, and we include fluid pumping as well as active electric effects. Considering a spherical geometry and a polarized tissue, our work shows that fluid pumping and tissue flexoelectricity play a crucial role in lumen formation. We furthermore explore the large variety of long-time states that are accessible for the cell aggregate and its lumen. Our work reveals a role of the coupling of mechanical, electrical, and hydraulic phenomena in tissue lumen formation.

Pattern formation of skin cancers: Effects of cancer proliferation and hydrodynamic interactions

We study pattern formation of skin cancers by means of numerical simulation of a binary system consisting of cancer and healthy cells. We extend the conventional model H for macrophase separations by considering a logistic growth of cancer cells and also a mechanical friction between dermis and epidermis. Importantly, our model exhibits a microphase separation due to the proliferation of cancer cells. By numerically solving the time evolution equations of the cancer composition and its velocity, we show that the phase separation kinetics strongly depends on the cell proliferation rate as well as on the strength of hydrodynamic interactions. A steady-state diagram of cancer patterns is established in terms of these two dynamical parameters and some of the patterns correspond to clinically observed cancer patterns. Furthermore, we examine in detail the time evolution of the average composition of cancer cells and the characteristic length of the microstructures. Our results demonstrate that different sequence of cancer patterns can be obtained by changing the proliferation rate and/or hydrodynamic interactions.

Quantitative cell-based model predicts mechanical stress response of growing tumor spheroids over various growth conditions and cell lines

Model simulations indicate that the response of growing cell populations on mechanical stress follows the same functional relationship and is predictable over different cell lines and growth conditions despite experimental response curves look largely different. We develop a hybrid model strategy in which cells are represented by coarse-grained individual units calibrated with a high resolution cell model and parameterized by measurable biophysical and cell-biological parameters. Cell cycle progression in our model is controlled by volumetric strain, the latter being derived from a bio-mechanical relation between applied pressure and cell compressibility. After parameter calibration from experiments with mouse colon carcinoma cells growing against the resistance of an elastic alginate capsule, the model adequately predicts the growth curve in i) soft and rigid capsules, ii) in different experimental conditions where the mechanical stress is generated by osmosis via a high molecular weight dextran solution, and iii) for other cell types with different growth kinetics from the growth kinetics in absence of external stress. Our model simulation results suggest a generic, even quantitatively same, growth response of cell populations upon externally applied mechanical stress, as it can be quantitatively predicted using the same growth progression function.

Colonic stent-induced mechanical compression may suppress cancer cell proliferation in malignant large bowel obstruction

BACKGROUND

The short-term safety and efficacy of insertion of a self-expandable metallic colonic stent (SEMS) followed by elective surgery, "bridge to surgery (BTS)", for malignant large bowel obstruction (MLBO) have been well described; however, the influence on long-term oncological outcomes is unclear. The aim of this study was to evaluate changes in oncological characteristics in colorectal cancer (CRC) tissues after SEMS insertion, focusing on growth factors, cell cycle and apoptosis.

METHODS

From January 2013 to September 2014, a total of 25 patients with MLBO who underwent BTS at our single institution were retrospectively included. Paired CRC tissue samples before (endoscopic biopsy) and after SEMS insertion (surgically resected) were collected from each patient. EGFR, VEGF, Ki-67, p27^kip1 and TUNEL expression were determined by immunohistochemistry.

RESULTS

No clinical or subclinical perforations evaluated by mechanical ulceration pathologically were observed. Epithelial exfoliation, tumour necrosis, infiltration of inflammatory cells and fibrosis were observed in SEMS-inserted surgically-resected specimens. Overall, 84% (21/25) and 60% (15/25) of patients exhibited no change or a decrease in staining category, respectively, for EGFR and VEGF expression after SEMS insertion. A significant decrease in Ki-67 expression was observed in surgically-resected specimens compared with endoscopic biopsy specimens (P < 0.01). The upstream cell cycle inhibitor, p27^kip1, was significantly increased after SEMS insertion (P = 0.049).

CONCLUSIONS

Although the long-term safety of BTS should be determined in a future clinical trial, mechanical compression by SEMS may suppress cancer cell proliferation and this result could provide some insights into the issue.

Permeability and viscoelastic fracture of a model tumor under interstitial flow

Interstitial flow in tumors is a key mechanism leading to cancer metastasis. Tumor growth is accompanied by the development of a leaky vasculature, which increases intratumoral pressure and generates an outward interstitial flow. This flow promotes tumor cell migration away from the tumor. The nature of such interstitial flow depends on the coupling between hydrodynamic conditions and material properties of the tumor, such as porosity and deformability. Here we investigate this coupling by means of a microfluidic model of interstitial flow through a tumor, which is represented by a tumor cell aggregate. For a weak intratumoral pressure, the model tumor behaves as a viscoelastic material of low permeability, which we estimate by means of a newly developed microfluidic device. As intratumoral pressure is raised, the model tumor deforms and its permeability increases. For a high enough pressure, localized intratumoral fracture occurs, which creates preferential flow paths and causes tumor cell detachment. The energy required to fracture depends on the rate of variation of intratumoral pressure, as explained here by a theoretical model originally derived to describe polymer adhesion. Besides the well-established picture of individual tumor cells migrating away under interstitial flow, our findings suggest that intratumoral pressures observed in tumors can suffice to detach tumor fragments, which may thus be an important mechanism to release cancer cells and initiate metastasis.

Predicting the role of microstructural and biomechanical cues in tumor growth and spreading

A continuum porous media model is developed for elucidating the role of the mechanical cues in regulating tumor growth and spreading. It is shown that stiffer matrices and higher cell-matrix adhesion limit tumor growth and spreading toward the surrounding tissue. Higher matrix porosities, conversely, favor the growth and the dissemination of tumor cells. This model could be used for predicting the response of malignant masses to novel therapeutic agents affecting directly the tumor microenvironment and its micromechanical cues.

Multicellular Tumor Spheroids (MCTS) as a 3D In Vitro Evaluation Tool of Nanoparticles

Multicellular tumor spheroid models (MCTS) are often coined as 3D in vitro models that can mimic the microenvironment of tissues. MCTS have gained increasing interest in the nano-biotechnology field as they can provide easily accessible information on the performance of nanoparticles without using animal models. Considering that many countries have put restrictions on animals testing, which will only tighten in the future as seen by the recent developments in the Netherlands, 3D models will become an even more valuable tool. Here, an overview on MCTS is provided, focusing on their use in cancer research as most nanoparticles are tested in MCTS for treatment of primary tumors. Thereafter, various types of nanoparticles-from self-assembled block copolymers to inorganic nanoparticles, are discussed. A range of physicochemical parameters including the size, shape, surface chemistry, ligands attachment, stability, and stiffness are found to influence nanoparticles in MCTS. Some of these studies are complemented by animal studies confirming that lessons from MCTS can in part predict the behaviour in vivo. In summary, MCTS are suitable models to gain additional information on nanoparticles. While not being able to replace in vivo studies, they can bridge the gap between traditional 2D in vitro studies and in vivo models.

Immotile Active Matter: Activity from Death and Reproduction

Unlike equilibrium atomic solids, biofilms-soft solids composed of bacterial cells-do not experience significant thermal fluctuations at the constituent level. However, living cells stochastically reproduce and die, provoking a mechanical response. We investigate the mechanical consequences of cellular death and reproduction by measuring surface-height fluctuations of biofilms containing two mutually antagonistic strains of Vibrio cholerae that kill one another on contact via the type VI secretion system. While studies of active matter typically focus on activity via constituent mobility, here, activity is mediated by reproduction and death events in otherwise immobilized cells. Biofilm surface topography is measured in the nearly homeostatic limit via white light interferometry. Although biofilms are far from equilibrium systems, measured surface-height fluctuation spectra resemble the spectra of thermal permeable membranes but with an activity-mediated effective temperature, as predicted by Risler, Peilloux, and Prost [Phys. Rev. Lett. 115, 258104 (2015)PRLTAO0031-900710.1103/PhysRevLett.115.258104]. By comparing the activity of killer strains of V. cholerae with that of genetically modified strains that cannot kill each other and validating with individual-based simulations, we demonstrate that extracted effective temperatures increase with the amount of death and reproduction and that death and reproduction can fluidize biofilms. Together, these observations demonstrate the unique physical consequences of activity mediated by death and reproduction events.

Solid tumors are poroelastic solids with a chemo-mechanical feedback on growth

The experimental evidence that a feedback exists between growth and stress in tumors poses challenging questions. First, the rheological properties (the "constitutive equations") of aggregates of malignant cells are still a matter of debate. Secondly, the feedback law (the "growth law") that relates stress and mitotic-apoptotic rate is far to be identified. We address these questions on the basis of a theoretical analysis of in vitro and in vivo experiments that involve the growth of tumor spheroids. We show that solid tumors exhibit several mechanical features of a poroelastic material, where the cellular component behaves like an elastic solid. When the solid component of the spheroid is loaded at the boundary, the cellular aggregate grows up to an asymptotic volume that depends on the exerted compression. Residual stress shows up when solid tumors are radially cut, highlighting a peculiar tensional pattern. By a novel numerical approach we correlate the measured opening angle and the underlying residual stress in a sphere. The features of the mechanobiological system can be explained in terms of a feedback of mechanics on the cell proliferation rate as modulated by the availability of nutrient, that is radially damped by the balance between diffusion and consumption. The volumetric growth profiles and the pattern of residual stress can be theoretically reproduced assuming a dependence of the target stress on the concentration of nutrient which is specific of the malignant tissue.

Spatiotemporal variation of endogenous cell-generated stresses within 3D multicellular spheroids

Multicellular spheroids serve as an excellent platform to study tissue behavior and tumor growth in a controlled, three-dimensional (3D) environment. While molecular and cellular studies have long used this platform to study cell behavior in 3D, only recently have studies using multicellular spheroids shown an important role for the mechanics of the microenvironment in a wide range of cellular processes, including during tumor progression. Despite the well-established relevance of mechanical cues to cell behavior and the numerous studies on mechanics using 2D cell culture systems, the spatial and temporal variations in endogenous cellular forces within growing multicellular aggregates remain unknown. Using cell-sized oil droplets with controlled physicochemical properties as force transducers in mesenchymal cell aggregates, we show that the magnitude of cell-generated stresses varies only weakly with spatial location within the spherical aggregate, but it increases considerably over time during aggregate compaction and growth. Moreover, our results indicate that the temporal increase in cellular stresses is due to increasing cell pulling forces transmitted via integrin-mediated cell adhesion, consistent with the need for larger intercellular pulling forces to compact cell aggregates.

Topological defects in epithelia govern cell death and extrusion

Epithelial tissues (epithelia) remove excess cells through extrusion, preventing the accumulation of unnecessary or pathological cells. The extrusion process can be triggered by apoptotic signalling, oncogenic transformation and overcrowding of cells. Despite the important linkage of cell extrusion to developmental, homeostatic and pathological processes such as cancer metastasis, its underlying mechanism and connections to the intrinsic mechanics of the epithelium are largely unexplored. We approach this problem by modelling the epithelium as an active nematic liquid crystal (that has a long range directional order), and comparing numerical simulations to strain rate and stress measurements within monolayers of MDCK (Madin Darby canine kidney) cells. Here we show that apoptotic cell extrusion is provoked by singularities in cell alignments in the form of comet-shaped topological defects. We find a universal correlation between extrusion sites and positions of nematic defects in the cell orientation field in different epithelium types. The results confirm the active nematic nature of epithelia, and demonstrate that defect-induced isotropic stresses are the primary precursors of mechanotransductive responses in cells, including YAP (Yes-associated protein) transcription factor activity, caspase-3-mediated cell death, and extrusions. Importantly, the defect-driven extrusion mechanism depends on intercellular junctions, because the weakening of cell-cell interactions in an α-catenin knockdown monolayer reduces the defect size and increases both the number of defects and extrusion rates, as is also predicted by our model. We further demonstrate the ability to control extrusion hotspots by geometrically inducing defects through microcontact printing of patterned monolayers. On the basis of these results, we propose a mechanism for apoptotic cell extrusion: spontaneously formed topological defects in epithelia govern cell fate. This will be important in predicting extrusion hotspots and dynamics in vivo, with potential applications to tissue regeneration and the suppression of metastasis. Moreover, we anticipate that the analogy between the epithelium and active nematic liquid crystals will trigger further investigations of the link between cellular processes and the material properties of epithelia.

Intracellular dynamics measurements with full field optical coherence tomography suggest hindering effect of actomyosin contractility on organelle transport

Intracellular motion can be quantitatively monitored in tissues using coherence-gated microscopic techniques. With full-field optical coherence tomography (FFOCT), the use of high numerical aperture microscope objectives provides a high resolution mapping of intracellular dynamics that are probed with subwavelength sensitivity. In the upper temporal bandwidth that we have used (1-6 Hz) the main contribution to the dynamic signal arises from the overall dynamical, optically heterogeneous cytoplasm. We propose a method to specifically study the impact of actomyosin contractility on the intracellular dynamic signal by performing high throughput, comparative measurements of multicellular aggregates with and without blebbistatin action, a selective inhibitor of class-II myosins that disrupts actomyosin contractile activity. Our results indicate a significant increase in the fraction of the signal that decorrelates within 1 second after inhibition of contractility. This observation mitigates the anticipated importance of actomyosin contractile forces to directly move organelles, but highlights their role in hindering organelle transport via their stiffening effect of the viscoelastic cytoplasm.

On the morphological stability of multicellular tumour spheroids growing in porous media

Multicellular tumour spheroids (MCTSs) are extensively used as in vitro system models for investigating the avascular growth phase of solid tumours. In this work, we propose a continuous growth model of heterogeneous MCTSs within a porous material, taking into account a diffusing nutrient from the surrounding material directing both the proliferation rate and the mobility of tumour cells. At the time scale of interest, the MCTS behaves as an incompressible viscous fluid expanding inside a porous medium. The cell motion and proliferation rate are modelled using a non-convective chemotactic mass flux, driving the cell expansion in the direction of the external nutrients' source. At the early stages, the growth dynamics is derived by solving the quasi-stationary problem, obtaining an initial exponential growth followed by an almost linear regime, in accordance with experimental observations. We also perform a linear-stability analysis of the quasi-static solution in order to investigate the morphological stability of the radially symmetric growth pattern. We show that mechano-biological cues, as well as geometric effects related to the size of the MCTS subdomains with respect to the diffusion length of the nutrient, can drive a morphological transition to fingered structures, thus triggering the formation of complex shapes that might promote tumour invasiveness. The results also point out the formation of a retrograde flow in the MCTS close to the regions where protrusions form, that could describe the initial dynamics of metastasis detachment from the in vivo tumour mass. In conclusion, the results of the proposed model demonstrate that the integration of mathematical tools in biological research could be crucial for better understanding the tumour's ability to invade its host environment.

Tumor proliferation and diffusion on percolation clusters

We study in silico the influence of host tissue inhomogeneity on tumor cell proliferation and diffusion by simulating the mobility of a tumor on percolation clusters with different homogeneities of surrounding tissues. The proliferation and diffusion of a tumor in an inhomogeneous tissue could be characterized in the framework of the percolation theory, which displays similar thresholds (0.54, 0.44, and 0.37, respectively) for tumor proliferation and diffusion in three kinds of lattices with 4, 6, and 8 connecting near neighbors. Our study reveals the existence of a critical transition concerning the survival and diffusion of tumor cells with leaping metastatic diffusion movement in the host tissues. Tumor cells usually flow in the direction of greater pressure variation during their diffusing and infiltrating to a further location in the host tissue. Some specific sites suitable for tumor invasion were observed on the percolation cluster and around these specific sites a tumor can develop into scattered tumors linked by some advantage tunnels that facilitate tumor invasion. We also investigate the manner that tissue inhomogeneity surrounding a tumor may influence the velocity of tumor diffusion and invasion. Our simulation suggested that invasion of a tumor is controlled by the homogeneity of the tumor microenvironment, which is basically consistent with the experimental report by Riching et al. as well as our clinical observation of medical imaging. Both simulation and clinical observation proved that tumor diffusion and invasion into the surrounding host tissue is positively correlated with the homogeneity of the tissue.

Cancer therapeutic potential of combinatorial immuno- and vasomodulatory interventions

Currently, most of the basic mechanisms governing tumour-immune system interactions, in combination with modulations of tumour-associated vasculature, are far from being completely understood. Here, we propose a mathematical model of vascularized tumour growth, where the main novelty is the modelling of the interplay between functional tumour vasculature and effector cell recruitment dynamics. Parameters are calibrated on the basis of different in vivo immunocompromised Rag1(-/-) and wild-type (WT) BALB/c murine tumour growth experiments. The model analysis supports that tumour vasculature normalization can be a plausible and effective strategy to treat cancer when combined with appropriate immunostimulations. We find that improved levels of functional tumour vasculature, potentially mediated by normalization or stress alleviation strategies, can provide beneficial outcomes in terms of tumour burden reduction and growth control. Normalization of tumour blood vessels opens a therapeutic window of opportunity to augment the antitumour immune responses, as well as to reduce intratumoral immunosuppression and induced hypoxia due to vascular abnormalities. The potential success of normalizing tumour-associated vasculature closely depends on the effector cell recruitment dynamics and tumour sizes. Furthermore, an arbitrary increase in the initial effector cell concentration does not necessarily imply better tumour control. We evidence the existence of an optimal concentration range of effector cells for tumour shrinkage. Based on these findings, we suggest a theory-driven therapeutic proposal that optimally combines immuno- and vasomodulatory interventions.

Defect-Mediated Morphologies in Growing Cell Colonies

Morphological trends in growing colonies of living cells are at the core of physiological and evolutionary processes. Using active gel equations, which include cell division, we show that shape changes during the growth can be regulated by the dynamics of topological defects in the orientation of cells. The friction between the dividing cells and underlying substrate drives anisotropic colony shapes toward more isotropic morphologies, by mediating the number density and velocity of topological defects. We show that the defects interact with the interface at a specific interaction range, set by the vorticity length scale of flows within the colony, and that the cells predominantly reorient parallel to the interface due to division-induced active stresses.

A question of time: tissue adaptation to mechanical forces

While much attention has been focused on the force-generating mechanisms responsible for shaping developing embryos, less is known about the ways in which cells in animal tissues respond to mechanical stimuli. Forces will arise within a tissue as the result of processes such as local cell death, growth and division, but they can also be an indirect consequence of morphogenetic movements in neighbouring tissues or be imposed from the outside, for example, by gravity. If not dealt with, the accumulation of stress and the resulting tissue deformation can pose a threat to tissue integrity and structure. Here we follow the time-course of events by which cells and tissues return to their preferred state following a mechanical perturbation. In doing so, we discuss the spectrum of biological and physical mechanisms known to underlie mechanical homeostasis in animal tissues.

Celebrating Soft Matter's 10th Anniversary: Cell division: a source of active stress in cellular monolayers

We introduce the notion of cell division-induced activity and show that the cell division generates extensile forces and drives dynamical patterns in cell assemblies. Extending the hydrodynamic models of lyotropic active nematics we describe turbulent-like velocity fields that are generated by the cell division in a confluent monolayer of cells. We show that the experimentally measured flow field of dividing Madin-Darby Canine Kidney (MDCK) cells is reproduced by our modeling approach. Division-induced activity acts together with intrinsic activity of the cells in extensile and contractile cell assemblies to change the flow and director patterns and the density of topological defects. Finally we model the evolution of the boundary of a cellular colony and compare the fingering instabilities induced by cell division to experimental observations on the expansion of MDCK cell cultures.

In silico tumor control induced via alternating immunostimulating and immunosuppressive phases

Despite recent advances in the field of Oncoimmunology, the success potential of immunomodulatory therapies against cancer remains to be elucidated. One of the reasons is the lack of understanding on the complex interplay between tumor growth dynamics and the associated immune system responses. Toward this goal, we consider a mathematical model of vascularized tumor growth and the corresponding effector cell recruitment dynamics. Bifurcation analysis allows for the exploration of model's dynamic behavior and the determination of these parameter regimes that result in immune-mediated tumor control. In this work, we focus on a particular tumor evasion regime that involves tumor and effector cell concentration oscillations of slowly increasing and decreasing amplitude, respectively. Considering a temporal multiscale analysis, we derive an analytically tractable mapping of model solutions onto a weakly negatively damped harmonic oscillator. Based on our analysis, we propose a theory-driven intervention strategy involving immunostimulating and immunosuppressive phases to induce long-term tumor control.

Compressive stress inhibits proliferation in tumor spheroids through a volume limitation

In most instances, the growth of solid tumors occurs in constrained environments and requires a competition for space. A mechanical crosstalk can arise from this competition. In this article, we dissect the biomechanical sequence caused by a controlled compressive stress on multicellular spheroids (MCSs) used as a tumor model system. On timescales of minutes, we show that a compressive stress causes a reduction of the MCS volume, linked to a reduction of the cell volume in the core of the MCS. On timescales of hours, we observe a reversible induction of the proliferation inhibitor, p27Kip1, from the center to the periphery of the spheroid. On timescales of days, we observe that cells are blocked in the cell cycle at the late G1 checkpoint, the restriction point. We show that the effect of pressure on the proliferation can be antagonized by silencing p27Kip1. Finally, we quantify a clear correlation between the pressure-induced volume change and the growth rate of the spheroid. The compression-induced proliferation arrest that we studied is conserved for five cell lines, and is completely reversible. It demonstrates a generic crosstalk between mechanical stresses and the key players of cell cycle regulation. Our results suggest a role of volume change in the sensitivity to pressure, and that p27Kip1 is strongly influenced by this change.

Mechanical induction of the tumorigenic β-catenin pathway by tumour growth pressure

The tumour microenvironment may contribute to tumorigenesis owing to mechanical forces such as fibrotic stiffness or mechanical pressure caused by the expansion of hyper-proliferative cells. Here we explore the contribution of the mechanical pressure exerted by tumour growth onto non-tumorous adjacent epithelium. In the early stage of mouse colon tumour development in the Notch(+)Apc(+/1638N) mouse model, we observed mechanistic pressure stress in the non-tumorous epithelial cells caused by hyper-proliferative adjacent crypts overexpressing active Notch, which is associated with increased Ret and β-catenin signalling. We thus developed a method that allows the delivery of a defined mechanical pressure in vivo, by subcutaneously inserting a magnet close to the mouse colon. The implanted magnet generated a magnetic force on ultra-magnetic liposomes, stabilized in the mesenchymal cells of the connective tissue surrounding colonic crypts after intravenous injection. The magnetically induced pressure quantitatively mimicked the endogenous early tumour growth stress in the order of 1,200 Pa, without affecting tissue stiffness, as monitored by ultrasound strain imaging and shear wave elastography. The exertion of pressure mimicking that of tumour growth led to rapid Ret activation and downstream phosphorylation of β-catenin on Tyr654, imparing its interaction with the E-cadherin in adherens junctions, and which was followed by β-catenin nuclear translocation after 15 days. As a consequence, increased expression of β-catenin-target genes was observed at 1 month, together with crypt enlargement accompanying the formation of early tumorous aberrant crypt foci. Mechanical activation of the tumorigenic β-catenin pathway suggests unexplored modes of tumour propagation based on mechanical signalling pathways in healthy epithelial cells surrounding the tumour, which may contribute to tumour heterogeneity.

Mechanical induction of the tumorigenic β-catenin pathway by tumour growth pressure

The tumour microenvironment may contribute to tumorigenesis owing to mechanical forces such as fibrotic stiffness or mechanical pressure caused by the expansion of hyper-proliferative cells. Here we explore the contribution of the mechanical pressure exerted by tumour growth onto non-tumorous adjacent epithelium. In the early stage of mouse colon tumour development in the Notch(+)Apc(+/1638N) mouse model, we observed mechanistic pressure stress in the non-tumorous epithelial cells caused by hyper-proliferative adjacent crypts overexpressing active Notch, which is associated with increased Ret and β-catenin signalling. We thus developed a method that allows the delivery of a defined mechanical pressure in vivo, by subcutaneously inserting a magnet close to the mouse colon. The implanted magnet generated a magnetic force on ultra-magnetic liposomes, stabilized in the mesenchymal cells of the connective tissue surrounding colonic crypts after intravenous injection. The magnetically induced pressure quantitatively mimicked the endogenous early tumour growth stress in the order of 1,200 Pa, without affecting tissue stiffness, as monitored by ultrasound strain imaging and shear wave elastography. The exertion of pressure mimicking that of tumour growth led to rapid Ret activation and downstream phosphorylation of β-catenin on Tyr654, imparing its interaction with the E-cadherin in adherens junctions, and which was followed by β-catenin nuclear translocation after 15 days. As a consequence, increased expression of β-catenin-target genes was observed at 1 month, together with crypt enlargement accompanying the formation of early tumorous aberrant crypt foci. Mechanical activation of the tumorigenic β-catenin pathway suggests unexplored modes of tumour propagation based on mechanical signalling pathways in healthy epithelial cells surrounding the tumour, which may contribute to tumour heterogeneity.

Stress distributions and cell flows in a growing cell aggregate

We discuss the short-time response of a multicellular spheroid to an external pressure jump. Our experiments show that 5 min after the pressure jump, the cell density increases in the centre of the spheroid but does not change appreciably close to the surface of the spheroid. This result can be explained if the cells are polarized which we show to be the case. Motivated by the experimental results, we develop a theory for polarized spheroids where the cell polarity is radial (except in a thin shell close to the spheroid surface). The theory takes into account the dependence of cell division and apoptosis rates on the local stress, the cell polarity and active stress generated by the cells and the dependence of active stress on the local pressure. We find a short-time increase of the cell density after a pressure jump that decays as a power law from the spheroid centre, which is in reasonable agreement with the experimental results. By comparing our theory to experiments, we can estimate the isotropic compression modulus of the tissue.

Fluorescent correlation spectroscopy measurements with adaptive optics in the intercellular space of spheroids

In this study we demonstrate the use of adaptive optics to correct the biasing effects of optical aberrations when measuring the dynamics of molecules diffusing between cells in multicellular spheroids. Our results indicate that, on average, adaptive optics leads to a reduction of the 3D size of the point spread function that is statistically significant in terms of measured number of molecules and diffusion time. The sensorless approach, which uses the molecular brightness as optimization metric, is validated in a complex, highly heterogeneous, biological environment. This work paves the way towards the design of accurate diffusion measurements of molecules in thick biological specimens.