Article Text

Download PDFPDF

MR spectroscopy indicates diffuse multiple sclerosis activity during remission
  1. I I Kirov1,
  2. V Patil1,
  3. J S Babb1,
  4. H Rusinek1,
  5. J Herbert2,
  6. O Gonen1
  1. 1
    Department of Radiology, New York University School of Medicine, New York, USA
  2. 2
    Department of Neurology, New York University School of Medicine, New York, USA
  1. Correspondence to Professor O Gonen, Department of Radiology, New York University School of Medicine, 660 First Avenue, 4th Floor, New York, NY 10016, USA; oded.gonen{at}nyumc.org

Abstract

Objective: To test the hypothesis that diffuse abnormalities precede axonal damage and atrophy in the MRI normal-appearing tissue of relapsing-remitting (RR) multiple sclerosis (MS) patients, and that these processes continue during clinical remission.

Methods: Twenty-one recently diagnosed mildly disabled (mean disease duration 2.3 years, mean Expanded Disability Status Scale score of 1.4) RR MS patients and 15 healthy matched controls were scanned with MRI and proton MR spectroscopic imaging (1H-MRSI) at 3 T. Metabolite concentrations: N-acetylaspartate (NAA) for neuronal integrity; choline (Cho) for membrane turnover rate; creatine (Cr) and myo-inositol (mI) for glial status were obtained in a 360 cm3 volume of interest (VOI) with 3D multivoxel 1H-MRSI. They were converted into absolute amounts using phantom replacement and normalised into absolute concentrations by dividing by the VOI tissue volume fraction obtained from MRI segmentation.

Results: The patients’ mean VOI tissue volume fraction, 0.92 and NAA concentration, 9.6 mM, were not different from controls’ 0.94 and 9.6 mM. In contrast, the patients’ mean Cr, Cho and mI levels 7.7, 1.9 and 4.1 mM were 9%, 14% and 20%, higher than the controls’ 7.1, 1.6 and 3.4 mM (p = 0.0097, 0.003 and 0.0023).

Conclusions: The absence of early tissue atrophy and apparent axonal dysfunction (NAA loss) in these RR MS patients suggests that both are preceded by diffuse glial proliferation (astrogliosis), as well as possible inflammation, demyelination and remyelination reflected by elevated mI, Cho and Cr, even during clinical remission and despite immunomodulatory treatment.

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

The sensitivity of T1- and T2-weighted MRI to white matter (WM) lesions and to regional as well as global atrophy has made this modality central to the diagnosis and treatment monitoring of multiple sclerosis (MS).1 2 3 4 Paradoxically, these findings correlate only weakly with clinical disability,5 an incongruity due to MRI’s insensitivity to microscopic pathology and lack of specificity to distinguish inflammation from demyelination, axonal loss or gliosis.1 6 Since MS is now recognised to be a diffuse disease of both normal-appearing WM (NAWM) and grey matter (GM),7 8 these sensitivity and specificity issues limit the potential of clinical MRI to fully assess the total load of the disease and consequently, to monitor treatment response.6

Although diffuse pathology that leads to microscopic axonal damage is often MRI-occult, it is nevertheless detectable with quantitative MR methods.6 For example, proton MR spectroscopy (1H-MRS) can probe tissue metabolism,9 whereas magnetisation transfer and diffusion tensor imaging can assess myelin and WM integrity.1 These have shown that NAWM abnormalities may (1) be present in all MS subgroups, (2) correlate with disability and cognitive impairment, and (3) occur prior to lesion development, all suggesting that they play a role in lesion formation.6 7

1H-MRS offers the ability to monitor several pathological processes via their specific metabolic markers: N-acetylaspartate (NAA) for neuronal integrity, whose decline correlated better with disability than MRI metrics;10 myo-inositol (mI) for glial proliferation,6 whose elevation has also correlated with impairment;11 creatine (Cr) for energy status and cell density and choline (Cho) for inflammation, demyelination and remyelination.10 Unfortunately, most studies to date use single-voxel or 2D 1H-MRS that sample under 10% of the brain, subjecting them to the implicit assumption that such small volumes are globally representative. In addition, quantification with metabolite ratios may have confounded interpretation since Cr (the frequent denominator) levels are reported also to vary.12 13

Only one study addressed both issues with absolute quantification in a large (∼½ litre) WM volume of relapsing-remitting (RR) MS patients and controls.12 It reported lower NAA and elevated Cho and Cr in the patients but could not ascribe the latter to either astrogliosis or myelin repair, since Cr levels are higher in both astrocytes and oligodendrocytes;14 and the long echo-time 1H-MRS used precluded detection of the astrocyte-specific marker—mI.6 Furthermore, the NAA loss could not distinguish isolated event(s) from chronic inflammation due to the long, 6+ years, disease duration. Our goal, therefore, is to quantify the glial, demyelination and remyelination markers (mI, Cho and Cr) in NAWM of recently diagnosed RR MS patients with short, TE = 35 ms, 1H-MRS to test the hypothesis that these processes precede axonal damage (NAA).

Materials and methods

Human subjects

Twenty-one (15 women) patients meeting the Poser criteria for RR MS were enrolled.15 Their mean age was 33 (range 21–42) years; the mean disease duration was 2.3 (range 0.3–5.2) years; and the mean Expanded Disability Status Scale (EDSS16) score was 1.25 (range 0–4), as shown in table 1. All patients were considered to have mild MS (EDSS 4.0 or less), since none demonstrated any limitation in ambulation. One patient (11 in the table) scored EDSS 4.0 due to permanent visual loss in one eye. Of the remaining 20 patients, 13 had clinical signs on examination consisting of hyper-reflexia, tandem gait ataxia or minimal sensory abnormality, but no clinical symptoms or disability (EDSS 1.5 or less); five patients presented with mild disability in only one functional system (EDSS 2.0), manifested as mild ataxia, nystagmus or impaired visual acuity; and only two patients demonstrated mild disability in more than one functional system: patient 3, EDSS 3.5 (monoparesis and impaired visual acuity), and 17, EDSS 2.5 (paraparesis and sensory disturbance). All were on immunomodulatory medication and none experienced a relapse in the preceding 3 months.

Table 1

Demographic, clinical, volumetric and metabolic data for patients (1–21, sorted by disease duration) and controls (22–36)

Fifteen age- and gender-matched controls (12 women, mean age 30, range 20–44 years) were also enrolled. Their “healthy” status was based on negative answers to a list of disqualifying neurological conditions before the study and unremarkable MRI afterwards. All participants gave written Institutional Review Board-approved informed consent.

MR acquisition

All measurements were done in a 3 T Trio MRI scanner (Siemens AG, Erlangen, Germany) with a TEM3000 transmit-receive head coil (MRInstruments, Minneapolis, Minnesota). For tissue segmentation, sagittal MP-RAGE (TE/TI/TR = 2.6/1800/1360 ms) was acquired at a field of view (FOV) of 240×240 mm2, 256×256 matrix and reconstructed in 1 mm thick axial and coronal slices. To guide the spectroscopic volume of interest (VOI), T2-weighted fluid attenuated inversion recovery (FLAIR) (TE/TI/TR = 88/2500/9000 ms) images were obtained at the same FOV and matrix but 3.7 mm slice thickness.

Following our chemical-shift imaging (CSI)-based shim procedure, a 10 cm anterior–posterior (AP)×8 cm left–right (LR)×4.5 cm inferior–superior (IS) = 360 cm3 1H-MRSI VOI was centred on the corpus callosum (fig 1) and excited using TE/TR = 35/1800 ms PRESS. The VOI was partitioned using Hadamard spectroscopic imaging17 into six axial slices which were encoded with 16×16 2D-CSI over a 16×16 cm2 (LR×AP) FOV. The VOI was defined in the planes of these slices by the two PRESS 11.2 ms long numerically optimised 180° pulses yielding 8×10×6 = 480 voxels, 1.0×1.0×0.75 cm3 each. The signal was acquired for 256 ms at ±1 kHz bandwidth. At two averages, the 1H-MRSI took 34 min, and the entire protocol lasted 1 h.

Figure 1

Left: sagittal (a) and coronal (b) T1-weighted and axial T2-weighted fluid attenuated inversion recovery (c) MRI of patient 12 in table 1, with the volume of interest and field of view (solid and dashed white frames) superimposed. Note characteristic periventricular hyperintensities on (c) (arrows) and little or no atrophy. Right: real part of the 8×10 (left–right (LR)×anterior–posterior (AP)) 1H spectra matrix from the volume of interest on (c), on common frequency (1.4 to 3.8 ppm) and intensity scales. Note the spectral resolution and signal-to-noise ratio in these 0.75 cm3 voxels and elevated choline (Cho) and myo-inositol (mI) at the lesions (black arrows). CSI, chemical-shift imaging; HSI, Hadamard Spectroscopic Imaging; NAA, N-acetylaspartate.

Tissue fraction (Tf) and lesion volume in the VOI

Since the VOI includes ventricles and sulci (cf fig. 1) correction for their cerebrospinal fluid (CSF) partial volume is needed. To this end the MP-RAGE images were segmented using our MIDAS package.18 The process starts by automatic detection of the WM signal intensity, IWM, in a periventricular “seed” region. Following selection of all pixels at or above 55, but below 135% (to exclude the fat) of IWM, a tissue-mask is constructed per slice in three steps: morphological erosion, recursive region growth retaining pixels connected to the “seed”; and morphological inflation to reverse the effect of erosion. Pixels whose intensity was under 55% of IWM were defined as CSF. Tf for each subject was the product of the number of “tissue” pixels in the VOI and their volume divided by the nominal 360 cm3 VOI-volume. Lesion volumes were obtained from the FLAIR images using the same procedure, but applying a threshold of 120% of IWM.

Absolute metabolite quantification

Processed offline using in-house software, the MRS data were voxel-shifted to align the CSI grid with the NAA VOI, Fourier-transformed in the time, AP and LR dimensions and Hadamard reconstructed along the IS direction. Each of the VOI’s 480 spectra was frequency-aligned and zero-order phased in reference to the NAA peak and then summed. This retained the individual spectra (narrow) line width and improved the signal-to-noise-ratio (SNR) by 480½≈22.12

Relative levels of the ith (NAA, Cr, Cho, mI) metabolite in the jth subject were estimated from their peak area, Sij, using spectral modelling.19 Sij were scaled into concentrations, Cij, relative to signals from a 2 L sphere of Civitro = 12.5, 10.0, 3.0 and 7.5 mM NAA, Cr, Cho and mI in water:

Embedded Image

where SR is the sphere’s metabolite signal, Vj180° and VR180° are the RF voltages needed for a non-selective 1 ms 180° inversion pulse on subject and reference sphere. The Cijs were corrected for relaxation time differences in vivo (T1vivo, T2vivo) and in the phantom (T1vitro, T2vitro) with a factor:12

Embedded Image

using the T1vivo = 1.4, 1.3 and 1.2 s, T2vivo = 343, 172 and 248 ms reported at 3 T for NAA, Cr and Cho20 21 and T1vitro = 605, 336 and 235 ms, T2vitro = 483, 288 and 200 ms in the phantom. Since J-coupling modulation and low voxel SNR for mI prevent obtaining its T1 and T2, we assumed f = 1 for it in equation 2. Finally, the Cijs were divided by Tf to convert into tissue concentrations.

Statistical analysis

Analysis of variance (ANOVA) based on ranks was used to compare patients’ and controls’ NAA, Cho, Cr, mI and Tf. A separate analysis was conducted for each measure. In each case, the observed values were converted to ranks (to better satisfy underlying distributional assumptions) that were then used as the dependent variable. Each model included age and gender as subject level covariates and subject group (patient versus control) as the between-subjects factor. Pearson product moment and Spearman rank correlation coefficients were used to evaluate association of disease duration, EDSS and lesion volume in the VOI with each measure among patients. When a significant correlation was detected, least-squares regression was used to characterise the change in that measure as a function of disease duration. SAS version 9.0 (SAS Institute, Cary, North Carolina) was used for statistical calculations.

Results

Our shim procedure adjusted the whole-head water line width to 27 (4) Hz full width at half maximum (FWHM) among the 36 subjects, which improved to 21 (3) Hz in the VOI. Examples of the size and placement of the VOI and 1H spectra are shown in figs 1, 2. The SNRs, estimated as the peak height divided by the root-mean-square of the noise in 16 800 voxels (35 subjects ×480 spectra each), were: NAA = 24.0 (10.0), Cr = 13.8 (4.8), Cho = 13.3 (4.8) and mI = 3.8 (1.0) (mean (SD)). The metabolites’ FWHM line width, Δω, was 5.5 (1.2) Hz across the 16 800 voxels reflecting a T2* ( = 1/πΔω assuming Lorentzian lines) of 57.8 (13.2) ms.

Figure 2

Axial T2-weighted fluid attenuated inversion recovery image of the patient (13) with the heaviest lesion load, superimposed with the volume of interest (white frame) and 8×10 (left–right×anterior–posterior) 1H spectra on common frequency (1.4 to 3.8 ppm) and intensity scales. Note extensive hyperintensities and visible choline (Cho) and myo-inositol (mI) elevation in voxels within these lesions (arrows) and in the normal-appearing white matter (asterisks), that is, their overall elevation is truly diffuse, not merely multifocal. Cr, creatine; NAA, N-acetylaspartate.

The 480 summed VOI spectra from each of the 21 patients and 15 controls superimposed with their model functions are shown in fig 3. Their SNRs range from a low of ∼75 for mI to over 500 for the NAA. The average line width in the sums increased only slightly to 6.8 (0.6) Hz versus the individual voxels (vide supra). These demonstrate the advantages of the aligned summation strategy: preservation of individual voxel line width in the sum and substantial SNR boost.12

Figure 3

Real part of the 1H 480 summed spectra from the volume of interest (black lines) for each of the 36 subjects in table 1: patients (1–21) and controls (22–36, circled), superimposed with the fitted model functions (grey lines) on common scales. Note: (1) excellent 75 (mI)—500 (N-acetylaspartate (NAA)) signal-to-noise ratios and line widths of the sums compared with the individual voxels in fig 1 (benefit of summation and alignment); (2) the quality of the fit; and (3) the visual similarity in NAA levels between patients and controls versus elevated creatine (Cr), choline (Cho) and myo-inositol (mI).

Average absolute metabolic concentrations obtained from the summed spectra of fig 3 using equations 1, 2, Tf-s and lesion volumes are given in table 1. Coefficients of variation (CV = standard deviation/mean) in controls (NAA: 13%; Cr: 11%; Cho: 12%; mI: 15%) were all lower (better) than the average CV (NAA: 21%; Cr: 17%; Cho: 22%; mI: 17%) of the 30 reports of healthy adults reviewed by Kreis,22 possibly reflecting the improved precision that is consequent to the higher SNR. Box plots (fig 4) show that the patients’ means: Tf = 0.92 and NAA = 9.6 mM are not significantly different from the controls’ 0.94 and 9.6 mM. (Based on their standard deviations, the smallest difference detectable with 80% power at the 5% significance level is 0.04 for Tf and 1.25 mM for NAA; alternatively, to establish the significance of these observed differences with 80% power would require larger cohorts: 56 subjects in each for Tf and 1100 for NAA.) In contrast, the patients’ Cr, Cho and mI levels, 7.7, 1.9 and 4.1 mM, were significantly (9%, 14% and 20%, p = 0.0097, 0.003 and 0.0023) higher than the controls’ 7.1, 1.6 and 3.4 mM even after Bonferroni correction. No correlations were found between any metric and EDSS, lesion volume in the VOI, or time from diagnosis except a trend (p = 0.06) for Cr increase with time.

Figure 4

Box plots displaying the 25%, median and 75% (box), 95% (whiskers) and outliers (*) of the N-acetylaspartate (NAA), creatine (Cr), choline (Cho) and myo-inositol (mI) concentrations as well as tissue fraction (Tf), distributions in the volumes of interest of patients and controls. Note the absence of atrophy and NAA loss versus significant 9, 14 and 20% elevation of Cr, Cho and mI in the recently diagnosed patients.

Discussion

Lack of convincing correlation between disability and T2 lesion load in MS is often attributed to the inability of MRI to account for microscopic NAWM damage.5 Since lesions comprise on average less than 3% of the brain volume,5 it is possible that milder but widespread pathology influences disability. Indeed, quantitative metrics that assess all the tissue in a given volume, for example, the NAA level, correlate better with clinical impairment.23 Yet, despite its specificity, 1H-MRS sensitivity is limited by SNR and WM/GM/CSF partial volume. Furthermore, placement of a single VOI over a lesion or NAMW may miss distal damage propagated by Wallerian and retrograde degeneration.24 25 26

To address both SNR and limited coverage, we assessed the metabolite levels in a large (360 cm3) VOI of mostly NAWM. Summing all 480 frequency-aligned spectra in the VOI yielded excellent spectral resolution and sufficient SNR to detect subtle metabolic changes. This approach has previously shown decreased Tf and NAA reflecting axonal loss and dysfunction, as well as elevated Cr and Cho in 11 RR MS patients of 6+ years’ disease duration.12 Based on those findings, our goals were to (1) test for the presence of astrogliosis (via its mI marker27) and (2) ascertain the temporal relationship between demyelination, gliosis and NAA decline, in order to (3) test the hypothesis that widespread pathogenesis precedes axonal damage and atrophy.

The results support the hypothesis and suggest the following sequence of disease evolution: (1) astrogliosis, possibly accompanied by inflammation, demyelination and remyelination, preceding (2) axonal damage, that is perhaps followed by (3) atrophy.

Diffuse abnormalities

Since our cohort’s average VOI lesion load was under 2% (cf table 1), differences of 9%, 14% and 20% for Cr, Cho and mI from controls must reflect widespread NAWM pathology (fig 2). Moreover, the metabolic pattern, seen in the global VOI spectra of fig 3, is similar to that observed in the lesions (cf figs 1, 2), known foci of ongoing activity.9 28 This similarity may suggest that the same processes may also be widespread in the NAWM, early in the disease course, even during remission and despite immunomodulatory treatment. This conjecture is supported by strong correlations of MRS observed elevated mI and Cho with biopsy findings (astrogliosis and inflammation) in lesions,29 and the well-documented NAWM abnormalities—astrogliosis and inflammation in autopsied tissue.7

Glial pathology

While Cr and Cho are found in all brain cells, their levels are highest in mature and progenitor oligodendrocytes and astrocytes.14 Their elevation, therefore, may represent glial activity: combinations of inflammation, demyelination and remyelination as well as gliosis. Corroboration and specificity are provided by the mI, an osmolyte and constituent of phosphate- and lipid-containing compounds involved in signal transduction in the brain.30 Its 1H-MRS peak is thought to originate from intracellular stores in astrocytes31 and its increase, therefore, indicates astrogliosis,6 that is, cell hypertrophy and hyperplasia.32

Evidence of glial abnormalities at the earliest stages is offered by the largest (96 patients) 1H-MRS study of clinically isolated syndrome suggestive of MS (CIS).33 It reported elevated mI, Cr and Cho (but no NAA decline) most pronounced in patients who also satisfied the McDonald criteria. Indeed, elevated mI suggestive of astrogliosis is a consistent finding in all disease durations.11 34 35 While astroglial scars are known to form on severed axons, the cause of astrogliosis in the context of (MRS inferred) healthy neurons is unclear. It is known that astrocytes directly contact the endothelial cells comprising the blood–brain barrier (BBB) and enhance its impermeability.36 Both the subtle BBB breakdown from microscopic inflammation37 and hypoperfusion38 that were proposed as possible pathogenic mechanisms may conceivably cause astrogliosis: the former in response to increased leakage through the BBB and (or) microplaque formation; the latter due to oxygen deprivation, as shown in a rat model of chronic ischaemia.39

Axonal sparing

NAA’s almost exclusive confinement to neurons renders its 1H-MRS signal a marker of their integrity.6 Since axonal loss is thought to proceed mainly by Wallerian or retrograde degeneration from lesions,25 40 conflicting reports of normal34 35 and low11 12 13 41 NAA may reflect different disease course, duration and subjective VOI placement. The last issue was addressed with whole-brain NAA quantification showing an average of 22% decline already in CIS.42 However, since over 60% of that NAA signal originates from cortical GM (known to be affected early in MS11 41 43), that finding is not necessarily in conflict with the NAA levels reported here which reflect mainly WM.

Despite contradictory reports, the concept of eventual NAA decreases,9 axonal loss7 and atrophy2 in NAWM is now generally accepted. We report on glial abnormalities early in the disease course, in the absence of, that is, prior to overt axonal injury, reflected in patients’ NAA level and Tf not significantly different from the controls’. These findings support the notion that glial processes including astrogliosis precede axonal damage, which in turn precedes tissue atrophy—the final non-specific end-point of the various pathologies of this disease.

Based on these findings we conjecture that some axonal damage may arise from astrogliosis and/or inflammation. Indeed, disability’s modest correlation with lesion load7 and histological evidence of axonal loss in the absence of apparent demyelination suggest that not all damage is due to lesions.44 Indeed, both biopsy and autopsy show that axonal damage correlates better with inflammatory infiltrates than with demyelination.45 46 The concept of astrogliosis-driven damage is supported by the EAE animal model of MS that shows astrocytic hypertrophy that is spatially and temporally related to axonal damage already in the earliest preinflammatory stage.47

Caveats

The sensitivity to metabolic changes reported here comes at a price—regional differences, especially between lesions and NAWM, are averaged out. Although the original localised information is still available (see figs 1, 2), the proposed strategy reflects a postprocessing choice of maximum sensitivity at the expense of localisation. Monitoring the global effects of diffuse pathogenesis may make this trade-off attractive, especially for the much weaker mI signal.

A cautionary note raised by this study underscores the importance of absolute quantification of 1H-MRS in MS versus use of metabolite ratios to assumed stable Cr levels. This and other reports35 that show increased Cr and unchanged NAA suggest that lower NAA/Cr (interpreted in the past as “axonal dysfunction”) may in fact be due primarily to elevated Cr. Similarly, the 15% Cho increase observed in patients here and in other studies will attenuate to ∼5% in a Cho/Cr due to the concomitant Cr increase, reducing the ratio’s discriminatory power.

Conclusion

Elevated mI, Cr (glial proliferation) and Cho (membrane turnover) but normal NAA levels and tissue fraction (absence of axonal injury and atrophy) in RR MS patients suggest ongoing diffuse inflammation, demyelination and (or) remyelination, as well as astrogliosis early in the disease course, even during clinical remission and despite immunomodulatory treatment. Consequently, mI, Cho and Cr may serve as earlier MRS markers of disease activity than the traditional NAA decline and atrophy that represent in order a later and the final stage of MS pathogenesis.

Acknowledgments

We thank AA Maudsley and BJ Soher for the use of their SITools-FITT spectral modelling software.

REFERENCES

Footnotes

  • Funding This work was supported by NIH Grants NS050520, NS29029 and EB01015.

  • Competing interests None.

  • Ethics approval Ethics approval was provided by New York University School of Medicine Institutional Review Board.

  • Patient consent Obtained.

  • Provenance and Peer review Not commissioned; externally peer reviewed.