Synaptic Protein Loss in Extracellular Vesicles Reflects Brain and Retinal Atrophy in People With Multiple Sclerosis.

OBJECTIVES: To assess whether the rate of change in synaptic proteins isolated from neuronally enriched extracellular vesicles (NEVs) is associated with brain and retinal atrophy in people with multiple sclerosis (MS). METHODS: People with MS were followed with serial blood draws, MRI (MRI), and optical coherence tomography (OCT) scans. NEVs were immunocaptured from plasma, and synaptopodin and synaptophysin proteins were measured using ELISA. Subject-specific rates of change in synaptic proteins, as well as brain and retinal atrophy, were determined and correlated. RESULTS: A total of 50 people with MS were included, 46 of whom had MRI and 45 had OCT serially. The rate of change in NEV synaptopodin was associated with whole brain (rho = 0.31; p = 0.04), cortical gray matter (rho = 0.34; p = 0.03), peripapillary retinal nerve fiber layer (rho = 0.37; p = 0.01), and ganglion cell/inner plexiform layer (rho = 0.41; p = 0.006) atrophy. The rate of change in NEV synaptophysin was also correlated with whole brain (rho = 0.31; p = 0.04) and cortical gray matter (rho = 0.31; p = 0.049) atrophy. DISCUSSION: NEV-derived synaptic proteins likely reflect neurodegeneration and may provide additional circulating biomarkers for disease progression in MS.

Development of Medical Imaging Data Standardization for Imaging-Based Observational Research: OMOP Common Data Model Extension.

The rapid growth of artificial intelligence (AI) and deep learning techniques require access to large inter-institutional cohorts of data to enable the development of robust models, e.g., targeting the identification of disease biomarkers and quantifying disease progression and treatment efficacy. The Observational Medical Outcomes Partnership Common Data Model (OMOP CDM) has been designed to accommodate a harmonized representation of observational healthcare data. This study proposes the Medical Imaging CDM (MI-CDM) extension, adding two new tables and two vocabularies to the OMOP CDM to address the structural and semantic requirements to support imaging research. The tables provide the capabilities of linking DICOM data sources as well as tracking the provenance of imaging features derived from those images. The implementation of the extension enables phenotype definitions using imaging features and expanding standardized computable imaging biomarkers. This proposal offers a comprehensive and unified approach for conducting imaging research and outcome studies utilizing imaging features.

Image harmonization improves consistency of intra-rater delineations of MS lesions in heterogeneous MRI.

Clinical magnetic resonance images (MRIs) lack a standard intensity scale due to differences in scanner hardware and the pulse sequences used to acquire the images. When MRIs are used for quantification, as in the evaluation of white matter lesions (WMLs) in multiple sclerosis, this lack of intensity standardization becomes a critical problem affecting both the staging and tracking of the disease and its treatment. This paper presents a study of harmonization on WML segmentation consistency, which is evaluated using an object detection classification scheme that incorporates manual delineations from both the original and harmonized MRIs. A cohort of ten people scanned on two different imaging platforms was studied. An expert rater, blinded to the image source, manually delineated WMLs on images from both scanners before and after harmonization. It was found that there is closer agreement in both global and per-lesion WML volume and spatial distribution after harmonization, demonstrating the importance of image harmonization prior to the creation of manual delineations. These results could lead to better truth models in both the development and evaluation of automated lesion segmentation algorithms.

Bile acid metabolites predict multiple sclerosis progression and supplementation is safe in progressive disease.

Background: Bile acid metabolism is altered in multiple sclerosis (MS) and tauroursodeoxycholic acid (TUDCA) supplementation ameliorated disease in mouse models of MS. Methods: Global metabolomics was performed in an observational cohort of people with MS followed by pathway analysis to examine relationships between baseline metabolite levels and subsequent brain and retinal atrophy. A double-blind, placebo-controlled trial, was completed in people with progressive MS (PMS), randomized to receive either TUDCA (2g daily) or placebo for 16 weeks. Participants were followed with serial clinical and laboratory assessments. Primary outcomes were safety and tolerability of TUDCA, and exploratory outcomes included changes in clinical, laboratory and gut microbiome parameters. Results: In the observational cohort, higher primary bile acid levels at baseline predicted slower whole brain, brain substructure and specific retinal layer atrophy. In the clinical trial, 47 participants were included in our analyses (21 in placebo arm, 26 in TUDCA arm). Adverse events did not significantly differ between arms (p=0.77). The TUDCA arm demonstrated increased serum levels of multiple bile acids. No significant differences were noted in clinical or fluid biomarker outcomes. Central memory CD4+ and Th1/17 cells decreased, while CD4+ naïve cells increased in the TUDCA arm compared to placebo. Changes in the composition and function of gut microbiota were also noted in the TUDCA arm compared to placebo. Conclusion: Bile acid metabolism in MS is linked with brain and retinal atrophy. TUDCA supplementation in PMS is safe, tolerable and has measurable biological effects that warrant further evaluation in larger trials with a longer treatment duration.

HACA3: A unified approach for multi-site MR image harmonization.

The lack of standardization and consistency of acquisition is a prominent issue in magnetic resonance (MR) imaging. This often causes undesired contrast variations in the acquired images due to differences in hardware and acquisition parameters. In recent years, image synthesis-based MR harmonization with disentanglement has been proposed to compensate for the undesired contrast variations. The general idea is to disentangle anatomy and contrast information from MR images to achieve cross-site harmonization. Despite the success of existing methods, we argue that major improvements can be made from three aspects. First, most existing methods are built upon the assumption that multi-contrast MR images of the same subject share the same anatomy. This assumption is questionable, since different MR contrasts are specialized to highlight different anatomical features. Second, these methods often require a fixed set of MR contrasts for training (e.g., both T1-weighted and T2-weighted images), limiting their applicability. Lastly, existing methods are generally sensitive to imaging artifacts. In this paper, we present Harmonization with Attention-based Contrast, Anatomy, and Artifact Awareness (HACA3), a novel approach to address these three issues. HACA3 incorporates an anatomy fusion module that accounts for the inherent anatomical differences between MR contrasts. Furthermore, HACA3 can be trained and applied to any combination of MR contrasts and is robust to imaging artifacts. HACA3 is developed and evaluated on diverse MR datasets acquired from 21 sites with varying field strengths, scanner platforms, and acquisition protocols. Experiments show that HACA3 achieves state-of-the-art harmonization performance under multiple image quality metrics. We also demonstrate the versatility and potential clinical impact of HACA3 on downstream tasks including white matter lesion segmentation for people with multiple sclerosis and longitudinal volumetric analyses for normal aging subjects. Code is available at https://github.com/lianruizuo/haca3.

Intersite brain MRI volumetric biases persist even in a harmonized multisubject study of multiple sclerosis.

BACKGROUND AND PURPOSE: Multicenter study designs involving a variety of MRI scanners have become increasingly common. However, these present the issue of biases in image-based measures due to scanner or site differences. To assess these biases, we imaged 11 volunteers with multiple sclerosis (MS) with scan and rescan data at four sites. METHODS: Images were acquired on Siemens or Philips scanners at 3 Tesla. Automated white matter lesion detection and whole-brain, gray and white matter, and thalamic volumetry were performed, as well as expert manual delineations of T1 magnetization-prepared rapid acquisition gradient echo and T2 fluid-attenuated inversion recovery lesions. Random-effect and permutation-based nonparametric modeling was performed to assess differences in estimated volumes within and across sites. RESULTS: Random-effect modeling demonstrated model assumption violations for most comparisons of interest. Nonparametric modeling indicated that site explained >50% of the variation for most estimated volumes. This expanded to >75% when data from both Siemens and Philips scanners were included. Permutation tests revealed significant differences between average inter- and intrasite differences in most estimated brain volumes (P < .05). The automatic activation of spine coil elements during some acquisitions resulted in a shading artifact in these images. Permutation tests revealed significant differences between thalamic volume measurements from acquisitions with and without this artifact. CONCLUSION: Differences in brain volumetry persisted across MR scanners despite protocol harmonization. These differences were not well explained by variance component modeling; however, statistical innovations for mitigating intersite differences show promise in reducing biases in multicenter studies of MS.

Deep gray matter substructure volumes and depressive symptoms in a large multiple sclerosis cohort.

BACKGROUND: Consistent findings on underlying brain features or specific structural atrophy patterns contributing to depression in multiple sclerosis (MS) are limited. OBJECTIVE: To investigate how deep gray matter (DGM) features predict depressive symptom trajectories in MS patients. METHODS: We used data from the MS Partners Advancing Technology and Health Solutions (MS PATHS) network in which standardized patient information and outcomes are collected. We performed whole-brain segmentation using SLANT-CRUISE. We assessed if DGM structures were associated with elevated depressive symptoms over follow-up and with depressive symptom phenotypes. RESULTS: We included 3844 participants (average age: 46.05 ± 11.83 years; 72.7% female) of whom 1905 (49.5%) experienced ⩾1 periods of elevated depressive symptoms over 2.6 ± 0.9 years mean follow-up. Higher caudate, putamen, accumbens, ventral diencephalon, thalamus, and amygdala volumes were associated with lower odds of elevated depressive symptoms over follow-up (odds ratio (OR) range per 1 SD (standard deviation) increase in volume: 0.88-0.94). For example, a 1 SD increase in accumbens or caudate volume was associated with 12% or 10% respective lower odds of having a period of elevated depressive symptoms over follow-up (for accumbens: OR: 0.88; 95% confidence interval (CI): 0.83-0.93; p < 0.001; for caudate: OR: 0.90; 95% CI: 0.85-0.96; p = 0.003). CONCLUSION: Lower DGM volumes were associated with depressive symptom trajectories in MS.

Trans-Synaptic Degeneration Following Acute Optic Neuritis in Multiple Sclerosis.

OBJECTIVE: To explore longitudinal changes in brain volumetric measures and retinal layer thicknesses following acute optic neuritis (AON) in people with multiple sclerosis (PwMS), to investigate the process of trans-synaptic degeneration, and determine its clinical relevance. METHODS: PwMS were recruited within 40 days of AON onset (n = 49), and underwent baseline retinal optical coherence tomography and brain magnetic resonance imaging followed by longitudinal tracking for up to 5 years. A comparator cohort of PwMS without a recent episode of AON were similarly tracked (n = 73). Mixed-effects linear regression models were used. RESULTS: Accelerated atrophy of the occipital gray matter (GM), calcarine GM, and thalamus was seen in the AON cohort, as compared with the non-AON cohort (-0.76% vs -0.22% per year [p = 0.01] for occipital GM, -1.83% vs -0.32% per year [p = 0.008] for calcarine GM, -1.17% vs -0.67% per year [p = 0.02] for thalamus), whereas rates of whole-brain, cortical GM, non-occipital cortical GM atrophy, and T2 lesion accumulation did not differ significantly between the cohorts. In the AON cohort, greater AON-induced reduction in ganglion cell+inner plexiform layer thickness over the first year was associated with faster rates of whole-brain (r = 0.32, p = 0.04), white matter (r = 0.32, p = 0.04), and thalamic (r = 0.36, p = 0.02) atrophy over the study period. Significant relationships were identified between faster atrophy of the subcortical GM and thalamus, with worse visual function outcomes after AON. INTERPRETATION: These results provide in-vivo evidence for anterograde trans-synaptic degeneration following AON in PwMS, and suggest that trans-synaptic degeneration may be related to clinically-relevant visual outcomes. ANN NEUROL 2023;93:76-87.

Bi-directional Synthesis of Pre- and Post-contrast MRI via Guided Feature Disentanglement.

Magnetic resonance imaging (MRI) with gadolinium contrast is widely used for tissue enhancement and better identification of active lesions and tumors. Recent studies have shown that gadolinium deposition can accumulate in tissues including the brain, which raises safety concerns. Prior works have tried to synthesize post-contrast T1-weighted MRIs from pre-contrast MRIs to avoid the use of gadolinium. However, contrast and image representations are often entangled during the synthesis process, resulting in synthetic post-contrast MRIs with undesirable contrast enhancements. Moreover, the synthesis of pre-contrast MRIs from post-contrast MRIs which can be useful for volumetric analysis is rarely investigated in the literature. To tackle pre- and post- contrast MRI synthesis, we propose a BI-directional Contrast Enhancement Prediction and Synthesis (BICEPS) network that enables disentanglement of contrast and image representations via a bi-directional image-to-image translation(I2I)model. Our proposed model can perform both pre-to-post and post-to-pre contrast synthesis, and provides an interpretable synthesis process by predicting contrast enhancement maps from the learned contrast embedding. Extensive experiments on a multiple sclerosis dataset demonstrate the feasibility of applying our bidirectional synthesis and show that BICEPS outperforms current methods.

Unsupervised MR harmonization by learning disentangled representations using information bottleneck theory.

In magnetic resonance (MR) imaging, a lack of standardization in acquisition often causes pulse sequence-based contrast variations in MR images from site to site, which impedes consistent measurements in automatic analyses. In this paper, we propose an unsupervised MR image harmonization approach, CALAMITI (Contrast Anatomy Learning and Analysis for MR Intensity Translation and Integration), which aims to alleviate contrast variations in multi-site MR imaging. Designed using information bottleneck theory, CALAMITI learns a globally disentangled latent space containing both anatomical and contrast information, which permits harmonization. In contrast to supervised harmonization methods, our approach does not need a sample population to be imaged across sites. Unlike traditional unsupervised harmonization approaches which often suffer from geometry shifts, CALAMITI better preserves anatomy by design. The proposed method is also able to adapt to a new testing site with a straightforward fine-tuning process. Experiments on MR images acquired from ten sites show that CALAMITI achieves superior performance compared with other harmonization approaches.

Autoencoder based self-supervised test-time adaptation for medical image analysis.

Deep neural networks have been successfully applied to medical image analysis tasks like segmentation and synthesis. However, even if a network is trained on a large dataset from the source domain, its performance on unseen test domains is not guaranteed. The performance drop on data obtained differently from the network's training data is a major problem (known as domain shift) in deploying deep learning in clinical practice. Existing work focuses on retraining the model with data from the test domain, or harmonizing the test domain's data to the network training data. A common practice is to distribute a carefully-trained model to multiple users (e.g., clinical centers), and then each user uses the model to process their own data, which may have a domain shift (e.g., varying imaging parameters and machines). However, the lack of availability of the source training data and the cost of training a new model often prevents the use of known methods to solve user-specific domain shifts. Here, we ask whether we can design a model that, once distributed to users, can quickly adapt itself to each new site without expensive retraining or access to the source training data? In this paper, we propose a model that can adapt based on a single test subject during inference. The model consists of three parts, which are all neural networks: a task model (T) which performs the image analysis task like segmentation; a set of autoencoders (AEs); and a set of adaptors (As). The task model and autoencoders are trained on the source dataset and can be computationally expensive. In the deployment stage, the adaptors are trained to transform the test image and its features to minimize the domain shift as measured by the autoencoders' reconstruction loss. Only the adaptors are optimized during the testing stage with a single test subject thus is computationally efficient. The method was validated on both retinal optical coherence tomography (OCT) image segmentation and magnetic resonance imaging (MRI) T1-weighted to T2-weighted image synthesis. Our method, with its short optimization time for the adaptors (10 iterations on a single test subject) and its additional required disk space for the autoencoders (around 15 MB), can achieve significant performance improvement. Our code is publicly available at: https://github.com/YufanHe/self-domain-adapted-network.

Joint Image and Label Self-Super-Resolution.

We propose a method to jointly super-resolve an anisotropic image volume along with its corresponding voxel labels without external training data. Our method is inspired by internally trained superresolution, or self-super-resolution (SSR) techniques that target anisotropic, low-resolution (LR) magnetic resonance (MR) images. While resulting images from such methods are quite useful, their corresponding LR labels-derived from either automatic algorithms or human raters-are no longer in correspondence with the super-resolved volume. To address this, we develop an SSR deep network that takes both an anisotropic LR MR image and its corresponding LR labels as input and produces both a super-resolved MR image and its super-resolved labels as output. We evaluated our method with 50 T 1-weighted brain MR images 4× down-sampled with 10 automatically generated labels. In comparison to other methods, our method had superior Dice across all labels and competitive metrics on the MR image. Our approach is the first reported method for SSR of paired anisotropic image and label volumes.

Frequency and phase correction of J-difference edited MR spectra using deep learning.

PURPOSE: To investigate whether a deep learning-based (DL) approach can be used for frequency-and-phase correction (FPC) of MEGA-edited MRS data. METHODS: Two neural networks (1 for frequency, 1 for phase) consisting of fully connected layers were trained and validated using simulated MEGA-edited MRS data. This DL-FPC was subsequently tested and compared to a conventional approach (spectral registration [SR]) and to a model-based SR implementation (mSR) using in vivo MEGA-edited MRS datasets. Additional artificial offsets were added to these datasets to further investigate performance. RESULTS: The validation showed that DL-based FPC was capable of correcting within 0.03 Hz of frequency and 0.4°of phase offset for unseen simulated data. DL-based FPC performed similarly to SR for the unmanipulated in vivo test datasets. When additional offsets were added to these datasets, the networks still performed well. However, although SR accurately corrected for smaller offsets, it often failed for larger offsets. The mSR algorithm performed well for larger offsets, which was because the model was generated from the in vivo datasets. In addition, the computation times were much shorter using DL-based FPC or mSR compared to SR for heavily distorted spectra. CONCLUSION: These results represent a proof of principle for the use of DL for preprocessing MRS data.

SMORE: A Self-Supervised Anti-Aliasing and Super-Resolution Algorithm for MRI Using Deep Learning.

High resolution magnetic resonance (MR) images are desired in many clinical and research applications. Acquiring such images with high signal-to-noise (SNR), however, can require a long scan duration, which is difficult for patient comfort, is more costly, and makes the images susceptible to motion artifacts. A very common practical compromise for both 2D and 3D MR imaging protocols is to acquire volumetric MR images with high in-plane resolution, but lower through-plane resolution. In addition to having poor resolution in one orientation, 2D MRI acquisitions will also have aliasing artifacts, which further degrade the appearance of these images. This paper presents an approach SMORE1 based on convolutional neural networks (CNNs) that restores image quality by improving resolution and reducing aliasing in MR images.2 This approach is self-supervised, which requires no external training data because the high-resolution and low-resolution data that are present in the image itself are used for training. For 3D MRI, the method consists of only one self-supervised super-resolution (SSR) deep CNN that is trained from the volumetric image data. For 2D MRI, there is a self-supervised anti-aliasing (SAA) deep CNN that precedes the SSR CNN, also trained from the volumetric image data. Both methods were evaluated on a broad collection of MR data, including filtered and downsampled images so that quantitative metrics could be computed and compared, and actual acquired low resolution images for which visual and sharpness measures could be computed and compared. The super-resolution method is shown to be visually and quantitatively superior to previously reported methods.

Experimental design and sample size considerations in longitudinal magnetic resonance imaging-based biomarker detection for multiple sclerosis.

Several modeling approaches have been developed to quantify differences in multiple sclerosis lesion evolution on magnetic resonance imaging to identify the effect of treatment on disease progression. These studies have limited clinical applicability due to onerous scan frequency and lengthy study duration. Efficient methods are needed to reduce the required sample size, study duration, and sampling frequency in longitudinal magnetic resonance imaging studies. We develop a data-driven approach to identify parameters of study design for evaluation of longitudinal magnetic resonance imaging biomarkers of multiple sclerosis lesion evolution. Our design strategies are considerably shorter than those described in previous studies, thus having the potential to lower costs of clinical trials. From a dataset of 36 multiple sclerosis patients with at least six monthly magnetic resonance imagings, we extracted new lesions and performed principal component analysis to estimate a biomarker that recapitulated lesion recovery. We tested the effect of multiple sclerosis disease modifying therapy on the lesion evolution index in three experimental designs and calculated sample sizes needed to appropriately power studies. Our proposed methods can be used to calculate required sample size and scan frequency in observational studies of multiple sclerosis disease progression as well as in designing clinical trials to find effects of treatment on multiple sclerosis lesion evolution.

Effect of disease-modifying therapies on subcortical gray matter atrophy in multiple sclerosis.

BACKGROUND: The effects of disease-modifying therapies (DMTs) on region-specific brain atrophy in multiple sclerosis (MS) are unclear. OBJECTIVE: To determine the effects of higher versus lower efficacy DMTs on rates of brain substructure atrophy in MS. METHODS: A non-randomized, observational cohort of people with MS followed with annual brain magnetic resonance imaging (MRI) was evaluated retrospectively. Whole brain, subcortical gray matter (GM), cortical GM, and cerebral white matter (WM) volume fractions were obtained. DMTs were categorized as higher (DMT-H: natalizumab and rituximab) or lower (DMT-L: interferon-beta and glatiramer acetate) efficacy. Follow-up epochs were analyzed if participants had been on a DMT for ⩾6 months prior to baseline and had at least one follow-up MRI while on DMTs in the same category. RESULTS: A total of 86 DMT epochs (DMT-H: n = 32; DMT-L: n = 54) from 78 participants fulfilled the study inclusion criteria. Mean follow-up was 2.4 years. Annualized rates of thalamic (-0.15% vs -0.81%; p = 0.001) and putaminal (-0.27% vs -0.73%; p = 0.001) atrophy were slower during DMT-H compared to DMT-L epochs. These results remained significant in multivariate analyses including demographics, clinical characteristics, and T2 lesion volume. CONCLUSION: DMT-H treatment may be associated with slower rates of subcortical GM atrophy, especially of the thalamus and putamen. Thalamic and putaminal volumes are promising imaging biomarkers in MS.

Evaluating the Impact of Intensity Normalization on MR Image Synthesis.

Image synthesis learns a transformation from the intensity features of an input image to yield a different tissue contrast of the output image. This process has been shown to have application in many medical image analysis tasks including imputation, registration, and segmentation. To carry out synthesis, the intensities of the input images are typically scaled-i.e., normalized-both in training to learn the transformation and in testing when applying the transformation, but it is not presently known what type of input scaling is optimal. In this paper, we consider seven different intensity normalization algorithms and three different synthesis methods to evaluate the impact of normalization. Our experiments demonstrate that intensity normalization as a preprocessing step improves the synthesis results across all investigated synthesis algorithms. Furthermore, we show evidence that suggests intensity normalization is vital for successful deep learning-based MR image synthesis.

DeepHarmony: A deep learning approach to contrast harmonization across scanner changes.

Magnetic resonance imaging (MRI) is a flexible medical imaging modality that often lacks reproducibility between protocols and scanners. It has been shown that even when care is taken to standardize acquisitions, any changes in hardware, software, or protocol design can lead to differences in quantitative results. This greatly impacts the quantitative utility of MRI in multi-site or long-term studies, where consistency is often valued over image quality. We propose a method of contrast harmonization, called DeepHarmony, which uses a U-Net-based deep learning architecture to produce images with consistent contrast. To provide training data, a small overlap cohort (n = 8) was scanned using two different protocols. Images harmonized with DeepHarmony showed significant improvement in consistency of volume quantification between scanning protocols. A longitudinal MRI dataset of patients with multiple sclerosis was also used to evaluate the effect of a protocol change on atrophy calculations in a clinical research setting. The results show that atrophy calculations were substantially and significantly affected by protocol change, whereas such changes have a less significant effect and substantially reduced overall difference when using DeepHarmony. This establishes that DeepHarmony can be used with an overlap cohort to reduce inconsistencies in segmentation caused by changes in scanner protocol, allowing for modernization of hardware and protocol design in long-term studies without invalidating previously acquired data.

Applications of a deep learning method for anti-aliasing and super-resolution in MRI.

Magnetic resonance (MR) images with both high resolutions and high signal-to-noise ratios (SNRs) are desired in many clinical and research applications. However, acquiring such images takes a long time, which is both costly and susceptible to motion artifacts. Acquiring MR images with good in-plane resolution and poor through-plane resolution is a common strategy that saves imaging time, preserves SNR, and provides one viewpoint with good resolution in two directions. Unfortunately, this strategy also creates orthogonal viewpoints that have poor resolution in one direction and, for 2D MR acquisition protocols, also creates aliasing artifacts. A deep learning approach called SMORE that carries out both anti-aliasing and super-resolution on these types of acquisitions using no external atlas or exemplars has been previously reported but not extensively validated. This paper reviews the SMORE algorithm and then demonstrates its performance in four applications with the goal to demonstrate its potential for use in both research and clinical scenarios. It is first shown to improve the visualization of brain white matter lesions in FLAIR images acquired from multiple sclerosis patients. Then it is shown to improve the visualization of scarring in cardiac left ventricular remodeling after myocardial infarction. Third, its performance on multi-view images of the tongue is demonstrated and finally it is shown to improve performance in parcellation of the brain ventricular system. Both visual and selected quantitative metrics of resolution enhancement are demonstrated.

Quantitative vibratory sensation measurement is related to sensory cortical thickness in MS.

Objective: Vibratory sensation is a quantifiable measure of physical dysfunction and is often related to spinal cord pathology; however, its association with relevant brain areas has not been fully explored. Our objective was to establish a cortical structural substrate for vibration sensation. Methods: Eighty-four individuals with multiple sclerosis (MS) (n = 54 relapsing, n = 30 progressive) and 28 controls participated in vibratory sensation threshold quantification at the great toe and a 3T MRI evaluating volume of the thalamus and cortical thickness primary and secondary sensory cortices. Results: After controlling for age, sex, and disability level, vibratory sensation thresholds were significantly related to cortical thickness of the anterior cingulate (P = 0.041), parietal operculum (P = 0.022), and inferior frontal gyrus pars operculum (P = 0.044), pars orbitalis (P = 0.007), and pars triangularis (P = 0.029). Within the progressive disease subtype, there were significant relationships between vibratory sensation and thalamic volume (P = 0.039) as well as reduced inferior frontal gyrus pars operculum (P = 0.014) and pars orbitalis (P = 0.005) cortical thickness. Conclusions: The data show significant independent relationships between quantitative vibratory sensation and measures of primary and secondary sensory cortices. Quantitative clinical measurement of vibratory sensation reflects pathological changes in spatially distinct brain areas and may supplement information captured by brain atrophy measures. Without overt relapses, monitoring decline in progressive forms of MS has proved challenging; quantitative clinical assessment may provide a tool to examine pathological decline in this cohort. These data suggest that quantitative clinical assessment may be a reliable way to examine pathological decline and have broader relevance to progressive forms of MS.