Quantitative evaluation of degenerated tendon model using combined optical coherence elastography and acoustic radiation force method

Abstract. Damage of collagen fibers in tendons is often directly related to changes in a tendon’s mechanical properties. Direct quantitative elasticity measurement of tendons will provide important information in tendon dysfunction diagnosis and treatment assessment. A feasibility study of quantifying the mechanical properties of a degenerated tendon model by a nondestructive imaging modality, which combines optical coherence elastography and acoustic radiation force (ARF) method, is presented. The degenerated tendon model was produced by the partial degradation of chicken tendons through incubation with collagenase at different concentrations and incubation times. A 30-kHz longitudinal ultrasound transducer was used to provide an ARF signal, which was detected by an ultra-high sensitive phase sensitive optical coherence tomography (PhS-OCT) system. The experimental results demonstrate that the combination of ARF method and PhS-OCT can measure the elasticity of tendon quantitatively. The corresponding changes in tendon elasticity due to the application of collagenase have been revealed by this new imaging modality. This method can potentially be used in the assessment of tissue engineering products and in the diagnosis and treatment progression of tendon diseases.


Introduction
Tendons are composed of uniaxial orientated collagen fibril bundles in hydrated elastin and proteoglycans matrix. [1][2][3][4][5] Its role is to reduce the impedance mismatch of the different tissues when transmitting contractile forces from the muscle to the skeleton 1,2,5 at the osseotendenous junction. Physical trauma in tendons due to high rate loading, laceration, or degeneration due to ageing will rupture/damage collagen fibers, which will cause changes in tendon mechanical behavior. 6,7 Thus, the measurement of the mechanical properties of tendons is of great importance, e.g., in predicting the severity of the injury to the tendons involved; 8 for diagnosis of tendon diseases; 9 in monitoring treatment procedures, e.g., tendon surgical repair; 10,11 or for quantitative evaluation of tendon construction in tendon tissue engineering. 12 A number of techniques have been developed for evaluation of the tendon mechanical properties, especially the Young's modulus of tendon tissue. 13,14 Extensive studies on isolated animal or human tendon material have been done under quasi-static loading conditions, which undergo elongation to failure. [13][14][15] However, these methods are destructive. Recently, non-destructive imaging techniques became increasingly popular to investigate muscle and tendon function, including magnetic resonance imaging (MRI) and real time ultrasonography and the corresponding elastography technique using strain or shear wave methods. 16,17 These modalities offer in vivo measurement of muscle-tendon morphology and mechanical properties. However, the spatial resolutions of these techniques are limited by the resolution of ultrasound or MRI, which is a scale of tens to hundreds of micron scale. This limits their application in the musculoskeletal system and for characterization of early-phase tendon disease diagnosis and treatment.
Optical coherence elastography (OCE) [18][19][20][21][22][23][24][25][26] has several advantages over other imaging modalities. It is a noninvasive imaging technique with ultra-high resolution, which is at least an order higher than ultrasound or MRI elastography. It is capable of providing ultra-high microstructural and mechanical characterization of human soft tissues, and therefore has promising clinical applications. This technique monitors the passive tissue responses of the samples stimulated by speckle [19][20][21][22][23] and the Doppler signals. [24][25][26] Phase sensitive optical coherence tomography (PhS-OCT) combined with shear or surface acoustic wave (SAW) methods have become increasingly popular in tissue characterization. Liang and Boppart used PhS-OCT to measure the shaker induced continuous surface wave in vivo on human skin. This method was able to quantify the bulk elastic properties of skin tissues. 27 Li et al. [28][29][30] introduced SAW by external shaker or Q-switch laser pulse on sample surface. Using this method, bulk quantitative elasticity information of different layers of phantom and in vivo human skin has been estimated. Kennedy et al. 31 applied OCE for in vivo characterization of normal and hydrated human skins, relying on the dynamic OCT system in combination with a designed ring actuator, which induces longitudinal vibration to samples. More recently, this technique was developed in providing the three-dimensional OCE images of human skin. 32 However, they do not provide a quantitative measurement of elasticity in these studies. Furthermore, there has been no report for measurements on tendon tissues using the OCT-based elastography modality.
This paper thus explores the use of OCE to provide localized and quantitative mechanical properties of tendons. To achieve this goal, we applied a modulated acoustic radiation force (ARF) to chicken tendons ex vivo to induce phase change caused by the vibration. The tendon samples are treated with collagenase to mimic degenerated tissues and then imbedded into 2% agar-agar phantom. A PhS-OCT was utilized to record the phase change of the agar phantom as well as the ex vivo tendon, from which the quantitative Young's modulus was deduced.

System Configuration
The system consists of both the generation of longitudinal waves (induced by an in-house designed and fabricated 30-kHz ultrasound transducer) and detection by PhS-OCT for soft tissue. The setup is schematically shown in Fig. 1.
In order to generate acoustic signals for the elasticity measurement, a transducer with resonant frequency of 30 kHz, which was driven by a square-wave modulated radiation force (RF) signal (50% duty cycle amplifier modulation), was used. The modulation frequency was set to 30 kHz, and the maximum actuator displacement applied to the sample was ∼50 nm to ensure the generated strain was in the pure linear-elastic regime. In all experiments, samples were directly placed on the ultrasound transducer. The ARF could then be transmitted from the transducer to compress the samples and trigger vibration in the axial direction.
Detection of the vibration signal was performed using a PhS-OCT system, which has been described in detail previously. 31,32 Briefly, the PhS-OCT system employed a superluminescent diode as a light source, with a center wavelength of ∼1310 nm and a bandwidth of ∼83 nm, implemented by a spectral domain configuration. The sample arm used an objective lens of ∼50-mm focal length to deliver the detection light on the sample that coupled the vibration signals into the PhS-OCT for detection. The system provided an axial resolution of 8.9 μm and a transverse resolution of 18 μm in air. The acquisition rate was determined by the spectrometer employed in the system that had a maximum rate of ∼92;000 A-scans∕s. For signal detection, the system was configured as M-mode acquisition while the OCT probe beam stayed constantly in one location of the sample. The system provided a phase sensitivity to the sample displacement of ∼200 pm when working in M-scan mode. Two thousand and forty-eight A-scans were acquired to obtain one frame of M-mode scan which takes about 22.3 ms. The structural image of the samples generated by the PhS-OCT system is shown as a function of depth.

Sample Preparation
To prove the accuracy of the ARF method, three double layer phantoms with ∼2% agar over ∼1% agar, ∼3% agar over ∼2% agar and ∼3% agar over ∼1% agar and with ∼1 mm of upper layer thickness were created for testing.
To investigate the effectiveness of the OCE system in the detection of degenerative changes in tendinopathy, collagenase from Clostridium histolyticum (type collagen, Sigma, China) was used to degrade the collagen in the tendon to form a degenerated tendon model similar to pathological diseases. The collagenase powder was dissolved in deionised water to different concentrations. Tendons dissected from fresh nonfrozen chicken (bought from supermarket) were treated with different collagenase concentrations (1 mg∕mL, 3.3 mg∕mL, and 10 mg∕mL) and with different incubation times (0, 30, 60, and 90 min) in the incubator at 37°C and 5% CO 2 . The collagenase was removed from the tendons in order to stop further digestion after the defined culture time.
After the treatments, tendons were imbedded in 2% agar-agar with ∼0.5 to 1 mm thickness from the agar phantom surface. Agar is known to have very low scattering properties; thus, milk was added to each phantom at a concentration of 1% to produce a scattering background similar to what is found in biological tissues. Agar also keeps the tendon hydrated. For each concentration and incubation time span, five samples were prepared for experiment.

Data Acquisition and Processing
The application of ARF generates vibration within the sample, which can be observed and detected as phase shift on PhS-OCT. We calculated the phase difference between the first and all the other A-lines (ϕ i − ϕ 1 ), where i is a value between 1 and 2048 and ϕ is the phase angle of the A-line at a specific depth. It has been observed that the phase changed at each depth in Fig. 2(a) as dark and light bands through time. These bands represent a frequency response at 30 kHz, as shown in Fig 2(c), which is the modulation of the ultrasound transducer. Figure 2(b) shows the phase differences, i.e., vibration through time at a depth of 700 μm from the sample surface [the depth location was identified with an arrow in Fig. 2(a)]. Vibration amplitude, which is directly related to the strain, was highly sensitive to the change of material elasticity. Under the same stress, softer material has higher strain and stiffer material has smaller strain. With respect to the Hooke's law of elasticity, Young's modulus is the ratio of stress to strain within the elastic region, which can be expressed as where E is the Young's modulus, σ is stress, and ε is strain. In the application of OCE to a two-layered sample, constant stress fields can be assumed to apply to both layers. Thus, under the same stress (σ 1 ≈ σ 2 ), in a heterogeneous material, the ratio of Young's modulus in different layers 1 and 2 is equal to the inverse ratio of the strain Here, the local strain is determined from the detected tissue displacement caused by the vibration amplitude: 31 where d 1 is the displacement amplitude estimated at a depth of z from the top of the sample and d 2 is the displacement amplitude estimate at a depth of z þ Δz. Figure 2(c) presents the results of fast Fourier transform (FFT) of the phase difference [FFT (Δϕ)] and therefore, the amplitude from the data shown in Fig 2(b). A peak at 30 kHz is clearly observed. The amplitude of this peak at each depth, after subtracting the background noise and averaging all the A-lines, is used to create the phase FFT amplitude signal as a function of depth [ Fig. 2(d)]. The phase difference caused by the vibration and the phase FFT amplitude are linearly related. 33 Combining Eqs. (2) and (3), if we take the Δz to be long enough, the ratio between the slope of FFT amplitude from agar phantom and tendon gives the ratio of strain, i.e., the inverse ratio of Young's modulus.
A linear polynomial model has been chosen to fit the averaged phase FFT amplitude of both 2% agar layer and tendon layer in one sample [ Fig. 2(d) green line and blue line] with 95% confidence bounds. The beginning and ending point of the fitting is decided by the surface of the agar and the interface of the two layers, however, we need to sacrifice some fitting length due to the inaccuracy of data, which may be caused by some noise at the surface or sharp attenuation of backscattering light intensity at the deep part. Then, the slopes of the linear fit can be determined. Taking previous measurement values of the Young's modulus for 2% agar-agar phantom, 193 kPa, 29 and using the above mentioned inversion method, the quantitative Young's modulus of tendon can be evaluated. Figure 3 shows the amplitude of the FFT of the phase changes caused by the ARF at different depths in both ∼2% agar over ∼1% agar and ∼3% agar over ∼2% agar phantoms. From Eqs. (2) and (3), if we make the Δz long enough to cover the whole layer thickness, the ratio of slope of the FFT amplitude from different layers of the sample is equal to the ratio of shear wave strain. Thus, we apply a linear polynomial model fitted to the averaged phase FFT amplitude of different layers of phantoms (Fig. 3, blue line and red dash line) with 95% confidence bounds. From Figs. 3(a) and 3(b), we find that at higher agar concentration, the slope of FFT amplitude curve is smaller, as under the same stress the stiffer material will have a smaller strain.   29 In Fig. 3(a), the estimated slopes of FFT amplitude curves for 1% agar layer and 2% agar layer are −67.51 and −25.57, respectively. The slope ratio of 2% and 1% is 2.64, which approximates the inverse ratio of Young's modulus of the two concentrations (2. 41 ¼ 193 kPa∕80 kPa). Similarly in Fig. 3(b), the estimated slope of FFT amplitude curves for 2% agar layer and 3% agar layer is −122.19 and −47.36. The ratio is 2.58, which approximates the inverse ratio of Young's modulus of 3% and 2% agar (2.66 ¼ 515 kPa∕193 kPa). We then try the phantom of ∼3% agar over ∼1% agar, which shows a more significant difference in Young's modulus (results shown in Fig. 4). As expected, the slope ratio of 1% agar and 3% agar vibration amplitude is 6.64 (−109.6∕ − 16.5), which is similar to the inverse ration of Young's modulus of 3% and 1% agar (6.43 ¼ 515 kPa∕80 kPa).

Results and Discussion
From the above phantom experiments, we have proven that the strain ratio, which is equal to the ratio of FFT vibration amplitude curve slopes, is approximate to the inverse ratio of its corresponded Young's modulus. Knowing the ratio and Young's modulus of any one material, the quantitative elasticity of connected other materials can be calculated easily.
There is an amplitude increase at the boundary between stiffer material and softer material, as seen in Figs. 3(b) and Fig. 4. This is due to the fact that the shear wave travels from stiffer material to softer material and so the transmitted wave amplitude will be greater than the incident shear wave. The percentage is dependent on the materials' acoustic impedances, which is the result of material density and acoustic wave velocity: where A t ðA i Þ is the transmitted (incident) wave amplitude, and z, ρ, and c are the acoustics impedance, density, and the acoustic wave velocity of the material. A stiffer material has higher density and shear wave velocity, thus, phenomena of amplitude increase can be observed. Figure 5(a) shows the real-time cross-sectional microstructures image (B-scan) of control tendon in a 2% agar phantom captured by the PhS-OCT system. The OCT image clearly demarcates the tendon tissue at ∼0.8 mm depth below the phantom surface in this sample. It provides an indicator and crossvalidation of the tendon boundary which can be used for the linear fit of the agar and the tendon layers. The vertical dash line indicates the depth range that can be detected by the OCT system. At the depth of 2 mm, the system also stops the phase shift detection. Figure 5(b) shows FFT amplitude of the phase changes caused by the ARF-induced vibration at different depths. Straight lines show linear model fits of the 2% agar area and control tendon without treatment. The estimated FFT amplitude curve slopes for the agar phantom layer and the control tendon are −9.367e þ 004 and −1377, respectively. This gives a phantom-tendon slope ratio of 68. We assume that both the agar layer and tendon layer are homogeneous, then this value can be treated as the ratio of strain for  these two layers. According to Eq. (2), the Young's modulus of control tendon can be estimated as 13124 kPa (∼193 kPa × 68). Figure 6 shows typical averaged FFT amplitudes and the linear fits of the phase differences caused by ARF-induced vibration as the function of depth of tendon-agar samples that were treated with collagenase solution (10 mg∕mL and 30 min). As seen in Fig. 6(a), the slope difference between ∼2% agar phantom and tendon sample was greatly reduced compared to the control group in Fig. 5(b). This result indicates the Young's modulus of the tendon sample has been reduced, as was expected because collagen fibers have been damaged by collagenase. With the treatment time increased to 60 and 90 min, the averaged FFT amplitude slopes of tendons became nearly equal to [ Fig. 6(b)] and larger [ Fig. 6(c)] than that of the agar phantom, which indicates that the Young's modulus was equal and smaller than that of agar phantom with 60 and 90 min treatment, respectively. Figure 7 shows the overall Young's modulus changes of collagenase treated tendon samples with different incubation times and concentrations, and the quantitative Young's modulus is summarized in Table 1. As expected, the Young's modulus   of tendon tissue is decreased as a function of collagenase digestion. Higher concentration of collagenase and longer treatment time led to a larger Young's modulus decrease. Figure 8 shows the OCT structure images, OCE images, and the overlaid images of tendon samples treated with collagenase for 1, 3.3, and 10 mg∕mL for 60 min, respectively. The twodimensional strain image was calculated based on the method developed from our previous study. 33 Applying Eq. (2) and using the Young's modulus value of 193 kPa for the agar phantom, the strain image was converted to the OCT elastography image with absolute Young's modulus values, as shown in Figs. 8(d)-8(f), instead of relative strain value as shown in Kennedy et al. 31 It shows that using the presented elasticity reconstruction approach based on combined ARF and OCE methods, quantitative elasticity images of an agar-tendon sample can be reconstructed with high resolution. A clear  Journal of Biomedical Optics 111417-6 November 2013 • Vol. 18 (11) tendon-agar boundary can be found in the elasticity images although the signal decreases at deeper locations. This is expected because the accuracy of the phase measurement relies on the intensity level from structure image. 33

Conclusions
A novel concept of combining an ARF method and PhS-OCT to characterize and monitor the quantitative information of chicken tendons ex vivo is reported in this paper. Experiments show that the method can efficiently monitor the elasticity change of the tendon due to different collagen fiber degradation. PhS-OCT is also an ideal tool to detect the ARF vibration as it provides high signal to noise ratio (SNR) on tendon samples. It is able to provide the B-frame image of the sample, which may validate the thickness estimation in elastography, while providing the elastography of the tendon tissues. The Young's modulus measured from our method corresponded well with the predicted results. A benefit of PhS-OCT based elastography, in addition to being noninvasive and noncontact, is real-time simultaneous tomographic imaging of the tendon samples. The concept reported in this paper can be developed further into a truly quantitative elastographic imaging system, which has great application value in tissue engineering to evaluate the growth of artificial tendon tissue in real-time as nondisruptive feedback for bioreactors. 34 With further development of the PhS-OCT system, this method may have a great potential in real clinical diagnosis to assess the physiological conditions of human tendons. In future work, this method can be combined with the needle probe OCT. The insertion of an OCT probe into deep tissue will overcome the penetration limits of OCT light beam.