Sage Journals: Discover world-class research

Abstract

Background:

While disability progression in multiple sclerosis (MS) is age-dependent, the exact timing of neuroaxonal degeneration throughout the patient lifespan remains unclear.

Objectives:

To investigate the influence of age on retinal neurodegeneration after the first unilateral episode of optic neuritis (ON).

Methods:

We measured peripapillary retinal nerve fibre layer (pRNFL) and macular ganglion cell inner plexiform layer (mGCIPL) absolute thickness and intereye differences (IEDs) in 61 MS patients with a single unilateral ON 6–24 months prior (MS-ON) and 61 MS patients without ON history (MS-NON). Participants were 1:1 matched based on age, sex, and Expanded Disability Status Scale (EDSS). Statistical significance was evaluated using mixed linear effects models.

Results:

Age and ON status of eyes significantly interacted in predicting absolute pRNFL and mGCIPL thickness (pRNFL: B = 0.28, 95% confidence interval (CI) = [0.03, 0.54], p = 0.03, mGCIPL: B = 0.24, 95% CI = [0.01, 0.46], p = 0.04).

Discussion:

Our data show evidence of increased relapse-associated neuroaxonal damage in older individuals. This suggests less neuronal resilience following attacks and implicates age-effects beyond insidious progression. The growing population of older people with MS (pwMS) may be at increased risk of neurodegeneration and permanent disability worsening after acute MS attacks. Longitudinal studies should confirm this.

Keywords

Aging neurodegeneration disability multiple sclerosis optic neuritis axonal loss

Introduction

Multiple sclerosis (MS) is an autoimmune disease of the CNS characterised by demyelination and neuroaxonal degeneration.¹ While autoimmune-mediated demyelination associated with acute relapses is effectively targeted by disease-modifying therapies (DMTs), our understanding of the causes and timing of neurodegeneration in MS remains limited.^1,2 Disability progression in MS is age-dependent.³ While younger people with MS (pwMS) are more prone to experiencing relapse-associated worsening (RAW), aging is typically accompanied by a gradual shift towards progression independent of relapse (PIRA).⁴ There is additional evidence that older patients with MS may exhibit incomplete recovery from relapses.⁵ However, the direct contribution of age to neurodegeneration after an acute attack remains unclear. Whether age-associated disability in MS is solely the result of cumulative injury or whether disability depends on the timing of injury, is currently unknown.

The afferent visual system provides a valuable model in MS research by uniquely enabling pathway-specific longitudinal assessment of neurodegeneration, while also being a frequent site of inflammatory activity.^6,7 Damage to retinal ganglion cells and their axons resulting from optic neuritis (ON) is linked to a decrease in macular ganglion cell inner plexiform layer (mGCIPL) and peripapillary retinal nerve fibre layer (pRNFL) thickness as measured by optical coherence tomography (OCT).⁷ OCT uses coherent light to produce high-resolution, non-invasive, cross-sectional images of the retina. Besides a reduction in absolute pRNFL and mGCIPL thickness, intereye differences (IEDs) of pRNFL and mGCIPL can be used as indicators of previous ON attacks and to quantify the severity of ON-related damage.^8–10

Regardless of the overall rate of neural decline, it is unclear how individual attacks contribute to neurodegeneration over the course of MS. An increasing vulnerability of neurons towards stressors or a chronic increase in neuronal stressor load might promote greater attack-related neurodegeneration in aging pwMS.¹¹ This study aimed to understand, whether the neuroaxonal loss resulting from a single MS attack is age-dependent. We hypothesised that higher age at ON is associated with greater subsequent retinal atrophy.

Methods

Participants

This cross-sectional analysis comprised 122 pwMS: Sixty-one pwMS with previous ON (MS-ON) and 61 pwMS without ON history (MS-NON) (Figure 1). Participants of both groups were derived from large observational cohorts at UCSF (UCSF EPIC study,¹² n = 50) and Charité (CIS, VIMS and BerlImmun¹³ studies, n = 72), with standardised protocols for assessments and trained personnel. Inclusion criteria for the MS-ON group were 1) a diagnosis of relapsing-remitting MS (RRMS) according to 2017 criteria,¹⁴ 2) age > 12 years and 3) a single episode of unilateral ON within 6-24 months before OCT. ON episodes were defined based on clinical history. Participants of the MS-NON group were included by manual 1:1 matching for age at OCT, sex, and Expanded Disability Status Scale (EDSS)¹⁵ score. EDSS scoring was performed under the supervision of an experienced MS neurologist. Participants with multiple or bilateral ON events, relapse activity < 3 months before OCT and eyes with retinal comorbidities (e.g. diabetic or hypertensive retinopathy, posterior vitreous detachment) were excluded.

Figure 1.

Cohort selection flow chart.

The study was approved by the local ethics committees in San Francisco and Berlin (Charité Universitätsmedizin Berlin: EA1/362/20, EA1/163/12, EA1/182/10) and performed in line with the current version of the Declaration of Helsinki. All participants gave written informed consent before study inclusion.

OCT acquisition and measurements

Experienced personnel performed OCT acquisition, quality control, and image data reading at UCSF (DC, SCM) and Charité (CB, GK). Images of macular and peripapillary retinal layers were acquired with spectral domain OCT devices (Spectralis, Heidelberg Engineering, Heidelberg, Germany). Pupil dilation was not performed. Quality control was executed in compliance with OSCAR-IB criteria,¹⁶ and OCT data were reported according to the APOSTEL guidelines.¹⁷ Semi-automatic layer segmentation was performed using the OCT manufacturer’s software (eye explorer, Heidelberg Engineering; versions 1.10.12.0 (UCSF) and 1.10.2.0 (Charité)). pRNFL thickness was quantified from 12° ring scans around the optic nerve head with activated eye tracker. mGCIPL thickness was calculated as average thickness from macular volume scans by placing 3.45 mm diameter cylinders around the fovea.

Visual function

Visual function was measured for each eye separately. High-contrast visual acuity (HCVA) was acquired using retro-illuminated Early Treatment in Diabetic Retinopathy (ETDRS) charts at 4 m distance. Decimal acuity (Charité) or letter acuity (UCSF) scores were converted to the logarithm of the Minimum Angle of Resolution (logMAR). Low-contrast visual acuity (LCVA) was tested with Sloan 2.5% contrast charts, with the letter acuity (range 0–60) as the outcome of interest. Correction was either carried out using optimal optical (Charité) or habitual correction (UCSF).

Statistical analysis

All parameters are reported as mean ± SD unless specified otherwise. The relationship between age and pRNFL/mGCIPL thickness was examined using linear mixed-effects models including an interaction term between ON status of eyes and age at OCT, as well as subject-specific intercepts. Contralateral, ON-unaffected eyes of MS-ON participants (CON eyes) and eyes from MS-NON participants (NON eyes) were combined to increase power. Age at OCT was centred around the mean to allow interpretation of ON effects at mean age. Estimated marginal means were derived from these models to predict pRNFL/mGCIPL mean thickness differences between ON-affected and ON-unaffected eyes at ages 20, 30, 40, and 50 years. Effect size of interaction terms was further assessed by obtaining marginal ΔR² and standardised coefficients. Marginal ΔR² was calculated as the difference in marginal R² between models with and without the interaction term. The association between age and pRNFL/mGCIPL IEDs per MS-ON participant was analysed with linear regression models. IEDs were calculated as the absolute value of the difference in pRNFL/mGCIPL thickness between both eyes. Percentage IEDs were obtained by normalising absolute IEDs to the eye with greater pRNFL/mGCIPL thickness. Furthermore, the MS-ON group was split based on median age (35 years). The difference between pRNFL/mGCIPL IEDs in these two subgroups was analysed with the Mann–Whitney U-test to allow for non-parametric testing of non-normally distributed data. Normality of linear model residuals was assessed graphically with density curves. Given the collinearity of the structural measures and the limited number of comparisons, we did not apply formal multiple-comparison corrections. The impact of centre, sex, disease duration, DMT treatment (yes/no), steroid treatment (yes/no, only for IED comparisons within the MS-ON group) and time between ON and OCT (only for IED comparisons within the MS-ON group) as additional predictor variables was investigated by comparing hierarchically nested models. Goodness of fit was assessed through likelihood ratio testing, implemented via the anova function in R. A sensitivity analysis was performed after outlier exclusion. Outliers were defined as datapoints outside the interval [Q1 – 1.5*interquartile range (IQR); Q3 + 1.5*IQR]. Finally, age-dependency of ON-related changes in HCVA and LCVA was analysed analogously to pRNFL/mGCIPL absolute thickness modelling. Due to centre-specific acquisition protocols, visual acuity analyses were separated by centre and coefficients were standardised to create comparable estimates. Statistical significance was defined as p < 0.05. All analyses were performed with R (Version 4.2.1) including the packages lme4, emmeans, MuMIn and parameters. Figures were created using R and Adobe Illustrator. Statistical information is reported in accordance with the SAMPL guidelines.¹⁸

Results

Patient population

Sixty-one participants with RRMS were included in both the MS-ON and MS-NON group (UCSF: n = 25, Charité: n = 36) (Table 1). The two groups were matched for sex, age and EDSS at OCT (sex (male, N (%)): MS-ON: 19 (31.1), MS-NON: 19 (31.1), p > 0.99; age at OCT (years, mean ± SD): MS-ON: 34.2 ± 8.2, MS-NON: 34.9 ± 7.7, p = 0.62; EDSS (median (IQR)): MS-ON: 1.5 (1.0–1.5), MS-NON: 1.5 (1.0–2.0), p = 0.36). Fifty-two of the 61 MS-ON participants received steroid treatment. Therapy escalation to plasma exchange was confirmed for one participant; for seven participants, this could not be determined. pRNFL thickness was analysed in 239 eyes, 242 eyes were included in mGCIPL thickness analyses. For pRNFL IED analyses, 119 participants, and for mGCIPL IED analyses, 121 participants were included. Three MS-ON participants were excluded from pRNFL analyses because of missing ringscans (n = 2) or unilateral posterior vitreous detachment (n = 1; unaffected eye included in pRNFL absolute thickness modelling). One MS-NON participant was excluded from mGCIPL analyses because of missing macular scans.

Table 1.

Cohort description.

	MS-NON	MS-ON
Participants (N)	61	61
Sex (male, N (%))	19 (31.1)	19 (31.1)
Age at OCT (years, mean ± SD)	34.9 ± 7.7	34.2 ± 8.2
Age at ON (years, mean ± SD)	–	33.4 ± 8.2
EDSS (median (IQR))	1.5 (1.0–2.0)	1.5 (1.0–1.5)
Ethnicity (N (%))
White	52 (85.2)	51 (83.6)
Multiple/Other	9 (14.8)	10 (16.4)
Disease duration at OCT (years, mean ± SD)	3.9 ± 4.6	1.3 ± 2.3
Time between ON and OCT (months, mean ± SD)	–	11.0 ± 4.1
Treatment efficacy at OCT (N (%))¹⁹
high	8 (13.1)	11 (18.0)
moderate	16 (26.2)	14 (23.0)
modest	12 (19.7)	15 (24.6)
none	25 (41.0)	20 (32.8)
unknown	0 (0.0)	1 (1.6)

Participants of the MS-ON group experienced a single unilateral ON episode, MS-NON participants had no history of ON. Abbreviations: EDSS = Expanded Disability Status Scale, IQR = interquartile range, N = number, OCT = optical coherence tomography, ON = optic neuritis.

ON leads to retinal atrophy

pRNFL/mGCIPL thickness was significantly lower in ON eyes compared to CON and NON eyes (ON eyes: pRNFL (median (IQR)): 92.00 (81.65–99.50) µm, mGCIPL (median (IQR)): 77.02 (69.53–83.55) µm; CON eyes: pRNFL (median (IQR)): 103.00 (94.17–108.27) µm, mGCIPL (median (IQR)): 87.72 (83.44–89.90) µm; NON eyes: pRNFL (median (IQR)): 100.83 (94.73–105.00) µm, mGCIPL (median (IQR)): 88.11 (82.81–93.76) µm; p value for all comparisons of ON eyes with CON or NON eyes < 0.001) (Table 2). In addition, pRNFL/mGCIPL IEDs were increased in MS-ON compared to MS-NON participants (MS-ON: pRNFL IED (median (IQR)): 9.65 (5.46–13.51) µm, mGCIPL IED (median (IQR)): 10.70 (5.55–17.13) µm; MS-NON: pRNFL IED (median (IQR)): 2.02 (1.00–4.00) µm, mGCIPLIED (median (IQR)): 1.38 (1.07–2.92) µm; p value for both comparisons < 0.001) (Table 2).

Table 2.

Description of OCT-derived outcome parameters in ON-affected and -unaffected eyes/participants.

		Absolute thickness (µm)		IEDs (µm)
		pRNFL median (IQR)	mGCIPL median (IQR)	pRNFL median (IQR)	mGCIPL median (IQR)
MS-ON	ON eyes	92.00 (81.65–99.50)	77.02 (69.53–83.55)	9.65 (5.46–13.51)	10.70 (5.55–17.13)
MS-ON	CON eyes	103.00 (94.17–108.27)	87.72 (83.44–89.90)	9.65 (5.46–13.51)	10.70 (5.55–17.13)
MS-NON	NON eyes	100.83 (94.73–105.00)	88.11 (82.81–93.76)	2.02 (1.00–4.00)	1.38 (1.07–2.92)

Retinal atrophy after ON is associated with age

As expected, older age was associated with lower pRNFL/mGCIPL thickness in ON-affected and -unaffected eyes (ON eyes: pRNFL: B = −0.56, 95% confidence interval (CI) = [−0.97, −0.15], p = 0.01; mGCIPL: −0.42, 95% CI = [−0.80, −0.04], p = 0.03; CON + NON eyes: pRNFL: B = −0.25, 95% CI = [−0.44, −0.06], p = 0.01; mGCIPL: −0.21, 95% CI = [−0.35, −0.07], p = 0.003) (Figure 2). Linear mixed-effects models revealed a significant interacting effect of age at OCT and ON status on both pRNFL and mGCIPL thickness (pRNFL: B = 0.28, 95% CI = [0.03, 0.54], p = 0.03, mGCIPL: B = 0.24, 95% CI = [0.01, 0.46], p = 0.04; Table 3). Standardised coefficients indicated that, with every 1 SD increase in age, the difference in pRNFL/mGCIPL between ON-unaffected eyes and their ON-affected counterparts increased by an additional 0.2 SD/0.19 SD. However, as revealed by ΔR², the interaction term only accounted for 0.8% of the variance in the pRNFL model and 0.6% in the mGCIPL model.

Figure 2.

Age dependency of pRNFL/mGCIPL thickness in ON-affected (orange) and -unaffected eyes (blue). a: pRNFL absolute thickness model. Fifty-nine ON-affected eyes (UCSF: n = 23, Charité: n = 36) and 180 ON-unaffected eyes (UCSF: n = 72, Charité: n = 108) were included. b: mGCIPL absolute thickness model. Sixty-one ON-affected eyes (UCSF: n = 25, Charité: n = 36) and 181 ON-unaffected eyes (UCSF: 75, Charité: 106) were included. Linear regression lines were based on mixed linear effects models with uncentred variables, including age at OCT (Model formula: pRNFL/mGCIPL thickness ~ Age at OCT*ON status).

Table 3.

Linear mixed-effects modelling of pRNFL/mGCIPL thickness (µm) in ON-affected and -unaffected eyes, depending on age.

	pRNFL ON vs. NON/CON eyes			mGCIPL ON vs. NON/CON eyes
	Estimate	95% CI	p value	Estimate	95% CI	p value
Intercept	89.31	[87.12, 91.51]	<0.001	74.66	[72.78, 76.54]	<0.001
ON status^a	11.60	[9.51, 13.67]	<0.001	13.21	[11.37, 15.06]	<0.001
Age at OCT^b	–0.54	[−0.81, −0.28]	<0.001	−0.46	[−0.69, −0.23]	<0.001
Age:ON^a,b,c	0.28	[0.03, 0.54]	0.03	0.24	[0.01, 0.46]	0.04

Model formula: pRNFL/mGCIPL thickness ~ Age at OCT*ON-status. a: ON status = ON-unaffected. b: Age at OCT, mean-centred (years). c: interaction of age at OCT and ON status of eyes. Abbreviations: CI = confidence interval, mGCIPL = macular ganglion cell inner plexiform layer, OCT = optical coherence tomography, ON = optic neuritis, pRNFL = peripapillary retinal nerve fibre layer.

Estimated marginal mean differences in pRNFL were predicted to increase from −7.49 µm (95% CI = [−11.72, −3.27]) at 20 years to −15.97 µm (95% CI = [−20.56, −11.39]) at 50 years of age (Table 4). For mGCIPL, the predicted mean difference between ON-affected and -unaffected eyes increased from −9.74 µm (95% CI = [−13.5, −5.99]) at 20 years to −16.89 µm (95% CI = [−20.96, −12.82]) at 50 years of age (Table 5).

Table 4.

Changes in estimated marginal mean pRNFL thickness differences between ON-affected and -unaffected eyes from 20 to 50 years of age.

Age^a years	pRNFL ON-affected^bµm, [95% CI]	pRNFL ON-unaffected^c µm, [95% CI]	Δ pRNFL µm, [95% CI]	Change in Δ pRNFL/year^d µm, [95% CI]	p value
20	97.21 [92.78, 101.65]	104.71 [101.25, 108.17]	−7.49 [−11.72, −3.27]	0.28 [0.03, 0.54]	0.03
30	91.77 [89.28, 94.26]	102.09 [100.18, 104.00]	−10.32 [−12.69, −7.96]
40	86.33 [83.6, 89.05]	99.47 [97.47, 101.48]	−13.15 [−15.73, −10.57]
50	80.88 [76.05, 85.72]	96.86 [93.24, 100.48]	−15.97 [−20.56, −11.39]

Estimated marginal means were derived from linear mixed-effects modelling of pRNFL absolute thickness (model formula: pRNFL thickness ~ Age at OCT*ON status + (1|ID)). a: Age = Age at OCT. b: ON-affected eyes = ON-eyes. c: ON-unaffected eyes = NON- and CON-eyes. d: year = year of age at OCT; interaction term and p value taken from original model. Abbreviations: CI = confidence interval, OCT = optical coherence tomography, ON = optic neuritis, pRNFL = peripapillary retinal nerve fibre layer.

Table 5.

Changes in estimated marginal mean mGCIPL thickness differences between ON-affected and -unaffected eyes from 20 to 50 years of age.

Age^a years	mGCIPL ON-affected^b µm, [95% CI]	mGCIPL ON-unaffected^c µm, [95% CI]	Δ mGCIPL µm, [95% CI]	Change in Δ mGCIPL/year^d µm, [95% CI]	p value
20	81.29 [77.47, 85.11]	91.03 [88.16, 93.91]	–9.74 [–13.5, –5.99]	0.24 [0.01, 0.46]	0.04
30	76.73 [74.59, 78.87]	88.86 [87.27, 90.45]	–12.13 [–14.23, –10.02]
40	72.18 [69.85, 74.50]	86.69 [85.02, 88.35]	–14.51 [–16.80, –12.22]
50	67.62 [63.48, 71.75]	84.51 [81.51, 87.51]	–16.89 [–20.96, –12.82]

Estimated marginal means were derived from linear mixed-effects modelling of mGCIPL absolute thickness (model formula: mGCIPL thickness ~ Age at OCT*ON status + (1|ID). a: Age = age at OCT. b: ON-affected eyes = ON-eyes. c: ON-unaffected eyes = NON- and CON-eyes. d: year = year of age at OCT; interaction term and p value taken from original model. Abbreviations: CI = confidence interval, mGCIPL = macular ganglion cell inner plexiform layer, OCT = optical coherence tomography, ON = optic neuritis.

Consistently, there was a numerical increase in pRNFL and mGCIPL IEDs of MS-ON participants with greater age, not meeting the predefined threshold of statistical significance (pRNFL: B = 0.28, 95% CI = [−0.09, 0.65], p = 0.13; mGCIPL: B = 0.26, 95% CI = [−0.07, 0.59], p = 0.12; Figures 3a and b). Only a small part of the variance of the dataset was explained by our models (multiple r² = 0.04 for both pRNFL and mGCIPL IEDs). No statistical difference in pRNFL (p = 0.71) or mGCIPL IEDs was observed (p = 0.26) when comparing MS-ON participants younger and older than the age median (35 years; Figure 3c). Participants from both age groups did not significantly differ in EDSS or disease duration (EDSS (median (IQR)): older participants: 1.5 (1.0-1.5), younger participants 1.5 (1.0-1.5), p = 0.67; disease duration at OCT (years; mean ± SD): older participants: 1.5 ± 2.4, younger participants: 1.1 ± 2.2, p = 0.52). Analyses of percentage IEDs yielded similar results (data not shown).

Figure 3.

Age dependency of pRNFL/mGCIPL IEDs in MS-ON and MS-NON participants. Participants of the MS-ON group (yellow) experienced a single unilateral ON episode, MS-NON participants (green) had no history of ON. a, b: Linear regression of age at OCT and pRNFL/mGCIPL IEDs. Model formula: pRNFL/mGCIPL IED ~ Age at OCT. pRNFL data from 119 participants (MS-ON: n = 58 (UCSF: n = 22, Charité: n = 36); MS-NON: n = 61 (UCSF: n = 25, Charité: n = 36)) was included. mGCIPL data stemmed from 121 participants (MS-ON: n = 61 (UCSF: n = 25, Charité: n = 36); MS-NON: n = 60 (UCSF: n = 25, Charité: n = 35)). c: Distribution of pRNFL/mGCIPL IEDs in older compared to younger MS-ON participants after median split. Median age: 35 years. Older group: n = 28/29, younger group: n = 30/32 (pRNFL/mGCIPL). The box depicts the IQR, the horizontal line indicates the median. Whiskers extend to Q1 – 1.5 × IQR and Q3 + 1.5 × IQR. Statistical testing was performed using the Mann–Whitney U-test.

Age (in IED models) or the interaction of age and ON status (in absolute thickness modelling) presented the best predictors of OCT-derived outcome parameters in all statistical models, when adding additional predictor variables (Supplemental Tables 1 and 2). In pRNFL/mGCIPL thickness models, age and ON status remained the only independent variables with statistical significance. Model estimates remained relatively stable (pRNFL: Age:ON interaction term = 0.28–0.30, mGCIPL: Age:ON interaction term = 0.23–0.24; pRNFL absolute IED: Age at OCT = 0.28–0.34, mGCIPL absolute IED: Age at OCT = 0.23–0.26). In addition, the studied outcome parameters were not significantly influenced by differing treatment efficacies (data not shown).

After outlier exclusion, the significant interaction of age and ON status persisted in modulating pRNFL thickness (Age: B = −0.52, 95% CI = [−0.77, −0.27], p < 0.001; ON status: B = 10.67, 95% CI = [8.72, 12.58], p < 0.001; Age:ON interaction: B = 0.25, 95% CI = [0.02, 0.48], p = 0.04). In mGCIPL models, age and ON status remained significant predictors of absolute thickness, but the interaction of these variables did not seem to majorly impact this parameter (Age: B = −0.34, 95% CI = [−0.55, −0.13], p = 0.002; ON status: B = 11.67, 95% CI = [10.09, 13.27], p < 0.001; Age:ON interaction: B = 0.13, 95% CI = [−0.07, 0.32], p = 0.21) (Supplemental Table 3).

Visual recovery after ON depends on age

HCVA data were available for 41 ON eyes (UCSF: n = 13, Charité: n = 28) and for 136 NON/CON eyes (UCSF: n = 51, Charité: n = 85). For LCVA, data were available for 41 ON eyes (UCSF: n = 13, Charité: n = 28) and for 138 NON/CON eyes (UCSF: n = 53, Charité: n = 85). Interestingly, no consistent effect of age and ON status on ETDRS logMAR was observed 6–24 months after ON. Prior ON was associated with higher ETDRS logMAR in UCSF participants, while age and the interaction of age and ON status did not seem to affect HCVA outcomes. Contrastingly, there was a trend towards an age-dependent increase of ETDRS logMAR in Charité participants, as well as a numerical interaction of age and ON status, which did not meet the predefined threshold of statistical significance (Charité: Age: B = 0.35, 95% CI = [0.00, 0.71], p = 0.05; ON status: B = -0.29, 95% CI = [−0.64, 0.06], p = 0.11; Age:ON interaction: B = −0.29, 95% CI [−0.65, 0.07], p = 0.12; UCSF: Age: B = 0.06, 95% CI = [−0.48, 0.60], p = 0.83; ON status: B = −0.60, 95% CI = [−1.11, −0.08], p = 0.02; Age:ON interaction: B = 0.01, 95% CI [–0.52, 0.54], p = 0.96) (Tables 6 and 7). In LCVA testing, age and prior ON were significantly associated with decreased letter counts in participants from both centres. Moreover, there was a significant interaction between age and ON status affecting the letter counts of Charité participants. Modelling of UCSF participant data showed a numerically similar effect, albeit not statistically significant (Charité: Age: B = −0.45, 95% CI = [–0.80, –0.10], p = 0.01; ON status: B = 0.82, 95% CI = [0.47, 1.18], p < 0.001; Age:ON interaction: B = 0.48, 95% CI = [0.11, 0.85], p = 0.01; UCSF: Age: B = -0.71, 95% CI = [−1.15, −0.28], p = 0.002; ON status: B = 0.86, 95% CI = [0.44, 1.28], p < 0.001; Age:ON interaction: B = 0.34, 95% CI = [−0.09, 0.77], p = 0.12) (Table 8).

Table 6.

Description of visual acuity parameters in ON-affected and -unaffected eyes by centre.

		MS-ON participants		MS-NON
	Centre	ON eyes	CON eyes	NON eyes
ETDRS (logMAR)	UCSF	0.02 (–0.10–0.02)	–0.10 (–0.12–0.06)	–0.07 (–0.16–0.03)
ETDRS (logMAR)	Charité	–0.20 (–0.20–0.10)	–0.20 (–0.23–0.10)	–0.20 (–0.20–0.10)
LCVA (Letter count)	UCSF	22.00 (19.00–31.00)	29.00 (24.00–34.00)	29.50 (26.75–35.00)
LCVA (Letter count)	Charité	44.50 (40.00–48.00)	49.00 (45.75–52.00)	49.00 (45.00–53.00)

All data presented as median (IQR) due to non-normality. Participants of the MS-ON group experienced a single unilateral ON episode, MS-NON participants had no history of ON. ON eyes: ON-affected eyes of MS-ON participants (UCSF: n = 13, Charité: n = 28), CON eyes: contralateral, ON-unaffected eyes of MS-ON participants (UCSF: n = 13, Charité: n = 28). NON eyes: eyes of MS-NON participants (UCSF: n = 38/40 (ETDRS/LCVA), Charité: n = 57). Abbreviations: ETDRS = Early Treatment of Diabetic Retinopathy Study, IQR = interquartile range, LCVA = low-contrast visual acuity, logMAR = logarithmic minimal angle of resolution, ON = optic neuritis, UCSF = University of California San Francisco.

Table 7.

Linear mixed-effects modelling of HCVA (ETDRS, logMAR) in ON-affected and -unaffected eyes, depending on age.

ETDRS	ON vs. NON/CON eyes Charité(logMAR)			ON vs. NON/CON eyes UCSF(logMAR)
	Estimate	95% CI	p value	Estimate	95% CI	p value
Intercept	0.22	[–0.13, 0.57]	0.21	0.48	[–0.05, 1.00]	0.07
ON status^a	−0.29	[–0.64, 0.06]	0.11	–0.60	[–1.11, –0.08]	0.02
Age at OCT^b	0.35	[0.00, 0.71]	0.05	0.06	[–0.48, 0.60]	0.83
Age:ON^a,^b,^c	−0.29	[–0.65, 0.07]	0.12	0.01	[–0.52, 0.54]	0.96

Model formula: ETDRS ~ Age at OCT*ON status; coefficients are standardised. a: ON status = ON-unaffected. b: Age at OCT, mean-centred (years). c: interaction of age at OCT and ON status of eyes. Abbreviations: CI = confidence interval, ETDRS = early treatment in diabetic retinopathy study, HCVA = high-contrast visual acuity, logMAR = logarithmic minimum angle of resolution, OCT = optical coherence tomography, ON = optic neuritis, UCSF = University of California San Francisco.

Table 8.

Linear mixed-effects modelling of LCVA (letter count) in ON-affected and -unaffected eyes, depending on age.

LCVA	ON vs. NON/CON eyes Charité(letter count)			ON vs. NON/CON eyes UCSF(letter count)
	Estimate	95% CI	p value	Estimate	95% CI	p value
Intercept	–0.63	[–0.97, –0.29]	0.003	–0.71	[–1.15, –0.28]	0.002
ON status^a	0.82	[0.47, 1.18]	<0.001	0.86	[0.44, 1.28]	<0.001
Age at OCT^b	–0.45	[–0.80, –0.10]	0.01	–0.72	[–1.16, –0.28]	0.002
Age:ON^a,^b,^c	0.48	[0.11, 0.85]	0.01	0.34	[–0.09, 0.77]	0.12

Model formula: LCVA ~ age at OCT*ON status; coefficients are standardised. a: ON status = ON-unaffected. b: Age at OCT, mean-centred (years). c: interaction of age at OCT and ON status of eyes. Abbreviations: CI = confidence interval, LCVA = low-contrast visual acuity, OCT = optical coherence tomography, ON = optic neuritis, UCSF = University of California San Francisco.

Discussion

While neurodegeneration is recognised as the cause of chronic disability in MS, drivers of age-related disability accrual are only poorly understood.³ We aimed to analyse the relationship between age and ON-related retinal atrophy. Previous studies have suggested greater overall neurodegeneration rates in younger pwMS,^20–22 but age differences in the contribution of individual attacks to the neurodegenerative load were, to our knowledge, previously unknown. A better understanding of age-related differences in vulnerability towards relapse-associated damage could help inform ongoing discussions regarding discontinuation of treatment in older pwMS.³⁵

Our data suggest that retinal neuroaxonal loss during acute events is age-dependent. This was reflected in the interaction of age and ON status of eyes in linear mixed-effects models of pRNFL and mGCIPL thickness, adding to the known age-dependent neuroaxonal loss in MS, which is also observed in the retina.^23,24 Older age was also associated with a tendency towards increased pRNFL/mGCIPL IEDs after ON. Moreover, the interaction between age and ON status also influenced LCVA 6–24 months after ON. This suggests that older age at ON is associated with impaired functional recovery. Multiple other studies have previously demonstrated an association of age with declining functional recovery from MS relapses.^5,23–26 However, the biological basis of this has not been determined. Our work suggests that early neuroaxonal loss is a likely significant contributor to this failure to recover functionally.

A previous study reported no predictive effect of age on HCVA 12 months after a first episode of ON in pwMS, while pRNFL thickness at 12 months was significantly associated with age at baseline in univariate regression analysis.²⁷ Due to the lack of a control group of pwMS without ON, this study could not separate the attack-specific contribution of age from the results of natural and MS-related aging.²⁷

Major strengths of our study include the close 1:1 matching of an MS-NON group to control for age- and attack-unrelated neurodegeneration, as well as a focus on patients after their first clinical ON attack. Since ON events often occur early in the disease, most participants were DMT-naïve at the ON event, mitigating possible effects of differing treatment efficacies. Older and younger MS-ON participants did not significantly differ with respect to EDSS and disease duration. This argues against the assumption that ON activity in older participants could be indicative of a more severe MS course.

Nevertheless, limitations have to be considered: The lack of baseline OCT data before ON events and the relatively narrow age span present limiting factors of our study. Owing to the cross-sectional design, it cannot be excluded that retinal asymmetry or thinning of retinal layers may be due to factors unrelated to age and retinal inflammation. As ON frequently presents the heralding event in MS, obtaining baseline data remains a challenge. However, the fact that IEDs remained low in MS-NON participants throughout the studied age span supports the underlying assumption that the increase in IEDs observed in MS-ON participants can be substantially attributed to the prior ON episode. ON-unaffected eyes of MS-ON participants can thus serve as baseline surrogates for the affected eyes. Considering that only 26% of participants were older than 40 years, the reported effects might be stronger when studied in cohorts with larger age ranges, although recruitment of older participants is limited by lower attack frequencies^28,29 and the requirement of no prior ON attack. The average difference in pRNFL/mGCIPL thickness of 0.24–0.28 µm per year of age at ON estimated by our models however amounts to an average difference of ~2–3 µm per decade of age. This difference could already be considered clinically meaningful: previously published cross-sectional data estimated an average reduction in pRNFL thickness of 2.2–2.9 µm to be associated with a 1-line change in HCVA testing.^30,31 In addition, our study suggests impairment of LCVA recovery in older patients with ON, further underlining the clinical implications of our structural findings.

The outlier-dependency of the mGCIPL data and the marginal statistical significance in IED analyses could be explained by the relatively small sample size. Generally, our variables of interest only explained small parts of the overall variance of our outcome data. Differences in attack severity might explain some of this additional data variability. Data on visual acuity at nadir was unfortunately not available to be included as a measure of attack severity. Given that older pwMS typically show decreased adaptive immune activity, it is plausible that relapse-related inflammation is less severe in this group.³² Differences in the time between ON and OCT could also have added to the data variability. Furthermore, overall disability of both MS-ON and MS-NON participants was relatively low. While this might be due to the short disease duration associated with first-time ON attacks, we cannot exclude that our study population might not be reflective of all disease courses.

Larger, prospectively gathered datasets from multiple centres and in more diverse populations are needed to replicate and expand upon these findings. The inclusion of additional measures of neurodegeneration is crucial, since age-related neurodegeneration has been suggested to be region- and assessment-dependent.²⁰ Attack definitions may however be more challenging to standardise when studying other areas of the CNS, for example, with magnetic resonance imaging (MRI).

Preclinical studies should further elucidate possible mechanisms behind age-related changes in neuroaxonal loss, such as impaired remyelination and changes in neuronal vulnerability. The degree of demyelination of the optic nerve has been shown to correspond with retinal neurodegeneration and remyelinating treatment has been associated with increased neuroaxonal stability.^33,34 Since impaired remyelination is a known pathophysiological feature associated with aging in MS, age-related remyelination failure and increased neuroaxonal loss might be linked.³² In addition, neuronal senescence might reduce the intrinsic capacity for self-protection. Unravelling the contribution of these complex mechanisms to age-associated neuronal loss could provide essential insights for future therapeutic strategies to prevent neurodegeneration and disability progression in MS.

Despite the relatively young age of our participants, our data suggest an association between age at ON and subsequent retinal neurodegeneration. In the clinical setting, an enhanced understanding of age-effects in MS attacks could aid both as a prognostic factor and help guide therapeutic decision-making. Although less frequent, individual MS attacks might have more detrimental effects in older patients. This conclusion is also relevant in the context of emerging discussions regarding discontinuation of treatment in aging pwMS.³⁵

Supplemental Material

sj-docx-1-msj-10.1177_13524585251374507 – Supplemental material for Investigating the impact of age on neurodegeneration after optic neuritis in multiple sclerosis using optical coherence tomography

Supplemental material, sj-docx-1-msj-10.1177_13524585251374507 for Investigating the impact of age on neurodegeneration after optic neuritis in multiple sclerosis using optical coherence tomography by Henriette KM Reinsberg, Frederike C Oertel, Daniel Casillas, Gelsomina Kaufhold, Shivany Condor Montes, Charlotte Bereuter, Makenna C Chapman, Pia S Sperber, Bruce A Cree, Judith Bellmann-Strobl, Stephen L Hauser, Tanja Schmitz-Hübsch, Christian Cordano, Hanna G Zimmermann, Ahmed Abdelhak, Friedemann Paul and Ari J Green in Multiple Sclerosis Journal

Footnotes

Acknowledgements

University of California, San Francisco MS-EPIC Team: Refujia Gomez, Meagan Harms, Tiffany Cooper, Anna Sindalovsky, Stacy Callier, Harkee Halait, Adam Santaniello, Adam Renschen, Simone Sacco MD, Nico Papinutto PhD, Ariele Greenfield MD, Joanne Guo MD, Sasha Gupta MD, Joseph Sabatino MD PhD, Richard Cuneo MD, Jill Hollenbach PhD, Jeffrey Gelfand MD MAS, Riley Bove MD MSc, Samuel Pleasure MD PhD, Michael Wilson MD, Roland G Henry PhD, Sergio Baranzini PhD, Jorge R Oksenberg PhD, Stephen L. Hauser MD.

HKMR thanks the DAAD PROMOS program and the German Academic Scholarship Foundation for fellowship support

Declaration of Conflicting Interests

The authors declared the following potential conflicts of interest with respect to the research, authorship, and/or publication of this article: F.C. Oertel serves as a board member of IMSVISUAL, received research grants by the Hertie foundation, Novartis, Deutsche Forschungsgemeinschaft and the Deutsche Gesellschaft für Neurologie, received speaker honoraria by UCB and fellowship support by NMSS and AAN.; B.A.C. Cree has received personal compensation for consulting from Alexion, Autobahn, Avotres, Biogen, Boston Pharma, EMD Serono, Hexal/Sandoz, Immunic AG, Kyverna, Neuron23, Novartis, Sanofi, and TG Therapeutics and received research support from Genentech and Kyverna.; J.Bellmann-Strobl received institutional research support from NEMOS e.V., Alexion and Bayer AG, personal compensation from Alexion, speaking honoraria and travel grants from Bayer Healthcare, Horizon/Amgen, Novartis and Sanofi-Aventis/Genzyme, in addition received compensation for serving on a scientific advisory board of Alexion, Roche and Merck, all unrelated to the presented work.; S.L. Hauser currently serves on the scientific advisory board of Accure, Alector, Annexon, Hinge Therapeutics. He has previously consulted for BD, Moderna, NGM Bio, Nurix Therapeutics, Pheno Therapeutics and previously served on the Board of Directors of Neurona. Dr. Hauser also has received travel reimbursement and writing support from F. Hoffmann-La Roche and Novartis AG for anti-CD20-therapy-related meetings and presentations. He has received grants by NIH/NINDS (R35NS111644), NMSS (SI-2001-35701) and Valhalla Foundation.; T. Schmitz-Hübsch received speaking honoraria from AbbVie, Bayer, Roche and research grants from bms/Celgene, Roche. C. Cordano reports funding from the Italian Multiple Sclerosis Foundation (FISM) Rita-Levi Montalcini Grant code 2019/BC/002.; H.G. Zimmermann reports grants from Novartis, unrelated to this study.; A.Abdelhak received research grants from the Department of Defence, the German Multiple Sclerosis Society, and the Weill Institute for Neurosciences, he received an advisory fee from OctaveBio and speaking honoraria from Roche, all unrelated to this study.; F.Paul received research support, grants and personal fees by Deutsche Forschungsgemeinschaft, Bundesministerium für Bildung und Forschung, Einstein Foundation, Guthy Jackson Charitable Foundation, EU FP7 Framework Programme, Biogen, Sanofi Genzyme, Merck Serono, Novartis, Bayer, Roche, Alexion, Horizon, Neuraxpharm, Almirall, GlaxoSmithKline, Hexal, Viela Bio, UCB, Mitsubishi Tanabe and Celgene. He is an advisory board member of Celgene, Roche, UCB and Merck, academic editor at PLos One and associate editor of Neurology^® Neuroimmunology & Neuroinflammation.; A. J. Green reports research support, grants and personal fees from NINDS (R01NS105741, R01AG062562, R01AG038791), NMSS (RG-1707-28564), SBIR, All May See, Rachleff Family, Westridge Foundation, JAMA Neurology, Anonymous philanthropy, Roche, Pipeline Pharmaceuticals, the National MS Society, UCSF Innovation Programme for Remyelination and Repair: IP/Patents Bazodoxefine and other miscellaneous outside the submitted work. In addition, Dr. Green has a patent Small Molecule drug for Remyelination pending and has worked on testing off label compounds for remyelination. The remaining authors have no conflicts to disclose.

Funding

The authors received no financial support for the research, authorship, and/or publication of this article.

Ethical Considerations

Consent to Participate

All participants gave written informed consent before study inclusion.

Consent for Publication

Not applicable.

Data Availability

Anonymised data used for this study are available from the corresponding author upon reasonable request.

ORCID iDs

Henriette KM Reinsberg

Frederike C Oertel

Shivany Condor Montes

Bruce A Cree

Stephen L Hauser

Christian Cordano

Ahmed Abdelhak

Ari J Green

Supplemental Material

Supplemental material for this article is available online.

References

Thompson

Baranzini

Geurts

, et al. Multiple sclerosis. Lancet 2018; 391: 1622–1636.

Vaughn

Jakimovski

Kavak

, et al. Epidemiology and treatment of multiple sclerosis in elderly populations. Nat Rev Neurol 2019; 15(6): 329–342.

Zeydan

Kantarci

. Impact of age on multiple sclerosis disease activity and progression. Curr Neurol Neurosci 2020; 20: 24.

Müller

Cagol

Lorscheider

, et al. Harmonizing definitions for progression independent of relapse activity in multiple sclerosis. JAMA Neurol 2023; 80: 1232–1245.

Conway

Zeydan

Uygunoglu

, et al. Age is a critical determinant in recovery from multiple sclerosis relapses. Mult Scler J 2018; 25: 1754–1763.

Balcer

. Optic neuritis. New Engl J Medicine 2006; 354: 1273–1280.

Graves

Oertel

Walt

, et al. Leveraging visual outcome measures to advance therapy development in neuroimmunologic disorders. Neurology Neuroimmunol Neuroinflammation 2021; 9: e1126.

Coric

Balk

Uitdehaag

BMJ

, et al. Diagnostic accuracy of optical coherence tomography inter-eye percentage difference for optic neuritis in multiple sclerosis. Eur J Neurol 2017; 24(12): 1479–1484.

Petzold

Fraser

Abegg

, et al. Diagnosis and classification of optic neuritis. Lancet Neurology 2022; 21: 1120–1134.

10.

Oertel

Zimmermann

Motamedi

, et al. Diagnostic value of intereye difference metrics for optic neuritis in aquaporin-4 antibody seropositive neuromyelitis optica spectrum disorders. J Neurol Neurosurg Psychiatry 2023; 94(7): 560–566.

11.

Sanai

Saini

Benedict

, et al. Aging and multiple sclerosis. Mult Scler J 2015; 22: 717–725.

12.

University of California San Francisco MS-EPIC Team, Cree

BAC

Hollenbach

, et al. Silent progression in disease activity–free relapsing multiple sclerosis. Ann Neurol 2019; 85(5): 653–666.

13.

Sperber

Brandt

Zimmermann

, et al. Berlin registry of neuroimmunological entities (BERLimmun): Protocol of a prospective observational study. BMC Neurol 2022; 22: 479.

14.

Thompson

Banwell

Barkhof

, et al. Diagnosis of multiple sclerosis: 2017 revisions of the McDonald criteria. Lancet Neurology 2018; 17: 162–173.

15.

Kurtzke

. Rating neurologic impairment in multiple sclerosis An Expanded Disability Status Scale (EDSS). Neurology 1983; 33(11): 1444–1452.

16.

Schippling

Balk

Costello

, et al. Quality control for retinal OCT in multiple sclerosis: Validation of the OSCAR-IB criteria. Mult Scler J 2014; 21: 163–170.

17.

Aytulun

Cruz- Herranz

Aktas

, et al APOSTEL 2.0 recommendations for reporting quantitative optical coherence tomography studies. Neurology 2021; 97: 68–79.

18.

Lang

Altman

. Basic statistical reporting for articles published in Biomedical journals: the ‘Statistical analyses and methods in the published literature’ or the SAMPL guidelines. Int J Nurs Stud 2015; 52: 5–9.

19.

Hauser

Cree

BAC

. Treatment of multiple sclerosis: A review. Am J Med 2020; 133: 1380–1390e2.

20.

Cordano

Nourbakhsh

Yiu

, et al. Differences in age-related retinal and cortical atrophy rates in multiple sclerosis. Neurology 2022; 99: e1685–e1693.

21.

Pulido-Valdeolivas

Andorrà

Gómez -Andrés

, et al. Retinal and brain damage during multiple sclerosis course: Inflammatory activity is a key factor in the first 5 years. Sci Rep 2020; 10: 13333.

22.

Balk

Cruz-Herranz

Albrecht

, et al. Timing of retinal neuronal and axonal loss in MS: A longitudinal OCT study. J Neurol 2016; 263(7): 1323–1331.

23.

Zanghì

Galgani

Bellantonio

, et al. Relapse-associated worsening in a real-life multiple sclerosis cohort: The role of age and pyramidal phenotype. Eur J Neurol 2023; 30(9): 2736–2744.

24.

Sotiropoulos

Lokhande

Healy

, et al. Relapse recovery in multiple sclerosis: effect of treatment and contribution to long-term disability. Mult Scler J Exp Transl Clin 2021; 7(2).

25.

Hosny

Shehata

Ahmed

, et al. Predictors of severity and outcome of multiple sclerosis relapses. BMC Neurol 2023; 23: 67.

26.

Cossburn

Ingram

Hirst

, et al. Age at onset as a determinant of presenting phenotype and initial relapse recovery in multiple sclerosis. Mult Scler J 2011; 18: 45–54.

27.

Deschamps

Shor

Vignal

, et al. Prospective longitudinal study on prognostic factors of visual recovery and structural change after a first episode of optic neuritis. Eur J Neurol 2022; 29(9): 2781–2791.

28.

Tremlett

Zhao

Joseph

, et al. Relapses in multiple sclerosis are age- and time-dependent. J Neurol Neurosurg Psychiatry 2008; 79(12): 1368–1374.

29.

Jakimovski

Eckert

Zivadinov

, et al. Considering patient age when treating multiple sclerosis across the adult lifespan. Expert Rev Neurother 2021; 21(3): 353–364.

30.

Walter

Ishikawa

Galetta

, et al. Ganglion cell loss in relation to visual disability in multiple sclerosis. Ophthalmology 2012; 119(6): 1250–1257.

31.

Fisher

Jacobs

Markowitz

, et al. Relation of visual function to retinal nerve fiber layer thickness in multiple sclerosis. Ophthalmology 2006; 113(2): 324–332.

32.

Graves

Krysko

Hua

, et al. Ageing and multiple sclerosis. Lancet Neurology 2023; 22(1): 66-77.

33.

Oertel

Krämer

Motamedi

, et al. Visually evoked potential as prognostic biomarker for neuroaxonal damage in multiple sclerosis from a multicenter longitudinal cohort. Neurol Neuroimmunol Neuroinflamm 2023; 10(3).

34.

Abdelhak

Cordano

Boscardin

, et al. Plasma neurofilament light chain levels suggest neuroaxonal stability following therapeutic remyelination in people with multiple sclerosis. J Neurol Neurosurg Psychiatry 2022; 93: 972–977.

35.

Corboy

Fox

Kister

, et al. Risk of new disease activity in patients with multiple sclerosis who continue or discontinue disease-modifying therapies (DISCOMS): A multicentre, randomised, single-blind, phase 4, non-inferiority trial. Lancet Neurol 2023; 22(7): 568–577.

Supplementary Material

Please find the following supplemental material available below.

For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.

For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.

0.00 MB

0.04 MB