Viscoelastic properties of suspended cells measured with shear flow deformation cytometry
Abstract
Numerous cell functions are accompanied by phenotypic changes in viscoelastic properties, and measuring them can help elucidate higher level cellular functions in health and disease. We present a highthroughput, simple and lowcost microfluidic method for quantitatively measuring the elastic (storage) and viscous (loss) modulus of individual cells. Cells are suspended in a highviscosity fluid and are pumped with high pressure through a 5.8 cm long and 200 µm wide microfluidic channel. The fluid shear stress induces large, ear ellipsoidal cell deformations. In addition, the flow profile in the channel causes the cells to rotate in a tanktreading manner. From the cell deformation and tank treading frequency, we extract the frequencydependent viscoelastic cell properties based on a theoretical framework developed by R. Roscoe [1] that describes the deformation of a viscoelastic sphere in a viscous fluid under steady laminar flow. We confirm the accuracy of the method using atomic force microscopycalibrated polyacrylamide beads and cells. Our measurements demonstrate that suspended cells exhibit powerlaw, soft glassy rheological behavior that is cellcycledependent and mediated by the physical interplay between the actin filament and intermediate filament networks.
Editor's evaluation
This paper describes an inexpensive but very powerful microfluidic approach to quantitatively determine the viscoelastic properties of living cells from their deformation in a flow. Its implementation seems simple so that even people not specialized in cell mechanics can use it, and the method offers the possibility to perform measurements on a large number of cells (up to 50100 per second). The data are compelling and this technique should set a new standard in the field.
https://doi.org/10.7554/eLife.78823.sa0eLife digest
Cells in the human body are viscoelastic: they have some of the properties of an elastic solid, like rubber, as well as properties of a viscous fluid, like oil. To carry out mechanical tasks – such as, migrating through tissues to heal a wound or to fight inflammation – cells need the right balance of viscosity and elasticity. Measuring these two properties can therefore help researchers to understand important cell tasks and how they are impacted by disease.
However, quantifying these viscous and elastic properties is tricky, as both depend on the timescale they are measured: when pressed slowly, cells appear soft and liquid, but they turn hard and thick when rapidly pressed. Here, Gerum et al. have developed a new system for measuring the viscosity and elasticity of individual cells that is fast, simple, and inexpensive.
In this new method, cells are suspended in a specialized solution with a consistency similar to machine oil which is then pushed with high pressure through channels less than half a millimeter wide. The resulting flow of fluid shears the cells, causing them to elongate and rotate, which is captured using a fast camera that takes 500 images per second. Gerum et al. then used artificial intelligence to extract each cell’s shape and rotation speed from these images, and calculated their viscosity and elasticity based on existing theories of how viscoelastic objects behave in fluids.
Gerum et al. also investigated how the elasticity and viscosity of cells changed with higher rotation frequencies, which corresponds to shorter timescales. This revealed that while higher frequencies made the cells appear more viscous and elastic, the ratio between these two properties remained the same. This means that researchers can compare results obtained from different experimental techniques, even if the measurements were carried out at completely different frequencies or timescales.
The method developed by Gerum et al. provides a fast an inexpensive way for analyzing the viscosity and elasticity of cells. It could also be a useful tool for screening the effects of drugs, or as a diagnostic tool to detect diseases that affect the mechanical properties of cells.
Introduction
Eukariotic cells can carry out complex mechanical tasks such as cell division, adhesion, migration, invasion, and force generation. These mechanical activities in turn are essential for higher order cell functions including differentiation, morphogenesis, wound healing, or inflammatory responses. Since cell mechanical activities are accompanied by phenotypic changes in the cell’s viscoelastic properties, measuring them can help elucidate higher order cell functions in health and disease (Urbanska et al., 2020). For example, the activation of neutrophils in response to proinflammatory agents is typically accompanied by a marked increase in cell stiffness (Frank, 1990; Fabry et al., 2001), which can subsequently lead to a sequestration of the stiffened cells in small capillaries for example of the lung (Doerschuk et al., 1993). This process may be relevant for the progression and exacerbation of inflammatory diseases such as coronavirus disease 2019.
In this report, we describe a quantitative, lowcost, highthroughput, and simple method to measure the viscoelastic properties of cells, specifically the storage modulus ${G}^{\prime}$, and the loss modulus ${G}^{\prime \prime}$. The cells are suspended in a highviscosity (0.5–10 Pa s) fluid (e.g. a 2% alginate solution) and are pumped at pressures of typically between 50 and 300 kPa through a several centimeter long microfluidic channel with a square cross section (200x200 µm in our setup). The fluid shear stress induces large cell deformations that are imaged using a complementary metaloxidesemiconductor (CMOS) camera at frame rates of up to 500 frames/s to achieve a measurement throughput of up to 100 cells/s. Images are stored and analyzed offline at a speed of around 50 frames/s on a standard desktop PC equipped with a graphics card.
The method takes advantage of two physical principles: First, the shear stress profile inside a long microfluidic channel depends only on the pressure gradient along the channel, which can be precisely controlled, and the channel geometry, which is fixed. Importantly, the shear stress profile does not depend on the viscosity of the cell suspension medium and smoothly increases from zero at the channel center to a maximum value at the channel walls. Accordingly, cells appear circular near the channel center and become increasingly elongated near the channel walls. As the width of the channel is significantly larger than the cell diameter, fluid shear stresses remain approximately constant across the cell surface, which considerably simplifies the fluid dynamics computations compared to existing microfluidic methods. From the stressstrain relationship, we estimate the storage modulus of the cell, which characterizes its elastic behavior.
Second, depending on the flow speed profile inside the channel, the cells rotate in a tanktreading manner, similar to a ball that is compressed between two countermoving parallel plates. Shearflow induced tanktreading was first theoretically explored by Einstein, 1906, and was later experimentally observed by H. SchmidSchönbein et al. in sheared red blood cell suspensions (SchmidSchöenbein and Wells, 1969; Fischer et al., 1978). Tanktreading arises as the flow speed of the suspension fluid in contact with the cell surface facing the channel center is larger compared to the flow speed at the opposite side. Hence, the rotational speed of this tanktreading motion increases with increasing shear rate near the channel walls. Tanktreading in combination with the cell’s viscous properties leads to energy dissipation, which limits the increase of cell strain at higher stresses near the channel walls. From this behavior, we extract the loss modulus of the cell, which characterizes its viscous behavior. Since the microfluidic channel is several centimeters long, most cells, with the exception of those in the center of the channel, have already undergone several full rotations before reaching the field of view. Therefore, the cell deformations are in a steady state, which is another major difference compared to existing microfluidic approaches and greatly simplifies the calculation of viscoelastic cell parameters.
For the calculation of viscoelastic cell parameter, we use a theoretical framework developed by Roscoe, 1967 that describes the deformation of a viscoelastic sphere in a viscous fluid under steady shear flow. This theory allows us to compute the stiffness (shear modulus) and viscosity of a cell from 5 measurable parameters. First, the fluid shear stress acting on the cell must be known, which we compute based on the extension of Poiseuille’s equation to channels with square crosssection (Delplace, 2018). Second, we measure the cell deformation (cell strain) from brightfield microscopy images. Third, we measure the alignment angle of the deformed cell with respect to the flow direction. This alignment angle depends on the ratio between cell viscosity and the viscosity of the suspension fluid. Fourth, we compute the local viscosity of the suspension fluid based on measurements of the radial flow speed profile in the channel, which we obtain from multiple images of the same cell during its passage through the channel. Fifth, since cell stiffness and cell viscosity are frequencydependent, we measure the tanktreading frequency of each cell.
The Roscoe model assumes that cells behave as a KelvinVoigt body consisting of an elastic spring in parallel with a resistive (or viscous) dashpot element. This then gives rise to a complex shear modulus with storage modulus ${G}^{\prime}$ and loss modulus ${G}^{\prime \prime}$, measured at twice the tank treading frequency (because a given volume element inside the cell is compressed and elongated twice during a full rotation). Roscoe theory, however, makes no assumptions about how ${G}^{\prime}$ and ${G}^{\prime \prime}$ might change as a function of frequency. A commonly used simplified assumption is that the elastic and viscous elements of the KelvinVoigt body are constant (Fregin et al., 2019). Hence, ${G}^{\prime}$ plotted versus frequency would be flat, and ${G}^{\prime \prime}$ would increase proportional with frequency. An alternative and, as we will show in this report, a more accurate model, known as the structural damping formalism, predicts that both ${G}^{\prime}$ and ${G}^{\prime \prime}$ increase with frequency according to a powerlaw (Fabry et al., 2001). In either case, to compare the stiffness and viscosity of cells that have experienced different tanktreading frequencies, it is important to scale the stiffness and fluidity of each cell to a reference frequency, for example of 1 Hz.
Using cell lines and calibrated polyacrylamide beads, we verify that our method provides accurate quantitative measurements of viscoelastic properties. Measurement results are not or only marginally influenced by experimental details such as the viscosity of the suspension fluid or the time point after suspending the cells. We demonstrate that the cell’s viscoelastic properties measured with our method conform to soft glassy powerlaw rheology that has been reported for a wide range of cells measured with different methods. We also show that our method can be used for doseresponse measurements of drugs that induce actin cytoskeleton disassembly, and that these responses are modulated by the cell cycle and the intermediate filament network of the cells.
Results
Measurement setup
We image the cells in brightfield mode while they are moving through the microchannel (Figure 1a–c). Using a neural network, we detect cells that are in focus at the midplane of the microchannel (Figure 1b), and segment their shapes (Figure 1d). We then quantify the cell position and cell shape by fitting an ellipse to the segmented cell image, from which we obtain the centroid coordinate (x_{0}, y_{0}), the length of the semimajor axis $a$ and the semiminor axis $b$, and the angular orientation β of the major axis with respect to the $x$(flow) direction (Figure 1e). From $a$ and $b$, we compute the cell strain $\u03f5$ using Equation 10 (Figure 2a). We also compute the local fluid shear stress $\sigma ({y}_{0})$ for a cellfree fluid at the cell’s centroid position using Equation 4 (Figure 1f).
Cell deformations under fluid shear stress
Cells are nearly circular in the center, and they elongate and align in flow direction near the channel walls (Figure 1c, Figure 2a, b) where they are exposed to higher fluid shear stress (Figure 1f). Cells imaged at the same position within the channel also tend to become more elongated with increasing pressure (Figure 1c). When we plot cell strain $\u03f5$ versus shear stress σ across the microfluidic channel (Figure 2c), we find that the cell strain increases nonlinearly with increasing fluid shear stress. In particular, the slope of the strain versus stress relationship decreases for higher stress values. This behavior is predominantly due to a dissipative process caused by the tank treadlike motion of the cells.
Tanktreading
The radial velocity gradient of the flow field (the shear rate $\dot{\gamma}$) creates a torque on the sheared and elongated cells and causes them to align in flow direction (Figures 1e and 2b) and to rotate in a tanktreading manner (Video 1): the cell’s elongated shape and alignment angle β remain stationary, but internally, the cell is constantly rotating as if being kneaded between two plates (SchmidSchöenbein and Wells, 1969; Fischer et al., 1978).
From a series of images that show the same cells as they flow through the channel, we compute the radial velocity profile $v(y)$ of the fluid flow (Equation 9, Figure 2d), the shear rate profile $\dot{\gamma}(y)$ (Equation 7, Figure 2e), and the tanktreading frequency f_{tt} of each cell (Figure 2h and i). We find that the tanktreading frequency of a cell is zero at the channel center and increases towards the channel walls (Figure 2j and k). At low shear rates (low driving pressure or near the channel center), the rotation rate ${\omega}_{\mathrm{tt}}/\dot{\gamma}$ of individual cells is close to the Einsteinlimit of 1/2, as theoretically predicted for spheres that are tanktreading in a Newtonian fluid (Einstein, 1906; Snijkers et al., 2011; Roscoe, 1967). Tanktreading dissipates energy in proportion to the cell’s internal viscosity, rotation frequency, and strain. This energy dissipation therefore limits the cell strain in regions of high shear rate and hence shear stress (Figure 2c).
Viscoelastic model
We can quantitatively explain the nonlinear strainstress relationship (Figure 2c) and its pressuredependency by a theoretical framework describing the deformation and alignment of viscoelastic spheres in a viscous fluid under steady shear flow (Roscoe, 1967). This theoretical framework (in the following referred to as Roscoetheory) predicts that the cell strain $\u03f5$ increases proportional with the shear stress σ and the sine of the alignment angle β, and inversely proportional with the elastic modulus ${G}^{\prime}$ of the cell (Equation 16). The alignment angle β in turn depends on the cell’s loss modulus ${G}^{\prime \prime}$, the local shear rate $\dot{\gamma}$ and the local sheardependent viscosity η of the suspension fluid (Equation 17). With increasing elastic modulus, cells are predicted to deform less (smaller strain $\u03f5$) and to align less in flow direction (larger alignment angle β) when exposed to a fixed shear stress and shear rate. With increasing loss modulus, cells are also predicted to deform less but to align more in flow direction. Thus, from the measurements of cell strain, alignment angle, local shear stress, local shear rate, and local viscosity, Roscoetheory allows us to compute the viscoelastic properties (${G}^{\prime}(\omega )$ and ${G}^{\prime \prime}(\omega )$) of individual cells at twice their specific angular tanktreading frequency, $\omega =2\cdot 2\pi {f}_{\mathrm{tt}}$.
Powerlaw behavior of cells
When we plot ${G}^{\prime}$ and $G}^{\mathrm{\prime}\mathrm{\prime}$ of individual cells versus twice their tanktreading frequency f_{tt} (Figure 3a), we find that the complex shear modulus $\stackrel{~}{G}={G}^{\mathrm{\prime}}+i{G}^{\mathrm{\prime}\mathrm{\prime}}$ of a cell population approximately follow a powerlaw relationship of the form
where $\mathrm{\Gamma}$ is the Gammafunction, $k$ is the elastic shear modulus (cell stiffness) referenced to an arbitrarily chosen frequency of 1 Hz by setting ${\omega}_{0}=2\pi \mathrm{rad}/\mathrm{s}$, α is the powerlaw exponent that characterizes the fluidity of the cell (zero indicating purely Hookean elastic behavior, unity indicating Newtonian viscous behavior), and $i=\sqrt{1}$ (Fabry et al., 2001). Such a behavior of a cell population emerges if the rheology of individual cells also follows a powerlaw relationship. Thus, using Equation 1, we can compare the mechanical behavior of cells measured at different tank treading frequencies by computing their stiffness $k$ (using Equation 21) and fluidity α (using Equation 22).
We find in agreement with previous reports (Desprat et al., 2005; Balland et al., 2006; Cai et al., 2013; Hecht et al., 2015; Bonakdar et al., 2016) that the individual stiffness values $k$ are typically lognormal distributed, and the fluidity values α are normal distributed (Figure 3b). Moreover, also in agreement with previous reports, we find an inverse relationship between stiffness and fluidity, whereby stiffer cells tend to be less fluidlike (Fabry et al., 2001; Smith et al., 2005; Lange et al., 2015). Due to this coupling, the mode of the twodimensional distribution of α and $k$ (the most common combination of α and $k$ among all cells, as estimated from the maximum of the Gaussian kerneldensity, Figure 3b), provides a robust measure for the mechanical behavior of a cell population.
Stress stiffening
To test if suspended cells exhibit stress stiffening, as previously reported (Lange et al., 2017), we increase the driving pressure from 100 kPa to 300 kPa, which increases the maximum shear stress at the channel wall from 116 Pa to 349 Pa (Figure 1f). Cell fluidity remains constant over this pressure range, but the median stiffness of the cell population increases with increasing pressure by 33% (Figure 3c and d). To explore to which extent this stiffness increase is caused by a higher shear stress as opposed to a higher shear rate, we keep the pressure constant at 200 kPa but increase the alginate concentration from 1.5% to 2.5% and therefore the viscosity of the suspension medium from 2.2 Pa⋅s to 9.2 Pa⋅s (zeroshear viscosity ${\eta}_{0}$ as determined with Equation 6). This causes the shear rate to decrease and leads to a slight but not statistically significant increase in stiffness and fluidity (Figure 3e and f). Hence, the increase of cell stiffness at a higher driving pressure is induced by stressstiffening and not by a higher shear rate. We also verify that cell stiffness and fluidity remain stable over a period of up to 60 min after suspending the cells in a 2% alginate solution (Figure 3g and h).
Validation with polyacrylamide beads
To evaluate the accuracy of our method, we measure 16 µm diameter polyacrylamide (PAAm) beads with three different nominal stiffnesses, in a range similar to living cells (Figure 4a–c). The frequencydependency of ${G}^{\prime}$ and ${G}^{\prime \prime}$ of the beads are calibrated using oscillatory atomic force microscopy (AFM), and conform to a powerlaw relationship with an additional Newtonian viscosity μ according to
with ${\omega}_{0}=2\pi $ rad/s (Figure 4b). Using shear flow deformation cytometry, we also find a powerlaw behavior (Figure 4c and d). As the maximum frequency remains below 10 Hz in these measurements, however, effect of the Newtonian viscosity term μ is less pronounced (Figure 4c), and we therefore perform a global fit of Equation 2 to the data using a constant μ for all conditions. The values of $k$ and α for beads with different acralymidebisacrylamid concentrations are comparable between AFM and shear flow deformation cytometry measurements (Figure 4—figure supplement 1). Moreover, $k$ and α are largely pressureindependent (from 0.2 to 2 bar; Figure 4e), as expected for a linear material such as PAAm. Fluidity is close to zero for strains below unity (α=0.092 for 5.9% ${C}_{\mathrm{AAmBis}}$, and α=0.074 for 3.9% ${C}_{\mathrm{AAmBis}}$), indicating predominantly elastic behavior as expected. Fluidity increases slightly at higher strains (Figure 4f), likely due to fluidinduced (poroelastic) relaxation processes (Kalcioglu et al., 2012). Together, these results demonstrate that our method provides quantitatively accurate estimates for the elastic and dissipative properties of soft spherical particles.
We next compare the viscoelastic properties of monocytic THP1 cells probed by shear flow cytometry and atomic force microscopy (AFM). We acquire forceindentation curves at rates of ~1 /s (Figure 5c), which is within the range of strain rates that cells experience in our shear flow cytometry setup. AFM measurements show that THP1 cells conform to powerlaw rheology with an additional Newtonian viscosity term according to Equation 2, from which we extract the shear modulus $k$ and fluidity α (Figure 5b). THP1 cells appear stiffer (at 1 Hz) and more fluidlike when measured with shear flow cytometry ($k$=82 Pa, α=0.44) compared to AFM ($k$=52 Pa, α=0.25). Despite these differences, AFM measurements confirm the applicability of powerlaw rheology, and they also show a lognormal distribution of cell stiffness $k$ as well as an inverse relationship between $k$ and fluidity α as seen in our shear flow measurements (Figure 5a and b).
Doseresponse measurements
We perform doseresponse measurements using latrunculin B (LatB), which prevents the polymerization of monomeric actin and leads to a depolymerization of the actin cytoskeleton (Urbanska et al., 2020). NIH3T3 fibroblasts soften with increasing doses of LatB (1–1000 nM) according to a sigmoidal (HillLangmuir) relationship, with a maximum response of 1.47fold and a halfmaximum dose of EC50=35.2 nM (Figure 6a). These responses agree with published data obtained using realtime deformability cytometry (RTDC) measurements on HL60 cells (maximum response 1.46fold, EC50=26.5 nM) (Urbanska et al., 2020). When we measure promyoblast HL60 suspension cells with our setup, EC50 is similar to published data (26.4 nM), but the maximum response is much higher (5.0fold) (Figure 6b).
Role of intermediate filaments
To explore the attenuated LatB responsiveness of NIH3T3 fibroblasts compared to HL60 leukemia cells, we reasoned that NIH3T3 cells express high levels of the intermediate filament protein vimentin (Figure 7a) that may protect the cells from excessive deformations when filamentous actin is depolymerized. To test this idea, we measure the stiffness of NIH3T3 and vimentinknockout (vim(/)) fibroblasts in response to 30 min treatment with cytochalasin D (2 µM), which binds to the barbed end of filamentous actin and—similar to LatB—leads to a net depolymerization of the actin cytoskeleton (Figure 7a). The NIH3T3 cell line has been established from mouse embryonic fibroblasts (MEFs) by spontaneous immortalization (Todaro and Green, 1963). We followed the corresponding protocol for MEFs obtained from vimentinknockout mouse embryos (ColucciGuyon et al., 1994). Thus, the three cell lines investigated here are of the same cell type. We find that cytochalasin D treated vim(/) cells soften by a considerably greater extent (2.16fold) compared to wildtype cells (1.22 fold) (Figure 7b and c), in support of the notion that vimentin stabilizes the cytoskeleton.
To explore if the cytoskeletonstabilizing effect of vimentin is a general feature also of other intermediate filament networks, we measure the cytochalasin D response of desmintransfected vimentin knockout MEFs (vim(/)des). Desmin, which is the dominant intermediate filament in skeletal muscle, forms an intermediate filament network in fibroblasts that is structurally similar to the vimentin network in wildtype cells (Figure 7a). Similar to vimentinexpressing MEFs, vim(/) desminexpressing MEFs also display an attenuated cytochalasin D response (1.37fold), confirming that both the vimentin and desmin intermediate filament network can protect cells from excessive deformations when filamentous actin is depolymerized (Figure 7c).
Cell cycle dependence
In our measurements, we observe that larger NIH3T3 cells tend to be softer compared to smaller cells (Figure 8a). We hypothesized that this weak sizedependence of cell stiffness might be attributable to cell cycle progression, which leads to changes in chromatin compaction and cell volume. To test this hypothesis, we extend our setup to acquire green fluorescent images alongside bright field images of cells transfected with a twocolor fluorescent Fucci cell cycle indicator (SakaueSawano et al., 2008). Fuccitransfected cells display high red and low green fluorescence when they are in G1 phase, and low red but increasing levels of green fluorescence as they progress into S, G2, and early Mphase (SakaueSawano et al., 2008). We measure the cell cycle distribution of NIHT3T cells before harvesting using epifluorescence microscopy (Figure 8b), and map the distribution to the green fluorescent intensities measured in our shear flow cytometry setup (Figure 8c).
We find as expected that cell radius increases with cell cycle progression (Figure 8d). In addition, cell stiffness steadily increases towards the end of the G1 and the beginning of the Sphase, and then rapidly decreases as the cell cycle progresses (Figure 8e). When we bin the cells according to their green fluorescent intensities (i.e. according to their cell cycle progression) and plot stiffness versus cell radius (Figure 8f), we find substantially larger and nonmonotonic fluctuations of cell stiffness versus cell radius, compared to the smaller, monotonic decrease of cell stiffness in the radiusbinned data (Figure 8c). These differences arise because changes in cell stiffness and cell radius occur at different stages of the cell cycle.
Discussion
Viscoelastic cell properties can be measured with established methods such as atomic force microscopy (Cordes et al., 2020), micropipette aspirations (Zhelev et al., 1994), or magnetic tweezer microrheology (Bonakdar et al., 2016). These methods have a relatively low throughput of typically below 10–100 cells/hr. The need to measure cell mechanical properties with substantially higher throughput led to the recent development of various microfluidic techniques (Urbanska et al., 2020) including hydrodynamic stretching (Gossett et al., 2012), realtime deformability cytometry (Otto et al., 2015; Fregin et al., 2019), microfiltration (Rowat et al., 2012), and microconstriction systems (Lange et al., 2015; Lange et al., 2017).
Our method builds on previously established highthroughput microfluidic approaches, with several modifications: We suspend cells in a medium that is pumped with high pressure (typically 50–300 kPa) through a long, parallel microfluidic channel with one inlet and outlet (no flowfocussing geometry is needed). Such simple microfluidic channels are commercially available at low cost, which we expect will facilitate a widespread adoption of the technique. The large driving pressure gives rise to sufficiently large (>50 Pa) shear stresses to induce measurable cell deformations. The high pressure can be controlled with a simple pressure regulator, without the need for a precise microfluidic controller—another advantage compared to existing methods that typically operate under lower pressure. The width and height of the channel (200 µm) are much larger than the cell diameter, which prevents clogging due to debris that is often encountered in microfluidic constrictionbased systems (Lange et al., 2015; Lange et al., 2017). Most importantly, the large channel diameter ensures that fluid shear stresses do not vary appreciably across the cell, which simplifies the analysis of cell mechanical properties as the cells do not deform into complex bullet or hourglasslike shapes as seen in other methods (Lange et al., 2015; Lange et al., 2017; Urbanska et al., 2020). By suspending the cells in a fluid with high viscosity (typically >1 Pa⋅s), we achieve a flow speed that is sufficiently low (<20 mm/s) so that the cells’ speed, position, and shape can be captured without motion blur at a typical exposure time of 30 µs using a standard CMOScamera mounted to a routine laboratory microscope.
The lateral flow profile in the channel causes a tanktreadinglike cell motion, which imposes periodic cell deformations with frequencies on the order of 10 Hz. At such low frequencies and strain rates, elastic cell properties dominate over viscous cell properties (Fabry et al., 2001; Mietke et al., 2015). The cell transit through the microfluidic channel lasts for several seconds, which is much longer than the period time of the cells’ tank treading rotation, implying that the measured cell deformations can safely be assumed to have reached a steadystate. Measuring viscoelastic parameters from steadystate cell deformations has a major advantage over existing microfluidic techniques in that no viscoelastic models (e.g. Maxwell, Voigt etc.) or preconceived response functions (e.g. exponential, biexponential, powerlaw etc.) are needed to fit the transient cell deformation responses during the cells’ passage through narrow constrictions and channels (Lange et al., 2015; Fregin et al., 2019). If for example an exponential function is fitted to a powerlaw creep response, the resulting viscoelastic cell properties would largely depend on the time scale of the experiment (e.g. the passage time of the cell through the microfluidic constriction or channel) and less so on the cell’s intrinsic mechanical properties (Fabry et al., 2001; Lange et al., 2015; Fregin et al., 2019) By contrast, the values obtained with our method are not influenced by the time scale of the experiment.
From images of the same cell as it is flowing through the channel, we estimate the tank treading frequency and the flow velocity; from the flow velocity profile across the channel, we compute the local shear rate (Equation 5) and the local sheardependent viscosity of the suspension fluid (Equation 6); from the radial cell position, we compute the local shear stress (Equation 4); from the cell shape, we compute the strain (Equation 10) and the alignment angle in flow direction. From these measurements, we finally compute the cell’s viscoelastic properties (stiffness and fluidity, Equation 21 and 22). Hence, once the flow velocity profile is known, we can determine the viscoelastic properties from a single image because all cell deformations are in a steadystate.
We provide userfriendly software for image acquisition and data analysis on a standard PC, which can be downloaded at https://github.com/fabrylab/shear_flow_deformation_cytometer (copy archived at swh:1:rev:2d35a697243c432cddd52e10d2e3c5121f922adf; Gerum, 2022). Currently, the method stores the acquired uncompressed images on a hard drive, which in the case of typically 10,000 images for a single experiment lasting 20 s amounts to a storage space of nearly 4 GB. The image data are analyzed afterwards, which at a rate of around 50 images per second can take several minutes. Future software developments and faster computer hardware will enable image analysis on the fly for realtime shear flow deformation cytometry.
The computation of viscoelastic properties is based on a theoretical model proposed by R. Roscoe that describes the deformation of homogeneous, isotropic, incompressible neoHookean viscoelastic spherical particles under fluid shear stress (Roscoe, 1967). Cells in suspensions, however, are known to deform nonlinearly (Lange et al., 2017), with stress or strain stiffening that is more pronounced than the stiffening predicted for a neoHookean material. Therefore, our measurements represent an effective secant modulus and not a smallstrain tangential modulus. Moreover, cells do not consist of a homogeneous material but of different components (e.g. the cell cortex and the nucleus) with different mechanical properties (Zhelev et al., 1994; Rowat et al., 2012; Mietke et al., 2015; Cordes et al., 2020). As a consequence, cells do not always deform into ellipsoidal shapes but occasionally deform into sigmoidal shapes, which becomes more pronounced in response to larger shear stresses or drugs that soften the cytoskeleton, such as cytochalasin D or latrunculin B.
Despite the simplified assumptions of the Roscoe theory, however, our cell rheological measurements agree with previously published findings that were obtained using a range of different methods and models, namely that suspended cells show a behavior that is consistent with powerlaw rheology, that the elasticity of individual cells is lognormal distributed, that the fluidity of individual cells is normaldistributed, and that stiffness and fluidity scale inversely (Fabry et al., 2001; Alcaraz et al., 2003; Desprat et al., 2005; Lange et al., 2015). These experimental findings are in agreement with predictions from soft glassy rheology (Sollich, 1998; Fabry et al., 2001). Moreover, we show that stiffness and fluidity values of polyacrylamide beads and cells measured with shear flow deformation cytometry agree quantitatively with AFM measurements.
Our measurements are insensitive to changes in the viscosity of the suspension medium, demonstrating that the fluidmechanical assumptions of the Roscoe theory hold in the case of living cells in a shearthinning suspension fluid. We find that cells appear stiffer when measured at higher driving pressures, likely due to stress or strainstiffening of the cells (Lange et al., 2017). When we measure linearly elastic polyacrylamide beads over a 10fold pressure range (from 20 to 200 kPa), we see a constant, pressureindependent shear modulus and agreement with the stiffness and fluidity values measured using AFM, demonstrating that the Roscoe theory gives quantitatively accurate estimates, regardless of driving pressure and suspension fluid viscosity.
Roscoe theory estimates the cell viscosity relative to the viscosity of the suspension fluid, which for a shear thinning fluid such as alginate can be difficult to measure. However, since we know the fluid profile in the microfluidic channel (from the flow speed of hundreds of cells), we can estimate the rheological properties of the suspension fluid, including its shear thinning behavior. This ability is more than a byproduct of our method and could be valuable for example for researchers interested in bioink development for applications in biofabrication. Moreover, we measure the complex rheology of the suspension fluid at the very same time and under the same conditions (temperature, range of shear rates) as the cells. Therefore, unlike other microfluidic cytometry methods (Otto et al., 2015; Fregin et al., 2019), our quantification of viscoelastic cell properties does not rely on separate measurements using coneplate or other types of rheometers. The rheological parameters of alginate solutions measured with our method closely agree with coneplate rheometer measurements, with relative deviations of 31% over a shear rate spanning 5 orders in magnitude (from 0.01 to 1000 s^{1}).
Our method measures each cell at a single tanktreading frequency that depends on the cell’s lateral position in the channel. Thus, with our method we sample the frequencydependent mechanical properties of a cell population simply by observing cells at different channel positions. By contrast, with existing methods, time or frequencydependent cell responses can only be explored by choosing different strain rates, for example via adjusting the driving pressure (Lange et al., 2015). The tanktreading frequency can be directly measured using particle flow analysis methods in a subset of the cells that shows small features with high contrast (Fischer et al., 1978). For the remaining cells, it is possible to estimate the tank treading frequency from the local shear rate according to an empirical equation (Equation 20). This equation holds for the cell types and suspension fluids used in our study, but we do not claim that it holds universally for other cell types or suspension fluids. For example, red blood cells exhibit a tumbling instead of a tanktreading motion at low shear rates (SchmidSchöenbein and Wells, 1969), and Equation 20 underestimates their tanktreading frequency at shear rates beyond $100{s}^{1}$(Fischer et al., 1978).
To demonstrate its practical applicability, we apply our method to measure the stiffness of HL60 cells in response to different doses of the actindepolymerizing agent latrunculin B. We find in agreement with previous observations a halfmaximum dose (EC50) of around 30 nM, but a considerably larger softening of the cells by a factor of 5.4fold at the highest dose of 1 µM, compared to a softening of only 1.5 fold that is seen with other microfluidic techniques constriction microfluidic constrictionbased deformability cytometry (cDC), and realtime deformability cytometry (RTDC) (Urbanska et al., 2020). This higher responsiveness is likely attributable to the relatively low cellular strain rates in our method, which are on the order of $10{s}^{1}$, compared to strain rates of around $100{s}^{1}$ in the case of RTDC. At these high strain rates, viscous cell behavior starts to dominate over cytoskeletonassociated elastic behavior (Fabry et al., 2001; Mietke et al., 2015). Accordingly, when cells are measured with extensional flow deformability, a method that operates at even higher strain rates in the kHzrange, they do not appreciably soften in response to LatB (Urbanska et al., 2020; Gossett et al., 2012).
We also demonstrate that the cell softening induced by cytochalasin D, another actindepolymerizing drug, is attenuated in the presence of intermediate filaments (vimentin or desmin), and becomes more pronounced when intermediate filaments are absent. This finding is in line with earlier reports that intermediate filaments protect cells against excessive strain (Patteson et al., 2020), and that the absence of vimentin in fibroblasts leaves the cells vulnerable to mechanical stress (Eckes et al., 1998). A physical interaction between vimentin intermediate filaments and Factin bundles as mediated by plectin has been demonstrated by high resolution immunoelectron microscopic methods (Svitkina et al., 1996). The importance of vimentinactin interactions has furthermore been corroborated by high resolution structured illumination microscopy in combination with cryoelectron tomography, revealing the intimate association and functional crosstalk between vimentin intermediate filaments and Factin stress fibers (Wu et al., 2022). Here, we have directly demonstrated that the stable introduction of a cytoplasmic intermediate filament protein into intermediate filamentfree cells restores their cytoskeletal functionality and mechanical stability.
Shear stress deformability cytometry can be combined with fluorescent imaging. Here, we image the viscoelastic properties of NIH3T3 cells together with the cell cycle using the fluorescent Fucci indicator. Our data demonstrate that NH3T3 cells stiffen during the course of cell cycle progression in G1 phase, with a maximum stiffness during late G1 – early Sphase, and then soften before they enter the G2 and MPhase. Since cell volume also increases during the transition from G1 to S phase, we find a slight overall dependence of cell stiffness on cell size in the case of NIH3T3 cells (Figure 8c). This cell size dependence is also detectable in HL60 and THP1cells (Figure 3—figure supplement 1).
In summary, shear flow deformation cytometry provides accurate quantitative measurements of elastic and dissipative cell properties at high throughput. The method can be easily and inexpensively implemented on standard or research grade microscopes. Unlike other highthroughput microfluidic methods, the cells are measured under near steadystate conditions at low to moderate strain rates where elastic responses dominate over viscous responses.
Methods
The measurement setup is depicted in Figure 1 a. Video 2 explains the measurement procedure. Cells are suspended in a highviscosity medium (e.g. a 2% alginate solution), and are pressed via a 10 cm long, 1mm inner diameter silicone tube through a 5.8 cm long microfluidic channel with a square cross section of 200x200 µm (CS10000090; Darwin Microfluidics, Paris, France). The driving air pressure of typically 1–3 bar is regulated with a pressure regulator (KPRG114/10, Knocks FluidTechnik, Selm, Germany) and can be switched on or off with a threeway valve (VHK204F04F; SMC, Egelsbach, Germany). The air pressure is measured with a digital pressure gauge (Digi04 0.4%, Empeo, Germany). Cells flowing through the channel are imaged in brightfield mode at 50–500 Hz (depending on the flow speed) with a CMOS camera (acA720520um, Basler, Germany) using a 40x0.4 NA objective (Leica) in combination with a 0.5 x video coupler attached to an inverted microscope. After passing the microchannel, the cells are collected in a waste reservoir.
Cell culture
Cells are cultured at 37 °C, 5% CO_{2} and 95% humidity and are split every 2–3 days for up to 20 passages.
Preparing cells for rheological measurements
Our method for measuring viscoelastic cell properties requires that the cells, if they are adherent to a cell culture dish (NIH3T3, vim(/), vim(/)des), are brought into suspension. For cells grown in 75 cm^{2} flasks, we remove the medium and wash the cells three times with 10 ml of 37 °C PBS. After removing the PBS, 5 ml of 0.05% trypsin/EDTA in PBS are added and distributed over the cells, and after 10 s, 4 ml of the supernatant are removed. Cells are then incubated for 3–5 min at 37 °C, 5% CO_{2}. 5 ml of 37 °C cell culture medium (Table 1) are added to the flask, and the cells are counted. If cells are already in suspension (THP1 and HL60 cells), the above steps are omitted. 10^{6} cells are taken out of the flask, centrifuged for 5 min at 25 rcf (NIH3T3, vim(/) and vim(/)des) or 290 rcf (HL60 and THP1) to remove the supernatant, gently mixed in 1 ml of equilibrated suspension fluid (see below), transferred to a 2 ml screwcup test tube, and centrifuged at 150 rcf for 30 s to remove air bubbles.
Suspension fluid preparation
Alginate solution is prepared freshly for the next day. Sodium alginate powder (Vivapharm alginate PH176, batch nr. 4503283839, JRS Pharma GmbH, Rosenberg, Germany, or alginic acid sodium salt from brown algae, A0682, Sigma Aldrich, for THP1 cells) is dispersed at a concentration of 1.5%, 2%, or 2.5% (w/v) in serumfree cell culture medium (Table 1). The alginate solution is mixed overnight with a magnetic stirrer at room temperature until all powder has been dissolved. The suspension fluid is then equilibrated by incubating for 6 hr at 37 °C, 5% CO_{2}. When prepared with RPMI media (but not when prepared with DMEM nor Sigma Aldrich alginate), the alginate solution is filtered with a 0.45 µm filter before use. 1 ml of alginate solution are then added to the cell pellet of 10^{6} cells in the Falcon tube and mixed using a positive displacement pipette (15314274, Gilson/Fisher Scientific) by slowly (~2 s cycle time) and repeatedly (10 x) sucking the liquid in and out. The alginatecell suspension is then transferred into a 2 ml screwcup test tube and centrifuged for 30 s at 150 rcf to remove air bubbles.
Drug treatment
Drugs are mixed in the alginate for at least 15 min with a magnetic stirrer at 350 rpm inside an incubator (37 °C, 5% CO_{2}, 95% relative humidity) prior to mixingin the cells. Cells are prepared as described above and mixed with the alginatedrug mixture using a positive displacement pipette by slowly (~2 s cycle time) and repeatedly (10 x) sucking the liquid in and out. The alginatedrugcell suspension is transferred into a 2 ml screwcup test tube and incubated for a prescribed time at 37 °C, 95% rH. Prior to measurements, the alginatedrugcell suspension is centrifuged at 150 rcf for 30 s to remove air bubbles.
Inhibition of actin polymerization on NIH3T3, vimentinknockout and desminknockin MEFs is performed with cytochalasin D (Cat. No. C8273; SigmaAldrich, St. Louis, MO). Cytochalasin D is dissolved in DMSO at a stock concentration of 20 mM. The equilibrated alginate (3 ml) is either mixed with cytochalasin D to a final concentration of 2 µM, or mixed with DMSO to a final concentration of 0.01% (DMSO control), or mixed with 3 µl of DMEM (negative control). Cells harvested from a single cell culture flask are split into three groups of 10^{6} cells, each group is suspended in one of the alginate solutions as described above, stored in an incubator for 15 min (alternating between either negative control of DMSO control), 30 min (drugtreated), and 45 min (alternating between either DMSO control or negative control), and measured.
Inhibition of actin polymerization on NIH3T3 cells is performed with latrunculin B (LatB, Cat. No. L5288; SigmaAldrich, St. Louis, MO, dissolved in DMSO at a stock concentration of 2 mM). We add 2 µl of LatB (stock) or 2 µl of DMSO to 4 ml of alginate (final concentration 1000 nM LatB, 0.2% DMSO), and mix with a magnetic stirrer at 350 rpm for 15 minutes. 1850 µl of the alginatedrug mixture is then added to 4 ml of alginate, mixed for 15 min, and the process is repeated to obtain a dilution series with LatB concentrations of 1000, 316, 100, 32, 20, 3.2, and 1 nM. The alginateDMSO mixture is diluted in the same way. Cells are prepared and mixed into the alginate as described above and stored at room temperature for 10 min (LatB) or 20 min (DMSO control) prior to measurements.
Image acquisition
Typically, 10,000 images per measurement are recorded with a CMOS camera (acA720520um, Basler, Germany) at a frame rate of 50–500Hz with an exposure time of 30 µs. To measure the flow speed, each cell has to be recorded in at least 2 consecutive images. Therefore, the frame rate $fr$ is chosen depending on the maximum flow speed v_{max} and the width of the region of interest (ROIx): $fr>{v}_{max}$ / (0.5 ROIx). In our setup, the ROIx is 248µm, resulting in a maximum flow speed of 41mm/s for a frame rate of 500Hz. To prevent motion blur, however, we keep the maximum flow speed to about 20mm/s.
Fluorescent images can be acquired in parallel with the bright field images. A 300 mW diodepumped solidstate laser (wavelength 473nm, VAIN473; Viasho, Beijing, China) serves as an epifluorescent light source, and a beam splitter projects the bright field and fluorescent images onto two synchronized cameras. To separate the light paths, the brightfield illumination is longpass filtered (>590nm), and a bandpass filter (500–550nm) is placed in front of the camera for the fluorescent channel.
We provide software for image acquisition (see below under Software flow chart), which includes a liveviewer and userfriendly interface for entering meta information (e.g. applied pressure, suspension medium, drug treatments) and configuration settings (e.g. frame rate, total number of images to be stored). The software is based on the pypylon library to record the images, and Python (Van Rossum and Drake, 2009) and Qt to provide the user interface.
Cell shape analysis
We normalize the brightfield images by subtracting the mean and dividing by the standard deviation of the pixel intensities. A neural network (UNet Ronneberger et al., 2015, tensorflow Abadi, 2016) trained on labeled images of different cell types and suspension media detects the cell outline and generates a binary mask, to which an ellipse is fitted ($x$,$y$ position of the ellipse center, its semimajor ($a$) and semiminor axis ($b$), and the angle of orientation β of the major axis with respect to the flow (x) direction, see Figure 1d and e, van der Walt et al., 2014). Binary masks that do not conform to an elliptical shape based on circumference or solidity criteria (e.g. due to cell doublets or erroneous cell outlines due to poor image contrast) are discarded.
Finding the channel mid plane and center line
Prior to recording the images, the microscope must be precisely focused to the mid plane (z=0, see Figure 1b) of the channel. To do so, we apply a small pressure (50–100 Pa) to the suspended cells and focus the microscope in phase contrast mode to the bottom of the microchannel, which can be unambiguously identified by stationary or very slowly flowing small debris. We then move the objective up by 75 µm, which corresponds to half the microchannel’s height (100 µm) divided by the refractive index of the suspension medium. We confirmed that the reproducibility of the method is within ±1.7 µm (rms) when a 40x0.6 NA objective is used.
The channel center line ($y=0$, see Figure 1b) is identified from the flow speed profile as a function of the radial ($y$) position. Flow speed is computed by tracking cells over subsequent images and dividing the distance they have moved in $x$direction by the time difference between images. A polynomial of the form.
is then fitted to the velocity profile to identify the center position of the channel (y_{c}), with the maximum flow speed v_{max} at the channel center as the second fit parameter, and the exponent ζ as the third fit parameter. $W$ is the channel width. The fit parameter y_{c} is then used to shift the image ycoordinate origin to the channel center. This procedure ensures that the channel does not need to be precisely centered in the camera’s field of view during the measurements. However, the channel should be aligned as precisely as possible with the field of view. To ensure alignment, we recommend to rotate the camera, as opposed to the slide that holds the channels.
Shear stress profile inside a channel with a square crosssection
The fluid shear stress σ in the mid plane of a channel (blue shading in Figure 1b) with length $L$ and square cross section of height $H$ and width $W$ only depends on the radial position $y$ and the total applied pressure $\mathrm{\Delta}P$ according to an infiniteseries expression (Delplace, 2018).
For all practical purposes, it is sufficient to compute the infinite series for the first 100 terms.
Equation 4 assumes laminar uniaxial parallel flow and neglects entrance and exit effects, which is justified for a long and narrow channel as used in this study ($L$=5.8 cm, $W$ = $H$ = 200 µm). Note that for a given channel geometry and pressure gradient $\mathrm{\Delta}P/L$, the shear stress profile $\sigma (y)$ does not depend on the viscosity of the fluid. Equation 4 remains approximately valid also for nonNewtonian e.g. shearthinning fluids. Equation 4 predicts that the shear stress is zero in the center of the channel and monotonically increases towards the channel wall (Figure 1f).
We take the shear stress $\sigma (y)$ at the cell center $y$ as the average stress acting on the cell. For cells that overstep the channel center, however, the nonmonotonic stress profile implies that the average stress can be larger than the stress at the cell center. Therefore, and because cells near the channel center deform and align only marginally, which makes the computation of mechanical properties errorprone, we exclude all cells from further analysis that are closer than one cell radius to the channel center.
Velocity profile, shear rate profile, and viscosity
The fit function (Equation 3) only approximates the true velocity profile, which is sufficient to efficiently and robustly find the channel center. For subsequent computations that require higher precision, we determine the velocity profile by integrating the shear rate. We compute the shear rate $\dot{\gamma}(y)$ as the shear stress σ (Equation 4) divided by the viscosity η.
For shear thinning fluids such as alginate solutions, the viscosity η is not constant but depends on the shear rate $\dot{\gamma}$. We describe the shear thinning behaviour of the viscosity by the Cross model (Cross, 1965)
with zeroshear viscosity ${\eta}_{0}$, relaxation time τ and powerlaw shear shearthinning exponent δ (Figure 2—figure supplement 2).
When Equation 6 is inserted into Equation 5, we obtain.
This equation can be written as.
and numerically solved for $\dot{\gamma}(y)$ by root finding using the NewtonRaphson method.
Finally, to obtain the velocity profile $v(y)$, we integrate the numerically obtained shear rate $\dot{\gamma}(y)$ over the channel, using 5 point Gaussian quadrature
with the boundary condition ${v}_{y=W/2}=0$. The viscosity parameters (${\eta}_{0}$, τ, δ) that best match the velocity profile are determined as follows. We choose five Gaussian quadrature points ${y}^{\prime}$ between (0,$W/2$) and numerically compute $\dot{\gamma}$ at the quadrature point ${y}^{\prime}$ using Equation 8. To ensure convergence, we start iterating with a value of $\dot{\gamma}$ that yields the maximum of the righthand side of Equation 8 plus a small number $\u03f5$. The weighted sum of $\dot{\gamma}$ at the Gaussian quadrature points ${y}^{\prime}$ is then the velocity at the radial position $y$. This procedure is repeated for different values of (${\eta}_{0}$, τ, δ) until a minimum of the squared differences between the measured and fitted velocity profile is found.
We find that the rheological parameters (${\eta}_{0}$, τ, δ) of the suspension medium obtained this way closely agree with coneplate rheology measurements (Müller et al., 2007). Moreover, the velocity profile for different pressure values can be accurately predicted (Figure 2—figure supplement 1), demonstrating that Equation 6 accurately describes the shear thinning behavior of the suspension fluid.
Computing the shear strain from the cell shape
Suspended cells under zero shear stress have an approximately circular shape with radius r_{0}. When exposed to constant shear stress, the cell deforms to an elliptical shape with semimajor axis $\stackrel{~}{a}=a/{r}_{0}$ and semiminor axes $\stackrel{~}{b}=b/{r}_{0}$ (in $x,y$direction) and $\stackrel{~}{c}=c/{r}_{0}$ (in zdirection), normalized to the radius r_{0} of the undeformed cell, so that $1=\stackrel{~}{a}\cdot \stackrel{~}{b}\cdot \stackrel{~}{c}$. Assuming the cell consists of an incompressible material and the stress inside the deformed cell is uniform, the strain $\u03f5$ can be computed from $\stackrel{~}{a}$, $\stackrel{~}{b}$ and $\stackrel{~}{c}$ using (Equation 10; Roscoe, 1967).
(corresponding to the righthand side of Equation 79 in Roscoe, 1967 without the sign error).
This requires solving a set of shape integrals that depend on the semimajor axis a and semiminor axis b.
With the integration variable λ. ${\mathrm{\Delta}}^{\prime}$ is defined as.
(Equation 11 corresponds to Equation 39 in Roscoe, 1967, and Equation 12 corresponds to Equation 18 in Roscoe, 1967.)
The shape integral $I$ is precomputed for different ratios of $\stackrel{~}{a}$ and $\stackrel{~}{b}$ and then taken from a lookup table.
Computing the cells’ storage and loss modulus
We calculate ${G}^{\prime}$ from σ, β, $a$, $b$ according to Roscoe, 1967.
(corresponding to the lefthand side of Equation 79 in Roscoe, 1967).
We calculate ${G}^{\prime \prime}$ from β, $\stackrel{~}{a}$, $\stackrel{~}{b}$, η, ω according to Roscoe, 1967.
(corresponding to Equation 80 in Roscoe, 1967).
with
(corresponding to Eqation 43 in Roscoe, 1967).
(corresponding to Equation 21 in Roscoe, 1967).
A given volume element inside the cell is compressed and elongated twice during a full rotation. Hence, the frequency ω at which ${G}^{\prime}$ and ${G}^{\prime \prime}$ is obtained using Equation 16 and Equation 17 is twice the angular tanktreading frequency $2{\omega}_{\mathrm{tt}}$.
Tank treading
We measure the tanktreading frequency as follows. We observe each cell as it travels through the fieldofview and cutout small image frames with the cell at its center (Figure 2h). We then track the movement of characteristic small features using optical flow estimated by the TVL1 algorithm (Zach et al., 2007; van der Walt et al., 2014), and calculate their speed and distance during their rotation around the cell’s center. The speed versus the ellipsecorrected radius is fitted with a linear relationship to determine the average angular speed (Figure 2i). The slope of this relationship is taken as the rotation frequency of the cell.
In cases where the tanktreading frequency cannot be measured (e.g. due to poor contrast or absence of cellinternal features that can be tracked), we estimate the tanktreading frequency following the approach outlined in Snijkers et al., 2011. Data shown in Figure 2j and k demonstrate that the measured rotation rate ${\omega}_{\mathrm{tt}}/\dot{\gamma}$ (angular frequency divided by the local shear rate) collapses onto a master relationship when plotted against the shear rate. The angular tanktreading frequency ${\omega}_{\mathrm{tt}}=2\pi {f}_{\mathrm{tt}}$ of the cells can then be predicted with an empirical relationship according to
when $\dot{\gamma}$ is given in units of 1 /s (Snijkers et al., 2011).
Scaling the rheology
Cells show powerlaw rheology according to Equation 1, which implies that the cell stiffness $k$ and the powerlaw exponent α (cell fluidity) fully describe the cell rheological properties. Cell stiffness $k$ and cell fluidity α can be obtained from ${G}^{\prime}$ and ${G}^{\prime \prime}$ by rearranging Equation 1 as follows
with $\omega =2{\omega}_{\mathrm{tt}}$ and ${\omega}_{0}=2\pi $ rad/s. We use a Gaussian kernel density estimation (Silverman, 1965; Virtanen et al., 2020) to compute the mode of the 2D distribution for stiffness $k$ and fluidity α, which corresponds to the “most representative” cell with the highest joint probability for stiffness $k$ and fluidity α.
Software flowchart
In the following, we summarize the sequence of steps and procedures for measuring cell mechanical properties with our method (Figure 1—figure supplement 1).
First, typically 10,000 image frames of cells flowing through the channel are recorded with an image acquisition program (recording.py, https://github.com/fabrylab/shear_flow_deformation_cytometer). Second, the images are analyzed offline with an evaluation pipeline (evaluate.py, https://github.com/fabrylab/shear_flow_deformation_cytometer). The pipeline loads the images and finds and segments cells at the focal plane using a neural network (Ronneberger et al., 2015). From the segmented cell shape, morphological properties ($x$,$y$ position, half major and minor axes $a$ and $b$, orientation β, solidity, circumference) are extracted using the regionprops method of the skimage library (van der Walt et al., 2014). Poorly or erroneously segmented cells that deviate from an elliptical shape are filtered out based on circumference and solidity criteria. From a measurement with 10,000 image frames, typically 5000–10,000 cells are identified for subsequent analysis.
Next, the program identifies cells that are detected across multiple subsequent frames, based on shape and position, computes the flow speed, and applies an particle image velocimetry algorithm to extract the tank treading frequency f_{tt}. Equation 3 is then fitted to the speed versus yposition relationship of all cells, yielding the channel center y_{c} and the maximum flow speed v_{max}.
The shear stress acting at the center position of each cell is computed using Equation 4.
The shear rate at the center position of each cell is computed using a set of equations as described above (Equations 5–9). This procedure also yields the parameters that describe the viscosity and shearthinning rheology of the suspension fluid (Equation 6).
The cell strain is computed from the half major and minor axis $a$ and $b$ using Equation 10. Subsequently, ${G}^{\prime}$ and ${G}^{\prime \prime}$ of each cell at twice its angular tank treading frequency is computed using Equation 16 and Equation 17.
To compare the mechanical properties of cells that have experienced different tanktreading frequencies, we scale ${G}^{\prime}$ and ${G}^{\prime \prime}$ to a frequency of 1 Hz using Equation 22 and Equation 21, yielding the stiffness $k$ and fluidity α of individual cells. The average stiffness $k$ and fluidity α of the cell population is determined from the maximum of the twodimensional Gaussian kernel density computed using the scipy.stats.gaussian_kde method of the scipy library (Silverman, 1965; Virtanen et al., 2020).
PAAm reference bead preparation
Polyacrylamide hydrogel microparticles (PAAm beads) are produced using a flowfocusing PDMSbased microfluidic chip described in Girardo et al., 2018. Briefly, a stream of a polyacrylamide pregel mixture is squeezed by two counterflowing streams of an oil solution to form droplets with a mean diameter in the range of 11.5–12.5 µm. The oil solution is prepared by dissolving ammonium Krytox surfactant (1.5% w/w), N,N,N’,N’tetramethylethylenediamine (0.4% v/v), and acrylic acid Nhydroxysuccinimide ester (0.1% w/v) in hydrofluoroether HFE 7500 (Ionic Liquid Technology, Germany). The pregel mixture is obtained by dissolving and mixing acrylamide (40% w/w), bisacrylamide (2% w/w) and ammonium persulfate (0.05% w/v) (all from Merck, Germany) in 10 mM Trisbuffer (pH 7.48). Particles with three different elasticities are obtained by diluting the pregel mixture in Trisbuffer to final acrylamidebisacrylamide concentrations of 3.9%, 5.9%, 6.9% respectively. Alexa Fluor 488 Hydrazide (ThermoFisher Scientific, Germany) is dissolved in D.I. water (stock solution 3 mg/ml) and added to the mixture for a final concentration of 55 µg/ml to make the particles fluorescent. Droplet gelation is carried out at 65 °C for 12 hr. The droplets are washed and resuspended in 1 x PBS.
Atomic force microscopy (AFM) of cells and PAAm beads
AFMbased microrheology measurements for PAAm beads are performed using a Nanowizard 4 (JPK BioAFM, Bruker Nano GmbH, Berlin). The measurements are carried out using a wedged cantilever with a flat surface parallel to the measurement dish. The cantilever is prepared by applying a UV curing glue to a tipless cantilever (PNPTRTL, nominal spring constant $k$ = 0.08 N/m used for the stiff (6.9% ${C}_{\mathrm{AAmBis}}$) beads, or Nanoworld or ArrowTL1, nominal spring constant $k$ = 0.03 N/m used for the medium (5.9% ${C}_{\mathrm{AAmBis}}$) and soft (3.9% ${C}_{\mathrm{AAmBis}}$) beads) as described in Stewart et al., 2013. Prior to each experiment, the optical lever sensitivity is measured from the forcedistance relationship of a polystyrene bead attached to a glass surface, and the cantilever spring constant is measured using the thermal noise method (Hutter and Bechhoefer, 1993). Measured spring constants are 0.09 N/m for PNPTRTL cantilevers, and 0.018 N/m for ArrowTL1cantilevers.
To perform the AFM microrheology measurements, the cantilever is lowered with a speed of 10 µm/s until a force between 1–3 nN is reached, corresponding to an indentation depth ${\delta}_{0}$ between 1.5–3 µm. The cantilever is then sinusoidally oscillated with an amplitude of 30 nm for a period of 10 cycles. This procedure is repeated for different oscillation frequencies in the range between 0.1–150 Hz. To extract the complex shear modulus ${G}^{*}$ of the PAAm beads, the forceindentation curves are analyzed as described in Alcaraz et al., 2003 using the Hertz model that describes the deformation of a soft sphere between two flat surfaces in the limit of small deformations. The complex shear modulus is then computed according to
where ν is the Poisson ratio of the PAAm bead (assumed to be 0.5), ω is the angular frequency of the oscillations, $F(\omega )$ and $d(\omega )$ are the Fourier transforms of the force and indentation signal, $R$ is the radius of the PAAm bead, ${\delta}_{0}$ is the initial indentation, and $b(0)$ is the hydrodynamic drag coefficient of the cantilever with the surrounding liquid. The hydrodynamic drag coefficient is measured as described in Alcaraz et al., 2002 and estimated to be $b(0)$ = 5.28 Ns/m for PNPTRTL cantilevers and $b(0)$ = 29.7 Ns/m for Arrow TL1 cantilevers.
AFMbased measurements for THP1 cells are performed with foursided regular pyramidaltipped MLCTbioDC(D) cantilevers (Bruker). The spring constant of the cantilever is measured from the thermal noise spectrum in air, and the optical lever sensitivity is measured from the thermal noise spectrum in liquid (Sumbul et al., 2020). The cells are immobilized to plastic petri dishes coated with polyLlysine at a concentration of 0.01 mg/mL for 10 min. Force curves are measured at 3 or more positions around the cell center for a constant indentation speed of 5 µm/s up to a maximum force of 0.8 Nn. At each position, at least 3 forcedistance curves are obtained. We determine the viscoelastic stepresponse stress relaxation function $E(t)$ of the cell by leastsquare fitting the theoretical force response to the measured force curve during indentation with a pyramidal tip (Efremov et al., 2017)
where $F$ is the force acting on the cantilever tip, δ is the indentation depth, $t$ is the time since initial contact, t_{m} is the duration of approach phase, t_{ind} is the duration of complete indentation cycle, and t_{1} is the auxiliary function determined by the equation
The viscoelastic step response function $E(t)$ is assumed to follow the relationship
where the reference time t_{0} is set to 1 s so that $k$ is the cell’s shear modulus measured at time $t$ = 0.159 s (corresponding to ω=1 rad/s as in the flow deformability measurements). The cell’s Poisson ratio ν is assumed to be 0.5, and α is the cell’s fluidity.
Rheology of alginate solutions
We measure the viscosity of the alginate solution at a temperature of 25 °C at shear rates between 0.01 s^{1} and 1000 s^{1} using a temperaturecontrolled rheometer (DHR3, TAInstruments, USA) with stainless steel cone and plate (diameter of 40 mm with a cone angle of 2° and a 65 µm truncation gap). Temperature is controlled with a Peltierelement. Equilibration time and measurement time are set to 30 seconds for every measurement point (logarithmic sweep, 5 points per decade). Every sample is rested for three minutes inside the rheometer to ensure temperature equilibration. A solvent trap with deionized water is used to prevent drying of the alginate samples.
Cell cycle measurement with Fucci
We use NIH3T3 cells that display the fluorescent ubiquitinationbased cell cycle indicator (FastFUCCI) reporter system after lentiviral transduction. The lentivirus is generated by transfection of LentiX 293T cells (Takara, #632180) with pBOBEF1FastFUCCIPuro (Addgene, #86849), a packaging plasmid psPAX2 (Addgene, #12260), and an envelope plasmid pCMVVSVG (Addgene, #8454), using Lipofectamine 2000 reagent (Invitrogen, #11668–019). 48 hr after transfection, infectious lentiviruscontaining supernatant is harvested, centrifuged (500 x g, 10 min), and 10fold concentrated using the LentiXconcentrator reagent (Takara, #631232). NIH3T3 cells are seeded 24 hr prior to transduction at a density of 10 000 per cm^{2}. Three days after transduction, cells are cultured for at least 5 additional days in medium containing puromycin (5 µg/ml) to select successfully transduced cells.
In our shear flow deformation cytometry setup, we measure only the green fluorescence signal, indicating cells in S, G2 and early Mphase (SakaueSawano et al., 2008), and deduce that cells with a green fluorescence intensity below a certain threshold are in G1 phase. To set this threshold, we measure both the red fluorescence signal (indicating cells in G1 phase SakaueSawano et al., 2008) and the green fluorescence signal of individual cells prior to harvesting, using an epifluorescence microscope. We then compute the greenfluorescence intensity threshold, normalized to the median intensity that best separates the cells in G1 phase from the cells in S, G2 and early Mphase. Because some cells fluoresce green and red at the same time, 22.6% of cells in G1 phase and 2.4% of the cells in S, G2 and early Mphase are erroneously classified when the classification is based on the green fluorescence signal alone. After harvesting and suspending the cells in alginate, they are measured in the shear flow setup. Brightfield images are analyzed as described above to segment cells that are in focus, and the fluorescence intensities are averaged over the segmented cell area.
Data availability
Software is made available at GitHub, https://github.com/fabrylab/shear_flow_deformation_cytometer, (copy archived at swh:1:rev:2d35a697243c432cddd52e10d2e3c5121f922adf). CSV files containing the data of all individual cells used for the study have been made available on Dryad (https://doi.org/10.5061/dryad.5hqbzkh8p).

Dryad Digital RepositoryViscoelastic properties of suspended cells measured with shear flow deformation cytometry.https://doi.org/10.5061/dryad.5hqbzkh8p
References

ConferenceTensorflow: A system for largescale machine learningIn 12th USENIX Symposium on Operating Systems Design and Implementation (OSDI. pp. 265–283.

Microrheology of human lung epithelial cells measured by atomic force microscopyBiophysical Journal 84:2071–2079.https://doi.org/10.1016/S00063495(03)750140

Power laws in microrheology experiments on living cells: comparative analysis and modelingPhysical Review. E, Statistical, Nonlinear, and Soft Matter Physics 74:021911.https://doi.org/10.1103/PhysRevE.74.021911

Quantifying celltocell variation in powerlaw rheologyBiophysical Journal 105:1093–1102.https://doi.org/10.1016/j.bpj.2013.07.035

Rheology of nonnewtonian fluids: a new flow equation for pseudoplastic systemsJournal of Colloid Science 20:417–437.https://doi.org/10.1016/00958522(65)90022X

Creep function of a single living cellBiophysical Journal 88:2224–2233.https://doi.org/10.1529/biophysj.104.050278

Comparison of neutrophil and capillary diameters and their relation to neutrophil sequestration in the lungJournal of Applied Physiology 74:3040–3045.https://doi.org/10.1152/jappl.1993.74.6.3040

Impaired mechanical stability, migration and contractile capacity in vimentindeficient fibroblastsJournal of Cell Science 111 (Pt 13):1897–1907.https://doi.org/10.1242/jcs.111.13.1897

Eine neue bestimmung der moleküldimensionenAnnalen Der Physik 324:289–306.https://doi.org/10.1002/andp.19063240204

Scaling the microrheology of living cellsPhysical Review Letters 87:148102.https://doi.org/10.1103/PhysRevLett.87.148102

SoftwareShear flow deformation cytometer, version swh:1:rev:2d35a697243c432cddd52e10d2e3c5121f922adfSoftware Heritage.

Standardized microgel beads as elastic cell mechanical probesJournal of Materials Chemistry. B 6:6245–6261.https://doi.org/10.1039/c8tb01421c

Calibration of atomic‐force microscope tipsThe Review of Scientific Instruments 64:1868–1873.https://doi.org/10.1063/1.1143970

Unbiased highprecision cell mechanical measurements with microconstrictionsBiophysical Journal 112:1472–1480.https://doi.org/10.1016/j.bpj.2017.02.018

Extracting cell stiffness from realtime deformability cytometry: theory and experimentBiophysical Journal 109:2023–2036.https://doi.org/10.1016/j.bpj.2015.09.006

ConferenceUnet: Convolutional networks for biomedical image segmentationIn International Conference on Medical image computing and computerassisted intervention. pp. 234–241.https://doi.org/10.1007/9783319245744

On the rheology of a suspension of viscoelastic spheres in a viscous liquidJournal of Fluid Mechanics 28:273.https://doi.org/10.1017/S002211206700206X

Nuclear envelope composition determines the ability of neutrophiltype cells to passage through micronscale constrictionsJournal of Biological Chemistry 288:8610–8618.https://doi.org/10.1074/jbc.M112.441535

Effect of viscoelasticity on the rotation of a sphere in shear flowJournal of NonNewtonian Fluid Mechanics 166:363–372.https://doi.org/10.1016/j.jnnfm.2011.01.004

Rheological constitutive equation for a model of soft glassy materialsPhysical Review E 58:738–759.https://doi.org/10.1103/PhysRevE.58.738

Onestep calibration of afm in liquidFrontiers in Physics 8:301.https://doi.org/10.3389/fphy.2020.00301

Plectin sidearms mediate interaction of intermediate filaments with microtubules and other components of the cytoskeletonThe Journal of Cell Biology 135:991–1007.https://doi.org/10.1083/jcb.135.4.991

Quantitative studies of the growth of mouse embryo cells in culture and their development into established linesThe Journal of Cell Biology 17:299–313.https://doi.org/10.1083/jcb.17.2.299

BookA duality based approach for realtime tvl1 optical flowIn: Hamprecht FA, Schnörr C, Jähne B, editors. Pattern Recognition. Berlin Heidelberg: Springer. pp. 214–233.https://doi.org/10.1007/9783540749363

Role of the membrane cortex in neutrophil deformation in small pipetsBiophysical Journal 67:696–705.https://doi.org/10.1016/S00063495(94)805296
Decision letter

Alphee MichelotReviewing Editor; Institut de Biologie du Développement, France

Anna AkhmanovaSenior Editor; Utrecht University, Netherlands

Clément CampilloReviewer; Université d'Evry, France

Timo BetzReviewer; Universtity of Göttingen, Germany
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 "Viscoelastic properties of suspended cells measured with shear flow deformation cytometry" for consideration by eLife. Your article has been reviewed by 2 peer reviewers, and the evaluation has been overseen by a Reviewing Editor and Anna Akhmanova as the Senior Editor. The following individuals involved in review of your submission have agreed to reveal their identity: Clément Campillo (Reviewer #1); Timo Betz (Reviewer #2).
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:
You will see from the comments of the two Reviewers that they agree on the usefulness of your method. However, they also point out some weaknesses in this study that need to be addressed in the revision. Essential revisions should include:
1) A much more indepth introduction and discussion of existing methods (see Reviewer 1's comment). The interest of a method paper is certainly to describe new protocols/tools, but also to discuss their interest (and also their disadvantages) compared to existing tools. This is particularly important as papers using microfluidics to characterize the mechanical properties of cells have recently been published (including the paper by Oliver Otto which is mentioned). This discussion would be useful for less specialized readers (e.g. cell biologists wanting to characterize their cells but not necessarily having a strong biophysics background). Why would it make a difference in their experiments to have a frequencydependent response? In what specific cases would having frequency dependent values allow them to discriminate cells better than static parameters?
These modifications to the manuscript should be particularly easy to make.
2) Further analysis of the experimental data to justify or refute the validity of a single powerlaw approach versus a two powerlaw approach, or a modification of the experimental setup to acquire and present data only in the linear regime (see reviewer comment 2). This would answer two questions:
(a) whether the discrepancy is truly due to strain stiffening, and
b) whether one can reliably use the data in the high frequency level to obtain the correct stiffness and powerlaw exponent. This point will certainly give you much more work, but it is essential to address it to convince us that the final values and model are reliable.
3) Presenting experiments investigating the role of intermediate filaments from the same cell line so that the results are easier to interpret.
Reviewer #1 (Recommendations for the authors):
I am not able to evaluate the validity of the Roscoe model and the equations 10 19, nevertheless the agreement between the measurements obtained by the two techniques used are very convincing. This is even surprising in the case of THP1 cells (Figure 5a and b), because SFDC and AFM probe cell mechanics at very different scales.
The novelty of the technique compared to Fregin NatComm 2019 has to be discussed in detail as the main claim of the article is the approach for cell mechanics measurements. From a physics point of view, the technique presented here gives for instance G' and G' as a function of frequency, whereas Fregin et al., give only effective cell elasticity and viscosity. Therefore, the measurement of cell's mechanics is more complete with this assay, as described in the manuscript l.217. Note that the experimental setup presented in the manuscript is simpler than the one presented by Fregin et al., which might favor its use by other labs. This justifies the publication of the article in eLife, as physical methods are in the scope of the journal. On the other hand, such a detailed physical description might not be required to only discriminate between cell populations, as it is the main application of this type of experiment.
The introduction is very short and the authors start presenting their technique immediately. I think the readers need more context on why these type of measurements are useful, which experimental techniques are used and the mechanical models on which these techniques rely. Some of these aspects are in the Discussion section, to highlight the interest of their technique but the introduction should be more detailed.
On the tank treading movement of cells in shear flows, several articles are cited but the article comports no discussion at all on this mechanism, and how it has been characterized for red blood cells for instance.
Similarly, some readers may not be familiar with viscoelastic models, the article is not very pedagogic on this.
Reviewer #2 (Recommendations for the authors):
While I really like the method, and think it should be published in eLife after successful revision, I am worried about the reliability of some parts of the paper. Mainly, I think the strain stiffening explanation needs to be better nailed. Why is it not possible to change the conditions, so that the deformation remains in the linear regime throughout the measurement. Even if correct, all the data that is in the strain stiffening regime would then lead to wrong stiffness and powerlaw exponent. Would it not make more sense, to only focus on the results as close to the measurements as possible? We have a G' and G' for a given frequency, why do we need to translate this into other parameters that are even less reliable?
Besides this main concern, I have a couple of other points that I would like to transmit to the authors to consider for a further improvement of the paper.
– Why is in figure 2a/b no data given for the position close to 0? I did not realize an explanation for this when I read the paper. Is it problematic in this regime to determine the deformation?
– In figure 2 g it is shown that for a large enough shear rate the rheometer and the flow cytometer values of viscosity are the same. If I understand this correctly, it also means that only right at the center the shear rate is so low that this matters… however in the distances that correspond to this shear no measurements are given. Is this right? If so, maybe it is worthwhile to mention this in more detail.
– Why seem the stiffness values of the PAA beads to be so pressure dependent as seen in figure 4e. And in figure S3, the values of the AFM and the flow cytometer should be plotted on the same y scale, or ideally even plotted over each other. What keeps you from adding the values a in b,c, and the same for d and g?
– In Equation 2 µ seems not to be defined.
– For the cell experiments with vimentin it is very confusing why you compare WT 3T3 with MEF desminknockouts and knockins of vimentin. It looks like you don't have WT MEF. Why comparing different cell types. This makes it hard to believe that your conclusion is correct. At least all differences between 3T3 and MEFs might be because of the different cell type. Use also MEFs for the 'normal' situation.
– Also, there is a statistical test in figure 7 b and c missing. How many cells did you measure for each of the datapoints in 7b,c?
– In the discussion you mention a local viscosity of the fluid. This is a bit misleading as it suggests that the viscosity depends on the position, but it depends on the shear rate (which depends on the position for a given pressure).
– It should be more focused on the actual measurement of G' G' at a single frequency. I think we don't know sufficiently well if a single powerlaw can explain the measurements, to use this hypothesis in a way that makes the reader think the method provides information about frequency dependence from a single snapshot.
– It would be good to get a reference for the statement that you find similar values for the E50 as other authors (line 269).
– Please tone down the statement of line 284286. There is no statistical test, and you compare apples with oranges… don't say that you have shown the effect of vimentin. All you see is the increase in stiffness in the knockin situation.
– It would be interesting to provide a measure of variance in the text describing the different \σ values. (line 502503)
– Why not using a nonshear thinning medium to do the measurements? Everything should be much more simple there.
– You treat \eta_0, \tau, and \δ as independent parameters. Is this really the case, or could some depend on each other? Also please give the values you measured for these parameters.
– You say that stress stiffening happens, but then you use a model (equations 21,22) that breaks down under these conditions. Here you basically say that the measurements you obtain in the high frequency regime do not work. Why are these not excluded, and how do you determine the stress at which we cannot take the G' and G' as unaffected.
– In Line 580 you say the AFM goes from 0.1 to 150Hz, but the measurements you show are up to 10000 rad/s. What is right?
https://doi.org/10.7554/eLife.78823.sa1Author response
Essential revisions:
You will see from the comments of the two Reviewers that they agree on the usefulness of your method. However, they also point out some weaknesses in this study that need to be addressed in the revision. Essential revisions should include:
1) A much more indepth introduction and discussion of existing methods (see Reviewer 1's comment). The interest of a method paper is certainly to describe new protocols/tools, but also to discuss their interest (and also their disadvantages) compared to existing tools. This is particularly important as papers using microfluidics to characterize the mechanical properties of cells have recently been published (including the paper by Oliver Otto which is mentioned). This discussion would be useful for less specialized readers (e.g. cell biologists wanting to characterize their cells but not necessarily having a strong biophysics background). Why would it make a difference in their experiments to have a frequencydependent response? In what specific cases would having frequency dependent values allow them to discriminate cells better than static parameters?
These modifications to the manuscript should be particularly easy to make.
We have followed these recommendations and have now explained in the Introduction why it is important to measure the frequencydependent properties (namely to compare measurements across different platforms and time/frequency scales).
2) Further analysis of the experimental data to justify or refute the validity of a single powerlaw approach versus a two powerlaw approach, or a modification of the experimental setup to acquire and present data only in the linear regime (see reviewer comment 2). This would answer two questions:
(a) whether the discrepancy is truly due to strain stiffening, and
b) whether one can reliably use the data in the high frequency level to obtain the correct stiffness and powerlaw exponent. This point will certainly give you much more work, but it is essential to address it to convince us that the final values and model are reliable.
As we explain in more detail in our response to the reviewers, we erroneously assumed that the probing frequency equals the tank treading frequency. However, the probing frequency is twice the tanktreading frequency. With this correction, a single powerlaw now describes our data much better. We have reanalyzed all data with this correction, and moreover we have changed the reference frequency (for the cell stiffness k) from 1 rad/s now to 1 Hz, which we think is more intuitive for most readers.
Further, we actually have used a twopower law description of cell mechanics already in our original submission (Equation 2), and we use a single power law description only where the frequency or shear rate of the measurements is so low that the second power law does not yet come into play, which is the case for all our shear flow deformation cytometry measurements on cells.
Our approach of describing the frequency response with a power law (Equations 1 and 2) has been validated by numerous different techniques, in vastly different cells, and by many different laboratories. It has been shown to be valid also for strainstiffening cells, it correctly models the frequency response of cells over a very large range of frequencies (spanning at least 5 orders in magnitude), and it allows for the comparison of data obtained at different frequencies or shear rates (which most other microfluidic techniques do not). We have provided more detailed explanations below.
3) Presenting experiments investigating the role of intermediate filaments from the same cell line so that the results are easier to interpret.
We actually have used the same cell type, a fibroblast line derived from mouse embryo fibroblasts (MEF) by spontaneous immortalization. It is important to note that the NIH3T3 cell line is derived from MEFs just as our Vim (/) cell line. We had stated this previously only in Methods, now explicitly state this also in the main text in Results.
Reviewer #1 (Recommendations for the authors):
I am not able to evaluate the validity of the Roscoe model and the equations 10 19, nevertheless the agreement between the measurements obtained by the two techniques used are very convincing. This is even surprising in the case of THP1 cells (Figure 5a and b), because SFDC and AFM probe cell mechanics at very different scales.
As we explain below in more detail in our answer to Reviewer #2, we had assumed that the probing frequency is the tanktreading frequency. However, the cells are actually probed at twice the tanktreading frequency. After correcting this, we find that measurements with our shear flow deformation cytometer show a higher fluidity compared to AFM measurements and also give a lower stiffness. With other words, the two techniques do not (any longer) give nearly identical absolute values (although they are reasonably similar), but we still find common cell behavior with both techniques, namely (1) powerlaw rheology, (2) a lognormal distribution of G’ and a normal distribution for α, and (3) an inverse relationship between G’ and α for individual cells. Moreover, in our AFM measurements we apply powerlaw rheology to individual cells, whereas with our microfluidic technique, we probe each cell only at a single frequency and therefore cannot directly demonstrate that individual cells show that same powerlaw behavior as the cell population. For these reasons, we still think that the comparison with AFM supports the validity of our approach and method.
The novelty of the technique compared to Fregin NatComm 2019 has to be discussed in detail as the main claim of the article is the approach for cell mechanics measurements.
This was also asked by the editor when we first submitted our manuscript, and we explained the differences to the method of Fregin and the novelty/advantages of our method at various places in Introduction and Discussion. These additions are highlighted in the manuscript in red.
From a physics point of view, the technique presented here gives for instance G' and G' as a function of frequency, whereas Fregin et al., give only effective cell elasticity and viscosity. Therefore, the measurement of cell's mechanics is more complete with this assay, as described in the manuscript l.217. Note that the experimental setup presented in the manuscript is simpler than the one presented by Fregin et al., which might favor its use by other labs. This justifies the publication of the article in eLife, as physical methods are in the scope of the journal. On the other hand, such a detailed physical description might not be required to only discriminate between cell populations, as it is the main application of this type of experiment.
The introduction is very short and the authors start presenting their technique immediately. I think the readers need more context on why these type of measurements are useful, which experimental techniques are used and the mechanical models on which these techniques rely. Some of these aspects are in the Discussion section, to highlight the interest of their technique but the introduction should be more detailed.
We followed the suggestion, and in Introduction we now briefly discuss two relevant models (the KelvinVoight model with constant parameters, as used by Fregin et al., and the structural damping model that we favor). While this is not meant to represent a comprehensive discussion of models, we hope it helps the reader to better understand what follows. We also added an example that we hope convincingly illustrates how cell rheology measurements might be interesting and useful even for potential readers outside the mechanobiology community.
On the tank treading movement of cells in shear flows, several articles are cited but the article comports no discussion at all on this mechanism, and how it has been characterized for red blood cells for instance.
We expanded our Introduction to explain the mechanism, included two more references, and also expanded our discussion on tanktreading.
Similarly, some readers may not be familiar with viscoelastic models, the article is not very pedagogic on this.
We also expanded the Introduction in this regard and have introduced the KelvinVoigt model and the structural damping model, both of which are relevant for what follows.
Reviewer #2 (Recommendations for the authors):
While I really like the method, and think it should be published in eLife after successful revision, I am worried about the reliability of some parts of the paper. Mainly, I think the strain stiffening explanation needs to be better nailed. Why is it not possible to change the conditions, so that the deformation remains in the linear regime throughout the measurement. Even if correct, all the data that is in the strain stiffening regime would then lead to wrong stiffness and powerlaw exponent.
Please see our more detailed explanations below where we argue that our method provides a secant modulus and not a small amplitude differential modulus, but that does not make the measurements or numbers “wrong”. Regarding the existence of a “linear regime”, we have addressed this in several previous studies (Kollmannsberger et al., Nonlinear viscoelasticity of adherent cells is controlled by cytoskeletal tension. Soft Matter, 2011 and Lange et al., Unbiased highprecision cell mechanical measurements with microconstrictions. Biophys J, 2017, Kah et al., HighForce Magnetic Tweezers with HysteresisFree Force Feedback, Biophys J 2020). The insight from these studies is that nonlinearities can start to appear even at small amplitudes of strain or stress.
Would it not make more sense, to only focus on the results as close to the measurements as possible? We have a G' and G' for a given frequency, why do we need to translate this into other parameters that are even less reliable?
It is possible to only report the “raw” G’ and G” of each cell, which is in essence what most other microfluidic cell deformability measurements provide. Our data do not necessarily need to be translated into other parameters (such as stiffness k and fluidity α) if one wants to compare different measurements obtained at similar tanktreading frequencies or time scales. However, as explained above, converting G’ and G” into frequencyindependent stiffness and fluidity values makes it possible to compare our measurements to other data obtained with other techniques that operate at different, nonoverlapping time or frequencyscales.
Besides this main concern, I have a couple of other points that I would like to transmit to the authors to consider for a further improvement of the paper.
– Why is in figure 2a/b no data given for the position close to 0? I did not realize an explanation for this when I read the paper. Is it problematic in this regime to determine the deformation?
At low shear stress, when cell deformations are small and cells appear nearly circular, it is difficult to reliably measure the alignment angle. Errors in the alignment angle affect the accuracy of our G” estimate. Furthermore, when cells in the center of the microfluidic channel are exposed to a nonmonotonic shear stress profile, they can deform (depending on their stiffness) into a bulletlike shape, which Roscoe theory does not consider. Therefore, we only include cells that experience a monotonic shear stress, i.e., that are completely on either side of the channel center. We have included this explanation in Methods.
– In figure 2 g it is shown that for a large enough shear rate the rheometer and the flow cytometer values of viscosity are the same. If I understand this correctly, it also means that only right at the center the shear rate is so low that this matters…
Actually, the data shown in Figure 2g demonstrates a good correspondence of the viscosity values between a coneplate rheometer and our device not only for high shear rates but also for low shear rates (e.g., below 10/s). Regarding the position within the microfluidic channel where shear thinning sets in: this depends on the applied pressure. At lower pressures (e.g., 0.5 bar, see Figure 2e), the shear rate of the suspension fluid remains below 100/s, and hence most of the cells do not experience the shear thinning regime. At 3 bar, by contrast, almost all cells are well in the shearthinning regime of the suspension medium except for those close to the center of the channel. For Roscoe theory, however, it does not matter if the suspension fluid is shear thinning or not, as long as we know the local value of the viscosity, which we do.
however in the distances that correspond to this shear no measurements are given. Is this right? If so, maybe it is worthwhile to mention this in more detail.
Correct, we do not evaluate cells at or near the center of the channel. The reason, however, is not that we cannot accurately determine the shear rate and viscosity at this location. Rather, the reason is the nonmonotonic shear stress profile, as explained above.
– Why seem the stiffness values of the PAA beads to be so pressure dependent as seen in figure 4e.
Actually, the stiffnessvalues we find for PAA beads are nearly pressure (and strain) independent, as expected for a linear material such as PAA (Figure 4e). However, we find that the power law exponent, which is very low for a nearly perfect elastic material such as PAA, tends to increase at large strains, presumably due to poroelastic effects (Figure 4f). Put simply, the PAA polymer network interacts with the background fluid of the PAA hydrogel, similar to the cytoskeleton that is dragged through the cytoplasm when cells are deformed. Hence the viscosity term in Equation 2 captures the high frequency rheology of both, cells and PAA beads.
And in figure S3, the values of the AFM and the flow cytometer should be plotted on the same y scale, or ideally even plotted over each other. What keeps you from adding the values a in b,c, and the same for d and g?
We followed this suggestion and now plot the AFM and flow cytometer data with the same yscale.
– In Equation 2 µ seems not to be defined.
– For the cell experiments with vimentin it is very confusing why you compare WT 3T3 with MEF desminknockouts and knockins of vimentin. It looks like you don't have WT MEF. Why comparing different cell types. This makes it hard to believe that your conclusion is correct. At least all differences between 3T3 and MEFs might be because of the different cell type. Use also MEFs for the 'normal' situation.
The NIH3T3 cell line has been established from mouse embryonic fibroblasts (MEFs) by spontaneous immortalization (Todaro and Green (1963) Quantitative studies of the growth of mouse embryo cells in culture and their development into established lines. J Cell Biol 17, 299313). We followed the corresponding protocol for MEFs obtained from vimentinknockout mouse embryos (GolucciGuyon et al., (1994) Mice lacking vimentin develop and reproduce without an obvious phenotype. Cell 79, 679694). The Vim(/)Des cell line was generated by stable transfection of the Vim(/) cell lines with desmin. Thus, the three cell lines are of the same cell type. This information was given in Methods; we now have also added an explanation in Results.
– Also, there is a statistical test in figure 7 b and c missing. How many cells did you measure for each of the datapoints in 7b,c?
The minimum number of cells was 313, the maximum number was 4884, with an average of 1813 cells per data point and a standard deviation of 1191. The average number is now specified in the figure legend. We also performed a statistical analysis and now show (Figure 7c) that the cytochalasinD response for the vin(/) cells is significantly increased compared to control, whereas the cytochalasinD responses between control and the desmin knockin cells are not significantly different.
– In the discussion you mention a local viscosity of the fluid. This is a bit misleading as it suggests that the viscosity depends on the position, but it depends on the shear rate (which depends on the position for a given pressure).
To clarify, we added “shear dependent” (local shear dependent viscosity).
– It should be more focused on the actual measurement of G' G' at a single frequency. I think we don't know sufficiently well if a single powerlaw can explain the measurements, to use this hypothesis in a way that makes the reader think the method provides information about frequency dependence from a single snapshot.
We agree. In the revised Introduction, we clarify that Roscoe theory makes no predictions about the frequencydependence of G’ and G” for that particular cell. However, we find a powerlaw in the G’ and G” vs frequency relationship of the ensemble, and this suggests that the G’ and G” vs frequency relationship of (most of the) individual cells also follow a powerlaw.
– It would be good to get a reference for the statement that you find similar values for the E50 as other authors (line 269).
We refer to Urbanska et al., Nat Meth 2020.
– Please tone down the statement of line 284286. There is no statistical test, and you compare apples with oranges… don't say that you have shown the effect of vimentin. All you see is the increase in stiffness in the knockin situation.
Regarding the apples and oranges, see our explanations above (all measurements, including those for wildtype cells, were performed with mouse embryonic fibroblastderived cell lines). Besides, we see more than just an increase in stiffness in the knockin situation, namely we find a decrease in the relative effects of actin cytoskeleton depolymerization. In our manuscript, we mainly focus on the latter observation, because such relative changes are thought to be less sensitive towards the different life histories of the WT MEF, Vim ko MEF and Des ki MEF cell lines. Nonetheless, also the fact that both the WT and Des ki MEFs are stiffer than the Vim ko MEFs supports our assertion that “the stable introduction of a cytoplasmic intermediate filament protein into intermediate filamentfree cells restores their cytoskeletal functionality and mechanical stability”. We also added a statistical test Mann WhitneyU, which supports this statement (see Figure 7c).
– It would be interesting to provide a measure of variance in the text describing the different \σ values. (line 502503)
The rather convoluted way of computing the acting shear stress (from the average stress at the extreme points) originated from the time when we still included cells that overlapped with the channel center. There, the shear stress at the cell center can reach zero, but a cell will still deform as the average shear stress is nonzero and approximately the average stress at the extreme points. But since we now exclude all cells that overstep the channel center, we directly use the shear stress at the cell center instead (the deviation between the two ways to calculate the acting stress, however, is less than 1% for most cells). This is now explained in Methods, together with an explanation why we exclude cells near the channel center. We have reevaluated all data accordingly.
– Why not using a nonshear thinning medium to do the measurements? Everything should be much more simple there.
The shearthinning behavior complicates matters only moderately, since we obtain the local viscosity (which is important only for the calculation of G") as a side result of our method. Nonetheless, we tried glycerin, but to achieve sufficiently high viscosities, the water content needs to be low, and the cells appeared to suffer. Methyl cellulose or xanthan gum solutions, two other commonly used cell suspension media which we tried and which can be used, are also shear thinning.
– You treat \eta_0, \tau, and \δ as independent parameters. Is this really the case, or could some depend on each other? Also please give the values you measured for these parameters.
We now show the parameters of eta_0, tau, and δ for different alginate concentrations in Figure 2supplementary figure 2. Regarding the question of the independence or covariance of the parameters – this is difficult to determine in a mathematically clean way as we extract these parameters from a fitting routine that iteratively reduces the error between the measured and modeled flow profil, but there are numerous intermediate calculations involved such as root finding and numerical integration. But this is what we can say:
eta_0 is the low shear stress viscosity and as such determines most sensitively the flow speed at the channel center. By contrast, δ describes the shear thinning behavior at large shear rates and determines the overall shape of the flow curve (which approaches a nearly quadratic function for a Newtonian fluid (δ = 0) and a plugflowlike profile for higher δ values). Because of their distinct effects, eta_0 and δ show a low covariance
Tau sets the shear rate at which the viscosity behavior transits from a Newtonian regime to a shearthinning regime. The value of tau affects both the magnitude and the shape of the velocity profile and therefore shows a higher covariance with both eta and δ. Between consecutive measurements of the same fluid, tau shows the largest relative fluctuations, whereas eta_0 and δ are more stable (the data points in Figure 2supplementary figure 2 represent individual measurements, some of it performed on the same day, i.e. with the same alginate preparation, but also on different days, i.e. each time with new alginate preparations of the same concentration).
– You say that stress stiffening happens, but then you use a model (equations 21,22) that breaks down under these conditions.
In fact, Roscoe theory considers a neoHookean material, which is strain stiffening. Hence, Roscoe theory does not break down in the case of a stress or strainstiffening material (and neither does power law rheology, see for example Kollmannsberger Soft Matter 2011 or Lange BJ 2017). However, the strain stiffening of cells appears to be more pronounced than the stiffening predicted for a neoHookean material. Therefore, our measurements represent an effective secant modulus and not a smallstrain tangential modulus, but nonetheless, also the G’ and G” data for large shear stresses (and higher tanktreading frequencies and larger strains) are meaningful. The cell stiffening behavior that we report in the manuscript (Figure 3c) agrees qualitatively with earlier findings (see for example Kollmannsberger Soft Matter 2011 or Lange BJ 2017) and demonstrates that the method can in principle be used to investigate the nonlinear behavior of cells. We added this information to Discussion. Also, we wish to point out that existing microfluidic cell deformability methods (with the exception of a method that we had developed some years ago, see Lange BJ 2015 and Lange BJ 2017) ignore mechanical nonlinearites, moreover they mostly ignore frequencydependent or strainrate dependent cell behavior, and typically they even ignore dissipate (viscous) properties and only report elastic cell properties. Therefore, we believe that our method represents an advance over these previous methods.
Here you basically say that the measurements you obtain in the high frequency regime do not work. Why are these not excluded, and how do you determine the stress at which we cannot take the G' and G' as unaffected.
We are saying the opposite, namely that measurements obtained in the high frequency regime are valid, but we exclude cells near the channel center that experience a low shear stress and low tanktreading frequency. Regarding the question if there is a shear stress threshold that separates a linear from a nonlinear regime, we have investigated this issue in two previous studies (Kollmannsberger2011, Lange2017) and did not find such a threshold. In other words, cells are nonlinear, and the complex shear modulus is always affected to some degree by the stress magnitude.
– In Line 580 you say the AFM goes from 0.1 to 150Hz, but the measurements you show are up to 10000 rad/s. What is right?
Thank you, we made a mistake and multiplied the frequency twice with 2 pi. This has now been corrected.
https://doi.org/10.7554/eLife.78823.sa2Article and author information
Author details
Funding
Deutsche Forschungsgemeinschaft (TRRSFB 225 subprojects A01)
 Elham Mirzahossein
Horizon 2020 (No 812772)
 Mar Eroles
Horizon 2020 (No 953121)
 Mar Eroles
Deutsche Forschungsgemeinschaft (B07)
 Stephan Gekle
Deutsche Forschungsgemeinschaft (A07)
 Stefan Schrüfer
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
This study was supported by the Deutsche Forschungsgemeinschaft (TRRSFB 225 project 326998133 subprojects A01, A07 and B07), and the European Union’s Horizon 2020 research and innovation programmes No 812772 (project Phys2BioMed, Marie SkłodowskaCurie grant) and No 953121 (project FLAMINGO). We thank Jonas Hazur and Aldo Boccaccini for helpful discussions and for providing the alginate.
Senior Editor
 Anna Akhmanova, Utrecht University, Netherlands
Reviewing Editor
 Alphee Michelot, Institut de Biologie du Développement, France
Reviewers
 Clément Campillo, Université d'Evry, France
 Timo Betz, Universtity of Göttingen, Germany
Publication history
 Preprint posted: January 12, 2022 (view preprint)
 Received: March 21, 2022
 Accepted: August 30, 2022
 Accepted Manuscript published: September 2, 2022 (version 1)
 Version of Record published: October 14, 2022 (version 2)
Copyright
© 2022, Gerum 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

 1,299
 Page views

 373
 Downloads

 0
 Citations
Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.
Download links
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)
Further reading

 Cell Biology
 Chromosomes and Gene Expression
Meiotic chromosome segregation relies on synapsis and crossover recombination between homologous chromosomes. These processes require multiple steps that are coordinated by the meiotic cell cycle and monitored by surveillance mechanisms. In diverse species, failures in chromosome synapsis can trigger a cell cycle delay and/or lead to apoptosis. How this key step in 'homolog engagement' is sensed and transduced by meiotic cells is unknown. Here we report that in C. elegans, recruitment of the Pololike kinase PLK2 to the synaptonemal complex triggers phosphorylation and inactivation of CHK2, an early meiotic kinase required for pairing, synapsis, and doublestrand break induction. Inactivation of CHK2 terminates doublestrand break formation and enables crossover designation and cell cycle progression. These findings illuminate how meiotic cells ensure crossover formation and accurate chromosome segregation.

 Cell Biology
 Physics of Living Systems
Motile cilia are hairlike cell extensions that beat periodically to generate fluid flow along various epithelial tissues within the body. In dense multiciliated carpets, cilia were shown to exhibit a remarkable coordination of their beat in the form of traveling metachronal waves, a phenomenon which supposedly enhances fluid transport. Yet, how cilia coordinate their regular beat in multiciliated epithelia to move fluids remains insufficiently understood, particularly due to lack of rigorous quantification. We combine experiments, novel analysis tools, and theory to address this knowledge gap. To investigate collective dynamics of cilia, we studied zebrafish multiciliated epithelia in the nose and the brain. We focused mainly on the zebrafish nose, due to its conserved properties with other ciliated tissues and its superior accessibility for noninvasive imaging. We revealed that cilia are synchronized only locally and that the size of local synchronization domains increases with the viscosity of the surrounding medium. Even though synchronization is local only, we observed global patterns of traveling metachronal waves across the zebrafish multiciliated epithelium. Intriguingly, these global wave direction patterns are conserved across individual fish, but different for left and right nose, unveiling a chiral asymmetry of metachronal coordination. To understand the implications of synchronization for fluid pumping, we used a computational model of a regular array of cilia. We found that local metachronal synchronization prevents steric collisions, cilia colliding with each other, and improves fluid pumping in dense cilia carpets, but hardly affects the direction of fluid flow. In conclusion, we show that local synchronization together with tissuescale cilia alignment coincide and generate metachronal wave patterns in multiciliated epithelia, which enhance their physiological function of fluid pumping.