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

Abstract
Introduction
Results
Discussion and conclusion
Materials and methods
Appendix 1
Data availability
References
Article and author information
Metrics

Abstract

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.

Introduction

Aside from biochemical signaling, cellular function and fate also depend on the mechanical state of the surrounding extracellular matrix (ECM) (Humphrey et al., 2014). The ECM is a non-cellular component of tissues providing a scaffold for cellular adhesion and triggering numerous mechanotransduction pathways, involved in morphogenesis and homeostasis (Vogel, 2018). An increasing number of studies in vivo and in vitro shows that changing the mechanical properties of the ECM by re-implanting tissues or changing the stiffness of the adherent substrate is sufficient to reverse aging (Segel et al., 2019), accelerate developmental processes (Barriga et al., 2018) or modulate tumor malignancy (Paszek et al., 2005; Tanner et al., 2012).

The importance of the mechanical context in cancer has been highlighted for a long time by experiments altering the composition and stiffness of the ECM (Levental et al., 2009). It has also been shown that the tumor growth is modulated by the mechanical compression caused by the tumor itself, as it expands in a confined environment (Fernandez-Sanchez et al., 2010; Nia et al., 2017). Such patho-physiological growth under pressure has also been studied in vitro. When multicellular aggregates are confined by soft gels (Helmlinger et al., 1997; Alessandri et al., 2013; Taubenberger et al., 2019) or submitted to a gentle osmotic compression (Montel et al., 2011; Dolega et al., 2017), their growth is substantially reduced. It has been demonstrated that the cell cytoskeleton is involved in the response to compression and can trigger the growth impediment through a cell-cycle inhibition (Taubenberger et al., 2019; Delarue et al., 2014). In addition, the cellular volume has been recently proposed to be a key parameter in the mechanosensitive pathway (Delarue et al., 2014; Han et al., 2020). Nevertheless, it is not known how such mild global compression is transduced to the individual cells of the aggregate to alter their proliferation.

Here, we posit that cells mainly respond to the mechanical stress transmitted by the ECM, when the aggregate is under compression. This hypothesis is motivated by two evidences. First, an aggregate is a composite material made of cells, extracellular matrix and interstitial fluid. The presence of hydrated extracellular matrix is evidenced by the abundance of fibronectin in the interstitial space (Figure 1a and Appendix 8). As the ECM is 100- to 1000-fold more compressible than the cells, it absorbs most of the deformation, but still transmits the mechanical stress to the cells. Second, whereas an osmotic pressure of a few kPa strongly reduces the cell proliferation within multicellular aggregates, an identical pressure has no effect on individual cells cultured on a Petri dish, in the absence of ECM (Montel et al., 2011). In addition, the use of drugs affecting the cytoskeleton organization has a negligible effect on the effective compressibility of multicellular aggregates (Appendix 12). This indicates that the volume loss under compression is mainly due to ECM dehydration (Dolega et al., 2021).

Figure 1

Download asset Open asset

Selective compression method.

(a) Immunofluorescent staining of fibronectin in the interstitial space of a multicellular spheroid made of CT26 cells (b) Schematic view of a cell (gray) embedded in extracellular matrix (filaments), permeated by interstitial fluid (light pink). (c) Big osmolytes (green) do not penetrate through the ECM and induce a global compression. Being much more compressible than the cells, the extracellular matrix absorbs most of the deformation and exert a positive stress on the cell. (d) Small osmolytes (blue) enter the ECM without exerting any osmotic pressure on it. Conversely, they compress the cell which, in turn, exerts a tension on the ECM.

To test the hypothesis that cells respond to the ECM deformation, we introduce an experimental method that uncouples the cell volume change from the mechanical stress transmitted to the cells through the ECM. We apply this method for both multicellular aggregates and individual cells embedded in a gelified ECM. In parallel, we present a theoretical framework to estimate both the displacement and the stress at the ECM/cell interface in response to an osmotic compression, and verify experimentally its qualitative prediction. At a longer timescale, we probe the effect of the ECM compression on the cellular response. In particular, we demonstrate that, even in the absence of cell deformation, the ECM alone regulates cell proliferation and motility.

Results

Selective-compression method

We developed a simple method to either selectively compress cells embedded in ECM or the whole aggregate composed of ECM and cells. This method is based on the use of osmolytes of different sizes. When big enough, the osmolytes do not infiltrate the ECM and thus compress the whole aggregate by dehydrating the ECM, which in turn mechanically compresses the cells (Monnier et al., 2016). When smaller than the exclusion size of the ECM, the osmolytes percolate through the ECM meshwork and compress the cells which can then pull on the ECM (see schematic in Figure 1b–d and Appendix 7). We already proved (Montel et al., 2012) that a gentle osmotic pressure $Π_{d}$ exerted using large dextran molecules considerably reduces the proliferation of cells inside multicellular spheroids. The effect was visible starting from $Π_{d} = 500 Pa$ and saturated at $Π_{d} ≃ 5 kPa$ . Unless explicitly stated, the experiments described in this article were performed at $Π_{d} ≃ 5 kPa$ , a value that minimizes the pressure, but exacerbates the biological effects.

We validated our approach by compressing ECM, cells and multicellular spheroids (MCS) using osmolytes with gyration radii $R_{g}$ respectively larger and smaller that the ECM pore sizes (Figure 2). As osmolytes, we used dextran molecules ranging from 10 to 2000 kDa. As a proxy of ECM, we used Matrigel (MG), a commercially available matrix secreted by cancer cells (Kleinman and Martin, 2005). To visualize the effect of the compression on the ECM, we prepared microbeads composed of matrigel, with a diameter of 100 μm (Figure 2a ). As shown in Figure 2a (top panel), fluorescent dextran molecules with a gyration radius below 5 nm (MW <70 kDa, hereafter called 'Small’; Granath, 1958) equally color the MG beads and the surrounding solution (left). Conversely, dextran molecules larger than 15 nm (MW >500 kDa, Big’) do not penetrate inside the MG beads, which appear darker than the surroundings (right). By following the evolution of the bead diameter subjected to $Π_{d} = 5 kPa$ (measurements taken before the compression and 45 min after the compression), we observed that small dextran molecules compress the matrigel beads by 2.5 ± 0.7% of their initial volume (Figure 2a, middle and bottom panels). Conversely, the same pressure caused by big dextran molecules occasions a much larger compression of 63 ± 5% (Figure 2b). The relatively minor compression occasioned by small dextran can be explained by thermodynamic theories involving chemical interaction between the matrix and the permeating polymer (Brochard, 1981; Bastide et al., 1981), an aspect that we neglect in this article.

Figure 2

Download asset Open asset

Cell and matrigel compression.

(a) Fluorescently labeled dextran molecules only permeate the beads (top-left panel) if their gyration radius is smaller than 5 nm. Otherwise (top-right panel) they are larger than the exclusion size of the matrigel network and are excluded from the bead. Compression of MG beads, occasioned by dextran molecules of two different sizes (Small: 70 kDa; Big: 500 kDa). Phase contrast images taken before and after the addition of pressure. (b) Beads lose 63 ± 5% of their initial volume when compressed using big dextran, and 2.5 ± 0.7% with small Dextran. N = 10. (c) Compression of individual cells using dextran of different sizes, with $Π_{d} =$ 15 kPa. At $Π_{d} =$ 5 kPa the compressibility of individual cells is not measurable. Cell compressibility is thus negligible in comparison to that of Matrigel. (d) MCS compression under $Π_{d} =$ 5 kPa, exerted using small (blue) and big (green) dextran molecules. In control experiments (0 kPa), the culture medium is replaced by fresh medium without dextran. (box : ±SEOM; error bars: ± SD • : single realizations).

Figure 2—source data 1 Data for Figure 2b,c,d.: https://cdn.elifesciences.org/articles/63258/elife-63258-fig2-data1-v2.xlsx
Download elife-63258-fig2-data1-v2.xlsx

Analogous experiments were performed using individual CT26 cells (murine colon carcinoma cells) and multicellular spheroids made with the same cell line. As the volume loss of individual cells is not measurable at $Π_{d} = 5$ kPa (Monnier et al., 2016), individual CT26 cells are submitted to $Π_{d} = 15$ kPa. At this pressure, we measured a relative compression $Δ V_{c} / V_{c} = 3.8 \pm 0.8 %$ (Figure 2c) , where $V_{c}$ is the cell volume and $Δ V_{c}$ the volume loss upon the application of $Π_{d}$ . This compression indicates that CT26 cells have an effective osmotic modulus $K_{c} =$ 400 ± 100 kPa. In contrast to single cells, MCS are much more compressible, as they lose up to 15% of their volume under an omostic pressure with big dextran of $Π_{d} = 5$ kPa (Figure 2d; See also Dolega et al., 2021 for a detailed mechanical analysis). Furthermore, these measurements indicate that MCS have a typical effective osmotic modulus of $K_{s} ≃$ 30 kPa, 15-folds smaller than that of individual cells (Dolega et al., 2021). In contrast, small dextran molecules have no measurable effect on the volume of MCS, for moderate osmotic pressures (up to $Π_{d} = 10$ kPa). However, larger pressures with these small osmolytes can lead to a cell compression within the MCS associated with a swelling of the interstitial space as we show in Section 'Selective compression of ECM in multicellular spheroids'.

These results confirm the ability of our method to discriminate between the effects occasioned by the compression of the whole MCS, and those due to the compression of the cells alone within the aggregate.

Theory: the effect of a selective compression applied to a cell nested in extracellular matrix

For simplicity, we consider the case of a single cell nested in a large -compared to the cell size- ball of ECM and subjected to the osmotic pressure $Π_{d}$ obtained by supplementing the culture medium with either small or big dextran. We assume that the small dextran can freely permeate in the ECM meshwork while the big one is excluded. Our aim is to compute the displacement of the cell boundary as well as the stress applied on the cell upon application of $Π_{d}$ in both conditions. Our model, detailed in Appendix 6.5, essentially couples a classical active pump-and-leak model (Hoppensteadt and Peskin, 2012) for the cell volume regulation through ion pumping and the constitutive behaviour of the ECM, which is assumed to be poro-elastic at a short timescale where remodeling is negligible. The cell cortex mechanics plays a negligible role in setting the cell volume since it involves stresses that are small compared to the osmotic forces. For simplicity, we neglect the mechano-sensitive nature of ion channels.

We show in Appendix A.5 that, for realistic estimates of the model parameters, the application of $Π_{d}$ with both small or big dextran leads to the same cell volume loss which does not involve the mechanical properties of the ECM but only the cell volume regulation system:

\frac{Δ V_{c}}{V_{c}} = \frac{Π_{d}}{(1 - β) Π_{e}},

where $Π_{e}$ is the osmotic pressure of ions in the culture medium and $β ≃ 0.1$ is a non-dimensional parameter representing the active pumping of ions (see Appendix A.4). The relation (Equation 1) shows that the reduction of the cell volume under compression is mainly resisted by the active osmotic equilibration of ions through the cell membrane. For relatively low pressures ( $Π_{d} ≪ Π_{e} ≃ 500$ kPa), the relative change of volume $Δ V_{c} / V_{c}$ is negligible. More quantitatively, formula (Equation 1) provides the estimate of the osmotic modulus of a cell $K_{c} = (1 - β) Π_{e} ≃ 450$ kPa which is in agreement with the value measured for CT26 cells.

However, the mechanical stress applied by the ECM to the cell is qualitatively and quantitatively different in the two situations. For big dextran, this stress is compressive as the dominating effect of the dextran is to compress the ECM which in turn compresses the cells. Within some realistic approximations the amount of this compressive stress (the traction force applied by the matrix on the cell) can be approximated as the applied osmotic pressure:

T_{big} = - Π_{d} < 0 .

In sharp contrast with the previous situation, for small dextran, the stress applied by the ECM on the cell is tensile. In fact, the dominating effect is that small dextran compresses the cells but not the ECM. Thus, cell compression is balanced by a tensile force in the ECM. This tension is given by

T_{small} = \frac{G Π_{d}}{3 (1 - β) Π_{e}} > 0,

where $G$ is the ECM shear modulus. Formulas (Equations 1, 2 and 3) hold in the ideal case, where osmolytes do not interact with the matrix and the axisymmetric system has stress free boundaries at infinity (the ECM ball radius is much larger than the cell radius).

In practice, for a moderate osmotic shock $Π_{d} ≃ 5$ kPa, the dextran concentration is much smaller than the characteristic ion concentration of the external medium (few hundreds millimolars) and the tension can be considered negligible: $T_{small} ≃ 20 Pa ≪ | T_{big} | ≃ 5 kPa$ because the ECM is soft. Therefore, in this condition, the presence of ECM makes the cell mechanically sensitive to a moderate osmotic compression using big dextran molecules, but not when using small dextran molecules. In both cases the cell volume is affected in the same negligible way, but the mechanical stress applied by small dextran on the cell is negligible compared to that exerted by big dextran.

If the osmotic pressure is further increased, the compression with small or big dextran can induce a measurable effect on the cell volume. However, the mechanical stress applied by the ECM to the cell remains fundamentally different in both situations: tensile for the small dextran and compressive for the big one.

Selective compression of ECM in multicellular spheroids

To test our theoretical predictions that the interstitial space is compressed under dextran pressure, we injected individual MCS (4–5 days old) into a 2D confiner microsystem and let them relax for few hours (Figure 3a). The MCS were thus immobilized and partially flattened inside the 2D confiner. In order to follow the evolution of the interstitial space under an osmotic compression, the culture medium was supplemented with a fluorescent tracer. The interstitial fluorescence was measured using two-photon microscopy (Figure 3b). The images of the confined multicellular aggregates were normalized to the fluorescence of the external medium and segmented with a thresholding procedure, and the signal exceeding the threshold value was integrated over the whole aggregate to quantify the total fluorescence of the interstitial space (Figure 3c). Due to optical limitations, we emphasized the effect by increasing the applied osmotic pressure to $Π_{d} =$ 40 kPa for small dextran and to $Π_{d} =$ 15 kPa for the big ones.

Figure 3

Download asset Open asset

Effect of small versus big dextran on tissue intercellular space.

(a) Schematic of the 2D confiner micro-device. The tissue is confined between the glass coverslip and the PDMS and does not move during medium exchange. (b) Two-photon images of the tissue before and after (20 min) osmotic shocks for dextran chains of 6 kDa (small) and 2MDa (big), for a given mass concentration of 100 g/L. Images were taken in the equatorial plane of the tissue, meaning 35 μm above the glass slide. Scale Bar: 50 μm (c) Mean fluorescence of the intercellular space averaged over the whole aggregate shown in panel b. (d) volume loss of spheroids submitted to $Π_{d}^{S m a l l}$ = 40 kPa (small dextran) and $Π_{d}^{B i g}$ = 15 kPa (big dextran).

Figure 3—source data 1 Data for Figure 3c,d.: https://cdn.elifesciences.org/articles/63258/elife-63258-fig3-data1-v2.xlsx
Download elife-63258-fig3-data1-v2.xlsx

In accordance with our theoretical predictions, we obtained two opposite behaviors, depending on the dextran size. Small dextran molecules induced a ∼35 ± 10% increase in the fluorescence intensity in the interstitial space (Figure 2c) while the total volume of the aggregate was reduced by ∼10% (Figure 2d). Simultaneously, the cell volume decreased (Appendix 9), thus stretching the ECM into occupying more interstitial space. In contrast, for big dextran we measured a loss of half the fluorescence, meaning that a large amount of interstitial liquid had left the intercellular space of the aggregate. The extracellular matrix was thus compressed as predicted by Equation (2) and the overall MCS volume of the whole aggregate was reduced by ∼17% (Figure 2d). We argued in Dolega et al., 2021 that the total volume reduction of the aggregate obtained with big dextran could be due mostly to the compressibility of the ECM, while the cells are quasi-incompressible. The volume reduction of the aggregate induced by a 15 kPa pressure did not differ much from the one obtained with a $Π_{d} = 5$ kPa compression, as the ECM was already fully squeezed at 5 kPa.

These results are consistent with our theoretical prediction that big and small dextran have an opposite effect on the matrix. The first puts the ECM under compression, while the latter puts the ECM under tension. Remarkably, in both cases, the cells within the aggregate undergo almost the same deformation.

ECM compression controls cell proliferation and motility

To understand the role of the ECM on the cell fate at longer timescale, we assessed the proliferation and the motility of cells within MCS cultured in the presence of small and big dextran. Figure 4a represents the equatorial cryosections of spheroids in the three mechanical states ( $Π_{d}$ = 0 kPa, $Π_{d}$ = 5 kPa small dextran, and $Π_{d}$ = 5 kPa big dextran). Proliferating cells were immuno-stained for Ki-67, a nuclear antigen present during the cell cycle, but absent in G0 phase (Gerdes et al., 1984). Whereas cells in control MCS ( $Π_{d} = 0$ Pa) present a rather uniform proliferation pattern, a global compression of MCS (big Dextran) stops cell division in the core and alters the overall MCS growth, as previously reported (Helmlinger et al., 1997; Alessandri et al., 2013; Montel et al., 2011). The density of Ki67-positive cells is reported in panel Figure 4b, as a function of the distance from the spheroid center and for the three conditions represented in panel Figure 4a. Remarkably, under pressure the density of Ki67-positive cells uniformly decreases across the MCS. Consequently, the ratio between the proliferating cells in the periphery of the MCS and those in its core increases under pressure: 2 without pressure, 2.5 under 5 kPa exerted by small dextran and 5 when the same pressure is exerted by big dextran. To quantify the change of cell division rate, we monitored the volumetric growth of the spheroid for three conditions (control, small and big dextran) and for several days (Figure 4c–d). In all cases, the spheroids initially grew exponentially (continuous lines). However, the MCS growth rate (time to double its volume) almost doubled under the big dextran compression, increasing from 36 ± 1 hr for the control and small dextran conditions (gray circles and blue squares, respectively) to 68 ± 4 hr for the compression with big dextran (green triangles).

Figure 4

Download asset Open asset

Growth of spheroids under pressure.

(a) Proliferating cells inside MCS revealed by immunostaining of KI67 with no pressure, under global compression of $Π_{d} = 5$ kPa (big dextran) and under selective compression of the cells by the same amount (small dextran). Scale bar: 100 $μ m$ (b) Density of Ki67-positive cells with respect to the distance from the center of the aggregate. Three conditions: No pressure (9 MCS), $Π_{d} = 5$ kPa with small dextran (26 MCS) and $Π_{d} = 5$ kPa with big dextran (19 MCS). Error bars = standard error of the mean. (c) Time evolution of the spheroid sizes (Full images are 700 × 700 μm) and (d) quantification of the volume increase, in the three reference conditions. (e) Cell migration speed within MCS also significantly depends on ECM compression. N = 5 independent experiments per condition. Error bars represent ± SEM. Experiments were repeated at least on three independent samples. (f) Division time of CT26 cells in 2D (Petri dish), respectively with no pressure (11.4 ± 0.5 hr), with $Π_{d} = 5$ kPa/small dextran (11.0 ± 0.3 hr) and with $Π_{d} = 5$ kPa/big dextran (10.9 ± 0.4 hr). (g) Mean velocity of individual cells on a Petri dish, before (N = 9) and after compression (N = 7).

Figure 4—source data 1 Data for Figure 4b.: https://cdn.elifesciences.org/articles/63258/elife-63258-fig4-data1-v2.mat
Download elife-63258-fig4-data1-v2.mat
Figure 4—source data 2 Data for Figure 4d.: https://cdn.elifesciences.org/articles/63258/elife-63258-fig4-data2-v2.mat
Download elife-63258-fig4-data2-v2.mat
Figure 4—source data 3 Data for Figure 4e.: https://cdn.elifesciences.org/articles/63258/elife-63258-fig4-data3-v2.mat
Download elife-63258-fig4-data3-v2.mat
Figure 4—source data 4 Data for Figure 4f,g.: https://cdn.elifesciences.org/articles/63258/elife-63258-fig4-data4-v2.xlsx
Download elife-63258-fig4-data4-v2.xlsx

Because experiments with MCS are typically performed in solution, where a metastatic behavior is not possible, we evaluated the cell motility within the aggregate, using the Dynamic Light Scattering technique introduced by Brunel et al., 2020 (see details in Appendix 11). The mean migration velocity of cells was reduced by 50% at $Π_{d}$ = 5 kPa with big dextran, as compared to the unstressed case (Figure 4e). Strikingly, both proliferation and motility remained almost unaltered when the MCS were exposed to an equivalent pressure ( $Π_{d} = 5$ kPa) applied by small dextran to selectively compress the cells while leaving the native ECM unstrained (small Dextran, blue).

To verify that neither proliferation nor motility are directly modified by the direct action of dextran in contact with the cells, we measured the proliferation and the velocity of individual cells plated in a Petri dish. Measurements were performed at low density to permit cell proliferation and migration. The results (Figure 4f-g) show that both proliferation and motility remained similar, before and after the addition of dextran at a final pressure $Π_{d} ≃$ 5 kPa.

Since the interstitial space is dehydrated under osmotic compression, cells may get in contact with each other, occasioning contact inhibition of proliferation and locomotion. However, it is also possible that cells sense and react to the stress in the ECM. To discriminate between these two hypotheses, we embeded individual cells in a MG matrix, before compressing the whole system with an osmotic pressure $Π_{d} = 5$ kPa using either small or big dextran. After a few days, we observed two clearly different phenotypes. Cells grown without pressure or in the presence of small dextran were sparse in the MG (Figure 5a, left panel). Conversely, cells cultured with big dextran proliferated locally (Figure 5a, right panel). Therefore, MG compression appears to inhibit cell motility and to promote the formation of mini-spheroids, which suggests that ECM compression has a direct effect on the cell-ECM biophysical signaling. The different cell morphology is particularly clear in the organization of the actin cytoskeleton. Cytoplasmic actin labeling revealed the presence of numerous protrusions, associated with high cell anisotropy, in cells cultured in a relaxed MG matrix (Figure 5b, left and middle panels), whereas cells appeared smooth and formed round structures, when the MG was compressed (Figure 5b, right panel). Of note, cells at the MG surface often extended outside the MG. Those cells not fully embedded in MG were excluded from our analysis.

Figure 5

Download asset Open asset

Individual cells in Matrigel.

(a) Hoechst-labeled cell nuclei superimposed to phase image. Images are taken after 2 days of proliferation in MG, either with small (left panel) or big (right panel) dextran molecules. Maximal projection from epifluorescence stacks. (b) Cell morphology and anisotropy revealed by labeling of cytoplasmic actin. Maximal projection of 50 μ m confocal Z-stack. In relaxed MG, the cells appear more elongated and with long protrusions. (c) Cell motility in MG under different compression states. Starting points of trajectories are translated to the origin, to highlight the typical distance over which cells move in the three compressive states. (d) Quantification of in-plane velocity extracted from mean square displacements, under different compression conditions. With no pressure or with small dextran (5 kPa), the average velocities are respectively 5.8 ± 0.8 $μ m$ / hr and 5.2 ± 0.5 $μ m$ / hr. Under 5 kPa exerted by big dextran, the cells are immobile (v = 0.5 ± 0.4 $μ m$ ), where the error is due to tracking uncertainties. (e) Temporal evolution of nuclear fluorescence intergrated over the whole sample. No pressure (ο), 5 kPa with small dextran ( $□$ ) and 5 kPa with big dextran ( $△$ ). (f) Cell proliferation rate in the three conditions. n = 15, from eight independent experiments. (g) Cell proliferation rate at different initial matrigel concentration, with no pressure. Boxes represent the mean values ± SEM, error bars correspond to the standard deviation, small markers are individual experiments and large markers the median.

Figure 5—source data 1 Data for Figure 5c,d.: https://cdn.elifesciences.org/articles/63258/elife-63258-fig5-data1-v2.mat
Download elife-63258-fig5-data1-v2.mat
Figure 5—source data 2 Data for Figure 5e,f,g.: https://cdn.elifesciences.org/articles/63258/elife-63258-fig5-data2-v2.xlsx
Download elife-63258-fig5-data2-v2.xlsx

These different morphologies also correlates with different motilities. Cells embedded in a compressed MG were nearly immobile, while they migrated through relaxed MG with a velocity comparable to that measured on flat surfaces. The results are summarized in Figure 5c, where we report $≃ 40$ trajectories per condition. To highlight differences and similarities between the three compression conditions, the starting points of all trajectories are translated to the origin and, although isotropic, they are divided in three quadrants. Quantification is reported in Figure 5d. From this experiment we conclude that whereas no appreciable differences are observable between control and the small dextran condition, the cell motility dramatically drops under MG compression with big dextran.

To quantify the effect of ECM compression on proliferation, we prepared several samples with the same number of hoechst-stained cells embedded in the MG and measured the overall fluorescence over time. Figure 5e shows the typical time evolution of the Hoechst signal for the three conditions: proliferation rate drops considerably when the MG is compressed (Big Dextran, $△$ ), compared to the case without pressure (ο), but also compared to the case where the pressure is selectively exerted on the cells with no MG compression (Small Dextran, $□$ ). Figure 5f quantifies the mean growth rate, measured on at least 15 samples for each condition, collected on eight independent experiments (different days and cell passages). Under pressure, the matrigel get dehydrated and compacted, which may directly influence cell proliferation. On the one hand a denser matrigel is less compressible and, thus, less favorable to cell proliferation (Baker et al., 2015). On the other hand, matrigel compression concentrates matrix-bound growth factors, which may promote cell division. To determine which effect dominates, we measured the proliferation rate at different initial MG concentrations, between 2 g/l and 8 g/l (experiments reported in panels a-f were performed with MG at 4.5 g/l). Our experiments show that within this range, the matrigel density has little to no effect on cell proliferation (see Figure 5g). These experiments confirmed that it is the compressive stress transmitted to the cells by the surrounding ECM, rather than a direct osmotic pressure on the cells, that strongly impacts cell motility and proliferation.

Discussion and conclusion

Large osmotic and mechanical pressures (of the order of 100 kPa) can cause a decrease in cell volume and consequently a deformation of the cell nucleus (Zhou et al., 2009; Kim et al., 2015) which may ultimately feedback on the cell proliferation. It has been recently proposed that the volume of the cell or its nucleus can be key to crucial processes such as proliferation, invasion, and differentiation Guo et al., 2017; Han et al., 2020. However the weak osmotic pressures (of the order of 1 kPa) that we apply have no measurable effect on the cell volume. In addition, it is well-known from a biological standpoint that such small volume perturbations are buffered by active regulatory processes in the cell (Hoffmann et al., 2009; Cadart et al., 2019). Yet, both cell proliferation and motility decrease in MCS submitted to weak osmotic compression. Our results show that, for such weak compressions, the cell volume is unchanged while the ECM located in between the cells is directly impacted. This mechano-sensitive role of the ECM could explain the reported evidences that osmotic pressures applied by big dextran and mechanical pressures similarly affect the growth of MCS. (Helmlinger et al., 1997; Alessandri et al., 2013; Montel et al., 2011). Indeed, in this case the osmotic pressure induces a mechanical one applied on the cells through the ECM drainage. Thanks to its bulk modulus $K_{E C M} ≃$ 1 kPa, the ECM behaves as a pressure sensor for the cell in the kPa range. Of note, stress relaxation in the ECM could occur through cleavage and remodeling of its components and such active processes should be quantified in the future.

Several mechanisms may explain how the dehydration of the extracellular matrix can result in an inhibition of proliferation and motility. First, the reduction of the interstitial space promotes interactions between neighbouring cells, which may activate contact inhibition signals of both proliferation and locomotion (Roycroft and Mayor, 2016). Second, the ECM porosity and tortuosity change within a compressed MCS, such that its effective permeability to oxygen, nutrients, growth factors and cytokines is reduced and might activate inhibition signals without cell-cell contact. However, both options are incompatible with the results we obtained with single cells embedded in MG (Figure 5). During the first 2–3 days after seeding, cells are either isolated or grouped in aggregates of two to four cells with limited cell-cell contacts. Additionally, a key factor limiting the diffusion of oxygen and nutrients in MCS is the tortuosity of the interstitial space (Bläßle et al., 2018). This constraint is simply absent in experiments with single cells embedded in MG, suggesting that the cell proliferation inhibition is most probably not related to hypoxia and starvation.

The present work therefore points at a direct mechanosensitive response of cells to the ECM deformation. The microscopic structure of the ECM is modified under compression (e.g. density increase and reduction of porosity), with consequences on the ECM rheology. Compression of the ECM is clearly accompanied by an increase in its bulk modulus and, due to its fibrillar structure, to a non-trivial and non-linear evolution of its stiffness (Sopher et al., 2018; Kurniawan et al., 2016). For example, the rheological properties of synthetic ECM have been shown to affect growth of aggregates and single cells through the regulation of streched-activated channels (Nam et al., 2019). As integrin-dependent signals and focal adhesion assembly are regulated by the stress and strain between the cell and the ECM, the osmotic compression may steer the fate of cells in terms of morphology, migration, and differentiation (Pelham and Wang, 1997; Choquet et al., 1997; Sunyer et al., 2016; Isenberg et al., 2009; Butcher et al., 2009; Engler et al., 2006; Staunton et al., 2019; Panzetta et al., 2019). This aspect is also relevant from an oncological point of view. Indeed the ECM is strongly modified in tumour tissues and the solid stress within tumors can reach several kPa, which is in accordance with the pressure applied here (Nia et al., 2017). For example, in tumors, there is a decrease in the ratio collagen/hyaluronan (Voutouri et al., 2016). The latter, more hydrophilic than the first one, tends to swell and stiffen the ECM. Whether a corrupted matrix is a contributing cause or the consequence of the neoplasia remains an open question, but the correlation between matrix mechanics and uncontrolled proliferation is more and more widely accepted (Bissell et al., 2002; Lelièvre and Bissell, 2006; Broders-Bondon et al., 2018).

In future experiments it will be crucial to identify whether the ECM compression and the associated changes in stiffness play a dominant role, or if – as we suggest – the mechanical stress applied on the cell through the ECM is the key ingredient directly triggering the cell biological adaptation in term of proliferation and motility.

Materials and methods

Key resources table

Reagent type (species) or resource	Designation	Source or reference	Identifiers	Additional information
Cell Line (mouse)	CT26	ATCC CRL-2638	RRID:CVCL_7256
Chemical compound,drug	Matrigel/MG	Corning; 354234
Chemical compound, drug	Small Dextran	Sigma Aldrich; D9260
Chemical compound, drug	Big Dextran	Sigma-Aldrich; D5376
Antibody	Anti-Fibronectin (Monoclonal mouse)	Sigma Aldrich; F7387	RRID:AB_476988	(1:200)
Antibody	Anti-Ki-67 (Polyclonal rabbit)	Sigma Aldrich; AB9260	RRID:AB_2142366	(1:200)

Cell culture, MCSs formation, and growth under mechanical stress

Request a detailed protocol

CT26 (mouse colon adenocarcinoma cells, ATCC CRL-2638); American Type Culture Collection were cultured under 37°C, 5% CO₂ in DMEM supplemented with 10% calf serum and 1% antibiotic/antimycotic (culture medium). Cells are texted every month for mycoplasma. None of the experiments was made using cells with mycoplasma. Spheroid were prepared on agarose cushion in 96 well plates at the concentration of 500 cell/well and centrifuged initially for 5 min at 800 rpm to accelerate aggregation. After 2 days, Dextran (molecular mass 1, 10, 40, 70, 100, 200, 500, and 2000 kDa; Sigma-Aldrich, St. Louis, MO) was added to the culture medium to exert mechanical stress, as previously described (Monnier et al., 2016). To follow spheroid growth over the time, phase contrast images were taken daily. Spheroid were kept under constant pressure over observation period. Images were analysed manually using Imagej. Each experiment was repeated three times, with 32 individual spheroids per condition.

Measurement of MCSs volume

Request a detailed protocol

The area of the MCS equatorial section was measured before and after addition of dextran, then converted to volume assuming that the MCS is spherical. To induce compression, half of the culture medium was removed and replaced with fresh medium containing dextran 2X (two-fold the target concentration). The error affecting this measurement mainly comes from the fact that spheroids rotate during buffer exchange. As they are not perfectly spherical, the area of the equatorial section may change by up to few percent. To homogenize the experiments, control spheroids (no pressure) were also measured before and after buffer exchange. In the latter case, 50% of the culture medium was simply aspirated and replaced by fresh medium not supplemented with dextran.

Fabrication of Matrigel beads

Request a detailed protocol

Matrigel beads (Matrigel Corning, Ref: 354234) were prepared using vortex method (Dolega et al., 2017). Oil phase of HFE-7500/PFPE-PEG (1.5% w/v) was cooled down to 4°C. For 400 μL of oil, 100 μL of Matrigel were added. Solution was vortexed at full speed for 20 s and subsequently kept at 37°C for 20 min for polymerization. Beads were eventually transferred to PBS phase by washing out the surfactant phase.

Fluorescence eXclusion method (single-cell volume measurements)

Request a detailed protocol

Cell volume was obtained using Fluorescence Exclusion microscopy (Cadart et al., 2017; Zlotek-Zlotkiewicz et al., 2015). Briefly, cells were incubated in PDMS chips, with medium supplemented with a fluorescent dye that does not enter the cells. Cells thus excluded fluorescence, and one extracted cellular volume by integrating the fluorescence intensity over the whole cell . Chips for volume measurements of single cells were made by pouring a mixture (1:10) of PMDS elastomer and curing agent (Sylgard 184) onto a brass master and cured at 80 °C for at least of 2 hr. Inlet and outlets were punched with a 3 mm biopsy puncher. Chips were prepared few days before, bounded with oxygen plasma for 30 s, warmed up at 80°C for 3 min then incubated with Poly-l-lysine (sigma) for 30 min to 1 hr, washed with PBS, then washed with dH₂O, dried and stored sealed with a paraffin film. The chambers were washed with PBS before cell injection. Imaging started within 10 min after cell injection in order to prevent adhesion and thus cells response to the shear stress generated by the medium exchange. Acquisition was performed at 37°C in CO₂ independent medium (Life Technologies) supplemented with 1 g/L FITC dextran (10 kDa, from Sigma Aldrich) on an epifluorescence microscope (Leica DMi8) with a 10x objective (NA. 0.3 from LEICA). Master molds were fabricated on a brass substrate with a micromilling machine (MiniMill/3; Minitech) using a 100-μm-diameter milling cutter (Minitech). Height profiles and surface roughness were measured with a vertical scanning interferometric profilometer (Brucker). 3D mold design and tool paths were generated using Autodesk Inventor Professional software (Autodesk). Molds for spheroid confinement were made with classical soft lithography techniques.

Tissue compression experiments

Request a detailed protocol

Spheroids were harvested 4 or 5 days after cell seeding and injected in the 2D confiner microsystem (Figure 3a) using a MFCS pressure controller (Fluigent). Spheroid were partially flattened between two parallel surfaces, perpendicular to the optical axis of the microscope, and rested for two to 5 hr to relax in the microsystem at 37°C in CO₂ independent medium. Before acquisition, medium supplemented with 2 g/L FITC-dextran (10 kDa from Sigma Aldrich) was injected to label the intercellular space. Medium exchange was performed manually using large inlets (<1 mm) during two-photon acquisition. Acquisitions were performed at 37°C on a Nikon C1 two-photon microscope coupled with a femtosecond laser at 780 nm with a 40x water-immersion (NA. 1.10) objective (Nikon). The 2D confiner chip was made by pouring PDMS elastomere and curing agent (1:10) into the mold and cured for at least 2 hr. The chips were bounded to glass coverslips with 30 s oxygen plasma, immediately after bounding. A solution of PLL-g-PEG (Surface Solutions) at 1 g/l was injected and incubated for 30 min in humid atmosphere to prevent cell surface adhesion during the experiment. The chips were washed with $d H_{2} O$ and dried and sealed with a paraffin film. Fluorescence of the Intercellular space (ICS) was measured using MatLab software. As control and dextran solutions have different levels of fluorescence, the fluorescence in the ICS was normalized by the one outside the ICS in order to compensate these variations. Then, first the tissue (cells and ICS) was segmented with a thresholding procedure. The threshold was determined in order to obtained accurate segmentation of the ICS before application of the osmotic stress. The surface of the ICS was computed as the ratio of pixels in the ICS to the number of pixels of the tissue. For each spheroid, 50 planes - 13 µm above and below the equatorial plane - were taken into account to compute the change of the ICS surface.

Cell culture in Matrigel

Request a detailed protocol

Experiments have been conceived to start the culture from individual cells embedded in Matrigel. At day 1, the cells were resuspended, then dispersed in a solution containing matrigel at the final concentration of 4.5 g/l. The cells were diluted to 10,000–50,000 cells/ml, a concentration at which the average distance between neighboring cells is about 250–400 μm. We therefore consider them as isolated entities. The MG/cell ensemble was gelified in 200 μl wells, at 37°C, for 30 min. To avoid cell sedimentation, we gently flipped the sample over, every two minutes. The samples were then redeposited in the incubator under three pressure conditions: no pressure and 5 kPa exerted by small and 5 kPa exerted big dextran.

Cells migration in Matrigel

Request a detailed protocol

To quantify cell migration in Matrigel, individual cells were observed by phase contrast microscopy. Z-stacks were collected every 20 min and for several days, with slices spaced by 50 μm. Then the full stack was projected to one single layer (maximum intensity projection). Cells were tracked manually in the 2D plane, using the ImageJ MTrackJ plugin (https://imagescience.org/meijering/software/mtrackj/).

Cryosectioning and immunostaining

Request a detailed protocol

Spheroids were fixed with 5% formalin (Sigma Aldrich, HT501128) in PBS for 30 min and washed once with PBS. For cryopreservation spheroids were exposed to sucrose at 10% (w/v) for 1 hr, 20% (w/v) for 1 hr and 30% (w/v) overnight at 4°C. Subsequently spheroids were transferred to a plastic reservoir and covered with Tisse TEK OCT (Sakura) in an isopropanol/dry ice bath. Solidified samples were brought to the cryotome (Leica CM3000) and sectioned into 15 μm slices. Cut layers were deposited onto poly-L-lysine coated glass slides (Sigma) and the region of interest was delineated with DAKO pen. Samples were stored at −20°C prior immunolabelling. For fibronectin and Ki67 staining samples were permeabilized with Triton X 0.5% in TBS (Sigma T8787) for 15 min at RT. Nonspecific sites were blocked with 3% BSA (Bovine serum Albumin) for 1 hr. Then, samples were incubated with first antibody (Fibronectin, Sigma F7387, 1/200 and Ki67; Millipore ab9260, 1/500) overnight at 4°C. Subsequently samples were thoroughly washed with TBS three times, for 15 min each. A second fluorescent antibody (goat anti-mouse Cy3, Invitrogen; 1/1000) was incubated for 40 min along with phalloidin (1/500, Alexa Fluor 488, Thermo Fisher Scientific). After extensive washing with TBS (four washes of 15 min) glass cover slides were mounted on the glass slides with a Progold mounting medium overnight (Life Technologies P36965) and stored at 4°C before imaging.

Statistical analysis

Request a detailed protocol

Student’s t-test (unpaired, two tailed, equal variances) was used to calculate statistical significance as appropriate by using the ttest2() function of Matlab (MathWorks). Statistical significance is given by *p<0.05; **p<0.01; ***p<0.001; ****p<0.0001.

Appendix 1

A Theoretical model of the osmotic compression of a single cell nested in matrigel

Our aim is to qualitatively understand the nature of the steady state mechanical stress and displacement of a cell nested in a matrix in two paradigmatic situations:

When some small osmolytes (typically dextran) that can permeate the matrix pores are introduced in the solution,
When some big osmolytes that are excluded from the matrix are introduced in the solution.

The matrix is a meshwork of biopolymers permeated by an aquaeous solution containing ions. These ions can also permeate the cell cytoplasm via specific channels and pumps integrated in the plasmic membrane (Hoffmann et al., 2009; Lang et al., 1998). For simplicity, we restrict our theoretical description to Na⁺, K⁺, and Cl^– ions which have specific channels and a well studied pump (Therien and Blostein, 2000) which actively pumps out three sodium ions in exchange of having two potassium ions in. Attached right under the cell membrane via some specific cross-linkers (Diz-Muñoz et al., 2010), the cell cortex is a thin ’muscle-like’ actin network cross-linked by passive and contractile cross-linkers such as myosin II. The cortex has been shown to be an important regulator of the cell surface tension (Clark and Paluch, 2011; Salbreux et al., 2012) as exemplified during motility (Hawkins et al., 2011; Farutin et al., 2019) and cell morphogenesis (Turlier et al., 2014; Sedzinski et al., 2011; Tinevez et al., 2009; Charras et al., 2008). The cell membrane and cortex enclose the cytoplasm a meshwork of macromolecules permeated by water and containing the aforementioned ions. See Appendix 1—figure 6 for a scheme of the model.

Appendix 1—figure 1

Download asset Open asset

Scheme of a cell nested in a porous matrix.

For simplicity we assume a spherical geometry with a cell of radius r_c inside a matrix ball of radius r_m. Each point in the space $x$ can therefore be localized by its radial position $x = r e_{r}$ where $e_{r}$ is radial unit vector. We assume a spherical symmetry of the problem such that all the introduced physical fields are independent of the angular coordinates θ and φ. Throughout this text, we restrict ourselves to a linear theory which typically holds when the deformation in the matrix is assumed to remain sufficiently small. A more quantitative theory would require to take into account both the non-linear aspects of the matrix deformation and the osmotic pressure created by the polymer.

A.1 Conservation laws at the cell-matrix interface

Water conservation

From Kedem-Katchalsky theory (Staverman, 1952; Kedem and Katchalsky, 1958; Kedem and Katchalsky, 1963; Baranowski, 1991; Elmoazzen et al., 2009), assuming that the aquaeous solvent moves through specific and passive channels, the aquaporins (Day et al., 2014), we can express the incoming water flux $𝐣_{w}$ in the cell at $r = r_{c}$ as (Yi et al., 2003; Hui et al., 2014; Strange, 1993; Hoffmann et al., 2009; Mori, 2012; Cadart et al., 2019):

𝐣_{w} . 𝐞_{r} = L_{p} [p_{m} - p_{c} - (Π_{m} - Π_{c})],

where $Π_{m, c}$ denote the osmotic pressures in the matrix phase and the cell while $p_{m, c}$ are the hydrostatic pressures defined with respect to the external (i.e. atmospheric) pressure. The so-called filtration coefficient $L_{p}$ is related to the permeability of aquaporins. In a dilute approximation which we again assume for simplicity, the osmotic pressure is dominated by the small molecules in solution and thus reads

Π_{m} = k_{B} T (N_{m} + K_{m} + C_{m} + D_{m}) a n d Π_{c} = k_{B} T (N_{c} + K_{c} + C_{c}),

where $k_{B}$ is the Boltzmann constant, $T$ the temperature, $N_{c, m}$ , $K_{c, m}$ and $C_{c, m}$ are the (number) concentrations of sodium, potassium and chloride in the cytoplasm and the extra-cellular medium and $D_{m}$ is the extra-cellular Dextran (necessary small as big are excluded) concentration in the matrix phase. We neglect in (Equation 5) the osmotic contribution associated with the large macromolecules composing the cell organelles and the cytoskeleton compared to the ionic contributions. In a similar manner, the osmotic contribution of the matrix polymer is also neglected. At steady state, the water flux vanishes ( $𝐣_{w} = 0$ ) leading to the relation at $r = r_{c}$ ,

p_{m} - p_{c} = Π_{m} - Π_{c} .

Ions conservation

As each ion travels through the plasma membrane via specific channels and pumps, the intensities of each ionic current at $r = r_{c}$ is given by Nernst-Planck laws (Mori, 2012),

\begin{matrix} i_{N} = g_{N} [v_{c} - \frac{k_{B} T}{q} \log (\frac{N_{m}}{N_{c}})] + 3 q j_{p} \\ i_{K} = g_{K} [v_{c} - \frac{k_{B} T}{q} \log (\frac{K_{m}}{K_{c}})] - 2 q j_{p} \\ i_{C} = g_{C} [v_{c} + \frac{k_{B} T}{q} \log (\frac{C_{m}}{C_{c}})], \end{matrix}

where $g_{N, K, C}$ are the respective conductivities of ions, v_c is the cell membrane potential, $q$ is the elementary charge and j_p is the pumping rate associated to the Na-K pump on the membrane which is playing a fundamental role for cellular volume control (Hoffmann et al., 2009). The factors 3 and 2 are related to the stochiometry of the sodium potassium pump. Again, in steady state, currents $i_{N, K, C} = 0$ , leading to the Gibbs-Donnan equilibrium:

N_{c} = N_{m} e^{- \frac{q (v_{c} - v_{N})}{k_{B} T}}, K_{c} = K_{m} e^{- \frac{q (v_{c} - v_{K})}{k_{B} T}} a n d C_{c} = C_{m} e^{\frac{q v_{c}}{k_{B} T}},

where the active potentials related to the pumping activity $v_{N, K}$ are $v_{N} = - 3 q j_{p} / g_{N}$ and $v_{K} = 2 q j_{p} / g_{K}$ .

Supposing that the cell membrane capacitance is vanishingly small (Mori, 2012), we can neglect the presence of surface charges and impose an electro-neutrality constraint for the intra-cellular medium:

\begin{matrix} N_{c} - C_{c} + K_{c} - ρ_{c} z_{c} = 0, \end{matrix}

where z_c is the average number of (negative in the physiological $pH = 7.4$ conditions) electric charges carried by macromolecules inside the cell and $ρ_{c}$ is their density. As macromolecules are trapped inside the cell membrane, we can express $ρ_{c} = X_{c} / (4 π r_{c}^{3} / 3)$ where $X_{c}$ is the number of macro-molecules which is fixed at short timescale and only increases slowly through synthesis as the amount of dry mass doubles during the cell cycle (Cadart et al., 2019).

Force balance

At the interface between the cell and the matrix $(r = r_{c})$ , we can express the mechanical balance as

𝚺_{c}^{bulk} 𝐞_{r} + σ_{c}^{surf} 𝐞_{r} = 𝚺_{m} 𝐞_{r} .

In (Equation 10), $𝚺_{c}^{bulk}$ is the Cauchy stress in the cytoplasm which we decompose into $𝚺_{c}^{bulk} = 𝚺_{c}^{skel} - p_{c} 𝐈$ , with a first contribution due to the cytoskeleton and a second contribution due to the hydrostatic pressure in the cytosol. The identity matrix is denoted I. The contribution due to the mechanical resistance of the cortex and membrane is denoted $σ_{c}^{surf}$ . In our spherical geometry, we can express $σ_{c}^{surf} = 2 γ_{c} / r_{c}$ where $γ_{c}$ is a surface tension in the cell contour. Finally $𝚺_{m}$ is the stress in the matrix phase for which we postulate a poro-elastic behavior such that, $𝚺_{m} = 𝚺_{m}^{el} (ϵ_{m}) - p_{m} 𝐈$ (the Biot coefficient [Coussy, 2004] is assumed to be one). where

𝚺_{m}^{el} = 2 G 𝐄_{m} + (K_{d} - \frac{2 G}{3}) tr (𝐄_{m}) 𝐈,

is the Hooke’s law with $𝐄_{m}$ the (small) elastic strain in the matrix, $G$ the shear modulus and $K_{d}$ the drained bulk modulus.

In the absence of cytoskeleton and external matrix (Equation 10) reduces to Laplace law:

\frac{2 γ_{c}}{r_{c}} = p_{c} - p_{m}

and more generally reads,

(𝚺_{c}^{skel} - 𝚺_{m}^{el}) 𝐞_{r} . 𝐞_{r} + \frac{2 γ_{c}}{r_{c}} = p_{c} - p_{m} .

Such relation provides the hydrostatic pressure jump at the cell membrane $(r = r_{c})$ entering in the osmotic balance (Equation 6) and, combining (Equation 6) and (Equation 12), we obtain

(𝚺_{c}^{skel} - 𝚺_{m}^{el}) 𝐞_{r} . 𝐞_{r} + \frac{2 γ_{c}}{r_{c}} = Π_{c} - Π_{m} .

Conservation laws in the extracellular matrix

Water conservation

Assuming that the extracellular fluid follows a Darcy law, mass conservation of the incompressible water permeating the matrix can be expressed as

\frac{\partial n}{\partial t} - \frac{κ}{μ} \frac{1}{r} \frac{\partial}{\partial r} (r \frac{\partial p_{m}}{\partial r}) = 0,

where $n$ is the matrix porosity, κ the matrix permeability and μ the fluid viscosity. At steady state, $\partial_{t} n = 0$ and (Equation 14) is associated with no flux boundary conditions at r_c and r_m given by

{\frac{\partial p_{m}}{\partial r} |}_{r_{m}, r_{c}} = 0 .

It follows that p_m is homogeneous in the matrix and its value is imposed by a relation similar to (Equation 6) with an infinitely permeable membrane at r_m :

p_{m} (r) = Π_{m} - Π_{e} .

In (Equation 15), $Π_{e}$ is the external osmotic pressure which reads

Π_{e} = k_{B} T (N_{e} + K_{e} + C_{e} + D_{e})

where $N_{e}$ , $K_{e}$ and $C_{e}$ denote the ions concentrations in the external solution and $D_{e}$ the concentration of Dextran added to the external solution.

Ions conservation

As we are interested in the steady-state only, the Poisson-Nernst fluxes of ions concentrations in the matrix locally vanish leading to:

\frac{\partial N_{m}}{\partial r} + \frac{N_{m} q}{k_{B} T} \frac{\partial v_{m}}{\partial r} = \frac{\partial K_{m}}{\partial r} + \frac{K_{m} q}{k_{B} T} \frac{\partial v_{m}}{\partial r} = \frac{\partial C_{m}}{\partial r} - \frac{C_{m} q}{k_{B} T} \frac{\partial v_{m}}{\partial r} = 0,

where $v_{m} (r)$ is the electro-static potential in the matrix.

As v_m is defined up to an additive constant, we chose that $v_{m} (r_{m}) = 0$ and, imposing the continuity of ions concentrations at the transition between the matrix and the external solution ${N_{m} |}_{r_{m}} = N_{e}$ , ${K_{m} |}_{r_{m}} = K_{e}$ and ${C_{m} |}_{r_{m}} = C_{e}$ , we obtain

N_{m} = N_{e} e^{- \frac{q v_{m}}{k_{B} T}}, K_{m} = K_{e} e^{- \frac{q v_{m}}{k_{B} T}} a n d C_{m} = C_{e} e^{\frac{q v_{m}}{k_{B} T}} .

Next, we again suppose for simplicity that the capacitance of both the porous matrix and the external media are vanishingly small leading to the electro-neutrality constraints

\begin{matrix} N_{m} + K_{m} - C_{m} - z_{m} ρ_{m} = 0 \\ N_{e} + K_{e} - C_{e} = 0, \end{matrix}

where z_m is the number of negative charges carried by the biopolymer chains forming the matrix and $ρ_{m}$ is their density. As we use uncharged Dextran, its concentration does not enter in expressions (Equation 18). Using, (Equation 17) in tandem with (Equation 18), we obtain

v_{m} = - \frac{k_{B} T}{q} \sinh^{- 1} (\frac{z_{m} ρ_{m}}{2 C_{e}}) .

Re-injecting this expression into (Equation 17), we obtain the steady state concentrations of ions in the matrix phase:

N_{m} = N_{e} e^{\sinh^{- 1} (\frac{z_{m} ρ_{m}}{2 C_{e}})}, K_{m} = K_{e} e^{\sinh^{- 1} (\frac{z_{m} ρ_{m}}{2 C_{e}})} a n d C_{m} = C_{e} e^{- \sinh^{- 1} (\frac{z_{m} ρ_{m}}{2 C_{e}})} .

Next, we make the realistic assumption that the chloride concentration (number of ions per unit volume) is much larger than the density of fixed charges carried by the polymer chains (number charges per unit volume): $z_{m} ρ_{m} / C_{e} ≪ 1$ . Indeed using the rough estimates of Section 'Cell volume in the reference situation', the average number of charge carried per amino acid is 0.06 and the typical concentration of matrix is 5 g/l. As the molar mass of an amino acid is roughly 150 g/mol, we can estimate in moles that $z_{m} ρ_{m} ≃ 2$ mM while $C_{e} ≃ 100$ mM. We can thus simplify (Equation 20) up to first order to obtain,

N_{m} = N_{e} (1 + \frac{z_{m} ρ_{m}}{2 C_{e}}), K_{m} = K_{e} (1 + \frac{z_{m} ρ_{m}}{2 C_{e}}) a n d C_{m} = C_{e} (1 - \frac{z_{m} ρ_{m}}{2 C_{e}}) .

As a result, we obtain that the only steady state contribution of

Π \overset{def}{=} Π_{e} - Π_{m} = k_{B} T (D_{e} - D_{m}) = {\begin{array}{lr} 0 & f o r s m a l l D e x t r a n \\ k_{B} T D_{e} & f o r b i g D e x t r a n, \end{array}

is imposed by Dextran since the ions only start to contribute to this difference at second order in the small parameter $z_{m} ρ_{m} / C_{e}$ . We therefore conclude that, in good approximation, $Π$ vanishes for small Dextran molecules that can permeate the matrix and equates to the imposed and known quantity $k_{B} T D_{e}$ for big Dextran molecules that cannot enter the matrix pores.

It then follows from Equation 15 that the hydrostatic pressure equilibrates with the imposed osmotic pressure,

p_{m} (r) = - Π .

Force balance

Using the spherical symmetry of the problem, the only non vanishing components of the stress tensor are $𝚺_{m}^{r r}$ and $𝚺_{m}^{θ θ} = 𝚺_{m}^{φ φ}$ . Therefore, the local stress balance reads

\frac{\partial 𝚺_{m}^{r r}}{\partial r} + \frac{2}{r} (𝚺_{m}^{r r} - 𝚺_{m}^{θ θ}) = 0,

Assuming a small enough displacement, the non-vanishing components of the strain tensor are given by, $𝐄_{m}^{r r} = \partial u_{r} / \partial r$ and $𝐄_{m}^{θ θ} = 𝐄_{m}^{φ φ} = u_{r} / r$ where u_r is the radial (and only non-vanishing) displacement component from an homogeneous reference configuration corresponding to a situation where the matrix is not subjected to any external loading and $r_{c, m} = R_{c, m}$ . Using the poro-elastic constitutive behavior (Equation 11), u_r satisfies

(K_{d} + \frac{4}{3} G) (\frac{\partial^{2} u_{r}}{\partial r^{2}} + \frac{2}{r} \frac{\partial u_{r}}{\partial r} - \frac{2}{r^{2}} u_{r}) = \frac{\partial p_{m}}{\partial r} .

This equation is supplemented with the traction free boundary condition at $r = r_{m}$

𝚺_{m} 𝐞_{r} = 0 .

Combined with (Equation 23), the two above Equations (24) and (25) lead to the solution

u_{r} (r) = ϵ^{0} r + \frac{r_{m}^{3} (Π + 3 K_{d} ϵ^{0})}{4 G r^{2}},

where the introduced constants $ϵ^{0}$ is found using the displacement continuity at the cell matrix- interface:

u_{r} (r_{c}) = u \overset{def}{=} r_{c} - R_{c},

with $u$ given by the change of the cell radius from a reference configuration with radius $R_{c}$ . The general expression of u_r therefore reads,

u_{r} (r) = \frac{u r_{c}^{2} (4 G r^{3} + 3 K_{d} r_{m}^{3}) + Π r_{m}^{3} (r_{c}^{3} - r^{3})}{r^{2} (4 G r_{c}^{3} + 3 K_{d} r_{m}^{3})},

leading to the following form of the total mechanical stress in the surrounding matrix:

Σ_{m} (r) = \frac{2 G r_{c}^{2} (3 K_{d} u + r_{c} Π)}{r^{3} (4 G r_{c}^{3} + 3 K_{d} r_{m}^{3})} \times (\begin{array}{ccc} 2 (r^{3} - r_{m}^{3}) & 0 & 0 \\ 0 & (2 r^{3} + r_{m}^{3}) & 0 \\ 0 & 0 & (2 r^{3} + r_{m}^{3}) \end{array}) .

A.3 Formulation of the model

Combining (Equation 5) with (Equation 13) and taking into account (Equation 21), we obtain the relation linking the cell mechanics and the osmotic pressures inside the cell and outside the matrix:

(𝚺_{c}^{skel} - 𝚺_{m}^{el}) 𝐞_{r} . 𝐞_{r} + \frac{2 γ_{c}}{r_{c}} = k_{B} T (N_{c} + K_{c} + C_{c} - N_{e} - K_{e} - C_{e} - D_{m}) .

We suppose that the stress in the cytoskeleton is regulated at a homeostatic tension such that $𝚺_{c}^{skel} 𝐞_{r} . 𝐞_{r} \overset{def}{=} Σ_{a}$ is a fixed given constant modeling the spontaneous cell contractility. We can then linearize the cell mechanical contributions close to $r_{c} = R_{c}$ to obtain

𝚺_{c}^{skel} 𝐞_{r} . 𝐞_{r} + \frac{2 γ_{c}}{r_{c}} = {\tilde{Σ}}_{a} - k_{c} u,

where ${\tilde{Σ}}_{a} = Σ_{a} + 2 γ_{c} / R_{c}$ and the effective cell mechanical stiffness is $k_{c} = 2 γ_{c} / R_{c}^{2}$ .

Using (Equation 23) and (Equation 29) close to $r_{c, m} = R_{c, m}$ we can express,

- 𝚺_{m}^{el} 𝐞_{r} . 𝐞_{r} = \frac{12 G K_{d} (R_{m}^{3} - R_{c}^{3}) u + (4 G + 3 K_{d}) R_{m}^{3} R_{c} Π}{4 G R_{c}^{4} + 3 K_{d} R_{c} R_{m}^{3}} .

We therefore finally get the linear relation,

{\tilde{Σ}}_{a} + \tilde{k} u + \tilde{α} Π = k_{B} T (N_{c} + K_{c} + C_{c} - N_{e} - K_{e} - C_{e} - D_{m}),

where,

\tilde{k} = - k_{c} + \frac{12 G K_{d} (R_{m}^{3} - R_{c}^{3})}{4 G R_{c}^{4} + 3 K_{d} R_{c} R_{m}^{3}} and \tilde{α} = \frac{(4 G + 3 K_{d}) R_{m}^{3}}{4 G R_{c}^{3} + 3 K_{d} R_{m}^{3}} .

In the limit where $R_{m} ≫ R_{c}$ ,

\tilde{k} = - k_{c} + \frac{4 G}{3 R_{c}} and \tilde{α} = 1 + \frac{4 G}{3 K_{d}} .

Next, using (Equation 8) and (Equation 21) and neglecting $z_{m} ρ_{m} / C_{e} ≪ 1$ we obtain the relation linking the externally controlled osmolarity with the cell and matrix mechanics:

\frac{{\tilde{Σ}}_{a} + \tilde{k} u + (\tilde{α} - 1) Π}{k_{B} T} = N_{e} (e^{- \frac{q (v - v_{N})}{k_{B} T}} - 1) + K_{e} (e^{- \frac{q (v - v_{K})}{k_{B} T}} - 1) + C_{e} (e^{\frac{q v}{k_{B} T}} - 1) - D_{e} .

In a similar way, we combine (Equation 8) with (Equation 9) with again (Equation 21) in the limit where $z_{m} ρ_{m} / C_{e} ≪ 1$ to express the electro-neutrality condition

N_{e} e^{- \frac{q (v - v_{N})}{k_{B} T}} + K_{e} e^{- \frac{q (v - v_{K})}{k_{B} T}} - C_{e} e^{\frac{q v}{k_{B} T}} = \frac{3 z_{c} X_{c}}{4 π R_{c}^{3}} (1 - \frac{3 u}{R_{c}}),

where we have additionally linearized the right handside close to $r_{c} = R_{c}$ .

The two Equations (31) and (32) constitute our final model.

A.4 Cell volume in the reference situation

We begin by computing the cell radius and the cell membrane potential in the reference configuration where by definition $u = 0$ and $Π = D_{e} = 0$ as no Dextran is present at all. In this case, we solve for the membrane potential $v_{c} \overset{def}{=} V_{c}$ and radius $R_{c}$ in (Equation 31) and (Equation 32) to find their reference values. This computation strictly follows Hoppensteadt and Peskin, 2012.

Defining the non-dimensional parameters,

β = \frac{N_{e} e^{q v_{N} / (k_{B} T)} + K_{e} e^{q v_{K} / (k_{B} T)}}{C_{e}} a n d σ = \frac{\tilde{Σ_{a}}}{k_{B} T C_{e}}

we find the reference radius and membrane potential,

R_{c} = {(\frac{3 z_{c} X_{c}}{4 π C_{e} \sqrt{{(σ + 2)}^{2} - 4 β}})}^{1 / 3}

and

V_{c} = \frac{k_{B} T}{2 q} \log (- \sqrt{{(σ + 2)}^{2} - 4 β} + σ + 2) .

Given that the typical concentration of chloride ions outside the cell is of the order of 100 milimolar, the osmotic pressure $k_{B} T C_{e}$ is of the order $10^{5} Pa$ (i.e. an atmosphere). In sharp contrast, the typical mechanical stresses in the cytoskeleton and the cortex are of the order of $10^{2} - 10^{3} Pa$ (Julicher et al., 2007). Therefore, the non-dimensional parameter σ is of the order of $σ \sim 10^{- 3}$ and will be neglected in the following. We then finally obtain the reference values,

R_{c} = {(\frac{3 z_{c} X_{c}}{8 π C_{e} \sqrt{1 - β}})}^{1 / 3}, V_{c} = \frac{k_{B} T}{q} \log (1 - \sqrt{1 - β}) .

The pumping rate enables the cell to maintain a finite a volume. When $j_{p} \to 0$ , $β \to 1$ and the cell swells to infinity because nothing balances the osmotic pressure due to the macromolecules trapped inside. So it is expected that dead cells will swell and lyse. The same happens if the pumping rate is to high. Indeed, as the membrane permeability of potassium is higher than the one of sodium, if the pumping rate is very high, a lot of potassium ions will be brought in (more than sodium ions will be expelled out) and to equilibrate osmolarity with the exterior, water will swell the cell until it bursts. Between these to unphysiological situations, computing the variation of volume with respect to the pumping rate, one gets that this variation vanishes when,

j_{p}^{o p t} = \frac{k_{B} T}{q^{2}} \frac{g_{N} g_{K}}{g_{N} + g_{K}} \log (\frac{N_{e} g_{K}}{K_{e} g_{N}}) .

At such pumping rate, the volume is less sensitive to small variations in the pumping rate that may occur.

Rough estimates

The computation of the effective charge carried by macromolecules is complex. The folding of proteins and the electrostatic screening of charges between them (Manning effect) plays a role. See (Barrat and Joanny, 1997) for a review. We can still make a rough estimate in the following way. We assume that macromolecules are mostly proteins. At physiological $pH = 7.4$ , three types of amino-acids carry a positive charge, Lysine (7%), Arginine (5.3%), Histidine (0.7%) while two others Aspartate (9.9%) and Glutamate (10.8%) carry a negative charge. Added to this, Histidine has a $pKa = 6$ smaller than the pH so the ratio of [histidine neutral base]/[histidine charged acid] is $10^{pH - pKa} = 25$ . Hence the contribution of histidine may be neglected. The occurrence of the aforementioned amino acids in the formation of proteins is also known. The average length of proteins is roughly 400 amino acids. We subsequently obtain the average effective number of negative charges as,

z_{c} = 400 (9.9 + 10.8 - 7 - 5.3) / 100 = 25 .

Such estimate needs to be refined and account for sugars and other macromolecules which carry more negative charges per chain but a interval from $z_{c} = 10$ to $z_{c} = 100$ charges is a plausible estimate.

The estimate of β requires the knowledge of the physiological external concentration of ions $C_{e} = 150 mM$ , $N_{e} = 140 mM$ and $K_{e} = 10 mM$ as well as conductances of sodium and potassium ions through the plasmic membrane. Here, again the situation is complicated since the dynamical opening of channels due to some change in the membrane potential (Hodgkin and Huxley, 1952) as well as the mechanical opening mediated by membrane stretching can play a role and affect these quantities. Nevertheless a rough estimate can be given (Yi et al., 2003)

g_{N} = 2 \times 10^{- 6} C . V^{- 1} . s^{- 1} and g_{K} = 4.5 \times 10^{- 5} C . V^{- 1} . s^{- 1}

Also the pump rate is estimated in Luo and Rudy, 1991,

j_{p} = 2.78 \times 10^{- 12} {mol.s}^{- 1} .

This pump rate is in good agreement with the optimal pump rate predicted by the model,

j_{p}^{o p t} = 3 \times 10^{- 12} {mol.s}^{- 1} .

This leads to an estimate of

β = 0.1 .

The density of macromolecules inside the cell is then found to be $ρ_{c} = 3 \times 10^{6}$ macromolecules per $μ m^{- 3}$ which is a correct order of magnitude (Milo, 2013). To further check the soundness of the above theory we can also compute the membrane potential and obtain $V_{c} = - 73$ mV in good agreement with classical values .

A.5 Osmotic compression of the cell

We now consider the case where, from the reference configuration, we impose an additional osmotic pressure in the external solution with Dextran polymers $Π_{d} = k_{B} T D_{e}$ . We recall that according to formula (Equation 22), $Π = 0$ for small Dextran molecules while $Π = Π_{d}$ for big Dextran molecules.

We use (Equation 31) and (Equation 32) to compute the ensuing small displacement $u$ . Assuming in good approximation that the osmotic pressure imposed by chloride ions is much larger (10⁵ Pa) than the mechanical resistance of the cell cortex and the external matrix (10³ Pa) $k_{B} T C_{e} ≫ \tilde{k} R_{c}$ we find that,

u = - \frac{{(z_{c} X_{c})}^{1 / 3} ((\tilde{α} - 1) Π + Π_{d})}{{4 3}^{2 / 3} π^{1 / 3} {(1 - β)}^{7 / 6} C_{e}^{4 / 3} k_{B} T} .

Strinkingly, making the realistic simplifying assumptions that $K_{d} ≫ G$ and $R_{m} ≫ R_{c}$ , leads to the same displacement of the cell membrane in the two situations of small and big Dextran:

u = - \frac{D_{e} {(z_{c} X_{c})}^{1 / 3}}{{4 3}^{2 / 3} π^{1 / 3} {(1 - β)}^{7 / 6} C_{e}^{4 / 3}},

showing that the two different osmotic loading are not distinguishable at that level. The main text relation (Equation 1) is obtained by assuming that the osmotic pressure of negatively charged ions is half the osmotic pressure of all ions.

However, the mechanical stress applied of the cell is completely different in both situations. For small Dextran, the mechanical stress confining the cell reads,

T^{small} = 𝚺_{m} 𝐞_{r} . {𝐞_{r} |}_{R_{c}} = \frac{2 D_{e} G}{3 C_{e} (1 - β)}

while for big Dextran it reads,

T^{big} = 𝚺_{m} 𝐞_{r} . {𝐞_{r} |}_{R_{c}} = - Π_{d} + T^{small} .

Since $T^{small} ≪ Π_{d}$ by at least one order of magnitude, the most important feature that changes between small and big Dextran is that $T^{small} > 0$ while $T^{big} < 0$ . The physical picture behind this is that small Dextran compresses the cell without draining the water out of the matrix. Therefore, the cell behaves as a small inclusion which volume is reduced by the osmotic compression. In response, the matrix is elastically pulling back to balance the stress at the interface. In contrast, for big Dextran, the water is drained out of the matrix which therefore compresses the cell.

Notice that, like the membrane displacement, the variation of the membrane potential $v \overset{def}{=} v_{c} - V_{c}$ is the same in the two situations:

v = - \frac{k_{B} T}{q} \frac{D_{e}}{2 \sqrt{1 - β} C_{e}},

where we have made the same previous simplifying assumptions that $k_{B} T C_{e} ≫ \tilde{k} R_{c}$ , $K_{d} ≫ G$ and $R_{m} ≫ R_{c}$ . Again such variation is negligibly small in our conditions where $D_{e} ≪ C_{e}$ by several order of magnitudes. This further indicates that the biological response of the cell in response to a big Dextran compression has a mechanical rather than an electro-static origin.

B Mechanical transmission of the stress through Matrigel

The purpose of this appendix is to verify that an extracellular matrix effectively transmits the mechanical stresses to the cells within it. To this end, we embedded soft polyacrylamide (PA) beads in a matrigel scaffold (see Appendix 1—figure 7a), before compressing the whole scaffold either with small or with big dextran. The polyacrylamide beads were fabricated as detailed in Dolega et al., 2017, they have a typical size of 20–50 μm, a bulk modulus modulus $K_{P A} ≃$ 15 kPa and are fluorescently labelled. We evaluated the volume reduction of the PA beads by imaging them before and after compression (Appendix 1—figure 7b). The bead volume was estimated by measuring the surface of its equatorial section, and assuming that the compression is isotropic. This assumption is valid only for beads located at the top of the dropled (z > 1 mm), where the interaction with the substrate is negligible. The bead volume change $Δ V / V$ (Appendix 1—figure 7c) is measured under three experimental conditions:

The culture medium is replaced by fresh medium with identical osmotic pressure (Control, n = 20 beads),
The medium is replaced by fresh medium supplemented with small Dextran at $Π_{d}$ = 5 kPa (Small, n = 20),
The fresh medium is supplemented with big Dextran at $Π_{d}$ = 5 kPa (Big, n = 18),

We observe that, when the MG is compressed by big dextran moelecules, the PA beads are also lose 25% of their volume, even though they are not directly in contact with the osmolytes (green). Such a volume loss is compatible with a mechanical pressure applied of the beads of few kPa (Dolega et al., 2017) and shows that the externally applied osmotic pressure results in a similar mechanical pressure applied on the beads through the drainage and compression of the MG meshwork. In contrast, if the dextran molecules are small enough to penetrate the MG and the PA beads, no mechanical stress is exerted on the inclusions (blue).

Appendix 1—figure 2

Download asset Open asset

Compression of polyacrylamide beads embedded in Matrigel.

(a) Schematic view of the experiment: compressible polyacrylamide beads (in red) are embedded in a matrigel drop and imaged before and after the addition of Dextran at $Π_{d}$ = 5 kPa. (b) Images of a fluorescent bead, respectively before and after compression. The bead volume is deduced from its equatorial section, assuming that the beads remain spherical after compression. (c) Volume decrease of polyacrylamide beads under compression occasioned by small (blue) and big (green) dextran.

C Interstitial space

C.1 Volume fraction of the intestitial space

To evaluate the volume fraction of the interstitial space in multi-cellular spheroids (MCS), we supplement the culture medium with sulforhodamine-B, a hydrophilic fluorophore that stains the extracellular space without penetrating the cells. From confocal sections of MCS (Appendix 1—figure 3a) we determine the thickness of the thin layer between two adjacent cells. By fitting the intensity profile to a Gaussian distribution (Appendix 1—figure 3b), and taking into account that the instrumental function (resolution 270 nm) broadens the profile, we estimate the extracellular layer to 0.9 ± 0.1 μm (histogram in Appendix 1—figure 3c; N = 132). With an average cell diameter of 20 μm, we evaluate that the fraction of extracellular space is approximately $n_{m} = V_{m} / V_{0} = 14 \pm 5 %$ .

Appendix 1—figure 3

Download asset Open asset

Volume fraction estimation.

(a) Confocal section of a MCS, the extracellular space of which is filled with sulforhodamine-B. (b) Intensity profile across two extracellular layers. The width of intercellular space is computed by fitting the intensity to a gaussian profile. (c) Distribution of the intercellular layer thicknesses.

C.2 Compression of ECM in the interstitial space of MCS

The evolution of the rheological properties of ECM filling the interstitial space of MCS is very difficult to evaluate, as the interstitial layer is extremely thin (see previous section). However, we can empirically define an exclusion-size (i.e. porosity), above which globular molecules do not penetrate the gel. To evaluate this exclusion-size, we dip the MCS in a solution containing fluorescent tracers with different radii. As shown in Appendix 1—figure 4a, tracers with $R_{S}$ = 4.4 nm and $R_{S}$ = 5.8 nm permeate the extracellular space of the MCS but not those larger than 14.8 nm. In order to quantify the relative amount of tracers inside the MCS, we compare the average fluorescence measured inside the MCS, $⟨ I_{I n} ⟩$ and in the surrounding solution $⟨ I_{O u t} ⟩$ . Appendix 1—figure 4b reports the relative intensities $⟨ I_{I n} ⟩$ / $⟨ I_{O u t} ⟩$ , obtained respectively at an external osmotic pressure $Π_{d} = 0$ Pa and at $Π_{d} = 5$ kPa. In both cases, the fluorescence level lowers with large tracers. From the results presented in this section, we deduce that:

The intercellular space is rich in fibronectin, a protein constitutive of the extracellular matrix;
The intercellular space has a porosity comparable to that of matrigel gelified in vitro;
The exclusion size of the intercellular space slightly decreases under compression, indicating a moderate compaction of the extracellular matrix.

Appendix 1—figure 4

Download asset Open asset

Exclusion size of the ECM in MCS.

(a) Confocal sections of three MCS dipped in culture media supplemented with Dextran of increasing molecular weights.To avoid saturation of $⟨ I_{O u t} ⟩$ , photomultiplier gain is kept low. This reduces the visibility of extracellular space inside the MCS. In the middle stripe of each image, the brightness is increased of the same amount to make the fluorescence of Dextran visible in the extracellular space. (b) Relative intensity $⟨ I_{I n} ⟩$ / $⟨ I_{O u t} ⟩$ for different Dextran sizes, respectively, without pressure (gray boxes) and with 5 kPa (green boxes). Box sizes and error bars represent respectively the standard error of the mean and the standard deviation. Controls results (no pressure) are obtained from 58 MCS per condition. Experiments under pressure are averaged over 16 (for 4.4 nm), 14 (for 5.8 nm), and 13 (for 14.8 nm) MCS.

D Cell volume

In order to estimate how cells react to an overall compression of the whole MCS, we measure the cell volume change within the aggregate. Cell contours are manually segmented from stack piles obtained with two-photon imaging. Manual segmentation is performed with Amira software. Cell volume is extracted before and after application of osmotic pressure for the same cell in order the measure its compression. Spherical cells are excluded of the analysis as they may undergo cell division and display rapid volume changes. Cells larger than 7000 μm³ are also discarded from the analysis as they may be two cells rather than one. The results are reported in Appendix 1—figure 5 and show that cells appear to be compressed both by small (40 kPa) and big (15 kPa) dextran molecules. The pressures are chosen to match with the experiments presented in Figure 3 of the main article. Notice that this method is much less accurate than the fluorescence exclusion’ method used to determine the volume of individual cells. The results have to be taken as qualitative.

Appendix 1—figure 5

Download asset Open asset

Cell volume change.

Volume loss of cells inside the spheroids, under $Π_{d}^{S m a l l}$ = 40 kPa and $Π_{d}^{B i g}$ = 15 kPa, respectively.

Appendix 1—figure 6

Download asset Open asset

Matrigel bulk modulus.

Compression of the Matrigel beads as a function of the osmotic pressure. Small circles correspond to individual measurements on different MG beads, large circles to the mean value at a given pressure and error bars to the standard deviation. Data are fitted to a hyperelastic model to determine the bulk modulus of matrigel (continuous line).

E Matrigel bulk modulus

In this section, we estimate the bulk modulus of the extracellular matrix. As interstitial ECM is difficult to characterize in-situ, we use matrigel (MG) beads to roughly estimate the rheological properties of ECM. Consistently with native ECM, large Dextran molecules were also excluded from microbeads made of MG suggesting an equivalent effective permeability (Dolega et al., 2021). To determine the bulk modulus of MG beads, we follow their compression at different dextran concentration. To facilitate the measurement, the beads are doped with fluorescent nanoparticles. In Appendix 1—figure 6, we display the volume decrease $V / V_{0}$ of MG beads as a function of the osmotic stresses, between 15 and 500 Pa (V₀ being the bead volume before compression). The continuous line represents the best fit to a Mooney-Rivlin model, the derivative of which represent the bulk modulus $K_{m}$ (Rivlin and Saunders, 1951). For small deformations, the best fit is obtained for a bulk modulus $K_{m} = 450 \pm 100$ Pa.

F Dynamic light scattering and motile activity inside MCS

The scope of this section is to illustrate how to determine cell motility inside an opaque multicellular aggregate. Previous works already indicate that pressure affects cell motility in multicellular spheroids, but the observations are limited to either the superficial layer (Alessandri et al., 2013) or to the long-term centripetal motion (Delarue et al., 2013). Recently, we developed a method to measure the cell velocity in the deep layers of MCS without using confocal microscopy, which is limited in terms of sample thickness and observation time. In our setup (Brunel et al., 2017; Brunel et al., 2020; Appendix 1—figure 7a), the MCS is observed by phase contrast (Appendix 1—figure 7c) and is simultaneously illuminated with an infrared laser (850 nm). The light scattered by the MCS in the forward direction produces an interference pattern, which is collected by a camera (Appendix 1—figure 7b). From the of temporal fluctuations of this pattern (signal shown in Appendix 1—figure 7d and its autocorrelation function in Appendix 1—figure 7e), one computes the average velocity of cells, moving inside the MCS. It has to be noticed that this technique, an evolution of the Dynamic Light Scattering, provides information on the 3D motility, and not only on the 2D motion as previously measured by Alessandri et al., 2013 at the surface of MCS. With this method, we measure the average speed in the three cases of interest: without pressure, when a pressure is selectively applied on the cells, but not on the ECM (small Dextran), and when the pressure is applied to the whole MCS (big Dextran). The results are shown in Appendix 1—figure 7f and g: whereas the average speed is comparable in the first two cases (10 ± 1 μm/hr; magenta and cyan), it is reduced by a factor of two when the compression is exerted on the entire MCS (4.8 ± 0.5 μm/hr, blue).

Appendix 1—figure 7

Download asset Open asset

Motile activity measured by Dynamic Light Scattering (DLS) under different pressure conditions.

(a) The experimental setup combines two counter-propagating optical pathways, with two different wavelengths. This allows us to observe the MCS simultaneously by DLS ( $λ_{D L S}$ = 850 nm, dark red in the sketch) and by phase contrast ( $λ_{p h a s e}$ = 530 nm, green). The DLS signal and phase contrast image are illustrated respectively in panels (b) and (c). DLS signals are acquired at different scattering vectors q and averaged over rings of equal q=|q| (colored sectors in panel (b)). (d) Time evolution of diffraction intensity at three different q; colors correspond to that of sectors in panel (b). (e) Intensity-Intensity autocorrelation functions for the three different scattering vectors. In the single-scattering regime, the intensity signal decorrelates in a typical timescale $τ = 1 / q v^{0}$ , where $v^{0}$ is the mean cell velocity inside the MCS. (**f–g**) respectively the decorrelation time as a function of q and the resulting mean cell velocity, measured in three different conditions: with no pressure, $v^{0}$ = 10 ± 1 μm/hr (magenta), with 5 kPa exerted by small Dextran, $v^{0}$ = 9.8 ± 0.8 μm/hr (cyan), and with $Π_{e} = 5$ kPa exerted by large Dextran, $v^{0}$ = 4.8 ± 0.5 μm/hr (blue). Box sizes correspond to the standard error of the mean (N = 5 MCS) and the error bars to the standard deviation.

G Cytoskeleton and compressibility

Aggregate compressibility depends primarily on three elements: the compressibility of the cells, that of the extracellular matrix, and the volumetric ratio between the two. In this section, we want to evaluate how the rheology of the cytoskeleton contributes to the apparent compressibility of MCS. To do so, we prepare six 96-wells plates of identical spheroids. To avoid the formation of a necrotic core, the MCS initial radius does not exceed 200 µm. Subsequently, five drugs and beg dextran molecules (2 MDa) are added to the different plates:

Y-27632 ROCK inhibitor, to reduce cell contractility,
Blebbistatin, to inhibit acto-myosin activity,
Cytochlasin D, to inhibit actin polymerization,
Nocodazole, to promote microtubule depolymerization,
Paclitaxel, to impede microtubule depolymerization.

As a control, a subset of spheroids are exposed either to the drug alone, without dextran in solution or to Dimethyl Sulfoxide (DMSO). The MCS volumes are obtained by measuring the surface of their equatorial planes and considering them as a perfectly spherical object. The volume is measured before adding the drug, then 45 min after exposure to the drug alone (gray in Appendix 1—figure 8) or to the drug supplemented with dextran (green in Appendix 1—figure 8) and normalized to the initial volume of each MCS.

We observe that drugs modify the MCS volume in different manners (see the figure below) as compared to the control (DMSO). This is in agreement with the fact such pharmacological perturbations are kwown to impact the single-cell volume in different maners (Stewart et al., 2011). When dextran is added to the solution, the spheroids get compressed from their initial state (with the drug); such final compression (Dextran+Drug) is comparable to that obtained with dextran alone but the amout of compression with respect to the intial state varies depending on the drug.

This result is compatible with our idea the 5 kPa Dextran compression reduces almost to the maximum the inter-cellular space and that cells are then almost fully connective in the final state. Thus, depeding on the amount of compression that the drug first creates, the ensuing compression with Dextran will change depending on the available inter-cellular space that remains. For instance, in the presence of cytochalasin, the extra-cellular space is already largely reduced compared to DMSO so when the osmotic compression follows, their is hardly no inter-cellular space which can still be compressed. This seems to be an additional indication that the MCS compressibility in response to a gentle osmotic pressure is more related to the rheology of the extracellular space than to the internal organization and contractility of the cytoskeleton as we argue in Dolega et al., 2021.

Appendix 1—figure 8

Download asset Open asset

Influence of cytoskeleton structure and contractility on MCS compressibility.

In gray, the relative change of the MCS volume after injection of drug (or DMSO alone in control experiments). In green, the vovlume change after addition of drug+dextran at 5 kPa. BoxPlot represent the standard error of the mean, error bars the standard deviation, and circles/squares the volume change of individual spheroiods.

Data availability

Data relating to figures 2, 3, 4, 5 and appendix B are available at https://osf.io/n6ra2/?view_only=059da2ebcd064b75bd12c0c2008b9a6a.

The following data sets were generated

1. Cappello G
(2021) Open Science Framework
ID n6ra2. Extracellular matrix in multicellular aggregates.

https://osf.io/n6ra2/?view_only=059da2ebcd064b75bd12c0c2008b9a6a

References

1. Alessandri K
2. Sarangi BR
3. Gurchenkov VV
4. Sinha B
5. Kießling TR
6. Fetler L
7. Rico F
8. Scheuring S
9. Lamaze C
10. Simon A
11. Geraldo S
12. Vignjevic D
13. Doméjean H
14. Rolland L
15. Funfak A
16. Bibette J
17. Bremond N
18. Nassoy P
(2013) Cellular capsules as a tool for multicellular spheroid production and for investigating the mechanics of tumor progression in vitro
PNAS 110:14843–14848.

https://doi.org/10.1073/pnas.1309482110
- PubMed
- Google Scholar
1. Baker BM
2. Trappmann B
3. Wang WY
4. Sakar MS
5. Kim IL
6. Shenoy VB
7. Burdick JA
8. Chen CS
(2015) Cell-mediated fibre recruitment drives extracellular matrix mechanosensing in engineered fibrillar microenvironments
Nature Materials 14:1262–1268.

https://doi.org/10.1038/nmat4444
- PubMed
- Google Scholar
1. Baranowski B
(1991) Non-equilibrium thermodynamics as applied to membrane transport
Journal of Membrane Science 57:119–159.

https://doi.org/10.1016/S0376-7388(00)80675-4
- Google Scholar
Conference
1. Barrat J-L
2. Joanny J-F
(1997)
Theory of polyelectrolyte solutions

Advances in Chemical Physics, Polymeric Systems.
- Google Scholar
1. Barriga EH
2. Franze K
3. Charras G
4. Mayor R
(2018) Tissue stiffening coordinates morphogenesis by triggering collective cell migration in vivo
Nature 554:523–527.

https://doi.org/10.1038/nature25742
- PubMed
- Google Scholar
(1981) Osmotic deswelling of gels by polymer solutions
Macromolecules 14:719–726.

https://doi.org/10.1021/ma50004a050
- Google Scholar
(2002) The organizing principle: microenvironmental influences in the normal and malignant breast
Differentiation 70:537–546.

https://doi.org/10.1046/j.1432-0436.2002.700907.x
- PubMed
- Google Scholar
1. Bläßle A
2. Soh G
3. Braun T
4. Mörsdorf D
5. Preiß H
6. Jordan BM
7. Müller P
(2018) Quantitative diffusion measurements using the open-source software PyFRAP
Nature Communications 9:1582.

https://doi.org/10.1038/s41467-018-03975-6
- PubMed
- Google Scholar
1. Brochard F
(1981) Polymer networks swollen by a homopolymer solution
Journal De Physique 42:505–511.

https://doi.org/10.1051/jphys:01981004203050500
- Google Scholar
(2018) Mechanotransduction in tumor progression: the dark side of the force
Journal of Cell Biology 217:1571–1587.

https://doi.org/10.1083/jcb.201701039
- PubMed
- Google Scholar
1. Brunel B
2. Blanch C
3. Gourrier A
4. Petrolli V
5. Delon A
6. Joanny J-F
7. Carminati R
8. Pierrat R
9. Cappello G
(2017) Structure and dynamics of multicellular assemblies measured by coherent light scattering
New Journal of Physics 19:073033.

https://doi.org/10.1088/1367-2630/aa7b0f
- Google Scholar
1. Brunel B
2. Levy V
3. Millet A
4. Dolega ME
5. Delon A
6. Pierrat R
7. Cappello G
(2020) Measuring cell displacements in opaque tissues: dynamic light scattering in the multiple scattering regime
Biomedical Optics Express 11:2277–2297.

https://doi.org/10.1364/BOE.388360
- PubMed
- Google Scholar
(2009) A tense situation: forcing tumour progression
Nature Reviews Cancer 9:108–122.

https://doi.org/10.1038/nrc2544
- PubMed
- Google Scholar
Book
(2017) Fluorescence eXclusion Measurement of Volume in Live Cells
Elsevier Ltd.

https://doi.org/10.1016/bs.mcb.2016.11.009
- Google Scholar
1. Cadart C
2. Venkova L
3. Recho P
4. Lagomarsino MC
5. Piel M
(2019) The physics of cell-size regulation across timescales
Nature Physics 15:993–1004.

https://doi.org/10.1038/s41567-019-0629-y
- Google Scholar
(2008) Life and times of a cellular bleb
Biophysical Journal 94:1836–1853.

https://doi.org/10.1529/biophysj.107.113605
- PubMed
- Google Scholar
(1997) Extracellular matrix rigidity causes strengthening of integrin-cytoskeleton linkages
Cell 88:39–48.

https://doi.org/10.1016/S0092-8674(00)81856-5
- PubMed
- Google Scholar
Book
1. Clark AG
2. Paluch E
(2011) Mechanics and regulation of cell shape during the cell cycle
In: Kubiak J, editors. Cell Cycle in Development. Springer. pp. 31–73.

https://doi.org/10.1007/978-3-642-19065-0_3
- Google Scholar
Book
1. Coussy O
(2004) Poromechanics
John Wiley & Sons.

https://doi.org/10.1002/0470092718
- Google Scholar
1. Day RE
2. Kitchen P
3. Owen DS
4. Bland C
5. Marshall L
6. Conner AC
7. Bill RM
8. Conner MT
(2014) Human aquaporins: regulators of transcellular water flow
Biochimica Et Biophysica Acta (BBA) - General Subjects 1840:1492–1506.

https://doi.org/10.1016/j.bbagen.2013.09.033
- PubMed
- Google Scholar
1. Delarue M
2. Montel F
3. Caen O
4. Elgeti J
5. Siaugue JM
6. Vignjevic D
7. Prost J
8. Joanny JF
9. Cappello G
(2013) Mechanical control of cell flow in multicellular spheroids
Physical Review Letters 110:138103.

https://doi.org/10.1103/PhysRevLett.110.138103
- PubMed
- Google Scholar
1. Delarue M
2. Montel F
3. Vignjevic D
4. Prost J
5. Joanny JF
6. Cappello G
(2014) Compressive stress inhibits proliferation in tumor spheroids through a volume limitation
Biophysical Journal 107:1821–1828.

https://doi.org/10.1016/j.bpj.2014.08.031
- PubMed
- Google Scholar
(2010) Control of directed cell migration in vivo by membrane-to-cortex attachment
PLOS Biology 8:e1000544.

https://doi.org/10.1371/journal.pbio.1000544
- PubMed
- Google Scholar
1. Dolega ME
2. Delarue M
3. Ingremeau F
4. Prost J
5. Delon A
6. Cappello G
(2017) Cell-like pressure sensors reveal increase of mechanical stress towards the core of multicellular spheroids under compression
Nature Communications 8:14056.

https://doi.org/10.1038/ncomms14056
- PubMed
- Google Scholar
1. Dolega M
2. Zurlo G
3. Goff ML
4. Greda M
5. Verdier C
6. Joanny J-F
7. Cappello G
8. Recho P
(2021) Mechanical behavior of multi-cellular spheroids under osmotic compression
Journal of the Mechanics and Physics of Solids 147:104205.

https://doi.org/10.1016/j.jmps.2020.104205
- Google Scholar
(2009) Osmotic transport across cell membranes in nondilute solutions: a new nondilute solute transport equation
Biophysical Journal 96:2559–2571.

https://doi.org/10.1016/j.bpj.2008.12.3929
- PubMed
- Google Scholar
1. Engler AJ
2. Sen S
3. Sweeney HL
4. Discher DE
(2006) Matrix elasticity directs stem cell lineage specification
Cell 126:677–689.

https://doi.org/10.1016/j.cell.2006.06.044
- PubMed
- Google Scholar
1. Farutin A
2. Étienne J
3. Misbah C
4. Recho P
(2019) Crawling in a fluid
Physical Review Letters 123:118101.

https://doi.org/10.1103/PhysRevLett.123.118101
- PubMed
- Google Scholar
Book
(2010) Mechanical induction in embryonic development and tumor growth: integrative cues through molecular to multicellular interplay and evolutionary perspectives
In: Serman F, editors. Methods in Cell Biology. Elsevier. pp. 295–321.

https://doi.org/10.1016/j.mod.2016.11.001
- Google Scholar
1. Gerdes J
2. Lemke H
3. Baisch H
4. Wacker HH
5. Schwab U
6. Stein H
(1984)
Cell cycle analysis of a cell proliferation-associated human nuclear antigen defined by the monoclonal antibody ki-67

Journal of Immunology 133:1710–1715.
- PubMed
- Google Scholar
1. Granath KA
(1958) Solution properties of branched dextrans
Journal of Colloid Science 13:308–328.

https://doi.org/10.1016/0095-8522(58)90041-2
- Google Scholar
1. Guo M
2. Pegoraro AF
3. Mao A
4. Zhou EH
5. Arany PR
6. Han Y
7. Burnette DT
8. Jensen MH
9. Kasza KE
10. Moore JR
11. Mackintosh FC
12. Fredberg JJ
13. Mooney DJ
14. Lippincott-Schwartz J
15. Weitz DA
(2017) Cell volume change through water efflux impacts cell stiffness and stem cell fate
PNAS 114:E8618–E8627.

https://doi.org/10.1073/pnas.1705179114
- PubMed
- Google Scholar
1. Han YL
2. Pegoraro AF
3. Li H
4. Li K
5. Yuan Y
6. Xu G
7. Gu Z
8. Sun J
9. Hao Y
10. Gupta SK
11. Li Y
12. Tang W
13. Tang X
14. Teng L
15. Fredberg JJ
16. Guo M
(2020) Cell swelling, softening and invasion in a three-dimensional breast cancer model
Nature Physics 16:101–108.

https://doi.org/10.1038/s41567-019-0680-8
- PubMed
- Google Scholar
(2011) Spontaneous contractility-mediated cortical flow generates cell migration in three-dimensional environments
Biophysical Journal 101:1041–1045.

https://doi.org/10.1016/j.bpj.2011.07.038
- PubMed
- Google Scholar
(1997) Solid stress inhibits the growth of multicellular tumor spheroids
Nature Biotechnology 15:778–783.

https://doi.org/10.1038/nbt0897-778
- PubMed
- Google Scholar
1. Hodgkin AL
2. Huxley AF
(1952) Propagation of electrical signals along giant nerve fibres
Proceedings of the Royal Society of London. Series B, Biological Sciences 140:177–183.

https://doi.org/10.1098/rspb.1952.0054
- PubMed
- Google Scholar
(2009) Physiology of cell volume regulation in vertebrates
Physiological Reviews 89:193–277.

https://doi.org/10.1152/physrev.00037.2007
- PubMed
- Google Scholar
Book
1. Hoppensteadt FC
2. Peskin CS
(2012) Modeling and Simulation in Medicine and the Life Sciences
Springer Science & Business Media.

https://doi.org/10.1007/978-0-387-21571-6
- Google Scholar
1. Hui TH
2. Zhou ZL
3. Qian J
4. Lin Y
5. Ngan AH
6. Gao H
(2014) Volumetric deformation of live cells induced by pressure-activated cross-membrane ion transport
Physical Review Letters 113:118101.

https://doi.org/10.1103/PhysRevLett.113.118101
- PubMed
- Google Scholar
(2014) Mechanotransduction and extracellular matrix homeostasis
Nature Reviews Molecular Cell Biology 15:802–812.

https://doi.org/10.1038/nrm3896
- PubMed
- Google Scholar
1. Isenberg BC
2. Dimilla PA
3. Walker M
4. Kim S
5. Wong JY
(2009) Vascular smooth muscle cell durotaxis depends on substrate stiffness gradient strength
Biophysical Journal 97:1313–1322.

https://doi.org/10.1016/j.bpj.2009.06.021
- PubMed
- Google Scholar
1. Julicher F
2. Kruse K
3. Prost J
4. Joanny J
(2007) Active behavior of the cytoskeleton
Physics Reports 449:3–28.

https://doi.org/10.1016/j.physrep.2007.02.018
- Google Scholar
1. Kedem O
2. Katchalsky A
(1958) Thermodynamic analysis of the permeability of biological membranes to non-electrolytes
Biochimica Et Biophysica Acta 27:229–246.

https://doi.org/10.1016/0006-3002(58)90330-5
- PubMed
- Google Scholar
1. Kedem O
2. Katchalsky A
(1963) Permeability of composite membranes. part 1.—Electric current, volume flow and flow of solute through membranes
Trans. Faraday Soc. 59:1918–1930.

https://doi.org/10.1039/TF9635901918
- Google Scholar
1. Kim DH
2. Li B
3. Si F
4. Phillip JM
5. Wirtz D
6. Sun SX
(2015) Volume regulation and shape bifurcation in the cell nucleus
Journal of Cell Science 128:3375–3385.

https://doi.org/10.1242/jcs.166330
- PubMed
- Google Scholar
Book
1. Kleinman HK
2. Martin GR
(2005) Matrigel: basement membrane matrix with biological activity
In: Martin G. R, editors. Seminars in Cancer Biology. Elsevier. pp. 378–386.

https://doi.org/10.1016/j.semcancer.2005.05.004
- Google Scholar
(2016) Mechanobiology of cell migration in the context of dynamic two-way cell-matrix interactions
Journal of Biomechanics 49:1355–1368.

https://doi.org/10.1016/j.jbiomech.2015.12.023
- PubMed
- Google Scholar
1. Lang F
2. Busch GL
3. Ritter M
4. Völkl H
5. Waldegger S
6. Gulbins E
7. Häussinger D
(1998) Functional significance of cell volume regulatory mechanisms
Physiological Reviews 78:247–306.

https://doi.org/10.1152/physrev.1998.78.1.247
- PubMed
- Google Scholar
Conference
1. Lelièvre SA
2. Bissell MJ
(2006)
Three dimensional cell culture: the importance of microenvironment in regulation of function

Reviews in Cell Biology and Molecular Medicine.
- Google Scholar
1. Levental KR
2. Yu H
3. Kass L
4. Lakins JN
5. Egeblad M
6. Erler JT
7. Fong SF
8. Csiszar K
9. Giaccia A
10. Weninger W
11. Yamauchi M
12. Gasser DL
13. Weaver VM
(2009) Matrix crosslinking forces tumor progression by enhancing integrin signaling
Cell 139:891–906.

https://doi.org/10.1016/j.cell.2009.10.027
- PubMed
- Google Scholar
1. Luo CH
2. Rudy Y
(1991) A model of the ventricular cardiac action potential. depolarization, Repolarization, and their interaction
Circulation Research 68:1501–1526.

https://doi.org/10.1161/01.res.68.6.1501
- PubMed
- Google Scholar
1. Milo R
(2013) What is the total number of protein molecules per cell volume? A call to rethink some published values
BioEssays 35:1050–1055.

https://doi.org/10.1002/bies.201300066
- PubMed
- Google Scholar
1. Monnier S
2. Delarue M
3. Brunel B
4. Dolega ME
5. Delon A
6. Cappello G
(2016) Effect of an osmotic stress on multicellular aggregates
Methods 94:114–119.

https://doi.org/10.1016/j.ymeth.2015.07.009
- PubMed
- Google Scholar
1. Montel F
2. Delarue M
3. Elgeti J
4. Malaquin L
5. Basan M
6. Risler T
7. Cabane B
8. Vignjevic D
9. Prost J
10. Cappello G
11. Joanny JF
(2011) Stress clamp experiments on multicellular tumor spheroids
Physical Review Letters 107:188102.

https://doi.org/10.1103/PhysRevLett.107.188102
- PubMed
- Google Scholar
1. Montel F
2. Delarue M
3. Elgeti J
4. Vignjevic D
5. Cappello G
6. Prost J
(2012) Isotropic stress reduces cell proliferation in tumor spheroids
New Journal of Physics 14:055008.

https://doi.org/10.1088/1367-2630/14/5/055008
- Google Scholar
1. Mori Y
(2012) Mathematical properties of pump-leak models of cell volume control and electrolyte balance
Journal of Mathematical Biology 65:875–918.

https://doi.org/10.1007/s00285-011-0483-8
- PubMed
- Google Scholar
1. Nam S
2. Gupta VK
3. Lee HP
4. Lee JY
5. Wisdom KM
6. Varma S
7. Flaum EM
8. Davis C
9. West RB
10. Chaudhuri O
(2019) Cell cycle progression in confining microenvironments is regulated by a growth-responsive TRPV4-PI3K/Akt-p27^Kip1 signaling Axis
Science Advances 5:eaaw6171.

https://doi.org/10.1126/sciadv.aaw6171
- PubMed
- Google Scholar
1. Nia HT
2. Liu H
3. Seano G
4. Datta M
5. Jones D
6. Rahbari N
7. Incio J
8. Chauhan VP
9. Jung K
10. Martin JD
11. Askoxylakis V
12. Padera TP
13. Fukumura D
14. Boucher Y
15. Hornicek FJ
16. Grodzinsky AJ
17. Baish JW
18. Munn LL
19. Jain RK
(2017) Solid stress and elastic energy as measures of tumour mechanopathology
Nature Biomedical Engineering 1:00004.

https://doi.org/10.1038/s41551-016-0004
- Google Scholar
(2019) Cell mechanosensing is regulated by substrate strain energy rather than stiffness
PNAS 116:22004–22013.

https://doi.org/10.1073/pnas.1904660116
- PubMed
- Google Scholar
1. Paszek MJ
2. Zahir N
3. Johnson KR
4. Lakins JN
5. Rozenberg GI
6. Gefen A
7. Reinhart-King CA
8. Margulies SS
9. Dembo M
10. Boettiger D
11. Hammer DA
12. Weaver VM
(2005) Tensional homeostasis and the malignant phenotype
Cancer Cell 8:241–254.

https://doi.org/10.1016/j.ccr.2005.08.010
- PubMed
- Google Scholar
1. Pelham RJ
2. Wang Y
(1997) Cell locomotion and focal adhesions are regulated by substrate flexibility
PNAS 94:13661–13665.

https://doi.org/10.1073/pnas.94.25.13661
- PubMed
- Google Scholar
1. Rivlin RS
2. Saunders D
(1951) Large elastic deformations of isotropic materials vii. experiments on the deformation of rubber
Philosophical Transactions of the Royal Society of London Series A, Mathematical and Physical Sciences 243:251–288.

https://doi.org/10.1098/rsta.1951.0004
- Google Scholar
1. Roycroft A
2. Mayor R
(2016) Molecular basis of contact inhibition of locomotion
Cellular and Molecular Life Sciences 73:1119–1130.

https://doi.org/10.1007/s00018-015-2090-0
- PubMed
- Google Scholar
(2012) Actin cortex mechanics and cellular morphogenesis
Trends in Cell Biology 22:536–545.

https://doi.org/10.1016/j.tcb.2012.07.001
- PubMed
- Google Scholar
1. Sedzinski J
2. Biro M
3. Oswald A
4. Tinevez JY
5. Salbreux G
6. Paluch E
(2011) Polar actomyosin contractility destabilizes the position of the cytokinetic furrow
Nature 476:462–466.

https://doi.org/10.1038/nature10286
- PubMed
- Google Scholar
1. Segel M
2. Neumann B
3. Hill MFE
4. Weber IP
5. Viscomi C
6. Zhao C
7. Young A
8. Agley CC
9. Thompson AJ
10. Gonzalez GA
11. Sharma A
12. Holmqvist S
13. Rowitch DH
14. Franze K
15. Franklin RJM
16. Chalut KJ
(2019) Niche stiffness underlies the ageing of central nervous system progenitor cells
Nature 573:130–134.

https://doi.org/10.1038/s41586-019-1484-9
- PubMed
- Google Scholar
1. Sopher RS
2. Tokash H
3. Natan S
4. Sharabi M
5. Shelah O
6. Tchaicheeyan O
7. Lesman A
(2018) Nonlinear elasticity of the ECM fibers facilitates efficient intercellular communication
Biophysical Journal 115:1357–1370.

https://doi.org/10.1016/j.bpj.2018.07.036
- PubMed
- Google Scholar
1. Staunton JR
2. So WY
3. Paul CD
4. Tanner K
(2019) High-frequency microrheology in 3D reveals mismatch between cytoskeletal and extracellular matrix mechanics
PNAS 116:14448–14455.

https://doi.org/10.1073/pnas.1814271116
- PubMed
- Google Scholar
1. Staverman AJ
(1952) Non-equilibrium thermodynamics of membrane processes
Transactions of the Faraday Society 48:176–185.

https://doi.org/10.1039/tf9524800176
- Google Scholar
(2011) Hydrostatic pressure and the actomyosin cortex drive mitotic cell rounding
Nature 469:226–230.

https://doi.org/10.1038/nature09642
- Google Scholar
Book
1. Strange K
(1993) Cellular and Molecular Physiology of Cell Volume Regulation
CRC Press.

https://doi.org/10.1201/9780367812140
- Google Scholar
(2016) Collective cell durotaxis emerges from long-range intercellular force transmission
Science 353:1157–1161.

https://doi.org/10.1126/science.aaf7119
- PubMed
- Google Scholar
(2012) Coherent angular motion in the establishment of multicellular architecture of glandular tissues
PNAS 109:1973–1978.

https://doi.org/10.1073/pnas.1119578109
- PubMed
- Google Scholar
1. Taubenberger AV
2. Girardo S
3. Träber N
4. Fischer-Friedrich E
5. Kräter M
6. Wagner K
7. Kurth T
8. Richter I
9. Haller B
10. Binner M
11. Hahn D
12. Freudenberg U
13. Werner C
14. Guck J
(2019) 3d microenvironment stiffness regulates tumor spheroid growth and mechanics via p21 and ROCK
Advanced Biosystems 3:e1900128.

https://doi.org/10.1002/adbi.201900128
- PubMed
- Google Scholar
1. Therien AG
2. Blostein R
(2000) Mechanisms of sodium pump regulation
American Journal of Physiology-Cell Physiology 279:C541–C566.

https://doi.org/10.1152/ajpcell.2000.279.3.C541
- PubMed
- Google Scholar
1. Tinevez JY
2. Schulze U
3. Salbreux G
4. Roensch J
5. Joanny JF
6. Paluch E
(2009) Role of cortical tension in bleb growth
PNAS 106:18581–18586.

https://doi.org/10.1073/pnas.0903353106
- PubMed
- Google Scholar
1. Turlier H
2. Audoly B
3. Prost J
4. Joanny JF
(2014) Furrow constriction in animal cell cytokinesis
Biophysical Journal 106:114–123.

https://doi.org/10.1016/j.bpj.2013.11.014
- PubMed
- Google Scholar
1. Vogel V
(2018) Unraveling the mechanobiology of extracellular matrix
Annual Review of Physiology 80:353–387.

https://doi.org/10.1146/annurev-physiol-021317-121312
- PubMed
- Google Scholar
(2016) Hyaluronan-Derived swelling of solid tumors, the contribution of collagen and Cancer cells, and implications for Cancer therapy
Neoplasia 18:732–741.

https://doi.org/10.1016/j.neo.2016.10.001
- PubMed
- Google Scholar
1. Yi CS
2. Fogelson AL
3. Keener JP
4. Peskin CS
(2003) A mathematical study of volume shifts and ionic concentration changes during ischemia and hypoxia
Journal of Theoretical Biology 220:83–106.

https://doi.org/10.1006/jtbi.2003.3154
- PubMed
- Google Scholar
1. Zhou EH
2. Trepat X
3. Park CY
4. Lenormand G
5. Oliver MN
6. Mijailovich SM
7. Hardin C
8. Weitz DA
9. Butler JP
10. Fredberg JJ
(2009) Universal behavior of the osmotically compressed cell and its analogy to the colloidal glass transition
PNAS 106:10632–10637.

https://doi.org/10.1073/pnas.0901462106
- PubMed
- Google Scholar
(2015) Optical volume and mass measurements show that mammalian cells swell during mitosis
Journal of Cell Biology 211:765–774.

https://doi.org/10.1083/jcb.201505056
- PubMed
- Google Scholar

Article and author information

Author details

Monika E Dolega

Université Grenoble Alpes, Laboratoire Interdisciplinaire de Physique, CNRS, Grenoble, France

Contribution
Conceptualization, Data curation, Validation, Investigation, Writing - original draft

Competing interests
No competing interests declared
Sylvain Monnier

Université de Lyon, Université Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, VILLEURBANNE, France

Contribution
Data curation, Investigation, Writing - review and editing

Competing interests
No competing interests declared
Benjamin Brunel

Université Grenoble Alpes, Laboratoire Interdisciplinaire de Physique, CNRS, Grenoble, France

Contribution
Investigation

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0003-2858-5074
Jean-François Joanny

Collège de France, PSL Research University, Paris, France

Contribution
Formal analysis, Writing - review and editing

Competing interests
No competing interests declared
Pierre Recho

Université Grenoble Alpes, Laboratoire Interdisciplinaire de Physique, CNRS, Grenoble, France

Contribution
Conceptualization, Formal analysis, Methodology, Writing - original draft

For correspondence
pierre.recho@univ-grenoble-alpes.fr

Competing interests
No competing interests declared
Giovanni Cappello

Université Grenoble Alpes, Laboratoire Interdisciplinaire de Physique, CNRS, Grenoble, France

Contribution
Conceptualization, Data curation, Formal analysis, Funding acquisition, Validation, Investigation, Methodology, Writing - original draft, Project administration

For correspondence
giovanni.cappello@univ-grenoble-alpes.fr

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-5012-367X

Funding

Agence Nationale de la Recherche (ANR-13-BSV5-0008-01)

Giovanni Cappello

Centre National de la Recherche Scientifique (MechanoBio 2018)

Giovanni Cappello

Ligue Contre le Cancer

Sylvain Monnier

Institut National de la Santé et de la Recherche Médicale (PC201407)

Giovanni Cappello

Centre National de la Recherche Scientifique

Pierre Recho

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We warmly thank J Prost and F Jülicher for drawing our attention to the potential impact of the poroelasticity in MCS, A Dawid and J Revilloud for the valuable suggestion to set up the proliferation assay in MG, C Verdier for the valuable exchanges about the evolution of the ECM rheology under stress, and T Boudou for the many fruitful discussion and careful proofreading of the manuscript. This work was supported by the Agence Nationale pour la Recherche (Grant ANR-13-BSV5-0008-01), by the Institut National de la Santé et de la Recherche Médicale (Grant PC201407), by the Centre National de la Recherche Scientifique (Grant MechanoBio 2018), by the Comité de Haute-Savoie de la Ligue contre le Cancer, and by a CNRS Momentum grant (PR).

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.