Integrated and efficient diffusion-relaxometry using ZEBRA

Time-division multiplexing (TDM) sequence removes bias in T2 estimation and relaxation-diffusion measurements

PURPOSE

To compare the performance of multi-echo (ME) and time-division multiplexing (TDM) sequences for accelerated relaxation-diffusion MRI (rdMRI) acquisition and to examine their reliability in estimating accurate rdMRI microstructure measures.

METHOD

The ME, TDM, and the reference single-echo (SE) sequences with six TEs were implemented using Pulseq with single-band (SB) and multi-band 2 (MB2) acceleration factors. On a diffusion phantom, the image intensities of the three sequences were compared, and the differences were quantified using the normalized RMS error (NRMSE). Shinnar-Le Roux (SLR) pulses were implemented for the SB-ME and SB-SE sequences to investigate the impact of slice profiles on ME sequences. For the in-vivo brain scan, besides the image intensity comparison and T2-estimates, different methods were used to assess sequence-related effects on microstructure estimation, including the relaxation diffusion imaging moment (REDIM) and the maximum-entropy relaxation diffusion distribution (MaxEnt-RDD).

RESULTS

TDM performance was similar to the gold standard SE acquisition, whereas ME showed greater biases (3-4× larger NRMSEs for phantom, 2× for in-vivo). T2 values obtained from TDM closely matched SE, whereas ME sequences underestimated the T2 relaxation time. TDM provided similar diffusion and relaxation parameters as SE using REDIM, whereas SB-ME exhibited a 60% larger bias in the  map and on average 3.5× larger bias in the covariance between relaxation-diffusion coefficients.

CONCLUSION

Our analysis demonstrates that TDM provides a more accurate estimation of relaxation-diffusion measurements while accelerating the acquisitions by a factor of 2 to 3.

DeepEMC-T2 mapping: Deep learning-enabled T2 mapping based on echo modulation curve modeling

PURPOSE

Echo modulation curve (EMC) modeling enables accurate quantification of T2 relaxation times in multi-echo spin-echo (MESE) imaging. The standard EMC-T2 mapping framework, however, requires sufficient echoes and cumbersome pixel-wise dictionary-matching steps. This work proposes a deep learning version of EMC-T2 mapping, called DeepEMC-T2 mapping, to efficiently estimate accurate T2 maps from fewer echoes.

METHODS

DeepEMC-T2 mapping was developed using a modified U-Net to estimate both T2 and proton density (PD) maps directly from MESE images. The network implements several new features to improve the accuracy of T2/PD estimation. A total of 67 MESE datasets acquired in axial orientation were used for network training and evaluation. An additional 57 datasets acquired in coronal orientation with different scan parameters were used to evaluate the generalizability of the framework. The performance of DeepEMC-T2 mapping was evaluated in seven experiments.

RESULTS

Compared to the reference, DeepEMC-T2 mapping achieved T2 estimation errors from 1% to 11% and PD estimation errors from 0.4% to 1.5% with ten/seven/five/three echoes, which are more accurate than standard EMC-T2 mapping. By incorporating datasets acquired with different scan parameters and orientations for joint training, DeepEMC-T2 exhibits robust generalizability across varying imaging protocols. Increasing the echo spacing and including longer echoes improve the accuracy of parameter estimation. The new features proposed in DeepEMC-T2 mapping all enabled more accurate T2 estimation.

CONCLUSIONS

DeepEMC-T2 mapping enables simplified, efficient, and accurate T2 quantification directly from MESE images without dictionary matching. Accurate T2 estimation from fewer echoes allows for increased volumetric coverage and/or higher slice resolution without prolonging total scan times.

Placental multimodal MRI prior to spontaneous preterm birth <32 weeks' gestation: An observational study

OBJECTIVE

To utilise combined diffusion-relaxation MRI techniques to interrogate antenatal changes in the placenta prior to extreme preterm birth among both women with PPROM and membranes intact, and compare this to a control group who subsequently delivered at term.

DESIGN

Observational study.

SETTING

Tertiary Obstetric Unit, London, UK.

POPULATION

Cases: pregnant women who subsequently spontaneously delivered a singleton pregnancy prior to 32 weeks' gestation without any other obstetric complications.

CONTROLS

pregnant women who delivered an uncomplicated pregnancy at term.

METHODS

All women consented to an MRI examination. A combined diffusion-relaxation MRI of the placenta was undertaken and analysed using fractional anisotropy, a combined T2*-apparent diffusion coefficient model and a combined T2*-intravoxel incoherent motion model, in order to provide a detailed placental phenotype associated with preterm birth. Subgroup analyses based on whether women in the case group had PPROM or intact membranes at time of scan, and on latency to delivery were performed.

MAIN OUTCOME MEASURES

Fractional anisotropy, apparent diffusion coefficients and T2* placental values, from two models including a combined T2*-IVIM model separating fast- and slow-flowing (perfusing and diffusing) compartments.

RESULTS

This study included 23 women who delivered preterm and 52 women who delivered at term. Placental T2* was lower in the T2*-apparent diffusion coefficient model (p < 0.001) and in the fast- and slow-flowing compartments (p = 0.001 and p < 0.001) of the T2*-IVIM model. This reached a higher level of significance in the preterm prelabour rupture of the membranes group than in the membranes intact group. There was a reduced perfusion fraction among the cases with impending delivery.

CONCLUSIONS

Placental diffusion-relaxation reveals significant changes in the placenta prior to preterm birth with greater effect noted in cases of preterm prelabour rupture of the membranes. Application of this technique may allow clinically valuable interrogation of histopathological changes before preterm birth. In turn, this could facilitate more accurate antenatal prediction of preterm chorioamnionitis and so aid decisions around the safest time of delivery. Furthermore, this technique provides a research tool to improve understanding of the pathological mechanisms associated with preterm birth in vivo.

T1 mapping from routine 3D T1-weighted inversion recovery sequences in clinical practice: comparison against reference inversion recovery fast field echo T1 scans and feasibility in multiple sclerosis

BACKGROUND AND PURPOSE

Quantitative T1 mapping can be an essential tool for assessing tissue injury in multiple sclerosis (MS). We introduce T1-REQUIRE, a method that converts a single high-resolution anatomical 3D T1-weighted Turbo Field Echo (3DT1TFE) scan into a parametric T1 map that could be used for quantitative assessment of tissue damage. We present the accuracy and feasibility of this method in MS.

METHODS

14 subjects with relapsing-remitting MS and 10 healthy subjects were examined. T1 maps were generated from 3DT1TFE images using T1-REQUIRE, which estimates T1 values using MR signal equations and internal tissue reference T1 values. Estimated T1 of lesions, white, and gray matter regions were compared with reference Inversion-Recovery Fast Field Echo T1 values and analyzed via correlation and Bland-Altman (BA) statistics.

RESULTS

159 T1-weighted (T1W) hypointense MS lesions and 288 gray matter regions were examined. T1 values for MS lesions showed a Pearson's correlation of r = 0.81 (p < 0.000), R2 = 0.65, and Bias = 4.18%. BA statistics showed a mean difference of -53.95 ms and limits of agreement (LOA) of -344.20 and 236.30 ms. Non-lesional normal-appearing white matter had a correlation coefficient of r = 0.82 (p < 0.000), R2 = 0.67, Bias = 8.78%, mean difference of 73.87 ms, and LOA of -55.67 and 203.41 ms.

CONCLUSIONS

We demonstrate the feasibility of retroactively derived high-resolution T1 maps from routinely acquired anatomical images, which could be used to quantify tissue pathology in MS. The results of this study will set the stage for testing this method in larger clinical studies for examining MS disease activity and progression.

Time-division multiplexing (TDM) sequence removes bias in T2 estimation and relaxation-diffusion measurements

Purpose

To compare the performance of multi-echo (ME) and time-division multiplexing (TDM) sequences for accelerated relaxation-diffusion MRI (rdMRI) acquisition and to examine their reliability in estimating accurate rdMRI microstructure measures.

Method

The ME, TDM, and the reference single-echo (SE) sequences with six echo times (TE) were implemented using Pulseq with single-band (SB-) and multi-band 2 (MB2-) acceleration factors. On a diffusion phantom, the image intensities of the three sequences were compared, and the differences were quantified using the normalized root mean squared error (NRMSE). For the in-vivo brain scan, besides the image intensity comparison and T2-estimates, different methods were used to assess sequence-related effects on microstructure estimation, including the relaxation diffusion imaging moment (REDIM) and the maximum-entropy relaxation diffusion distribution (MaxEnt-RDD).

Results

TDM performance was similar to the gold standard SE acquisition, whereas ME showed greater biases (3-4× larger NRMSEs for phantom, 2× for in-vivo). T2 values obtained from TDM closely matched SE, whereas ME sequences underestimated the T2 relaxation time. TDM provided similar diffusion and relaxation parameters as SE using REDIM, whereas SB-ME exhibited a 60% larger bias in the map and on average 3.5× larger bias in the covariance between relaxation-diffusion coefficients.

Conclusion

Our analysis demonstrates that TDM provides a more accurate estimation of relaxation-diffusion measurements while accelerating the acquisitions by a factor of 2 to 3.

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.

Advanced magnetic resonance imaging detects altered placental development in pregnancies affected by congenital heart disease

Congenital heart disease (CHD) is the most common congenital malformation and is associated with adverse neurodevelopmental outcomes. The placenta is crucial for healthy fetal development and placental development is altered in pregnancy when the fetus has CHD. This study utilized advanced combined diffusion-relaxation MRI and a data-driven analysis technique to test the hypothesis that placental microstructure and perfusion are altered in CHD-affected pregnancies. 48 participants (36 controls, 12 CHD) underwent 67 MRI scans (50 control, 17 CHD). Significant differences in the weighting of two independent placental and uterine-wall tissue components were identified between the CHD and control groups (both pFDR < 0.001), with changes most evident after 30 weeks gestation. A significant trend over gestation in weighting for a third independent tissue component was also observed in the CHD cohort (R = 0.50, pFDR = 0.04), but not in controls. These findings add to existing evidence that placental development is altered in CHD. The results may reflect alterations in placental perfusion or the changes in fetal-placental flow, villous structure and maturation that occur in CHD. Further research is needed to validate and better understand these findings and to understand the relationship between placental development, CHD, and its neurodevelopmental implications.

Optimisation of quantitative brain diffusion-relaxation MRI acquisition protocols with physics-informed machine learning

Diffusion-relaxation MRI aims to extract quantitative measures that characterise microstructural tissue properties such as orientation, size, and shape, but long acquisition times are typically required. This work proposes a physics-informed learning framework to extract an optimal subset of diffusion-relaxation MRI measurements for enabling shorter acquisition times, predict non-measured signals, and estimate quantitative parameters. In vivo and synthetic brain 5D-Diffusion-T1-T2∗-weighted MRI data obtained from five healthy subjects were used for training and validation, and from a sixth participant for testing. One fully data-driven and two physics-informed machine learning methods were implemented and compared to two manual selection procedures and Cramér-Rao lower bound optimisation. The physics-informed approaches could identify measurement-subsets that yielded more consistently accurate parameter estimates in simulations than other approaches, with similar signal prediction error. Five-fold shorter protocols yielded error distributions of estimated quantitative parameters with very small effect sizes compared to estimates from the full protocol. Selected subsets commonly included a denser sampling of the shortest and longest inversion time, lowest echo time, and high b-value. The proposed framework combining machine learning and MRI physics offers a promising approach to develop shorter imaging protocols without compromising the quality of parameter estimates and signal predictions.

What contributes to disability in progressive MS? A brain and cervical cord-matched quantitative MRI study

BACKGROUND

We assessed the ability of a brain-and-cord-matched quantitative magnetic resonance imaging (qMRI) protocol to differentiate patients with progressive multiple sclerosis (PMS) from controls, in terms of normal-appearing (NA) tissue abnormalities, and explain disability.

METHODS

A total of 27 patients and 16 controls were assessed on the Expanded Disability Status Scale (EDSS), 25-foot timed walk (TWT), 9-hole peg (9HPT) and symbol digit modalities (SDMT) tests. All underwent 3T brain and (C2-C3) cord structural imaging and qMRI (relaxometry, quantitative magnetisation transfer, multi-shell diffusion-weighted imaging), using a fast brain-and-cord-matched protocol with brain-and-cord-unified imaging readouts. Lesion and NA-tissue volumes and qMRI metrics reflecting demyelination and axonal loss were obtained. Random forest analyses identified the most relevant volumetric/qMRI measures to clinical outcomes. Confounder-adjusted linear regression estimated the actual MRI-clinical associations.

RESULTS

Several qMRI/volumetric differences between patients and controls were observed (p < 0.01). Higher NA-deep grey matter quantitative-T1 (EDSS: beta = 7.96, p = 0.006; 9HPT: beta = -0.09, p = 0.004), higher NA-white matter orientation dispersion index (TWT: beta = -3.21, p = 0.005; SDMT: beta = -847.10, p < 0.001), lower whole-cord bound pool fraction (9HPT: beta = 0.79, p = 0.001) and higher NA-cortical grey matter quantitative-T1 (SDMT = -94.31, p < 0.001) emerged as particularly relevant predictors of greater disability.

CONCLUSION

Fast brain-and-cord-matched qMRI protocols are feasible and identify demyelination - combined with other mechanisms - as key for disability accumulation in PMS.

Advanced magnetic resonance imaging detects altered placental development in pregnancies affected by congenital heart disease

Congenital heart disease (CHD) is the most common congenital malformation and is associated with adverse neurodevelopmental outcomes. The placenta is crucial for healthy fetal development and placental development is altered in pregnancy when the fetus has CHD. This study utilized advanced combined diffusion-relaxation MRI and a data-driven analysis technique to test the hypothesis that placental microstructure and perfusion are altered in CHD-affected pregnancies. 48 participants (36 controls, 12 CHD) underwent 67 MRI scans (50 control, 17 CHD). Significant differences in the weighting of two independent placental and uterine-wall tissue components were identified between the CHD and control groups (both pFDR<0.001), with changes most evident after 30 weeks gestation. A Significant trend over gestation in weighting for a third independent tissue component was also observed in the CHD cohort (R = 0.50, pFDR=0.04), but not in controls. These findings add to existing evidence that placental development is altered in CHD. The results may reflect alterations in placental perfusion or the changes in fetal-placental flow, villous structure and maturation that occur in CHD. Further research is needed to validate and better understand these findings and to understand the relationship between placental development, CHD, and its neurodevelopmental implications.

Non-invasive mapping of human placenta microenvironments throughout pregnancy with diffusion-relaxation MRI

INTRODUCTION

In-vivo measurements of placental structure and function have the potential to improve prediction, diagnosis, and treatment planning for a wide range of pregnancy complications, such as fetal growth restriction and pre-eclampsia, and hence inform clinical decision making, ultimately improving patient outcomes. MRI is emerging as a technique with increased sensitivity to placental structure and function compared to the current clinical standard, ultrasound.

METHODS

We demonstrate and evaluate a combined diffusion-relaxation MRI acquisition and analysis pipeline on a sizable cohort of 78 normal pregnancies with gestational ages ranging from 15 + 5 to 38 + 4 weeks. Our acquisition comprises a combined T2*-diffusion MRI acquisition sequence - which is simultaneously sensitive to oxygenation, microstructure and microcirculation. We analyse our scans with a data-driven unsupervised machine learning technique, InSpect, that parsimoniously identifies distinct components in the data.

RESULTS

We identify and map seven potential placental microenvironments and reveal detailed insights into multiple microstructural and microcirculatory features of the placenta, and assess their trends across gestation.

DISCUSSION

By demonstrating direct observation of micro-scale placental structure and function, and revealing clear trends across pregnancy, our work contributes towards the development of robust imaging biomarkers for pregnancy complications and the ultimate goal of a normative model of placental development.

Probing brain tissue microstructure with MRI: principles, challenges, and the role of multidimensional diffusion-relaxation encoding

Diffusion MRI uses the random displacement of water molecules to sensitize the signal to brain microstructure and to properties such as the density and shape of cells. Microstructure modeling techniques aim to estimate these properties from acquired data by separating the signal between virtual tissue 'compartments' such as the intra-neurite and the extra-cellular space. A key challenge is that the diffusion MRI signal is relatively featureless compared with the complexity of brain tissue. Another challenge is that the tissue microstructure is wildly different within the gray and white matter of the brain. In this review, we use results from multidimensional diffusion encoding techniques to discuss these challenges and their tentative solutions. Multidimensional encoding increases the information content of the data by varying not only the b-value and the encoding direction but also additional experimental parameters such as the shape of the b-tensor and the echo time. Three main insights have emerged from such encoding. First, multidimensional data contradict common model assumptions on diffusion and T2 relaxation, and illustrates how the use of these assumptions cause erroneous interpretations in both healthy brain and pathology. Second, many model assumptions can be dispensed with if data are acquired with multidimensional encoding. The necessary data can be easily acquired in vivo using protocols optimized to minimize Cramér-Rao lower bounds. Third, microscopic diffusion anisotropy reflects the presence of axons but not dendrites. This insight stands in contrast to current 'neurite models' of brain tissue, which assume that axons in white matter and dendrites in gray matter feature highly similar diffusion. Nevertheless, as an axon-based contrast, microscopic anisotropy can differentiate gray and white matter when myelin alterations confound conventional MRI contrasts.

Mapping the individual human cortex using multidimensional MRI and unsupervised learning

Human evolution has seen the development of higher-order cognitive and social capabilities in conjunction with the unique laminar cytoarchitecture of the human cortex. Moreover, early-life cortical maldevelopment has been associated with various neurodevelopmental diseases. Despite these connections, there is currently no noninvasive technique available for imaging the detailed cortical laminar structure. This study aims to address this scientific and clinical gap by introducing an approach for imaging human cortical lamina. This method combines diffusion-relaxation multidimensional MRI with a tailored unsupervised machine learning approach that introduces enhanced microstructural sensitivity. This new imaging method simultaneously encodes the microstructure, the local chemical composition and importantly their correlation within complex and heterogenous tissue. To validate our approach, we compared the intra-cortical layers obtained using our ex vivo MRI-based method with those derived from Nissl staining of postmortem human brain specimens. The integration of unsupervised learning with diffusion-relaxation correlation MRI generated maps that demonstrate sensitivity to areal differences in cytoarchitectonic features observed in histology. Significantly, our observations revealed layer-specific diffusion-relaxation signatures, showing reductions in both relaxation times and diffusivities at the deeper cortical levels. These findings suggest a radial decrease in myelin content and changes in cell size and anisotropy, reflecting variations in both cytoarchitecture and myeloarchitecture. Additionally, we demonstrated that 1D relaxation and high-order diffusion MRI scalar indices, even when aggregated and used jointly in a multimodal fashion, cannot disentangle the cortical layers. Looking ahead, our technique holds the potential to open new avenues of research in human neurodevelopment and the vast array of disorders caused by disruptions in neurodevelopment.

The impact of head orientation with respect to B0 on diffusion tensor MRI measures

Diffusion tensor MRI (DT-MRI) remains the most commonly used approach to characterise white matter (WM) anisotropy. However, DT estimates may be affected by tissue orientation w.r.t. B → 0 due to local gradients and intrinsic T 2 orientation dependence induced by the microstructure. This work aimed to investigate whether and how diffusion tensor MRI-derived measures depend on the orientation of the head with respect to the static magnetic field, B → 0 . By simulating WM as two compartments, we demonstrated that compartmental T 2 anisotropy can induce the dependence of diffusion tensor measures on the angle between WM fibres and the magnetic field. In in vivo experiments, reduced radial diffusivity and increased axial diffusivity were observed in white matter fibres perpendicular to B → 0 compared to those parallel to B → 0 . Fractional anisotropy varied by up to 20 % as a function of the angle between WM fibres and the orientation of the main magnetic field. To conclude, fibre orientation w.r.t. B → 0 is responsible for up to 7 % variance in diffusion tensor measures across the whole brain white matter from all subjects and head orientations. Fibre orientation w.r.t. B → 0 may introduce additional variance in clinical research studies using diffusion tensor imaging, particularly when it is difficult to control for (e.g., fetal or neonatal imaging, or when the trajectories of fibres change due to, e.g., space occupying lesions).

Assessing within-subject rates of change of placental MRI diffusion metrics in normal pregnancy

PURPOSE

Studying placental development informs when development is abnormal. Most placental MRI studies are cross-sectional and do not study the extent of individual variability throughout pregnancy. We aimed to explore how diffusion MRI measures of placental function and microstructure vary in individual healthy pregnancies throughout gestation.

METHODS

Seventy-nine pregnant, low-risk participants (17 scanned twice and 62 scanned once) were included. T2 -weighted anatomical imaging and a combined multi-echo spin-echo diffusion-weighted sequence were acquired at 3 T. Combined diffusion-relaxometry models were performed using both aT 2 * $$ {\mathrm{T}}_2^{\ast } $$ -ADC and a bicompartmentalT 2 * $$ {\mathrm{T}}_2^{\ast } $$ -intravoxel-incoherent-motion (T 2 * IVIM $$ {\mathrm{T}}_2^{\ast}\;\mathrm{IVIM} $$ ) model fit.

RESULTS

There was a significant decline in placentalT 2 * $$ {\mathrm{T}}_2^{\ast } $$ and ADC (both P < 0.01) over gestation. These declines are consistent in individuals forT 2 * $$ {\mathrm{T}}_2^{\ast } $$ (covariance = -0.47), but not ADC (covariance = -1.04). TheT 2 * IVIM $$ {\mathrm{T}}_2^{\ast}\;\mathrm{IVIM} $$ model identified a consistent decline in individuals over gestation inT 2 * $$ {\mathrm{T}}_2^{\ast } $$ from both the perfusing and diffusing placental compartments, but not in ADC values from either. The placental perfusing compartment fraction increased over gestation (P = 0.0017), but this increase was not consistent in individuals (covariance = 2.57).

CONCLUSION

Whole placentalT 2 * $$ {\mathrm{T}}_2^{\ast } $$ and ADC values decrease over gestation, although onlyT 2 * $$ {\mathrm{T}}_2^{\ast } $$ values showed consistent trends within subjects. There was minimal individual variation in rates of change ofT 2 * $$ {\mathrm{T}}_2^{\ast } $$ values from perfusing and diffusing placental compartments, whereas trends in ADC values from these compartments were less consistent. These findings probably relate to the increased complexity of the bicompartmentalT 2 * IVIM $$ {\mathrm{T}}_2^{\ast}\;\mathrm{IVIM} $$ model, and differences in how different placental regions evolve at a microstructural level. These placental MRI metrics from low-risk pregnancies provide a useful benchmark for clinical cohorts.

Low-Field Combined Diffusion-Relaxation MRI for Mapping Placenta Structure and Function

Purpose

Demonstrating quantitative multi-parametric mapping in the placenta with combined T 2 ∗ -diffusion MRI at low-field (0.55T).

Methods

We present 57 placental MRI scans performed on a commercially available 0.55T scanner. We acquired the images using a combined T 2 ∗ -diffusion technique scan that simultaneously acquires multiple diffusion preparations and echo times. We processed the data to produce quantitative T 2 ∗ and diffusivity maps using a combined T 2 ∗ -ADC model. We compared the derived quantitative parameters across gestation in healthy controls and a cohort of clinical cases.

Results

Quantitative parameter maps closely resemble those from previous experiments at higher field strength, with similar trends in T 2 ∗ and ADC against gestational age observed.

Conclusion

Combined T 2 ∗ -diffusion placental MRI is reliably achievable at 0.55T. The advantages of lower field strength - such as cost, ease of deployment, increased accessibility and patient comfort due to the wider bore, and increased T 2 ∗ for larger dynamic ranges - can support the widespread roll out of placental MRI as an adjunct to ultrasound during pregnancy.

Spatiotemporal tissue maturation of thalamocortical pathways in the human fetal brain

The development of connectivity between the thalamus and maturing cortex is a fundamental process in the second half of human gestation, establishing the neural circuits that are the basis for several important brain functions. In this study, we acquired high-resolution in utero diffusion magnetic resonance imaging (MRI) from 140 fetuses as part of the Developing Human Connectome Project, to examine the emergence of thalamocortical white matter over the second to third trimester. We delineate developing thalamocortical pathways and parcellate the fetal thalamus according to its cortical connectivity using diffusion tractography. We then quantify microstructural tissue components along the tracts in fetal compartments that are critical substrates for white matter maturation, such as the subplate and intermediate zone. We identify patterns of change in the diffusion metrics that reflect critical neurobiological transitions occurring in the second to third trimester, such as the disassembly of radial glial scaffolding and the lamination of the cortical plate. These maturational trajectories of MR signal in transient fetal compartments provide a normative reference to complement histological knowledge, facilitating future studies to establish how developmental disruptions in these regions contribute to pathophysiology.

Multi-readout DWI with a reduced FOV for studying the coupling between diffusion and T 2 * relaxation in the prostate

PURPOSE

To develop a DWI sequence with multiple readout echo-trains in a single shot (multi-readout DWI) over a reduced FOV, and to demonstrate its ability to achieve high data acquisition efficiency in the study of coupling between diffusion and relaxation in the human prostate.

METHODS

The proposed multi-readout DWI sequence plays out multiple EPI readout echo-trains after a Stejskal-Tanner diffusion preparation module. Each EPI readout echo-train corresponded to a distinct effective TE. To maintain a high spatial resolution with a relatively short echo-train for each readout, a 2D RF pulse was used to limit the FOV. Experiments were performed on the prostate of six healthy subjects to acquire a set of images with three b values (0, 500, and 1000 s/mm² ) and three TEs (63.0, 78.8, and 94.6 ms), producing three ADC maps at different TEs and three T 2 * $$ {T}_2^{\ast } $$ maps at different b values.

RESULTS

Multi-readout DWI enabled a threefold acceleration without compromising the spatial resolution when compared with a conventional single-readout sequence. Images with three b values and three TEs were obtained in 3 min 40 s with an adequate SNR (≥ 26.9). The ADC values (1.45 ± 0.13, 1.52 ± 0.14, and 1.58 ± 0.15  μm 2 / ms $$ {\upmu \mathrm{m}}^2/\mathrm{ms} $$ ; P < 0.01) exhibited an increasing trend as TEs increased (63.0 ms, 78.8 ms, and 94.6 ms), whereas T 2 * $$ {T}_2^{\ast } $$ values (74.78 ± 13.21, 63.21 ± 7.84, and 56.61 ± 5.05 ms; P < 0.01) decreases as the b values increased (0, 500, and 1000 s/mm² ).

CONCLUSION

The multi-readout DWI sequence over a reduced FOV provides a time-efficient technique to study the coupling between diffusion and relaxation times.

Multi-modal MRI reveals changes in placental function following preterm premature rupture of membranes

PURPOSE

Preterm premature rupture of membranes complicates up to 40% of premature deliveries. Fetal infection may occur in the absence of maternal symptoms, delaying diagnosis and increasing morbidity and mortality. A noninvasive antenatal assessment of early signs of placental inflammation is therefore urgently required.

METHODS

Sixteen women with preterm premature rupture of membranes < 34 weeks gestation and 60 women with uncomplicated pregnancies were prospectively recruited. A modified diffusion-weighted spin-echo single shot EPI sequence with a diffusion preparation acquiring 264 unique parameter combinations in < 9 min was obtained on a clinical 3 Tesla MRI scanner. The data was fitted to a 2-compartment T 2 * $$ {\mathrm{T}}_2^{\ast } $$ -intravoxel incoherent motion model comprising fast and slowly circulating fluid pools to obtain quantitative information on perfusion, density, and tissue composition. Z values were calculated, and correlation with time from between the rupture of membranes and the scan, gestational age at delivery, and time between scan and delivery assessed.

RESULTS

Placental T 2 * $$ {\mathrm{T}}_2^{\ast } $$ was significantly reduced in preterm premature rupture of membranes, and the 2-compartmental model demonstrated that this decline is mainly linked to the perfusion component observed in the placental parenchyma. Multi-modal MRI measurement of placental function is linked to gestational age at delivery and time from membrane rupture.

CONCLUSION

More complex models and data acquisition can potentially improve fitting of the underlying etiology of preterm birth compared with individual single-contrast models and contribute to additional insights in the future. This will need validation in larger cohorts. A multi-modal MRI acquisition between rupture of the membranes and delivery can be used to measure placental function and is linked to gestational age at delivery.

MR Fingerprinting with b-Tensor Encoding for Simultaneous Quantification of Relaxation and Diffusion in a Single Scan

PURPOSE

Although both relaxation and diffusion imaging are sensitive to tissue microstructure, studies have reported limited sensitivity and robustness of using relaxation or conventional diffusion alone to characterize tissue microstructure. Recently, it has been shown that tensor-valued diffusion encoding and joint relaxation-diffusion quantification enable more reliable quantification of compartment-specific microstructural properties. However, scan times to acquire such data can be prohibitive. Here, we aim to simultaneously quantify relaxation and diffusion using MR fingerprinting (MRF) and b-tensor encoding in a clinically feasible time.

METHODS

We developed multidimensional MRF scans (mdMRF) with linear and spherical b-tensor encoding (LTE and STE) to simultaneously quantify T1, T2, and ADC maps from a single scan. The image quality, accuracy, and scan efficiency were compared between the mdMRF using LTE and STE. Moreover, we investigated the robustness of different sequence designs to signal errors and their impact on the maps.

RESULTS

T1 and T2 maps derived from the mdMRF scans have consistently high image quality, while ADC maps are sensitive to different sequence designs. Notably, the fast imaging steady state precession (FISP)-based mdMRF scan with peripheral pulse gating provides the best ADC maps that are free of image distortion and shading artifacts.

CONCLUSION

We demonstrated the feasibility of quantifying T1, T2, and ADC maps simultaneously from a single mdMRF scan in around 24 s/slice. The map quality and quantitative values are consistent with the reference scans.

Primary Multiparametric Quantitative Brain MRI: State-of-the-Art Relaxometric and Proton Density Mapping Techniques

This review on brain multiparametric quantitative MRI (MP-qMRI) focuses on the primary subset of quantitative MRI (qMRI) parameters that represent the mobile ("free") and bound ("motion-restricted") proton pools. Such primary parameters are the proton densities, relaxation times, and magnetization transfer parameters. Diffusion qMRI is also included because of its wide implementation in complete clinical MP-qMRI application. MP-qMRI advances were reviewed over the past 2 decades, with substantial progress observed toward accelerating image acquisition and increasing mapping accuracy. Areas that need further investigation and refinement are identified as follows: (a) the biologic underpinnings of qMRI parameter values and their changes with age and/or disease and (b) the theoretical limitations implicitly built into most qMRI mapping algorithms that do not distinguish between the different spatial scales of voxels versus spin packets, the central physical object of the Bloch theory. With rapidly improving image processing techniques and continuous advances in computer hardware, MP-qMRI has the potential for implementation in a wide range of clinical applications. Currently, three emerging MP-qMRI applications are synthetic MRI, macrostructural qMRI, and microstructural tissue modeling.

Assessment of the fetal lungs in utero

Antenatal diagnosis of abnormal pulmonary development has improved significantly over recent years because of progress in imaging techniques. Two-dimensional ultrasound is the mainstay of investigation of pulmonary pathology during pregnancy, providing good prognostication in conditions such as congenital diaphragmatic hernia; however, it is less validated in other high-risk groups such as those with congenital pulmonary airway malformation or preterm premature rupture of membranes. Three-dimensional assessment of lung volume and size is now possible using ultrasound or magnetic resonance imaging; however, the use of these techniques is still limited because of unpredictable fetal motion, and such tools have also been inadequately validated in high-risk populations other than those with congenital diaphragmatic hernia. The advent of advanced, functional magnetic resonance imaging techniques such as diffusion and T2* imaging, and the development of postprocessing pipelines that facilitate motion correction, have enabled not only more accurate evaluation of pulmonary size, but also assessment of tissue microstructure and perfusion. In the future, fetal magnetic resonance imaging may have an increasing role in the prognostication of pulmonary abnormalities and in monitoring current and future antenatal therapies to enhance lung development. This review aims to examine the current imaging methods available for assessment of antenatal lung development and to outline possible future directions.

Mapping astrogliosis in the individual human brain using multidimensional MRI

There are currently no non-invasive imaging methods available for astrogliosis assessment or mapping in the central nervous system despite its essential role in the response to many disease states, such as infarcts, neurodegenerative conditions, traumatic brain injury and infection. Multidimensional MRI is an increasingly employed imaging modality that maximizes the amount of encoded chemical and microstructural information by probing relaxation (T1 and T2) and diffusion mechanisms simultaneously. Here, we harness the exquisite sensitivity of this imagining modality to derive a signature of astrogliosis and disentangle it from normative brain at the individual level using machine learning. We investigated ex vivo cerebral cortical tissue specimens derived from seven subjects who sustained blast-induced injuries, which resulted in scar-border forming astrogliosis without being accompanied by other types of neuropathological abnormality, and from seven control brain donors. By performing a combined post-mortem radiology and histopathology correlation study we found that astrogliosis induces microstructural and chemical changes that are robustly detected with multidimensional MRI, and which can be attributed to astrogliosis because no axonal damage, demyelination or tauopathy were histologically observed in any of the cases in the study. Importantly, we showed that no one-dimensional T1, T2 or diffusion MRI measurement can disentangle the microscopic alterations caused by this neuropathology. Based on these findings, we developed a within-subject anomaly detection procedure that generates MRI-based astrogliosis biomarker maps ex vivo, which were significantly and strongly correlated with co-registered histological images of increased glial fibrillary acidic protein deposition (r = 0.856, P < 0.0001; r = 0.789, P < 0.0001; r = 0.793, P < 0.0001, for diffusion-T2, diffusion-T1 and T1-T2 multidimensional data sets, respectively). Our findings elucidate the underpinning of MRI signal response from astrogliosis, and the demonstrated high spatial sensitivity and specificity in detecting reactive astrocytes at the individual level, and if reproduced in vivo, will significantly impact neuroimaging studies of injury, disease, repair and aging, in which astrogliosis has so far been an invisible process radiologically.

Neurite Exchange Imaging (NEXI): A minimal model of diffusion in gray matter with inter-compartment water exchange

Biophysical models of diffusion in white matter have been center-stage over the past two decades and are essentially based on what is now commonly referred to as the "Standard Model" (SM) of non-exchanging anisotropic compartments with Gaussian diffusion. In this work, we focus on diffusion MRI in gray matter, which requires rethinking basic microstructure modeling blocks. In particular, at least three contributions beyond the SM need to be considered for gray matter: water exchange across the cell membrane - between neurites and the extracellular space; non-Gaussian diffusion along neuronal and glial processes - resulting from structural disorder; and signal contribution from soma. For the first contribution, we propose Neurite Exchange Imaging (NEXI) as an extension of the SM of diffusion, which builds on the anisotropic Kärger model of two exchanging compartments. Using datasets acquired at multiple diffusion weightings (b) and diffusion times (t) in the rat brain in vivo, we investigate the suitability of NEXI to describe the diffusion signal in the gray matter, compared to the other two possible contributions. Our results for the diffusion time window 20-45 ms show minimal diffusivity time-dependence and more pronounced kurtosis decay with time, which is well fit by the exchange model. Moreover, we observe lower signal for longer diffusion times at high b. In light of these observations, we identify exchange as the mechanism that best explains these signal signatures in both low-b and high-b regime, and thereby propose NEXI as the minimal model for gray matter microstructure mapping. We finally highlight multi-b multi-t acquisition protocols as being best suited to estimate NEXI model parameters reliably. Using this approach, we estimate the inter-compartment water exchange time to be 15 - 60 ms in the rat cortex and hippocampus in vivo, which is of the same order or shorter than the diffusion time in typical diffusion MRI acquisitions. This suggests water exchange as an essential component for interpreting diffusion MRI measurements in gray matter.

Diffusion-relaxation scattered MR signal representation in a multi-parametric sequence

This work focuses on obtaining a more general diffusion magnetic resonance imaging (MRI) signal representation that accounts for a longitudinal T₁ and transverse T₂^⋆ relaxations while at the same time integrating directional diffusion in the context of scattered multi-parametric acquisitions, where only a few diffusion gradient directions and b-values are available for each pair of echo and inversion times. The method is based on the three-dimensional simple harmonic oscillator-based reconstruction and estimation (SHORE) representation of the diffusion signal, which enables the estimation of the orientation distribution function and the retrieval of various quantitative indices such as the generalized fractional anisotropy or the return-to-the-origin probability while simultaneously resolving for T₁ and T₂^⋆ relaxation times. Our technique, the Relax-SHORE, has been tested on both in silico and in vivo diffusion-relaxation scattered MR data. The results show that Relax-SHORE is accurate in the context of scattered acquisitions while guaranteeing flexibility in the diffusion signal representation from multi-parametric sequences.

What's new and what's next in diffusion MRI preprocessing

Diffusion MRI (dMRI) provides invaluable information for the study of tissue microstructure and brain connectivity, but suffers from a range of imaging artifacts that greatly challenge the analysis of results and their interpretability if not appropriately accounted for. This review will cover dMRI artifacts and preprocessing steps, some of which have not typically been considered in existing pipelines or reviews, or have only gained attention in recent years: brain/skull extraction, B-matrix incompatibilities w.r.t the imaging data, signal drift, Gibbs ringing, noise distribution bias, denoising, between- and within-volumes motion, eddy currents, outliers, susceptibility distortions, EPI Nyquist ghosts, gradient deviations, B1 bias fields, and spatial normalization. The focus will be on "what's new" since the notable advances prior to and brought by the Human Connectome Project (HCP), as presented in the predecessing issue on "Mapping the Connectome" in 2013. In addition to the development of novel strategies for dMRI preprocessing, exciting progress has been made in the availability of open source tools and reproducible pipelines, databases and simulation tools for the evaluation of preprocessing steps, and automated quality control frameworks, amongst others. Finally, this review will consider practical considerations and our view on "what's next" in dMRI preprocessing.

Quantitative mapping of the brain's structural connectivity using diffusion MRI tractography: A review

Diffusion magnetic resonance imaging (dMRI) tractography is an advanced imaging technique that enables in vivo reconstruction of the brain's white matter connections at macro scale. It provides an important tool for quantitative mapping of the brain's structural connectivity using measures of connectivity or tissue microstructure. Over the last two decades, the study of brain connectivity using dMRI tractography has played a prominent role in the neuroimaging research landscape. In this paper, we provide a high-level overview of how tractography is used to enable quantitative analysis of the brain's structural connectivity in health and disease. We focus on two types of quantitative analyses of tractography, including: 1) tract-specific analysis that refers to research that is typically hypothesis-driven and studies particular anatomical fiber tracts, and 2) connectome-based analysis that refers to research that is more data-driven and generally studies the structural connectivity of the entire brain. We first provide a review of methodology involved in three main processing steps that are common across most approaches for quantitative analysis of tractography, including methods for tractography correction, segmentation and quantification. For each step, we aim to describe methodological choices, their popularity, and potential pros and cons. We then review studies that have used quantitative tractography approaches to study the brain's white matter, focusing on applications in neurodevelopment, aging, neurological disorders, mental disorders, and neurosurgery. We conclude that, while there have been considerable advancements in methodological technologies and breadth of applications, there nevertheless remains no consensus about the "best" methodology in quantitative analysis of tractography, and researchers should remain cautious when interpreting results in research and clinical applications.

Accelerating joint relaxation-diffusion MRI by integrating time division multiplexing and simultaneous multi-slice (TDM-SMS) strategies

PURPOSE

To accelerate the acquisition of relaxation-diffusion imaging by integrating time-division multiplexing (TDM) with simultaneous multi-slice (SMS) for EPI and evaluate imaging quality and diffusion measures.

METHODS

The time-division multiplexing (TDM) technique and SMS method were integrated to achieve a high slice-acceleration (e.g., 6×) factor for acquiring relaxation-diffusion MRI. Two variants of the sequence, referred to as TDM3e-SMS and TDM2s-SMS, were developed to simultaneously acquire slice groups with three distinct TEs and two slice groups with the same TE, respectively. Both sequences were evaluated on a 3T scanner with in vivo human brains and compared with standard single-band (SB) -EPI and SMS-EPI using diffusion measures and tractography results.

RESULTS

Experimental results showed that the TDM3e-SMS sequence with total slice acceleration of 6 (multiplexing factor (MP) = 3 × multi-band factor (MB) = 2) provided similar image intensity and microstructure measures compared to standard SMS-EPI with MB = 2, and yielded less bias in intensity compared to standard SMS-EPI with MB = 4. The three sequences showed a similar positive correlation between TE and mean kurtosis (MK) and a negative correlation between TE and mean diffusivity (MD) in white matter. Multi-fiber tractography also shows consistency of results in TE-dependent measures between different sequences. The TDM2s-SMS sequence (MP = 2, MB = 2) also provided imaging measures similar to standard SMS-EPI sequences (MB = 2) for single-TE diffusion imaging.

CONCLUSIONS

The TDM-SMS sequence can provide additional 2× to 3× acceleration to SMS without degrading imaging quality. With the significant reduction in scan time, TDM-SMS makes joint relaxation-diffusion MRI a feasible technique in neuroimaging research to investigate new markers of brain disorders.

Distortion-Free Diffusion Imaging Using Self-Navigated Cartesian Echo-Planar Time Resolved Acquisition and Joint Magnitude and Phase Constrained Reconstruction

Echo-planar time resolved imaging (EPTI) is an effective approach for acquiring high-quality distortion-free images with a multi-shot EPI (ms-EPI) readout. As with traditional ms-EPI acquisitions, inter-shot phase variations present a main challenge when incorporating EPTI into a diffusion-prepared pulse sequence. The aim of this study is to develop a self-navigated Cartesian EPTI-based (scEPTI) acquisition together with a magnitude and phase constrained reconstruction for distortion-free diffusion imaging. A self-navigated Cartesian EPTI-based diffusion-prepared pulse sequence is designed. The different phase components in EPTI diffusion signal are analyzed and an approach to synthesize a fully phase-matched navigator for the inter-shot phase correction is demonstrated. Lastly, EPTI contains richer magnitude and phase information than conventional ms-EPI, such as the magnitude and phase correlations along the temporal dimension. The potential of these magnitude and phase correlations to enhance the reconstruction is explored. The reconstruction results with and without phase matching and with and without phase or magnitude constraints are compared. Compared with reconstruction without phase matching, the proposed phase matching method can improve the accuracy of inter-shot phase correction and reduce signal corruption in the final diffusion images. Magnitude constraints further improve image quality by suppressing the background noise and thereby increasing SNR, while phase constraints can mitigate possible image blurring from adding magnitude constraints. The high-quality distortion-free diffusion images and simultaneous diffusion-relaxometry imaging capacity provided by the proposed EPTI design represent a highly valuable tool for both clinical and neuroscientific assessments of tissue microstructure.

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.

Quantifying myelin in crossing fibers using diffusion-prepared phase imaging: Theory and simulations

PURPOSE

Myelin has long been the target of neuroimaging research. However, most available techniques can only provide a voxel-averaged estimate of myelin content. In the human brain, white matter fiber pathways connecting different brain areas and carrying different functions often cross each other in the same voxel. A measure that can differentiate the degree of myelination of crossing fibers would provide a more specific marker of myelination.

THEORY AND METHODS

One MRI signal property that is sensitive to myelin is the phase accumulation. This sensitivity is used by measuring the phase accumulation of the signal remaining after diffusion-weighting, which is called diffusion-prepared phase imaging (DIPPI). Including diffusion-weighting before estimating the phase accumulation has two distinct advantages for estimating the degree of myelination: (1) It increases the relative contribution of intra-axonal water, whose phase is related linearly to the thickness of the surrounding myelin (in particular the log g-ratio); and (2) it gives directional information, which can be used to distinguish between crossing fibers. Here the DIPPI sequence is described, an approach is proposed to estimate the log g-ratio, and simulations are used and DIPPI data acquired in an isotropic phantom to quantify other sources of phase accumulation.

RESULTS

The expected bias is estimated in the log g-ratio for reasonable in vivo acquisition parameters caused by eddy currents (~4%-10%), remaining extra-axonal signal (~15%), and gradients in the bulk off-resonance field (<10% for most of the brain).

CONCLUSION

This new sequence may provide a g-ratio estimate per fiber population crossing within a voxel.

Model-based multi-parameter mapping

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

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.

Measuring compartmental T2-orientational dependence in human brain white matter using a tiltable RF coil and diffusion-T2 correlation MRI

The anisotropy of brain white matter microstructure manifests itself in orientational-dependence of various MRI contrasts, and can result in significant quantification biases if ignored. Understanding the origins of this orientation-dependence could enhance the interpretation of MRI signal changes in development, ageing and disease and ultimately improve clinical diagnosis. Using a novel experimental setup, this work studies the contributions of the intra- and extra-axonal water to the orientation-dependence of one of the most clinically-studied parameters, apparent transverse relaxation T2. Specifically, a tiltable receive coil is interfaced with an ultra-strong gradient MRI scanner to acquire multidimensional MRI data with an unprecedented range of acquisition parameters. Using this setup, compartmental T2 can be disentangled based on differences in diffusional-anisotropy, and its orientation-dependence further elucidated by re-orienting the head with respect to the main magnetic field B→0. A dependence of (compartmental) T2 on the fibre orientation w.r.t. B→0 was observed, and further quantified using characteristic representations for susceptibility- and magic angle effects. Across white matter, anisotropy effects were dominated by the extra-axonal water signal, while the intra-axonal water signal decay varied less with fibre-orientation. Moreover, the results suggest that the stronger extra-axonal T2 orientation-dependence is dominated by magnetic susceptibility effects (presumably from the myelin sheath) while the weaker intra-axonal T2 orientation-dependence may be driven by a combination of microstructural effects. Even though the current design of the tiltable coil only offers a modest range of angles, the results demonstrate an overall effect of tilt and serve as a proof-of-concept motivating further hardware development to facilitate experiments that explore orientational anisotropy. These observations have the potential to lead to white matter microstructural models with increased compartmental sensitivity to disease, and can have direct consequences for longitudinal and group-wise T2- and diffusion-MRI data analysis, where the effect of head-orientation in the scanner is commonly ignored.

An efficient and combined placental T 1 -ADC acquisition in pregnancies with and without pre-eclampsia

PURPOSE

To provide a new approach to jointly assess microstructural and molecular properties of the human placenta in vivo fast and efficiently and to present initial evidence in cohorts of healthy pregnancies and those affected by pre-eclampsia.

METHODS

Slice and diffusion preparation shuffling, built on the previously proposed ZEBRA method, is presented as a robust and fast way to obtain T 1 and apparent diffusivity coefficient (ADC) values. Joint modeling and evaluation is performed on a cohort of healthy and pre-eclamptic participants at 3T.

RESULTS

The datasets show the ability to obtain robust and fast T 1 -ADC measurements. Significant decay over gestation in T 1 (-11 ms/week, P < . 05 ) and a trend toward significance in ADC (-0.23 mm/ s 2 /week, P = .08) values can be observed in a control cohort. Values for the pre-eclamptic pregnancies show a negative trend for both ADC and T 1 .

CONCLUSIONS

The presented sequence allows the simultaneous acquisition of 2 of the most promising quantitative parameters to study placental insufficiency-identified individually as relevant in previous studies-in under 2 minutes. This allows dynamic assessment of physiological processes, reduced inconsistency in spatial comparisons due to reduced motion artefacts and opens novel avenues for analysis. Initial results in pre-eclamptic placentas, with depicted changes in both ADC and T 1 , illustrate its potential to identify cases of placental insufficiency. Future work will focus on expanding the field-of-view using multi-band acceleration techniques and the expansion to larger and more diverse patient groups.

Accelerated diffusion and relaxation-diffusion MRI using time-division multiplexing EPI

PURPOSE

To develop a time-division multiplexing echo-planar imaging (TDM-EPI) sequence for approximately two- to threefold acceleration when acquiring joint relaxation-diffusion MRI data with multiple TEs.

METHODS

The proposed TDM-EPI sequence interleaves excitation and data collection for up to 3 separate slices at different TEs and uses echo-shifting gradients to disentangle the overlapping echo signals during the readout period. By properly arranging the sequence event blocks for each slice and adjusting the echo-shifting gradients, diffusion-weighted images from separate slices can be acquired. Therefore, we present 2 variants of the sequence. A single-TE TDM-EPI is presented to demonstrate the concept. Next, a multi-TE TDM-EPI is presented to highlight the advantages of the TDM approach for relaxation-diffusion imaging. These sequences were evaluated on a 3 Tesla scanner with a water phantom and in vivo human brain data.

RESULTS

The single-TE TDM-EPI sequence can simultaneously acquire 2 slices with a maximum b value of 3000 s/mm² and 2.5 mm isotropic resolution using interleaved readout windows with TE ≈ 78 ms. With the same b value and resolution, the multi-TE TDM-EPI sequence can simultaneously acquire 2 or 3 separate slices using interleaved readout sections with shortest TE ≈ 70 ms and ΔTE ≈ 30 ms. Phantom and in vivo experiments have shown that the proposed TDM-EPI sequences can provide similar image quality and diffusion measures as conventional EPI readouts with multiple echoes but can reduce the overall relaxation-diffusion protocol scan time by approximately two- to threefold.

CONCLUSION

TDM-EPI is a novel approach to acquire diffusion imaging data at multiple TEs. This enables a significant reduction in acquisition time for relaxation-diffusion MRI experiments but without compromising image quality and diffusion measurements, thus removing a significant barrier to the adoption of relaxation-diffusion MRI in clinical research studies of neurological and mental disorders.

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

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

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.

Three-dimensional self-attention conditional GAN with spectral normalization for multimodal neuroimaging synthesis

PURPOSE

To develop a new 3D generative adversarial network that is designed and optimized for the application of multimodal 3D neuroimaging synthesis.

METHODS

We present a 3D conditional generative adversarial network (GAN) that uses spectral normalization and feature matching to stabilize the training process and ensure optimization convergence (called SC-GAN). A self-attention module was also added to model the relationships between widely separated image voxels. The performance of the network was evaluated on the data set from ADNI-3, in which the proposed network was used to predict PET images, fractional anisotropy, and mean diffusivity maps from multimodal MRI. Then, SC-GAN was applied on a multidimensional diffusion MRI experiment for superresolution application. Experiment results were evaluated by normalized RMS error, peak SNR, and structural similarity.

RESULTS

In general, SC-GAN outperformed other state-of-the-art GAN networks including 3D conditional GAN in all three tasks across all evaluation metrics. Prediction error of the SC-GAN was 18%, 24% and 29% lower compared to 2D conditional GAN for fractional anisotropy, PET and mean diffusivity tasks, respectively. The ablation experiment showed that the major contributors to the improved performance of SC-GAN are the adversarial learning and the self-attention module, followed by the spectral normalization module. In the superresolution multidimensional diffusion experiment, SC-GAN provided superior predication in comparison to 3D Unet and 3D conditional GAN.

CONCLUSION

In this work, an efficient end-to-end framework for multimodal 3D medical image synthesis (SC-GAN) is presented. The source code is also made available at https://github.com/Haoyulance/SC-GAN.

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.

Data-Driven multi-Contrast spectral microstructure imaging with InSpect: INtegrated SPECTral component estimation and mapping

•

Unsupervised learning technique for spectroscopic analysis of quantitative MRI.

•

Shares information across voxels to improve estimation of multi-dimensional or single-dimensional spectra.

•

Spectral maps are dramatically improved compared to existing approaches.

•

Can potentially identify and map tissue environments; in placental diffusion-relaxometry MRI we demonstrate that it identifies components that correspond to distinct tissue types.

We introduce and demonstrate an unsupervised machine learning technique for spectroscopic analysis of quantitative MRI experiments. Our algorithm supports estimation of one-dimensional spectra from single-contrast data, and multidimensional correlation spectra from simultaneous multi-contrast data. These spectrum-based approaches allow model-free investigation of tissue properties, but require regularised inversion of a Laplace transform or Fredholm integral, which is an ill-posed calculation. Here we present a method that addresses this limitation in a data-driven way. The algorithm simultaneously estimates a canonical basis of spectral components and voxelwise maps of their weightings, thereby pooling information across whole images to regularise the ill-posed problem. We show in simulations that our algorithm substantially outperforms current voxelwise spectral approaches. We demonstrate the method on multi-contrast diffusion-relaxometry placental MRI scans, revealing anatomically-relevant sub-structures, and identifying dysfunctional placentas. Our algorithm vastly reduces the data required to reliably estimate spectra, opening up the possibility of quantitative MRI spectroscopy in a wide range of new applications. Our InSpect code is available at github.com/paddyslator/inspect.

Efficient whole-brain tract-specific T1 mapping at 3T with slice-shuffled inversion-recovery diffusion-weighted imaging

PURPOSE

Most voxels in white matter contain multiple fiber populations with different orientations and levels of myelination. Conventional T₁ mapping measures 1 T₁ value per voxel, representing a weighted average of the multiple tract T₁ times. Inversion-recovery diffusion-weighted imaging (IR-DWI) allows the T₁ times of multiple tracts in a voxel to be disentangled, but the scan time is prohibitively long. Recently, slice-shuffled IR-DWI implementations have been proposed to significantly reduce scan time. In this work, we demonstrate that we can measure tract-specific T₁ values in the whole brain using simultaneous multi-slice slice-shuffled IR-DWI at 3T.

METHODS

We perform simulations to evaluate the accuracy and precision of our crossing fiber IR-DWI signal model for various fiber parameters. The proposed sequence and signal model are tested in a phantom consisting of crossing asparagus pieces doped with gadolinium to vary T₁ , and in 2 human subjects.

RESULTS

Our simulations show that tract-specific T₁ times can be estimated within 5% of the nominal fiber T₁ values. Tract-specific T₁ values were resolved in subvoxel 2 fiber crossings in the asparagus phantom. Tract-specific T₁ times were resolved in 2 different tract crossings in the human brain where myelination differences have previously been reported; the crossing of the cingulum and genu of the corpus callosum and the crossing of the corticospinal tract and pontine fibers.

CONCLUSION

Whole-brain tract-specific T₁ mapping is feasible using slice-shuffled IR-DWI at 3T. This technique has the potential to improve the microstructural characterization of specific tracts implicated in neurodevelopment, aging, and demyelinating disorders.

Resolving bundle-specific intra-axonal T2 values within a voxel using diffusion-relaxation tract-based estimation

At the typical spatial resolution of MRI in the human brain, approximately 60-90% of voxels contain multiple fiber populations. Quantifying microstructural properties of distinct fiber populations within a voxel is therefore challenging but necessary. While progress has been made for diffusion and T1-relaxation properties, how to resolve intra-voxel T2 heterogeneity remains an open question. Here a novel framework, named COMMIT-T2, is proposed that uses tractography-based spatial regularization with diffusion-relaxometry data to estimate multiple intra-axonal T2 values within a voxel. Unlike previously-proposed voxel-based T2 estimation methods, which (when applied in white matter) implicitly assume just one fiber bundle in the voxel or the same T2 for all bundles in the voxel, COMMIT-T2 can recover specific T2 values for each unique fiber population passing through the voxel. In this approach, the number of recovered unique T2 values is not determined by a number of model parameters set a priori, but rather by the number of tractography-reconstructed streamlines passing through the voxel. Proof-of-concept is provided in silico and in vivo, including a demonstration that distinct tract-specific T2 profiles can be recovered even in the three-way crossing of the corpus callosum, arcuate fasciculus, and corticospinal tract. We demonstrate the favourable performance of COMMIT-T2 compared to that of voxelwise approaches for mapping intra-axonal T2 exploiting diffusion, including a direction-averaged method and AMICO-T2, a new extension to the previously-proposed Accelerated Microstructure Imaging via Convex Optimization (AMICO) framework.

Placental magnetic resonance imaging in chronic hypertension: A case-control study

INTRODUCTION

We aimed to explore the use of magnetic resonance imaging (MRI) in vivo as a tool to elucidate the placental phenotype in women with chronic hypertension.

METHODS

In case-control study, women with chronic hypertension and those with uncomplicated pregnancies were imaged using either a 3T Achieva or 1.5T Ingenia scanner. T2-weighted images, diffusion weighted and T1/T2* relaxometry data was acquired. Placental T2*, T1 and apparent diffusion coefficient (ADC) maps were calculated.

RESULTS

129 women (43 with chronic hypertension and 86 uncomplicated pregnancies) were imaged at a median of 27.7 weeks' gestation (interquartile range (IQR) 23.9-32.1) and 28.9 (IQR 26.1-32.9) respectively. Visual analysis of T2-weighted imaging demonstrated placentae to be either appropriate for gestation or to have advanced lobulation in women with chronic hypertension, resulting in a greater range of placental mean T2* values for a given gestation, compared to gestation-matched controls. Both skew and kurtosis (derived from histograms of T2* values across the whole placenta) increased with advancing gestational age at imaging in healthy pregnancies; women with chronic hypertension had values overlapping those in the control group range. Upon visual assessment, the mean ADC declined in the third trimester, with a corresponding decline in placental mean T2* values and showed an overlap of values between women with chronic hypertension and the control group.

DISCUSSION

A combined placental MR examination including T2 weighted imaging, T2*, T1 mapping and diffusion imaging demonstrates varying placental phenotypes in a cohort of women with chronic hypertension, showing overlap with the control group.

Toward Quantification: Microstructure and Magnetic Resonance Fingerprinting

Quantitative magnetic resonance imaging (MRI) is a long-standing challenge. We advocate that the origin of the problem is the simplification applied in commonly used models of the MRI signal relation to the target parameters of biological tissues. Two research fields are briefly reviewed as ways to respond to the challenge of quantitative MRI, both experiencing an exponential growth right now. Microstructure MRI strives to build physiology-based models from cells to signal and, given the signal, back to the cells again. Magnetic resonance fingerprinting aims at efficient simultaneous determination of multiple signal parameters. The synergy of these yet disjoined approaches promises truly quantitative MRI with specific target-oriented diagnostic tools rather than universal imaging methods.

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.