3D magnetic resonance fingerprinting on a low-field 50 mT point-of-care system prototype: evaluation of muscle and lipid relaxation time mapping and comparison with standard techniques

OBJECTIVE

To implement magnetic resonance fingerprinting (MRF) on a permanent magnet 50 mT low-field system deployable as a future point-of-care (POC) unit and explore the quality of the parameter maps.

MATERIALS AND METHODS

3D MRF was implemented on a custom-built Halbach array using a slab-selective spoiled steady-state free precession sequence with 3D Cartesian readout. Undersampled scans were acquired with different MRF flip angle patterns and reconstructed using matrix completion and matched to the simulated dictionary, taking excitation profile and coil ringing into account. MRF relaxation times were compared to that of inversion recovery (IR) and multi-echo spin echo (MESE) experiments in phantom and in vivo. Furthermore, B₀ inhomogeneities were encoded in the MRF sequence using an alternating TE pattern, and the estimated map was used to correct for image distortions in the MRF images using a model-based reconstruction.

RESULTS

Phantom relaxation times measured with an optimized MRF sequence for low field were in better agreement with reference techniques than for a standard MRF sequence. In vivo muscle relaxation times measured with MRF were longer than those obtained with an IR sequence (T₁: 182 ± 21.5 vs 168 ± 9.89 ms) and with an MESE sequence (T₂: 69.8 ± 19.7 vs 46.1 ± 9.65 ms). In vivo lipid MRF relaxation times were also longer compared with IR (T₁: 165 ± 15.1 ms vs 127 ± 8.28 ms) and with MESE (T₂: 160 ± 15.0 ms vs 124 ± 4.27 ms). Integrated ΔB₀ estimation and correction resulted in parameter maps with reduced distortions.

DISCUSSION

It is possible to measure volumetric relaxation times with MRF at 2.5 × 2.5 × 3.0 mm³ resolution in a 13 min scan time on a 50 mT permanent magnet system. The measured MRF relaxation times are longer compared to those measured with reference techniques, especially for T₂. This discrepancy can potentially be addressed by hardware, reconstruction and sequence design, but long-term reproducibility needs to be further improved.

Quantitative MRI in leukodystrophies

Leukodystrophies constitute a large and heterogeneous group of genetic diseases primarily affecting the white matter of the central nervous system. Different disorders target different white matter structural components. Leukodystrophies are most often progressive and fatal. In recent years, novel therapies are emerging and for an increasing number of leukodystrophies trials are being developed. Objective and quantitative metrics are needed to serve as outcome measures in trials. Quantitative MRI yields information on microstructural properties, such as myelin or axonal content and condition, and on the chemical composition of white matter, in a noninvasive fashion. By providing information on white matter microstructural involvement, quantitative MRI may contribute to the evaluation and monitoring of leukodystrophies. Many distinct MR techniques are available at different stages of development. While some are already clinically applicable, others are less far developed and have only or mainly been applied in healthy subjects. In this review, we explore the background, current status, potential and challenges of available quantitative MR techniques in the context of leukodystrophies.

Fast and accurate T2 mapping using Bloch simulations and low-rank plus sparse matrix decomposition

PURPOSE

MRI's T2 relaxation time is one of the key contrast mechanisms for clinical diagnosis and prognosis of pathologies. Mapping this relaxation time, however, involves extensive scan times, which are needed to collect quantitative data, thereby impeding its integration into clinical routine. This study employs a low-rank plus sparse (L + S) signal decomposition approach in order to reconstruct accurate T2-maps from highly undersampled multi-echo spin-echo (MESE) MRI data.

METHODS

Two new algorithms are presented: the first uses standard L + S approach, where both L and S are iteratively updated. The second technique, dubbed SPArse and fixed RanK (SPARK), uses a fixed-rank L, under the assumption that most MESE information is found in the L component and that this rank can be pre-calculated. The utility of these new techniques is demonstrated on in vivo brain and calf data at x2 to x6 acceleration factors.

RESULTS

Accelerated T2 maps showed improved accuracy compared to fully sampled ground truth maps, when using L + S and SPARK techniques vis-à-vis standard GRAPPA acceleration.

CONCLUSION

SPARK provides accurate T2 maps with increased robustness to the selection of reconstruction parameters making it suitable to a wide range of applications and facilitating the use of quantitative T2 information in clinical settings.

Assessing post-traumatic changes in cartilage using T1ρ dispersion parameters

Degeneration of cartilage can be studied non-invasively with quantitative MRI. A promising parameter for detecting early osteoarthritis in articular cartilage is T_1ρ, which can be tuned via the amplitude of the spin-lock pulse. By measuring T_1ρ at several spin-lock amplitudes, the dispersion of T_1ρ is obtained. The aim of this study is to find out if the dispersion contains diagnostically relevant information complementary to a T_1ρ measurement at a single spin-lock amplitude. To this end, five differently acquired dispersion parameters are utilized; A, B, τ_c, T_1ρ/T₂, and R₂ - R_1ρ. An open dataset of an equine model of post-traumatic cartilage was utilized to assess the T_1ρ dispersion parameters for the evaluation of cartilage degeneration. Firstly, the parameters were compared for their sensitivity in detecting degenerative changes. Secondly, the relationship of the dispersion parameters to histological and biomechanical reference parameters was studied. Parameters A, T_1ρ/T₂, and R₂ - R_1ρ were found to be sensitive to lesion-induced changes in the cartilage within sample. Strong correlations of several dispersion parameters with optical density, as well as with collagen fibril angle were found. Most of the dispersion parameters correlated strongly with individual T_1ρ values. The results suggest that dispersion parameters can in some cases provide a more accurate description of the biochemical composition of cartilage as compared to conventional MRI parameters. However, in most cases the information given by the dispersion parameters is more of a refinement than complementary to conventional quantitative MRI.

Reliability of spinal cord measures based on synthetic T1-weighted MRI derived from multiparametric mapping (MPM)

Short MRI acquisition time, high signal-to-noise ratio, and high reliability are crucial for image quality when scanning healthy volunteers and patients. Cross-sectional cervical cord area (CSA) has been suggested as a marker of neurodegeneration and potential outcome measure in clinical trials and is conventionally measured on T₁-weigthed 3D Magnetization Prepared Rapid Acquisition Gradient-Echo (MPRAGE) images. This study aims to reduce the acquisition time for the comprehensive assessment of the spinal cord, which is typically based on MPRAGE for morphometry and multi-parameter mapping (MPM) for microstructure. The MPRAGE is replaced by a synthetic T₁-w MRI (synT₁-w) estimated from the MPM, in order to measure CSA. SynT₁-w images were reconstructed using the MPRAGE signal equation based on quantitative maps of proton density (PD), longitudinal (R₁) and effective transverse (R₂*) relaxation rates. The reliability of CSA measurements from synT₁-w images was determined within a multi-center test-retest study format and validated against acquired MPRAGE scans by assessing the agreement between both methods. The response to pathological changes was tested by longitudinally measuring spinal cord atrophy following spinal cord injury (SCI) for synT₁-w and MPRAGE using linear mixed effect models. CSA measurements based on the synT₁-w MRI showed high intra-site (Coefficient of variation [CoV]: 1.43% to 2.71%) and inter-site repeatability (CoV: 2.90% to 5.76%), and only a minor deviation of -1.65 mm² compared to MPRAGE. Crucially, by assessing atrophy rates and by comparing SCI patients with healthy controls longitudinally, differences between synT₁-w and MPRAGE were negligible. These results demonstrate that reliable estimates of CSA can be obtained from synT₁-w images, thereby reducing scan time significantly.

Quantitative muscle MRI displays clinically relevant myostructural abnormalities in long-term ICU-survivors: a case-control study

BACKGROUND

Long-term data on ICU-survivors reveal persisting sequalae and a reduced quality-of-life even after years. Major complaints are neuromuscular dysfunction due to Intensive care unit acquired weakness (ICUAW). Quantitative MRI (qMRI) protocols can quantify muscle alterations in contrast to standard qualitative MRI-protocols.

METHODS

Using qMRI, the aim of this study was to analyse persisting myostructural abnormalities in former ICU patients compared to controls and relate them to clinical assessments. The study was conducted as a cohort/case-control study. Nine former ICU-patients and matched controls were recruited (7 males; 54.8y ± 16.9; controls: 54.3y ± 11.1). MRI scans were performed on a 3T-MRI including a mDTI, T2 mapping and a mDixonquant sequence. Water T2 times, fat-fraction and mean values of the eigenvalue (λ₁), mean diffusivity (MD), radial diffusivity (RD) and fractional anisotropy (FA) were obtained for six thigh and seven calf muscles bilaterally. Clinical assessment included strength testing, electrophysiologic studies and a questionnaire on quality-of-life (QoL). Study groups were compared using a multivariate general linear model. qMRI parameters were correlated to clinical assessments and QoL questionnaire using Pearson´s correlation.

RESULTS

qMRI parameters were significantly higher in the patients for fat-fraction (p < 0.001), water T2 time (p < 0.001), FA (p = 0.047), MD (p < 0.001) and RD (p < 0.001). Thighs and calves showed a different pattern with significantly higher water T2 times only in the calves. Correlation analysis showed a significant negative correlation of muscle strength (MRC sum score) with FA and T2-time. The results were related to impairment seen in QoL-questionnaires, clinical testing and electrophysiologic studies.

CONCLUSION

qMRI parameters show chronic next to active muscle degeneration in ICU survivors even years after ICU therapy with ongoing clinical relevance. Therefore, qMRI opens new doors to characterize and monitor muscle changes of patients with ICUAW. Further, better understanding on the underlying mechanisms of the persisting complaints could contribute the development of personalized rehabilitation programs.

Cardiac Magnetic Resonance Fingerprinting: Potential Clinical Applications

PURPOSE OF REVIEW

Cardiac magnetic resonance fingerprinting (cMRF) has developed as a technique for rapid, multi-parametric tissue property mapping that has potential to both improve cardiac MRI exam efficiency and expand the information captured. In this review, we describe the cMRF technique, summarize technical developments and in vivo reports, and highlight potential clinical applications.

RECENT FINDINGS

Technical developments in cMRF continue to progress rapidly, including motion compensated reconstruction, additional tissue property quantification, signal time course analysis, and synthetic LGE image generation. Such technical developments can enable simplified CMR protocols by combining multiple evaluations into a single protocol and reducing the number of breath-held scans. cMRF continues to be reported for use in a range of pathologies; however barriers to clinical implementation remain. Technical developments are described in this review, followed by a focus on potential clinical applications that they may support. Clinical translation of cMRF could shorten protocols, improve CMR accessibility, and provide additional information as compared to conventional cardiac parametric mapping methods. Current needs for clinical implementation are discussed, as well as how those needs may be met in order to bring cMRF from its current research setting to become a viable tool for patient care.

Joint k-TE Space Image Reconstruction and Data Fitting for T2 Mapping

Objectives

To develop a joint k-TE reconstruction algorithm to reconstruct the T2-weighted (T2W) images and T2 map simultaneously.

Materials and Methods

The joint k-TE reconstruction model was formulated as an optimization problem subject to a self-consistency condition of the exponential decay relationship between the T2W images and T2 map. The objective function included a data fidelity term enforcing the agreement between the solution and the measured k-space data, together with a spatial regularization term on image properties of the T2W images. The optimization problem was solved using Alternating-Direction Method of Multipliers (ADMM). We tested the joint k-TE method in phantom data and healthy volunteer scans with fully-sampled and under-sampled k-space lines. Image quality of the reconstructed T2W images and T2 map, and the accuracy of T2 measurements derived by the joint k- TE and the conventional signal fitting method were compared.

Results

The proposed method improved image quality with reduced noise and less artifacts on both T2W images and T2 map, and increased measurement consistency in T2 relaxation time measurements compared with the conventional method in all data sets.

Conclusions

The proposed reconstruction method outperformed the conventional magnitude image-based signal fitting method in image quality and stability of quantitative T2 measurements.

Dynamics of progressive degeneration of major spinal pathways following spinal cord injury: A longitudinal study

BACKGROUND

Following spinal cord injury (SCI), disease processes spread gradually along the spinal cord forming a spatial gradient with most pronounced changes located at the lesion site. However, the dynamics of this gradient in SCI patients is not established.

OBJECTIVE

This study tracks the spatiotemporal dynamics of remote anterograde and retrograde spinal tract degeneration in the upper cervical cord following SCI over two years utilizing quantitative MRI.

METHODS

Twenty-three acute SCI patients (11 paraplegics, 12 tetraplegics) and 21 healthy controls were scanned with a T1-weighted sequence for volumetry and a FLASH sequence for myelin-sensitive magnetization transfer saturation (MTsat) of the upper cervical cord. We estimated myelin content from MTsat maps within the corticospinal tracts (CST) and dorsal columns (DC) and measured spinal cord atrophy by means of left-right width (LRW) and anterior-posterior width (APW) on the T1-weighted images across cervical levels C1-C3. MTsat in the CST and LRW were considered proxies for retrograde degeneration, while MTsat in the DC and APW provided evidence for anterograde degeneration, respectively. Using regression models, we compared the temporal and spatial trajectories of these MRI readouts between tetraplegics, paraplegics, and controls over a 2-year period and assessed their associations with clinical improvement.

RESULTS

Linear rates and absolute differences in myelin-sensitive MTsat indicated retrograde and anterograde neurodegeneration in the CST and DC, respectively. Changes in MTsat within the CST and in LRW progressively developed over time forming a gradient towards lower cervical levels by 2 years after injury, especially in tetraplegics (change per cervical level in MTsat: -0.247 p.u./level, p = 0.034; in LRW: -0.323 mm/level, p = 0.024). MTsat within the DC was already decreased at cervical levels C1-C3 at baseline (1.5 months after injury) in both tetra- and paraplegics, while linear decreases in APW over time were similar across C1-C3, preserving the spatial gradient. The relative improvement in light touch score was associated with MTsat within the DC at baseline (r_s = 0.575, p = 0.014).

CONCLUSION

Rostral and remote to the injury, the CST and DC show ongoing structural changes, indicative of myelin reductions and atrophy within 2 years after SCI. While anterograde degeneration in the DC was already detectable uniformly at C1-C3 early following SCI, retrograde degeneration in the CST developed over time revealing specific spatial and temporal neurodegenerative gradients. Disentangling and quantifying such dynamic pathological processes may provide biomarkers for regenerative and remyelinating therapies along entire spinal pathways.

Self-labelled encoder-decoder (SLED) for multi-echo gradient echo-based myelin water imaging

PURPOSE

Reconstruction of high quality myelin water imaging (MWI) maps is challenging, particularly for data acquired using multi-echo gradient echo (mGRE) sequences. A non-linear least squares fitting (NLLS) approach has often been applied for MWI. However, this approach may produce maps with limited detail and, in some cases, sub-optimal signal to noise ratio (SNR), due to the nature of the voxel-wise fitting. In this study, we developed a novel, unsupervised learning method called self-labelled encoder-decoder (SLED) to improve gradient echo-based MWI data fitting.

METHODS

Ultra-high resolution, MWI data was collected from five mouse brains with variable levels of myelination, using a mGRE sequence. Imaging data was acquired using a 7T preclinical MRI system. A self-labelled, encoder-decoder network was implemented in TensorFlow for calculation of myelin water fraction (MWF) based on the mGRE signal decay. A simulated MWI phantom was also created to evaluate the performance of MWF estimation.

RESULTS

Compared to NLLS, SLED demonstrated improved MWF estimation, in terms of both stability and accuracy in phantom tests. In addition, SLED produced less noisy MWF maps from high resolution MR microscopy images of mouse brain tissue. It specifically resulted in lower noise amplification for all mouse genotypes that were imaged and yielded mean MWF values in white matter ROIs that were highly correlated with those derived from standard NLLS fitting. Lastly, SLED also exhibited higher tolerance to low SNR data.

CONCLUSION

Due to its unsupervised and self-labeling nature, SLED offers a unique alternative to analyze gradient echo-based MWI data, providing accurate and stable MWF estimations.

A unified model for reconstruction and R2* mapping of accelerated 7T data using the quantitative recurrent inference machine

Quantitative MRI (qMRI) acquired at the ultra-high field of 7 Tesla has been used in visualizing and analyzing subcortical structures. qMRI relies on the acquisition of multiple images with different scan settings, leading to extended scanning times. Data redundancy and prior information from the relaxometry model can be exploited by deep learning to accelerate the imaging process. We propose the quantitative Recurrent Inference Machine (qRIM), with a unified forward model for joint reconstruction and R₂^*-mapping from sparse data, embedded in a Recurrent Inference Machine (RIM), an iterative inverse problem-solving network. To study the dependency of the proposed extension of the unified forward model to network architecture, we implemented and compared a quantitative End-to-End Variational Network (qE2EVN). Experiments were performed with high-resolution multi-echo gradient echo data of the brain at 7T of a cohort study covering the entire adult life span. The error in reconstructed R₂^* from undersampled data relative to reference data significantly decreased for the unified model compared to sequential image reconstruction and parameter fitting using the RIM. With increasing acceleration factor, an increasing reduction in the reconstruction error was observed, pointing to a larger benefit for sparser data. Qualitatively, this was following an observed reduction of image blurriness in R₂^*-maps. In contrast, when using the U-Net as network architecture, a negative bias in R₂^* in selected regions of interest was observed. Compressed Sensing rendered accurate, but less precise estimates of R₂^*. The qE2EVN showed slightly inferior reconstruction quality compared to the qRIM but better quality than the U-Net and Compressed Sensing. Subcortical maturation over age measured by a linearly increasing interquartile range of R₂^* in the striatum was preserved up to an acceleration factor of 9. With the integrated prior of the unified forward model, the proposed qRIM can exploit the redundancy among repeated measurements and shared information between tasks, facilitating relaxometry in accelerated MRI.

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.

CLIPPERS: Multiparametric and quantitative MRI features

Chronic lymphocytic inflammation with pontine perivascular enhancement responsive to steroids (CLIPPERS) is a rare chronic central-nervous-system inflammatory disorder that became known only recently, and the pathogenesis of CLIPPERS remains poorly understood. This report presents clinical and radiological features of a rare case: a young female patient who rapidly died of suspected CLIPPERS. Helpful multiparametric MRI diagnostic criteria are proposed that can help discriminate CLIPPERS from non-CLIPPERS pathologies. We reviewed clinical history, symptoms, quantitative data from brain multiparametric MRI before and after treatment, and histopathological data. Perfusion-weighted imaging revealed a decrease in regional cerebral blood flow by 31% and in cerebral blood volume by 64%, with a moderate increase in transit time and in time to peak by up to 23% in affected pontine and cerebral white matter. As estimated by diffusion tensor imaging, there was elevated density of tracts (n/mm²) and a decrease of fraction anisotropy (×10^-3 mm/s²) in the patient's pons as compared to a healthy control: density of tracts = 13.5 vs 12.4 and fraction anisotropy = 0.32 vs 0.45, respectively. Macromolecular proton fraction values proved to be reduced (15.8% and 14.5% in the control, respectively) in the patient's cerebral peduncles by 3% and in the pons by 4.1% and in a periventricular white matter lesion by 6.4% (11.3% in the normal-looking contralateral hemisphere). Based on our findings, we argue that quantitative MRI techniques may be a valuable source of biomarkers and reliable diagnostic criteria and can shed light on the pathogenesis and exact nosological position of this disorder.

Clinical Utility of Magnetic Resonance Imaging Biomarkers for Identifying Nonalcoholic Steatohepatitis Patients at High Risk of Progression: A Multicenter Pooled Data and Meta-Analysis

BACKGROUND & AIMS

Nonalcoholic fatty liver disease (NAFLD) is increasing in prevalence worldwide. NAFLD is associated with excess risk of all-cause mortality, and its progression to nonalcoholic steatohepatitis (NASH) and fibrosis accounts for a growing proportion of cirrhosis and hepatocellular cancer and thus is a leading cause of liver transplant worldwide. Noninvasive precise methods to identify patients with NASH and NASH with significant disease activity and fibrosis are crucial when the disease is still modifiable. The aim of this study was to examine the clinical utility of corrected T1 (cT1) vs magnetic resonance imaging (MRI) liver fat for identification of NASH participants with nonalcoholic fatty liver disease activity score ≥4 and fibrosis stage (F) ≥2 (high-risk NASH).

METHODS

Data from five clinical studies (n = 543) with participants suspected of NAFLD were pooled or used for individual participant data meta-analysis. The diagnostic accuracy of the MRI biomarkers to stratify NASH patients was determined using the area under the receiver operating characteristic curve (AUROC).

RESULTS

A stepwise increase in cT1 and MRI liver fat with increased NAFLD severity was shown, and cT1 was significantly higher in participants with high-risk NASH. The diagnostic accuracy (AUROC) of cT1 to identify patients with NASH was 0.78 (95% CI, 0.74-0.82), for liver fat was 0.78 (95% CI, 0.73-0.82), and when combined with MRI liver fat was 0.82 (95% CI, 0.78-0.85). The diagnostic accuracy of cT1 to identify patients with high-risk NASH was good (AUROC = 0.78; 95% CI, 0.74-0.82), was superior to MRI liver fat (AUROC = 0.69; 95% CI, 0.64-0.74), and was not substantially improved by combining it with MRI liver fat (AUROC = 0.79; 95% CI, 0.75-0.83). The meta-analysis showed similar performance to the pooled analysis for these biomarkers.

CONCLUSIONS

This study shows that quantitative MRI-derived biomarkers cT1 and liver fat are suitable for identifying patients with NASH, and cT1 is a better noninvasive technology than liver fat to identify NASH patients at greatest risk of disease progression. Therefore, MRI cT1 and liver fat have important clinical utility to help guide the appropriate use of interventions in NAFLD and NASH clinical care pathways.

Determination of oxygen relaxivity in oxygen nanobubbles at 3 and 7 Tesla

OBJECTIVE

Oxygen-loaded nanobubbles have shown potential for reducing tumour hypoxia and improving treatment outcomes, however, it remains difficult to noninvasively measure the changes in partial pressure of oxygen (PO₂) in vivo. The linear relationship between PO₂ and longitudinal relaxation rate (R₁) has been used to noninvasively infer PO₂ in vitreous and cerebrospinal fluid, and therefore, this experiment aimed to investigate whether R₁ is a suitable measurement to study oxygen delivery from such oxygen carriers.

METHODS

T₁ mapping was used to measure R₁ in phantoms containing nanobubbles with varied PO₂ to measure the relaxivity of oxygen (r_1Ox) in the phantoms at 7 and 3 T. These measurements were used to estimate the limit of detection (LOD) in two experimental settings: preclinical 7 T and clinical 3 T MRI.

RESULTS

The r_1Ox in the nanobubble solution was 0.00057 and 0.000235 s^-1/mmHg, corresponding to a LOD of 111 and 103 mmHg with 95% confidence at 7 and 3 T, respectively.

CONCLUSION

This suggests that T₁ mapping could provide a noninvasive method of measuring a > 100 mmHg oxygen delivery from therapeutic nanobubbles.

Inconsistencies in atlas-based volumetric measures of the human nucleus basalis of Meynert: A need for high-resolution alternatives

The nucleus basalis of Meynert (nbM) is the major source of cortical acetylcholine (ACh) and has been related to cognitive processes and to neurological disorders. However, spatially delineating the human nbM in MRI studies remains challenging. Due to the absence of a functional localiser for the human nbM, studies to date have localised it using nearby neuroanatomical landmarks or using probabilistic atlases. To understand the feasibility of MRI of the nbM we set our four goals; our first goal was to review current human nbM region-of-interest (ROI) selection protocols used in MRI studies, which we found have reported highly variable nbM volume estimates. Our next goal was to quantify and discuss the limitations of existing atlas-based volumetry of nbM. We found that the identified ROI volume depends heavily on the atlas used and on the probabilistic threshold set. In addition, we found large disparities even for data/studies using the same atlas and threshold. To test whether spatial resolution contributes to volume variability, as our third goal, we developed a novel nbM mask based on the normalized BigBrain dataset. We found that as long as the spatial resolution of the target data was 1.3 mm isotropic or above, our novel nbM mask offered realistic and stable volume estimates. Finally, as our last goal we tried to discern nbM using publicly available and novel high resolution structural MRI ex vivo MRI datasets. We find that, using an optimised 9.4T quantitative T₂^⁎ ex vivo dataset, the nbM can be visualised using MRI. We conclude caution is needed when applying the current methods of mapping nbM, especially for high resolution MRI data. Direct imaging of the nbM appears feasible and would eliminate the problems we identify, although further development is required to allow such imaging using standard (f)MRI scanning.

Single-shot multi-parametric mapping based on multiple overlapping-echo detachment (MOLED) imaging

Multi-parametric quantitative magnetic resonance imaging (mqMRI) allows the characterization of multiple tissue properties non-invasively and has shown great potential to enhance the sensitivity of MRI measurements. However, real-time mqMRI during dynamic physiological processes or general motions remains challenging. To overcome this bottleneck, we propose a novel mqMRI technique based on multiple overlapping-echo detachment (MOLED) imaging, termed MQMOLED, to enable mqMRI in a single shot. In the data acquisition of MQMOLED, multiple MR echo signals with different multi-parametric weightings and phase modulations are generated and acquired in the same k-space. The k-space data is Fourier transformed and fed into a well-trained neural network for the reconstruction of multi-parametric maps. We demonstrated the accuracy and repeatability of MQMOLED in simultaneous mapping apparent proton density (APD) and any two parameters among T₂, T₂*, and apparent diffusion coefficient (ADC) in 130-170 ms. The abundant information delivered by the multiple overlapping-echo signals in MQMOLED makes the technique potentially robust to system imperfections, such as inhomogeneity of static magnetic field or radiofrequency field. Benefitting from the single-shot feature, MQMOLED exhibits a strong motion tolerance to the continuous movements of subjects. For the first time, it captured the synchronous changes of ADC, T₂, and T₁-weighted APD in contrast-enhanced perfusion imaging on patients with brain tumors, providing additional information about vascular density to the hemodynamic parametric maps. We expect that MQMOLED would promote the development of mqMRI technology and greatly benefit the applications of mqMRI, including therapeutics and analysis of metabolic/functional processes.

Accuracy and repeatability of joint sparsity multi-component estimation in MR Fingerprinting

MR fingerprinting (MRF) is a promising method for quantitative characterization of tissues. Often, voxel-wise measurements are made, assuming a single tissue-type per voxel. Alternatively, the Sparsity Promoting Iterative Joint Non-negative least squares Multi-Component MRF method (SPIJN-MRF) facilitates tissue parameter estimation for identified components as well as partial volume segmentations. The aim of this paper was to evaluate the accuracy and repeatability of the SPIJN-MRF parameter estimations and partial volume segmentations. This was done (1) through numerical simulations based on the BrainWeb phantoms and (2) using in vivo acquired MRF data from 5 subjects that were scanned on the same week-day for 8 consecutive weeks. The partial volume segmentations of the SPIJN-MRF method were compared to those obtained by two conventional methods: SPM12 and FSL. SPIJN-MRF showed higher accuracy in simulations in comparison to FSL- and SPM12-based segmentations: Fuzzy Tanimoto Coefficients (FTC) comparing these segmentations and Brainweb references were higher than 0.95 for SPIJN-MRF in all the tissues and between 0.6 and 0.7 for SPM12 and FSL in white and gray matter and between 0.5 and 0.6 in CSF. For the in vivo MRF data, the estimated relaxation times were in line with literature and minimal variation was observed. Furthermore, the coefficient of variation (CoV) for estimated tissue volumes with SPIJN-MRF were 10.5% for the myelin water, 6.0% for the white matter, 5.6% for the gray matter, 4.6% for the CSF and 1.1% for the total brain volume. CoVs for CSF and total brain volume measured on the scanned data for SPIJN-MRF were in line with those obtained with SPM12 and FSL. The CoVs for white and gray matter volumes were distinctively higher for SPIJN-MRF than those measured with SPM12 and FSL. In conclusion, the use of SPIJN-MRF provides accurate and precise tissue relaxation parameter estimations taking into account intrinsic partial volume effects. It facilitates obtaining tissue fraction maps of prevalent tissues including myelin water which can be relevant for evaluating diseases affecting the white matter.

Clinical translation of quantitative magnetic resonance imaging biomarkers - An overview and gap analysis of current practice

PURPOSE

This overview of the current landscape of quantitative magnetic resonance imaging biomarkers (qMR IBs) aims to support the standardisation of academic IBs to assist their translation to clinical practice.

METHODS

We used three complementary approaches to investigate qMR IB use and quality management practices within the UK: 1) a literature search of qMR and quality management terms during 2011-2015 and 2016-2020; 2) a database search for clinical research studies using qMR IBs during 2016-2020; and 3) a survey to ascertain the current availability and quality management practices for clinical MRI scanners and associated equipment at research institutions across the UK.

RESULTS

The analysis showed increased use of all qMR methods between the periods 2011-2015 and 2016-2020 and diffusion-tensor MRI and volumetry to be popular methods. However, the "translation ratio" of journal articles to clinical research studies was higher for qMR methods that have evidence of clinical translation via a commercial route, such as fat fraction and T2 mapping. The number of journal articles citing quality management terms doubled between the periods 2011-2015 and 2016-2020; although, its proportion relative to all journal articles only increased by 3.0%. The survey suggested that quality assurance (QA) and quality control (QC) of data acquisition procedures are under-reported in the literature and that QA/QC of acquired data/data analysis are under-developed and lack consistency between institutions.

CONCLUSIONS

We summarise current attempts to standardise and translate qMR IBs, and conclude by outlining the ideal quality management practices and providing a gap analysis between current practice and a metrological standard.

A comparative assessment of myelin-sensitive measures in multiple sclerosis patients and healthy subjects

INTRODUCTION

Multiple Sclerosis (MS) is a common neurological disease primarily characterized by myelin damage in lesions and in normal - appearing white and gray matter (NAWM, NAGM). Several quantitative MRI (qMRI) methods are sensitive to myelin characteristics by measuring specific tissue biophysical properties. However, there are currently few studies assessing the relative reproducibility and sensitivity of qMRI measures to MS pathology in vivo in patients.

METHODS

We performed two studies. The first study assessed of the sensitivity of qMRI measures to MS pathology: in this work, we recruited 150 MS and 100 healthy subjects, who underwent brain MRI at 3 T including quantitative T1 mapping (qT1), quantitative susceptibility mapping (QSM), magnetization transfer saturation imaging (MTsat) and myelin water imaging for myelin water fraction (MWF). The sensitivity of qMRIs to MS focal pathology (MS lesions vs peri-plaque white/gray matter (PPWM/PPGM)) was studied lesion-wise; the sensitivity to diffuse normal appearing (NA) pathology was measured using voxel-wise threshold-free cluster enhancement (TFCE) in NAWM and vertex-wise inflated cortex analysis in NAGM. Furthermore, the sensitivity of qMRI to the identification of lesion tissue was investigated using a voxel-wise logistic regression analysis to distinguish MS lesion and PP voxels. The second study assessed the reproducibility of myelin-sensitive qMRI measures in a single scanner. To evaluate the intra-session and inter-session reproducibility of qMRI measures, we have investigated 10 healthy subjects, who underwent two brain 3 T MRIs within the same day (without repositioning), and one after 1-week interval. Five region of interest (ROIs) in white and deep grey matter areas were segmented, and inter- and intra- session reproducibility was studied using the intra-class correlation coefficient (ICC). Further, we also investigated the voxel-wise reproducibility of qMRI measures in NAWM and NAGM.

RESULTS

qT1 and QSM showed the highest sensitivity to distinguish MS focal WM and cortical pathology from peri-plaque WM (P < 0.0001), although QSM also showed the highest variance when applied to lesions. MWF and MTsat exhibited the highest sensitivity to NAWM pathology (P < 0.01). On the other hand, qT1 appeared to be the most sensitive measure to NAGM pathology (P < 0.01). All myelin-sensitive qMRI measures exhibited high inter/intra sessional ICCs in various WM and deep GM ROIs, in NAWM and in NAGM (ICC 0.82 ± 0.12).

CONCLUSION

This work shows that the applied qT1, MWF, MTsat and QSM are highly reproducible and exhibit differential sensitivity to focal and diffuse WM and GM pathology in MS patients.

Using quantitative MRI to study the association of isocitrate dehydrogenase (IDH) status with oxygen metabolism and cellular structure changes in glioma

OBJECTIVE

To investigate the characteristics of oxygen metabolism and the cellular structure of glioma using quantitative MRI to predict the isocitrate dehydrogenase 1 (IDH1) status and to further understand the biological characteristics of gliomas.

METHODS

In this retrospective study, 94 patients with gliomas eventually received quantitative MRI measures to study oxygen metabolism. The oxygen metabolism biomarker maps (oxygen extraction fraction [OEF] and cerebral metabolic rate of oxygen [CMRO₂]) and the tissue-cellular-specific (R₂t*) MRI relaxation parameter were evaluated in different regions of glioma.

RESULTS

MRI results showed differences in oxygen metabolism measures in all patients with gliomas of different IDH1 statuses. Compared to patients with IDH1 mutant gliomas, patients with IDH1 wild type gliomas showed increased (P < 0.01) CMRO₂, OEF, cerebral blood volume [CBF], and R₂t* measures in tumor regions, while only OEF, CBF and R₂t* were found to be increased (P < 0.05) in the peritumoral area. OEF achieved the best performance for distinguishing IDH1 wild type and mutant gliomas in the tumor area (AUC = 0.732, P < 0.001). R₂t* values correlated with Ki-67(R = 0.35, P < 0.001) in the tumor area, while no significant correlations between Ki-67 and R₂t* were found in the peritumoral area (R = 0.19, P = 0.072).

CONCLUSION

Quantitative MRI has potential applications in studying the tumor and peritumoral areas of glioma, and it has the ability to predict and reveal the characteristics of oxygen metabolism and cellular structure in different regions of gliomas.

Determining the internal orientation, degree of ordering, and volume of elongated nanocavities by NMR: Application to studies of plant stem

This study investigates the fibril nanostructure of fresh celery samples by modeling the anisotropic behavior of the transverse relaxation time (T₂) in nuclear magnetic resonance (NMR). Experimental results are interpreted within the framework of a previously developed theory, which was successfully used to model the nanostructures of several biological tissues as a set of water filled nanocavities, hence explaining the anisotropy the T₂ relaxation time in vivo. An important feature of this theory is to determine the degree of orientational ordering of the nanocavities, their characteristic volume, and their average direction with respect to the macroscopic sample. Results exhibit good agreement between theory and experimental data, which are, moreover, supported by optical microscopic resolution. The quantitative NMR approach presented herein can be potentially used to determine the internal ordering of biological tissues noninvasively.

Repeatability of diffusion-weighted magnetic resonance imaging in head and neck cancer at a 1.5 T MR-Linac

BACKGROUND AND PURPOSE

Functional information acquired through diffusion-weighted magnetic resonance imaging (DW-MRI) may be beneficial for personalized head and neck cancer (HNC) radiotherapy. Technical validation is required before DW-MRI based radiotherapy interventions can be realized clinically. The aim of this study was to assess the repeatability of apparent diffusion coefficients (ADC) derived from DW-MRI in HNC using echo-planar imaging (EPI) on a 1.5 T MR-Linac.

MATERIAL AND METHODS

A total of eleven HNC patients underwent test/retest DW-MRI scans at least once per week during fractionated radiotherapy at the MR-Linac. An EPI DW-MRI test scan (b=0, 150, 500 s/mm²) was acquired before the start of adaptive MR-guided radiotherapy in addition to an identical retest scan after irradiation. Volumes-of-interest (VOI) were defined manually for parotid (PTs) and submandibular glands (SMs), gross tumor volume (GTV) and lymph nodes (LNs). Mean ADC was calculated for all VOI in all test/retest scans. Absolute/relative repeatability coefficients (RCs/relRCs) as well as intraclass correlation coefficients (ICCs) were determined for all VOI.

RESULTS

A total of 81 datasets were analyzed. Mean test ADC values were 1380/1416, 950/1010, 1520 and 1344·10^-6 mm²/s for left/right SM and PT, GTV and LNs, respectively. Accordingly, RC (relRC) values were determined as 271/281 (19.4/21.8%) and 138/155 (13.3/15.2%), 457 (31.3%) and 310·10^-6 mm²/s (23.5%). ICC resulted in 0.80/0.87, 0.97/0.94, 0.75 and 0.83 for left/right SM and PT, GTV and LNs, respectively.

CONCLUSION

The repeatability of ADC derived from EPI DW-MRI at the 1.5 T MR-Linac appears reasonable to be used for future biologically adapted MR-guided radiotherapy.

Reproducibility and repeatability of magnetic resonance imaging in dementia

PURPOSE

Individualised predictive models of cognitive decline require disease-monitoring markers that are repeatable. For wide-spread adoption, such markers also need to be reproducible at different locations. This study assessed the repeatability and reproducibility of MRI markers derived from a dementia protocol.

METHODS

Six participants were scanned at three different sites with a 3T MRI scanner. The protocol employed: T1-weighted (T1w) imaging, resting state functional MRI (rsfMRI), arterial spin labelling (ASL), diffusion-weighted imaging (DWI), T2-weighted fluid attenuation inversion recovery (FLAIR), T2-weighted (T2w) imaging, and susceptibility weighted imaging (SWI). Participants were scanned repeatedly, up to six times over a maximum period of five years. One participant was also scanned a further three times on sequential days on one scanner. Fifteen derived metrics were computed from the seven different modalities.

RESULTS

Reproducibility (coefficient of variation; CoV, across sites) was best for T1w derived grey matter, white matter and hippocampal volume (CoV < 1.5%), compared to rsfMRI and SWI derived metrics (CoV, 19% and 21%). For a given metric, long-term repeatability (CoV across time) was comparable to reproducibility, with short-term repeatability considerably better.

CONCLUSIONS

Reproducibility and repeatability were assessed for a suite of markers calculated from a dementia MRI protocol. In general, structural markers were less variable than functional MRI markers. Variability over time on the same scanner was comparable to variability measured across different scanners. Overall, the results support the viability of multi-site longitudinal studies for monitoring cognitive decline.