High-resolution T1 mapping of the brain at 3T with driven equilibrium single pulse observation of T1 with high-speed incorporation of RF field inhomogeneities (DESPOT1-HIFI)

Presence of slowly expanding lesions in multiple sclerosis predicts progressive demyelination within lesions and normal-appearing tissue over time

BACKGROUND

Multiple sclerosis (MS) slowly expanding lesions (SELs) are defined on magnetic resonance imaging (MRI) as contiguous regions of pre-existing focal non-contrast-enhancing T2 lesions with constant and concentric local expansion on conventional T1-weighted and T2-weighted images. SELs are associated with an increased risk of disability progression.

METHODS

Myelin-related changes detected using myelin water fraction (MWF) and magnetisation transfer ratio (MTR) in SELs and T2 lesions were measured over 192 weeks in participants with relapsing MS.

RESULTS

In participants with SELs (SEL+), SELs (MWF: 0.12 ± 0.03, MTR: 33.1 ± 3.6 pu) showed reduced myelin measures at baseline compared to T2 lesions (MWF: 0.13 ± 0.02, MTR: 35.1 ± 2.4 pu). In participants without SELs (SEL-), T2 lesions had higher myelin measures (MWF: 0.15 ± 0.02, MTR: 36.2 ± 2.0 pu) compared to T2 lesions in SEL+. Over 4 years, only SELs showed decreases in MWF (-11.4%). The percentage of abnormal voxels within normal-appearing white matter was higher in SEL+ and increased over time (SEL+ MWF Week 0: 0.56%, Week 192: 0.98%; SEL- MWF Week 0: 0.13%, Week 192: 0.25%).

CONCLUSION

Our results indicate progressive focal and global demyelination in SEL+ participants and that the presence of SELs might be a biomarker for participants with ongoing diffuse or smouldering inflammation within the whole brain.

Blood biomarkers of vascular dysfunction in small vessel disease progression: Insights from a longitudinal neuroimaging study

INTRODUCTION

This study explored the relationship between blood biomarkers of cerebrovascular function and small vessel disease (SVD) neuroimaging markers and cognitive outcomes in highly-phenotyped participants.

METHODS

We conducted cross-sectional and 1-year longitudinal analyses on 181 patients with mild ischemic stroke, enriched for SVD features. We examined relationships between a panel of 13 blood biomarkers and magnetic resonance imaging (MRI) markers of SVD (structural lesions, diffusion-weighted imaging [DWI]-positive lesions, blood-brain barrier (BBB) permeability, and cerebrovascular reactivity (CVR), and cognition.

RESULTS

In linear mixed models, vascular endothelial growth factor was significantly associated with incident DWI-positive lesions over 1 year. Intercellular adhesion molecule-1 was linked with lower CVR while platelet-derived growth factor-subunit B and Endothelin-1 were associated with higher CVR. Platelet-Selectin levels were associated with mild cognitive impairment at 1 year.

DISCUSSION

Our results support the role of endothelial and pericyte dysfunction in SVD burden and progression and suggest that specific biomarkers relate to distinct SVD manifestations.

HIGHLIGHTS

Small vessel disease (SVD) lacks specific or predictive biomarker signatures. Vascular endothelial growth factor levels were linked to incident lesions detected over 1 year. Circulating intercellular adhesion molecule-1 related to lower cerebrovascular reactivity. Platelet-selectin levels were associated with mild cognitive impairment longitudinally. These findings could help stratify patients at high-risk of rapid-progression SVD.

A Water Relaxation Atlas for Age- and Region-Specific Metabolite Concentration Correction at 3 T

Metabolite concentration estimates from magnetic resonance spectroscopy (MRS) are typically quantified using water referencing, correcting for relaxation-time differences between metabolites and water. One common approach is to correct the reference signal for differential relaxation within three tissue compartments (gray matter, white matter, and cerebrospinal fluid) using fixed literature values. However, water relaxation times (T1 and T2) vary between brain locations and with age. MRS studies, even those measuring metabolite levels across the lifespan, often ignore these effects, because of a lack of reference data. The purpose of this study is to develop a water relaxometry atlas and to integrate location- and age-appropriate relaxation values into the MRS analysis workflow. One hundred one volunteers (51 men, 50 women; ~10 male, and 10 female participants per decade from the 20s to 60s) were recruited. T1-weighted MPRAGE images ((1-mm)3 isotropic resolution) were acquired. Whole-brain water T1 and T2 measurements were made with DESPOT ((1.4 mm)3 isotropic resolution) at 3T. T1 and T2 maps were registered to the JHU MNI-SS/EVE atlas using affine and LDDMM transformation. The atlas's 268 parcels were reduced to 130 by combining homologous parcels. Mean T1 and T2 values were calculated for each parcel in each subject. Linear models of T1 and T2 as functions of age were computed, using age - 30 as the predictor. Reference atlases of "age-30-intercept" and age-slope for T1 and T2 were generated. The atlas-based workflow was integrated into Osprey, which co-registers MRS voxels to the atlas and calculates location- and age-appropriate water relaxation parameters for quantification. The water relaxation aging atlas revealed significant regional and tissue differences in water relaxation behavior across adulthood. Using location- and subject-appropriate reference values in the MRS analysis workflow removes a current methodological limitation and is expected to reduce quantification biases associated with water-referenced tissue correction, especially for studies of aging.

Heterogeneity and Penumbra of White Matter Hyperintensities in Small Vessel Diseases Determined by Quantitative MRI

BACKGROUND

White matter hyperintensities (WMHs) are established structural imaging markers of cerebral small vessel disease. The pathophysiologic condition of brain tissue varies over the core, the vicinity, and the subtypes of WMH and cannot be interpreted from conventional magnetic resonance imaging. We aim to improve our pathophysiologic understanding of WMHs and the adjacently injured normal-appearing white matter in terms of microstructural and microvascular alterations using quantitative magnetic resonance imaging in patients with sporadic and genetic cerebral small vessel disease.

METHODS

Structural T2-weighted imaging, multishell diffusion imaging, and dynamic contrast-enhanced magnetic resonance imaging were performed at 3T in 44 participants with sporadic cerebral small vessel disease and 32 participants with monogenic cerebral small vessel disease (cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy; 59±12 years, 41 males) between June 2017 and May 2020 as part of the prospective, multicenter (Edinburgh, the United Kingdom; Maastricht, the Netherlands; and Munich, Germany), observational INVESTIGATE-SVDs study (Imaging Neurovascular, Endothelial and Structural Integrity in Preparation to Treat Small Vessel Diseases). The mean diffusivity, free water content, and perfusion (all derived from multishell diffusion imaging), as well as the blood-brain barrier leakage and plasma volume fraction (derived from dynamic contrast-enhanced magnetic resonance imaging), were compared between deep and periventricular WMH types using paired t tests. Additional spatial analyses were performed inside and outside the WMH types to determine the internal heterogeneity and the extent of the penumbras, that is, adjacent white matter at risk for conversion to WMH.

RESULTS

Periventricular WMH had higher mean diffusivity, higher free water content, and more plasma volume compared with deep WMH (P<0.001, P=0.01, and P<0.001, respectively). No differences were observed in perfusion (P=0.94) and blood-brain barrier leakage (P=0.65) between periventricular and deep WMHs. The spatial analyses inside WMH and the adjacent white matter revealed a gradual gradient in white matter microstructure, free water content, perfusion, and plasma volume but not in blood-brain barrier leakage.

CONCLUSIONS

We showed different pathophysiological heterogeneity of the 2 WMH types. Periventricular WMHs display more severe damage and fluid accumulation compared with deep WMH, whereas deep WMHs reflect stronger hypoperfusion in the lesion's core.

REGISTRATION

URL: https://www.isrctn.com; Unique identifier: ISRCTN10514229.

T1 Thermometry for Deep Brain Stimulation Applications: A Comparison between Rapid Gradient Echo Sequences

Background

T1 thermometry is considered a straight method for the safety monitoring of patients with deep brain stimulation (DBS) electrodes against radiofrequency-induced heating during Magnetic Resonance Imaging (MRI), requiring different sequences and methods.

Objective

This study aimed to compare two T1 thermometry methods and two low specific absorption rate (SAR) imaging sequences in terms of the output image quality.

Material and Methods

In this experimental study, a gel phantom was prepared, resembling the brain tissue properties with a copper wire inside. Two types of rapid gradient echo sequences, namely radiofrequency-spoiled and balanced steady-state free precession (bSSFP) sequences, were used. T1 thermometry was performed by either T1-weighted images with a high SAR sequence to increase heating around the wire or T1 mapping methods.

Results

The balanced steady-state free precession (bSSFP) sequence provided higher image quality in terms of spatial resolution (1×1×1.5 mm3 compared with 1×1×3 mm3) at a shorter acquisition time. The susceptibility artifact was also less pronounced for the bSSFP sequence compared with the radiofrequency-spoiled sequence. A temperature increase, of up to 8 ℃, was estimated using a high SAR sequence. The estimated change in temperature was reduced when using the T1 mapping method.

Conclusion

Heating induced during MRI of implanted electrodes could be estimated using high-resolution T1 maps obtained from inversion recovery bSSFP sequence. Such a method gives a direct estimation of heating during the imaging sequence, which is highly desirable for safe MRI of DBS patients.

Black-Blood Magnetization Prepared 2 Rapid Acquisition Gradient Echoes: A Fast and Three-Dimensional MR Black-Blood T1 Mapping Technique for Quantitative Assessment of Atherosclerosis and Venous Thrombosis

BACKGROUND

Blood flow signals may be a confounder in quantifying T1 values of plaque or thrombus and how to realize black-blood T1 mapping remains a challenge task.

PURPOSE

To develop a fast and three-dimensional black-blood T1 mapping technique for quantitative assessment of atherosclerosis and venous thrombosis.

STUDY TYPE

Sequence development and optimization via phantoms and volunteers as well as pilot prospective.

PHANTOM AND SUBJECTS

Numerical simulations, a standard phantom, 8 healthy volunteers (mean age, 22 ± 1 years; 5 males), and 19 patients (mean age, 57 ± 14 years; 13 males) with atherosclerosis or venous thrombosis.

FIELD STRENGTH/SEQUENCE

3T/inversion recovery spin-echo sequence (IR-SE), magnetization prepared 2 rapid acquisition gradient echoes (MP2RAGE), and black-blood prepared MP2RAGE (BB-MP2RAGE).

ASSESSMENT

The black-blood preparation (i.e., delay alternating with nutation for tailored excitation, DANTE) was incorporated into MP2RAGE for black-blood T1 mapping. The BB-MP2RAGE was optimized numerically based on the Bloch equation, and then the phantom study was performed to verify the accuracy of T1 mapping by BB-MP2RAGE against IR-SE and MP2RAGE. Preliminary clinical validation was prospectively performed to assess the flow suppression effect and its potential application in plaque and thrombosis identification.

STATISTICAL TESTS

Pearson correlation test, Bland-Altman analysis, paired t-test, and intraclass correlation coefficient. A P value <0.05 indicates a statistically significant difference.

RESULTS

Phantom experiments showed comparable accuracy of T1 maps by BB-MP2RAGE with IR-SE and MP2RAGE (all r2 > 0.99); Compared to MP2RAGE, BB-MP2RAGE effectively nulled the blood flow signals, and had a significant improvement in contrast-to-noise ratio between static tissue and blood (250.5 ± 66.6 vs. 91.9 ± 35.9). BB-MP2RAGE can quantify plaque or thrombus T1 relaxation time with blood flow signal suppression.

DATA CONCLUSION

Accurate T1 mapping with sufficient blood flow suppression was achieved by BB-MP2RAGE. BB-MP2RAGE has the potential to quantitatively characterize atherosclerosis and venous thrombosis.

LEVEL OF EVIDENCE

1 TECHNICAL EFFICACY: Stage 1.

Atlas-based assessment of hypomyelination: Quantitative MRI in Pelizaeus-Merzbacher disease

Pelizaeus-Merzbacher disease (PMD) is a rare childhood hypomyelinating leukodystrophy. Quantification of the pronounced myelin deficit and delineation of subtle myelination processes are of high clinical interest. Quantitative magnetic resonance imaging (qMRI) techniques can provide in vivo insights into myelination status, its spatial distribution, and dynamics during brain maturation. They may serve as potential biomarkers to assess the efficacy of myelin-modulating therapies. However, registration techniques for image quantification and statistical comparison of affected pediatric brains, especially those of low or deviant image tissue contrast, with healthy controls are not yet established. This study aimed first to develop and compare postprocessing pipelines for atlas-based quantification of qMRI data in pediatric patients with PMD and evaluate their registration accuracy. Second, to apply an optimized pipeline to investigate spatial myelin deficiency using myelin water imaging (MWI) data from patients with PMD and healthy controls. This retrospective single-center study included five patients with PMD (mean age, 6 years ± 3.8) who underwent conventional brain MRI and diffusion tensor imaging (DTI), with MWI data available for a subset of patients. Three methods of registering PMD images to a pediatric template were investigated. These were based on (a) T1-weighted (T1w) images, (b) fractional anisotropy (FA) maps, and (c) a combination of T1w, T2-weighted, and FA images in a multimodal approach. Registration accuracy was determined by visual inspection and calculated using the structural similarity index method (SSIM). SSIM values for the registration approaches were compared using a t test. Myelin water fraction (MWF) was quantified from MWI data as an assessment of relative myelination. Mean MWF was obtained from two PMDs (mean age, 3.1 years ± 0.3) within four major white matter (WM) pathways of a pediatric atlas and compared to seven healthy controls (mean age, 3 years ± 0.2) using a Mann-Whitney U test. Our results show that visual registration accuracy estimation and computed SSIM were highest for FA-based registration, followed by multimodal, and T1w-based registration (SSIMFA = 0.67 ± 0.04 vs. SSIMmultimodal = 0.60 ± 0.03 vs. SSIMT1 = 0.40 ± 0.14). Mean MWF of patients with PMD within the WM pathways was significantly lower than in healthy controls MWFPMD = 0.0267 ± 0.021 vs. MWFcontrols = 0.1299 ± 0.039. Specifically, MWF was measurable in brain structures known to be myelinated at birth (brainstem) or postnatally (projection fibers) but was scarcely detectable in other brain regions (commissural and association fibers). Taken together, our results indicate that registration accuracy was highest with an FA-based registration pipeline, providing an alternative to conventional T1w-based registration approaches in the case of hypomyelinating leukodystrophies missing normative intrinsic tissue contrasts. The applied atlas-based analysis of MWF data revealed that the extent of spatial myelin deficiency in patients with PMD was most pronounced in commissural and association and to a lesser degree in brainstem and projection pathways.

Comparison of test-retest reproducibility of DESPOT and 3D-QALAS for water T 1 and T 2 mapping

Purpose

Relaxometry, specifically T 1 and T 2 mapping, has become an essential technique for assessing the properties of biological tissues related to various physiological and pathological conditions. Many techniques are being used to estimate T 1 and T 2 relaxation times, ranging from the traditional inversion or saturation recovery and spin-echo sequences to more advanced methods. Choosing the appropriate method for a specific application is critical since the precision and accuracy of T 1 and T 2 measurements are influenced by a variety of factors including the pulse sequence and its parameters, the inherent properties of the tissue being examined, the MRI hardware, and the image reconstruction. The aim of this study is to evaluate and compare the test-retest reproducibility of two advanced MRI relaxometry techniques (Driven Equilibrium Single Pulse Observation of T 1 and T 2, DESPOT, and 3D Quantification using an interleaved Look-Locker acquisition Sequence with a T 2 preparation pulse, QALAS), for T 1 and T 2 mapping in a healthy volunteer cohort.

Methods

10 healthy volunteers underwent brain MRI at 1.3 mm3 isotropic resolution, acquiring DESPOT and QALAS data (~11.8 and ~5 minutes duration, including field maps, respectively), test-retest with subject repositioning, on a 3.0 Tesla Philips Ingenia Elition scanner. To reconstruct the T 1 and T 2 maps, we used an equation-based algorithm for DESPOT and a dictionary-based algorithm that incorporates inversion efficiency and B 1 -field inhomogeneity for QALAS. The test-retest reproducibility was assessed using the coefficient of variation (CoV), intraclass correlation coefficient (ICC) and Bland-Altman plots.

Results

Our results indicate that both the DESPOT and QALAS techniques demonstrate good levels of test-retest reproducibility for T 1 and T 2 mapping across the brain. Higher whole-brain voxel-to-voxel ICCs are observed in QALAS for T 1 (0.84 ± 0.039) and in DESPOT for T 2 (0.897 ± 0.029). The Bland-Altman plots show smaller bias and variability of T 1 estimates for QALAS (mean of -0.02 s, and upper and lower limits of -0.14 and 0.11 s, 95% CI) than for DESPOT (mean of -0.02 s, and limits of -0.31 and 0.27 s). QALAS also showed less variability (mean 1.08 ms, limits -1.88 to 4.04 ms) for T 2 compared to DESPOT (mean of 2.56 ms, and limits -17.29 to 22.41 ms). The within-subject CoVs for QALAS range from 0.6% (T 2 in CSF) to 5.8% (T 2 in GM), while for DESPOT they range from 2.1% (T 2 in CSF) to 6.7% (T 2 in GM). The between-subject CoVs for QALAS range from 2.5% (T 2 in GM) to 12% (T 2 in CSF), and for DESPOT they range from 3.7% (T 2 in WM) to 9.3% (T 2 in CSF).

Conclusion

Overall, QALAS demonstrated better reproducibility for T 1 and T 2 measurements than DESPOT, in addition to reduced acquisition time.

Intra-scan RF power amplifier drift correction

PURPOSE

The drift in radiofrequency (RF) power amplifiers (RFPAs) is assessed and several contributing factors are investigated. Two approaches for prospective correction of drift are proposed and their effectiveness is evaluated.

METHODS

RFPA drift assessment encompasses both intra-pulse and inter-pulse drift analyses. Scan protocols with varying flip angle (FA), RF length, and pulse repetition time (TR) are used to gauge the influence of these parameters on drift. Directional couplers (DICOs) monitor the forward waveforms of the RFPA outputs. DICOs data is stored for evaluation, allowing calculation of correction factors to adjust RFPAs' transmit voltage. Two correction methods, predictive and run-time, are employed: predictive correction necessitates a calibration scan, while run-time correction calculates factors during the ongoing scan.

RESULTS

RFPA drift is indeed influenced by the RF duty-cycle, and in the cases examined with a maximum duty-cycle of 66%, the potential drift is approximately 41% or 15%, depending on the specific RFPA revision. Notably, in low transmit voltage scenarios, FA has minimal impact on RFPA drift. The application of predictive and run-time drift correction techniques effectively reduces the average drift from 10.0% to less than 1%, resulting in enhanced MR signal stability.

CONCLUSION

Utilizing DICO recordings and implementing a feedback mechanism enable the prospective correction of RFPA drift. Having a calibration scan, predictive correction can be utilized with fewer complexity; for enhanced performance, a run-time approach can be employed.

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.

Patient-Specific Vascular Flow Phantom for MRI- and Doppler Ultrasound Imaging

OBJECTIVE

Intraoperative Doppler ultrasound imaging of human brain vasculature is an emerging neuro-imaging modality that offers vascular brain mapping with unprecedented spatiotemporal resolution. At present, however, access to the human brain using Doppler Ultrasound is only possible in this intraoperative context, posing a significant challenge for validation of imaging techniques. This challenge necessitates the development of realistic flow phantoms outside of the neurosurgical operating room as external platforms for testing hardware and software. An ideal ultrasound flow phantom should provide reference-like values in standardized topologies such as a slanted pipe, and allow for measurements in structures closely resembling vascular morphology of actual patients. Additionally, the phantom should be compatible with other clinical cerebrovascular imaging modalities. To meet these criteria, we developed and validated a versatile, multimodal MRI- and ultrasound Doppler phantom.

METHODS

Our approach incorporates the latest advancements in phantom research using tissue-mimicking material and 3D-printing with water-soluble resin to create wall-less patient-specific lumens, compatible for ultrasound and MRI.

RESULTS

We successfully produced three distinct phantoms: a slanted pipe, a y-shape phantom representing a bifurcating vessel and an arteriovenous malformation (AVM) derived from clinical Digital Subtraction Angiography (DSA)-data of the brain. We present 3D ultrafast power Doppler imaging results from these phantoms, demonstrating their ability to mimic complex flow patterns as observed in the human brain. Furthermore, we showcase the compatibility of our phantom with Magnetic Resonance Imaging (MRI).

CONCLUSION

We developed an MRI- and Doppler Ultrasound-compatible flow-phantom using customizable, water-soluble resin prints ranging from geometrical forms to patient-specific vasculature.

Contrast-agent-based perfusion MRI code repository and testing framework: ISMRM Open Science Initiative for Perfusion Imaging (OSIPI)

PURPOSE

Software has a substantial impact on quantitative perfusion MRI values. The lack of generally accepted implementations, code sharing and transparent testing reduces reproducibility, hindering the use of perfusion MRI in clinical trials. To address these issues, the ISMRM Open Science Initiative for Perfusion Imaging (OSIPI) aimed to establish a community-led, centralized repository for sharing open-source code for processing contrast-based perfusion imaging, incorporating an open-source testing framework.

METHODS

A repository was established on the OSIPI GitHub website. Python was chosen as the target software language. Calls for code contributions were made to OSIPI members, the ISMRM Perfusion Study Group, and publicly via OSIPI websites. An automated unit-testing framework was implemented to evaluate the output of code contributions, including visual representation of the results.

RESULTS

The repository hosts 86 implementations of perfusion processing steps contributed by 12 individuals or teams. These cover all core aspects of DCE- and DSC-MRI processing, including multiple implementations of the same functionality. Tests were developed for 52 implementations, covering five analysis steps. For T1 mapping, signal-to-concentration conversion and population AIF functions, different implementations resulted in near-identical output values. For the five pharmacokinetic models tested (Tofts, extended Tofts-Kety, Patlak, two-compartment exchange, and two-compartment uptake), differences in output parameters were observed between contributions.

CONCLUSIONS

The OSIPI DCE-DSC code repository represents a novel community-led model for code sharing and testing. The repository facilitates the re-use of existing code and the benchmarking of new code, promoting enhanced reproducibility in quantitative perfusion imaging.

Accelerated liver water T1 mapping using single-shot continuous inversion-recovery spiral imaging

PURPOSE

Liver T1 mapping techniques typically require long breath holds or long scan time in free-breathing, need correction for B 1 + inhomogeneities and process composite (water and fat) signals. The purpose of this work is to accelerate the multi-slice acquisition of liver water selective T1 (wT1) mapping in a single breath hold, improving the k-space sampling efficiency.

METHODS

The proposed continuous inversion-recovery (IR) Look-Locker methodology combines a single-shot gradient echo spiral readout, Dixon processing and a dictionary-based analysis for liver wT1 mapping at 3 T. The sequence parameters were adapted to obtain short scan times. The influence of fat, B 1 + inhomogeneities and TE on the estimation of T1 was first assessed using simulations. The proposed method was then validated in a phantom and in 10 volunteers, comparing it with MRS and the modified Look-Locker inversion-recovery (MOLLI) method. Finally, the clinical feasibility was investigated by comparing wT1 maps with clinical scans in nine patients.

RESULTS

The phantom results are in good agreement with MRS. The proposed method encodes the IR-curve for the liver wT1 estimation, is minimally sensitive to B 1 + inhomogeneities and acquires one slice in 1.2 s. The volunteer results confirmed the multi-slice capability of the proposed method, acquiring nine slices in a breath hold of 11 s. The present work shows robustness to B 1 + inhomogeneities (wT 1 , No B 1 + = 1.07 wT 1 , B 1 + - 45.63 , R 2 = 0.99 ) , good repeatability (wT 1 , 2 ° = 1 . 0 wT 1 , 1 ° - 2.14 , R 2 = 0.96 ) and is in better agreement with MRS (wT 1 = 0.92 wT 1 MRS + 103.28 , R 2 = 0.38 ) than is MOLLI (wT 1 MOLLI = 0.76 wT 1 MRS + 254.43 , R 2 = 0.44 ) . The wT1 maps in patients captured diverse lesions, thus showing their clinical feasibility.

CONCLUSION

A single-shot spiral acquisition can be combined with a continuous IR Look-Locker method to perform rapid repeatable multi-slice liver water T1 mapping at a rate of 1.2 s per slice without a B 1 + map. The proposed method is suitable for nine-slice liver clinical applications acquired in a single breath hold of 11 s.

Phosphorus-31: A table-top method for 3D B1-field amplitude and phase measurements

A novel method of high-spatial-resolution, 3D B1-field distribution measurements is presented. The method is independent of the MR-scanner, and it allows for automated acquisitions of complete maps of all magnetic field vector components for both proton and heteronuclear MR coils of arbitrary geometrical shapes. The advantage of the method proposed here, compared with methods based on measurements with an MR-scanner, is that a complete image of both receive and transmit B1-fields, including the phase of the B1-field, can be acquired. The B1 field maps obtained in this manner can be used for absolute quantification of metabolites in MRS experiments, as well as for intensity compensations in imaging experiments, both of which are important concepts in biological and medical MR applications. Another use might be in coil development and testing. A comparison with B1 field magnitude maps obtained with an MR-scanner was included to validate the accuracy of the proposed method.

Magnetic Resonance Imaging Tissue Signatures Associated With White Matter Changes Due to Sporadic Cerebral Small Vessel Disease Indicate That White Matter Hyperintensities Can Regress

BACKGROUND

White matter hyperintensities (WMHs) might regress and progress contemporaneously, but we know little about underlying mechanisms. We examined WMH change and underlying quantitative magnetic resonance imaging tissue measures over 1 year in patients with minor ischemic stroke with sporadic cerebral small vessel disease.

METHODS AND RESULTS

We defined areas of stable normal-appearing white matter, stable WMHs, progressing and regressing WMHs based on baseline and 1-year brain magnetic resonance imaging. In these areas we assessed tissue characteristics with quantitative T1, fractional anisotropy (FA), mean diffusivity (MD), and neurite orientation dispersion and density imaging (baseline only). We compared tissue signatures cross-sectionally between areas, and longitudinally within each area. WMH change masks were available for N=197. Participants' mean age was 65.61 years (SD, 11.10), 59% had a lacunar infarct, and 68% were men. FA and MD were available for N=195, quantitative T1 for N=182, and neurite orientation dispersion and density imaging for N=174. Cross-sectionally, all 4 tissue classes differed for FA, MD, T1, and Neurite Density Index. Longitudinally, in regressing WMHs, FA increased with little change in MD and T1 (difference estimate, 0.011 [95% CI, 0.006-0.017]; -0.002 [95% CI, -0.008 to 0.003] and -0.003 [95% CI, -0.009 to 0.004]); in progressing and stable WMHs, FA decreased (-0.022 [95% CI, -0.027 to -0.017] and -0.009 [95% CI, -0.011 to -0.006]), whereas MD and T1 increased (progressing WMHs, 0.057 [95% CI, 0.050-0.063], 0.058 [95% CI, 0.050 -0.066]; stable WMHs, 0.054 [95% CI, 0.045-0.063], 0.049 [95% CI, 0.039-0.058]); and in stable normal-appearing white matter, MD increased (0.004 [95% CI, 0.003-0.005]), whereas FA and T1 slightly decreased and increased (-0.002 [95% CI, -0.004 to -0.000] and 0.005 [95% CI, 0.001-0.009]).

CONCLUSIONS

Quantitative magnetic resonance imaging shows that WMHs that regress have less abnormal microstructure at baseline than stable WMHs and follow trajectories indicating tissue improvement compared with stable and progressing WMHs.

Lower fractional anisotropy without evidence for neuro-inflammation in patients with early-phase schizophrenia spectrum disorders

Various lines of research suggest immune dysregulation as a potential therapeutic target for negative and cognitive symptoms in schizophrenia spectrum disorders (SSD). Immune dysregulation would lead to higher extracellular free-water (EFW) in cerebral white matter (WM), which may partially underlie the frequently reported lower fractional anisotropy (FA) in SSD. We aim to investigate differences in EFW concentrations - a presumed proxy for neuro-inflammation - between early-phase SSD patients (n = 55) and healthy controls (HC; n = 37), and to explore immunological and cognitive correlates. To increase specificity for EFW, we study several complementary magnetic resonance imaging contrasts that are sensitive to EFW. FA, mean diffusivity (MD), magnetization transfer ratio (MTR), myelin water fraction (MWF) and quantitative T1 and T2 were calculated from diffusion-weighted imaging (DWI), magnetization transfer imaging (MTI) and multicomponent driven equilibrium single-pulse observation of T1/T2 (mcDESPOT). For each measure, WM skeletons were constructed with tract-based spatial statistics. Multivariate SSD-HC comparisons with WM skeletons and their average values (i.e. global WM) were not statistically significant. In voxel-wise analyses, FA was significantly lower in SSD in the genu of the corpus callosum and in the left superior longitudinal fasciculus (p < 0.04). Global WM measures did not correlate with immunological markers (i.e. IL1-RA, IL-6, IL-8, IL-10 and CRP) or cognition in HC and SSD after corrections for multiple comparisons. We confirmed lower FA in early-phase SSD patients. However, nonFA measures did not provide additional evidence for immune dysregulation or for higher EFW as the primary mechanism underlying the reported lower FA values in SSD.

B1 inhomogeneity-corrected T1 mapping and quantitative magnetization transfer imaging via simultaneously estimating Bloch-Siegert shift and magnetization transfer effects

PURPOSE

To introduce a method of inducing Bloch-Siegert shift and magnetization Transfer Simultaneously (BTS) and demonstrate its utilization for measuring binary spin-bath model parameters free pool spin-lattice relaxation (T 1 F $$ {T}_1^{\mathrm{F}} $$ ), macromolecular fraction (f $$ f $$ ), magnetization exchange rate (k F $$ {k}_{\mathrm{F}} $$ ) and local transmit field (B 1 + $$ {B}_1^{+} $$ ).

THEORY AND METHODS

Bloch-Siegert shift and magnetization transfer is simultaneously induced through the application of off-resonance irradiation in between excitation and acquisition of an RF-spoiled gradient-echo scheme. Applying the binary spin-bath model, an analytical signal equation is derived and verified through Bloch simulations. Monte Carlo simulations were performed to analyze the method's performance. The estimation of the binary spin-bath parameters withB 1 + $$ {B}_1^{+} $$ compensation was further investigated through experiments, both ex vivo and in vivo.

RESULTS

Comparing BTS with existing methods, simulations showed that existing methods can significantly biasT 1 $$ {T}_1 $$ estimation when not accounting for transmitB 1 $$ {B}_1 $$ heterogeneity and MT effects that are present. Phantom experiments further showed that the degree of this bias increases with increasing macromolecular proton fraction. Multi-parameter fit results from an in vivo brain study generated values in agreement with previous literature. Based on these studies, we confirmed that BTS is a robust method for estimating the binary spin-bath parameters in macromolecule-rich environments, even in the presence ofB 1 + $$ {B}_1^{+} $$ inhomogeneity.

CONCLUSION

A method of estimating Bloch-Siegert shift and magnetization transfer effect has been developed and validated. Both simulations and experiments confirmed that BTS can estimate spin-bath parameters (T 1 F $$ {T}_1^{\mathrm{F}} $$ ,f $$ f $$ ,k F $$ {k}_{\mathrm{F}} $$ ) that are free fromB 1 + $$ {B}_1^{+} $$ bias.

Cortical Structure Differences in Relation to Age, Sexual Attractions, and Gender Dysphoria in Adolescents: An Examination of Mean Diffusivity and T1 Relaxation Time

Recent research found that the combination of masculine gender identity and gynephilia was associated with cortical T1 relaxation time, which is considered to reflect gray matter density. We hypothesized that mean diffusivity (MD), a diffusion tensor imaging metric that reflects the degree to which water movement is free versus constrained, in combination with T1 relaxation time would provide further insight regarding cortical tissue characteristics. MD and T1 relaxation time were measured in 76 cortical regions in 15 adolescents assigned female at birth who experience gender dysphoria (GD AFAB) and were not receiving hormone therapy, 17 cisgender girls, and 14 cisgender boys (ages 12-17 years). Sexual orientation was represented by the degree of androphilia-gynephilia and the strength of sexual attraction. In multivariate analyses, cortical T1 relaxation time showed a weak but statistically significant positive association with MD across the cortex, suggesting that macromolecule-rich cortical tissue also tends to show water movement that is somewhat more constrained. In further multivariate analyses, in several left frontal, parietal, and temporal regions, the combination of shorter T1 relaxation time and faster MD was associated with older age and greater gynephilia in GD AFAB individuals and cisgender boys and with stronger attractions in cisgender boys only. Thus, for these cortical regions in these groups, older age, gynephilia, and stronger attractions (cisgender boys only) were associated with macromolecule-rich tissue in which water movement was freer-a pattern that some prior research suggests is associated with greater cell density and size. Overall, this study indicates that investigating T1 relaxation time and MD together can further inform how cortical gray matter tissue characteristics relate to age and psychosexuality.

Ktrans Calculation Using Reference Method Corrected Native T10 for Breast Cancer Diagnosis

Purpose

The objective of the study is to use multiple tube phantoms to generate correction factor at different spatial locations for each breast coil cuff to correct the native T10 value in the corresponding spatial location of the breast lesion. The corrected T10 value was used to compute Ktrans and analyze its diagnostic accuracy in the classification of target condition, i.e., breast tumors into malignant and benign.

Materials and Methods

Both in vitro phantom study (external reference) and patient's studies were acquired on simultaneous positron emission tomography/magnetic resonance imaging (PET/MRI) Biograph molecular magnetic resonance (mMR) system using 4 channel mMR breast coil. The spatial correction factors derived using multiple tube phantom were used for a retrospective analysis of dynamic contrast-enhanced (DCE) MRI data of 39 patients with a mean age of 50 years (31-77 years) having 51 enhancing breast lesions.

Results

Corrected and non-corrected receiver operating characteristic (ROC) curve analysis revealed a mean Ktrans value of 0.64 min-1 and 0.60 min-1, respectively. The sensitivity, specificity, positive predictive value (PPV), negative predictive value (NPV), and overall accuracy for non-corrected data were 86.21%, 81.82%, 86.20%, 81.81%, and 84.31%, respectively, and for corrected data were 93.10%, 86.36%, 90%, 90.47%, and 90.20% respectively. The area under curve (AUC) of corrected data was improved to 0.959 (95% confidence interval [CI] 0.862-0.994) from 0.824 (95% CI 0.694-0.918) of non-corrected data, and for NPV, it was improved to 90.47% from 81.81%, respectively.

Conclusion

T10 values were normalized using multiple tube phantom which was used for computation of Ktrans. We found significant improvement in the diagnostic accuracy of corrected Ktrans values that results in better characterization of breast lesions.

Rapid fat-water separated T1 mapping using a single-shot radial inversion-recovery spoiled gradient recalled pulse sequence

T1 mapping is increasingly used in clinical practice and research studies. With limited scan time, existing techniques often have limited spatial resolution, contrast resolution and slice coverage. High fat concentrations yield complex errors in Look-Locker T1 methods. In this study, a dual-echo 2D radial inversion-recovery T1 (DEradIR-T1) technique was developed for fast fat-water separated T1 mapping. The DEradIR-T1 technique was tested in phantoms, 5 volunteers and 28 patients using a 3 T clinical MRI scanner. In our study, simulations were performed to analyze the composite (fat + water) and water-only T1 under different echo times (TE). In standardized phantoms, an inversion-recovery spin echo (IR-SE) sequence with and without fat saturation pulses served as a T1 reference. Parameter mapping with DEradIR-T1 was also assessed in vivo, and values were compared with modified Look-Locker inversion recovery (MOLLI). Bland-Altman analysis and two-tailed paired t-tests were used to compare the parameter maps from DEradIR-T1 with the references. Simulations of the composite and water-only T1 under different TE values and levels of fat matched the in vivo studies. T1 maps from DEradIR-T1 on a NIST phantom (Pcomp  = 0.97) and a Calimetrix fat-water phantom (Pwater  = 0.56) matched with the references. In vivo T1 was compared with that of MOLLI: R comp 2 = 0.77 ; R water 2 = 0.72 . In this work, intravoxel fat is found to have a variable, echo-time-dependent effect on measured T1 values, and this effect may be mitigated using the proposed DRradIR-T1.

In vivo measurements of lamination patterns in the human cortex

The laminar composition of the cerebral cortex is tightly connected to the development and connectivity of the brain, as well as to function and pathology. Although most of the research on the cortical layers is done with the aid of ex vivo histology, there have been recent attempts to use magnetic resonance imaging (MRI) with potential in vivo applications. However, the high-resolution MRI technology and protocols required for such studies are neither common nor practical. In this article, we present a clinically feasible method for assessing the laminar properties of the human cortex using standard pulse sequence available on any common MRI scanner. Using a series of low-resolution inversion recovery (IR) MRI scans allows us to calculate multiple T1 relaxation time constants for each voxel. Based on the whole-brain T1 -distribution, we identify six different gray matter T1 populations and their variation across the cortex. Based on this, we show age-related differences in these population and demonstrate that this method is able to capture the difference in laminar composition across varying brain areas. We also provide comparison to ex vivo high-resolution MRI scans. We show that this method is feasible for the estimation of layer variability across large population cohorts, which can lead to research into the links between the cortical layers and function, behavior and pathologies that was heretofore unexplorable.

Synthetic MR: Physical principles, clinical implementation, and new developments

Current clinical MR imaging practices rely on the qualitative assessment of images for diagnosis and treatment planning. While contrast in MR images is dependent on the spin parameters of the imaged tissue, pixel values on MR images are relative and are not scaled to represent any tissue properties. Synthetic MR is a fully featured imaging workflow consisting of efficient multiparameter mapping acquisition, synthetic image generation, and volume quantitation of brain tissues. As the application becomes more widely available on multiple vendors and scanner platforms, it has also gained widespread adoption as clinicians begin to recognize the benefits of rapid quantitation. This review will provide details about the sequence with a focus on the physical principles behind its relaxometry mechanisms. It will present an overview of the products in their current form and some potential issues when implementing it in the clinic. It will conclude by highlighting some recent advances of the technique, including a 3D mapping method and its associated applications. This article is protected by copyright. All rights reserved.

Rapid MR relaxometry using deep learning: An overview of current techniques and emerging trends

Quantitative mapping of MR tissue parameters such as the spin-lattice relaxation time (T1 ), the spin-spin relaxation time (T2 ), and the spin-lattice relaxation in the rotating frame (T1ρ ), referred to as MR relaxometry in general, has demonstrated improved assessment in a wide range of clinical applications. Compared with conventional contrast-weighted (eg T1 -, T2 -, or T1ρ -weighted) MRI, MR relaxometry provides increased sensitivity to pathologies and delivers important information that can be more specific to tissue composition and microenvironment. The rise of deep learning in the past several years has been revolutionizing many aspects of MRI research, including image reconstruction, image analysis, and disease diagnosis and prognosis. Although deep learning has also shown great potential for MR relaxometry and quantitative MRI in general, this research direction has been much less explored to date. The goal of this paper is to discuss the applications of deep learning for rapid MR relaxometry and to review emerging deep-learning-based techniques that can be applied to improve MR relaxometry in terms of imaging speed, image quality, and quantification robustness. The paper is comprised of an introduction and four more sections. Section 2 describes a summary of the imaging models of quantitative MR relaxometry. In Section 3, we review existing "classical" methods for accelerating MR relaxometry, including state-of-the-art spatiotemporal acceleration techniques, model-based reconstruction methods, and efficient parameter generation approaches. Section 4 then presents how deep learning can be used to improve MR relaxometry and how it is linked to conventional techniques. The final section concludes the review by discussing the promise and existing challenges of deep learning for rapid MR relaxometry and potential solutions to address these challenges.

MULTI-parametric MR imaging with fLEXible design (MULTIPLEX)

PURPOSE

To introduce a gradient echo (GRE) -based method, namely MULTIPLEX, for single-scan 3D multi-parametric MRI with high resolution, signal-to-noise ratio (SNR), accuracy, efficiency, and acquisition flexibility.

THEORY

With a comprehensive design with dual-repetition time (TR), dual flip angle (FA), multi-echo, and optional flow modulation features, the MULTIPLEX signals contain information on radiofrequency (RF) B_1t fields, proton density, T₁ , susceptibility and blood flows, facilitating multiple qualitative images and parametric maps.

METHODS

MULTIPLEX was evaluated on system phantom and human brains, via visual inspection for image contrasts and quality or quantitative evaluation via simulation, phantom scans and literature comparison. Region-of-interest (ROI) analysis was performed on parametric maps of the system phantom and brain scans, extracting the mean and SD of the T₁ , T2∗ , proton density (PD), and/or quantitative susceptibility mapping (QSM) values for comparison with reference values or literature.

RESULTS

One MULTIPLEX scan offers multiple sets of images, including but not limited to: composited PDW/T₁ W/ T2∗ W, aT₁ W, SWI, MRA (optional), B_1t map, T₁ map, T2∗ / R2∗ maps, PD map, and QSM. The quantitative error of phantom T₁ , T2∗ and PD mapping were <5%, and those in brain scans were in good agreement with literature. MULTIPLEX scan times for high resolution (0.68 × 0.68 × 2 mm³ ) whole brain coverage were about 7.5 min, while processing times were <1 min. With flow modulation, additional MRA images can be obtained without affecting the quality or accuracy of other images.

CONCLUSION

The proposed MUTLIPLEX method possesses great potential for multi-parametric MR imaging.

Myelin water imaging in relapsing multiple sclerosis treated with ocrelizumab and interferon beta-1a

•

2-Year change in MS myelin water fraction favored ocrelizumab over interferon.

•

Matched healthy controls showed no change in myelin water fraction over 2 years.

•

Ocrelizumab appears to protect against demyelination in MS white matter and lesions.

Background

Myelin water imaging is a magnetic resonance imaging (MRI) technique that quantifies myelin damage and repair in multiple sclerosis (MS) via the myelin water fraction (MWF).

Objective

In this substudy of a phase 3 therapeutic trial, OPERA II, MWF was assessed in relapsing MS participants assigned to interferon beta-1a (IFNb-1a) or ocrelizumab (OCR) during a two-year double-blind period (DBP) followed by a two-year open label extension (OLE) with ocrelizumab treatment.

Methods

MWF in normal appearing white matter (NAWM), including both whole brain NAWM and 5 white matter structures, and chronic lesions, was assessed in 29 OCR and 26 IFNb-1a treated participants at weeks 0, 24, 48 and 96 (DBP), and weeks 144 and 192 (OLE), and in white matter for 23 healthy control participants at weeks 0, 48 and 96.

Results

Linear mixed-effects models of data from baseline to week 96 showed a difference in the change in MWF over time favouring ocrelizumab in all NAWM regions. At week 192, lesion MWF was lower for participants originally randomised to IFNb-1a compared to those originally randomised to OCR. Controls showed no change in MWF over 96 weeks in any region.

Conclusion

Ocrelizumab appears to protect against demyelination in MS NAWM and chronic lesions and may allow for a more permissive micro environment for remyelination to occur in focal and diffusely damaged tissue.

Using quantitative MRI to study brain responses to immune challenge with interferon-α

Inflammatory processes in the Central Nervous System (CNS) have been proposed to mediate the association between peripheral inflammation and the development of psychiatric disorders, but we currently lack sensitive measures of CNS inflammation for human studies in vivo. Here we used quantitative MRI (qMRI) to explore the in vivo central response to a peripheral immune challenge in healthy humans, and we assessed whether changes in quantitative relaxometry MRI parameters were associated with changes in peripheral inflammation. Quantitative relaxation times (T1 & T2) and Proton Density (PD) were measured in n ​= ​6 healthy males (mean age ​= ​30.5 ​± ​6.8 years) in two MRI assessments collected before and 24 ​hours after a subcutaneous injection of the proinflammatory cytokine and immune activator, interferon-alpha (IFN-α). Serum levels of immune markers and markers of blood-brain barrier integrity (S100B) were also measured before and after the injection. Region of interest and histogram-based analyses (optimized for the small sample size) showed a statistically significant increase of both T1 (t₍₅₎ ​= ​3.78, p ​= 0.013) and PD (t₍₅₎ ​= ​2.91, p ​= ​0.033) parameters in the bilateral hippocampus after IFN-α administration. T1 peak values in bilateral hippocampus were positively correlated with serum Tumour Necrosis Factor-alpha levels at 24 ​h after the injection, when this cytokine peaked (Spearman's rho ​= ​0.67, p ​= ​0.018) and negatively correlated with S100B levels (Spearman's rho ​= ​-0.826, p ​= ​0.001). Our data suggest that peripheral administration of IFN-α produces acute changes in brain qMRI which might indicate a brain immune response. This is supported by the association of such changes with low-grade peripheral inflammation.

A deep learning approach for synthetic MRI based on two routine sequences and training with synthetic data

BACKGROUND AND OBJECTIVE

Synthetic magnetic resonance imaging (MRI) is a low cost procedure that serves as a bridge between qualitative and quantitative MRI. However, the proposed methods require very specific sequences or private protocols which have scarcely found integration in clinical scanners. We propose a learning-based approach to compute T1, T2, and PD parametric maps from only a pair of T1- and T2-weighted images customarily acquired in the clinical routine.

METHODS

Our approach is based on a convolutional neural network (CNN) trained with synthetic data; specifically, a synthetic dataset with 120 volumes was constructed from the anatomical brain model of the BrainWeb tool and served as the training set. The CNN learns an end-to-end mapping function to transform the input T1- and T2-weighted images to their underlying T1, T2, and PD parametric maps. Then, conventional weighted images unseen by the network are analytically synthesized from the parametric maps. The network can be fine tuned with a small database of actual weighted images and maps for better performance.

RESULTS

This approach is able to accurately compute parametric maps from synthetic brain data achieving normalized squared error values predominantly below 1%. It also yields realistic parametric maps from actual MR brain acquisitions with T1, T2, and PD values in the range of the literature and with correlation values above 0.95 compared to the T1 and T2 maps obtained from relaxometry sequences. Further, the synthesized weighted images are visually realistic; the mean square error values are always below 9% and the structural similarity index is usually above 0.90. Network fine tuning with actual maps improves performance, while training exclusively with a small database of actual maps shows a performance degradation.

CONCLUSIONS

These results show that our approach is able to provide realistic parametric maps and weighted images out of a CNN that (a) is trained with a synthetic dataset and (b) needs only two inputs, which are in turn obtained from a common full-brain acquisition that takes less than 8 min of scan time. Although a fine tuning with actual maps improves performance, synthetic data is crucial to reach acceptable performance levels. Hence, we show the utility of our approach for both quantitative MRI in clinical viable times and for the synthesis of additional weighted images to those actually acquired.

Whole brain and deep gray matter structure segmentation: Quantitative comparison between MPRAGE and MP2RAGE sequences

Objective

T1-weighted MRI images are commonly used for volumetric assessment of brain structures. Magnetization prepared 2 rapid gradient echo (MP2RAGE) sequence offers superior gray (GM) and white matter (WM) contrast. This study aimed to quantitatively assess the agreement of whole brain tissue and deep GM (DGM) volumes obtained from MP2RAGE compared to the widely used MP-RAGE sequence.

Methods

Twenty-nine healthy participants were included in this study. All subjects underwent a 3T MRI scan acquiring high-resolution 3D MP-RAGE and MP2RAGE images. Twelve participants were re-scanned after one year. The whole brain, as well as DGM segmentation, was performed using CAT12, volBrain, and FSL-FAST automatic segmentation tools based on the acquired images. Finally, contrast-to-noise ratio between WM and GM (CNR_WG), the agreement between the obtained tissue volumes, as well as scan-rescan variability of both sequences were explored.

Results

Significantly higher CNR_WG was detected in MP2RAGE vs. MP-RAGE (Mean ± SD = 0.97 ± 0.04 vs. 0.8 ± 0.1 respectively; p<0.0001). Significantly higher total brain GM, and lower cerebrospinal fluid‏ volumes were obtained from MP2RAGE vs. MP-RAGE based on all segmentation methods (p<0.05 in all cases). Whole-brain voxel-wise comparisons revealed higher GM tissue probability in the thalamus, putamen, caudate, lingual gyrus, and precentral gyrus based on MP2RAGE compared with MP-RAGE. Moreover, significantly higher WM probability was observed in the cerebellum, corpus callosum, and frontal-and-temporal regions in MP2RAGE vs. MP-RAGE. Finally, MP2RAGE showed a higher mean percentage of change in total brain GM compared to MP-RAGE. On the other hand, MP-RAGE demonstrated a higher overtime percentage of change in WM and DGM volumes compared to MP2RAGE.

Conclusions

Due to its higher CNR, MP2RAGE resulted in reproducible brain tissue segmentation, and thus is a recommended method for volumetric imaging biomarkers for the monitoring of neurological diseases.

An asymmetrical whole-body birdcage RF coil without RF shield for hyperpolarized 129 Xe lung MR imaging at 1.5 T

PURPOSE

This study describes the development and testing of an asymmetrical xenon-129 (¹²⁹ Xe) birdcage radiofrequency (RF) coil for ¹²⁹ Xe lung ventilation imaging at 1.5 Tesla, which allows proton (¹ H) system body coil transmit-receive functionality.

METHODS

The ¹²⁹ Xe RF coil is a whole-body asymmetrical elliptical birdcage constructed without an outer RF shield to enable ¹ H imaging. B 1 + field homogeneity and flip angle mapping of the ¹²⁹ Xe birdcage RF coil and ¹ H system body RF coil with the ¹²⁹ Xe RF coil in situ were evaluated in the MR scanner. The functionality of the ¹²⁹ Xe birdcage RF coil was demonstrated through hyperpolarized ¹²⁹ Xe lung ventilation imaging with the birdcage in both transceiver configuration and transmit-only configuration when combined with an 8-channel ¹²⁹ Xe receive-only RF coil array. The functionality of ¹ H system body coil with the ¹²⁹ Xe RF coil in situ was demonstrated by acquiring coregistered ¹ H lung anatomical MR images.

RESULTS

The asymmetrical birdcage produced a homogeneous B 1 + field (±10%) in agreement with electromagnetic simulations. Simulations indicated an optimal detuning configuration with 4 diodes. The obtained g-factor of 1.4 for acceleration factor of R = 2 indicates optimal array configuration. Coregistered ¹ H anatomical images from the system body coil along with ¹²⁹ Xe lung images demonstrated concurrent and compatible arrangement of the RF coils.

CONCLUSION

A large asymmetrical birdcage for homogenous B 1 + transmission with high sensitivity reception for ¹²⁹ Xe lung MRI at 1.5 Tesla has been demonstrated. The unshielded asymmetrical birdcage design enables ¹ H structural lung MR imaging in the same exam.

Imaging neurovascular, endothelial and structural integrity in preparation to treat small vessel diseases. The INVESTIGATE-SVDs study protocol. Part of the SVDs@Target project

Background

Sporadic cerebral small vessel disease (SVD) and cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy (CADASIL) share clinical and neuroimaging features and possibly vascular dysfunction(s). However few studies have included both conditions, assessed more than one vascular dysfunction simultaneously, or included more than one centre. The INVESTIGATE-SVDs study will assess several cerebrovascular dysfunctions with MRI in participants with sporadic SVD or CADASIL at three European centres.

Methods

We will recruit participants with sporadic SVDs (ischaemic stroke or vascular cognitive impairment) and CADASIL in Edinburgh, Maastricht and Munich. We will perform detailed clinical and neuropsychological phenotyping of the participants, and neuroimaging including structural MRI, cerebrovascular reactivity MRI (CVR: using carbon dioxide challenge), phase contrast MRI (arterial, venous and CSF flow and pulsatility), dynamic contrast-enhanced MRI (blood brain barrier (BBB) leakage) and multishell diffusion imaging. Participants will measure their blood pressure (BP) and its variability over seven days using a telemetric device.

Discussion

INVESTIGATE-SVDs will assess the relationships of BBB integrity, CVR, pulsatility and CSF flow in sporadic SVD and CADASIL using a multisite, multimodal MRI protocol. We aim to establish associations between these measures of vascular function, risk factors particularly BP and its variability, and brain parenchymal lesions in these two SVD phenotypes. Additionally we will test feasibility of complex multisite MRI, provide reliable intermediary outcome measures and sample size estimates for future trials.

Diffusion-Weighted and Dynamic Contrast-Enhanced MRI Derived Imaging Metrics for Stereotactic Body Radiotherapy of Pancreatic Ductal Adenocarcinoma: Preliminary Findings

We aimed to assess longitudinal changes in quantitative imaging metric values obtained from diffusion-weighted (DW-) and dynamic contrast-enhanced magnetic resonance imaging (DCE)-MRI at pre-treatment (TX[0]), immediately after the first fraction of stereotactic body radiotherapy (D1-TX[1]), and 6 weeks post-TX (Post-TX[2]) in patients with pancreatic ductal adenocarcinoma. Ten enrolled patients (n = 10) underwent DW- and DCE-MRI examinations on a 3.0 T scanner. The apparent diffusion coefficient, ADC (mm²/s), was derived from DW imaging data using a monoexponential model. The tissue relaxation rate, R _1t, time-course data were fitted with a shutter-speed model, which provides estimates of the volume transfer constant, K ^trans (min^-1), extravascular extracellular volume fraction, v_e , and mean lifetime of intracellular water protons, τ _i (seconds). Wilcoxon rank-sum test compared the mean values, standard deviation, skewness, kurtosis, and relative percentage (r, %) changes (Δ) in ADC, K ^trans, v_e , and τ _i values between the magnetic resonance examinations. rADC_Δ2-0 values were significantly greater than rADC_Δ1-0 values (P = .009). rK ^trans _Δ2-0 values were significantly lower than rK ^trans _Δ1-0 values (P = .048). rv_{e Δ2-1} and rv_eΔ2-0 values were significantly different (P = .016). rτ _{i Δ2-1} values were significantly lower than rτ _{i Δ2-0} values (P = .008). For group comparison, the pre-TX mean and kurtosis of ADC (P = .18 and P = .14), skewness and kurtosis of K ^trans values (P = .14 for both) showed a leaning toward significant difference between patients who experienced local control (n = 2) and failed early (n = 4). DW- and DCE-MRI-derived quantitative metrics could be useful biomarkers to evaluate longitudinal changes to stereotactic body radiotherapy in patients with pancreatic ductal adenocarcinoma.

MESMERISED: Super-accelerating T1 relaxometry and diffusion MRI with STEAM at 7 T for quantitative multi-contrast and diffusion imaging

There is an increasing interest in quantitative imaging of T₁, T₂ and diffusion contrast in the brain due to greater robustness against bias fields and artifacts, as well as better biophysical interpretability in terms of microstructure. However, acquisition time constraints are a challenge, particularly when multiple quantitative contrasts are desired and when extensive sampling of diffusion directions, high b-values or long diffusion times are needed for multi-compartment microstructure modeling. Although ultra-high fields of 7 T and above have desirable properties for many MR modalities, the shortening T₂ and the high specific absorption rate (SAR) of inversion and refocusing pulses bring great challenges to quantitative T₁, T₂ and diffusion imaging. Here, we present the MESMERISED sequence (Multiplexed Echo Shifted Multiband Excited and Recalled Imaging of STEAM Encoded Diffusion). MESMERISED removes the dead time in Stimulated Echo Acquisition Mode (STEAM) imaging by an echo-shifting mechanism. The echo-shift (ES) factor is independent of multiband (MB) acceleration and allows for very high multiplicative (ESxMB) acceleration factors, particularly under moderate and long mixing times. This results in super-acceleration and high time efficiency at 7 T for quantitative T₁ and diffusion imaging, while also retaining the capacity to perform quantitative T₂ and B₁ mapping. We demonstrate the super-acceleration of MESMERISED for whole-brain T₁ relaxometry with total acceleration factors up to 36 at 1.8 mm isotropic resolution, and up to 54 at 1.25 mm resolution qT₁ imaging, corresponding to a 6x and 9x speedup, respectively, compared to MB-only accelerated acquisitions. We then demonstrate highly efficient diffusion MRI with high b-values and long diffusion times in two separate cases. First, we show that super-accelerated multi-shell diffusion acquisitions with 370 whole-brain diffusion volumes over 8 b-value shells up to b = 7000 s/mm² can be generated at 2 mm isotropic in under 8 minutes, a data rate of almost a volume per second, or at 1.8 mm isotropic in under 11 minutes, achieving up to 3.4x speedup compared to MB-only. A comparison of b = 7000 s/mm² MESMERISED against standard MB pulsed gradient spin echo (PGSE) diffusion imaging shows 70% higher SNR efficiency and greater effectiveness in supporting complex diffusion signal modeling. Second, we demonstrate time-efficient sampling of different diffusion times with 1.8 mm isotropic diffusion data acquired at four diffusion times up to 290 ms, which supports both Diffusion Tensor Imaging (DTI) and Diffusion Kurtosis Imaging (DKI) at each diffusion time. Finally, we demonstrate how adding quantitative T₂ and B₁⁺ mapping to super-accelerated qT₁ and diffusion imaging enables efficient quantitative multi-contrast mapping with the same MESMERISED sequence and the same readout train. MESMERISED extends possibilities to efficiently probe T₁, T₂ and diffusion contrast for multi-component modeling of tissue microstructure.

Magnetic resonance parameter mapping using model-guided self-supervised deep learning

PURPOSE

To develop a model-guided self-supervised deep learning MRI reconstruction framework called reference-free latent map extraction (RELAX) for rapid quantitative MR parameter mapping.

METHODS

Two physical models are incorporated for network training in RELAX, including the inherent MR imaging model and a quantitative model that is used to fit parameters in quantitative MRI. By enforcing these physical model constraints, RELAX eliminates the need for full sampled reference data sets that are required in standard supervised learning. Meanwhile, RELAX also enables direct reconstruction of corresponding MR parameter maps from undersampled k-space. Generic sparsity constraints used in conventional iterative reconstruction, such as the total variation constraint, can be additionally included in the RELAX framework to improve reconstruction quality. The performance of RELAX was tested for accelerated T1 and T2 mapping in both simulated and actually acquired MRI data sets and was compared with supervised learning and conventional constrained reconstruction for suppressing noise and/or undersampling-induced artifacts.

RESULTS

In the simulated data sets, RELAX generated good T1 /T2 maps in the presence of noise and/or undersampling artifacts, comparable to artifact/noise-free ground truth. The inclusion of a spatial total variation constraint helps improve image quality. For the in vivo T1 /T2 mapping data sets, RELAX achieved superior reconstruction quality compared with conventional iterative reconstruction, and similar reconstruction performance to supervised deep learning reconstruction.

CONCLUSION

This work has demonstrated the initial feasibility of rapid quantitative MR parameter mapping based on self-supervised deep learning. The RELAX framework may also be further extended to other quantitative MRI applications by incorporating corresponding quantitative imaging models.

Sparse precontrast T1 mapping for high-resolution whole-brain DCE-MRI

PURPOSE

To develop and evaluate an efficient precontrast T₁ mapping technique suitable for quantitative high-resolution whole-brain dynamic contrast-enhanced-magnetic resonance imaging (DCE-MRI).

METHODS

Variable flip angle (VFA) T₁ mapping was considered that provides 1 × 1 × 2 mm³ resolution to match a recent high-resolution whole-brain DCE-MRI protocol. Seven FAs were logarithmically spaced from 1.5° to 15°. T₁ and M₀ maps were estimated using model-based reconstruction. This approach was evaluated using an anatomically realistic brain tumor digital reference object (DRO) with noise-mimicking 3T neuroimaging and fully sampled data acquired from one healthy volunteer. Methods were also applied on fourfold prospectively undersampled VFA data from 13 patients with high-grade gliomas.

RESULTS

T₁ -mapping precision decreased with undersampling factor R, althoughwhereas bias remained small before a critical R. In the noiseless DRO, T₁ bias was <25 ms in white matter (WM) and <11 ms in brain tumor (BT). T₁ standard deviation (SD) was <119.5 ms in WM (coefficient of variation [COV] ~11.0%) and <253.2 ms in BT (COV ~12.7%). In the noisy DRO, T₁ bias was <50 ms in WM and <30 ms in BT. For R ≤ 10, T₁ SD was <107.1 ms in WM (COV ~9.9%) and <240.9 ms in BT (COV ~12.1%). In the healthy subject, T₁ bias was <30 ms for R ≤ 16. At R = 4, T₁ SD was 171.4 ms (COV ~13.0%). In the prospective brain tumor study, T₁ values were consistent with literature values in WM and BT.

CONCLUSION

High-resolution whole-brain VFA T₁ mapping is feasible with sparse sampling, supporting its use for quantitative DCE-MRI.

Sources of systematic error in DCE-MRI estimation of low-level blood-brain barrier leakage

PURPOSE

Dynamic contrast-enhanced (DCE) -MRI with Patlak model analysis is increasingly used to quantify low-level blood-brain barrier (BBB) leakage in studies of pathophysiology. We aimed to investigate systematic errors due to physiological, experimental, and modeling factors influencing quantification of the permeability-surface area product PS and blood plasma volume v_p , and to propose modifications to reduce the errors so that subtle differences in BBB permeability can be accurately measured.

METHODS

Simulations were performed to predict the effects of potential sources of systematic error on conventional PS and v_p quantification: restricted BBB water exchange, reduced cerebral blood flow, arterial input function (AIF) delay and B 1 + error. The impact of targeted modifications to the acquisition and processing were evaluated, including: assumption of fast versus no BBB water exchange, bolus versus slow injection of contrast agent, exclusion of early data from model fitting and B 1 + correction. The optimal protocol was applied in a cohort of recent mild ischaemic stroke patients.

RESULTS

Simulation results demonstrated substantial systematic errors due to the factors investigated (absolute PS error ≤ 4.48 × 10^-4 min^-1 ). However, these were reduced (≤0.56 × 10^-4 min^-1 ) by applying modifications to the acquisition and processing pipeline. Processing modifications also had substantial effects on in-vivo normal-appearing white matter PS estimation (absolute change ≤ 0.45 × 10^-4 min^-1 ).

CONCLUSION

Measuring subtle BBB leakage with DCE-MRI presents unique challenges and is affected by several confounds that should be considered when acquiring or interpreting such data. The evaluated modifications should improve accuracy in studies of neurodegenerative diseases involving subtle BBB breakdown.

A model-based framework for correcting B 1 + inhomogeneity effects in magnetization transfer saturation and inhomogeneous magnetization transfer saturation maps

PURPOSE

In this work, we propose that Δ B 1 + -induced errors in magnetization transfer (MT) saturation (MT_sat ) maps can be corrected with use of an R₁ and B 1 + map and through numerical simulations of the sequence.

THEORY AND METHODS

One healthy subject was scanned at 3.0T using a partial quantitative MT protocol to estimate the relationship between observed R₁ (R_1,obs ) and apparent bound pool size ( M 0 , a p p B ) in the brain. MT_sat values were simulated for a range of B 1 + , R_1,obs , and M 0 , a p p B . An equation was fit to the simulated MT_sat , then a linear relationship between R_1,obs and M 0 , a p p B was generated. These results were used to generate correction factor maps for the MT_sat acquired from single-point data. The proposed correction was compared to an empirical correction factor with different MT-preparation schemes.

RESULTS

M 0 , a p p B was highly correlated with R_1,obs (r > 0.96), permitting the use of R_1,obs to estimate M 0 , a p p B for B 1 + correction. All B 1 + corrected MT_sat maps displayed a decreased correlation with B 1 + compared to uncorrected MT_sat and MT_sat corrected with an empirical factor in the corpus callosum. There was good agreement between the proposed approach and the empirical correction with radiofrequency saturation at 2 kHz, with larger deviations seen when using saturation pulses further off-resonance and in inhomogeneous (ih) MT_sat maps.

CONCLUSION

The proposed correction decreases the dependence of MT_sat on B 1 + inhomogeneities. Furthermore, this flexible framework permits the use of different saturation protocols, making it useful for correcting B 1 + inhomogeneities in ihMT.

A Multi-Modal MRI Analysis of Cortical Structure in Relation to Gender Dysphoria, Sexual Orientation, and Age in Adolescents

Gender dysphoria (GD) is characterized by distress due to an incongruence between experienced gender and sex assigned at birth. Sex-differentiated brain regions are hypothesized to reflect the experienced gender in GD and may play a role in sexual orientation development. Magnetic resonance brain images were acquired from 16 GD adolescents assigned female at birth (AFAB) not receiving hormone therapy, 17 cisgender girls, and 14 cisgender boys (ages 12–17 years) to examine three morphological and microstructural gray matter features in 76 brain regions: surface area (SA), cortical thickness (CT), and T1 relaxation time. Sexual orientation was represented by degree of androphilia-gynephilia and sexual attraction strength. Multivariate analyses found that cisgender boys had larger SA than cisgender girls and GD AFAB. Shorter T1, reflecting denser, macromolecule-rich tissue, correlated with older age and stronger gynephilia in cisgender boys and GD AFAB, and with stronger attractions in cisgender boys. Thus, cortical morphometry (mainly SA) was related to sex assigned at birth, but not experienced gender. Effects of experienced gender were found as similarities in correlation patterns in GD AFAB and cisgender boys in age and sexual orientation (mainly T1), indicating the need to consider developmental trajectories and sexual orientation in brain studies of GD.

Comparison of multi echo T2 relaxation and steady state approaches for myelin imaging in the central nervous system

The traditional approach for measuring myelin-associated water with quantitative magnetic resonance imaging (MRI) uses multi-echo T₂ relaxation data to calculate the myelin water fraction (MWF). A fundamentally different approach, abbreviated “mcDESPOT”, uses a more efficient steady-state acquisition to generate an equivalent metric (f_M). Although previous studies have demonstrated inherent instability and bias in the complex mcDESPOT analysis procedure, f_M has often been used as a surrogate for MWF. We produced and compared multivariate atlases of MWF and f_M in healthy human brain and cervical spinal cord (available online) and compared their ability to detect multiple sclerosis pathology. A significant bias was found in all regions (p < 10^–5), albeit reversed for spinal cord (f_M-MWF =  − 3.4%) compared to brain (+ 6.2%). MWF and f_M followed an approximately linear relationship for regions with MWF <  ~ 10%. For MWF >  ~ 10%, the relationship broke down and f_M no longer increased in tandem with MWF. For multiple sclerosis patients, MWF and f_M Z score maps showed overlapping areas of low Z score and similar trends between patients and brain regions, although those of f_M generally had greater spatial extent and magnitude of severity. These results will guide future choice of myelin-sensitive quantitative MRI and improve interpretation of studies using either myelin imaging approach.

Bi-component T2 mapping correlates with articular cartilage material properties

Non-invasive estimation of cartilage material properties is useful for understanding cartilage health and creating subject-specific computational models. Bi-component T2 mapping measured using Multi-Component Driven Equilibrium Single Shot Observation of T₁ and T₂ (mcDESPOT) is sensitive for detecting cartilage degeneration within the human knee joint, but has not been correlated with cartilage composition and mechanical properties. Therefore, the purpose of this study was to investigate the relationship between bi-component T₂ parameters measured using mcDESPOT at 3.0 T and cartilage composition and mechanical properties. Ex-vivo patellar cartilage specimens harvested from five human cadaveric knees were imaged using mcDESPOT at 3.0 T. Cartilage samples were removed from the patellae, mechanically tested to determine linear modulus and dissipated energy, and chemically tested to determine proteoglycan and collagen content. Parameter maps of single-component T₂ relaxation time (T₂), the T₂ relaxation times of the fast relaxing macromolecular bound water component (T_2F) and slow relaxing bulk water component (T_2S), and the fraction of the fast relaxing macromolecular bound water component (F_F) were compared to mechanical and chemical measures using linear regression. F_F was significantly (p < 0.05) correlated with energy dissipation and linear modulus. T₂ was significantly (p ≤ 0.05) correlated with elastic modulus at 1 Hz and energy dissipated at all frequencies. There were no other significant (p = 0.13-0.97) correlations between mcDESPOT parameters and mechanical properties. F_F was significantly (p = 0.04) correlated with proteoglycan content. There were no other significant (p = 0.19-0.92) correlations between mcDESPOT parameters and proteoglycan or collagen content. This study suggests that F_F measured using mcDESPOT at 3.0 T could be used to non-invasively estimate cartilage proteoglycan content, elastic modulus, and energy dissipation.

A practical method for post-acquisition reduction of bias in fast, whole-brain B1-maps

Large consistent differences have been observed between maps of the flip angle correction factor (commonly called "B₁-maps") produced with different fast methods in the human brain. We present an empirical procedure for first-order multiplicative bias correction that can be applied when more than one B₁-mapping method is available. We use a B₁-map measurement in a calibration phantom as a reference and the voxel-wise histogram mode between ratios of B₁-maps produced from different methods to calculate determine the bias as a multiplicative correcting scale factor. Institutional implementations of four common methods of B₁-mapping were assessed: Method of Slopes, FSE and EPI double angle methods (DAM), and Bloch-Siegert. In human subjects, the multiplicative bias used to correct for each of the four methods was: Method of Slopes = 1.005, FSE-DAM = 0.956, EPI-DAM = 1.080, and Bloch-Siegert = 1.128. Scaling to remove this bias between methods produces more consistent B₁-maps which enable more consistent values for any computations requiring flip angle correction. In addition, we present evidence that the corrected B₁ maps, using our calibration method, are also more accurate.

Gray and white matter are both affected in classical galactosemia: An explorative study on the association between neuroimaging and clinical outcome

BACKGROUND

Classical Galactosemia (CG) is an inherited disorder of galactose metabolism caused by a deficiency of the galactose-1-phosphate uridylyltransferase (GALT) enzyme resulting in neurocognitive complications. As in many Inborn Errors of Metabolism, the metabolic pathway of CG is well-defined, but the pathophysiology and high variability in clinical outcome are poorly understood. The aim of this study was to investigate structural changes of the brain of CG patients on MRI and their association with clinical outcome.

METHODS

In this prospective cohort study an MRI protocol was developed to evaluate gray matter (GM) and white matter (WM) volume of the cerebrum and cerebellum, WM hyperintensity volume, WM microstructure and myelin content with the use of conventional MRI techniques, diffusion tensor imaging (DTI) and quantitative T1 mapping. The association between several neuroimaging parameters and both neurological and intellectual outcome was investigated.

RESULTS

Twenty-one patients with CG (median age 22 years, range 8-47) and 24 controls (median age 30, range 16-52) were included. Compared to controls, the WM of CG patients was lower in volume and the microstructure of WM was impaired both in the whole brain and corticospinal tract (CST) and the lower R1 values of WM, GM and the CST were indicative of less myelin. The volume of WM lesions were comparable between patients and controls. The 9/16 patients with a poor neurological outcome (defined as the presence of a tremor and/or dystonia), demonstrated a lower WM volume, an impaired WM microstructure and lower R1 values of the WM indicative of less myelin content compared to 7/16 patients without movement disorders. In 15/21 patients with a poor intellectual outcome (defined as an IQ < 85) both GM and WM were affected with a lower cerebral and cerebellar WM and GM volume compared to 6/21 patients with an IQ ≥ 85. Both the severity of the tremor (as indicated by the Tremor Rating Scale) and IQ (as continuous measure) were associated with several neuroimaging parameters such as GM volume, WM volume, CSF volume, WM microstructure parameters and R1 values of GM and WM.

CONCLUSION

In this explorative study performed in patients with Classical Galactosemia, not only WM but also GM pathology was found, with more severe brain abnormalities on MRI in patients with a poor neurological and intellectual outcome. The finding that structural changes of the brain were associated with the severity of long-term complications indicates that quantitative MRI techniques could be of use to explain neurological and cognitive dysfunction as part of the disease spectrum. Based on the clinical outcome of patients, the absence of widespread WM lesions and the finding that both GM and WM are affected, CG could be primarily a GM disease with secondary damage to the WM as a result of neuronal degeneration. To investigate this further the course of GM and WM should be evaluated in longitudinal research, which could also clarify if CG is a neurodegenerative disease.

Methodological Considerations for Neuroimaging in Deep Brain Stimulation of the Subthalamic Nucleus in Parkinson's Disease Patients

Deep brain stimulation (DBS) of the subthalamic nucleus is a neurosurgical intervention for Parkinson's disease patients who no longer appropriately respond to drug treatments. A small fraction of patients will fail to respond to DBS, develop psychiatric and cognitive side-effects, or incur surgery-related complications such as infections and hemorrhagic events. In these cases, DBS may require recalibration, reimplantation, or removal. These negative responses to treatment can partly be attributed to suboptimal pre-operative planning procedures via direct targeting through low-field and low-resolution magnetic resonance imaging (MRI). One solution for increasing the success and efficacy of DBS is to optimize preoperative planning procedures via sophisticated neuroimaging techniques such as high-resolution MRI and higher field strengths to improve visualization of DBS targets and vasculature. We discuss targeting approaches, MRI acquisition, parameters, and post-acquisition analyses. Additionally, we highlight a number of approaches including the use of ultra-high field (UHF) MRI to overcome limitations of standard settings. There is a trade-off between spatial resolution, motion artifacts, and acquisition time, which could potentially be dissolved through the use of UHF-MRI. Image registration, correction, and post-processing techniques may require combined expertise of traditional radiologists, clinicians, and fundamental researchers. The optimization of pre-operative planning with MRI can therefore be best achieved through direct collaboration between researchers and clinicians.

Whole-brain R1 predicts manganese exposure and biological effects in welders

Manganese (Mn) is a neurotoxicant that, due to its paramagnetic property, also functions as a magnetic resonance imaging (MRI) T1 contrast agent. Previous studies in Mn toxicity have shown that Mn accumulates in the brain, which may lead to parkinsonian symptoms. In this article, we trained support vector machines (SVM) using whole-brain R1 (R1 = 1/T1) maps from 57 welders and 32 controls to classify subjects based on their air Mn concentration ([Mn]_Air), Mn brain accumulation (ExMn_Brain), gross motor dysfunction (UPDRS), thalamic GABA concentration (GABA_Thal), and total years welding. R1 was highly predictive of [Mn]_Air above a threshold of 0.20 mg/m³ with an accuracy of 88.8% and recall of 88.9%. R1 was also predictive of subjects with GABA_Thal having less than or equal to 2.6 mM with an accuracy of 82% and recall of 78.9%. Finally, we used an SVM to predict age as a method of verifying that the results could be attributed to Mn exposure. We found that R1 was predictive of age below 48 years of age with accuracies ranging between 75 and 82% with recall between 94.7% and 76.9% but was not predictive above 48 years of age. Together, this suggests that lower levels of exposure (< 0.20 mg/m³ and < 18 years of welding on the job) do not produce discernable signatures, whereas higher air exposures and subjects with more total years welding produce signatures in the brain that are readily identifiable using SVM.

An 8-element Tx/Rx array utilizing MEMS detuning combined with 6 Rx loops for 19 F and 1 H lung imaging at 1.5T

PURPOSE

To firstly improve the attainable image SNR of ¹⁹ F and ¹ H C₃ F₈ lung imaging at 1.5 tesla using an 8-element transmit/receive (Tx/Rx) flexible vest array combined with a 6-element Rx-only array, and to secondly evaluate microelectromechanical systems for switching the array elements between the 2 resonant frequencies.

METHODS

The Tx efficiency and homogeneity of the 8-element array were measured and simulated for ¹ H imaging in a cylindrical phantom and then evaluated for in vivo ¹⁹ F/¹ H imaging. The added improvement provided by the 6-element Rx-only array was quantified through simulation and measurement and compared to the ultimate SNR. It was verified through the measurement of isolation that microelectromechanical systems switches provided broadband isolation of Tx/Rx circuitry such that the ¹⁹ F tuned Tx/Rx array could be effectively used for both ¹⁹ F and ¹ H nuclei.

RESULTS

For ¹ H imaging, the measured Tx efficiency/homogeneity (mean ± percent SD; 6.79 μ T / kW ± 26 % ) was comparable to that simulated ( 7.57 μ T / kW ± 20 % ). The 6 additional Rx-only loops increased the mean Rx sensitivity when compared to the 8-element array by a factor of 1.41× and 1.45× in simulation and measurement, respectively. In regions central to the thorax, the simulated SNR of the 14-element array achieves ≥70% of the ultimate SNR when including noise from the matching circuits and preamplifiers. A measured microelectromechanical systems switching speed of 12 µs and added minimum 22 dB of isolation between Tx and Rx were sufficient for Tx/Rx switching in this application.

CONCLUSION

The described single-tuned array driven at ¹⁹ F and ¹ H, utilizing microelectromechanical systems technology, provides excellent results for ¹⁹ F and ¹ H dual-nuclear lung ventilation imaging.

Microstructural imaging in temporal lobe epilepsy: Diffusion imaging changes relate to reduced neurite density

PURPOSE

Previous imaging studies in patients with refractory temporal lobe epilepsy (TLE) have examined the spatial distribution of changes in imaging parameters such as diffusion tensor imaging (DTI) metrics and cortical thickness. Multi-compartment models offer greater specificity with parameters more directly related to known changes in TLE such as altered neuronal density and myelination. We studied the spatial distribution of conventional and novel metrics including neurite density derived from NODDI (Neurite Orientation Dispersion and Density Imaging) and myelin water fraction (MWF) derived from mcDESPOT (Multi-Compartment Driven Equilibrium Single Pulse Observation of T1/T2)] to infer the underlying neurobiology of changes in conventional metrics.

METHODS

20 patients with TLE and 20 matched controls underwent magnetic resonance imaging including a volumetric T1-weighted sequence, multi-shell diffusion from which DTI and NODDI metrics were derived and a protocol suitable for mcDESPOT fitting. Models of the grey matter-white matter and grey matter-CSF surfaces were automatically generated from the T1-weighted MRI. Conventional diffusion and novel metrics of neurite density and MWF were sampled from intracortical grey matter and subcortical white matter surfaces and cortical thickness was measured.

RESULTS

In intracortical grey matter, diffusivity was increased in the ipsilateral temporal and frontopolar cortices with more restricted areas of reduced neurite density. Diffusivity increases were largely related to reductions in neurite density, and to a lesser extent CSF partial volume effects, but not MWF. In subcortical white matter, widespread bilateral reductions in fractional anisotropy and increases in radial diffusivity were seen. These were primarily related to reduced neurite density, with an additional relationship to reduced MWF in the temporal pole and anterolateral temporal neocortex. Changes were greater with increasing epilepsy duration. Bilaterally reduced cortical thickness in the mesial temporal lobe and centroparietal cortices was unrelated to neurite density and MWF.

CONCLUSIONS

Diffusivity changes in grey and white matter are primarily related to reduced neurite density with an additional relationship to reduced MWF in the temporal pole. Neurite density may represent a more sensitive and specific biomarker of progressive neuronal damage in refractory TLE that deserves further study.

Prenatal exposure to maternal cigarette smoking and structural properties of the human corpus callosum

Alterations induced by prenatal exposure to nicotine have been observed in experimental (rodent) studies. While numerous developmental outcomes have been associated with prenatal exposure to maternal cigarette smoking (PEMCS) in humans, the possible relation with brain structure is less clear. Here we sought to elucidate the relation between PEMCS and structural properties of human corpus callosum in adolescence and early adulthood in a total of 1,747 youth. We deployed three community-based cohorts of 446 (age 25-27 years, 46% exposed), 934 (age 12-18 years, 47% exposed) and 367 individuals (age 18-21 years, 9% exposed). A mega-analysis revealed lower mean diffusivity in the callosal segments of exposed males. We speculate that prenatal exposure to maternal cigarette smoking disrupts the early programming of callosal structure and increases the relative portion of small-diameter fibres.

Estimating axon conduction velocity in vivo from microstructural MRI

The conduction velocity (CV) of action potentials along axons is a key neurophysiological property central to neural communication. The ability to estimate CV in humans in vivo from non-invasive MRI methods would therefore represent a significant advance in neuroscience. However, there are two major challenges that this paper aims to address: (1) Much of the complexity of the neurophysiology of action potentials cannot be captured with currently available MRI techniques. Therefore, we seek to establish the variability in CV that can be captured when predicting CV purely from parameters that have been reported to be estimatable from MRI: inner axon diameter (AD) and g-ratio. (2) errors inherent in existing MRI-based biophysical models of tissue will propagate through to estimates of CV, the extent to which is currently unknown. Issue (1) is investigated by performing a sensitivity analysis on a comprehensive model of axon electrophysiology and determining the relative sensitivity to various morphological and electrical parameters. The investigations suggest that 85% of the variance in CV is accounted for by variation in AD and g-ratio. The observed dependency of CV on AD and g-ratio is well characterised by the previously reported model by Rushton. Issue (2) is investigated through simulation of diffusion and relaxometry MRI data for a range of axon morphologies, applying models of restricted diffusion and relaxation processes to derive estimates of axon volume fraction (AVF), AD and g-ratio and estimating CV from the derived parameters. The results show that errors in the AVF have the biggest detrimental impact on estimates of CV, particularly for sparse fibre populations (AVF<0.3). For our equipment set-up and acquisition protocol, CV estimates are most accurate (below 5% error) where AVF is above 0.3, g-ratio is between 0.6 and 0.85 and AD is high (above 4μm). CV estimates are robust to errors in g-ratio estimation but are highly sensitive to errors in AD estimation, particularly where ADs are small. We additionally show CV estimates in human corpus callosum in a small number of subjects. In conclusion, we demonstrate accurate CV estimates are possible in regions of the brain where AD is sufficiently large. Problems with estimating ADs for smaller axons presents a problem for estimating CV across the whole CNS and should be the focus of further study.

Feasibility of Quantitative Magnetic Resonance Fingerprinting in Ovarian Tumors for T1 and T2 Mapping in a PET/MR Setting

Multiparametric magnetic resonance imaging (MRI) can be used to characterize many cancer subtypes including ovarian cancer. Quantitative mapping of MRI relaxation values, such as T 1 and T 2 mapping, is promising for improving tumor assessment beyond conventional qualitative T 1- and T 2-weighted images. However, quantitative MRI relaxation mapping methods often involve long scan times due to sequentially measuring many parameters. Magnetic resonance fingerprinting (MRF) is a new method that enables fast quantitative MRI by exploiting the transient signals caused by the variation of pseudorandom sequence parameters. These transient signals are then matched to a simulated dictionary of T 1 and T 2 values to create quantitative maps. The ability of MRF to simultaneously measure multiple parameters, could represent a new approach to characterizing cancer and assessing treatment response. This feasibility study investigates MRF for simultaneous T 1, T 2, and relative proton density (rPD) mapping using ovarian cancer as a model system.