|
The cornea is a transparent connective tissue that provides approximately two-thirds of the optical refracting power of the eye. Its shape and structure are critical for normal vision. Alterations in the biomechanical characteristics of the cornea can have a profound influence on its health, structural integrity, and normal function, and can lead to significant visual losses and even to permanent blindness.1 Therefore, it is critically important to understand the interplay between the mechanical properties of the cornea and its physiological function. Quantitative characterization of the biomechanical properties of the cornea could provide useful diagnostic information for early detection of corneal pathologies and for planning refractive surgery.2,3 Several studies have demonstrated that significant changes in the corneal microstructure take place with age, including increasing fibril diameter4 and stiffness.5 Most approaches to measure the biomechanical properties of the biological tissue have relied on inducing a stimulus and measuring the tissue’s response to that stimulus. The stimulus can be induced in many ways, e.g., exciting with mechanical force,6,7 acoustic radiation force,8 or by using a laser pulse.9 Combining different stimuli and methods to measure the mechanical response to these stimuli, several elasticity imaging modalities have been proposed such as magnetic resonance elastography,7,10–12 ultrasound elastography,13,14 acoustic radiation force imaging,15,16 and supersonic shear imaging.8,17,18 However, all of these methods require significant amplitude of tissue stimulation in order to produce a measurable signal. The only clinically accepted instrument for measurement of corneal biomechanics is ocular response analyzer (ORA), which uses a controlled air pulse to induce mechanical deformation of the eye.19 The ORA utilizes an electro-optical system to record applanation (flattening) of the corneal surface during both inward and outward deflection in response to an air pulse. Biomechanical properties are then estimated from the time taken for the corneal applanations to occur.19,20 Therefore, ORA also requires a large displacement of the corneal surface for evaluating biomechanical properties. However, the predictability of this system is still under investigation.21 Our method utilizes optical coherence tomography (OCT), a noninvasive in-depth imaging technology, to detect the vibrations on anterior and posterior corneal surfaces. OCT was first applied in 1998 by Schmitt to measure microscopic deformations due to compressive stress.22 Since then, several groups have been utilizing OCT to obtain the biomechanical properties of tissues.9,23–26 De la Torre et al. have measured displacement fields with a sensitivity of 10 nm in porcine cornea ex vivo.27 Corneal biomechanics using OCT have been demonstrated by Li et al. and Ford et al. in ex vivo28 and in situ29 conditions, respectively. Recently, Alonso-Caneiro et al. measured corneal dynamics in vivo by observing the corneal response to an air pulse from a series of depth-wise images generated from swept source OCT.23 However, they had to displace the cornea to the order of 1 mm in order to perform measurements. No work related to estimation of corneal biomechanics in vivo with minute corneal displacement has been reported yet, mainly due to the difficulty of inducing vibrations and detecting small amplitudes simultaneously in life conditions. In this work we utilize a previously developed method to detect vibrations on corneal surface in vivo with high sensitivity.30,31 The method measures the parameters required for estimating biomechanical properties of the cornea from a very low amplitude ( a few microns) excitation. To the best of our knowledge, this is the first time the wave parameters have been quantified in in vivo mice cornea using very low amplitude vibrations. In proof-of-principle studies, we utilized mechanical stimulus for inducing waves in the sample and phase-stabilized OCT for sensing the tissue response (Fig. 1). The details of the phase-stabilized swept source optical coherence elastography (PhS-SSOCE) and the method of obtaining phase response of the sample can be found in our previous publications.32–34 Briefly, the PhS-SSOCE consists of an electrically driven excitation unit and PhS-swept source optical coherence tomography (SSOCT) for sensing the vibrations. The excitation unit consists of a thin wire attached to a speaker diaphragm that responds to a rectangular pulse generated by a pulse generator (duration of 0.1 ms). The tip of the wire was rounded to remove sharp edges (contact area of ) to insure no damage was introduced to the tissues. With the onset of the adjustable input pulse, the tip of the wire comes in contact with the sample and introduces elastic flexural waves in the cornea. The PhS-SSOCT consists of a Mach-Zehnder interferometer for imaging and a fiber Bragg grating (FBG) for triggering. The phase-stability of the system is measured as 0.09 rad that corresponds to 1.9 nm of displacement sensitivity. We chose PhS-SSOCT for these studies due to its ability to extract phase information from extended tissue depths. PhS-SSOCE readings were performed at different spatial locations as shown in Fig. 1(c). From these measurements of the phase responses, the wave attenuation and wave velocity were quantified. Wave attenuation was directly obtained from the amplitude of the phase response, whereas the wave velocity was measured from the delay observed in the phase response between two points of known distances. Therefore, as shown in the insert c of Fig. 1, phase responses were measured at increasing 0.5 mm incremental distances away from the location where mechanical stimulus was applied. A typical temporal response of the corneal surface to the mechanical excitation measured at an arbitrary spatial location is shown in Fig. 2(a). The corneal surface displacements measured at increasing distances relative to the excitation pulse are shown in Fig. 2(b). These data clearly depict that the amplitude of the wave decreases with an increase in the distance from the point of stimulus induction. These amplitudes were used to plot a wave propagation tomogram from both anterior and posterior surfaces of mice cornea as shown in Fig. 3(a) and 3(b), respectively, where the attenuation of the wave amplitude with the increasing distances is clearly seen. In these experiments, amplitudes of the mechanical waves recorded at different spatial locations were from approximately 10 μm to 2 μm, highest near the excitation point and gradually decreasing as they propagate in the cornea. Data shown in Fig. 2(b) also demonstrate that the measured amplitude is delayed for the points located farther away from the excitation point. This phase delay can be used to estimate the time taken by the wave to travel over that distance. Therefore, by taking the ratio between the distance and the time delay, velocity of the elastic wave can be calculated. As a proof-of-concept, these velocities were calculated from the anterior surface of the mouse cornea for 1, 6, and 13-month-old animals. In these pilot studies, only one mouse has been used for each age group; however, the measurements were made three times at different spatial locations of the excitation/recording. The velocity as a function of age is shown in Fig. 4 and demonstrates clear increase of the velocities with the animal age. Our previous phantom studies showed that higher velocities correspond to stiffer phantoms.30 Additionally, recent studies demonstrate clear dependence of ocular tissue elasticity as a function of age.35 Therefore, higher velocities measured in corneas of older mice suggest that the elasticity of the cornea is increasing with the age. The results shown in this letter demonstrate the capability of the PhS-SSOCE technique to measure tissue response to very minute stimulus. This method can be directly applied to any kind of stimulus to measure the wave parameters. Furthermore, using the theory of elasticity and reconstructive approaches, the wave parameters can be used to evaluate the biomechanical properties of tissue in vivo and completely noninvasively. AcknowledgmentsThis work was supported by NIH Grant 1R01EY022362-01 and, in part, by the Federal Target Program “Scientific and Scientific-Pedagogical Personnel of Innovative Russia” for 2009–2013 Grant 14.B37.21.1238. ReferencesM. C. ChiramboD. Benezra,
“Causes of blindness among students in blind school institutions in a developing country,”
Br. J. Ophthalmol., 60
(9), 665
–668
(1976). http://dx.doi.org/10.1136/bjo.60.9.665 BJOPAL 0007-1161 Google Scholar
I. F. ComaishM. A. Lawless,
“Progressive post-LASIK keratectasia: biomechanical instability or chronic disease process?,”
J. Cataract Refractive Surg., 28
(12), 2206
–2213
(2002). http://dx.doi.org/10.1016/S0886-3350(02)01698-X 0886-3350 Google Scholar
L. T. Nordan,
“Keratoconus: diagnosis and treatment,”
Int. Ophthalmol. Clin., 37
(1), 51
–63
(1997). http://dx.doi.org/10.1097/00004397-199703710-00005 IOPCAV 0020-8167 Google Scholar
A. Daxeret al.,
“Collagen fibrils in the human corneal stroma: structure and aging,”
Invest. Ophthalmol. Visual Sci., 39
(3), 644
–648
(1998). IOVSDA 0146-0404 Google Scholar
A. Elsheikhet al.,
“Assessment of corneal biomechanical properties and their variation with age,”
Curr. Eye Res., 32
(1), 11
–19
(2007). http://dx.doi.org/10.1080/02713680601077145 CEYRDM 0271-3683 Google Scholar
K. W. Hollmanet al.,
“Strain imaging of corneal tissue with an ultrasound elasticity microscope,”
Cornea, 21
(1), 68
–73
(2002). http://dx.doi.org/10.1097/00003226-200201000-00015 CORNDB 0277-3740 Google Scholar
D. V. Litwilleret al.,
“MR elastography of the ex vivo bovine globe,”
J. Magn. Reson. Imaging, 32
(1), 44
–51
(2010). http://dx.doi.org/10.1002/jmri.v32:1 1053-1807 Google Scholar
M. Tanteret al.,
“High-resolution quantitative imaging of cornea elasticity using supersonic shear imaging,”
IEEE Trans. Med. Imag., 28
(12), 1881
–1893
(2009). http://dx.doi.org/10.1109/TMI.2009.2021471 ITMID4 0278-0062 Google Scholar
S. Niet al.,
“Effect of glucose on the stress-strain behavior of ex-vivo rabbit cornea,”
Exper. Eye Res., 92
(5), 353
–360
(2011). http://dx.doi.org/10.1016/j.exer.2011.02.005 EXERA6 0014-4835 Google Scholar
A. Manducaet al.,
“Magnetic resonance elastography: non-invasive mapping of tissue elasticity,”
Med. Image Anal., 5
(4), 237
–254
(2001). http://dx.doi.org/10.1016/S1361-8415(00)00039-6 1361-8415 Google Scholar
P. J. McCrackenet al.,
“Mechanical transient-based magnetic resonance elastography,”
Magn. Reson. Med., 53
(3), 628
–639
(2005). http://dx.doi.org/10.1002/(ISSN)1522-2594 MRMEEN 0740-3194 Google Scholar
R. Muthupillaiet al.,
“Magnetic resonance elastography by direct visualization of propagating acoustic strain waves,”
Science, 269
(5232), 1854
–1857
(1995). http://dx.doi.org/10.1126/science.7569924 SCIEAS 0036-8075 Google Scholar
J. Ophiret al.,
“Elastography: ultrasonic estimation and imaging of the elastic properties of tissues,”
Proc. Institution of Mechanical Engineers—Part H—J. Eng. Med., 213
(3), 203
–233
(1999). http://dx.doi.org/10.1243/0954411991534933 0954-4119 Google Scholar
K. J. ParkerR. M. Lerner,
“Sonoelasticity of organs: shear waves ring a bell,”
J. Ultrasound Med., 11
(8), 387
–392
(1992). JUMEDA 0278-4297 Google Scholar
M. D’OnofrioA. GallottiR. P. Mucelli,
“Tissue quantification with acoustic radiation force impulse imaging: measurement repeatability and normal values in the healthy liver,”
Am. J. Roentgenol., 195
(1), 132
–136
(2010). http://dx.doi.org/10.2214/AJR.09.3923 AJROAM 0092-5381 Google Scholar
K. Nightingale,
“Acoustic radiation force impulse (ARFI) imaging: a review,”
Curr. Med. Imaging Rev., 7
(4), 328
–339
(2011). http://dx.doi.org/10.2174/157340511798038657 1573-4056 Google Scholar
J. BercoffM. TanterM. Fink,
“Supersonic shear imaging: a new technique for soft tissue elasticity mapping,”
IEEE Trans. Ultrason. Ferroelectr. Freq. Control, 51
(4), 396
–409
(2004). http://dx.doi.org/10.1109/TUFFC.2004.1295425 Google Scholar
M. Tanteret al.,
“Quantitative assessment of breast lesion viscoelasticity: initial clinical results using supersonic shear imaging,”
Ultrasound Med. Biol., 34
(9), 1373
–1386
(2008). http://dx.doi.org/10.1016/j.ultrasmedbio.2008.02.002 USMBA3 0301-5629 Google Scholar
D. A. Luce,
“Determining in vivo biomechanical properties of the cornea with an ocular response analyzer,”
J. Cataract Refractive Surg., 31
(1), 156
–162
(2005). http://dx.doi.org/10.1016/j.jcrs.2004.10.044 0886-3350 Google Scholar
G. J. OrssengoD. C. Pye,
“Determination of the true intraocular pressure and modulus of elasticity of the human cornea in vivo,”
Bull. Math. Biol., 61
(3), 551
–572
(1999). http://dx.doi.org/10.1006/bulm.1999.0102 BMTBAP 0092-8240 Google Scholar
A. Saadet al.,
“Biomechanical properties of keratoconus suspect eyes,”
Invest. Ophthalmol. Visual Sci., 51
(6), 2912
–2916
(2010). http://dx.doi.org/10.1167/iovs.09-4304 IOVSDA 0146-0404 Google Scholar
J. Schmitt,
“OCT elastography: imaging microscopic deformation andstrain of tissue,”
Opt. Express, 3
(6), 199
–211
(1998). http://dx.doi.org/10.1364/OE.3.000199 OPEXFF 1094-4087 Google Scholar
D. Alonso-Caneiroet al.,
“Assessment of corneal dynamics with high-speed swept source Optical Coherence Tomography combined with an air puff system,”
Opt. Express, 19
(15), 14188
–14199
(2011). http://dx.doi.org/10.1364/OE.19.014188 OPEXFF 1094-4087 Google Scholar
B. F. Kennedyet al.,
“In vivo dynamic optical coherence elastography using a ring actuator,”
Opt. Express, 17
(24), 21762
–21772
(2009). http://dx.doi.org/10.1364/OE.17.021762 OPEXFF 1094-4087 Google Scholar
X. Lianget al.,
“Dynamic spectral-domain optical coherence elastography for tissue characterization,”
Opt. Express, 18
(13), 14183
–14190
(2010). http://dx.doi.org/10.1364/OE.18.014183 OPEXFF 1094-4087 Google Scholar
X. Lianget al.,
“Optical micro-scale mapping of dynamic biomechanical tissue properties,”
Opt. Express, 16
(15), 11052
–11065
(2008). http://dx.doi.org/10.1364/OE.16.011052 OPEXFF 1094-4087 Google Scholar
M. H. De la Torre-IbarraP. D. RuizJ. M. Huntley,
“Double-shot depth-resolved displacement field measurement using phase-contrast spectral optical coherence tomography,”
Opt. Express, 14
(21), 9643
–9656
(2006). http://dx.doi.org/10.1364/OE.14.009643 OPEXFF 1094-4087 Google Scholar
C. Liet al.,
“Noncontact all-optical measurement of corneal elasticity,”
Opt. Lett., 37
(10), 1625
–1627
(2012). http://dx.doi.org/10.1364/OL.37.001625 OPLEDP 0146-9592 Google Scholar
M. R. Fordet al.,
“Method for optical coherence elastography of the cornea,”
J. Biomed. Opt., 16
(1), 016005
–016007
(2011). http://dx.doi.org/10.1117/1.3526701 JBOPFO 1083-3668 Google Scholar
R. K. Manapuramet al.,
“Estimation of shear wave velocity in gelatin phantoms utilizing PhS-SSOCT,”
Laser Phys., 22
(9), 1439
–1444
(2012). http://dx.doi.org/10.1134/S1054660X12090101 LAPHEJ 1054-660X Google Scholar
R. K. Manapuramet al.,
“Assessment of wave propagation on surfaces of crystalline lens with phase sensitive optical coherence tomography,”
Laser Phys. Lett., 8
(2), 164
–168
(2011). http://dx.doi.org/10.1002/lapl.201010102 1612-2011 Google Scholar
R. ManapuramV. ManneK. Larin,
“Development of phase-stabilized swept-source OCT for the ultrasensitive quantification of microbubbles,”
Laser Phys., 18
(9), 1080
–1086
(2008). http://dx.doi.org/10.1134/S1054660X08090144 LAPHEJ 1054-660X Google Scholar
R. K. ManapuramV. G. R. ManneK. V. Larin,
“Phase-sensitive swept source optical coherence tomography for imaging and quantifying of microbubbles in clear and scattering media,”
J. Appl. Phys., 105
(10), 102040
(2009). http://dx.doi.org/10.1063/1.3116614 JAPIAU 0021-8979 Google Scholar
R. K. Manapuramet al.,
“3D assessment of mechanical wave propagation in the crystalline eye lens using PhS-SSOCT,”
Proc. SPIE, 7885 78851V
(2011). http://dx.doi.org/10.1117/12.875064 PSISDG 0277-786X Google Scholar
N. E. K. CartwrightJ. R. TyrerJ. Marshall,
“Age-related differences in the elasticity of the human cornea,”
Invest. Ophthalmol. Visual Sci., 52
(7), 4324
–4329
(2011). http://dx.doi.org/10.1167/iovs.09-4798 IOVSDA 0146-0404 Google Scholar
|