STrategically Acquired Gradient Echo (STAGE) imaging, part II: Correcting for RF inhomogeneities in estimating T1 and proton density

Accelerated 3D metabolite T1 mapping of the brain using variable-flip-angle SPICE

PURPOSE

To develop a practical method to enable 3D T1 mapping of brain metabolites.

THEORY AND METHODS

Due to the high dimensionality of the imaging problem underlying metabolite T1 mapping, measurement of metabolite T1 values has been currently limited to a single voxel or slice. This work achieved 3D metabolite T1 mapping by leveraging a recent ultrafast MRSI technique called SPICE (spectroscopic imaging by exploiting spatiospectral correlation). The Ernst-angle FID MRSI data acquisition used in SPICE was extended to variable flip angles, with variable-density sparse sampling for efficient encoding of metabolite T1 information. In data processing, a novel generalized series model was used to remove water and subcutaneous lipid signals; a low-rank tensor model with prelearned subspaces was used to reconstruct the variable-flip-angle metabolite signals jointly from the noisy data.

RESULTS

The proposed method was evaluated using both phantom and healthy subject data. Phantom experimental results demonstrated that high-quality 3D metabolite T1 maps could be obtained and used for correction of T1 saturation effects. In vivo experimental results showed metabolite T1 maps with a large spatial coverage of 240 × 240 × 72 mm3 and good reproducibility coefficients (< 11%) in a 14.5-min scan. The metabolite T1 times obtained ranged from 0.99 to 1.44 s in gray matter and from 1.00 to 1.35 s in white matter.

CONCLUSION

We successfully demonstrated the feasibility of 3D metabolite T1 mapping within a clinically acceptable scan time. The proposed method may prove useful for both T1 mapping of brain metabolites and correcting the T1-weighting effects in quantitative metabolic imaging.

In Vivo visualization of white matter fiber tracts in the brainstem using low flip angle double echo 3D gradient echo imaging at 3T

BACKGROUND

White matter (WM) fiber tracts in the brainstem communicate with various regions in the cerebrum, cerebellum, and spinal cord. Clinically, small lesions, malformations, or histopathological changes in the brainstem can cause severe neurological disorders. A direct and non-invasive assessment approach could bring valuable information about the intricate anatomical variations of the white matter fiber tracts and nuclei. Although tractography from diffusion tensor imaging has been commonly used to map the WM fiber tracts connectivity, it is difficult to differentiate the complex WM tracts anatomically. Both high field MRI methods and ultrahigh-field MRI methods at 7T and 11.7 T have been used to enhance the contrast of WM fiber tracts. Despite their promising results, it is still challenging to achieve wide clinical adoption at 3T. In this study, we explored a clinically feasible method using a proton density weighted (PDW) 3D gradient echo (GRE) sequence to directly image the WM fiber tracts in the brainstem at 3T in vivo.

METHODS

We optimized a 3D high resolution, double echo, short TR, PDW GRE sequence on 5 healthy volunteers using a clinical 3T scanner to visualize the complicated anatomy of WM fiber tracts in the brain stem. Tissue properties including T1, proton density and T2* from in vivo quantitative MRI data were used for simulations to determine the optimal flip angle for the sequence. The visualization of multiple WM fiber tracts in the brainstem was assessed qualitatively and quantitatively using relative contrast and contrast-to-noise ratio (CNR). To improve the CNR, the final images were created by averaging over all echoes from two consecutive scans at the optimal flip angle. The results were compared to anatomical atlases and histology sections to identify the major fiber tracts. All the identified major fiber tracts were labeled on axial, sagittal and coronal slices.

RESULTS

The WM fiber tracts were found to have distinct hypointense signal throughout the brainstem and most of the major WM fiber tracts, such as the corticospinal tract, medial lemniscus, medial longitudinal fasciculus, and central tegmental tract, in the brainstem up to and including the thalamus were identified in all subjects. Both qualitative and quantitative evaluations showed that the 3° scan offered the best contrast for WM fiber tracts for a TR of 20 ms. The average over the first two echo times and two consecutive 3° scans gave a CNR of 47.8 ± 6.2 for the pyramidal tracts in particular and CNRs values greater than 6.5 ± 2.4 for the rest of the fiber tracts.

CONCLUSIONS

All the major fiber tracts in the brainstem could be visualized. Given the reasonably short scan time of 10 min at 3T, double echo PDW GRE sequence is a very practical approach for clinical adoption.

The STRAT-PARK cohort: A personalized initiative to stratify Parkinson's disease

The STRAT-PARK initiative aims to provide a platform for stratifying Parkinson's disease (PD) into biological subtypes, using a bottom-up, multidisciplinary biomarker-based and data-driven approach. PD is a heterogeneous entity, exhibiting high interindividual clinicopathological variability. This diversity suggests that PD may encompass multiple distinct biological entities, each driven by different molecular mechanisms. Molecular stratification and identification of disease subtypes is therefore a key priority for understanding and treating PD. STRAT-PARK is a multi-center longitudinal cohort aiming to recruit a total of 2000 individuals with PD and neurologically healthy controls from Norway and Canada, for the purpose of identifying molecular disease subtypes. Clinical assessment is performed annually, whereas biosampling, imaging, and digital and neurophysiological phenotyping occur every second year. The unique feature of STRAT-PARK is the diversity of collected biological material, including muscle biopsies and platelets, tissues particularly useful for mitochondrial biomarker research. Recruitment rate is ∼150 participants per year. By March 2023, 252 participants were included, comprising 204 cases and 48 controls. STRAT-PARK is a powerful stratification initiative anticipated to become a global research resource, contributing to personalized care in PD.

PENCIL imaging: A novel approach for neuromelanin sensitive MRI in Parkinson's disease

BACKGROUND

Parkinson's disease (PD) is associated with the loss of neuromelanin (NM) and increased iron in the substantia nigra (SN). Magnetization transfer contrast (MTC) is widely used for NM visualization but has limitations in brain coverage and scan time. This study aimed to develop a new approach called Proton-density Enhanced Neuromelanin Contrast in Low flip angle gradient echo (PENCIL) imaging to visualize NM in the SN.

METHODS

This study included 30 PD subjects and 50 healthy controls (HCs) scanned at 3T. PENCIL and MTC images were acquired. NM volume in the SN pars compacta (SNpc), normalized image contrast (Cnorm), and contrast-to-noise ratio (CNR) were calculated. The change of NM volume in the SNpc with age was analyzed using the HC data. A group analysis compared differences between PD subjects and HCs. Receiver operating characteristic (ROC) analysis and area under the curve (AUC) calculations were used to evaluate the diagnostic performance of NM volume and CNR in the SNpc.

RESULTS

PENCIL provided similar visualization and structural information of NM compared to MTC. In HCs, PENCIL showed higher NM volume in the SNpc than MTC, but this difference was not observed in PD subjects. PENCIL had higher CNR, while MTC had higher Cnorm. Both methods revealed a similar pattern of NM volume in SNpc changes with age. There were no significant differences in AUCs between NM volume in SNpc measured by PENCIL and MTC. Both methods exhibited comparable diagnostic performance in this regard.

CONCLUSIONS

PENCIL imaging provided improved CNR compared to MTC and showed similar diagnostic performance for differentiating PD subjects from HCs. The major advantage is PENCIL has rapid whole-brain coverage and, when using STAGE imaging, offers a one-stop quantitative assessment of tissue properties.

Strategically Acquired Gradient Echo (STAGE) Imaging, part IV: Constrained Reconstruction of White Noise (CROWN) Processing as a Means to Improve Signal-to-Noise in STAGE Imaging at 3 Tesla

Increasing the signal-to-noise ratio (SNR) has always been of critical importance for magnetic resonance imaging. Although increasing field strength provides a linear increase in SNR, it is more and more costly as field strength increases. Therefore, there is a major effort today to use signal processing methods to improve SNR since it is more efficient and economical. There are a variety of methods to improve SNR such as averaging the data at the expense of imaging time, or collecting the data with a lower resolution, all of these methods, including imaging processing methods, usually come at the expense of loss of image detail or image blurring. Therefore, we developed a new mathematical approach called CROWN (Constrained Reconstruction of White Noise) to enhance SNR without loss of structural detail and without affecting scanning time. In this study, we introduced and tested the concept behind CROWN specifically for STAGE (strategically acquired gradient echo) imaging. The concept itself is presented first, followed by simulations to demonstrate its theoretical effectiveness. Then the SNR improvement on proton spin density (PSD) and R2⁎ maps was investigated using brain STAGE data acquired from 10 healthy controls (HCs) and 10 patients with Parkinson's disease (PD). For the PSD and R2* maps, the SNR and CNR between white matter and gray matter were improved by a factor of 1.87 ± 0.50 and 1.72 ± 0.88, respectively. The white matter hyperintensity lesions in PD patients were more clearly defined after CROWN processing. Using these improved maps, simulated images for any repeat time, echo time or flip angle can be created with improved SNR. The potential applications of this technology are to trade off the increased SNR for higher resolution images and/or faster imaging.

Imaging Biomarkers of Peripheral Nerves: Focus on Magnetic Resonance Neurography and Ultrasonography

Peripheral neuropathy is a prevalent and debilitating condition affecting millions of individuals globally. Magnetic resonance neurography (MRN) and ultrasonography (US) are noninvasive methods offering comprehensive visualization of peripheral nerves, using anatomical and functional imaging biomarkers to ensure accurate evaluation. For optimized MRN, superior and high-resolution two-dimensional and three-dimensional imaging protocols are essential. The anatomical MRN and US imaging markers include quantitative measures of nerve and fascicular size and signal, and qualitative markers of course and morphology. Among them, quantitative markers of T2-signal intensity ratio are sensitive to nerve edema-like signal changes, and the T1-mapping technique reveals nerve and muscle tissue fatty and fibrous compositional alterations.The functional markers are derived from physiologic properties of nerves, such as diffusion characteristics or blood flow. They include apparent diffusion coefficient from diffusion-weighted imaging and fractional anisotropy and tractography from diffusion tensor imaging to delve into peripheral nerve microstructure and integrity. Peripheral nerve perfusion using dynamic contrast-enhanced magnetic resonance imaging estimates perfusion parameters, offering insights into nerve health and neuropathies involving edema, inflammation, demyelination, and microvascular alterations in conditions like type 2 diabetes, linking nerve conduction pathophysiology to vascular permeability alterations.Imaging biomarkers thus play a pivotal role in the diagnosis, prognosis, and monitoring of nerve pathologies, thereby ensuring comprehensive assessment and elevating patient care. These biomarkers provide valuable insights into nerve structure, function, and pathophysiology, contributing to the accurate diagnosis and management planning for peripheral neuropathy.

Characterization of white matter lesions in multiple sclerosis using proton density and T1-relaxation measures

PURPOSE

Although lesion dissemination in time is a defining characteristic of multiple sclerosis (MS), there is a limited understanding of lesion heterogeneity. Currently, conventional sequences such as fluid attenuated inversion recovery (FLAIR) and T1-weighted (T1W) data are used to assess MS lesions qualitatively. Estimating water content could provide a measure of local tissue rarefaction, or reduced tissue density, resulting from chronic inflammation. Our goal was to utilize the proton spin density (PD), derived from a rapid, multi-contrast STAGE (strategically acquired gradient echo) protocol to characterize white matter (WM) lesions seen on T2W, FLAIR and T1W data.

MATERIALS AND METHODS

Twenty (20) subjects with relapsing-remitting MS were scanned at 3 T using T1W, T2-weighted, FLAIR and strategically acquired gradient echo (STAGE) sequences. PD and T1 maps were derived from the STAGE data. Disease severity scores, including Extended Disability Status Scale (EDSS) and Multiple Sclerosis Functional Composite (MSFC), were correlated with total, high PD and high T1 lesion volumes. A probability map of high PD regions and all lesions across all subjects was generated. Five perilesional normal appearing WM (NAWM) bands surrounding the lesions were generated to compare the median PD and T1 values in each band with the lesional values and the global WM.

RESULTS

T1W intensity was negatively correlated with PD as expected (R = -0.87, p < 0.01, R2 = 0.756) and the FLAIR signal was suppressed for high PD volumes within the lesions, roughly for PD ≥ 0.85. The threshold for high PD and T1 regions was set to 0.909 and 1953.6 ms, respectively. High PD regions showed a high probability of occurrence near the boundary of the lateral ventricles. EDSS score and nine-hole peg test (dominant and non-dominant hand) were significantly correlated with the total lesion volume and the volumes of high PD and T1 regions (p < 0.05). There was a significant difference in PD/T1 values between the high PD/T1 regions within the lesions and the remaining lesional tissue (p < 0.001). In addition, the PD values of the first NAWM perilesional band directly adjacent to the lesional boundary displayed a significant difference (p < 0.05) compared to the global WM.

CONCLUSION

Lesions with high PD and T1s had the highest probability of occurrence at the boundary of the lateral ventricles and likely represent chronic lesions with significant local tissue rarefaction. Moreover, the perilesional NAWM exhibited subtly increasing PD and T1 values from the NAWM up to the lesion boundary. Unlike on the T1 maps, the perilesional band adjacent to the lesion boundary possessed a significantly higher PD value than the global WM PD values. This shows that PD maps were sensitive to the subtle changes in NAWM surrounding the lesions.

Enhancing Nigrosome-1 Sign Identification via Interpretable AI using True Susceptibility Weighted Imaging

BACKGROUND

Nigrosome 1 (N1), the largest nigrosome region in the ventrolateral area of the substantia nigra pars compacta, is identifiable by the "N1 sign" in long echo time gradient echo MRI. The N1 sign's absence is a vital Parkinson's disease (PD) diagnostic marker. However, it is challenging to visualize and assess the N1 sign in clinical practice.

PURPOSE

To automatically detect the presence or absence of the N1 sign from true susceptibility weighted imaging by using deep-learning method.

STUDY TYPE

Prospective.

POPULATION/SUBJECTS

453 subjects, including 225 PD patients, 120 healthy controls (HCs), and 108 patients with other movement disorders, were prospectively recruited including 227 males and 226 females. They were divided into training, validation, and test cohorts of 289, 73, and 91 cases, respectively.

FIELD STRENGTH/SEQUENCE

3D gradient echo SWI sequence at 3T; 3D multiecho strategically acquired gradient echo imaging at 3T; NM-sensitive 3D gradient echo sequence with MTC pulse at 3T.

ASSESSMENT

A neuroradiologist with 5 years of experience manually delineated substantia nigra regions. Two raters with 2 and 36 years of experience assessed the N1 sign on true susceptibility weighted imaging (tSWI), QSM with high-pass filter, and magnitude data combined with MTC data. We proposed NINet, a neural model, for automatic N1 sign identification in tSWI images.

STATISTICAL TESTS

We compared the performance of NINet to the subjective reference standard using Receiver Operating Characteristic analyses, and a decision curve analysis assessed identification accuracy.

RESULTS

NINet achieved an area under the curve (AUC) of 0.87 (CI: 0.76-0.89) in N1 sign identification, surpassing other models and neuroradiologists. NINet localized the putative N1 sign within tSWI images with 67.3% accuracy.

DATA CONCLUSION

Our proposed NINet model's capability to determine the presence or absence of the N1 sign, along with its localization, holds promise for enhancing diagnostic accuracy when evaluating PD using MR images.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY: Stage 1.

Impaired T1 mapping and Tmax during the first 7 days after ischemic stroke. A retrospective observational study

OBJECTIVE

To measure the relative T1 (rT1) value in different hypo-perfused regions after ischemic stroke using T1 mapping derived by Strategically Acquired Gradient Echo (STAGE) and assess its relationship with onset time and severity of ischemia.

MATERIALS AND METHODS

Sixty-three patients with acute anterior circulation ischemic stroke from 2017 to 2022 who underwent STAGE, diffusion weighted imaging (DWI) and dynamic susceptibility contrast perfusion weighted imaging (DSC-PWI) within 7 days were retrospectively enrolled. The areas with reduced diffusion and hypo-perfusion were segmented based on apparent diffusion coefficient (ADC) value < 0.62 × 10-3mm2/s and time-to-maximum (Tmax) thresholds (4, 6, 8, and 10 seconds). We measured the T1 value in the diffusion reduced and every 2 s Tmax strata regions and calculated rT1 (T1ipsi/T1contra) to explore the relationship between rT1 value, Tmax, and onset time.

RESULTS

rT1 value was increased in diffusion reduced (1.42) and hypo-perfused regions (1.02, 1.06, 1.12, 1.27, Tmax 4-6 s, 6-8 s, 8-10 s, > 10 s, respectively; all different from 1, P < 0.001). rT1 value was positively correlated with Tmax (rs = 0.61, P < 0.001) and onset time in area with reduced diffusion (rs = 0.39, P = 0.014).

CONCLUSIONS

Increased rT1 value in different hypo-perfused brain regions using T1 mapping derived by STAGE may reflect the edema; it was associated with the severity of Tmax and showed a weak correlation with the onset time in diffusion reduced areas.

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.

Multi-parametric MRI for radiotherapy simulation

Magnetic resonance imaging (MRI) has become an important imaging modality in the field of radiotherapy (RT) in the past decade, especially with the development of various novel MRI and image-guidance techniques. In this review article, we will describe recent developments and discuss the applications of multi-parametric MRI (mpMRI) in RT simulation. In this review, mpMRI refers to a general and loose definition which includes various multi-contrast MRI techniques. Specifically, we will focus on the implementation, challenges, and future directions of mpMRI techniques for RT simulation.

Computed cancer magnetic susceptibility imaging (canχ): Computational inverse mappings of cancer MRI

PURPOSE

We report a new cancer imaging modality in the contrast of tissue intrinsic susceptibility property by computed inverse magnetic resonance imaging (CIMRI).

METHODS

In MRI physics, an MRI signal is formed from tissue magnetism source (primarily magnetic susceptibility χ) through a cascade of MRI-introduced transformations (e.g. dipole-convolved magnetization) involving MRI setting parameters (e.g. echo time). In two-step computational inverse mappings (from phase image to internal fieldmap to susceptibility source), we could remove the MRI transformations and imaging parameters, thereby obtaining χ-depicted cancer images (canχ) from MRI phase images. Canχ is computationally implemented from clinical cancer MRI phase image by CIMRI.

RESULTS

As a result of MRI effect removal through computational inverse mappings, the reconstructed χ map (canχ) could provide a new cancerous tissue depiction in contrast of tissue intrinsic magnetism property (i.e. diamagnetism vs paramagnetism) as in an off-scanner state (e.g. in absence of main field B₀).

CONCLUSION

Through retrospective clinical cancer MRI data analysis, we reported on the canχ method in technical details and demonstrated its feasibility of innovating cancer imaging in the contrast of tissue intrinsic paramagnetism/diamagnetism property (in a cancer tissue state free from MRI effect).

Synthesizing Complex-Valued Multicoil MRI Data from Magnitude-Only Images

Despite the proliferation of deep learning techniques for accelerated MRI acquisition and enhanced image reconstruction, the construction of large and diverse MRI datasets continues to pose a barrier to effective clinical translation of these technologies. One major challenge is in collecting the MRI raw data (required for image reconstruction) from clinical scanning, as only magnitude images are typically saved and used for clinical assessment and diagnosis. The image phase and multi-channel RF coil information are not retained when magnitude-only images are saved in clinical imaging archives. Additionally, preprocessing used for data in clinical imaging can lead to biased results. While several groups have begun concerted efforts to collect large amounts of MRI raw data, current databases are limited in the diversity of anatomy, pathology, annotations, and acquisition types they contain. To address this, we present a method for synthesizing realistic MR data from magnitude-only data, allowing for the use of diverse data from clinical imaging archives in advanced MRI reconstruction development. Our method uses a conditional GAN-based framework to generate synthetic phase images from input magnitude images. We then applied ESPIRiT to derive RF coil sensitivity maps from fully sampled real data to generate multi-coil data. The synthetic data generation method was evaluated by comparing image reconstruction results from training Variational Networks either with real data or synthetic data. We demonstrate that the Variational Network trained on synthetic MRI data from our method, consisting of GAN-derived synthetic phase and multi-coil information, outperformed Variational Networks trained on data with synthetic phase generated using current state-of-the-art methods. Additionally, we demonstrate that the Variational Networks trained with synthetic k-space data from our method perform comparably to image reconstruction networks trained on undersampled real k-space data.

Visualizing the deep cerebellar nuclei using quantitative susceptibility mapping: An application in healthy controls, Parkinson's disease patients and essential tremor patients

The visualization and identification of the deep cerebellar nuclei (DCN) (dentate [DN], interposed [IN] and fastigial nuclei [FN]) are particularly challenging. We aimed to visualize the DCN using quantitative susceptibility mapping (QSM), predict the contrast differences between QSM and T2* weighted imaging, and compare the DCN volume and susceptibility in movement disorder populations and healthy controls (HCs). Seventy-one Parkinson's disease (PD) patients, 39 essential tremor patients, and 80 HCs were enrolled. The PD patients were subdivided into tremor dominant (TD) and postural instability/gait difficulty (PIGD) groups. A 3D strategically acquired gradient echo MR imaging protocol was used for each subject to obtain the QSM data. Regions of interest were drawn manually on the QSM data to calculate the volume and susceptibility. Correlation analysis between the susceptibility and either age or volume was performed and the intergroup differences of the volume and magnetic susceptibility in all the DCN structures were evaluated. For the most part, all the DCN structures were clearly visualized on the QSM data. The susceptibility increased as a function of volume for both the HC group and disease groups in the DN and IN (p < .001) but not the FN (p = .74). Only the volume of the FN in the TD-PD group was higher than that in the HCs (p = .012), otherwise, the volume and susceptibility among these four groups did not differ significantly. In conclusion, QSM provides clear visualization of the DCN structures. The results for the volume and susceptibility of the DCN can be used as baseline references in future studies of movement disorders.

Quantifying calcium changes in the fetal spine using quantitative susceptibility mapping as extracted from STAGE imaging

OBJECTIVES

To evaluate calcium deposition in the fetal spine in vivo during the second and third trimesters using quantitative susceptibility mapping (QSM).

METHODS

Fifty-four pregnant women in their second and third trimesters underwent a 2D multi-echo STrategically Acquired Gradient Echo (STAGE) MR imaging protocol at 3T covering the fetal spine. The first echo data was used for QSM processing. A linear regression model was used to assess the correlation between magnetic susceptibility and gestational age (GA). A paired sample t-test was used to compare the consistency of QSM measurements from each sequence.

RESULTS

The magnetic susceptibility of the fetal spine decreased linearly with advancing GA, with a slope of -52.3 parts per billion (ppb)/week and a Pearson correlation coefficient (r) of 0.83 (p < 0.001). In 37 subjects for whom the STAGE local QSM data were available from both flip angles, the average magnetic susceptibility values were -1111 ± 278 ppb and -1081 ± 262 ppb for FA = 8° and FA = 40°, respectively. These means were not statistically different according to a paired sample t-test (p = 0.156).

CONCLUSIONS

QSM is a reliable technique for evaluating calcium deposition and bone mineral density of fetal vertebrae. Our results demonstrate an increase in fetal calcium levels as a function of GA. These measures might be able to provide reference values for calcium content in the fetal spine during the second and third trimesters.

KEY POINTS

• Calcium deposition and mineralization in the fetal spine, evaluated by vertebral magnetic susceptibility, increased with advancing gestational age. • Our results provide reference values for calcium content in the fetal spine during the second and third trimesters.

Augmented T1 -weighted steady state magnetic resonance imaging

T₁ contrasts obtained using short-TR incoherent steady state gradient echo (GRE) methods are generally suboptimal, to which non-T₁ factors in the signals play a major part. In this work, we proposed an augmented T₁ -weighted (aT₁ W) method to extract the signal ratio between routine GRE T₁ W and proton density-weighted signals that effectively removes the non-T₁ effects from the original T₁ W signals, including proton density, T₂ * decay, and coil sensitivity. A recently proposed multidimensional integration (MDI) technique was incorporated in the aT₁ W calculation for better signal-to-noise ratio (SNR) performance. For comparison between aT₁ W and T₁ W results, Monte Carlo noise analysis was performed via simulation and on scanned data, and region-of-interest (ROI) analysis and comparison was performed on the system phantom. For brain scans, the image contrast, noise behavior, and SNR of aT₁ W images were compared with routine GRE and inversion-recovery-based T₁ W images. The proposed aT₁ W method yielded saliently improved T₁ contrasts (potentially > 30% higher contrast-to-noise ratio [CNR]) than routine GRE T₁ W images. Good spatial homogeneity and signal consistency as well as high SNR/CNR were achieved in aT₁ W images using the MDI technique. For contrast-enhanced (CE) imaging, aT₁ W offered stronger post-CE contrast and better boundary delineation than T₁ MPRAGE images while using a shorter scan time.

Quantitative susceptibility mapping as a measure of cerebral oxygenation in neonatal piglets

Prominence of cerebral veins using susceptibility weighted magnetic resonance imaging (SWI) has been used as a qualitative indicator of cerebral venous oxygenation (CvO2). Quantitative susceptibility mapping (QSM) adds more precision to the assessment of CvO2, but has not been applied to neonatal hypoxic ischemic injury (HII). We proposed to study QSM measures of venous susceptibility and their correlation with direct measures of brain oxygenation and cerebral blood flow (CBF) in the neonatal piglet. The association of QSM intravascular cerebral venous susceptibility, with brain tissue O2 tension, CBF, cortical tissue oxyhemoglobin saturation, and the partial pressure of oxygen in arterial blood measurement during various oxygenation states was determined by linear regression. Compared to normoxia, venous susceptibility in the straight sinus increased 56.8 ± 25.4% during hypoxia, while decreasing during hyperoxia (23.5 ± 32.9%) and hypercapnia (23.3 ± 73.1%), which was highly correlated to all other measures of oxygenation (p < 0.0001) but did not correlate to CBF (p = 0.82). These findings demonstrate a strong relationship between venous susceptibility and brain tissue O2 tension. Our results suggest that QSM-derived venous susceptibility is sensitive to cerebral oxygenation status across various oxygenation states.

Three-dimensional simultaneous brain mapping of T1, T2, T2∗ and magnetic susceptibility with MR Multitasking

PURPOSE

To develop a new technique that enables simultaneous quantification of whole-brain T1 , T2 , T 2 ∗ , as well as susceptibility and synthesis of six contrast-weighted images in a single 9.1-minute scan.

METHODS

The technique uses hybrid T2 -prepared inversion-recovery pulse modules and multi-echo gradient-echo readouts to collect k-space data with various T1, T2, and T 2 ∗ weightings. The underlying image is represented as a six-dimensional low-rank tensor consisting of three spatial dimensions and three temporal dimensions corresponding to T1 recovery, T2 decay, and multi-echo behaviors, respectively. Multiparametric maps were fitted from reconstructed image series. The proposed method was validated on phantoms and healthy volunteers, by comparing quantitative measurements against corresponding reference methods. The feasibility of generating six contrast-weighted images was also examined.

RESULTS

High quality, co-registered T1 , T2 , and T 2 ∗ susceptibility maps were generated that closely resembled the reference maps. Phantom measurements showed substantial consistency (R2 > 0.98) with the reference measurements. Despite the significant differences of T1 (p < .001), T2 (p = .002), and T 2 ∗ (p = 0.008) between our method and the references for in vivo studies, excellent agreement was achieved with all intraclass correlation coefficients greater than 0.75. No significant difference was found for susceptibility (p = .900). The framework is also capable of synthesizing six contrast-weighted images.

CONCLUSION

The MR Multitasking-based 3D brain mapping of T1 , T2 , T 2 ∗ , and susceptibility agrees well with the reference and is a promising technique for multicontrast and quantitative imaging.

MULTI-parametric MR imaging with fLEXible design (MULTIPLEX)

PURPOSE

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

THEORY

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

METHODS

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

RESULTS

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

CONCLUSION

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

A Comparison of Magnetic Resonance Imaging Methods to Assess Multiple Sclerosis Lesions: Implications for Patient Characterization and Clinical Trial Design

Magnetic resonance imaging (MRI) is a sensitive imaging modality for identifying inflammatory and/or demyelinating lesions, which is critical for a clinical diagnosis of MS and evaluating drug responses. There are many unique means of probing brain tissue status, including conventional T1 and T2 weighted imaging (T1WI, T2WI), T2 fluid attenuated inversion recovery (FLAIR), magnetization transfer, myelin water fraction, diffusion tensor imaging (DTI), phase-sensitive inversion recovery and susceptibility weighted imaging (SWI), but no study has combined all of these modalities into a single well-controlled investigation. The goals of this study were to: compare different MRI measures for lesion visualization and quantification; evaluate the repeatability of various imaging methods in healthy controls; compare quantitative susceptibility mapping (QSM) with myelin water fraction; measure short-term longitudinal changes in the white matter of MS patients and map out the tissue properties of the white matter hyperintensities using STAGE (strategically acquired gradient echo imaging). Additionally, the outcomes of this study were anticipated to aid in the choice of an efficient imaging protocol reducing redundancy of information and alleviating patient burden. Of all the sequences used, T2 FLAIR and T2WI showed the most lesions. To differentiate the putative demyelinating lesions from inflammatory lesions, the fusion of SWI and T2 FLAIR was used. Our study suggests that a practical and efficient imaging protocol combining T2 FLAIR, T1WI and STAGE (with SWI and QSM) can be used to rapidly image MS patients to both find lesions and study the demyelinating and inflammatory characteristics of the lesions.

All Central Nervous System Neuro- and Vascular-Communication Channels Are Surrounded With Cerebrospinal Fluid

Background: Recent emerging evidence has highlighted the potential critical role of cerebrospinal fluid (CSF) in cerebral waste clearance and immunomodulation. It is already very well-established that the central nervous system (CNS) is completely submerged in CSF on a macro-level; but to what extent is this true on a micro-level? Specifically, within the peri-neural and peri-vascular spaces within the CNS parenchyma. Therefore, the objective of this study was to use magnetic resonance imaging (MRI) to simultaneously map the presence of CSF within all peri-neural (cranial and spinal nerves) and peri-vascular spaces in vivo in humans. Four MRI protocols each with five participants were used to image the CSF in the brain and spinal cord. Our findings indicated that all CNS neuro- and vascular-communication channels are surrounded with CSF. In other words, all peri-neural spaces surrounding the cranial and spinal nerves as well as all peri-vascular spaces surrounding MRI-visible vasculature were filled with CSF. These findings suggest that anatomically, substance exchange between the brain parenchyma and outside tissues including lymphatic ones can only occur through CSF pathways and/or vascular pathways, warranting further investigation into its implications in cerebral waste clearance and immunity.

Automation of Quantifying Axonal Loss in Patients with Peripheral Neuropathies through Deep Learning Derived Muscle Fat Fraction

BACKGROUND

Axonal loss denervates muscle, leading to an increase of fat accumulation in the muscle. Therefore, fat fraction (FF) in whole limb muscle using MRI has emerged as a monitoring biomarker for axonal loss in patients with peripheral neuropathies. In this study, we are testing whether deep learning-based model can automate quantification of the FF in individual muscles. While individual muscle is smaller with irregular shape, manually segmented muscle MRI images have been accumulated in this lab; and make the deep learning feasible.

PURPOSE

To automate segmentation on muscle MRI images through deep learning for quantifying individual muscle FF in patients with peripheral neuropathies.

STUDY TYPE

Retrospective.

SUBJECTS

24 patients and 19 healthy controls.

FIELD STRENGTH/SEQUENCES

3T; Interleaved 3D GRE.

ASSESSMENT

A 3D U-Net model was implemented in segmenting muscle MRI images. This was enabled by leveraging a large set of manually segmented muscle MRI images. B1+ and B1- maps were used to correct image inhomogeneity. Accuracy of the automation was evaluated using Pixel Accuracy (PA), Dice Coefficient (DC) in binary masks; and Bland-Altman and Pearson correlation by comparing FF values between manual and automated methods.

STATISTICAL TESTS

PA and DC were reported with their median value and standard deviation. Two methods were compared using the ± 95% confidence intervals (CI) of Bland-Altman analysis and the Pearson's coefficient (r2 ).

RESULTS

DC values were from 0.83 ± 0.17 to 0.98 ± 0.02 in thigh and from 0.63 ± 0.18 to 0.96 ± 0.02 in calf muscles. For FF values, the overall ± 95% CI and r2 were [0.49, -0.56] and 0.989 in thigh and [0.84, -0.71] and 0.971 in the calf.

DATA CONCLUSION

Automated results well agreed with the manual results in quantifying FF for individual muscles. This method mitigates the formidable time consumption and intense labor in manual segmentations; and enables the use of individual muscle FF as outcome measures in upcoming longitudinal studies.

LEVEL OF EVIDENCE

3 TECHNICAL EFFICACY STAGE: 1.

Fetal brain tissue characterization at 1.5 T using STrategically Acquired Gradient Echo (STAGE) imaging

OBJECTIVES

To estimate human fetal brain MRI tissue properties including apparent T₁ (T_1app) and apparent proton density (PD_app) by using a rapid multi-contrast acquisition protocol called STrategically Acquired Gradient Echo (STAGE) imaging.

METHODS

STAGE data were collected using two flip angles (15° and 60°, with a TR = 600 ms) for 30 pregnant women at 1.5 T (15 healthy controls: gestational age (GA) range 19 + 1/7 weeks to 34 + 5/7 weeks; 11 abnormal subjects with ventriculomegaly: GA range 21 + 5/7 weeks to 31 + 5/7 weeks; 4 subjects with other abnormalities). Both T_1app and PD_app maps of the fetal brain were calculated from the STAGE data. A region-of-interest-based approach was used to measure T_1app and PD_app in the subplate/intermediate zone (SP/IZ), cortical plate (CP), and cerebrospinal fluid (CSF) in the fetal brain.

RESULTS

The ratios of T_1appSP/IZ/T_1appCP and PD_appSP/IZ/PD_appCP were larger than unity while T_1appSP/IZ/T_1appCSF and PD_appSP/IZ/PD_appCSF were both less than unity.

CONCLUSIONS

STAGE imaging provides a potential practical approach to estimate multi-parametric properties of the human fetal brain.

KEY POINTS

• STAGE is feasible in measuring fetal brain tissue properties. • Water content in cortical plate and subplate/intermediate zone approaches that of cerebrospinal fluid in early gestational ages.

Imaging iron and neuromelanin simultaneously using a single 3D gradient echo magnetization transfer sequence: Combining neuromelanin, iron and the nigrosome-1 sign as complementary imaging biomarkers in early stage Parkinson's disease

Diagnosing early stage Parkinson's disease (PD) is still a clinical challenge. Previous studies using iron, neuromelanin (NM) or the Nigrosome-1 (N1) sign in the substantia nigra (SN) by themselves have been unable to provide sufficiently high diagnostic performance for these methods to be adopted clinically. Our goal in this study was to extract the NM complex volume, iron content and volume representing the entire SN, and the N1 sign as potential complementary imaging biomarkers using a single 3D magnetization transfer contrast (MTC) gradient echo sequence and to evaluate their diagnostic performance and clinical correlations in early stage PD. A total of 40 early stage idiopathic PD subjects and 40 age- and sex-matched healthy controls (HCs) were imaged at 3T. NM boundaries (representing the SN pars compacta (SNpc) and parabrachial pigmented nucleus) and iron boundaries representing the total SN (SNpc and SN pars reticulata) were determined semi-automatically using a dynamic programming (DP) boundary detection algorithm. Receiver operating characteristic analyses were performed to evaluate the utility of these imaging biomarkers in diagnosing early stage PD. A correlation analysis was used to study the relationship between these imaging measures and the clinical scales. We also introduced the concept of NM and total iron overlap volumes to demonstrate the loss of NM relative to the iron containing SN. Furthermore, all 80 cases were evaluated for the N1 sign independently. The NM and SN volumes were lower while the iron content was higher in the SN for PD subjects compared to HCs. Interestingly, the PD subjects with bilateral loss of the N1 sign had the highest iron content. The area under the curve (AUC) values for the average of both hemispheres for single measures were: .960 for NM complex volume; .788 for total SN volume; .740 for SN iron content and .891 for the N1 sign. Combining NM complex volume with each of the following measures through binary logistic regression led to AUC values for the averaged right and left sides of: .976 for total iron content; .969 for total SN volume, .965 for overlap volume and .983 for the N1 sign. We found a negative correlation between SN volume and UPDRS-III (R² = .22, p = .002). While the N1 sign performed well, it does not contain any information about iron content or NM quantitatively, therefore, marrying this sign with the NM and iron measures provides a better physiological explanation of what is happening when the N1 sign disappears in PD subjects. In summary, the combination of NM complex volume, SN volume, iron content and the N1 sign as derived from a single MTC sequence provides complementary information for understanding and diagnosing early stage PD.

Towards in vivo g-ratio mapping using MRI: Unifying myelin and diffusion imaging

BACKGROUND

The g-ratio, quantifying the comparative thickness of the myelin sheath encasing an axon, is a geometrical invariant that has high functional relevance because of its importance in determining neuronal conduction velocity. Advances in MRI data acquisition and signal modelling have put in vivo mapping of the g-ratio, across the entire white matter, within our reach. This capacity would greatly increase our knowledge of the nervous system: how it functions, and how it is impacted by disease.

NEW METHOD

This is the second review on the topic of g-ratio mapping using MRI.

RESULTS

This review summarizes the most recent developments in the field, while also providing methodological background pertinent to aggregate g-ratio weighted mapping, and discussing pitfalls associated with these approaches.

COMPARISON WITH EXISTING METHODS

Using simulations based on recently published data, this review reveals caveats to the state-of-the-art calibration methods that have been used for in vivo g-ratio mapping. It highlights the need to estimate both the slope and offset of the relationship between these MRI-based markers and the true myelin volume fraction if we are really to achieve the goal of precise, high sensitivity g-ratio mapping in vivo. Other challenges discussed in this review further evidence the need for gold standard measurements of human brain tissue from ex vivo histology.

CONCLUSIONS

We conclude that the quest to find the most appropriate MRI biomarkers to enable in vivo g-ratio mapping is ongoing, with the full potential of many novel techniques yet to be investigated.

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

OBJECTIVES

The strategically acquired gradient echo (STAGE) protocol, developed for 3T scanners, allows one to derive quantitative maps such as T1, T2*, proton density, and quantitative susceptibility mapping in about 5 min. Our aim was to adapt the STAGE sequences for 1.5T scanners which are still commonly used in clinical practice. Furthermore, the accuracy and repeatability of the STAGE-derived T1 estimate were tested.

METHODS

Flip angle (FA) optimization was performed using a theoretical simulation by maximizing signal-to-noise ratio, contrast-to-noise ratio, and T1 precision. The FA choice was further refined with the ISMRM/NIST phantom and in vivo acquisitions. The accuracy of the T1 estimate was assessed by comparing STAGE-derived T1 values with T1 maps obtained with an inversion recovery sequence. T1 accuracy was investigated for both the phantom and in vivo data. Finally, one subject was acquired 10 times once a week and a group of 27 subjects was scanned once. The T1 coefficient of variation (COV) was computed to assess scan-rescan and physiological variability, respectively.

RESULTS

The FA_1,2 = 7°,38° were identified as the optimal FA pair at 1.5T. The T1 estimate errors were below 3% and 5% for phantom and in vivo measurements, respectively. COV for different tissues ranged from 1.8 to 4.8% for physiological variability, and between 0.8 and 2% for scan-rescan repeatability.

CONCLUSION

The optimized STAGE protocol can provide accurate and repeatable T1 mapping along with other qualitative images and quantitative maps in about 7 min on 1.5T scanners. This study provides the groundwork to assess the role of STAGE in clinical settings.

KEY POINTS

• The STAGE imaging protocol was optimized for use on 1.5T field strength scanners. • A practical STAGE protocol makes it possible to derive quantitative maps (i.e., T1, T2*, PD, and QSM) in about 7 min at 1.5T. • The T1 estimate derived from the STAGE protocol showed good accuracy and repeatability.

Multi-Echo Quantitative Susceptibility Mapping for Strategically Acquired Gradient Echo (STAGE) Imaging

Purpose

To develop a method to reconstruct quantitative susceptibility mapping (QSM) from multi-echo, multi-flip angle data collected using strategically acquired gradient echo (STAGE) imaging.

Methods

The proposed QSM reconstruction algorithm, referred to as “structurally constrained Susceptibility Weighted Imaging and Mapping” scSWIM, performs an ℓ₁ and ℓ₂ regularization-based reconstruction in a single step. The unique contrast of the T1 weighted enhanced (T1WE) image derived from STAGE imaging was used to extract reliable geometry constraints to protect the basal ganglia from over-smoothing. The multi-echo multi-flip angle data were used for improving the contrast-to-noise ratio in QSM through a weighted averaging scheme. The measured susceptibility values from scSWIM for both simulated and in vivo data were compared to the: original susceptibility model (for simulated data only), the multi orientation COSMOS (for in vivo data only), truncated k-space division (TKD), iterative susceptibility weighted imaging and mapping (iSWIM), and morphology enabled dipole inversion (MEDI) algorithms. Goodness of fit was quantified by measuring the root mean squared error (RMSE) and structural similarity index (SSIM). Additionally, scSWIM was assessed in ten healthy subjects.

Results

The unique contrast and tissue boundaries from T1WE and iSWIM enable the accurate definition of edges of high susceptibility regions. For the simulated brain model without the addition of microbleeds and calcium, the RMSE was best at 5.21ppb for scSWIM and 8.74ppb for MEDI thanks to the reduced streaking artifacts. However, by adding the microbleeds and calcium, MEDI’s performance dropped to 47.53ppb while scSWIM performance remained the same. The SSIM was highest for scSWIM (0.90) and then MEDI (0.80). The deviation from the expected susceptibility in deep gray matter structures for simulated data relative to the model (and for the in vivo data relative to COSMOS) as measured by the slope was lowest for scSWIM + 1%(−1%); MEDI + 2%(−11%) and then iSWIM −5%(−10%). Finally, scSWIM measurements in the basal ganglia of healthy subjects were in agreement with literature.

Conclusion

This study shows that using a data fidelity term and structural constraints results in reduced noise and streaking artifacts while preserving structural details. Furthermore, the use of STAGE imaging with multi-echo and multi-flip data helps to improve the signal-to-noise ratio in QSM data and yields less artifacts.

Synthetic MRI: Technologies and Applications in Neuroradiology

Synthetic MRI is a technique that synthesizes contrast-weighted images from multicontrast MRI data. There have been advances in synthetic MRI since the technique was introduced. Although a number of synthetic MRI methods have been developed for quantifying one or more relaxometric parameters and for generating multiple contrast-weighted images, this review focuses on several methods that quantify all three relaxometric parameters (T₁ , T₂ , and proton density) and produce multiple contrast-weighted images. Acquisition, quantification, and image synthesis techniques are discussed for each method. We discuss the image quality and diagnostic accuracy of synthetic MRI methods and their clinical applications in neuroradiology. Based on this analysis, we highlight areas that need to be addressed for synthetic MRI to be widely implemented in the clinic. LEVEL OF EVIDENCE: 5 TECHNICAL EFFICACY STAGE: 1.

Plaque characteristics of middle cerebral artery assessed using strategically acquired gradient echo (STAGE) and vessel wall MR contribute to misery downstream perfusion in patients with intracranial atherosclerosis

OBJECTIVES

To assess plaque vulnerability of the middle cerebral artery (MCA) using strategically acquired gradient echo (STAGE) versus high-resolution vessel wall MRI (hr-vwMRI), and explore the relationship between plaque characteristics and misery downstream perfusion.

METHODS

Ninety-one patients with single MCA atherosclerotic plaques underwent STAGE and hr-vwMRI were categorized into a group with misery perfusion and a group without based on the Alberta Stroke Program Early CT score (MTT-ASPECTS) with a threshold of 6. Plaque characteristics including inner lumen area (IWA), susceptibility, presence of hyperintensity within plaque (HIP), surface irregularity, stenosis degree, remodeling index, lipid ratio, and enhancement grade were compared between the two groups. The vulnerability of each plaque was retrospectively assessed on both STAGE and hr-vwMRI according to the combination of plaque features. Logistic regression analysis and ROC curve were performed to evaluate the effect of plaque characteristics on the presence of misery perfusion.

RESULTS

Taking hr-vwMRI as the reference, STAGE showed good efficiency in detecting vulnerable plaques. Patients with misery perfusion had less IWA, higher stenosis degree, more irregular surface and HIP, higher enhancement grade, and susceptibility (p < 0.01 for all). Higher susceptibility and stenosis degree were independent predictors for the occurrence of misery perfusion (p = 0.025, p = 0.048). The AUC was 0.900 for the combination of the two variables.

CONCLUSION

STAGE shows good efficiency to assess MCA plaque vulnerability versus hr-vwMRI. Plaque susceptibility evaluated using STAGE provides incremental value to predict misery perfusion combined with hr-vwMRI plaque features.

KEY POINTS

• STAGE has good efficiency in evaluating MCA plaque vulnerability versus hr-vwMRI. • Higher plaque susceptibility assessed using STAGE and higher grade luminal stenosis based on hr-vwMRI attribute to misery downstream perfusion. • STAGE provides incremental value on the understanding of plaque vulnerability in addition to conventional hr-vwMRI.

Rapid multicontrast brain imaging on a 0.35T MR-linac

PURPOSE

Magnetic resonance-guided radiation therapy (MRgRT) has shown great promise for localization and real-time tumor monitoring. However, to date, quantitative imaging has been limited for low field MRgRT. This work benchmarks quantitative T1, R2*, and Proton Density (PD)mapping in a phantom on a 0.35 T MR-linac and implements a novel acquisition method, STrategically Acquired Gradient Echo (STAGE). To further validate STAGE in a clinical setting, a pilot study was undertaken in a cohort of brain tumor patients to elucidate opportunities for longitudinal functional imaging with an MR-linac in the brain.

METHODS

STAGE (two triple-echo gradient echo (GRE) acquisitions) was optimized for a 0.35T low-field MR-linac. Simulations were performed to choose two flip angles to optimize signal-to-noise ratio (SNR) and T1-mapping precision. Tradeoffs between SNR, scan time, and spatial resolution for whole-brain coverage were evaluated in healthy volunteers. Data were inputted into a STAGE processing pipeline to yield four qualitative images (T1-weighted, enhanced T1-weighted, proton-density (PD) weighted, and simulated FLuid-Attenuated Inversion Recovery (sFLAIR)), and three quantitative datasets (T1, PD, and R2*). A benchmarking ISMRM/NIST phantom consisting of vials with variable NiCl₂ and MnCl₂ concentrations was scanned using variable flip angles (VFA) (2-60 degrees) and inversion recovery (IR) methods at 0.35 T. STAGE and VFA T1 values of vials were compared to IR T1 values. As measures of agreement with reference values and repeatability, relative error (RE) and coefficient of variability (CV) were calculated, respectively, for quantitative MR values within the phantom vials (spheres). To demonstrate feasibility, longitudinal STAGE data (pretreatment, weekly, and ~ 2 months post-treatment) were acquired in an IRB-approved pilot study of brain tumor cases via the generation of temporal and differential quantitative MRI maps.

RESULTS

In the phantom, RE of measured VFA T1 and STAGE relative to IR reference values were 7.0 ± 2.5% and 9.5 ± 2.2% respectively. RE for the PD vials was 8.1 ± 6.8% and CV for phantom R2* measurements was 10.1 ± 9.9%. Simulations and volunteer experiments yielded final STAGE parameters of FA = 50°/10°, 1 × 1 × 3 mm³ resolution, TR = 40 ms, TE = 5/20/34 ms in 10 min (64 slices). In the pilot study of brain tumor patients, differential maps for R2* and T1 maps were sensitive to local tumor changes and appeared similar to 3 T follow-up MRI datasets.

CONCLUSION

Quantitative T1, R2*, and PD mapping are promising at 0.35 T agreeing well with reference data. STAGE phantom data offer quantitative representations comparable to traditional methods in a fraction of the acquisition time. Initial feasibility of implementing STAGE at 0.35 T in a patient brain tumor cohort suggests that detectable changes can be observed over time. With confirmation in a larger cohort, results may be implemented to identify areas of recurrence and facilitate adaptive radiation therapy.

Peripheral nerve magnetic resonance imaging

Magnetic resonance imaging (MRI) has been used extensively in revealing pathological changes in the central nervous system. However, to date, MRI is very much underutilized in evaluating the peripheral nervous system (PNS). This underutilization is generally due to two perceived weaknesses in MRI: first, the need for very high resolution to image the small structures within the peripheral nerves to visualize morphological changes; second, the lack of normative data in MRI of the PNS and this makes reliable interpretation of the data difficult. This article reviews current state-of-the-art capabilities in in vivo MRI of human peripheral nerves. It aims to identify areas where progress has been made and those that still require further improvement. In particular, with many new therapies on the horizon, this review addresses how MRI can be used to provide non-invasive and objective biomarkers in the evaluation of peripheral neuropathies. Although a number of techniques are available in diagnosing and tracking pathologies in the PNS, those techniques typically target the distal peripheral nerves, and distal nerves may be completely degenerated during the patient’s first clinic visit. These techniques may also not be able to access the proximal nerves deeply embedded in the tissue. Peripheral nerve MRI would be an alternative to circumvent these problems. In order to address the pressing clinical needs, this review closes with a clinical protocol at 3T that will allow high-resolution, high-contrast, quantitative MRI of the proximal peripheral nerves.

Visualizing the lateral habenula using susceptibility weighted imaging and quantitative susceptibility mapping

The habenulae consist of a pair of small nuclei which bridge the limbic forebrain and midbrain monoaminergic centers. They are implicated in major depressive disorders due to abnormal phasic response when provoked by a conditioned stimulus. The lateral habenula (Lhb) is believed to be involved in dopamine metabolism and is now a target for deep brain stimulation, a treatment which has shown promising anti-depression effects. We imaged the habenulae with susceptibility weighted imaging (SWI) and quantitative susceptibility mapping (QSM) in order to localize the lateral habenula. Fifty-six healthy controls were recruited for this study. For the quantitative assessment, we traced the structure to compute volume from magnitude images and mean susceptibility bilaterally for the habenula on QSM. Thresholding methods were used to delineate the Lhb habenula on QSM. SWI, true SWI (tSWI), and QSM data were subjectively reviewed for increased Lhb contrast. SWI, QSM, and tSWI showed bilateral signal changes in the posterior location of the habenulae relative to the anterior location, which may indicate increased putative iron content within the Lhb. This signal behavior was shown in 41/44 (93%) subjects. In summary, it is possible to localize the lateral component of the habenula using SWI and QSM at 3 T.

STrategically Acquired Gradient Echo (STAGE) imaging, part III: Technical advances and clinical applications of a rapid multi-contrast multi-parametric brain imaging method

One major thrust in radiology today is image standardization with a focus on rapidly acquired quantitative multi-contrast information. This is critical for multi-center trials, for the collection of big data and for the use of artificial intelligence in evaluating the data. Strategically acquired gradient echo (STAGE) imaging is one such method that can provide 8 qualitative and 7 quantitative pieces of information in 5 min or less at 3 T. STAGE provides qualitative images in the form of proton density weighted images, T1 weighted images, T2* weighted images and simulated double inversion recovery (DIR) images. STAGE also provides quantitative data in the form of proton spin density, T1, T2* and susceptibility maps as well as segmentation of white matter, gray matter and cerebrospinal fluid. STAGE uses vendors' product gradient echo sequences. It can be applied from 0.35 T to 7 T across all manufacturers producing similar results in contrast and quantification of the data. In this paper, we discuss the strengths and weaknesses of STAGE, demonstrate its contrast-to-noise (CNR) behavior relative to a large clinical data set and introduce a few new image contrasts derived from STAGE, including DIR images and a new concept referred to as true susceptibility weighted imaging (tSWI) linked to fluid attenuated inversion recovery (FLAIR) or tSWI-FLAIR for the evaluation of multiple sclerosis lesions. The robustness of STAGE T1 mapping was tested using the NIST/NIH phantom, while the reproducibility was tested by scanning a given individual ten times in one session and the same subject scanned once a week over a 12-week period. Assessment of the CNR for the enhanced T1W image (T1WE) showed a significantly better contrast between gray matter and white matter than conventional T1W images in both patients with Parkinson's disease and healthy controls. We also present some clinical cases using STAGE imaging in patients with stroke, metastasis, multiple sclerosis and a fetus with ventriculomegaly. Overall, STAGE is a comprehensive protocol that provides the clinician with numerous qualitative and quantitative images.

Brain vascular changes in adults with congenital heart disease: A systematic review

Less information is available on brain integrity in adults with congenital heart disease than on brain changes in newborns and children with heart defects. Nevertheless, the number of adults with congenital heart disease is increasing rapidly and it has been shown that adults with congenital heart disease develop dementia almost twice as frequently as adults in the general population.

In the context of a rapidly growing congenital heart disease population, neuroradiological-oriented investigations of biomarkers distinctive for vascular damage, brain aging, and possible cognitive impairment is a crucial challenge. We provide an overview of the existing literature on neuroimaging studies in adults with congenital heart disease and discuss methodology issues to further investigate this subject. Overall, we aim to raise awareness of the importance of brain health studies in adults with congenital heart disease given the likely increasing impact on social and healthcare systems.

•

The number of adults with congenital heart disease is increasing rapidly.

•

There is considerable concern regarding global health outcomes in this population.

•

Evidence shows that patients with congenital heart disease are more prone to cerebrovascular accidents.

•

Neuroimaging tools may aid the detection of subtle vascular brain changes before cognitive impairment occurs.

Flexible proton density (PD) mapping using multi-contrast variable flip angle (VFA) data

Quantitative proton density (PD) maps measure the amount of free water, which is important for non-invasive tissue characterization in pathology and across lifespan. PD mapping requires the estimation and subsequent removal of factors influencing the signal intensity other than PD. These factors include the T1, T2* relaxation effects, transmit field inhomogeneities, receiver coil sensitivity profile (RP) and the spatially invariant factor that is required to scale the data. While the transmit field can be reliably measured, the RP estimation is usually based on image post-processing techniques due to limitations of its measurement at magnetic fields higher than 1.5 T. The post-processing methods are based on unified bias-field/tissue segmentation, fitting the sensitivity profile from images obtained with different coils, or on the linear relationship between T1 and PD. The scaling factor is derived from the signal within a specific tissue compartment or reference object. However, these approaches for calculating the RP and scaling factor have limitations particularly in severe pathology or over a wide age range, restricting their application. We propose a new approach for PD mapping based on a multi-contrast variable flip angle acquisition protocol and a data-driven estimation method for the RP correction and map scaling. By combining all the multi-contrast data acquired at different echo times, we are able to fully correct the MRI signal for T2* relaxation effects and to decrease the variance and the entropy of PD values within tissue class of the final map. The RP is determined from the corrected data applying a non-parametric bias estimation, and the scaling factor is based on the median intensity of an external calibration object. Finally, we compare the signal intensity and homogeneity of the multi-contrast PD map with the well-established effective PD (PD*) mapping, for which the RP is based on concurrent bias field estimation and tissue classification, and the scaling factor is estimated from the mean white matter signal. The multi-contrast PD values homogeneity and accuracy within the cerebrospinal fluid (CSF) and deep brain structures are increased beyond that obtained using PD* maps. We demonstrate that the multi-contrast RP approach is insensitive to anatomical or a priori tissue information by applying it in a patient with extensive brain abnormalities and for whole body PD mapping in post-mortem foetal imaging.

Quantifying iron content in magnetic resonance imaging

Measuring iron content has practical clinical indications in the study of diseases such as Parkinson's disease, Huntington's disease, ferritinopathies and multiple sclerosis as well as in the quantification of iron content in microbleeds and oxygen saturation in veins. In this work, we review the basic concepts behind imaging iron using T2, T2*, T2', phase and quantitative susceptibility mapping in the human brain, liver and heart, followed by the applications of in vivo iron quantification in neurodegenerative diseases, iron tagged cells and ultra-small superparamagnetic iron oxide (USPIO) nanoparticles.