ABSTRACT
Objectives To find predictive factors and to construct predictive models using epidemiological and clinical preoperative factors for the visual acuity change after intracorneal ring segment (ICRS) implantation in keratoconus patients.
Methods Medical records of 287 keratoconic eyes implanted with ICRS at Chula Refractive Surgery Centre of a tertiary university hospital (Bangkok, Thailand) between January 2012 and March 2022 were retrospectively reviewed for epidemiological and clinical preoperative variables, including those derived from Scheimpflug tomography. The correlation between these variables and postoperative uncorrected and corrected distance visual acuity change (ΔUDVA and ΔCDVA; logMAR scale) at 6 months was explored. Two hundred forty eight eyes (excluding cases with unmeasurable refraction) were randomized into 2 groups: the equation group (198 eyes) and the validation group (50 eyes). Multiple linear regression analysis was used to develop the predictive model from the equation group.
Results Twenty preoperative variables were statistically correlated with ΔCDVA. Only the preoperative corrected distance visual acuity (CDVAp) strongly correlated with ΔCDVA (R = - 0.746). ΔUDVA correlated with only preoperative visual acuity, corneal astigmatism, and maximum Ambrósio relational thickness (ARTmax). CDVAp, spherical power (SPH), and anterior minimum sagittal curvature (Rmin) were the best predictors of ΔCDVA. The proposed model, ΔCDVA = 0.589 - 0.713×CDVAp - 0.010×SPH - 0.082×Rmin, had an acceptable predictability (R2 = 53.4%). The prediction was correct in 82% of the eyes within 0.22 logMAR.
Conclusions Potential predictive factors and models for ICRS-induced changes in visual acuity were proposed as adjunctive tools for clinicians. Such tools could be used for case selection and during counseling before ICRS implantation to maximize surgical outcomes.
1. INTRODUCTION
Keratoconus is an ectatic corneal disorder characterized by progressive biomechanical corneal instability. This can lead to corneal apex thinning, corneal protrusion, irregular astigmatism, and—sometimes—central corneal scarring [1]. Keratoconus is generally bilateral and asymmetrical. Typical onset is during puberty, with progression for around 40 years. Spectacles and rigid contact lenses were the initial treatments used for keratoconus. With progression, this condition can cause irreversible visual loss, leading to corneal transplantation in 10–20% of cases [2]. A reversible and less-invasive surgical option than corneal transplantation is intrastromal corneal ring segment (ICRS) implantation. The ICRS exerts an arc-shortening effect on the corneal architecture, flattening the cornea while improving astigmatism and contact lens tolerance [3–5]. Several studies have emphasized the advantages of ICRS, including its removability, stability, and security, by eliminating the need for an intraocular procedure [6–9].
The corneal response to ICRS in keratoconus remains unpredictable. ICRS improves both uncorrected (UDVA) and corrected visual acuity (CDVA) [9–12]. However, some patients, particularly those whose preoperative corrected visual acuity (CDVAp) was not compromised, experienced worse vision post-implantation [13]. Furthermore, postoperative changes in vision do not regularly correlate with geometric improvements in the cornea [14]. Numerous studies have attempted to determine the relationship between various factors and surgical outcomes to avoid unsatisfactory results [4, 7, 11, 13, 15–26]. However, the impact of these factors on ICRS-induced visual outcomes remains inconclusive, and some factors reportedly increase the risk of post-implantation complications. For example, atopic dermatitis has been associated with increased ICRS extrusion [27].
In this study, we aimed to identify predictors of visual outcomes at 6 months (or within 5 – 12 months) following ICRS implantation. We also created a mathematical model using the abovementioned factors to predict postoperative visual acuities quantitatively. In the future, this model may facilitate patient selection and aid preoperative patient education efforts.
2. MATERIALS AND METHODS
2.1 Study subjects
This study protocol was approved by the Institutional Review Board (IRB) of the Faculty of Medicine, Chulalongkorn University, Thailand (COA no. 1517/2022 and IRB no. 525/63) and adhered to the tenets of the Declaration of Helsinki. The trial was registered with the Thai Clinical Trial Registry (TCTR20200929001).
This retrospective study was conducted during November 2022 by collecting data from medical records of patients who underwent ICRS implantation at the Chula Refractive Surgery Center, King Chulalongkorn Memorial Hospital, Bangkok, Thailand. Although the authors had access to information that could identify individual patients, the data in the record forms were anonymized to maintain patients’ confidentiality. The primary objective was to identify predictors of postoperative visual acuity change 6 months (or within 5-12 months) following ICRS implantation in patients with keratoconus. The medical records of all eligible cases (351 eyes) were included in this study. The inclusion criteria were as follows: patients with keratoconus implanted with ICRS (Ferrara Ring, AJL, Boecillo, Spain) at the Chula Refractive Surgery Centre between January 2012 and March 2022. Keratoconus diagnosis was based on characteristic signs noted on slit-lamp examination, corneal topography, and tomography. The patient with contact lens intolerance and without contraindication of the ICRS implantation (such as central corneal scar) became a candidate for the surgery. The exclusion criteria were as follows: (i) lack of visual acuity measurement during 5 - 12 months after ICRS implantation (39 eyes), (ii) absence of preoperative data (14 eyes), and (iii) vision-affecting complication and/or ICRS removal within the follow-up period (11 eyes). Ultimately, 64 eyes were excluded from the study, yielding a study cohort of 287 eyes.
2.2 Surgical technique
All included subjects underwent ICRS implantation by one of five corneal specialists using femtosecond laser-assisted tunnel creation under topical anesthesia. The arc length, thickness, and number of ring segments were selected according to a previously described nomogram based on the type of keratoconus, the steepest axis, Q-value, and anterior corneal topographic astigmatism. A disposable suction ring was placed and centered after marking the visual axis. Next, an intrastromal tunnel was created by a laser beam with a spot size of 3 µm while focused on a predetermined depth (75% of corneal thickness) from the anterior corneal surface at a 5.0 optical zone. After creating the main incision, the ICRS was implanted using the complete aseptic technique. Postoperative medications included topical steroids and antibiotics.
2.3 Data collection and outcome measurement
The epidemiological, pre-implantation, and post-implantation data were obtained from patients’ medical records. The epidemiological data included each patient’s sex, the age at onset, age at ICRS implantation, and history of atopy and eye rubbing. Pre-implantation assessments included uncorrected and corrected distance visual acuities (UDVAp and CDVAp; logMAR), intraocular pressure (IOP; mmHg) by applanation/noncontact tonometer, spherical power (SPH; D) and cylindrical power (CYL; D) from autorefractor, and corneal data evaluated using Scheimpflug tomography (Pentacam HR, Oculus, Wetzlar, Germany). The Scheimpflug tomography-derived variables included the mean keratometry of the front (Km-front; D) and back surfaces (Km-back; D), maximum keratometry (Kmax; D), anterior corneal astigmatism (CA), Q-value of the front (Q-front) and the back surface (Q-back), pachymetry (microns), topographic indices, pachymetric indices, and corneal aberration. Six distances between each pair of four topographic landmarks (apex, pupil center, thinnest point, and Kmax point) calculated by the Pythagorean Theorem were also analyzed. The aberrometry (from Zernike analysis) comprised the total root mean square (total rms), root mean square of lower-order (LOA), and higher-order aberration (HOA). Moreover, the K-factor (KF), which was previously reported as a potential predictor of ΔCDVA, was also determined in this study [21]. The KF was indirectly derived from the Scheimpflug tomography by calculating the multiple products of flattest keratometry and corneal astigmatism of the front corneal surface. The post-implantation visual outcome was assessed by determining changes in the UDVA and CDVA at six months (or within 5-12 months) post-surgery and comparing them to preoperative values (ΔUDVA and ΔCDVA). Holladay’s technique was applied to transform visual acuity into the LogMAR scale for cases where visual acuity was reported as counting fingers, hand motion, light perception, or no light perception. The complications from ICRS implantation were retrieved as secondary outcomes.
2.4 Statistical analyses
Statistical analyses were performed using SPSS Statistics for Windows, version 21 (IBM Corporation, Armonk, NY, USA). Considering the large sample size, the data distribution was assumed to be normal. In comparison between the equation group and validation group, Fisher’s exact test and unpaired t-test were used for categorical data and continuous data, respectively. Paired t-test was used to compare pre- and postoperative visual acuity.
Univariate analysis was initially performed to evaluate the correlation between each preoperative variable and the visual outcomes (ΔUDVA and ΔCDVA) using the Pearson correlation coefficient (R). Dummy variables were applied for categorical variables. Female was coded as 0 while male was coded as 1. The history of atopy and eye rubbing was coded as 0 for absence and 1 for presence. The variable with p < 0.2 in this step was considered as a potential predictor and selected for the multivariate regression analysis. The eyes which missed the value of those potential predictors would be excluded from the regression study. In this step, 39 eyes with unmeasurable refractive errors were excluded. The remaining 248 eyes were randomly divided into 2 groups: 198 eyes were assigned to a group for developing the predictive equation (equation group) and the other 50 eyes to a group for validating the proposed equation (validation group).
A stepwise method was applied to construct the predictive model. Model assumptions were evaluated by analyzing the Durbin-Watson test (to confirm the lack of correlation between errors), mean Cook’s distance (to detect influential points or outliers), collinearity tolerance, and the variance inflation factor (VIF). The level of significance was set at p < 0.05.
3. RESULTS
Among 287 eyes of 230 patients was comprised of 154 (53.7%) right and 133 (46.3%) left eyes. There were 161 males (70.0%) and 69 females (30.0%). The mean age at onset was 20.11 ± 7.86 years (range: 7–50 years), while the mean age at ring implantation time was 26.97 ± 8.56 years (range: 13–59 years). There were 149 (51.9%) eyes with atopy and 177 (61.7%) with frequent eye rubbing. Table 1 summarizes the preoperative numerical data’s means, standard deviations, and ranges. The mean duration of postoperative measurement was 7.0 ± 1.4 months (5.0–12.0 months). We observed significant improvement in the UDVA and CDVA from 1.08 ± 0.55 logMAR to 0.73 ± 0.44 logMAR (Snellen equivalent from 20/239 to 20/109) and 0.51 ± 0.42 logMAR to 0.35 ± 0.28 logMAR (Snellen equivalent from 20/65 to 20/44), respectively. The means of ΔUDVA and ΔCDVA were -0.34 ± 0.48 (95% CI: - 0.40 to -0.29, p < 0.001) and -0.17 ± 0.35 (95% CI: -0.21 to -0.13, p < 0.001), respectively.
The patient demographic data of the equation group and validation group were presented in Table 2. showing the similarity between these two groups.
3.1 Univariate linear regression analysis
Table 3 summarizes the Pearson correlation coefficients between the postoperative outcomes (ΔUDVA and ΔCDVA) and each preoperative variable. Neither ΔUDVA nor ΔCDVA significantly correlated with the epidemiological variables. Four variables were significantly correlated with ΔUDVA: UDVAp (R = -0.643, p < 0.001), CDVAp (R = -0.337, p < 0.001), CA (R = 0.119, p = 0.043) and ARTmax (R = 0.118, p = 0.046).
Twenty preoperative variables were statistically correlated with ΔCDVA. Only CDVAp demonstrated a strong correlation with ΔCDVA (R = -0.746). In contrast, the others had a weak correlation (|R| < 0.40), including UDVAp, Km-front, Km-back, KF, Kmax, Q-front, Q-back, thickness at pupil center (CPT), thickness at the apex (CAT), thinnest pachymetry (CTT), Belin/Ambrósio Enhanced Ectasia Display (BADD), average pachymetric progression index (PIavg), ARTmax, the index of surface variance (ISV), central keratoconus index (CKI), Rmin, total rms, LOA, and HOA.
3.2 Multivariable linear regression analysis
3.2.1 Predictive model for ΔUDVA after six months
According to univariate analysis, twelve preoperative variables with p < 0.20 were introduced to the multiple regression analysis. By applying the stepwise method, UDVAp, and PIavg were selected as explanatory variables of ΔUDVA. (summary of the result is shown in Table 4). The proposed equation was
The predictability (R2) was 0.408, and the adjusted R2 was 0.402, with F = 67.101 (p < 0.001). No influential points or outliers were detected (mean Cook’s distance, 0.006 ± 0.014). The independence of the residuals (Durbin-Watson test = 2.079) and the lack of multicollinearity were confirmed.
3.2.2 Predictive model for ΔCDVA after six months
Twenty-eight preoperative variables with p < 0.20 were recruited for the multiple regression analysis. Using the stepwise method, CDVAp, SPH, and Rmin were selected as explanatory variables of ΔCDVA. (summary of the result is shown in Table 5). The proposed equation was
The predictability (R2) was 0.541, and the adjusted R2 was 0.534, with F = 76.194 (p < 0.001). No influential points or outliers were detected (mean Cook’s distance, 0.010 ± 0.040). The independence of the residuals (Durbin-Watson test = 2.224) and the lack of multicollinearity were confirmed. Patients with poorer preoperative visual acuity or higher CDVAp logMAR value demonstrated greater improvement in visual acuity after ring implantation. The model indicated that patients with more negative SPH (more myopic) and smaller Rmin would likely demonstrate poor visual gain when adjusted for CDVAp.
3.3 Validation of the predictive models
For the predictive model of ΔUDVA, the average absolute residual between predicted UDVA and actual UDVA was 0.294 ± 0.238, and the correct prediction was achieved only in 44% of cases within 0.22 logMAR.
For the predictive model of ΔCDVA, the average absolute residual between predicted CDVA and actual CDVA was 0.159 ± 0.159, and the correct prediction was achieved in 82% of cases within 0.22 logMAR
3.4 Complications of ICRS implantation
Both intraoperative and postoperative complications among the 351 eyes were recorded as secondary outcomes. There was one intraoperative complication in which a corneal perforation at the incision site using a stromal spreader was documented. This eye was successfully re-operated using a manual technique, approximately seven months following the first surgery. Postoperative complications included segment migration (17 eyes, 4.84%), ICRS-related infection (14 eyes, 3.99%), corneal melting (3 eyes, 0.85%), peri-annular deposits (2 eyes, 0.56%) and spontaneous fragmentation of the ring segment (1 eye, 0.28%). Twelve eyes underwent ring removal, four of which were because of infection-induced corneal melting. One patient requested ring explantation from his right eye owing to concerns about corneal irritation and a mild visual acuity drop post-surgery from CDVA 20/16 to 20/25.
4. DISCUSSION
According to our literature review, this may be the largest retrospective study of the relationship between ICRS-induced visual acuity changes and various demographic and preoperative factors of patients with keratoconus. The diversity of studied variable types was another strength of our study. This study was also the first attempt to create an equation to predict ΔUDVA following ICRS implantation. In the multiple linear regression analysis, UDVAp and PIavg seemed to be good predictors for ΔUDVA. However, this model had low predictability (Adjust R2 = 0.402) and low validity (accuracy rate less than 50%).
Regarding our univariate analysis, patients with poorer baseline visual acuity (larger value on the logMAR scale) were anticipated to demonstrate greater UDVA and CDVA improvement. This finding was consistent with those of previous studies [13, 15, 18, 21]. More myopic spherical power weakly correlated with more considerable CDVA gain (R = 0.123). This was in contrast to a study by Alio et al. [16], which found that patients with better visual acuity gain had less myopia than other groups. Several parameters derived from Scheimpflug-based tomography demonstrated significant correlations with ΔCDVA; however, only corneal astigmatism significantly correlated with ΔUDVA. Based on the Pearson coefficients, patients with more advanced baseline topography (steeper Km, K factor, and Kmax; and more prolate Q-value) tended to gain more lines of postoperative CDVA. However, our findings differed from previous studies, which found that patients with lower mean keratometry experienced more visual improvement [7, 16, 24].
The corneal thickness at the apex, thinnest point, and pupil center were significantly correlated with ΔCDVA in a similar direction, i.e., thinner corneal thickness tended to demonstrate more CDVA improvement. This finding contrasts that of Zare et al. [7], who stated that a patient with a thicker cornea (thinnest corneal thickness > 400 µm) would demonstrate greater improvements in UDVA and CDVA six months post-surgery. The relationships between ΔCDVA and all pachymetric indices were also shown in our study. ICRS-induced visual gain increased in patients with poorer pachymetric indices (higher preoperative BADD, higher PIavg, and lower ARTmax).
For topographic indices, higher ISV and CKI values and smaller Rmin significantly correlated with a more remarkable improvement in CDVA. The ISV is calculated from the standard deviation of individual sagittal radii from mean curvature. Higher ISV is an indicator of increased anterior corneal surface irregularity. CKI is the ratio between mean radius of curvature in the peripheral zone and mean radius of curvature in central zone, while Rmin denotes the maximum steepness of the cone [28]. As the CKI is higher, the cone of keratoconic cornea become far steeper, and Rmin usually become smaller. It can be inferred that the cornea with a more irregular surface and the steeper cone will gain a more favorable visual outcome from ICRS implantation. Nonetheless, a previous study had reported in the opposite direction that lower preoperative ISV and larger Rmin had a strong correlation with a better gain of UDVA and CDVA [24].
All aberrometry variables had a weak correlation with CDVA change in our study. The highly aberrated cornea was related to more gain of CDVA. However, many previous studies failed to demonstrate a significant correlation between aberrometry and visual acuity change after ICRS implantation [15, 22, 23].
Distances between the corneal apex, thinnest point, pupil center, and Kmax point were also studied for the hypothesis that longer distances might reflect more cone eccentricity and have more effect on visual acuity. However, this study failed to show any significant correlation between these distances and visual acuity change. Gatzioufas et al. also studied the association between visual outcome after ring implantation with distance from apex to thinnest point and from apex to maximum keratometry. They found that there was no statistically significant association [17].
In this study, the age at onset, age at surgery, history of rubbing or atopy, or intraocular pressure did not significantly correlate with ICRS-induced visual acuity change. Despite females experiencing more gain of UDVA than males in the univariate analysis, this association was disproven in the multivariate analysis, similar to previous results [19].
According to the result of univariate analysis in our study, CDVA improvement was anticipated to be better in more advanced keratoconus because it is usually associated with poorer baseline visual acuity, higher myopia, steeper keratometry, more prolateness, thinner pachymetry, higher surface irregularity, steeper cone, and higher aberration.
Our proposed model for predicting ΔCDVA included different predictors from those previously reported [21, 22]. Three variables, including CDVAp, SPH, and Rmin, were selected by the multivariate analysis in the predictive equation for ΔCDVA. This model corresponded to the findings of Pena-Garcia et al. [21] where CDVAp emerged as the most critical predictor of ΔCDVA at six months. While no keratometric parameters appeared in our model, the previously proposed model included at least one keratometric parameter [21, 22]. This may be attributed to different sets of potential variables in the multiple regression analysis and differences in ICRS manufacturers in each study. It could be noted that the sign of coefficients of SPH and Rmin had reversed from plus (in univariate analysis) to minus. This revealed that if two patients have equal CDVAp, one with less myopia and larger Rmin will be predicted to gain more CDVA improvement.
This study was limited by its sole focus on visual acuity at about six months post-surgery, as this time point appeared to be associated with short-term outcomes. Thus, further research with a longer follow-up is necessary to identify factors that predict the long-term effects of ICRS. Another limitation was a high rate of missing data due to the retrospective nature of this study. We excluded some cases because their medical records lacked pre-defined parameters; therefore, our results may not reflect real-world outcomes for all patients undergoing ICRS. The utility of the proposed models might be limited because its variables contained the spherical power from autorefraction, which is frequently unmeasurable in advanced keratoconus. Researchers are also interested in conducting studies on corneal biomechanical properties evaluated by the dynamic Scheimpflug analyzer (Corvis ST, Oculus, Wetzlar, Germany). This is based on the hypothesis that biomechanical properties of the cornea can be potential predictors for ΔCDVA following ICRS implantation [22]. Unfortunately, we could not perform the aforementioned study because a few patients underwent the dynamic Scheimpflug analysis before the ICRS surgery at our center. In conclusion, the baseline visual acuity, refractive error, and numerous parameters from Scheimpflug tomography can predict ICRS-induced changes in visual acuity six months post-surgery. These predictors, alongside the proposed model, could assist with appropriate case selection and correct timing for ICRS implantation. Considering these parameters may maximize visual gains for patients with keratoconus.
Data Availability
All relevant data are within the manuscript and its Supporting Information files.
SUPPORTING INFORMATION
S1 File. The data spreadsheet (XLSX).
ACKNOWLEDGEMENTS
We would like to thank Dollapas Punpanich, M.S., a consultant for statistical analysis, and Editage (www.editage.com) for English language editing.
Footnotes
Declarations of interest: None of the authors had any financial or proprietary interest to declare.
Financial support None to disclose.