Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Corneal pulsation and biomechanics during induced ocular pulse. An ex-vivo pilot study

  • Maja M. Rogala,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Validation, Visualization, Writing – original draft, Writing – review & editing

    Affiliation Department of Mechanics, Materials Science and Engineering, Wrocław University of Science and Technology, Wrocław, Poland

  • Daniel Lewandowski,

    Roles Conceptualization, Methodology, Resources, Validation, Writing – review & editing

    Affiliation Department of Mechanics, Materials Science and Engineering, Wrocław University of Science and Technology, Wrocław, Poland

  • Jerzy Detyna,

    Roles Formal analysis, Resources, Validation, Writing – review & editing

    Affiliation Department of Mechanics, Materials Science and Engineering, Wrocław University of Science and Technology, Wrocław, Poland

  • Agnieszka Antończyk,

    Roles Resources, Writing – review & editing

    Affiliation Department of Surgery, Wrocław University of Environmental and Life Sciences, Wrocław, Poland

  • Monika E. Danielewska

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Resources, Supervision, Validation, Writing – original draft, Writing – review & editing

    monika.danielewska@pwr.edu.pl

    Affiliation Department of Biomedical Engineering, Wrocław University of Science and Technology, Wrocław, Poland

Abstract

The purpose of this study was to ascertain the relationships between the amplitude of the corneal pulse (CP) signal and the parameters of corneal biomechanics during ex-vivo intraocular pressure (IOP) elevation experiments on porcine eyes with artificially induced ocular pulse cycles. Two experiments were carried out using porcine eyes. In the first one, a selected eye globe was subjected to three IOP levels (15, 30 and 45 mmHg), where changes in physical ocular pulse amplitude were controlled by infusion/withdrawal volumes (ΔV). In the second experiment, six eyes were subjected to IOP from 15 mmHg to 45 mmHg in steps of 5 mmHg with a constant ΔV, where corneal deformation parameters were measured using Corvis ST. In both experiments, at each IOP, the CP and IOP signals were acquired synchronically using a non-contact ultrasonic distance sensor and a pressure transmitter, respectively. Based on the amplitudes of the CP and IOP signals ocular pulse based corneal rigidity index (OPCRI) was calculated. Results indicate positive correlations between ΔV and the physical ocular pulse amplitude, and between ΔV and the corneal pulse amplitude (both p < 0.001). OPCRI was found to increase with elevated IOP. Furthermore, IOP statistically significantly differentiated changes in OPCRI, the amplitudes of CP and IOP signals and in most of the corneal deformation parameters (p < 0.05). The partial correlation analysis, with IOP as a control variable, revealed a significant correlation between the length of the flattened cornea during the first applanation (A1L) and the corneal pulse amplitude (p = 0.002), and between A1L and OPCRI (p = 0.003). In conclusion, this study proved that natural corneal pulsations, detected with a non-contact ultrasonic technique, reflect pressure-volume dynamics and can potentially be utilized to assess stiffness of the cornea. The proposed new rigidity index could be a simple approach to estimating corneal rigidity.

Introduction

In recent years, corneal biomechanics has become of particular interest for improving refractive surgeries [1,2], better understanding pathogenesis of corneal diseases such as keratoconus [35] and eye diseases linked to intraocular pressure (IOP) elevation experienced in glaucoma [68]. Cornea, like many other soft tissues, represents viscoelastic model of the material [9,10]. The dynamic deformation response of the cornea is described by its viscoelastic properties and originates mainly from molecular rearrangement as a response to the mechanical load application [11,12]. Non-linear stress-strain relationship, described in the literature, indicates that the mechanical properties of the cornea change with the applied stress [13,14]. The cornea, as an integral part of the outer ocular coat, is directly subjected to the internal load of the IOP [13,15]. Thus, as it was reported earlier [7,1618], IOP itself affects stiffness of the ocular shell. In turn, the corneal mechanical resistance influences IOP fluctuations, e.g. Ocular Pulse Amplitude (OPA) and IOP spikes, as it has recently been proven in both in-vivo [7,19] and ex-vivo studies [20].

Current efforts are notably focused on developing an accurate and at the same time non-invasive in-vivo methodology for determining biomechanical properties of the cornea, which have a considerable clinical relevance. Here, one can specify air-puff systems for the dynamic corneal response evaluation with ultra-high-speed Scheimpflug camera [2124] or OCT [2528], Brillouin microscopy [29,30] and ocular pulse elastography [31]. The first is the only commercially available solution, which makes it the one most commonly used, however, the last one is particularly interesting because it aims at involving natural corneal deformation due to the ocular pulse (OP) phenomenon [32]. In general, OP characterizes pulse-synchronous eye volume changes resulting from the pulsatile variations in IOP [33], closely linked to the blood circulation [7,34] and aqueous humor dynamics [35].

Numerous efforts have been made in human in-vivo studies to investigate the influence of corneal biomechanical properties on IOP measurements [21,36,37], however, without any knowledge about the true IOP. Studies that considered comparing IOP measurements with the manometrically determined IOP are limited, mainly because of a limited number of patients and ethical restrictions [38,39]. On the other hand, the effect of corneal biomechanics on the accurate IOP readings can be estimated using an appropriate mathematical model, such as that proposed by Liu and Roberts [40]. Ex-vivo studies, however, provide an opportunity to link estimated IOP values to true IOP by taking into account controlled corneal biomechanical properties [16,41,42]. Including OP cycles simulations to such experiments can be a valuable source of additional knowledge about the corneal biomechanics intrinsic to the ocular dynamics. To the best of our knowledge, the first ex-vivo test with OP cycles simulation was performed very recently by Liu [32]. The applied methodology combined the advantage of in-vivo tests (a possibility to closely reflect the real pulsatile IOP conditions) with the advantage of ex-vivo tests (the ability to fully control those conditions). Those ex-vivo OP simulations gave the opportunity to thoroughly investigate the IOP dynamics simultaneously with the biomechanical changes of cornea.

Similarly, parameters of ocular dynamic due to OP have been assessed via examining the natural corneal pulse (CP) defined as a superposition of slight semi-periodic corneal surface expansion and longitudinal eye globe movements [4345]. Measurements of the CP signal, registered with non-contact and non-invasive ultrasonic technique [46], were performed both in humans [45,4750], and in animals [51]. The advantage of this technique is the ability to measure the natural pulsation of the cornea without using any external mechanical excitation system that would affect natural tissue dynamics. In our earlier study, in anesthetized rabbits [51], the CP signal’s parameters were referred to true IOP, however, without controlled corneal biomechanical properties at each IOP value. Combining elements of either the CP cycles and corneal biomechanics for different IOP conditions can bring new insight into their mutual dependencies.

The aim of this study was to ascertain the relationships between the corneal pulse amplitude and the parameters of corneal biomechanics during the ex-vivo IOP elevation experiments on porcine eyes with artificially induced OP cycles.

Materials and methods

Specimens preparation

Porcine eye globes were obtained from the local abattoir Dworecki’s Meat Processing Plant (Golejowo, Poland) at most 6 h post-mortem and tested within the maximum of 12 h post-mortem, similarly as in the work of Kling et al [52]. The residues of the eye muscles were mechanically separated from each eye, then the eye globes were placed in a storage medium of phosphate-buffered saline (PBS) solution at 4°C before the experiment began. Before qualifying porcine eye globes for the measurement, the corneas were carefully checked using a slit lamp in order to assess existing edema, endothelial damage, mechanical damage and corneal transparency. The inclusion criteria were: lack of mechanical damages of the cornea, conjunctiva and sclera, corneal transparency remaining throughout the entire procedure, optic nerve intact, and central corneal thickness (CCT) ranging between 950 μm and 1200 μm. This narrow range was used to ensure that eye globes have corneas of similar thickness, because it can influence the biomechanical properties of the tissue [53], and hence the ability of cornea to deform during the ocular pulse simulation.

Each eye globe underwent the same procedure. First, an eye was carefully placed in a custom-made holder internally padded with PBS-moistened cotton. Restraint of the whole eye movement was obtained by gently binding the optic nerve to the holder with the polyamide sewing thread. The thread was slightly tensed during the measurements so that the eye globe was lightly pressed to the rear wall of the holder. Then, optical biometry was performed using the IOL Master 700 (Carl Zeiss Meditec AG, Jena, Germany) to ensure uniformity of the samples in terms of their geometrical parameters and to exclude outliers. In order to maintain proper hydration of the cornea during the measurement, a few drops of Eusol-C (Alchimia, Ponte San Nicolo, Padova, Italy) were applied to corneal surface.

Dynamic inflation experiment with induced OP

Dynamic behavior of the corneal tissue was investigated during inflation test. From that moment, a mineral oil was applied to the eye globe surface to prevent loss of hydration. This procedure was also aimed at prevention of the tissue swelling [54]. The 19G injection needle was inserted into the anterior chamber of the eye from the corneoscleral limbus area in order to control the value of the intraocular pressure (IOP) and to provide the pressure fluctuations simulating the OP. This relatively large size of the needle was selected to avoid pressure loss and prevent delays in the experimental setup. The needle was attached to the pressure transmitter P-30 (WIKA, Germany) of relative pressure range 0–187.5 mmHg (0.05% accuracy class) and to the custom-made syringe pump via non-deformable polytetrafluoroethylene tubing and the set of quick couplings, all filled with PBS and free of air bubbles (Fig 1). The software consisted of two parts. The first one directly controlled the work of the pump motor. It operated on a microprocessor (Teensy 3.5, MK64FX512VMD12, 120 MHz, ARM Cortex-M4a and FPU) with codes written in C language. The second part, supervising the pumping process and communicating microprocessor with the PC, was written using Keysight VEE (Pro 9.3, Keysight Technologies Inc.).

thumbnail
Fig 1. A scheme of the inflation test rig for the CP signal measurement during induced ocular pulse (OP) cycles and for corneal biomechanical measurements.

https://doi.org/10.1371/journal.pone.0228920.g001

Two experiments were considered. In the first one, a selected eye globe was subjected to three levels of IOP (15, 30, and 45 mmHg), where the OP was generated as a sine wave with a frequency of 1.2 Hz and an amplitude, OPAp, corresponding to controlled infusion/withdrawal PBS volume forced by the pump. The test was performed for the five different volumes (ΔV): 40, 60, 80, 90 and 100 μl (Fig 2) set at each IOP value. OPAp stands for the physical ocular pulse amplitude value [16] and is the difference between peak and trough in the pressure signal.

thumbnail
Fig 2. A scheme of the CP and IOP signal measurement procedure during ocular pulse (OP) simulations with controlled infusion/withdrawal PBS volume (ΔV) in the three-step inflation experiment.

https://doi.org/10.1371/journal.pone.0228920.g002

During all the enabled OP simulations, the CP signal was acquired in a non-contact manner using ultrasonic distance sensor (UltraLab, Wrocław, Poland) [46]. Synchronically with the CP signal, the IOP signal was registered with the pressure transmitter. The three 10-second length continuous data acquisitions were performed at each setting.

The second experiment involved a number of eye globes which were subjected to IOP increased from 15 mmHg to 45 mmHg in steps of 5 mmHg. The infusion/withdrawal volume was adjusted to obtain pressure amplitude value at about 4 mmHg (mean ± SD, 4.25 ± 0.25 mmHg) for the IOP equal to 15 mmHg, which gives from about 70 μl to 90 μl of PBS. The determined volume was maintained at each IOP value for the given eye that, in turn, should lead to a natural increase of OPAp as a result of the increase of IOP according to the results obtained from previous in-vivo animal and human studies [7,55]. This effect is due to the decreasing ability of the ocular tissue to deform under higher internal pressure.

At each IOP, when the OP simulation was enabled, the CP and IOP signals were registered as in the first experiment. Also, in this experiment, static cases of IOP (without OP) were investigated (see Fig 3). When induction of OP was disabled, three corneal biomechanical measurements were performed using a dynamic Scheimpflug analyzer, Corvis ST (Oculus, Optikgeräte GmbH, Wetzlar, Germany). The device provided 10 corneal deformation parameters that were included in the analysis: the time from the start until the first/second applanation (A1T/A2T), the length of the flattened cornea during the first/second applanation (A1L/A2L), the corneal velocity at the two applanation events (A1V, A2V), and the four parameters describing the moment of the maximum deformation of the cornea: the amplitude (DA), the distance between corneal bending points (PD), the time from the start until the highest concavity of cornea (HCT), and the radius of the curvature at the corneal apex (R).

thumbnail
Fig 3. A scheme of the CP and IOP signal measurement procedure during ocular pulse (OP) simulations with constant infusion/withdrawal PBS volume (ΔV) together with corneal biomechanical measurements in the seven-step inflation experiment.

https://doi.org/10.1371/journal.pone.0228920.g003

Signal analysis

The CP and IOP signals were acquired at a sampling frequency of 400 Hz. CP was determined with the accuracy below 1 μm whereas IOP with a precision of 0.1 mmHg. Both signals were numerically processed using custom written program in MATLAB (MathWorks, Inc. Natick, MA, USA). In the case of the CP signal linear trend was removed. Then, as in the case of OPAp, the corneal pulse amplitude (CPA) was determined as the difference between the median value of successive peaks and the median value of the successive troughs in the signal.

In the second experiment, the considered parameters from the three measurements assigned to the specific IOP value were averaged and then normalized to their mean value calculated for the baseline IOP of 15 mmHg for each eye.

All experimental data are included in the S1 File of Supporting information.

Biomechanics

In order to estimate the stiffness of the corneal tissue in the particular experimental setup considered here, a new index named ocular pulse based corneal rigidity index (OPCRI) is introduced for the OP cycles. The methodology is analogous to the conventional calculation of the Young modulus E, however, it solves the problem of the inability to assess the exact value of E without extracting corneal tissue from the eye globe. The CPA and OPAp were taken as surrogate data for the stress-strain relationship, where OPAp, as an internal pressure change, substitutes the stress value, while CPA, as a corneal surface expansion, is the indirect measure of the resulting strain. Consequently, the higher the OPCRI value is obtained the more pressure is needed to yield a change in the ocular volume. Here, OPCRI was used to determine the corneal tissue stiffness at different IOP.

Statistical analysis

Statistical analyses were conducted using Statistica (StatSoft, Inc., USA). Firstly, in order to find if the experimental data meet the assumptions for parametric testing the Shapiro–Wilk W test of normality and Levene’s test for equality of variances were applied. Because not all of the parameters fulfilled the prerequisites of parametric test procedures the Friedman ANOVA test was employed to ascertain if the IOP statistically significantly differentiate the CP signal parameters and those describing the corneal biomechanics. In addition, depending on the Fisher’s exact test results, either linear or exponential modelling was used for those relationships. Furthermore, partial correlation was performed to investigate the dependencies between CPA or OPCRI and corneal deformation parameters, with IOP set as a control variable. The significance level α was set to 0.05 for all tests used in this study.

Results

Six porcine eye globes passed the inclusion criteria. The ocular biometric parameters with its mean and median values are gathered in the Table 1.

The first experiment showed that the elevation of IOP results in an increase of physical ocular pulse amplitude (OPAp) for the same infusion/withdrawal volume. Fig 4 shows the linear models of OPAp versus ΔV for each IOP value. The regression coefficients obtained for IOP = 15, 30 and 45 mmHg are 0.03, 0.05, and 0.07 mmHg/μl, respectively. Similar effect was observed for corneal pulse amplitude (CPA) (Fig 4). However, in this case corneal tissue response to changes in ΔV is weaker after raising IOP (the subsequent regression coefficients are: 0.23, 0.12 and 0.11 μm/μl).

thumbnail
Fig 4. Regression analysis of physical ocular pulse amplitude (OPAp), corneal pulse amplitude (CPA), and ocular pulse based corneal rigidity index (OPCRI), with infusion/withdrawal volume ΔV, for three different IOP values.

https://doi.org/10.1371/journal.pone.0228920.g004

An increase in OPCRI can be observed with the increase in IOP. However, OPCRI remains stable at the particular IOP (regression coefficients: −0.00004, −0.00088, and 0.00007 mmHg/(μm·μl) for IOP = 15, 30, and 45 mmHg, respectively) (see Fig 4).

In the second experiment, significant differences were present among the IOP in the changes of almost all the corneal deformation parameters (Friedman ANOVA, p < 0.05). The only exceptions were the A1L (Friedman ANOVA, p = 0.416), A2L (Friedman ANOVA, p = 0.087), and HCT (Friedman ANOVA, p = 0.208). IOP significantly differentiate also OP parameters (CPA and OPAp) and OPCRI (Friedman ANOVA, p = 0.013, p < 0.001, p < 0.001, respectively) (Fig 5). Table 2 summarizes these results.

thumbnail
Fig 5. Boxplots of deformation amplitude (DA), corneal pulse amplitude (CPA), ocular pulse based corneal rigidity index (OPCRI), and physical ocular pulse amplitude (OPAp) normalized to their baseline values obtained at 15 mmHg, for different IOP values.

https://doi.org/10.1371/journal.pone.0228920.g005

thumbnail
Table 2. Results of the Friedman ANOVA test presenting IOP impact on the changes in the corneal deformation parameters and the ocular pulse signal parameters.

In addition, depending on the Fisher’s exact test outcome, either linear or exponential model was applied.

https://doi.org/10.1371/journal.pone.0228920.t002

The results of partial correlation are shown in Table 3. Both CPA and OPCRI correlate with A1L, whereas only CPA correlates with A2T.

thumbnail
Table 3. Results of partial correlation between normalized corneal pulse amplitude (CPA) or ocular pulse based corneal rigidity index (OPCRI) and normalized corneal deformation parameters, with IOP used as a control variable.

https://doi.org/10.1371/journal.pone.0228920.t003

Discussion

The pressure-volume relationship in the anterior chamber is of particular interest for finding the effective method for in-vivo measuring ocular stiffness and for accurate IOP reading, which is crucial in diagnostics of ophthalmic diseases [18,41]. The pressure-volume dynamics is strictly related to IOP fluctuations and biomechanics of ocular tissues [7]. The pulsatile component of IOP comes from the pulsatile ocular blood flow induced by the arterial blood pulse [34]. The difference between the systolic and diastolic IOP values in the IOP characteristic is named the ocular pulse amplitude (OPA), whose average values range from 2 mmHg to 4 mmHg in healthy subjects [5658]. In clinical practice, the dynamic contour tonometer (DCT; Pascal tonometer) is the only currently commercially available instrument to in-vivo measure the IOP wave in a contact manner after the cornea anesthesia [59]. This measurement procedure prevents the natural expansion of the cornea and influences corneal biomechanical properties which are crucial to accurately estimate the IOP. This study contributes to this knowledge by providing information on the relationship between corneal pulsation and the IOP and its fluctuations during the ex-vivo IOP elevation experiment with porcine eyes with artificially induced OP cycles using the ultrasonic technique which enables registering natural corneal pulsation without any external stimuli and tissue contact. The use of porcine eyes was dictated by the fact that they share many anatomical and physiological similarities to those of the human eye [6062]. In many ex-vivo studies related to biomechanical properties of eye tissues, the porcine cornea can be used as a substitute model for human cornea study, according to the stress-strain pattern [42,63,64].

For the first ex-vivo inflation experiment, the results show that both physical ocular pulse amplitude (OPAp) and corneal pulse amplitude (CPA) positively correlate with infusion/withdrawal volume. However, the IOP elevation magnifies the relationship between changes in the internal pressure and those in the ocular volume (OPA versus ΔV) and, simultaneously, weakens the corneal tissue response to the changes in ocular volume (CPA versus ΔV). An increase of OPCRI—an index that was found to be uncorrelated with ΔV—at elevated IOP confirmed that biomechanical properties of the cornea change with the applied stress leading to rigidity increase. Hence, the corneal stiffness growth, concurrent with the IOP raise, prevents tissue ability to deform (represented by CPA) but also causes an increase in the amplitude of IOP changes (represented by OPAp). These findings are in agreement with a study by Clayson et al. [20] on porcine eyes, which showed that the induced corneal stiffening, obtained by crosslinking treatment, and the IOP rise significantly impact the magnitude of IOP spikes. In human studies, the positive correlation between IOP and OPA has been also reported in healthy subjects [58] as well as in patients with ophthalmic diseases [7,19] showing increased mechanical resistance of the ocular wall at high IOP.

In the second ex-vivo experiment, the IOP impact on corneal deformation parameters and OP signal parameters was assessed. The results of corneal deformation obtained using Corvis ST are consistent with those of Bao and colleagues [16] showing that the majority of biomechanical metrics are correlated with IOP and their relationship can be well described using either linear or exponential model (see Table 2).

In the current study, a new approach, which additionally takes OP signal parameters into account, has shown that OP signal parameters are correlated with IOP. When comparing CPA and deformation amplitude changes with IOP rise, one can clearly see that increasing the internal load reduces corneal ability to expand naturally during OP cycles and, at the same time, to deform in response to the air puff. This conclusion is consistent with OPCRI behavior under incremental rise of IOP, highlighting the increase in corneal rigidity with the load, as it has also been shown in the first experiment.

The OPAp has been found to increase with IOP rise, even though it was triggered with the constant ΔV at all IOP values. This observation confirms the results obtained in the first experiment of this study and also corresponds to the previous in-vivo animal [55] and human [7] studies demonstrating the increased OPA with artificial increase in IOP. Dastiridou and colleagues [7] suggested that a higher OPA may appear in eyes with increased ocular rigidity.

Ocular rigidity is usually obtained for the whole pressure-volume characteristic of the eye, which is approximated as a linear [65] or exponential model [7]. In this paper, the rigidity index, OPCRI, was introduced as a measure of the corneal response to the induced internal pressure pulsation imitating OP. This means that OPCRI calculation is focused on the small recurring changes of the load rather than a single inflation event. It is worth noting that OPCRI is based on corneal pulsation measured without any external stimulation, in contrast to Corvis ST, where complete applanation cycle of the cornea is required to acquire this information. OPCRI enables simple and indirect estimation of the tissue rigidity at different IOP. The study shows that the relationship between IOP and corneal rigidity in the form of OPCRI, as well as between IOP and CPA, can be modelled in a linear fashion. It is worth emphasizing that OPCRI is not the only alternative (with respect to Corvis ST) approach to estimating corneal stiffness. Another measure, termed ocular pulse stiffness index (OPSI), based on induced IOP pulsation was proposed by Pavlatos and colleagues [31]. The OPSI was introduced to verify effectiveness of the ocular pulse ultrasound elastography aiming at evaluating corneal biomechanics on the basis of naturally occurring OP and without the necessity of exerting external force needed to induce tissue deformation. However, the main difference between OPSI and OPCRI lies in the applied ultrasound techniques, which, in the case of OPCRI, allows registering corneal surface pulsation in an air-coupled manner without disturbing natural corneal dynamics. Worth noting is that both techniques are complementary as one provides high spatial resolution whereas the other provides high temporal resolution.

In this study, the relationships between corneal deformation parameters and OP signal parameters have been examined taking into account the influence of IOP on these parameters (please refer to Table 3). The results indicate that A1L correlates with both CPA and OPCRI, underlining the role of corneal biomechanical properties in changes of natural corneal pulsation. Recent studies have revealed that the applanation length is mostly affected by CCT and age of the subjects rather than by a change of IOP [66]. Hence, it was suggested that A1L could be used to determine the stiffness of the cornea. In the light of the above, the current study demonstrates that CPA values could be utilized for estimating, in an indirect way, corneal biomechanical changes linked to the higher corneal stiffness.

Limitations of this study include potential postmortem reactions of the ex-vivo specimens, e.g., corneal swelling, on the resulting analysis. Special effort was taken in order to minimize those effects during the experimental procedure. Nevertheless, the amount of any bias resulted from corneal swelling is likely to be the same for all considered cornea. Eusol-C and mineral oil were applied to maintain proper hydration of the tissue. The study relies on ex-vivo porcine eye globes substituting human eyes, which are difficult to obtain from human donors. Even so, engaging porcine model proved to be quite problematic. Some elements of the procedures used at abattoir may result in optical or mechanical defects of the eye globes, e.g., corneal opacity or impaired corneal epithelium. Therefore, only very limited number of porcine eye globes satisfied the inclusion criteria. In that sense, our study has a pilot character. In the first experiment, increasing infusion/withdrawal PBS volume at the particular IOP value, with a frequency maintained constant, results in both raise of the OPAp and acceleration of the pressure change during OP cycles. This small variations of the velocity of the load were neglected in this study, however, we are conscious that this aspect is important in the context of viscoelastic properties of the cornea [10,67]. Constant frequency of the OP simulation was set to reflect typical heart rhythm in a healthy human.

Conclusions

Summarizing, this study revealed that the corneal biomechanical changes related to increase in IOP can be detected by observing natural corneal pulsations with a non-contact ultrasonic technique. Specifically, it was shown that the increase in IOP fluctuations and corneal stiffness co-occurring with the IOP growth could be estimated indirectly based on the corneal pulse amplitude. The newly proposed rigidity index, calculated on the basis of the CP signal, could be an approximate estimator of the corneal stiffness corresponding to both IOP variations and changes in biomechanical properties of the cornea. More importantly, unlike in elastography, this knowledge can be acquired without disturbing natural corneal dynamic.

Acknowledgments

The authors would like to thank Prof. D. Robert Iskander for his support and fruitful comments on an earlier version of the manuscript, and also Dr Cezary Sielużycki and Dr Małgorzata Kostyszak for proofreading.

References

  1. 1. Sinha Roy A, Dupps WJ, Roberts CJ. Comparison of biomechanical effects of small-incision lenticule extraction and laser in situ keratomileusis: Finite-element analysis. J Cataract Refract Surg. 2014;40: 971–980. pmid:24857440
  2. 2. Frings A, Linke SJ, Bauer EL, Katz T. Effects of laser in situ keratomileusis (LASIK) on corneal biomechanical measurements with the Corvis ST tonometer. Clin Ophthalmol. 2015; 305–311. pmid:25709393
  3. 3. Shah S, Laiquzzaman M, Bhojwani R, Mantry S, Cunliffe I. Assessment of the biomechanical properties of the cornea with the ocular response analyzer in normal and keratoconic eyes. Invest Ophthalmol Vis Sci. 2007;48: 3026–3031. pmid:17591868
  4. 4. Wolffsohn JS, Safeen S, Shah S, Laiquzzaman M. Changes of corneal biomechanics with keratoconus. Cornea. 2012;31: 849–854. pmid:22495031
  5. 5. Johnson RD, Nguyen MT, Lee N, Hamilton DR. Corneal biomechanical properties in normal, forme fruste keratoconus, and manifest keratoconus after statistical correction for potentially confounding factors. Cornea. 2011;30: 516–523. pmid:21045653
  6. 6. Coudrillier B, Tian J, Alexander S, Myers KM, Quigley HA, Nguyen TD. Biomechanics of the human posterior sclera: Age- and glaucoma-related changes measured using inflation testing. Invest Ophthalmol Vis Sci. 2012;53: 1714–1728. pmid:22395883
  7. 7. Dastiridou AI, Ginis HS, de Brouwere D, Tsilimbaris MK, Pallikaris IG. Ocular rigidity, ocular pulse amplitude, and pulsatile ocular blood flow: The effect of intraocular pressure. Invest Ophthalmol Vis Sci. 2009;50: 5718–5722. pmid:19608534
  8. 8. Pérez-Merino P, Bekesi N, Fernández-López A. Intraocular pressure and three-dimensional corneal biomechanics. In: Koprowski R, editor. Air-Puff Tonometers Challenges and insights. IOP Publishing Ltd; 2019.
  9. 9. Fung Y. C. Biomechanics: Mechanical Properties of Living Tissues. New York: Springer Science & Business Media; 1993.
  10. 10. Hatami-Marbini H. Viscoelastic shear properties of the corneal stroma. J Biomech. 2014;47: 723–728. pmid:24368145
  11. 11. Bell JS, Hayes S, Whitford C, Sanchez-Weatherby J, Shebanova O, Vergari C, et al. The hierarchical response of human corneal collagen to load. Acta Biomater. 2018;65: 216–225. pmid:29128531
  12. 12. Holzapfel GA. Biomechanics of soft tissue. In: Lemaitre J, editor. The handbook of materials behavior models. Boston: Academic Press; 2001. pp. 1049–1063.
  13. 13. Anderson K, El-Sheikh A, Newson T. Application of structural analysis to the mechanical behaviour of the cornea. J R Soc Interface. 2004;1: 3–15. pmid:16849148
  14. 14. Kling S, Hafezi F. Corneal biomechanics–a review. Ophthalmic Physiol Opt. 2017;37: 240–252. pmid:28125860
  15. 15. Śródka W. Biomechanical model of human eyeball and its applications. Opt Appl. 2009;39: 401–413.
  16. 16. Bao F, Deng M, Wang Q, Huang J, Yang J, Whitford C, et al. Evaluation of the relationship of corneal biomechanical metrics with physical intraocular pressure and central corneal thickness in exvivo rabbit eye globes. Exp Eye Res. 2015;137: 11–17. pmid:26026878
  17. 17. Morris HJ, Tang J, Perez BC, Pan X, Hart RT, Weber PA, et al. Correlation between biomechanical responses of posterior sclera and IOP elevations during micro intraocular volume change. Investig Ophthalmol Vis Sci. 2013;54: 7215–7222.
  18. 18. Girard MJA, Dupps WJ, Baskaran M, Scarcelli G, Yun SH, Quigley HA, et al. Translating ocular biomechanics into clinical practice: Current state and future prospects. Curr Eye Res. 2015;40: 1–18. pmid:24832392
  19. 19. Knecht PB, Bosch MM, Michels S, Mannhardt S, Schmid U, Bosch MA, et al. The ocular pulse amplitude at different intraocular pressure: a prospective study. Acta Ophthalmol. 2011;89: e466–e471. pmid:21401909
  20. 20. Clayson K, Pan X, Pavlatos E, Short R, Morris H, Hart RT, et al. Corneoscleral stiffening increases IOP spike magnitudes during rapid microvolumetric change in the eye. Exp Eye Res. 2017;165: 29–34. pmid:28864177
  21. 21. Leung CKS, Ye C, Weinreb RN. An ultra-high-speed Scheimpflug camera for evaluation of corneal deformation response and its impact on IOP measurement. Investig Ophthalmol Vis Sci. 2013;54: 2885–2892.
  22. 22. Hon Y, Lam AKC. Corneal deformation measurement using Scheimpflug noncontact tonometry. Optom Vis Sci. 2013;90: e1–e8. pmid:23238261
  23. 23. Ambrósio RJ, Ramos I, Luz A, Faria FC, Steinmueller A, Krug M, et al. Dynamic ultra high speed Scheimpflug imaging for assessing corneal biomechanical properties. Rev Bras Oftalmol. 2013;72: 99–102.
  24. 24. Jędzierowska M. Corvis ST tonometer and the possibility of analysing corneal deformation dynamics during intraocular pressure measurement. In: Koprowski R, editor. Air-Puff Tonometers Challenges and insights. IOP Publishing Ltd; 2019.
  25. 25. Gora M, Karnowski K, Szkulmowski M, Kaluzny BJ, Huber R, Kowalczyk A, et al. Ultra high-speed swept source OCT imaging of the anterior segment of human eye at 200 kHz with adjustable imaging range. Opt Express. 2009;17: 14880–14894. pmid:19687967
  26. 26. Alonso-Caneiro D, Karnowski K, Kaluzny BJ, Kowalczyk A, Wojtkowski M. Assessment of corneal dynamics with high-speed swept source Optical Coherence Tomography combined with an air puff system. Opt Express. 2011;19: 14188–14199. pmid:21934782
  27. 27. Dorronsoro C, Pascual D, Pérez-Merino P, Kling S, Marcos S. Dynamic OCT measurement of corneal deformation by an air puff in normal and cross-linked corneas. Biomed Opt Express. 2012;3: 473–487. pmid:22435096
  28. 28. Wang S, Larin KV, Li J, Vantipalli S, Manapuram RK, Aglyamov S, et al. A focused air-pulse system for optical-coherence-tomography-based measurements of tissue elasticity. Laser Phys Lett. 2013;10: 075605. pmid:29805349
  29. 29. Scarcelli G, Pineda R, Yun SH. Brillouin optical microscopy for corneal biomechanics. Investig Ophthalmol Vis Sci. 2012;53: 185–190.
  30. 30. Scarcelli G, Kling S, Quijano E, Pineda R, Marcos S, Yun SH. Brillouin microscopy of collagen crosslinking: Noncontact depth-dependent analysis of corneal elastic modulus. Investig Ophthalmol Vis Sci. 2013;54: 1418–1425.
  31. 31. Pavlatos E, Chen H, Clayson K, Pan X, Liu J. Imaging corneal biomechanical responses to ocular pulse using high-frequency ultrasound. IEEE Trans Med Imaging. 2018;37: 663–670. pmid:29408793
  32. 32. Liu J, Pavlatos E, Pan X, Clayson K. Influence of physiological parameters on corneal deformation in response to ocular pulse. IOVS. 2017; ARVO E-Abstract 4334.
  33. 33. Galin MA, Baras I, Best M. The nature of the ocular pulse. J Head Face Pain. 1969;9: 112–118.
  34. 34. Silver DM, Farrell RA, Langham ME, O’Brien V, Schilder P. Estimation of pulsatile ocular blood flow from intraocular pressure. Acta Ophthalmol Suppl (Oxf). 1989;191: 25–29.
  35. 35. Ramos RF, Stamer WD. Effects of cyclic intraocular pressure on conventional outflow facility. Investig Ophthalmol Vis Sci. 2008;49: 275–281.
  36. 36. Medeiros FA, Weinreb RN. Evaluation of the influence of corneal biomechanical properties on intraocular pressure measurements using the ocular response analyzer. J Glaucoma. 2006;15: 364–370. pmid:16988597
  37. 37. Touboul D, Roberts C, Kérautret J, Garra C, Maurice-Tison S, Saubusse E, et al. Correlations between corneal hysteresis, intraocular pressure, and corneal central pachymetry. J Cataract Refract Surg. 2008;34: 616–622. pmid:18361984
  38. 38. Eisenberg DL, Sherman BG, McKeown CA, Schuman JS. Tonometry in adults and children: A manometric evaluation of pneumatonometry, applanation, and TonoPen in vitro and in vivo. Ophthalmology. 1998;105: 1173–1181. pmid:9663218
  39. 39. Whitacre MM, Stein RA, Hassanein K. The effect of corneal thickness on applanation tonometry. Am J Ophthalmol. 1993;115: 592–596. pmid:8488910
  40. 40. Liu J, Roberts CJ. Influence of corneal biomechanical properties on intraocular pressure measurement. J Cataract Refract Surg. 2005;31: 146–155. pmid:15721707
  41. 41. Eliasy A, Chen KJ, Vinciguerra R, Maklad O, Vinciguerra P, Ambrósio R, et al. Ex-vivo experimental validation of biomechanically-corrected intraocular pressure measurements on human eyes using the CorVis ST. Exp Eye Res. 2018;175: 98–102. pmid:29908883
  42. 42. Elsheikh A, Alhasso D, Rama P. Biomechanical properties of human and porcine corneas. Exp Eye Res. 2008;86: 783–790. pmid:18396276
  43. 43. Kowalska MA, Kasprzak HT. Ultrasonic measurement of binocular longitudinal corneal apex movements and their correlation to cardiopulmonary system. Biocybern Biomed Eng. 2008;28: 35–43.
  44. 44. Northrop RB, Nilakhe SS. A no-touch ocular pulse measurement system for the diagnosis of carotid occlusions. IEEE Trans Biomed Eng. 1977;24: 139–148. pmid:892817
  45. 45. Iskander RD, Kasprzak HT. Dynamics in longitudinal eye movements and corneal shape. Ophthalmic Physiol Opt. 2006;26: 572–579. pmid:17040421
  46. 46. Licznerski TJ, Jaroński J, Kosz D. Ultrasonic system for accurate distance measurement in the air. Ultrasonics. 2011;51: 960–965. pmid:21669450
  47. 47. Kowalska MA, Kasprzak HT, Iskander DR, Danielewska M, Mas D. Ultrasonic in vivo measurement of ocular surface expansion. IEEE Trans Biomed Eng. 2011;58: 674–680. pmid:21177153
  48. 48. Danielewska ME, Iskander DR, Kowalska M, Kasprzak HT. Phase dependencies between longitudinal corneal apex displacement and cardiovascular signals: Is the ocular pulse influenced by the electrical activity of the heart? Clin Exp Optom. 2012;95: 631–637. pmid:22827844
  49. 49. Danielewska ME, Iskander DR. Age and ocular dicrotism—spectral representations of corneal indentation pulse and cardiovascular signals in healthy and glaucomatous eyes. IOVS. 2015;56: ARVO E-Abstract 3369.
  50. 50. Kasprzak HT, Iskander DR. Spectral characteristics of longitudinal corneal apex velocities and their relation to the cardiopulmonary system. Eye. 2007;21: 1212–1219. pmid:16936640
  51. 51. Rogala MM, Danielewska ME, Antończyk A, Kiełbowicz Z, Rogowska ME, Kozuń M, et al. In-vivo corneal pulsation in relation to in-vivo intraocular pressure and corneal biomechanics assessed in-vitro. An animal pilot study. Exp Eye Res. 2017;162: 27–36. pmid:28689748
  52. 52. Kling S, Bekesi N, Dorronsoro C, Pascual D, Marcos S. Corneal viscoelastic properties from finite-element analysis of in vivo air-puff deformation. PLoS One. 2014;9: e104904. pmid:25121496
  53. 53. Kotecha A, Elsheikh A, Roberts CR, Zhu H, Garway-Heath DF. Corneal thickness- and age-related biomechanical properties of the cornea measured with the Ocular Response Analyzer. Invest Ophthalmol Vis Sci. 2006;47: 5337–5347. pmid:17122122
  54. 54. Hatami-Marbini H, Rahimi A. Effects of bathing solution on tensile properties of the cornea. Exp Eye Res. 2014;120: 103–108. pmid:24333541
  55. 55. Lawrence C, Schlegel WA. Ophthalmic pulse studies. I. Influence of intraocular pressure. Invest Ophthalmol. 1966;5: 515–525. pmid:5929287
  56. 56. Schwenn O, Troost R, Vogel A, Grus F, Beck S, Pfeiffer N. Ocular pulse amplitude in patients with open angle glaucoma, normal tension glaucoma, and ocular hypertension. Br J Ophthalmol. 2002;86: 981–984. pmid:12185121
  57. 57. Hoffmann EM, Grus F-H, Pfeiffer N. Intraocular pressure and ocular pulse amplitude using dynamic contour tonometry and contact lens tonometry. BMC Ophthalmol. 2004;4: 4. pmid:15038831
  58. 58. Kaufmann C, Bachmann LM, Robert YC, Thiel MA. Ocular pulse amplitude in healthy subjects as measured by dynamic contour tonometry. Arch Ophthalmol. 2006;124: 1104–1108. pmid:16908812
  59. 59. Kaufmann C, Bachmann LM, Thiel MA, The PU. Comparison of dynamic contour tonometry with Goldmann applanation tonometry. Invest Ophthalmol. 2004; 3118–3121.
  60. 60. Sanchez I, Martin R, Ussa F, Fernandez-Bueno I. The parameters of the porcine eyeball. Graefe’s Arch Clin Exp Ophthalmol. 2011;249: 475–482.
  61. 61. Gelatt KN, Gilger BC, Kern TJ. Veterinary ophthalmology. No.Ed.5. John Wiley & Sons; 2012.
  62. 62. Ruiz-Ederra J, García M, Hernández M, Urcola H, Hernández-Barbáchano E, Araiz J, et al. The pig eye as a novel model of glaucoma. Exp Eye Res. 2005;81: 561–569. pmid:15949799
  63. 63. Zeng Y, Yang J, Huang K, Lee Z, Lee X. A comparison of biomechanical properties between human and porcine cornea. J Biomech. 2001;34: 533–537. pmid:11266678
  64. 64. Bao F, Jiang L, Wang X, Zhang D, Wang Q, Zeng Y. Assessment of the ex vivo biomechanical properties of porcine cornea with inflation test for corneal xenotransplantation. J Med Eng Technol. 2012;36: 17–21. pmid:22085017
  65. 65. Pallikaris IG, Kymionis GD, Ginis HS, Kounis GA, Tsilimbaris MK. Ocular rigidity in living human eyes. Invest Ophthalmol Vis Sci. 2005;46: 409–414. pmid:15671262
  66. 66. Reisdorf S. Factors contributing to corneal deformation in ex vivo and in vivo measurements. In: Roberts CJ, Liu J, editors. Corneal Biomechanics: from theory to practice. Amsterdam: Kugler Publications; 2016. pp. 106–107.
  67. 67. Elsheikh A, Kassem W, Jones SW. Strain-rate sensitivity of porcine and ovine corneas. Acta Bioeng Biomech. 2011;13: 25–36. pmid:21761808