A viscoelastic biomechanical model of the cornea describing the effect of viscosity and elasticity on hysteresis

Corneal biomechanics in early diagnosis of keratoconus using artificial intelligence

Keratoconus is a blinding eye disease that affects activities of daily living; therefore, early diagnosis is crucial. Great efforts have been made toward an early diagnosis of keratoconus. Recent studies have shown that corneal biomechanics is associated with the occurrence and progression of keratoconus. Hence, detecting changes in corneal biomechanics may provide a novel strategy for early diagnosis. However, an early keratoconus diagnosis remains challenging due to the subtle and localized nature of its lesions. Artificial intelligence has been used to help address this problem. Herein, we reviewed the literature regarding three aspects of keratoconus (keratoconus, early keratoconus, and keratoconus grading) based on corneal biomechanical properties using artificial intelligence. Furthermore, we summarized the current research progress, limitations, and possible prospects.

Corneal biomechanics and their association with severity of lens dislocation in Marfan syndrome

PURPOSE

To investigate corneal biomechanical properties and its associations with the severity of lens dislocation in patients with Marfan syndrome.

METHODS

A total of 30 patients with Marfan syndrome and 30 age-, sex- and axial length (AL)-matched controls were recruited. Corneal biomechanical parameters of both groups were measured by CorVis ST and were compared between groups. Potential associations between corneal biomechanical parameters and severity of lens dislocation were also investigated.

RESULTS

Lower applanation 1 velocity (A1V) (0.13 ± 0.004 vs. 0.15 ± 0.003, P = 0.016), shorter applanation 2 time (A2T)(22.64 ± 0.11 vs. 22.94 ± 0.11, P = 0.013), longer peak distance (PD) (5.03 ± 0.07 vs. 4.81 ± 0.05, P = 0.008), longer radius (R) of highest concavity (7.44 ± 0.16 vs. 6.93 ± 0.14, P = 0.012), greater Ambrosio relational thickness horizontal (ARTh) (603 ± 20 vs. 498 ± 12, P < 0.001), and integrated radius (IR) (8.32 ± 0.25 vs. 8.95 ± 0.21, P = 0.033) were detected among Marfan eyes compared with controls (all P < 0.05). Marfan individuals with more severe lens dislocation tended to have increased stiffness parameter as longer A1T, slower A1V, shorter A2T, slower application 2 velocity (A2V), smaller PD and smaller Distance Amplitude (DA) (P < 0.05).

CONCLUSION

Marfan patients were detected to have increased corneal stiffness compared with normal subjects. Corneal biomechanical parameters were significantly associated with the severity of lens dislocation in Marfan patients.

Corneal biomechanics and diagnostics: a review

PURPOSE

Corneal biomechanics is an emerging field and the interest into physical and biological interrelations in the anterior part of the eye has significantly increased during the past years. There are many factors that determine corneal biomechanics such as hormonal fluctuations, hydration and environmental factors. Other factors that can affect the corneas are the age, the intraocular pressure and the central corneal thickness. The purpose of this review is to evaluate the factors affecting corneal biomechanics and the recent advancements in non-destructive, in vivo measurement techniques for early detection and improved management of corneal diseases.

METHODS

Until recently, corneal biomechanics could not be directly assessed in humans and were instead inferred from geometrical cornea analysis and ex vivo biomechanical testing. The current research has made strides in studying and creating non-destructive and contactless techniques to measure the biomechanical properties of the cornea in vivo.

RESULTS

Research has indicated that altered corneal biomechanics contribute to diseases such as keratoconus and glaucoma. The identification of pathological corneas through the new measurement techniques is imperative for preventing postoperative complications.

CONCLUSIONS

Identification of pathological corneas is crucial for the prevention of postoperative complications. Therefore, a better understanding of corneal biomechanics will lead to earlier diagnosis of ectatic disorders, improve current refractive surgeries and allow for a better postoperative treatment.

Accuracy of an Air-Puff Dynamic Tonometry Biomarker to Discriminate the Corneal Biomechanical Response in Patients With Keratoconus

PURPOSE

The aim of this study was to assess accuracy of the mean corneal stiffness ( kc , N/m) parameter to discriminate between patients with keratoconus and age-matched healthy subjects.

METHODS

Dynamic Scheimpflug imaging tonometry was performed with Corvis ST (Oculus Optikgeräte GmbH, Germany) in patients with keratoconus (n = 24; study group) and age-matched healthy subjects (n = 32; control). An image processing algorithm was developed to analyze the video sequence of the Corvis ST air-puff event and to determine the geometric and temporal parameters that correlated with the corneal tissue biomechanical properties. A modified 3-element viscoelastic model was used to derive the kc parameter, which represented the corneal tissue resistance to deformation under load. Receiver operating characteristic curves were used to assess the overall diagnostic performance for determining the area under the curve, sensitivity, and specificity of the kc in assessing the corneal tissue deformation to the Corvis ST air-puff event in keratoconus and control eyes. The Corvis Biomechanical Index ( CBI ) was analyzed for external validation.

RESULTS

The kc parameter was significantly different between keratoconus and controls ( P < 0.001), ranging from 24.9 ±3.0 to 34.2 ±3.5 N/m, respectively. It was highly correlated with CBI (r = -0.69; P < 0.001); however, the kc parameter had greater specificity (94%) than CBI (75%), whereas the 2 biomarkers had similar area under the curve (0.98 vs. 0.94) and sensitivity (96% vs. 92%) in predicting the occurrence of keratoconus.

CONCLUSIONS

The kc parameter extracted by video processing analysis of dynamic Scheimpflug tonometry data was highly accurate in discriminating patients with clinically manifest keratoconus compared with controls.

Corneal retardation time as an ocular hypertension disease indicator

Objective.A detailed analysis of the corneal retardation timeτas a highly related parameter to the intraocular pressure (IOP), and its plausible role as an indicator of ocular hypertension disease.Approach.A simple theoretical expression forτis derived within the corneal viscoelastic model of Kelvin-Voigt with 3 elements. This retardation time can be easily calculated from the well-known signal and pressure amplitudes of non-contact tonometers like the Ocular Response Analyzer (ORA). Then, a population-based study was performed where 100 subjects aged from 18 to 30 were analyzed (within this group, about 10% had an elevated IOP with more than 21 mmHg).Main results.A clear relationship between the corneal retardation time and the corneal-compensated intraocular pressure (IOPcc) was found, underlying the risk for ocular hypertensive (OHT) subjects with lowerτvalues to develop hypertension illnesses (due to the inability of poorly viscoelastic corneas to absorb IOP fluctuations, resulting in probable optic nerve damage).Significance.Our results might provide an useful tool to systematically discern which OHT patients (and even those with normal IOP values) are more likely to suffer glaucoma progression and, consequently, ensure an early diagnosis.

Localized Refractive Changes Induced by Symmetric and Progressive Asymmetric Intracorneal Ring Segments Assessed with a 3D Finite-Element Model

To build a representative 3D finite element model (FEM) for intracorneal ring segment (ICRS) implantation and to investigate localized optical changes induced by different ICRS geometries, a hyperelastic shell FEM was developed to compare the effect of symmetric and progressive asymmetric ICRS designs in a generic healthy and asymmetric keratoconic (KC) cornea. The resulting deformed geometry was assessed in terms of average curvature via a biconic fit, sagittal curvature (K), and optical aberrations via Zernike polynomials. The sagittal curvature map showed a locally restricted flattening interior to the ring (Kmax -11 to -25 dpt) and, in the KC cornea, an additional local steepening on the opposite half of the cornea (Kmax up to +1.9 dpt). Considering the optical aberrations present in the model of the KC cornea, the progressive ICRS corrected vertical coma (-3.42 vs. -3.13 µm); horizontal coma (-0.67 vs. 0.36 µm); and defocus (2.90 vs. 2.75 µm), oblique trefoil (-0.54 vs. -0.08 µm), and oblique secondary astigmatism (0.48 vs. -0.09 µm) aberrations stronger than the symmetric ICRS. Customized ICRS designs inspired by the underlying KC phenotype have the potential to achieve more tailored refractive corrections, particularly in asymmetric keratoconus patterns.

Human cornea thermo-viscoelastic behavior modelling using standard linear solid model

BACKGROUND

Corneal biomechanics is of great interest to researchers recently. Clinical findings relate them to corneal diseases and to outcomes of refractive surgery. To have a solid understanding of corneal diseases' progression, it is important to understand corneal biomechanics. Also, they are essential for better explaining outcomes of refractive surgeries and their undesired consequences. There is a difficulty for studying corneal biomechanics in-vivo and multiple limitations arise for ex-vivo studies. Hence mathematical modelling is considered as a proper solution to overcome such obstacles. Mathematical modelling of cornea in-vivo allows studying corneal viscoelasticity with taking into consideration all boundary conditions existing in real in-vivo situation.

METHODS

Three mathematical models are used to simulate corneal viscoelasticity and thermal behavior in two different loading situations: constant and transient loading. Two models of the three are used for viscoelasticity simulation which are Kelvin-Voigt and standard linear solid models. Also, temperature rise due to the ultrasound pressure push is calculated using bioheat transfer model for both the axial direction and as a 2D spatial map using the third model (standard linear solid model).

RESULTS

Viscoelasticity simulation results show that standard linear solid model is efficient for describing the viscoelastic behavior of human cornea in both loading conditions. Results show also that the deformation amplitude obtained from standard linear solid model is more reasonable for corneal soft-tissue deformation with respect to corresponding clinical findings than that obtained from Kelvin-Voigt model. Thermal behavior results estimated corneal temperature rise to be roughly 0.2 °C, which conforms with FDA regulations for soft tissue safety.

CONCLUSION

Standard Linear Solid (SLS) model is better describing the human corneal behavior in response to constant and transient load more efficiently. Temperature rise (TR) for the corneal tissue of about 0.2 °C is conforming with FDA regulations and even less than the FDA regulations for soft tissue safety.

Predicting corneal cross-linking treatment efficacy with real-time assessment of corneal riboflavin concentration

PURPOSE

To assess predictability of tissue biomechanical stiffening induced by UV-A light-mediated real-time assessment of riboflavin concentration during corneal crosslinking (CXL) of human donor tissues.

SETTING

Studio Italiano di Oftalmologia, Rome, Italy.

DESIGN

Laboratory study.

METHODS

20 sclerocorneal tissues were randomly stratified to undergo CXL with either the epithelium intact (n = 12) or removed (n = 8). Samples underwent corneal soaking with 0.22% riboflavin formulation (RitSight) with dosing time of t = 10 minutes and t = 20 minutes in epithelium-off and epithelium-on protocols, respectively. All tissues underwent 9-minute UV-A irradiance at 10 mW/cm 2 using theranostic device (C4V CHROMO4VIS). The device used controlled UV-A light irradiation to induce both imaging and treatment of the cornea, providing a real-time measure of corneal riboflavin concentration and treatment efficacy (ie, theranostic score) during surgery. Tissue biomechanics were assessed with an air-puff device (Corvis), which was performed before and after treatment. A 3-element viscoelastic model was developed to fit the corneal deformation response to air-puff excitation and to calculate the mean corneal stiffness parameter (k c ).

RESULTS

Significant corneal tissue stiffening ( P < .05) was induced by the theranostic UV-A device in either CXL treatment protocol. Significant correlation was found between the theranostic score and the increase in k c ( R = 0.75; P = .003). The score showed high accuracy (94%) and precision (94%) to predict correctly samples that had improved tissue biomechanical strengthening.

CONCLUSIONS

Real-time assessment of corneal riboflavin concentration provided a predictive and precise approach for significant improvement of tissue strength on individual corneas, regardless of CXL treatment protocol.

In vivo corneal elastography: A topical review of challenges and opportunities

Clinical measurement of corneal biomechanics can aid in the early diagnosis, progression tracking, and treatment evaluation of ocular diseases. Over the past two decades, interdisciplinary collaborations between investigators in optical engineering, analytical biomechanical modeling, and clinical research has expanded our knowledge of corneal biomechanics. These advances have led to innovations in testing methods (ex vivo, and recently, in vivo) across multiple spatial and strain scales. However, in vivo measurement of corneal biomechanics remains a long-standing challenge and is currently an active area of research. Here, we review the existing and emerging approaches for in vivo corneal biomechanics evaluation, which include corneal applanation methods, such as ocular response analyzer (ORA) and corneal visualization Scheimpflug technology (Corvis ST), Brillouin microscopy, and elastography methods, and the emerging field of optical coherence elastography (OCE). We describe the fundamental concepts, analytical methods, and current clinical status for each of these methods. Finally, we discuss open questions for the current state of in vivo biomechanics assessment techniques and requirements for wider use that will further broaden our understanding of corneal biomechanics for the detection and management of ocular diseases, and improve the safety and efficacy of future clinical practice.

Biaxial hyperelastic and anisotropic behaviors of the corneal anterior central stroma along the preferential fibril orientations. Part II: Quantitative computational analysis of mechanical response of stromal components

To study the hyperelastic and anisotropic behaviors of the central anterior stroma for patients with myopia, 40 corneal stromal specimens extracted after small incision lenticule extraction (SMILE) surgery were used in the biaxial extension test along two preferential fibril orientations. An improved collagen fibril crimping constitutive model with a specific physical meaning was proposed to analyze the hyperelasticity and anisotropy of the stroma. The effective elastic modulus of the two families of preferentially oriented collagen fibrils and the stiffness of the non-collagenous matrix along all three directions were compared according to the specific physical meaning of the parameters. Anisotropic behavior was found in the hyperelastic properties of the corneal anterior central stroma in the preferential fibril orientations. The stiffness of non-collagenous matrix is significantly larger in the optical axis direction than in the nasal-temporal (NT) and superior-inferior (SI) directions. Moreover, individual differences between males and females slightly impact on hyperelastic and anisotropic behaviors. The differences of these behaviors were significant in the comparison of the left and right eyes. These results have a guiding significance for the accurate design of surgical plans for refractive surgery according to a patient's condition and have a driving value for the further exploration of the biomechanical properties of the whole cornea.

Clinical Ocular Biomechanics: Where Are We after 20 Years of Progress?

Purpose: Ocular biomechanics is an assessment of the response of the structures of the eye to forces that may lead to disease development and progression, or influence the response to surgical intervention. The goals of this review are (1) to introduce basic biomechanical principles and terminology, (2) to provide perspective on the progress made in the clinical study and assessment of ocular biomechanics, and (3) to highlight critical studies conducted in keratoconus, laser refractive surgery, and glaucoma in order to aid interpretation of biomechanical parameters in the laboratory and in the clinic.Methods: A literature review was first conducted of basic biomechanical studies related to ocular tissue. The subsequent review of ocular biomechanical studies was limited to those focusing on keratoconus, laser refractive surgery, or glaucoma using the only two commercially available devices that allow rapid assessment of biomechanical response in the clinic.Results: Foundational studies on ocular biomechanics used a combination of computer modeling and destructive forces on ex-vivo tissues. The knowledge gained from these studies could not be directly translated to clinical research and practice until the introduction of non-contact tonometers that quantified the deformation response of the cornea to an air puff, which represents a non-destructive, clinically appropriate load. The corneal response includes a contribution from the sclera which may limit corneal deformation. Two commercial devices are available, the Ocular Response Analyzer which produces viscoelastic parameters with a customized load for each eye, and the Corvis ST which produces elastic parameters with a consistent load for every eye. Neither device produces the classic biomechanical properties reported in basic studies, but rather biomechanical deformation response parameters which require careful interpretation.Conclusions: Research using clinical tools has enriched our understanding of how ocular disease alters ocular biomechanics, as well as how ocular biomechanics may influence the pathophysiology of ocular disease and response to surgical intervention.

Repeatability and correlation of corneal biomechanical measurements obtained by Corvis ST in orthokeratology patients

PURPOSE

To evaluate the repeatability of the corneal biomechanical measurements obtained by Corvis ST in post-orthokeratology patients and analyze the correlation between the biomechanical and ocular parameters.

METHODS

Fifty-one eyes of 51 myopic subjects were included in this study. The biomechanical parameters were assessed using Corvis ST. Repeatability was assessed using one-way ANOVA based on within-subject standard deviation (Sw), repeatability coefficient (RC), intraclass correlation coefficient (ICC) and correlation of variation (CoV). The correlation was evaluated using Pearson correlation analysis.

RESULTS

All parameters measured by Corvis ST, except length of flattened cornea at the first and second applanations (A1L and A2L), showed a good intraobserver repeatability after a 3-month follow-up period. The ICC values for A1L and A2L were 0.444 and 0.654, whereas the other parameters were higher than 0.8. Similar trends were obtained for CoV, wherein the CoV values for A1L and A2L were greater than 13 %. The corneal biomechanical parameters were correlated with age, refraction, axial length (AL), steep and flat keratometry before and after orthokeratology, and central corneal thickness (CCT). Following orthokeratology treatment, post-keratometry demonstrated a higher correlation with stiffness parameter at first applanation (SP-A1), velocity of corneal apex at the first applanation (A1V), and radius than pre-keratometry, which showed a weak correlation with SP-A1.

CONCLUSION

Corneal biomechanical parameters assessed using Corvis ST demonstrated a good repeatability, except A1L and A2L. The corneal biomechanical parameters were correlated with age, refraction, AL and pre- and post-keratometry. Thus, Corvis ST is a suitable device for investigating biomechanical parameter.

Theranostic-guided corneal cross-linking: Preclinical evidence on a new treatment paradigm for keratoconus

Theranostics is an emerging therapeutic paradigm of personalized medicine; the term refers to the simultaneous integration of therapy and diagnostics. In this work, theranostic-guided corneal cross-linking was performed on 10 human sclero-corneal tissues. The samples were soaked with 0.22% riboflavin formulation and underwent 9 minutes UV-A irradiance at 10 mW/cm² using theranostic device, which provided both a measure of corneal riboflavin concentration and a theranostic score estimating treatment efficacy in real time. A three-element viscoelastic model was developed to fit the deformation response of the cornea to air-puff excitation of dynamic tonometry and to calculate the mean corneal stiffness parameter before and after treatment. Significant correlation was found between the theranostic score and the increase in mean corneal stiffness (R = 0.80; P < .001). Accuracy and precision of the theranostic score in predicting the induced corneal tissue stiffening were both 90%. The riboflavin concentration prior to starting the UV-A photo-therapy phase was the most important variable to allow corneal cross-linking to be effective. Theranostic UV-A light mediated imaging and therapy enables the operator to adopt a precise approach for achieving highly predictable biomechanical strengthening on individual corneas.

In Vivo Determination of the Human Corneal Elastic Modulus Using Vibrational Optical Coherence Tomography

Purpose

To determine the in vivo elastic modulus of the human cornea using vibrational optical coherence tomography (VOCT).

Methods

Vibrational analysis coupled with optical coherence tomography (OCT) was used to obtain the resonant frequency (RF) and elastic modulus of corneal structural components. VOCT corneal thickness values were measured using OCT images and correlated with corneal thickness determined with Pentacam (Oculus, Wetzlar, Germany). Moduli were obtained at two locations: central cornea (CC) and inferior cornea (IC). Measurements were obtained with and without anesthetic eye drops to assess their effect on the modulus measurements.

Results

VOCT thickness values correlated positively (R² = 0.97) and linearly (y = 1.039x–16.89) with those of Pentacam. Five RF peaks (1–5) were present, although their presence was variable across eyes. The RF for peaks 1 to 5 in the CC and IC ranged from 73.5 ± 4.9 to 239 ± 3 Hz and 72.1 ± 6.3 to 238 ± 4 Hz, respectively. CC and IC moduli for peaks 1 to 5 ranged from 1.023 ± 0.104 to 6.87 ± 0.33 MPa and 0.98 ± 0.15 to 6.52 ± 0.79 MPa, respectively. Topical anesthesia did not significantly alter the modulus (P > 0.05 for all), except for peak 2 in the CC (P < 0.05).

Conclusions

This pilot study demonstrates the utility of VOCT as an in vivo, noninvasive technology to measure the elastic modulus in human corneas. The structural origin of these moduli is hypothesized based on previous reports, and further analyses are necessary for confirmation.

Translational Relevance

This work presents VOCT as a novel approach to assess the in vivo elastic modulus of the cornea, an indicator of corneal structural integrity and health.

Determine Corneal Biomechanical Parameters by Finite Element Simulation and Parametric Analysis Based on ORA Measurements

Purpose: The Ocular Response Analyzer (ORA) is one of the most commonly used devices to measure corneal biomechanics in vivo. Until now, the relationship between the output parameters and corneal typical biomechanical parameters was not clear. Hence, we defined the output parameters of ORA as ORA output parameters. This study aims to propose a method to determine corneal biomechanical parameters based on ORA measurements by finite element simulation and parametric analysis.

Methods: Finite element analysis was used to simulate the mechanics process of ORA measurements with different intraocular pressure (IOP), corneal geometrical parameters and corneal biomechanical parameters. A simplified geometrical optics model was built to simulate the optical process of the measurements to extract ORA output parameters. After that, 70% of the simulated data was used to establish the quantitative relationship between corneal biomechanical parameters and ORA output parameters by parametric analysis and 30% of the simulated data was used to validate the established model. Besides, ten normal subjects were included to evaluate the normal range of corneal biomechanical parameters calculated from ORA.

Results: The quantitative relationship between corneal biomechanical parameters and ORA output parameters is established by combining parametric analysis with finite element simulation. The elastic modulus (E) and relaxation limit (G_∞) of the ten normal subjects were 0.65 ± 0.07 MPa and 0.26 ± 0.15, respectively.

Conclusions: A method was proposed to determine corneal biomechanical parameters based on the results of ORA measurements. The magnitude of the corneal biomechanical parameters calculated according to our method was reasonable.

Risk Factors for Corneal Striae in Eyes After Glaucoma Surgery

PRCIS

Eyes with corneal striae had steeper cornea, induced astigmatism, and higher corneal hysteresis (CH), which implies a relationship between striae, corneal shape, and the cornea's resistance to deformation at low intraocular pressures (IOPs).

BACKGROUND

Anterior corneal striae (ACS) are associated with low IOP. However, the clinical significance of ACS is unclear. Here, we aim to evaluate differences in eyes with striae compared with eyes without striae.

METHODS

Adults with ACS (cases) and without ACS (controls) ≥8 weeks after glaucoma surgery with an IOP ≤10 mm Hg were enrolled. Optical coherence tomography and optical biometry were performed. CH, defined as the difference in pressure between corneal indentation and reformation in response to an air jet, was obtained by the ocular response analyzer. Hypotony maculopathy (HM) was defined as optic disc swelling, vascular tortuosity attributed to hypotony, or clinical presence of chorioretinal folds confirmed on OCT.

RESULTS

One hundred sixteen eyes (76 cases, 40 controls) were included. Cases had a lower IOP compared with controls (6.5±2.3 vs. 8.5±1, P<0.0001). A 1 mm Hg increase in CH increased ACS odds [odds ratio (OR)=1.51, P=0.01]. A 1 D increase in the flattest presurgical and postsurgical corneal power increased ACS odds by 1.83 (P=0.01) and 1.41 (P=0.02), respectively. Astigmatism increased in eyes with ACS by 1.11 D (P<0.001). ACS odds were increased with every 1 minute increase in mitomycin-C duration (OR=1.58, P=0.047) and decreased with the use of topical glaucoma medication (OR=0.62, P=0.03). Visual acuity decreased from logarithm of the minimal angle of resolution 0.22 (20/33 Snellen) presurgery to 0.28 (20/38) postsurgery (P=0.008), independent of ACS. HM occurred in 19% of cases (P=0.05). A higher postsurgical CH increased HM odds (OR=1.8, P=0.003). HM predicted a 0.41 mm decrease in axial length (P<0.0001), independent of IOP.

CONCLUSION

ACS were associated with a steeper cornea, induced astigmatism, and higher CH, suggesting a relationship between striae, corneal shape, and the cornea's ability to resist deformations at lower IOP. CH, HM, and axial length shortening were associated independently of IOP.

Viscoelastic properties' characterization of corneal stromal models using non-contact surface acoustic wave optical coherence elastography (SAW-OCE)

Viscoelastic characterization of the tissue-engineered corneal stromal model is important for our understanding of the cell behaviors in the pathophysiologic altered corneal extracellular matrix (ECM). The effects of the interactions between stromal cells and different ECM characteristics on the viscoelastic properties during an 11-day culture period were explored. Collagen-based hydrogels seeded with keratocytes were used to replicate human corneal stroma. Keratocytes were seeded at 8 × 10³ cells per hydrogel and with collagen concentrations of 3, 5 and 7 mg/ml. Air-pulse-based surface acoustic wave optical coherence elastography (SAW-OCE) was employed to monitor the changes in the hydrogels' dimensions and viscoelasticity over the culture period. The results showed the elastic modulus increased by 111%, 56% and 6%, and viscosity increased by 357%, 210% and 25% in the 3, 5 and 7 mg/ml hydrogels, respectively. To explain the SAW-OCE results, scanning electron microscope was also performed. The results confirmed the increase in elastic modulus and viscosity of the hydrogels, respectively, arose from increased fiber density and force-dependent unbinding of bonds between collagen fibers. This study reveals the influence of cell-matrix interactions on the viscoelastic properties of corneal stromal models and can provide quantitative guidance for mechanobiological investigations which require collagen ECM with tuneable viscoelastic properties.

Refractive changes of a new asymmetric intracorneal ring segment with variable thickness and base width: A 2D finite-element model

Purpose

To evaluate the local geometric effects of a unilateral intrastromal ring segment with a combined variation of ring thickness and base width in a finite element simulation, and to compare it against the isolated effect of thickness or base width variation alone.

Methods

A two-dimensional finite-element model of a transversely isotropic cornea was created assuming either axisymmetric stress or plane strain condition. The model geometry was composed of a three-layered corneal tissue (epithelium, anterior and posterior stroma) fixed at the limbus. The implantation of a triangular-shape asymmetric ring segment with varying ring thickness (150 to 300 μm) and base width (600 to 800 μm) was simulated. Also, changes induced by thickness or base width alone were studied and compared their combined effect in the asymmetric ring segment. Geometrical deformation of the simulated cornea and sagittal curvature were the main parameters of study.

Results

Increasing ring thickness and base width along the arc of the asymmetric ring segment produced a more pronounced flattening in this part of the ring. The asymmetric design did find a good balance between maximizing corneal flattening at one end and minimizing it at the other end, compared to the isolated effect of ring thickness and width. Ring thickness was the most robust parameter in flattening both, the central and peripheral cornea.

Conclusion

The finite-element model permitted a theoretical study of corneal deformation undergoing implantation of realistic and hypothetical ring geometries. Intracorneal asymmetric ring segments with varying thickness and base width can be a good alternative in corneas with asymmetric keratoconus phenotypes.

The effect of corneal crosslinking on the rigidity of the cornea estimated using a modified algorithm for the Schiøtz tonometer

Purpose:

The aim of this study was to test a method for estimating corneal rigidity before and after cross-linking (CXL) using a Schiøtz tonometer.

Methods:

The study was performed in the Kyiv City Clinical Ophthalmological Hospital “Eye Microsurgical Center”, Ukraine. This was a prospective, consecutive, randomized, masked, case-by-case, clinical study. Corneal rigidity, indicated by the gradient (G) between lg applied weight and corresponding lg scale reading during Schiøtz tonometry, were obtained by increasing (A-mode) then reducing (D-mode) weights by two operators [A] in keratoconus, post-CXL and control subjects for estimation of (i) interoperator and (ii) intersessional errors, (iii) intergroup differences; [B] before and after CXL. Central corneal thickness CCT was measured by scanning slit pachymetry. ANOVA, t tests, linear regression were the statistical tools used.

Results:

Average interoperator difference (ΔG) was –0.120 (SD = ±0.294, 95%CI = –0.175 to –0.066). A significant correlation between ΔG and the mean of each pair of G values was found (r = –0.196, n = 112, P = 0.038). Intersessional differences in mean G values were insignificant (P > 0.05). There was a significant correlation between G at first session (X₁) and difference between sessions (ΔG) [Operator 1, ΔG = 0.598x₁–0.461, r = 0.601, n = 27, P = 0.009]. Significant intergroup differences in G were found (Operator 1, one-way ANOVA, F = 4.489, P = 0.014). The difference (Δ) between the pre-(X₂) and post-CXL treatment G values was significantly associated with the pre-CXL treatment value (Operator 1, Δ = 1.970x₂-1.622, r = 0.642, n = 18, P = <.001). G values were correlated with CCT in keratoconus and post-CXL.

Conclusion:

Corneal rigidity (G) estimated using the Schiøtz tonometer can be useful for detecting changes after CXL. However, G values are linked to CCT, can vary from time-to-time and the procedure is operator dependent.

Corneal Biomechanical Assessment with Ultra-High-Speed Scheimpflug Imaging During Non-Contact Tonometry: A Prospective Review

BACKGROUND

In recent years, increasing interest has arisen in the application of data from corneal biomechanics in many areas of ophthalmology, particularly to assist in the detection of early corneal ectasia or ectasia susceptibility, to predict corneal response to surgical or therapeutic interventions and in glaucoma management. Technology has evolved and, recently, the Scheimpflug principle was associated with a non-contact air-puff tonometer, allowing a thorough analysis of corneal biomechanics and a biomechanically corrected intraocular pressure assessment, opening up new perspectives both in ophthalmology and in other medical areas. Data from corneal biomechanics assessment are being integrated in artificial intelligence models in order to increase its value in clinical practice.

OBJECTIVE

To review the state of the art in the field of corneal biomechanics assessment with special emphasis to the technology based on ultra-high-speed Scheimpflug imaging during non-contact tonometry.

SUMMARY

A meticulous literature review was performed until the present day. We used 136 published manuscripts as our references. Both information from healthy individuals and descriptions of possible associations with systemic diseases are described. Additionally, it exposed information regarding several fields of ocular pathology, from cornea and ocular surface through areas of refractive surgery and glaucoma until vascular and structural diseases of the chorioretinal unit.

Direct Evidence of Symmetry between Bilateral Human Corneas in Biomechanical Properties: A Comparison Study with Fresh Corneal Tissue

Purpose. To investigate the difference between the eyes from the same human with respect to the biomechanical properties of fresh corneal tissues and investigate the assumption of similarity of the corneal biomechanical properties between the eyes. Methods. Strip specimens extracted through a small incision lenticule extraction (SMILE) surgery were tested using a uniaxial tensile test. The specimens were extracted vertically. Low-strain tangent modulus (LSTM), high-strain tangent modulus (HSTM), and tensile strength (${\sigma }_b$) were the biomechanical parameters used in the comparison of the eyes from the same human. Results. Ninety corneal specimens from 45 persons were included in this study. The LSTM of the left and right eyes were 1.34 ± 0.52 and 1.37 ± 0.46 MPa, while the HSTM were 50.53 ± 7.51 and 49.41 ± 7.01 MPa, respectively. There was no significant difference between the eyes in terms of LSTM, HSTM, and${\sigma }_b\left(P=0.813,0.335,\text{and\hspace{0.17em}}0.605,\text{resp}.\right)$. The LSTM and HSTM were significantly correlated with the spherical equivalent (SE) ($P\le 0.01,P=0.001$, resp.). Conclusions. The assumption that the corneal biomechanical properties of the eyes from the same human are similar has been confirmed for the first time using fresh human corneal tissue. This finding may be useful in further biomechanical studies.

Dynamic Scheimpflug Ocular Biomechanical Parameters in Untreated Primary Open Angle Glaucoma Eyes

Purpose

To characterize the corneal biomechanical properties of glaucoma eyes by comparing the dynamic Scheimpflug biomechanical parameters between untreated glaucoma and control eyes.

Methods

Cross-sectional observational data of dynamic Scheimpflug analyzer (Corvis ST) examinations were retrospectively collected from 35 eyes of 35 consecutive patients with untreated normal tension glaucoma and 35 eyes of 35 healthy patients matched on age and IOP. Ten biomechanical parameters were compared between the two groups using multivariable models adjusting for IOP, central corneal thickness, age, and axial length. The Benjamini-Hochberg method was used to correct for multiple comparison.

Results

In multivariable models, glaucoma was associated with smaller applanation 1 time (P < 0.001, coefficient = −0.5865), applanation 2 time (P = 0.012, coefficient = −0.1702), radius (P = 0.006, coefficient = −0.5447), larger peak distance (P = 0.011, coefficient = 0.1023), deformation amplitude ratio at 1 mm (P < 0.001, coefficient = 0.072), and integrated radius (P < 0.001, coefficient = 1.094). These associations consistently indicate greater compliance of the cornea in glaucoma eyes.

Conclusions

Untreated normal tension glaucoma eyes were more compliant than healthy eyes. The greater compliance (smaller stiffness) of normal tension glaucoma eyes may increase the risk of optic nerve damage. These results suggest the relevance of measuring biomechanical properties of glaucoma eyes.

Plastic modification of human cornea in vivo: applications to clinical refractive procedures

Purpose:

The purpose of this report is to quantify how pressure applied to the human cornea, either physiological or intentional, affects its curvature. In particular, how pneumatic procedures flatten the central cornea and keep it flat over time, thereby decreasing the patient’s myopia.

Methods:

A viscoelastic model is developed for plastic deformation which gives us the basic governing equations of the elastic and plastic strain of corneal stroma. The model is applied to data from corneas of six patients who underwent pneumatic keratology (NEumatica Keratologia) to reduce their myopia.

Results:

The model shows corneal dimensional stability for long periods of time after NEumatica Keratologia that decay with an exponential time constant. Separate equations are developed that relate corneal plastic strain to the pressure applied and its duration ε = σ₀t₁/η₁, to change in refraction ε = 2 × ΔRefr, to keratometry radius increase ε = ΔR/R, and to corneal thinning ε = sqr (Δh/h). The average values obtained for ε from the patients’ data are 3%, 3.2%, 3%, and 2.6%, respectively, all in remarkable agreement. The average refraction change is found to remain stable at ΔRefr = +1.67D ± 5.2%. Clinical data yield good agreement of theory and treatment results.

Conclusions:

The model proposed is a good description of NEumatica Keratologia outcomes. Practical applications include the long-term stable correction of myopia with refractive procedures. High myopia subjects can benefit from this procedure because NEumatica Keratologia corrects and protects the central cornea radius by stretching the peripheral cornea.

Corneal Biomechanical Changes After Uneventful Phacoemulsification in Patients With Type 2 Diabetes Mellitus and Patients Without Diabetes

PURPOSE

To compare corneal biomechanical changes after uneventful phacoemulsification cataract surgery between type 2 diabetic (DM) and nondiabetic patients.

METHODS

Forty-four diabetic (44 eyes) and 44 (44 eyes) age and sex-matched non-DM controls with age-related cataract were enrolled in this prospective observational study. Corneal hysteresis (CH), corneal resistance factor (CRF), and corneal-compensated intraocular pressure (IOPcc) were evaluated by using the ocular response analyzer; central corneal thickness was recorded by using the Pentacam HR. Patients were evaluated preoperatively and 1 and 6 months after surgery.

RESULTS

In the DM group, CH was observed to be significantly lower than preoperative value (9.8 ± 1.5 mm Hg) at 1 month (9.4 ± 1.2 mm Hg, P = 0.040), but not at the 6-month follow-up (9.6 ± 1.6 mm Hg, P = 0.437),whereas it did not change significantly in the non-DM group (preoperative 9.8 ± 1.3 mm Hg vs. 1 month 9.6 ± 1.1 mm Hg vs. 6 months 10.1 ± 1.1 mm Hg, P > 0.05). CRF was significantly lower than the preoperative values at 1 and 6 months in both groups (P ≤ 0.001). Postoperative CH change was significantly associated with preoperative CH (P < 0.001), preoperative IOPcc (P = 0.004), and IOPcc change (P < 0.001), whereas CRF change was only correlated with preoperative CRF (P < 0.001). There was a significant postoperative IOPcc reduction 6 months after surgery (P < 0.001) in both DM and non-DM groups; however, central corneal thickness was not found to significantly change in the postoperative period (P > 0.05).

CONCLUSIONS

This study showed that phacoemulsification causes a significant and persistent decrease in intraocular pressure and CRF in both groups, whereas CH recovered to preoperative values, although more slowly in patients with diabetes.

Comparison of Corneal Biomechanical Properties among Axial Myopic, Nonaxial Myopic, and Nonmyopic Eyes

Purpose

To compare corneal deformation characteristics using ultra-high-speed Scheimpflug camera (Corvis ST) in patients with nonmyopic (NM), mild-to-moderate nonaxial myopic (MM), and high axial myopic (HM) eyes.

Methods

In this cross-sectional study, normal subjects aged >40 years with no history of ocular laser/surgery were classified according to axial length (AL) and spherical equivalence (SE) into three groups: (1) NM (SE > −0.50 D and AL < 26 mm), (2) MM (SE −6.00 D to −0.50 D and AL < 26 mm), and (3) HM (SE ≤ −6.00 D and AL ≥ 26 mm). Seven parameters including corneal deformation amplitude (CDA), inward/outward corneal applanation length, inward/outward corneal velocity (ICV and OCV), peak distance, and radius were measured. Pearson correlation and linear mixed-effects model were done.

Results

A total of 180 eyes were recruited. 98 eyes were NM, 30 eyes were MM, and 52 eyes were HM. There were significant correlations of OCV to the degree of refractive error (r = 0.203, p < 0.001) and AL (r = −0.242, p < 0.001). After adjusting for age, sex, intraocular pressure, and corneal thickness, there was significantly higher CDA (β = 0.07, p < 0.001), faster OCV (β = −0.08, p < 0.001), and smaller radius (β = −0.39, p=0.01) in the HM group compared to the NM group.

Conclusion

The higher CDA, faster OCV, and smaller radius found in the HM may suggest that these eyes have reduced ocular stiffness and may be less stable and more prone to stress.

The Association Between Keratoconus and Mitral Valve Prolapse: A Meta-Analysis

Objective

The debate pertaining to the association between Keratoconus (KC) and Mitral Valve Prolapse (MVP) continues to occur among physicians. The results of cross-sectional studies attempting to present the co-existing prevalence of these two diseases remain indeterminate. We compiled the first meta-analysis to determine the pattern of prevalence between the two diseases.

Methods

Two separate literature searches for cross-sectional studies were performed for this meta-analysis. The first search encompassed finding literature comparing the prevalence of KC between patients with MVP and a control group. The second search pertained to finding studies comparing the prevalence of MVP patients with KC and a control group.

Results

Six studies reported the prevalence of MVP in patients with KC and a control group. The prevalence was 41.6% in patients with KC and 11.5% in patients without KC (OR = 7.06 [95% CI = 2.41-20.64]). There was a significant heterogeneity among the studies (I² = 84%). Two studies showed the prevalence of KC in patients with MVP and a control group. The prevalence was 17.0% in patients with KC and 2.9% in the control group (OR = 5.07 [95% CI = 1.08-23.83]). There was no heterogeneity within the analysis (I² = 0%).

Conclusion

There is a statistically significant co-existing prevalence between MVP and KC. Patients with KC are more likely to present with MVP, and patients with MVP are more likely to present with KC.

Effect of congenital blepharoptosis on corneal biomechanical properties and changes after ptosis surgery

OBJECTIVES

We studied the difference in the corneal biomechanical parameters of ptotic and fellow eyes in patients with congenital blepharoptosis. The correlations between corneal biomechanical parameters and demographic or ocular parameters, and the changes after surgery were also researched.

METHODS

The corneal biomechanical parameters were measured by Corvis ST tonometry. The central corneal thickness (CCT), axial length (AL) and keratometry measurements were performed with LenStar LS900, and intraocular pressure (IOP) by non-contact applanation tonometry. The parameters were evaluated for the effect of ptosis and the relationship of corneal biomechanical parameters. These examinations were repeated 6 months after blepharoptosis surgery.

RESULTS

Twenty-nine patients were enroled. The Corvis ST parameters (Deformation amplitude [DA], A1 times, and A1 velocity), CCT, IOP with NCT, IOP with corrected, differed significantly between ptotic and fellow eyes. CCT was significantly positively correlated with Length A1 and IOP with Corvis, and negatively correlated with IOP corrected by Corvis of the ptotic eyes. The same tendency was found in the fellow eyes. Six months after the ptosis surgery, the differences in corneal biomechanics parameters between ptotic eyes and fellow eyes were not significantly changed.

CONCLUSIONS

Congenital blepharoptosis causes significant corneal biomechanical changes measured by Corvis ST. The ptotic eyes had thicker and less deformable corneas. The differences in corneal biomechanics between ptotic eyes and fellow eyes were mostly related to CCT changes. Six months after surgery, these differences in corneal biomechanics did not change significantly.

Corneal biomechanical properties after SMILE versus FLEX, LASIK, LASEK, or PRK: a systematic review and meta-analysis

BACKGROUND

The aim of this study was to compare the postoperative corneal biomechanical properties between small incision lenticule extraction (SMILE) and other corneal refractive surgeries.

METHODS

A systematic review and meta-analysis were conducted. Articles from January 2005, to April 2019, were identified searching PubMed, EMBASE, Web of Science, and International Clinical Trials Registry Platform. Studies that compared SMILE with other corneal refractive surgeries on adult myopia patients and evaluated corneal biomechanics were included. Multiple effect sizes in each study were combined. Random-effects model was conducted in the meta-analysis.

RESULTS

Twenty-two studies were included: 5 randomized controlled trials (RCTs), 9 prospective and 6 retrospective cohort studies, and 2 cross-sectional studies. Using the combined effect of corneal hysteresis (CH) and corneal resistance factor (CRF), which were obtained from ocular response analyzer (ORA), the pooled Hedges' g of SMILE versus femtosecond laser-assisted in situ keratomileusis (FS-LASIK) was 0.41 (95% CI, 0.00 to 0.81; p = 0.049; I² = 78%), versus LASIK was 1.31 (95% CI, 0.54 to 2.08; p < 0.001; I² = 77%), versus femtosecond lenticule extraction (FLEX) was - 0.01 (95% CI, - 0.31 to 0.30; p = 0.972; I² = 20%), and versus the group of photorefractive keratectomy (PRK) and laser-assisted sub-epithelial keratectomy (LASEK) was - 0.26 (95% CI, - 0.67 to 0.16; p = 0.230; I² = 54%). The summary score of Corvis ST (CST) after SMILE was comparable to FS-LASIK/LASIK with the pooled Hedges' g = - 0.05 (95% CI, - 0.24 to 0.14; p = 0.612, I² = 55%).

CONCLUSIONS

In terms of preserving corneal biomechanical strength after surgeries, SMILE was superior to either FS-LASIK or LASIK, while comparable to FLEX or PRK/LASEK group based on the results from ORA. More studies are needed to apply CST on evaluating corneal biomechanics after refractive surgeries.

Structural and Biomechanical Corneal Differences between Type 2 Diabetic and Nondiabetic Patients

Purpose

To analyze and compare corneal structural and biomechanical properties, characterized by corneal hysteresis (CH) and resistance factor (CRF), between patients with and without type 2 diabetes mellitus (DM), and determine the main ocular variables that influence them.

Methods

Sixty diabetic and 48 age- and sex-matched non-DM patients were enrolled in this cross-sectional study. The DM group was analyzed according to DM duration (<or ≥ 10 years), HbA1c levels (<or ≥ 7%), and presence of retinopathy. CH and CRF were evaluated using the Ocular Response Analyzer® (ORA). Central corneal thickness (CCT) was determined by Scheimpflug tomography (Pentacam® HR). Intraocular pressure was obtained with ORA (IOPcc) and Goldmann applanation tonometry (IOP-GAT). Univariate and multivariate linear regression analyses were performed to evaluate the relationship between demographical, clinical, and ocular variables with the biomechanical properties.

Results

There were no statistically significant differences in the CH and the CRF between DM and non-DM groups (p=0.637 and p=0.439, respectively). Also, there was no statistical difference between groups for the CCT, IOPcc, or IOP-GAT. Multivariate linear regression analysis showed that CH was positively associated with CCT (p < 0.001) and negatively associated with IOPcc (p < 0.001), while CRF was positively associated with CCT (p < 0.001) and IOPcc (p=0.014).

Conclusion

The CCT and IOPcc were found to be the main parameters that affect corneal biomechanical properties both in diabetics and controls. In this study, there was no significant effect of DM type 2 on corneal biomechanics.

Evaluation of corneal elastic modulus based on Corneal Visualization Scheimpflug Technology

Background

Corneal biomechanical properties are important for the diagnosis of corneal diseases, individualized design and prognosis of corneal surgery. Clinical available devices such as Ocular Response Analyzer (ORA) and Corneal Visualization Scheimpflug Technology (Corvis ST) can provide corneal biomechanics related parameters, while corneal elastic modulus cannot be extracted directly from them at present. The aim of this study is to suggest a method to determine corneal elastic modulus based on the results of Corvis ST test according to Reissner’s theory on the relation between stress and small displacement in shallow spherical shell.

Results

Five rabbits (10 eyes) and 10 healthy humans (20 eyes) were measured with Corvis ST to obtain the normal range of corneal elastic modulus. Results showed Corneal elastic modulus of rabbit was 0.16 MPa to 0.35 MPa, human corneal elastic modulus was 0.16–0.30 MPa. Rabbit corneas were also measured at different intraocular pressures (IOP), and results showed corneal elastic modulus, first applanation time (A1T) and stiffness parameter (SP-A1) were positively correlated with IOP. Deformation amplitude (DA), the second applanations time (A2T), and peak distance (PD) were negatively correlated with IOP. Finite element method was used to simulate the Corvis measurements according to the calculated elastic modulus and the simulated corneal apical displacements were agreement with experimental results in general.

Conclusions

The method to determine corneal elastic modulus based on Corvis test according to the relationship between force and displacements of shallow spherical shell is convenient and effective.

A Review of Structural and Biomechanical Changes in the Cornea in Aging, Disease, and Photochemical Crosslinking

The study of corneal biomechanics is motivated by the tight relationship between biomechanical properties and visual function within the ocular system. For instance, variation in collagen fibril alignment and non-enzymatic crosslinks rank high among structural factors which give rise to the cornea's particular shape and ability to properly focus light. Gradation in these and other factors engender biomechanical changes which can be quantified by a wide variety of techniques. This review summarizes what is known about both the changes in corneal structure and associated changes in corneal biomechanical properties in aging, keratoconic, and photochemically crosslinked corneas. In addition, methods for measuring corneal biomechanics are discussed and the topics are related to both clinical studies and biomechanical modeling simulations.

Corneal Biomechanical Properties of Keratoconic Eyes Following Penetrating Keratoplasty

Objectives

To investigate the corneal biomechanical properties of keratoconic eyes following penetrating keratoplasty (PKP).

Materials and Methods

Thirty-five patients (70 eyes) were enrolled to this prospective study. Operated and contralateral keratoconic eyes were defined as Group 1 and 2, respectively. All patients underwent ophthalmological examination and measurements of corneal biomechanical properties by Ocular Response Analyzer (ORA), intraocular pressure (IOP) by Goldmann applanation tonometry, and central corneal thickness (CCT) by Pentacam. Shapiro-Wilk W test was performed to test normality of the data. The statistical significance was evaluated with the paired t-test and Wilcoxon signed ranks test. Pearson correlation and Spearman rho tests were used for correlation analysis.

Results

The average age and male/female ratio were 31.34±11.65 (15-60) years and 21/14, respectively. The mean values of the data obtained from Group 1 and 2 respectively were: corneal hysteresis (CH): 9.35±1.66, 8.18±1.84 mmHg (p=0.013), corneal resistance factor (CRF): 9.48±1.96, 7.14±2.05 mmHg (p<0.001), IOPcc: 16.90±4.32, 14.26±3.69 mmHg (p=0.004), IOPg: 15.45±4.61, 10.91±3.97 mmHg (p<0.001), IOPapl: 14.26±3.11, 13.09±2.54 mmHg (p=0.046), and central corneal thickness (CCT): 545.64±60.82, 442.60±68.14 μM (p<0.001). The positive correlation between CH and CRF was moderate (r=0.444) in Group 1 and strong (r=0.770) in Group 2. There was a moderate negative correlation between CH and IOPcc in both groups (r=-0.426, r=-0.423), but CH was not correlated with IOPg or IOPapl in either group. There were weak to strong positive correlations between CRF and all IOP values in both groups. There was no correlation between CRF and CCT in Group 1 (r=0.075) and a very weak correlation in Group 2 (r=0.237). Only IOPcc and IOPg were strongly correlated in both groups.

Conclusion

Better understanding of corneal biomechanical properties is essential for elucidating the pathophysiology and diagnosis of several corneal pathologies such as keratoconus. The biomechanical properties of keratoconic eyes seem to be closer to normal values after PKP.

Association between Rates of Visual Field Progression and Intraocular Pressure Measurements Obtained by Different Tonometers

PURPOSE

To investigate the associations between intraocular pressure (IOP) measurements obtained by different tonometric methods and rates of visual field loss in a cohort of patients with glaucoma followed over time.

DESIGN

Prospective, observational cohort study.

PARTICIPANTS

This study included 213 eyes of 125 glaucomatous patients who were followed for an average of 2.4±0.6 years.

METHODS

At each visit, IOP measurements were obtained using Goldmann applanation tonometry (GAT), the Ocular Response Analyzer (ORA) (Reichert, Inc., Depew, NY), corneal-compensated IOP (IOP_cc), and the ICare Rebound Tonometer (RBT) (Tiolat, Oy, Helsinki, Finland). Rates of visual field loss were assessed by standard automated perimetry (SAP) mean deviation (MD). Linear mixed models were used to investigate the relationship between mean IOP by each tonometer and rates of visual field loss over time, while adjusting for age, race, central corneal thickness, and corneal hysteresis.

MAIN OUTCOME MEASURES

Strength of associations (R²) between IOP measurements from each tonometer and rates of SAP MD change over time.

RESULTS

Average values for mean IOP over time measured by GAT, ORA, and RBT were 14.4±3.3, 15.2±4.2, and 13.4±4.2 mmHg, respectively. Mean IOP_cc had the strongest relationship with SAP MD loss over time (R² = 24.5%) and was significantly different from the models using mean GAT IOP (R² = 11.1%; 95% confidence interval [CI] of the difference, 6.6-19.6) and mean RBT IOP (R²= 5.8%; 95% CI of the difference, 11.1-25.0).

CONCLUSIONS

Mean ORA IOP_cc was more predictive of rates of visual field loss than mean IOP obtained by GAT or RBT. By correcting for corneal-induced artifacts, IOPcc measurements may present significant advantages for predicting clinically relevant outcomes in patients with glaucoma.

Review of Corneal Biomechanical Properties Following LASIK and SMILE for Myopia and Myopic Astigmatism

Worldwide, femtosecond Laser Assisted In-situ Keratomileusis (LASIK) is a well known and commonly used refractive technique, although Small Incision Lenticule Extraction (SMILE) has become increasingly popular since it was introduced in 2011. In LASIK, a corneal flap is cut with a microkeratome or femtosecond laser, followed by thinning of the stromal bed with excimer laser ablation. In SMILE, a minor intrastromal lenticule is cut with a femtosecond laser and subsequently removed through a small incision, leaving the anterior and strongest part of the cornea almost intact. Both LASIK and SMILE require cutting of corneal lamellae that may reduce the biomechanical stability of the cornea, with the potential risk of corneal iatrogenic ectasia as a severe complication. However, SMILE preserves the anterior corneal integrity and may, in theory, better preserve the corneal biomechanical strength than LASIK after surgery. A review aimed to examine the current literature that describes and compares the corneal biomechanical properties after Laser Assisted In-situ Keratomileusis (LASIK) and Small Incision Lenticule Extraction (SMILE). A comprehensive search was performed in Pubmed.gov using the following search queries: Corneal biomechanical properties, corneal biomechanics, ocular response analyser, ocular response analyzer, ORA, ex vivo, in vitro, Corvis, Corvis ST, LASIK, and SMILE.

Corneal Stiffness and Its Relationship With Other Corneal Biomechanical and Nonbiomechanical Parameters in Myopic Eyes of Chinese Patients

PURPOSE

To investigate corneal stiffness and analyze its possible influence on other corneal biomechanical and nonbiomechanical parameters in myopic eyes of Chinese patients.

METHODS

A total of 387 healthy (right) myopic eyes were included in this study. Corneal visualization Scheimpflug technology was used to measure the corneal stiffness parameter at the first applanation (SP-A1), deformation amplitude ratio at 1.0 and at 2.0 mm, biomechanically corrected intraocular pressure (bIOP), and corneal deformation parameters during the first applanation (A1: A1-time, A1-length, and A1-velocity), second applanation (A2; A2-time, A2-length, and A2-velocity), and highest concavity (HC; HC-time, HC-radius, HC deformation amplitude, and HC peak distance). The Pentacam was used to evaluate central corneal thickness, mean corneal curvature (Km), anterior corneal central elevation, and corneal asphericity (Q value) of the anterior cornea.

RESULTS

The mean SP-A1 and bIOP values were 90.46 ± 15.39 mm Hg/mm and 13.5 ± 1.85 mm Hg, respectively. The SP-A1 increased with age (β = 0.41 [95% CI, 0.28%-0.54%]; P < 0.0001). Corneal stiffness was positively correlated with central corneal thickness, uncorrected intraocular pressure, and bIOP and was also significantly correlated with all corneal deformation parameters (P < 0.05), except A1-length, HC-time, and HC peak distance. Corneal stiffness was negatively correlated with central corneal elevation of the anterior surface (r = -0.124, P = 0.014) and mean corneal curvature (r = -0.114, P = 0.025) and positively correlated with the Q value of the anterior surface (r = 0.109, P = 0.032).

CONCLUSIONS

Corneal stiffness increased with age. Stronger corneal stiffness may be related to flattening of the cornea and higher intraocular pressure.

In vivo assessment of changes in corneal hysteresis and lamina cribrosa position during acute intraocular pressure elevation in eyes with markedly asymmetrical glaucoma

Purpose

To investigate the biomechanical response of the cornea, lamina cribrosa (LC), and prelaminar tissue (PT) to an acute intraocular pressure (IOP) increase in patients with markedly asymmetrical glaucoma and in healthy controls.

Patients and methods

A total of 24 eyes of 12 patients with markedly asymmetrical primary open-angle glaucoma (POAG) and 12 eyes of 12 healthy patients were examined with spectral-domain optical coherence tomography (SD-OCT) and ocular response analyzer (ORA) at baseline and during acute IOP elevation by means of an ophthalmodynamometer. The displacement of the LC and PT and the change in corneal hysteresis (CH) and corneal resistance factor (CRF) were evaluated.

Results

Following a mean IOP increase of 12.3±2.4 mmHg, eyes with severe glaucoma demonstrated an overall mean anterior displacement of the LC (−6.58±26.09 µm) as opposed to the posterior laminar displacement in eyes with mild glaucoma (29.08±19.28 µm) and in healthy eyes (30.3±10.9; p≤0.001 and p=0.001, respectively). The PT displaced posteriorly during IOP elevation in all eyes. The CH decreased in eyes with severe glaucoma during IOP elevation (from 9.30±3.65 to 6.92±3.04 mmHg; p=0.012), whereas the CRF increased markedly in eyes with mild glaucoma (from 8.61±2.30 to 12.38±3.64; p=0.002) and in eyes with severe glaucoma (from 9.02±1.48 to 15.20±2.06; p=0.002). The increase in CRF correlated with the anterior displacement of the LC in eyes with severe glaucoma.

Conclusion

Eyes with severe glaucoma exhibited a mean overall anterior displacement of the anterior laminar surface, while eyes with mild glaucoma and healthy eyes showed a posterior displacement of the LC during IOP elevation. The CH decreased significantly from baseline only in eyes with severe glaucoma, but the CRF increased significantly in all glaucomatous eyes. The CRF increase correlated with the anterior displacement of the LC in eyes with severe glaucoma.

THz imaging system for in vivo human cornea

Terahertz (THz) imaging of corneal tissue water content (CTWC) is a proposed method for early, accurate detection and study of corneal diseases. Despite promising results from ex vivo and in vivo cornea studies, interpretation of the reflectivity data is confounded by the contact between corneal tissue and rigid dielectric window used to flatten the imaging field. This work develops a novel imaging system and image reconstruction methods specifically for nearly spherical targets such as human cornea. A prototype system was constructed using a 650 GHz multiplier source and Schottky diode detector. Resolution and imaging field strength measurement from characterization targets correlate well with those predicted by the quasioptical theory and physical optics analysis. Imaging experiments with corneal phantoms and ex vivo corneas demonstrate the hydration sensitivity of the imaging system and reliable measurement of CTWC. We present successful acquisition of non-contact THz images of in vivo human cornea, and discuss strategies for optimizing the imaging system design for clinical use.

In Vivo Noninvasive Measurement of Young's Modulus of Elasticity in Human Eyes: A Feasibility Study

PURPOSE

Abnormal ocular biomechanical properties may be important for understanding the risk of glaucoma. However, there are no clinical methods for measuring standard material properties in patients. In this feasibility study we demonstrated proof-of-principle for a novel method, ultrasound surface wave elastography (USWE), to determine the in vivo Young's modulus of elasticity of corneas in normal human eyes.

METHODS

In total, 20 eyes of 10 healthy subjects (mean age, 51.4±7.2; ±SD; range, 43 to 64 y) were studied. A spherical-tipped probe (3-mm diameter) was placed on closed eyelids and generated a gentle harmonic vibration at 100 Hz for 0.1 second. Wave speed propagation in the cornea was measured by USWE, and Young's modulus was calculated from the wave speed. Associations between Young's modulus and intraocular pressure (IOP), age, central corneal thickness, and axial length were explored by the Pearson correlation. Statistical significance was determined by using generalized estimating equation models to account for possible correlation between fellow eyes.

RESULTS

Mean IOP was 12.8±2.7 mm Hg. Mean wave speed in the cornea was 1.82±0.10 m/s. Young's modulus of elasticity was 696±113 kPa and was correlated with IOP (r=0.57; P=0.004), but none of the other variables (P>0.1).

CONCLUSIONS

USWE is a novel noninvasive technique for measuring ocular biomechanical properties. Corneal Young's modulus in normal eyes is associated with IOP, consistent with measurements in cadaver eyes. Further work is needed to determine elasticity in other ocular tissues, particularly the sclera, and if elasticity is altered in glaucoma patients.

Estimation of the Corneal Young's Modulus In Vivo Based on a Fluid-Filled Spherical-Shell Model with Scheimpflug Imaging

Current intraocular pressure (IOP) measurement using air puff could be erroneous without applying proper corrections. Although noncontact tonometry is not considered to be accurate, it is still popularly used by eye clinics. It is thus necessary to extract the correct information from their results. This study proposes a practical approach to correctly measure IOP in vivo. By embedding a new model-based correction to the Corvis® ST, we can extract the corneal Young's modulus from the patient data. This Young's modulus can be used to correct the IOP readings. The tests were applied to 536 right eyes of 536 healthy subjects (228 male and 308 female) between March of 2012 and April of 2016. The tests were applied to patients at the Department of Ophthalmology, National Taiwan University Hospital and the Hung-Chuo Eye Clinics. The statistical analysis showed that the value for the Young's modulus was independent of all the other parameters collected from the Corvis ST, including the corneal thickness and the intraocular pressure. Therefore, it is important to independently measure the Young's modulus instead of depending on the correlation with the other parameters. This study adds the methodology of measuring corneal stiffness in vivo for ophthalmologists' reference in diagnosis.

Characteristics of corneal biomechanical responses detected by a non-contact scheimpflug-based tonometer in eyes with glaucoma

PURPOSE

To determine the corneal biomechanical properties in eyes with glaucoma using a non-contact Scheimpflug-based tonometer.

METHODS

Corneal biomechanical responses were examined using a non-contact Scheimpflug-based tonometer. The tonometer parameters of the normal control group (n = 75) were compared with those of the glaucoma group (n = 136), including an analysis of glaucoma subgroups categorized by visual field loss.

RESULTS

After adjusting for potential confounding factors, including the intraocular pressure (IOP), central corneal thickness (CCT), age and axial length, the deformation amplitude was smaller in the glaucoma group (1.09 ± 0.02 mm) than in the normal control group (1.12 ± 0.02 mm; p value = 0.031). The deformation amplitude and the deflection amplitude of the severe glaucoma group (1.12 ± 0.02 mm and 0.92 ± 0.01 mm) were significantly greater than that of the early glaucoma group (1.07 ± 0.01  mm and 0.88 ± 0.11 mm, p = 0.006 and p = 0.031), whereas that of the moderate glaucoma group (1.09 ± 0.02 mm and 0.90 ± 0.02 mm) was greater than that of the early glaucoma group, but this difference was not statistically significant. The deformation amplitude showed a negative correlation with the CCT in the normal control group (r = -0.235), with a weaker negative relationship observed in the early glaucoma group (r = -0.099). However, in the moderate and severe glaucoma groups, the deformation amplitude showed a positive relationship with the CCT, showing an inverse relationship. The duration and number of antiglaucomatous eyedrops used had negative correlations with the CCT in eyes with moderate and severe glaucoma.

CONCLUSION

Overall, the glaucoma group showed significantly less deformable corneas than did the normal controls, even after adjusting for the IOP, CCT, age and axial length. However, there were also differences according to the severity of glaucoma, where the corneal deformation amplitude was greater in the severe glaucoma group compared to the early glaucoma group. The combined effects of stiffening due to glaucoma and increased viscoelastic properties caused by the chronic use of antiglaucomatous eyedrops may have resulted in the present findings.

Stromal striae: a new insight into corneal physiology and mechanics

We uncover the significance of a previously unappreciated structural feature in corneal stroma, important to its biomechanics. Vogt striae are a known clinical indicator of keratoconus, and consist of dark, vertical lines crossing the corneal depth. However we detected stromal striae in most corneas, not only keratoconus. We observed striae with multiple imaging modalities in 82% of 118 human corneas, with pathology-specific differences. Striae generally depart from anchor points at Descemet’s membrane in the posterior stroma obliquely in a V-shape, whereas in keratoconus, striae depart vertically from posterior toward anterior stroma. Optical coherence tomography shear wave elastography showed discontinuity of rigidity, and second harmonic generation and scanning electron microscopies showed undulation of lamellae at striae locations. Striae visibility decreased beyond physiological pressure and increased beyond physiological hydration. Immunohistology revealed striae to predominantly contain collagen VI, lumican and keratocan. The role of these regions of collagen VI linking sets of lamellae may be to absorb increases in intraocular pressure and external shocks.

Measurement repeatability of the dynamic Scheimpflug analyzer

PURPOSE

To evaluate the repeatability of corneal deformation parameters measured using a dynamic Scheimpflug analyzer and the impact of baseline clinical factors on the repeatability of each parameter.

STUDY DESIGN

Retrospective, cross-sectional study.

METHODS

Forty-eight eyes (48 healthy subjects; mean age, 49.0 ± 19.5 years) underwent repeated examinations with the Scheimpflug analyzer to evaluate the test-retest variability. The intraclass correlation coefficient (ICC) and repeatability coefficient as indicators of variability were computed for 35 parameters measured with the Scheimpflug analyzer. The associations between the magnitude of the test-retest variability and baseline factors, such as age, axial length (AL), intraocular pressure (IOP), and central corneal thickness (CCT), were analyzed.

RESULTS

The test-retest repeatability was excellent for 22 (62.9%) of 35 parameters (ICC ≥ 0.75), good for seven (20%), (ICC ≥ 0.6), fair for four (11.4%), (ICC ≥ 0.4), and poor for two (5.7%) parameters (ICC < 0.4). Age was associated positively with the magnitude of variability in 13 (37.1%) parameters; measurement variability was affected significantly by AL (5 parameters, 14.3%) and CCT (7 parameters, 20%) but, except for one parameter not by IOP.

CONCLUSION

Most parameters of the dynamic Scheimpflug analyzer showed favorable measurement reliability in healthy subjects. However, six parameters showed poor-to-fair repeatability. Age, AL, and CCT significantly affected the repeatability of several parameters. These results should be considered when clinicians use this device in clinical practice.

A viscoelastic anisotropic hyperelastic constitutive model of the human cornea

A constitutive model based on the continuum mechanics theory has been developed which represents interlamellar cohesion, regional variation of collagen fibril density, 3D anisotropy and both age-related viscoelastic and hyperelastic stiffening behaviour of the human cornea. Experimental data gathered from a number of previous studies on 48 ex vivo human cornea (inflation and shear tests) enabled calibration of the constitutive model by numerical analysis. Wide-angle X-ray scattering and electron microscopy provided measured data which quantify microstructural arrangements associated with stiffness. The present study measures stiffness parallel to the lamellae of the cornea which approximately doubles with an increase in strain rate from 0.5 to 5%/min, while the underlying stromal matrix provides a stiffness 2–3 orders of magnitude lower than the lamellae. The model has been simultaneously calibrated to within 3% error across three age groups ranging from 50 to 95 years and three strain rates across the two loading scenarios. Age and strain-rate-dependent material coefficients allow numerical simulation under varying loading scenarios for an individual patient with material stiffness approximated by their age. This present study addresses a significant gap in numerical representation of the cornea and has great potential in daily clinical practice for the planning and optimisation of corrective procedures and in preclinical optimisation of diagnostic procedures.

Corneal biomechanical properties in healthy children measured by corneal visualization scheimpflug technology

Background

The aim of this study was to evaluate corneal biomechanical properties in a population of healthy children in China using corneal visualization Scheimpflug technology (CST).

Methods

All children underwent complete bi-ocular examinations. CST provided intraocular pressure (IOP) and corneal biomechanical parameters, including time, velocity, length and deformation amplitude at first applanation (A1T, A1V, A1L, A1DA), at second applanation (A2T, A2V, A2L, A2DA), highest concavity time (HCT), maximum deformation amplitude (MDA), peak distance (PD), and radius of curvature (RoC). Pearson correlation analysis was used to assess the impacts of demographic factors, central corneal thickness (CCT), spherical equivalent (SE), and IOP on corneal biomechanics.

Results

One hundred eight subjects (32 girls and 76 boys) with the mean age of 10.80 ± 4.13 years (range 4 to18 years) were included in the final analyses. The right and left eyes were highly symmetrical in SE (p = 0.082), IOP (p = 0.235), or CCT (p = 0.210). Mean A1T of the right eyes was 7.424 ± 0.340 ms; the left eyes 7.451 ± 0.365 ms. MDA was 0.993 ± 0.102 mm in the right eyes and 0.982 ± 0.100 mm in the left eyes. Mean HCT of the right eyes was 16.675 ± 0.502 ms; the left eyes 16.735 ± 0.555 ms. All CST parameters of both eye were remarkably symmetrical with the exception of A2L (p = 0.006), A1DA (p = 0.025). The majority of CST parameters of both eyes were significantly correlated with CCT and IOP (p < 0.05). However, age, SE, and sex exert little influence on the CST measurements.

Conclusions

This study found interocular symmetry in corneal biomechanics in healthy children eyes. Several CST biomechanical parameters in children are modified by CCT and IOP.