A mechano-osmotic feedback couples cell volume to the rate of cell deformation

  1. Larisa Venkova
  2. Amit Singh Vishen
  3. Sergio Lembo
  4. Nishit Srivastava
  5. Baptiste Duchamp
  6. Artur Ruppel
  7. Alice Williart
  8. Stéphane Vassilopoulos
  9. Alexandre Deslys
  10. Juan Manuel Garcia Arcos
  11. Alba Diz-Muñoz
  12. Martial Balland
  13. Jean-François Joanny
  14. Damien Cuvelier
  15. Pierre Sens  Is a corresponding author
  16. Matthieu Piel  Is a corresponding author
  1. Institut Curie, PSL Research University, CNRS, UMR 144, France
  2. Institut Pierre Gilles de Gennes, PSL Research University, France
  3. Institut Curie, PSL Research University, CNRS, UMR 168, France
  4. Cell Biology and Biophysics Unit, European Molecular Biology Laboratory, Germany
  5. Laboratoire Interdisciplinaire de Physique, France
  6. Sorbonne Université, Institut National de la Santé et de la Recherche Médicale, Institut de Myologie, Centre de Recherche en Myologie, France
  7. Sorbonne Université, INSERM, France

Abstract

Mechanics has been a central focus of physical biology in the past decade. In comparison, how cells manage their size is less understood. Here, we show that a parameter central to both the physics and the physiology of the cell, its volume, depends on a mechano-osmotic coupling. We found that cells change their volume depending on the rate at which they change shape, when they spontaneously spread or when they are externally deformed. Cells undergo slow deformation at constant volume, while fast deformation leads to volume loss. We propose a mechanosensitive pump and leak model to explain this phenomenon. Our model and experiments suggest that volume modulation depends on the state of the actin cortex and the coupling of ion fluxes to membrane tension. This mechano-osmotic coupling defines a membrane tension homeostasis module constantly at work in cells, causing volume fluctuations associated with fast cell shape changes, with potential consequences on cellular physiology.

Editor's evaluation

The paper by Venkova et al. is a comprehensive study of mammalian cell volume dynamics during the common cellular process of adhesion and spreading on a flat substrate, osmotic changes, and mechanical confinement. The paper reveals a complex interplay between cell water/ion regulation, cytoskeletal processes, and mechanical deformation of the cell. The topic is important in cell physiology and should be of considerable interest to cell biologists, mechanobiologists and biophysicists.

https://doi.org/10.7554/eLife.72381.sa0

Introduction

In recent years, in vivo imaging has revealed that, in a variety of physiological and pathological contexts, cells undergo large deformations (Weigelin et al., 2014), sometimes being squeezed to a tenth of their resting diameter. Migrating cells, in particular fast-moving immune or cancer cells, can deform to a large extent in only a few minutes (Thiam et al., 2016; Vargas et al., 2017; Beunk et al., 2019), for example, when they cross an endothelial barrier (Subramanian et al., 2020). Even faster deformations, below the second timescale, can be observed in circulating cells pushed through small blood capillaries. Altogether, these examples show that large cell deformations are physiological and occur across a large range of timescales. Large cell shape changes must involve significant changes in volume, surface area, or both. But the number of studies on cell volume modulation upon cell deformation is still very small (Guilak, 1995; Liu et al., 2020). It is still not clear whether the material that cells are made of is rather poroelastic (Moeendarbary et al., 2013), losing volume when pressed, or behaves like a liquid droplet, extending its surface area at constant volume. Two articles, measuring volume using 3D reconstruction from confocal slices, report that cells that are more spread are smaller in volume (Guo et al., 2017; Xie et al., 2018), leading to a higher density and potential long-term effects on cell fate (Guo et al., 2017). On the other hand, another article, using volume measurements by fluorescence exclusion (FXm), reports no or slightly positive correlation between spreading area and cell volume (Perez Gonzalez et al., 2018), reflecting the fact that as cells grow larger, their spreading area also increases.

Different models have recently been proposed to explain a coupling between cells shape changes and cell volume modulation (Guo et al., 2017; Xie et al., 2018; Adar and Safran, 2020; McEvoy et al., 2020). Most of them are based on the same type of scenario: depending on the timescale and extent of the deformation, cell shape changes can stress the cell surface, including the membrane and the actin cortex (Chugh and Paluch, 2018). This stress can be relaxed due to cortex turnover, unfolding of membrane reservoirs (Pietuch and Janshoff, 2013), and detachment of the membrane from the cortex with the formation of blebs (Tinevez et al., 2009). Stress in these structures can also lead to the modulation of ion fluxes (Jiang and Sun, 2013) resulting in cell volume changes. Despite its broad relevance for cell mechanics and cell physiology, the consequences of this type of scenario have not been explored in depth experimentally.

Using FXm to accurately measure volume in live cells (Zlotek-Zlotkiewicz et al., 2015; Cadart et al., 2017; Cadart et al., 2018), we found that when cells deform as they spread on a flat adhesive surface, the degree of volume changes depends on the speed of spreading. To explain this observation, we propose an extension of the classical pump and leak model (PLM; Cadart et al., 2019), including a mechano-osmotic coupling activated upon cell deformation occurring faster than the membrane tension/actin cortex relaxation timescale. We further probe our model assumptions and predictions experimentally by characterizing the cell volume response during osmotic shocks and during ultra-fast (ms) mechanical cell deformation, as well as by performing tether pulling experiments on spreading cells to assess their membrane tension. We believe that our observations, together with this novel physical model, constitute strong evidence for the existence of a mechano-osmotic coupling constantly at work in animal cells and modulating their volume as they deform.

Results

Cell volume is not correlated to the final steady-state spreading area but significantly decreases during cell spreading

We first asked whether in a population growing and dividing at steady state, cells display a correlation between their spreading area and their volume. We used HeLa cells expressing hGeminin-mCherry, which accumulates in the nucleus during the S phase. Cell spreading area was measured using phase contrast and cell volume using FXm (Cadart et al., 2017, Figure 1A images). We did not find any strong correlation between spreading area and volume for HeLa hGeminin-mCherry, as well as for HeLa EMBL (Kyoto) and RPE-1 cells measuring their spreading area with reflection interference contrast microscopy (RICM) (Rädler and Sackmann, 1993; Cuvelier et al., 2007), larger cells in volume being also slightly more spread (Figure 1A graph and Figure 1—figure supplement 1A and B). A clearer positive correlation was observed for 3T3 fibroblasts, which were also generally more spread for a given volume (Figure 1—figure supplement 1C). Using the hGeminin cell cycle marker, we observed that S/G2 cells tend to be larger and more spread than G1 cells (Figure 1A graph), suggesting that the positive correlation is simply due to cell growth, with cells increasing their spreading area as they grow. Using live-cell recording of phase, volume, and hGeminin, we also considered cells at given windows of time following cell division, to examine the correlation between volume and surface area at a given cell cycle stage and thus independently of cell growth. Considering the same group of cells at various times after mitosis, or after the G1/S transition, we could not observe any correlation between cell volume and spreading area at any given cell cycle stage (Figure 1A). Finally, to extend the range of spreading areas considered, we used adhesive micro-patterns with areas smaller than the average spontaneous steady-state spreading area of HeLa EMBL (Kyoto) cells (Figure 1B images). We found that the distribution of volumes did not change when cells were plated on smaller adhesive patterns (Figure 1B graph). Overall, these experiments suggest that, as reported before (Perez Gonzalez et al., 2018), there is no strong correlation, at the cell population level, between spreading area and cell volume, independently of the cell cycle stage.

Figure 1 with 3 supplements see all
Cell volume in spreading cells.

(A) Left: Composite of FXm in GFP channel and fluorescent image in mCherry channel of HeLa hgem-mCherry cells. Scale bar: 100 µm. Right: Relation between volume and spreading area of HeLa hgem-mCherry cells at the different cell cycle stages, N=3: M+1h (n=131) correlation coefficient R=0.11, M+4h (n=131) R=0.23, M+6h (n=131) R=0.26, G1/S (n=99) R=0.20, G1/S+4h (n=92) R=0.22. Error bars represent standard deviation. (B) Left: Typical images of micropatterns and typical images of cells plated on micropatterns. Scale bar: 10 µm. Right: Average volume of HeLa Kyoto cells plated on the patterns (measurements are done 4 hr after cell plating) of different shape and size in comparison with non-patterned cells. Blue: rectangle 30×13 µm2 (n=131, N=2); red: rectangle 40×7 µm2 (n=214, N=2); purple: circle, r=15 µm (n=338, N=4); green: circle, r=10 µm (n=242, N=3); gray: non-patterned cells (n=286, N=3). Error bars represent standard deviation. There is nos statistical difference between patterned cells and non-patterned cells: rectangle 30×13 p=0.15, rectangle 40×7 µm2 p=0.96, r=15 µm p=0.63, r=10 µm p=0.97. (C) Top: Side view of a HeLa-Lifeact (black) cell spreading on fibronectin-coated glass. Scale bar: 20 µm. Bottom: Scheme of shape transition during cell spreading. (D) Left: FXm and RICM imaging of a HeLa Kyoto cell spreading on fibronectin-coated glass. Scale bar: 20 µm. Right: Volume (red) and spreading area (blue) of cell represented on the left panel. (E) Average normalized volume of control HeLa Kyoto cells (blue, n=127, N=3) spreading on fibronectin-coated glass, or plated on PLL-PEG-coated glass (cyan, n=493, N=5). Error bars represent standard deviation. (F) Average spreading area of control HeLa Kyoto cells (n=125, N=3), spreading on fibronectin-coated glass. Error bars represent standard deviation. (G) Linear derivatives dA/dt (blue) and dV/dt (red) for average spreading area and volume represented on (F) and (E) for sliding window 10 min. (H) Average normalized volume (left) and spreading area (right) of control HeLa cells divided in three categories based on their initial spreading speed, N=3: slow (n=42), moderate (n=43), fast (n=42). Error bars represent standard deviation. (I) Volume flux (dV/dt) of single control HeLa Kyoto cells (n=194, N=3) plotted versus their spreading speed (dA/dt) at the first 10 min of spreading. The data are fitted with linear regression y=−0.31x−0.71, R2=0.19. Error bars represent standard deviation. Color code indicate three groups of cells represented on (H). (J) Median volume flux (dV/dt) of HeLa Kyoto cells plotted versus median spreading speed (dA/dt) at the different time intervals of spreading (n=194, N=3). Error bars represent standard deviation.

Figure 1—source data 1

Data tables related to quantifications in Figure 1.

https://cdn.elifesciences.org/articles/72381/elife-72381-fig1-data1-v2.zip

Previous studies also reported volume loss during cell spreading (Guo et al., 2017). When plated on a fibronectin-coated substrate, HeLa EMBL cells showed a transition from a sphere to a half-sphere in about 15 min, then continued spreading by extending lamellipodial protrusions (Figure 1C and Figure 1—video 1). We recorded spreading cells, combining FXm to measure volume and RICM to measure spreading area accurately (Figure 1D and Figure 1—video 2). RICM images showed an initial spreading phase of about 15±10 min until the radius of the contact region equaled that of the cell, corresponding to a hemispheric cap cell shape, which was followed by an extension of lamellipodial protrusions. Cell spreading was accompanied by a small (5% on average) but significant loss of volume, typically occurring during the first 20 min of spreading and followed by a volume increase of about 5%/hr, in the range of the expected cell growth for a doubling time of about 20 hr (Figure 1E, F and G). The same was observed for cells that had been synchronized by serum starvation, but with a smaller standard deviation (Figure 1—figure supplement 1D). The absolute precision of our measurements is within 10%, evaluated by comparing the initial average volume of cell populations in different chambers and on different days (Figure 1—figure supplement 1E). But the accuracy of the measurement when following the same individual cell is rather 1% (Figure 1—figure supplement 1F, evaluated by measuring the volume of the same cell multiple times at a 30 ms frame rate). Combining quantitative phase and volume measurement, we found that only cell volume decreased while dry mass remained constant over the few tens of minutes of initial cell spreading (Figure 1—figure supplement 1G), causing a transient density increase (Figure 1—figure supplement 1H and G). This suggests a loss of water (and probably small osmolites like ions) from the cell, similar to volume regulatory decrease following a hypo-osmotic shock (Pedersen et al., 2013). Cells plated on PLL-PEG, instead of fibronectin, did not spread and displayed an increase in the volume of about 7%/hr (Figure 1E). This result, together with our observation on steady-state spread cells, suggests that the spreading dynamics rather than the final spreading area might be coupled to the loss of volume.

The loss of volume during spreading depends on the spreading speed

Taking advantage of the intrinsic variability in the cell spreading dynamics, we considered single cell volume and spreading trajectories. We observed that individual spreading cells could display a large range of volume loss (Figure 1—figure supplement 1I). Pooling cells together according to their spreading speed, we observed that faster-spreading cells were losing more volume, whereas slow-spreading cells lost less volume or did not lose volume at all (Figure 1H). To further validate this correlation, we measured the initial spreading speed for the first 10 min of spreading and plotted it against the rate of volume loss, for individual cells (Figure 1I). The graph clearly shows that faster spreading cells also lose volume faster in this initial spreading phase. Spreading speed and volume loss are both slowing down with time (Figure 1J), whereas absolute spreading area increases (Figure 1—figure supplement 1J). Conversely, the amount of volume lost was not correlated with the initial cell volume (Figure 1—figure supplement 1K). Overall, these data show that volume loss in spreading cells is a transient phenomenon correlated with the spreading kinetics and not the absolute spreading area.

Early spreading dynamics were shown to strongly depend on the properties of the actomyosin cortex (Cuvelier et al., 2007). Hence, we affected F-actin with a low dose (100 nM) of Latrunculin A (Lat A) which still allowed cell spreading, and myosin with the 100 µM ROCK inhibitor Y-27632 (Y-27, Figure 2A). As expected, we found that Lat A-treated cells spread slower, while Y-27-treated cells spread faster than control cells (Figure 2B). Accordingly, Lat A treated cells lost less volume (2–3%) than control cells, while Y-27 treated cells lost more (10–15%, Figure 2C and D, Figure 2—figure supplement 1A). Y-27 treated cells plated on PLL-PEG substrate on which they could not spread, increased their volume like control cells (Figure 2—figure supplement 1B), thereby showing that larger volume loss was not due to the drug treatment itself but was a result of the spreading kinetics in the presence of the drug. This coupling between spreading speed and volume loss was also found to be very similar for other cell types, RPE-1 and 3T3 (Figure 2—figure supplement 1C-I), although RPE-1 cells displayed an initial phase of volume increase before eventually losing volume (Figure 2—figure supplement 1E, F, H, I), a phenomenon that we have not investigated further in this article. This initial phase of volume increase was lost upon Y-27 treatment (Figure 2—figure supplement 1F, Figure 2—figure supplement 1I), suggesting that it was due to induction of contractility through mechanotransduction pathways (Burridge et al., 2019). In order to increase spreading speed without changing cell contractility, we added 1 mM MnCl2 during the spreading experiment. This treatment increased spreading speed, as previously reported (Edwards et al., 1988) and also increased the volume loss (Figure 2—figure supplement 1J and Figure 2D). These data together suggest a general effect of spreading speed on volume modulation, with a loss of volume reaching up to 20% for fast-spreading cells.

Figure 2 with 1 supplement see all
Cell volume depends on spreading speed.

(A) RICM imaging of control HeLa Kyoto cell or cell treated with 100 nM Latrunculin A or with 100 µM Y-27632 spreading on fibronectin-coated glass. Scale bar: 20 µm. (B) Average spreading area of control HeLa Kyoto cells (gray, n=125, N=3), cells treated with 100 nM Latrunculin A (orange, n=30, N=2) or 100 µM Y-27632 (red, n=98, N=3). Error bars represent standard error. (C) Average normalized volume of control HeLa Kyoto cells (gray, n=125, N=3), cells treated with 100 nM Latrunculin A (orange, n=30, N=2) or 100 µM Y-27632 (red, n=98, N=3). Error bars represent standard error. (D) Median volume flux (dV/dt) of control (blue, n=194, N=3), 100 µM Y-27632 (red, n=121, N=4), 100 nM Latrunculin A (orange, n=41, N=3) or 1 mM MnCl2 (N=3, n=57) treated HeLa Kyoto cells plotted versus their spreading speed (dA/dt) at the first 10 min of spreading. Average dV/dt(dA/dt) for three groups of control cells from Figure 1I are shown on the graph. The dashed line is a linear regression for control cells from panel 1I. Error bars represent standard error.

Figure 2—source data 1

Data tables related to quantifications in Figure 2.

https://cdn.elifesciences.org/articles/72381/elife-72381-fig2-data1-v2.zip

The classical pump and leak model describes properly the osmotic response of non-spreading cells but cannot explain the water loss during spreading

Because the observation of a coupling between volume modulation and spreading speed was not reported before, we asked whether classical volume modulation models are sufficient to explain it. Our results show that during fast spreading, cells lose more than 10% of their water content. Water loss exceeding 1% is considered to be dominated by osmotic volume regulation (see Appendix 1 and Cadart et al., 2019; Yellin et al., 2018). Volume set point and large volume modulation such as volume regulatory response following osmotic shocks can be accounted for by the general theoretical framework of the PLM (see Appendix 1 and Cadart et al., 2019; Tosteson and Hoffman, 1960; Kay, 2017). Briefly, the cell volume is determined by an osmotic balance involving the active pumping of specific ions (mainly sodium and potassium) to compensate for the pressure from impermeant solutes in the cell (Figure 3A). The PLM has been verified experimentally on several occasions, mostly with indirect methods for cell volume measurements Kay and Blaustein, 2019. We thus decided to check that we could reproduce these results with our cells. We performed series of osmotic shock experiments using PEG400 or distilled water, while recording cell volume by FXm (Figure 3B and Figure 3—figure supplement 1A, B and Figure 3—video 1). Cells showed the expected response to both hypo and hyperosmotic shocks, with a fast change in volume (less than a minute timescale) followed by a slower adaptation (timescale of minutes) (Figure 3C and Figure 3—figure supplement 1C). We also checked, using quantitative phase measure of dry mass, that these fast changes in volume were not accompanied by any change in dry mass and thus corresponded to water (and ion) fluxes, as expected (Figure 3—figure supplement 1D). Because of timescale separation between water flux in the seconds timescales and active ion transport, which takes minutes, upon an osmotic shock, cells first display a passive response corresponding to water fluxes, followed by a slower response due to ion exchanges. The Ponder relation (Ponder and Saslow, 1931), which relates the relative change in cell volume right after the shock (at timescale of seconds), to the relative difference of osmotic pressure imposed experimentally, corresponds to the passive cell response. Ponder’s plot showed a very good agreement with previous reports (Pritchard and Guilak, 2004; Zhou et al., 2009), with a linear relation between the change in volume and the change in osmotic pressure, over a large range of imposed external osmolality (Figure 3D) and corresponds to about 30% of osmotically inactive volume (volume occupied by large molecules or solid components). As shown by others (Pritchard and Guilak, 2004), we find that the Ponder relation does not depend on the integrity of the actin cytoskeleton, as cells treated with 2 µM Lat A show the same relation. Ponder relation does not depend either on the inhibition of sodium/proton exchanger NHE1 by 50 µm EIPA (Figure 3D). These experiments also allowed us to estimate the bulk modulus of the cells defined as B=VisoΔPΔV, where Viso is the volume in isosmotic state, ΔV is the volume change, induced by the osmotic pressure difference ΔP (order of GPa, Figure 3—figure supplement 1E), which is in good agreement with previous measurements (Guo et al., 2017; Monnier et al., 2016). These results show both that our cell volume measurements are accurate, even for small volume changes, and that our cells display the expected response to osmotic shocks, explained by the classical PLM, in agreement with previously published results (Ponder and Saslow, 1931; Pritchard and Guilak, 2004; Zhou et al., 2009; Roffay et al., 2021).

Figure 3 with 2 supplements see all
Verifying the Pump and Leak Model and its mecano-sensitive extension.

(A) Schematic of “pump-leak” model (PLM). In brief, the plasma membrane let ions and water pass, with specific channels which increase their permeation coefficient (the ‘leak’). Ions can also be pumped out of the cell (the ‘pump’). Outside the cell, the concentration of ions is about 300 mM, while it is only of about 150 mM inside the cell. The quantity of ion species also differs, with more anions outside the cell (because proteins are on average negatively charged), more sodium outside the cell, and more potassium inside. To achieve the equilibrium, sodium and chloride need to be constantly pumped outside of the cell. The cell also contains nonpermeant osmolytes: charged osmolytes (‘zY’ at the scheme) and neutral osmolytes (‘Y’ at the scheme), for example, proteins, amino-acids, and sugars (which cannot pass the membrane at the same rate as water or ions). Considering these equilibria and the resulting osmotic balance, the PLM predicts the volume of the cell (see the model in Appendix 1 and Cadart et al., 2019). Numbers at the panels are taken from the model in Appendix 1. (B) Left: FXm images of HeLa Kyoto cells exposed to media exchange of same osmolarity, hypertonic and hypotonic. Right: Volume of the cells represented on the left panel. (C) Average normalized HeLa Kyoto cells volume response to osmotic shocks of different magnitudes acquired every 30 s. Number of cells in the experiments: control Piso/P=1 (n=51, N=1), Piso/P=1.11 (n=30, N=1), Piso/P=2 (n=17, N=1), Piso/P=0.82 (n=33, N=1), Piso/P=0.6 (n=67, N=1). (D) Ponder’s relation for control HeLa Kyoto cells (blue), treated with 2 µM Lat A (orange), or 50 µM EIPA (green). Each point represents average value of single experiment. Average number of cells in each experiment n~58. Error bars represent standard deviation. Dashed line represents linear regression fit for control cells y=0.67x+0.33, R2=0.98. Coefficient 0.67 refers to the ratio of osmotically active volume to total volume named ‘R’ in the Appendix 1. (E) Tether force measurements of control HeLa Kyoto cells (blue) and treated with 100 µM Y-27632 (red) at the ‘spreading phase’ (measurements are performed within 30–90 min after cell plating), ‘steady state, spread’ (measurements are performed within 4–5 hr after plating) or ‘steady state, non-spread’ (measurements are performed within 4–5 hr after plating on 20 µm diameter micropatterns). For ‘control, spreading phase’ n=50, N=6; for ‘Y-27, spreading phase’ n=27, N=3; for ‘control, steady state, spread’ n=55, N=9; for ‘Y-27, steady state’ n=21, N=3, for ‘control, steady state, non-spread’ n=18, N=3, for ‘Y-27, steady state, non-spread’ n=20, N=3. Error bars represent standard deviation. The results of statistical tests are shown at the graph.

Figure 3—source data 1

Data tables related to quantifications in Figure 3.

https://cdn.elifesciences.org/articles/72381/elife-72381-fig3-data1-v2.zip

The classical PLM does not account for the cell shape changes and mechanics. Several additional mechanisms have been proposed to account for the coupling between cell shape and cell volume. A recent model proposed a direct extension of the PLM to account for cell spreading, by including the assumption that channels and pumps are working differently on the adhered and the free surface of the cell (Adar and Safran, 2020). Nevertheless, such a model does not predict an effect of the spreading speed on volume, but rather an effect of the spreading area itself, while our data suggest that the opposite is true in our experiments. The correlation of volume loss with spreading speed suggests that the effect on volume could be due to a change in the mechanical state of the cell surface (membrane and/or cortex). A contribution of contractility to volume regulation has been proposed before (Tao and Sun, 2015). To test whether cell contractility could directly affect the cell volume, we recorded the volume of the same non-spreading cells, plated on PLL, before and after treatment with the contractility inhibitor (Rho-associated protein kinase ROCK inhibitor) Y-27632 (Figure 3—figure supplement 1F). We did not observe any significant difference between a change of medium with a control medium and with a medium containing the drug. This shows that Y-27632 has an effect on volume only during cell spreading. To confirm that Y-27632 treatment during cell spreading reduces traction forces (and thus cortical tension), we also performed traction force microscopy (Sabass et al., 2008) during cells spreading (Figure 3—figure supplement 1G). We found that, in control cells, the total traction energy, and thus the contractility, increased during spreading as expected, but that it was not the case in Y-27632 treated cells in which it remained very low throughout spreading (similarly to cells treated for the Arp2/3 inhibitor CK-666 which are spreading slower). Our results (Figure 2) show that the inhibition of contractility during spreading increased the volume loss, while if contractility had a direct effect on volume via the force balance at the cell surface, its inhibition should lead to a larger volume.

These observations are consistent with order-of-magnitude estimates. The osmolarity of the cell is the sum of contribution of impermeant osmolytes and different species of permeable ions. Well accepted orders of magnitudes for the ion concentrations inside the cells (about 200 mM) and outside (about 300 mM) suggest that the concentration of the impermeant osmolytes should be about 100 mM, corresponding to an osmotic pressure of about 105 Pa. Taking the cortical tension to be of order of 1 mN/m and cell radius to be about 10 µm we get the hydrostatic pressure difference (Laplace pressure) produced by the cell cortex to be about 100 Pa. It means that cortical tension would be able to impose a volume change that would ‘concentrate’ the impermeant osmolytes by a maximum of only about 1%. For small variation around this steady state, a 10% decrease in cell volume will increase the osmotic pressure of the trapped osmolytes by 104 Pa. For this osmotic pressure to be balanced by an increase in the hydrostatic pressure, the cortical tension would need to increase by a factor of 100, which we consider to be too large to be realistic. This reasoning leads to the generally accepted result that even a large increase of contractility (e.g., by a factor of 10) would only change the volume by 1%. Overall, these considerations demonstrate that the current versions of the PLM or its extensions cannot explain our observations.

These results also confirm our interpretation that decreasing contractility leads to a larger loss of volume during spreading mostly indirectly, because it increases the speed of cell spreading.

A mechano-sensitive PLM including a mechano-osmotic coupling predicts the observed relation between spreading speed and volume loss

We thus engaged in proposing a modified model (see the full model in Appendix 1), to combine PLM with cell mechanics and shape. To account for the observation of an increase in volume loss for faster spreading cells, we made the assumption, like the other models discussed above, that an element coupling cell mechanics (which is directly affected by the spreading speed) to the ion fluxes need to be added to the classical PLM. Ion channels and pumps can be affected by membrane tension, as demonstrated multiple times by others (Cox et al., 2019), we thus chose to implement this mechano-osmotic coupling, similarly to the model discussed above (Xie et al., 2018; McEvoy et al., 2020; Jiang and Sun, 2013; Tao and Sun, 2015), but with a full PLM, including permeant and impermeant solutes (see also Yellin et al., 2018; Li et al., 2021). We also chose to implement cell growth in the model, using the experimentally measured rate (about 5% volume increase per hour), because it plays a significant role at timescales of a few tens of minutes, thus overlapping with the timescale of the latest part of the spreading process.

An important assumption of such a model is that faster spreading cells display a higher membrane tension specifically during the phase of fast spreading. To test this hypothesis, we pulled membrane tubes with an atomic force spectrometer tip following a well-established protocol (see Materials and methods for further details; Diz-Muñoz et al., 2016). This allows to measure the tether force, which varies with the square root of the membrane tension in the absence of a cell cortex. It is a common readout of an apparent membrane tension, even if it is not a direct measure, because it also depends on the membrane interaction with the cell cortex (Sitarska and Diz-Muñoz, 2020). Because tether force measures take a significant time and require adhered cells, we started measures 30 min after cell plating and performed them for the following hour. These measures correspond to the spreading phase, which is comparable to the phase of volume loss in the single-cell spreading experiments (in which time zero is taken as the time when each single cell starts to spread, and not the initial seeding time, which gives on average a delay of about 10 min). This early time point designated as ‘spreading phase’ thus overlaps with the phase in which cells are still spreading and losing volume. We then performed a second measure within 4–5 hr after seeding, which is long enough after the spreading to consider that cells have recovered their steady-state growing phase and the effect of their initial spreading on volume is lost (see more details in Materials and methods). We found that cells treated with Y-27, which spread faster than control cells, displayed a higher tether force during the spreading phase, while the force was similar to control cells at steady state (Figure 3E), showing that the increase was not due to the drug treatment itself. Importantly, we also showed that cells that stayed rounded because they had been plated on small (20 µm in diameter) micropatterns, measured 4–5 hr after seeding, had a low membrane tension, treated or not with Y-27632. This experiment suggests an effect of spreading speed on membrane tension. This is consistent with the hypothesis that membrane tension might modulate cell volume upon fast cell shape changes.

Another important requirement to propose that an extended PLM could explain our observations on spreading cells is that the timescales involved in the osmotic and in the spreading phenomena match. To estimate the typical timescales of water and ion fluxes, we performed a detailed characterization of the cell response to osmotic shocks. We first made high time resolution recordings of cell swelling and shrinking upon a change in the external osmolarity (Figure 4A and Figure 4—figure supplement 1A, B and Figure 4—video 1). The change of volume occurred in a timescale of seconds, as expected. These experiments provided the rate of cell water entry and exit as a function of the difference in osmotic pressure (Figure 4B). This allowed us to estimate the permeability (Figure 4—figure supplement 1), which appeared smaller for hyper-osmotic shocks than for hypo-osmotic shocks, as reported previously (Chara et al., 2005; Peckys et al., 2011), although the reason for this difference is not understood. We next characterized the longer, minutes timescale of volume adaptation (Figure 4C). It showed that volume adapted faster for larger shocks. At the level of individual cells, the response was quite homogenous for the recovery from hypertonic shocks, while there was a higher cell-cell variability during recovery from hypotonic shocks (Figure 4—figure supplement 1D), with cells showing only partial recovery, especially for large hypo-osmotic shocks. Despite these complex single-cell behaviors, these experiments provide clear evidence, as well known from decades of studies of this phenomenon, of an active volume regulation mechanism on the timescale of minutes, setting the typical timescale for ion fluxes. These two timescales, seconds for water flows through the membrane and minutes for ion fluxes, are basic assumptions of the PLM model verified by our experiments. Importantly, the rate of volume change observed for small shocks is similar to the rate of volume change during cell spreading experiments (about 10 μm3/min). This justifies the use of a mechanosensitive PLM to explain the cell spreading data.

Figure 4 with 3 supplements see all
The mecano-sensitive PLM explains the volume loss in spreading cells.

(A) Average normalized volume of HeLa Kyoto cells during initial response to osmotic shocks of different magnitudes measured with high time resolution, 100 ms. Number of cells in the experiments: control Piso/P=1.25 (n=13, N=1), Piso/P=2.01 (n=19, N=1), Piso/P=0.78 (n=15, N=1), Piso/P=0.4 (n=17, N=1). Error bars represent standard deviation. (B) Average volume flux in HeLa Kyoto cells during initial response to osmotic shocks of different magnitudes. Each point represents average value of single experiment, average number of cells in each experiment n~12. Error bars represent standard deviation. (C) Average volume flux in HeLa Kyoto cells during regulatory volume adaptation. Each point represents average value of single experiment, average number of cells in each experiment n~48. Error bars represent standard deviation. (D) Fits of the spreading data from the model using best fit parameters on average normalized volume (right) and spreading area (left) of control cells divided into three categories represented in Figure 1H. (E) Z-plane of HeLa LifeAct-mcherry (red) cell before and after addition of distilled water, cell membrane is stained with CellMask Green (cyan). Scale bar 10 µm. (F) Left top: Volume (black) and propidium iodide (PI) intensity of single HeLa Kyoto cell exposed to distilled water. Dashed line indicates the time of distilled water addition. Reaching of maximum cell volume is followed by cell membrane rupture, volume decrease, and PI entry into the cell. Left bottom: Corresponding FXm images and PI staining. Right: Distribution of ratio between maximum volume cells reach before bursting induced by exposure of distilled water and their initial volume (n=63, N=3).

Figure 4—source data 1

Data tables related to quantifications in Figure 4.

https://cdn.elifesciences.org/articles/72381/elife-72381-fig4-data1-v2.zip

In brief, in this new model (see details in Appendix 1), the dependence of the volume on membrane tension is through the mechanosensitivity of the ion channels and pumps. In the linear regime, we assume that small changes in tension lead to a small change in ion transport rates so that the volume change is proportional to the change in tension. The equation for change in volume reads

(1) dδVViso=αdγγiso,

where δV=Vt-Viso1+rgrowtht, with rgrowth and being the growth rate of the cell and time, respectively. If the mechanosensitivity parameter is negative, then the volume will decrease upon an increase in the tension. We note that Equation 1 is valid for small changes in volume corresponding to small changes in tension. The other assumption we make in writing Equation 1 is that over the timescale of cell-spreading the ion and water transport have equilibrated. Since we are dealing with linearized equations, Equation 1 includes the complex properties of the ion transport like voltage and concentration dependence of ion channels and pumps. We explicitly evaluate the mechanosensitivity parameter, by analyzing a model of ion transport with three ion species—sodium, potassium, and chloride. We find that the sign and magnitude depend on the mechanosensitivity of the potassium and sodium channels, and on the ion concentrations before spreading. For physiological values of parameters found in the literature (see values in Appendix 1), we expect to be negative. Note that other effects not accounted for in our simple estimate, such as the voltage and concentration dependence of ion channels and the existence of co-transporters, could affect the value and sign of. To relate tension variations to the rate of cell spreading, we model surface tension using a Maxwell fluid model, with a relaxation timescale and elastic modulus, driven by the rate change of total surface area. The tension dynamics reads

(2) 1+τddtδγγiso=kτγisodAtotdt

The elastic modulus characterizes the short-time elastic response, while the relaxation timescale accounts for the existence of tension homeostasis mechanisms that have a longer response time. During cell spreading, the total surface area will increase leading to a spreading rate-dependent increase in tension, which will relax back to the homeostatic value, in agreement with the tether pulling experiments reported in Figure 3E. To estimate the total surface area, we take the cell shape to be that of a spherical cap (Figure 4—video 2; we discuss in more details the possible shape approximations and their relation to the measured spreading area in the Appendix 1). Combining the tension dynamics in Equation 2 with tension/volume coupling in Equation 1 leads to the following effective viscoelastic model for volume dynamics driven by a change in the total area,

(3) 1+τddtδVViso=-ξτ1AtotdAtotdt,

where the effective elastic modulus ξ=-Atotkα/γiso is proportional to the effective elasticity of the membrane and to the magnitude of the mechanosensitivity parameter relating volume loss to the tension increase and is also inversely proportional to the surface tension. The total area itself depends on the volume, we can write Equation 3 as

(4) 1+τeffddtVViso=1+rgrowtht+τ-ξτf2V,AcdAcdt,

where τeff=τ1+ξf1V,AcViso , f1V,Ac, and f2V,Ac are functions that are given by the geometry of the cell, which relate the change in total area to the change in volume, and change in contact area respectively. The model parameter τeff is an effective relaxation timescale for volume that depends on τ, volume, and contact area. The difference between τ and τeff stems from the complicated geometrical relationship between contact area and total area in a spreading cell. The contact area dynamics is fitted to an exponentially saturating function. The two input parameters for the models are τa : the timescale of cell spreading and A0 : the saturating value, obtained by fitting the cell spreading curve. Fitting the volume dynamics yields the two model parameters ξ and τ for the mechanosensitivity of ion transport and membrane mechanics.

These parameters allowed us to fit the various experimental data and their values are discussed in more details below. Importantly, this simple extension of the PLM predicts the observed proportionality between volume loss and the speed of spreading, and no dependency on the absolute cell spreading area (Equation 3, Figures 1I, 2C and D). We conclude that this new model is able to explain our main observation and constitutes a robust implementation of a membrane tension homeostasis mechanism within the PLM framework. We propose to call it the mechanosensitive PLM.

Fitting the spreading and volume data with the mechano-sensitive PLM

To further test the capacity of the model to explain our observations, we performed a fit of our experimental data using the mechanosensitive PLM model and analyzed the parameters obtained. We used the three groups of control cells defined in Figure 1H, sorted based on spreading speed during the first 10 min. The spreading parameters were extracted from the experimental spreading data, and the model allowed a satisfactory fit of the experimental volume data (Figure 4D).

Because the mechanosensitive PLM assumes a coupling between membrane tension and ion fluxes, how much volume is lost by a cell during spreading depends on whether the cell deforms in a rather elastic or viscous regime. The transition between these regimes is defined by the relative values of the spreading rate and the effective tension relaxation time—if the spreading rate is faster than the relaxation time, the cell deforms in a rather elastic regime, and as a result, the membrane gets tensed and the cell loses volume. The effective relaxation time depends on the two main fitting parameters, the bare tension relaxation timescale (which varies in the minutes to tens of minutes timescale) and the stiffness (which varies around one). When fitting the three classes of fast, intermediate, and slow-spreading cells, we found that the values of the fitting parameters (Appendix 1—table 3 in Appendix 1), do vary significantly for the three classes. However, this variation could not explain the difference in volume loss (see Appendix 1), which must therefore be attributed to the difference in spreading speed. We conclude that our mechanosensitive PLM not only captures properly the coupling of spreading kinetics on volume modulation but that the parameter fitting suggests that the key ingredient of the model, the finite response time of the mechano-osmotic feedback, might be the cause of the volume loss in fast-spreading cells.

Volume loss upon fast cell deformation depends on branched actin and on changes in ion fluxes

We then asked what could be the origin of the increase in surface tension during fast cell spreading. We first evaluated the total amount of cell membrane available. We exposed cells to distilled water and first imaged actin and membrane staining. It showed a rapid full unfolding of membrane reservoirs (Figure 4E) before the cell exploded. We then used propidium iodide (PI) to identify the timing of plasma membrane rupture (Figure 4F). We found that on average, the plasma membrane ruptured when cells reached 5.7 times their initial volume, which corresponds to an excess of membrane surface area of about 3.3 times, in accordance with previous measures (Ting-Beall et al., 1993; Guilak et al., 2002).

It means that cells have a very large excess of membrane surface area and that membrane tension could not arise from a limitation in the total amount of membrane. Nevertheless, the plasma membrane being bound to the underlying cytoskeleton, its restricted unfolding could generate an increase in tension depending on the rate of cell deformation. This would explain why the volume loss depends on the spreading speed. We thus further tested the role of the actin cytoskeleton in volume loss during cell spreading.

Because branched actin was shown to more specifically interact with the plasma membrane (Diz-Muñoz et al., 2016; Lieber et al., 2013) and modulate membrane tension, we used HeLa EMBL cells treated with the Arp2/3 complex inhibitor CK-666, and combined the treatment with Y-27 to induce fast spreading. We found that CK-666 treatment alone induced both a slower spreading and lower volume loss (2–3%, Figure 5A and B), similar to the low Lat A treatment (Figure 2B and C), which was well fitted by the mechanosensitive PLM (fits on Figure 5A and B). Treatment with Y-27 increased the spreading speed of CK-666 treated cells, but the volume did not decrease in this fast-spreading condition (Figure 5C and D). The mechanosensitive PLM could fit these data by adjusting the parameter coupling surface tension to the change of activity of ion pumps (Appendix 1—table 4 in Appendix 1). To directly test the role of ion fluxes in the volume loss, we targeted two main players: first, stretch-activated calcium channels (including Piezo), using gadolinium chloride (GdCl3) and second NHE1, the sodium/proton exchanger, using EIPA. Neither of them have a direct role in volume regulation (Figure 5—figure supplement 1A), because there are too few calcium ions in the cell (for Piezo) or protons (for NHE1). Affecting direct volume regulation by targeting the transport of sodium or potassium ions would be too detrimental for the cell and change dramatically the initial cell volume. By contrast, inhibition of these two channels did not affect the initial average cell volume (Figure 5—figure supplement 1 A). Piezo was chosen because it is known to be mechanosensitive and because calcium acts upstream of many fast cell response pathways, including regulating cell contractility (Clapham, 2007) and it was proposed before to play a role in cell volume regulation (Hua et al., 2010). We also chose to target NHE1 because its inhibition was previously shown to affect fast cell volume changes, at mitotic entry (Byun et al., 2015; Miettinen et al., 2021) and during cell spreading (Xie et al., 2018), as we confirmed in our experiments (Figure 5—figure supplement 1B). Treatment with GdCl3 led to an increase in volume loss (from 5% to 8%), which could be fully accounted for by the increase in cell spreading speed (Figure 5E and F) and well fitted by the mechanosensitive PLM. Because the effect of GdCl3 can be purely accounted for by the increase in spreading speed and thus does not disrupt the mechanosensitive PLM, we speculate that it could be explained by a reduction of cell contractility, since preventing calcium entry might reduce actomyosin motors activation (Clapham, 2007), and thus increase the spreading speed (Wakatsuki et al., 2003). On the other hand, HeLa EMBL cells treated with EIPA, while spreading slightly faster than control cells, lost less volume (Figure 5G and H). Combining EIPA with Y-27 showed that, despite a fast spreading speed comparable to Y-27 treated cells, NHE1 inhibition fully prevented volume loss (Figure 5G and H), an effect that we also observed for RPE-1 cells (Figure 5—figure supplement 1C). Inhibition of NHE1, which is known to affect ion transport, is thus fully preventing volume loss during fast spreading. Importantly none of the drug treatments performed significantly affected the initial volume of cells prior to spreading (Figure 5—figure supplement 1A). This is consistent with the importance of changes in ion fluxes in the mechanosensitive PLM.

Figure 5 with 1 supplement see all
The mecano-sensitive PLM depends on branched actin and modulation of ion fluxes and constitutes a surface tension homeostasis mecanism.

(A) Two parameter fits for the spreading kinetics using the exponential saturation anzatz (see text) on average area of control cells (blue, n=73, N=1), 100 nM Latrunculin A (orange, n=30, N=2) or 100 µM CK-666 (magenta, n=37, N=2) treated. Error bars represent standard deviation. (B) Fits from the model using best fit parameters on average normalized volume of control cells (blue, n=73, N=1), 100 nM Latrunculin A (orange, n=30, N=2) or 100 µM CK-666 (magenta, n=37, N=2) treated. Error bars represent standard deviation. (C) Two parameter fits for the spreading kinetics using the exponential saturation anzatz (see text) on average area of control cells (gray, n=73, N=1) or combination of 100 µM CK-666 and 100 µM Y-27632 (violet, n=24, N=1) treated. Error bars represent standard deviation. (D) Fits from the model using best fit parameters on average normalized volume of control cells (gray, n=73, N=1) or combination of 100 µM CK-666 and 100 µM Y-27632 (violet, n=24, N=1) treated. Error bars represent standard deviation. (E) Two parameter fits for the spreading kinetics using the exponential saturation anzatz (see text) on average area of control cells (gray, n=73, N=1) or 100 µM GdCl3 (cyan, n=30, N=2) treated. Error bars represent standard deviation. (F) Fits from the model using best fit parameters on average normalized volume of control cells (gray, n=73, N=1) or 100 µM GdCl3 (cyan, n=30, N=2) treated. Error bars represent standard deviation. (G) Two parameter fits for the spreading kinetics using the exponential saturation anzatz (see text) on average area of control cells (gray, n=73, N=1), 100 µM Y-27632 (red, n=21, N=1), 50 µM EIPA (green, n=73, N=1), or combination of 50 µM EIPA and 100 µM Y-27632 (dark cyan, n=30, N=2) treated. Error bars represent standard deviation. (H) Fits from the model using best fit parameters on average normalized volume of control cells (gray, n=73, N=1), 100 µM Y-27632 (red, n=21, N=1), 50 µM EIPA (green, n=73, N=1), or combination of 50 µM EIPA and 100 µM Y-27632 (dark cyan, n=30, N=2) treated. Error bars represent standard deviation. (I) Predicted by model, plots for difference between tension without mechano-osmotic coupling (for α = 0 and kγiso = 100) and tension with mechano-osmotic coupling (for fitted and kγiso = 100). (J) Tether force measurements of control HeLa Kyoto cells (gray, for ‘spreading phase’ n=50, N=6; for ‘steady state’ n=55, N=9), treated with Y-27632 (gray, for ‘spreading phase’ n=27, N=3; for ‘steady state’ n=21, N=3), CK-666+Y-27 (purple, for ‘spreading phase’ n=25, N=3; for ‘steady state’ n=19, N=3), EIPA (green, for ‘spreading phase’ n=23, N=3; for ‘steady state’ n=18, N=3), EIPA+Y-27 (dark cyan, for ‘spreading phase’ n=23, N=3; for ‘steady state’ n=15, N=3) during the first 30–90 min after plating or 4–5 hr after plating. Error bars represent standard deviation. The results of statistical tests are shown at the graph. (K) Volume flux (dV/dt) of single control HeLa Kyoto cells (n=194, N=3), treated with Lat A (n=41, N=3), CK-666 (n=54, N=3), Y-27 (n=121, N=4), EIPA (n=117, N=3), GdCl3 (n=53, N=3), CK-666+Y-27 (n=74, N=3), EIPA+Y-27 (n=50, N=3), MnCl2 (N=3, n=57) plotted versus their spreading speed (dA/dt) at the first 10 min of spreading. Error bars represent standard error. (L) Scheme of mechanosensistive “pump-leak” model. (M) Scheme representing cell volume regulation in response to deformations.

Figure 5—source data 1

Data tables related to quantifications in Figure 5.

https://cdn.elifesciences.org/articles/72381/elife-72381-fig5-data1-v2.zip

The mechano-osmotic coupling moderates the membrane tension increase in fast-spreading cells, acting as a membrane tension homeostasis mechanism

The mechanosensitive PLM predicts that inhibition of the mechano-osmotic coupling in fast-spreading cells, would prevent the associated volume loss, and lead to an increase in the membrane tension during the spreading phase. To give a qualitative prediction of this effect, the membrane tension value was extracted from the model, for the various conditions tested experimentally, using the fit on the experimental data, and the tension values were compared for the case of a model with or without a mechano-osmotic coupling (the difference between the two values is given in Figure 5I). This model prediction shows that the largest increase in tension, in case the mechano-osmotic coupling was absent, is expected in the case of the fastest spreading cells (Y-27 or GdCl3 treatments). On the other hand, in case the mechano-osmotic coupling is already disrupted by the treatment (e.g., in the case of EIPA treatment), or if cells spread very slowly (e.g., Lat A treatment), the predicted difference is small. To test this prediction experimentally, we performed tether pulling experiments (Figure 5J). These experiments showed that, as predicted by the model, disrupting the mechano-osmotic coupling in fast-spreading cells (Y-27 plus CK-666 or Y-27 plus EIPA) leads to the highest tether force values in the spreading phase, while the steady-state values did not change. Similarly, combined EIPA and Y-27 treated cells showed higher tension than Y-27 or EIPA alone. Tension was highest during early spreading compared to steady-state spread cells, suggesting that the increase was due to spreading and not to the drug treatments alone, even though EIPA alone also add an effect on steady-state tension. This shows that, in these cells, the coupling between membrane tension and volume regulation is lost, and that fast spreading in the absence of volume loss induces higher tension increase (Figure 5J and K), as predicted by the model (Figure 5I). Taken together, these experiments confirm the validity of our mechanosensitive PLM. They also support the existence of a membrane tension homeostasis mechanism that reduces the extent of changes in membrane tension upon fast cell shape changes by modulating the relative contribution of surface expansion and volume loss.

Discussion

A mechano-osmotic coupling leads to volume loss in fast-spreading cells

Our detailed characterization of cell volume during cell spreading revealed that, while cell volume is not related to the steady-state shape of the cell, it is modulated by the rate of cell shape change. We propose that this is due to a coupling between cell membrane tension and rates of ion fluxes (Figure 5L). An extension of the classical PLM including this coupling can account for our observations of cell volume during cell spreading. Measures of membrane tension during spreading and at steady state under a variety of conditions confirmed that fast spreading is associated with a transient increase in membrane tension and that preventing volume modulation leads to even higher membrane tension, as predicted by the model. Taken together, these experiments and this model suggest the existence of a mechano-osmotic coupling at the level of the cell membrane, which acts as a membrane tension homeostasis mechanism by reducing membrane tension changes upon fast cell deformation (Figure 5M).

The role of membrane binding to the actin cortex in inducing membrane tension and volume loss in fast-spreading cells

A central hypothesis in the model is that the physical coupling between the actin cortex and the cell membrane leads to an increase in membrane tension when the rate of deformation is faster than the relaxation time of the actin cortex and membrane ensemble. To verify this hypothesis, we performed membrane tether experiments in various conditions, during spreading and at steady state (Figure 5J). Membrane-to-cortex attachment is at least partly mediated by proteins of the ERM family Chugh and Paluch, 2018. Thus, we performed an additional experiment using an inhibitor of Ezrin, membrane-cytoskeleton linker (20 µM NSC668394) and monitored cell volume during spreading. We found that while spreading was similar or even slightly faster during the initial phase, treated cells lost less volume than control cells (Figure 6A), consistent with a role of cortex/membrane coupling in mediating the effect of spreading kinetics on volume loss. To further investigate the ultrastructure of the cell cortex during spreading, we unroofed (Vassilopoulos et al., 2019) Hela cells after 30 min spreading on a fibronectin-coated substrate. Following electron microscopy imaging confirmed the different extent of spreading in the various conditions assayed, and the perturbation of branched actin in the CK-666 treated cells (Figure 6B). Membrane folds and structures such as clathrin-coated pits and caveolae were present in all conditions. Although their number and degree of curvature did not change significantly, the different populations of caveolae appeared hard to quantify and compare on spreading cells without underestimating the amount of flat caveolae. We conclude that, while our biophysical measures and fluorescence imaging gave a clear indication of changes in membrane tension and folding state during large cell deformations, further investigations are needed to precisely describe the change in the state of the membrane and its association to the actin cortex in this context. In particular, branched actin is known to be important for endo- and exocytosis (Toshima et al., 2005; Tran et al., 2015) and, on the timescale of several minutes to tens of minutes, it could regulate membrane tension, especially during cell spreading (Gauthier et al., 2011) and thus contribute to cell volume modulation.

Role of membrane to actin attachment on volume loss and state of the plasma membrane in spreading cells.

(A) Left: Average spreading area of control HeLa Kyoto cells (gray, n=125, N=1) or 20 µM NSC (light blue, n=101, N=3). Middle: Average normalized volume of control HeLa Kyoto cells (gray, n=125, N=1), or 20 µM NSC (light blue, n=101, N=3). Error bars represent standard deviations. Right: Volume flux (dV/dt) plotted versus their spreading speed (dA/dt) of single control HeLa Kyoto cells and treated with various drugs and represented in Figure 5K and treated with 20 µM NSC (light blue, n=101, N=3). Error bars represent standard error. (B) Platinum replica electron microscopy survey views of the cytoplasmic surface in control, Y-27632, CK-666, or CK-666+Y-27632-treated unroofed Hela cells spread on glass coverslips for 30 min. Extracellular substrate is pseudo-colored in red. For each panel, high magnification views corresponding to the boxed regions are shown on the right.

Figure 6—source data 1

Data tables related to quantifications in Figure 6.

https://cdn.elifesciences.org/articles/72381/elife-72381-fig6-data1-v2.zip

The sign of the volume change upon fast cell deformation

Although the effect of mechanosensitive ion channels on volume has been discussed before for simplified systems considering one or two solutes (Guo et al., 2017; Xie et al., 2018; Adar and Safran, 2020; Jiang and Sun, 2013), the relation between mechanosensitivity of ion channels and pumps and volume change is far from obvious. The sign of the volume change upon an increase in tension depends on whether the contribution of ions to the osmotic pressure increases or decreases. Since the cell is always osmotically balanced, if the concentration of the ions inside the cell decreases, the concentration of the trapped molecules should increase by decreasing the volume such that the cell osmotic pressure stays constant and vice-versa. For instance, an increase in sodium channel conductance upon an increase in tension leads to an increase in volume, whereas an increase in potassium conductance upon an increase in tension leads to a decrease in volume. Using a detailed model of ion transport, we show that for physiological values of parameters, as observed in experiments, the volume is indeed expected to decrease upon an increase in tension.

Values of fitting parameters for the mechanosensitive PLM suggest a role for branched actin in modulating ion fluxes

This new mechanosensitive PLM gives a Maxwell viscoelastic model for the volume with two fitting parameters, the effective stiffness, and the bare relaxation timescale. In most cases, we get a good quantitative fit for cells treated with different drugs that perturb the cytoskeleton and the ion channels (Appendix 1). For the Y-27 and Lat A treated cells, the fits only qualitatively capture the temporal dynamics of the volume. One of the reasons for an imperfect fit for these two drugs could be the failure of the spherical cap approximation used to estimate the surface area (see more discussion on the shape estimates and the parameters used for cell surface area in the model, Appendix 1). Control cells, and cells treated with Y-27, EIPA, and CK-666 show less than 30% variation in the value of ξ, implying that most of the volume loss is explained by differences in spreading speed. Cells treated with GdCl3 show a larger decrease of 60% but stay in the same range of parameters (and they are close to the same line in the dV/dt versus dA/dt summary graph shown in Figure 5K). However, for the cells treated with Y-27+EIPA and Y-27+CK-666, ξ decreases by an order of magnitude, leading to low volume loss even though the cells are spreading fast (Appendix 1—table 4 in Appendix 1). This decrease of ξ could be either due to a decrease in the elasticity of the membrane or due to the decrease in the value of the mechanosensitivity parameter. Spreading experiments show that, for both Y-27+CK-666 and Y-27+EIPA, membrane tension reaches the highest values. This means that, in both cases, spreading is still inducing an increase in membrane tension, and the absence of volume loss reinforces the effect on membrane tension. It suggests that the elasticity parameter is not affected but rather the volume-tension electromechanical coupling. This could mean that, unexpectedly, branched actin networks are specifically required for this coupling. This could be due to a direct association of branched actin with ion channels and pumps (Mazzochi et al., 2006; Shaw and Koleske, 2021).

Volume loss in ultra-fast deforming cells

The PML and our mechanosensitive extension are meant to explain volume changes at the minutes timescale, which correspond to the time needed for large enough ion fluxes to take place. To test the limit of validity of the model, we imposed fast (less than a second timescale) deformation on the cell, we used our previously developed cell confiner (Le Berre et al., 2014). This device can impose a precise height on cells and thus gives access to a large range of deformations (Figure 7A and Figure 7—figure supplement 1A and Figure 7—video 1). RICM measure of the cell contact area showed a range of spreading similar to what was observed during spontaneous cell spreading (Figure 7B and Figure 7—figure supplement 1B). In addition, imaging of the plasma membrane showed that confinement below 10 μm induced a clear loss of membrane folds and reservoirs (Figure 7C images and Figure 7—video 2), while treatment with Lat A induced the formation of large membrane blebs and less extension of the cell diameter upon confinement (Figure 7C graph and Figure 5—figure supplement 1C and Figure 7—video 2). This suggests that cell confinement, like hypo-osmotic shocks, induces membrane reservoir unfolding, and that Lat A treatment, by reducing the membrane anchorage and causing bleb formation, reduces the surface expansion following confinement. FXm volume measurement combined with confinement showed a strong loss of volume of confined HeLa EMBL control cells, while Lat A treated cells kept a constant volume (Figure 7D and Figure 7—figure supplement 1D showing for both treatments the decrease in FXm background intensity corresponding to the confiner height; and Figure 7—video 3). In control cells, stronger confinement led to larger volume loss, while Lat A treated cells showed no significant volume loss except for the lowest confinement height (Figure 7E). Conversely, treatment with Y-27632 to reduce contractility had no effect on volume loss upon confinement (Figure 7E), confirming that contractility has no direct effect on cell volume and that the effect of Y-27632 on cell volume during cell spreading is mostly via an increase of the spreading speed. The loss of volume in control cells corresponded to a deformation at an almost constant surface area (Figure 7F, calculated from the volume, see Appendix 1). Below 5 μm height, the cell surface significantly increased, which also corresponded to the formation of large blebs (Figure 7—figure supplement 1E). This loss of volume induced by fast confinement was also found in other cell types, RPE-1 and HEK-293 (Figure 7—figure supplement 1F, G) and was also previously observed in confined Dictyostelium cells (Srivastava et al., 2020). Overall, these experiments show that fast imposed cell deformation induces an actin-dependent loss of volume (up to 30%), at almost constant surface area.

Figure 7 with 6 supplements see all
Volume modulation during ultra-fast cell flattening and during cell migration through collagen matrices.

(A) 3D-membrane reconstruction of HeLa expressing MyrPalm-GFP (black) cells cell shape under different confinement heights, side view. Scale bar: 10 µm. (B) Contact area with bottom glass substrate of Hela Kyoto cells under different confinement heights. Average number of cells in each experiment n~79, for each condition N=1. Error bars represent standard deviation. (C) Left: Z-plane of control and 2 µM Lat A treated HeLa cells under 20 µm and 7.6 µm confinement heights. Cell membrane is stained with CellMask Far Red (white). Scale bar: 10 µm. Right: Average CellMask intensity plotted versus distance from cell body contour on the middle Z-plane of HeLa-MYH9-GFP-LifeAct-mcherry cells. Number of cells in each condition n=10, N=1 for each condition. Error bars represent standard deviation. (D) Average normalized volume of control (blue, n=48, N=1) and 2 µM Lat A treated (orange, n=32, N=1) HeLa Kyoto cells during dynamic confinement experiment. Dashed line indicates the moment of confinement. Error bars represent standard deviation. (E) Average normalized volume of HeLa Kyoto cells (blue) and cells treated with Lat A 2 µM (orange) or 5 µM (yellow) or 100 µM Y-27 (red) under different confinement heights. Each data point represents an average of N~10 experiments; each experiment contains n~160 individual cells. Error bars represent standard deviation. There is statistically significant difference between control and Lat A for the heights ~12 µm and ~7 µm (p=0.02 and p<10–8). There is no statistically significant difference between control and Y-27 for the heights ~12 µm and ~7 µm (p=0.4 and p=0.3). (F) Projected surface (computed from volume represented in panel Figure 5G) of HeLa Kyoto cells under different confinement heights. Dashed line indicates the confinement height that corresponds to blebs appearance. Error bars represent standard deviation. (G) Volume of single HeLa Kyoto cells during dynamic confinement experiment. Dashed line indicates the moment of confinement. (H) FXm images of HeLa Kyoto cells during dynamic confinement experiment taken with high NA objective. (I) Top: Volume of single DC migrating in collagen. Bottom: Corresponding FXm images. (J) Volume of single DCs migrating in collagen with the different speeds, N=1. Left: <0.5 µm/min (n=14), middle: >0.5 µm/min (n=10) and <1 µm/min, right: >1 µm/min (n=19). (K) Coefficient of variation of volume flux dV/dt computed for 10 min intervals during single DCs migration in collagen plotted versus their average speed (n=43, N=1). The line is the fit for coefficient of variation of DCs volume using best-fit parameters (see Appendix 1).

To better estimate the speed of deformation imposed by the confiner, we imaged at high frame rate during the confinement process. It showed that, even with a time lapse of 30 ms, the volume loss happened between two consecutive frames (Figure 7G a Figure 7—figure supplement 1H, I and Figure 7—video 4). Only volume is lost and not dry mass (Figure 7—figure supplement 1J), which suggests that only water and probably small solutes are lost. Nevertheless, the speed of volume change is not compatible with our mechanosensitive PLM, as in this model, volume loss occurs in the minutes timescale due to changes in ion transport rates. Fast imaging of the fluorescent medium surrounding the cells used for FXm indeed showed a transient appearance of streams of darker fluid (non labeled, thus coming from the cells) emanating from confined groups of cells (Figure 7H and Figure 7—video 4), likely corresponding to the expelled water and osmolites. Overall, these confinement experiments suggest that, although at this timescale of milliseconds, the mechanism of volume loss very likely differs from the context of spontaneous cell spreading, it is also induced by an increase in membrane tension, and requires the presence of the actin cytoskeleton.

Within the PLM framework, for a given osmolarity of an external medium, the cell volume may change either due to a change in hydrostatic or osmotic pressure. Fast compression can increase the cortical tension, which can cause an increase in hydrostatic pressure of the cell. However, the maximum hydrostatic pressure in the cell before the membrane detaches from the cortex is of the order of 102 Pa, thus producing no direct effect on the cell volume, as discussed before (Cadart et al., 2019; Salbreux et al., 2012). Hence, the observed volume loss of 10–15% can only be due to a change in the osmolarity of the cell, and not to a change in hydrostatic pressure. For ion transport to take place at timescales of milliseconds, the transport rates of channels and pumps would need to increase by 4 orders of magnitude. Such an increase can be easily attained if the high tension upon compression leads to transient formation of pores in the plasma membrane (observed in spreading GUVs; Karatekin et al., 2003). If these pores are small enough to allow for free ion transport but do not let the larger molecules trapped in the cell pass through (which should be the case since the dry mass was found to remain constant), the cell volume will increase rather than decrease (a consequence of the Donnan effect; Sperelakis, 2012). The formation of pores thus cannot explain our observations. Another mechanism that may lead to volume decrease upon compression without losing the trapped osmolites requires a selective increase of the ion conductance upon compression, but by orders of magnitude. Whether the ion conductance can increase by 4 orders of magnitude by mechanical stretching requires further investigation. Finally, it is also possible that due to its poroelastic nature (Moeendarbary et al., 2013; Charras et al., 2009), the cytoplasm behaves as a gel-like structure, and that water and osmolites are pressed out of the cell upon confinement, without changing the osmotic balance nor the dry mass (Sachs and Sivaselvan, 2015). In conclusion, confinement experiments confirm that fast deformation is associated with volume loss in an actin-dependent way, also suggesting a coupling between cell mechanics and volume regulation. However, they are hard to fully interpret in physical terms. This means that such a simple experiment as squeezing a cell cannot yet be understood with the current general knowledge on cell biophysics, pointing to a need for further investigations of the physics of large cell deformations. Such deformations are likely to occur in physiological contexts such as circulation of white blood cells and circulating cancer cells through small capillaries and may lead to volume change as was shown in vitro (Liu et al., 2020).

Volume fluctuations in fast migrating immune cells can be explained by the mechanosensitive PLM

While our mechanosensitive PLM might be limited in the interpretation of cell deformations occurring below the second timescale, it captures well the larger timescales, based on a modulation of ion fluxes by membrane tension. Such timescales correspond to deformations that cells experience, for example, as they migrate through dense tissues. This implies that migrating cells might display volume fluctuations. To test this prediction, we used a classical cell migration assay with fast-moving bone marrow-derived dendritic cells (DCs) from mice embedded in a collagen gel (Vargas et al., 2016). The collagen gel mixed with fluorescent dextran was assembled inside a cell volume measurement chamber (Figure 7—figure supplement 1K and Figure 7—video 5). Because of the low fraction of collagen in the solution and the homogeneity of the fluorescent background, regular FXm measurements could be performed. We observed that the cell volume changed by a few percent as single cells moved through the collagen gel (Figure 7I), with periods of cell protrusion corresponding to a decrease in cell volume. To assess whether these fluctuations in volume were related to the migration of cells, we split individual cells into three groups according to their average speed and plotted their volume (in %) as a function of time (Figure 7J). This clearly showed that faster moving cells displayed larger volume fluctuations. Finally, to get a more quantitative assessment of the correlation, we plotted the coefficient of variation of the volume against the speed (Figure 7K), for single cells shown in (Figure 7J). Faster cells displayed more volume fluctuations. Interestingly, this relation was well fitted by an extension of the model to cells moving through a meshwork (fit in Figure 7K, and see Appendix 1 for the model extension). This experiment suggests that the mechano-osmotic coupling that we describe in our study is at work in migrating cells, inducing larger volume fluctuations (and thus larger density changes) in faster migrating cells. These volume and density fluctuations could thus be present in a large range of cells in physiological conditions, with yet unknown consequences on cell physiology and behavior.

Beyond the potential functional significance of volume and density fluctuations associated with cell shape changes, our observations and our model demonstrate that a membrane tension homeostasis mechanism is constantly at work in mammalian cells. This mechanism is most likely due to crosstalk between mechanical, osmotic, and electrical properties of the cell pointing to the importance of taking into account complex coupling between various physical parameters to understand cellular physics and physiology.

Materials and methods

Key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
Cell line (Homo sapiens)HeLa EMBL (Kyoto)Gift from Valérie Doye
Cell line (H. sapiens)RPE1ATCC
Cell line (H. sapiens)3T3ATCC – from Alba Diez-Munoz Lab, EMBL, Heidelberg, Germany
Cell line (H. sapiens)HEK-293Gift from Liam Holt lab, NYU, NewYork
Chemical compound, drugFetal bovine serumPAN-BiotechP30-193306Use at 10%
Chemical compound, drugDextran, Alexa Fluor 647; 10,000 MW, Anionic, FixableSigma-AldrichD22914Stock at 10 mg/ml in PBS
Chemical compound, drugFluorescein isothiocyanate-dextran; 10,000 MWSigma-AldrichFD10SStock at 50 mg/ml in PBS
Chemical compound, drugFibronectinSigma-AldrichF1141-1MG50 µg/ml in PBS
Chemical compound, drugPoly-L-lysineSigma-AldrichP8920Use at 0.01%
Chemical compound, drugPLL-PEGSuSoS0.1 mg/ml solution in HEPES
Software, algorithmSoftware for FXm image analysis and volume calculationAvailable upon request to the authorsRRID:SCR_001622

Cell culture and drug treatment

Request a detailed protocol

HeLa EMBL (Kyoto), and derived cell lines HeLa LifeAct, HeLa Myrpalm-GFP-LiFeact mCherry, Hela hgem-mCherry, RPE-1, 3T3-ATCC, and HEK-293 cells were maintained in Dulbecco’s modified Eagle medium with Glutamax (DMEM/Glutamax; Gibco) supplemented with 10% fetal bovine serum (FBS; Biowest or PAN-Biotech) and 1% penicillin-streptomycin solution Thermo Fisher Scientific, and stored at 37°C and 5% CO2. All cell lines were regularly tested for mycoplasma contamination.

Bone marrow-derived DCs were obtained by differentiation of bone marrow precursors for 10 days in DCs medium (IMDM-Glutamax, FCS 10%, pen-strep 100 U/ml, and 2-ME 50 μM) supplemented with granulocyte-macrophage colony-stimulating factor (GM-CSF)-containing supernatant (50 ng/ml) obtained from transfected J558 cell line, as previously described (Barbier et al., 2019).

Latrunculin A (Sigma-Aldrich, used in the final concentrations: 100 nM, 2 or 5 µm), CK-666 (Sigma-Aldrich, used in the final concentrations: 100 µM), EIPA (Tocris Bioscience, used in the final concentrations: 50 µM), NSC668394 (Sigma-Aldrich, used in the final concentrations: 20 µM) dissolved in DMSO (Sigma-Aldrich), Y-27632 (Tocris Bioscience, used in the final concentrations: 100 µM), and GdCl3 (Sigma-Aldrich, used in the final concentrations: 100 µM) dissolved in H2O. Manganese(II) chloride solution (Sigma-Aldrich) was used in the final concertation of 1 mM. Incubations with drugs were done for suspended cells 30 min prior to experiments.

For volume measurements, 10 kDa dextran conjugated with different fluorophores were used in the final concentration of 1 mg/ml: fluorescein isothiocyanate-dextran (Sigma-Aldrich) or Alexa Fluor 647 (Thermo Fisher Scientific).

For serum starvation experiments, plated cells were incubated overnight in DMEM without FBS. Prior to the experiments, cells were detached with EDTA and resuspended in the DMEM without FBS collected from cells or in the fresh DMEM supplemented with 10% FBS and incubated for 30 min in suspension.

Cell cycle stage detection

Request a detailed protocol

The cell cycle state of the cells is indicated by the expression of h-Geminin protein which is expressed by cells from the start of S phase until mitosis (Sakaue-Sawano et al., 2013) in HeLa hgem-mCherry cell line. To quantify the fluorescence of geminin in the nucleus, first, a background subtraction is performed on the images using the ImageJ software. An ROI is used to define an area containing the background fluorescence in the image. An average value of the ROI is then subtracted from all the frames. Subsequently, an ROI is drawn to drawn as close to the cell, as possible, and then the mean gray value is measured across all the frames.

Monitoring of cell volume and contact area while spreading

Request a detailed protocol

PDMS chambers were prepared as described in Cadart et al., 2017. The typical height of PDMS chambers for volume measurements was 20 µm. PDMS chambers were incubated with 50 µg/ml fibronectin (Sigma-Aldrich) in phosphate-buffered saline (PBS) for 1 hr, washed and incubated overnight with culture medium. Cells were detached with warm Versen (Gibco) and resuspended in medium collected from cells to facilitate spreading.

In case of measurements of non-adherent cells, we used chambers incubated with PLL-PEG coating (0.1 mg/ml solution in HEPES, SuSoS), washed and incubated overnight with culture medium without FBS. Cells were detached with Trypsin and resuspended in a fresh culture medium.

The cell volume measurement explained in details in Cadart et al., 2017 and used in were coupled with spreading area measurement performed by RICM (Rädler and Sackmann, 1993; Cuvelier et al., 2007). Microscopy was performed at 37°C with 5% CO2 atmosphere. Imaging was started immediately after cell injection into the chamber with 1-min time interval. Imaging was performed using a ZEISS Z1 Observer epi-fluorescence microscope equipped with an Orca-Flash 4 Camera (Hamamatsu), 20× Plan-Apochromat objective, NA0.8 and the software Metamorph (Molecular Devices).

The volume extraction was performed with a MatLab software as described in Cadart et al., 2017.

The analysis of spreading and contact area was performed manually using the ImageJ software. The borders of the cell were delimited manually and then the area, and different shape descriptors were extracted. For the volume and spreading area data, first experimental point was taken not as experimental time point t=0, but at the first point where contact area was detectable by RICM and not exceeding 100 µm2. For average volume and spreading area during spreading values of different cells were averaged at every time point. As the measurements of spreading area and volume were done not at every time point, start time point differs from cell to cell, and because of the different duration of the experiments, averaging leads to the appearance of outliers that were deleted manually. For better visualization of the experiments, we used continuous standard deviation at the graphs instead of error bars by using matplotlib.axes.axes.fill_between. For better visualization to avoid the gaps in the continuous standard deviations, we used the average of closest points in case a value was missed by the reasons explained above.

For HeLa and 3T3 cells, initial speed of spreading dAdt and volume flux dVdt was calculated as linear slope in the first 10 min after measurable cell to substrate contact. For RPE-1 cells, initial dAdt and dVdt were calculated as linear slopes in the 10 min prior to the time point when spreading area is equal to cross-section area of cell in initial non-spread state, and in the first 10 min after that time point.

Measurements of cell volume and spreading area at the level of population (and for micropatterns experiments) were done 4 hr after cell seeding.

Micropatterning

Request a detailed protocol

Cells were patterned using the existed technique (Azioune et al., 2011) or PRIMO (Alveole) in case of tether pulling experiments.

Side-view microscopy

Request a detailed protocol

Glass slide was attached to glass bottom dish by UV-glue, the position of glass was slightly tilted from perpendicular to the dish bottom. Glass was coated with fibronectin and washed with medium. Cells were detached with Versen and resuspended in warm medium collected from cells and incubated for 30 min. Then drop of cell was added to the dish, close to the angle between dish bottom and attached glass. Dish was placed to the incubator for 2 min to allow cell initial attachment to the tilted glass. Then 2 ml of medium collected from cells was added to the dish and microscopy started with a time frame of 1 min. Imaging was performed using a ZEISS Z1 Observer epi-fluorescence microscope 20× NA0.4.

Monitoring of cell volume during cell migration in the collagen

Request a detailed protocol

Collagen mix was prepared on ice to delay polymerization: 25 µl 10× PBS + 25 µl culture medium + 55 µl collagen + 140 µl culture medium with DCs (2*106 /ml) + 5 µl FITC-dextran + 1.3 µl NaOH.

Immediately after mixing, suspension was added into PDMS chamber for volume measurements with height of 12 µm. Microscopy was started ~10 min after injection. Imaging was performed using a ZEISS Z1 Observer epi-fluorescence microscope equipped with 20× NA0.8.

Cell velocity during migration in collagen gel was calculated for 10 min intervals. Cell position was defined as a center of mass of a binary mask applied on FXm images of cells.

Monitoring of cell volume during osmotic shock

Request a detailed protocol

PDMS chambers were coated with 0.01% PLL (Sigma-Aldrich) to prevent cell detachment during changing medium and maintaining cell round shape during experiment, then washed and incubated overnight with culture medium without FBS. Cells were detached with Trypsin. Isoosmotic medium was exchanged to the medium with known osmolarity typically for 2.5 min after beginning of acquisition. Full medium exchange in the chamber takes less than 1 s. Imaging was performed using a ZEISS Z1 Observer epi-fluorescence microscope equipped with 20× NA0.4, and 20× NA0.8 in case of the stream movies. Hypoosmotic solutions were made by water addition to culture medium, hyperosmotic by addition PEG400. Osmolarity of working solutions was measured by osmometer Type 15 M (Löser Messtechnik).

Cell rupture in response to distilled water exposure was monitored by PI (1 µg/ml) (Sigma-Aldrich) intensity inside the cell.

Volume flux for passive response to osmotic shock was defined as a linear slope at the linear region of volume curves defined manually.

Adaptation speed for osmotic shock recovery was calculated as a linear slope starting from the minimum or maximum volume value achieved during passive response (for hyper or hypoosmotic shock) at 5-min interval.

Monitoring of cell volume under confinement

Request a detailed protocol

Cells were detached with Trypsin and resuspended in fresh culture medium. Both static six-well confiner and dynamic confiner were used according to experimental procedure described in Le Berre et al., 2014. Imaging was performed using a ZEISS Z1 equipped with 20× long-distance objective NA0.4.

For volume measurements performed with dynamic confiner experiments, bottom glass was coated with 0.01% PLL, that prevented cell escape from the field of view and allowed following the same cells before and after confinement.

Calculation of surface area of non-confined cell was done with the assumption of spherical cell shape, and of confined cells with the assumption of cylindrical cell shape, based on measured cell volume.

Spinning disk microscopy

Request a detailed protocol

Qualitative imaging for osmotic shock and confinement experiments was performed with spinning disk set-up (Leica DMi8). 63× and 100× oil objectives were used. CellMask (Invitrogen) staining was performed in warm PBS solution (1 µl of dye to 1000 µl PBS).

Filopodia were manually segmented. Filopodia density is plotted as number of filopodia per µm of cell body diameter. Bleb was manually segmented from middle plane images. For membrane density measurements on cell contour, cells were background substrated, and resliced by their contour, where most of the membrane marker accumulates. An average projection was plotted for 3 µm around the cell edge.

Dry mass measurements

Request a detailed protocol

Mass measurement was performed by quantitative phase microscopy using Phasics camera (Aknoun et al., 2015). Images were acquired by Phasics camera every 15 min for 35 hr during the duration of the experiment. To get the reference image, 32 empty fields were acquired on the PDMS chips and a median image was calculated. Custom MATLAB scripts were written by Quantacell for analysis of interferograms (images acquired by phasics). The interferograms were associated with reference images to measure the optical path difference and then separated into phase, intensity and phase cleaned images (background set to 1000 and field is cropped to remove edges). Background was then cleaned using gridfit method and a watershed algorithm was used to separate cells that touch each other. Mass was then calculated by integrating the intensity of the whole cell.

Tether pulling

Request a detailed protocol

For apparent membrane tension measurements, tether force was measured with single-cell atomic force spectroscopy by extruding tethers from the plasma membrane on top of the nucleus of HeLa EMBL cells. Cellview glass bottom dishes (Greiner) were coated for 1 hr with fibronectin (50 μg/ml; Sigma-Aldrich). Cells were incubated for 30 min with in the presence of drugs or vehicle, then plated, and probed either during spreading (from 30 to 90 min after plating) or at steady state (fully spread; from 4 to 5 hr after plating). To perform experiments on non-spread cells, fibronectin-coated circles (Ø 20 um) were micropatterned onto Cellview glass bottom dishes (Greiner) using PRIMO (Alveole) following the manufacturer’s recommendations.

Tether extrusion was performed on a CellHesion 200 BioAFM (Bruker) integrated into an Eclipse Ti inverted light microscope (Nikon). OBL-10 Cantilevers (spring constant ~60 pN/nm; Bruker) were mounted on the spectrometer, calibrated using the thermal noise method (reviewed in Houk et al., 2012) and coated for 1 hr at 37°C with 2.5 mg/ml Concanavalin A (Sigma-Aldrich), which binds polysaccharides expressed on the surface of the cell (Goldstein and So, 1965). Before the measurements, cantilevers were rinsed in PBS and cells were washed and probed in DMEM/Glutamax (Gibco) supplemented with 2% FBS (Life Technologies) and 1% penicillin-streptomycin solution (Life Technologies). Measurements were run at 37°C with 5% CO2 and samples were used no longer than 1 hr for data acquisition.

Tether force was measured at 0 velocity, which is linearly proportional to apparent membrane tension, assuming constant membrane bending rigidity (Hochmuth et al., 1996). In brief, approach velocity was set to 0.5 µm/s while contact force and contact time ranged between 100 and 200 pN and 100 ms to 10 s, respectively. The latter two parameters were experimentally tuned before every tether pulling attempt, aiming to reach a tradeoff between the maximization of the probability to extrude single tethers, and the reduction of experimental stress on the cells. The larger the contact time and force, the higher is the probability of formation of bonds between the molecules of Concanavalin A on the surface of the cantilever and the polysaccharides on the surface of the cell. On the other side, the lower those two parameters, the lower is the stress experienced by the cell during the contact with the cantilever. As a general trend, contact force and time must be increased over the course of the experiment, owing to the depletion of Concanavalin A from the cantilever (Krieg et al., 2008).

To ensure tether force measurement at 0 velocity, after contacting the cell surface, the cantilever was retracted for 10 µm at a velocity of 10 µm/s. The position was then kept constant for 30 s and tether force was recorded at the moment of tether breakage at a sampling rate of 2000 Hz. Each tether extrusion attempt lasted about 3 min and each cell was probed until three single tethers were successfully extruded or for a maximum of 10 min.

Force-time curves resulting from successful tether extrusions were analyzed using the JPK Data Processing Software. Tether force values from tethers extracted from the cell were then averaged, and each cell was accounted as a single data point.

Traction force measurements

Request a detailed protocol

Force measurements were conducted directly after seeding the cells on the sample and spreading was observed for 90 min on an inverted microscope (Nikon Ti-E2) with a Orca Flash 4.0 sCMOS camera (Hamamatsu) and a temperature control system set at 37°C. To avoid shaking the cells during stage movement, a POC-R2 sample holder in closed perfusion configuration was used and cells were seeded with a syringe right before image acquisition. The medium was supplemented with 20 mM of HEPES in order to buffer the pH during the experiment. Force measurements were performed using a method described previously (Tseng et al., 2011). In short, fluorescent beads were embedded in a polyacrylamide substrate with 20 kPa rigidity and images of those beads were taken during cell spreading. The first frame, before cells started attaching to the substrate, served as unstressed reference image. The displacement field analysis was done using a homemade algorithm based on the combination of particle image velocimetry and single-particle tracking. After correcting for experimental drift, bead images were divided into smaller subimages of 20.7 µm width. By cross correlating the subimages of the stressed and the unstressed state, mean displacement of the subimage can be measured. After correcting for this displacement, the window size is divided by 2 and the procedure is repeated twice. On the final subimages, single-particle tracking was performed to obtain a subpixel resolution displacement measurement. From the bead displacement measurements, a displacement field was then interpolated on a regular grid with 1.3 µm spacing. Cellular traction forces were calculated using Fourier transform traction cytometry with zero-order regularization (Sabass et al., 2008; Milloud et al., 2017), under the assumption that the substrate is a linear elastic half-space and considering only displacement and stress tangential to the substrate. To calculate the strain energy stored in the substrate, stress and displacement field were multiplied with each other and with the grid pixel area and then summed up over the whole cell. All calculations and image processing were performed with MATLAB.

Electron microscopy

Request a detailed protocol

Hela cells plated on fibronectin-coated glass coverslips for 30 min were disrupted by scanning the coverslip with rapid sonicator pulses in KHMgE buffer (70 mM KCl, 30 mM HEPES, 5 mM MgCl2, 3 mM EGTA, and pH 7.2). Paraformaldehyde 2%/glutaraldehyde 2%-fixed cells were further sequentially treated with 0.5% OsO4, 1% tannic acid, and 1% uranyl acetate prior to graded ethanol dehydration and Hexamethyldisilazane substitution (HMDS, Sigma-Aldrich). Dried samples were then rotary-shadowed with 2 nm of platinum and 5–8 nm of carbon using an ACE600 high vacuum metal coater (Leica Microsystems). Platinum replicas were floated off the glass with 5% hydrofluoric acid, washed several times by floatation on distilled water, and picked up on 200 mesh formvar/carbon-coated EM grids. The grids were mounted in a eucentric side-entry goniometer stage of a transmission electron microscope operated at 80 kV (Philips, model CM120) and images were recorded with a Morada digital camera (Olympus). Images were processed in Adobe Photoshop to adjust brightness and contrast and presented in inverted contrast.

Statistical analysis

Request a detailed protocol

Error bars represent standard deviation and in some cases standard error, which is specified in the legend. When applicable Shapiro-Wilk test was used to test for normality of data. Student’s t-test was chosen for statistical testing of normal distributed data, while Mann-Whitney U-test was performed on non-normal distributed data.

Appendix 1

Extended PLM for cell volume

The volume of the cell changes due to water flux driven by the difference of the osmotic and hydrostatic pressure. The volume dynamics is given by

(A1) dVdt=LpA((p(V)pext)(P(V)Pext))

where Lp is the cell’s hydraulic conductivity, is the total surface area, pext and Pext are the hydrostatic and osmotic pressure of the external medium, respectively. The volume dependence of the hydrostatic pV and osmotic pressure PV of the cell is explicitly shown. The volume dependence of hydrostatic pressure is through the force balance at the cell surface (Cadart et al., 2019) in which both tension and curvature are functions of cell size. In the dilute limit, the osmotic pressure is PV=kBTi=1NXi+i=1MYj , where T is the temperature, Xi is the concentration of solute species that are either actively or passively transported across the cell, and Yj is the concentration of impermeant solutes. The concentration of the trapped solute is Yj=yj/V-Vsolid where yj is the total number of molecules of type j in the cell, and Vsolid is the solid volume of the cell that is inaccessible to the solute molecules. The solid volume is essentially the sum of the volume taken up by the proteins and the DNA. The osmotic pressure of the external medium isPext=kBTi=1NX0i+j=1MY0j.

In Equation A1, the contribution of the hydrostatic pressure to the volume (or the volume change) is negligible in comparison to the osmotic pressure, and this even for small variation near the steady-state volume. Here is how we reach this conclusion: the concentration of the trapped osmolytes and its counterions is estimated to be about 100 mM, the osmotic pressure corresponding to this is about 105 Pa. Taking the cortical tension to be of order of 1 mN/m and cell size to be about 10 µm, we get the hydrostatic pressure difference of about 100 Pa. A 10% decrease in cell volume will increase the osmotic pressure of the trapped osmolytes by 104 Pa. For this osmotic pressure to be balanced by an increase in the hydrostatic pressure, the cortical tension would need to increase by a factor of 100, which we consider to be too large for this term to be relevant. Therefore, we find it reasonable to ignore the contribution of the hydrostatic pressure difference in the water flux equation. In this limit the volume dynamics is given by

(A2) dVdt=LpAkBT(i=1N(XiX0i)+j=1M(YjY0j)).

The main solutes that are transported across the cell membrane are sodium, potassium, and chloride (Kay, 2017). The sodium and potassium are actively transported through sodium-potassium (Na+/K+) pumps. The pump transports two potassium ions into the cell in exchange for three sodium ions transported out of the cell. This leads to the enrichment of potassium inside the cell and sodium outside the cell. Along with the Na+/K+ pump, there are various cotransporters and channels that passively transport the other ions - chloride, hydrogen ion, carbonates, and so on. A general expression of the ion transport including different cotransporters reads

(A3) d[(VVsolid)Xi]dt=jΛij(kBTlogXjX0j+zjΔΦ)+ Si,

where the first term is a passive flux of ion i due to the electrochemical potential difference of ion j (summed over all j), Λij is element of a symmetric matrix that couples the flux of the ions i to the chemical potential difference of ion j , ΔΦ is the electric potential energy difference between the inside and the outside of the cell ,zj is the electric charge of species j. The second term ,Si is the source term due to active pumping. The non-diagonal terms of the transport matrix Λ are due to the co-transporters (like NKCC, NHE, KCl, etc.) that couples the flux of ion j to that of i. In general, the permeability Λij and Si depend on the concentration of ions, electric potential difference, and membrane tension. Their volume dependence may be due to the transporters’ mechanosensitivity or due to feedback from signaling molecules that are sensitive to volume change.

The membrane potential is determined using the electroneutrality condition in the cell: i=1NziXi+j=1MzjYj=0. Using the electroneutrality condition and Equation A3, the potential is given by

(A4) ΔΦ=i=1Nzi(k=1NΛiklog(XkX0k)+Si)i=1Nk=1NΛikzizk.

Grouping the positive and negative ions and taking the trapped solute to be negatively charged, we get

(A5) i=1N|zi|Xi+j=1M|zj|Yj=k=1N+|zk|Xk+,

where N- is the number of negative ionic species and N+ is the number of positive ion species.

For monovalent ions, substituting Equation A5 in Equation A2 we get the steady-state volume as

(A6) V(t)Vsolid=j=1M(|zj|+1)yj(t)(βPext2i=1NXi(t)),

where Pext=kBTi=1NX0i+j=1MY0j , and β=1/kBT.

Fast timescales: cell as an osmometer

At the isotonic steady-state condition, the cell volume as given by Equation A6 is

(A7) VisoVsolid=j=1M(zj+1)yj(βPiso2i=1NXi),

where Pext=Piso is the osmotic pressure of the external medium at the isotonic condition. For hypotonic shock the medium is diluted leading to osmotic pressure Pext=rhypoPiso where rhypo is the dilution factor, and for hyperosmotic shock, PEG400 is added to the external medium leading to osmotic pressure Pext=rhyperPiso, where rhyper=Piso+PPEG/Piso. After the osmotic shock, the cell volume changes in response to the new extracellular osmolarity. There is a clear separation in the timescale between water transport and ion transport (Cadart et al., 2019). The ion concentration changes over minutes whereas the volume change due to water flux is on the timescale of seconds. This timescale separation divides the volume dynamics into ‘fast’ passive response, in which the water flows in and out of the cell with constant number of ions within the cell, and ‘slow’ response, in which the ions are transported across the cell. Over the timescale of seconds, the number of ions inside the cells is constant, that is, XiV-Vsolid=XiisoViso-Vsolid . Substituting this in Equation A2, we get

(A8) dVdt=PisoLpA(RVisoV(1R)VisoPextPiso),

where R=1-Vsolid/Viso. Thus, we see that the volume dynamics is well approximated by the Van't-Hoff relation with a fixed number of solutes in the cell. This equation at steady state gives the maximum (minimum) volume (Vm) after the fast hypoosmotic (hyperosmotic) shock. At steady state, we get the Ponder’s relation Vm/Viso=RPiso/Pext+(1R). From Equation A8, we get the rate of volume change just after the shock as

(A9) dVdt|t0=PisoLpAiso(1PextPiso).

Comparing Equation (A9) with the experimentally measured rate of volume increase just after the shock, we calculated Lp (Figure 4—figure supplement 1C). The cell volume over the minute’s timescale can be changed from this osmotic shock value by tuning the ion channels and pumps. We first consider the case when the ion transport does not change before and after the osmotic shock. The hypo-osmotic shock in the experiments is attained by dilution, we can see from Equation A4 that the membrane potential is constant, and from Equation A3, we see that the right hand side remains constant. This implies that after the fast increase in volume, the ions reach a new concentration, which is their steady state value, hence, there is no further ion flow at long timescale. However, in the experiments, it has been observed that the cell goes through a volume decrease. This necessarily implies a feedback mechanism for regulatory volume decrease.

In hyperosmotic shock condition, the membrane potential does change; hence, at long times, the volume of the cell increases to a value larger than the maximum decrease. This value is still less than the isotonic volume. The experiments show almost perfect adaptation, implying a regulatory volume increase.

Mechano-osmotic mechanism for cell volume regulation

There is a clear timescale separation between fluid flow, which is of the order of seconds, and spreading kinetics, which is of the order of minutes. Hence, over the timescale of spreading the cell is in osmotic balance with the external medium. The volume of the cell changes quasi-statically with the change in ion concentration according to Equation (A6). The change in volume can be due to change in the number of impermeant ions, or due to a change in the concentration of negatively charged ions. The rate change of the volume, obtained by taking the time derivative of Equation A6 and linearizing is

(A10) 1VisodVdt=1 VisodVsoliddt+Rj=1M(zj+1)yj(iso)j=1M(zj+1)dyjdt+2 R(βPiso2i=1NXi(iso))i=1NdXidt

where the first and second term on the right are due to a change in the number of impermeable molecules in the cell due to growth and the third term is due to the change in ion concentration due to change of ion transport rates. The volume-dependent feedback could affect either of the terms. We assume that the change in volume due to feedback is on the ion transport parameters, and the change in impermeable ions is only due to growth. We take the growth rate to be a rgrowth for all the trapped molecules, the growth in the number of trapped molecules also contributes to the growth in the volume of the solid fraction. For simplicity we take

(A11) dyjdt=rgrowth yj   and   dVsoliddt=rgrowthVsolid .

In the following we take rgrowth=0.05/hr, which gives a doubling time of 20 hr.

Our hypothesis is that the cell shape changes upon spreading generate a transient increase in membrane tension; this increase in tension leads to the activation of mechanosensitive ion channels and induces ion flux leading to volume change.

Slow spreading would then induce a lower transient tension increase that would lead to a smaller volume loss. We use a simplified model that includes the mechanosensitivity of ion transport. Assuming that the ion transport parameters vary quasi-statically, we use the following phenomenological expression for the change in ion concentration:

(A12) 2 R(βPiso2i=1NXi(iso))i=1NdXidt=α1γisodγmdt,

where the term on the right accounts for the change in ion concentration due to the feedback from membrane tension γm on mechanosensitive ion transporters. This defines the mechanosensitivity parameter α. Substituting Equation A11 and Equation A12 into Equation A10 we get

(A13) dδVViso=αdδγmγiso.

where δV=V-Viso1+rgrowtht and δγm=γ-γiso. Thus, we see that the change in volume is proportional to change in tension. For volume to decrease upon an increase of tension, the coefficient α should be negative. For a simplified model of transport of three ions—chloride, sodium, and potassium, we later show that α , in general, can take both positive and negative values. For the physiological value of the parameters, we find that α is negative if the increase in potassium permeability is much larger than that of sodium.

Effect of mechano-osmotic coupling of volume recovery upon osmotic shock

We can ask if a purely mechano-osmotic coupling could be sufficient to give a volume recovery upon osmotic shock. The answer to that question is that the mechano-osmotic coupling will affect the volume dynamics after the shock, but it is not sufficient to give a volume recovery at long timescales. This is because at long times the membrane tension goes back to its homeostatic value (γiso) and hence the ion transport parameters Λij and Si in Equation A3 also recover their pre-shock values. However, since the osmolarity of the external medium is different even at long time the volume will not go back to the value before the shock. The mechano-osmotic coupling will thus affect the volume recovery dynamics, but not the asymptotic value of the volume.

Sign of tension and volume coupling

We now solve for volume, electric potential, and α, for a minimal model which only accounts for three permeant ions: sodium, potassium, and chloride. The external medium only includes these ions and the cytoplasm also includes an impermeant negatively charged species. Following references Adar and Safran, 2020; Tosteson and Hoffman, 1960; Kay, 2017; Kay and Blaustein, 2019, we consider a simplified model of ion transport, that is, we ignore the contribution of co-transporters to the ion flux and take the matrix Λ in Equation A3 to be diagonal. The dynamics of ion transport through the cell membrane as given by Equation A3 now reads

(A14) d[(VVsolid)Na+]dt=ΛNa(kBTlogNa+Na0++ΔΦ)3Sp,
(A15) d[(VVsolid)K+]dt=ΛK(kBTlogK+K0++ΔΦ)+2Sp,
(A16) d[(VVsolid)Cl]dt=ΛCl(kBTlogClCl0ΔΦ),

where ΔΦ=ΦΦ0 is the electric potential energy difference between inside and outside of the cell, Sp is the effective permeability of the membrane to ion X due to ion channels, and Sp is the activity of the sodium-potassium pump – Na+/K+ ATPase, the prefactors—minus three and two—are due to the fact that the pump exchanges two potassium for three sodium. Cells pump sodium out of the cell, implying a positive value of Sp. Since the chloride ions are not actively transported there is no active pump contribution to the chloride flux. At steady state, Equation A14 to Equation A16 gives

(A17) Na+=Na0+eβ(ΔΦ3Sp/ΛNa),K+=K0+eβ(ΔΦ2Sp/ΛK),andCl=Cl0eβΔΦ

Note that if ΛNa,ΛK and Sp depend on the ion concentration and the electric potential difference, then Equation A17 is a set of implicit equations that need to be solved self-consistently.

In the following, we consider a single impermeant species (M=1 in Equation A2) of number in the cell Y, and effective charge z, which represents the proteins and the small impermeant charged molecules like phosphate ions. At steady state the osmotic balance reads

(A18) Na++K++Cl+yVVsolid=Na0++K0++Cl0.

The electroneutrality condition inside the cell and in the external medium is given by Na++K+=Clzy/V-Vsolid and Cl0-=Na0++K0+ , respectively. From the electroneutrality condition, Equation A17, and Equation A18 we get the potential difference across the cell membrane, and the cell volume. The volume thus obtained reads

(A19) V=Vsolid+(z+1)y2Cl0(1eβΔΦ).

For the volume to be finite we need ΔΦ<0, consistent with different experimental measurements. Substituting Equation A19 in the electroneutrality condition inside the cell, we obtain a quadratic equation for eβΔΦ . However, only one of the two roots leads to ΔΦ<0. The electric potential difference thus obtained reads

(A20) eβΔΦ=zz2(z21)((1δ)e3βSp/ΛNa+δ e2βSp/ΛK)(z1),

where δ=K0+/CI0. We can now evaluate the change in concentration of chloride ions when tension is changed. Identifying Xi with CI the l.h.s of Equation A12 reads

(A21) 1(Cl0Cliso)dCldt=Cl0(Cl0Cliso)deβΔΦdγmdγmdt,

where we have used Piso=2CI0. The electric potential difference depends on tension due to mechanosensitivity of the ion channels and pumps. Varying Equation A20 we get

(A22) deβΔΦ=(z+1)Naiso+2Cl0(zeβΔΦ(z1))3βSpΛNa((dΛNaΛNadSpSp)Kiso+Naiso+2ΛNa3ΛK(dΛKΛKdSpSp)).

For a small change in tension the channels and pumps change by a small value given by the relation

(A23) dΛNaΛNa=αNadγmγiso,dΛKΛK=αKdγmγiso,anddSpSp=αSdγmγiso.

We expect an increase in the channel values and decrease in the pump values due to an increase in tension, implying that the proportionality factors should be positive. Substituting this we get

(A24) α=R (z+1)Naiso+(αNa+αS)2 (Cl0Cliso-)(zeβΔΦz+1)3βSpΛNa(12ΛNa3ΛKKiso+Naiso+αK+αSαNa+αS).

Substituting the parameter values from Appendix 1—table 5 into Equation A24 we get

(A25) α0.0084(1.7αK+0.7αSαNa).

Thus we see that we expect to get α<0 if the values of αK,αNa, and αS are comparable. For slightly different parameters values, it is indeed possible that α>0. Fitting the experimentally measured cell volume and the spreading data, we get ξ=A0kα/γ1, taking α10-3 , this gives A0k/γ10 . In other words, a 10% increase of contact area leads to a doubling of tension in the elastic regime.

Membrane tension as a function of rate change of contact area

We model membrane as a Maxwell viscoelastic element (elastic at short time and viscous at long time), taking the change in membrane tension γm to be proportional to the contact area Ac , that is,

(A26) (1+τddt)(δγmγiso)=kτγisodAcdt,

where k is the elastic modulus and τ is the tension relaxation timescale that depends on various factors related to cortex organization as well on the membrane turnover. Note that k is not the elastic response of the lipid bilayer, it is an effective parameter that is related to the cells ability to access its membrane reservoirs upon stretching. Substituting Equation A26 in Equation A13, we get

(A27) 1VisoτdδVdt+1VisoδV=ξτ1A0dAcdt,

where ξ=-A0kα/γ0 . For a time much less than the volume relaxation timescale τ the rate of volume change is proportional to the rate of spreading

(A28) 1VisodδVdt=ξA0dAcdt.

As observed in the experiments, faster the spreading rate larger the initial rate of volume loss. We take the time-series of the contact area as input in Equation A27. We first fit the contact area times series to the following equation Act=A01-exp-t/τa, and thus obtain the best fit values of A0 and τa . Solving Equation A27 with this expression for Act, we get

(A29) VViso=1+rgrowtht+ξτet/τaet/τττa.

We then obtain the best fit value of τ and ξ by numerically fitting the volume dynamics to the measured volume time series using the ‘NonlinearModelFit’ function in Mathematica.

Fast and slow-spreading control cells

As seen from the estimates the drug treatment affects multiple parameters. To confirm the relation between the spreading speed and volume loss we sorted the control cells (n=127) into three equal groups based on their average spreading speed in the first 10 min. We find that, indeed, the fast-spreading cells lose more volume than the slow-spreading cell. We fit the three groups—slow, intermediate, and fast-spreading cells—to the model and obtain the best-fit parameters. We find that the parameters characterizing the mechano-osmotic feedback (ξ andτ) are similar for the three classes. Compared to the cells spreading at intermediate speed, the rate of volume loss of the slowest spreading cell was 50% less than that of the fast-spreading cell was 50% more. The average of all the control cells was close to the group of cells with intermediate spreading speed and the fast-spreading cells behaved similarly to the GdCl3 treated cells. The value of the parameters is listed in Appendix 1—table 1.

Appendix 1—table 1
List of fitted parameter values for the control cells are grouped based on spreading speed when the tension depends on change in contact area.
Condition/parametersA0 (µm2)τa(min)τ(min)ξ
Slow spreading62063.398.50.21
Moderate spreading48825.664.20.16
Fast spreading65617.299.40.16

For different drug treatments

The parameter τ varies over a wide range of values, whereas ξ is about one tenth for all cases except Y-27 and about half for the Y-27 treated cells. For fast-spreading cells treated with Y-27 and GdCl3, the spreading rate A0/τa is about 1.8 times that of control. However, the initial rate of volume loss from the model is about four times for Y-27 treated cells and about 1.5 times for the GdCl3 treated cells. The parameter ξ/A0 increases by a factor of 1.5 for the Y27 treated cells and decreases by a factor of 0.75 for the GdCl3 treated cells. For EIPA treated cells, the spreading rate is the same as control but the parameter ξ decreases by 20% as seen by a smaller volume loss. The initial rate of volume loss and rate of spreading is similar for the Lat A and CK-666 treated cells.

The volume recovery timescale is quite variable for the different drug treatments. Over the measured timescale of an hour, the volume recovery of the GdCl3 and CK-666 treated cells is mainly due to growth. The value of the parameters is listed in Appendix 1—table 2.

Appendix 1—table 2
List of fitted parameter values for different drug treatment when the tension depends on change in contact area.
Condition/parametersA0 (µm2)τa(min)τ(min)ξ
control47925.953.50.16
Y-2783324.416.40.67
EIPA54933.730.70.18
Lat A152221213.80.88
GdCl362818.3171.20.17
CK-66628514.90.09
CK-666+Y-2773825.6360.1
EIPA+Y-271035290.02

Membrane tension as a function of total area

The membrane tension dynamics is given by following Maxwell viscoelastic model:

(A30) (1+τddt)δγmγiso=kγisoτAtdAtdt,

where γiso is the homeostatic value of membrane tension, τ is the relaxation timescale, k is the membrane elasticity, and At is the total cell area. Substituting Equation A13 in Equation A30 the equation for volume dynamics reads

(A31) (1+τddt)δVViso=ξτAtdAtdt,

where ξ=-kα/γiso . To compute the total cell area we need to make assumptions about the shape of the cell. At early times, the shape can be well approximated by a spherical cap. However, at a later stage of spreading it is not a reasonable assumption. Since we do not know more about the three dimensional shape of the cell, for simplicity, we model the cell as a spherical cap through the duration of spreading. With this assumption, we can calculate the total area in terms of the instantaneous volume and the contact area. The volume and the total surface area of the spherical cap are given by

(A32) V=16h(3Ac+πh2)andAt=(2Ac+πh2),

where h is the height of the cell and Ac is the contact area. Eliminating the height from Equation A32, we get the total area in terms of the volume and contact area, which reads

(A33) At=(2Ac+π(Ac+(3πV+Ac3+9πV2)2/3π(3πV+Ac3+9πV2)1/3)2).

Taking the time derivative of Equation A33 and using Equation A32, we get the following expression for rate change of the total area as function of rate change of the volume and the contact area:

(A34) At˙=24πVV˙(At+Ac)(AtAc)+2AcAc˙AtAc.

Substituting Equation A33 into Equation A31, we get

(A35) (1+τddt)VViso=1+(τ+t)gξτ(f1(V,Ac)dVdt+f2(V,Ac)dAcdt),

where f1V,Ac=24πV/AtAt2-Ac2 and f2V,Ac=2Ac/AtAt-Ac .

Rearranging the terms we get

(A36) (1+τeffddt)VViso=1+(τ+t)gξτf2(V,Ac)dAcdt,

where τeffτ1+ξf1V,AcViso . Thus, we see that the volume relaxation timescale is normalized by the volume dependent term in total area, and this effective timescale is larger than the bare tension relaxation timescale.

The effective tension dynamics obtained by using Equation A13 in Equation A36 is

(A37) (1+τeffddt)δγγiso=kττeffγisoτeff(rgrowth+f2(V,Ac)dAcdt),

We now fit Equation A36 with the experimentally measured volume to obtain the best-fit parameters. For this, we fit the cell spreading data with an exponentially saturating function of the form

(A38) Ac(t)=A0(1A0Ac(0)A0et/ta),

where Ac0 is equal to the initial contact area which we obtain from the data. We solve Equation A36 numerically, using the NDSolve function in Mathematica, for a set of parameter values ξ and τ. We then select the parameters which minimize the error between the numerically calculated volume and the experimentally measured volume using L2 norm.

The best-fit parameter values for control cells grouped in three groups based on spreading speeds are listed in Appendix 1—table 3 and that for cells treated with different drugs are listed in Appendix 1—table 4. We see that the value of τ is much smaller compared to the case in which we use contact area as a proxy for total area. This is due to the fact that the volume relaxation time is normalized by the dependence of the total area of the spherical cap on the volume. We also find that parameter ξ is an order of magnitude larger compared to the fit using contact area. Although the parameters for the cells in three control groups (Appendix 1—table 3) vary significantly, this variation cannot explain the difference in volume loss. If we use the same spreading speed for the three sets of parameters, the resulting volume loss is actually maximum for the parameter corresponding to the slowest spreading cell. This implies that the observed difference in volume loss must be attributed to the difference in spreading speed itself (Appendix 1—figure 1).

Appendix 1—figure 1
Model Figure 1.

Fits for volume taking the contact area of the fast-spreading cells for the three sets of parameters. Model predictions for the cell volume, based on fits from data, for slow (green), moderate (cyan), and fast (pink) spreading.

We now compare the change in tension for two cases, one with finite volume tension coupling, that is, α0, and compare it with the case α=0, which implies ξ=0. For ξ=0, the volume is given byV=Viso1+rgrowtht.

Substituting this volume in Equation (A30), we can compute the difference between the two tensions. We find that during the initial spreading the change in tension when there is volume loss due to spreading is always lower than the case when there is not volume loss. In the later part of the spreading, this is not the case. However, since the spherical cap model is more reasonable at the start of the spreading rather than later. This supports the hypothesis that the functional role of the volume loss may be to prevent rapid increase in tension due to fast spreading.

Fast confinement

From Equation A1, we see that, within the PLM framework, for a fixed external medium, the cell volume can change either due to change in hydrostatic or osmotic pressure. The Membrane rupture tension is 20mN/m, for radius of curvature of about 5 this gives ΔP4×103Pa. The change in volume due to this pressure increase as given by Equation A1 for constant ion concentration is

(A39) VisoV=ΔPPiso+1.

For external osmolarity of 300mM,,iso7×105Pa, which gives ΔP/Piso0.5×10-2. Thus, we see that even at the rupture tension the hydrostatic pressure can change the volume only by about 1%. Hence, the volume loss upon fast confinement that is of the order of 10% cannot be explained by just the increase of the hydrostatic pressure due to compression.

Within the framework of PLM, the volume change of the order of 10% can only be due to change in the osmolarity of the cell, which requires transport of ions. For ions transport to take place at timescales of milliseconds, the rates need to increase by 4 orders of magnitude. Such increase can be easily attained due to pore formation. However, formation of small pores that allow the ions to leak through but does not discriminate between the different ions will lead to an increase rather than a decrease in volume. This is due to the fact that the concentration of ions outside the cell is larger than that inside, hence once the pores open the ion flux is into the cell and the water flux follows the ions flux.

Volume fluctuation during cell migration

The cell area fluctuates as the cell passes through the collagen matrix.

We take the total area to be At=A0sinωt , where ω~vcell/lmesh . As before assuming a viscoelastic model for tension driven by the area At , we can compute the volume fluctuation due to mechano-osmotic coupling with membrane tension. The volume dynamics is given by

(A40) (1+τddt)VViso=ξτAisodAtdt

From Equation A40, the standard deviation of the volume is

(A41) V~2Viso=ξτω1+ω2τ2A0Aiso

If the distance cell moves in time τ is much smaller than the mesh size vcellτlmesh then the standard deviation of volume increases linearly with the cell velocity. As the speed increases, standard deviation of the volume saturates to the value ξA0/Aiso . We fit Equation A41 to the experimentally measured values (Figure 7K). Taking the standard deviation of the volume at zero velocity to be 0.01 gives the best fit parameter values to be τ/lmesh=1.8min/µm and ξA0/Aiso~0.04. For a mesh size of about 5 µm and area change of the order of few percents we get τ~10min, and ξ~1, which is in the same range as that obtained when fitting the volume change upon spreading (see Appendix 1—table 4).

Appendix 1—table 3
List of fitted parameter values for the control cells are grouped based on spreading speed when the tension depends on change in total surface area.
Condition/parametersA0 (µm2)τa(min)τ(min)ξ
Slow spreading62063.352.5
Moderate spreading48825.6130.9
Fast spreading65617.2510.3
Appendix 1—table 4
List of fitted parameter values for different drug treatment when the tension depends on change in total surface area.
Condition/parametersA0 (µm2)τa(min)τ(min)ξ
Control47925.9121
Y-2783324.48.51.3
EIPA54933.78.50.8
Lat A152221223.4
GdCl362818.334.50.4
CK-66628514.929.50.9
CK-666+Y-2773825.6180.2
EIPA+Y-271,035294.50.1
Appendix 1—table 5
List of parameter values.
ParameterDefinitionEstimate
Concentration of chloride ions in the medium150 mM (Kay, 2017)
K0+Concentration of potassium ions in the medium3 mM (Kay, 2017)
Na0+Concentration of sodium ions in the medium147 mM (Kay, 2017)
ClisoConcentration of chloride ions in the cell4 mM (Equation (A17))
Kiso+Concentration of potassium ions in the cell130 mM (Equation (A17))
Naiso+Concentration of sodium ions in the cell5 mM (Equation (A17))
YConcentration of trapped particles in the cell30 mM (osmotic balance)
zYConcentration of trapped charges in the cell131 mM (electro-neutrality)
ΔΦPotential difference across the plasma membrane–90 mV (Equation (A5))
kBTTemperature (25°C)4.1*10–21 J
ΛNaPermeability of Na+2*10–8 moles/(m2*s) (Kay, 2017)
ΛKPermeability of K+6*10–7 moles/(m2*s) (Kay, 2017)
ΛClPermeability of CI4*10–7 moles/(m2*s) (Kay, 2017)
βSpNaK Atpase pumping rate1.4*10–7 moles/(m2*s)
LpHydraulic conductivity10–(12–13) m/(Pa*s) (Equation (A9))
RRatio of osmotically active volume to total volume0.7

Data availability

All data generated or analysed during this study are included in the manuscript and supporting file; all the raw analysed data shown in the figure panels in the article are available in the accompanying Source Data files.

References

    1. Burridge K
    2. Monaghan-Benson E
    3. Graham DM
    (2019) Mechanotransduction: from the cell surface to the nucleus via RhoA
    Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences 374:20180229.
    https://doi.org/10.1098/rstb.2018.0229
    1. Ting-Beall HP
    2. Needham D
    3. Hochmuth RM
    (1993)
    Volume and osmotic properties of human neutrophils
    Blood 81:2774–2780.

Decision letter

  1. Felix Campelo
    Reviewing Editor; Institute of Photonic Sciences, Spain
  2. Anna Akhmanova
    Senior Editor; Utrecht University, Netherlands
  3. Sean X Sun
    Reviewer; Johns Hopkins University, United States

Our editorial process produces two outputs: i) public reviews designed to be posted alongside the preprint for the benefit of readers; ii) feedback on the manuscript for the authors, including requests for revisions, shown below. We also include an acceptance summary that explains what the editors found interesting or important about the work.

Decision letter after peer review:

Thank you for submitting your article "A mechano-osmotic feedback couples cell volume to the rate of cell deformation" for consideration by eLife. Your article has been reviewed by 3 peer reviewers, and the evaluation has been overseen by a Reviewing Editor and Anna Akhmanova as the Senior Editor. The following individual involved in review of your submission has agreed to reveal their identity: Sean X. Sun (Reviewer #1).

The reviewers have discussed their reviews with one another, and the Reviewing Editor has drafted this to help you prepare a revised submission.

Essential revisions:

- Discuss the accuracy/validity of the FXm method (Rev. 1, point 1).

- Is there an effect on acto-myosin contraction? Include the effect of hydraulic pressure in the model or further discuss why it can be neglected (Rev. 1, point 2 and recommendation 2,3). Related to this, consider also tje effect of myosin II mechanosensitive activation (Rev. 2 point 5).

- Similarly, the constant term in the model for the active pumping needs to be either reformulated or the further justified (Rev. 1, point 3).

- How is membrane tension at early times after seeding (e.g. 10min condition)? (Rev. 2, point 1).

- Somewhat related, it is not clear the role membrane tension plays in volume regulation, in particular with respect to the comparison between drug treatments and the theoretical model (Figure 6I,J) (Rev. 3 main comment). Would it be possible to monitor tension dynamics as suggested by Rev. 3?

- Add details and discuss further on the pulling experiments (Rev. 2, point 2 + Rev. 1 recommendation 4).

- Improve on the statistical analyses of the results (Rev. 2 point 9).

- I agree this is a very complex paper with lots of data. Probably the authors should consider re-writing some parts of the paper more clearly and expose better the different parts of the paper (Rev. 1, point 4 and recommendation 1). Although splitting the paper in two is of course a possibility, I think I would suggest that the distinction between the two scenarios (osmotic shock vs. spreading) is made more clear by re-sectioning and maybe adapting the model for each case, following the reviewer's suggestions (rev. 1 point 4 and also comment 3 of rev. 2).

- Discuss on the role (if any) of the initial cell volume in the observed effects (e.g. volume regulated by spreading speed) (Rev. 2, point 2).

- Direct measurements on calcium-mediated activation of myosin II (Rev. 2 point 6).

- Related to NHE1, is it also involved in volume regulation under osmotic stress? (Rev. 2 point 7).

- Discuss on alternative Arp2/3 functions that can lead to volume regulation (Rev. 2 point 8).

Reviewer #1 (Recommendations for the authors):

1) It has to be said that the paper is quite complex to read. There is almost too much data. 3 different scenarios, spreading, osmotic shock, and confinement, are considered. While I appreciate that these phenomena are inter-related, as it is it's difficult to read. Each scenario can be explored in more detail, with models specifically addressing the scenarios. I understand this would require extensive revision, so I leave it up to author's discretion.

Consider splitting the paper and concentrate on 1 or 2 phenomena only. Personally, I think the dynamics during cell spreading and osmotic shock are significantly different. Osmotic shock involves a very large change in ion concentration, therefore might be more appropriately described by the model. During spreading and confinement, there is a small volume change, and I think the correct approach is to consider active contraction and ion dynamics together.

2) It could be more appropriate to model spreading and confinement using a 1-solute phenomenological model, but include myosin contraction.

3) Some comments could be more rigorous. For example, the tether force is not just a function of membrane tension, but also active contraction. Comments about 1-solute model, etc.

4) What is the significance of measuring tether force 4hrs vs. 30min?

Reviewer #2 (Recommendations for the authors):

The authors find a correlation between dV/dt and dA/dt in spreading cells and argue that an increase in tension is associated with cell volume loss. It is however not fully clear how a fast increase in tension is build up during the initial phases of cell spreading to explain the highest change of dV/dt at short time scales.

Tension values were also measured via AFM at 30 min after spreading when dV/dt and dA/dt are zero. The authors should clarify how this time point is informative to argue about tension driving volume change. In this context, it remains to be tested how tension evolves at earlier time points during spreading. Given also the higher tension in cells for specific treatments at 30 min (e.g. Y-27632), why would this not lead to a longer period of cell volume decrease compared to control cells considering the authors' argument of tension homeostasis in cells? In this context, the authors could also extend their comments on the results of panel 6I in the main text and how far they match with the measured apparent membrane tension.

The presentation of AFM measurements and data would benefit from a more detailed description of the methodology. Tether pulling experiments were indicated at 30min/4h but given the time it takes to perform a single experiment and probe multiple cells the authors should comment on the specific time duration of experiments to facilitate the interpretation of results. In addition, parameters as tether pulling positions (cell body versus protrusions) and tether pulling statistics would enable to assess the presented data more rigorously. Can the authors further comment on the highly variable contact times used for AFM tether pulling experiments, is it associated with specific experimental conditions? Were drug treatments further performed 30min prior to measurements or was the incubation performed for 4h, to clarify long-term effects of drugs and potential cell adaptations.

Authors define categories of fast/intermediate and slow spreading cells. Can the authors clarify if initial cell volumes is relevant to consider and different in these categories? Does initial cell volume also change with drug treatments (e.g. Y-27632 treatment, a condition in which the cells spread faster)? The authors could further comment in the main text on potential parameters regulating spreading rate variability in cells.

The authors use confinement and osmotic shocks to test model assumptions and to study cell volume changes at short time scales. However, the authors so far do not consider applying the modified PLM to volume adaption in osmotic shock conditions or potential long-term volume changes in confinement. Can the mechano-sensitive PLM also capture volume adaption in osmotic shock conditions or what are the potential limitations? In addition, do cells in confinement show a volume adaptation over time that would follow similar principles as in spreading cells or osmotic shock experiments?

The authors argue for a role of myosin II in counteracting volume loss in spreading cells. Would this effect also be relevant in confinement as the authors previously showed a mechanosensitive activation of myosin II in confined cells? Does volume loss occur similarly when myosin II activity is inhibited? Do the authors expect a similar role of myosin II in regulating volume change in confinement and volume adaptation in osmotic shock conditions and spreading cells?

The authors explain the Gadolinium effect by referring to the fact that calcium increase might activate myosin II. Direct measurements such as MLCK treatment or different calcium interference approaches would allow to support this argument directly.

Does NHE1 play a similar role during adaptation to osmotic stress conditions as observed in spreading cells? It would further be interesting to show whether the function of NHE1 for mechano-sensitive volume regulation is observed in the other cell types tested in this study.

Arp2/3 is identified as a molecular regulator of cell volume in spreading cells and the authors speculate on its role in regulating ion channel/pump activity. The authors could further include a discussion on alternative functions of Arp2/3 such as regulating endo/exocytosis relevant for cell volume regulation.

The discussion could further benefit from adding further arguments on the physiological relevance of cell volume regulation, considering the observed range of volume change and the focus on the specific mechanical stress conditions studied in this work.

Statistical testing, the authors present N=1 in multiple panels of the presented data and should re-check the statistical significance of these results. Including p-values directly in the graphs would ease the reading of figure panels.

Propidium iodide staining could be used to monitor plasma membrane rupture in confinement, considering the high leakage/volume loss in confined cells. It also appears that the authors did not comment on the high increase in surface area in Figure 5H in the main text.

LatA treatments are often used at concentrations which only partially depolymerize the cell cortex but concentrations vary and could be more clearly commented in the main text. Would higher LatA concentrations increase cell spreading in confinement and alter the cell volume change?

Related to Figure 1G, the authors state "causing a transient density increase". Panel S1G however indicates an increased density with a stable plateau, can the authors clarify this statement?

The authors state that "with a loss of volume reaching up to 20% for fast-spreading cells" and "actin-dependent loss of volume (up to 30%)". The max volume loss of cells is not clearly visible in the current graphs and appears exaggerated in the text. The authors could include individual data points or change to a different data illustration to visualize the range of cell volume loss better.

The authors could list the specific cell types tested in the study in various statements in the main text that appear rather generic otherwise.

How does tau_a versus tau in the table relate to tau_eff and tau in the main text? Providing a clearer introduction and overview of model parameters would support the readability.

Comments on figure panels:

Figure 1A, cells are not well visible in this figure panel and a higher magnification or zoom in could be used.

Figure 1B, indicate channels in legend and check for consistency in figure panel.

Figure 1C, for the micropattern data indicate spreading time relevant for volume calculation. The authors also argue that they use micropatterns to reduce spreading area to a defined size compared to control conditions. Comparing Figure 1B and 1A, the cell areas are however similar and control cell spreading area appears to be higher in panel 1B versus 1A?

Figure 1I, J, indicate in the legend at which time point dV/dt was measured. How was the scale bar in panel 1I calculated?

Figure 2C, there appears to be a gap in the data of Y-27632?

Figure 3A, the explanation of the schematic of the PLM could be explained in more detail to ease the reading.

Figure 3D vs 4B, why is the Ponders equation in panel 3D independent of actin but not in panel 4B? LatA data in panel 4B appear to be not commented in main text.

Figure 4B, check sign of dV/dt compared to other figure panels as it appears neg/pos for hypo/hyper conditions when cells swell/shrink.

Figure 4C, over which duration was adaptation measured?

Figure 5E, the legend indicates that a double expression line of myosin II/actin was used but only one channel is shown (membrane?). The visualization of myosin II and actin would be informative to include.

Figure 6K, indicating measured tension values in Figure 6K or matching the color code between panel 6J and 6K would ease the reading of figures.

Figure S2G Why is the correlation before and after lamellopodia changing slope, appears to be not commented and the figure is not cited in the main text.

Figure S2F, what do the lines indicate?

Reviewer #3 (Recommendations for the authors):

It would be good to insert figure S4H in main figure, as it clearly how speed of deformation affects the volume change.

The absence of correlation between volume and spreading area at a given spreading rate is interesting and new. It is probably worth discuss the authors' results in light of Ref 11, that conclude on a weak correlation between spread area and volume, and ref 9 that observe a strong correlation between area/volume.

Fig1F

Could the authors consider including the spreading area change in time on PLL-PEG?

Fig1D-H

Color code is not consistent and make the figure harder to read.

Page 5, first paragraph

Please explain which criteria you used to pool the cells according to their spreading speed: do you have 3 clear populations, well separated by spreading speed, or you had to apply cutoffs sort cells in different populations?

Could the authors indicate the three corresponding populations ('slow spreading', 'moderate spreading', 'fast spreading') on Figure 1I? Consider clarifying the color code.

Can the authors plot Figure 1J for each population ('slow spreading', 'moderate spreading', 'fast spreading')? Is the same tendency still observed?

Page 5, second paragraph

In RPE-1, dynamics of spreading area (Figure S2D) and volume (Figure S2E right) for control cell and under Y-27 are very similar, while it is not the case in HeLa cells. Can the authors comment on it?

Figure 2B-C

Axes do not have the same minima/maxima so complicated to read and to compare with Figure 1H.

Figure 2D

Axis do not have the same minima/maxima. It would be easier to read and compare with Figure 1I if they had the same range. Are the linear fits in fig2D and fig1I giving the same results? They seemed to be different.

Page 5, third paragraph

The results of this part are very similar to the paper cited above (https://doi.org/10.1101/2021.01.22.427801), probably worth mentioning. In the main text, authors should specify that you used PEG400 to apply hyperosmotic shocks. Do the authors find differences if they apply osmotic shocks with salts or adding non-metabolized sugars such as sucrose or sorbitol?

Page 6, second paragraph

Under Y27 treatment, cell spread faster and their membrane tension measured using AFM is 30% higher than control 30 min after seeding. As LatA treated cells spread slower, wouldn't the authors expect a smaller membrane tension? Dai and Sheetz-1999 looked at the effect of cytochalasin D and indeed observed a reduced membrane tension.

Page 8

What are the differences, beside the timescale, between Figure 3C and Figure 4A? It would be worth mentioning the way these different were acquired in the figure legend.

Page 9 last paragraph + page 10 first paragraph

Is intensity of CellMask a proper tool to extrapolate on membrane folds and reservoirs (Figure 5E)? The spread around the contour line seems larger for 7.5um than for no confinement for control cells, and the reverse for Lat A treated cells. I am not sure how these intensity plot can be quantitatively linked to the size of the membrane reservoir.

Also, in the before to last paragraph in page 9, I kindly disagree with the statement that "membrane tension could not arise from a limitation in membrane availability". Isn't it that unfolding the membrane partially, not reaching the full size of the membrane reservoir, would still change membrane tension, because of opposing cytoskeletal and adhesive forces trying to fold the membrane?

Figure 5F

LatA treated cell keep a constant volume after confinement while control cell decrease volume. Since you have a condition (Y27 treated cell) that have opposing effect of LatA on spreading area and volume, I wonder what would be the effect of confinement on Y27 treated cells. Would cells 'overreact?'

Figure 6 A-H

Y-axis scale is different for each graph, would be better to fit one range and keep it so that the reader can compare from one to another.

Figure 6J

Do you affect absolute cell volume when adding those drugs?

[Editors' note: further revisions were suggested prior to acceptance, as described below.]

Thank you for submitting your article "A mechano-osmotic feedback couples cell volume to the rate of cell deformation" for consideration by eLife. Your article has been reviewed by 3 peer reviewers, and the evaluation has been overseen by a Reviewing Editor and Anna Akhmanova as the Senior Editor. The following individual involved in the review of your submission has agreed to reveal their identity: Sean X. Sun (Reviewer #1).

The reviewers have discussed their reviews with one another, and the Reviewing Editor has drafted this to help you prepare a revised submission.

Essential revisions:

We all agreed that the revised version of the manuscript is much improved and clearly presents a very interesting and compelling story, well done!

Before the formal acceptance of the manuscript, one of the reviewers (rev. #1) made two final comments regarding (1) the contractile stress's role in volume control; and (2) the relative importance of membrane tension and voltage for the effects presented here. I think it might be a good idea to consider these comments and decide whether you think you can adapt the discussion part of the text so your points are clearly laid down.

Reviewer #1 (Recommendations for the authors):

The authors have addressed my comments and the manuscript is much improved. I can recommend publication. I think there is generally no disagreement in our thinking of contractility and cell volume regulation. It is true that contractile force, or pressure, alone cannot really change the volume very much. This was already shown in Yellin et al., Biophys. J. 2018, and Li and Sun, Frontiers Cell Dev. Bio., 2022, and I agree with the authors' general thinking. However, the discussion in Pg. 7, I think, can be improved. Right now, it sounds like contractile stress is being compared with the total osmotic pressure of impermeable solutes inside the cell. But in reality, the contractile stress should be compared with the total osmotic pressure difference between inside and outside of the cell, which is much smaller than 100mM. During 10% volume change during cell spreading, I think since it's happening slowly (order of minutes), the ion flux and contractile stress should change simultaneously, possibly with contractile stress changing slowly, while ion flux happening very fast. But, the contractile stress must change, since at minimum, the cell is changing shape, and therefore cortical stress must change to maintain relative force balance. I'm still not completely sure whether is ok to neglect it, since the steady-state result of the proposed model is zero osmotic pressure difference between inside and outside, but in reality, there is about 1kPa of hydraulic pressure difference and 1kPa osmotic pressure (~0.5mM) difference. This difference may be small enough to be negligible in the current context. But it may not be small in other contexts, such as during cell motility.

The other issue I hope the authors can clarify in the discussion is the equal importance of membrane voltage and membrane tension in this story. Ion exchangers such as NHE and NaK are affected by both. The current dependence of exchanger flux with respect to voltage is based on just passing an ion through a constant electric field. This is generally not true for ion exchangers. I think most studies of ion exchangers show a voltage gating effect. This, after all, is the basis of action potentials. Therefore, these important facts should be clearly stated.

Reviewer #2 (Recommendations for the authors):

The authors have made an effort to provide a detailed response to the reviewer's comments and the manuscript is of high quality to be published in Elife. The work presents a relevant step forward to understand the coupling of cell mechanics during cell spreading and volume regulation and opens up interesting questions to be addressed in future studies.

Reviewer #3 (Recommendations for the authors):

The authors have completely addressed all my points, and I enthusiastically recommend publication.

https://doi.org/10.7554/eLife.72381.sa1

Author response

Essential revisions:

- Discuss the accuracy/validity of the FXm method (Rev. 1, point 1).

This has been done. The result is that, from one experiment to another, the absolute measure of the average volume of a population of cells varies by +/- 10% as mentioned by the referee. But in this article, most of the measures are temporal measures of the same single cell. In this case, the accuracy of the measure is +/- 1%, as we now show experimentally. See more details in the answer to Referee 1, point 1 below.

- Is there an effect on acto-myosin contraction? Include the effect of hydraulic pressure in the model or further discuss why it can be neglected (Rev. 1, point 2 and recommendation 2,3). Related to this, consider also tje effect of myosin II mechanosensitive activation (Rev. 2 point 5).

This has been done. The experimental results show that acto-myosin contraction has no direct effect on cell volume. This is also what the classical theory predicts (and thus others before us, see 1). The reason comes from the fact that while ions represent the most abundant osmolytes in the cell (about 200 mM), the amount of nonionic osmolytes (about 100 mM) cannot be neglected. The force needed to ‘concentrate’ these nonionic osmolytes even by a few percent is much larger than what the acto-myosin cortex can produce. The only way to change the volume of the cell by a few % or more is by changing the ion content, as predicted by the classical pump and leak model (or change the nonionic content, but that is even much slower and is rather on the timescale of cell growth). This means that the effect of myosin contractility on the volume of spreading cells is only indirect, via the modulation of the spreading speed. See also the more detailed answers to the referees below.

- Similarly, the constant term in the model for the active pumping needs to be either reformulated or the further justified (Rev. 1, point 3).

This has been addressed. See the detailed discussion on this modeling aspect in the answer to the referee.

- How is membrane tension at early times after seeding (e.g. 10min condition)? (Rev. 2, point 1).

This point has been addressed. It is not possible to measure membrane tension at early times for two reasons: the measure takes time, and the cells need to be adhered enough (they should not move). We thus instead provide a measure of the membrane tension for cells which are adhered but not spread, with a spreading area corresponding to the early times after seeding, using small adhesive micropatterns. See more details in the answer to the referee.

- Somewhat related, it is not clear the role membrane tension plays in volume regulation, in particular with respect to the comparison between drug treatments and the theoretical model (Figure 6I,J) (Rev. 3 main comment). Would it be possible to monitor tension dynamics as suggested by Rev. 3?

This point was mostly due to a confusing labelling of the axis of Figure 6I, which made it hard to understand what the graph shows, despite details given in the figure legend. We changed the labelling. What this graph shows is the difference in tension, as predicted by the model, comparing the case with or without the mecano-osmotic coupling effect, for each of the condition, using the fits from the experimental data. The largest difference is for the Y-27 treatment, because in this case the cells are spreading fastest and not having the mechano-osmotic feedback would lead to a large increase in tension. On the other hand, for conditions in which the mechano-osmotic feedback is already altered (EIPA treatment for example), the difference of removing it in the model is small. Having that in mind, it is then clear that the predictions of the model from Figure 6I actually correspond perfectly to the experimental observations for the membrane tension – affecting the mechano-osmotic coupling in Y-27 treated cells (for example the Y-27 plus EIPA or the 7-27 plus CK 666 treatments) leads to a further increase in the membrane tension assessed from tether pulling, mostly during the fast spreading phase. We understand that it is not a very easy point to understand, but we hope that it is now more clear. We also modified the text describing these results (p. 12).

As for the suggestion of Referee 3 to use Flipper probes, we have tried and provide the results that we obtained, but we failed to get good enough data to lead to significant results. We also discussed with the lab which developed these probes (Aurélien Roux) and Aurélien confirmed that such experiments are hard to do (live cell imaging with the probe) and would require us to perform a series of experiments in his lab, which was not possible to organize given the current situation and the short time for the revisions. Moreover, comparing strictly the results from tether pulling and the Flipper probe is a project per se. These experiments will be performed in the future, because they are indeed very interesting, but it is out of the scope of this article.

- Add details and discuss further on the pulling experiments (Rev. 2, point 2 + Rev. 1 recommendation 4).

This point has been addressed. See the details in the answer to the referees below.

- Improve on the statistical analyses of the results (Rev. 2 point 9).

This point has been addressed.

- I agree this is a very complex paper with lots of data. Probably the authors should consider re-writing some parts of the paper more clearly and expose better the different parts of the paper (Rev. 1, point 4 and recommendation 1). Although splitting the paper in two is of course a possibility, I think I would suggest that the distinction between the two scenarios (osmotic shock vs. spreading) is made more clear by re-sectioning and maybe adapting the model for each case, following the reviewer's suggestions (rev. 1 point 4 and also comment 3 of rev. 2).

This has been addressed be rewriting extensively the introduction and part of the discussion of the article. We feel that the main issue was that we did not make it clear enough that our article is solely centered on understanding the volume modulation during cell spreading and that the central result is that volume modulation depends on spreading speed. The other aspects, osmotic shocks and confinement, were introduced to test some aspects of the model, and to characterize the (almost) purely osmotic and (almost) purely mechanical volume responses of the cell, while spreading corresponds to a regime where these two aspects overlap due to the modulation of ion fluxes by cell membrane tension.

- Discuss on the role (if any) of the initial cell volume in the observed effects (e.g. volume regulated by spreading speed) (Rev. 2, point 2).

This has been addressed by showing that there is no effect of the initial volume of the cells, and that the various drug treatments do not affect the initial volume distributions. See more details in the answer to the referee (in fact Referee 2 point 3, as well as Referee 3 point 21).

- Direct measurements on calcium-mediated activation of myosin II (Rev. 2 point 6).

This point has not been addressed, because we think it is out of the scope of the study. We have modified the text to avoid speculation on that point. We also performed a few experiments but decided not to include them because they were not conclusive. See more details in the answer to the referee.

- Related to NHE1, is it also involved in volume regulation under osmotic stress? (Rev. 2 point 7).

This point has been addressed. The answer is that NHE1 inhibition has in general an effect on the speed of cell volume increase, so it affects the adaptation to osmotic shocks, but, as expected, it has not effect on the immediate response to the shock (the Ponder’s relation).

- Discuss on alternative Arp2/3 functions that can lead to volume regulation (Rev. 2 point 8).

This has been implemented.

Reviewer #1 (Recommendations for the authors):

1) It has to be said that the paper is quite complex to read. There is almost too much data. 3 different scenarios, spreading, osmotic shock, and confinement, are considered. While I appreciate that these phenomena are inter-related, as it is it's difficult to read. Each scenario can be explored in more detail, with models specifically addressing the scenarios. I understand this would require extensive revision, so I leave it up to author's discretion.

Consider splitting the paper and concentrate on 1 or 2 phenomena only. Personally, I think the dynamics during cell spreading and osmotic shock are significantly different. Osmotic shock involves a very large change in ion concentration, therefore might be more appropriately described by the model. During spreading and confinement, there is a small volume change, and I think the correct approach is to consider active contraction and ion dynamics together.

We apologize for the confusion introduced by the way we wrote the article. In fact, we completely agree with the referee, and we do not think that the dynamics of volume during spreading and osmotic shock are the same at all. The goal of our article is really centered on understanding volume loss during cell spreading, and our central result is that it depends on the rate of spreading, which was not reported before and is the basis for the model we propose. We added experiments on osmotic shocks and on confinement to better characterize our cells in the context of a (almost) purely osmotic response (osmotic shocks) and an (almost) purely mechanical response (fast confinement), because we believe that in the context of cell spreading, the timescale is such that both osmotic and mechanical responses overlap and are coupled, leading to our observation of a deformation rate dependent volume loss.

We have tried to clarify this point by rewriting the end of the introduction and also largely rewriting the first part of the article and moving some elements to supplementary figures or to discussion. We hope that the focus of the paper is now more clear.

2) It could be more appropriate to model spreading and confinement using a 1-solute phenomenological model, but include myosin contraction.

We respectfully disagree with the referee on that point. Because the referee revealed his name, we know that we have a disagreement on this precise point, and we have tried in our answer to his main point 2 above to explain our view. We hope that we convinced him that our view of the question is correct. In brief, on the theoretical side, we do not believe that a one solute model is an appropriate model to use in our case, because the volume changes of several % to more than 10% cannot be accounted by changes in contraction (because of the 100 mM of nonionic osmolytes present inside the cell, which cannot be ‘concentrated’ by the force produced by the acto-myosin cortex). Experimentally, we now added results showing that modulating contractility has no direct effect on cell volume, which supports our interpretation that the effect we observe of decreasing contractility during cell spreading is indirectly due to a modulation of the cell spreading speed. We agree that it could also have other indirect effects, which we now mentioned. But in any case, what we observed is that reducing contraction leads to larger volume loss during cell spreading, which goes in the opposite way of a direct effect of contraction on volume. If forces produced by contraction were large enough to have an effect on volume, myosin inhibition should lead to an increase in cell volume and not a reduction as we observed.

3) Some comments could be more rigorous. For example, the tether force is not just a function of membrane tension, but also active contraction. Comments about 1-solute model, etc.

We apologize for these approximations, and we have modified the text accordingly (e.g. top of p 8 for the tether force).

4) What is the significance of measuring tether force 4hrs vs. 30min?

This is an important point, and we apologize if it was not clearly explained in the text. These timings are the time at which the tether force is measured after cells have been plated to spread. 30 min after plating correspond to the spreading phase and the phase of volume loss, while 4 hours after platting, cells have finished spreading for several hours and have returned to steady-state growth. Accordingly, inducing faster spreading with Y27632 leads to an increase in tether force only during the spreading phase (30 min) and not at steady-state (4h). In general, for all treatments, when there is an effect leading to a higher tether force, it is observed mostly at 30 min and goes down at 4h, meaning that it is a transient increase in tether force related to the spreading phase.

We have explained that more clearly in the text when we present these results (p 8).

Reviewer #2 (Recommendations for the authors):

The authors find a correlation between dV/dt and dA/dt in spreading cells and argue that an increase in tension is associated with cell volume loss. It is however not fully clear how a fast increase in tension is build up during the initial phases of cell spreading to explain the highest change of dV/dt at short time scales.

We thank the reviewer for this point. We emphasize that one should make the distinction between the actual value of the tension and the rate at which the tension changes. In our model, the rate of volume change is directly proportional to the rate of tension change (Equation1). Although the tension needs some time to increase above its resting value (it needs the cell to spread), the rate at which the tension changes is (at short time, in the elastic regime) proportional to the rate at which the cell area increases upon spreading, which is typically fastest at the onset of spreading. This naturally leads to the conclusion that the rate of volume loss will be highest at the onset of spreading (as observed). Importantly, in our analysis, we took time zero for spreading based on the first time point at which we are able to detect spreading area with 20X magnification, and not the plating time. Once the cells start to spread, they do so at various speeds, and the spreading speed can vary also in time for each single cell, getting faster, or slower (see individual cell curves Figure S1I), but, on average, the spreading speed is faster in the first 10 to 20 minutes after the beginning of spreading, then it slows down (Figure 1G). Our analysis shows that faster spreading cells lose more volume, and that this happens mostly during the first 10 to 20 minutes after the start of spreading (Figure 1J).

Tension values were also measured via AFM at 30 min after spreading when dV/dt and dA/dt are zero. The authors should clarify how this time point is informative to argue about tension driving volume change. In this context, it remains to be tested how tension evolves at earlier time points during spreading. Given also the higher tension in cells for specific treatments at 30 min (e.g. Y-27632), why would this not lead to a longer period of cell volume decrease compared to control cells considering the authors' argument of tension homeostasis in cells? In this context, the authors could also extend their comments on the results of panel 6I in the main text and how far they match with the measured apparent membrane tension.

We agree with the reviewer that these points need to be clarified. Unfortunately, due to technical limitations, it is not possible to combine tether forces measurements and cell volume measure by FXm (one needs a cantilever, the other is performed in a closed chamber). The two measures were thus performed separately. In addition, because of the time it takes to measure tether force and the difficulty to do it on very early spreading cells, it was not possible to follow the tether force from the time zero of spreading on single cells. We thus showed two measures at two time points, one as early as possible, which corresponds to the 30-90 min time after cell plating and the second 4-5 h after plating. It is important to mention that time zero is not exactly the same for the curves showing volume measures and for the tether force. For the volume and spreading curves, time zero was taken when a contact area between the cell and the coverslip became visible by RICM, which could also vary from cell to cell, and occurred after flow equilibration in the FXm chamber. On the other hand, the time 0 for the tether force was taken as the time when cells were seeded on the dish. Also, for volume, the time 0 is set for each single cell, while for the tether force, it is a global time 0 for the dish and the measures are done across a window of time which in overall corresponds to the phase of spreading area increase, and we call it now “spreading phase” This is the best we could achieve technically. Overall, we estimate that there is about 10 minutes of delay on average between the time 0 of these two measures. We have now made it more clear in the text (p 8) and in the Methods (Supplementary file, p. 19-20).

The comment of the referee also made us realize that we do not provide a value for the initial tether force prior to spreading. Because, as explained above we cannot obtain a very early measure (when cells are not adhered at all, it is not possible to perform tether force measure), we decided to measure the tether force for adhered non-spread cells. To achieve rounded spherical cells, we used small round micropatterns of 20 µm in diameter. We found that the tether force values, for both control and Y-27632 treated cells, are the same for adhered non-spread cells on small micropatterns and for spread cells at steady-state (4h after seeding).

This control has now been added to the Figure 3E and commented in the main text (p. 8).

From the modeling point of view, our theory is based on the assumption that the membrane behaves as a viscoelastic material sensitive to the rate of change of total cell area. During the initial phase of cell spreading (timescale below the relaxation timescale \tau) the membrane is elastic, and the increase in tension is proportional to the increase in the cell area. This explains why the rate of tension increase is fastest during the initial stage of spreading, when the spreading rate is fastest. At longer timescales, the tension starts to decrease due to the viscous relaxation. Since the volume change is proportional to the negative of tension change, volume follows the same trend.

At longer timescales, the tension relaxes to its homeostatic value. At some intermediate time, which depends both on the spreading rate and the viscoelastic relaxation time, but not on the initial tension value, the tension will reach its maximum values, and the volume will be at its minimum. The time of this minimum volume varies between ten to twenty minutes for most of the experiments.

We apologize for the mislabeling of the former Figure 6I (now Figure 5I). This plot shows the theoretical estimate for the difference in tension (in the units of homeostatic tension) between the case when the cell loses its volume upon spreading (as observed in experiments) compared to the hypothetical situation when the cell does not lose volume upon spreading (α = 0). The positive value of the tension difference predicts that the cell tension would have been higher if the cell were not losing volume upon spreading, which suggests a possible physiological role for the mechano-osmotic feedback we uncovered. It also matches our experimental observations for drug treatments which reduce or abolish the volume loss during spreading and correspond to higher tether force only at short time.

The presentation of AFM measurements and data would benefit from a more detailed description of the methodology. Tether pulling experiments were indicated at 30min/4h but given the time it takes to perform a single experiment and probe multiple cells the authors should comment on the specific time duration of experiments to facilitate the interpretation of results. In addition, parameters as tether pulling positions (cell body versus protrusions) and tether pulling statistics would enable to assess the presented data more rigorously. Can the authors further comment on the highly variable contact times used for AFM tether pulling experiments, is it associated with specific experimental conditions? Were drug treatments further performed 30min prior to measurements or was the incubation performed for 4h, to clarify long-term effects of drugs and potential cell adaptations.

We apologize for this and, as explained in the point above, indeed the tether pulling experiments take some time, which is why we cannot have a better time resolution for these. We have now added a more detailed description, including the time 0 issue, as described in the answer to the point 1 of the referee (p. 8). We hope it clarifies the aspects raised by the referee. This is the updated section in Materials and methods:

“Tether pulling

For apparent membrane tension measurements, tether force was measured with single cell atomic force spectroscopy by extruding tethers from the plasma membrane on top of the nucleus of HeLa EMBL cells. Cellview glass bottom dishes (Greiner) were coated for 1 h with fibronectin (50 μg/ml; Σ-Aldrich). Cells were incubated for 30 min with in the presence of drugs or vehicle, then plated, and probed either during spreading (from 30 to 90 min after plating) or at steady state (fully spread; from 4 to 5 h after plating). To perform experiments on non-spread cells, Fibronectin-coated circles (Ø 20 um) were micropatterned onto Cellview glass bottom dishes (Greiner) using PRIMO (Alveole) following the manufacturer’s recommendations.

Tether extrusion was performed on a CellHesion 200 BioAFM (Bruker) integrated into an Eclipse Ti inverted light microscope (Nikon). OBL-10 Cantilevers (spring constant ~60 pN/nm; Bruker) were mounted on the spectrometer, calibrated using the thermal noise method (reviewed in19) and coated for 1 h at 37° C with 2.5 mg/ml Concanavalin A (Σ-Aldrich), which binds polysaccharides expressed on the surface of the cell20. Before the measurements, cantilevers were rinsed in PBS and cells were washed and probed in Dulbecco’s Modified Eagle Medium with Glutamax (DMEM/Glutamax; Gibco) supplemented with 2% FBS (Life Technologies) and 1% penicillin-streptomycin solution (Life Technologies). Measurements were run at 37° C with 5% CO2 and samples were used no longer than 1 h for data acquisition.

Tether force was measured at 0 velocity at which is linearly proportional to apparent membrane tension, assuming constant membrane bending rigidity21. In brief, approach velocity was set to 0.5 µm/s while contact force and contact time ranged between 100 to 200 pN and 100 ms to 10 s respectively. The latter two parameters were experimentally tuned before every tether pulling attempt, aiming to reach a tradeoff between the maximization of the probability to extrude single tethers, and the reduction of experimental stress on the cells. The larger the contact time and force, the higher is the probability of formation of bonds between the molecules of Concanavalin A on the surface of the cantilever and the polysaccharides on the surface of the cell. On the other side, the lower those two parameters, the lower is the stress experienced by the cell during the contact with the cantilever. As a general trend, contact force and time must be increased over the course of the experiment, owing to the depletion of Concanvalin A from the cantilever22.

To ensure tether force measurement at 0 velocity, after contacting the cell surface the cantilever was retracted for 10 µm at a velocity of 10 µm/s. The position was then kept constant for 30 s and tether force was recorded at the moment of tether breakage at a sampling rate of 2000 Hz. Each tether extrusion attempt last about 3 min and each cell was probed until 3 single tethers were successfully extruded or for a maximum of 10 min.

Force-time curves resulting from successful tether extrusions were analyzed using the JPK Data Processing Software. Tether force values from tethers extracted from the cell were then averaged, and each cell was accounted as a single data point. For every experimental condition the Shapiro–Wilk test was used to test for normality of data. Student t-test was chosen for statistical testing of normal distributed data, while Mann–Whitney U-test was performed on non-normal distributed data.”

Authors define categories of fast/intermediate and slow spreading cells. Can the authors clarify if initial cell volumes is relevant to consider and different in these categories? Does initial cell volume also change with drug treatments (e.g. Y-27632 treatment, a condition in which the cells spread faster)? The authors could further comment in the main text on potential parameters regulating spreading rate variability in cells.

We thank the referee for pointing to this initial volume parameter. In Author response image 1 we show that, for most drug treatments, as well as for the three categories of cells, the average initial volume differs by less than 10% (within the accuracy of the measure for different experiments, see Referee 1 point 1). The only exception is for treatment with NSC (the ezrin inhibitor), for which cells start with a larger volume by about 25%, and we do not have an explanation for that.

Author response image 1
Graph showing difference in % of average initial volume of 3 groups of control cells and cells treated with drugs comparing with average initial volume of control cells.

We also plotted the degree of volume loss (in %) during spreading versus the initial volume for control cells. It showed that there is no correlation between these two measures, even for differences in volume over a factor of more than 2. Which makes the differences in initial volume that we observed between conditions (Author response image 1) irrelevant for the degree of volume loss during spreading.

We have now added these controls as Figure S1K and commented them on p 5

The authors use confinement and osmotic shocks to test model assumptions and to study cell volume changes at short time scales. However, the authors so far do not consider applying the modified PLM to volume adaption in osmotic shock conditions or potential long-term volume changes in confinement. Can the mechano-sensitive PLM also capture volume adaption in osmotic shock conditions or what are the potential limitations? In addition, do cells in confinement show a volume adaptation over time that would follow similar principles as in spreading cells or osmotic shock experiments?

We thank the referee for these suggestions. We performed the osmotic shock and confinement experiments to characterize the (almost) purely osmotic and mechanical responses and we used parameters that we extracted from the osmotic shock experiments to compute membrane permeability. Our model aims at understanding, and is limited to, the timescales at which there is an overlap between these two types of responses. The timescale of volume loss due to confinement deformation is extremely fast, less than 30 ms, while the adaptation to osmotic shocks take several tens of minutes (the PLM operates in the range of minutes, which corresponds to the typical rates of ion transport). For this reason, we believe that the volume change during the confinement involves different physics than the one involved in the current model as its timescale is too fast to only involve the standard ion transport. Our model does not include fast timescales. It is a limitation of the model that is now explicitly mentioned in the text (p.14-15).

On the other hand, the referee is right that a timescale at which the model could become relevant is the adaptation/recovery from the confinement experiment. The model would suggest that, following confinement, cells could recover their volume as membrane tension decreases, in addition to steady-state growth – this would be true only if there is no loss of non-permeable osmolytes through transient pores in the membrane during confinement (e.g. amino-acids, which are a large part of these and can play a significant role in cell osmolarity).

We thus checked what was the cell increase in volume in the 20 minutes following confinement, for cells which had been confined at 20 microns (control) and at 10 microns (which induces a volume loss). These experiments are very difficult, because the height of the confinement chamber has to remain stable enough over the measurement time to allow it to serve for FXm, in the % range. Because this did not work, we used the background level of fluorescence as a proxy for the height to correct for changes in height of the confiner over the duration of the experiment. This allowed us to obtain three independent experiments for both height (in many cases the height was too unstable and experiments had to be discarded). As shown in Author response image 2, the results are consistent with the model, with more volume increase for the cells which had been confined to 10 microns. While this experiment is interesting and opens new perspectives for our model, we believe it is out of the scope of this article and will need to be repeated more times and better controlled. We thus show it just for the referee’s benefit.

Author response image 2
Graph showing the rate of volume increase after initial volume loss induced by confinement.

Each point represents the average rate of individual experiment. Average amount of cells in each experiment n=73.

Similarly, cells display a transient phase of faster growth following the spreading phase, which, according to our model, corresponds to the combination of growth and relaxation of membrane tension leading to a volume recovery phase. This can be seen when recording cells after spreading for a longer time. However, this belongs to another ongoing study in the lab, that we do not wish to include here.

From the theoretical point of view, while our model could potentially capture some aspects of the osmotic and confinement experiments (the recovery phase for example), extending it would require to introduce elements which have not been investigated in this study and is thus beyond the scope of this work. The osmotic shock experiments do however provide us with valuable information needed to support various assumptions we make in the model, which we detail below.

From the osmotic shock experiment, we find that there are two timescales in the volume dynamics. At fast timescales of seconds, there is a passive influx of water, this allows us to take the water transport at a steady state during the slow-spreading regime. At a slower timescale, there is volume recovery due to regulatory volume increase or decrease. This has been studied previously in the literature and these regulatory volume changes are associated with the regulation of different ion co-transporters like NKCC, KCl, and NHE. The model presented in this work could be easily extended to include these co-transporters, but that would involve many more parameters (see our response to Reviewer#1 point 3). We now provide a formal description of what such an extended model would look like in the Supplementary file (p 3). From our model, we can ask if a purely mechano-osmotic coupling could be sufficient to give a volume recovery upon osmotic shock. The answer to that question is that the mechano-osmotic coupling will affect the volume dynamics after the shock, but it is not sufficient to give a volume recovery at long timescales. This is because at long times the membrane tension goes back to its homeostatic value and hence the ion transport parameters are the same before the shock and a long time after the shock. However, since the osmolarity of the external medium is different even at long time the volume will not go back to the value before the shock. The mechano-osmotic coupling will indeed affect the volume recovery dynamics, but not the asymptotic value of the volume.

In summary, the volume change at different timescale is possibly regulated by different mechanisms, with some overlapping regimes, and our model focuses on the case of cell spreading.

The authors argue for a role of myosin II in counteracting volume loss in spreading cells. Would this effect also be relevant in confinement as the authors previously showed a mechanosensitive activation of myosin II in confined cells? Does volume loss occur similarly when myosin II activity is inhibited? Do the authors expect a similar role of myosin II in regulating volume change in confinement and volume adaptation in osmotic shock conditions and spreading cells?

We first want to clarify that we do not think that myosin II (or actin cortex contractility) regulates cell volume directly, as explained in the answer to referee 1 point 2. What our experiments suggest is that reducing contractility has an effect on volume during cell spreading (and not on non-spreading cells) because it accelerates the spreading speed. So myosin II is effectively counteracting volume loss in spreading cells via its effect on reducing the spreading speed. As a consequence, there is no reason for myosin II to have an effect on volume modulation under confinement, as it has no direct effect on volume (see answer to referee 1, point 2). To test it experimentally, we performed confinement with control cells and cells incubated with Y-27632 30 min prior to the experiments. We found that there is no statistically significant difference between control and Y-27632 treated cells.

Author response image 3
Graph showing average relative volume of control and Y-27 treated cells under different confinement heights.

Each point is average relative volume of N=14-30 (average amount of cells n=165) experiments and N=4-5 (average amount of cells n=166) for Y-27 treated cells. The results of t-test for each height: h=20 µm p=0.13, h=11 µm p=0.41, h=5 µm p=0.29.

It is now discussed on p 14-15.

Similarly, there is no reason (from our model perspective) for Y-27 treated cells to display a different volume adaptation following osmotic shocks. We already showed that even disrupting the actin cortex with Lat A does not significantly affect the response to osmotic shocks. Because we had already a lot of experiments to perform, we chose to not address this point experimentally (which would have required performing an entire set of osmotic shock experiments in the presence of Y-27632).

The authors explain the Gadolinium effect by referring to the fact that calcium increase might activate myosin II. Direct measurements such as MLCK treatment or different calcium interference approaches would allow to support this argument directly.

As suggested by the referee, we performed cell spreading experiments in the presence of ML-7 treatment (MLCK inhibitor). We 25 µM as was suggested in 23.We did not see any effect on average cell volume and spreading upon this treatment:

Author response image 4
Graph showing average relative volume and average spreading area of control and treated with 25 µM ML-7 cells (N=3, n=25).

Even when using 100 mM, we were not able to mimic the effect on spreading speed of Y27632 (or Gadolinium). Cell showed actually a slightly slower spreading speed, maybe due to a toxic effect of the drug at that dose. Consistently, initial volume loss was not increased compared to controls, but the cell did not recover and kept loosing volume, again maybe due to a toxic effect at this dose, since we did not observe that at lower doses.

Because these results are hard to interpret, since we have no positive control for the action of ML-7 at low doses, and high doses have little effect on spreading speed and potentially some off target effects, we cannot really interpret them and we chose to leave them out of the article.

What we wanted to test with the Gadolinium treatment was the involvement of Piezo channels, because they are classical mechanosensors. Our results show clearly that the inhibition of these channels lead to both a faster spreading speed and more volume loss, with an effect which remains on the trend defined by control cells and other treatments affecting spreading speed like Y-27632. Because we think that showing that this effect goes through the effect of Calcium on contractility (which is the most obvious from the literature), we chose to tone down our interpretation of this particular point by rephrasing our conclusion (p 11).

To complement that point, we added a novel experiment, to vary the spreading speed independently of a drug treatment. For that, we chose to add MnCl2, which is known to increase cell spreading24. The result gave a clear confirmation, with a faster spreading speed and more volume loss. This result is now added in Figure 2D and S2J and in the text (p 5).

Does NHE1 play a similar role during adaptation to osmotic stress conditions as observed in spreading cells? It would further be interesting to show whether the function of NHE1 for mechano-sensitive volume regulation is observed in the other cell types tested in this study.

We thank the reviewer for this comment. Indeed, we have not well justified our use of the NHE inhibitor EIPA. Our aim was not to directly affect the major ion pumps involved in volume regulation (which would indeed rather be the NaK exchanger), because that would likely strongly impact the initial volume of the cell and not only the volume response to spreading, making the interpretation more difficult. We based our choice on previous publication, e.g.13,showing that EIPA inhibited the main fast volume changes previously reported for cultured cells: it was shown to inhibit volume loss in spreading cells, as well as mitotic cell swelling 15,25. Using EIPA, we also found that, while the initial volume was only slightly affected, the volume loss was completely abolished even in fast spreading cells (Y-27 and EIPA combined treatment, Figure 5H). This clearly proves that the volume loss behavior can be abolished, without changing the speed of spreading, which was our main aim with this experiment.

The most direct effect of inhibiting NHE exchangers is to change the cell pH16,17, which, given the low number of H protons in the cell (negligible contribution to cells osmotic pressure), cannot affect the cell volume directly. A well-studied mechanism through which proton transport can have indirect effect on cell volume is through the effect of pH on ion transporters or due to the coupling between NHE and HCO3/Cl exchanger. The latter case is well studied in the literature18. In brief, the flux of proton out of the cell through the NHE due to Na gradient leads to an outflux of HC03 and an influx of Cl. The change in Cl concentration will have an effect on the osmolarity and cell volume.

We thus performed hyperosmotic shocks with this drug and we found that, as expected, it had no effect on the immediate volume change (the Ponder’s relation), but affected the rate of volume recovery (combined with cell growth). Overall, the cells treated with EIPA showed a faster volume increase, which is what is expected if active pumping rate is reduced. This is in contrast with the above mentioned mechanism of volume regulation which will to lead to a reduced volume recovery of EIPA treated cells. This leads us to conclude that there is potentially another effect of NHE perturbation. Changing the pH will have a large impact on the functioning of many other processes, in particular, it can have an effect on ion transport16. Overall, the cells treated with EIPA showed a faster volume increase, which is what is expected if active pumping rate is reduced.

A full mechanistic explanation of the effect of this drug is beyond the scope of this work. Because of this we are not analyzing the effect of EIPA on the model parameter α in detail. We now clarified our interpretation of these results in the main text of the article (p 11, Figure 3D and S5B).

Author response image 5
Graphs showing Ponder’s relation of control cells (the line indicated linear regression for control cells); only control cells, linear regression y=0.

67+0.33, R2=0.98 (left); only EIPA treated cells, linear regression y=0.71+0.32, R2=0.96 (right).

Author response image 6
Graph showing the effect on RVI and RVD.

We plotted average volume at 30 min after osmotic shock for control and EIPA treated cells (where P/Piso=0.61 and 0.72 are hypoosmotic shock). EIPA treatment, in a majority of cases, increased the cell volume gained during adaptation..

In conclusion, EIPA treatment tends to increase the rate of cell volume increase in all the conditions tested, which is consistent with a reduction of the pumping rate of ion pumps (Now shown in Figure S5B and p 11).

To answer the second part of the referee’s comment, we tested the role of EIPA and EIPA+Y-27 treatment on RPE-1 cells during spreading. We found that both EIPA and EIPA+Y-27 treated cells lost less volume than control, while their spreading speed was not significantly affected, which is coherent with our results for HeLa cells.

Author response image 7
Graph the relation between dV/dt and dA/dt during lamellipodia formation.

These results are now shown as Figure S5C; p 11.

Arp2/3 is identified as a molecular regulator of cell volume in spreading cells and the authors speculate on its role in regulating ion channel/pump activity. The authors could further include a discussion on alternative functions of Arp2/3 such as regulating endo/exocytosis relevant for cell volume regulation.

The discussion could further benefit from adding further arguments on the physiological relevance of cell volume regulation, considering the observed range of volume change and the focus on the specific mechanical stress conditions studied in this work.

We thank the reviewer for this suggestion. We have added a sentence on endo/exocytosis and Arp2/3 (p 13). Cell spreading was shown to be accompanied by an increase in exocytosis, providing more surface area to the plasma membrane and reducing membrane tension26. From previous literature, this is supposed to be triggered when the cell is already significantly spread, so probably later in the spreading process than the cell volume loss we are reporting here, which happens readily in the first minutes of spreading. In addition, the volumes related to endo-exocytosis events are very small compared to the volume changes considered here, so these processes are unlikely to be contributing significantly to the volume change in fast spreading cells. We also now suggest in the discussion that it might be interesting to investigate how affecting endo/exocytosis changes volume modulation during spreading, as it might contribute indirectly by modulating membrane tension (p 13).

In terms of physiological relevance, we propose a potential contribution to membrane tension homeostasis. Nevertheless, more generally, what we think our work reveals is the existence of a fundamental mechanism, which we call mechano-osmotic coupling, probably constantly at work in cells and coupling their volume change to their mechanical state. This means that it could introduce noise in the volume growth rate27 while not affecting their growth in mass (see Figure S1G). What our work suggests is thus that cell shape changes might be a source of density fluctuations. It is known that cells maintain a tight regulation of their global density, because it is a central parameter that can affect all the cell biochemistry (there is an ongoing work in our lab proving that point and showing how cells cope with that, and maintain density homeostasis despite fluctuations produced by cell shape changes). We discuss this aspect in the discussion of our article (p 14).

Statistical testing, the authors present N=1 in multiple panels of the presented data and should re-check the statistical significance of these results. Including p-values directly in the graphs would ease the reading of figure panels.

We apologize for the missing information. In most cases, experiments shown with N=1 in the main figure were representative examples, with further quantification. This was maybe not well explained and has been corrected (for most of the experiments with N=1 additional experiments were in supplementary dataset). In cases in which the number of shown experiments was really N=1, repeats are now added to the article. The changes are shown in Author response table 1. There are only few experiments left with N=1, which are more there for illustration purposes but have results not directly relevant for the core of the article. They could be removed if the referee feels that they are superfluous.

Author response table 1
Old panelNew panelStatistics improvement
1A1AIt was an error in the legend, actual N=3
1E1EThere were 2 more experiments in the attached datasets. During the revision 2 additional experiment were done. N=5
3C-D,4A-C3C-D,4A-CThe osmolarity of control medium and osmotic shock solution differ from experiment to experiment. Thus instead of averaging of experiments with similar osmolarity prior and after shock we plotted them as individual points representing the average of each individual experiment at the panels 3D, 4B and 4C. These panels also represent the total amount of osmotic shock experiment performed in this study. For the experiment with 30 s time lapse the total amount of experiments is: control N=28, Lat A N=4, EIPA N=8. For the experiment with 100 ms time lapse the total amount of experiments is: control N=11, N=5.Panels 3C and 4A are the representative examples.
5D7BThe actual confinement height differs from experiment to experiment. Thus instead of averaging of experiments with similar heights we plotted them as individual points representing the average of each individual experiment. The total amount of experiments is N=9.
5E7CNo improvements of statistics due to the limited revision time.
5F, 5I7D, 7GNo improvements of statistics due to the limited revision time.
S1AS1BThere were 2 more experiments in the attached datasets. We merged it with the experiment showed on the panel. N=3.
S1BS1AThere were 2 more experiments in the attached datasets. We merged it with the experiment showed on the panel. N=3.
S1ES1DNo improvements of statistics due to the limited revision time.
S1FS1GNo improvements of statistics due to the limited revision time.
S2AS2BDuring the revision 2 additional experiment were done. N=3
S2DThere were 2 experiments for CK-666. An additional experiment was performed. N=3.
S2ES2D-FThere were 2 more experiments for PLL-PEG in the attached datasets. We merged it with the experiment showed on the panel. N=3 for PLL-PEG.There were 2 experiments for CK-666. An additional experiment was performed. N=3.
S3BWe removed this panel.
S3CS3BNo improvements of statistics due to the limited revision time.
S3ES3DDuring the revision 2 additional experiment were done. N=3
S5JS7JThere was 1 more experiment in the attached datasets. We merged it with the experiment showed on the panel. N=2.
S6AS3GIt was an error in the legend, actual N=2
S6B-C6AIt was an error in the legend, actual N=2. But one more experiment was analyzed and added. N=3.
S6I7KNo improvements of statistics due to the limited revision time

We also included more results of statistical tests in the figure panels or legends when it was possible.

Propidium iodide staining could be used to monitor plasma membrane rupture in confinement, considering the high leakage/volume loss in confined cells. It also appears that the authors did not comment on the high increase in surface area in Figure 5H in the main text.

We thank the referee for this suggestion. We had in fact also thought of the possibility that large holes form in the plasma membrane, and we had performed experiments with Propidium iodide, but we could not observe any cell displaying a positive staining, except for a small percentage of dying cells which served as positive control that the labelling worked. It is possible that membrane ruptures (if it is indeed the mechanism of volume loss) are too transient and close fast enough to prevent a significant π staining. The confining and volume loss time in these experiments is below 100 ms, which means that the holes (if there are any) probably do not remain open longer than that.

Author response image 8
Graph showing percentage of cells with positive propidium iodide staining for different confinement heights (for each height N=2-3, average number of cells in each experiment n=196).</Author response image 8 title/legend>.

There might have been a misunderstanding on the meaning of the surface area increase shown in Figure 5H (now Figure 7F). It is written in the main text: “Below 5 μm height, the cell surface significantly increased, which also corresponded to the formation of large blebs”, but it might be confusing and we changed it into “Below 5 μm height, the cell apparent surface area significantly increased due to the formation of large blebs”. The surface area increase is mostly the surface of blebs. Blebs are forming when the plasma membrane detaches from the underlying cell cortex, this releases a lot of membrane wrinkles which constitute most of the bleb surface area, as shown by others28. It is thus not surprising that the apparent surface area increases when blebs form. This is why we did not comment it more.

LatA treatments are often used at concentrations which only partially depolymerize the cell cortex but concentrations vary and could be more clearly commented in the main text. Would higher LatA concentrations increase cell spreading in confinement and alter the cell volume change?

For cell spreading experiments, we used a low concentration (100 nM), otherwise cells would not spread at all. The data show that for this concentration, the cells spread more slowly and accordingly loose less volume. For the confinement, we used much higher doses (2 and 5 µM), which disrupted strongly the cell cortex, leading to formation of large blebs even prior to confinement. In this case the cells did not lose any volume at all even when strongly confined. Our interpretation is that the loss of cortical actin, because of the formation of blebs and the release of membrane folds, made a large amount of free membrane available and thus membrane was not tensed and the volume remained constant. We do not think that a further increase in Lat A concentration would give a different result.

Related to Figure 1G, the authors state "causing a transient density increase". Panel S1G however indicates an increased density with a stable plateau, can the authors clarify this statement?

This is because the data presented here do not run for a long enough time. The density increases because the volume drops but mass keeps increasing at a normal rate. This density increase is eventually resorbed, but over a very long time (several hours), thanks to an adaptation of the mass growth rate. This is a process that we are currently studying in details, but it is out of the scope of this study.

The authors state that "with a loss of volume reaching up to 20% for fast-spreading cells" and "actin-dependent loss of volume (up to 30%)". The max volume loss of cells is not clearly visible in the current graphs and appears exaggerated in the text. The authors could include individual data points or change to a different data illustration to visualize the range of cell volume loss better.

We apologize for this and we now show graphs which display more clearly the end of the distribution (now Figure S2A).

The authors could list the specific cell types tested in the study in various statements in the main text that appear rather generic otherwise.

We apologize if it was not clear. We went through the text and tried to fix it when it was unclear which cells had been used. The figure legends have also been checked.

How does tau_a versus tau in the table relate to tau_eff and tau in the main text? Providing a clearer introduction and overview of model parameters would support the readability.

The model parameter \tau is defined in the main text above Equation 2, it is the viscoelastic relaxation timescale of tension. The experimental parameter \tau_a is the timescale of cell spreading, which is obtained by fitting the cell speeding curve with an exponentially saturating function. The model parameter \tau_eff is defined in the main text below Equation 4. This is an effective relaxation timescale for volume that depends on \tau, volume, and contact area. The difference between \tau and \tau_eff stems from the complicated geometrical relationship between contact area and total area in a spreading cell. The two fitting parameters A0 and tau_a are obtained by fitting the contact area. And the two parameters \xi and \tau are obtained from the volume dynamics.

Comments on figure panels:

Figure 1A, cells are not well visible in this figure panel and a higher magnification or zoom in could be used.

Done.

Figure 1B, indicate channels in legend and check for consistency in figure panel.

Done.

Figure 1C, for the micropattern data indicate spreading time relevant for volume calculation.

The time for volume measurements on micropatterns was 4h after plating, the same as we used for measuring volume in cell population. It is indicated in the legend and Materials and methods now.

The authors also argue that they use micropatterns to reduce spreading area to a defined size compared to control conditions. Comparing Figure 1B and 1A, the cell areas are however similar and control cell spreading area appears to be higher in panel 1B versus 1A?

We apologize for the confusion. We now say more clearly in the text that these are different cell lines: HeLa EMBL control cells (also called HeLa Kyoto) for patterns and most of the experiments in the paper and a sub clone of HeLa cells expressing hgem-mcherry which turn out to have slightly different characteristics although they also are HeLa cells.

Figure 1I, J, indicate in the legend at which time point dV/dt was measured.

Done.

How was the scale bar in panel 1I calculated?

Color bar represents Kernel density. Kernel density estimation (KDE) is a non-parametric way to estimate the probability density function of a random variable. It was calculated with Python scipy.stats.gaussian_kde(). https://github.com/scipy/scipy/blob/v1.7.1/scipy/stats/kde.py#L41-L637

We replaced this color code now by 3 groups of control cells as was suggested by the Reviewer 3.

Figure 2C, there appears to be a gap in the data of Y-27632?

The referee is right. This is due to the averaging of several experiments with an outlier which pulled the average down and the corresponding data point was removed. The full graph is shown in Author response image 9. We also described the outlier removal procedure in the Material and Method and the Transparent report.

Author response image 9

Figure 3A, the explanation of the schematic of the PLM could be explained in more detail to ease the reading.

We tried to improve it.

Figure 3D vs 4B, why is the Ponders equation in panel 3D independent of actin but not in panel 4B? LatA data in panel 4B appear to be not commented in main text.

We are not sure we understand this comment of the referee. Figure 4B does not correspond to the Ponder’s relation, but to the speed of volume change during the initial passive response to the osmotic shock. The referee probably thinks that there is an effect of Lat A because the points for Lat A do not cover all the P/Piso, but only a small range around 1. On the contrary, they show that there is no effect. To make it clearer, Lat A points have been removed from the panel 4B.

Figure 4B, check sign of dV/dt compared to other figure panels as it appears neg/pos for hypo/hyper conditions when cells swell/shrink.

Done. We checked but we have not detected an issue. Note that the x axis in Figure 4B is P/Piso meaning above 1 when there is an hyper osmotic shock (thus volume decreases and dv/dt is negative, as can be seen on the graph).

Figure 4C, over which duration was adaptation measured?

We apologize for the unclear sentence in Materials and methods: “Adaptation speed for osmotic shock recovery was calculated as a linear slope at 5 min intervals.” We changed it to make it clearer: “Adaptation speed for osmotic shock recovery was calculated as a linear slope starting from the minimum or maximum volume value achieved during passive response (for hyper or hypoosmotic shock) on 5 min interval”

Figure 5E, the legend indicates that a double expression line of myosin II/actin was used but only one channel is shown (membrane?). The visualization of myosin II and actin would be informative to include.

We used this cell line, but our aim was mostly to visualize the membrane in this particular experiment, which was labelled with CellMask Far Red (this is what was imaged and what is shown). We have removed the mention that the cell actually expresses the other markers, since we do not use them.

Figure 6K, indicating measured tension values in Figure 6K or matching the color code between panel 6J and 6K would ease the reading of figures.

Done – note that when data were shown before in another graph, we made them appear in grey, to make sure it is clear that they are the same dataset and not a new one for the particular panel

Figure S2G Why is the correlation before and after lamellopodia changing slope, appears to be not commented and the figure is not cited in the main text.

We observed that RPE1 cells display an initial short phase during which volume increases, which is dependent on Myosin II contractility (it disappears in Y-27632 treated cells). It seems that the extent of this initial small increase in volume also depends on the speed of the very initial spreading (in the first ~4 minutes). This is what is shown on the graph. Once the cell starts to spread and extend a lamellipodia, the volume decreases like in HeLa cells. We have no explanation for this initial bump in volume in RPE1 cells. We now say that more explicitly in the text (p 5).

Figure S2F, what do the lines indicate?

These are linear regressions. It is now mentioned in the legend

Reviewer #3 (Recommendations for the authors):

It would be good to insert figure S4H in main figure, as it clearly how speed of deformation affects the volume change.

There is no Figure S4H, the referee is probably mentioning another panel but we could not identify which one it is. In general, we have tried to rethink which panels should be shown in the main figures and which should be placed in supplementary figures, to make it more clear that the focus of the article is on cell spreading. If the referee mentioned the panel S6I we now replaced it with the Model Figure 2 and it is the Figure 7K now.

The absence of correlation between volume and spreading area at a given spreading rate is interesting and new. It is probably worth discuss the authors' results in light of Ref 11, that conclude on a weak correlation between spread area and volume, and ref 9 that observe a strong correlation between area/volume.

We thank the reviewer for this comment, and we agree that this is actually the central result and message of our article. We apologize if it was not clear enough; we have now better highlighted this result by rewriting extensively some parts of the article.

Fig1F

Could the authors consider including the spreading area change in time on PLL-PEG?

The spreading area that we measure in the article is more precisely the contact area measured by RICM. In the case of PLL-PEG there is no detectable contact area in RICM. The cells ‘float’ over the substrate, they cannot attach on it, at least in the duration of the experiments.

Fig1D-H

Color code is not consistent and make the figure harder to read.

We tried our best to improve that point.

Page 5, first paragraph

Please explain which criteria you used to pool the cells according to their spreading speed: do you have 3 clear populations, well separated by spreading speed, or you had to apply cutoffs sort cells in different populations?

The cells were simply separated in three groups of equal sizes. This is now better illustrated showing the groups on the speed distribution. With this grouping, we still had enough cells in each group (~ 40 cells) and enough groups to argue some correlation between spreading speed and volume.

Could the authors indicate the three corresponding populations ('slow spreading', 'moderate spreading', 'fast spreading') on Figure 1I?

We thank the referee for the suggestion, it is actually more useful than the color code we had used originally, which was not very informative.

Consider clarifying the color code.

We have replaced the color code (which was giving the local density of points) and was not very useful, by the three groups of cells, as suggested by the referee.

Can the authors plot Figure 1J for each population ('slow spreading', 'moderate spreading', 'fast spreading')? Is the same tendency still observed?

The tendency is observed within each group, as shown on the graphs in Author response image 10, with of course less range of spreading speed for slower spreading cells. In all cases cells tend to spread faster at the beginning and loose more volume, then slow down and stop losing volume. We do not feel the graphs are necessary to add to the article, but they can be added if the referee thinks they are useful.

Author response image 10
Graphs showing relation between average dV/dt and dA/dt for different time intervals of spreading for 3 groups of cells: fast, moderate and slow spreading.

Page 5, second paragraph

In RPE-1, dynamics of spreading area (Figure S2D) and volume (Figure S2E right) for control cell and under Y-27 are very similar, while it is not the case in HeLa cells. Can the authors comment on it?

We do not have a good explanation for that. Our simplest explanation is that RPE1 cells already have a very low basal contractility, as we reported before many times. So inhibiting myosin II activity has a mild effect on spreading speed. It seems to mostly affect the very early spreading dynamics – our interpretation of that is that the initial contact of the RPE1 cells induces the classical cascade directed by integrin binding and leading to a transient initial contraction followed by actual spreading and lamellipodia extension. On the other hand, HeLa cells have a high basal contractility level and do not display this initial contraction. They immediately spread, but with a high basal contractility which limits their spreading speed. To test these hypothesis, we would need to study in details the process of spreading in these two cells lines and the state of the acto-myosin cytoskeleton. It would be very interesting but is out of the scope of this study, which focuses on the coupling between spreading and volume changes.

Figure 2B-C

Axes do not have the same minima/maxima so complicated to read and to compare with Figure 1H.

Note that in 1H, the goal is the comparison, in control cells, between different spreading speeds, while 2B shows the effect of Y-27 which produces faster spreading than the fastest spreading cells in the controls. All these are then placed together on a single graph for comparison in Figure 5K.

Figure 2D

Axis do not have the same minima/maxima. It would be easier to read and compare with Figure 1I if they had the same range.

This is the same point as above. If we use the same range, all the points in figure 1J will be pushed together, it will not help the reading.

Are the linear fits in fig2D and fig1I giving the same results? They seemed to be different.

This is a good point from the referee. They were not the same. We have now changed that in all figures reporting this dV/dt versus dA/dt to always use the fit taken from the control cells. It shows that the Y27632 treated cells loose a little bit more volume than what control cells would at that speed of spreading if the relation was purely linear on the whole range of speeds. This could be due to an additional effect of Y27/contractility, but it is a small one and the point falls quite close to the control trend.

Page 5, third paragraph

The results of this part are very similar to the paper cited above (https://doi.org/10.1101/2021.01.22.427801), probably worth mentioning.

This is indeed true, since these are almost the same experiments (volume is just measured in a different way). We have now added the citation and mentioned the similarity (p6).

In the main text, authors should specify that you used PEG400 to apply hyperosmotic shocks.

Done (p 6)

Do the authors find differences if they apply osmotic shocks with salts or adding non-metabolized sugars such as sucrose or sorbitol?

We chose not to do this set of experiments, for the sake of time, because the central message of the article is not on the osmotic response. These experiments were done to characterize the response of the cells, to support some hypothesis of the model. However, it is well accepted in the literature, that Ponder’s relation is not affected by the choice of osmotic agent29.

Page 6, second paragraph

Under Y27 treatment, cell spread faster and their membrane tension measured using AFM is 30% higher than control 30 min after seeding. As LatA treated cells spread slower, wouldn't the authors expect a smaller membrane tension? Dai and Sheetz-1999 looked at the effect of cytochalasin D and indeed observed a reduced membrane tension.

The referee is right. However, we have not included that condition. As the referee states, this is an expected result, and given the number of conditions already included, we thought that it would not add much.

Page 8

What are the differences, beside the timescale, between Figure 3C and Figure 4A? It would be worth mentioning the way these different were acquired in the figure legend.

The only difference is the time-lapse of the acquisition. In 3C, images were acquired every 30 second over 60 minutes, while in 4A they were acquired every 100 ms over ~7 seconds. This is now mentioned in the figure legends.

Page 9 last paragraph + page 10 first paragraph

Is intensity of CellMask a proper tool to extrapolate on membrane folds and reservoirs (Figure 5E)? The spread around the contour line seems larger for 7.5um than for no confinement for control cells, and the reverse for Lat A treated cells. I am not sure how these intensity plot can be quantitatively linked to the size of the membrane reservoir.

We believe that CellMask dyes are a good proxy to quantify the amount of membrane. The spread around the contour line is actually shown quantitatively on the right graph next to the images. And it is very clear from the graph that the width at mid-height is not more for 7.5 than controls – it is mostly given by the optical resolution. We now give the non-confined. What matters is the integrated intensity. These dies contain a lipophilic moiety for membrane loading and a negatively charged hydrophilic dye for “anchoring” of the probe in the plasma membrane. These amphipathic dyes provide a uniform staining of the plasma membrane. While they cannot be used to assess the presence of nanofolds or domains, they can provide bulk information on the surface area of the plasma membrane in contact with the extracellular medium, and they stain large folds and membrane ruffles. Depletion of membrane ruffles and folds would result in less CellMask 30–32. In this article the drop of CellMask intensity at the plasma membrane serves an an indication of depletion of fold, ruffles and small filopodia. In agreement with this, in cells stained with CellMask, blebs display a thinner and dimmer CellMask signal compared to the cortex-bound membrane – this is because blebs have a bare flat membrane and thus contain less plasma membrane reservoirs than the cortex-bound28,33. We thus think that our measure, although maybe not perfectly quantitative for the amount of membrane, can be used to assess the difference in the amount of membrane folds between conditions.

Also, in the before to last paragraph in page 9, I kindly disagree with the statement that "membrane tension could not arise from a limitation in membrane availability". Isn't it that unfolding the membrane partially, not reaching the full size of the membrane reservoir, would still change membrane tension, because of opposing cytoskeletal and adhesive forces trying to fold the membrane?

We agree with the referee and we apologize for the sentence which was not clear enough. What we mean is that the increase in membrane tension could not arise from a limitation in the total amount of membrane, since there is a very large reservoir. And indeed, it increases due to the phenomenon mentioned by the referee. We modified the sentence accordingly (p 14)

Figure 5F

LatA treated cell keep a constant volume after confinement while control cell decrease volume. Since you have a condition (Y27 treated cell) that have opposing effect of LatA on spreading area and volume, I wonder what would be the effect of confinement on Y27 treated cells. Would cells 'overreact?'

The condition of Y27 treatment does not have a direct effect on volume in the cell spreading experiments, but it has an effect only via the increase in spreading speed. We have now added experiments to prove that more directly (see Referee 1 point 2). As a consequence, it should not have any effect in the confinement experiments, and this is what we find experimentally (see Referee 1 point 2)

Figure 6 A-H

Y-axis scale is different for each graph, would be better to fit one range and keep it so that the reader can compare from one to another.

Because the effects of the various drugs have different amplitude, it is not really possible to use the same scales for all these graphs. The comparison can be made on the panel K, where all the values are placed together for that purpose.

Figure 6J

Do you affect absolute cell volume when adding those drugs?

Because the effects of the various drugs have different amplitude, it is not really possible to use the same scales for all these graphs. The comparison can be made on the panel K, where all the values are placed together for that purpose.

[Editors' note: further revisions were suggested prior to acceptance, as described below.]

Essential revisions:

We all agreed that the revised version of the manuscript is much improved and clearly presents a very interesting and compelling story, well done!

Before the formal acceptance of the manuscript, one of the reviewers (rev. #1) made two final comments regarding (1) the contractile stress's role in volume control; and (2) the relative importance of membrane tension and voltage for the effects presented here. I think it might be a good idea to consider these comments and decide whether you think you can adapt the discussion part of the text so your points are clearly laid down.

Reviewer #1 (Recommendations for the authors):

The authors have addressed my comments and the manuscript is much improved. I can recommend publication. I think there is generally no disagreement in our thinking of contractility and cell volume regulation. It is true that contractile force, or pressure, alone cannot really change the volume very much. This was already shown in Yellin et al., Biophys. J. 2018, and Li and Sun, Frontiers Cell Dev. Bio., 2022, and I agree with the authors' general thinking. However, the discussion in Pg. 7, I think, can be improved. Right now, it sounds like contractile stress is being compared with the total osmotic pressure of impermeable solutes inside the cell. But in reality, the contractile stress should be compared with the total osmotic pressure difference between inside and outside of the cell, which is much smaller than 100mM. During 10% volume change during cell spreading, I think since it's happening slowly (order of minutes), the ion flux and contractile stress should change simultaneously, possibly with contractile stress changing slowly, while ion flux happening very fast. But, the contractile stress must change, since at minimum, the cell is changing shape, and therefore cortical stress must change to maintain relative force balance. I'm still not completely sure whether is ok to neglect it, since the steady-state result of the proposed model is zero osmotic pressure difference between inside and outside, but in reality, there is about 1kPa of hydraulic pressure difference and 1kPa osmotic pressure (~0.5mM) difference. This difference may be small enough to be negligible in the current context. But it may not be small in other contexts, such as during cell motility.

We have added the references that the referee mentioned, which support our argument that hydrostatic pressure has a negligible role in volume changes that are more than 1%. We have also emphasised the importance of membrane and cortical tension in other cellular phenomena such as shape change, which are not considered here since the spreading rate is a measured input of the model (p.7 2nd paragraph). We are not sure to understand the referee’s point about the significance of a 1kPa osmotic pressure difference in the context of cell motility. The motility data reported in our paper (for a particular cell type – dendritic cells – migrating in collagen) shows volume changes of order 10-20% (Figure 7), which cannot be explained by a 1% modulation of the osmotic pressure difference.

The other issue I hope the authors can clarify in the discussion is the equal importance of membrane voltage and membrane tension in this story. Ion exchangers such as NHE and NaK are affected by both. The current dependence of exchanger flux with respect to voltage is based on just passing an ion through a constant electric field. This is generally not true for ion exchangers. I think most studies of ion exchangers show a voltage gating effect. This, after all, is the basis of action potentials. Therefore, these important facts should be clearly stated.

We certainly agree with the referee that voltage-gating is very important for many type of channels and exchanger. We now more clearly state the assumptions behind our estimate of the coupling parameter \α, and their limitations (p.9 below Equation1). However, we are unsure about what the referee means when he mentions the "the equal importance of membrane voltage and membrane tension in this story”. In our mode, the ion concentrations and the potential difference are dynamical variables of the PLM. The only coupling between the PLM and cell spreading is through membrane tension. Changes of tension during spreading will certainly affect the transmembrane potential, which in turn affects the rates of ion transport, but this is capture by Equation1 at the linear level. we have now made this more explicit in the text (p.9 below Equation1).

https://doi.org/10.7554/eLife.72381.sa2

Article and author information

Author details

  1. Larisa Venkova

    1. Institut Curie, PSL Research University, CNRS, UMR 144, Paris, France
    2. Institut Pierre Gilles de Gennes, PSL Research University, Paris, France
    Present address
    Institute of Biochemistry and Cellular Genetics, University of Bordeaux, Paris, France
    Contribution
    Conceptualization, Formal analysis, Investigation, LV performed the imaging, volume measurements, analysis for cell spreading experiments with help of DC, BD, NS and AD. LV performed the imaging, volume measurements, analysis for cell confinement and osmotic shock experiments. LV. and MP. conceived the study and designed experiments. LV and MP wrote the manuscript., Methodology, Software, Visualization, Writing – original draft, Writing – review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-5721-7962
  2. Amit Singh Vishen

    Institut Curie, PSL Research University, CNRS, UMR 168, Paris, France
    Contribution
    ASV developed the mechano-sensitive PLM model under the supervision of J-FJ and PS and fitted the model predictions to the experimental data, Conceptualization, Formal analysis, Investigation, Methodology, Software, Supervision, Visualization, Writing – original draft, Writing – review and editing, Lead author for the physical model
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-8938-0852
  3. Sergio Lembo

    Cell Biology and Biophysics Unit, European Molecular Biology Laboratory, Heidelberg, Germany
    Contribution
    Formal analysis, Investigation, Methodology, SL performed the tether force measurements under the supervision of AD-M, Supervision, Visualization, Writing – review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-2253-8771
  4. Nishit Srivastava

    1. Institut Curie, PSL Research University, CNRS, UMR 144, Paris, France
    2. Institut Pierre Gilles de Gennes, PSL Research University, Paris, France
    Contribution
    Formal analysis, Investigation, Methodology, NS performed mass measurements and cell cycle experiments
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-4177-6123
  5. Baptiste Duchamp

    1. Institut Curie, PSL Research University, CNRS, UMR 144, Paris, France
    2. Institut Pierre Gilles de Gennes, PSL Research University, Paris, France
    Contribution
    BD performed imaging, volume measurements and analysis of cell plated on micropatterns under the supervision of DC, Formal analysis, Investigation, Methodology, Supervision, Writing – review and editing
    Competing interests
    No competing interests declared
  6. Artur Ruppel

    Laboratoire Interdisciplinaire de Physique, Grenoble, France
    Contribution
    AR performed the traction force experiments under the supervision of MB, Conceptualization, Formal analysis, Investigation, Methodology, Supervision, Visualization, Writing – review and editing
    Competing interests
    No competing interests declared
  7. Alice Williart

    1. Institut Curie, PSL Research University, CNRS, UMR 144, Paris, France
    2. Institut Pierre Gilles de Gennes, PSL Research University, Paris, France
    Contribution
    AW and LV performed the experiments with the FLIPPER probe, and AW performed the image analysis, Formal analysis, Investigation
    Competing interests
    No competing interests declared
  8. Stéphane Vassilopoulos

    Sorbonne Université, Institut National de la Santé et de la Recherche Médicale, Institut de Myologie, Centre de Recherche en Myologie, Paris, France
    Contribution
    Conceptualization, Formal analysis, Investigation, Methodology, SV performed electron microscopy experiments, Visualization
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-0172-330X
  9. Alexandre Deslys

    1. Institut Curie, PSL Research University, CNRS, UMR 144, Paris, France
    2. Institut Pierre Gilles de Gennes, PSL Research University, Paris, France
    Contribution
    Formal analysis, Investigation, Methodology
    Competing interests
    No competing interests declared
  10. Juan Manuel Garcia Arcos

    1. Institut Curie, PSL Research University, CNRS, UMR 144, Paris, France
    2. Institut Pierre Gilles de Gennes, PSL Research University, Paris, France
    Contribution
    Conceptualization, Formal analysis, Investigation, JMGA performed the analysis of cell membranes under confinement, Methodology, Supervision, Validation, Writing – review and editing
    Competing interests
    No competing interests declared
  11. Alba Diz-Muñoz

    Cell Biology and Biophysics Unit, European Molecular Biology Laboratory, Heidelberg, Germany
    Contribution
    ADM supervised the tether pulling experiments, Conceptualization, Supervision, Validation, Writing – review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-6864-8901
  12. Martial Balland

    Laboratoire Interdisciplinaire de Physique, Grenoble, France
    Contribution
    Conceptualization, MB supervised the traction force microscopy experiments, Methodology, Supervision, Validation, Writing – review and editing
    Competing interests
    No competing interests declared
  13. Jean-François Joanny

    Institut Curie, PSL Research University, CNRS, UMR 168, Paris, France
    Contribution
    Conceptualization, Formal analysis, Investigation, JFJ participated in establishing and supervising the physical modeling, Methodology, Supervision, Validation, Writing – review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-6966-3222
  14. Damien Cuvelier

    1. Institut Curie, PSL Research University, CNRS, UMR 144, Paris, France
    2. Institut Pierre Gilles de Gennes, PSL Research University, Paris, France
    3. Sorbonne Université, INSERM, Paris, France
    Contribution
    Conceptualization, DC supervised the micropatterning experiments, Formal analysis, Investigation, Methodology, Validation, Writing – original draft, Writing – review and editing
    Competing interests
    No competing interests declared
  15. Pierre Sens

    Institut Curie, PSL Research University, CNRS, UMR 168, Paris, France
    Contribution
    Conceptualization, Investigation, Methodology, PS supervised and coordinated the physical modeling, Supervision, Validation, Writing – original draft, Writing – review and editing
    For correspondence
    pierre.sens@curie.fr
    Competing interests
    Reviewing editor, eLife
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-4523-3791
  16. Matthieu Piel

    1. Institut Curie, PSL Research University, CNRS, UMR 144, Paris, France
    2. Institut Pierre Gilles de Gennes, PSL Research University, Paris, France
    Contribution
    Conceptualization, MP supervised and coordinated the study, Methodology, Project administration, Supervision, Validation, Writing – original draft, Writing – review and editing
    For correspondence
    matthieu.piel@curie.fr
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-2848-177X

Funding

Agence Nationale de la Recherche (ANR-19-CE13-0030)

  • Matthieu Piel

Agence Nationale de la Recherche (ANR-10-EQPX-34)

  • Matthieu Piel

Agence Nationale de la Recherche (ANR-10-IDEX-0001-02 PSL)

  • Matthieu Piel

Agence Nationale de la Recherche (ANR-10-LABX-31)

  • Matthieu Piel

Fondation pour la Recherche Médicale (FDT201805005592)

  • Larisa Venkova

Human Frontier Science Program (LT000305/2018-L)

  • Nishit Srivastava

Agence Nationale de la Recherche (ANR‐17‐CE13‐0020‐02)

  • Amit Singh Vishen

European Union's Horizon 2020 research and innovation programme (Marie Sklodowska-Curie grant agreement no. 641639)

  • Larisa Venkova

Human Frontier Science Program

  • Nishit Srivastava

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

Acknowledgements

The authors thank Jian Shi, Rafaële Attia, and Li Wang for help with photolithography; Pierre Recho and Romain Rollin for fruitful discussions; Pablo Vargas for help with DCs experiment; Olivier Thouvenin and Matthieu Maurin for assistance with optical details; Vincent Frasier for help with 3D-reconstruction of confined cells. This study was supported by a French Agence Nationale de la Recherche (ANR) grant to MP (ANR-19-CE13-0030). This study has also received the support of Institut Pierre-Gilles de Gennes-IPGG (Equipement d’Excellence, 'Investissements d’avenir,' program ANR-10-EQPX-34) and laboratoire d’excellence, 'Investissements d’avenir' program ANR-10-IDEX-0001-02 PSL and ANR-10-LABX-31. The authors also thank IPGG technical platform for providing equipment and technical assistance. LV has received funding from the European Union’s Horizon 2020 Research and Innovation Program under the Marie Sklodowska-Curie Grant agreement no. 641639, and Fondation pour la Recherche Médicale (FDT201805005592). NS acknowledges support from the Human Frontier Science Program (LT000305/2018 L). ASV acknowledges the support of ANR under reference ANR‐17‐CE13‐0020‐02.

Senior Editor

  1. Anna Akhmanova, Utrecht University, Netherlands

Reviewing Editor

  1. Felix Campelo, Institute of Photonic Sciences, Spain

Reviewer

  1. Sean X Sun, Johns Hopkins University, United States

Version history

  1. Preprint posted: June 8, 2021 (view preprint)
  2. Received: July 21, 2021
  3. Accepted: April 8, 2022
  4. Accepted Manuscript published: April 13, 2022 (version 1)
  5. Version of Record published: May 10, 2022 (version 2)

Copyright

© 2022, Venkova et al.

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.

Metrics

  • 5,585
    Page views
  • 966
    Downloads
  • 14
    Citations

Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Larisa Venkova
  2. Amit Singh Vishen
  3. Sergio Lembo
  4. Nishit Srivastava
  5. Baptiste Duchamp
  6. Artur Ruppel
  7. Alice Williart
  8. Stéphane Vassilopoulos
  9. Alexandre Deslys
  10. Juan Manuel Garcia Arcos
  11. Alba Diz-Muñoz
  12. Martial Balland
  13. Jean-François Joanny
  14. Damien Cuvelier
  15. Pierre Sens
  16. Matthieu Piel
(2022)
A mechano-osmotic feedback couples cell volume to the rate of cell deformation
eLife 11:e72381.
https://doi.org/10.7554/eLife.72381

Further reading

    1. Cell Biology
    Herschel S Dhekne, Francesca Tonelli ... Suzanne R Pfeffer
    Research Advance Updated

    Activating mutations in the leucine-rich repeat kinase 2 (LRRK2) cause Parkinson’s disease. LRRK2 phosphorylates a subset of Rab GTPases, particularly Rab10 and Rab8A, and we showed previously that these phosphoRabs play an important role in LRRK2 membrane recruitment and activation (Vides et al., 2022). To learn more about LRRK2 pathway regulation, we carried out an unbiased, CRISPR-based genome-wide screen to identify modifiers of cellular phosphoRab10 levels. A flow cytometry assay was developed to detect changes in phosphoRab10 levels in pools of mouse NIH-3T3 cells harboring unique CRISPR guide sequences. Multiple negative and positive regulators were identified; surprisingly, knockout of the Rab12 gene was especially effective in decreasing phosphoRab10 levels in multiple cell types and knockout mouse tissues. Rab-driven increases in phosphoRab10 were specific for Rab12, LRRK2-dependent and PPM1H phosphatase-reversible, and did not require Rab12 phosphorylation; they were seen with wild type and pathogenic G2019S and R1441C LRRK2. As expected for a protein that regulates LRRK2 activity, Rab12 also influenced primary cilia formation. AlphaFold modeling revealed a novel Rab12 binding site in the LRRK2 Armadillo domain, and we show that residues predicted to be essential for Rab12 interaction at this site influence phosphoRab10 and phosphoRab12 levels in a manner distinct from Rab29 activation of LRRK2. Our data show that Rab12 binding to a new site in the LRRK2 Armadillo domain activates LRRK2 kinase for Rab phosphorylation and could serve as a new therapeutic target for a novel class of LRRK2 inhibitors that do not target the kinase domain.

    1. Cell Biology
    Ling-Yun Zhou, Chen-Xi Jin ... Hao Wu
    Research Article Updated

    The MRTF–SRF pathway has been extensively studied for its crucial role in driving the expression of a large number of genes involved in actin cytoskeleton of various cell types. However, the specific contribution of MRTF–SRF in hair cells remains unknown. In this study, we showed that hair cell-specific deletion of Srf or Mrtfb, but not Mrtfa, leads to similar defects in the development of stereocilia dimensions and the maintenance of cuticular plate integrity. We used fluorescence-activated cell sorting-based hair cell RNA-Seq analysis to investigate the mechanistic underpinnings of the changes observed in Srf and Mrtfb mutants, respectively. Interestingly, the transcriptome analysis revealed distinct profiles of genes regulated by Srf and Mrtfb, suggesting different transcriptional regulation mechanisms of actin cytoskeleton activities mediated by Srf and Mrtfb. Exogenous delivery of calponin 2 using Adeno-associated virus transduction in Srf mutants partially rescued the impairments of stereocilia dimensions and the F-actin intensity of cuticular plate, suggesting the involvement of Cnn2, as an Srf downstream target, in regulating the hair bundle morphology and cuticular plate actin cytoskeleton organization. Our study uncovers, for the first time, the unexpected differential transcriptional regulation of actin cytoskeleton mediated by Srf and Mrtfb in hair cells, and also demonstrates the critical role of SRF–CNN2 in modulating actin dynamics of the stereocilia and cuticular plate, providing new insights into the molecular mechanism underlying hair cell development and maintenance.