Correlation of IOP with Corneal Acoustic Impedance in Porcine Eye Model

Zhang, Jun; Zhang, Yi; Li, Yang; Chen, Ruimin; Shung, K. Kirk; Richter, Grace; Zhou, Qifa

doi:https://doi.org/10.1155/2017/2959717

BioMed Research International

On this page

Abstract Introduction Results Discussion Disclosure Conflicts of Interest Authors’ Contributions Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2017 | Article ID 2959717 | https://doi.org/10.1155/2017/2959717

Correlation of IOP with Corneal Acoustic Impedance in Porcine Eye Model

Jun Zhang,^1,2Yi Zhang,³Yang Li,²Ruimin Chen,²K. Kirk Shung,²Grace Richter,³and Qifa Zhou^2,3

Academic Editor: Sergio Claudio Saccà

Received23 Feb 2017

Revised24 Apr 2017

Accepted14 May 2017

Published31 May 2017

Abstract

Purpose. The aim of this study is to correlate the intraocular pressure (IOP) change with the acoustic impedance of the cornea, in order to propose a noncontact and noninvasive method for IOP monitoring. Methods and Materials. A highly focused transducer (frequency 47-MHz; bandwidth 62%) was made to measure the echo from the anterior and posterior surfaces of intact porcine eyes, respectively. A multilayered transmission and reflection model was used to calculate the acoustic impedance. The linear relationship between acoustic impedance and intraocular pressure was analyzed by statistical method. Result. During pressure elevation from 10 mm Hg to 50 mm Hg, the mean acoustic impedance of the posterior cornea increased from 1.5393 to 1.5698 MRayl, which showed a strong linear correlation (; ). Meanwhile, the mean value of the anterior cornea increased from 1.5399 to 1.5519 MRayl, and a less significant correlation was observed (; ). Conclusion. This study revealed a linear correlation between intraocular pressure and acoustic impedance of the cornea, thus demonstrating a potentially important method to noninvasively measure the intraocular pressure in vivo.

1. Introduction

Elevated intraocular pressure (IOP) is the main risk factor for glaucoma, which is a leading cause of irreversible optic nerve damage for blindness [1, 2]. Nowadays, several methods of IOP measurement have been developed for clinical use, including the gold-standard Goldmann applanation tonometer (GAT), the dynamic contour tonometer (DCT), Tono-Pen, and Icare rebound tonometer. A contact lens sensor (CLS) has also been developed for IOP monitoring [3–6], capable of providing data needed for the clinician to characterize certain types of glaucoma based on fluctuations in IOP for one 24-hour period. However, IOP measurements by some of the currently developed devices may be affected by ocular properties such as corneal curvature, hysteresis, and thickness [6, 7]. There is an unmet need for a noninvasive and accurate method for IOP measurement.

Previous works in ocular biomechanics have focused on studying the influence of intraocular pressure on the biomechanical properties of tissue such as cornea, sclera, and optic nerve. For example, strip extensiometry and inflation tests were always used for measurement of the cornea elastic modulus according to the nonlinear strain-stress curve [8]. However, due to different ocular conditions, there is inconsistency of the measured results in the prior studies [9–11]. Mathematical models were also established to improve the understanding of the complex biomechanical behavior of the eyeball. Anderson et al. used nonlinear finite-element modelling to study the behavior of the cornea under different loading states [12]. Elsheikh et al. used shell theory for pressure-deformation results analysis and found that the correlation between the modulus of elasticity in the cornea and intraocular pressure had a positive linear correlation, although the geometry deformation and stress change in the cornea were complex [13, 14]. This gives us a clue that we could use elasticity of the cornea or related parameters to characterize the IOP level.

Ultrasonic waves have been used to measure elastic properties of the cornea and sclera. Tanter et al. used supersonic shear waves to get three-dimensional quantitative maps of the corneal elasticity at two IOP levels (10 mmHg and 20 mmHg) [15–17]. Liu et al. used an ultrasound speckle tracking method to map 3D strains of the porcine sclera inflated from 15 mmHg to 19 mmHg [18–20]. However, these methods did not provide a direct relationship between the elastic modulus and intraocular pressure. Ultrasonic waves have also been used to measure the pressure in tissue directly. Eisenbrey et al. used a subharmonic aided pressure estimation (SHAPE) method to monitor portal hypertension in patients [21, 22]. Although this is a relatively noninvasive technique, a microbubble is still needed to generate a harmonic wave, and for the eye, this would require a surgical procedure to inject a bubble and thus would be rather invasive. Zhang et al. used a magnetic shaker to generate low frequency shear waves to estimate carpal tunnel pressure by measuring the shear wave speed in the tendon [23, 24]. The obstacle to apply this technique in the eye is that the shaker would vibrate the eye and potentially damage the delicate ocular structures. Recently, Liu et al. used a quantitative ultrasound spectroscopy method and the reflection amplitude method to measure the acoustic impedance of the canine cornea at a fixed intraocular pressure and found that the acoustic impedance has a linear relationship with the elastic modulus of the cornea, which increased nonlinearly with stress [25–27]. As mentioned above, the relationship of IOP and elasticity was also found to be linearly correlated. Theoretically, the correlation of IOP and corneal acoustic impedance could be achieved, and this relationship is expected to be linear. It is simple to characterize the IOP level directly based on the change of acoustic impedance, which is easily obtained by measuring the pulse echo of ultrasonic waves. So far, to the best of our knowledge, no study has reported on the relationship of IOP and corneal acoustic impedance.

In this study, we used a custom-made highly focused transducer (47 MHz, F# 0.9) to acquire the anteriorly and posteriorly reflected signal of cadaver porcine eyes at intraocular pressure ranging from 10 mmHg to 50 mmHg with an increasing interval of 10 mmHg. The reflected signals were used to calculate the acoustic impedance of the anterior and posterior segment of the cornea, which were based on reflection and transmission coefficient models established by continuous monitoring of stress and strain. Statistical analysis was used to correlate the acoustic impedance and the intraocular pressure.

2. Material and Method

2.1. Sample Preparation

Six fresh porcine eyes were collected from a local slaughterhouse and stored in a box at 0°C. The eyeball was place on a rubber pad with a hole in the center to maintain its position during inflation due to pressure elevation. Four pins were inserted into the extraocular muscles and the rubber pad to fix the eyeball. A 30-gauge needle connected with an infusion line was inserted though the limbus into the anterior chamber for IOP elevation. Balanced salt solution was used for both insertion into the chamber and immersion of the eyeball during ultrasound reflection testing. The IOP was decided by the height of saline bottle. The experimental temperature is 20°C. The same eye was measured three times with the same condition, and the mean values were used for analysis.

2.2. Transducers and System

A highly focused transducer was made to measure the echo from the anterior and posterior surfaces of the cornea, respectively. The frequency of the transducer is 47-MHz with a bandwidth of 62%. The element is 4 mm diameter and made of LiNbO₃ (Boston Piezo-Optics, Bellingham, MA). Matching layers were used to acquire the large bandwidth. A steel ball with diameter of 7.2 mm was used to make foci of 3.6 mm. The pulse echo and frequency spectrum of the transducer are shown in Figure 1.

The ultrasound system and pressure elevation setup are shown in Figure 2. The eyeball was fixed on a rubber holder. The transducer was aligned to the center of the eyeball by motor controller (ESP301, Newport) and scanned vertically with a step of 1 nm. The pulser/receiver (Panametrics 5900PR, Olympus, Waltham, MA) was set with a bandwidth from 1 MHz to 100 MHz to transmit and receive ultrasound signals. A high-speed acquisition card (GAGE) was used to acquire the signal with a sampling frequency of 1 GHz and 10 A-line scans and transmit it to personal computer via custom software. An infusion bottle with a needle was used to inject the balanced salt solution into the anterior chamber by the height difference between the bottleneck and eyeball. The height of the bottle could be adjusted to change the level of IOP.

2.3. Data Processing

To accurately adjust the transducer position to get the maximum pulse echo, one-dimensional scanning data was used to get the maximum amplitude and related time of flight by fitting of the acquired signals to a parabolic envelope curve.

Reflection coefficient models were used to calculate the acoustic impedance of the anterior and posterior cornea, respectively, and the highly focused transducer could completely separate the two reflected signals. For the anterior reflection, the reflection coefficient was only dependent on the acoustic impedance ratio of the immersion saline and anterior cornea. So the cornea impedance could be calculated by the following formula [27] when the impedance of saline is known:

is the max amplitude of the reflected signal, while is the amplitude of the incident signal. The echo of a quartz with known density and longitudinal wave velocity was used as the reference signal to calculate by a similar formula as (1). As the reflection interface was not planar, the curvature coefficient was used for modification as mentioned in [28].

For posterior reflection, as Young’s modulus of the anterior and posterior cornea proved to be different [29], the acoustic impedance was also considered separately. The saline, cornea, and anterior chamber could be considered as a three-layered medium. So the incident signal and the reflected signal has the following relationship:

, , , and are the acoustic impedance of saline, anterior cornea, posterior cornea, and aqueous humor, respectively.

The Pearson correlation analysis [27] was used to correlate the corneal acoustic impedance and intraocular pressure. The correlation coefficients were calculated by functions provided by MATLAB Software.

3. Results

Figure 3 shows the time of flight difference at different intraocular pressures. For the anterior reflected signals, the difference of time that received the focused echo was limited by ±5 ns, which indicated that the experiment and data processing method exhibited good repeatability. For the posterior reflected signals, the difference of time of flight increased with IOP elevation, according to the fitted curve. There is an inconsistence that the thickness of the cornea would change or not with the IOP elevation [30, 31]. If the thickness of the cornea was not affected by changing IOP from 0 to 50 mmHg [30], the time of flight change is dominated by the velocity change of the cornea, which was induced by the elasticity change during pressure elevation. However, if the cornea thickness changed [31], the time of flight will affected both by velocity and thickness changing, which may lead to a nonlinear relationship as described in Figure 3(b). We cannot make a direct relationship between time of flight and IOP change due to complex deformation of eyeball. So we propose an acoustic impedance method which is not affected by the thickness change of cornea.

(a)

(b)

Figure 4 is the amplitude of the echo from the anterior and posterior interface. The dots represent the experimental amplitude of the echo as a function of the IOP level, while the solid line represents a linear fitting line. During the pressure elevation from 10 mm Hg to 50 mm Hg, the anterior amplitude increased from to and the posterior amplitude increased from to . Strong correlation was observed (; ) for the posterior amplitude with IOP, while less significant correlation was observed (; ) for anterior amplitude.

(a)

(b)

Figure 5 is the acoustic impedance of the anterior and posterior cornea calculated by reflection model using the amplitude data in Figure 4. The dots represent the experimental acoustic impedance as a function of the IOP level, while the solid line represents a linearly fitted line. During the pressure elevation from 10 mm Hg to 50 mm Hg, the mean acoustic impedance of the anterior cornea increased from 1.5399 to 1.5519 MRayl. The mean value of posterior acoustic impedance increased from 1.5393 to 1.5698 MRayl. Strong correlation was observed (; ) for the posterior acoustic impedance with IOP, while less significant correlation was observed (; ) for anterior acoustic impedance.

(a)

(b)

4. Discussion

The measured corneal acoustic impedance showed a significant correlation with the intraocular pressure. However, data dispersion occurred due to the complex deformation behavior and elastic property change during pressure elevation.

For biomechanics analysis, the globe will get inflated as the intraocular pressure increases. Due to the irregular structure of the segments and their different properties, the deformation of cornea geometry, compression, and elasticity changes with IOP changes. According to Laplace law, the stress applied to the cornea is calculated using intraocular pressure , cornea radius , and cornea thickness :The corneal thickness would vary little or become thinner if the IOP changes in a reasonable range, as mentioned before [30, 31]. The corneal radius should increase during the inflation process, so according to (3), stress should increase during IOP elevation. The relationship between stress and strain for the cornea is nonlinear, which means that the elastic modulus (secant or tangent) should increase as stress increases. The linear relationship of stress and elastic modulus has been demonstrated by shell theory and prior experiments [14]. The matrix of the cornea is a viscoelastic and nearly incompressible material; this means that the density of the cornea should not vary during inflation [32]. According to the above discussion, IOP should have a direct effect on the acoustic impedance, which is the product of density and acoustic velocity. However, due to the complexity of the deformation process, an explicit expression of the relationship could not be found. Fortunately, the linear relationship between the acoustic impedance and elastic modulus of the cornea has been demonstrated experimentally in [27]. Thus, we can safely conclude that the IOP and cornea impedance has a linear relationship, which is also demonstrated in our experiments. However, data dispersion was observed. This may due to the assumptions and approximations we made, in spite of experimental error. For example, the change in cornea curvature, which will affect the reflection coefficient of incident ultrasound wave, was ignored.

It is important to accurately measure the acoustic impedance because a large difference in intraocular pressure could only induce a small change in impedance, as shown in Figure 5. Several technical methods should be taken to improve accuracy. Firstly, a precision device should be used to align the center of the ultrasonic transducer and the cornea for vertical scanning in an in vitro measurement. It is also possible to design a device for in vivo measurement of acoustic impedance, since the cornea could be reached directly by ultrasound without anatomic barriers. Secondly, both the anterior and posterior reflected signal of the cornea could be acquired by a highly focused transducer to calculate the impedance. Young’s modulus of the anterior and posterior cornea showed a large difference in [29]. However, a smaller difference was observed in acoustic impedance in our experiment. Although cornea is composed of five layers, the innermost and outermost layers are too thin to affect ultrasound propagation, so the acoustic property of the corneal stroma is dominant. The similar anterior and posterior corneal acoustic impedance could be mutual verification to each other. Thirdly, steps should be taken to avoid or minimize the hydration of cornea during IOP elevation in vitro studies. The 20% dextran-saline solution could be used for insertion solution and immersion solution. The immersion time should be controlled to a few minutes to avoid the change of central corneal thickness.

Further studies include the following. The cornea may reach a high strain level (4%) during intraocular elevation, which could be calculated by Laplace law with the experimental parameters. The nonlinear relationship of stress and strain should be considered in the correlation between acoustic impedance and IOP. An inflation test combined with ultrasonic test should be done together, so that the elastic modulus could be calculated by shell theory model, and cornea geometry change could be recorded by laser device [14]. The strain could be calculated from ultrasonic shear wave imaging [15] and ultrasound speckle tracking [18]. The abundant data could be used to rebuild and validate the inflation model. Animal studies are expected to demonstrate more accurate correlation of corneal acoustic impedance and intraocular pressure since it is more physiologic. A low intensity level ultrasound test could be used for in vivo study to validate the possibility of intraocular pressure measurement using acoustic impedance of the cornea.

In summary, a strong correlation between intraocular pressure and acoustic impedance of the cornea was demonstrated in the present study. This correlation provided a potentially important method to noninvasively measure the intraocular pressure in vivo. This method of intraocular pressure measurement may prove useful for corneal biomechanics research as well as glaucoma clinical studies in the future.

Disclosure

The authors alone are responsible for the content and writing of the paper.

Conflicts of Interest

The authors report no conflicts of interest.

Authors’ Contributions

Jun Zhang and Yi Zhang contributed equally.

Acknowledgments

This work was supported by the NIH Grants R01EY026091 and P41-EB2002182. The work was also funded by China Scholarship Council.

References

H. Quigley and A. T. Broman, “The number of people with glaucoma worldwide in 2010 and 2020,” British Journal of Ophthalmology, vol. 90, no. 3, pp. 262–267, 2006.
View at: Publisher Site | Google Scholar
D. S. Friedman, R. C. Wolfs, B. J. O'Colmain et al., “Prevalence of open-angle glaucoma among adults in the United States,” Archives of Ophthalmology, vol. 122, no. 4, pp. 532–538, 2004.
View at: Publisher Site | Google Scholar
D. Touboul, C. Roberts, J. Kérautret et al., “Correlations between corneal hysteresis, intraocular pressure, and corneal central pachymetry,” Journal of Cataract and Refractive Surgery, vol. 34, no. 4, pp. 616–622, 2008.
View at: Publisher Site | Google Scholar
A. S. Wang, L. M. Alencar, R. N. Weinreb et al., “Repeatability and reproducibility of goldmann applanation, dynamic contour, and ocular response analyzer tonometry,” Journal of Glaucoma, vol. 22, no. 2, pp. 127–132, 2013.
View at: Publisher Site | Google Scholar
D. A. Luce, “Determining in vivo biomechanical properties of the cornea with an ocular response analyzer,” Journal of Cataract and Refractive Surgery, vol. 31, no. 1, pp. 156–162, 2005.
View at: Publisher Site | Google Scholar
K. Mansouri, J. H. K. Liu, R. N. Weinreb, A. Tafreshi, and F. A. Medeiros, “Analysis of continuous 24-hour intraocular pressure patterns in glaucoma,” Investigative Ophthalmology and Visual Science, vol. 53, no. 13, pp. 8050–8056, 2012.
View at: Publisher Site | Google Scholar
J. M. Martinez-De-La-Casa, J. Garcia-Feijoo, A. Fernandez-Vidal, C. Mendez-Hernandez, and J. Garcia-Sanchez, “Ocular response analyzer versus Goldmann applanation tonometry for intraocular pressure measurements,” Investigative Ophthalmology and Visual Science, vol. 47, no. 10, pp. 4410–4414, 2006.
View at: Publisher Site | Google Scholar
S. Kling, H. Ginis, and S. Marcos, “Corneal biomechanical properties from two-dimensional corneal flap extensiometry: Application to UV-Riboflavin cross-linking,” Investigative Ophthalmology and Visual Science, vol. 53, no. 8, pp. 5010–5015, 2012.
View at: Publisher Site | Google Scholar
N. Garcia-Porta, P. Fernandes, A. Queiros, J. Salgado-Borges, M. Parafita-Mato, and J. M. González-Méijome, “Corneal biomechanical properties in different ocular conditions and new measurement techniques,” ISRN Ophthalmology, vol. 2014, Article ID 724546, 19 pages, 2014.
View at: Publisher Site | Google Scholar
K. E. Hamilton and D. C. Pye, “Young's modulus in normal corneas and the effect on applanation tonometry,” Optometry and Vision Science, vol. 85, no. 6, pp. 445–450, 2008.
View at: Publisher Site | Google Scholar
E. Sjøntoft and C. Edmund, “In vivo determination of young's modulus for the human cornea,” Bulletin of Mathematical Biology, vol. 49, no. 2, pp. 217–232, 1987.
View at: Publisher Site | Google Scholar
K. Anderson, A. El-Sheikh, and T. Newson, “Application of structural analysis to the mechanical behaviour of the cornea,” Journal of the Royal Society Interface, vol. 1, no. 1, pp. 3–15, 2004.
View at: Publisher Site | Google Scholar
A. Elsheikh, D. Alhasso, and P. Rama, “Biomechanical properties of human and porcine corneas,” Experimental Eye Research, vol. 86, no. 5, pp. 783–790, 2008.
View at: Publisher Site | Google Scholar
A. Elsheikh, D. Wang, and D. Pye, “Determination of the modulus of elasticity of the human cornea,” Journal of Refractive Surgery, vol. 23, no. 8, pp. 808–818, 2007.
View at: Google Scholar
T.-M. Nguyen, J.-F. Aubry, D. Touboul et al., “Monitoring of cornea elastic properties changes during UV-A/riboflavin-induced corneal collagen cross-linking using supersonic shear wave imaging: A pilot study,” Investigative Ophthalmology and Visual Science, vol. 53, no. 9, pp. 5948–5954, 2012.
View at: Publisher Site | Google Scholar
T.-M. Nguyen, J.-F. Aubry, M. Fink, J. Bercoff, and M. Tanter, “In vivo evidence of porcine cornea anisotropy using supersonic shear wave imaging,” Investigative Ophthalmology and Visual Science, vol. 55, no. 11, pp. 7545–7552, 2014.
View at: Publisher Site | Google Scholar
M. Tanter, D. Touboul, J.-L. Gennisson, J. Bercoff, and M. Fink, “High-resolution quantitative imaging of cornea elasticity using supersonic shear imaging,” IEEE Transactions on Medical Imaging, vol. 28, no. 12, pp. 1881–1893, 2009.
View at: Publisher Site | Google Scholar
B. Cruz Perez, E. Pavlatos, H. J. Morris et al., “Mapping 3D Strains with Ultrasound Speckle Tracking: Method Validation and Initial Results in Porcine Scleral Inflation,” Annals of Biomedical Engineering, vol. 44, no. 7, pp. 2302–2312, 2016.
View at: Publisher Site | Google Scholar
J. Tang and J. Liu, “Ultrasonic measurement of scleral cross-sectional strains during elevations of intraocular pressure: Method validation and initial results in posterior porcine sclera,” Journal of Biomechanical Engineering, vol. 134, no. 9, Article ID 91007, 2012.
View at: Publisher Site | Google Scholar
K. Gijsbertse, A. M. Sprengers, M. M. Nillesen, H. H. G. Hansen, N. Verdonschot, and C. L. De Korte, “Three-dimensional ultrasound strain imaging of skeletal muscles,” in Proceedings of the IEEE International Ultrasonics Symposium, IUS 2015, twn, October 2015.
View at: Publisher Site | Google Scholar
J. R. Eisenbrey, J. K. Dave, V. G. Halldorsdottir et al., “Subharmonic aided pressure estimation in patients with suspected portal hypertension,” in Proceedings of the 2011 IEEE International Ultrasonics Symposium, IUS 2011, pp. 620–623, October 2011.
View at: Publisher Site | Google Scholar
V. G. Halldorsdottir, J. K. Dave, L. M. Leodore et al., “Subharmonic contrast microbubble signals for noninvasive pressure estimation under static and dynamic flow conditions,” Ultrasonic Imaging, vol. 33, no. 3, pp. 153–164, 2011.
View at: Publisher Site | Google Scholar
Y. Wang, B. Qiang, X. Zhang et al., “A non-invasive technique for estimating carpal tunnel pressure by measuring shear wave speed in tendon: A feasibility study,” Journal of Biomechanics, vol. 45, no. 16, pp. 2927–2930, 2012.
View at: Publisher Site | Google Scholar
X. Zhang, B. Qiang, and J. Greenleaf, “Comparison of the surface wave method and the indentation method for measuring the elasticity of gelatin phantoms of different concentrations,” Ultrasonics, vol. 51, no. 2, pp. 157–164, 2011.
View at: Publisher Site | Google Scholar
J. Liu, X. He, X. Pan, and C. J. Roberts, “Ultrasonic model and system for measurement of corneal biomechanical properties and validation on phantoms,” Journal of Biomechanics, vol. 40, no. 5, pp. 1177–1182, 2007.
View at: Publisher Site | Google Scholar
X. He and J. Liu, “A quantitative ultrasonic spectroscopy method for noninvasive determination of corneal biomechanical properties,” Investigative Ophthalmology and Visual Science, vol. 50, no. 11, pp. 5148–5154, 2009.
View at: Publisher Site | Google Scholar
X. He and J. Liu, “Correlation of corneal acoustic and elastic properties in a canine eye model,” Investigative Ophthalmology and Visual Science, vol. 52, no. 2, pp. 731–736, 2011.
View at: Publisher Site | Google Scholar
J. Tang and J. Liu, “Variance of Speed of Sound and Correlation with Acoustic Impedance in Canine Corneas,” Ultrasound in Medicine and Biology, vol. 37, no. 10, pp. 1714–1721, 2011.
View at: Publisher Site | Google Scholar
G. Lombardo, S. Serrao, M. Rosati, and M. Lombardo, “Analysis of the viscoelastic properties of the human cornea using scheimpflug imaging in inflation experiment of eye globes,” PLoS ONE, vol. 9, no. 11, Article ID 0112169, 2014.
View at: Publisher Site | Google Scholar
M. J. Doughty and M. L. Zaman, “Human corneal thickness and its impact on intraocular pressure measures: a review and meta-analysis approach,” Survey of Ophthalmology, vol. 44, no. 5, pp. 367–408, 2000.
View at: Publisher Site | Google Scholar
Y.-W. Park, M.-B. Jeong, E. R. Lee et al., “Acute changes in central corneal thickness according to experimental adjustment of intraocular pressure in normal canine eyes,” Journal of Veterinary Medical Science, vol. 75, no. 11, pp. 1479–1483, 2013.
View at: Publisher Site | Google Scholar
T. M. Nejad, C. Foster, and D. Gongal, “Finite element modelling of cornea mechanics: A review,” Arquivos Brasileiros de Oftalmologia, vol. 77, no. 1, pp. 60–65, 2014.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2017 Jun Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

1273

Downloads

810

Citations