Corrigendum: Distribution of corneal geometric landmarks and relationship between their distances and biomechanical parameters in the development of keratoconus
- 1Beijing Institute of Ophthalmology, Beijing Tongren Eye Center, Beijing Tongren Hospital, Capital Medical University and Beijing Ophthalmology and Visual Sciences Key Laboratory, Beijing, China
- 2Beijing Advanced Innovation Center for Big Data-Based Precision Medicine, Beihang University & Capital Medical University, Beijing, China
- 3School of Biomedical Engineering, Capital Medical University, Beijing, China
- 4Department of Medical Engineering, Peking Union Medical College Hospital, Chinese Academy of Medical Sciences, Beijing, China
- 5The First People's Hospital of Xuzhou, Jiangsu, China
- 6Beijing Key Laboratory of Fundamental Research on Biomechanics in Clinical Application, Capital Medical University, Beijing, China
- 7Department of Ophthalmology, Chinese People's Liberation Army General Hospital, Beijing, China
Purpose: To analyze the changes in coordinates and distances among three typical geometric landmarks of the cornea, namely, the thinnest point (TP), maximum curvature (Kmax), and corneal apex (AP) during the development of keratoconus, and explore the potential relationship between these changes and the abnormalities of corneal biomechanics.
Methods: Normal eyes (n = 127), clinical keratoconic eyes (CKC, n = 290), and the eyes of forme fruste keratoconus (FFKC, n = 85) were included; among them, the CKC group was classified into four grades based on the Topographic Keratoconus Classification (TKC) provided by Pentacam. A total of 38 Corvis ST output parameters and three distance parameters of three typical landmarks (DKmax-AP, DTP-AP, and DKmax-TP) based on Pentacam were included. The differences of parameters among the abovementioned six groups (Normal, FFKC, and CKC stage I to CKC stage IV) were analyzed. Spearman’s rank correlation test was performed to choose several dynamic corneal response (DCR) parameters that could best reflect the changes of corneal biomechanical characteristics during the progression of the disease, and the Pearson’s or Spearman’s correlation test was conducted to determine the association between the three distances and the selected DCR parameters in each grade. In addition, by flipping the X coordinate of the left eye on the vertical axis to reflect the direction of the right eye, the coordinates of TP and Kmax in different developmental grades were highlighted.
Results: From CKC stage II, the three geometric landmark distances commenced to correlate with the corneal DCR parameters (CBI, SPA1, IR, DA Ratio 2, ARTh, MIR, Radius, Pachy, and DA Ratio 1), which could better represent the changes of biomechanical properties from normal cornea to keratoconus. From normal cornea to CKC stage IV, the coordinates of Kmax were gradually tended to the inferior temporal region from dispersion, while TP was always concentrated in the inferior temporal region. Although DKmax-AP, DKmax-TP, and DTP-AP all showed a gradual decreasing trend with the progress of the disease, the first two did not change significantly, and only DTP-AP significantly approached AP in the later stage of disease development. In addition, from the FFKC group, the corresponding values of DKmax-TP in each disease development group were smaller than DKmax-AP.
Conclusions: In the later stage of keratoconus, the relationship between the three typical landmark distance parameters and DCR parameters is stronger, and even the weakening of corneal biomechanical properties may be accompanied by the merger of typical landmark positions.
Introduction
With respect to in-depth understanding of the biomechanical mechanism of corneal diseases, evaluation of the biomechanical properties of cornea has greatly attracted scholars’ attention in terms of prevention and treatment of corneal dilatation diseases, especially keratoconus (Herber et al., 2019). It has been gradually found that the biomechanical properties of the cornea depend on the collagen fiber, collagen fiber bundle, and their spatial structure composition (Oxlund and Simonsen, 1985), and studies (Meek et al., 2005; Catalán-López et al., 2018) have shown that the mechanical strength of keratoconus is often remarkably weaker than that of normal cornea. In addition, a number of scholars have pointed out that the morphological changes of keratoconus may be secondary to the changes in corneal biomechanics (Roberts et al., 2017; Sedaghat et al., 2018).
At present, the Ocular Response Analyzer (ORA) and Corneal Visualization Scheimpflug Technology (Corvis ST) are the two most recognized devices for the measurement of cornea biomechanics in vivo. Among them, ORA cannot display the process of corneal deformation dynamically in real time, and its main biomechanical parameters, corneal hysteresis (CH) and corneal resistance factor (CRF), are derived by analyzing the measured waveforms (McMonnies, 2012), while Corvis ST can dynamically record the whole process under impulse pressure and generate DCR parameters to reflect the biomechanical characteristics of cornea (Elham et al., 2017). For now, the combined diagnosis of keratoconus with Corvis ST and three-dimensional anterior segment analysis and diagnosis system Pentacam, which characterizes the morphological features of cornea, has been clinically recognized.
That is to say, although morphological changes are not the primary cause of keratoconus, its diagnostic value cannot be easily ignored. As we all know, the most typical morphological change of keratoconus is the thinning and protruding of cornea (Hashemi et al., 2019). To date, numerous studies (Galletti et al., 2015; Chan et al., 2018) have concentrated on the analysis of the numerical changes of the thinnest thickness and the maximum curvature caused by the gradual corneal protrusion, while few people have analyzed the coordinates of these two typical geometric landmarks and the distance between them and the central apex of the cornea. Only several studies simply compared the distance between the cornea apex and the thinnest point of cornea in different groups, and found that there were remarkable differences among normal cornea, subclinical keratoconus (Muftuoglu et al., 2013), and clinical keratoconus (Muftuoglu et al., 2015).
As mentioned above, DCR parameters and distance parameters of typical geometric landmarks representing morphological characteristics are both significantly different between normal cornea and keratoconus. Thus, we speculated that there may be a certain correlation between the DCR parameters and distance parameters with the development of keratoconus disease.
The main purpose of the present study was to evaluate the relationship between the distances among the typical geometric landmarks of cornea and the DCR parameters output by Corvis ST in the assumed different grades of keratoconus development, and explore more potential patterns of disease development.
Patients and Methods
Study Subjects
This prospective comparative study included patients with clinical keratoconus (CKC group), forme fruste keratoconus (FFKC group), and candidates undergoing refractive surgery with normal corneas (Normal group).
A diagnosis of keratoconus was made if the eye met the following conditions (Tian et al., 2014; Cui et al., 2016): (1) an irregular cornea determined by distorted keratometry mires, distortion of the retinoscopic or ophthalmoscopic red reflex (or a combination of the two); (2) with at least one of the following biomicroscopic signs: Vogt’s striae, Fleischer’s ring of >2 mm arc, or corneal scarring consistent with keratoconus. In the CKC group, the seven classes of topographic keratoconus classification (TKC) provided by Pentacam (Oculus, Wetzlar, Germany) were included (i.e., 1, 1-2, 2, 2-3, 3, 3-4, and 4). Those classes could be divided into four four subgroups based on their TKC number: TKC 1 was classified as CKC stage I, TKC 1-2 and 2 were classified as CKC stage II, TKC 2-3 and 3 were classified as CKC stage III, and TKC 3-4 and 4 were classified as CKC stage IV.
An eye was diagnosed with FFKC if it was the fellow eye of a patient with keratoconus and showed the following features (Peña-García et al., 2016): (1) a normal-appearing cornea on slit-lamp examination, retinoscopy, and ophthalmoscopy; (2) topography was normal with no asymmetric bowtie and no focal or inferior steepening pattern; (3) the level of TKC was normal, namely, it was “-”; and (4) the patient had no history of contact lens use, ocular surgery, or trauma. For participants undergoing refractive surgery, only one eye from each individual was chosen using a random numbers table. In addition, the TKC level of eyes in the Normal group and FFKC group was normal.
Exclusion criteria were a history of undergoing ocular surgery and cases with eye diseases that may potentially affect the outcomes. For contact lens-wearing patients, they were asked to remove soft contact lenses at least 2 weeks and rigid contact lenses at least 1 month before assessment. Data were collected from May 2013 to January 2020 at the Beijing Tongren Hospital affiliated to Capital Medical University (Beijing, China). All participants signed the written informed consent form prior to commencing the study. The study was carried out in accordance with the Declaration of Helsinki, and it was approved by the Ethics Committee of the Beijing Tongren Hospital affiliated to Capital Medical University.
Ocular Examination
A comprehensive ocular examination was performed on all eyes, including uncorrected visual acuity, slit-lamp and fundoscopic examinations, Pentacam tomographic examination, and Corvis ST (Oculus; Wetzlar, Germany) measurements. All measurements were undertaken between 9:00 a.m. and 5:00 p.m. by the same trained ophthalmologists during the same visit.
Pentacam Measurement
The Pentacam software (ver. 1.20r134) reconstructs a three-dimensional (3D) image of the entire anterior segment from the anterior surface of the cornea to the posterior surface of the lens by utilizing the high-speed rotating Scheimpflug system. Details and principles of the Pentacam are described elsewhere (Cui et al., 2016). Only scans that the Pentacam’s “quality specification” (QS) function determined as “OK” were included for analysis.
In the study, we focused on three points on the cornea: corneal apex (AP), thinnest point (TP), and maximum curvature (Kmax). Extracted parameters from Pentacam data for analysis included the coordinates of TP (TP X, TP Y) and Kmax (Kmax X, Kmax Y); then, we calculated the absolute distances from the cornea apex (geometric center of the examination [x = 0; y = 0]): DTP-AP and DKmax-AP, and the distance between Kmax and TP (DKmax-TP). The formula of distance was as follows:
Corvis ST Measurement
The Corvis ST software (ver. 1.5r1902) measures dynamic corneal deformation response to an air-puff pulse. Details and principles of the Corvis ST are described elsewhere (Elham et al., 2017). The measurements were checked under the QS window; only correct measurements were accepted (comment box reading “OK”). If the comment box was marked yellow or red, the examination was repeated. The following parameters were detected by Corvis ST: intraocular pressure (IOP), Pachymetry (Pachy), time from the initiation of air puff until the first applanation (A1T), second applanation (A2T) and maximum deformation (HCT), corneal velocity at the first (A1V) and second applanation (A2V), peak distance (PD) and radius of curvature (Radius), deformation amplitude at the first applanation (A1DA), second applanation (A2DA) and highest concavity (HCDA), deflection length at the first applanation (A1DLL), second applanation (A2DLL) and highest concavity (HCDLL), deflection amplitude at the first applanation (A1DLA), second applanation (A2DLA) and highest concavity (HCDLA), deflection area at the first applanation (A1DLAr), second applanation (A2DLAr) and highest concavity (HCDLAr), delta arc length at the first applanation (A1dArcL), second applanation (A2dArcL) and highest concavity (HCdArcL), max time and length at deflection amplitude (DLAMT, DLAML), max time and amplitude of whole eye movement (WEMT, WEMA), delta arc length max (dArcLM) and PachySlope, the maximal value of the ratio between deformation amplitude at the apex and that at 1 (DA Ratio 1) and 2 mm (DA Ratio 2) from the corneal apex, max inverse radius (MIR) and integrated radius (IR), Ambrósio relational thickness to the horizontal profile (ARTh), Biomechanical-corrected intraocular pressure (bIOP), stiffness parameter at first applanation (SPA1), and Corvis biomechanical index (CBI).
Statistical Analysis
Statistical analysis was performed using SPSS 20.0 (IBM, Armonk, NY, USA). Drawing was completed by R Core Team (version 3.6.1; R Foundation for Statistical Computing, Vienna, Austria; https://www.R-project.org/) software and GraphPad Prism software version 8.0, respectively.
The Shapiro–Wilk test was used to assess normal distribution of quantitative data. The normally distributed data were expressed as the mean ± standard deviation (SD), while abnormally distributed data were presented as median and range of variation.
One-way analysis of variance (ANOVA) test or non-parametric Kruskal–Wallis test was used to analyze the differences among the four subgroups of CKC group, Normal group, and FFKC group. Spearman’s rank correlation test was performed to assess correlation among all parameters measured by Corvis ST and the developmental grades of keratoconus (Rank-group). Then, the Pearson’s or Spearman’s correlation test was applied to determine the association between the distance of three geometric landmarks and the above selected DCR parameters in each grade, and Bonferroni correction was performed (p < 0.0056). Moreover, we plotted the variation trend of three typical landmark distance parameters and the selected DCR parameters with the progress of disease stage. The differences of parameters in any disease stage and its adjacent previous disease stage were compared by independent sample t-test or nonparametric Mann–Whitney test, and Bonferroni correction was performed (p < 0.01).
In addition, with flipping the X coordinate of the left eye on the vertical axis to reflect the direction of the right eye, the coordinates of Kmax and TP in different developmental grades of keratoconus were drawn. p < 0.05 was considered statistically significant, except for the Bonferroni correction.
Results
Herein, 290 eyes of 223 patients (mean age, 23.19 ± 7.39 years old; range of age, 9–53 years old) were assigned to the CKC group, of whom both eyes of 59 patients were included, one eye of 85 patients was involved because of unilateral keratoconus, and one eye of 79 patients was included because the fellow eye had undergone eye surgery or the quality of the examination did not meet the predefined criteria. The normal contralateral eye of the unilateral keratoconus constituted the FFKC group (mean age, 23.61 ± 7.73 years old; range of age, 10–49 years old). The Normal group consisted of 127 normal individuals (mean age, 24.39 ± 4.38 years old; range of age, 15–37 years old), and only one eye per person was randomly evaluated. There were no statistically significant differences between the groups in age distribution (p = 0.065, the Kruskal–Wallis test).
With the exception of DKmax-AP, HCT, A2DA, A2DLL, DLAMT, and WEMA, statistically significant differences in other parameters were found among four subgroups of the CKC group, Normal group, and FFKC group (p < 0.05) (Table 1).
Figure 1 shows the correlation among 38 parameters measured by Corvis ST and developmental grades of keratoconus. The results unveiled that there were nine DCR parameters that were strongly correlated with Rank-group (|r| > 0.6), and they were CBI, SPA1, IR, DA Ratio 2, ARTh, MIR, Radius, and Pachy in the order of correlation from high to low. Figure 2 shows the change trend of the nine DCR parameters screened above with the increase of disease stage. Obviously, these nine parameters would change significantly in the process of disease progression (normal to CKC stage IV), that is, gradually increased or decreased.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/766163/fbioe-09-766163-HTML-r3/image_m/fbioe-09-766163-g001.jpg)
FIGURE 1. Spearman correlations of DCR parameters and the developmental grades of keratoconus (strong correlation: 0.6 ≤ |r| < 1 and p < 0.05, weak correlation: |r| < 0.6 and p < 0.05, irrelevant: p ≥ 0.05).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/766163/fbioe-09-766163-HTML-r3/image_m/fbioe-09-766163-g002.jpg)
FIGURE 2. The change tendency of DCR parameters, which are highly correlated as disease progression with the disease stage increasing. a, b, c, d, and e indicate that the parameters in a certain stage are statistically different from the corresponding parameters in the previous adjacent stage (p < 0.01); 0 = Normal group, 1 = FFKC group, 2 = CKC Stage I group, 3 = CKC Stage II group, 4 = CKC Stage III group, 5 = CKC Stage IV group.
Figure 3 illustrates the coordinates of Kmax and TP in different developmental grades of keratoconus, and Figure 4 depicts the variation trend of three typical landmark distance parameters of geometric landmarks with the progress of disease stage. As shown in Figure 3, from normal to CKC stage IV, the coordinates of Kmax were gradually concentrated from dispersion, and finally tended to the inferior temporal region, while the coordinates of TP were always concentrated in the inferior temporal region. However, it can be seen from Figure 4 that although DKmax-AP, DKmax-TP, and DTP-AP all showed a gradual decreasing trend with the progress of the disease, the first two (DKmax-AP and DKmax-TP) did not change significantly, and only DTP-AP significantly approached AP in the later stage of disease development.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/766163/fbioe-09-766163-HTML-r3/image_m/fbioe-09-766163-g003.jpg)
FIGURE 3. The coordinate distribution of TP and Kmax in different developmental grades of keratoconus, where the X coordinate of the left eye was flipped on the vertical axis to reflect the direction of the right eye (TP = Thinnest point; Kmax = Maximum keratometry of the front surface).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/766163/fbioe-09-766163-HTML-r3/image_m/fbioe-09-766163-g004.jpg)
FIGURE 4. The variation trend of three typical landmark distance parameters with the progress of disease stage (0 = Normal group, 1 = FFKC group, 2 = CKC Stage I group, 3 = CKC Stage II group, 4 = CKC Stage III group, 5 = CKC Stage IV group; DKmax-AP: the absolute distances from the maximum curvature of the anterior surface to the apex of the cornea; DKmax-TP: the absolute distances from the maximum curvature of the anterior surface to the thinnest point; DTP-AP: the absolute distances from the cornea apex (geometric center of the examination [x = 0; y = 0]) to the thinnest point).
Figure 5 represents the correlation between the three geometric landmark distance parameters and the abovementioned nine DCR parameters in each grade from normal to CKC stage IV. From the results, it could be seen that starting from CKC stage II, the distance parameters of the three geometric landmarks were correlated with the DCR parameters (p < 0.0056), and with the aggravation of disease, the correlation between the two parameters of DKmax-AP and DKmax-TP with the DCR parameters was basically strengthened.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/766163/fbioe-09-766163-HTML-r3/image_m/fbioe-09-766163-g005.jpg)
FIGURE 5. Correlation between the three geometric landmark distance parameters and the nine DCR parameters in each grade from normal to CKC Stage IV (the small lattice filled with color indicates that the correlation is statistically significant, in which red indicates that the correlation coefficient is positive, blue indicates that the correlation coefficient is negative, and the depth of color indicates the absolute value of the correlation coefficient).
Discussion
The purpose of the present study was to analyze the changes in coordinates and distances between the three typical geometric landmarks, namely, the thinnest point (TP), the maximum curvature (Kmax), and corneal apex (AP) during the development of keratoconus (Normal, FFKC, and CKC stage I to CKC stage IV), as well as to explore the potential relationship among these changes and the abnormalities of corneal biomechanics.
Our results showed that from normal cornea to CKC stage IV, the coordinates of Kmax gradually tended to focus from the scattered distribution, and finally concentrated on the region under the temporal cornea. However, although the coordinates of Kmax were converging gradually, in fact, its variation range was basically stable within the range of 2 mm in diameter. For TP, as mentioned in other studies (Muftuoglu et al., 2013; Huseynli et al., 2018), it is always concentrated in the inferior temporal region of the cornea, indicating inferior decentration of the thinnest point of the cornea in eyes with early keratoconus. In addition, although there was no significant change in DKmax-AP and DKmax-TP in the process of disease progression, it should be pointed out that DKmax-TP values of each disease grade from FFKC were basically smaller than DKmax-AP, which also illustrated that Kmax might be more inclined to corresponding TP than AP with the development of keratoconus to a certain extent.
Previous studies (Kaya et al., 2007; Ashwin et al., 2009) found that the average distance between the TP and AP in normal eyes was 0.52–1.01 mm, and that in keratoconus eyes was 0.78 mm, which was basically consistent with the results of the current study: normal group (0.72 mm), FFKC group (0.79 mm), and CKC group (0.70 mm). It can be seen from the results of this study that from CKC stage II, the value of DTP-AP would gradually decrease with significant changes, that is, TP would be relatively close to AP in the late stage of disease development. Combined with the discussion in the previous paragraph, it seems that in the late development of keratoconus disease, there will be a trend of proximity or merging between the typical landmark positions of the cornea.
In the present study, there were nine DCR parameters (CBI, SPA1, IR, DA Ratio 2, ARTh, MIR, Radius, Pachy, and DA Ratio 1) that had strong correlation with Rank-group; that is to say, they could well represent the changes in morphological characteristics and biomechanical behavior of cornea during the progress of keratoconus. Studies (Kataria et al., 2019; Yang et al., 2019) pointed out the significant role of these parameters in the diagnosis of keratoconus, which was confirmed by the area under the receiver operating characteristic (ROC) curve. Among the above parameters, in addition to Pachy, ARTh is another DCR parameter that more characterizes the changes of corneal geometric characteristics, which is calculated by the division between corneal thickness at the thinnest point and pachymetric progression index, and a lower value means a faster increase of thickness toward the periphery or a thinner cornea (Vinciguerra et al., 2016). The DA Ratio 1 and DA Ratio 2 are measured at 1 or 2 mm from the center (Wang et al., 2017), and a greater value indicates less resistance to the cornea or a softer cornea. For SPA1, it is developed by using displacement of the apex from the undeformed state to first applanation in the deformation process, and it more characterizes the stiffness behavior of cornea to resist deformation (Vinciguerra et al., 2016). That is, the smaller the SPA1, the smaller the overall stiffness of the cornea. The MIR is the maximum value of radius of curvature during concave phase of the deformation (Yang et al., 2019), while IR is defined as the area under the inverse radius curve with respect to time, and with the progress of keratoconus, both gradually increases, that is, the cornea gradually softens. Of course, the significant trend of the above parameters that can characterize the biomechanical properties of cornea with the increase of disease developmental grades further shows that the cornea will gradually soften during the process of disease progression, that is, the biomechanical properties will gradually weaken (Scarcelli et al., 2014).
According to the results, we found that starting from CKC stage II, the three geometric landmark distance parameters, namely, DKmax-TP, DKmax-AP, and DTP-AP, began to be correlated with the selected biomechanical parameters, and with the aggravation of disease, the correlation between the two parameters (DKmax-TP and DKmax-AP) and the corneal biomechanical parameters was basically strengthened. On this basis, if we synthesize the above discussion, that is, there will be certain proximity or merger between the typical landmark positions of cornea in the later stage of disease development, it is not difficult to believe that the weakening of corneal biomechanical properties may be accompanied by the merger between the typical landmark positions of cornea, especially in the later stage of disease progression.
Of course, our study also had certain limitations. First of all, we did not consider the cone center and the distances in this study because the cone centers of keratoconus were not available in the Pentacam data. Thus, if we can obtain the cone centers of anterior and posterior corneal surfaces by using the Pentacam data, further relevant analysis on the coordinate of the cone center should be carried out in the future research. In addition, this study was not a longitudinal study in the strict sense. The reason why the cross-sectional data were used to study some changes in the development of keratoconus was that the course of keratoconus disease itself was irreversible and would continue to deteriorate over time.
In conclusion, in the later stage of keratoconus, the relationship between the three typical landmark distance parameters and DCR parameters is stronger, and even the weakening of corneal biomechanical properties may be accompanied by the merger of typical landmark positions. It is believed that these findings have certain value for us to further understand the significance of corneal response parameters under external force. Moreover, the coordinate of Kmax may also be used as a reference parameter to judge the disease stage in the future; that is, it will gradually tend to focus from dispersion with the progress of the disease.
Data Availability Statement
The original contributions presented in the study are included in the article/supplementary material. Further inquiries can be directed to the corresponding authors.
Ethics Statement
The studies involving human participants were reviewed and approved by the Ethics Committee of the Beijing Tongren Hospital. Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin.
Author Contributions
LT and HZ contributed equally to this work. HZ, L-LG, XQ, and DZ collected data that met the inclusion criteria. LT conceived the original idea, provided guidance for clinical knowledge, and revised the article. HZ completed the statistical analysis and produced the first draft of the manuscript. LL and YW helped analyze the data and revised the article. YJ and HX-Z helped supervise the project and gave critical revision of article. All authors contributed to manuscript revision, read, and approved the submitted version.
Funding
This work was financially supported by grants from the National Natural Science Foundation of China (32171304 and 82171101); Beijing Nova Program (Z181100006218099); and The Open Research Fund from Beijing Advanced Innovation Center for Big Data-Based Precision Medicine, Beijing Tongren Hospital, Beihang University & Capital Medical University (BHTR-KFJJ-202001).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Ashwin, P. T., Shah, S., Pushpoth, S., Wehbeh, L., and Ilango, B. (2009). The Relationship of Central Corneal Thickness (CCT) to Thinnest Central Cornea (TCC) in Healthy Adults. Contact Lens and Anterior Eye 32 (2), 64–67. doi:10.1016/j.clae.2008.07.006
Catalán-López, S., Cadarso-Suárez, L., López-Ratón, M., and Cadarso-Suárez, C. (2018). Corneal Biomechanics in Unilateral Keratoconus and Fellow Eyes with a Scheimpflug-Based Tonometer. Optom. Vis. Sci. 95 (7), 608–615. doi:10.1097/OPX.0000000000001241
Chan, T. C. Y., Wang, Y. M., Yu, M., and Jhanji, V. (2018). Comparison of Corneal Tomography and a New Combined Tomographic Biomechanical Index in Subclinical Keratoconus. J. Refract Surg. 34 (9), 616–621. doi:10.3928/1081597X-20180705-02
Cui, J., Zhang, X., Hu, Q., Zhou, W.-Y., and Yang, F. (2016). Evaluation of Corneal Thickness and Volume Parameters of Subclinical Keratoconus Using a Pentacam Scheimflug System. Curr. Eye Res. 41 (7), 923–926. doi:10.3109/02713683.2015.1082188
Elham, R., Jafarzadehpur, E., Hashemi, H., Amanzadeh, K., Shokrollahzadeh, F., Yekta, A., et al. (2017). Keratoconus Diagnosis Using Corvis ST Measured Biomechanical Parameters. J. Curr. Ophthalmol. 29 (3), 175–181. doi:10.1016/j.joco.2017.05.002
Galletti, J. D., Ruiseñor Vázquez, P. R., Minguez, N., Delrivo, M., Bonthoux, F. F., Pförtner, T., et al. (2015). Corneal Asymmetry Analysis by Pentacam Scheimpflug Tomography for Keratoconus Diagnosis. J. Refract Surg. 31 (2), 116–123. doi:10.3928/1081597X-20150122-07
Hashemi, H., Khabazkhoob, M., Pakzad, R., Bakhshi, S., Ostadimoghaddam, H., Asaharlous, A., et al. (2019). Pentacam Accuracy in Discriminating Keratoconus from Normal Corneas: A Diagnostic Evaluation Study. Eye Contact Lens 45 (1), 46–50. doi:10.1097/ICL.0000000000000531
Herber, R., Ramm, L., Spoerl, E., Raiskup, F., Pillunat, L. E., and Terai, N. (2019). Assessment of Corneal Biomechanical Parameters in Healthy and Keratoconic Eyes Using Dynamic Bidirectional Applanation Device and Dynamic Scheimpflug Analyzer. J. Cataract Refract Surg. 45 (6), 778–788. doi:10.1016/j.jcrs.2018.12.015
Huseynli, S., Salgado-Borges, J., and Alio, J. L. (2018). Comparative Evaluation of Scheimpflug Tomography Parameters between Thin Non-keratoconic, Subclinical Keratoconic, and Mild Keratoconic Corneas. Eur. J. Ophthalmol. 28 (5), 521–534. doi:10.1177/1120672118760146
Kataria, P., Padmanabhan, P., Gopalakrishnan, A., Padmanaban, V., Mahadik, S., and Ambrósio, R. (2019). Accuracy of Scheimpflug-Derived Corneal Biomechanical and Tomographic Indices for Detecting Subclinical and Mild Keratectasia in a South Asian Population. J. Cataract Refract Surg. 45 (3), 328–336. doi:10.1016/j.jcrs.2018.10.030
Kaya, V., Karakaya, M., Utine, C. A., Albayrak, S., Oge, O. F., and Yilmaz, O. F. (2007). Evaluation of the Corneal Topographic Characteristics of Keratoconus with Orbscan II in Patients with and without Atopy. Cornea 26 (8), 945–948. doi:10.1097/ICO.0b013e3180de1e04
McMonnies, C. W. (2012). Assessing Corneal Hysteresis Using the Ocular Response Analyzer. Optom. Vis. Sci. 89 (3), E343–E349. doi:10.1097/OPX.0b013e3182417223
Meek, K. M., Tuft, S. J., Huang, Y., Gill, P. S., Hayes, S., Newton, R. H., et al. (2005). Changes in Collagen Orientation and Distribution in Keratoconus Corneas. Invest. Ophthalmol. Vis. Sci. 46 (6), 1948–1956. doi:10.1167/iovs.04-1253
Muftuoglu, O., Ayar, O., Hurmeric, V., Orucoglu, F., and Kılıc, I. (2015). Comparison of Multimetric D index with Keratometric, Pachymetric, and Posterior Elevation Parameters in Diagnosing Subclinical Keratoconus in Fellow Eyes of Asymmetric Keratoconus Patients. J. Cataract Refract Surg. 41 (3), 557–565. doi:10.1016/j.jcrs.2014.05.052
Muftuoglu, O., Ayar, O., Ozulken, K., Ozyol, E., and Akıncı, A. (2013). Posterior Corneal Elevation and Back Difference Corneal Elevation in Diagnosing Forme Fruste Keratoconus in the Fellow Eyes of Unilateral Keratoconus Patients. J. Cataract Refract Surg. 39 (9), 1348–1357. doi:10.1016/j.jcrs.2013.03.023
Oxlund, H., and Simonsen, A. H. (1985). Biochemical Studies of normal and Keratoconus Corneas. Acta Ophthalmol. (Copenh) 63 (6), 666–669. doi:10.1111/j.1755-3768.1985.tb01578.x
Peña-García, P., Peris-Martínez, C., Abbouda, A., and Ruiz-Moreno, J. M. (2016). Detection of Subclinical Keratoconus Through Non-Contact Tonometry and the Use of Discriminant Biomechanical Functions. J. Biomech. 49 (3), 353–363. doi:10.1016/j.jbiomech.2015.12.031
Roberts, C. J., Mahmoud, A. M., Bons, J. P., Hossain, A., Elsheikh, A., Vinciguerra, R., et al. (2017). Introduction of Two Novel Stiffness Parameters and Interpretation of Air Puff-Induced Biomechanical Deformation Parameters with a Dynamic Scheimpflug Analyzer. J. Refract Surg. 33 (4), 266–273. doi:10.3928/1081597X-20161221-03
Scarcelli, G., Besner, S., Pineda, R., and Yun, S. H. (2014). Biomechanical Characterization of Keratoconus Corneas Ex Vivo with Brillouin Microscopy. Invest. Ophthalmol. Vis. Sci. 55 (7), 4490–4495. doi:10.1167/iovs.14-14450
Sedaghat, M.-R., Momeni-Moghaddam, H., Ambrósio, R., Heidari, H.-R., Maddah, N., Danesh, Z., et al. (2018). Diagnostic Ability of Corneal Shape and Biomechanical Parameters for Detecting Frank Keratoconus. Cornea 37 (8), 1025–1034. doi:10.1097/ICO.0000000000001639
Tian, L., Ko, M. W. L., Wang, L.-k., Zhang, J.-y., Li, T.-j., Huang, Y.-f., et al. (2014). Assessment of Ocular Biomechanics Using Dynamic Ultra High-Speed Scheimpflug Imaging in Keratoconic and normal Eyes. J. Refract Surg. 30 (11), 785–791. doi:10.3928/1081597X-20140930-01
Vinciguerra, R., Ambrósio, R., Elsheikh, A., Roberts, C. J., Lopes, B., Morenghi, E., et al. (2016). Detection of Keratoconus with a New Biomechanical Index. J. Refract Surg. 32 (12), 803–810. doi:10.3928/1081597X-20160629-01
Wang, Y. M., Chan, T. C. Y., Yu, M., and Jhanji, V. (2017). Comparison of Corneal Dynamic and Tomographic Analysis in Normal, Forme Fruste Keratoconic, and Keratoconic Eyes. J. Refract Surg. 33 (9), 632–638. doi:10.3928/1081597X-20170621-09
Yang, K., Xu, L., Fan, Q., Zhao*, D., and Ren*, S. (2019). Repeatability and Comparison of New Corvis ST Parameters in normal and Keratoconus Eyes. Sci. Rep. 9 (1), 15379. doi:10.1038/s41598-019-51502-4
Glossary
A1DA Deformation amplitude at the first applanation
A1dArcL Chang in arc length (length along the curved line) at first applanation
A1DLA Deflection amplitude at the first applanation
A1DLAr Deflection area at the first applanation
A1DLL Deflection length at the first applanation
A1T Time from the initiation of air puff until the first applanation
A1V The first velocity of applanation
A2DA Deformation amplitude at the second applanation
A2dArcL Chang in arc length (length along the curved line) at second applanation
A2DLA Deflection amplitude at the second applanation
A2DLAr Deflection area at the second applanation
A2DLL Deflection length at the second applanation
A2T Time from the initiation of air puff until the second applanation
A2V The second velocity of applanation
ANOVA One-way analysis of variance
AP Corneal thickness at the apex
ARTh Ambrósio relational thickness to the horizontal profile
bIOP Biomechanical-corrected intraocular pressure
CBI Corvis Biomechanical Index
CH Corneal hysteresis
CKC Clinical keratoconus
Corvis ST Corneal Visualization Scheimpflug Technology
CRF Corneal resistance factor
DA Ratio 1 The maximal value of the ratio between deformation amplitude at the apex and that at 1 mm from the corneal apex
DA Ratio 2 The maximal value of the ratio between deformation amplitude at the apex and that at 2 mm from the corneal apex
dArcLM Delta arc length max
DCR Dynamic corneal response
DKmax-AP The absolute distances from the maximum curvature of the anterior surface to the apex of the cornea
DKmax-TP The absolute distances from the maximum curvature of the anterior surface to the thinnest point
DLAML Max length at deflection amplitude
DLAMT Max time at deflection amplitude
DTP-AP The absolute distances from the cornea apex to the thinnest point
FFKC Forme fruste keratoconus
HCDA Deformation amplitude at the highest concavity
HCdArcL Chang in arc length (length along the curved line) at the highest concavity
HCDLAr Deflection area at the highest concavity
HCDLL Deflection length at the highest concavity
HCT Time from the start until the highest concavity
IOP Intraocular pressure
IR Integrated radius
Kmax Maximum keratometry from the anterior corneal surface
Kmax X The position of maximum keratometry in X direction relative to apex
Kmax Y The position of maximum keratometry in Y direction relative to apex
MIR Max inverse radius
ORA Ocular Response Analyzer
Pachy Pachymetry
PD Peak distance
QS Quality specification
Radius Radius of curvature
SD Standard deviation
SPA1 Stiffness parameter at first applanation
TKC Topographic keratoconus classification
TP Corneal thickness at the thinnest point
TP X The X position of thinnest corneal thickness related to apex
TP Y The Y position of thinnest corneal thickness related to apex
WEMA Max amplitude of whole eye movement
WEMT Max time of whole eye movement
Keywords: keratoconus, forme fruste keratoconus, morphology, biomechanics, geometric landmark
Citation: Tian L, Zhang H, Guo L-L, Qin X, Zhang D, Li L, Wu Y, Jie Y and Zhang H (2021) Distribution of Corneal Geometric Landmarks and Relationship Between Their Distances and Biomechanical Parameters in the Development of Keratoconus. Front. Bioeng. Biotechnol. 9:766163. doi: 10.3389/fbioe.2021.766163
Received: 28 August 2021; Accepted: 21 October 2021;
Published: 22 December 2021.
Edited by:
Yan Wang, Tianjin Eye Hospital, ChinaReviewed by:
Hongli Yang, Devers Eye Institute, United StatesAike Qiao, Beijing University of Technology, China
Copyright © 2021 Tian, Zhang, Guo, Qin, Zhang, Li, Wu, Jie and Zhang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Haixia Zhang, Zhanghx@ccmu.edu.cn; Ying Jie, jie_yingcn@aliyun.com
†These authors have contributed equally to this work and share first authorship.