Diffusion-relaxation correlation spectroscopic imaging: A multidimensional approach for probing microstructure

Data Processing in Multidimensional MRI For Biomarker Identification: Is It Necessary?

Multidimensional MRI (MD-MRI) is an emerging technique that holds promise for identifying tissue characteristics that could be indicative of pathologies. Before these characteristics can be interpreted, MD-MRI measurements are converted into an spectrum. These spectra are then utilized to obtain some understanding of the underlying tissue microstructure, often through the use of statistical, machine learning, and mathematical modeling methods. The aim of this study was to compare outcomes of using unprocessed MDMRI signals for statistical regression in comparison to the corresponding spectra. Backed by a theoretical argument, we described an experimental procedure regressing both MDMRI signals and spectra to histological outcomes intrasubject. Through using multiple conventional ML methods, and a proposed method using convex sets, we aimed to see which yielded the highest accuracy. Both theory and experimental evidence suggest that, without a priori information, statistical regression was best performed on the MDMRI signal. We conclude, barring any a priori information regarding tissue changes, there is no significant advantage to performing regression analysis on reconstructed spectra in the process of biomarker identification.

Ex vivo massively multidimensional diffusion-relaxation correlation MRI: scan-rescan reproducibility and caveats

Massively multidimensional diffusion-relaxation correlation MRI (MMD-MRI) provides information beyond the traditional voxel-averaged metric that may better characterize the microstructural characteristics of biological tissues. MMD-MRI reproducibility has been established in clinical settings, but has yet to be thoroughly evaluated under preclinical conditions, where superior hardware and modulated gradient waveforms enhance its performance. In this study, we investigate the reproducibility of MMD-MRI on a micro-imaging system using ex vivo mouse brains. Notably, the estimated signal fractions of intra-voxel spectral components in the MD-MRI distribution, corresponding to white and gray matter, along with the frequency-dependent parameters, demonstrated high reproducibility. We identified bias between scan and rescan in some of the metrics, which we attribute to the time gap between repeated scans pointing to a long-time progressive fixation effect. We compare our results with in vivo results from clinical scanners and show the reproducibility of diffusion frequency-dependent metrics to benefit from the improved gradient hardware on our preclinical setup. Our results inform future micro-imaging ex vivo MMD-MRI studies of the reproducibility of MMD-MRI metrics and their dependence on fixation time.

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.

Advanced Diffusion-Weighted MRI for Cancer Microstructure Assessment in Body Imaging, and Its Relationship With Histology

Diffusion-weighted magnetic resonance imaging (DW-MRI) aims to disentangle multiple biological signal sources in each imaging voxel, enabling the computation of innovative maps of tissue microstructure. DW-MRI model development has been dominated by brain applications. More recently, advanced methods with high fidelity to histology are gaining momentum in other contexts, for example, in oncological applications of body imaging, where new biomarkers are urgently needed. The objective of this article is to review the state-of-the-art of DW-MRI in body imaging (ie, not including the nervous system) in oncology, and to analyze its value as compared to reference colocalized histology measurements, given that demonstrating the histological validity of any new DW-MRI method is essential. In this article, we review the current landscape of DW-MRI techniques that extend standard apparent diffusion coefficient (ADC), describing their acquisition protocols, signal models, fitting settings, microstructural parameters, and relationship with histology. Preclinical, clinical, and in/ex vivo studies were included. The most used techniques were intravoxel incoherent motion (IVIM; 36.3% of used techniques), diffusion kurtosis imaging (DKI; 16.7%), vascular, extracellular, and restricted diffusion for cytometry in tumors (VERDICT; 13.3%), and imaging microstructural parameters using limited spectrally edited diffusion (IMPULSED; 11.7%). Another notable category of techniques relates to innovative b-tensor diffusion encoding or joint diffusion-relaxometry. The reviewed approaches provide histologically meaningful indices of cancer microstructure (eg, vascularization/cellularity) which, while not necessarily accurate numerically, may still provide useful sensitivity to microscopic pathological processes. Future work of the community should focus on improving the inter-/intra-scanner robustness, and on assessing histological validity in broader contexts. LEVEL OF EVIDENCE: NA TECHNICAL EFFICACY: Stage 2.

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.

Reproducibility of spatial penalty-based methodologies for intravoxel incoherent motion analysis with diffusion MRI

Objective was to assess the precision and reproducibility of spatial penalty-based intravoxel incoherent motion (IVIM) methods in comparison to the conventional bi-exponential (BE) model-based IVIM methods. IVIM-MRI (11 b-values; 0-800 s/mm2) of forty patients (N = 40; Age = 17.7 ± 5.9 years; Male:Female = 30:10) with biopsy-proven osteosarcoma were acquired on a 1.5 Tesla scanner at 3 time-points: (i) baseline, (ii) after 1-cycle and (iii) after 3-cycles of neoadjuvant chemotherapy. Diffusion coefficient (D), Perfusion coefficient (D*) and Perfusion fraction (f) were estimated at three time-points in whole tumor and healthy muscle tissue using five methodologies (1) BE with three-parameter-fitting (BE), (2) Segmented-BE with two-parameter-fitting (BESeg-2), (3) Segmented-BE with one-parameter-fitting (BESeg-1), (4) BE with adaptive Total-Variation-penalty (BE + TV) and (5) BE with adaptive Huber-penalty (BE + HPF). Within-subject coefficient-of-variation (wCV) and between-subject coefficient-of-variation (bCV) of IVIM parameters were measured in healthy and tumor tissue. For precision and reproducibility, intra-scan comparison of wCV and bCV among five IVIM methods were performed using Friedman test followed by Wilcoxon-signed-ranks (WSR) post-hoc test. Experimental results demonstrated that BE + TV and BE + HPF showed significantly (p < 10-3) lower wCV and bCV for D (wCV: 24-32%; bCV: 22-31%) than BE method (wCV: 38-49%; bCV: 36-46%) across three time-points in healthy muscle and tumor. BE + TV and BE + HPF also demonstrated significantly (p < 10-3) lower wCV and bCV for estimating D* (wCV: 89-108%; bCV: 83-102%) and f (wCV: 55-60%; bCV: 56-60%) than BE, BESeg-2 and BESeg-1 methods (D*-wCV: 102-122%; D*-bCV: 98-114% and f-wCV: 96-130%; f-bCV: 94-125%) in both tumor and healthy tissue across three time-points. Spatial penalty based IVIM analysis methods BE + TV and BE + HPF demonstrated lower variability and improved precision and reproducibility in the current clinical settings.

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.

To characterize small renal cell carcinoma using diffusion relaxation correlation spectroscopic imaging and apparent diffusion coefficient based histogram analysis: a preliminary study

PURPOSE

To study the capability of diffusion-relaxation correlation spectroscopic imaging (DR-CSI) on subtype classification and grade differentiation for small renal cell carcinoma (RCC). Histogram analysis for apparent diffusion coefficient (ADC) was studied for comparison.

MATERIALS AND METHODS

A total of 61 patients with small RCC (< 4 cm) were included in the retrospective study. MRI data were reviewed, including a multi-b (0-1500 s/mm2) multi-TE (51-200 ms) diffusion weighted imaging (DWI) sequence. Region of interest (ROI) was delineated manually on DWI to include solid tumor. For each patient, a D-T2 spectrum was fitted and segmented into 5 compartments, and the volume fractions VA, VB, VC, VD, VE were obtained. ADC mapping was calculated, and histogram parameters ADC 90th, 10th, median, standard deviation, skewness and kurtosis were obtained. All MRI metrices were compared between clear cell RCC (ccRCC) and non-ccRCC group, and between high-grade and low-grade group. Receiver operator curve analysis was used to assess the corresponding diagnostic performance.

RESULTS

Significantly higher ADC 90th, ADC 10th and ADC median, and significantly lower DR-CSI VB was found for ccRCC compared to non-ccRCC. Significantly lower ADC 90th, ADC median and significantly higher VB was found for high-grade RCC compared to low-grade. For identifying ccRCC from non-ccRCC, VB showed the highest area under curve (AUC, 0.861) and specificity (0.882). For differentiating high- from low-grade, ADC 90th showed the highest AUC (0.726) and specificity (0.786), while VB also displayed a moderate AUC (0.715).

CONCLUSION

DR-CSI may offer improved accuracy in subtype identification for small RCC, while do not show better performance for small RCC grading compared to ADC histogram.

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.

Characterization of parotid gland tumors using diffusion-relaxation correlation spectrum imaging: a preliminary study

AIM

To assess the performance of diffusion-relaxation correlation spectrum imaging (DR-CSI) in the characterization of parotid gland tumors.

MATERIALS AND METHODS

Twenty-five pleomorphic adenomas (PA) patients, 9 Warthin's tumors (WT) patients and 7 malignant tumors (MT) patients were prospectively recruited. DR-CSI (7 b-values combined with 5 TEs, totally 35 diffusion-weighted images) was scanned for pre-treatment assessment. Diffusion (D)-T2 signal spectrum summating all voxels were built for each patient, characterized by D-axis with range 0∼5 × 10-3 mm2/s, and T2-axis with range 0∼300ms. With boundaries of 0.5 and 2.5 × 10-3 mm2/s for D, all spectra were divided into three compartments labeled A (low D), B (mediate D) and C (high D). Volume fractions acquired from each compartment (VA, VB, VC) were compared among PA, WT and MT. Diagnostic performance was assessed using receiver operating characteristic analysis and area under the curve (AUC).

RESULTS

Each subtype of parotid tumors had their specific D-T2 spectrum. PA showed significantly lower VA (8.85 ± 4.77% vs 20.68 ± 10.85%), higher VB (63.40 ± 8.18% vs 43.05 ± 7.16%), and lower VC (27.75 ± 8.51% vs 36.27 ± 11.09) than WT (all p<0.05). VB showed optimal diagnostic performance (AUC 0.969, sensitivity 92.00%, specificity 100.00%). MT showed significantly higher VA (21.23 ± 12.36%), lower VB (37.09 ± 6.43%), and higher VC (41.68 ± 13.72%) than PA (all p<0.05). Similarly, VB showed optimal diagnostic performance (AUC 0.994, sensitivity 96.00%, specificity 100.00%). No significant difference of VA, VB and VC was found between WT and MT.

CONCLUSIONS

DR-CSI might be a promising and non-invasive way for characterizing parotid gland tumors.

Resiliency to Alzheimer's disease neuropathology can be distinguished from dementia using cortical astrogliosis imaging

Despite the presence of significant Alzheimer's disease (AD) pathology, characterized by amyloid β (Aβ) plaques and phosphorylated tau (pTau) tangles, some cognitively normal elderly individuals do not inevitably develop dementia. These findings give rise to the notion of cognitive 'resilience', suggesting maintained cognitive function despite the presence of AD neuropathology, highlighting the influence of factors beyond classical pathology. Cortical astroglial inflammation, a ubiquitous feature of symptomatic AD, shows a strong correlation with cognitive impairment severity, potentially contributing to the diversity of clinical presentations. However, noninvasively imaging neuroinflammation, particularly astrogliosis, using MRI remains a significant challenge. Here we sought to address this challenge and to leverage multidimensional (MD) MRI, a powerful approach that combines relaxation with diffusion MR contrasts, to map cortical astrogliosis in the human brain by accessing sub-voxel information. Our goal was to test whether MD-MRI can map astroglial pathology in the cerebral cortex, and if so, whether it can distinguish cognitive resiliency from dementia in the presence of hallmark AD neuropathological changes. We adopted a multimodal approach by integrating histological and MRI analyses using human postmortem brain samples. Ex vivo cerebral cortical tissue specimens derived from three groups comprised of non-demented individuals with significant AD pathology postmortem, individuals with both AD pathology and dementia, and non-demented individuals with minimal AD pathology postmortem as controls, underwent MRI at 7 T. We acquired and processed MD-MRI, diffusion tensor, and quantitative T 1 and T 2 MRI data, followed by histopathological processing on slices from the same tissue. By carefully co-registering MRI and microscopy data, we performed quantitative multimodal analyses, leveraging targeted immunostaining to assess MD-MRI sensitivity and specificity towards Aβ, pTau, and glial fibrillary acidic protein (GFAP), a marker for astrogliosis. Our findings reveal a distinct MD-MRI signature of cortical astrogliosis, enabling the creation of predictive maps for cognitive resilience amid AD neuropathological changes. Multiple linear regression linked histological values to MRI changes, revealing that the MD-MRI cortical astrogliosis biomarker was significantly associated with GFAP burden (standardized β=0.658, pFDR<0.0001), but not with Aβ (standardized β=0.009, p FDR =0.913) or pTau (standardized β=-0.196, p FDR =0.051). Conversely, none of the conventional MRI parameters showed significant associations with GFAP burden in the cortex. While the extent to which pathological glial activation contributes to neuronal damage and cognitive impairment in AD is uncertain, developing a noninvasive imaging method to see its affects holds promise from a mechanistic perspective and as a potential predictor of cognitive outcomes.

In vivo disentanglement of diffusion frequency-dependence, tensor shape, and relaxation using multidimensional MRI

Diffusion MRI with free gradient waveforms, combined with simultaneous relaxation encoding, referred to as multidimensional MRI (MD-MRI), offers microstructural specificity in complex biological tissue. This approach delivers intravoxel information about the microstructure, local chemical composition, and importantly, how these properties are coupled within heterogeneous tissue containing multiple microenvironments. Recent theoretical advances incorporated diffusion time dependency and integrated MD-MRI with concepts from oscillating gradients. This framework probes the diffusion frequency,ω $$ \omega $$ , in addition to the diffusion tensor,D $$ \mathbf{D} $$ , and relaxation,R 1 $$ {R}_1 $$ ,R 2 $$ {R}_2 $$ , correlations. AD ω - R 1 - R 2 $$ \mathbf{D}\left(\omega \right)-{R}_1-{R}_2 $$ clinical imaging protocol was then introduced, with limited brain coverage and 3 mm3 voxel size, which hinder brain segmentation and future cohort studies. In this study, we introduce an efficient, sparse in vivo MD-MRI acquisition protocol providing whole brain coverage at 2 mm3 voxel size. We demonstrate its feasibility and robustness using a well-defined phantom and repeated scans of five healthy individuals. Additionally, we test different denoising strategies to address the sparse nature of this protocol, and show that efficient MD-MRI encoding design demands a nuanced denoising approach. The MD-MRI framework provides rich information that allows resolving the diffusion frequency dependence into intravoxel components based on theirD ω - R 1 - R 2 $$ \mathbf{D}\left(\omega \right)-{R}_1-{R}_2 $$ distribution, enabling the creation of microstructure-specific maps in the human brain. Our results encourage the broader adoption and use of this new imaging approach for characterizing healthy and pathological tissues.

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.

Grading Clear Cell Renal Cell Carcinoma Grade Using Diffusion Relaxation Correlated MR Spectroscopic Imaging

BACKGROUND

Clear cell renal cell carcinoma (ccRCC) is the most common subtype of RCC, and accurate grading is crucial for prognosis and treatment selection. Biopsy is the reference standard for grading, but MRI methods can improve and complement the grading procedure.

PURPOSE

Assess the performance of diffusion relaxation correlation spectroscopic imaging (DR-CSI) in grading ccRCC.

STUDY TYPE

Prospective.

SUBJECTS

79 patients (age: 58.1 +/- 11.5 years; 55 male) with ccRCC confirmed by histopathology (grade 1, 7; grade 2, 45; grade 3, 18; grade 4, 9) following surgery.

FIELD STRENGTH/SEQUENCE

3.0 T MRI scanner. DR-CSI with a diffusion-weighted echo-planar imaging sequence and T2-mapping with a multi-echo spin echo sequence.

ASSESSMENT

DR-CSI results were analyzed for the solid tumor regions of interest using spectrum segmentation with five sub-region volume fraction metrics (VA , VB , VC , VD , and VE ). The regulations for spectrum segmentation were determined based on the D-T2 spectra of distinct macro-components. Tumor size, voxel-wise T2, and apparent diffusion coefficient (ADC) values were obtained. Histopathology assessed tumor grade (G1-G4) for each case.

STATISTICAL TESTS

One-way ANOVA or Kruskal-Wallis test, Spearman's correlation (coefficient, rho), multivariable logistic regression analysis, receiver operating characteristic curve analysis, and DeLong's test. Significance criteria: P < 0.05.

RESULTS

Significant differences were found in ADC, T2, DR-CSI VB , and VD among the ccRCC grades. Correlations were found for ccRCC grade to tumor size (rho = 0.419), age (rho = 0.253), VB (rho = 0.553) and VD (rho = -0.378). AUC of VB was slightly larger than ADC in distinguishing low-grade (G1-G2) from high-grade (G3-G4) ccRCC (0.801 vs. 0.762, P = 0.406) and G1 from G2 to G4 (0.796 vs. 0.647, P = 0.175), although not significant. Combining VB , VD , and VE had better diagnostic performance than combining ADC and T2 for differentiating G1 from G2-G4 (AUC: 0.814 vs 0.643).

DATA CONCLUSION

DR-CSI parameters are correlated with ccRCC grades, and may help to differentiate ccRCC grades.

EVIDENCE LEVEL

2 TECHNICAL EFFICACY STAGE: 2.

Investigation of clear cell renal cell carcinoma grades using diffusion-relaxation correlation spectroscopic imaging with optimized spatial-spectrum analysis

OBJECTIVES

To differentiate high-grade from low-grade clear cell renal cell carcinoma (ccRCC) using diffusion-relaxation correlation spectroscopic imaging (DR-CSI) spectra in an equal separating analysis.

METHODS

Eighty patients with 86 pathologically confirmed ccRCCs who underwent DR-CSI were enrolled. Two radiologists delineated the region of interest. The spectrum was derived based on DR-CSI and was further segmented into multiple equal subregions from 2*2 to 9*9. The agreement between the 2 radiologists was assessed by the intraclass correlation coefficient (ICC). Logistic regression was used to establish the regression model for differentiation, and 5-fold cross-validation was used to evaluate its accuracy. McNemar's test was used to compare the diagnostic performance between equipartition models and the traditional parameters, including the apparent diffusion coefficient (ADC) and T2 value.

RESULTS

The inter-reader agreement decreased as the divisions in the equipartition model increased (overall ICC ranged from 0.859 to 0.920). The accuracy increased from the 2*2 to 9*9 equipartition model (0.68 for 2*2, 0.69 for 3*3 and 4*4, 0.70 for 5*5, 0.71 for 6*6, 0.78 for 7*7, and 0.75 for 8*8 and 9*9). The equipartition models with divisions >7*7 were significantly better than ADC and T2 (vs ADC: P = .002-.008; vs T2: P = .001-.004).

CONCLUSIONS

The equipartition method has the potential to analyse the DR-CSI spectrum and discriminate between low-grade and high-grade ccRCC.

ADVANCES IN KNOWLEDGE

The evaluation of DR-CSI relies on prior knowledge, and how to assess the spectrum derived from DR-CSI without prior knowledge has not been well studied.

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.

In vivo disentanglement of diffusion frequency-dependence, tensor shape, and relaxation using multidimensional MRI

Diffusion MRI with free gradient waveforms, combined with simultaneous relaxation encoding, referred to as multidimensional MRI (MD-MRI), offers microstructural specificity in complex biological tissue. This approach delivers intravoxel information about the microstructure, local chemical composition, and importantly, how these properties are coupled within heterogeneous tissue containing multiple microenvironments. Recent theoretical advances incorporated diffusion time dependency and integrated MD-MRI with concepts from oscillating gradients. This framework probes the diffusion frequency, ω , in addition to the diffusion tensor, D , and relaxation, R 1 , R 2 , correlations. A D ( ω ) - R 1 - R 2 clinical imaging protocol was then introduced, with limited brain coverage and 3 mm3 voxel size, which hinder brain segmentation and future cohort studies. In this study, we introduce an efficient, sparse in vivo MD-MRI acquisition protocol providing whole brain coverage at 2 mm3 voxel size. We demonstrate its feasibility and robustness using a well-defined phantom and repeated scans of five healthy individuals. Additionally, we test different denoising strategies to address the sparse nature of this protocol, and show that efficient MD-MRI encoding design demands a nuanced denoising approach. The MD-MRI framework provides rich information that allows resolving the diffusion frequency dependence into intravoxel components based on their D ( ω ) - R 1 - R 2 distribution, enabling the creation of microstructure-specific maps in the human brain. Our results encourage the broader adoption and use of this new imaging approach for characterizing healthy and pathological tissues.

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.

Diffusion-relaxation correlation spectrum imaging for predicting tumor consistency and gross total resection in patients with pituitary adenomas: a preliminary study

OBJECTIVE

To evaluate the ability of diffusion-relaxation correlation spectrum imaging (DR-CSI) to predict the consistency and extent of resection (EOR) of pituitary adenomas (PAs).

METHODS

Forty-four patients with PAs were prospectively enrolled. Tumor consistency was evaluated at surgery as either soft or hard, followed by histological assessment. In vivo DR-CSI was performed and spectra were segmented following to a peak-based strategy into four compartments, designated A (low ADC), B (mediate ADC, short T2), C (mediate ADC, long T2), and D (high ADC). The corresponding volume fractions ([Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text]) along with the ADC and T2 values were calculated and assessed using univariable analysis for discrimination between hard and soft PAs. Predictors of EOR > 95% were analyzed using logistic regression model and receiver-operating-characteristic analysis.

RESULTS

Tumor consistency was classified as soft (n = 28) or hard (n = 16). Hard PAs presented higher [Formula: see text] (p = 0.001) and lower [Formula: see text] (p = 0.013) than soft PAs, while no significant difference was found in other parameters. [Formula: see text] significantly correlated with the level of collagen content (r = 0.448, p = 0.002). Knosp grade (odds ratio [OR], 0.299; 95% confidence interval [CI], 0.124-0.716; p = 0.007) and [Formula: see text] (OR, 0.834, per 1% increase; 95% CI, 0.731-0.951; p = 0.007) were independently associated with EOR > 95%. A prediction model based on these variables yielded an AUC of 0.934 (sensitivity, 90.9%; specificity, 90.9%), outperforming the Knosp grade alone (AUC, 0.785; p < 0.05).

CONCLUSION

DR-CSI may serve as a promising tool to predict the consistency and EOR of PAs.

CLINICAL RELEVANCE STATEMENT

DR-CSI provides an imaging dimension for characterizing tissue microstructure of PAs and may serve as a promising tool to predict the tumor consistency and extent of resection in patients with PAs.

KEY POINTS

• DR-CSI provides an imaging dimension for characterizing tissue microstructure of PAs by visualizing the volume fraction and corresponding spatial distribution of four compartments ([Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text]). • [Formula: see text] correlated with the level of collagen content and may be the best DR-CSI parameter for discrimination between hard and soft PAs. • The combination of Knosp grade and [Formula: see text] achieved an AUC of 0.934 for predicting the total or near-total resection, outperforming the Knosp grade alone (AUC, 0.785).

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.

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.

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.

A Multicompartmental Diffusion Model for Improved Assessment of Whole-Body Diffusion-weighted Imaging Data and Evaluation of Prostate Cancer Bone Metastases

Purpose To develop a multicompartmental signal model for whole-body diffusion-weighted imaging (DWI) and apply it to study the diffusion properties of normal tissue and metastatic prostate cancer bone lesions in vivo. Materials and Methods This prospective study (ClinicalTrials.gov: NCT03440554) included 139 men with prostate cancer (mean age, 70 years ± 9 [SD]). Multicompartmental models with two to four tissue compartments were fit to DWI data from whole-body scans to determine optimal compartmental diffusion coefficients. Bayesian information criterion (BIC) and model-fitting residuals were calculated to quantify model complexity and goodness of fit. Diffusion coefficients for the optimal model (having lowest BIC) were used to compute compartmental signal-contribution maps. The signal intensity ratio (SIR) of bone lesions to normal-appearing bone was measured on these signal-contribution maps and on conventional DWI scans and compared using paired t tests (α = .05). Two-sample t tests (α = .05) were used to compare compartmental signal fractions between lesions and normal-appearing bone. Results Lowest BIC was observed from the four-compartment model, with optimal compartmental diffusion coefficients of 0, 1.1 × 10^-3, 2.8 × 10^-3, and >3.0 ×10^-2 mm²/sec. Fitting residuals from this model were significantly lower than from conventional apparent diffusion coefficient mapping (P < .001). Bone lesion SIR was significantly higher on signal-contribution maps of model compartments 1 and 2 than on conventional DWI scans (P < .008). The fraction of signal from compartments 2, 3, and 4 was also significantly different between metastatic bone lesions and normal-appearing bone tissue (P ≤ .02). Conclusion The four-compartment model best described whole-body diffusion properties. Compartmental signal contributions from this model can be used to examine prostate cancer bone involvement. Keywords: Whole-Body MRI, Diffusion-weighted Imaging, Restriction Spectrum Imaging, Diffusion Signal Model, Bone Metastases, Prostate Cancer Clinical trial registration no. NCT03440554 Supplemental material is available for this article. © RSNA, 2023 See also commentary by Margolis in this issue.

Detection of prostate cancer using diffusion-relaxation correlation spectrum imaging with support vector machine model - a feasibility study

BACKGROUND

To evaluate the performance of diffusion-relaxation correlation spectrum imaging (DR-CSI) with support vector machine (SVM) in detecting prostate cancer (PCa).

METHODS

In total, 114 patients (mean age, 66 years, range, 48-87 years) who received a prostate MRI and underwent biopsy were enrolled in three stages. Thirty-nine were assigned for the exploration stage to establish the model, 18 for the validation stage to choose the appropriate scale for mapping and 57 for the test stage to compare the diagnostic performance of the DR-CSI and PI-RADS.

RESULTS

In the exploration stage, the DR-CSI model was established and performed better than the ADC and T₂ values (both P < 0.001). The validation result shows that at least 2 pixels were required for both the long-axis and short-axis in the mapping procedure. In the test stage, DR-CSI had higher accuracy than PI-RADS ≥ 3 as a positive finding based on patient (84.2% vs. 63.2%, P = 0.004) and lesion (78.8% vs. 57.6%, P = 0.001) as well as PI-RADS ≥ 4 on lesion (76.5% vs. 64.7%, P = 0.029), while there was no significant difference between DR-CSI and PI-RADS ≥ 4 based on patient (P = 0.508). For clinically significant PCa, DR-CSI had higher accuracy than PI-RADS ≥ 3 based on patients (84.2% vs. 63.2%, P = 0.004) and lesions (62.4% vs. 48.2%, P = 0.036). There was no significant difference between DR-CSI and PI-RADS ≥ 4 (P = 1.000 and 0.845 for the patient and lesion levels, respectively).

CONCLUSIONS

DR-CSI combined with the SVM model may improve the diagnostic accuracy of PCa.

TRIAL REGISTRATION

This study was approved by the Ethics Committee of our institute (Approval No. KY2018-213). Written informed consent was obtained from all participants.

COnstrained Reference frame diffusion TEnsor Correlation Spectroscopic (CORTECS) MRI: A practical framework for high-resolution diffusion tensor distribution imaging

High-resolution imaging studies have consistently shown that in cortical tissue water diffuses preferentially along radial and tangential orientations with respect to the cortical surface, in agreement with histology. These dominant orientations do not change significantly even if the relative contributions from microscopic water pools to the net voxel signal vary across experiments that use different diffusion times, b-values, TEs, and TRs. With this in mind, we propose a practical new framework for imaging non-parametric diffusion tensor distributions (DTDs) by constraining the microscopic diffusion tensors of the DTD to be diagonalized using the same orthonormal reference frame of the mesoscopic voxel. In each voxel, the constrained DTD (cDTD) is completely determined by the correlation spectrum of the microscopic principal diffusivities associated with the axes of the voxel reference frame. Consequently, all cDTDs are inherently limited to the domain of positive definite tensors and can be reconstructed efficiently using Inverse Laplace Transform methods. Moreover, the cDTD reconstruction can be performed using only data acquired efficiently with single diffusion encoding, although it also supports datasets with multiple diffusion encoding. In tissues with a well-defined architecture, such as the cortex, we can further constrain the cDTD to contain only cylindrically symmetric diffusion tensors and measure the 2D correlation spectra of principal diffusivities along the radial and tangential orientation with respect to the cortical surface. To demonstrate this framework, we perform numerical simulations and analyze high-resolution dMRI data from a fixed macaque monkey brain. We estimate 2D cDTDs in the cortex and derive, in each voxel, the marginal distributions of the microscopic principal diffusivities, the corresponding distributions of the microscopic fractional anisotropies and mean diffusivities along with their 2D correlation spectra to quantify the cDTD shape-size characteristics. Signal components corresponding to specific bands in these cDTD-derived spectra show high specificity to cortical laminar structures observed with histology. Our framework drastically simplifies the measurement of non-parametric DTDs in high-resolution datasets with mesoscopic voxel sizes much smaller than the radius of curvature of the underlying anatomy, e.g., cortical surface, and can be applied retrospectively to analyze existing diffusion MRI data from fixed cortical tissues.

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 B0 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 B0 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.

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.

Improved diffusion parameter estimation by incorporating T2 relaxation properties into the DKI-FWE model

The free water elimination (FWE) model and its kurtosis variant (DKI-FWE) can separate tissue and free water signal contributions, thus providing tissue-specific diffusional information. However, a downside of these models is that the associated parameter estimation problem is ill-conditioned, necessitating the use of advanced estimation techniques that can potentially bias the parameter estimates. In this work, we propose the T₂-DKI-FWE model that exploits the T₂ relaxation properties of both compartments, thereby better conditioning the parameter estimation problem and providing, at the same time, an additional potential biomarker (the T₂ of tissue). In our approach, the T₂ of tissue is estimated as an unknown parameter, whereas the T₂ of free water is assumed known a priori and fixed to a literature value (1573 ms). First, the error propagation of an erroneous assumption on the T₂ of free water is studied. Next, the improved conditioning of T₂-DKI-FWE compared to DKI-FWE is illustrated using the Cramér-Rao lower bound matrix. Finally, the performance of the T₂-DKI-FWE model is compared to that of the DKI-FWE and T₂-DKI models on both simulated and real datasets. The error due to a biased approximation of the T₂ of free water was found to be relatively small in various diffusion metrics and for a broad range of erroneous assumptions on its underlying ground truth value. Compared to DKI-FWE, using the T₂-DKI-FWE model is beneficial for the identifiability of the model parameters. Our results suggest that the T₂-DKI-FWE model can achieve precise and accurate diffusion parameter estimates, through effective reduction of free water partial volume effects and by using a standard nonlinear least squares approach. In conclusion, incorporating T₂ relaxation properties into the DKI-FWE model improves the conditioning of the model fitting, while only requiring an acquisition scheme with at least two different echo times.

Stabilization of parameter estimates from multiexponential decay through extension into higher dimensions

Analysis of multiexponential decay has remained a topic of active research for over 200 years. This attests to the widespread importance of this problem and to the profound difficulties in characterizing the underlying monoexponential decays. Here, we demonstrate the fundamental improvement in stability and conditioning of this classic problem through extension to a second dimension; we present statistical analysis, Monte-Carlo simulations, and experimental magnetic resonance relaxometry data to support this remarkable fact. Our results are readily generalizable to higher dimensions and provide a potential means of circumventing conventional limits on multiexponential parameter estimation.

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.

Data-driven algorithm for myelin water imaging: Probing subvoxel compartmentation based on identification of spatially global tissue features

PURPOSE

Multicomponent analysis of MRI T₂ relaxation time (mcT₂ ) is commonly used for estimating myelin content by separating the signal at each voxel into its underlying distribution of T₂ values. This voxel-based approach is challenging due to the large ambiguity in the multi-T₂ space and the low SNR of MRI signals. Herein, we present a data-driven mcT₂ analysis, which utilizes the statistical strength of identifying spatially global mcT₂ motifs in white matter segments before deconvolving the local signal at each voxel.

METHODS

Deconvolution is done using a tailored optimization scheme, which incorporates the global mcT₂ motifs without additional prior assumptions regarding the number of microscopic components. The end results of this process are voxel-wise myelin water fraction maps.

RESULTS

Validations are shown for computer-generated signals, uniquely designed subvoxel mcT₂ phantoms, and in vivo human brain. Results demonstrated excellent fitting accuracy, both for the numerical and the physical mcT₂ phantoms, exhibiting excellent agreement between calculated myelin water fraction and ground truth. Proof-of-concept in vivo validation is done by calculating myelin water fraction maps for white matter segments of the human brain. Interscan stability of myelin water fraction values was also estimated, showing good correlation between scans.

CONCLUSION

We conclude that studying global tissue motifs prior to performing voxel-wise mcT₂ analysis stabilizes the optimization scheme and efficiently overcomes the ambiguity in the T₂ space. This new approach can improve myelin water imaging and the investigation of microstructural compartmentation in general.

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.

Emerging placental biomarkers of health and disease through advanced magnetic resonance imaging (MRI)

Placental dysfunction is a major cause of fetal demise, fetal growth restriction, and preterm birth, as well as significant maternal morbidity and mortality. Infant survivors of placental dysfunction are at elevatedrisk for lifelong neuropsychiatric morbidity. However, despite the significant consequences of placental disease, there are no clinical tools to directly and non-invasively assess and measure placental function in pregnancy. In this work, we will review advanced MRI techniques applied to the study of the in vivo human placenta in order to better detail placental structure, architecture, and function. We will discuss the potential of these measures to serve as optimal biomarkers of placental dysfunction and review the evidence of these tools in the discrimination of health and disease in pregnancy. Efforts to advance our understanding of in vivo placental development are necessary if we are to optimize healthy pregnancy outcomes and prevent brain injury in successive generations. Current management of many high-risk pregnancies cannot address placental maldevelopment or injury, given the standard tools available to clinicians. Once accurate biomarkers of placental development and function are constructed, the subsequent steps will be to introduce maternal and fetal therapeutics targeting at optimizing placental function. Applying these biomarkers in future studies will allow for real-time assessments of safety and efficacy of novel interventions aimed at improving maternal-fetal well-being.

Quantitative diffusion MRI of the abdomen and pelvis

Diffusion MRI has enormous potential and utility in the evaluation of various abdominal and pelvic disease processes including cancer and noncancer imaging of the liver, prostate, and other organs. Quantitative diffusion MRI is based on acquisitions with multiple diffusion encodings followed by quantitative mapping of diffusion parameters that are sensitive to tissue microstructure. Compared to qualitative diffusion-weighted MRI, quantitative diffusion MRI can improve standardization of tissue characterization as needed for disease detection, staging, and treatment monitoring. However, similar to many other quantitative MRI methods, diffusion MRI faces multiple challenges including acquisition artifacts, signal modeling limitations, and biological variability. In abdominal and pelvic diffusion MRI, technical acquisition challenges include physiologic motion (respiratory, peristaltic, and pulsatile), image distortions, and low signal-to-noise ratio. If unaddressed, these challenges lead to poor technical performance (bias and precision) and clinical outcomes of quantitative diffusion MRI. Emerging and novel technical developments seek to address these challenges and may enable reliable quantitative diffusion MRI of the abdomen and pelvis. Through systematic validation in phantoms, volunteers, and patients, including multicenter studies to assess reproducibility, these emerging techniques may finally demonstrate the potential of quantitative diffusion MRI for abdominal and pelvic imaging applications.

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.

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.

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.

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.

SNR-enhanced diffusion MRI with structure-preserving low-rank denoising in reproducing kernel Hilbert spaces

PURPOSE

To introduce, develop, and evaluate a novel denoising technique for diffusion MRI that leverages nonlinear redundancy in the data to boost the SNR while preserving signal information.

METHODS

We exploit nonlinear redundancy of the dMRI data by means of kernel principal component analysis (KPCA), a nonlinear generalization of PCA to reproducing kernel Hilbert spaces. By mapping the signal to a high-dimensional space, a higher level of redundant information is exploited, thereby enabling better denoising than linear PCA. We implement KPCA with a Gaussian kernel, with parameters automatically selected from knowledge of the noise statistics, and validate it on realistic Monte Carlo simulations as well as with in vivo human brain submillimeter and low-resolution dMRI data. We also demonstrate KPCA denoising on multi-coil dMRI data.

RESULTS

SNR improvements up to 2.7 × were obtained in real in vivo datasets denoised with KPCA, in comparison to SNR gains of up to 1.8 × using a linear PCA denoising technique called Marchenko-Pastur PCA (MPPCA). Compared to gold-standard dataset references created from averaged data, we showed that lower normalized root mean squared error was achieved with KPCA compared to MPPCA. Statistical analysis of residuals shows that anatomical information is preserved and only noise is removed. Improvements in the estimation of diffusion model parameters such as fractional anisotropy, mean diffusivity, and fiber orientation distribution functions were also demonstrated.

CONCLUSION

Nonlinear redundancy of the dMRI signal can be exploited with KPCA, which allows superior noise reduction/SNR improvements than the MPPCA method, without loss of signal information.

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.

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.