Magnetomotive optical coherence elastography for microrheology of biological tissues

Abstract. Optical coherence elastography (OCE) is an established paradigm for measuring biomechanical properties of tissues and cells noninvasively, in real time, and with high resolution. We present a different development of a spectral domain OCE technique that enables simultaneous measurements of multiple biomechanical parameters of biological tissues. Our approach extends the capabilities of magnetomotive OCE (MM-OCE), which utilizes iron oxide magnetic nanoparticles (MNPs) distributed and embedded in the specimens as transducers for inducing motion. Step-wise application of an external magnetic field results in displacements in the tissue specimens that are deduced from sensitive phase measurements made with the MM-OCE system. We analyzed freshly excised rabbit lung and muscle tissues. We observe that while they present some similarities, rabbit lung and muscle tissue displacements display characteristic differentiating features. Both tissue types undergo a fast initial displacement followed by a rapidly damped oscillation and the onset of creep. However, the damping is faster in muscle compared to lung tissue, while the creep is steeper in muscle. This approach has the potential to become a novel way of performing real-time measurements of biomechanical properties of tissues and to enable the development of different diagnostic and monitoring tools in biology and medicine.


Introduction
The biomechanical properties of tissues and cells are known to contain information about the state of health of the organs or organisms where they reside.A well-known example is the physical palpation of the skin surface and underlying tissue structures to determine their gross biomechanical properties in relation to adjacent regions of tissue.A small local hardening in the breast is often found in a number of pathologies, from benign sebaceous cysts to malignant breast cancer.Thus, there is a great deal of interest in having access to versatile diagnostic tools that could differentiate pathological tissue from normal tissue in unique ways.The field of elastography involves the use of an imaging modality to provide data on the displacement, distortion, or dynamics of tissue under various forms of static or dynamic stress.2][3][4][5][6][7] However, there remains a need to develop diagnostic elastography tools with the capability of measuring biomechanical properties with high resolution on the scale of small tissues and cells.
A number of methods have been reported to assess the biomechanical properties of tissues using optical methods.The high spatial resolution possible with optical methods offers the ability to provide greater displacement sensitivity compared to other elastography techniques based on ultrasound and magnetic resonance elastography. 80][11][12][13][14][15][16] OCE measurements have been used to assess the biomechanical properties of a wide variety of tissues of interest, such as measuring the stiffness of arteries, 17 different skin layers, 18 tissue engineering constructs, 19 and the cornea. 20,21CE requires the application of a mechanical load (stress) on the sample, and by measuring the resulting displacements, the viscoelastic properties of the sample can be extracted.Mechanical stresses are induced in the samples of interest through either dynamic or static methods.In dynamic methods, the excitation source could be external, such as mechanical compression using piezoelectric (PZT) actuators, [22][23][24] or it could be an internal excitation, such as by using acoustic radiation force to displace the internal tissue structures. 25,26arious methods could be used to quantify the displacement of the tissues.Earlier methods for elastography using optical coherence tomography (OCT) relied on speckle-tracking algorithms to estimate tissue displacement, however, speckle de-correlation and spatial sampling limit the range of displacements that can be reliably measured. 27Phase-resolved measurements, where the axial displacements are tracked based on the phase difference between adjacent A-scans, have been more widely employed to estimate small scale displacements. 16The lower and upper limits on the displacement range using phase-resolved methods depend on the phase noise and phase wrapping, respectively.Some research groups have also used shear wave propagation 25 or surface acoustic waves 28 in tissue samples to measure the stiffness of the tissue.
Externally applied stresses using a load cell or PZT actuators tend to be larger in scale and coverage, which may limit the high-resolution mapping of the mechanical properties.We chose to employ magnetic nanoparticles (MNPs) in this study for dynamic mechanical actuation as these could provide localized mechanical perturbations enabling the high-resolution probing of the tissue viscoelastic properties.Their small size allows them to diffuse through the tissues of interest, thereby allowing these to be used as local probes to assess the local mechanical environment of the tissues.It has been previously shown that they can enhance the contrast in OCT images 29 and can be effective transducers for sensitive magnetomotive OCE (MM-OCE) measurements. 30By measuring the response of the tissue to a step excitation at high spatial and temporal resolution, we propose a new way of simultaneously measuring several biomechanical parameters of interest, such as the elastic modulus, tissue displacement rise time, creep factor, and the natural frequency of biological tissues.

Tissue Specimen Preparation
Under a protocol approved by the Institutional Animal Care and Use Committee at the University of Illinois at Urbana-Champaign, tissues were harvested from two New Zealand white rabbits (Covance, Princeton, New Jersey) and sectioned in cubes of approximately 1.5 × 1.5 × 1.5 cm.The fresh tissues were then completely immersed for 4 h in a room-temperature saline suspension of MNPs (#637106, average diameter 25 nm, Sigma-Aldrich, Inc., St. Louis, Missouri) with a concentration of 10 mg∕ml.Prior to imaging, the tissues were rinsed with fresh saline solution and set at rest for 5 min in order to avoid fluid flow transients in the specimen during imaging.

Imaging System
The MM-OCE imaging system has been previously described in detail. 30Figure 1 shows a schematic of the optical system that includes the solenoid coil.In brief, a titanium:sapphire laser (KMLabs, Boulder, Colorado) with a center wavelength of 800 nm and a bandwidth of 120 nm provided the probing light for a spectral-domain OCT system.The average power incident on the specimens was 10 mW.The axial resolution was 1.7 μm and the transverse resolution was 16 μm in the specimens, assuming an average tissue index of refraction of 1.4.A small, custom-built, computer-controlled solenoid coil situated 2 mm above the specimen provided the switchable magnetic field for triggering the motion of the MNPs.The strength of the magnetic field at the measurement locations was approximately 400 G.

Data Acquisition
Microrheology measurements were made from M-mode imaging data acquired at camera line rates of 5 and 20 kHz to allow high time-resolution displacement measurements of the initial dynamic response of the tissues as well as longer time measurements that would show the mechanical behavior of the tissue past the initial perturbation.Measurements were made from four different locations close to the center of each specimen.Each M-mode image consisted of 10,000 axial lines and captured four cycles of a constant step magnetic field. 30The phase content of the acquired signals was used to determine the displacements occurring in the tissue specimens, as previously demonstrated. 30The displacement sensitivity of the system was 11 nm, allowing us to detect very small displacements in real time and with high resolution.As we have shown previously, this approach to displacement measurements is preferable to amplitude measurements, whose displacement sensitivity is much lower on the order of 1 μm. 30

Data Analysis
The normalized displacement data was analyzed in MATLAB to extract the following biomechanical parameters of interest: the natural frequency of oscillation from the underdamped oscillations, the rise time, defined as the interval between the magnetic field being switched on and the moment the tissue reached the first peak of the oscillation, and the exponential creep term of the displacement following the initial oscillation.Figure 2 shows Fig. 1 Biological magnetomotive optical coherence elastography (MM-OCE) setup.A standard 800 nm spectral-domain OCT system is used with an electromagnetic coil placed in the sample arm.A magnetic field strength of approximately 400 G was used to perturb the MNPs inside the sample and changes in the magnitude and phase of the interference pattern were detected in synchrony with the A-line scan rate and activation of the magnetic field coil.
representative data for lung tissue and highlights the relevant features of the signals that are utilized to calculate the above parameters.The data shown in Fig. 2(a) was acquired with 20 kHz line scan rate and describes the fast response of the tissue while the 5 kHz data in Fig. 2(b) characterizes the slow response.In order to quantify the description of these regimes, we fit the 20 kHz data to the equation y ¼ A × sinð2πft þ θÞ × e −bt where y is the displacement, A is the amplitude of the oscillation, f is the natural frequency of oscillation, t is the time, θ is a phase factor, and b is the damping parameter.The creep parameter c is calculated by fitting the 5 kHz data to the equation y ¼ 1 − a × e ct where a is the amplitude factor.The creep parameter quantifies how fast the tissue approaches a certain deformation state.
Young's moduli of the tissue specimens were calculated based on a Voigt model. 31In brief, Young's modulus E is given by the equation E ¼ kðL∕SÞ, where k is the elastic constant of the tissue specimen, L is the height of the tissue, and S is the cross-sectional area of the specimen.k is calculated from the expression k ¼ ½ð2πfÞ 2 þ b 2 × m where m is the mass of the tissue specimen.

Results and Discussion
The tissues measured with MM-OCE contained MNPs that had diffused into the specimens while being immersed in a concentrated suspension of MNPs in saline.Similar to our previous study with silicone tissue phantoms, 30 the MNPs likely bind to the extracellular matrix of the tissue, or perhaps are internalized into some cells, and provide access to interrogate the biomechanical microenvironment of the tissue specimens when perturbed with an external magnetic field.We chose to utilize a magnetic field applied in the form of a step function that would be kept on for a period of time sufficiently long to allow the observation of the dynamic response of the tissue specimen immediately after the MNPs are set in motion.
In a preliminary set of experiments, several types of rabbit tissues (adipose, heart, kidney, muscle, lung, and liver) were harvested and MM-OCE was performed.The displacement data for this set of tissue specimens is shown in Fig. 3.It is evident that each type of tissue has a distinctive response.We chose to study the lung and muscle tissues in more depth because they illustrate distinctly different types of responses, appeared more uniform than the other types of tissue, and could be handled more easily, so that each could be sectioned into simple geometric shapes such as cubes.This is important because the boundary conditions can influence the measurements of the natural frequency of oscillation, based on which Young's moduli are calculated. 13,31n another set of experiments, we performed MM-OCE measurements on tissues acquired from a different rabbit.Magnetomotive OCT images [Figs.4(a   (normalized displacements) obtained from the lung and muscle specimens are shown in Figs.4(c) and 4(d).M-mode MM-OCE data was taken at four nearby locations in each specimen and the displacements from bright scatterers along with the depth of the axial scans were calculated from the phase data. 30Depending on the tissue specimen and location, the number of bright scatterers varied from 6 to 20.The displacement data from these scatterers was averaged and utilized for postprocessing.Figures 4(e) and 4(f) show the noise-level signals from the same specimens when no magnetic field was present.Similar noise-level signals were observed in tissue specimens without MNPs, both with and without the switched magnetic field (data not shown).We note that the noise fluctuations of the displacement are extremely small compared to the signals we record when the magnetic field is switched from off to on or vice versa, with a ratio of maximum displacement to noise level displacement of at least 100∶1.
Table 1 summarizes the averaged values of the natural frequency of oscillation from the underdamped oscillations, the rise time, the exponential creep term of the displacement following the initial oscillation, and Young's modulus.The natural frequencies of oscillation are comparable to values measured previously in rat lung and muscle. 32Young's modulus values are also consistent with the previous measurements, though we note that there is an acknowledged wide range of values in the literature, due to biological variability, physical constraints, chemical environment, and the geometry of the samples. 14,32n our previous study with silicone tissue phantoms, the displacements we observed contained an underdamped oscillation followed by an equilibrium position. 30Several cycles of the oscillation were present for these highly homogeneous samples.The more heterogeneous rabbit lung and muscle tissues exhibit a similar trend in that an underdamped oscillation is also present in their response.However, the damping in these tissue specimens is much faster than in tissue phantoms as only one (in the case of lung) or less (in the case of muscle) cycle of oscillation is observed.Moreover, the tissues undergo  creep, as opposed to rapidly settling at an equilibrium position.
Creep occurs when the fluid inside the tissue redistributes in response to the tissue deformation due to a constant force.
Creep tests are common in the study of mechanical properties of different materials, including soft biological tissue. 33It is apparent that muscle experiences a steeper creep than lung, while the amplitude of its underdamped oscillation is relatively smaller.Also, when the magnetic field is turned off and the MNPs (and thus the forces on the tissue) are released, the tissue does not revert to its initial position, and some residual strain is present, more so in muscle than in lung.
In conclusion, we have demonstrated the application of MM-OCE in biological tissues, and how phase-sensitive measurements of the tissue response from an applied external step-function magnetic field can yield several important biomechanical properties of the tissue under investigation.These differences observed demonstrate that the behaviors of different types of tissues carry a unique set of information that results from a multitude of factors, such as the physio-chemical constituents of the tissue, the micro-level bonds in the extracellular matrix, the types of cells that make up the tissue, etc. Future studies will begin to elucidate the dependency of each of these factors in the extracted MM-OCE biomechanical data.As a novel optical elastography technique, MM-OCE can potentially provide access to this biomechanical information with high spatial and temporal resolution.This method can be extended to map depth-resolved two-dimensional and three-dimensional tissue biomechanical properties and can be a promising technique for advancing the study of biomechanics on the micro-scale.
) and 4(b)] of rabbit lung and muscle tissues acquired from the center of the specimens show the presence of the MNPs, as demonstrated by the green channel signal.Representative MM-OCE signals

Fig. 2
Fig. 2 Representative lung tissue response.(a) The displacement data acquired with a 20 kHz camera line rate contains a fast underdamped oscillation from which the rise time is calculated.The natural frequency of the tissue can be calculated by fitting the data after the creep component is subtracted.(b) The displacement data acquired with a 5 kHz camera line rate is fitted to determine the creep parameter, which describes how fast the tissue displaces over a large time scale compared to the initial oscillation.

Fig. 3
Fig. 3 Representative MM-OCE signals for different types of rabbit tissues acquired at a line scan rate of 20 kHz.The first transition corresponds to switching the B field from off to on, while the second one corresponds to switching the field from on to off.

Fig. 4
Fig. 4 MM-OCT and MM-OCE of biological tissues.(a and b) B-mode MM-OCT images of rabbit lung and muscle tissues where the green channel represents the MM-OCT signal and the red channel is the structural OCT signal.(c and d) Representative normalized displacements calculated from Mmode phase MM-OCE data taken with a line rate of 20 kHz from rabbit lung and muscle while the magnetic field B is modulated in a step-wise manner (square wave) for four cycles (B off followed by B on). (e and f) Normalized displacement data of lung and muscle in the absence of magnetic fields.

Table 1
Biomechanical parameters of rabbit lung and muscle tissues measured with magnetomotive optical coherence elastography.