Christina J Azevedo1, Eve Overton1, Sankalpa Khadka1, Jessica Buckley1, Shuang Liu1, Mehul Sampat1, Orhun Kantarci1, Christine Lebrun Frenay1, Aksel Siva1, Darin T Okuda1, Daniel Pelletier1. 1. Department of Neurology (C.J.A., E.O., S.K., J.B., S.L., M.S., D.P.) and Department of Diagnostic Radiology (D.P.), Yale University, New Haven, CT; Department of Neurology (O.K.), Mayo Clinic College of Medicine, Rochester, MN; Service de Neurologie (C.L.F.), Hospital Pasteur, Nice, France; Department of Neurology (A.S.), Cerrahpasa School of Medicine, University of Istanbul, Turkey; and Department of Neurology (D.T.O.), University of Texas Southwestern Medical Center, Dallas, TX.
Abstract
OBJECTIVE: Increasing evidence indicates that the thalamus may be a location of early neurodegeneration in multiple sclerosis (MS). Our objective was to identify the presence of gray matter volume loss and thinning in patients with radiologically isolated syndrome (RIS). METHODS: Sixty-three participants were included in this case-control study. Twenty-one patients with RIS were age- and sex-matched to 42 healthy controls in a 1:2 ratio. All participants underwent brain MRIs on a single 3T scanner. After lesion segmentation and inpainting, 1 mm(3)-isometric T1-weighted images were submitted to FreeSurfer (v5.2). Normalized cortical and deep gray matter volumes were compared between patients with RIS and controls using t tests, and thalamic volumes were correlated with white matter lesion volumes using Pearson correlation. Exploratory cortical thickness maps were created. RESULTS: Although traditional normalized total gray and white matter volumes were not statistically different between patients with RIS and controls, normalized left (0.0046 ± 0.0005 vs 0.0049 ± 0.0004, p = 0.006), right (0.0045 ± 0.0005 vs 0.0048 ± 0.0004, p = 0.008), and mean (0.0045 ± 0.0005 vs 0.0049 ± 0.0004, p = 0.004) thalamic volumes were significantly lower in patients with RIS (n = 21, mean age 41.9 ± 12.7 years) than in controls (n = 42, mean age 41.4 ± 11.2 years). Thalamic volumes correlated modestly with white matter lesion volumes (range: r = -0.35 to -0.47). CONCLUSION: Our data provide novel evidence of thalamic atrophy in RIS and are consistent with previous reports in early MS stages. Thalamic volume loss is present early in CNS demyelinating disease and should be further investigated as a metric associated with neurodegeneration.
OBJECTIVE: Increasing evidence indicates that the thalamus may be a location of early neurodegeneration in multiple sclerosis (MS). Our objective was to identify the presence of gray matter volume loss and thinning in patients with radiologically isolated syndrome (RIS). METHODS: Sixty-three participants were included in this case-control study. Twenty-one patients with RIS were age- and sex-matched to 42 healthy controls in a 1:2 ratio. All participants underwent brain MRIs on a single 3T scanner. After lesion segmentation and inpainting, 1 mm(3)-isometric T1-weighted images were submitted to FreeSurfer (v5.2). Normalized cortical and deep gray matter volumes were compared between patients with RIS and controls using t tests, and thalamic volumes were correlated with white matter lesion volumes using Pearson correlation. Exploratory cortical thickness maps were created. RESULTS: Although traditional normalized total gray and white matter volumes were not statistically different between patients with RIS and controls, normalized left (0.0046 ± 0.0005 vs 0.0049 ± 0.0004, p = 0.006), right (0.0045 ± 0.0005 vs 0.0048 ± 0.0004, p = 0.008), and mean (0.0045 ± 0.0005 vs 0.0049 ± 0.0004, p = 0.004) thalamic volumes were significantly lower in patients with RIS (n = 21, mean age 41.9 ± 12.7 years) than in controls (n = 42, mean age 41.4 ± 11.2 years). Thalamic volumes correlated modestly with white matter lesion volumes (range: r = -0.35 to -0.47). CONCLUSION: Our data provide novel evidence of thalamic atrophy in RIS and are consistent with previous reports in early MS stages. Thalamic volume loss is present early in CNS demyelinating disease and should be further investigated as a metric associated with neurodegeneration.
The thalamus is a major relay center with extensive cortical and subcortical anatomic connections. Increasing evidence suggests that the thalamus may be a location of early neurodegeneration in multiple sclerosis (MS). Thalamic atrophy has been reported in early relapsing-remitting MS (RRMS),[1] clinically isolated syndrome (CIS) at presentation,[2] and pediatric MS.[3] Thalamic volume loss is likely due in part to disconnection caused by white matter lesions; using diffusion tensor imaging tractography, a relationship between thalamic volume and white matter lesion volume within thalamocortical projections has been observed.[4] In addition to the thalamus, cortical gray matter may exhibit early volume loss or thinning. Regional cortical thinning has been identified as early as CIS.[5]Radiologically isolated syndrome (RIS) is defined as incidentally identified white matter anomalies within the CNS that are suggestive of MS in healthy individuals who lack clinical symptomatology associated with CNS demyelination, after exclusion of other possible etiologies.[6] Approximately 30% of patients with RIS will have a sentinel clinical demyelinating event within 5 years of the identification of their abnormal MRI[6,7]; thus, for many patients, RIS represents the earliest visible manifestation of demyelinating disease, or a preclinical stage of MS.The objective of this work was to investigate the presence of gray matter volume loss and thinning in patients with RIS using the FreeSurfer image analysis suite, an unbiased segmentation method. In addition to very robust and fully automated subcortical anatomic segmentation, FreeSurfer is capable of detecting submillimeter differences in cortical thickness between groups.[8]
METHODS
Participants.
Healthy individuals not exposed to a disease-modifying therapy fulfilling 2009 RIS criteria[6] were retrospectively identified at the University of California, San Francisco (UCSF) MS Center. For each participant with RIS, the first 2 available age- and sex-matched healthy controls were selected from a larger dataset (i.e., 1:2 ratio).
Standard protocol approvals, registrations, and patient consents.
The protocol was approved by the Committee on Human Research at UCSF, and informed consent was obtained from all participants.
MRI acquisition.
All study subjects were scanned on a single 3T GE Signa scanner using a standardized protocol. The protocol included an inversion recovery spoiled gradient echo T1-weighted isotropic volumetric sequence (3D-IRSPGR, 1 × 1 × 1 mm3, 180 slices, echo time/repetition time/inversion time = 2/7/400 ms, flip angle = 15°, 256 × 256 × 180 matrix, 240 × 240 × 180 mm3 field of view, number of excitations = 1).
MRI postprocessing.
The 3D-IRSPGR served as the input image into the postprocessing pipeline. T1 lesion masks were derived from manual segmentation of T1-visible white matter lesions, as described previously.[6] None of the patients with RIS had visible thalamic lesions. MS lesions were inpainted[9] on the 3D-IRSPGR prior to submitting the image to FreeSurfer. Inpainting was critical in order to avoid voxel misclassification due to juxtacortical MS lesions, which are of similar intensity to cortical gray matter on gradient echo images and can result in parcellation errors in FreeSurfer's output. The inpainted 3D-IRSPGRs were submitted to FreeSurfer (v5.2; Desikan-Killiany atlas for surface parcellation and probabilistic atlas for subcortical segmentation). FreeSurfer output was reviewed by an experienced MRI reader/postprocessor (C.J.A.) and manually corrected using methods within FreeSurfer such as control points, white matter edits, and dura edits. Cases were rerun as needed until no further improvement was achieved by additional editing. On average, 3 runs per case were performed, but some cases required up to 7 runs. Particular attention was paid to error-prone regions such as the mesial temporal lobes.
Statistical analyses.
Volumes (i.e., total gray matter, total white matter, deep gray nuclei) were extracted directly from FreeSurfer's automated segmentation and normalized for head size using the total intracranial volume (also taken from FreeSurfer). These values (unitless) were compared between patients with RIS and controls using t tests. Among patients with RIS, thalamic volumes were correlated with log-transformed white matter lesion volumes using Pearson correlation. Statistical analyses were performed using SAS (v9.2). All hypothesis tests in this portion of the analysis were performed at a nominal α = 0.05 significance level without adjustment for multiple comparisons. Adjustment for multiple comparisons is generally performed when several independent hypothesis tests are conducted on one population. In this case, volumes for a given participant are likely to be correlated with each other and are thus not independent. As such, a correction would be too strict. Scientific judgment can be used to indicate where caution is warranted in interpreting results with p <0.05.An exploratory, group-level, vertex-wise analysis of cortical thickness was performed using QDEC, a tool within FreeSurfer. Briefly, as part of its cortical parcellation, FreeSurfer reconstructs each individual's cortical surface using a meshwork of ∼150,000 triangles. Each point of intersection between triangles is known as a vertex, and each vertex has an associated XYZ coordinate. QDEC creates an average cortical thickness map per group of participants and then generates a significant difference map by performing a generalized linear model at each vertex. Although we matched on age and sex, the final QDEC result was adjusted for both of these variables to avoid the possibility of residual confounding. QDEC's output was corrected for multiple comparisons using a Monte Carlo simulation procedure.
RESULTS
Participant demographics.
Baseline demographics for all participants are given in table 1.
Table 1
Participant demographics
Participant demographics
Comparison of volumes between patients with RIS and controls.
Normalized total parenchymal (0.782 ± 0.046 vs 0.798 ± 0.035, p = 0.12) and total white matter (0.33 ± 0.02 vs 0.33 ± 0.02, p = 0.6) volumes were not statistically different between patients with RIS and controls. However, a trend toward lower normalized total gray matter volumes was seen in patients with RIS compared with controls (p = 0.06, table 2). When total gray matter was divided into cortical and subcortical gray matter, normalized cortical gray matter volume did not differ significantly between patients with RIS and controls (p = 0.13), but normalized subcortical gray matter volume was significantly lower in patients with RIS than in controls (p = 0.04). Among all subcortical gray matter structures, the thalamus was the main driver of this difference. Normalized left (p = 0.006), right (p = 0.008), and mean (p = 0.004) thalamic volumes were significantly lower in patients with RIS than in controls.
Table 2
Normalized gray matter regions in patients with radiologically isolated syndrome (RIS) and controls
Normalized gray matter regions in patients with radiologically isolated syndrome (RIS) and controls
Correlation between thalamic volumes and white matter lesion volumes.
Among patients with RIS, mean thalamic (r = −0.42, p = 0.06), left thalamic (r = −0.35, p = 0.13), and right thalamic (r = −0.47, p = 0.03) volumes were modestly correlated with white matter lesion volumes. The correlations were in the expected direction and generally exhibited a trend toward statistical significance. The asymmetry in the strength of the correlation for right and left thalami was not expected a priori and may have been due to chance.
Cortical thinning.
The final result from QDEC analysis is shown in the figure. An area measuring approximately 3 cm in the right superior and inferior parietal gyri (Desikan-Killiany labels) was significantly thinner in patients with RIS than in controls. It was present at Monte Carlo thresholds of p <0.05 and p <0.01, and was stable up to a threshold of p <0.005.
Figure
Vertex-wise significant difference map (QDEC) of cortical thickness between patients with RIS and controls
An area measuring approximately 3 cm in the right superior and inferior parietal gyri (Desikan-Killiany labels, in red) was significantly thinner in patients with radiologically isolated syndrome (RIS) than in controls at a Monte Carlo threshold of p <0.005 and adjusted for age and sex.
Vertex-wise significant difference map (QDEC) of cortical thickness between patients with RIS and controls
An area measuring approximately 3 cm in the right superior and inferior parietal gyri (Desikan-Killiany labels, in red) was significantly thinner in patients with radiologically isolated syndrome (RIS) than in controls at a Monte Carlo threshold of p <0.005 and adjusted for age and sex.
DISCUSSION
These data provide novel in vivo evidence to suggest that tissue loss may precede clinical manifestations of CNS demyelination. These findings add to previous work reporting that thalamic atrophy is present in early stages of CNS demyelinating disease, including early RRMS,[1] CIS,[2] and pediatric MS.[3]Thalamic atrophy is likely due to several converging aspects of MS disease pathology. Although we did not see thalamic lesions in our RIS population, thalamic demyelinating lesions[10] and iron deposition[11] likely contribute to thalamic volume loss. In addition, neuronal density has been shown to be reduced in nonlesional thalamic tissue,[12] which may be driven by disconnection from white matter lesions in the various tracts that project from and into the thalamus[4] (e.g., thalamocortical, opticothalamic, and spinothalamic projections). Further studies at very high field strengths examining the interplay of these factors, including the contribution of thalamic demyelinating lesions and iron deposition to thalamic volume loss, will be of interest.Although others have reported decreased whole brain and cortical gray matter volumes in patients with RIS compared with healthy controls,[13] we did not see such differences in our data. Our findings are consistent with a CIS cohort[2] in which there was no difference in normalized whole brain or gray matter volumes but there was a reduction in thalamic volume using voxel-based morphometry. To our knowledge, there are no other published studies examining cortical thickness in patients with RIS. We report our finding of a ∼3 cm area of thinning in the right superior and inferior parietal gyri but regard it as a preliminary finding without a specific rationale for why these particular gyri are affected in our dataset. In a group of patients with CIS, this region was not significantly different from controls using different FreeSurfer methods and input images.[5] Further studies examining whole brain and regional gray matter volumes and cortical thickness in early demyelinating disease will be of interest in order to replicate these findings.Our study has certain limitations, including a small sample size and a single-center design. In addition, cortical lesions are of great interest in the field currently but are generally not visible on conventional MRIs, including the IRSPGRs we used; as such, this work is not able to address the impact of cortical lesions on cortical thinning. Despite these limitations, high-resolution images and a careful MRI postprocessing approach were used, and the results are felt to be robust.The reflection of multiple aspects of MS pathology, the presence of thalamic volume loss in early disease, the relative ease of measuring thalamic volume by various methods, and the observed correlation of thalamic volume with cognition[14] make thalamic volume an attractive potential MRI clinical trial metric associated with the neurodegenerative aspect of MS and potentially RIS. Further studies are needed to investigate the temporal evolution and the effect of disease-modifying therapies on thalamic volume loss. In the context of RIS, larger longitudinal studies are currently under way. Although there may be other areas of gray matter volume loss or thinning in RIS, thalamic volume will be of particular interest in these studies and will be examined as a potential predictor of conversion to CIS/MS.
Authors: D T Okuda; E M Mowry; A Beheshtian; E Waubant; S E Baranzini; D S Goodin; S L Hauser; D Pelletier Journal: Neurology Date: 2008-12-10 Impact factor: 9.910
Authors: Alberto Cifelli; Marzena Arridge; Peter Jezzard; Margaret M Esiri; Jacqueline Palace; Paul M Matthews Journal: Ann Neurol Date: 2002-11 Impact factor: 10.422
Authors: B Aubert-Broche; V Fonov; R Ghassemi; S Narayanan; D L Arnold; B Banwell; J G Sled; D L Collins Journal: Neuroimage Date: 2011-03-14 Impact factor: 6.556
Authors: Robert Zivadinov; Mari Heininen-Brown; Claudiu V Schirda; Guy U Poloni; Niels Bergsland; Christopher R Magnano; Jacqueline Durfee; Cheryl Kennedy; Ellen Carl; Jesper Hagemeier; Ralph H B Benedict; Bianca Weinstock-Guttman; Michael G Dwyer Journal: Neuroimage Date: 2011-07-27 Impact factor: 6.556
Authors: M K Houtchens; R H B Benedict; R Killiany; J Sharma; Z Jaisani; B Singh; B Weinstock-Guttman; C R G Guttmann; R Bakshi Journal: Neurology Date: 2007-09-18 Impact factor: 9.910
Authors: R G Henry; M Shieh; D T Okuda; A Evangelista; M L Gorno-Tempini; D Pelletier Journal: J Neurol Neurosurg Psychiatry Date: 2008-05-09 Impact factor: 10.154
Authors: Darin T Okuda; Aksel Siva; Orhun Kantarci; Matilde Inglese; Ilana Katz; Melih Tutuncu; B Mark Keegan; Stacy Donlon; Le H Hua; Angela Vidal-Jordana; Xavier Montalban; Alex Rovira; Mar Tintoré; Maria Pia Amato; Bruno Brochet; Jérôme de Seze; David Brassat; Patrick Vermersch; Nicola De Stefano; Maria Pia Sormani; Daniel Pelletier; Christine Lebrun Journal: PLoS One Date: 2014-03-05 Impact factor: 3.240
Authors: Burcu Zeydan; Val J Lowe; Christopher G Schwarz; Scott A Przybelski; Nirubol Tosakulwong; Samantha M Zuk; Matthew L Senjem; Jeffrey L Gunter; Rosebud O Roberts; Michelle M Mielke; Eduardo E Benarroch; Moses Rodriguez; Mary M Machulda; Timothy G Lesnick; David S Knopman; Ronald C Petersen; Clifford R Jack; Kejal Kantarci; Orhun H Kantarci Journal: Mult Scler Date: 2017-05-05 Impact factor: 6.312