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. Elastography methods have been developed using ultrasound, magnetic resonance, and X-ray computed tomography, among others.220.127.116.11.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.8 Optical coherence elastography (OCE) has been recognized as a versatile approach for measuring and mapping biomechanical properties of biological tissues on the micron scale.910.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,21 OCE 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,2223.–24 or it could be an internal excitation, such as by using acoustic radiation force to displace the internal tissue structures.25,26
Various 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.27 Phase-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.16 The 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 propagation25 or surface acoustic waves28 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 images29 and can be effective transducers for sensitive magnetomotive OCE (MM-OCE) measurements.30 By 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 . 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 . 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.
The MM-OCE imaging system has been previously described in detail.30 Figure 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.
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.30 The phase content of the acquired signals was used to determine the displacements occurring in the tissue specimens, as previously demonstrated.30 The 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
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 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 where is the displacement, is the amplitude of the oscillation, is the natural frequency of oscillation, is the time, is a phase factor, and is the damping parameter. The creep parameter is calculated by fitting the 5 kHz data to the equation where 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.31 In brief, Young’s modulus is given by the equation , where is the elastic constant of the tissue specimen, is the height of the tissue, and is the cross-sectional area of the specimen. is calculated from the expression where 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,31
In another set of experiments, we performed MM-OCE measurements on tissues acquired from a different rabbit. Magnetomotive OCT images [Figs. 4(a) 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 (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.30 Depending 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 .
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.32 Young’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,32
Biomechanical parameters of rabbit lung and muscle tissues measured with magnetomotive optical coherence elastography.
|Creep parameter (1/s)||−8.6±1.9||−6.0±2.4|
|Rise time (ms)||7.1±0.7||11.3±0.5|
|Young’s modulus (kPa)||94.6±15.2||28.5±2.8|
In our previous study with silicone tissue phantoms, the displacements we observed contained an underdamped oscillation followed by an equilibrium position.30 Several 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.33 It 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.
This work was supported in part by a grant from the National Institutes of Health (R01 EB009073, S.A.B.). We thank Eric Chaney and Dr. Marina Marjanovic from the Beckman Institute for Advanced Science and Technology for their help with harvesting the tissue specimens. We also thank Dr. Steven Adie and Yue Wang for helpful discussions on the theory of elasticity. Additional information can be found at http://biophotonics.illinois.edu.