ORIGINAL RESEARCH article

Front. Bioeng. Biotechnol., 22 December 2021

Sec. Biomechanics

Volume 9 - 2021 | https://doi.org/10.3389/fbioe.2021.766163

Distribution of Corneal Geometric Landmarks and Relationship Between Their Distances and Biomechanical Parameters in the Development of Keratoconus

  • 1. Beijing Institute of Ophthalmology, Beijing Tongren Eye Center, Beijing Tongren Hospital, Capital Medical University and Beijing Ophthalmology and Visual Sciences Key Laboratory, Beijing, China

  • 2. Beijing Advanced Innovation Center for Big Data-Based Precision Medicine, Beihang University & Capital Medical University, Beijing, China

  • 3. School of Biomedical Engineering, Capital Medical University, Beijing, China

  • 4. Department of Medical Engineering, Peking Union Medical College Hospital, Chinese Academy of Medical Sciences, Beijing, China

  • 5. The First People's Hospital of Xuzhou, Jiangsu, China

  • 6. Beijing Key Laboratory of Fundamental Research on Biomechanics in Clinical Application, Capital Medical University, Beijing, China

  • 7. Department of Ophthalmology, Chinese People's Liberation Army General Hospital, Beijing, China

Article metrics

View details

6

Citations

2,6k

Views

1,2k

Downloads

Abstract

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).

TABLE 1

Normal group (n = 142)FFKC group (n = 93)CKC group (n = 290)p-value
Stage I (n = 54)Stage II (n = 123)Stage III (n = 82)Stage IV (n = 31)
Geometric landmarks distances
 DKmax-AP [mm]1.02 (0.06–4.78)1.13 (0.06–5.13)1.12 (0.13–3.03)1.17 ± 0.660.89 (0.13–2.84)0.85 ± 0.590.058a
 DTP-AP [mm]0.72 ± 0.190.76 (0.29–1.92)0.77 ± 0.300.74 (0.06–2.96)0.58 (0.13–2.25)0.41 (0.18–1.26)<0.001a
 DKmax-TP [mm]1.04 (0.01–5.17)1.16 (0.14–5.00)0.98 (0.21–2.52)0.93 (0.08–2.45)0.83 ± 0.400.43 (0.13–1.91)<0.001a
DCR parameters (Corvis ST’s output parameters)
 IOP [mmHg]14.5 (9.0–30.5)13.6 ± 2.113.5 (8.0–30.0)11.5 (5.0–19.0)10.5 ± 2.310.4 ± 3.0<0.001a
 Pachy [μm]543 ± 34518 ± 33502 ± 26479 (402–591)461 ± 32449 ± 37<0.001a
 A1T [ms]7.319 (6.782–8.990)7.172 ± 0.2257.156 (6.590–8.934)6.992 ± 0.2406.857 ± 0.2226.873 ± 0.287<0.001a
 A1V [m/s]0.146 ± 0.0210.158 (0.095–0.196)0.160 ± 0.0210.170 ± 0.0200.181 ± 0.0230.185 ± 0.020<0.001a
 A2T [ms]21.892 ± 0.40021.961 ± 0.37821.902 ± 0.37822.132 ± 0.35022.249 ± 0.33822.260 ± 0.396<0.001b
 A2V [m/s]−0.277 (−0.363 to −0.118)−0.288 (−0.356 to −0.166)−0.297 ± 0.044−0.322 ± 0.046−0.345 (−0.562 to −0.149)−0.348 ± 0.049<0.001a
 HCT [ms]16.863 (15.477–18.108)16.863 (15.246–18.249)17.094 (15.708–17.787)16.863 (15,477–18.018)17.094 (14.784–17.787)16.772 ± 0.4090.098a
 PD [mm]5.194 (4.228–5.768)5.209 ± 0.2685.142 ± 0.2895.268 ± 0.2865.311 ± 0.2785.283 (4.569–5.559)0.009a
 Radius [mm]7.387 ± 0.7666.833 ± 0.7686.316 ± 0.6785.850 (4.152–9.040)5.391 ± 0.6705.090 ± 0.867<0.001a
 A1DA [mm]0.131 ± 0.0110.132 (0.098–0.150)0.136 ± 0.0100.137 (0.105–0.191)0.141 ± 0.0150.154 ± 0.017<0.001a
 HCDA [mm]1.074 ± 0.1101.116 ± 0.1101.138 ± 0.1161.192 (0.943–1.676)1.281 (1.037–1.877)1.298 (1.125–1.756)<0.001a
 A2DA [mm]0.373 (0.244–0.617)0.373 ± 0.0680.381 ± 0.0640.367 (0.226–0.567)0.398 ± 0.0740.394 (0.268–0.702)0.246a
 A1DLL [mm]2.313 (1.835–2.895)2.300 (1.891–2.684)2.319 ± 0.1292.339 ± 0.1432.362 ± 0.1842.394 ± 0.1730.008a
 HCDLL [mm]6.832 (5.614–7.788)6.708 ± 0.4986.575 ± 0.5186.674 ± 0.4606.610 ± 0.4726.380 ± 0.4960.003a
 A2DLL [mm]2.716 (1.660–4.234)2.652 (1.630–4.510)2.608 (1.582–6.441)2.866 ± 0.6422.910 ± 0.6802.877 ± 0.5100.910a
 A1DLA [mm]0.094 ± 0.0080.095 ± 0.0070.099 ± 0.0080.103 (0.079–0.166)0.110 ± 0.0140.120 ± 0.016<0.001a
 HCDLA [mm]0.923 ± 0.1040.965 ± 0.1050.967 ± 0.1161.061 ± 0.1341.122 (0.870–1.734)1.158 ± 0.122<0.001a
 A2DLA [mm]0.103 (0.080–0.142)0.107 (0.080–0.142)0.104 (0.078–0.162)0.115 (0.090–0.250)0.124 (0.083–0.209)0.138 ± 0.017<0.001a
 DLAML [mm]0.938 ± 0.1050.979 ± 0.1110.972 ± 0.1311.063 (0.825–1.558)1.140 (0.905–1.936)1.169 ± 0.120<0.001a
 DLAMT [ms]16.388 (14.637–17.455)16.334 (13.677–17.684)16.436 (14.350–17.396)16.438 (14.535–17.605)16.297 (8.419–17.397)16.219 ± 0.6090.978a
 WEMA [mm]0.270 (0.156–0.504)0.271 ± 0.0660.276 ± 0.0620.260 (0.143–0.445)0.278 (0.155–0.482)0.255 (0.181–0.562)0.943a
 WEMT [ms]21.742 (20.585–23.230)21.862 (20.536–23.701)21.454 ± 1.76021.823 (20.070–23.544)21.997 ± 0.46721.879 ± 0.460<0.001a
 A1DLAr [mm2]0.180 ± 0.0270.176 ± 0.0230.185 ± 0.0230.190 ± 0.0300.206 ± 0.0380.221 ± 0.045<0.001b
 HCDLAr [mm2]3.457 ± 0.5493.577 ± 0.5553.470 ± 0.5903.878 (2.433–6.218)3.981 (2.846–7.195)4.067 ± 0.520<0.001a
 A2DLAr [mm2]0.233 ± 0.0410.234 (0.171–0.369)0.232 (0.167–0.373)0.260 (0.152–0.736)0.280 (0.152–0.534)0.307 ± 0.056<0.001a
 A1dArcL [mm]−0.019 (−0.029 to −0.008)−0.019 (−0.022 to −0.005)−0.019 ± 0.003−0.020 ± 0.004−0.021 (−0.040 to −0.012)−0.025 ± 0.007<0.001a
 HCdArcL [mm]−0.135 ± 0.028−0.127 ± 0.024−0.115 ± 0.020−0.119 (-0.230 to −0.014)−0.115 (−0.231 to −0.019)−0.115 ± 0.032<0.001a
 A2dArcL [mm]−0.022 (−0.036 to −0.004)−0.023 (−0.040 to −0.016)−0.023 (-0.036 to -0.016)−0.025 (−0.057 to 0.006)−0.029 ± 0.008-0.031 (-0.051 to 0.002)<0.001a
 dArcLM [mm]−0.153 ± 0.035−0.148 ± 0.030−0.146 ± 0.064−0.142 (−0.452 to −0.034)−0.135 (−0.399 to −0.041)−0.137 (−0.270 to -0.050)0.023a
 MIR [mm−1]0.167 ± 0.0160.178 (0.140–0.259)0.199 ± 0.0240.209 ± 0.0260.223 (0.170–0.336)0.251 ± 0.040<0.001a
 DA Ratio 24.188 (3.081–5.946)4.646 ± 0.4404.880 ± 0.5965.409 ± 0.7866.208 ± 0.8616.615 ± 1.013<0.001a
 PachySlope [μm]50.435 (29.036–98.805)50.480 ± 10.24755.290 ± 13.33565.865 (−24.281–145.630)76.817 (38.351–154.705)98.906 ± 24.373<0.001a
 DA Ratio 11.558 (1.407–1.729)1.593 ± 0.0441.613 ± 0.0521.650 (1.509–1.918)1.710 (1.560–1.971)1.744 ± 0.074<0.001a
 ARTh415.883 ± 86.639392.170 ± 92.140322.116 (195.308–603.057)226.520 (28.989–630.625)174.476 (71.233–417.461)133.123 ± 56.762<0.001a
 bIOP [mmHg]14.7 (10.8–26.5)14.3 ± 1.914.6 (10.4–29.2)13.3 ± 2.312.4 ± 2.212.6 ± 2.9<0.001a
 IR [mm−1]7.918 (5.300–12.555)8.898 ± 1.0259.698 ± 1.19210.858 ± 1.71611.967 (8.085–17.603)13.401 ± 2.298<0.001a
 SPA1102.040 ± 19.50585.088 ± 16.69277.961 ± 18.80459.822 ± 15.63947.491 ± 12.51241.476 ± 14.298<0.001b
 CBI0.065 (0.000–1.000)0.404 (0.001–1.000)0.957 (0.000–1.000)1.000 (0.001–1.000)1.000 (0.943–1.000)1.000 (0.970–1.000)<0.001a

Comparison of distance of three geometric landmarks and DCR parameters by groups.

a

Non-parametric Kruskal–Wallis test.

b

One-way ANOVA test.

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.

FIGURE 1

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).

FIGURE 2

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.

FIGURE 3

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).

FIGURE 4

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.

FIGURE 5

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.

Statements

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.

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

References

  • 1

    AshwinP. T.ShahS.PushpothS.WehbehL.IlangoB. (2009). The Relationship of Central Corneal Thickness (CCT) to Thinnest Central Cornea (TCC) in Healthy Adults. Contact Lens and Anterior Eye32 (2), 6467. 10.1016/j.clae.2008.07.006

  • 2

    Catalán-LópezS.Cadarso-SuárezL.López-RatónM.Cadarso-SuárezC. (2018). Corneal Biomechanics in Unilateral Keratoconus and Fellow Eyes with a Scheimpflug-Based Tonometer. Optom. Vis. Sci.95 (7), 608615. 10.1097/OPX.0000000000001241

  • 3

    ChanT. C. Y.WangY. M.YuM.JhanjiV. (2018). Comparison of Corneal Tomography and a New Combined Tomographic Biomechanical Index in Subclinical Keratoconus. J. Refract Surg.34 (9), 616621. 10.3928/1081597X-20180705-02

  • 4

    CuiJ.ZhangX.HuQ.ZhouW.-Y.YangF. (2016). Evaluation of Corneal Thickness and Volume Parameters of Subclinical Keratoconus Using a Pentacam Scheimflug System. Curr. Eye Res.41 (7), 923926. 10.3109/02713683.2015.1082188

  • 5

    ElhamR.JafarzadehpurE.HashemiH.AmanzadehK.ShokrollahzadehF.YektaA.et al (2017). Keratoconus Diagnosis Using Corvis ST Measured Biomechanical Parameters. J. Curr. Ophthalmol.29 (3), 175181. 10.1016/j.joco.2017.05.002

  • 6

    GallettiJ. D.Ruiseñor VázquezP. R.MinguezN.DelrivoM.BonthouxF. F.PförtnerT.et al (2015). Corneal Asymmetry Analysis by Pentacam Scheimpflug Tomography for Keratoconus Diagnosis. J. Refract Surg.31 (2), 116123. 10.3928/1081597X-20150122-07

  • 7

    HashemiH.KhabazkhoobM.PakzadR.BakhshiS.OstadimoghaddamH.AsaharlousA.et al (2019). Pentacam Accuracy in Discriminating Keratoconus from Normal Corneas: A Diagnostic Evaluation Study. Eye Contact Lens45 (1), 4650. 10.1097/ICL.0000000000000531

  • 8

    HerberR.RammL.SpoerlE.RaiskupF.PillunatL. E.TeraiN. (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), 778788. 10.1016/j.jcrs.2018.12.015

  • 9

    HuseynliS.Salgado-BorgesJ.AlioJ. L. (2018). Comparative Evaluation of Scheimpflug Tomography Parameters between Thin Non-keratoconic, Subclinical Keratoconic, and Mild Keratoconic Corneas. Eur. J. Ophthalmol.28 (5), 521534. 10.1177/1120672118760146

  • 10

    KatariaP.PadmanabhanP.GopalakrishnanA.PadmanabanV.MahadikS.AmbrósioR.Jr (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), 328336. 10.1016/j.jcrs.2018.10.030

  • 11

    KayaV.KarakayaM.UtineC. A.AlbayrakS.OgeO. F.YilmazO. F. (2007). Evaluation of the Corneal Topographic Characteristics of Keratoconus with Orbscan II in Patients with and without Atopy. Cornea26 (8), 945948. 10.1097/ICO.0b013e3180de1e04

  • 12

    McMonniesC. W. (2012). Assessing Corneal Hysteresis Using the Ocular Response Analyzer. Optom. Vis. Sci.89 (3), E343E349. 10.1097/OPX.0b013e3182417223

  • 13

    MeekK. M.TuftS. J.HuangY.GillP. S.HayesS.NewtonR. H.et al (2005). Changes in Collagen Orientation and Distribution in Keratoconus Corneas. Invest. Ophthalmol. Vis. Sci.46 (6), 19481956. 10.1167/iovs.04-1253

  • 14

    MuftuogluO.AyarO.HurmericV.OrucogluF.KılıcI. (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), 557565. 10.1016/j.jcrs.2014.05.052

  • 15

    MuftuogluO.AyarO.OzulkenK.OzyolE.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), 13481357. 10.1016/j.jcrs.2013.03.023

  • 16

    OxlundH.SimonsenA. H. (1985). Biochemical Studies of normal and Keratoconus Corneas. Acta Ophthalmol. (Copenh)63 (6), 666669. 10.1111/j.1755-3768.1985.tb01578.x

  • 17

    Peña-GarcíaP.Peris-MartínezC.AbboudaA.Ruiz-MorenoJ. M. (2016). Detection of Subclinical Keratoconus Through Non-Contact Tonometry and the Use of Discriminant Biomechanical Functions. J. Biomech.49 (3), 353363. 10.1016/j.jbiomech.2015.12.031

  • 18

    RobertsC. J.MahmoudA. M.BonsJ. P.HossainA.ElsheikhA.VinciguerraR.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), 266273. 10.3928/1081597X-20161221-03

  • 19

    ScarcelliG.BesnerS.PinedaR.YunS. H. (2014). Biomechanical Characterization of Keratoconus Corneas Ex Vivo with Brillouin Microscopy. Invest. Ophthalmol. Vis. Sci.55 (7), 44904495. 10.1167/iovs.14-14450

  • 20

    SedaghatM.-R.Momeni-MoghaddamH.AmbrósioR.JrHeidariH.-R.MaddahN.DaneshZ.et al (2018). Diagnostic Ability of Corneal Shape and Biomechanical Parameters for Detecting Frank Keratoconus. Cornea37 (8), 10251034. 10.1097/ICO.0000000000001639

  • 21

    TianL.KoM. W. L.WangL.-k.ZhangJ.-y.LiT.-j.HuangY.-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), 785791. 10.3928/1081597X-20140930-01

  • 22

    VinciguerraR.AmbrósioR.JrElsheikhA.RobertsC. J.LopesB.MorenghiE.et al (2016). Detection of Keratoconus with a New Biomechanical Index. J. Refract Surg.32 (12), 803810. 10.3928/1081597X-20160629-01

  • 23

    WangY. M.ChanT. C. Y.YuM.JhanjiV. (2017). Comparison of Corneal Dynamic and Tomographic Analysis in Normal, Forme Fruste Keratoconic, and Keratoconic Eyes. J. Refract Surg.33 (9), 632638. 10.3928/1081597X-20170621-09

  • 24

    YangK.XuL.FanQ.Zhao*D.Ren*S. (2019). Repeatability and Comparison of New Corvis ST Parameters in normal and Keratoconus Eyes. Sci. Rep.9 (1), 15379. 10.1038/s41598-019-51502-4

Summary

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

Volume

9 - 2021

Edited by

Yan Wang, Tianjin Eye Hospital, China

Reviewed by

Hongli Yang, Devers Eye Institute, United States

Aike Qiao, Beijing University of Technology, China

Updates

Copyright

*Correspondence: Haixia Zhang, ; Ying Jie,

†These authors have contributed equally to this work and share first authorship.

This article was submitted to Biomechanics, a section of the journal Frontiers in Bioengineering and Biotechnology

Disclaimer

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics