Biomechanical parameters of the cornea measured with the Ocular Response Analyzer in normal eyes
© Kotecha et al.; licensee BioMed Central Ltd. 2014
Received: 5 July 2013
Accepted: 24 January 2014
Published: 30 January 2014
To evaluate the relationships between Reichert Ocular Response Analyzer (ORA) parameters corneal hysteresis (CH) and corneal response factor (CRF) and ocular dimensions, age and intraocular pressure.
Two hundred and twelve eyes of 212 participants with no ocular pathology had CH and CRF measured with the ORA. Intraocular pressure (IOP) was measured with the Dynamic Contour tonometer and central corneal thickness (CCT) was also evaluated. Partial least squares linear regression (PLSLR) analyses were performed to examine the relationships between each response variable, CH and CRF, and the predictor variables age, corneal curvature (CC), axial length (AL), CCT and IOP.
CH was positively associated with CCT and negatively associated with age (scaled coefficients: CCT 0.62, p < 0.0001; age -0.55, p <0.0001; r2 = 0.25). CRF was positively associated with CCT and DCT IOP and negatively associated with age and AL (scaled coefficients: CCT 0.89, p < 0.0001; DCT IOP 0.46, p < 0.01; age – 0.60, p < 0.0001; AL -0.37, p < 0.01; r2 = 0.43). There was no significant association between CC and CH or CRF.
The study suggests that age and CCT are strongly associated with CH and CRF, and that the latter is also influenced by AL and IOP. However, the variables studied could explain only 25% and 43% of the measured variation in CH and CRF, respectively, suggesting other factors also affect the values of these measurements.
Knowledge of the cornea’s biomechanical properties is important for a wide variety of applications within ophthalmology. Corneal biomechanical variations are known to affect the accuracy of intraocular pressure measurements [1, 2], may be used to identify early corneal disease [3, 4], and may assist with predicting refractive outcomes following corneal refractive surgery . It has also been suggested that corneal biomechanical properties may reflect globe biomechanics and thus give an indication of the susceptibility of developing glaucomatous damage [6, 7].
Until recently, most investigations evaluating corneal biomechanics were based on ex vivo tissue [8, 9]. However, the Reichert Ocular Response Analyzer (ORA; Reichert Ophthalmic Instruments, Buffalo, NY, USA) has facilitated an in vivo measurement of aspects of corneal biomechanical properties. The device measures the central corneal response to indentation by a rapid jet of air and provides two metrics of corneal biomechanics, corneal hysteresis (CH) and the corneal response factor (CRF). It is thought that CH predominantly reflects the viscous dampening properties of the cornea whilst CRF, a metric empirically derived to be strongly correlated with central corneal thickness (CCT), is thought to be most associated with the cornea’s elastic response . However, how these metrics relate to conventional biomechanical measurements is still relatively unclear. Studies examining the ORA biomechanical metrics have suggested associations with age, CCT and IOP [11, 12]; however, some suggest no association with these parameters.
The purpose of this study was to evaluate the relationships between ORA generated corneal biomechanical characteristics and ocular dimensions, age and IOP in a cohort of normal eyes. Understanding the relationships between ORA measures and ocular dimensions may help better elucidate the significance of, and applications for, the metrics it produces.
This was a prospective cross-sectional observational study. The study had the approval of the local research ethical committee (REC; London-City Road and Hampstead REC and Riverside REC, London, United Kingdom) and informed consent, according to the tenets of the Declaration of Helsinki, was obtained from each participant prior to examination. Study participants were recruited from staff, spouses and friends of patients attending Moorfields Eye Hospital, London, UK, or St. Luc University Hospital, Université Catholique de Louvain, Belgium. Data was collected between October 2008 and April 2011. Participants underwent a complete ophthalmic investigation, including visual field testing with the Humphrey Visual Field Analyser (HFA: Zeiss Humphrey Systems, Dublin, California, USA) SITA standard threshold 24–2 strategy. Participants were excluded from the study if they had any signs of corneal pathology, corneal astigmatism ≥2 dioptres or a history of incisional or intraocular surgery, a visual field defect, suspicious optic disc appearance, IOP ≥ 30 mmHg, a history of diabetes or a family history of glaucoma in a first degree relative. Soft-contact lens wearers were required to remove their lenses at least 24 hours prior to study participation; rigid contact lens wearers were excluded from the study. Of eligible participants, only one randomly chosen eye was measured for the study.
Study participants first had axial length (AL) and corneal curvature (CC) measurements made with the IOLMaster (version 3.01, Carl Zeiss Meditech AG, Jena, Germany). Prior to instillation of topical anaesthesia, participants underwent ORA measurements and three good quality waveform scans, defined as having symmetry in height between the two peaks of the waveform, were recorded and the mean value used in subsequent analysis.
Following instillation of topical corneal anaesthesia (proxymetacaine hydrochloride 0.5% with fluorescein sodium 0.25%), IOP measurements were made using both the Goldmann applanation tonometer (GAT) and the dynamic contour tonometer (DCT). Two GAT IOP and three DCT IOP measurements were made in a randomised order. Only DCT measurements with a 'quality’ reading of 1, 2 or 3 were accepted and the first DCT reading was discarded in accordance with manufacturer’s guidance. A minimum two-minute interval was left between IOP measurements with each device to minimise the tonographic effects of repeated tonometry measurements . The mean of two IOP readings, taken with each of the two instruments, was calculated for each participant and was used in the analyses. Measurements of CCT, using an ultrasound pachymeter (Altair, Optikron 2000, Roma, Italy), were made at the end of the visit; the average of 3 readings was recorded. All measurements were performed by one of four investigators (AK, AS, DS or SP).
Based on previous pilot data from an unrelated dataset, it was calculated that a sample of 194 eyes was required to achieve a correlation between CH and age of rho = -0.2, with 80% power at the p < 0.05 level.
Multivariate statistical analyses are effective tools to identify and explore the relationships between a response variable and several predictor variables. In this study, the response variable of interest was either CH or CRF and the predictor variables of interest were age, CC, CCT, IOP and AL. The theoretical assumptions of standard multiple linear regression (MLR) analysis limit its application to explore the relationships between CH, CRF and the predictors for this study as the predictor variables are correlated, a statistical phenomenon known as multicollinearity. In circumstances where there is multicollinearity in the predictors, partial least squares linear regression (PLSLR) offers a robust alternative to MLR . Like MLR, the objective of PLSLR is to describe the relationships between the response and predictor variables.
For this study, two PLSLR models were fitted to a subset of the data under investigation. The 'calibration’ dataset consisted of a proportion (80%) of participants randomly selected from the complete sample used to construct a predictive model for CH and CRF. The number of components in each PLSLR model was chosen using the 'leave-one-out method’ . The predictive performance of each model was then tested on the data excluded from the calibration data, which consisted of the remainder of participants. As the predictor variables have differing units of measurement, PLSLR models were also fitted using scaled predictor variables (where each predictor was weighted by dividing it by its standard deviation) in order to evaluate the relative impact of each variable on the response.
Demographics of study cohort
Total cohort n = 212
19.0 to 92.6
21.4 to 28.7
7.1 to 8.6
Corneal astigmatism (dioptres)
0.0 to 4.0
490 to 633
6.0 to 25.5
9.7 to 25.0
ORA IOPcc (mmHg)
8.7 to 29.0
Correlations between CH, CRF and other parameters
Correlation table showing Spearman’s rho and significance values for parameters
Spearman’s rho coefficient (p)
GAT IOP (mmHg)
DCT IOP (mmHg)
GAT IOP (mmHg)
DCT IOP (mmHg)
Previous work has shown the DCT to be less affected by variations in corneal biomechanical properties ; consequently, for this study, the DCT IOP measurement was used as a surrogate measure of 'true’ IOP. Thus, in total, the impact of five predictor variables on CH and CRF were considered.
PLSLR analysis: corneal hysteresis and corneal response factor models
For the CH model, four components explained 25% of the variance in CH. For the CRF model, four components explained 43% of the variance in CRF.
Coefficients of PLSLR predictive model for CH in the calibration dataset
PLSLR coefficient (non-scaled)
PLSLR coefficient (scaled)
DCT IOP (mmHg)
Coefficients of PLSLR predictive model for CRF in the calibration dataset
PLSLR coefficient (non-scaled)
Scaled PLSLR coefficient
DCT IOP (mmHg)
The results of this study suggest that ORA-generated metrics of corneal biomechanical properties in non-glaucomatous eyes are correlated with age, CCT and, in the case of CRF, AL. The relationships between these variables are complex and the predictor variables do not explain all the variation in CH or CRF measurements.
The data were analysed with PLSLR analysis since the predictor variables were correlated. In PLSLR, components (known as 'latent variables’) are extracted from the predictor variables by maximising the covariance between the predictor variables and the response variable.
It has been suggested that CH represents the viscoelastic capacity of the cornea, that is, the cornea’s ability to dampen and dissipate applied energy. Hysteresis is dependent on the relative contributions of both elasticity and viscosity, and it has been shown that alterations in either component will have very different and sometimes opposing effects on measured hysteresis . Our finding that CH reduces with age corroborates the findings of experimental ex vivo studies that show an increase in collagen cross-linking with age , which result in a reduction in the viscosity of the cornea and thus an increase in the 'stiffness’ of the structure [19, 20]. The data are also in agreement with previous clinical work evaluating the effect of age on ORA measured corneal biomechanical properties [11, 12, 21]. Taken together, these findings suggest that aging results in an overall reduction in the dampening capacity of the cornea.
Our study also agrees with previous findings that CH is positively associated with CCT [11, 21]. This is unsurprising as one might expect a thicker cornea to have a greater viscous dampening capacity. Recent work has also found that CH reduces with increasing axial length , and it has been suggested that these findings are indicative of altered biomechanics in axially myopic eyes. However, other reports have found no association between axial length and CH . In the present study, Spearman’s correlation test suggested that CH was negatively associated with both axial length and corneal curvature; however, in the PLSLR model, these parameters were not significant (although axial length approached significance). These findings may be a reflection of the ethnic differences in the cohorts studied; our participants were predominantly of Caucasian ethnicity, whilst those in the other studies were of Far East Asian origin. Surprisingly, Spearman’s correlation test also suggested that CH was positively associated with DCT IOP and GAT IOP; however, in the PLSLR model the association between DCT IOP and CH only approached significance. This suggests that the univariate association was either a spurious finding, or as a result of correlations with other predictor variables that were subsequently dealt with in the PLSLR model. Further work is required to establish the significance of relationships that may exist between axial length, IOP and CH.
CRF was also found to increase with CCT and reduce with age. CRF was intended to quantify the overall corneal viscoelastic resistance to indentation and was developed to be strongly associated with CCT . The fact that CRF reduces with age is counterintuitive, as it might be expected that the increase in corneal 'stiffness’ resultant from an age-related increase in corneal collagen cross-linking would result in an increased resistance to deformation; however, our results do agree with previous findings . The CRF represents a metric of corneal resistance to a near instantaneous indentation force applied axially. Our data suggest that the assumption that CRF reflects overall corneal rigidity may be an oversimplification and that other factors need to be considered when interpreting its value.
In the PLSLR model, CRF was also found to be positively associated with IOP. The cornea is hyperelastic; thus, the 'stretching’ of the cornea under high IOP conditions will render it a 'stiffer’ structure. Of interest is the finding that a lower CRF was associated with longer axial length, which was highlighted in the PLSLR model. This agrees with some  but not others’  findings. It is not clear what this relationship between CRF and axial length biomechanically represents. One may postulate that if CRF represents overall corneal rigidity, since longer eyes may have altered scleral biomechanics  a reduction in CRF may be an expected finding, in that a reduced corneal rigidity would be associated with a reduced scleral rigidity. However, further work is required to establish the significance of the relationship and how corneal biomechanics relate, if at all, to scleral biomechanics.
Of note is the proportion of variance in corneal biomechanical measurements predicted by the studied parameters when all the inter-relationships are considered. In the calibration PLSLR model, only 25% of the variation in measured CH could be explained by the variables studied, which reduced even further when applied to the test dataset. Whilst part of this may indicate true inter-individual variations in CH measurements, there may be other explanations for this observation. It is possible that this reflects a reduced signal-to-noise ratio in the CH measurement which masks the true effect of the predictors studied . Furthermore, there may be significant and as yet unmeasured variables that effect the CH measurement. These may include the degree of corneal indentation and the rate and maximal level of external air pressure application, all of which have been shown experimentally to affect hysteresis [17, 29]. In contrast, approximately 43% of CRF variation could be explained by the variables studied. This suggests that the CRF measurement can be characterised by CCT, age, IOP and AL. However, a large proportion of the measurement variation cannot be explained by these variables and may represent a combination of actual inter-individual CRF differences, measurement noise and other as yet undetermined factors.
In conclusion, this study of corneal biomechanics in normal eyes describes the complex interactions between ocular characteristics and ORA metrics, and finds that both age and CCT are significantly associated with CH and CRF. However, only a proportion of the variation in both these ORA metrics could be described by the ocular characteristics measured, implying that there are other elements contributing to the CH and CRF measurement.
Finally, condensing corneal biomechanical measures to a single summary metric will never completely describe the cornea’s properties. Recently, investigators evaluating the ORA applanation signal have found that variations in specific signal elements are better descriptors of the corneal response to indentation particularly following refractive surgery procedures [30, 31]. Further work is required to establish the validity of these new parameters and how they relate to more conventional biomechanical measures.
The authors would like to thank Reichert Corporation for the loan of the ORA.
Drs Kotecha and Russell and Professor Garway-Heath have received a proportion of their funding from the Department of Health’s National Institute for Health Research Biomedical Research Centre for Ophthalmology at Moorfields Eye Hospital NHS Foundation Trust and the UCL Institute of Ophthalmology.
Presented in part at the Association for Research in Vision and Ophthalmology meeting, Fort Lauderdale, USA, May 2010.
- Lau W, Pye D: Changes in corneal biomechanics and applanation tonometry with induced corneal swelling. Invest Ophthalmol Vis Sci. 2011, 52 (6): 3207-3214. 10.1167/iovs.10-6754.View ArticlePubMedGoogle Scholar
- Tonnu PA, Ho T, Newson T, El Sheikh A, Sharma K, White E, Bunce C, Garway-Heath D: The influence of central corneal thickness and age on intraocular pressure measured by pneumotonometry, non-contact tonometry, the Tono-Pen XL, and Goldmann applanation tonometry. Br J Ophthalmol. 2005, 89 (7): 851-854. 10.1136/bjo.2004.056622.View ArticlePubMedPubMed CentralGoogle Scholar
- Schweitzer C, Roberts CJ, Mahmoud AM, Colin J, Maurice-Tison S, Kerautret J: Screening of forme fruste keratoconus with the ocular response analyzer. Invest Ophthalmol Vis Sci. 2010, 51 (5): 2403-2410. 10.1167/iovs.09-3689.View ArticlePubMedGoogle Scholar
- del Buey MA, Cristobal JA, Ascaso FJ, Lavilla L, Lanchares E: Biomechanical properties of the cornea in Fuchs’ corneal dystrophy. Invest Ophthalmol Vis Sci. 2009, 50 (7): 3199-3202. 10.1167/iovs.08-3312.View ArticlePubMedGoogle Scholar
- Roberts C: Biomechanics of the cornea and wavefront-guided laser refractive surgery. J Refract Surg. 2002, 18 (5): S589-S592.PubMedGoogle Scholar
- Congdon NG, Broman AT, Bandeen-Roche K, Grover D, Quigley HA: Central corneal thickness and corneal hysteresis associated with glaucoma damage. Am J Ophthalmol. 2006, 141 (5): 868-875. 10.1016/j.ajo.2005.12.007.View ArticlePubMedGoogle Scholar
- Wells AP, Garway-Heath DT, Poostchi A, Wong T, Chan KC, Sachdev N: Corneal biomechanics and optic nerve head compliance: corneal hysteresis but Not corneal thickness correlates with optic nerve surface deformation in glaucoma patients. Invest Ophthalmol Vis Sci. 2008, 49 (8): 3262-3268. 10.1167/iovs.07-1556.View ArticlePubMedGoogle Scholar
- Dupps WJ, Wilson SE: Biomechanics and wound healing in the cornea. Exp Eye Res. 2006, 83 (4): 709-720. 10.1016/j.exer.2006.03.015.View ArticlePubMedPubMed CentralGoogle Scholar
- Dupps WJ: Hysteresis: new mechanospeak for the ophthalmologist. J Cataract Refract Surg. 2007, 33 (9): 1499-1501. 10.1016/j.jcrs.2007.07.008.View ArticlePubMedGoogle Scholar
- Kotecha A: What biomechanical properties of the cornea are relevant for the clinician?. Surv Ophthalmol. 2007, 52 (Suppl 2): S109-S114.View ArticlePubMedGoogle Scholar
- 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 (12): 5337-5347. 10.1167/iovs.06-0557.View ArticlePubMedGoogle Scholar
- Leite MT, Alencar LM, Gore C, Weinreb RN, Sample PA, Zangwill LM, Medeiros FA: Comparison of corneal biomechanical properties between healthy blacks and whites using the Ocular Response Analyzer. Am J Ophthalmol. 2010, 150 (2): 163-168. 10.1016/j.ajo.2010.02.024. e161View ArticlePubMedPubMed CentralGoogle Scholar
- Recep OF, Hasiripi H, Vayisoglu E, Kalayci D, Sarikatipoglu H: Accurate time interval in repeated tonometry. Acta Ophthalmol Scand. 1998, 76 (5): 603-605. 10.1034/j.1600-0420.1998.760518.x.View ArticlePubMedGoogle Scholar
- Wold H: Estimation of Principal Components and Related Models by Iterative Least Squares. Multivariate Analysis. Edited by: Krishnaiah PR. 1966, New York: AcademicGoogle Scholar
- Lachenbruch PA, Mickey MR: Estimation of error rates in discriminant analysis. Technometrics. 1968, 10 (1): 1-11. 10.1080/00401706.1968.10490530.View ArticleGoogle Scholar
- Boehm AG, Weber A, Pillunat LE, Koch R, Spoerl E: Dynamic contour tonometry in comparison to intracameral IOP measurements. Invest Ophthalmol Vis Sci. 2008, 49 (6): 2472-2477. 10.1167/iovs.07-1366.View ArticlePubMedGoogle Scholar
- Glass DH, Roberts CJ, Litsky AS, Weber PA: A viscoelastic biomechanical model of the cornea describing the effect of viscosity and elasticity on hysteresis. Invest Ophthalmol Vis Sci. 2008, 49 (9): 3919-3926. 10.1167/iovs.07-1321.View ArticlePubMedGoogle Scholar
- Daxer A, Misof K, Grabner B, Ettl A, Fratzl P: Collagen fibrils in the human corneal stroma: structure and aging. Invest Ophthalmol Vis Sci. 1998, 39 (3): 644-648.PubMedGoogle Scholar
- Elsheikh A, Wang D, Brown M, Rama P, Campanelli M, Pye D: Assessment of corneal biomechanical properties and their variation with age. Curr Eye Res. 2007, 32 (1): 11-19. 10.1080/02713680601077145.View ArticlePubMedGoogle Scholar
- Knox Cartwright NE, Tyrer JR, Marshall J: Age-related differences in the elasticity of the human cornea. Invest Ophthalmol Vis Sci. 2011, 52 (7): 4324-4329. 10.1167/iovs.09-4798.View ArticlePubMedGoogle Scholar
- Lau W, Pye D: A clinical description of Ocular Response Analyser measurements. Invest Ophthalmol Vis Sci. 2011, 52 (6): 2911-2916. 10.1167/iovs.10-6763.View ArticlePubMedGoogle Scholar
- Narayanaswamy A, Chung RS, Wu RY, Park J, Wong WL, Saw SM, Wong TY, Aung T: Determinants of corneal biomechanical properties in an adult chinese population. Ophthalmology. 2011, 118 (7): 1253-1259.PubMedGoogle Scholar
- Lim L, Gazzard G, Chan YH, Fong A, Kotecha A, Sim EL, Tan D, Tong L, Saw SM: Cornea biomechanical characteristics and their correlates with refractive error in Singaporean children. Invest Ophthalmol Vis Sci. 2008, 49 (9): 3852-3857. 10.1167/iovs.07-1670.View ArticlePubMedGoogle Scholar
- Methodology For Corneal Compensated Iop And Corneal Resistance Factor For An Ocular Response Analyzer. Invest Ophthalmol Vis Sci. 2006, 47: E-Abstract 2266-
- Fontes BM, Ambrosio R, Alonso RS, Jardim D, Velarde GC, Nose W: Corneal biomechanical metrics in eyes with refraction of -19.00 to +9.00 D in healthy Brazilian patients. J Refract Surg. 2008, 24 (9): 941-945.PubMedGoogle Scholar
- Chang PY, Chang SW, Wang JY: Assessment of corneal biomechanical properties and intraocular pressure with the Ocular Response Analyzer in childhood myopia. Br J Ophthalmol. 2010, 94 (7): 877-881. 10.1136/bjo.2009.158568.View ArticlePubMedGoogle Scholar
- McBrien NA, Jobling AI, Gentle A: Biomechanics of the sclera in myopia: extracellular and cellular factors. Optom Vis Sci. 2009, 86 (1): E23-E30. 10.1097/OPX.0b013e3181940669.View ArticlePubMedGoogle Scholar
- Kotecha A, Crabb DP, Spratt A, Garway-Heath DF: The relationship between diurnal variations in intraocular pressure measurements and central corneal thickness and corneal hysteresis. Invest Ophthalmol Vis Sci. 2009, 50 (9): 4229-4236. 10.1167/iovs.08-2955.View ArticlePubMedGoogle Scholar
- Elsheikh A, Wang D, Rama P, Campanelli M, Garway-Heath D: Experimental assessment of human corneal hysteresis. Curr Eye Res. 2008, 33 (3): 205-213. 10.1080/02713680701882519.View ArticlePubMedGoogle Scholar
- Mikielewicz M, Kotliar K, Barraquer RI, Michael R: Air-pulse corneal applanation signal curve parameters for the characterisation of keratoconus. Br J Ophthalmol. 2011, 95 (6): 793-798. 10.1136/bjo.2010.188300.View ArticlePubMedGoogle Scholar
- Kerautret J, Colin J, Touboul D, Roberts C: Biomechanical characteristics of the ectatic cornea. J Cataract Refract Surg. 2008, 34 (3): 510-513. 10.1016/j.jcrs.2007.11.018.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2415/14/11/prepub
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.