Predictive Modeling of Long-Term Glaucoma Progression Based on Initial Ophthalmic Data and Optic Nerve Head Characteristics

Purpose: The purpose of this study was to develop a model, based on initial optic nerve head (ONH) characteristics, predictive of long-term rapid retinal nerve fiber layer (RNFL) thinning in patients with open-angle glaucoma (OAG). Methods: This study evaluated 712 eyes with OAG that had been followed up for >5 years with annual evaluation of RNFL thickness. Baseline ophthalmic features were incorporated into the machine learning models for prediction of faster RNFL thinning. The model was trained and tested using a random forest (RF) method, and was interpreted using Shapley additive explanations. Factors associated with faster rate of RNFL thinning were statistically evaluated using a decision tree. Results: The RF model showed that greater lamina cribrosa (LC) curvature, higher intraocular pressure (IOP), visual field mean deviation converging towards −5 dB, and thinner peripapillary choroid at baseline were the four most significant features predicting faster RNFL thinning. Partial interaction between the features showed that larger LC curvature was a strong factor for faster RNFL thinning when it exceeded approximately 12.0. When the LC curvature was ≤12, higher initial IOP and thinner peripapillary choroid played a role in the rapid RNFL thinning. Based on the decision tree, higher IOP (>26.5 mm Hg), greater laminar curvature (>13.95), and thinner peripapillary choroid (≤117.5 μm) were the 3 most important determinants affecting the rate of RNFL thinning. Conclusions: Baseline ophthalmic data and ONH characteristics of patients with OAG were predictive of eyes at risk of faster progression. Combinations of important characteristics, such as IOP, LC curvature, and choroidal thickness, could stratify eyes into groups with different rates of RNFL thinning. Translational Relevance: This work lays the foundations for developing prediction models to estimate glaucoma prognosis based on initial ONH characteristics.