Breath-hold BOLD fMRI without CO2 sampling enables estimation of venous cerebral blood volume: potential use in normalization of stimulus-evoked BOLD fMRI data.

BOLD fMRI signal has been used in conjunction with vasodilatory stimulation as a marker of cerebrovascular reactivity (CVR): the relative change in cerebral blood flow (CBF) arising from a unit change in the vasodilatory stimulus. Using numerical simulations, we demonstrate that the variability in the relative BOLD signal change induced by vasodilation is strongly influenced by the variability in deoxyhemoglobin-containing cerebral blood volume (CBV), as this source of variability is likely to be more prominent than that of CVR. It may, therefore, be more appropriate to describe the relative BOLD signal change induced by an isometabolic vasodilation as a proxy of deoxygenated CBV (CBVdHb) rather than CVR. With this in mind, a new method was implemented to map a marker of CBVdHb, termed BOLD-CBV, based on the normalization of voxel-wise BOLD signal variation by an estimate of the intravascular venous BOLD signal from voxels filled with venous blood. The intravascular venous BOLD signal variation, recorded during repeated breath-holding, was extracted from the superior sagittal sinus in a cohort of 27 healthy volunteers and used as a regressor across the whole brain, yielding maps of BOLD-CBV. In the same cohort, we demonstrated the potential use of BOLD-CBV for the normalization of stimulus-evoked BOLD fMRI by comparing group-level BOLD fMRI responses to a visuomotor learning task with and without the inclusion of voxel-wise vascular covariates of BOLD-CBV and the BOLD signal change per mmHg variation in end-tidal carbon dioxide (BOLD-CVR). The empirical measure of BOLD-CBV accounted for more between-subject variability in the motor task-induced BOLD responses than BOLD-CVR estimated from end-tidal carbon dioxide recordings. The new method can potentially increase the power of group fMRI studies by including a measure of vascular characteristics and has the strong practical advantage of not requiring experimental measurement of end-tidal carbon dioxide, unlike traditional methods to estimate BOLD-CVR. It also more closely represents a specific physiological characteristic of brain vasculature than BOLD-CVR, namely blood volume.

Changes of brain parenchyma free water fraction reflect tissue damage and impaired processing speed in multiple sclerosis.

Free water fraction (FWF) represents the amount of water per unit volume of brain parenchyma, which is not bound to macromolecules. Its excess in multiple sclerosis (MS) is related to increased tissue loss. The use of mcDESPOT (multicomponent driven single pulse observation of T1 and T2), a 3D imaging method which exploits both the T1 and T2 contrasts, allows FWF to be derived in clinically feasible times. However, this method has not been used to quantify changes of FWF and their potential clinical impact in MS. The aim of this study is to investigate the changes in FWF in MS patients and their relationship with tissue damage and cognition, under the hypothesis that FWF is a proxy of clinically meaningful tissue loss. To this aim, we tested the relationship between FWF, MS lesion burden and information processing speed, evaluated via the Symbol Digit Modalities Test (SDMT). In addition to standard sequences, used for T1- and T2-weighted lesion delineation, the mcDESPOT sequence with 1.7 mm isotropic resolution and a diffusion weighted imaging protocol (b = 0, 1200 s/mm2, 40 diffusion directions) were employed at 3 T. The fractional anisotropy map derived from diffusion data was used to define a subject-specific white matter (WM) atlas. Brain parenchyma segmentation returned masks of gray matter (GM) and WM, and normal-appearing WM (NAWM), in addition to the T1 and T2 lesion masks (T1L and T2L, respectively). Ninety-nine relapsing-remitting MS patients (age = 43.3 ± 9.9 years, disease duration 12.3 ± 7.7 years) were studied, together with twenty-five healthy controls (HC, age = 38.8 ± 11.0 years). FWF was higher in GM and NAWM of MS patients, compared to GM and WM of HC (both p < .001). In MS patients, FWF was the highest in the T1L and GM, followed by T2L and NAWM, respectively. FWF increased significantly with T1L and T2L volume (ρ ranging from 0.40 to 0.58, p < .001). FWF in T2L was strongly related to both T1L volume and the volume ratio T1L/T2L (ρ = 0.73, p < .001). MS patients performed worse than HC in the processing speed test (mean ± SD: 54.1 ± 10.3 for MS, 63.8 ± 10.8 for HC). FWF in GM, T2L, perilesional tissue and NAWM increased with SDMT score reduction (ρ = -0.30, -0.29, -0.33 respectively and r = -.30 for T2L, all with p < .005). A regional analysis, conducted to determine which NAWM regions were of particular importance to explain the relationship between FWF and cognitive impairment, revealed that FWF spatial variance was negatively related to SDMT score in the corpus callosum and the superior longitudinal fasciculus, WM structures known to be associated with cognitive impairment, in addition to the left corticospinal tract, the sagittal stratum, the right anterior limb of internal capsule. In conclusion, we found excess free water in brain parenchyma of MS patients, an alteration that involved not only MS lesions, but also the GM and NAWM, impinging on brain function and negatively associated with cognitive processing speed. We suggest that the FWF metric, derived from noninvasive, rapid MRI acquisitions and bearing good biological interpretability, may prove valuable as an MRI biomarker of tissue damage and associated cognitive impairment in MS.

Quantitative myelin imaging with MRI and PET: an overview of techniques and their validation status.

Myelin is the protective sheath wrapped around axons, consisting of a phospholipid bilayer with water between the wraps. The measurement of damage to the myelin sheaths, the evaluation of the efficacy of therapies aiming to promote remyelination and monitoring the degree of brain maturation in children all require non-invasive quantitative myelin imaging methods. To date, various myelin imaging techniques have been developed. Five different MRI approaches can be distinguished based on their biophysical principles: (i) imaging of the water between the lipid bilayers directly (e.g. myelin water imaging); (ii) imaging the non-aqueous protons of the phospholipid bilayer directly with ultra-short echo-time techniques; (iii) indirect imaging of the macromolecular content (e.g. magnetization transfer; inhomogeneous magnetization transfer); (iv) mapping of the effects of the myelin sheath's magnetic susceptibility on the MRI signal (e.g. quantitative susceptibility mapping); and (v) mapping of the effects of the myelin sheath on water diffusion. Myelin imaging with PET uses radioactive molecules with high affinity to specific myelin components, in particular myelin basic protein. This review aims to give an overview of the various myelin imaging techniques, their biophysical principles, image acquisition, data analysis and their validation status.

Cerebral oxygen metabolism from MRI susceptibility.

This article provides an overview of MRI methods exploiting magnetic susceptibility properties of blood to assess cerebral oxygen metabolism, including the tissue oxygen extraction fraction (OEF) and the cerebral metabolic rate of oxygen (CMRO2). The first section is devoted to describing blood magnetic susceptibility and its effect on the MRI signal. Blood circulating in the vasculature can have diamagnetic (oxyhemoglobin) or paramagnetic properties (deoxyhemoglobin). The overall balance between oxygenated and deoxygenated hemoglobin determines the induced magnetic field which, in turn, modulates the transverse relaxation decay of the MRI signal via additional phase accumulation. The following sections of this review then illustrate the principles underpinning susceptibility-based techniques for quantifying OEF and CMRO2. Here, it is detailed whether these techniques provide global (OxFlow) or local (Quantitative Susceptibility Mapping - QSM, calibrated BOLD - cBOLD, quantitative BOLD - qBOLD, QSM+qBOLD) measurements of OEF or CMRO2, and what signal components (magnitude or phase) and tissue pools they consider (intravascular or extravascular). Validations studies and potential limitations of each method are also described. The latter include (but are not limited to) challenges in the experimental setup, the accuracy of signal modeling, and assumptions on the measured signal. The last section outlines the clinical uses of these techniques in healthy aging and neurodegenerative diseases and contextualizes these reports relative to results from gold-standard PET.

Quantitative susceptibility mapping identifies hippocampal and other subcortical grey matter tissue composition changes in temporal lobe epilepsy.

Temporal lobe epilepsy (TLE) is associated with widespread brain alterations. Using quantitative susceptibility mapping (QSM) alongside transverse relaxation rate ( R 2 * ), we investigated regional brain susceptibility changes in 36 patients with left-sided (LTLE) or right-sided TLE (RTLE) secondary to hippocampal sclerosis, and 27 healthy controls (HC). We compared three susceptibility calculation methods to ensure image quality. Correlations of susceptibility and R 2 * with age of epilepsy onset, frequency of focal-to-bilateral tonic-clonic seizures (FBTCS), and neuropsychological test scores were examined. Weak-harmonic QSM (WH-QSM) successfully reduced noise and removed residual background field artefacts. Significant susceptibility increases were identified in the left putamen in the RTLE group compared to the LTLE group, the right putamen and right thalamus in the RTLE group compared to HC, and a significant susceptibility decrease in the left hippocampus in LTLE versus HC. LTLE patients who underwent epilepsy surgery showed significantly lower left-versus-right hippocampal susceptibility. Significant R 2 * changes were found between TLE and HC groups in the amygdala, putamen, thalamus, and in the hippocampus. Specifically, decreased R2 * was found in the left and right hippocampus in LTLE and RTLE, respectively, compared to HC. Susceptibility and R 2 * were significantly correlated with cognitive test scores in the hippocampus, globus pallidus, and thalamus. FBTCS frequency correlated positively with ipsilateral thalamic and contralateral putamen susceptibility and with R 2 * in bilateral globi pallidi. Age of onset was correlated with susceptibility in the hippocampus and putamen, and with R 2 * in the caudate. Susceptibility and R 2 * changes observed in TLE groups suggest selective loss of low-myelinated neurons alongside iron redistribution in the hippocampi, predominantly ipsilaterally, indicating QSM's sensitivity to local pathology. Increased susceptibility and R 2 * in the thalamus and putamen suggest increased iron content and reflect disease severity.

Pathophysiology of multiple sclerosis damage and repair: Linking cerebral hypoperfusion to the development of irreversible tissue loss in multiple sclerosis using magnetic resonance imaging.

BACKGROUND AND PURPOSE: Reduced cerebral perfusion has been observed in multiple sclerosis (MS) and may contribute to tissue loss both acutely and chronically. Here, we test the hypothesis that hypoperfusion occurs in MS and relates to the presence of irreversible tissue damage. METHODS: In 91 patients with relapsing MS and 26 healthy controls (HC), gray matter (GM) cerebral blood flow (CBF) was assessed using pulsed arterial spin labeling. GM volume, T1 hypointense and T2 hyperintense lesion volumes (T1LV and T2LV, respectively), and the proportion of T2-hyperintense lesion volume that appears hypointense on T1-weighted magnetic resonance imaging (T1LV/T2LV) were quantified. GM CBF and GM volume were evaluated globally, as well as regionally, using an atlas-based approach. RESULTS: Global GM CBF was lower in patients (56.9 ± 12.3 mL/100 g/min) than in HC (67.7 ± 10.0 mL/100 g/min; p < 0.001), a difference that was widespread across brain regions. Although total GM volume was comparable between groups, significant reductions were observed in a subset of subcortical structures. GM CBF negatively correlated with T1LV (r = -0.43, p = 0.0002) and T1LV/T2LV (r = -0.37, p = 0.0004), but not with T2LV. CONCLUSIONS: GM hypoperfusion occurs in MS and is associated with irreversible white matter damage, thus suggesting that cerebral hypoperfusion may actively contribute and possibly precede neurodegeneration by hampering tissue repair abilities in MS.

Multi-echo quantitative susceptibility mapping: how to combine echoes for accuracy and precision at 3 Tesla.

PURPOSE: To compare different multi-echo combination methods for MRI QSM. Given the current lack of consensus, we aimed to elucidate how to optimally combine multi-echo gradient-recalled echo signal phase information, either before or after applying Laplacian-base methods (LBMs) for phase unwrapping or background field removal. METHODS: Multi-echo gradient-recalled echo data were simulated in a numerical head phantom, and multi-echo gradient-recalled echo images were acquired at 3 Tesla in 10 healthy volunteers. To enable image-based estimation of gradient-recalled echo signal noise, 5 volunteers were scanned twice in the same session without repositioning. Five QSM processing pipelines were designed: 1 applied nonlinear phase fitting over TEs before LBMs; 2 applied LBMs to the TE-dependent phase and then combined multiple TEs via either TE-weighted or SNR-weighted averaging; and 2 calculated TE-dependent susceptibility maps via either multi-step or single-step QSM and then combined multiple TEs via magnitude-weighted averaging. Results from different pipelines were compared using visual inspection; summary statistics of susceptibility in deep gray matter, white matter, and venous regions; phase noise maps (error propagation theory); and, in the healthy volunteers, regional fixed bias analysis (Bland-Altman) and regional differences between the means (nonparametric tests). RESULTS: Nonlinearly fitting the multi-echo phase over TEs before applying LBMs provided the highest regional accuracy of χ $$ \chi $$ and the lowest phase noise propagation compared to averaging the LBM-processed TE-dependent phase. This result was especially pertinent in high-susceptibility venous regions. CONCLUSION: For multi-echo QSM, we recommend combining the signal phase by nonlinear fitting before applying LBMs.

Development, validation, qualification, and dissemination of quantitative MR methods: Overview and recommendations by the ISMRM quantitative MR study group.

On behalf of the International Society for Magnetic Resonance in Medicine (ISMRM) Quantitative MR Study Group, this article provides an overview of considerations for the development, validation, qualification, and dissemination of quantitative MR (qMR) methods. This process is framed in terms of two central technical performance properties, i.e., bias and precision. Although qMR is confounded by undesired effects, methods with low bias and high precision can be iteratively developed and validated. For illustration, two distinct qMR methods are discussed throughout the manuscript: quantification of liver proton-density fat fraction, and cardiac T1 . These examples demonstrate the expansion of qMR methods from research centers toward widespread clinical dissemination. The overall goal of this article is to provide trainees, researchers, and clinicians with essential guidelines for the development and validation of qMR methods, as well as an understanding of necessary steps and potential pitfalls for the dissemination of quantitative MR in research and in the clinic.

Deep Learning-Based Neuromelanin MRI Changes of Isolated REM Sleep Behavior Disorder.

BACKGROUND: Isolated REM sleep behavior disorder (iRBD) is considered a prodromal stage of parkinsonism. Neurodegenerative changes in the substantia nigra pars compacta (SNc) in parkinsonism can be detected using neuromelanin-sensitive MRI. OBJECTIVE: To investigate SNc neuromelanin changes in iRBD patients using fully automatic segmentation. METHODS: We included 47 iRBD patients, 134 early Parkinson's disease (PD) patients and 55 healthy volunteers (HVs) scanned at 3 Tesla. SNc regions-of-interest were delineated automatically using convolutional neural network. SNc volumes, volumes corrected by total intracranial volume, signal-to-noise ratio (SNR) and contrast-to-noise ratio were computed. One-way general linear models (GLM) analysis of covariance (ANCOVA) was conducted while adjusting for age and sex. RESULTS: All SNc measurements differed significantly between the three groups (except SNR in iRBD). Changes in iRBD were intermediate between those in PD and HVs. CONCLUSIONS: Using fully automated SNc segmentation method and neuromelanin-sensitive imaging, iRBD patients showed neurodegenerative changes in the SNc at a lower level than in PD patients. © 2022 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Regional Selectivity of Neuromelanin Changes in the Substantia Nigra in Atypical Parkinsonism.

BACKGROUND: Neurodegeneration in the substantia nigra pars compacta (SNc) in parkinsonian syndromes may affect the nigral territories differently. OBJECTIVE: The objective of this study was to investigate the regional selectivity of neurodegenerative changes in the SNc in patients with Parkinson's disease (PD) and atypical parkinsonism using neuromelanin-sensitive magnetic resonance imaging (MRI). METHODS: A total of 22 healthy controls (HC), 38 patients with PD, 22 patients with progressive supranuclear palsy (PSP), 20 patients with multiple system atrophy (MSA, 13 with the parkinsonian variant, 7 with the cerebellar variant), 7 patients with dementia with Lewy body (DLB), and 4 patients with corticobasal syndrome were analyzed. volume and signal-to-noise ratio (SNR) values of the SNc were derived from neuromelanin-sensitive MRI in the whole SNc. Analysis of signal changes was performed in the sensorimotor, associative, and limbic territories of the SNc. RESULTS: SNc volume and corrected volume were significantly reduced in PD, PSP, and MSA versus HC. Patients with PSP had lower volume, corrected volume, SNR, and contrast-to-noise ratio than HC and patients with PD and MSA. Patients with PSP had greater SNR reduction in the associative region than HC and patients with PD and MSA. Patients with PD had reduced SNR in the sensorimotor territory, unlike patients with PSP. Patients with MSA did not differ from patients with PD. CONCLUSIONS: This study provides the first MRI comparison of the topography of neuromelanin changes in parkinsonism. The spatial pattern of changes differed between PSP and synucleinopathies. These nigral topographical differences are consistent with the topography of the extranigral involvement in parkinsonian syndromes. © 2022 International Parkinson and Movement Disorder Society.

The spatiotemporal changes in dopamine, neuromelanin and iron characterizing Parkinson's disease.

In Parkinson's disease, there is a progressive reduction in striatal dopaminergic function, and loss of neuromelanin-containing dopaminergic neurons and increased iron deposition in the substantia nigra. We tested the hypothesis of a relationship between impairment of the dopaminergic system and changes in the iron metabolism. Based on imaging data of patients with prodromal and early clinical Parkinson's disease, we assessed the spatiotemporal ordering of such changes and relationships in the sensorimotor, associative and limbic territories of the nigrostriatal system. Patients with Parkinson's disease (disease duration < 4 years) or idiopathic REM sleep behaviour disorder (a prodromal form of Parkinson's disease) and healthy controls underwent longitudinal examination (baseline and 2-year follow-up). Neuromelanin and iron sensitive MRI and dopamine transporter single-photon emission tomography were performed to assess nigrostriatal levels of neuromelanin, iron, and dopamine. For all three functional territories of the nigrostriatal system, in the clinically most and least affected hemispheres separately, the following was performed: cross-sectional and longitudinal intergroup difference analysis of striatal dopamine and iron, and nigral neuromelanin and iron; in Parkinson's disease patients, exponential fitting analysis to assess the duration of the prodromal phase and the temporal ordering of changes in dopamine, neuromelanin or iron relative to controls; and voxel-wise correlation analysis to investigate concomitant spatial changes in dopamine-iron, dopamine-neuromelanin and neuromelanin-iron in the substantia nigra pars compacta. The temporal ordering of dopaminergic changes followed the known spatial pattern of progression involving first the sensorimotor, then the associative and limbic striatal and nigral regions. Striatal dopaminergic denervation occurred first followed by abnormal iron metabolism and finally neuromelanin changes in the substantia nigra pars compacta, which followed the same spatial and temporal gradient observed in the striatum but shifted in time. In conclusion, dopaminergic striatal dysfunction and cell loss in the substantia nigra pars compacta are interrelated with increased nigral iron content.

Planning of gamma knife radiosurgery (GKR) for brain arteriovenous malformations using triple magnetic resonance angiography (triple-MRA).

PURPOSE: Intra-arterial Digital Subtraction Angiography (DSA) is the gold standard technique for radiosurgery target delineation in brain Arterio-Venous Malformations (AVMs). This study aims to evaluate whether a combination of three Magnetic Resonance Angiography sequences (triple-MRA) could be used for delineation of brain AVMs for Gamma Knife Radiosurgery (GKR). METHODS: Fifteen patients undergoing DSA for GKR targeting of brain AVMs also underwent triple-MRA: 4D Arterial Spin Labelling based angiography (ASL-MRA), Contrast-Enhanced Time-Resolved MRA (CE-MRA) and High Definition post-contrast Time-Of-Flight angiography (HD-TOF). The arterial phase of the AVM nidus was delineated on triple-MRA by an interventional neuroradiologist and a consultant neurosurgeon (triple-MRA volume). Triple-MRA volumes were compared to AVM targets delineated by the clinical team for delivery of GKR using the current planning paradigm, i.e., stereotactic DSA and volumetric MRI (DSA volume). Difference in size, degree of inclusion (DI) and concordance index (CcI) between DSA and triple-MRA volumes are reported. RESULTS: AVM target volumes delineated on triple-MRA were on average 9.8% smaller than DSA volumes (95%CI:5.6-13.9%; SD:7.14%; p = .003). DI of DSA volume in triple-MRA volume was on average 73.5% (95%CI:71.2-76; range: 65-80%). The mean percentage of triple-MRA volume not included on DSA volume was 18% (95%CI:14.7-21.3; range: 7-30%). CONCLUSION: The technical feasibility of using triple-MRA for visualisation and delineation of brain AVMs for GKR planning has been demonstrated. Tighter and more precise delineation of AVM target volumes could be achieved by using triple-MRA for radiosurgery targeting. However, further research is required to ascertain the impact this may have in obliteration rates and side effects.

Investigating the effect of flow compensation and quantitative susceptibility mapping method on the accuracy of venous susceptibility measurement.

Quantitative susceptibility mapping (QSM) is a promising non-invasive method for obtaining information relating to oxygen metabolism. However, the optimal acquisition sequence and QSM reconstruction method for reliable venous susceptibility measurements are unknown. Full flow compensation is generally recommended to correct for the influence of venous blood flow, although the effect of flow compensation on the accuracy of venous susceptibility values has not been systematically evaluated. In this study, we investigated the effect of different acquisition sequences, including different flow compensation schemes, and different QSM reconstruction methods on venous susceptibilities. Ten healthy subjects were scanned with five or six distinct QSM sequence designs using monopolar readout gradients and different flow compensation schemes. All data sets were processed using six different QSM pipelines and venous blood susceptibility was evaluated in whole-brain segmentations of the venous vasculature and single veins. The quality of vein segmentations and the accuracy of venous susceptibility values were analyzed and compared between all combinations of sequences and reconstruction methods. The influence of the QSM reconstruction method on average venous susceptibility values was found to be 2.7-11.6 times greater than the influence of the acquisition sequence, including flow compensation. The majority of the investigated QSM reconstruction methods tended to underestimate venous susceptibility values in the vein segmentations that were obtained. In summary, we found that multi-echo gradient-echo acquisition sequences without full flow compensation yielded venous susceptibility values comparable to sequences with full flow compensation. However, the QSM reconstruction method had a great influence on susceptibility values and thus needs to be selected carefully for accurate venous QSM.

Longitudinal Changes in Neuromelanin MRI Signal in Parkinson's Disease: A Progression Marker.

BACKGROUND: Development of reliable and accurate imaging biomarkers of dopaminergic cell neurodegeneration is necessary to facilitate therapeutic drug trials in Parkinson's disease (PD). Neuromelanin-sensitive MRI techniques have been effective in detecting neurodegeneration in the substantia nigra pars compacta (SNpc). The objective of the current study was to investigate longitudinal neuromelanin signal changes in the SNpc in PD patients. METHODS: In this prospective, longitudinal, observational case-control study, we included 140 PD patients and 64 healthy volunteers divided into 2 cohorts. Cohort I included 99 early PD patients (disease duration, 1.5 ± 1.0 years) and 41 healthy volunteers analyzed at baseline (V1), where 79 PD patients and 32 healthy volunteers were rescanned after 2.0 ± 0.2 years of follow-up (V2). Cohort II included 41 progressing PD patients (disease duration, 9.3 ± 3.7 years) and 23 healthy volunteers at V1, where 30 PD patients were rescanned after 2.4 ± 0.5 years of follow-up. Subjects were scanned at 3 T MRI using 3-dimensional T1-weighted and neuromelanin-sensitive imaging. Regions of interest were delineated manually to calculate SN volumes, volumes corrected by total intracranial volume, signal-to-noise ratio, and contrast-to-noise ratio. RESULTS: Results showed (1) significant reduction in volume and volume corrected by total intracranial volume between visits, greater in progressing PD than nonsignificant changes in healthy volunteers; (2) no significant effects of visit for signal intensity (signal-to-noise ratio); (3) significant interaction in volume between group and visit; (4) greater volume corrected by total intracranial volume at baseline in female patients and greater decrease in volume and increase in the contrast-to-noise ratio in progressing female PD patients compared with male patients; and (5) correlations between neuromelanin SN changes and disease severity and duration. CONCLUSIONS: We observed a progressive and measurable decrease in neuromelanin-based SN signal and volume in PD, which might allow a direct noninvasive assessment of progression of SN loss and could represent a target biomarker for disease-modifying treatments. © 2021 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Automated Categorization of Parkinsonian Syndromes Using Magnetic Resonance Imaging in a Clinical Setting.

BACKGROUND: Machine learning algorithms using magnetic resonance imaging (MRI) data can accurately discriminate parkinsonian syndromes. Validation in patients recruited in routine clinical practice is missing. OBJECTIVE: The aim of this study was to assess the accuracy of a machine learning algorithm trained on a research cohort and tested on an independent clinical replication cohort for the categorization of parkinsonian syndromes. METHODS: Three hundred twenty-two subjects, including 94 healthy control subjects, 119 patients with Parkinson's disease (PD), 51 patients with progressive supranuclear palsy (PSP) with Richardson's syndrome, 35 with multiple system atrophy (MSA) of the parkinsonian variant (MSA-P), and 23 with MSA of the cerebellar variant (MSA-C), were recruited. They were divided into a training cohort (n = 179) scanned in a research environment and a replication cohort (n = 143) examined in clinical practice on different MRI systems. Volumes and diffusion tensor imaging (DTI) metrics in 13 brain regions were used as input for a supervised machine learning algorithm. To harmonize data across scanners and reduce scanner-dependent effects, we tested two types of normalizations using patient data or healthy control data. RESULTS: In the replication cohort, high accuracies were achieved using volumetry in the classification of PD-PSP, PD-MSA-C, PSP-MSA-C, and PD-atypical parkinsonism (balanced accuracies: 0.840-0.983, area under the receiver operating characteristic curves: 0.907-0.995). Performances were lower for the classification of PD-MSA-P, MSA-C-MSA-P (balanced accuracies: 0.765-0.784, area under the receiver operating characteristic curve: 0.839-0.871) and PD-PSP-MSA (balanced accuracies: 0.773). Performance using DTI was improved when normalizing by controls, but remained lower than that using volumetry alone or combined with DTI. CONCLUSIONS: A machine learning approach based on volumetry enabled accurate classification of subjects with early-stage parkinsonism, examined on different MRI systems, as part of their clinical assessment. © 2020 International Parkinson and Movement Disorder Society.

Spatiotemporal changes in substantia nigra neuromelanin content in Parkinson's disease.

This study aimed to investigate the spatiotemporal changes in neuromelanin-sensitive MRI signal in the substantia nigra and their relation to clinical scores of disease severity in patients with early or progressing Parkinson's disease and patients with idiopathic rapid eye movement sleep behaviour disorder (iRBD) exempt of Parkinsonian signs compared to healthy control subjects. Longitudinal T1-weighted anatomical and neuromelanin-sensitive MRI was performed in two cohorts, including patients with iRBD, patients with early or progressing Parkinson's disease, and control subjects. Based on the aligned substantia nigra segmentations using a study-specific brain anatomical template, parametric maps of the probability of a voxel belonging to the substantia nigra were calculated for patients with various degrees of disease severity and controls. For each voxel in the substantia nigra, probability map of controls, correlations between signal-to-noise ratios on neuromelanin-sensitive MRI in patients with iRBD and Parkinson's disease and clinical scores of motor disability, cognition and mood/behaviour were calculated. Our results showed that in patients, compared to the healthy control subjects, the volume of the substantia nigra was progressively reduced for increasing disease severity. The neuromelanin signal changes appeared to start in the posterolateral motor areas of the substantia nigra and then progressed to more medial areas of this region. The ratio between the volume of the substantia nigra in patients with Parkinson's disease relative to the controls was best fitted by a mono-exponential decay. Based on this model, the pre-symptomatic phase of the disease started at 5.3 years before disease diagnosis, and 23.1% of the substantia nigra volume was lost at the time of diagnosis, which was in line with previous findings using post-mortem histology of the human substantia nigra and radiotracer studies of the human striatum. Voxel-wise patterns of correlation between neuromelanin-sensitive MRI signal-to-noise ratio and motor, cognitive and mood/behavioural clinical scores were localized in distinct regions of the substantia nigra. This localization reflected the functional organization of the nigrostriatal system observed in histological and electrophysiological studies in non-human primates (motor, cognitive and mood/behavioural domains). In conclusion, neuromelanin-sensitive MRI enabled us to assess voxel-wise modifications of substantia nigra's morphology in vivo in humans, including healthy controls, patients with iRBD and patients with Parkinson's disease, and identify their correlation with nigral function across all motor, cognitive and behavioural domains. This insight could help assess disease progression in drug trials of disease modification.

Investigating the accuracy and precision of TE-dependent versus multi-echo QSM using Laplacian-based methods at 3 T.

PURPOSE: Multi-echo gradient-recalled echo acquisitions for QSM enable optimizing the SNR for several tissue types through multi-echo (TE) combination or investigating temporal variations in the susceptibility (potentially reflecting tissue microstructure) by calculating one QSM image at each TE (TE-dependent QSM). In contrast with multi-echo QSM, applying Laplacian-based methods (LBMs) for phase unwrapping and background field removal to single TEs could introduce nonlinear temporal variations (independent of tissue microstructure) into the measured susceptibility. Here, we aimed to compare the effect of LBMs on the QSM susceptibilities in TE-dependent versus multi-echo QSM. METHODS: TE-dependent recalled echo data simulated in a numerical head phantom and gradient-recalled echo images acquired at 3 T in 10 healthy volunteers. Several QSM pipelines were tested, including four distinct LBMs: sophisticated harmonic artifact reduction for phase data (SHARP), variable-radius sophisticated harmonic artifact reduction for phase data (V-SHARP), Laplacian boundary value background field removal (LBV), and one-step total generalized variation (TGV). Results from distinct pipelines were compared using visual inspection, summary statistics of susceptibility in deep gray matter/white matter/venous regions of interest, and, in the healthy volunteers, regional susceptibility bias analysis and nonparametric tests. RESULTS: Multi-echo versus TE-dependent QSM had higher regional accuracy, especially in high-susceptibility regions and at shorter TEs. Everywhere except in the veins, a processing pipeline incorporating TGV provided the most temporally stable TE-dependent QSM results with an accuracy similar to multi-echo QSM. CONCLUSIONS: For TE-dependent QSM, carefully choosing LBMs can minimize the introduction of LBM-related nonlinear temporal susceptibility variations.

Investigating the oxygenation of brain arteriovenous malformations using quantitative susceptibility mapping.

Brain arteriovenous malformations (AVMs) are congenital vascular anomalies characterized by arteriovenous shunting through a network of coiled and tortuous vessels. Because of this anatomy, the venous drainage of an AVM is hypothesized to contain more oxygenated, arterialized blood than healthy veins. By exploiting the paramagnetic properties of deoxygenated hemoglobin in venous blood using magnetic resonance imaging (MRI) quantitative susceptibility mapping (QSM), we aimed to explore venous density and oxygen saturation (SvO2) in patients with a brain AVM. We considered three groups of subjects: patients with a brain AVM before treatment using gamma knife radiosurgery (GKR); patients three or more years post-GKR treatment; and healthy volunteers. First, we investigated the appearance of AVMs on QSM images. Then, we investigated whether QSM could detect increased SvO2 in the veins draining the malformations. In patients before GKR, venous density, but not SvO2, was significantly larger in the hemisphere containing the AVM compared to the contralateral hemisphere (p = 0.03). Such asymmetry was not observed in patients after GKR or in healthy volunteers. Moreover, in all patients before GKR, the vein immediately draining the AVM nidus had a higher SvO2 than healthy veins. Therefore, QSM can be used to detect SvO2 alterations in brain AVMs. However, since factors such as flow-induced signal dephasing or the presence of hemosiderin deposits also strongly affect QSM image contrast, AVM vein segmentation must be performed based on alternative MRI acquisitions, e.g., time of flight magnetic resonance angiography or T1-weighted images. This is the first study to show, non-invasively, that AVM draining veins have a significantly larger SvO2 than healthy veins, which is a finding congruent with arteriovenous shunting.

A flow-diffusion model of oxygen transport for quantitative mapping of cerebral metabolic rate of oxygen (CMRO2) with single gas calibrated fMRI.

One promising approach for mapping CMRO2 is dual-calibrated functional MRI (dc-fMRI). This method exploits the Fick Principle to combine estimates of CBF from ASL, and OEF derived from BOLD-ASL measurements during arterial O2 and CO2 modulations. Multiple gas modulations are required to decouple OEF and deoxyhemoglobin-sensitive blood volume. We propose an alternative single gas calibrated fMRI framework, integrating a model of oxygen transport, that links blood volume and CBF to OEF and creates a mapping between the maximum BOLD signal, CBF and OEF (and CMRO2). Simulations demonstrated the method's viability within physiological ranges of mitochondrial oxygen pressure, PmO2, and mean capillary transit time. A dc-fMRI experiment, performed on 20 healthy subjects using O2 and CO2 challenges, was used to validate the approach. The validation conveyed expected estimates of model parameters (e.g., low PmO2), with spatially uniform OEF maps (grey matter, GM, OEF spatial standard deviation ≈ 0.13). GM OEF estimates obtained with hypercapnia calibrated fMRI correlated with dc-fMRI (r = 0.65, p = 2·10-3). For 12 subjects, OEF measured with dc-fMRI and the single gas calibration method were correlated with whole-brain OEF derived from phase measures in the superior sagittal sinus (r = 0.58, p = 0.048; r = 0.64, p = 0.025 respectively). Simplified calibrated fMRI using hypercapnia holds promise for clinical application.

Cerebrovascular reactivity in multiple sclerosis is restored with reduced inflammation during immunomodulation.

Cerebrovascular reactivity (CVR) reflects the capacity of the brain's vasculature to increase blood flow following a vasodilatory stimulus. Reactivity is an essential property of the brain's blood vessels that maintains nutrient supplies in the face of changing demand. In Multiple Sclerosis (MS), CVR may be diminished with brain inflammation and this may contribute to neurodegeneration. We test the hypothesis that CVR is altered with MS neuroinflammation and that it is restored when inflammation is reduced. Using a breath-hold task during functional Magnetic Resonance Imaging (MRI), we mapped grey matter and white matter CVRs (CVRGM and CVRWM, respectively) in 23 young MS patients, eligible for disease modifying therapy, before and during Interferon beta treatment. Inflammatory activity was inferred from the presence of Gadolinium enhancing lesions at MRI. Eighteen age and gender-matched healthy controls (HC) were also assessed. Enhancing lesions were observed in 12 patients at the start of the study and in 3 patients during treatment. Patients had lower pre-treatment CVRGM (p = 0.04) and CVRWM (p = 0.02) compared to HC. In patients, a lower pre-treatment CVRGM was associated with a lower GM volume (r = 0.60, p = 0.003). On-treatment, there was an increase in CVRGM (p = 0.02) and CVRWM (p = 0.03) that negatively correlated with pre-treatment CVR (GM: r = - 0.58, p = 0.005; WM: r = - 0.60, p = 0.003). CVR increased when enhancing lesions reduced in number (GM: r = - 0.48, p = 0.02, WM: r = - 0.62, p = 0.003). Resolution of inflammation may restore altered cerebrovascular function limiting neurodegeneration in MS. Imaging of cerebrovascular function may thereby inform tissue physiology and improve treatment monitoring.