Feasibility of high-frequency percussions in people with severe acquired brain injury and tracheostomy: an observational study.

People with severe acquired brain injury (pwSABI) frequently experience pulmonary complications. Among these, atelectasis can occur as a result of pneumonia, thus increasing the chance of developing acute respiratory failure. Respiratory physiotherapy contribution to the management of atelectasis in pwSABI is yet poorly understood. We conducted a retrospective analysis on 15 non-cooperative pwSABI with tracheostomy and spontaneously breathing, hospitalized and treated with high-frequency percussion physiotherapy between September 2018 and February 2021 at the Neurological Rehabilitation Unit of the IRCCS "S.Maria Nascente - Fondazione Don Gnocchi", Milan. Our primary aim was to investigate the feasibility of such a physiotherapy intervention method. Then, we assessed changes in respiratory measures (arterial blood gas analysis and peripheral night-time oxygen saturation) and high-resolution computed tomography lung images, evaluated before and after the physiotherapy treatment. The radiological measures were a modified radiological atelectasis score (mRAS) assigned by two radiologists, and an opacity score automatically provided by the software CT Pneumonia Analysis® that identifies the regions of abnormal lung patterns. Treatment diaries showed that all treatments were completed, and no adverse events during treatment were registered. Among the 15 pwSABI analyzed, 8 were treated with IPV® and 7 with MetaNeb®. After a median of 14 (I-III quartile=12.5-14.5) days of treatment, we observed a statistical improvement in various arterial blood gas measures and peripheral night-time oxygen saturation measures. We also found radiological improvement or stability in more than 80% of pwSABI. In conclusion, our physiotherapy approach was feasible, and we observed respiratory parameters and radiological improvements. Using technology to assess abnormal tomographic patterns could be of interest to disentangle the short-term effects of respiratory physiotherapy on non-collaborating people.

Comparison of multicenter MRI protocols for visualizing the spinal cord gray matter.

PURPOSE: Spinal cord gray-matter imaging is valuable for a number of applications, but remains challenging. The purpose of this work was to compare various MRI protocols at 1.5 T, 3 T, and 7 T for visualizing the gray matter. METHODS: In vivo data of the cervical spinal cord were collected from nine different imaging centers. Data processing consisted of automatically segmenting the spinal cord and its gray matter and co-registering back-to-back scans. We computed the SNR using two methods (SNR_single using a single scan and SNR_diff using the difference between back-to-back scans) and the white/gray matter contrast-to-noise ratio per unit time. Synthetic phantom data were generated to evaluate the metrics performance. Experienced radiologists qualitatively scored the images. We ran the same processing on an open-access multicenter data set of the spinal cord MRI (N = 267 participants). RESULTS: Qualitative assessments indicated comparable image quality for 3T and 7T scans. Spatial resolution was higher at higher field strength, and image quality at 1.5 T was found to be moderate to low. The proposed quantitative metrics were found to be robust to underlying changes to the SNR and contrast; however, the SNR_single method lacked accuracy when there were excessive partial-volume effects. CONCLUSION: We propose quality assessment criteria and metrics for gray-matter visualization and apply them to different protocols. The proposed criteria and metrics, the analyzed protocols, and our open-source code can serve as a benchmark for future optimization of spinal cord gray-matter imaging protocols.

Integrating large-scale neuroimaging research datasets: Harmonisation of white matter hyperintensity measurements across Whitehall and UK Biobank datasets.

Large scale neuroimaging datasets present the possibility of providing normative distributions for a wide variety of neuroimaging markers, which would vastly improve the clinical utility of these measures. However, a major challenge is our current poor ability to integrate measures across different large-scale datasets, due to inconsistencies in imaging and non-imaging measures across the different protocols and populations. Here we explore the harmonisation of white matter hyperintensity (WMH) measures across two major studies of healthy elderly populations, the Whitehall II imaging sub-study and the UK Biobank. We identify pre-processing strategies that maximise the consistency across datasets and utilise multivariate regression to characterise study sample differences contributing to differences in WMH variations across studies. We also present a parser to harmonise WMH-relevant non-imaging variables across the two datasets. We show that we can provide highly calibrated WMH measures from these datasets with: (1) the inclusion of a number of specific standardised processing steps; and (2) appropriate modelling of sample differences through the alignment of demographic, cognitive and physiological variables. These results open up a wide range of applications for the study of WMHs and other neuroimaging markers across extensive databases of clinical data.

Quantitative MRI using STrategically Acquired Gradient Echo (STAGE): optimization for 1.5 T scanners and T1 relaxation map validation.

OBJECTIVES: The strategically acquired gradient echo (STAGE) protocol, developed for 3T scanners, allows one to derive quantitative maps such as T1, T2*, proton density, and quantitative susceptibility mapping in about 5 min. Our aim was to adapt the STAGE sequences for 1.5T scanners which are still commonly used in clinical practice. Furthermore, the accuracy and repeatability of the STAGE-derived T1 estimate were tested. METHODS: Flip angle (FA) optimization was performed using a theoretical simulation by maximizing signal-to-noise ratio, contrast-to-noise ratio, and T1 precision. The FA choice was further refined with the ISMRM/NIST phantom and in vivo acquisitions. The accuracy of the T1 estimate was assessed by comparing STAGE-derived T1 values with T1 maps obtained with an inversion recovery sequence. T1 accuracy was investigated for both the phantom and in vivo data. Finally, one subject was acquired 10 times once a week and a group of 27 subjects was scanned once. The T1 coefficient of variation (COV) was computed to assess scan-rescan and physiological variability, respectively. RESULTS: The FA1,2 = 7°,38° were identified as the optimal FA pair at 1.5T. The T1 estimate errors were below 3% and 5% for phantom and in vivo measurements, respectively. COV for different tissues ranged from 1.8 to 4.8% for physiological variability, and between 0.8 and 2% for scan-rescan repeatability. CONCLUSION: The optimized STAGE protocol can provide accurate and repeatable T1 mapping along with other qualitative images and quantitative maps in about 7 min on 1.5T scanners. This study provides the groundwork to assess the role of STAGE in clinical settings. KEY POINTS: • The STAGE imaging protocol was optimized for use on 1.5T field strength scanners. • A practical STAGE protocol makes it possible to derive quantitative maps (i.e., T1, T2*, PD, and QSM) in about 7 min at 1.5T. • The T1 estimate derived from the STAGE protocol showed good accuracy and repeatability.

Combined DTI-fMRI Analysis for a Quantitative Assessment of Connections Between WM Bundles and Their Peripheral Cortical Fields in Verbal Fluency.

Diffusion tensor imaging (DTI) tractography and functional magnetic resonance imaging (fMRI) are powerful techniques to elucidate the anatomical and functional aspects of brain connectivity. However, integrating these approaches to describe the precise link between structure and function within specific brain circuits remains challenging. In this study, a novel DTI-fMRI integration method is proposed, to provide the topographical characterization and the volumetric assessment of the functional and anatomical connections within the language circuit. In a group of 21 healthy elderly subjects (mean age 68.5 ± 5.8 years), the volume of connection between the cortical activity elicited by a verbal fluency task and the cortico-cortical fiber tracts associated with this function are mapped and quantified. An application of the method to a case study in neuro-rehabilitation context is also presented. Integrating structural and functional data within the same framework, this approach provides an overall view of white and gray matter when studying specific brain circuits.

Corticospinal tract integrity is related to primary motor cortex thinning in relapsing-remitting multiple sclerosis.

BACKGROUND: The relationship between white matter injury and cortical atrophy development in relapsing-remitting multiple sclerosis (RRMS) remains unclear. OBJECTIVES: To investigate the associations between corticospinal tract integrity and cortical morphology measures of the primary motor cortex in RRMS patients and healthy controls. METHODS: 51 RRMS patients and 30 healthy controls underwent MRI examination for cortical reconstruction and assessment of corticospinal tract integrity. Partial correlation and multiple linear regression analyses were used to investigate the associations of focal and normal appearing white matter (NAWM) injury of the corticospinal tract with thickness and surface area measures of the primary motor cortex. Relationships between MRI measures and clinical disability as assessed by the Expanded Disability Status Scale and disease duration were also investigated. RESULTS: In patients only, decreased cortical thickness was related to increased corticospinal tract NAWM mean, axial and radial diffusivities in addition to corticospinal tract lesion volume. The final multiple linear regression model for PMC thickness retained only NAWM axial diffusivity as a significant predictor (adjusted R(2)= 0.270, p= 0.001). Clinical measures were associated with NAWM corticospinal tract integrity measures. CONCLUSIONS: Primary motor cortex thinning in RRMS is related to alterations in connected white matter and is best explained by decreased NAWM integrity.

An anatomy-based lumped parameter model of cerebrospinal venous circulation: can an extracranial anatomical change impact intracranial hemodynamics?

BACKGROUND: The relationship between extracranial venous system abnormalities and central nervous system disorders has been recently theorized. In this paper we delve into this hypothesis by modeling the venous drainage in brain and spinal column areas and simulating the intracranial flow changes due to extracranial morphological stenoses. METHODS: A lumped parameter model of the cerebro-spinal venous drainage was created based on anatomical knowledge and vessels diameters and lengths taken from literature. Each vein was modeled as a hydraulic resistance, calculated through Poiseuille's law. The inputs of the model were arterial flow rates of the intracranial, vertebral and lumbar districts. The effects of the obstruction of the main venous outflows were simulated. A database comprising 112 Multiple Sclerosis patients (Male/Female = 42/70; median age ± standard deviation = 43.7 ± 10.5 years) was retrospectively analyzed. RESULTS: The flow rate of the main veins estimated with the model was similar to the measures of 21 healthy controls (Male/Female = 10/11; mean age ± standard deviation = 31 ± 11 years), obtained with a 1.5 T Magnetic Resonance scanner. The intracranial reflux topography predicted with the model in cases of internal jugular vein diameter reduction was similar to those observed in the patients with internal jugular vein obstacles. CONCLUSIONS: The proposed model can predict physiological and pathological behaviors with good fidelity. Despite the simplifications introduced in cerebrospinal venous circulation modeling, the key anatomical feature of the lumped parameter model allowed for a detailed analysis of the consequences of extracranial venous impairments on intracranial pressure and hemodynamics.

Grey matter damage in progressive multiple sclerosis versus amyotrophic lateral sclerosis: a voxel-based morphometry MRI study.

Primary progressive multiple sclerosis (PPMS) and amyotrophic lateral sclerosis (ALS) seem to share some clinical and pathological features. MRI studies revealed the presence of grey matter (GM) atrophy in both diseases, but no comparative data are available. The objective was to compare the regional patterns of GM tissue loss in PPMS and ALS with voxel-based morphometry (VBM). Eighteen PPMS patients, 20 ALS patients, and 31 healthy controls (HC) were studied with a 1.5 Tesla scanner. VBM was performed to assess volumetric GM differences with age and sex as covariates. Threshold-free cluster enhancement analysis was used to obtain significant clusters. Group comparisons were tested with family-wise error correction for multiple comparisons (p < 0.05) except for HC versus MND which was tested at a level of p < 0.001 uncorrected and a cluster threshold of 20 contiguous voxels. Compared to HC, ALS patients showed GM tissue reduction in selected frontal and temporal areas, while PPMS patients showed a widespread bilateral GM volume decrease, involving both deep and cortical regions. Compared to ALS, PPMS patients showed tissue volume reductions in both deep and cortical GM areas. This preliminary study confirms that PPMS is characterized by a more diffuse cortical and subcortical GM atrophy than ALS and that, in the latter condition, brain damage is present outside the motor system. These results suggest that PPMS and ALS may share pathological features leading to GM tissue loss.

Real-Time Phase-Contrast MRI to Monitor Cervical Blood and Cerebrospinal Fluid Flow Beat-by-Beat Variability.

Beat-by-beat variability (BBV) rhythms are observed in both cardiovascular (CV) and intracranial (IC) compartments, yet interactions between the two are not fully understood. Real-Time Phase-Contrast (RT-PC) MRI sequence was acquired for 30 healthy volunteers at 1st cervical level on a 3T scanner. The arterial (AF), venous (VF), and cerebrospinal fluid (CSF) flow (CSFF) were computed as velocity integrals over the internal carotid artery, internal jugular vein, and CSF. AF, VF, and CSFF signals were segmented in inspiration and expiration beats, to assess the respiration influence. Systolic and diastolic BBV, and heart period series underwent autoregressive power spectral density analysis, to evaluate the low-frequency (LF, Mayer waves) and high frequency (HF, respiratory waves) components. The diastolic VF had the largest BBV. LF power was high in the diastolic AF series, poor in all CSFF series. The pulse wave analyses revealed higher mean amplitude during inspiration. Findings suggests a possible role of LF modulation of IC resistances and propagation of HF waves from VF to AF and CCSF. PC-RT-MRI could provide new insight into the interaction between CV and IC regulation and pave the way for a detailed analysis of the cerebrovascular effects of varied respiration patterns due to exercise and rehabilitation.

Development and Testing of SPIDER-NET: An Interactive Tool for Brain Connectogram Visualization, Sub-Network Exploration and Graph Metrics Quantification.

Brain connectomics consists in the modeling of human brain as networks, mathematically represented as numerical connectivity matrices. However, this representation may result in difficult interpretation of the data. To overcome this limitation, graphical representation by connectograms is currently used via open-source tools, which, however, lack user-friendly interfaces and options to explore specific sub-networks. In this context, we developed SPIDER-NET (Software Package Ideal for Deriving Enhanced Representations of brain NETworks), an easy-to-use, flexible, and interactive tool for connectograms generation and sub-network exploration. This study aims to present SPIDER-NET and to test its potential impact on pilot cases. As a working example, structural connectivity (SC) was investigated with SPIDER-NET in a group of 17 healthy controls (HCs) and in two subjects with stroke injury (Case 1 and Case 2, both with a focal lesion affecting part of the right frontal lobe, insular cortex and subcortical structures). 165 parcels were determined from individual structural magnetic resonance imaging data by using the Destrieux atlas, and defined as nodes. SC matrices were derived with Diffusion Tensor Imaging tractography. SC matrices of HCs were averaged to obtain a single group matrix. SC matrices were then used as input for SPIDER-NET. First, SPIDER-NET was used to derive the connectogram of the right hemisphere of Case 1 and Case 2. Then, a sub-network of interest (i.e., including gray matter regions affected by the stroke lesions) was interactively selected and the associated connectograms were derived for Case 1, Case 2 and HCs. Finally, graph-based metrics were derived for whole-brain SC matrices of Case 1, Case 2 and HCs. The software resulted effective in representing the expected (dis) connectivity pattern in the hemisphere affected by the stroke lesion in Cases 1 and 2. Furthermore, SPIDER-NET allowed to test an a priori hypothesis by interactively extracting a sub-network of interest: Case 1 showed a sub-network connectivity pattern different from Case 2, reflecting the different clinical severity. Global and local graph-based metrics derived with SPIDER-NET were different between cases with stroke injury and HCs. The tool proved to be accessible, intuitive, and interactive in brain connectivity investigation and provided both qualitative and quantitative evidence.

Mid-term MRI evaluation reveals microstructural white matter alterations in COVID-19 fully recovered subjects with anosmia presentation.

Background: Little is still known about the mid/long-term effects of coronavirus disease 2019 (COVID-19) on the brain, especially in subjects who have never been hospitalized due to the infection. In this neuroimaging exploratory study, we analyzed the medium-term effect of COVID-19 on the brain of people who recovered from COVID-19, experienced anosmia during the acute phase of the disease, and have never been hospitalized due to SARS-Co-V-2 infection. Methods: Forty-three individuals who had (COV+, n = 22) or had not (COV-, n = 21) been infected with SARS-Co-V-2 were included in the study; the two groups were age- and sex-matched and were investigated using 3T magnetic resonance imaging (MRI). Gray matter (GM) volume, white matter (WM) hyperintensity volume, WM microstrutural integrity (i.e. fractional anisotropy [FA], mean diffusivity [MD], axial diffusivity [AD], radial diffusivity [RD]) and cerebral blood flow (CBF) differences between the two groups were tested with either analysis of covariance or voxel-wise analyses. Results were family wise error (FWE) corrected. Results: No significant differences between COV+ and COV- groups were observed in terms of GM volume, WM hyperintensity volume, and CBF. Conversely, local WM microstructural alterations were detected in COV+ when compared with COV- with tract-based spatial statistics. Specifically, COV+ showed lower FA (pFWE-peak = 0.035) and higher RD (pFWE-peak = 0.038) than COV- in several WM regions. Conclusion: COVID-19 may produce mid/long-term microstructural effect on the brain, even in case of mild-to-moderate disease not requiring hospitalization. Further investigation and additional follow-ups are warranted to assess if the alterations reported in this study totally recover over time. As brain alterations could increase the risk of cognitive decline, greater knowledge of their trajectories is crucial to aid neurorehabilitation treatments.

Blood and cerebrospinal fluid flow oscillations measured with real-time phase-contrast MRI: breathing mode matters.

BACKGROUND: Cervical blood and cerebrospinal fluid (CSF) flow rates can be quantified with Phase-contrast (PC) MRI, which is routinely used for clinical studies. Previous MRI studies showed that venous and CSF flow alterations are linked to various pathological conditions. Since it is well known that, besides the heart beating, the thoracic pump influences the blood and CSF dynamics, we studied the effect of different respiration modes on blood and CSF flow rates using a real-time (RT)-PC prototype. METHODS: Thirty healthy volunteers were examined with a 3 T scanner. A RT-PC sequence was acquired at the first cervical level to quantify the flow rates of internal carotid arteries, internal jugular veins (IJVs) and CSF. Each RT-PC acquisition was repeated three times, while the subjects were asked to breathe in three different ways for 60 s each: freely (F), with a constant rate (PN) and with deep and constant respiration rate (PD). The average flow rates were computed, they were removed from the respective signals and integrated in the inspiratory and expiratory phases (differential volumes). Finally, the power spectral density was computed for each detrended flow rate. High- and very-high frequency peaks were identified on the spectra while their frequencies were compared to the respiratory and cardiac frequencies estimated using a thoracic belt and a pulse oximeter. The area under the spectra was computed in four 0.5 Hz-wide ranges, centered on the high-frequency peak, on very-high frequency peak and its 2nd and 3rd harmonics, and then they were normalized by the flow rate variance. The effect of breathing patterns on average flow rates, on systolic and diastolic peaks, and on the normalized power was tested. Finally, the differential volumes of inspiration were compared to those of expiration. RESULTS: The frequencies of the high- and very-high spectral peaks corresponded to the respiratory and cardiac frequencies. The average flow rate progressively decreased from F to PN to PD breathing, and the cardiac modulations were less predominant especially for the IJVs. The respiratory modulation increased with PD breathing. The average volumes displaced in the inspiratory phases were not significantly different from those of the expiratory one. CONCLUSIONS: The spectral analyses demonstrated higher respiratory modulations in PD compared to free breathing, even prevailing the cardiac modulation in the IJVs, showing an increment of the thoracic pump affecting the flow rate shape.

Cardiac and Respiratory Influences on Intracranial and Neck Venous Flow, Estimated Using Real-Time Phase-Contrast MRI.

The study of brain venous drainage has gained attention due to its hypothesized link with various neurological conditions. Intracranial and neck venous flow rate may be estimated using cardiac-gated cine phase-contrast (PC)-MRI. Although previous studies showed that breathing influences the neck's venous flow, this aspect could not be studied using the conventional segmented PC-MRI since it reconstructs a single cardiac cycle. The advent of real-time PC-MRI has overcome these limitations. Using this technique, we measured the internal jugular veins and superior sagittal sinus flow rates in a group of 16 healthy subjects (12 females, median age of 23 years). Comparing forced-breathing and free-breathing, the average flow rate decreased and the respiratory modulation increased. The flow rate decrement may be due to a vasoreactive response to deep breathing. The respiratory modulation increment is due to the thoracic pump's greater effect during forced breathing compared to free breathing. These results showed that the breathing mode influences the average blood flow and its pulsations. Since effective drainage is fundamental for brain health, rehabilitative studies might use the current setup to investigate if respiratory exercises positively affect clinical variables and venous drainage.

Quantitative MRI Harmonization to Maximize Clinical Impact: The RIN-Neuroimaging Network.

Neuroimaging studies often lack reproducibility, one of the cardinal features of the scientific method. Multisite collaboration initiatives increase sample size and limit methodological flexibility, therefore providing the foundation for increased statistical power and generalizable results. However, multisite collaborative initiatives are inherently limited by hardware, software, and pulse and sequence design heterogeneities of both clinical and preclinical MRI scanners and the lack of benchmark for acquisition protocols, data analysis, and data sharing. We present the overarching vision that yielded to the constitution of RIN-Neuroimaging Network, a national consortium dedicated to identifying disease and subject-specific in-vivo neuroimaging biomarkers of diverse neurological and neuropsychiatric conditions. This ambitious goal needs efforts toward increasing the diagnostic and prognostic power of advanced MRI data. To this aim, 23 Italian Scientific Institutes of Hospitalization and Care (IRCCS), with technological and clinical specialization in the neurological and neuroimaging field, have gathered together. Each IRCCS is equipped with high- or ultra-high field MRI scanners (i.e., ≥3T) for clinical or preclinical research or has established expertise in MRI data analysis and infrastructure. The actions of this Network were defined across several work packages (WP). A clinical work package (WP1) defined the guidelines for a minimum standard clinical qualitative MRI assessment for the main neurological diseases. Two neuroimaging technical work packages (WP2 and WP3, for clinical and preclinical scanners) established Standard Operative Procedures for quality controls on phantoms as well as advanced harmonized quantitative MRI protocols for studying the brain of healthy human participants and wild type mice. Under FAIR principles, a web-based e-infrastructure to store and share data across sites was also implemented (WP4). Finally, the RIN translated all these efforts into a large-scale multimodal data collection in patients and animal models with dementia (i.e., case study). The RIN-Neuroimaging Network can maximize the impact of public investments in research and clinical practice acquiring data across institutes and pathologies with high-quality and highly-consistent acquisition protocols, optimizing the analysis pipeline and data sharing procedures.

Measuring Aqueduct of Sylvius Cerebrospinal Fluid Flow in Multiple Sclerosis Using Different Software.

Aqueduct of Sylvius (AoS) cerebrospinal fluid flow can be quantified using phase-contrast (PC) Magnetic Resonance Imaging. The software used for AoS segmentation might affect the PC-derived measures. We analyzed AoS PC data of 30 people with multiple sclerosis and 19 normal controls using three software packages, and estimated cross-sectional area (CSA), average and highest AoS velocity (Vmean and Vmax), flow rate and volume. Our aims were to assess the repeatability and reproducibility of each PC-derived measure obtained with the various software packages, including in terms of group differentiation. All the variables had good repeatability, except the average Vmean, flow rate and volume obtained with one software package. Substantial to perfect agreement was seen when evaluating the overlap between the AoS segmentations obtained with different software packages. No variable was significantly different between software packages, with the exception of Vmean diastolic peak and CSA. Vmax diastolic peak differentiated groups, regardless of the software package. In conclusion, a clinical study should preliminarily evaluate the repeatability in order to interpret its findings. Vmax seemed to be a repeatable and reproducible measure, since the pixel with its value is usually located in the center of the AoS, and is thus unlikely be affected by ROI size.

Automated Assessment of the Substantia Nigra Pars Compacta in Parkinson's Disease: A Diffusion Tensor Imaging Study.

The substantia nigra (SN) pars compacta (SNpc) and pars reticulata (SNpr) are differentially affected in Parkinson's disease (PD). Separating the SNpc and SNpr is challenging with standard magnetic resonance imaging (MRI). Diffusion tensor imaging (DTI) allows for the characterization of SN microstructure in a non-invasive manner. In this study, 29 PD patients and 28 healthy controls (HCs) were imaged with 1.5T MRI for DTI. Images were nonlinearly registered to standard space and SNpc and SNpr DTI parameters were measured. ANCOVA and receiver operator characteristic (ROC) analyses were performed. Clinical associations were assessed with Spearman correlations. Multiple corrections were controlled for false discovery rate. PD patients presented with significantly increased SNpc axial diffusivity (AD) (1.207 ± 0.068 versus 1.156 ± 0.045, p = 0.024), with ROC analysis yielding an under the curve of 0.736. Trends with Unified Parkinson's Disease Rating Scale (UPDRS) III scores were identified for SNpc MD (rs = 0.449), AD (rs = 0.388), and radial diffusivity (rs = 0.391) (all p < 0.1). A trend between baseline SNpr MD and H&Y change (rs = 0.563, p = 0.081) over 2.9 years of follow-up was identified (n = 14). In conclusion, SN microstructure shows robust, clinically meaningful associations in PD.

Relationship between MRI perfusion and clinical severity in multiple sclerosis.

Perfusion alterations within several brain regions have been shown in multiple sclerosis patients using different magnetic resonance imaging (MRI) techniques. Furthermore, MRI-derived brain perfusion metrics have been investigated in association with multiple sclerosis phenotypes, physical disability, and cognitive impairment. However, a review focused on these aspects is still missing. Our aim was to review all the studies investigating the relationship between perfusion MRI and clinical severity during the last fifteen years to understand the clinical relevance of these findings. Perfusion differences among phenotypes were observed both with 1.5T and 3T scanners, with progressive multiple sclerosis presenting with lower perfusion values than relapsing-remitting multiple sclerosis patients. However, only 3T scanners showed a statistically significant distinction. Controversial results about the association between MRI-derived perfusion metrics and physical disability scores were found. However, the majority of the studies showed that lower brain perfusion and longer transit time are associated with more severe physical disability and worse cognitive performances.