Kalpana Mandal

Atef Asnacios

Bruno Goud

Jean-Baptiste Manneville
email: baptiste.manneville@curie.fr

Mapping intracellular mechanics on micropatterned substrates

Keywords: microrheology, optical tweezers, viscoelasticity, cytoskeleton, active forces, membrane, cancer, micropatterning, spatial analysis, power law, standard linear liquid

Significance statement

Cell mechanics is crucial for many, if not all, cell functions. Although mechanics at the scale of the whole cell is extensively documented, there are only few methods to measure intracellular mechanics at the local scale. Here we develop a technique that allows us 1) to map the spatial variations of intracellular mechanical parameters; 2) to study how the actin and microtubule cytoskeleton, intracellular membranes and ATP-dependent active forces contribute to intracellular mechanics; and 3) to differentiate normal and cancer cells. Because intracellular mechanical maps can detect subtle differences in the spatial distribution of mechanical parameters even in the absence of any change in their average values, our approach could provide a new diagnostic and prognostic tool for cancers.

Introduction

Cell mechanics play a crucial role in many cellular functions that are critical during development or are altered during pathologies such as cancers. For instance, cell migration, cell adhesion and cell division have all been shown to depend on cell mechanics [START_REF] Guillot | Mechanics of epithelial tissue homeostasis and morphogenesis[END_REF][START_REF] Lange | Cell and tissue mechanics in cell migration[END_REF][START_REF] Louveaux | The mechanics behind cell division[END_REF][START_REF] Vedula | Collective cell migration: a mechanistic perspective[END_REF][START_REF] Asnacios | The mechanics behind cell polarity[END_REF]. Most studies have focused on the mechanical crosstalk between the cell and its microenvironment at the whole cell or the tissue level. It has been shown recently that intracellular mechanics could also be used to differentiate cancer cells from normal cells [START_REF] Li | Intracellular particle tracking as a tool for tumor cell characterization[END_REF][START_REF] Baker | Cancer cell stiffness: integrated roles of threedimensional matrix stiffness and transforming potential[END_REF][START_REF] Guo | Probing the Stochastic, Motor-Driven Properties of the Cytoplasm Using Force Spectrum Microscopy[END_REF][START_REF] Agus | A physical sciences network characterization of non-tumorigenic and metastatic cells[END_REF]. Several intracellular elements participate in the mechanical response of the cytoplasm and most of them may be modified during cancer progression. The three types of fibers constituting the cytoskeleton, actin, microtubules and intermediate filaments, clearly contribute to cell mechanics [START_REF] Huber | Cytoskeletal crosstalk: when three different personalities team up[END_REF][START_REF] Wen | Polymer physics of the cytoskeleton[END_REF][START_REF] Coughlin | Changes in cytoskeletal dynamics and nonlinear rheology with metastatic ability in cancer cell lines[END_REF][START_REF] Laudadio | Rat airway smooth muscle cell during actin modulation: rheology and glassy dynamics[END_REF].

The role of internal membranes, molecular crowding or active force generators in the cytoplasm is less documented but could also be significant. During cancer development, signaling associated to the cytoskeleton as well as membrane trafficking, cell polarity and intracellular organization are deregulated [START_REF] Halaoui | Rewiring cell polarity signaling in cancer[END_REF][START_REF] Goldenring | A central role for vesicle trafficking in epithelial neoplasia: intracellular highways to carcinogenesis[END_REF][START_REF] Zegers | Rho GTPases in collective cell migration[END_REF], supporting the idea that the mechanics of the interior of cancer cells could differ from that of normal cells.

Passive or active microrheology techniques have been developed to measure intracellular mechanical properties, mostly based on the tracking of endogenous granules or internalized particles [START_REF] Crocker | Multiple-particle tracking and two-point microrheology in cells[END_REF][START_REF] Weihs | Bio-microrheology: a frontier in microrheology[END_REF][START_REF] Bausch | Measurement of local viscoelasticity and forces in living cells by magnetic tweezers[END_REF][START_REF] Broedersz | Measurement of nonlinear rheology of cross-linked biopolymer gels[END_REF][START_REF] Almonacid | Active diffusion positions the nucleus in mouse oocytes[END_REF][START_REF] Berret | Local viscoelasticity of living cells measured by rotational magnetic spectroscopy[END_REF]. The experimental results are generally analyzed using two main classes of models, viscoelastic models with a finite number of viscous (dashpots) and elastic (springs) elements and power-law models (see [START_REF] Kasza | The cell as a material[END_REF][START_REF] Stamenović | Rheological behavior of mammalian cells[END_REF][START_REF] Lim | Mechanical models for living cells--a review[END_REF][START_REF] Pritchard | Mechanics of biological networks: from the cell cytoskeleton to connective tissue[END_REF] for reviews). An increasing amount of experimental evidence point to weak power-law rheology as a common feature of cell mechanical responses [START_REF] Guo | Probing the Stochastic, Motor-Driven Properties of the Cytoplasm Using Force Spectrum Microscopy[END_REF][START_REF] Fabry | Scaling the microrheology of living cells[END_REF][START_REF] Cai | Quantifying cell-to-cell variation in power-law rheology[END_REF][START_REF] Balland | Power laws in microrheology experiments on living cells: Comparative analysis and modeling[END_REF][START_REF]A comparative mechanical analysis of plant and animal cells reveals convergence across kingdoms[END_REF][START_REF] Hoffman | The consensus mechanics of cultured mammalian cells[END_REF][START_REF] Guo | The role of vimentin intermediate filaments in cortical and cytoplasmic mechanics[END_REF]. Previous experiments and associated modelling have aimed at measuring intracellular mechanical parameters averaged over the whole cell interior. However, the cell interior is an extremely heterogeneous medium characterized by length scales varying on several orders of magnitude from the nanometer scale to the micrometer scale.

As a consequence, the intracellular viscoelasticity strongly depends on the position probed within the cell cytoplasm [START_REF] Tseng | Micromechanical mapping of live cells by multiple-particle-tracking microrheology[END_REF][START_REF] Hecht | Imaging viscoelastic properties of live cells by AFM: power-law rheology on the nanoscale[END_REF][START_REF] Daniels | Differences in the microrheology of human embryonic stem cells and human induced pluripotent stem cells[END_REF]. The spatial variations of intracellular mechanics are hidden when only cell-averaged parameters are measured while they clearly contain relevant information regarding the internal organization of the cell. In the context of cancer progression, it could thus be crucial to measure not only the averaged intracellular mechanical properties but also their spatial variations.

We have developed an experimental approach to map in a controlled manner, for the first time to our knowledge, the position-dependent rheological properties of the cell interior by combining optical tweezers-based active microrheology [START_REF] Guet | Mechanical role of actin dynamics in the rheology of the Golgi complex and in Golgiassociated trafficking events[END_REF] with soft lithography micropatterning. Recently experiments with atomic force microscopy (AFM) [START_REF] Hecht | Imaging viscoelastic properties of live cells by AFM: power-law rheology on the nanoscale[END_REF] and AFM or magnetic twisting cytometry combined with micropatterning [START_REF] Rigato | Atomic Force Microscopy Mechanical Mapping of Micropatterned Cells Shows Adhesion Geometry-Dependent Mechanical Response on Local and Global Scales[END_REF][START_REF] Park | Mapping the cytoskeletal prestress[END_REF] have been performed but did not discriminate between the mechanical contribution of the prestressed cell cortex and that of the cell interior. In contrast, our approach isolates the contribution of the cell interior. Plating cells on polarized crossbow-shaped adhesive micropatterns allows us to standardize the intracellular organization and reduce cell-to-cell variability [START_REF] Mandal | Cell dipole behaviour revealed by ECM subcellular geometry[END_REF][START_REF] Schauer | Probabilistic density maps to study global endomembrane organization[END_REF] and to quantify the spatial variations of the mechanical properties. Using power-law rheology, we analyze the relaxation of internalized microspheres following an imposed displacement in an optical trap and measure the local intracellular complex shear modulus G. The real part of G, the storage modulus G', is related to the elasticity of the bead microenvironment, while the imaginary part of G, the loss modulus G'', is related to its viscous-like behavior. We also use the more classical standard linear liquid (SLL) viscoelastic model [START_REF] Bausch | Measurement of local viscoelasticity and forces in living cells by magnetic tweezers[END_REF][START_REF] Chan | Myosin II Activity Softens Cells in Suspension[END_REF] to analyze the relaxation experiments in terms of viscoelastic components and compare the results with our power-law based model. The resulting maps characterize the spatial variations of the intracellular complex shear modulus and allow us to compare its distribution upon treatments with various drugs or inhibitors. We first map the shear modulus in retinal pigment epithelial (RPE1) cells and quantify the contribution to intracellular mechanics of the actin and microtubule cytoskeletons, ATP-dependent active processes and internal membranes. We further show that breast metastatic cells can be clearly distinguished from non-tumorigenic cells based on the spatial distribution of their mechanical properties.

Results

Microrheological measurements combined with micropatterning to map intracellular mechanics

Cells plated on adhesive micropatterned substrates are known to exhibit a standardized shape and intracellular organization [START_REF] Schauer | Probabilistic density maps to study global endomembrane organization[END_REF]. We took advantage of this property to map intracellular mechanics in cells plated on crossbow-shaped micropatterns (Fig. 1A). Optical tweezers-based microrheology was used to measure intracellular mechanics. Beads were first incubated with cells so that typically one or two beads were internalized by endocytosis and transported to various locations in the cell cytoplasm. After plating and full spreading of the cells on the micropatterns, a step stress was applied to an optically trapped bead by displacing the microscope stage (Fig. 1B). The viscoelastic relaxation of the bead towards the trap center was followed by tracking the position of the bead with time. The initial position of the bead relative to the crossbow micropattern was recorded. This procedure was repeated for a large number of cells in order to sample a maximal area on the micropattern and build the spatial maps (see Materials and Methods and Figs. 1C andS1).

In order to test the robustness of our mechanical maps, the relaxation experiments were analyzed by three independent methods: i) a phenomenological and model-independent approach, ii) a standard viscoelastic model, and iii) a specific analysis based on power law rheology. First we used a phenomenological and model-independent approach as described in [START_REF] Guet | Mechanical role of actin dynamics in the rheology of the Golgi complex and in Golgiassociated trafficking events[END_REF]. Qualitatively, in a rigid microenvironment, the initial bead displacement is large and its relaxation is weak. We measured two parameters to characterize the local intracellular rigidity of the bead microenvironment, the initial bead displacement (X b ) following the step displacement of the microscope stage (X s =0.5 µm) and the rigidity index (RI) defined as the ratio between the bead displacement (x b (t)) and the stage displacement integrated along the bead relaxation curve 𝑅𝐼 = ∫ 𝑥 𝑏 (𝑡) 𝑑𝑡 𝑇 0 𝑋 𝑠 𝑇 (Fig. S2A). The value of the bead step amplitude X b is always smaller than X s and the value of the rigidity index RI falls between 0 and 1. A rigid microenvironment is characterized by X b ≲ X s and RI≃1, while a soft microenvironment is characterized by X b ≪ X s and RI≃0. Although the two measured parameters do not directly relate to usual rheological parameters, they allow to compare the intracellular rigidity in various cellular contexts in a modelindependent manner [START_REF] Guet | Mechanical role of actin dynamics in the rheology of the Golgi complex and in Golgiassociated trafficking events[END_REF].

Second, we considered a simple three elements viscoelastic analytical model [START_REF] Bausch | Measurement of local viscoelasticity and forces in living cells by magnetic tweezers[END_REF] composed of a Kelvin-Voigt body (a dashpot and a spring in parallel) and a dashpot in series, known as the standard linear liquid (SLL) model (Fig. S2B and Supplementary Materials and Methods). The relaxation curve of the bead towards the trap center was fitted using this model to determine the local spring constant  and dashpot viscosity  of the bead microenvironment (Fig. S3 and Supplementary Materials and Methods) in the reference frame of the micropattern.

Third, and most importantly, we developed a theoretical framework based on power-law rheology.

We first calculated the Laplace transform 𝑥 ̃𝑏(𝑠) of the bead displacement x b (t) as in [START_REF] Desprat | Creep function of a single living cell[END_REF]. We found (see Materials and Methods): 𝑥 ̃𝑏(𝑠) = 𝑋 𝑠 𝜏𝑠 2 1-𝑒 -𝑠𝜏 1+𝜎 * 𝑠𝐽 ̃ for 𝑡 > 𝜏, where the stage step displacement is modelled by a linear increase from 0 µ𝑚 to 𝑋 𝑠 = 0.5 µ𝑚 between 𝑡 = 0 and 𝑡 = 𝜏 seconds, 𝐽 ̃ is the Laplace transform of the creep function 𝐽(𝑡) and 𝜎 * = 𝑘 𝑡𝑟𝑎𝑝 𝐿 𝑐 4𝜋𝑅 2 is a geometrical factor involving the trap stiffness 𝑘 𝑡𝑟𝑎𝑝 , the bead radius 𝑅 and a characteristic length 𝐿 𝑐 . Assuming that the creep function behaves as a power law with an exponent  and a prefactor A, 𝐽(𝑡) = 𝐴𝑡 𝛼 , we then used inverse Laplace transforms and gamma functions to obtain the bead displacement:

𝑥 𝑏 (𝑡) 𝑋 𝑠 = 𝑡 𝜏 ∑ [-Γ(1+𝛼)𝜎 * 𝐴𝑡 𝛼 ] 𝑛 Γ(2+𝑛𝛼) ∞ 𝑛=0 - 𝑡-𝜏 𝜏 ∑ [-Γ(1+𝛼)𝜎 * 𝐴(𝑡-𝜏) 𝛼 ] 𝑛 Γ(2+𝑛𝛼) ∞ 𝑛=0
. We used this equation, that we term power-law (PL) model in the following, to fit the relaxation of the bead position 𝑥 𝑏 (𝑡) as a function of time and measure the prefactor A and the exponent . The complex shear modulus G= G'+i G'' was then deduced from A and  (see Materials and Methods). The models described above both assume a linear mechanical response of the system. Therefore we checked that the experiments were performed in a linear regime by varying the step size linearly around 𝑋 𝑠 = 0.5 µ𝑚. We found that the bead step amplitude X b increased linearly while none of the mechanical parameters derived from either the SLL model or the power law model significantly differed (Fig. S2 C-E).

We first applied our active microrheology technique to map the complex shear modulus in retinal pigment epithelial (RPE1) cells. When RPE1 cells are plated on crossbow-shaped micropatterns, they adopt a polarized intracellular organization, with the microtubule organizing center and the Golgi apparatus positioned in front of the nucleus which is found on average in the center of the micropattern (Fig. S1B) [START_REF] Schauer | Probabilistic density maps to study global endomembrane organization[END_REF]. Three-dimensional imaging of beads internalized in micropatterned cells shows that the bead z-position (5.0±.4 µm) is well distributed within the maximum cell height (10.6 ±0.4 µm) (Fig. S2F-H). A striking feature of the relaxation curves is that they strongly depend on the bead position in the reference frame of the micropattern (Fig. 1D). Both the amplitude of the bead displacement and the dynamics of the relaxation varied with the bead position. To visualize these variations, we generated maps of the storage modulus G', of the loss modulus G'', and of the shear modulus G 0 defined as the 𝐺 0 = |𝐺| = |𝐺 ′ + 𝑖𝐺′′| obtained using the PL model (Fig. 1E). Although the three maps exhibited strong spatial fluctuations, they clearly showed that the values of the three parameters decrease from the cell center towards the cell periphery. Furthermore, the spatial distribution of shear modulus is polarized along the polarity axis of the crossbow-shaped micropattern (Fig. 1F), with moduli increasing along the cell rearto-front axis. This is consistent with the distribution of intracellular organelles and cytoskeletal elements reported previously in the same cell type [START_REF] Schauer | Probabilistic density maps to study global endomembrane organization[END_REF][START_REF] Théry | Anisotropy of cell adhesive microenvironment governs cell internal organization and orientation of polarity[END_REF]. We also note a slight left/right asymmetry in the distribution of the shear modulus (Figs. 1E-F). The map of the ratio G' /G'' (Fig. 1E) shows that elasticity contributes more to cell rigidity than viscosity in agreement with previous results [START_REF] Guo | Probing the Stochastic, Motor-Driven Properties of the Cytoplasm Using Force Spectrum Microscopy[END_REF][START_REF] Crocker | Multiple-particle tracking and two-point microrheology in cells[END_REF][START_REF] Weihs | Bio-microrheology: a frontier in microrheology[END_REF][START_REF] Cai | Quantifying cell-to-cell variation in power-law rheology[END_REF] and that the contribution of elasticity relative to viscosity increases towards the cell periphery. The results of the PL model were confirmed by generating maps of the mechanical parameters obtained with the modelindependent approach and the SLL model: the rigidity index RI, the bead step amplitude X b , the spring constant  and the viscosity  (Fig. S3). Finally, to validate the use of micropatterns, we show that the local values of the rigidity index are less dispersed in cells plated on crossbow-shaped micropatterns than in non-patterned cells (Fig. S3G).

Role of the actin and microtubule cytoskeletons in the spatial distribution of intracellular viscoelasticity

The cytoskeleton maintains cell architecture, generates forces and contractility and plays an active role in many cellular processes such as cell division or cell migration. Actin and microtubules have been shown to contribute to intra and extracellular mechanics [START_REF] Wen | Polymer physics of the cytoskeleton[END_REF][START_REF] Mandal | Cell dipole behaviour revealed by ECM subcellular geometry[END_REF][START_REF] Van Citters | The role of F-actin and myosin in epithelial cell rheology[END_REF][START_REF] Bursac | Cytoskeletal remodelling and slow dynamics in the living cell[END_REF][START_REF] Wakatsuki | Effects of cytochalasin D and latrunculin B on mechanical properties of cells[END_REF][START_REF] Wilhelm | Rotational magnetic endosome microrheology: Viscoelastic architecture inside living cells[END_REF]. To determine the contribution of the cytoskeletal networks to the local intracellular viscoelastic properties, we have focused on actin and microtubules and used drug treatments that affect their dynamics (Fig. S4).

When microtubules were depolymerized with nocodazole, the intracellular shear modulus slightly but significantly decreased compared to non-treated cells. Both the storage and the loss moduli decreased. In contrast, microtubule stabilization by taxol induced a strong increase in the shear modulus (Fig. 2A,B). To investigate the role of actin dynamics, we used two drugs that disrupt microfilament organization with different mechanisms of action, latrunculin A (LatA) which forms 1:1 complexes with monomeric G-actin and cytochalasin D (CytoD) which caps growing ends of microfilaments and disrupts actin organization by increasing the number of actin filament ends and forming aggregates of dense filamentous actin [START_REF] Cooper | New insights into mechanism and regulation of actin capping protein[END_REF]. In contrast with microtubules, we did not try stabilizing the actin cytoskeleton because jasplakinolide, the drug classically used to stabilize actin, blocks actin dynamics by disrupting microfilaments and inducing polymerization of monomeric actin into amorphous aggregates and its effects on intracellular mechanics are difficult to interpret [START_REF] Guet | Mechanical role of actin dynamics in the rheology of the Golgi complex and in Golgiassociated trafficking events[END_REF]. While LatA decreased the shear modulus through a strong decrease in both storage and loss moduli, as expected from previous reports [START_REF] Van Citters | The role of F-actin and myosin in epithelial cell rheology[END_REF][START_REF] Wakatsuki | Effects of cytochalasin D and latrunculin B on mechanical properties of cells[END_REF][START_REF] Wilhelm | Rotational magnetic endosome microrheology: Viscoelastic architecture inside living cells[END_REF], CytoD had no significant effect (Fig. 2A,B). Similar results were obtained with the two model-independent parameters, rigidity index and bead step amplitude, or with the spring constant and viscosity derived from the SLL model (Fig. S5A). To investigate further the differences between the effects of LatA and CytoD on intracellular mechanics, we have measured the diffusion coefficient of beads trapped by the optical tweezers in the presence or absence of LatA or CytoD (Fig. 2C). Consistently with the viscoelastic relaxation experiments, we found that beads diffused faster in cells treated with LatA while they diffused slightly slower in the presence of CytoD than in control cells.

We next built the spatial maps of the shear modulus (Fig. 3A) and the storage and loss moduli (Fig. S6), as well as the maps of the rigidity index and the bead step amplitude (Fig. S5B) and of the spring constant and viscosity of the SLL model (Fig. S5C). In the presence of the drugs, the maps showed a less regular pattern with more secondary maxima at the periphery of the micropattern. That all drug treatments induced mechanical spatial heterogeneities was further indicated by the shorter autocorrelation length of the drug-treated maps of the rigidity index compared to the control map (Fig. S5D).

Because our mapping method uses the center of the micropattern as a reference frame, we also asked how the intracellular rigidity spatially distributes relative to the nucleus. We plotted the variations of the shear modulus as a function of the distance separating the bead and the nuclear envelope (Fig. 3B, see Materials and Methods). In agreement with previous results showing that on average the nucleus is located in the center of the micropattern [START_REF] Schauer | Probabilistic density maps to study global endomembrane organization[END_REF] and consistent with the maps shown in Fig. 1E, we found that the shear modulus decreases with the distance to the nucleus in control conditions (Fig. 3B, left panel) yielding a negative shear modulus gradient (Fig. 3C). Both microtubule depolymerization by nocodazole and stabilization by taxol increased the shear modulus in close proximity to the nucleus, intensifying the negative gradient (Fig. 3B,C). In contrast, destabilizing actin with LatA induced a decrease in the shear modulus in the perinuclear region, which abolished the gradient. In the presence of CytoD, the shear modulus was more heterogeneously distributed between 2 µm and 6 µm from the nucleus, yielding a positive shear modulus gradient (Fig. 3B,C). As in the control situation, the G'/G'' ratio was greater than 1 for all tested perturbations of the actin or microtubule cytoskeleton, indicating that elastic storage contributes more to the shear modulus than viscous loss. However their relative contribution changed depending on the drug treatment as shown by the variations of the G'/G'' ratio (Fig. S5E, left panel) or by cross-correlation of the rigidity maps with the elasticity or viscosity maps (Fig. S5E, right panel). Nocodazole and CytoD increased the predominance of elastic storage over viscous loss, while taxol and LatA decreased it slightly. Together these results not only confirm the contribution of microfilaments and microtubules to cell mechanics, but also show that they control the spatial distribution of intracellular mechanical parameters.

Impact of intracellular membranes

Surprisingly, the contribution of intracellular membranes to cell mechanics has never been quantified to our knowledge. We have used brefeldin A (BFA), an inhibitor of the Arf1 exchange factor GBF1, to investigate the mechanical role of intracellular membranes. By blocking COPI-dependent membrane trafficking, BFA disperses Golgi membranes which first form long tubules before being absorbed in the endoplasmic reticulum [START_REF] Lippincott-Schwartz | Rapid redistribution of Golgi proteins into the ER in cells treated with brefeldin A: Evidence for membrane cycling from Golgi to ER[END_REF]. BFA was also shown to induce membrane tubulation of endosomes and lysosomes [START_REF] Lippincott-Schwartz | Brefeldin A's effects on endosomes, lysosomes, and the TGN suggest a general mechanism for regulating organelle structure and membrane traffic[END_REF]. The effect of BFA is reversible and when BFA is washed-out the Golgi apparatus reassembles.

We first asked whether perturbing intracellular membranes using BFA could affect viscoelasticity of the cytoplasm. The live Golgi lipid marker BodipyTR-ceramide was used to visualize Golgi membranes.

We found that BFA had minor effects on the relaxation curves (Fig. 4A). The bead displacement was shifted to larger values, indicating a small increase in intracellular rigidity. Accordingly, the shear modulus deduced from the PL model slightly but significantly increased (Fig. 4B). The mechanical parameters deduced from the model-independent approach (Fig. S7A) or the SLL model (Fig. S7B) all increased but the only statistically significant effect we measured was an increase in viscosity. After wash-out and recovery, the values decreased back to levels close to the control. Similarly, BFA did not induce any major difference on the spatial maps (Fig. 4C and Fig. S7C-E) except for an increase in G 0 3-6 µm away from the nucleus leading to an increase in the shear modulus gradient (Fig. 4D). These results point to a small but significant contribution of internal membranes to intracellular mechanics.

Energy depletion stiffens the cell cytoplasm

Several energy-dependent cellular processes, including cytoskeletal dynamics and molecular motor activity, can affect intracellular mechanics. It was shown recently using force spectrum microscopy that non-equilibrium ATP-driven fluctuating forces dominate thermal forces on time-scales greater than one second and enhance the movement of intracellular particles [START_REF] Guo | Probing the Stochastic, Motor-Driven Properties of the Cytoplasm Using Force Spectrum Microscopy[END_REF]. We have used ATP depletion to study the effects of active forces on intracellular viscoelasticity with our microrheology technique. The viscoelastic relaxation was much slower after ATP depletion (Fig. 5A) indicating a stiffening of the cytoplasm. Consistently, the shear modulus was significantly higher in ATP depleted cells, mostly due to an increase in the storage modulus G' (Fig. 5B). The model-independent rigidity index and bead step amplitude (Fig. S8A) and the spring constant and viscosity from the SLL model experienced a similar increase (Fig. S8B). No significant changes in the microtubule and actin organization were observed upon ATP depletion (Fig. S8E) suggesting that the effects on intracellular rigidity were not indirectly due to effects on the cytoskeleton. Spatial mapping evidenced a striking increase in the shear modulus at the center of the micropattern due to ATP depletion (Fig. 5C and Fig. S8C) that yielded a sharper shear modulus gradient than in control cells (Fig. 5D). Accordingly we also detected an increase in the perinuclear region of the four other mechanical parameters we measured (Fig. S8D). These results show that active ATP-dependent forces soften the cell interior, especially at the cell center.

Spatially resolved microrheology to differentiate normal and cancer cells

Several recent reports indicate that the cytoplasm of metastatic cells is softer than the cytoplasm of nontumorigenic cells [START_REF] Li | Intracellular particle tracking as a tool for tumor cell characterization[END_REF][START_REF] Guo | Probing the Stochastic, Motor-Driven Properties of the Cytoplasm Using Force Spectrum Microscopy[END_REF][START_REF] Agus | A physical sciences network characterization of non-tumorigenic and metastatic cells[END_REF][START_REF] Coughlin | Changes in cytoskeletal dynamics and nonlinear rheology with metastatic ability in cancer cell lines[END_REF][START_REF] Guck | Optical deformability as an inherent cell marker for testing malignant transformation and metastatic competence[END_REF]. We have used our technique to differentiate between non-tumorigenic breast cells (MCF-10A) and highly metastatic breast cancer cells (MDA-MB-231). Averaged data showed that the shear modulus (Fig. 6A,B) and all the other measured mechanical parameters (Fig. S9A,B) were significantly lower in the case of metastatic cells, with up to a three-fold decrease depending on the parameter. The spatial maps obtained from metastatic cells further showed a global reduction in the shear modulus (Fig. 6C) and the other parameters (Fig. S9C,D) throughout the micropattern surface. This difference was not due to geometrical effects as the averaged cell volume was not significantly different between the two cell lines (1905±185 µm 3 and 2015±125 µm 3 for MCF-10A and MDA-MB-231 respectively with a p-value equals to 0.3, as measured from bright field images of cells in suspension; see also [START_REF] Agus | A physical sciences network characterization of non-tumorigenic and metastatic cells[END_REF]). The spatial maps also exhibited a similar smoother spatial distribution than in non-tumorigenic cells (Figs. 6C and S9C,D). Accordingly, the auto-correlation lengths of the spring constant and the viscosity form the SLL model were larger in metastatic cells (Fig. S9E). Measurements of the G'/G'' ratio and cross-correlation of the rigidity map with the spring constant and viscosity maps showed that the relative contribution of elasticity decreases in metastatic cells (Fig. S9F). Finally, in the reference frame of the nucleus, the shear modulus was always smaller in metastatic cells, but the difference was more pronounced towards the cell periphery due to opposite shear modulus gradients (Fig. 6D). These results suggest that, while cancer progression decreases both elastic storage and viscous loss, the balance shifts towards viscosity and the spatial distribution of the mechanical parameters gets more uniform.

We next asked whether the differences in intracellular mechanics between non-tumorigenic and metastatic cells could correlate with the spatial organization of the actin and/or microtubule cytoskeletons.

In MCF-10A non-tumorigenic cells, a higher microtubule density was observed in the cell center in a perinuclear area and actin stress fibers were visible throughout the cell cytoplasm. In contrast, in MDA-MB-231 metastatic cells, the microtubule cytoskeleton appeared more uniformly spread and less dense and microfilaments mostly localized at the cell periphery (Fig. 6E andS9G). Accordingly, and as expected from the results in RPE1 cells, depolymerizing actin with latrunculin A in MCF-10A cells decreased intracellular rigidity, while stabilizing microtubules with taxol in MDA-MB-231 cells increased rigidity (Fig. S9H). These findings suggest that the mechanical differences observed between non-tumorigenic and metastatic cells could arise from their different organization of the cytoskeleton and its associated signaling.

Discussion

We have developed an approach to map the spatial distribution of the intracellular mechanical properties of the cytoplasm of living cells. Our technique allows to measure cell-averaged parameters that can be compared with those obtained with more classical approaches, but also reveals the spatial variations of intracellular mechanics. Previous work has either studied local mechanics without controlling the cell geometry and thus with large cell-to-cell variability [START_REF] Berret | Local viscoelasticity of living cells measured by rotational magnetic spectroscopy[END_REF][START_REF] Tseng | Micromechanical mapping of live cells by multiple-particle-tracking microrheology[END_REF][START_REF] Hecht | Imaging viscoelastic properties of live cells by AFM: power-law rheology on the nanoscale[END_REF] or measured regional rigidity variations of the cell surface in controlled geometries using atomic force microscopy or magnetic twisting cytometry [START_REF] Cai | Quantifying cell-to-cell variation in power-law rheology[END_REF][START_REF] Rigato | Atomic Force Microscopy Mechanical Mapping of Micropatterned Cells Shows Adhesion Geometry-Dependent Mechanical Response on Local and Global Scales[END_REF][START_REF] Park | Mapping the cytoskeletal prestress[END_REF]. Central to our method is the use of adhesive micropatterns to standardize the intracellular organization. By combining micropatterning with optical tweezers active microrheology, we are able to spatially resolve the local mechanical properties of the cell interior. As expected, plating cells on micropatterns decreases the variability in the local rigidity measurements (Fig. S3G). In our viscoelastic relaxation experiments, neither the stress nor the strain is constant, in contrast to classical creep experiments [START_REF] Desprat | Creep function of a single living cell[END_REF]. The analysis of the experimental data thus required a unique theoretical framework. The simple analytical standard linear liquid (SLL) viscoelastic model we use to describe the intracellular microenvironment fits the experimental data remarkably well. However models based on a finite number of viscoelastic elements and corresponding relaxation times are no longer considered to be well adapted to describe cell mechanics and many more elaborate theoretical models have been proposed, including glassy rheology [START_REF] Bursac | Cytoskeletal remodelling and slow dynamics in the living cell[END_REF][START_REF] Mandadapu | On the cytoskeleton and soft glassy rheology[END_REF], poroelasticity [START_REF] Moeendarbary | The cytoplasm of living cells behaves as a poroelastic material[END_REF] or tensegrity [START_REF] Ingber | Tensegrity, cellular biophysics, and the mechanics of living systems[END_REF]. A growing number of reports point to models based on power-law rheology as most relevant [START_REF] Guo | Probing the Stochastic, Motor-Driven Properties of the Cytoplasm Using Force Spectrum Microscopy[END_REF][START_REF] Fabry | Scaling the microrheology of living cells[END_REF][START_REF] Cai | Quantifying cell-to-cell variation in power-law rheology[END_REF][START_REF] Balland | Power laws in microrheology experiments on living cells: Comparative analysis and modeling[END_REF][START_REF]A comparative mechanical analysis of plant and animal cells reveals convergence across kingdoms[END_REF][START_REF] Hoffman | The consensus mechanics of cultured mammalian cells[END_REF][START_REF] Guo | The role of vimentin intermediate filaments in cortical and cytoplasmic mechanics[END_REF]. Consequently we also developed a specific theoretical analysis based on power-law (PL) rheology (see Materials and Methods) to interpret our experiments. In our PL-based analysis, only two fit parameters are used, the exponent and the prefactor of the underlying power-law. In contrast three fit parameters are used in the SLL model.

Despite this, the PL model fits the experimental data equally well, strengthening the idea that power-law rheology may also govern intracellular mechanics. The order of magnitude of the mechanical parameters derived from both models are consistent since the spring constant 𝜇 from the SLL model can be related to the elastic storage modulus

𝐺 ′ using 𝐺 ′ = 𝜇 6𝜋𝑅
, where 𝑅=1 µm is the bead radius. With the typical values we measured for the spring constant 𝜇~200 pN/µm, we obtain 𝐺′~10 Pa in good agreement with the values of G' obtained with the PL model and with values reported previously for the cytoplasm [START_REF] Guo | Probing the Stochastic, Motor-Driven Properties of the Cytoplasm Using Force Spectrum Microscopy[END_REF][START_REF] Berret | Local viscoelasticity of living cells measured by rotational magnetic spectroscopy[END_REF][START_REF] Cai | Quantifying cell-to-cell variation in power-law rheology[END_REF].

Besides quantifying the complex shear modulus with the PL model and the spring constant and the viscosity with the SLL model, we also introduce the rigidity index, a local indicator of intracellular rigidity independent of any underlying model. The rigidity maps provide a convenient tool to visualize the local variations in intracellular mechanics.

Spatial mapping of the shear modulus clearly depicts a decrease in rigidity from the center of the micropattern towards its periphery in RPE1 cells and a polarization along the axis of the crossbow micropattern (Fig. 1) reflecting the polarized organization of the cytoskeleton and intracellular compartments. The central region of the micropattern statistically corresponds to the perinuclear region, a zone particularly rich in internal membrane compartments such as the endoplasmic reticulum (ER), the Golgi apparatus and endosomes and also dense in cytoskeletal fibers. The shear modulus decreases between 1 µm and 2.5 µm away from the nucleus edge then stabilizes (Fig. 3B). Interestingly this decrease correlates with a decrease in the density of vimentin intermediate filaments, but not with the density of ER membranes (Fig. S10), pointing to a role of intermediate filaments in regulating the local rigidity in close vicinity to the nucleus. We have shown previously that intracellular rigidity increases close to the Golgi apparatus [START_REF] Guet | Mechanical role of actin dynamics in the rheology of the Golgi complex and in Golgiassociated trafficking events[END_REF] due to the actin cytoskeleton. In good agreement, we find here that actin depolymerization by LatA induces a local decrease in the shear modulus 1-2 µm away from the nucleus (Fig. 3). Surprisingly, dispersing Golgi membranes with BFA had only a weak effect at the center (Fig. 4),

suggesting that Golgi associated actin may remain present on BFA-insensitive Golgi remnants or Golgi matrix components [START_REF] Nakamura | Characterization of a cis-Golgi matrix protein, GM130[END_REF][START_REF] Hendricks | Golgi proteins persist in the tubulovesicular remnants found in brefeldin A-treated pancreatic acinar cells[END_REF]. BFA most significant effect was to increase the shear modulus, and especially viscous loss, closer to the cell periphery, about 3-6 µm away from the nucleus (Fig. 4C,D and S7). This may be due to endosome tubulation as reported upon treatment with BFA [START_REF] Lippincott-Schwartz | Brefeldin A's effects on endosomes, lysosomes, and the TGN suggest a general mechanism for regulating organelle structure and membrane traffic[END_REF].

Interestingly, we note that the mechanical maps of non-treated RPE1 cells (Fig. 1E) and the distribution of the shear modulus (Fig. 1F) are not strictly symmetric with respect to the polarity axis of the crossbow micropattern. The averaged images of actin, microtubules and vimentin intermediate filaments (Figs. S4 andS10B) show that the cytoskeleton too lacks perfect symmetry. Similarly, when we examined the position of the nucleus carefully, we found that on average the nucleus position is slightly shifted towards the upper right quadrant of the micropattern (Fig. S1B). This left/right asymmetry may be related to an intrinsic actin-driven chirality reported recently in micropatterned cells [START_REF] Tee | Cellular chirality arising from the self-organization of the actin cytoskeleton[END_REF].

To further take advantage of the spatial information provided by our mapping technique, we have used image correlation tools to quantify the spatial variations of the maps and the spatial differences between maps. This approach is relevant because each data point in the map actually probes a relatively large volume in the cytoplasm, since the 2 µm diameter bead is displaced by 0.5 µm. Using large beads allows to integrate the mechanics at all length scales below the bead size, while smaller submicrometric probes may omit essential mechanical features originating from micron-scale cytoplasmic structures such as intracellular membrane compartments or cytoskeletal fibers. One striking example is given by the antagonistic effects of depolymerizing or stabilizing microtubules. The averaged viscoelastic measurements show that depolymerizing microtubules with nocodazole decreases rigidity, as previously reported ( 47) while stabilizing microtubules with taxol increases rigidity (Figs. 2 andS5). The analysis of the spatial maps gives insights into the local effects of both types of perturbations (Figs. 3 andS5). When microtubules are depolymerized, rigidity increases in the perinuclear region and drops everywhere else, yielding the lowest auto-correlation length and the steepest shear modulus gradient of all the conditions we tested. Since no obvious effect of nocodazole on the actin cytoskeleton was observed (Fig. S4), the increase in perinuclear rigidity may be attributed to intermediate filaments which are known to contribute to cell rigidity (for a recent review see ( 58)) and to accumulate around the nucleus when microtubules are depolymerized [START_REF] Dupin | Cytoplasmic intermediate filaments mediate actin-driven positioning of the nucleus[END_REF][START_REF] Gurland | Stable, detyrosinated microtubules function to localize vimentin intermediate filaments in fibroblasts[END_REF]. In contrast, when microtubules are stabilized, the rigidity level increases, especially towards the pattern periphery, probably due to the presence of stable microtubule plus-ends at the cell periphery. Similarly, the local contribution of actin in intracellular rheology can be investigated in details using our method, as best illustrated by the differences we observed between the effects of two actin destabilizing drugs, LatA and CytoD. Although actin destabilization is mostly thought to soften the cell interior, its precise effect on intracellular mechanics is still not clear [START_REF] Van Citters | The role of F-actin and myosin in epithelial cell rheology[END_REF] and probably depends on the cell type [START_REF] Grady | Cell elasticity with altered cytoskeletal architectures across multiple cell types[END_REF]. Surprisingly, while in LatA treated cells the shear modulus dropped and bead diffusion was facilitated, we did not detect any effect upon treatment with CytoD (Fig. 2). Moreover, the two drugs had contrasting effects on the spatial distribution of the shear modulus. Maps of LatA treated cells were smoother, with a longer auto-correlation length and a very low shear modulus gradient (Figs. 3 andS5).

These differences probably reflect the differences in the action mechanisms of these two drugs [START_REF] Wakatsuki | Effects of cytochalasin D and latrunculin B on mechanical properties of cells[END_REF][START_REF] Cooper | New insights into mechanism and regulation of actin capping protein[END_REF].

A crosstalk between actin, microtubules and intermediate filaments could also be involved since actin destabilization affects the distribution of both microtubules (LatA in Fig. S4) and intermediate filaments in the perinuclear region [START_REF] Dupin | Cytoplasmic intermediate filaments mediate actin-driven positioning of the nucleus[END_REF][START_REF] Gruenbaum | Intermediate filaments: a dynamic network that controls cell mechanics[END_REF].

There has recently been a growing interest in the role of active cellular forces, such as forces generated by cytoskeletal dynamics, molecular motors or acto-myosin contractility, in intracellular mechanics and conflicting results have been reported [START_REF] Guo | Probing the Stochastic, Motor-Driven Properties of the Cytoplasm Using Force Spectrum Microscopy[END_REF][START_REF] Chan | Myosin II Activity Softens Cells in Suspension[END_REF][START_REF] Van Citters | The role of F-actin and myosin in epithelial cell rheology[END_REF][START_REF] Mizuno | Nonequilibrium mechanics of active cytoskeletal networks[END_REF][START_REF] Trepat | Universal physical responses to stretch in the living cell[END_REF][START_REF] Hale | Resolving the role of actoymyosin contractility in cell microrheology[END_REF]. Here we measure an unambiguous stiffening of the cytoplasm of RPE1 cells upon energy depletion, reflected by an elevated shear modulus localized at the cell center up to 5 µm away from the nucleus (Figs. 5 andS8). A likely explanation is that energy depletion inhibits cytoskeletal dynamic remodeling resulting in more crosslinked networks [START_REF] Mizuno | Nonequilibrium mechanics of active cytoskeletal networks[END_REF].

Finally, our method allows to clearly differentiate cancer cells from normal cells, based not only on their averaged intracellular mechanical properties, but also on the spatial distribution of their intracellular mechanical parameters. As previously reported for several types of cancers (6-9, 12, 51), the cytoplasm of metastatic breast cancer cells is softer than their non-tumorigenic counterparts. Our method, in contrast with previous studies, yields highly statistically significant decreases in both elasticity and viscosity (Fig. 6). The additional information provided by the spatial mapping on adhesive micropatterns further shows that maps of metastatic cells are smoother, more auto-correlated and comparatively more dominated by viscous loss than non-tumorigenic ones. The organization of both actin and microtubule cytoskeletons differs markedly in both cell types, which could at least in part explain the differences in intracellular mechanics (Fig. S9H). In cancer cells, microtubules appear more homogeneously distributed and less bundled and actin fibers are less abundant in the cell center and accumulate at the cell periphery (Figs. 6 andS9). Extrapolating our results on the mechanical effects of energy depletion, softening of cancer cells could also originate from a higher out-of-equilibrium activity. Our study emphasizes the need for spatial information in intracellular mechanics. Standardizing cell shape by micropatterning is a key aspect of our work. Because intracellular mechanical maps can detect subtle differences in the spatial distribution of mechanical parameters even in the absence of any change in their average values, our approach should provide insights into cancer cell mechanics.

Materials and Methods

Cell culture, micropatterning and reagents

Retinal pigment epithelial (RPE-1) cells were grown in DMEM/F12 medium supplemented with 10% fetal bovine serum (FBS) at 37ºC with 5% CO 2 . MCF-10A mammary epithelial cells were grown in DMEM/F12 medium supplemented with 5% horse serum, 20 ng/ml EGF, 100 ng/ml cholera toxin, and 10 µg/ml insulin. MDA-MB-231 metastatic breast cancer cells were grown in DMEM supplemented with 10% FBS (vol/vol). 900±30 µm 2 crossbow-shaped (18 µm radius circumscribing circle) micropatterns were printed on PEG-coated glass coverslips by deep UV photolithography then coated with 50 μg /ml fibronectin and 20 μg/ml Alexa546-fibrinogen (Sigma). For MCF-10A and MDA-MB-231 cells, 50 µg/ml concanavalin A was added to facilitate cell adhesion to the micropatterns. Cells were seeded on freshly prepared protein coated micropatterns and allowed to spread for at least 2-3hrs. Non-adherent cells were washed off by rinsing with culture medium. For rheological measurements cells were incubated overnight with 2 μm diameter fluorescently labeled latex beads (660/690 fluorescence, Bangs Laboratories, Inc) before seeding them on the micropatterns. 10 mM Hepes was added before microrheology experiments.

For immunofluorescence experiments, cells were fixed with 4% paraformaldehyde or with methanol. Primary antibodies were anti-β tubulin (Sigma #T4026), anti-GM130 (BD Transduction laboratories) and anti-vimentin (Santa Cruz). Actin was labelled using Alexa546-phalloidin (Invitrogen).

Secondary antibodies were from Jackson ImmunoResearch Laboratories. The live Golgi marker BodipyFL-ceramide (Invitrogen) was used to label the Golgi complex. Cells were incubated with 5 µM BodipyFL-ceramide for 15-30 min at 4ºC, washed with cold HBSS (Hank's Balanced Salt Solution, Invitrogen) then with warm culture medium and incubated at 37 ºC with 5% CO 2 for a further 30 mins. The endoplasmic reticulum (ER) was labeled using ER Tracker (Molecular Probes/ Thermo Fisher). Cellular membranes were stained in live cells using CellMask (Molecular Probes/Thermo Fisher).

To destabilize the cytoskeleton, cells were incubated for 30 min with 2 µM nocodazole, 2 µM taxol, 50 nM latrunculin A (LatA), or 200 nM cytochalasin D (CytoD). To perturb internal membranes, Brefeldin A (BFA) was incubated at 3 μg/ml for 15-30 mins before microrheological measurements in cells previously labeled using BodipyFL-ceramide (see above). After the rheological measurements in the presence of BFA, cells were washed with culture medium and were allowed to recover for 1 hr before performing rheological measurements again in the absence of BFA. Control cells were also labeled with bodipyFL-ceramide. To deplete ATP, cells were incubated with 10 mM deoxyglucose and 10 mM sodium azide (NaN 3 ) in Ca 2+ and Mg 2+ enriched PBS for 30 mins. Control cells were incubated with 5 mM glucose for 30 mins.

Optical tweezers based intracellular microrheology, image acquisition and data processing

The set-up combining optical tweezers and fast confocal microscopy was described in details previously [START_REF] Guet | Mechanical role of actin dynamics in the rheology of the Golgi complex and in Golgiassociated trafficking events[END_REF]. Briefly, a single fixed optical trap was built by coupling a 1060-1100 nm infra-red laser beam (2W maximal output power, IPG Photonics) to the back port of an inverted Eclipse microscope (Nikon) equipped with a resonant laser confocal A1R scanner (Nikon), a 37°C incubator and a nanometric piezostage (Mad City Labs). Coverslips with micropatterned cells containing typically 1 to 3 internalized beads were mounted in a Ludin chamber. A bead was first trapped by the optical tweezers (1 W laser output power, corresponding to 150 mW on the sample). An automated 0.5 µm step displacement was applied to the cell via the nanometric piezostage. Following the cell displacement, the bead was moved out of the trap center due to the viscoelastic nature of its microenvironment. The initial bead displacement, which we termed bead step amplitude and denoted with X b , was smaller than the stage displacement X s =0.5 µm. Because the optical tweezers position is fixed and acts as a spring on the bead, the bead position x b (t) relaxes from its maximal position X b towards the center of the optical trap. The bead relaxation was recorded by acquiring images at 120 frames/sec in the resonant mode of the confocal scanner for 40 s using the NIS Nikon software.

Images of fixed cells were acquired as z-stacks separated by 0.5 µm in the galvanometric mode of the confocal scanner with an 8 frames averaging using the NIS Nikon software. For single cell images, the z-stack was projected as a maximum intensity projection using Image J. For averaged cell images, zstacks were projected as average intensity projections using Image J then automatically aligned and averaged with a custom-written Image J macro. The averaged images are displayed using a fire look up table scaled between the minimum and the maximum intensity levels to highlight intensity variations.

Power law model

The relationships between the stress  and the strain  exerted on the bead microenvironment define the relaxation function G(t) and the creep function J(t) of the cytoplasm. If the initial strain is zero, then:

𝜎(𝑡) = ∫
𝑥 𝑏 (𝑡) = 𝑋 𝑠 [ 𝑡 𝜏 ∑ [-Γ(1+𝛼)𝜎 * 𝐴𝑡 𝛼 ] 𝑛 Γ(2+𝑛𝛼) ∞ 𝑛=0 - 𝑡-𝜏 𝜏 ∑ [-Γ(1+𝛼)𝜎 * 𝐴(𝑡-𝜏) 𝛼 ] 𝑛 Γ(2+𝑛𝛼) ∞ 𝑛=0
] (Eq. 3).

We use the first four terms of Eq. 3, that we term power-law (PL) model, to fit the first ten seconds of the relaxation of the bead position 𝑥 𝑏 (𝑡) as a function of time knowing 𝑋 𝑠 = 0.5 µ𝑚 and 𝜏 = 30 ms .The prefactor A and the exponent  are the only two fit parameters. Because at longer times bead tracking is more noisy due to intracellular dynamics even when cells are plated on micropatterns [START_REF] Tee | Cellular chirality arising from the self-organization of the actin cytoskeleton[END_REF], the data was fitted using the following weights: 0.9, 0.09 and 0.01 for 0<t<2s, 2<t<5s and 5<t<10s respectively).

The complex shear modulus G=G'+iG'' is then deduced from A and  as in [START_REF] Desprat | Creep function of a single living cell[END_REF] using

𝐺 ′ = (2𝜋) 𝛼 cos (𝛼

Spatial mapping

To generate the spatial maps, the two model-independent parameters (the rigidity index (RI) and the bead step amplitude (X b )), the spring constant  and the viscosity  deduced from the fits with the SLL model and the complex shear modulus (its value G 0 , its real part G' and its imaginary part G'') deduced from the fits with the PL model were measured at different bead locations within the cell. For each cell, one or two beads were used. A bright field image (differential interference contrast, DIC) and a confocal fluorescence image were taken simultaneously to visualize the crossbow micropattern labelled with Alexa546fibrinogen, the fluorescent bead and the nucleus. From these images, the X-Y position of the bead relative to center of the micropattern was measured. The distance between the bead and the nucleus was defined as the shortest distance between the center of the bead and the nuclear envelope. Each rheological parameter was plotted as a function of its X-Y position relative to the center of the micropattern using the 'contour' plot type in the OriginPro software (Origin Lab Corporation). Spatial maps were smoothed using 18 contour levels and the 'smoothing' function in OriginPro with the following parameter: 'total points increase factor'=500, 'smoothing parameter'=0.005). The convex hull of the set of bead X-Y coordinates was used as the map boundary. Data from N=19-166 beads, depending on the experimental conditions, were pooled to sample a maximum surface area on the micropattern. The coordinates of the most peripheral beads differed between experimental conditions or cell treatments, leading to different shapes of the map boundary. The polar plots were obtained by averaging the data along the polar radius in 30° intervals. The plots of the shear modulus as a function of the distance to the nucleus were obtained by averaging the data along all directions in 1 µm intervals from the nucleus.

Statistical analysis

Results are representative of at least three independent experiments. Data are expressed as means ± standard error of the mean. Statistical relevance was evaluated using Student's t-tests and the P value is indicated (n.s, nonsignificant; *P < 0.01, **P < 0.001, and ***P < 0.0001, except when indicated).

Supplementary Materials and Methods

Standard linear liquid (SLL) model

The microenvironment of the bead was described by a three-element viscoelastic model composed of a Kelvin-Voigt body (a spring of spring constant  and a dashpot of viscosity  in parallel) and a dashpot of viscosity  0 in series (Fig. S2b). The forces acting on the bead after the step displacement of the stage X s are due to the viscoelastic response of the microenvironment and to the optical trap that acts as a spring of stiffness k trap . The trap stiffness was calibrated using Stokes law in water and was 214.4±2.7 pN.µm.W -1

. Neglecting inertia in Newton's law yields the following equation verified by 𝑥 𝑏 (𝑡) :

𝑘 𝑡𝑟𝑎𝑝 𝑥 𝑏 + 6𝜋𝑅𝜂(𝑥̇𝑏 -𝑥̇1) + 𝜇(𝑥 𝑏 -𝑥 1 ) = 0 (Eq. 1),
where R=1 µm is the bead radius and x 1 is the position of the junction between the Kelvin-Voigt body and the dashpot in series relative to the trap center. Stress continuity between the Kelvin-Voigt body and the dashpot in series imposes:

6𝜋𝑅𝜂 0 (𝑥̇1 -𝑥̇𝑠) = 6𝜋𝑅𝜂(𝑥̇𝑏 -𝑥̇1) + 𝜇(𝑥 𝑏 -𝑥 1 ),
where x s is the position of the stage relative to the trap center. Taking 𝑥 𝑠 (𝑡) = 𝑋 𝑠 = 𝑐𝑠𝑡 for 𝑡 > 0 to model the stage step displacement, we find:

6𝜋𝑅𝜂 0 𝑥̇1 = -𝑘 𝑡𝑟𝑎𝑝 𝑥 𝑏 (Eq. 2).
Combining equations ( 1) and ( 2) gives the equation verified by 𝑥 𝑏 (𝑡) :

𝑥 𝑏 ̈+ (𝛼 + 𝛽 + 𝛾)𝑥 𝑏 ̇+ 𝛼𝛽𝑥 𝑏 = 0 (Eq. 3)

with 𝛼 = 𝑘 𝑡𝑟𝑎𝑝 /(6𝜋𝑅𝜂 0 ), 𝛽 = 𝜇 /(6𝜋𝑅𝜂) and 𝛾 = 𝑘 𝑡𝑟𝑎𝑝 /(6𝜋𝑅𝜂). Solving equation 3 with the initial conditions 𝑥 𝑏 (𝑡 = 0) = 𝑋 𝑏 and 𝑥 1 (𝑡 = 0) = 𝑋 𝑠 gives 𝑥 𝑏 (𝑡) = 𝑋 0 𝑒 -𝑡/𝜏 +𝑋 0 ′ 𝑒 -𝑡/𝜏′ (Eq. 4), with the two relaxation times  and ' satisfying (𝜏 -1 )(𝜏 ′ -1 ) = 𝛼𝛽 and 𝜏 -1 + 𝜏 ′ -1 = 𝛼 + 𝛽 + 𝛾.

The bead displacement relative to the trap center 𝑥 𝑏 (𝑡) was measured as a function of time using a homemade single particle tracking Matlab routine (kindly provided by Gil Toombes, National Institutes of Health, Bethesda, and Patricia Bassereau, UMR168 CNRS-Institut Curie, Paris) and fitted for T=10 s with equation 4 to yield the four parameters X 0 , X 0 ', and '. From these values, the initial bead step amplitude X b and the three parameters ,, and  were then calculated using the initial conditions:

𝑋 𝑏 = 𝑋 0 + 𝑋 0 ′ and 𝛽 = 1 𝑋 𝑠 ( 𝑋 𝑏 -𝑋 0 𝜏 + 𝑋 𝑏 -𝑋 0 ′ 𝜏′
). Knowing ,, and we deduced the rheological parameters 𝜇, 𝜂 and 𝜂 0 . Note that if the dashpot in series 𝜂 0 is removed (𝜂 0 → ∞), we obtain a single exponential (𝜏 ′-1 → 0) as in Ref. [START_REF] Guet | Mechanical role of actin dynamics in the rheology of the Golgi complex and in Golgiassociated trafficking events[END_REF] in the main text.

Correlation analysis

Contour maps were first converted to grey level images. Auto-correlation and cross-correlation were performed using a custom written Matlab routine based on the corr2 and xcorr2 functions available from the Signal Processing Toolbox™. The auto-correlation length of a given map A was defined as the width of the best Gaussian fit to the auto-correlation image obtained by cross-correlating A with itself. Crosscorrelation between two maps A and B was defined as the normalized ratio

Mean squared displacement analysis

Beads diffusing in the optical trap were tracked using single particle tracking. The two-dimensional mean squared displacement (MSD) was calculated from the X-Y coordinates of the bead and plotted as a function of the time interval ∆𝑡. A restricted diffusion in a 'cage' model 46

Figure S8

Energy depletion stiffens the cell interior and perturbs the spatial distribution of intracellular viscoelasticity.

Figure S10

(A) Correlation between the shear modulus (Fig. 3B) and the density of vimentin intermediate filaments (green) or the density of membranes of the endoplasmic reticulum (ER, red) as a function of the distance to the nucleus in RPE1 cells plated on crossbow micropatterns. Cells were fixed and stained for vimentin using an anti-vimentin antibody or stained for ER membranes using a live cell ER marker then fixed.

., 2 1 𝑠 2

 212 In the case of our optical tweezers relaxation experiment the stress is given by 𝜎(where x b and x s are the bead and the stage displacements relative to the trap center respectively, R is the bead radius, k trap is the trap stiffness and L c is length characteristic of the deformation of the bead microenvironment. Replacing in the expression of the stress, we find: is a geometrical factor. Modelling the stage step displacement by a linear increase from 0 µ𝑚 to 𝑋 𝑠 = 0.5 µ𝑚 between 𝑡 = 0 and 𝑡 = 𝜏 seconds and taking the Laplace transform (𝐿𝑇(𝑓(𝑡)) = 𝑓 ̃(𝑠) = ∫ 𝑓(𝑡)𝑒 -𝑠𝑡 ∞ 0 𝑑𝑡) of Eq.1 as in (42), we find 𝑥 ̃𝑏(𝑠) ̃ for 0 ≤ 𝑡 ≤ 𝜏 and 𝑥 ̃𝑏(𝑠) = 𝑋 𝑠 𝜏𝑠 2 1-𝑒 -𝑠𝜏 1+𝜎 * 𝑠𝐽 ̃ for 𝑡 > 𝜏, where 𝑥 ̃𝑏(𝑠) is the Laplace transform of the bead displacement x b (t) and 𝐽 ̃(𝑠) is the Laplace transform of the creep function 𝐽(𝑡) (𝐽 ̃(𝑠) = 𝐺 ̃(𝑠) with 𝐺 ̃(𝑠) the Laplace transform of the relaxation function 𝐺(𝑡)). Developing in powers of 𝜎 * 𝑠𝐽 ̃ gives: 𝑥 ̃𝑏(𝑠) = 𝑋 𝑠 𝜏𝑠 2 (1 -𝑒 -𝑠𝜏 ) ∑ [-𝜎 * 𝑠𝐽 ̃]𝑛 ∞ 𝑛=0 (Eq. 2). Next we assume that the creep function behaves as a power law with an exponent  and a prefactor A, 𝐽(𝑡) = 𝐴𝑡 𝛼 (42). The Laplace transform of 𝐽(𝑡) is thus 𝐽 ̃(𝑠) = 𝐴Γ(1 + 𝛼)/𝑠 1+𝛼 where Γ(1 + 𝛼) = ∫ 𝑒 -𝑡 𝑡 𝛼 ∞ 0 𝑑𝑡 is the gamma function. Using the inverse Laplace transform and the following properties: 𝐿𝑇(𝑓(𝑡 -𝜏)) = 𝑒 -𝑠𝜏 𝑓 ̃(𝑠), 𝑠 𝑛-1 (𝐿𝑇(𝐴𝑡 𝛼 )) 𝑛 = Γ 𝑛 (1+𝛼) Γ(1+𝑛𝛼) 𝐿𝑇((𝐴𝑡 𝛼 ) 𝑛 ) and (1 + 𝑛𝛼)Γ(1 + 𝑛𝛼) = Γ(2 + 𝑛𝛼), we obtain the bead displacement from Eq. 2:

𝜋 2 ) 2 )

 22 /(𝐴 Γ(1 + 𝛼)) 𝐺 ′′ = (2𝜋) 𝛼 sin (𝛼 𝜋 /(𝐴 Γ(1 + 𝛼)).

Figure 1 Spatially

 1 Figure 1 Spatially-resolved intracellular microrheology. (A) Bright-field image of an RPE1 cell plated on a crossbow micropattern. 2 µm diameter beads were internalized in the cell (black arrow). The dashed and dotted lines show the micropattern and cell contours respectively. Scale bar, 10 µm. (B) (Top) Schematic of the step relaxation experiment. The bead is initially trapped at the center of the optical tweezers. A time t=0 s, the cell is moved in a stepwise fashion to the left on a distance X s =0.5 μm. The bead is displaced on a distance X b <X s due to the viscoelastic nature of its microenvironment. (Bottom) Single and averaged bead displacement plots. Error bars represent standard error mean (N=152 beads from 147 cells). (C) Distribution of the beads locations relative to the micropattern. Position (x=0, y=0) denotes the center of the crossbow micropattern. 𝜃 denotes the angle between the bead position and the polarity axis (𝜃 = 0) of the crossbow micropattern. The dashed blue contour represents the average nucleus shape (see Fig. S1B). (D) The relaxation strongly depends upon the bead location within the cell: bead 1 located in the central region of the micropattern relaxes more slowly than bead 2 located at the periphery. The blue lines are fits to the experimental data using the power law (PL) model. (E) Maps showing the spatial distribution of the storage modulus (G'), the loss modulus (G''), the shear modulus G 0 and the G'/G'' ratio. The maps are color-coded from blue (low values) to red (high values). Data are from N=152 beads from 147 cells. The dashed white contour represents the average nucleus shape (see Fig. S1B). (F) Polar plot of the shear modulus (G 0 ) as a function of the angle 𝜃. Data were pooled and averaged in 30° intervals.

Figure 2

 2 Figure 2 Contribution of the actin and microtubule cytoskeleton to the spatially averaged mechanical properties of the cytoplasm. (A) Relaxation curves for control cells (black) and for cells treated with drugs destabilizing microtubules (nocodazole, taxol) or microfilaments (latrunculin A and cytochalasin D) (red). Note that the cytoD and control curves superimpose. (B) Average values of the shear modulus G 0 , the storage modulus G', and the loss modulus G'' in the presence of the indicated drug. Data shown in (A) and (B) are averages from N= 152, 44, 56, 24 and 34 beads in control cells, nocodazole-, taxol-, LatA-and CytoDtreated cells respectively. (C) Mean square displacement (MSD) analysis of the motion of 2 µm diameter beads diffusing inside the optical trap in control cells (black) and in cells treated with LatA (red) or CytoD (blue). The averaged MSD plots as a function of time (left) were fitted using a caged diffusion model (see Supplementary Materials and Methods) to obtain the diffusion coefficient (right). Data shown are averages from N=54, 27, and 38 beads for control cells, LatA and CytoD treated cells respectively. Error bars represent standard errors. p-values are determined from Student's t-test for unpaired samples with respect to control cells (*** p<0.0001; ** p<0.001; * p<0.01).

Figure 3

 3 Figure 3 Distinct roles of the actin and microtubule cytoskeletons in the spatial distribution of the mechanical properties of the cytoplasm. (A) Maps showing the distribution of the shear modulus (G 0 ) in the presence of nocodazole, taxol, LatA and CytoD. Data are from N=44, 56, 24 and 34 beads in cells treated with nocodazole, taxol, LatA and CytoD respectively. (B) Plots of the shear modulus as a function of the distance to the nucleus in control cells (left) or cells treated with drugs perturbing microtubules (middle) or actin (right). (C) Shear modulus gradient obtained by a linear fit of the plots shown in (B).

Figure 4

 4 Figure 4Contribution of internal membranes to intracellular rigidity. (A) Relaxation curves for control cells (black), for cells treated with Brefeldin A (BFA, red) and after washout of the BFA followed by a recovery period (blue). (B) Average values of the shear modulus G 0 , the storage modulus G', and the loss modulus G'' in RPE1 control cells, in cells treated with Brefeldin A (BFA) or in cells after BFA wash-out and recovery. (C) Maps showing the distribution of the shear modulus G 0 in control cells, in cells treated with BFA and after recovery from BFA treatment. (D) Plots of the shear modulus as a function of the distance to the nucleus in control cells, in cells treated with BFA and after recovery from BFA treatment (left and middle panels) and corresponding shear modulus gradients (right panel). Data were obtained from N=17, 24 and 21 beads in control, BFA and recovery conditions respectively. Error bars represent standard errors. pvalues are determined from Student's t-test for unpaired samples with respect to control cells (* p<0.01; n.s. non significant).

Figure 5

 5 Figure 5 Energy depletion stiffens the cell cytoplasm. (A) Relaxation curves for control cells (black) and for cells depleted in ATP (red). (B) Average values of the shear modulus G 0 , the storage modulus G', and the loss modulus G'' in control cells and in cells depleted with ATP. (C) Maps showing the distribution of the shear modulus G 0 in control cells and in cells depleted in ATP. (D) Plots of the shear modulus as a function of the distance to the nucleus in control cells and in cells depleted in ATP (left panel) and corresponding shear modulus gradients (right panel). Data were obtained from N=33 and 37 beads in control and ATP depletion conditions respectively. Error bars represent standard errors. p-values are determined from Student's t-test for unpaired samples with respect to control cells (** p<0.001; * p<0.01; n.s. non significant).

Figure 6

 6 Figure 6Mapping intracellular mechanics shows that metastatic breast cancer cells are softer than nontumorigenic cells and reveals differences in the distribution of intracellular viscoelasticity. (A) Relaxation curves for non-tumorigenic MCF-10A cells (black) and for metastatic MDA-MB-231 cells (red). (B) Average values of the shear modulus G 0 , the storage modulus G', and the loss modulus G'' in MCF-10A and MDA-MB-231 cells. (C) Maps showing the distribution of the shear modulus G 0 in MCF-10A and

  𝐴 ̅ and 𝐵 ̅ are the averages of A and B. The cross-correlation is such that -1<R 12 <1 with large positive values indicating strong correlation, low values indicating weak correlation and strong negative values indicating inverse correlation. Note that because of the definition of R 12 , no error bar is available for cross-correlation. To avoid edge effects when the maps were of different shapes (see above), the correlation analysis was restricted to a sub-region of interest (about 100 µm 2 ) at the center of the micropattern when required.

2 [

 2 was used to fit the MSD(∆𝑡) curve: 𝑀𝑆𝐷(∆𝑡) = 𝑅 𝑐𝑎𝑔𝑒 4𝐷 𝑐𝑎𝑔𝑒 ∆𝑡 + 𝑐𝑠𝑡, where 𝑅 𝑐𝑎𝑔𝑒 is the cage radius, 𝐷 𝑐𝑎𝑔𝑒 is the diffusion coefficient of the cage and 𝐷 is the diffusion coefficient of the bead. The diffusion coefficient 𝐷 was deduced from the fit. the area under the stage displacement curve (black dashed line). The bead step amplitude X b is defined as the initial displacement of the bead following the step displacement of the stage (X s =0.5 µm in all experiments except in Fig.S2 C-E). (B) Standard linear liquid (SLL) viscoelastic model composed of a dashpot of viscosity  0 in series with a Kelvin-Voigt body (spring of spring constant  and dashpot of viscosity  in parallel). x s , x b and x 1 denote respectively the positions of the stage, the bead and the junction between the spring and the Kelvin-Voigt body with respect to the center of the optical trap. The first ten seconds of the relaxation curves were fitted using the model and the rheological parameters were deduced from the fit (see Methods). (C-E) To check whether the experiments were performed in a linear regime, the step displacement of the stage was varied from 0.4 to 0.6 µm. The corresponding relaxation curves of the bead position were averaged (C) and fitted using the SLL model or the PL model. The model-independent rigidity index (RI) and bead step amplitude (X b ) were measured (D,E top panels) and the spring constant and viscosity (E, middle panels) and the complex shear modulus (E, bottom panel) were deduced from the SLL or PL fits. (N=16, 9 and 23 for control 0.4, 0.5 and 0.6 µm stage step sizes respectively). (F-H) 3D localization of 2 µm diameter beads internalized in RPE1 cells plated on crossbow-shaped micropatterns. Live cells were stained with the membrane marker CellMask and imaged using confocal microscopy. A typical example (F), the averaged x-z and y-z profiles (G, top two images) and the maximum intensity projections of the x-z and y-z profiles (G, bottom two images) are shown (from N=10 cells). The positions of the 16 beads internalized in these cells are marked in red on the maximum intensity projections. The cell maximum height and the bead z-position are quantified in (H) from N=30 beads in 22 cells. Scale bars, 5 µm.

Figure S3 (

 S3 Figure S3(A,B) Analysis of the relaxation of the same two beads as in Fig.1C,D using the viscoelastic SLL model. The blue lines are fits to the experimental data using the SLL model. (C,D) Maps showing the spatial distribution of the model-independent parameters, the rigidity index (RI) and the bead step amplitude X b (C), and of the viscoelastic parameters obtained with the SLL model, the spring constant and the viscosity  (D). Data are from N=166 beads from 147 cells. Average values of the mechanical parameters are SI=0.39±0.01, X b =0.349±0.005 µm, =224±15 pN/µm, =18.5±2.0 Pa.s. (E) Polar plot of the rigidity index (RI) as a function of the angle 𝜃. Data were pooled and averaged in 30° intervals. (F) Crosscorrelation of the rigidity index map with respectively the bead step amplitude (X b ) map, the spring constant ()map and the viscosity () map. Note that the scale is between 0.7 and 1.0 and that the crosscorrelation is high (values close to 1) in the three measurements. (G) Dispersion of the rigidity index (RI) in cells plated on crossbow-shaped micropatterns and in non-patterned cells. The averaged nucleus shape in micropatterned cells is taken from Fig. S1B. Data are from N=24 beads chosen inside the red rectangular region (-2.5 < x < 2.5 µm and 0 < y < 7 µm) in patterned cells and from N=25 beads located

(

  A,B) Average values of the model-independent rigidity index and bead step amplitude (A) and of the spring constant and the viscosity obtained with the SLL model (B) in control cells or in ATP-depleted cells.(C,D) Maps showing the distribution of the storage modulus G' and the loss modulus G'' obtained with the PL model (C), the distribution of the rigidity index and the bead step amplitude (D, top maps) and the distribution of the spring constant and the viscosity obtained with the SLL model (D, bottom maps) in control cells or in ATP-depleted cells. Data were obtained from N=33 and 37 beads in control and ATP depletion conditions respectively. (E) Immunofluorescence images of microtubules and actin in individual fixed cells (left panels) or averaged (right panels) in control cells or in ATP-depleted cells. Averaged images are from N= 10 cells for control cells and for ATP depleted cells respectively. Scale bar, 10 μm.

Figure S9 Metastatic

 S9 Figure S9 Metastatic breast cancer cells are softer than non-tumorigenic cells and exhibit a different distribution of intracellular viscoelasticity. (A,B) Average values of the model-independent rigidity index and bead step amplitude (A) and of the spring constant and the viscosity obtained with the SLL model (B) in nontumorigenic MCF-10A breast cells and in metastatic MDA-MB-231 breast cancer cells. (C,D) Maps showing the distribution of the rigidity index and the bead step amplitude (C, top maps), the distribution of the spring constant and the viscosity obtained with the SLL model (C, bottom maps) and the distribution of the storage modulus G' and the loss modulus G'' obtained with the PL model (D) in MCF-10A and in MDA-MB-231 cells. Data were obtained from N=43 and 54 beads in MCF-10A cells and in MDA-MB-231 cells respectively. (E) Auto-correlation lengths of the rigidity index, the spring constant and the viscosity maps. (F) Relative contributions of elastic storage and viscous loss to the global rigidity of the cytoplasm evaluated by calculating the averaged ratios between G' and G'' obtained by the PL analysis (left) or by cross-correlating the rigidity map with the spring constant map or the viscosity map in each cell type (right). (G) Averaged immunofluorescence images of the microtubule and actin cytoskeletons in MCF-10A and in MDA-MB-231 cells from N=10 cells for MCF-10A and for MDA-MB-231 cells. Scale bar, 10 μm. (H) Average values of the rigidity index in MCF-10A cells treated or not with latrunculin A (LatA) and in MDA-MB-231 cells treated or not with taxol (Tax). Data were normalized by their respective controls and were

Fluorescence intensities were quantified

  Fig. S2

  Ratio between RI µ cross-correlationand RI η cross-correlation -15 -10 -5

  Fig. S6

Acknowledgements

We thank Vera Costa for providing the MCF-10A and MDA-MB-231 cells. We thank David Guet, Mathieu Pinot and Charlotte Alibert for technical assistance with the optical tweezers set-up. We thank Timo Betz, Maxime Dahan, Mathieu Coppey, Paul Janmey, and Sébastien Manneville for fruitful discussions. K.M. was funded by the Plan Cancer 2012-2014 INSERM-Tecsan (grant number PC201125) and by a grant from the Labex CelTisPhyBio. A.A. acknowledges funding from Agence Nationale de la Recherche (ImmunoMeca, ANR-12-BSV5-0007-01; Initiatives d'excellence, Idex ANR-11-IDEX-0005-02; and Labex Who Am I?, ANR-11-LABX-0071).

Supplementary figure legends