Corneal axonal loss as an imaging biomarker of neurodegeneration in multiple sclerosis: a longitudinal study

Abstract

Background:

Resourceful endpoints of axonal loss are needed to predict the course of multiple sclerosis (MS). Corneal confocal microscopy (CCM) can detect axonal loss in patients with clinically isolated syndrome and established MS, which relates to neurological disability.

Objective:

To assess corneal axonal loss over time in relation to retinal atrophy, and neurological and radiological abnormalities in MS.

Methods:

Patients with relapsing-remitting (RRMS) (n = 68) or secondary progressive MS (SPMS) (n = 15) underwent CCM and optical coherence tomography. Corneal nerve fibre density (CNFD-fibres/mm²), corneal nerve branch density (CNBD-branches/mm²), corneal nerve fibre length (CNFL-mm/mm²) and retinal nerve fibre layer (RNFL-μm) thickness were quantified along with neurological and radiological assessments at baseline and after 2 years of follow-up. Age-matched, healthy controls (n = 20) were also assessed.

Results:

In patients with RRMS compared with controls at baseline, CNFD (p = 0.004) and RNFL thickness (p < 0.001) were lower, and CNBD (p = 0.003) was higher. In patients with SPMS compared with controls, CNFD (p < 0.001), CNFL (p = 0.04) and RNFL thickness (p < 0.001) were lower. For identifying RRMS, CNBD had the highest area under the receiver operating characteristic (AUROC) curve (0.99); and for SPMS, CNFD had the highest AUROC (0.95). At follow-up, there was a further significant decrease in CNFD (p = 0.04), CNBD (p = 0.001), CNFL (p = 0.008) and RNFL (p = 0.002) in RRMS; in CNFD (p = 0.04) and CNBD (p = 0.002) in SPMS; and in CNBD (p = 0.01) in SPMS compared with RRMS. Follow-up corneal nerve loss was greater in patients with new enhancing lesions and optic neuritis history.

Conclusion:

Progressive corneal and retinal axonal loss was identified in patients with MS, especially those with more active disease. CCM may serve as an imaging biomarker of axonal loss in MS.

Keywords

axonal loss biomarker corneal confocal microscopy multiple sclerosis neurodegeneration

Introduction

The assessment of clinical disability and neuroimaging constitute the mainstay for assessing disease progression and therapeutic decision making in multiple sclerosis (MS). While the burden of inflammation on magnetic resonance imaging (MRI) can predict clinical outcomes to some extent,¹ quantitative evaluation of lesion metrics has failed to account fully for disease progression.² Moreover, in more advanced disease, neurodegeneration may occur in the absence of neuroinflammation.³ Indeed, neuronal and axonal loss are major pathological features of irreversible and progressive MS⁴ accompanied by progressive brain atrophy and ventricular enlargement.⁵ However, the accuracy of automated brain volumetric assessment has been questioned,⁶ with the demonstration of hypertrophy of the putamen⁷ and an increase in brain volume in 26.3% of longitudinal scans.⁸

There has been considerable interest in evaluating surrogate markers of neurodegeneration to monitor and predict disease progression, and as endpoints in trials of neuroprotection in MS.⁹ Optical coherence tomography (OCT) has been used to show retinal nerve fibre layer (RNFL) thinning in MS and has been related to the severity of axonal damage in the visual pathways¹⁰ and brain atrophy.¹¹ Corneal confocal microscopy (CCM) is a rapid, non-invasive ophthalmic imaging technique, which enables in vivo assessment of the corneal subbasal nerves.¹² CCM has been used to demonstrate axonal loss in a range of peripheral neuropathies, including diabetic neuropathy,¹³ Friedreich’s ataxia,¹⁴ chronic inflammatory demyelinating polyneuropathy¹⁵ and chemotherapy-induced neuropathy.¹⁶ Longitudinal studies have also established that CCM can predict incident diabetic neuropathy,¹⁷ and clinical trials^18–20 have shown that CCM identifies early nerve regeneration which precedes an improvement in symptoms, neurophysiology and intra-epidermal nerve density.

We and others have also demonstrated corneal nerve loss in patients with MS,^21–26 Parkinson’s disease²⁷ and dementia.²⁸ In our recent study of a cohort of patients with relapsing-remitting MS (RRMS), followed over 2 years, we showed that progressive corneal nerve loss was associated with increasing disability.²⁹ We aimed to compare corneal and retinal axonal loss longitudinally in patients with MS in relation to neurological and radiological abnormalities.

Methods

Study participants and setting

The present study constitutes a pooled analysis of two prospective, observational studies conducted in Doha, Qatar (September 2016–March 2020), and Konya, Turkey (May 2016–July 2020). This research followed the tenets of the Declaration of Helsinki and obtained approvals from the institutional review boards of Weill Cornell Medicine–Qatar (1500064 with approval date 2 May 2016), Hamad Medical Corporation (15218/15 with approval date 5 September 2016) and the research ethics committee of Necmettin Erbakan University, Meram Faculty of Medicine (2016/570 with approval date 27 May 2016; and 2018/1134 with approval date 5 January 2018). Participants or their legally authorized representatives gave informed written consent to participate in this study. Patients with RRMS (n = 68) and secondary progressive MS (SPMS) (n = 15) underwent assessment of neurological disability, CCM and OCT at baseline and follow-up using a standardized protocol. Age-matched, healthy controls (n = 20) were also assessed at baseline. Four patients with MS (RRMS = 3, SPMS = 1) were lost to follow-up. Reporting of results in this study followed the STARD guidelines³⁰ and a participant flow chart is presented in Figure 1. Inclusion criteria were diagnosis of MS based on the 2010 McDonald criteria³¹ and age 18–75 years. Exclusion criteria were regular contact lens use (>3 times/week), presence of chronic ophthalmic disease, active optic neuritis (ON), or history of refractive or cataract surgery. Patients with metabolic, rheumatologic or neurologic co-morbidities associated with peripheral neuropathy were excluded based on glycated haemoglobin A1c, presence of anti-nuclear antibody, serum B₁₂/folate, immunoglobulins and medical history.

Figure 1.

Participant flow chart.

Clinical assessments

The expanded disability status scale (EDSS)³² was used to rate neurological impairment. The EDSS is a physician-administered composite ranging from 0 (normal) to 10 (death due to MS) in 0.5 increments (from EDSS > 1). Scores from 0 to 4 evaluate general neurological function; 4 to 6 focuses on walking ability; and scores > 6 indicate loss of independence. The MS severity score (MSSS) was calculated from the EDSS and MS duration.³³ The annualized relapse rate (ARR) was calculated as the ratio of the total number of relapses divided by the total number of days in the study for the group multiplied by 365. Clinical assessments preceded ophthalmic assessments. Brain MRI (Magnetom Skyra 3T, Siemens Medical Systems, Erlangen, Germany) was performed for all follow-up patients within 1 month of clinical and ophthalmic assessments and T1-weighted images were captured at axial, coronal and three-dimensional (3D) MP-RAGE sagittal planes. Enhancing lesions at follow-up were detected on post-contrast, T1-weighted images with fat saturation.

CCM

CCM scans (Heidelberg Retinal Tomograph III Rostock Cornea Module, Heidelberg Engineering GmbH, Heidelberg, Germany) were performed as per our established methodology^13,34 in the same room, under constant light conditions by four trained examiners. Scans were performed on the right eye of healthy controls and patients without ON, on the ON eye in patients with history of unilateral ON, or on the right eye in patients with history of bilateral ON. Longitudinal studies have established that RNFL loss is greater in ON versus non-ON eyes.^35,36 The same has not yet been established for CCM. Based on central location and clarity (i.e. visibility of subbasal nerves without pressure lines), three non-overlapping images/patient from the central subbasal nerve plexus were manually analysed using CCMetrics (Courtesy of Prof. RA Malik, Weill Cornell Medicine–Qatar of Cornell University, Doha, Qatar) as per our established protocol³⁷ by the same examiner. The specific parameters measured were corneal nerve fibre density (CNFD) (fibres/mm²), corneal nerve branch density (CNBD) (branches/mm²) and corneal nerve fibre length (CNFL) (mm/mm²), and the results were calculated as an average of analysed images per participant. Examiners were masked to the subtype of MS, clinical and MRI examination results during the ophthalmic examination and image analysis phase. Figure 2 illustrates CCM images from a healthy control (2a, b) and patients with MS (2c, d) at baseline.

Figure 2.

Representative CCM images from controls and patients with multiple sclerosis. Panel shows an image from the corneal subbasal nerve plexus of a healthy control (a). Same image analysed with CCMetrics (b). Corneal subbasal nerve plexus images from a patient with RRMS (c) and SPMS (d) showing corneal axonal loss compared with the healthy control.

OCT

Peripapillary RNFL thickness measurements were performed with a spectral-domain OCT (Spectralis OCT, Heidelberg Engineering GmbH, Heidelberg, Germany) on the same eye scanned by CCM, under the same light conditions without pupil dilation. RNFL measurements were performed by using circular scans with a scanning angle of 12°, which equates to a retinal diameter of 3.5 mm when assuming a standard corneal curvature of 7.7 mm and with the eye tracker activated to minimize motion artefacts. All RNFL scans in this study were performed in high-speed mode, and a signal strength of ⩾ 20/30 was set as the minimum acceptable quality along with optimal centering of the circular scan on the optic nerve head. Follow-up RNFL scans were performed with reference to their respective baseline scans using the built-in ‘progression’ function.

Statistical analysis

Statistical analysis was performed using Prism (version 9.1.0 for Mac, GraphPad Software Inc, San Diego, CA, US). Data were confirmed to follow a normal distribution by means of a Shapiro–Wilk test (p > 0.05) and by visual inspection of the respective Q-Q plots. Based on the variability seen in our pilot MS data^21,23,26 and previously published studies,^13,34,38 we estimated (paired t-test) that a sample of 80 participants will give 90% power to detect an effect size of at least 0.365 in CCM measures with a type 1 error of 0.05. One-way analysis of variance (post hoc Tukey) was used for between-group comparisons. Fisher’s exact test was used to calculate differences in categorical data. A paired t-test was used for within-group comparisons and an unpaired t-test was used for the subgroup analysis. Receiver operating characteristic (ROC) curve analysis was performed to assess the discriminating ability of corneal nerve parameters and RNFL based on the area under the ROC (AUROC) curve, and sensitivity and specificity identified using Youden’s index. Multiple linear regression modelling was used to assess the influence of MS-specific features on the change (Δ) in neurological disability (as expressed by MSSS, EDSS), corneal and retinal nerve parameters calculated as the difference between follow-up and baseline. ΔMSSS, ΔEDSS, ΔCNFD, ΔCNBD, ΔCNFL and ΔRNFL were used as the dependent variables and presence of new enhancing lesions, positive ON history, use of disease-modifying treatments and ARR as the independent variables. Two models were used: an unadjusted (model 1); and an adjusted one for age and sex (model 2). As a second step, we used multiple linear regression modelling to assess the influence of baseline and follow-up clinical parameters as independent variables (time since diagnosis, EDSS, MSSS) on ΔCNFD, ΔCNBD, ΔCNFL and ΔRNFL, adjusting for subject’s age and sex. Continuous data are expressed as [mean difference (standard error of mean difference), 95% confidence interval, p value]. The reported p values are two-sided and a p ⩽ 0.05 was considered significant.

Protocol and data availability statement

The study protocol and all anonymized, individual-level data used in this manuscript are available to qualified researchers by direct request to the corresponding author.

Results

Baseline clinical and demographic results

In healthy controls compared with RRMS and SPMS, there was no significant difference in age or sex. Compared with RRMS, patients with SPMS had a higher EDSS [2.4 (0.5), 1.4–3.4, p < 0.001], MSSS [2.5 (0.7), 1.1–4.0, p < 0.001] and longer time since diagnosis [2.6 (1.1), 0.5–4.7, p = 0.02). There was no significant difference in age [4.1 (2.4), –0.7 to 8.9, p = 0.09], sex [44 (65%) versus 10 (67%), p > 0.99], ON prevalence [37 (54) versus 8 (53)], the total number of relapses [0.2 (0.7), –1.2 to 1.6, p = 0.76] and ARR [0.7 (0.4), –0.04 to 1.5, p = 0.06] (Table 1).

Table 1.

Baseline clinical and demographic characteristics.

	Controls	RRMS	SPMS	p value
n (%)	20	68 (82)	15 (18)	–
Age (years)	34.1 (1.6)	35.5 (1.0)	39.6 (2.3)	0.13
Sex (% women)	11 (55)	44 (65)	10 (67)	0.77
Time since diagnosis (years)	–	6.9 (0.5)	9.5 (1.0)	0.02
ON history, n (%)	–	37 (54)	8 (53)	0.99
Relapses	–	2.8 (0.3)	3.0 (0.5)	0.8
EDSS	–	1.7 (0.2)	4.1 (0.7)	<0.001
MSSS	–	2.8 (0.3)	5.3 (0.2)	<0.001
ARR	–	1.5 (0.2)	2.2 (0.4)	0.06
DMT use, n (%)	–	55 (81)	15 (100)	–
Beta interferon, n (%)	–	19 (29)	2 (13.4)	–
Fingolimod, n (%)	–	15 (22)	6 (40)	–
Dimethyl fumarate, n (%)	–	6 (9)	6 (40)	–
Other,^a n (%)	–	15 (22)	1 (6.7)	–
CNFD (fibres/mm²)	41.3 (1.3)	34.3 (1.2)^A	30.3 (1.6)^A	<0.001
CNBD (branches/mm²)	101.0 (8.0)	151.0 (7.6)^A	122.0 (15.4)	0.002
CNFL (mm/mm²)	26.1 (0.7)	25.1 (0.8)	21.6 (1.5)^A	0.06
RNFL (μm)	98.9 (1.7)	85.0 (1.7)^A	75.4 (4.6)^A,B	<0.001

ARR, annualized relapse rate; CNBD, corneal nerve branch density; CNFD, corneal nerve fibre density; CNFL, corneal nerve fibre length; DMT, disease-modifying therapy; EDSS, expanded disability status scale; MSSS, multiple sclerosis severity score; ON, optic neuritis; RNFL, retinal nerve fibre layer; RRMS, relapsing-remitting multiple sclerosis; SPMS, secondary progressive multiple sclerosis.

Data are expressed as mean (standard error of mean) or as n (%) where applicable. For multiple comparisons, p value corresponds to the overall p value of the test statistic (ANOVA, Fisher’s). Superscripted letters indicate: ^Asignificantly different from controls; ^Bsignificantly different from RRMS. CNFD: significantly different in RRMS versus controls (p = 0.004); and in SPMS versus controls (p < 0.001). CNBD: significantly different in RRMS versus controls (p = 0.002). CNFL significantly different in SPMS versus controls (p = 0.04). RNFL significantly different in RRMS versus controls (p < 0.001); in SPMS versus controls (p < 0.001) and versus RRMS (p = 0.05).

Teriflunomide (n = 6), Azathioprine (n = 6), Glatiramer Acetate (n = 2), Natalizumab (n = 1) and Alemtuzumab (n = 1).

Baseline CCM and OCT

The detailed results are presented in Table 1 and Figure 3(a)–(d).

Figure 3.

Graphic illustration of baseline CCM and OCT data and ROC curves. Upper half: graphs indicate mean ± range (1%–99%) for baseline CNFD (a), CNFL (b), CNBD (c) and RNFL (d) in controls (left) compared with RRMS (middle) and SPMS (right). ROC curves illustrating the capacity of CCM and OCT measures to discriminate between controls and patients with RRMS (e); controls and patients with SPMS (f); and between MS subtypes (g). The diagonal line indicates where sensitivity = 1 – specificity.

RRMS

CNFD [7.1 (2.2), 1.9 to 12.2, p = 0.004] and RNFL [14.1 (3.5), 5.8 to 22.4, p < 0.001] were significantly lower, and CNBD [−49.7 (14.5), –84.3 to −15.1, p = 0.003] was significantly higher compared with healthy controls.

SPMS

CNFD [11.1 (3.0), 4.0 to 18.1, p < 0.001], CNFL [4.5 (2.0), 0.04 to 9.0, p = 0.04] and RNFL [23.5 (4.7), 12.3 to 34.8, p < 0.001] were lower compared with healthy controls. Compared with RRMS, RNFL [9.6 (4.0), 0.2 to 19.0, p = 0.05] was lower.

ROC curve analysis

The detailed results are presented in Table 2 and Figure 3(e)–(g).

Table 2.

Receiver operating characteristic curve analysis results.

	Cut-off point	AUROC	95% CI	p value	Sensitivity	Specificity	LR
Controls versus RRMS
CNFD	<39.6	0.77	0.67–0.87	<0.001	0.78	0.70	2.6
CNBD	<53.1	0.99	0.98–1.00	<0.001	0.96	0.95	19.1
CNFL	<27.4	0.80	0.70–0.89	<0.001	0.72	0.70	2.4
RNFL	<94.5	0.82	0.73–0.91	<0.001	0.79	0.74	3.0
Controls versus SPMS
CNFD	<36.7	0.94	0.87–1.00	<0.001	0.87	0.75	3.5
CNBD	>102.6	0.64	0.44–0.84	0.16	0.67	0.55	1.5
CNFL	<24.2	0.82	0.66–0.99	0.001	0.80	0.70	2.7
RNFL	<92.5	0.91	0.77–1.0	<0.001	0.87	0.79	4.1
RRMS versus SPMS
CNFD	<29.9	0.63	0.50–0.76	0.12	0.60	0.65	1.7
CNBD	>45.8	0.91	0.78–1.00	<0.001	0.87	0.87	6.5
CNFL	<25.8	0.90	0.83–0.97	<0.001	0.87	0.82	4.9
RNFL	<81.5	0.70	0.54–0.85	0.02	0.67	0.62	1.7

AUROC, area under the receiver operating characteristic curve; CI, confidence interval; CNBD, corneal nerve branch density; CNFD, corneal nerve fibre density; CNFL, corneal nerve fibre length; LR, likelihood ratio; RNFL, retinal nerve fibre layer; RRMS, relapsing-remitting multiple sclerosis; SPMS, secondary progressive multiple sclerosis.

Healthy controls versus RRMS

CNBD [AUROC = 0.99, sensitivity/specificity = 0.96/0.95, 95% confidence interval = 0.98–1.0, likelihood ratio (LR) = 19.1] showed the highest discriminative performance followed by RNFL (0.82, 0.79/0.74, 0.73–0.91, 3.0).

Healthy controls versus SPMS

CNFD (0.94, 0.87/0.75, 0.87–1.0, 3.5) showed the highest discriminative performance followed by RNFL (0.91, 0.87/0.79, 0.77–1.0, 4.1).

RRMS versus SPMS

CNBD (0.91, 0.87/0.87, 0.78–1.0, 6.6) showed the highest discriminative performance followed by CNFL (0.9, 0.87/0.82, 0.83–0.97, 4.9).

Follow-up clinical and demographic results

The average follow-up time was [20.7 (4.1), 19.7 to 21.7] months. New gadolinium-enhancing lesions were clinically confirmed in n = 19 (23%) of patients. Although there was no significant change in the number of patients on DMT, n = 27 (33%) patients switched or discontinued DMT by their follow-up visit. At follow-up, differences between RRMS and SPMS in EDSS, MSSS, ARR and relapses were comparable with baseline. Two additional patients (3%) in the RRMS group had a clinically confirmed episode of ON since their baseline visit, which had resolved by the time of their follow-up visit. In the RRMS group, there was a significant increase in relapses [0.2 (0.06), 0.1 to 0.3, p < 0.001]; ARR [0.1 (0.03), 0.06 to 0.2, p < 0.001] and MSSS [−0.3 (0.1), –0.6 to −0.05, p = 0.02) at follow-up compared with baseline. In the SPMS group, there was a significant increase in relapses [0.4 (0.1), 0.1–0.7, p = 0.009] and ARR [0.2 (0.1), 0.04–0.5, p = 0.02].

Follow-up CCM and OCT

RRMS

There was a decrease in CNFD [−2.3 (1.2), –4.6 to −0.01, p = 0.04], CNBD [−28.5 (8.2), –44.9 to −12.0, p = 0.001], CNFL [−1.8 (0.7), –3.2 to −0.5, p = 0.008] and RNFL [−1.2 (0.4), –2.0 to −0.5, p = 0.002] at follow-up compared with baseline; see Figure 4(a)–(d).

Figure 4.

Trajectories of MS subtypes at follow-up. Graphs illustrating the mean difference between baseline and follow-up visits for patients with RRMS (dashed line) and SPMS (continuous line) for CNFD (a), CNFL (b), CNBD (c) and RNFL (d).

SPMS

There was a decrease in CNFD [−3.1 (2.6), –8.7 to 2.5, p = 0.04] and CNBD [−47.4 (11.8), –73.0 to −21.9, p = 0.002] at follow-up compared with baseline. In SPMS compared with RRMS, there was a significant difference in CNBD [−47.9 (18.3), –84.4 to −11.5, p = 0.01], see Figure 4(a)–(d).

Subgroup analysis

Clinical parameters, CCM and RNFL in relation to new enhancing lesions at follow-up

In patients with new enhancing lesions, there was a decrease in CNFD [−2.8 (1.4), to −5.7 to 0.08, p = 0.05], CNBD [−44.0 (14.3), –73.9 to −14.1, p = 0.006], CNFL [−2.9 (0.7), –4.4 to −1.5, p < 0.001] and RNFL [−2.2 (0.8), –3.9 to −0.5, p = 0.01] with a higher number of relapses [0.5 (0.1), 0.2 to 0.7, p < 0.001], EDSS [0.5 (0.1), 0.2 to 0.8, p < 0.001] and ARR [0.3 (0.06), 0.1–0.4, p = 0.001] at follow-up compared with baseline. In patients without new enhancing lesions, there was a decrease in CNBD [−23.2 (8.2), –39.6 to −6.7, p = 0.007] and RNFL [−0.9 (0.4), –1.6 to −0.2, p = 0.02] with a higher number of relapses [0.2 (0.05), 0.1 to 0.3, p < 0.001] and ARR [0.1 (0.03), 0.05 to 0.2, p = 0.001] and lower MSSS [−0.6 (0.2), –0.9 to −0.3, p < 0.001].

Clinical parameters, CCM and RNFL in relation to ON

In patients with a history of ON, there was a decrease in CNBD [−32.0 (9.5), to −51.0 to −12.9, p = 0.002], CNFL [−1.9 (0.8), −3.6 to −0.3, p = 0.02] and RNFL [−1.0 (0.5), −2.0 to −0.09, p = 0.03] with a higher number of relapses [0.3 (0.08), 0.2 to 0.5, p < 0.001] and ARR [0.2 (0.04), 0.1 to 0.3, p < 0.001]. In patients without ON history, there was a decrease in RNFL [−1.3 (0.4), –2.1 to −0.4, p = 0.005] with a higher number of relapses [0.2 (0.06), 0.04 to 0.3, p = 0.01] and ARR [0.1 (0.04), 0.02 to 0.2, p = 0.01]. In patients with left eye ON versus non-left eye ON (i.e. right eye or bilateral ON), there was no significant difference at baseline in CNFD [4.8 (2.5), −0.3 to 10.0, p = 0.07], CNBD [−9.9 (19.6), –49.4 to 29.6, p = 0.61], CNFL [1.5 (1.9), –2.3 to 5.4, p = 0.4] and RNFL [2.6 (4.5), –6.4 to 11.6, p = 0.56]. Similarly, there was no significant difference at follow-up in CNFD [4.7 (2.8), –1.0 to 10.3, p = 0.1], CNBD [27.7 (19.0), –10.6 to 65.9, p = 0.15], CNFL [2.8 (1.8), –0.9 to 6.4, p = 0.13] and RNFL [0.3 (4.5), –8.8 to 9.5, p = 0.94]. There was no difference in age [1.4 (2.7), –4.0 to 6.7, p = 0.61] and time since diagnosis [−0.9 (1.2), –3.4 to 1.6, p = 0.45].

Clinical parameters, CCM and RNFL in relation to disease-modifying therapy status

In patients with MS who switched or discontinued DMT at follow-up, there was a decrease in CNBD [−23.8 (8.8), –41.9 to −5.7, p = 0.01] with a higher number of relapses [0.3 (0.08), 0.1 to 0.5, p = 0.003] and ARR [0.1 (0.05), 0.04 to 0.2, p = 0.007]. In patients with MS with no change in DMT status, there was a decrease in CNBD [−27.0 (10.7), –48.4 to −5.7, p = 0.01], RNFL [−1.4 (0.3), –2.1 to −0.8, p < 0.001] and MSSS [−0.4 (0.2), –0.7 to −0.08, p = 0.01] with a higher number of relapses [0.3 (0.07), 0.1 to 0.4, p < 0.001] and ARR [0.2 (0.04), 0.08 to 0.2, p < 0.001].

Multiple linear regression analysis

The detailed results of the multiple linear regression modelling are presented in Tables 3 and 4. A model with ΔMSSS as the dependent variable and new enhancing lesions, ON history, DMT use and ARR as the independent variables predicted 8.8% of the variability in ΔMSSS in the unadjusted model and 6.5% in the adjusted model. Presence of new enhancing lesions was associated with an unadjusted mean difference of [0.959 (0.305), p = 0.002] and adjusted mean difference of [0.96 (0.309), p = 0.003] (Table 3). A model with ΔEDSS as the dependent variable predicted 8.1% of the variability in ΔEDSS in the unadjusted model and 5.9% in the adjusted model. Presence of new enhancing lesions was associated with an unadjusted mean difference of [0.637 (0.204), p = 0.003] and adjusted mean difference of [0.631 (0.207), p = 0.003] (Table 3). A model with ΔCNBD as the dependent variable predicted 8.3% of the variability in ΔCNBD in the unadjusted model and 8.9% in the adjusted model. Higher ARR was associated with an unadjusted mean difference of [14.95 (4.84), p = 0.003] and adjusted mean difference of [14.35 (4.85), p = 0.004] (Table 3). For ΔCNFL, the unadjusted model predicted 4.4% of the variability and higher ARR was associated with an unadjusted mean difference of [0.84 (0.42), p = 0.04] (Table 3). A regression model with ΔCNBD as the dependent variable and MSSS at baseline and follow-up as the independent variable adjusted for age and sex predicted 7.3% of the variability in ΔCNBD [7.05 (2.74), p = 0.01] at baseline and 4.7% [5.89 (2.8), p = 0.04] at follow-up (Table 4).

Table 3.

Multiple linear regression models: unadjusted (model 1) and adjusted for age and sex (model 2).

Dependent variable	Independent variables	Model 1			Model 2
Dependent variable	Independent variables	R² adjusted	Unadjusted mean difference (SEM)	p value	R² adjusted	Adjusted mean difference (SEM)	p value
ΔMSSS	New Gd lesions	0.09	0.96 (0.31)	0.002	0.07	0.96 (0.31)	0.003
	ON history		−0.18 (0.26)	0.48		−0.19 (0.26)	0.47
	DMT use		0.21 (0.38)	0.59		0.21 (0.39)	0.60
	ARR		−0.05 (0.09)	0.55		−0.05 (0.09)	0.56
ΔEDSS	New Gd lesions	0.08	0.64 (0.2)	0.003	0.06	0.63 (0.21)	0.003
	ON history		−0.05 (0.17)	0.78		−0.05 (0.18)	0.77
	DMT use		0.12 (0.26)	0.64		0.11 (0.26)	0.66
	ARR		0.02 (0.06)	0.72		0.02 (0.06)	0.75
ΔCNFD	New Gd lesions	−0.03	−0.68 (2.53)	0.80	−0.03	−0.78 (2.54)	0.76
	ON history		−0.85 (2.14)	0.69		–1.14 (2.15)	0.60
	DMT use		–1.624 (3.179)	0.61		–1.85 (3.18)	0.56
	ARR		0.87 (0.73)	0.24		0.82 (0.73)	0.27
ΔCNBD	New Gd lesions	0.08	–14.12 (16.81)	0.40	0.09	–15.31 (16.79)	0.37
	ON history		–15.81 (14.21)	0.27		–17.49 (14.23)	0.22
	DMT use		–9.37 (21.11)	0.66		11.20 (21.06)	0.60
	ARR		14.95 (4.84)	0.003		14.35 (4.85)	0.004
ΔCNFL	New Gd lesions	0.04	–1.37 (1.46)	0.35	0.03	–1.40 (1.47)	0.35
	ON history		–2.04 (1.23)	0.10		–2.14 (1.25)	0.09
	DMT use		−0.75 (1.83)	0.68		−0.82 (1.85)	0.66
	ARR		0.84 (0.42)	0.04		0.83 (0.43)	0.06
ΔRNFL	New Gd lesions	0.002	–1.35 (0.76)	0.08	0.02	–1.34 (0.75)	0.08
	ON history		0.41 (0.64)	0.62		0.54 (0.64)	0.40
	DMT use		−0.47 (0.95)	0.52		−0.40 (0.94)	0.67
	ARR		−0.02 (0.22)	0.94		−0.09 (0.22)	0.97

ARR, annualized relapse rate; CNBD, corneal nerve branch density; CNFD, corneal nerve fibre density; CNFL, corneal nerve fibre length; DMT, disease-modifying therapy; EDSS, expanded disability status scale; Gd, gadolinium; MSSS, multiple sclerosis severity score; ON, optic neuritis; RNFL, retinal nerve fibre layer; SEM, standard error of the mean.

Table 4.

Multiple linear regression models adjusted for age and sex.

Dependent variable	Independent variables	Baseline			Follow-up
Dependent variable	Independent variables	R² adjusted	Adjusted mean difference (SEM)	p value	R2 adjusted	Adjusted mean difference (SEM)	p value
ΔCNFD	Time since diagnosis	−0.01	0.10 (0.27)	0.72	−0.01	0.13 (0.28)	0.64
	EDSS	−0.01	−0.14 (0.53)	0.79	−0.01	−0.78 (1.38)	0.56
	MSSS	−0.01	−0.01 (0.41)	0.97	−0.01	−0.08 (0.41)	0.84
ΔCNBD	Time since diagnosis	0.001	–1.32 (1.89)	0.49	−0.003	−0.74 (1.94)	0.70
	EDSS	0.03	6.31 (3.65)	0.09	−0.003	–4.03 (9.37)	0.67
	MSSS	0.07	7.05 (2.74)	0.01	0.05	5.88 (2.83)	0.04
ΔCNFL	Time since diagnosis	−0.03	−0.01 (0.16)	0.93	−0.03	−0.02 (0.17)	0.91
	EDSS	−0.03	0.05 (0.32)	0.87	−0.001	–1.26 (0.80)	0.12
	MSSS	−0.02	0.29 (0.24)	0.24	−0.03	0.15 (0.25)	0.56
ΔRNFL	Time since diagnosis	0.008	0.03 (0.08)	0.75	0.01	0.02 (0.09)	0.84
	EDSS	0.02	−0.13 (0.17)	0.43	0.01	−0.25 (0.40)	0.54
	MSSS	0.01	−0.06 (0.12)	0.63	0.02	−0.12 (0.12)	0.32

CNBD, corneal nerve branch density; CNFD, corneal nerve fibre density; CNFL, corneal nerve fibre length; EDSS, expanded disability status scale; MSSS, multiple sclerosis severity score; RNFL, retinal nerve fibre layer; SEM, standard error of the mean.

Discussion

We demonstrate corneal and retinal axonal loss in patients with RRMS and SPMS. Corneal and retinal axonal measures show excellent capacity to discriminate healthy controls from patients with RRMS and SPMS and patients with RRMS from SPMS. Furthermore, both CCM and OCT show progressive corneal and retinal axonal loss. Subgroup analysis showed greater corneal and retinal axonal loss in patients with new gadolinium-enhancing lesions and history of ON. Furthermore, regression modelling revealed a significant association between the change in CNBD and ARR at follow-up, and change in CNBD was also moderately associated with neurological disability.

Axonal degeneration is the major driver of progressive disability in MS, and a key finding from a pathological study was that smaller axons are more susceptible to injury.³⁹ CCM is a non-invasive imaging modality for rapid quantification of small fibres in diabetic¹³ and other neuropathies,^14,15 and predicts incident diabetic neuropathy.¹⁷ We and others have also shown significant corneal axonal loss in RRMS,^21–26 which was associated with neurological disability. In the present study, we show significant progressive corneal nerve loss and RNFL thinning in both RRMS and SPMS and higher CNBD in patients with RRMS. CNBD at follow-up differed significantly not only between patients with MS and controls but also between RRMS and SPMS suggesting cumulative axonal loss with more advanced MS course. We have previously shown higher CNBD in patients with Parkinson’s disease, which was related to neurological disability.²⁷ Clinical trials^18,40 of DMTs in diabetic neuropathy have demonstrated that early nerve regeneration is characterized by an increase in branch density. Studies have shown that acute axonal damage is partially reversible in the initial stages of MS,^4,41 and in this context increased CNBD may indicate nerve regeneration.

The ROC curve analysis revealed that baseline corneal nerve and RNFL thinning, a well-established measure of axonal degeneration in MS, have comparable diagnostic utility for both RRMS and SPMS. Furthermore, CNBD has excellent sensitivity and specificity for RRMS, while CNFD and CNFL are comparable with RNFL in discriminating patients with SPMS from healthy controls. Overall, our results are in keeping with previous findings on the diagnostic capacity of CCM for diabetic neuropathy. Ahmed et al.⁴² reported a sensitivity/specificity of 0.85/0.84 for a diagnosis of diabetic neuropathy with a positive LR of 5.3. Petropoulos et al.¹³ reported that reduced CNFD was associated with a sensitivity/specificity of 0.79/0.78 and a positive LR of 4.6 for diabetic neuropathy. Chen et al.³⁸ showed a 0.82/0.71 sensitivity/specificity for diabetic neuropathy. In the present study, CNBD had the highest discriminating capacity for RRMS, indicating a more dynamic role for branches, particularly during the early, highly active stage of MS.

RNFL thinning is progressive over time and predicts disability worsening.^43,44 However, attenuation of RNFL thinning has been described with longer disease duration in MS.⁴⁵ In a population with diabetes, there is a continuing decline in CCM measures over 8 years which is greater in patients who subsequently develop neuropathy.¹⁷ Bitirgen et al.²⁹ recently reported a progressive reduction in corneal nerve fibre area and width over 2 years in a small cohort of patients with RRMS. In the present study, we demonstrate a significant reduction in CNFD, CNBD and CNFL at follow-up in both RRMS and SPMS.

Subgroup analysis showed that both corneal and retinal axonal loss were greater in patients with clinically confirmed new gadolinium-enhancing lesions at follow-up, and that patients with a history of ON had greater loss of corneal nerve fibres and branches at follow-up compared with patients without ON. However, regression modelling revealed that ARR and MSSS were only associated with ΔCNBD at follow-up, and the presence of new enhancing lesions was associated with ΔMSSS and ΔEDSS. A potential explanation for these findings is that corneal axonal alterations are more susceptible to cumulative inflammation reflected by a higher relapse frequency over time as opposed to acute inflammation at a single timepoint. Previous studies have reported accelerated RNFL thinning after ON⁴⁶ paralleled by trans-synaptic changes in the central nervous system. No association has been previously reported between ON history and corneal nerve density.^21–25

Previous studies have reported an association between CCM measures and neurological disability^21–24,29 with the strength of the relationship increasing for greater EDSS.²⁵ In the present study, ΔCNBD was associated with the number of relapses and disability at baseline and follow-up, and ΔCNFL was associated with ΔEDSS and ΔMSSS. It has been postulated that neurodegeneration in MS may initially be detectable only on highly sensitive techniques, before it becomes diffuse and irreversible with progressive disability.⁴⁷ In this context, our findings imply a varying axonal burden as detected by CCM on patients with MS, particularly during the earlier stage, which increases in relation to relapses and worsening of disability.

We acknowledge several limitations in our study. The small sample size limits the generalizability of our findings. In addition, the relatively short follow-up time does not allow interpretation of our findings in relation to meaningful worsening of neurological disability and brain atrophy. We also acknowledge a lack of follow-up in the controls. However, it has been previously established that corneal nerves of healthy persons remain stable over a 3-year period.⁴⁸ Finally, our ophthalmic imaging protocol targeting the ON eye rather than the non-ON eye may have skewed our results towards greater alterations at both timepoints. Nolan-Kenney et al.⁴⁶ have described inter-eye difference thresholds in relation to the ON status, including subclinical ON. In the present study, we did not assess for subclinical ON in patients with MS, and the extent to which it could have influenced our results is unknown. The magnitude of differences between ON and non-ON eyes and in relation to subclinical ON should be established.

In conclusion, this is the largest study to date showing progressive corneal axonal loss in relation to radiological and clinical severity in patients with MS. CCM may add to the imaging toolbox as a surrogate marker of axonal degeneration and regeneration in MS. Larger studies with longer follow-up are required to assess the relationship between corneal nerve loss, ganglion cell inner plexiform thinning and disease progression.

Footnotes

Acknowledgements

The authors would like to thank Dr Zehra Akpinar from the Department of Neurology, Meram Faculty of Medicine, Necmettin Erbakan University, for her contribution of study participants.

Declarations

ORCID iDs

Saadat Kamran

Rayaz A. Malik

References

Brex

Ciccarelli

O’Riordan

, et al. A longitudinal study of abnormalities on MRI and disability from multiple sclerosis. N Engl J Med 2002; 346: 158–164.

Barkhof

Calabresi

Miller

, et al. Imaging outcomes for neuroprotection and repair in multiple sclerosis trials. Nat Rev Neurol 2009; 5: 256–266.

Pérez-Cerdá

Sánchez-Gómez

Matute

The link of inflammation and neurodegeneration in progressive multiple sclerosis. Mult Scler Demyelinating Disord 2016; 1: 1–8.

Ferguson

Matyszak

Esiri

, et al. Axonal damage in acute multiple sclerosis lesions. Brain 1997; 120: 393–399.

Zhuo

Duan

, et al. Subtyping relapsing–remitting multiple sclerosis using structural MRI. J Neurol 2021; 268: 1808–1817.

de Sitter

Verhoeven

Burggraaff

, et al. Reduced accuracy of MRI deep grey matter segmentation in multiple sclerosis: an evaluation of four automated methods against manual reference segmentations in a multi-center cohort. J Neurol 2020; 267: 3541–3554.

Morelli

Baldini

Sartori

, et al. Early putamen hypertrophy and ongoing hippocampus atrophy predict cognitive performance in the first ten years of relapsing-remitting multiple sclerosis. Neurol Sci 2020; 41: 2893–2904.

Uher

Bergsland

Krasensky

, et al. Interpretation of brain volume increase in multiple sclerosis. J Neuroimaging 2021; 31: 401–407.

Tur

Moccia

Barkhof

, et al. Assessing treatment outcomes in multiple sclerosis trials and in the clinical setting. Nat Rev Neurol 2018; 14: 75–93.

10.

Kolbe

Marriott

van der Walt

, et al. Diffusion tensor imaging correlates of visual impairment in multiple sclerosis and chronic optic neuritis. Invest Ophthalmol Vis Sci 2012; 53: 825–832.

11.

Saidha

Al-Louzi

Ratchford

, et al. Optical coherence tomography reflects brain atrophy in multiple sclerosis: a four-year study. Ann Neurol 2015; 78: 801–813.

12.

Müller

Marfurt

Kruse

, et al. Corneal nerves: structure, contents and function. Exp Eye Res 2003; 76: 521–542.

13.

Petropoulos

Alam

Fadavi

, et al. Rapid automated diagnosis of diabetic peripheral neuropathy with in vivo corneal confocal microscopyautomated detection of diabetic neuropathy. Invest Ophthalmol Vis Sci 2014; 55: 2071–2078.

14.

Pagovich

Zhao

, et al. Corneal confocal microscopy: Neurologic disease biomarker in Friedreich ataxia. Ann Neurol 2018; 84: 893–904.

15.

Stettner

Hinrichs

Guthoff

, et al. Corneal confocal microscopy in chronic inflammatory demyelinating polyneuropathy. Ann Clin Transl Neurol 2015; 3: 88–100.

16.

Ferdousi

Azmi

Petropoulos

, et al. Corneal confocal microscopy detects small fibre neuropathy in patients with upper gastrointestinal cancer and nerve regeneration in chemotherapy induced peripheral neuropathy. PLoS ONE 2015; 10: e0139394.

17.

Lewis

EJH

Lovblom

Ferdousi

, et al. Rapid corneal nerve fiber loss: a marker of diabetic neuropathy onset and progression. Diabetes Care 2020; 43: 1829–1835.

18.

Azmi

Jeziorska

Ferdousi

, et al. Early nerve fibre regeneration in individuals with type 1 diabetes after simultaneous pancreas and kidney transplantation. Diabetologia 2019; 62: 1478–1487.

19.

Culver

Dahan

Bajorunas

, et al. Cibinetide improves corneal nerve fiber abundance in patients with sarcoidosis-associated small nerve fiber loss and neuropathic pain. Invest Ophthalmol Vis Sci 2017; 58: BIO52–BIO60.

20.

Lewis

Perkins

Lovblom

, et al. Effect of omega-3 supplementation on neuropathy in type 1 diabetes: A 12-month pilot trial. Neurology 2017; 88: 2294–2301.

21.

Bitirgen

Akpinar

Malik

, et al. Use of corneal confocal microscopy to detect corneal nerve loss and increased dendritic cells in patients with multiple sclerosis. JAMA Ophthalmol 2017; 135: 777–782.

22.

Mikolajczak

Zimmermann

Kheirkhah

, et al. Patients with multiple sclerosis demonstrate reduced subbasal corneal nerve fibre density. Mult Scler 2017; 23: 1847–1853.

23.

Petropoulos

Kamran

, et al. Corneal confocal microscopy: an imaging endpoint for axonal degeneration in multiple sclerosis. Invest Ophthalmol Vis Sci 2017; 58: 3677–3681.

24.

Testa

De Santis

Scotto

, et al. Neuroaxonal degeneration in patients with multiple sclerosis: an optical coherence tomography and in vivo corneal confocal microscopy study. Cornea 2020; 39: 1221–1226.

25.

Fernandes

Luís

Cardigos

, et al. Corneal subbasal nerve plexus evaluation by in vivo corneal confocal microscopy in multiple sclerosis: a potential new biomarker. Curr Eye Res 2021; 46: 1452–1459.

26.

Petropoulos

Fitzgerald

Oakley

, et al. Corneal confocal microscopy demonstrates axonal loss in different courses of multiple sclerosis. Sci Rep 2021; 11: 21688.

27.

Kass-Iliyya

Javed

Gosal

, et al. Small fiber neuropathy in Parkinson’s disease: A clinical, pathological and corneal confocal microscopy study. Parkinsonism Relat Disord 2015; 21: 1454–1460.

28.

Ponirakis

Al Hamad

Sankaranarayanan

, et al. Association of corneal nerve fiber measures with cognitive function in dementia. Ann Clin Transl Neurol 2019; 6: 689–697.

29.

Bitirgen

Akpinar

Uca

, et al. Progressive loss of corneal and retinal nerve fibers in patients with multiple sclerosis: A 2-year follow-up study. Transl Vis Sci Technol 2020; 9: 37–37.

30.

Bossuyt

Reitsma

Bruns

, et al. STARD 2015: an updated list of essential items for reporting diagnostic accuracy studies. Clin Chem 2015; 61: 1446–1452.

31.

Polman

Reingold

Banwell

, et al. Diagnostic criteria for multiple sclerosis: 2010 revisions to the McDonald criteria. Ann Neurol 2011; 69: 292–302.

32.

Kurtzke

JF.

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

33.

Roxburgh

Seaman

Masterman

, et al. Multiple Sclerosis Severity Score Using disability and disease duration to rate disease severity. Neurology 2005; 64: 1144–1151.

34.

Petropoulos

Alam

Fadavi

, et al. Corneal nerve loss detected with corneal confocal microscopy is symmetrical and related to the severity of diabetic polyneuropathy. Diabetes Care 2013; 36: 3646–3651.

35.

Abalo-Lojo

Treus

Arias

, et al. Longitudinal study of retinal nerve fiber layer thickness changes in a multiple sclerosis patients cohort: a long-term 5 year follow-up. Mult Scler Relat Disord 2018; 19: 124–128.

36.

Talman

Bisker

Sackel

, et al. Longitudinal study of vision and retinal nerve fiber layer thickness in multiple sclerosis. Ann Neurol 2010; 67: 749–760.

37.

Petropoulos

Manzoor

Morgan

, et al. Repeatability of in vivo corneal confocal microscopy to quantify corneal nerve morphology. Cornea 2013; 32(5): e83–e89.

38.

Chen

Graham

Dabbah

, et al. Small nerve fiber quantification in the diagnosis of diabetic sensorimotor polyneuropathy: comparing corneal confocal microscopy with intraepidermal nerve fiber density. Diabetes Care 2015; 38: 1138–1144.

39.

Evangelou

Konz

Esiri

, et al. Size-selective neuronal changes in the anterior optic pathways suggest a differential susceptibility to injury in multiple sclerosis. Brain 2001; 124: 1813–1820.

40.

Ponirakis

Abdul-Ghani

Jayyousi

, et al. Effect of treatment with exenatide and pioglitazone or basal-bolus insulin on diabetic neuropathy: a substudy of the Qatar Study. BMJ Open Diabetes Res Care 2020; 8: e001420.

41.

Schirmer

Merkler

König

, et al. Neuroaxonal regeneration is more pronounced in early multiple sclerosis than in traumatic brain injury lesions. Brain Pathol 2013; 23: 2–12.

42.

Ahmed

Bril

Orszag

, et al. Detection of diabetic sensorimotor polyneuropathy by corneal confocal microscopy in type 1 diabetes. Diabetes Care 2012; 35: 821–828.

43.

Rothman

Murphy

Fitzgerald

, et al. Retinal measurements predict 10-year disability in multiple sclerosis. Ann Clin Transl Neurol 2019; 6: 222–232.

44.

Martinez-Lapiscina

Arnow

Wilson

, et al. Retinal thickness measured with optical coherence tomography and risk of disability worsening in multiple sclerosis: a cohort study. Lancet Neurol 2016; 15: 574–584.

45.

Balk

Cruz-Herranz

Albrecht

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

46.

Nolan-Kenney

Liu

Akhand

, et al. Optimal intereye difference thresholds by optical coherence tomography in multiple sclerosis: An international study. Ann Neurol 2019; 85: 618–629.

47.

Confavreux

Vukusic

Natural history of multiple sclerosis: a unifying concept. Brain 2006; 129: 606–616.

48.

Dehghani

Pritchard

Edwards

, et al. Morphometric stability of the corneal subbasal nerve plexus in healthy individuals: a 3-year longitudinal study using corneal confocal microscopy. Invest Ophthalmol Vis Sci 2014; 55: 3195–3199.