Multidimensional correlation spectroscopic imaging of exponential decays: From theoretical principles to in vivo human applications

A data-driven intravoxel mean diffusivities distribution approach for molecular classifications and MIB-1 prediction of gliomas

BACKGROUND

Measuring non-parametric intravoxel mean diffusivity distributions (MDDs) using magnetic resonance imaging (MRI) is a sensitive method for detecting intracellular diffusivity changes during physiological alterations. Histological and molecular glioma classifications are essential for prognosis and treatment, with distinct water diffusion dynamics among subtypes.

PURPOSE

We developed a data-driven approach using a fully connected network (FCN) to enhance the speed and stability of calculating MDDs across varying SNRs, enable tumor microstructural mapping, and test its reliability in identifying MIB-1 labeling index (LI) levels and molecular status of gliomas.

METHODS

An FCN was trained to learn the mapping between the simulated diffusion decay curves and the ground truth MDDs. We performed 5 000 000 simulation curves with various diffusivity components and random SNR∈ [ 30 , 300 ] $ \in [ {30,\ 300} ]$ . Eighty percent of simulation curves were used for the FCN training, 10% for validation, and the others were external tests for the FCN performance evaluation. In vivo data were collected to evaluate its clinical reliability. One hundred one patients (44 years ± $ \pm $ 14, 67 men) with gliomas and six healthy controls underwent a 3.0 T MRI examination with a spin echo-echo planar imaging (SE-EPI) diffusion-weighted imaging (DWI) sequence. The trained FCN was employed to calculate MDDs of each brain voxel by voxel. We used the Fuzzy C-means algorithm to cluster the MDDs of tumor voxels, facilitating the characterization of distinct glioma tissues. Quantitative assessments were conducted through sectional integrals of the MDDs, demarcated by six bands to derive signal fractions (f n , n = 1 - 6 ${{f}_n},\ n = 1 -6$ ) and diffusivities of the maximum peaks (D p e a k ${{D}_{peak}}$ ). Cosine similarity scores (CSS) were used for MDD similarity. ANOVA and Mann-Whitney U test were used for difference analysis. Logistic regression and area under the receiver operator characteristic curve (AUC) were used for classification evaluation.

RESULTS

The simulation results showed that the FCN-based MDD approach (FCN-MDD) achieved higher CSS than non-negative least squares-based MDD (NNLS-MDD). For in vivo data, the spectra of ET and NET obtained by FCN-MDD are more distinguishable than NNLS-MDD. Fraction maps delineate the characteristics of different tumor tissues (enhancing and non-enhancing tumor, edema, and necrosis).f 3 , f 4 , D p e a k ${{f}_3},\ {{f}_4},{{D}_{peak}}$ showed a positive and negative correlation with MIB-1 respectively (r = 0.568 , r = - 0.521 , r = - 0.654 $r = 0.568,\ r = - 0.521,\ r = - 0.654$ , allp < 0.001 $p < 0.001$ ). The AUC ofD p e a k ${{D}_{peak}}$ for predicting MIB-1 LI levels was 0.900 (95% CI, 0.826-0.974), versus 0.781 (0.677-0.886) of ADC. The highest AUC of isocitrate dehydrogenase (IDH) mutation status, assessed by a logistic regression model (f 1 + f 3 ${{f}_1} + {{f}_3}$ ) was 0.873 (95% CI, 0.802-0.944).

CONCLUSION

The proposed FCN-MDD method was more robust to variations in SNR and less reliant on empirically set regularization values than the NNLS-MDD method. FCN-MDD also enabled qualitative and quantitative evaluation of the composition of gliomas.

Assessing renal interstitial fibrosis using compartmental, non-compartmental, and model-free diffusion MRI approaches

OBJECTIVE

To assess renal interstitial fibrosis (IF) using diffusion MRI approaches, and explore whether corticomedullary difference (CMD) of diffusion parameters, combination among MRI parameters, or combination with estimated glomerular filtration rate (eGFR) benefit IF evaluation.

METHODS

Forty-two patients with chronic kidney disease were included, undergoing MRI examinations. MRI parameters from apparent diffusion coefficient (ADC), intra-voxel incoherent motion (IVIM), diffusion kurtosis imaging (DKI), and diffusion-relaxation correlated spectrum imaging (DR-CSI) were obtained both for renal cortex and medulla. CMD of these parameters was calculated. Pathological IF scores (1-3) were obtained by biopsy. Patients were divided into mild (IF = 1, n = 23) and moderate-severe fibrosis (IF = 2-3, n = 19) groups. Group comparisons for MRI parameters were performed. Diagnostic performances were assessed by the receiver operator's curve analysis for discriminating mild from moderate-severe IF patients.

RESULTS

Significant inter-group differences existed for cortical ADC, IVIM-D, IVIM-f, DKI-MD, DR-CSI VB, and DR-CSI VC. Significant inter-group differences existed in ΔADC, ΔMD, ΔVB, ΔVC, ΔQB, and ΔQC. Among the cortical MRI parameters, VB displayed the highest AUC = 0.849, while ADC, f, and MD also showed AUC > 0.8. After combining cortical value and CMD, the diagnostic performances of the MRI parameters were slightly improved except for IVIM-D. Combining VB with f brings the best performance (AUC = 0.903) among MRI bi-variant models. A combination of cortical VB, ΔADC, and eGFR brought obvious improvement in diagnostic performance (AUC 0.963 vs 0.879, specificity 0.826 vs 0.896, and sensitivity 1.000 vs 0.842) than eGFR alone.

CONCLUSION

Our study shows promising results for the assessment of renal IF using diffusion MRI approaches.

CRITICAL RELEVANCE STATEMENT

Our study explores the non-invasive assessment of renal IF, an independent and effective predictor of renal outcomes, by comparing and combining diffusion MRI approaches including compartmental, non-compartmental, and model-free approaches.

KEY POINTS

Significant difference exists for diffusion parameters between mild and moderate-severe IF. Generally, cortical parameters show better performance than corresponding CMD. Bi-variant model lifts the diagnostic performance for assessing IF.

A genetic optimisation and iterative reconstruction framework for sparse multi-dimensional diffusion-relaxation correlation MRI

Multi-dimensional diffusion-relaxation correlation (DRC) magnetic resonance imaging (MRI) techniques have recently been developed to investigate tissue microstructures. Sub-voxel tissue heterogeneity is resolved from the local correlation distributions of relaxation time and molecular diffusivity. However, the implementation of these techniques considerably increases the total acquisition time, and simply reducing the scan time may be at the expense of detailed structural resolution. To overcome these limitations, an optimised framework was proposed for acquiring microstructural maps of the human brain on a clinically feasible timescale. First, the acquisition parameters of the multi-dimensional DRC MRI method were sparsely optimised using a genetic algorithm with a fitness function according to the spectral resolution of the correlation map, hardware requirements, and total scan time. Next, the acquired DRC MRI data were processed using a proposed numerical algorithm based on the dynamic inverse Laplace transform (ILT). Prior knowledge from one-dimensional data was then utilised in the iterative procedure to improve the spectral resolution. Finally, the proposed framework was validated using Monte Carlo simulations and experimental data acquired from healthy participants on an MRI scanner. The results demonstrated that the suggested approach is feasible for offering high-resolution DRC maps that correspond to distinct microstructures with a limited amount of optimised acquisition data from two-dimensional DRC sampling space. By significantly reducing scan time while retaining structural resolution, this approach may enable multi-dimensional DRC MRI to be more widely used for quantitative evaluation in biological and medical settings.

In Vivo Myelin Water Quantification Using Diffusion-Relaxation Correlation MRI: A Comparison of 1D and 2D Methods

Multidimensional Magnetic Resonance Imaging (MRI) is a versatile tool for microstructure mapping. We use a diffusion weighted inversion recovery spin echo (DW-IR-SE) sequence with spiral readouts at ultra-strong gradients to acquire a rich diffusion-relaxation data set with sensitivity to myelin water. We reconstruct 1D and 2D spectra with a two-step convex optimization approach and investigate a variety of multidimensional MRI methods, including 1D multi-component relaxometry, 1D multi-component diffusometry, 2D relaxation correlation imaging, and 2D diffusion-relaxation correlation spectroscopic imaging (DR-CSI), in terms of their potential to quantify tissue microstructure, including the myelin water fraction (MWF). We observe a distinct spectral peak that we attribute to myelin water in multi-component T1 relaxometry, T1-T2 correlation, T1-D correlation, and T2-D correlation imaging. Due to lower achievable echo times compared to diffusometry, MWF maps from relaxometry have higher quality. Whilst 1D multi-component T1 data allows much faster myelin mapping, 2D approaches could offer unique insights into tissue microstructure and especially myelin diffusion.

Enabling early detection of knee osteoarthritis using diffusion-relaxation correlation spectrum imaging

AIM

To present a technique that enables detection of early stage OA of the knee using diffusion-relaxation correlation spectrum imaging (DR-CSI).

MATERIALS AND METHODS

Fifty-five early osteoarthritis patients (OA, Kellgren-Lawrence [KL] score 1 to 2; mean age, 56.4 years) and 49 healthy volunteers (mean age, 56.7 years) were underwent magnetic resonance imaging (MRI) with T2-mapping and DR-CSI techniques. Maps of mean apparent diffusion coefficient (ADC), T2 relaxation time and volume fraction Vi for DR-CSI compartment i (A, B, C, D) sensitivity, specificity, and positive and negative likelihood ratio (PLR, NLR) were assessed to determine the diagnostic accuracy for detection of early-stage degeneration of knee articular cartilage. The structural abnormalities of articular cartilage were evaluated using modified Whole-Organ MR Imaging Scores (WORMS).

RESULTS

All intra- and interobserver agreements for DR-CSI compartment volume fractions and modified WORMS of cartilage were excellent. Early OA versus the controls had higher V_C, lower V_A and V_B (p<0.001), but comparable V_D (p>0.05). V_A, V_B and V_C had a moderate association with WORMS. No significant correlation was identified between V_D and WORMS. V_C had better ability than V_A,V_B, V_D, T2 and ADC to discriminate early OA patients from healthy controls (area under the curve, 0.898). Sensitivity, specificity, PLR, and NLR of V_C with a cut-off value of 29.9% were 81.8% (95% confidence interval [CI], 69.1-90.9%), 95.9% (86-99.5%), 20.05% (5.13-78.34%), and 0.19% (0.11-0.33%).

CONCLUSIONS

DR-CSI compartment volume fractions may be sensitive indicators for detecting early-stage degeneration in knee articular cartilage.

Deep Learning Initialized Compressed Sensing (Deli-CS) in Volumetric Spatio-Temporal Subspace Reconstruction

Introduction

Spatio-temporal MRI methods enable whole-brain multi-parametric mapping at ultra-fast acquisition times through efficient k-space encoding, but can have very long reconstruction times, which limit their integration into clinical practice. Deep learning (DL) is a promising approach to accelerate reconstruction, but can be computationally intensive to train and deploy due to the large dimensionality of spatio-temporal MRI. DL methods also need large training data sets and can produce results that don't match the acquired data if data consistency is not enforced. The aim of this project is to reduce reconstruction time using DL whilst simultaneously limiting the risk of deep learning induced hallucinations, all with modest hardware requirements.

Methods

Deep Learning Initialized Compressed Sensing (Deli-CS) is proposed to reduce the reconstruction time of iterative reconstructions by "kick-starting" the iterative reconstruction with a DL generated starting point. The proposed framework is applied to volumetric multi-axis spiral projection MRF that achieves whole-brain T1 and T2 mapping at 1-mm isotropic resolution for a 2-minute acquisition. First, the traditional reconstruction is optimized from over two hours to less than 40 minutes while using more than 90% less RAM and only 4.7 GB GPU memory, by using a memory-efficient GPU implementation. The Deli-CS framework is then implemented and evaluated against the above reconstruction.

Results

Deli-CS achieves comparable reconstruction quality with 50% fewer iterations bringing the full reconstruction time to 20 minutes.

Conclusion

Deli-CS reduces the reconstruction time of subspace reconstruction of volumetric spatio-temporal acquisitions by providing a warm start to the iterative reconstruction algorithm.

Nonparametric 5D D-R2 distribution imaging with single-shot EPI at 21.1 T: Initial results for in vivo rat brain

In vivo human diffusion MRI is by default performed using single-shot EPI with greater than 50-ms echo times and associated signal loss from transverse relaxation. The individual benefits of the current trends of increasing B₀ to boost SNR and employing more advanced signal preparation schemes to improve the specificity for selected microstructural properties eventually may be cancelled by increased relaxation rates at high B₀ and echo times with advanced encoding. Here, initial attempts to translate state-of-the-art diffusion-relaxation correlation methods from 3 T to 21.1 T are made to identify hurdles that need to be overcome to fulfill the promises of both high SNR and readily interpretable microstructural information.

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.

Nonparametric D-R1-R2 distribution MRI of the living human brain

Diffusion-relaxation correlation NMR can simultaneously characterize both the microstructure and the local chemical composition of complex samples that contain multiple populations of water. Recent developments on tensor-valued diffusion encoding and Monte Carlo inversion algorithms have made it possible to transfer diffusion-relaxation correlation NMR from small-bore scanners to clinical MRI systems. Initial studies on clinical MRI systems employed 5D D-R₁ and D-R₂ correlation to characterize healthy brain in vivo. However, these methods are subject to an inherent bias that originates from not including R₂ or R₁ in the analysis, respectively. This drawback can be remedied by extending the concept to 6D D-R₁-R₂ correlation. In this work, we present a sparse acquisition protocol that records all data necessary for in vivo 6D D-R₁-R₂ correlation MRI across 633 individual measurements within 25 min-a time frame comparable to previous lower-dimensional acquisition protocols. The data were processed with a Monte Carlo inversion algorithm to obtain nonparametric 6D D-R₁-R₂ distributions. We validated the reproducibility of the method in repeated measurements of healthy volunteers. For a post-therapy glioblastoma case featuring cysts, edema, and partially necrotic remains of tumor, we present representative single-voxel 6D distributions, parameter maps, and artificial contrasts over a wide range of diffusion-, R₁-, and R₂-weightings based on the rich information contained in the D-R₁-R₂ distributions.

Combined diffusion-relaxometry microstructure imaging: Current status and future prospects

Microstructure imaging seeks to noninvasively measure and map microscopic tissue features by pairing mathematical modeling with tailored MRI protocols. This article reviews an emerging paradigm that has the potential to provide a more detailed assessment of tissue microstructure-combined diffusion-relaxometry imaging. Combined diffusion-relaxometry acquisitions vary multiple MR contrast encodings-such as b-value, gradient direction, inversion time, and echo time-in a multidimensional acquisition space. When paired with suitable analysis techniques, this enables quantification of correlations and coupling between multiple MR parameters-such as diffusivity, T 1 , T 2 , and T 2 ∗ . This opens the possibility of disentangling multiple tissue compartments (within voxels) that are indistinguishable with single-contrast scans, enabling a new generation of microstructural maps with improved biological sensitivity and specificity.

Multidimensional MRI for Characterization of Subtle Axonal Injury Accelerated Using an Adaptive Nonlocal Multispectral Filter

Multidimensional MRI is an emerging approach that simultaneously encodes water relaxation (T₁ and T₂) and mobility (diffusion) and replaces voxel-averaged values with subvoxel distributions of those MR properties. While conventional (i.e., voxel-averaged) MRI methods cannot adequately quantify the microscopic heterogeneity of biological tissue, using subvoxel information allows to selectively map a specific T₁-T₂-diffusion spectral range that corresponds to a group of tissue elements. The major obstacle to the adoption of rich, multidimensional MRI protocols for diagnostic or monitoring purposes is the prolonged scan time. Our main goal in the present study is to evaluate the performance of a nonlocal estimation of multispectral magnitudes (NESMA) filter on reduced datasets to limit the total acquisition time required for reliable multidimensional MRI characterization of the brain. Here we focused and reprocessed results from a recent study that identified potential imaging biomarkers of axonal injury pathology from the joint analysis of multidimensional MRI, in particular voxelwise T₁-T₂ and diffusion-T₂ spectra in human Corpus Callosum, and histopathological data. We tested the performance of NESMA and its effect on the accuracy of the injury biomarker maps, relative to the co-registered histological reference. Noise reduction improved the accuracy of the resulting injury biomarker maps, while permitting data reduction of 35.7 and 59.6% from the full dataset for T₁-T₂ and diffusion-T₂ cases, respectively. As successful clinical proof-of-concept applications of multidimensional MRI are continuously being introduced, reliable and robust noise removal and consequent acquisition acceleration would advance the field towards clinically-feasible diagnostic multidimensional MRI protocols.

Diffuse axonal injury has a characteristic multidimensional MRI signature in the human brain

Axonal injury is a major contributor to the clinical symptomatology in patients with traumatic brain injury. Conventional neuroradiological tools, such as CT and MRI, are insensitive to diffuse axonal injury (DAI) caused by trauma. Diffusion tensor MRI parameters may change in DAI lesions; however, the nature of these changes is inconsistent. Multidimensional MRI is an emerging approach that combines T1, T2, and diffusion, and replaces voxel-averaged values with distributions, which allows selective isolation of specific potential abnormal components. By performing a combined post-mortem multidimensional MRI and histopathology study, we aimed to investigate T1-T2-diffusion changes linked to DAI and to define their histopathological correlates. Corpora callosa derived from eight subjects who had sustained traumatic brain injury, and three control brain donors underwent post-mortem ex vivo MRI at 7 T. Multidimensional, diffusion tensor, and quantitative T1 and T2 MRI data were acquired and processed. Following MRI acquisition, slices from the same tissue were tested for amyloid precursor protein (APP) immunoreactivity to define DAI severity. A robust image co-registration method was applied to accurately match MRI-derived parameters and histopathology, after which 12 regions of interest per tissue block were selected based on APP density, but blind to MRI. We identified abnormal multidimensional T1-T2, diffusion-T2, and diffusion-T1 components that are strongly associated with DAI and used them to generate axonal injury images. We found that compared to control white matter, mild and severe DAI lesions contained significantly larger abnormal T1-T2 component (P = 0.005 and P < 0.001, respectively), and significantly larger abnormal diffusion-T2 component (P = 0.005 and P < 0.001, respectively). Furthermore, within patients with traumatic brain injury the multidimensional MRI biomarkers differentiated normal-appearing white matter from mild and severe DAI lesions, with significantly larger abnormal T1-T2 and diffusion-T2 components (P = 0.003 and P < 0.001, respectively, for T1-T2; P = 0.022 and P < 0.001, respectively, for diffusion-T2). Conversely, none of the conventional quantitative MRI parameters were able to differentiate lesions and normal-appearing white matter. Lastly, we found that the abnormal T1-T2, diffusion-T1, and diffusion-T2 components and their axonal damage images were strongly correlated with quantitative APP staining (r = 0.876, P < 0.001; r = 0.727, P < 0.001; and r = 0.743, P < 0.001, respectively), while producing negligible intensities in grey matter and in normal-appearing white matter. These results suggest that multidimensional MRI may provide non-invasive biomarkers for detection of DAI, which is the pathological substrate for neurological disorders ranging from concussion to severe traumatic brain injury.

Whole-Brain Imaging of Subvoxel T1-Diffusion Correlation Spectra in Human Subjects

T1 relaxation and water mobility generate eloquent MRI tissue contrasts with great diagnostic value in many neuroradiological applications. However, conventional methods do not adequately quantify the microscopic heterogeneity of these important biophysical properties within a voxel, and therefore have limited biological specificity. We describe a new correlation spectroscopic (CS) MRI method for measuring how T1 and mean diffusivity (MD) co-vary in microscopic tissue environments. We develop a clinical pulse sequence that combines inversion recovery (IR) with single-shot isotropic diffusion encoding (IDE) to efficiently acquire whole-brain MRIs with a wide range of joint T1-MD weightings. Unlike conventional diffusion encoding, the IDE preparation ensures that all subvoxel water pools are weighted by their MDs regardless of the sizes, shapes, and orientations of their corresponding microscopic diffusion tensors. Accordingly, IR-IDE measurements are well-suited for model-free, quantitative spectroscopic analysis of microscopic water pools. Using numerical simulations, phantom experiments, and data from healthy volunteers we demonstrate how IR-IDE MRIs can be processed to reconstruct maps of two-dimensional joint probability density functions, i.e., correlation spectra, of subvoxel T1-MD values. In vivo T1-MD spectra show distinct cerebrospinal fluid and parenchymal tissue components specific to white matter, cortical gray matter, basal ganglia, and myelinated fiber pathways, suggesting the potential for improved biological specificity. The one-dimensional marginal distributions derived from the T1-MD correlation spectra agree well with results from other relaxation spectroscopic and quantitative MRI studies, validating the T1-MD contrast encoding and the spectral reconstruction. Mapping subvoxel T1-diffusion correlations in patient populations may provide a more nuanced, comprehensive, sensitive, and specific neuroradiological assessment of the non-specific changes seen on fluid-attenuated inversion recovery (FLAIR) and diffusion-weighted MRIs (DWIs) in cancer, ischemic stroke, or brain injury.

Toward nonparametric diffusion- T1 characterization of crossing fibers in the human brain

PURPOSE

To estimate T 1 for each distinct fiber population within voxels containing multiple brain tissue types.

METHODS

A diffusion- T 1 correlation experiment was carried out in an in vivo human brain using tensor-valued diffusion encoding and multiple repetition times. The acquired data were inverted using a Monte Carlo algorithm that retrieves nonparametric distributions P ( D , R 1 ) of diffusion tensors and longitudinal relaxation rates R 1 = 1 / T 1 . Orientation distribution functions (ODFs) of the highly anisotropic components of P ( D , R 1 ) were defined to visualize orientation-specific diffusion-relaxation properties. Finally, Monte Carlo density-peak clustering (MC-DPC) was performed to quantify fiber-specific features and investigate microstructural differences between white matter fiber bundles.

RESULTS

Parameter maps corresponding to P ( D , R 1 ) 's statistical descriptors were obtained, exhibiting the expected R 1 contrast between brain tissue types. Our ODFs recovered local orientations consistent with the known anatomy and indicated differences in R 1 between major crossing fiber bundles. These differences, confirmed by MC-DPC, were in qualitative agreement with previous model-based works but seem biased by the limitations of our current experimental setup.

CONCLUSIONS

Our Monte Carlo framework enables the nonparametric estimation of fiber-specific diffusion- T 1 features, thereby showing potential for characterizing developmental or pathological changes in T 1 within a given fiber bundle, and for investigating interbundle T 1 differences.

Region-optimized virtual (ROVir) coils: Localization and/or suppression of spatial regions using sensor-domain beamforming

PURPOSE

In many MRI scenarios, magnetization is often excited from spatial regions that are not of immediate interest. Excitation of uninteresting magnetization can complicate the design of efficient imaging methods, leading to either artifacts or acquisitions that are longer than necessary. While there are many hardware- and sequence-based approaches for suppressing unwanted magnetization, this paper approaches this longstanding problem from a different and complementary angle, using beamforming to suppress signals from unwanted regions without modifying the acquisition hardware or pulse sequence. Unlike existing beamforming approaches, we use a spatially invariant sensor-domain approach that can be applied directly to raw data to facilitate image reconstruction.

THEORY AND METHODS

We use beamforming to linearly mix a set of original coils into a set of region-optimized virtual (ROVir) coils. ROVir coils optimize a signal-to-interference ratio metric, are easily calculated using simple generalized eigenvalue decomposition methods, and provide coil compression.

RESULTS

ROVir coils were compared against existing coil-compression methods, and were demonstrated to have substantially better signal suppression capabilities. In addition, examples were provided in a variety of different application contexts (including brain MRI, vocal tract MRI, and cardiac MRI; accelerated Cartesian and non-Cartesian imaging; and outer volume suppression) that demonstrate the strong potential of this kind of approach.

CONCLUSION

The beamforming-based ROVir framework is simple to implement, has promising capabilities to suppress unwanted MRI signal, and is compatible with and complementary to existing signal suppression methods. We believe that this general approach could prove useful across a wide range of different MRI applications.

Diffusion-PEPTIDE: Distortion- and blurring-free diffusion imaging with self-navigated motion-correction and relaxometry capabilities

PURPOSE

To implement the time-resolved relaxometry PEPTIDE technique into a diffusion acquisition to provide self-navigated, distortion- and blurring-free diffusion imaging that is robust to motion, while simultaneously providing T₂ and T 2 ∗ mapping.

THEORY AND METHODS

The PEPTIDE readout was implemented into a spin-echo diffusion acquisition, enabling reconstruction of a time-series of T₂ - and T 2 ∗ -weighted images, free from conventional echo planar imaging (EPI) distortion and blurring, for each diffusion-encoding. Robustness of PEPTIDE to motion and shot-to-shot phase variation was examined through a deliberate motion-corrupted diffusion experiment. Two diffusion-relaxometry in vivo brain protocols were also examined: (1)1 × 1 × 3 mm³ across 32 diffusion directions in 20 min, (2)1.5 × 1.5 × 3.0 mm³ across 6 diffusion-weighted images in 3.4 min. T₂ , T 2 ∗ , and diffusion parameter maps were calculated from these data. As initial exploration of the rich diffusion-relaxometry data content for use in multi-compartment modeling, PEPTIDE data were acquired of a gadolinium-doped asparagus phantom. These datasets contained two compartments with different relaxation parameters and different diffusion orientation properties, and T₂ relaxation variations across these diffusion directions were explored.

RESULTS

Diffusion-PEPTIDE showed the capability to provide high quality diffusion images and T₂ and T 2 ∗ maps from both protocols. The reconstructions were distortion-free, avoided potential resolution losses exceeding 100% in equivalent EPI acquisitions, and showed tolerance to nearly 30° of rotational motion. Expected variation in T₂ values as a function of diffusion direction was observed in the two-compartment asparagus phantom (P < .01), demonstrating potential to explore diffusion-PEPTIDE data for multi-compartment modeling.

CONCLUSIONS

Diffusion-PEPTIDE provides highly robust diffusion and relaxometry data and offers potential for future applications in diffusion-relaxometry multi-compartment modeling.

Fast multicomponent 3D-T1ρ relaxometry

NMR relaxometry can provide information about the relaxation of the magnetization in different tissues, increasing our understanding of molecular dynamics and biochemical composition in biological systems. In general, tissues have complex and heterogeneous structures composed of multiple pools. As a result, bulk magnetization returns to its original state with different relaxation times, in a multicomponent relaxation. Recovering the distribution of relaxation times in each voxel is a difficult inverse problem; it is usually unstable and requires long acquisition time, especially on clinical scanners. MRI can also be viewed as an inverse problem, especially when compressed sensing (CS) is used. The solution of these two inverse problems, CS and relaxometry, can be obtained very efficiently in a synergistically combined manner, leading to a more stable multicomponent relaxometry obtained with short scan times. In this paper, we will discuss the details of this technique from the viewpoint of inverse problems.

Three-dimensional whole-brain simultaneous T1, T2, and T1ρ quantification using MR Multitasking: Method and initial clinical experience in tissue characterization of multiple sclerosis

PURPOSE

To develop a 3D whole-brain simultaneous T1/T2/T1ρ quantification method with MR Multitasking that provides high quality, co-registered multiparametric maps in 9 min.

METHODS

MR Multitasking conceptualizes T1/T2/T1ρ relaxations as different time dimensions, simultaneously resolving all three dimensions with a low-rank tensor image model. The proposed method was validated on a phantom and in healthy volunteers, comparing quantitative measurements against corresponding reference methods and evaluating the scan-rescan repeatability. Initial clinical validation was performed in age-matched relapsing-remitting multiple sclerosis (RRMS) patients to examine the feasibility of quantitative tissue characterization and to compare with the healthy control cohort. The feasibility of synthesizing six contrast-weighted images was also examined.

RESULTS

Our framework produced high quality, co-registered T1/T2/T1ρ maps that closely resemble the reference maps. Multitasking T1/T2/T1ρ measurements showed substantial agreement with reference measurements on the phantom and in healthy controls. Bland-Altman analysis indicated good in vivo repeatability of all three parameters. In RRMS patients, lesions were conspicuously delineated on all three maps and on four synthetic weighted images (T2-weighted, T2-FLAIR, double inversion recovery, and a novel "T1ρ-FLAIR" contrast). T1 and T2 showed significant differences for normal appearing white matter between patients and controls, while T1ρ showed significant differences for normal appearing white matter, cortical gray matter, and deep gray matter. The combination of three parameters significantly improved the differentiation between RRMS patients and healthy controls, compared to using any single parameter alone.

CONCLUSION

MR Multitasking simultaneously quantifies whole-brain T1/T2/T1ρ and is clinically promising for quantitative tissue characterization of neurological diseases, such as MS.