Significance: The ability to measure the micro-mechanical properties of biological tissues and biomaterials is crucial for numerous fields of cancer research, including tumor mechanobiology, tumor-targeting drug delivery, and therapeutic development. Aim: Our goal is to provide a renewed perspective on the mainstream techniques used for micro-mechanical evaluation of biological tissues and biomimetic scaffoldings. We specifically focus on portraying the outlook of laser speckle micro-rheology (LSM), a technology that quantifies the mechanical properties of biomaterials and tissues in a rapid, non-contact manner. Approach: First, we briefly explain the motivation and significance of evaluating the tissue micro-mechanics in various fields of basic and translational cancer research and introduce the key concepts and quantitative metrics used to explain the mechanical properties of tissue. This is followed by reviewing the general active and passive themes of measuring micro-mechanics. Next, we focus on LSM and elaborate on the theoretical grounds and working principles of this technique. Then, the perspective for measuring the micro-mechanical properties via LSM is outlined. Finally, we draw an overview picture of LSM in cancer mechanobiology research. Results: With the continued emergence of new approaches for measuring the mechanical attributes of biological tissues, the field of micro-mechanical imaging is at its boom. As one of these competent innovations, LSM presents a tremendous potential for both technical maturation and prospective applications in cancer biomechanics and mechanobiology research. Conclusion: By elaborating the current viewpoint of LSM, we expect to accelerate the expansion of this approach to new territories in both technological domains and applied fields. This renewed perspective on LSM may also serve as a road map for other micro-mechanical measurement concepts to be applied for answering mechanobiological questions. |
1.IntroductionIt is increasingly recognized that the micro-scale viscoelastic properties of tissue confer critical micro-mechanical cues that orchestrate nearly all aspects of cellular function, including growth and differentiation.1,2 This tightly controlled, viscoelastic micro-environment is essential to organ development, wound healing, and normal homeostasis in biological tissues.3–5 Irregular micro-mechanical remodeling of tissue is implicated in a broad spectrum of pathologies, including cardiovascular disease, fibro-proliferative disorders, hematological diseases, and cancer.1,6,7 Moreover, meticulously tuned deformability and compliance are imperative for the design of synthetic scaffolds and hydrogels, to mimic normal and pathological structure and function of the tissue microenvironment, and to adjust the release of nutrients, oxygen, reagents, and therapeutics in drug-delivery applications.8 The development of novel tools for measuring the viscoelastic properties of tissue at microscopic length scales is instrumental for understanding the mechanical attributes of cellular structures and their microenvironment and identifying their relationships to the progression of tumor malignancies. These tools are also imperative for staging of the diseases, based on significant mechanical contrast between normal and abnormal tissue states, guiding drug delivery, and devising new lines of therapeutic strategies that target micro-mechanical aberrations. The complex micro-mechanical milieu of tissues may be characterized by the shear viscoelastic modulus, , where and are the elastic and viscous moduli, respectively, and is the angular loading frequency.9 Mechanical rheometry is the traditional approach for evaluating the of materials. In this invasive process, the specimen is sandwiched between two parallel plates, and the sinusoidal shear strain, , is applied to the specimen. This induces an oscillatory stress to the sample, . The rheometer evaluates the twisting and displacements of the plates and calculates the , as the stress-to-strain ratio, i.e., , for a limited frequency range. Averaging the cumulative response over the whole sample volume in the rheometer obscures the mechanical heterogeneities, and this traditional tool returns only the bulk mechanical response of the specimen. Moreover, this approach calls for a relatively large sample volumes and is not conducive for rare and precious biomaterials and tissue specimens.10 The need for non-invasive micro-mechanical testing has motivated the development of several innovative technologies for evaluating the micro-scale viscoelastic properties of tissues and biomaterials. These techniques involve a variety of “active” approaches that measure tissue response to an extrinsic force, and “passive” approaches that rely on thermal fluctuations within tissues to probe micro-mechanical behavior. Here, we briefly review a variety of active mechanical testing approaches and focus our attention on evolution of passive techniques as novel means for non-invasive micro-mechanical assessment. In particular, we detail our perspectives on laser speckle micro-rheology (LSM), a paradigm for realizing micro-mechanical characterization and elaborate on new basic and translational opportunities in cancer mechanobiology research. 2.Measuring Micro-Mechanics2.1.Active MethodsActive micro-mechanical sensing methods frequently involve physical contact with the specimen to apply external forces to induce displacement and deformation. These active methods frequently involve sophisticated instrumentations, with the advantage of applying large forces to probe stiff materials beyond the linear viscoelastic regime. The conventional tool for measuring the micro-scale viscoelastic properties of materials is the atomic force microscopy (AFM)-based indentation. In AFM, a flexible cantilever harboring a micron-sized tip is mounted on a precision stage that tracks the tip displacement as it pokes the specimen. Optical position sensing is used to retrieve the cantilever deflection and calculate the force. Fitting an appropriate model to the force–distance curve yields the elastic indentation modulus () of the specimen.11 Due to its contact-based nature, AFM is inherently invasive. Because the measurements are exuberantly long, they are frequently limited to infinitesimal field of views of a few 10s of . Moreover, AFM only reflects the mechanical properties over depths of 10 s of . In addition, AFM measurements are subject to adhesion and other electrostatic forces. Finally, AFM-based indentation yields only the elastic modulus at a fixed indentation rate and does not probe the frequency dependence or the dynamic viscosity. To circumvent these limitations, several innovative tools are developed for evaluating the micro-mechanical properties of cells and their microenvironment. For example, magnetic twisting cytometry (MTC) evaluates the force transduction of cell membrane to measure the viscoelastic remodeling of cells.12,13 In MTC, an oscillatory magnetic torque is applied on a magnetic bead, bound to a transmembrane receptor. Bead rotation induces an oscillatory twisting stress within the cell. A magnetometer evaluates the of the cell as the ratio of rotary torque to bead rotation.12 MTC was initially used for demonstrating the integrin-mediated cellular mechano-transduction.12 Later, MTC measurements of bronchial epithelial cells yield the in range () and revealed a weak power-law behavior of , .13 The loss tangent () was 0.5 and varied minimally with frequency, indicating a predominantly elastic behavior.13 The frequency dependence and loss tangent may be readily evaluated in MTC. However, since the bead–cell contact area is not precisely known, this technique is mostly qualitative. In addition, frequency range is limited by sampling rate of magnetometers.13 Optical tweezer active micro-rheology (OTAM) is another active technique that employs a highly focused laser beam to trap small dielectric particles.14 By moving the beam, particles apply a local stress to the surrounding material. Sinusoidal oscillation of the beads in the optical trap and tracing the nano-meter scale motion of the bead yield local over a wide range of to ().14 OTAM has been applied for investigating the cell-extra-cellular matrix (ECM) micro-mechanical coupling and hemostasis and showed that at high frequencies (), and exhibits distinct power-law behaviors, i.e., , within cells () and the ECM ().14 Moreover, the frequency, , at which the mechanical properties transition from elastic to viscous (i.e., ) is drastically smaller in cells (2 kHz) compared to ECM (7 kHz).14 Despite its wide frequency range, the small magnitude of applied forces (pN range) limits this approach to compliant specimens. In addition, because of the high NA, local heating and phototoxic effects are inevitable. Quantitative micro-elastography (QME) is a novel variant of compression optical coherence elastography (OCE)15 The parent compression OCE that maps the tissue strain/deformation is response to a static applied force and provides a qualitative elastogram.16 By incorporating a compliant stress sensor, QME permits mapping the absolute Young’s modulus. QME exhibits up to 100-fold contrast improvement by distinguishing the tissue components that have similar deformation but distinct elasticities.15 In highly heterogeneous tissue, the stress varies with depth, violating the QME assumptions.17 Inverse methods may be used to back out the stress distribution in the tissue volume, based on the surface stress, depth-resolved strain, and micro-structural information afforded by the baseline OCT images.18 Because the gradient of displacement is calculated over multiple pixels along the depth, axial resolution of QME is reduced 5 to 10 times compared to the underlying OCT to .19 The lateral resolution is also reduced because the incompressible stress sensor translates a step change in elasticity within the sample (i.e., a sharp feature boundary) to a stress gradient. As such, the stress–strain ratio presents a blurred map of elastic modulus, particularly in deeper sections.16 2.2.Passive MethodsEndogenous motions and vibrations may also be used as a source of deformation, enabling all-optical, non-contact mapping of the viscoelastic properties within the biological tissue.10 For passive measurements, materials must be sufficiently soft to permit detectable motions in response to subtle internal fluctuations. The passive methods may be advantageous compared to active techniques in that the measurements are within the linear viscoelastic regime.20 In addition, because passive methods do not require a loading mechanism, the instrumentation is often simplified, and the hardware frequently has a smaller footprint compared to active techniques. Moreover, since no physical contact is needed, these methods are often non-invasive. Here, we review some of the prominent passive techniques. Brillouin microscopy is based on the inelastic scattering of light, due to its interaction with the local spontaneous acoustic waves within the tissue.21 Brillouin frequency shift is proportionate to the velocity of these intrinsic acoustic vibrations, which are in turn proportionate to the square root of the Young’s modulus at the GHz frequency range, with the exact relationship being specific to the tissue or cell type.22 Implemented using a confocal geometry, Brillouin microscopy provides a non-destructive, label- and contact-free method for probing the viscoelastic properties of biological samples with diffraction-limited axial and lateral resolutions (1.5 and ), as demonstrated in the context of in-vivo tissue and in-vitro cell levels.22,23 Particle tracking micro-rheology (PTM) is another passive approach that involves injecting fluorescently labeled sub-micron beads into the cytoplasm of live cells.24,25 High-frame rate fluorescent video microscopy is used to track the motion of particles. Statistical averaging of particle trajectories yields the mean-square displacement (MSD) .24 Replacing the MSD in the generalized Stokes–Einstein relation (GSER) returns the complex frequency-dependent of the cytoplasm over the frequency range of 0.1 to .26 The upper frequency is limited to the rate at which the particle position is recorded, i.e., the frame rate of the camera.25 Because PTM permits obtaining the entire particle trajectory, it is further possible to analyze the individual particles’ trajectories beyond the spatially averaged MSD and gain insight into the local micro-rheology of the cytoplasm.10,26 Dynamic light scattering (DLS) may also be used for evaluating the Brownian displacements of extrinsic particles.27 In DLS, time-dependent intensity fluctuations of light, that is singly scattered from extrinsic particles, are collected through a pinhole, captured by a photodiode, and analyzed by a digital correlator.20,27 Due to the single-scattered nature of collected light, the light intensity fluctuates very slowly and needs to be analyzed for long times ( or tens of thousands of seconds), to detect perceptible MSD of the tracer particles.10 In addition, the sample must be transparent, homogeneous, and ergodic, so that the extensive temporal averaging provides a statistically accurate ensemble averaging.10 Moreover, large particles’ displacements are needed to sufficiently change the optical path length and provoke a noticeable phase shift and a detectable speckle intensity variation, yielding DLS insensitive to smaller MSD of particles in stiffer samples.27 Development of the diffusing wave spectroscopy (DWS) enabled evaluating the opaque, highly multiple-scattering media, where light propagation is nearly diffusive.27 Nevertheless, DLS and DWS are limited to the extremes of transparent or highly opaque specimens. LSM is a novel optical approach that affords mapping the micro-scale viscoelastic properties of biomaterials and tissues without using extrinsic particles, in a passive, non-invasive manner.28–30 Speckle is a grainy intensity pattern that forms by the self-interference of coherent laser beam, as it backscatters from the turbid materials, including opaque tissue.31 Thermal Brownian motion of intrinsic light scattering particles continuously alter the relative optical phase shifts of backscattered rays and provoke temporal speckle intensity fluctuation. These fluctuations are exquisitely sensitive to particle displacements, and in-turn to the viscoelastic properties of the tissue particles’ microenvironment.6,28,30,32–38 Because in LSM light is multiply scattered from the tissue particles, even minute motions (fraction of a wavelength, in the order of nm) of particles, encountered within each lightpath, give rise to cumulative phase shifts that induce perceptible speckle intensity fluctuations. The acute susceptibility of speckle intensity to particle displacements as small as a few angstroms in LSM makes this technique capable of characterizing the mechanical properties in stiffer samples, in which particle displacements are reduced to subtle vibrations. 3.Laser Speckle Micro-Rheology: Technological Perspectives3.1.Principles of OperationA typical LSM optical setup is displayed in Fig. 1(a).6,28,32–35,38,39 In short, a polarized laser beam is focused on the sample surface. The backscattered speckle patterns are filtered by a linear polarizer and focused by a macro-lens on to a high-speed CMOS camera sensor.6,28,32–35,38,39 The choice of camera frame rate and acquisition time are motivated by the range of frequencies of interests, over which the is to be evaluated. They are further adjusted according to the sample dynamics and relaxation times, such that the speckle patterns exhibit a high contrast between dark and bright spots, and that the intensity autocorrelation curve fully decorrelates well before the end of acquisition.6,32–34 In LSM, quantitative analysis of the speckle dynamics returns the viscoelastic properties of the tissue microenvironment. The working principles of LSM are displayed in the flowchart of Fig. 1(b) as previously detailed in several publications. Briefly, the speckle intensity temporal autocorrelation curve, , is obtained by measuring the correlation between pixel intensities over the time series of speckle images according to the following equation: where and are the intensities of individual pixels at times and , respectively, and represents spatio-temporal averaging over the ensemble of all pixels in the frames and over all reference frame times, . From the curve, the MSD, i.e., , of scattering particles may be retrieved. For low viscosity liquids, scattering particles simply diffuse throughout the sample and this random motion corresponds to a linear trend for MSD. On the other hand, in highly elastic solids, the particle displacements reduce to sluggish vibrations around the equilibrium, which is evidenced by a saturated MSD.32,34,35 For viscoelastic materials, the MSD curve takes on more complex forms as a function of the lag time. Once the MSD is evaluated, the GSER is used to extract the complex shear modulus, as a function of frequency as follows: where is the Fourier transform and is the scattering particle radius. Moreover, stands for the Boltzmann constant (), is the temperature (degrees Kelvin), and corresponds to the frequency at which the is evaluated, and conceptually parallels the mechanical oscillation frequency in the conventional rheometry as detailed above.6,20,33–35 Due to the discrete and finite nature of the data points, and to avoid numerical errors at frequency limits, an algebraic approximation of Fourier transform is frequently used to calculate the . This is achieved by fitting the MSD to a power-law form, i.e., , where is the log–log slope of MSD. Subsequently, the complex may be expressed as The equation above returns both the elastic (storage) and viscous (loss) moduli, i.e., the real and imaginary parts of the modulus. Here, is the gamma function.6,20,33–354.4.Influence of Optical Properties and Scattering Particle Size DistributionIn biological tissues, evaluating the from the measured curve is not straightforward due to contributions of two main factors.6,33 The first factor is the influence of optical absorption and scattering, as speckle fluctuations are driven by both particles’ displacements and the number of particles encountered in photons’ paths. For example, Fig. 2(a) displays the curves, obtained using Eq. (1), for aqueous glycerol mixtures of 90% glycerol and 10% water, dispersed with varying volume fractions of scattering particles, ranging from 0.04% to 2%, with corresponding values of 1.3 to (). The theoretical DLS and DWS curves (dotted lines) are also displayed. As the concentration of scattering particles is increased, the curves vary between the DLS and DWS limits, for single- and rich multiple-scattering, demonstrtaing that the curve depends on both mechancial and optical properties of the sample. The latter is determined by the optical absorption and reduced scattering coefficients (, ) of the tissue.6,33 In other words, when increases, light rays experience larger number of scattering events. As such, the speckle fluctuations increase and the decays more rapidly. Conversely, when increases, the longer optical paths that involve larger number of particles interaction are pruned due to absorption. As such, the speckle fluctuations are reduced and the decays slowly. For biological tissues, and are not known beforehand. To address this, speckle frames may be temporally averaged to calculate the diffuse reflectance profile (DRP) of the specimen. From the DRP, the optical properties may be deduced in two ways. The first approach involves fitting a model, derived from light diffusion approximation, to the radial variations of DRP from the beam focus point, which yields both of and .6,33,40 Nevertheless, curve is in effect modulated by the ratio of . This ratio maybe conveniently extracted from the total reflectance of the specimen, which can be obtained by comparing the DRP of the specimen to that of a standard reflector.34,40 The second factor in calculating the from MSD via the GSER is that size of scattering particles is needed.34 Conceptually, in addition to optical properties, size of scattering particles also modifies the speckle fluctuations, with larger particles presenting slower dynamics.34 Similar to optical properties, the average size of scattering particles may also be derived from temporally averaged speckle frames, i.e., diffuse reflectance. However, in this case, the speckle patterns may be acquired by placing a linear polarizing filter in front of the camera that is oriented parallel to the polarizer in the path of the illumination laser beam. Under this condition, the DRP exhibits polar angle variations that reveal the average scattering size, .34,35 More specifically, when the average particle size increases in the 0.1- to range, DRP evolves from a bi-lobular shape to a clover-like pattern, with the size of the additional lobes increasing with the particle size, as seen in Fig. 2(b).34 The size-dependent changes of DRP patterns are likely due to transition between isotropic Rayleigh scattering to forwardly directed Mie scattering, as the particle size increases compared to the laser source wavelength. Figures 3(a) and 3(b) display representative examples of LSM measurements in viscoelastic biofluids and predominantly elastic hydrogels, calculated by following the steps outlined in the flowchart of Fig. 1(b) and through compensating for variations of optical properties and scattering particle size, as discussed above.6,34,35 The corresponding conventional mechanical rheometer measurements are also displayed, yet are only valid at lower frequencies, where the inertial effects are negligible. The onset of inertial effects is variable among different specimens and is usually lower in low viscosity liquid, as elaborated in the next section. These figures demonstrate the capability of LSM in evaluating specimens of assorted viscoelastic and elastoviscous behaviors over an expanded frequency range. Figures 3(c) and 3(d) display the strong and statistically significant correlation between LSM measurements versus rheometer and AFM in assorted hydrogel phantoms, which highlights the large dynamic range that exceeds to 40 kPa at .35 4.Perspectives for Micro-Mechanical Measurements4.1.Dynamic Range of Viscoelastic Modulus MeasurementsBecause in LSM light is multiply scattered from the tissue particles, even minute motions (fraction of a wavelength, in the order of nm) of particles, encountered within each lightpath, give rise to cumulative phase shifts that induce perceptible speckle intensity fluctuations. The rapidly fluctuating speckle patterns are captured via a high-speed camera, for a few seconds at a time to estimate the MSD. The acute susceptibility of speckle intensity to particle displacements as small as a few angstroms in LSM makes this technique capable of characterizing the mechanical properties in stiffer samples, in which particle displacements are reduced to subtle vibrations. In particular, the upper limit of viscoelastic modulus accessible to LSM is set by the size of scattering particles, , the ability to resolve their infinitesimal motions, , and the thermal energy, , to .10 The highly sensitive multi-speckle detection scheme of LSM enables resolving displacements of and permits probing the viscoelastic modulus over 5 to 6 decades of magnitude, with an upper limit in the order of to 100 kPa. In comparison, other passive techniques such as PTM resolve only considerably larger displacements via video microscopy and exhibit much smaller accessible range of in the order of a few Pa.24,26 These distinct features of LSM permit extending the passive micro-rheology to the context of stiffer biological tissues, in which intrinsic scattering particles exhibit arbitrary concentrations and size distributions. 4.2.Frequency LimitsAs observed in Figs. 3(a) and 3(b), the upper frequency limit is reached with the advent of inertial effects in the conventional rheolometry measurements. For the rheometer, the shear strain waves applied to the sample via a parallel plate decay exponentially as they penetrate the specimen because of inertial effects. The penetration depth of shear strain waves is given by , where is the density and is the upper frequency limit, suggesting that is inversely proportionate to the frequency. When < typical gap size or sample thickness, the inertial effects dominate, and the measurements become unreliable. For low viscosity fluids, this correspond to upper frequency limit of and 100 Hz, as evidenced by the curves of Fig. 3.10,41 Likewise, in LSM the shear waves propagated by the motion of the Brownian scattering particles decay exponentially from the surface of the bead through the sample microenvironment. Nevertheless, the inertial effects only dominate when the , where is the scattering particle size. Therefore, in LSM the small, sub-micron size of tissue scattering particles postpones the onset of inertial effect to -regime, opening a substantially larger window of timescales and frequencies, compared to both the rheometer and other micro-mechanical testing techniques.10,41 The upper frequency range of competing micro-mechanical measurement techniques, namely PTM, MTC, and OTAM, are 10 Hz, 16, and 15 kHz, respectively.12–14,24,25 While LSM may in principle probe such a wide range of sample dynamics, in practice the high-frequency limit is further influenced by the CMOS camera; when operated at high frame rate, the camera can provide finer temporal resolution extending the frequency range of LSM. In the typical studies, a frame rate of 750 fps for instance would limit the higher frequency boundary to within 100s of Hz. By employing acquisition speeds, for instance in the order of few 100 kHz, higher frequencies in the order of may be achieved, which are not accessible to standard mechanical rheometry. The low-frequency bound, , is determined by the time scale at which longitudinal modes become significant compared to the shear modes, excited in the system.24,26,42 In mechanical rheology, the top plate only applies a shear strain to the sample. Therefore, in principle, no lower frequency limit exists in the case of mechanical rheology. Nevertheless, in practice evaluating the at frequencies below 0.1 Hz may require tediously long times. In LSM, on the other hand, scattering particles respond to all the thermally excited modes, including the longitudinal modes of the elastic network. At high frequencies, the viscoelastic material may be modeled as an elastic network that is viscously coupled to the surrounding liquid. Due to the incompressible nature of the viscous liquid, no compressional strain is sustained in the microenvironment and the Brownian displacements are entirely due to excited shear modes. At low strain rates, the motion of network and liquid is decoupled and the liquid drains freely through the network. The onset of longitudinal modes is set by the ratio of sample elasticity and viscosity, as well as the ratio of network mesh to Brownian particle size, as . Here, is the elastic modulus of the material, is the pore size, is the viscosity of the specimen, and a is the scattering particle size. This translates into a frequency of about 0.159 Hz for a typical soft material where elastic modulus is about 100 times larger than the viscosity, and the mesh size is one-tenth of the radius of the embedded probe. 4.3.Viscous and Elastic Moduli and their Frequency DependenceThe LSM evaluates the complex and is capable of measuring both the elastic (storage) and viscous (loss) moduli as a function of frequency that form the real and imaginary parts of the complex modulus, i.e., . More specifically, once the is evaluated through the GSER equation, these two moduli may be expressed as where is the log–log slope of the MSD at . By affording the ability to measure both viscous and elastic moduli, the LSM surmounts the other micro-mechanical elastography techniques, such as QME, that only evaluate the elastic response of the tissue microenvironment. This is particularly important in the context of cellular mechanobiology studies, where the interplay between micro-scale elastic and viscous moduli of cellular microenvironment is hypothesized to modulate the oncogenic signaling pathways in terms that may not be predicted solely based on the elastic properties.43,44 Moreover, in the context of whole tissue, the complex modulus may not always be explained by alone, and the contribution of loss and storage moduli to the perceived “stiffness” depends on the specific tissue type and the time-scale/loading frequency of measurements. The LSM in principle permits evaluating the various aspects of viscoelasticity and frequency-dependence, including the power-law behavior of and , (i.e., ), and the frequency dependence of loss tangent (, which present invaluable diagnostic yields. For instance, the ability to measure both elastic and viscous moduli, together with the large frequency range of LSM, permits identifying the cross-over frequency, , at which the mechanical properties transition from elastic to viscous (sol–gel transition). The onset of , i.e., transition to viscous behavior, is expected to be smaller in cells compared to the ECM.14 In addition, a reduced is hypothesized to be a critical indicator of malignant behavior.14,45When using the algebraic approximation of GSER in LSM, we inherently assumed that MSD behavior in time domain and frequency-dependence is in symmetry. Therefore, the log-slope of MSD, i.e., α, directly returns both the frequency-dependent exponent of power-law behavior and the loss tangent of the , as . Therefore, the power-law exponent may be readily probed from the MSD and no additional curve-fitting steps are needed. For a purely elastic material, particle displacements reduce to subtle vibration, MSD is constant, and , which implies that . On the other hand, for a purely viscous liquid, particle displacements are diffusive, MSD grows linearly with time, , and . Therefore, when MSD is only due to passive Brownian motion, a physically realistic range for is 0 to 1. Nevertheless, bulk motion of tissue due to environmental vibrations could induce actively driven, super-diffusive displacements and increase beyond 1.46 For instance, in the extreme case of a laminarly flowing purely viscous biofluid reaches 2. Therefore, even minute drifts that are not perceptible in magnitude of MSD, and in turn, could significantly influence the and in turn the evaluated and . Thus, an effective vibration cancellation mechanism is imperative for accurately measuring and in LSM. 5.LSM in Cancer Mechanobiology ResearchThe stiffness of the extracellular matrix or tumor cell substrate has been frequently used to explain the mechano-regulation of oncogenic signaling pathways that drive the hallmarks of cancer, such as proliferation, morphogenesis, migration, and invasion, as well as vascular modulation and immune response alterations47 Cancer mechanobiology research has thus utilized mechanical testing methods such as conventional rheometry that measures the bulk or average stiffness of the ECM. The tumor microenvironment, however, is far from a homogeneous solid.2,48 It is speculated that cellular mechano-sensors perceive both the elastic stiffness and viscous dissipation at the microscale, and respond to the superposition of these cues, by either activating or suppressing their malignant traits, in terms that are not sufficiently explained by the baseline bulk stiffness.43,44 Micromechanical heterogeneity of the tumor microenvironment limits the uniform and effective penetration of the cytotoxic drugs around and within the tumor and thus affects the treatment outcome.49–51 It further restricts the even distribution of oxygen and nutrient within the tumor, and thus maintains a hypoxic niche for proliferation of drug-resistant cancer stem cells.52 In this context, the capability of the LSM for high-resolution, non-contact mapping of tumor micro-environment within in-vitro culture systems of cancer cells grown on engineered ECM substrates, under sterile condition is instrumental for providing a renewed and refined perspective of tumor mechanobiology. Furthermore, the susceptibility of LSM measurements to sub-wavelength scattering particle displacements, and in turn sensitivity to moduli changes of less than 1 Pa, makes this technology highly desirable for longitudinal studies of micro-mechanical modulation. Figures 4(a) and 4(b) demonstrate the spatial resolution of LSM micro-mechanical maps in a micro-fabricated polydimethylsiloxane-polyethylene glycol diacrylate (PDMS-PEGDA) phantoms, featuring highly stiff PDMS bars of 1-mm long, and 200, 150, 100, 80, and wide within a background of compliant PEGDA 5% gel.35 From this figure, it is evident that LSM resolves micro-mechanical features in the order of a few 10s . This is comparable with the length-scales of cell-ECM micro-mechanical interactions and speaks to utility of LSM for probing mechanobiological processes at the interface of cells and their microenvironment. We have previously conducted a longitudinal study of micro-mechanical remodeling in an in-vitro OvCar5 cells, embedded in a Matrigel ECM,53 as they aggressively migrated and assembled into large multicellular 3D tumor micro-spheroid clusters in the span of 21 days (unpublished results). Figures 4(c) and 4(d) displays the bright field image of a micro-metastatic ovarian cancer nodules at the seventh day of the study, along with the corresponding spatial map of , evaluated by LSM. Mechanically stiffer regions are apparent even at the cellular levels and a mechanical transformation from surrounding matrix is observed. An important consideration when conducting LSM in cell culture systems is the influence of active intra-cellular motion on the speckle intensity fluctuations. Prior studies by other groups have investigated the distinctions between the time scales of Brownian motions and these active processes. In particular, speckle fluctuations caused by diffusion of gold particles in the ECM and intracellular adenosine triphosphate (ATP) motors have been previously studied.54 Accordingly, the time scale of ATP-driven speckle decorrelation is in the order of 5 s and is captured only when the frame rate is . Given that LSM acquires speckle images at much higher frames rates in the order of 100s and 1000s of frames per second, the contribution of non-Brownian dynamics to the LSM measurements is minimal. We independently verified this by evaluating 3 plane 4% agarose gels and 3 others seeded with 3T3 fibroblast cells. T-test declared no significant difference in , evaluated at and 100 Hz between the two groups (). These studies open exciting opportunities for future applications of LSM to obtain insights on the mechanically mediated cross-talk between cells and their 3D microenvironment with high resolution and sensitivity. LSM may further assist with the development of new lines of therapeutics that focus on alleviating micro-mechanical abnormalities. The traditional understanding of tumor mechanobiology has suggested several mechano-therapies that aimed to regress the course of cancer by reducing the tumor stiffness. For instance, anti-fibrotic reagents and matrix metalloproteinase inhibitors are proposed to target the ECM micromechanical properties.55,56 Nevertheless, these approaches largely failed in clinical trials, likely due to extensive toxicity associated with off-target interruption of tissue functions in healthy organs.2,57 LSM affords the capability to evaluate the micro-mechanical properties of tumor cells cultured on engineered substrates, and permits investigating the interaction and interplay of cell–ECM viscoelastic properties. In this way, it could assist developing successful therapeutic strategies based on both restoration of the ECM micro-mechanical properties and disruption of the aberrant cellular responses to ideally strike all the modifiable targets at the same time.2 By quantifying the changes in micro-mechanical properties of the ECM, LSM could facilitate the development of such optimal interventions for clinical benefit across multiple diseases. 6.Summary and ConclusionWe recapitulated the ability of LSM for evaluating key quantitative metrics of elastic, viscous, and viscoelastic moduli as well as loss tangent of tissues and biomaterials. We further renumerated the unique features of this technology, including large dynamic-range, unparalleled frequency limits, and exquisite sensitivity, which permit investigating fundamental questions on micromechanical aspects of cancer etiology and drug resistance. Beyond cancer research applications, the LSM can also be used to study mechanosensitive processes pertinent in wound healing, atherosclerosis, as well as neurodegenerative and orthopedic diseases.49,58–61 We envision these insights will open new perspectives for rheology-informed disease prognosis and drug development for improved therapeutic efficacy in the future. DisclosuresDr. Nadkarni and Dr. Hajjarian are inventors of coagulation sensing technology to be licensed to Coalesenz, Inc., which is developing a point-of-care system to measure coagulation parameters. Dr. Nadkarni has equity and serves on the board of the Coalesenz, Inc. Dr. Nadkarni’s interests were reviewed and are managed by Massachusetts General Hospital and Partners Healthcare in accordance with their conflict-of-interest policies. AcknowledgmentsThe authors acknowledge research funding from the National Institutes of Health (No. 1R01HL119867-01). ReferencesT. R. Cox and J. T. Erler,
“Remodeling and homeostasis of the extracellular matrix: implications for fibrotic diseases and cancer,”
Disease Models Mech., 4 165
–178
(2011). https://doi.org/10.1242/dmm.004077 Google Scholar
M. C. Lampi and C. A. Reinhart-King,
“Targeting extracellular matrix stiffness to attenuate disease: from molecular mechanisms to clinical trials,”
Sci. Transl. Med., 10 eaao0475
(2018). https://doi.org/10.1126/scitranslmed.aao0475 STMCBQ 1946-6234 Google Scholar
D. E. Ingber,
“Mechanical control of tissue morphogenesis during embryological development,”
Int. J. Dev. Biol., 50 255
–266
(2006). https://doi.org/10.1387/ijdb.052044di IJDBE5 0214-6282 Google Scholar
E. Farge,
“Mechanotransduction in development,”
Curr. Top. Dev. Biol., 95 243
–265
(2011). https://doi.org/10.1016/B978-0-12-385065-2.00008-6 CTDBA5 0070-2153 Google Scholar
L. A. Barnes et al.,
“Mechanical forces in cutaneous wound healing: emerging therapies to minimize scar formation,”
Adv. Wound Care, 7 47
–56
(2018). https://doi.org/10.1089/wound.2016.0709 Google Scholar
Z. Hajjarian and S. K. Nadkarni,
“Evaluation and correction for optical scattering variations in laser speckle rheology of biological fluids,”
PLoS One, 8 e65014
(2013). https://doi.org/10.1371/journal.pone.0065014 POLNCL 1932-6203 Google Scholar
T. Saneyasu, R. Akhtar and T. Sakai,
“Molecular cues guiding matrix stiffness in liver fibrosis,”
Biomed. Res. Int., 2016 1
–11
(2016). https://doi.org/10.1155/2016/2646212 Google Scholar
T. R. Hoare and D. S. Kohane,
“Hydrogels in drug delivery: progress and challenges,”
Polymer, 49 1993
–2007
(2008). https://doi.org/10.1016/j.polymer.2008.01.027 POLMAG 0032-3861 Google Scholar
A. P. Sarvazyan et al.,
“Shear wave elasticity imaging: a new ultrasonic technology of medical diagnostics,”
Ultrasound Med. Biol., 24 1419
–1435
(1998). https://doi.org/10.1016/S0301-5629(98)00110-0 USMBA3 0301-5629 Google Scholar
M. V. Gardel, D. Weitz,
“Microrheology,”
Microscale Diagnostic Techniques, Springer, Berlin, Heidelberg
(2005). Google Scholar
A. J. Engler et al.,
“Microtissue elasticity: measurements by atomic force microscopy and its influence on cell differentiation,”
Methods Cell Biol., 83 521
–545
(2007). https://doi.org/10.1016/S0091-679X(07)83022-6 MCBLAG 0091-679X Google Scholar
N. Wang and D. E. Ingber,
“Probing transmembrane mechanical coupling and cytomechanics using magnetic twisting cytometry,”
Biochem. Cell Biol., 73 327
–335
(1995). https://doi.org/10.1139/o95-041 Google Scholar
M. Puig-De-Morales et al.,
“Measurement of cell microrheology by magnetic twisting cytometry with frequency domain demodulation,”
J. Appl. Physiol., 91 1152
–1159
(2001). https://doi.org/10.1152/jappl.2001.91.3.1152 Google Scholar
J. R. Staunton et al.,
“High-frequency microrheology in 3D reveals mismatch between cytoskeletal and extracellular matrix mechanics,”
Proc. Natl. Acad. Sci. U.S.A., 116 14448
–14455
(2019). https://doi.org/10.1073/pnas.1814271116 Google Scholar
M. S. Hepburn et al.,
“Three-dimensional imaging of cell and extracellular matrix elasticity using quantitative micro-elastography,”
Biomed. Opt. Express, 11 867
–884
(2020). https://doi.org/10.1364/BOE.383419 BOEICL 2156-7085 Google Scholar
B. F. Kennedy et al.,
“Investigation of optical coherence microelastography as a method to visualize cancers in human breast tissue,”
Cancer Res., 75 3236
–3245
(2015). https://doi.org/10.1158/0008-5472.CAN-14-3694 CNREA8 0008-5472 Google Scholar
L. Chin et al.,
“Simplifying the assessment of human breast cancer by mapping a micro-scale heterogeneity index in optical coherence elastography,”
J. Biophotonics, 10 690
–700
(2017). https://doi.org/10.1002/jbio.201600092 Google Scholar
M. M. Doyley,
“Model-based elastography: a survey of approaches to the inverse elasticity problem,”
Phys. Med. Biol., 57 R35
–R73
(2012). https://doi.org/10.1088/0031-9155/57/3/R35 PHMBA7 0031-9155 Google Scholar
B. F. Kennedy, K. M. Kennedy and D. D. Sampson,
“A review of optical coherence elastography: fundamentals, techniques and prospects,”
IEEE J. Sel. Top. Quantum Electron., 20
(2), 7101217
(2014). https://doi.org/10.1109/JSTQE.2013.2291445 IJSQEN 1077-260X Google Scholar
T. G. Mason and D. A. Weitz,
“Optical measurements of frequency-dependent linear viscoelasticity moduli of complex fluids,”
Phys. Rev. Lett., 74 1250
–1253
(1995). https://doi.org/10.1103/PhysRevLett.74.1250 PRLTAO 0031-9007 Google Scholar
L. Brillouin,
“Diffusion de la lumière et des rayons X par un corps transparent homogène,”
Ann. Phys., 9 88
–122
(1922). https://doi.org/10.1051/anphys/192209170088 Google Scholar
G. Scarcelli et al.,
“Noncontact three-dimensional mapping of intracellular hydromechanical properties by Brillouin microscopy,”
Nat. Methods, 12 1132
–1134
(2015). https://doi.org/10.1038/nmeth.3616 1548-7091 Google Scholar
G. Scarcelli and S. H. Yun,
“In vivo Brillouin optical microscopy of the human eye,”
Opt. Express, 20 9197
–9202
(2012). https://doi.org/10.1364/OE.20.009197 OPEXFF 1094-4087 Google Scholar
J. A. McGlynn, N. Wu and K. M. Schultz,
“Multiple particle tracking microrheological characterization: fundamentals, emerging techniques and applications,”
J. Appl. Phys., 127 201101
(2020). https://doi.org/10.1063/5.0006122 JAPIAU 0021-8979 Google Scholar
D. Wirtz,
““Particle-tracking microrheology of living cells: principles and applications,”
Annu. Rev. Biophys., 38 301
–326
(2009). https://doi.org/10.1146/annurev.biophys.050708.133724 ARBNCV 1936-122X Google Scholar
Y. Tseng, T. P. Kole and D. Wirtz,
“Micromechanical mapping of live cells by multiple-particle-tracking microrheology,”
Biophys. J., 83 3162
–3176
(2002). https://doi.org/10.1016/S0006-3495(02)75319-8 BIOJAU 0006-3495 Google Scholar
D. A. Weitz, D. J. Pine, and W. Brown,
“Diffusing wave spectroscopy,”
Dynamic Light Scattering: The Method and Some Applications, Oxford University Press, New York
(1993). Google Scholar
S. K. Nadkarni et al.,
“Characterization of atherosclerotic plaques by laser speckle imaging,”
Circulation, 112 885
–892
(2005). https://doi.org/10.1161/CIRCULATIONAHA.104.520098 CIRCAZ 0009-7322 Google Scholar
S. K. Nadkarni et al.,
“Laser speckle imaging of atherosclerotic plaques through optical fiber bundles,”
J. Biomed. Opt., 13 054016
(2008). https://doi.org/10.1117/1.2982529 JBOPFO 1083-3668 Google Scholar
Z. Hajjarian and S. Nadkarni,
“Tutorial on laser speckle rheology: technology, applications, and opportunities,”
J. Biomed. Opt., 25 050801
(2020). https://doi.org/10.1117/1.JBO.25.5.050801 Google Scholar
J. W. Goodman, Speckle Phenomena in Optics: Theory and Applications, 2nd edSPIE, Bellingham, Washington
(2020). Google Scholar
Z. Hajjarian and S. K. Nadkarni,
“Evaluating the viscoelastic properties of tissue from laser speckle fluctuations,”
Sci. Rep., 2 316
(2012). https://doi.org/10.1038/srep00316 SRCEC3 2045-2322 Google Scholar
Z. Hajjarian and S. K. Nadkarni,
“Correction of optical absorption and scattering variations in laser speckle rheology measurements,”
Opt. Express, 22 6349
–6361
(2014). https://doi.org/10.1364/OE.22.006349 OPEXFF 1094-4087 Google Scholar
Z. Hajjarian and S. K. Nadkarni,
“Estimation of particle size variations for laser speckle rheology of materials,”
Opt. Lett., 40 764
–767
(2015). https://doi.org/10.1364/OL.40.000764 OPLEDP 0146-9592 Google Scholar
Z. Hajjarian et al.,
“Laser speckle rheology for evaluating the viscoelastic properties of hydrogel scaffolds,”
Sci. Rep., 6 37949
(2016). https://doi.org/10.1038/srep37949 SRCEC3 2045-2322 Google Scholar
Z. Hajjarian et al.,
“In-vivo mechanical characterization of coronary atherosclerotic plaques in living swine using intravascular laser speckle imaging,”
Biomed. Opt. Express, 12 2064
–2078
(2021). https://doi.org/10.1364/BOE.418939 BOEICL 2156-7085 Google Scholar
Z. Hajjarian et al.,
“Intravascular laser speckle imaging catheter for the mechanical evaluation of the arterial wall,”
J. Biomed. Opt., 16 026005
(2011). https://doi.org/10.1117/1.3533322 JBOPFO 1083-3668 Google Scholar
S. K. Nadkarni et al.,
“Measurement of fibrous cap thickness in atherosclerotic plaques by spatiotemporal analysis of laser speckle images,”
J. Biomed. Opt., 11 021006
(2006). https://doi.org/10.1117/1.2186046 JBOPFO 1083-3668 Google Scholar
M. M. Tripathi et al.,
“Assessing blood coagulation status with laser speckle rheology,”
Biomed. Opt. Express, 5 817
–831
(2014). https://doi.org/10.1364/BOE.5.000817 BOEICL 2156-7085 Google Scholar
T. J. Farrell, M. S. Patterson and B. Wilson,
“A diffusion theory model of spatially resolved, steady-state diffuse reflectance for the noninvasive determination of tissue optical properties in vivo,”
Med. Phys., 19 879
–888
(1992). https://doi.org/10.1118/1.596777 MPHYA6 0094-2405 Google Scholar
V. Breedveld and D. J. Pine,
“Microrheology as a tool for high-throughput screening,”
J. Mater. Sci., 38 4461
–4470
(2003). https://doi.org/10.1023/A:1027321232318 JMTSAS 0022-2461 Google Scholar
B. Schnurr et al.,
“Determining microscopic viscoelasticity in flexible and semiflexible polymer networks from thermal fluctuations,”
Macromolecules, 30 7781
–7792
(1997). https://doi.org/10.1021/ma970555n MAMOBX 0024-9297 Google Scholar
E. E. Charrier et al.,
“Elasticity-dependent response of malignant cells to viscous dissipation,”
Biomech. Model. Mechanobiol., 20 145
–154
(2021). https://doi.org/10.1007/s10237-020-01374-9 BMMICD 1617-7959 Google Scholar
O. Chaudhuri et al.,
“Effects of extracellular matrix viscoelasticity on cellular behaviour,”
Nature, 584 535
–546
(2020). https://doi.org/10.1038/s41586-020-2612-2 Google Scholar
R. Sinkus et al.,
“MR elastography of breast lesions: understanding the solid/liquid duality can improve the specificity of contrast-enhanced MR mammography,”
Magn. Reson. Med., 58 1135
–1144
(2007). https://doi.org/10.1002/mrm.21404 MRMEEN 0740-3194 Google Scholar
N. Gal, D. Lechtman-Goldstein and D. Weihs,
“Particle tracking in living cells: a review of the mean square displacement method and beyond,”
Rheol. Acta, 52 425
–443
(2013). https://doi.org/10.1007/s00397-013-0694-6 RHEAAK 0035-4511 Google Scholar
M. W. Pickup, J. K. Mouw and V. M. Weaver,
“The extracellular matrix modulates the hallmarks of cancer,”
EMBO Rep., 15 1243
–1253
(2014). https://doi.org/10.15252/embr.201439246 Google Scholar
D. T. Butcher, T. Alliston and V. M. Weaver,
“A tense situation: forcing tumour progression,”
Nat. Rev. Cancer, 9 108
–122
(2009). https://doi.org/10.1038/nrc2544 NRCAC4 1474-175X Google Scholar
L. E. Barney et al.,
“The predictive link between matrix and metastasis,”
Curr. Opin. Chem. Eng., 11 85
–93
(2016). https://doi.org/10.1016/j.coche.2016.01.001 Google Scholar
J. Dittmer and B. Leyh,
“The impact of tumor stroma on drug response in breast cancer,”
Semin. Cancer Biol., 31 3
–15
(2015). https://doi.org/10.1016/j.semcancer.2014.05.006 SECBE7 1044-579X Google Scholar
R. K. Jain, J. D. Martin and T. Stylianopoulos,
“The role of mechanical forces in tumor growth and therapy,”
Annu. Rev. Biomed. Eng., 16 321
–346
(2014). https://doi.org/10.1146/annurev-bioeng-071813-105259 ARBEF7 1523-9829 Google Scholar
D. M. Gilkes, G. L. Semenza and D. Wirtz,
“Hypoxia and the extracellular matrix: drivers of tumour metastasis,”
Nat. Rev. Cancer, 14 430
–439
(2014). https://doi.org/10.1038/nrc3726 NRCAC4 1474-175X Google Scholar
A. O. Abu-Yousif et al.,
“PuraMatrix encapsulation of cancer cells,”
J. Vis. Exp.,
(34), e1692
(2009). https://doi.org/10.3791/1692 Google Scholar
A. L. Oldenburg et al.,
“Motility-, autocorrelation-, and polarization-sensitive optical coherence tomography discriminates cells and gold nanorods within 3D tissue cultures,”
Opt. Lett., 38 2923
–2926
(2013). https://doi.org/10.1364/OL.38.002923 OPLEDP 0146-9592 Google Scholar
D. E. Ingber,
“Can cancer be reversed by engineering the tumor microenvironment?,”
Semin. Cancer Biol., 18 356
–364
(2008). https://doi.org/10.1016/j.semcancer.2008.03.016 SECBE7 1044-579X Google Scholar
B. Diop-Frimpong et al.,
“Losartan inhibits collagen I synthesis and improves the distribution and efficacy of nanotherapeutics in tumors,”
Proc. Natl. Acad. Sci. U.S.A., 108 2909
–2914
(2011). https://doi.org/10.1073/pnas.1018892108 PNASA6 0027-8424 Google Scholar
A. Winer, S. Adams and P. Mignatti,
“Matrix metalloproteinase inhibitors in cancer therapy: turning past failures into future successes,”
Mol. Cancer Ther., 17 1147
–1155
(2018). https://doi.org/10.1158/1535-7163.MCT-17-0646 Google Scholar
N. D. Evans et al.,
“Substrate stiffness affects early differentiation events in embryonic stem cells,”
Eur. Cell Mater., 18 1
–14
(2009). https://doi.org/10.22203/eCM.v018a01 Google Scholar
C. Y. Yip et al.,
“Calcification by valve interstitial cells is regulated by the stiffness of the extracellular matrix,”
Arterioscler. Thromb. Vasc. Biol., 29 936
–942
(2009). https://doi.org/10.1161/ATVBAHA.108.182394 Google Scholar
F. J. Byfield et al.,
“Endothelial actin and cell stiffness is modulated by substrate stiffness in 2D and 3D,”
J. Biomech., 42 1114
–1119
(2009). https://doi.org/10.1016/j.jbiomech.2009.02.012 JBMCB5 0021-9290 Google Scholar
N. Yamamura et al.,
“Effects of the mechanical properties of collagen gel on the in vitro formation of microvessel networks by endothelial cells,”
Tissue Eng., 13 1443
–1453
(2007). https://doi.org/10.1089/ten.2006.0333 1937-3341 Google Scholar
BiographyZeinab Hajjarian received her PhD in electrical engineering from the Pennsylvania State University in 2009. Following her doctoral work, she completed her post-doctoral fellowship in the lab of Dr. Seemantini Nadkarni at the Wellman Center for Photomedicine. In 2013, she was promoted to an assistant of electrical engineering at the Wellman Center for Photomedicine, Massachusetts General Hospital, and an instructor of dermatology at Harvard Medical School. Her research is focused on developing non-destructive optical micro-rheology techniques for mechanical assessment of biological tissue and studying biomechanical cues in normal development and disease. Seemantini K. Nadkarni is an associate professor at Harvard Medical School and directs her laboratory at the Wellman Center for Photomedicine, Massachusetts General Hospital. The research initiatives in the Nadkarni Lab are focused on inventing optical in-vitro diagnostic devices for whole blood testing at the point-of-care, developing new imaging strategies for quantifying the micromechanical properties of the tumor extracellular matrix, and engineering endoscopic probes for intra-luminal imaging of soft tissue biomechanics for a variety of clinical applications. |