Impairment of odor discrimination and identification is associated with disability progression and gray matter atrophy of the olfactory system in MS.

BACKGROUND: Impairment of odor discrimination (D), identification (I), and threshold (T) are characteristic features of multiple sclerosis (MS).
OBJECTIVE: To identify patterns of gray matter concentration (GMC) associated with different qualities of olfactory function.
METHODS: Olfactory function (T and combined DI score) was measured by Sniffin' Sticks-Test over 2 years longitudinally, and T1-weighted 3-T magnetic resonance imaging (MRI) was performed in 37 MS patients and 18 matched healthy controls (HCs). Statistical parametric mapping (SPM) was applied to objectively identify changes of voxel-wise-GMC throughout the entire brain volume and to correlate image parameters with odor function.
RESULTS: SPM localized significant GMC decreases in the anterior cingulum as well as temporomesial and frontobasal brain areas of the MS group compared with HCs, and revealed significant correlations between lower DI scores and GMC decreases in the olfactory gyrus, anterior cingulum, temporal regions including the parahippocampus, and putamen. Contrarily, no correlations were found between T and GMC. Patients with disability progression had significantly lower mean temporomesial/putamen GMC (0.782 vs 0.804, p = 0.004) compared to patients without Expanded Disability Status Scale (EDSS) progression.
CONCLUSION: Impairment of DI, but not T is associated with GM atrophy in brain regions related to olfactory function. Further studies are warranted to investigate DI scores and temporomesial/putamen GMC as biomarkers for disability progression.

Introduction

Dysfunction of various qualities of the olfactory sense has been increasingly recognized in multiple sclerosis (MS) with prevalences ranging from 30% to 75%, mostly depending on different MS populations and testing methods.[1,2] The current paradigm of MS pathology comprises a disease process involving both inflammation and neurodegeneration with the former being the prominent feature of early and relapsing disease phases while the latter is deemed the key contributor to accumulation of disability and development of progressive disease stages.[3]

The capacity to correctly identify odors (identification) and to discriminate them (discrimination) is affected in MS patients suffering from disability progression, whereas olfactory threshold is transiently altered during acute relapses in early and active MS resolving in the absence of a relapse.[4,5] In this regard, the sum score of discrimination and identification (DI score) was shown to be correlated with disease duration, physical disability, reduced cognitive function, and reduced retinal thickness of MS patients.[6] Hence, odor identification and discrimination were suggested to serve as parameters reflecting neurodegeneration in MS, while decreased odor thresholds indicate short-term inflammatory disease activity.[4,5]

Structural assessment of olfactory dysfunction using magnetic resonance imaging (MRI) was shown to be associated with the number and activity of MS-related plaques in frontal and temporal brain regions, atrophy of the olfactory bulb and tract, a decrease in white matter integrity, and recently atrophy of olfactory-related gray matter volume.[2,7-12]

However, these studies have predominantly considered olfactory dysfunction as one domain and hence did not investigate associations of different qualities of olfactory function such as odor identification, discrimination, and threshold, and the corresponding brain structure. Also, there is a lack of longitudinal studies investigating gray matter volume and olfactory function.

Thus, in this study we aimed to identify patterns of gray matter alterations associated with different qualities of olfactory function applying voxel-based morphometry of gray matter concentration (GMC), an early marker of structural brain alteration, to a cohort of early-staged MS patients in a longitudinal design.

Methods

Between May 2015 and February 2016, we recruited 37 MS patients and 18 healthy controls (HCs) from the MS clinic of the Department of Neurology at the Medical University Innsbruck. HCs did not report any health-associated complaints. The study was approved by the ethics committee of the Medical University Innsbruck (ethical approval number: AM3743-281/4.3) and all participants gave written informed consent before inclusion. MS patients were diagnosed according to the 2010 McDonald’s criteria.[7] Patients and HCs had to be aged between 18 and 55 years. Exclusion criteria were cognitive impairment defined as a score of 25 or lower in the Mini-Mental State Examination, or a history of chronic oto-rhino-laryngeal disease (such as chronic rhinitis, nasal polyposis, or sinus disease), head trauma, toxic exposures, upper respiratory infections at the time of assessment, previous radiation, nasal and/or oral surgery, and symptoms of nasal obstruction or other neurological diseases known to be associated with olfactory disturbances.

Study visits were conducted at baseline and at 2 years of follow-up. At every visit, demographic data, smoking habits, neurological and pharmacological history including occurrence and date of relapses and disease-modifying therapy, Expanded Disability Status Scale (EDSS), subjective olfactory function, and use of drugs and hormonal contraceptives were obtained from each participant.[8] EDSS progression was defined as a confirmed EDSS increase of ⩾1.0 point in patients with a baseline score of ⩽5.5, or an increase of ⩾0.5 points in patients with a baseline score of >5.5 sustained for at least 12 months as compared to baseline.[9] If a relapse, verified by the treating neurologist, had occurred within 6 months, EDSS was only considered when confirmed after 6 months following the visit.

The Beck Depression Inventory (BDI) was performed at each visit to screen for depression. Depression was defined as a score of 18 or higher on the BDI.[10]

Olfactory function was assessed at each visit using the extended version of the Sniffin’ Sticks test (Burghart Medizintechnik, Wedel, Germany) according to the manufacturer’s instruction including change of testing sticks every 6 months.[11] The Sniffin’ Sticks is a three-stage test, based on pen-like odor-dispensing devices. In stage 1, odor threshold is assessed using n-butanol as a single odorant. Three sticks are presented to each subject in a randomized order, two containing solvent and the third containing the odorant at a certain dilution. The subject is asked to identify the stick with the odorant. In stage 2, the subject is presented with three pens and asked to discriminate between two pens containing the same odorant and one containing a different one. In stage 3, the ability of the subject to identify an odorant is assessed by testing 16 different odors from a single pen by a forced choice from four options. The threshold, discrimination, and identification scores were rated for each patient. The maximum score at each stage is 16 points and reflects optimal olfactory function. Lower scores are associated with impaired ability to identify and discriminate odors and an increased threshold for odor perception. The overall TDI score is not ideal for quantification of olfactory dysfunction in MS as it displays less robust correlations to clinical variables compared to its subscores.[4] Alterations in threshold are masked by better performance in discrimination/identification or vice versa. Therefore, we regarded threshold as a separate domain to discrimination and identification summed in a composite score (DI score) providing better correlation to clinical variables than each subscore alone.[4]

The normative values are based on data from 3000 healthy subjects.[12] Olfactory testing was postponed for 4 weeks, if the patient had a relapse or received corticosteroids within 4 weeks or if upper respiratory tract infections were present at the time of assessment.

MRI data acquisition

All MRI measurements were performed on a 3.0 T whole-body MR scanner (Magnetom Skyra, Siemens, Erlangen, Germany) equipped with a 20-channel head coil. The MRI protocol included 3D-T1-weighted and T2-SPACE sequences. T1-MPRAGE sequence parameters were 160 slices TR = 1850 milliseconds, TE = 3.01 milliseconds, inversion time = 900 milliseconds, slice thickness = 1 mm, matrix = 256 × 204 pixels, number of excitations = 1, flip angle = 8, and FOV = 220 mm × 178.75 mm. T2-SPACE parameters were set to the following: 160 slices, slice thickness = 1 mm, TR = 4800 milliseconds, TE = 388 milliseconds, flip angle = 120, and FOV = 230 × 222.81 mm.

Image postprocessing

T2-SPACE and T1-MPRAGE images were brain extracted. Registration of T2-SPACE and T1-MPRAGE was conducted by the FMRIB’s linear image registration tool FLIRT implemented into FSL software (FMRIB library 5.0, Oxford, UK).[13] An automated lesion prediction algorithm as implemented in the lesion segmentation toolbox for statistical parametric mapping (SPM) was used to segment T2 hyperintense and gray matter hypointense lesions in SPACE and T1-weighted images. The tool gives total lesion volume and lesion number. The segmented lesions were reviewed by an expert (C.S.).[14-16]

To avoid a priori assumptions through region of interest (ROI) analysis on brain areas of potential interests, gray matter density measures were subjected to SPM (SPM12, Wellcome Department of Cognitive Neurology, London, UK) implemented in MATLAB 9.2 (MathsWorks Inc., Sherborn, MA, USA), a technique that objectively localizes focal changes of voxel values throughout the entire brain volume.[14] Voxel-based morphometry of the gray matter compartment was performed by applying the standard version of the diffeomorphic anatomical registration using exponentiated lie algebra toolbox implemented in SPM12 to have a high-dimensional normalization protocol.[16]

Segmented and unmodulated images were transformed from the study-specific diffeomorphic anatomical registration space into Montreal Neurological Institute (MNI) space and smoothed by a Gaussian kernel of 8 mm × 8 mm × 8 mm in order to accommodate interindividual anatomic variability and to improve signal-to-noise ratios for the statistical analysis. Brodmann areas (BAs) were obtained by the MNI atlas implemented into the software package MRICRON.[17] A masking threshold of 10% of the lower image signal was applied to reduce signal noise.

Statistical analysis

Statistical analysis was performed using SPSS 25.0 (SPSS Inc, Chicago, IL, USA). Categorical variables were expressed in frequencies and percentages. Continuous variables were tested for normal distribution by Kolmogorov–Smirnov test. In case of normal distribution, variables were displayed as mean and standard deviation or 95% confidence interval. If variables were not normally distributed, they were displayed as median and range. Univariate comparisons were done by chi-square-test, Fisher’s exact test, Mann–Whitney U test, or independent t-test as appropriate. Effect sizes have been calculated using Cohen’s d values where appropriate.

The obtained datasets allowed for comparisons of GMC in analogous voxel regions between healthy volunteers and MS patients. A general linear model was set up to compare GMC values between groups using a two-tailed t-contrast. MRI parameters and scores of olfactory testing (threshold and DI score) were entered into the design matrix of the general linear model and examined by t-contrast. All results were corrected for multiple comparisons at family-wise-error I (FWE) p < 0.05 with the statistical threshold set to p < 0.001 for group comparison and correlation analysis. Age and total intracranial volume were entered as covariates.

A post hoc t-test analysis of the mean GMC of the voxel-cluster displaying significant positive correlations of GMC decreases and DI score reductions was performed. Since the olfactory cortex refers to a constellation of processing modules distributed predominantly throughout the mediotemporal and inferofrontal cortices, the respective voxel-cluster was separated into a frontal and a temporomesial/putamen region with the regional boundaries defined by the MNI atlas.[18-20]

Subsequently, GMC values of these two ROIs were compared between patients with and without EDSS progression within 2 years from baseline using t-test analysis.

In an explorative approach, we performed receiver-operating characteristic (ROC) curve analysis to identify the best possible cutoff values discriminating between patients with and without EDSS progression in the 2-year span after the baseline visit, respectively.

Results

Characteristics of MS patients are displayed in Table 1. The patient and HC groups were matched for age and sex. Five MS patients were excluded from the SPM analyses for insufficient MRI quality. The cumulative white matter lesion volume of the entire MS cohort was 164.16 ccm. The filled-in gray matter lesion volume ranged from 0.01 to 3.74 ccm (mean 1.16 ± 1.09 ccm).

Table 1.

Characteristics of MS patients and HC.

	MS (n = 37)	HC (n = 18)	p-Value
Female[a]	29 (78.4%)	14 (77.7%)	0.999[d]
Age at baseline (years)[b]	33.8 (6.9)	30.6 (6.6)	0.107[e]
Disease duration (years)[b]	5.5 (4.4)	na
RRMS[a]	33 (89.2%)	na
SPMS[a]	3 (8.1%)	na
PPMS[a]	1 (2.7%)
Number of relapses in last year[b]	0.73 (0.93)	na
EDSS at baseline[c]	2.0 (0–4.5)	na
EDSS at year 2[c]	2.5 (0–5.5)	na
EDSS progression baseline—year 2[a]	15 (40.5%)	na
Treated with DMT[a]	32 (86.5%)	na
Interferon beta[a]	13 (35.1%)	na
Glatiramer acetate[a]	10 (27.0%)	na
Fingolimod[a]	1 (2.7%)	na
Natalizumab[a]	8 (21.6%)	na
Beck Depression Index[b]	5.1 (5.6)	na
Depression (BDI ⩾18)[a]	2 (5.4%)	na
MMSE[c]	30 (28–30)	na
Smokers[a]	11 (29.7%)	na
Baseline threshold[c]	7.0 (4.0–10.0)	na
Threshold age-adjusted hyposmia[a]	15 (40.5%)	na
Baseline discrimination[c]	13 (11–16)	na
Discrimination age-adjusted hyposmia[a]	4 (10.8%)	na
Baseline identification[c]	13 (9–16)	na
Identification age-adjusted hyposmia[a]	5 (13.5%)	na
2-year threshold[c]	7.25 (4.25–10.0)	na
Threshold age-adjusted hyposmia[a]	15 (40.5%)	na
2-year discrimination[c]	12 (11–16)	na
Discrimination age-adjusted hyposmia[a]	7 (18.9%)	na
2-year identification[c]	11 (8–16)	na
Identification age-adjusted hyposmia[a]	9 (24.3%)	na

BDI: Beck Depression Index; DMT: disease-modifying therapy; EDSS: Expanded Disability Status Scale; HC: healthy controls; MMSE: Mini-Mental Status Examination; MS: multiple sclerosis; Na: not applicable; PMS: progressive multiple sclerosis; RRMS: relapsing-remitting multiple sclerosis.

Number (percentage).

Mean and standard deviation.

Median and range.

Calculated by Fisher’s exact test.

Calculated by t-test.

While DI score significantly decreased during the observation period, olfactory threshold did not change significantly (Table 1, Figure 1).

Figure 1.

Development of olfactory threshold (panel a), odor discrimination (b), identification (c), and combined DI score (d).

p-Values calculated by Mann–Whitney U test. Cyan lines indicate median scores.

SPM localized significant decreases in GMC between MS patients and HCs, most prominent bilaterally in the cingulum (BA 23 and 24), parahippocampal gyri (BA 28), and in the right straight gyrus, frontal superior and medial orbital gyrus, and olfactory gyrus (BA 11), at FWE corrected p < 0.001 (Supplemental Table 1, Supplemental Figure 1).

Decreases of the DI score correlated significantly with decreased GMC in the right olfactory gyrus and anterior cingulate (BA 11), the middle temporal gyrus and temporal pole (BA 20, 21), the left parahippocampus (BA 35), and bilaterally in the amygdala (BA 34) and the putamen (p < 0.001; FWE corrected; Table 2, Figure 2). No significant associations were found between olfactory threshold scores and GMC.

Table 2.

Brain regions of significant positive correlations of gray matter concentration decreases and combined olfactory function discrimination and identification score in patients with multiple sclerosis.

Cerebral region	Cluster size (mm3)	MNI coordinates			t-value	Cohen’s d	Height threshold
		x	y	z
Right middle temporal gyrus and temporal pole (BA 20, 21)	7780	−50	−8	−18	4.6	1.91	0.001
Right putamen		−34	−8	−5	4.45	1.85
Right amygdala		22	−2	−15	4.31	1.79
Right frontal superior and medial orbital gyrus		−15	35	−20	4.57	1.9
Straight gyrus (BA 11)		−9	33	−18	4.2	1.74
Right olfactory gyrus (BA 25)		−10	12	−15	3.77	1.57
Right anterior cingulate (BA 11)		−8	39	−4	4.2	1.76
Left amygdala (BA 34)	3490	29	2	−18	4.41	1.84
Left putamen		29	−2	−4	4.08	1.70
Left parahippocampus (BA 35)		19	−4	−27	3.75	1.56

BA: Brodmann area; MNI: Montreal Neurological Institute coordinates.

Cerebral regions are family-wise error corrected at the threshold of 0.001 of the cluster level.

Figure 2.

Correlations of gray matter concentration decreases and the combined olfactory function discrimination and identification score. Statistical parametric mapping (t-value) axial intensity projection maps showing areas of significant positive correlations of gray matter concentration decreases and the combined olfactory function discrimination and identification score (upper panel) and its effect size representing Cohen’s d values (lower panel) in a cohort of patients with multiple sclerosis rendered on to a stereotactically normalized magnetic resonance imaging scan. The number at the bottom right corner of each MRI scan corresponds to the z coordinate in Montreal Neurological Institute space.

In both the frontal region and the temporomesial/putamen ROIs, GMC correlated significantly with DI score but not with threshold (Figure 3). Patients with EDSS progression had significantly lower mean GMC in the temporomesial/putamen ROI (0.782 vs 0.804, p = 0.004, Cohen’s d = 2.103) compared to patients without EDSS progression, while the frontal ROI was not statistically significant (0.688 vs 0.692, p = 0.074, Cohen’s d = 1.478). Mean DI score was significantly lower in patients suffering an EDSS progression (23.2 vs 27.0, p = 0.004, Cohen’s d = 2.354), while threshold scores did not differ (6.8 vs 7.4, p = 0.345, Cohen’s d = 0.921).

Figure 3.

Gray matter concentration in relation to olfactory threshold (a) and (c) and DI score (b) and (d).

DI score: combined score of discrimination and identification; GMC: gray matter concentration.

p-Values calculated by general linear models.

The DI score displayed the best discrimination ability at a cutoff lower than 24 points (specificity 94%, sensitivity 63%, area under the ROC curve (AUC) 0.86) followed by temporomesial/putamen GMC (specificity 82%, sensitivity 58%, AUC 0.81) with the cutoff at <0.79 (Table 3).

Table 3.

ROC analysis for prediction of EDSS progression within 2 years.

	Cutoff value	Specificity (%)	Sensitivity (%)	AUC (SD)
Frontal ROI GMC	<0.69	74	34	0.68 (0.53–0.78)
Temporomesial/putamen ROI GMC	<0.79	82	58	0.81 (0.65–0.93)
DI-score	<24	94	63	0.86 (0.72–0.96)

AUC: area under the ROC curve; DI: discrimination and identification sum score; EDSS: expanded Disability Status Scale; GMC: gray matter concentration; ROC: receiver-operating characteristic; ROI: region of interest; SD: standard deviation.

Discussion

In this study, we identified a robust correlation between impairment of DI score and decreased GMC in the putamen and temporomesial brain regions in MS patients. On the contrary, we did not find any significant associations between olfactory threshold and GMC alterations.

In recent years, detailed assessments of odor function in MS revealed the DI score as strongly correlated with clinical measures of neurodegeneration (i.e. EDSS, cognitive function) mostly affected in patients suffering from disability progression, but independent of inflammatory disease activity (i.e. occurrence and number of relapses).[1,4,6,21] On the contrary, impairment of olfactory threshold is correlated with short-term inflammatory relapse activity, transiently resolving in the absence of relapse activity but independent of parameters of neurodegeneration.[4-6,22]

However, the underlying pathophysiology is unclear. It has been hypothesized that olfactory threshold is a function of more peripheral parts of the olfactory system (“cable function”) affected by a bystander inflammation of the olfactory tract during phases of clinical disease activity, while discrimination and identification are complex cortical functions (“processing functions”) likely to be mainly affected by neurodegenerative processes.[4,22] The plausibility of this concept is strengthened by observations regarding the transient nature of threshold impairment compared to the irreversible affection of discrimination and identification.[4,5] Also, retinal nerve fiber layer thickness, which is a reliable and robust surrogate biomarker of neuroaxonal degeneration, is correlated with DI score, but not with threshold.[6,23]

Gray matter atrophy quantified by means of MRI is a well-established and widely accessible imaging marker of neurodegeneration in MS.[24] It was shown to occur in all MS phenotypes and to be associated with disability accumulation.[25] In line with our findings, brain areas displaying early atrophy include the cingulate cortex, insular and temporal cortical GM, and the deep gray matter (i.e. putamen and caudate).[25,26] Interestingly, histopathological examinations revealed the neostriatum (putamen and caudate) and the deep gray matter as brain regions exhibiting the greatest extent of demyelination in early MS.[27] Since the striatum receives significant inputs from the motor cortex and the association cortices, it was speculated that enhanced vulnerability of these structures may arise from retrograde neurodegeneration secondary to higher lesion load.[25]

While former studies applied gray matter volume measures, a recent study investigated GMC and successfully revealed subtle changes in subcortical and frontoparietal regions of MS patients displaying little signs of physical disability (EDSS score <3.0) but cognitive impairment. Hence, reduction of GMC seems to precede overt brain atrophy measured by volumetric methods.[28]

In accordance with previously reported morphometric investigations in MS, categorical voxel-based analysis identified significant GMC decrease in the right olfactory gyrus, the anterior cingulate (BA 11), the middle temporal gyrus and temporal pole (BA 20, 21), the left parahippocampus (BA 35), and bilaterally in the amygdala (BA 34) and the putamen in MS patients compared to HCs.[25,26]

Linking GM atrophy and olfactory dysfunction in MS, the volume of “olfactory-related GM” (predominantly in the parahippocampal gyrus, PCG) was previously reported to be decreased in patients with olfactory dysfunction.[22] However, this study did not differentiate between odor identification, discrimination, and threshold. In our study, we identified a highly significant association between PCG density and DI scores, but not olfactory threshold scores. The PCG along with the amygdala and parts of the anterior and posterior cingulate is involved in olfactory memory, semantic memory and recollection, and odor familiarity, features that are assessed by DI but not threshold.[29,30]

Current evidence suggests that olfactory dysfunction in MS has to be considered as a result of a dual pathology: (1) alteration of discrimination and identification reflects GM atrophy in regions directly related to olfactory function and disease progression, whereas (2) impairment of olfactory threshold does not (and might instead reflect inflammatory disease activity). Our study is supporting this concept by providing for the first time a direct link between impairment of DI score, disability progression, and GM atrophy in a longitudinal design.

In this context, dysfunction of distinct olfactory qualities is emerging as a biomarker of both inflammation and subsequent neurodegeneration. The spectrum of disease-modifying treatment options has significantly broadened over the last decade and, as a consequence, treatment decisions have become much more complex. To enable evidence-based decision making, it is critical to have biomarkers for the assessment of individual risk for short- to mid-term disability progression. Clinical disability assessment by the EDSS reflects neurodegenerative damage incompletely, with substantial delay and low sensitivity and does not provide a linear marker, as greater rates of change are observed for lower EDSS scores.[31]

Based on our results, both DI score and temporomesial/putamen GMC are suggested candidates for biomarkers for the prediction of short- to mid-term disability progression. However, as this is an explorative study with small sample size, our results certainly have to be interpreted cautiously and require validation in an independent larger cohort.

In conclusion, impairment of odor discrimination and identification, but not threshold, is associated with gray matter atrophy in regions related to olfactory function. This finding suggests to validate the assessment of DI score and temporomesial/putamen GMC as biomarkers to predict disability progression and monitor neurodegeneration in MS.

Click here for additional data file.

Supplemental material, MSJ838205_supplementary_figures for Impairment of odor discrimination and identification is associated with disability progression and gray matter atrophy of the olfactory system in MS by Gabriel Bsteh, Ruth Steiger, Noora Tuovinen, Harald Hegen, Klaus Berek, Sebastian Wurth, Michael Auer, Franziska Di Pauli, Elke R Gizewski, Florian Deisenhammer, Thomas Berger and Christoph Scherfler in Multiple Sclerosis Journal

Click here for additional data file.

Supplemental material, MSJ838205_supplementary_table_1 for Impairment of odor discrimination and identification is associated with disability progression and gray matter atrophy of the olfactory system in MS by Gabriel Bsteh, Ruth Steiger, Noora Tuovinen, Harald Hegen, Klaus Berek, Sebastian Wurth, Michael Auer, Franziska Di Pauli, Elke R Gizewski, Florian Deisenhammer, Thomas Berger and Christoph Scherfler in Multiple Sclerosis Journal