Influence of RF spoiling on the stability and accuracy of T1 mapping based on spoiled FLASH with varying flip angles

Improved gradient echo magnitude- and phase-based mapping of T 2 $$ {\mathrm{T}}_2 $$ using multiple RF spoiling increments at 3T and 7T

PURPOSE

The transverse relaxation time T  2 $$ {}_2 $$ holds significant relevance in clinical applications and research studies. Conventional T  2 $$ {}_2 $$ mapping approaches rely on spin-echo sequences, which require lengthy acquisition times and involve high radiofrequency (RF) power deposition. An alternative gradient echo (GRE) phase-based T  2 $$ {}_2 $$ mapping method, utilizing steady-state acquisitions at one small RF spoil phase increment, was recently demonstrated. Here, a modified magnitude- and phase-based T  2 $$ {}_2 $$ mapping approach is proposed, which improvesT 2 $$ {\mathrm{T}}_2 $$ estimations by simultaneous fitting ofT 1 $$ {\mathrm{T}}_1 $$ and signal amplitude (A ∝ P D $$ A\propto PD $$ ) at three or more RF spoiling phase increments, instead of assuming a fixedT 1 $$ {\mathrm{T}}_1 $$ value.

METHODS

The feasibility of the magnitude-phase-based method was assessed by simulations, in phantom and in vivo experiments using skipped-CAIPI three-dimensional-echo-planar imaging (3D-EPI) for rapid GRE imaging.T 2 $$ {\mathrm{T}}_2 $$ ,T 1 $$ {\mathrm{T}}_1 $$ and PD estimations obtained by our method were compared toT 2 $$ {\mathrm{T}}_2 $$ of the phase-based method andT 1 $$ {\mathrm{T}}_1 $$ and PD of spoiled GRE-based multi-parameter mapping using a multi-echo version of the same sequence.

RESULTS

The agreement of the proposedT 2 $$ {\mathrm{T}}_2 $$ with ground truth and referenceT 2 $$ {\mathrm{T}}_2 $$ values was higher than that of phase-basedT 2 $$ {\mathrm{T}}_2 $$ in simulations and in phantom data. While phase-basedT 2 $$ {\mathrm{T}}_2 $$ overestimation increases with actualT 2 $$ {\mathrm{T}}_2 $$ andT 1 $$ {\mathrm{T}}_1 $$ , the proposed method is accurate over a large range of physiologically meaningfulT 2 $$ {\mathrm{T}}_2 $$ andT 1 $$ {\mathrm{T}}_1 $$ values. At the same time, precision is improved. In vivo results were in line with these observations.

CONCLUSION

Accurate magnitude-phase-based T  2 $$ {}_2 $$ mapping is feasible in less than 5 min scan time for 1 mm nominal isotropic whole-head coverage at 3T and 7T.

Characterization of the orientation dependence of magnetization transfer measures in single and crossing-fiber white matter

PURPOSE

To fully characterize the orientation dependence of magnetization transfer (MT) and inhomogeneous MT (ihMT) measures in the whole white matter (WM), for both single-fiber and crossing-fiber voxels.

METHODS

A characterization method was developed using the fiber orientation obtained from diffusion MRI (dMRI) with diffusion tensor imaging (DTI) and constrained spherical deconvolution. This allowed for characterization of the orientation dependence of measures in all of WM, regardless of the number of fiber orientation in a voxel. Furthermore, the orientation dependence inside 31 different WM bundles was characterized to evaluate the homogeneity of the effect. Variation of the results within and between-subject was assessed from a 12-subject dataset.

RESULTS

Previous results for single-fiber voxels were reproduced and a novel characterization was produced in voxels of crossing fibers, which seems to follow trends consistent with single-fiber results. Heterogeneity of the orientation dependence across bundles was observed, but homogeneity within similar bundles was also highlighted. Differences in behavior between MT and ihMT measures, as well as the ratio and saturation versions of these, were noted.

CONCLUSION

Orientation dependence characterization was proven possible over the entirety of WM. The vast range of effects and subtleties of the orientation dependence on MT measures showed the need for, but also the challenges of, a correction method.

Fast magnetization transfer saturation imaging of the brain using MP2RAGE T1 mapping

PURPOSE

Magnetization transfer saturation (MTsat) mapping is commonly used to examine the macromolecular content of brain tissue. This study compared variable flip angle (VFA) T1 mapping against compressed-sensing MP2RAGE (csMP2RAGE) T1 mapping for accelerating MTsat imaging.

METHODS

VFA, MP2RAGE, and csMP2RAGE were compared against inversion-recovery T1 in an aqueous phantom at 3 T. The same 1-mm VFA, MP2RAGE, and csMP2RAGE protocols were acquired in 4 healthy subjects to compare T1 and MTsat. Bloch-McConnell simulations were used to investigate differences between the phantom and in vivo T1 results. Ten healthy controls were imaged twice with the csMP2RAGE MTsat protocol to quantify repeatability.

RESULTS

The MP2RAGE and csMP2RAGE protocols were 13.7% and 32.4% faster than the VFA protocol, respectively. At these scan times, all approaches provided strong repeatability and accurate T1 times (< 5% difference) in the phantom, but T1 accuracy was more impacted by T2 for VFA than for MP2RAGE. In vivo, VFA estimated longer T1 times than MP2RAGE and csMP2RAGE. Simulations suggest that the differences in the T1 measured using VFA, MP2RAGE, and inversion recovery could be explained by the magnetization-transfer effects. In the test-retest experiment, we found that the csMP2RAGE has a minimum detectable change of 2.3% for T1 mapping and 7.8% for MTsat imaging.

CONCLUSIONS

We demonstrated that MP2RAGE can be used in place of VFA T1 mapping in an MTsat protocol. Furthermore, a shorter scan time and high repeatability can be achieved using the csMP2RAGE sequence.

Harnessing myelin water fraction as an imaging biomarker of human cerebral aging, neurodegenerative diseases, and risk factors influencing myelination: A review

Myelin water fraction (MWF) imaging has emerged as a promising magnetic resonance imaging (MRI) biomarker for investigating brain function and composition. This comprehensive review synthesizes the current state of knowledge on MWF as a biomarker of human cerebral aging, neurodegenerative diseases, and risk factors influencing myelination. The databases used include Web of Science, Scopus, Science Direct, and PubMed. We begin with a brief discussion of the theoretical foundations of MWF imaging, including its basis in MR physics and the mathematical modeling underlying its calculation, with an overview of the most adopted MRI methods of MWF imaging. Next, we delve into the clinical and research applications that have been explored to date, highlighting its advantages and limitations. Finally, we explore the potential of MWF to serve as a predictive biomarker for neurological disorders and identify future research directions for optimizing MWF imaging protocols and interpreting MWF in various contexts. By harnessing the power of MWF imaging, we may gain new insights into brain health and disease across the human lifespan, ultimately informing novel diagnostic and therapeutic strategies.

Multicompartment imaging of the brain using a comprehensive MR imaging protocol

In this study, we describe a comprehensive 3D magnetic resonance imaging (MRI) protocol designed to assess major tissue and fluid components in the brain. The protocol comprises four different sequences: 1) magnetization transfer prepared Cones (MT-Cones) for two-pool MT modeling to quantify macromolecular content; 2) short-TR adiabatic inversion-recovery prepared Cones (STAIR-Cones) for myelin water imaging; 3) proton-density weighted Cones (PDw-Cones) for total water imaging; and 4) highly T2 weighted Cones (T2w-Cones) for free water imaging. By integrating these techniques, we successfully mapped key brain components-namely macromolecules, myelin water, intra/extracellular water, and free water-in ten healthy volunteers and five patients with multiple sclerosis (MS) using a 3T clinical scanner. Brain macromolecular proton fraction (MMPF), myelin water proton fraction (MWPF), intra/extracellular water proton fraction (IEWPF), and free water proton fraction (FWPF) values were generated in white matter (WM), grey matter (GM), and MS lesions. Excellent repeatability of the protocol was demonstrated with high intra-class correlation coefficient (ICC) values. In MS patients, the MMPF and MWPF values of the lesions and normal-appearing WM (NAWM) were significantly lower than those in normal WM (NWM) in healthy volunteers. Moreover, we observed significantly higher FWPF values in MS lesions compared to those in NWM and NAWM regions. This study demonstrates the capability of our technique to volumetrically map major brain components. The technique may have particular value in providing a comprehensive assessment of neuroinflammatory and neurodegenerative diseases of the brain.

Whole brain multiparametric mapping in two minutes using a dual-flip-angle stack-of-stars blipped multi-gradient-echo acquisition

A new MRI technique is presented for three-dimensional fast simultaneous whole brain mapping of myelin water fraction (MWF), T1, proton density (PD), R2*, magnetic susceptibility (QSM), and B1 transmit field (B1+). Phantom and human (N = 9) datasets were acquired using a dual-flip-angle blipped multi-gradient-echo (DFA-mGRE) sequence with a stack-of-stars (SOS) trajectory. Images were reconstructed using a subspace-based algorithm with a locally low-rank constraint. A novel joint-sparsity-constrained multicomponent T2*-T1 spectrum estimation (JMSE) algorithm is proposed to correct for the T1 saturation effect and B1+/B1- inhomogeneities in the quantification of MWF. A tissue-prior-based B1+ estimation algorithm was adapted for B1 correction in the mapping of T1 and PD. In the phantom study, measurements obtained at an acceleration factor (R) of 12 using prospectively under-sampled SOS showed good consistency (R2 > 0.997) with Cartesian reference for R2*/T1app/M0app. In the in vivo study, results of retrospectively under-sampled SOS with R = 6, 12, 18, showed good quality (structure similarity index measure > 0.95) compared with those of fully-sampled SOS. Besides, results of prospectively under-sampled SOS with R = 12 showed good consistency (intraclass correlation coefficient > 0.91) with Cartesian reference for T1/PD/B1+/MWF/QSM/R2*, and good reproducibility (coefficient of variation < 7.0 %) in the test-retest analysis for T1/PD/B1+/MWF/R2*. This study has demonstrated the feasibility of simultaneous whole brain multiparametric mapping with a two-minute scan using the DFA-mGRE SOS sequence, which may overcome a major obstacle for neurological applications of multiparametric MRI.

Improving quantitative MRI using self-supervised deep learning with model reinforcement: Demonstration for rapid T1 mapping

PURPOSE

This paper proposes a novel self-supervised learning framework that uses model reinforcement, REference-free LAtent map eXtraction with MOdel REinforcement (RELAX-MORE), for accelerated quantitative MRI (qMRI) reconstruction. The proposed method uses an optimization algorithm to unroll an iterative model-based qMRI reconstruction into a deep learning framework, enabling accelerated MR parameter maps that are highly accurate and robust.

METHODS

Unlike conventional deep learning methods which require large amounts of training data, RELAX-MORE is a subject-specific method that can be trained on single-subject data through self-supervised learning, making it accessible and practically applicable to many qMRI studies. Using quantitativeT 1 $$ {\mathrm{T}}_1 $$ mapping as an example, the proposed method was applied to the brain, knee and phantom data.

RESULTS

The proposed method generates high-quality MR parameter maps that correct for image artifacts, removes noise, and recovers image features in regions of imperfect image conditions. Compared with other state-of-the-art conventional and deep learning methods, RELAX-MORE significantly improves efficiency, accuracy, robustness, and generalizability for rapid MR parameter mapping.

CONCLUSION

This work demonstrates the feasibility of a new self-supervised learning method for rapid MR parameter mapping, that is readily adaptable to the clinical translation of qMRI.

Fast and accessible T2 mapping using off-resonance corrected DESPOT2 with application to 3D prostate

PURPOSE

Most T1 and T2 mapping take long acquisitions or needs specialized sequences not widely accessible on clinical scanners. An available solution is DESPOT1/T2 (Driven equilibrium single pulse observation of T1/T2). DESPOT1/T2 uses Spoiled gradient-echo (SPGR) and balanced Steady-State Free Precession (bSSFP) sequences, offering an accessible and reliable way for 3D accelerated T1/T2 mapping. However, bSSFP is prone to off-resonance artifacts, limiting the application of DESPOT2 in regions with high susceptibility contrasts, like the prostate. Our proposal, DESPO+, employs the full bSSFP and SPGR models with a dictionary-based method to reconstruct 3D T1/T2 maps in the prostate region without off-resonance banding.

METHODS

DESPO+ modifies the bSSFP acquisition of the original variable flip angle DESPOT2. DESPO+ uses variable repetition and echo times, employing a dictionary-based method of the full bSSFP and SPGR models to reconstruct T1, T2, and Proton Density (PD) simultaneously. The proposed DESPO+ method underwent testing through simulations, T1/T2 phantoms, and on fourteen healthy subjects.

RESULTS

The results reveal a significant reduction in T2 map banding artifacts compared to the original DESPOT2 method. DESPO+ approach reduced T2 errors by up to seven times compared to DESPOT2 in simulations and phantom experiments. We also synthesized in-vivo T1-weighted/T2-weighted images from the acquired maps using a spin-echo model to verify the map's quality when lacking a reference. For in-vivo imaging, the synthesized images closely resemble those from the clinical MRI protocol, reducing scan time by around 50% compared to traditional spin-echo T1-weighted/T2-weighted acquisitions.

CONCLUSION

DESPO+ provides an off-resonance insensitive and clinically available solution, enabling high-resolution 3D T1/T2 mapping and synthesized T1-weighted/T2-weighted images for the entire prostate, all achieved within a short scan time of 3.6 min, similar to DESPOT1/T2.

Brain microstructure is linked to cognitive fatigue in early multiple sclerosis

Cognitive fatigue is a major symptom of Multiple Sclerosis (MS), from the early stages of the disease. This study aims to detect if brain microstructure is altered early in the disease course and is associated with cognitive fatigue in people with MS (pwMS) compared to matched healthy controls (HC). Recently diagnosed pwMS (N = 18, age < 45 years old) with either a Relapsing-Remitting or a Clinically Isolated Syndrome course of the disease, and HC (N = 19) matched for sex, age and education were analyzed. Quantitative multiparameter maps (MTsat, PD, R1 and R2*) of pwMS and HC were calculated. Parameters were extracted within the normal appearing white matter, cortical grey matter and deep grey matter (NAWM, NACGM and NADGM, respectively). Bayesian T-test for independent samples assessed between-group differences in brain microstructure while associations between score at a cognitive fatigue scale and each parameter in each tissue class were investigated with Generalized Linear Mixed Models. Patients exhibited lower MTsat and R1 values within NAWM and NACGM, and higher R1 values in NADGM compared to HC. Cognitive fatigue was associated with PD measured in every tissue class and to MTsat in NAWM, regardless of group. Disease-specific negative correlations were found in pwMS in NAWM (R1, R2*) and NACGM (R1). These findings suggest that brain microstructure within normal appearing tissues is already altered in the very early stages of the disease. Moreover, additional microstructure alterations (e.g. diffuse and widespread demyelination or axonal degeneration) in pwMS may lead to disease-specific complaint of cognitive fatigue.

Evidence for interictal blood-brain barrier dysfunction in people with epilepsy

OBJECTIVE

Interictal blood-brain barrier dysfunction in chronic epilepsy has been demonstrated in animal models and pathological specimens. Ictal blood-brain barrier dysfunction has been shown in humans in vivo using an experimental quantitative magnetic resonance imaging (MRI) protocol. Here, we hypothesized that interictal blood-brain barrier dysfunction is also present in people with drug-resistant epilepsy.

METHODS

Thirty-nine people (21 females, mean age at MRI ± SD = 30 ± 8 years) with drug-resistant epilepsy were prospectively recruited and underwent interictal T1-relaxometry before and after administration of a paramagnetic contrast agent. Likewise, quantitative T1 was acquired in 29 people without epilepsy (12 females, age at MRI = 48 ± 18 years). Quantitative T1 difference maps were calculated and served as a surrogate imaging marker for blood-brain barrier dysfunction. Values of quantitative T1 difference maps inside hemispheres ipsilateral to the presumed seizure onset zone were then compared, on a voxelwise level and within presumed seizure onset zones, to the contralateral side of people with epilepsy and to people without epilepsy.

RESULTS

Compared to the contralateral side, ipsilateral T1 difference values were significantly higher in white matter (corrected p < .05), gray matter (uncorrected p < .05), and presumed seizure onset zones (p = .04) in people with epilepsy. Compared to people without epilepsy, significantly higher T1 difference values were found in the anatomical vicinity of presumed seizure onset zones (p = .004). A subgroup of people with hippocampal sclerosis demonstrated significantly higher T1 difference values in the ipsilateral hippocampus and in regions strongly interconnected with the hippocampus compared to people without epilepsy (corrected p < .01). Finally, z-scores reflecting the deviation of T1 difference values within the presumed seizure onset zone were associated with verbal memory performance (p = .02) in people with temporal lobe epilepsy.

SIGNIFICANCE

Our results indicate a blood-brain barrier dysfunction in drug-resistant epilepsy that is detectable interictally in vivo, anatomically related to the presumed seizure onset zone, and associated with cognitive deficits.

Validation of single reference variable flip angle (SR-VFA) dynamic T1 mapping with T2 * correction using a novel rotating phantom

PURPOSE

To validate single reference variable flip angle (SR-VFA) dynamic T1 mapping with and without T2 * correction against inversion recovery (IR) T1 measurements.

METHODS

A custom cylindrical phantom with three concentric compartments was filled with variably doped agar to produce a smooth spatial gradient of the T1 relaxation rate as a function of angle across each compartment. IR T1 , VFA T1 , and B1 + measurements were made on the phantom before rotation, and multi-echo stack-of-radial dynamic images were acquired during rotation via an MRI-compatible motor. B1 + -corrected SR-VFA and SR-VFA-T2 * T1 maps were computed from the sliding window reconstructed images and compared against rotationally registered IR and VFA T1 maps to determine the percentage error.

RESULTS

Both VFA and SR-VFA-T2 * T1 maps fell within 10% of IR T1 measurements for a low rotational speed, with a mean accuracy of 2.3% ± 2.6% and 2.8% ± 2.6%, respectively. Increasing rotational speed was found to decrease the accuracy due to increasing temporal smoothing over ranges where the T1 change had a nonconstant slope. SR-VFA T1 mapping was found to have similar accuracy as the SR-VFA-T2 * and VFA methods at low TEs (˜<2 ms), whereas accuracy degraded strongly with later TEs. T2 * correction of the SR-VFA T1 maps was found to consistently improve accuracy and precision, especially at later TEs.

CONCLUSION

SR-VFA-T2 * dynamic T1 mapping was found to be accurate against reference IR T1 measurements within 10% in an agar phantom. Further validation is needed in mixed fat-water phantoms and in vivo.

Rapid quantitative magnetization transfer imaging: Utilizing the hybrid state and the generalized Bloch model

PURPOSE

To explore efficient encoding schemes for quantitative magnetization transfer (qMT) imaging with few constraints on model parameters.

THEORY AND METHODS

We combine two recently proposed models in a Bloch-McConnell equation: the dynamics of the free spin pool are confined to the hybrid state, and the dynamics of the semi-solid spin pool are described by the generalized Bloch model. We numerically optimize the flip angles and durations of a train of radio frequency pulses to enhance the encoding of three qMT parameters while accounting for all eight parameters of the two-pool model. We sparsely sample each time frame along this spin dynamics with a three-dimensional radial koosh-ball trajectory, reconstruct the data with subspace modeling, and fit the qMT model with a neural network for computational efficiency.

RESULTS

We extracted qMT parameter maps of the whole brain with an effective resolution of 1.24 mm from a 12.6-min scan. In lesions of multiple sclerosis subjects, we observe a decreased size of the semi-solid spin pool and longer relaxation times, consistent with previous reports.

CONCLUSION

The encoding power of the hybrid state, combined with regularized image reconstruction, and the accuracy of the generalized Bloch model provide an excellent basis for efficient quantitative magnetization transfer imaging with few constraints on model parameters.

Multimodal imaging of microstructural cerebral alterations and loss of synaptic density in Alzheimer's disease

Multiple neuropathological events are involved in Alzheimer's disease (AD). The current study investigated the concurrence of neurodegeneration, increased iron content, atrophy, and demyelination in AD. Quantitative multiparameter magnetic resonance imaging (MRI) maps providing neuroimaging biomarkers for myelination and iron content along with synaptic density measurements using [18F] UCB-H PET were acquired in 24 AD and 19 Healthy controls (19 males). The whole brain voxel-wise group comparison revealed demyelination in the right hippocampus, while no significant iron content difference was detected. Bilateral atrophy and synaptic density loss were observed in the hippocampus and amygdala. The multivariate GLM (mGLM) analysis shows a bilateral difference in the hippocampus and amygdala, right pallidum, left fusiform and temporal lobe suggesting that these regions are the most affected despite the diverse differences in brain tissue properties in AD. Demyelination was identified as the most affecting factor in the observed differences. Here, the mGLM is introduced as an alternative for multiple comparisons between different modalities, reducing the risk of false positives while informing about the co-occurrence of neuropathological processes in AD.

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.

B1 mapping using pre-learned subspaces for quantitative brain imaging

PURPOSE

To develop a machine learning-based method for estimation of both transmitter and receiver B1 fields desired for correction of the B1 inhomogeneity effects in quantitative brain imaging.

THEORY AND METHODS

A subspace model-based machine learning method was proposed for estimation of B1t and B1r fields. Probabilistic subspace models were used to capture scan-dependent variations in the B1 fields; the subspace basis and coefficient distributions were learned from pre-scanned training data. Estimation of the B1 fields for new experimental data was achieved by solving a linear optimization problem with prior distribution constraints. We evaluated the performance of the proposed method for B1 inhomogeneity correction in quantitative brain imaging scenarios, including T1 and proton density (PD) mapping from variable-flip-angle spoiled gradient-echo (SPGR) data as well as neurometabolic mapping from MRSI data, using phantom, healthy subject and brain tumor patient data.

RESULTS

In both phantom and healthy subject data, the proposed method produced high-quality B1 maps. B1 correction on SPGR data using the estimated B1 maps produced significantly improved T1 and PD maps. In brain tumor patients, the proposed method produced more accurate and robust B1 estimation and correction results than conventional methods. The B1 maps were also applied to MRSI data from tumor patients and produced improved neurometabolite maps, with better separation between pathological and normal tissues.

CONCLUSION

This work presents a novel method to estimate B1 variations using probabilistic subspace models and machine learning. The proposed method may make correction of B1 inhomogeneity effects more robust in practical applications.

Review of oxygen-enhanced lung mri: Pulse sequences for image acquisition and T1 measurement

Oxygen-enhanced MR imaging (OE-MRI) is a special proton imaging technique that can be performed without modifying the scanner hardware. Many fundamental studies have been conducted following the initial reporting of this technique in 1996, illustrating the high potential for its clinical application. This review aims to summarise and analyse current pulse sequences and T1 measurement methods for OE-MRI, including fundamental theories, existing pulse sequences applied to OE-MRI acquisition and T1 mapping. Wash-in and wash-out time identify lung function and are sensitive to ventilation; thus, dynamic OE-MRI is also discussed in this review. We compare OE-MRI with the primary competitive technique, hyperpolarised gas MRI. Finally, an overview of lower-field applications of OE-MRI is highlighted, as relatively recent publications demonstrated positive results. Lower-field OE-MRI, which is lower than 1.5 T, could be an alternative modality for detecting lung diseases. This educational review is aimed at researchers who want a quick summary of the steps needed to perform pulmonary OE-MRI with a particular focus on sequence design, settings, and quantification methods.

Optimal flip angles for in vivo liver 3D T1 mapping and B1+ mapping at 3T

PURPOSE

The spoiled gradient recalled echo (SPGR) sequence with variable flip angles (FAs) enables whole liverT 1 $$ {T}_1 $$ mapping at high spatial resolutions but is strongly affected byB 1 + $$ {B}_1^{+} $$ inhomogeneities. The aim of this work was to study how the precision of acquiredT 1 $$ {T}_1 $$ maps is affected by theT 1 $$ {T}_1 $$ andB 1 + $$ {B}_1^{+} $$ ranges observed in the liver at 3T, as well as how noise propagates from the acquired signals into the resultingT 1 $$ {T}_1 $$ map.

THEORY

TheT 1 $$ {T}_1 $$ variance was estimated through the Fisher information matrix with a total noise variance including, for the first time, theB 1 + $$ {B}_1^{+} $$ map noise as well as contributions from the SPGR noise.

METHODS

Simulations were used to find the optimal FAs for both theB 1 + $$ {B}_1^{+} $$ mapping andT 1 $$ {T}_1 $$ mapping. The simulations results were validated in 10 volunteers.

RESULTS

Four optimized SPGR FAs of 2°, 2°, 15°, and 15° (TR = 4.1 ms) andB 1 + $$ {B}_1^{+} $$ map FAs of 65° and 130° achieved aT 1 $$ {T}_1 $$ coefficient of variation of 6.2 ± 1.7% across 10 volunteers and validated our theoretical model. Four optimal FAs outperformed five uniformly spaced FAs, saving the patient one breath-hold. For the liverB 1 + $$ {B}_1^{+} $$ andT 1 $$ {T}_1 $$ parameter space at 3T, a higher return inT 1 $$ {T}_1 $$ precision was obtained by investing FAs in the SPGR acquisition rather than in theB 1 + $$ {B}_1^{+} $$ map.

CONCLUSION

A novel framework was developed and validated to calculate the SPGRT 1 $$ {T}_1 $$ variance. This framework efficiently identifies optimal FA values and determines the total number of SPGR andB 1 + $$ {B}_1^{+} $$ measurements needed to achieve a desiredT 1 $$ {T}_1 $$ precision.

Potential cognitive and neural benefits of a computerised cognitive training programme based on Structure Learning in healthy adults: study protocol for a randomised controlled trial

BACKGROUND

Cognitive flexibility refers to the capacity to shift between conceptual representations particularly in response to changes in instruction and feedback. It enables individuals to swiftly adapt to changes in their environment and has significant implications for learning. The present study focuses on investigating changes in cognitive flexibility following an intervention programme-Structure Learning training.

METHODS

Participants are pseudo-randomised to either the Training or Control group, while matched on age, sex, intelligence and cognitive flexibility performance. In the Training group, participants undergo around 2 weeks of training (at least 13 sessions) on Structure Learning. In the Control group, participants do not have to undergo any training and are never exposed to the Structure Learning task. The effects of Structure Learning training are investigated at both the behavioural and neural level. We measured covariates that can influence an individual's training performance before the training phase and outcome measures that can potentially show training benefits after the training phase. At the behavioural level, we investigated outcomes in both cognitive and social aspects with a primary focus on executive functions. At the neural level, we employed a multimodality approach and investigated potential changes to functional connectivity patterns, neurometabolite concentration in the frontal brain regions, and brain microstructure and myelination.

DISCUSSION

We reported the development of a novel training programme based on Structure Learning that aims to hone a general learning ability to potentially achieve extensive transfer benefits across various cognitive constructs. Potential transfer benefits can be exhibited through better performance in outcome measures between Training and Control participants, and positive associations between training performance and outcomes after the training in Training participants. Moreover, we attempt to substantiate behavioural findings with evidence of neural changes across different imaging modalities by the Structure Learning training.

TRIAL REGISTRATION

National Institutes of Health U.S. National Library of Medicine ClinicalTrials.gov NCT05611788. Registered on 7 November 2022.

PROTOCOL VERSION

11 May 2023.

Focal epilepsy without overt epileptogenic lesions: no evidence of microstructural brain tissue damage in multi-parametric quantitative MRI

Background and purpose

In patients with epilepsies of structural origin, brain atrophy and pathological alterations of the tissue microstructure extending beyond the putative epileptogenic lesion have been reported. However, in patients without any evidence of epileptogenic lesions on diagnostic magnetic resonance imaging (MRI), impairment of the brain microstructure has been scarcely elucidated. Using multiparametric quantitative (q) magnetic resonance imaging MRI, we aimed to investigate diffuse impairment of the microstructural tissue integrity in MRI-negative focal epilepsy patients.

Methods

27 MRI-negative patients with focal epilepsy (mean age 33.1 ± 14.2 years) and 27 matched healthy control subjects underwent multiparametric qMRI including T1, T2, and PD mapping at 3 T. After tissue segmentation based on synthetic anatomies, mean qMRI parameter values were extracted from the cerebral cortex, the white matter (WM) and the deep gray matter (GM) and compared between patients and control subjects. Apart from calculating mean values for the qMRI parameters across the respective compartments, voxel-wise analyses were performed for each tissue class.

Results

There were no significant differences for mean values of quantitative T1, T2, and PD obtained from the cortex, the WM and the deep GM between the groups. Furthermore, the voxel-wise analyses did not reveal any clusters indicating significant differences between patients and control subjects for the qMRI parameters in the respective compartments.

Conclusions

Based on the employed methodology, no indication for an impairment of the cerebral microstructural tissue integrity in MRI-negative patients with focal epilepsy was found in this study. Further research will be necessary to identify relevant factors and mechanisms contributing to microstructural brain tissue damage in various subgroups of patients with epilepsy.

Quantitative magnetization transfer MRI unbiased by on-resonance saturation and dipolar order contributions

PURPOSE

To demonstrate the bias in quantitative MT (qMT) measures introduced by the presence of dipolar order and on-resonance saturation (ONRS) effects using magnetization transfer (MT) spoiled gradient-recalled (SPGR) acquisitions, and propose changes to the acquisition and analysis strategies to remove these biases.

METHODS

The proposed framework consists of SPGR sequences prepared with simultaneous dual-offset frequency-saturation pulses to cancel out dipolar order and associated relaxation (T_1D ) effects in Z-spectrum acquisitions, and a matched quantitative MT (qMT) mathematical model that includes ONRS effects of readout pulses. Variable flip angle and MT data were fitted jointly to simultaneously estimate qMT parameters (macromolecular proton fraction [MPF], T_2,f , T_2,b , R, and free pool T₁ ). This framework is compared with standard qMT and investigated in terms of reproducibility, and then further developed to follow a joint single-point qMT methodology for combined estimation of MPF and T₁ .

RESULTS

Bland-Altman analyses demonstrated a systematic underestimation of MPF (-2.5% and -1.3%, on average, in white and gray matter, respectively) and overestimation of T₁ (47.1 ms and 38.6 ms, on average, in white and gray matter, respectively) if both ONRS and dipolar order effects are ignored. Reproducibility of the proposed framework is excellent (ΔMPF = -0.03% and ΔT₁  = -19.0 ms). The single-point methodology yielded consistent MPF and T₁ values with respective maximum relative average bias of -0.15% and -3.5 ms found in white matter.

CONCLUSION

The influence of acquisition strategy and matched mathematical model with regard to ONRS and dipolar order effects in qMT-SPGR frameworks has been investigated. The proposed framework holds promise for improved accuracy with reproducibility.

In vivo multi-parameter mapping of the habenula using MRI

The habenula is a small, epithalamic brain structure situated between the mediodorsal thalamus and the third ventricle. It plays an important role in the reward circuitry of the brain and is implicated in psychiatric conditions, such as depression. The importance of the habenula for human cognition and mental health make it a key structure of interest for neuroimaging studies. However, few studies have characterised the physical properties of the human habenula using magnetic resonance imaging because its challenging visualisation in vivo, primarily due to its subcortical location and small size. To date, microstructural characterization of the habenula has focused on quantitative susceptibility mapping. In this work, we complement this previous characterisation with measures of longitudinal and effective transverse relaxation rates, proton density and magnetisation transfer saturation using a high-resolution quantitative multi-parametric mapping protocol at 3T, in a cohort of 26 healthy participants. The habenula had consistent boundaries across the various parameter maps and was most clearly visualised on the longitudinal relaxation rate maps. We have provided a quantitative multi-parametric characterisation that may be useful for future sequence optimisation to enhance visualisation of the habenula, and additionally provides reference values for future studies investigating pathological differences in habenula microstructure.

Reliability of multi-parameter mapping (MPM) in the cervical cord: A multi-center multi-vendor quantitative MRI study

MRI based multicenter studies which target neurological pathologies affecting the spinal cord and brain - including spinal cord injury (SCI) - require standardized acquisition protocols and image processing methods. We have optimized and applied a multi-parameter mapping (MPM) protocol that simultaneously covers the brain and the cervical cord within a traveling heads study across six clinical centers (Leutritz et al., 2020). The MPM protocol includes quantitative maps (magnetization transfer saturation (MT), proton density (PD), longitudinal (R1), and effective transverse (R2*) relaxation rates) sensitive to myelination, water content, iron concentration, and morphometric measures, such as cross-sectional cord area. Previously, we assessed the repeatability and reproducibility of the brain MPM data acquired in the five healthy participants who underwent two scan-rescans (Leutritz et al., 2020). This study focuses on the cervical cord MPM data derived from the same acquisitions to determine its repeatability and reproducibility in the cervical cord. MPM matrices of the cervical cord were generated and processed using the hMRI and the spinal cord toolbox. To determine reliability of the cervical MPM data, the intra-site (i.e., scan-rescan) coefficient of variation (CoV), inter-site CoV, and bias within region of interests (C1, C2 and C3 levels) were determined. The range of the mean intra- and inter-site CoV of MT, R1 and PD was between 2.5% and 12%, and between 1.1% and 4.0% for the morphometric measures. In conclusion, the cervical MPM data showed a high repeatability and reproducibility for key imaging biomarkers and hence can be employed as a standardized tool in multi-center studies, including clinical trials.

First clinical application of a novel T1 mapping of the whole brain

BACKGROUND

The aim of this study was to evaluate the reproducibility and clinical value of the novel single-shot T1 mapping method for rapid and accurate multi-slice coverage of the whole brain, described by Wang et al. 2015.

METHODS

At a field strength of 3 Tesla, T1 mappings of 139 patients (51 of them without pathologic findings) and two repeats of five volunteers were performed at 0.5 mm in-plane resolution. Mean T1 values were determined in 18 manually segmented regions-of-interest without pathologic findings. Reproducibility of the repeated scans was calculated using mean coefficient of variations. Pathologies were grouped and separately evaluated.

RESULTS

The mean age of the cohort was 49 (range 1-95 years). T1 relaxation times for ordinary brain and pathologies were in accordance with the literature values. Intra- and inter-subject reproducibility was excellent, and mean coefficient of variations were 2.4% and 3.8%, respectively.

DISCUSSION

The novel rapid T1 mapping method is a reliable magnetic resonance imaging technique for identifying and quantifying normal brain structures and may thus serve as a basis for assessing pathologies. The fast and parallel online calculation enables a comfortable use in everyday clinical practice. We see a possible clinical value in a large spectrum of diseases, which should be investigated in further studies.

Reproducibility of rapid multi-parameter mapping at 3T and 7T with highly segmented and accelerated 3D-EPI

PURPOSE

Quantitative multi-parameter mapping (MPM) has been shown to provide good longitudinal and cross-sectional reproducibility for clinical research. Unfortunately, acquisition times (TAs) are typically infeasible for routine scanning at high resolutions.

METHODS

A fast whole-brain MPM protocol based on interleaved multi-shot 3D-EPI with controlled aliasing (SC-EPI) at 3T and 7T is proposed and compared with MPM using a standard spoiled gradient echo (FLASH) sequence. Four parameters (R₁ , PD, R 2 * $$ {R}_2^{\ast } $$ , and MTsat) were measured in less than 3 min at 1 mm isotropic resolution. Five subjects went through the same scanning sessions twice at each scanner. The intra-subject coefficient of variation (scan-rescan) (CoV) was estimated for each protocol and scanner to assess the longitudinal reproducibility.

RESULTS

At 3T, the CoV of SC-EPI ranged between 1.2%-4.8% for PD and R₁ , 2.8%-10.6% for R 2 * $$ {R}_2^{\ast } $$ and MTsat, which was comparable with FLASH (0.6%-4.9% for PD and R₁ , 2.6%-11.3% for R 2 * $$ {R}_2^{\ast } $$ and MTsat). At 7T, where the SC-EPI TA was reduced to ∼2 min, the CoV of SC-EPI (1.4%-10.6% for PD, R₁ , and R 2 * $$ {R}_2^{\ast } $$ ) was 1.2-2.4 times larger than the CoV of FLASH (1.0%-15%) and MTsat showed much higher variability across subjects. The SC-EPI-MPM protocol at 3T showed high reproducibility and yielded stable quantitative maps at a clinically feasible resolution and scan time, whereas at 7T, MT saturation homogeneity needs to be improved.

CONCLUSION

SC-EPI-based MPM is feasible as an additional MRI modality in clinical or population studies where the parameters offer great potential as biomarkers.

Vendor-neutral sequences and fully transparent workflows improve inter-vendor reproducibility of quantitative MRI

PURPOSE

We developed an end-to-end workflow that starts with a vendor-neutral acquisition and tested the hypothesis that vendor-neutral sequences decrease inter-vendor variability of T1, magnetization transfer ratio (MTR), and magnetization transfer saturation-index (MTsat) measurements.

METHODS

We developed and deployed a vendor-neutral 3D spoiled gradient-echo (SPGR) sequence on three clinical scanners by two MRI vendors. We then acquired T1 maps on the ISMRM-NIST system phantom, as well as T1, MTR, and MTsat maps in three healthy participants. We performed hierarchical shift function analysis in vivo to characterize the differences between scanners when the vendor-neutral sequence is used instead of commercial vendor implementations. Inter-vendor deviations were compared for statistical significance to test the hypothesis.

RESULTS

In the phantom, the vendor-neutral sequence reduced inter-vendor differences from 8% to 19.4% to 0.2% to 5% with an overall accuracy improvement, reducing ground truth T1 deviations from 7% to 11% to 0.2% to 4%. In vivo, we found that the variability between vendors is significantly reduced (p = 0.015) for all maps (T1, MTR, and MTsat) using the vendor-neutral sequence.

CONCLUSION

We conclude that vendor-neutral workflows are feasible and compatible with clinical MRI scanners. The significant reduction of inter-vendor variability using vendor-neutral sequences has important implications for qMRI research and for the reliability of multicenter clinical trials.

The Compressed Sensing MP2RAGE as a Surrogate to the MPRAGE for Neuroimaging at 3 T

OBJECTIVES

The magnetization-prepared 2 rapid acquisition gradient echo (MP2RAGE) sequence provides quantitative T1 maps in addition to high-contrast morphological images. Advanced acceleration techniques such as compressed sensing (CS) allow its acquisition time to be compatible with clinical applications. To consider its routine use in future neuroimaging protocols, the repeatability of the segmented brain structures was evaluated and compared with the standard morphological sequence (magnetization-prepared rapid gradient echo [MPRAGE]). The repeatability of the T1 measurements was also assessed.

MATERIALS AND METHODS

Thirteen healthy volunteers were scanned either 3 or 4 times at several days of interval, on a 3 T clinical scanner, with the 2 sequences (CS-MP2RAGE and MPRAGE), set with the same spatial resolution (0.8-mm isotropic) and scan duration (6 minutes 21 seconds). The reconstruction time of the CS-MP2RAGE outputs (including the 2 echo images, the MP2RAGE image, and the T1 map) was 3 minutes 33 seconds, using an open-source in-house algorithm implemented in the Gadgetron framework.Both precision and variability of volume measurements obtained from CAT12 and VolBrain were assessed. The T1 accuracy and repeatability were measured on phantoms and on humans and were compared with literature.Volumes obtained from the CS-MP2RAGE and the MPRAGE images were compared using Student t tests (P < 0.05 was considered significant).

RESULTS

The CS-MP2RAGE acquisition provided morphological images of the same quality and higher contrasts than the standard MPRAGE images. Similar intravolunteer variabilities were obtained with the CS-MP2RAGE and the MPRAGE segmentations. In addition, high-resolution T1 maps were obtained from the CS-MP2RAGE. T1 times of white and gray matters and several deep gray nuclei are consistent with the literature and show very low variability (<1%).

CONCLUSIONS

The CS-MP2RAGE can be used in future protocols to rapidly obtain morphological images and quantitative T1 maps in 3-dimensions while maintaining high repeatability in volumetry and relaxation times.

In vivo Estimation of Axonal Morphology From Magnetic Resonance Imaging and Electroencephalography Data

Purpose

We present a novel approach that allows the estimation of morphological features of axonal fibers from data acquired in vivo in humans. This approach allows the assessment of white matter microscopic properties non-invasively with improved specificity.

Theory

The proposed approach is based on a biophysical model of Magnetic Resonance Imaging (MRI) data and of axonal conduction velocity estimates obtained with Electroencephalography (EEG). In a white matter tract of interest, these data depend on (1) the distribution of axonal radius [P(r)] and (2) the g-ratio of the individual axons that compose this tract [g(r)]. P(r) is assumed to follow a Gamma distribution with mode and scale parameters, M and θ, and g(r) is described by a power law with parameters α and β.

Methods

MRI and EEG data were recorded from 14 healthy volunteers. MRI data were collected with a 3T scanner. MRI-measured g-ratio maps were computed and sampled along the visual transcallosal tract. EEG data were recorded using a 128-lead system with a visual Poffenberg paradigm. The interhemispheric transfer time and axonal conduction velocity were computed from the EEG current density at the group level. Using the MRI and EEG measures and the proposed model, we estimated morphological properties of axons in the visual transcallosal tract.

Results

The estimated interhemispheric transfer time was 11.72 ± 2.87 ms, leading to an average conduction velocity across subjects of 13.22 ± 1.18 m/s. Out of the 4 free parameters of the proposed model, we estimated θ – the width of the right tail of the axonal radius distribution – and β – the scaling factor of the axonal g-ratio, a measure of fiber myelination. Across subjects, the parameter θ was 0.40 ± 0.07 μm and the parameter β was 0.67 ± 0.02 μm^−α.

Conclusion

The estimates of axonal radius and myelination are consistent with histological findings, illustrating the feasibility of this approach. The proposed method allows the measurement of the distribution of axonal radius and myelination within a white matter tract, opening new avenues for the combined study of brain structure and function, and for in vivo histological studies of the human brain.

Reliability of quantitative multiparameter maps is high for magnetization transfer and proton density but attenuated for R1 and R2 * in healthy young adults

We investigate the reliability of individual differences of four quantities measured by magnetic resonance imaging‐based multiparameter mapping (MPM): magnetization transfer saturation (MT), proton density (PD), longitudinal relaxation rate (R₁), and effective transverse relaxation rate (R₂*). Four MPM datasets, two on each of two consecutive days, were acquired in healthy young adults. On Day 1, no repositioning occurred and on Day 2, participants were repositioned between MPM datasets. Using intraclass correlation effect decomposition (ICED), we assessed the contributions of session‐specific, day‐specific, and residual sources of measurement error. For whole‐brain gray and white matter, all four MPM parameters showed high reproducibility and high reliability, as indexed by the coefficient of variation (CoV) and the intraclass correlation (ICC). However, MT, PD, R₁, and R₂* differed markedly in the extent to which reliability varied across brain regions. MT and PD showed high reliability in almost all regions. In contrast, R₁ and R₂* showed low reliability in some regions outside the basal ganglia, such that the sum of the measurement error estimates in our structural equation model was higher than estimates of between‐person differences. In addition, in this sample of healthy young adults, the four MPM parameters showed very little variability over four measurements but differed in how well they could assess between‐person differences. We conclude that R₁ and R₂* might carry only limited person‐specific information in some regions of the brain in healthy young adults, and, by implication, might be of restricted utility for studying associations to between‐person differences in behavior in those regions.

Correcting inter-scan motion artifacts in quantitative R1 mapping at 7T

PURPOSE

Inter-scan motion is a substantial source of error in R 1 estimation methods based on multiple volumes, for example, variable flip angle (VFA), and can be expected to increase at 7T where B 1 fields are more inhomogeneous. The established correction scheme does not translate to 7T since it requires a body coil reference. Here we introduce two alternatives that outperform the established method. Since they compute relative sensitivities they do not require body coil images.

THEORY

The proposed methods use coil-combined magnitude images to obtain the relative coil sensitivities. The first method efficiently computes the relative sensitivities via a simple ratio; the second by fitting a more sophisticated generative model.

METHODS

R 1 maps were computed using the VFA approach. Multiple datasets were acquired at 3T and 7T, with and without motion between the acquisition of the VFA volumes. R 1 maps were constructed without correction, with the proposed corrections, and (at 3T) with the previously established correction scheme. The effect of the greater inhomogeneity in the transmit field at 7T was also explored by acquiring B 1 + maps at each position.

RESULTS

At 3T, the proposed methods outperform the baseline method. Inter-scan motion artifacts were also reduced at 7T. However, at 7T reproducibility only converged on that of the no motion condition if position-specific transmit field effects were also incorporated.

CONCLUSION

The proposed methods simplify inter-scan motion correction of R 1 maps and are applicable at both 3T and 7T, where a body coil is typically not available. The open-source code for all methods is made publicly available.

Multi-parameter quantitative mapping of R1, R2*, PD, and MTsat is reproducible when accelerated with Compressed SENSE

Multi-parameter mapping (MPM) magnetic resonance imaging (MRI) provides quantitative estimates of the longitudinal and effective transverse relaxation rates R1 and R2*, proton density (PD), and magnetization transfer saturation (MTsat). Thereby, MPM enables better comparability across sites and time than conventional weighted MRI. However, for MPM, several contrasts must be acquired, resulting in prolonged measurement durations and thus preventing MPM's application in clinical routines. State-of-the-art imaging acceleration techniques such as Compressed SENSE (CS), a combination of compressed sensing and sensitivity encoding, can be used to reduce the scan time of MPM. However, the accuracy and precision of the resulting quantitative parameter maps have not been systematically evaluated. In this study, we therefore investigated the effect of CS acceleration on the fidelity and reproducibility of MPM acquisitions. In five healthy volunteers and in a phantom, we compared MPM metrics acquired without imaging acceleration, with the standard acceleration (SENSE factor 2.5), and with Compressed SENSE with acceleration factors 4 and 6 using a 32-channel head coil. We evaluated the reproducibility and repeatability of accelerated MPM using data from three scan sessions in gray and white matter volumes-of-interest (VOIs). Accelerated MPM provided precise and accurate quantitative parameter maps. For most parameters, the results of the CS-accelerated protocols correlated more strongly with the non-accelerated protocol than the standard SENSE-accelerated protocols. Furthermore, for most VOIs and contrasts, coefficients of variation were lower when calculated from data acquired with different imaging accelerations within a single scan session than from data acquired in different scan sessions. These results suggest that MPM with Compressed SENSE acceleration factors up to at least 6 yields reproducible quantitative parameter maps that are highly comparable to those acquired without imaging acceleration. Compressed SENSE can thus be used to considerably reduce the scan duration of R1, R2*, PD, and MTsat mapping, and is highly promising for clinical applications of MPM.

MRI of focal cortical dysplasia

Focal cortical dysplasia (FCD) are histopathologically categorized in ILAE type I to III. Mild malformations of cortical development (mMCD) including those with oligodendroglial hyperplasia (MOGHE) are to be integrated into this classification yet. Only FCD type II have distinctive MRI and molecular genetics alterations so far. Subtle FCD including FCD type II located in the depth of a sulcus are often overlooked requiring the use of dedicated sequences (MP2RAGE, FLAWS, EDGE) and/or voxel (VBM)- or surface-based (SBM) postprocessing. The added value of 7 Tesla MRI has to be proven yet.

Relaxation time of brain tissue in the elderly assessed by synthetic MRI

BACKGROUND

Synthetic MRI (SyMRI) is a quantitative technique that allows measurements of T1 and T2 relaxation times (RTs). Brain RT evolution across lifespan is well described for the younger population. The aim was to study RTs of brain parenchyma in a healthy geriatric population in order to define the normal value of structures in this group population. Normal values for geriatric population could help find biomarker for age-related brain disease.

MATERIALS AND METHODS

Fifty-four normal-functioning individuals (22 females, 32 males) with mean age of 83 years (range 56-98) underwent SyMRI. RT values in manually defined ROIs (centrum semiovale, middle cerebellar peduncles, thalamus, and insular cortex) and in segmented whole-brain components (brain parenchyma, gray matter, white matter, myelin, CSF, and stromal structures) were extracted from the SyMRI segmentation software. Patients' results were combined into the group age. Main ROI-based and whole-brain results were compared for the all dataset and for age group results as well.

RESULTS

For white matter, RTs between ROI-based analyses and whole-brain results for T2 and for T1 were statistically different and a trend of increasing T1 in centrum semiovale and cerebellar peduncle was observed. For gray matter, thalamic T1 was statistically different from insular T1. A difference was also found between left and right insula (p < .0001). T1 RTs of ROI-based and whole-brain-based analyses were statistically different (p < .0001). No significant difference in T1 and T2 was found between age groups on ROI-based analysis, but T1 in centrum semiovale and thalamus increased with age. No statistical difference between age groups was found for the various segmented volumes except for myelin between 65-74 years of age and the 95-105 years of age groups (p = .038).

CONCLUSIONS

SyMRI is a new tool that allows faster imaging and permits to obtain quantitative T1 and T2. By defining RT values of different brain components of normal-functioning elderly individuals, this technique may be used as a biomarker for clinical disorders like dementia.

Post mortem brain temperature and its influence on quantitative MRI of the brain

Objective

MRI temperature sensitivity presents a major issue in in situ post mortem MRI (PMMRI), as the tissue temperatures differ from living persons due to passive cooling of the deceased. This study aims at computing brain temperature effects on the MRI parameters to correct for temperature in PMMRI, laying the foundation for future projects on post mortem validation of in vivo MRI techniques.

Materials and methods

Brain MRI parameters were assessed in vivo and in situ post mortem using a 3 T MRI scanner. Post mortem brain temperature was measured in situ transethmoidally. The temperature effect was computed by fitting a linear model to the MRI parameters and the corresponding brain temperature.

Results

Linear positive temperature correlations were observed for T₁, T₂* and mean diffusivity in all tissue types. A significant negative correlation was observed for T₂ in white matter. Fractional anisotropy revealed significant correlations in all gray matter regions except for the thalamus.

Discussion

The linear models will allow to correct for temperature in post mortem MRI. Comparing in vivo to post mortem conditions, the mean diffusivity, in contrast to T₁ and T₂, revealed additional effects besides temperature, such as cessation of perfusion and active diffusion.

Model-based multi-parameter mapping

Quantitative MR imaging is increasingly favoured for its richer information content and standardised measures. However, computing quantitative parameter maps, such as those encoding longitudinal relaxation rate (R1), apparent transverse relaxation rate (R2*) or magnetisation-transfer saturation (MTsat), involves inverting a highly non-linear function. Many methods for deriving parameter maps assume perfect measurements and do not consider how noise is propagated through the estimation procedure, resulting in needlessly noisy maps. Instead, we propose a probabilistic generative (forward) model of the entire dataset, which is formulated and inverted to jointly recover (log) parameter maps with a well-defined probabilistic interpretation (e.g., maximum likelihood or maximum a posteriori). The second order optimisation we propose for model fitting achieves rapid and stable convergence thanks to a novel approximate Hessian. We demonstrate the utility of our flexible framework in the context of recovering more accurate maps from data acquired using the popular multi-parameter mapping protocol. We also show how to incorporate a joint total variation prior to further decrease the noise in the maps, noting that the probabilistic formulation allows the uncertainty on the recovered parameter maps to be estimated. Our implementation uses a PyTorch backend and benefits from GPU acceleration. It is available at https://github.com/balbasty/nitorch.

Sensitivity and resolution improvement for in vivo magnetic resonance current-density imaging of the human brain

PURPOSE

Magnetic resonance current-density imaging (MRCDI) combines MRI with low-intensity transcranial electrical stimulation (TES; 1-2 mA) to map current flow in the brain. However, usage of MRCDI is still hampered by low measurement sensitivity and image quality.

METHODS

Recently, a multigradient-echo-based MRCDI approach has been introduced that presently has the best-documented efficiency. This MRCDI approach has now been advanced in three directions and has been validated by phantom and in vivo experiments. First, the importance of optimum spoiling for brain imaging was verified. Second, the sensitivity and spatial resolution were improved by using acquisition weighting. Third, navigators were added as a quality control measure for tracking physiological noise. Combining these advancements, the optimized MRCDI method was tested by using 1 mA TES for two different injection profiles.

RESULTS

For a session duration of 4:20 min, the new MRCDI method was able to detect TES-induced magnetic fields at a sensitivity level of 84 picotesla, representing a twofold efficiency increase against our original method. A comparison between measurements and simulations based on personalized head models showed a consistent increase in the coefficient of determination of ΔR² = 0.12 for the current-induced magnetic fields and ΔR² = 0.22 for the current flow reconstructions. Interestingly, some of the simulations still clearly deviated from the measurements despite the strongly improved measurement quality. This highlights the utility of MRCDI to improve head models for TES simulations.

CONCLUSION

The achieved sensitivity improvement is an important step from proof-of-concept studies toward a broader application of MRCDI in clinical and basic neuroscience research.

Age-Dependent Changes in Total and Free Water Content of In Vivo Human Lenses Measured by Magnetic Resonance Imaging

Purpose

To use magnetic resonance imaging (MRI) to measure age-dependent changes in total and free water in human lenses in vivo.

Methods

Sixty-four healthy adults aged 18 to 86 years were recruited, fitted with a 32-channel head receiver coil, and placed in a 3 Tesla clinical MR scanner. Scans of the crystalline lens were obtained using a volumetric interpolated breath-hold examination sequence with dual flip angles, which were corrected for field inhomogeneity post-acquisition using a B1-map obtained using a turbo-FLASH sequence. The spatial distribution and content of corrected total (ρlens) and free (T1) water along the lens optical axis were extracted using custom-written code.

Results

Lens total water distribution and content did not change with age (all P > 0.05). In contrast to total water, a gradient in free water content that was highest in the periphery relative to the center was present in lenses across all ages. However, this initially parabolic free water gradient gradually developed an enhanced central plateau, as indicated by increasing profile shape parameter values (anterior: 0.067/y, P = 0.004; posterior: 0.050/y, P = 0.020) and central free water content (1.932 ms/y, P = 0.022) with age.

Conclusions

MRI can obtain repeatable total and free water measurements of in vivo human lenses. The observation that the lens steady-state free, but not total, water gradient is abolished with age raises the possibility that alterations in protein-water interactions are an underlying cause of the degradation in lens optics and overall vision observed with aging.

A fast and novel method for amide proton transfer-chemical exchange saturation transfer multislice imaging

Amide proton transfer-chemical exchange saturation transfer (APT-CEST) imaging provides important information for the diagnosis and monitoring of tumors. For such analysis, complete coverage of the brain is advantageous, especially when registration is performed with other magnetic resonance (MR) modalities, such as MR spectroscopy (MRS). However, the acquisition of Z-spectra across several slices via multislice imaging may be time-consuming. Therefore, in this paper, we present a new approach for fast multislice imaging, allowing us to acquire 16 slices per frequency offset within 8 s. The proposed fast CEST-EPI sequence employs a presaturation module, which drives the magnetization into the steady-state equilibrium for the first frequency offset. A second module, consisting of a single CEST pulse (for maintaining the steady-state) followed by an EPI acquisition, passes through a loop to acquire multiple slices and adjacent frequency offsets. Thus, the whole Z-spectrum can be recorded much faster than the conventional saturation scheme, which employs a presaturation for each single frequency offset. The validation of the CEST sequence parameters was performed by using the conventional saturation scheme. Subsequently, the proposed and a modified version of the conventional CEST sequence were compared in vitro on a phantom with different T1 times and in vivo on a brain tumor patient. No significant differences between both sequences could be found in vitro. The in vivo data yielded almost identical MTR_asym contrasts for the white and gray matter as well as for tumor tissue. Our results show that the proposed fast CEST-EPI sequence allows for rapid data acquisition and provides similar CEST contrasts as the modified conventional scheme while reducing the scanning time by approximately 50%.

RESUMEN: A flexible class of multi-parameter qMRI protocols

PURPOSE

To introduce a class of fast 3D quantitative MRI (qMRI) schemes (RESUME^N, for N=1,…,4) that allow for a thorough characterization of microstructural properties of brain tissues.

METHODS

An arbitrary multi-echo GRE acquisition optimized for quantitative susceptibility mapping (QSM) is complemented with an appropriate low flip-angle GRE sequence drawn from four possible choices. The acquired signals are processed to analytically derive the longitudinal relaxation (R₁) and free induction decay (R₂^∗) rates, as well as the proton density (PD) and QSM. A comprehensive modeling of the excitation and B₁^- profiles and of the RF-spoiling is included in the acquisition and processing pipeline.

RESULTS

The RESUME^N maps appear homogeneous throughout the field-of-view and exhibit comparable values and high SNR across the considered range of N values.

CONCLUSIONS

The introduced schemes represent a class of robust and flexible strategies to derive a thorough and fast qMRI study, suitable for a whole-brain acquisition with isotropic voxel resolution of 700 μm in less than 15 min.

System-specific evaluation of the dual flip angle MRI technique for quantitative T1 measurement

PURPOSE

To investigate if the accuracy of the dual flip angle (DFA) technique for T₁ measurement is affected by the system-specific RF excitation performance.

METHODS

A T₁ phantom, made of 12 vials of unique T₁ value ranging approximately from 200 ms to 2000 ms, was built and tested on seven different clinical scanners. For each experiment, the reference T₁ of each vial was obtained by the inversion recovery-based technique, and the DFA technique was applied repeatedly with several flip angle (FA) pairs conventionally proposed as optimal. The accuracy of the DFA technique for each FA pair was then evaluated by comparing the measured T₁ values for the vials to the references. Any variation of the accuracy was then evaluated across different FA pairs, and across different MRI systems. To improve accuracy with a selected FA pair, the signal ratio (SR) curve, obtained from the phantom, was utilized in a calibration strategy of the DFA technique.

RESULTS

When combined for all the vials, the average ratio of the measured T₁ to the reference generally increased as the FA pair window gradually slid from the smaller to the larger FA values. Furthermore, among several optimal FA pairs, the pair of the best accuracy varied slightly by the MRI system. The accuracy for any FA pair could be improved when the calibration strategy was utilized.

CONCLUSIONS

The RF excitation performance may vary by the specific FA pair and by the specific MRI system, influencing the accuracy of the DFA technique. The system-specific evaluation, and, if needed, its calibration, would help improve the accuracy of the DFA technique.

Model-based reconstruction for simultaneous multi-slice T1 mapping using single-shot inversion-recovery radial FLASH

PURPOSE

To develop a single-shot multi-slice T 1 mapping method by combing simultaneous multi-slice (SMS) excitations, single-shot inversion-recovery (IR) radial fast low-angle shot (FLASH), and a nonlinear model-based reconstruction method.

METHODS

SMS excitations are combined with a single-shot IR radial FLASH sequence for data acquisition. A previously developed single-slice calibrationless model-based reconstruction is extended to SMS, formulating the estimation of parameter maps and coil sensitivities from all slices as a single nonlinear inverse problem. Joint-sparsity constraints are further applied to the parameter maps to improve T 1 precision. Validations of the proposed method are performed for a phantom and for the human brain and liver in 6 healthy adult subjects.

RESULTS

Phantom results confirm good T 1 accuracy and precision of the simultaneously acquired multi-slice T 1 maps in comparison to single-slice references. In vivo human brain studies demonstrate the better performance of SMS acquisitions compared to the conventional spoke-interleaved multi-slice acquisition using model-based reconstruction. Aside from good accuracy and precision, the results of 6 healthy subjects in both brain and abdominal studies confirm good repeatability between scan and re-scans. The proposed method can simultaneously acquire T 1 maps for 5 slices of a human brain ( 0.75 × 0.75 × 5 mm 3 ) or 3 slices of the abdomen ( 1.25 × 1.25 × 6 mm 3 ) within 4 seconds.

CONCLUSIONS

The IR SMS radial FLASH acquisition together with a nonlinear model-based reconstruction enable rapid high-resolution multi-slice T 1 mapping with good accuracy, precision, and repeatability.

Imperfect spoiling in variable flip angle T1 mapping at 7T: Quantifying and minimizing impact

Purpose

The variable flip angle (VFA) approach to T₁ mapping assumes perfectly spoiled transverse magnetisation at the end of each repetition time (TR). Despite radiofrequency (RF) and gradient spoiling, this condition is rarely met, leading to erroneous T₁ estimates (T1app). Theoretical corrections can be applied but make assumptions about tissue properties, for example, a global T₂ time. Here, we investigate the effect of imperfect spoiling at 7T and the interaction between the RF and gradient spoiling conditions, additionally accounting for diffusion. We provide guidance on the optimal approach to maximise the accuracy of the T₁ estimate in the context of 3D multi‐echo acquisitions.

Methods

The impact of the spoiling regime was investigated through numerical simulations, phantom and invivo experiments.

Results

The predicted dependence of T1app on tissue properties, system settings, and spoiling conditions was observed in both phantom and in vivo experiments. Diffusion effects modulated the dependence of T1app on both B1+ efficiency and T₂ times.

Conclusion

Error in T1app can be minimized by using an RF spoiling increment and gradient spoiler moment combination that minimizes T₂‐dependence and safeguards image quality. Although the diffusion effect was comparatively small at 7T, correction factors accounting for this effect are recommended.

Click here for author‐reader discussions

The traveling heads 2.0: Multicenter reproducibility of quantitative imaging methods at 7 Tesla

OBJECT

This study evaluates inter-site and intra-site reproducibility at ten different 7 T sites for quantitative brain imaging.

MATERIAL AND METHODS

Two subjects - termed the "traveling heads" - were imaged at ten different 7 T sites with a harmonized quantitative brain MR imaging protocol. In conjunction with the system calibration, MP2RAGE, QSM, CEST and multi-parametric mapping/relaxometry were examined.

RESULTS

Quantitative measurements with MP2RAGE showed very high reproducibility across sites and subjects, and errors were in concordance with previous results and other field strengths. QSM had high inter-site reproducibility for relevant subcortical volumes. CEST imaging revealed systematic differences between the sites, but reproducibility was comparable to results in the literature. Relaxometry had also very high agreement between sites, but due to the high sensitivity, differences caused by different applications of the B1 calibration of the two RF coil types used were observed.

CONCLUSION

Our results show that quantitative brain imaging can be performed with high reproducibility at 7 T and with similar reliability as found at 3 T for multicenter studies of the supratentorial brain.

Mapping grip force to motor networks

AIM

There is ongoing debate about the role of cortical and subcortical brain areas in force modulation. In a whole-brain approach, we sought to investigate the anatomical basis of grip force whilst acknowledging interindividual differences in connectivity patterns. We tested if brain lesion mapping in patients with unilateral motor deficits can inform whole-brain structural connectivity analysis in healthy controls to uncover the networks underlying grip force.

METHODS

Using magnetic resonance imaging (MRI) and whole-brain voxel-based morphometry in chronic stroke patients (n=55) and healthy controls (n=67), we identified the brain regions in both grey and white matter significantly associated with grip force strength. The resulting statistical parametric maps (SPMs) provided seed areas for whole-brain structural covariance analysis in a large-scale community dwelling cohort (n=977) that included beyond volume estimates, parameter maps sensitive to myelin, iron and tissue water content.

RESULTS

The SPMs showed symmetrical bilateral clusters of correlation between upper limb motor performance, basal ganglia, posterior insula and cortico-spinal tract. The covariance analysis with the seed areas derived from the SPMs demonstrated a widespread anatomical pattern of brain volume and tissue properties, including both cortical, subcortical nodes of motor networks and sensorimotor areas projections.

CONCLUSION

We interpret our covariance findings as a biological signature of brain networks implicated in grip force. The data-driven definition of seed areas obtained from chronic stroke patients showed overlapping structural covariance patterns within cortico-subcortical motor networks across different tissue property estimates. This cumulative evidence lends face validity of our findings and their biological plausibility.

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.

Brain plasticity dynamics during tactile Braille learning in sighted subjects: Multi-contrast MRI approach

A growing body of empirical evidence supports the notion of diverse neurobiological processes underlying learning-induced plasticity changes in the human brain. There are still open questions about how brain plasticity depends on cognitive task complexity, how it supports interactions between brain systems and with what temporal and spatial trajectory. We investigated brain and behavioural changes in sighted adults during 8-months training of tactile Braille reading whilst monitoring brain structure and function at 5 different time points. We adopted a novel multivariate approach that includes behavioural data and specific MRI protocols sensitive to tissue properties to assess local functional and structural and myelin changes over time. Our results show that while the reading network, located in the ventral occipitotemporal cortex, rapidly adapts to tactile input, sensory areas show changes in grey matter volume and intra-cortical myelin at different times. This approach has allowed us to examine and describe neuroplastic mechanisms underlying complex cognitive systems and their (sensory) inputs and (motor) outputs differentially, at a mesoscopic level.

Voxel-Based quantitative MRI reveals spatial patterns of grey matter alteration in multiple sclerosis

Despite robust postmortem evidence and potential clinical importance of gray matter (GM) pathology in multiple sclerosis (MS), assessing GM damage by conventional magnetic resonance imaging (MRI) remains challenging. This prospective cross‐sectional study aimed at characterizing the topography of GM microstructural and volumetric alteration in MS using, in addition to brain atrophy measures, three quantitative MRI (qMRI) parameters—magnetization transfer (MT) saturation, longitudinal (R1), and effective transverse (R2*) relaxation rates, derived from data acquired during a single scanning session. Our study involved 35 MS patients (14 relapsing–remitting MS; 21 primary or secondary progressive MS) and 36 age‐matched healthy controls (HC). The qMRI maps were computed and segmented in different tissue classes. Voxel‐based quantification (VBQ) and voxel‐based morphometry (VBM) statistical analyses were carried out using multiple linear regression models. In MS patients compared with HC, three configurations of GM microstructural/volumetric alterations were identified. (a) Co‐localization of GM atrophy with significant reduction of MT, R1, and/or R2*, usually observed in primary cortices. (b) Microstructural modifications without significant GM loss: hippocampus and paralimbic cortices, showing reduced MT and/or R1 values without significant atrophy. (c) Atrophy without significant change in microstructure, identified in deep GM nuclei. In conclusion, this quantitative multiparametric voxel‐based approach reveals three different spatially‐segregated combinations of GM microstructural/volumetric alterations in MS that might be associated with different neuropathology.

Mapping Human Laryngeal Motor Cortex during Vocalization

The representations of the articulators involved in human speech production are organized somatotopically in primary motor cortex. The neural representation of the larynx, however, remains debated. Both a dorsal and a ventral larynx representation have been previously described. It is unknown, however, whether both representations are located in primary motor cortex. Here, we mapped the motor representations of the human larynx using functional magnetic resonance imaging and characterized the cortical microstructure underlying the activated regions. We isolated brain activity related to laryngeal activity during vocalization while controlling for breathing. We also mapped the articulators (the lips and tongue) and the hand area. We found two separate activations during vocalization-a dorsal and a ventral larynx representation. Structural and quantitative neuroimaging revealed that myelin content and cortical thickness underlying the dorsal, but not the ventral larynx representation, are similar to those of other primary motor representations. This finding confirms that the dorsal larynx representation is located in primary motor cortex and that the ventral one is not. We further speculate that the location of the ventral larynx representation is in premotor cortex, as seen in other primates. It remains unclear, however, whether and how these two representations differentially contribute to laryngeal motor control.

Characterization of 3-Dimensional Printing and Casting Materials for use in Magnetic Resonance Imaging Phantoms at 3 T

Imaging phantoms are used to calibrate and validate the performance of magnetic resonance imaging (MRI) systems. Many new materials have been developed for additive manufacturing (three-dimensional [3D] printing) processes that may be useful in the direct printing or casting of dimensionally accurate, anatomically accurate, patient-specific, and/or biomimetic MRI phantoms. The T1, T2, and T2* spin relaxation times of polymer samples were tested to discover materials for use as tissue mimics and structures in MRI phantoms. This study included a cohort of polymer compounds that was tested in cured form. The cohort consisted of 101 standardized polymer samples fabricated from: two-part silicones and polyurethanes used in commercial casting processes; one-part optically cured polyurethanes used in 3D printing; and fused deposition thermoplastics used in 3D printing. The testing was performed at 3 T using inversion recovery, spin echo, and gradient echo sequences for T1, T2, and T2*, respectively. T1, T2, and T2* values were plotted with error bars to allow the reader to assess how well a polymer matches a tissue for a specific application. A correlation was performed between T1, T2, T2* values and material density, elongation, tensile strength, and hardness. Two silicones, SI_XP-643 and SI_P-45, may be usable mimics for reported liver values; one silicone, SI_XP-643, may be a useful mimic for muscle; one silicone, SI_XP-738, may be a useful mimic for white matter; and four silicones, SI_P-15, SI_GI-1000, SI_GI-1040, and SI_GI-1110, may be usable mimics for spinal cord. Elongation correlated to T2 (p = 0.0007), tensile strength correlated to T1 (p = 0.002), T2 (p = 0.0003), and T2* (p = 0.003). The 80 samples not providing measurable signal with T1, T2, T2* relaxation values too short to measure with the standard sequences, may be useful for MRI-invisible fixturing and medical devices at 3 T.

Multiparameter quantitative histological MRI values in high-grade gliomas: a potential biomarker of tumor progression

Background

Conventional MRI poorly distinguishes brain parenchyma microscopically invaded by high-grade gliomas (HGGs) from the normal brain. By contrast, quantitative histological MRI (hMRI) measures brain microstructure in terms of physical MR parameters influenced by histochemical tissue composition. We aimed to determine the relationship between hMRI parameters in the area surrounding the surgical cavity and the presence of HGG recurrence.

Methods

Patients were scanned after surgery with an hMRI multiparameter protocol that allowed for estimations of longitudinal relaxation rate (R1) = 1/T1, effective transverse relaxation rate (R2)*=1/T2*, magnetization transfer saturation (MT_sat), and proton density. The initial perioperative zone (IPZ) was segmented on the postoperative MRI. Once recurrence appeared on conventional MRI, the area of relapsing disease was delineated (extension zone, EZ). Conventional MRI showing recurrence and hMRI were coregistered, allowing for the extraction of parameters R1, R2*, MT_sat, and PD in 3 areas: the overlap area between the IPZ and EZ (OZ), the peritumoral brain zone, PBZ (PBZ = IPZ - OZ), and the area of recurrence (RZ = EZ - OZ).

Results

Thirty-one patients with HGG who underwent gross-total resection were enrolled. MT_sat and R1 were the most strongly associated with tumor progression. MT_sat was significantly lower in the OZ and RZ, compared to PBZ. R1 was significantly lower in RZ compared to PBZ. PD was significantly higher in OZ compared to PBZ, and R2* was higher in OZ compared to PBZ or RZ. These changes were detected 4 to 120 weeks before recurrence recognition on conventional MRI.

Conclusions

HGG recurrence was associated with hMRI parameters' variation after initial surgery, weeks to months before overt recurrence.