Individual Prognostication of Disease Activity and Disability Worsening in Multiple Sclerosis With Retinal Layer Thickness z Scores

Retinal optical coherence tomography (OCT) provides promising prognostic imaging biomarkers for future disease activity in multiple sclerosis (MS). However, raw OCT-derived measures have multiple dependencies, supporting the need for establishing reference values adjusted for possible confounders. The purpose of this study was to investigate the capacity for age-adjusted scores of OCT-derived measures to prognosticate future disease activity and disability worsening in people with MS (PwMS). We established age-adjusted OCT reference data using generalized additive models for location, scale, and shape for peripapillary retinal nerve fiber layer (pRNFL) and ganglion cell-inner plexiform layer (GCIP) thicknesses, involving 910 and 423 healthy eyes, respectively. Next, we transformed the retinal layer thickness of PwMS from 3 published studies into age-adjusted scores (pRNFL-z and GCIP-z) based on the reference data. Finally, we investigated the association of pRNFL-z or GCIP-z as predictors with future confirmed disability worsening (Expanded Disability Status Scale score increase) or disease activity (failing of the no evidence of disease activity [NEDA-3] criteria) as outcomes. Cox proportional hazards models or logistic regression analyses were applied according to the original studies. Optimal cutoffs were identified using the Akaike information criterion as well as location with the log-rank and likelihood-ratio tests. In the first cohort (n = 863), 172 PwMS (24%) had disability worsening over a median observational period of 2.0 (interquartile range [IQR]:1.0-3.0) years. Low pRNFL-z (≤-2.04) were associated with an increased risk of disability worsening (adjusted hazard ratio (aHR) [95% CI] = 2.08 [1.47-2.95], 3.82e ). In the second cohort (n = 170), logistic regression analyses revealed that lower pRNFL-z showed a higher likelihood for disability accumulation at the two-year follow-up (reciprocal odds ratio [95% CI] = 1.51[1.06-2.15], = 0.03). In the third cohort (n = 78), 46 PwMS (59%) did not maintain the NEDA-3 status over a median follow-up of 2.0 (IQR: 1.9-2.1) years. PwMS with low GCIP-z (≤-1.03) had a higher risk of showing disease activity (aHR [95% CI] = 2.14 [1.03-4.43], = 0.04). Compared with raw values with arbitrary cutoffs, applying the score approach with optimal cutoffs showed better performance in discrimination and calibration (higher Harrell's concordance index and lower integrated Brier score). In conclusion, our work demonstrated