Multidimensional correlation MRI

Simple MATLAB and Python scripts for multi-exponential analysis

Multi-exponential decay is prevalent in magnetic resonance spectroscopy, relaxation, and imaging. This paper describes simple MATLAB and Python functions and scripts for regularized multi-exponential analysis methods for 1D and 2D data and example test problems and experiments. Regularized least-squares solutions provide production-quality outputs with robust stopping rules in ~5 and ~20 lines of code for 1D and 2D inversions, respectively. The software provides an open-architecture simple solution for transforming exponential decay data to the distribution of their decay lifetimes. Examples from magnetic resonance relaxation of a complex fluid, a Danish North Sea crude oil, and fluid mixtures in porous materials-brine/crude oil mixture in North Sea reservoir chalk-are presented. Developed codes may be incorporated in other software or directly used by other researchers, in magnetic resonance relaxation, diffusion, and imaging or other physical phenomena that require multi-exponential analysis.

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

OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSION

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

CRITICAL RELEVANCE STATEMENT

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

KEY POINTS

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

Excessive interstitial free-water in cortical gray matter preceding accelerated volume changes in individuals at clinical high risk for psychosis

Recent studies show that accelerated cortical gray matter (GM) volume reduction seen in anatomical MRI can help distinguish between individuals at clinical high risk (CHR) for psychosis who will develop psychosis and those who will not. This reduction is suggested to represent atypical developmental or degenerative changes accompanying an accumulation of microstructural changes, such as decreased spine density and dendritic arborization. Detecting the microstructural sources of these changes before they accumulate into volume loss is crucial. Our study aimed to detect these microstructural GM alterations using diffusion MRI (dMRI). We tested for baseline and longitudinal group differences in anatomical and dMRI data from 160 individuals at CHR and 96 healthy controls (HC) acquired in a single imaging site. Of the CHR individuals, 33 developed psychosis (CHR-P), while 127 did not (CHR-NP). Among all participants, longitudinal data was available for 45 HCs, 17 CHR-P, and 66 CHR-NP. Eight cortical lobes were examined for GM volume and GM microstructure. A novel dMRI measure, interstitial free water (iFW), was used to quantify GM microstructure by eliminating cerebrospinal fluid contribution. Additionally, we assessed whether these measures differentiated the CHR-P from the CHR-NP. In addition, for completeness, we also investigated changes in cortical thickness and in white matter (WM) microstructure. At baseline the CHR group had significantly higher iFW than HC in the prefrontal, temporal, parietal, and occipital lobes, while volume was reduced only in the temporal lobe. Neither iFW nor volume differentiated between the CHR-P and CHR-NP groups at baseline. However, in many brain areas, the CHR-P group demonstrated significantly accelerated changes (iFW increase and volume reduction) with time than the CHR-NP group. Cortical thickness provided similar results as volume, and there were no significant changes in WM microstructure. Our results demonstrate that microstructural GM changes in individuals at CHR have a wider extent than volumetric changes or microstructural WM changes, and they predate the acceleration of brain changes that occur around psychosis onset. Microstructural GM changes, as reflected by the increased iFW, are thus an early pathology at the prodromal stage of psychosis that may be useful for a better mechanistic understanding of psychosis development.

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.

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.

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.

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.

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.

Correspondence of mean apparent propagator MRI metrics with phosphorylated tau and astrogliosis in chronic traumatic encephalopathy

Chronic traumatic encephalopathy is a neurodegenerative disease that is diagnosed and staged based on the localization and extent of phosphorylated tau pathology. Although its identification remains the primary diagnostic criteria to distinguish chronic traumatic encephalopathy from other tauopathies, the hyperphosphorylated tau that accumulates in neurofibrillary tangles in cortical grey matter and perivascular regions is often accompanied by concomitant pathology such as astrogliosis. Mean apparent propagator MRI is a clinically feasible diffusion MRI method that is suitable to characterize microstructure of complex biological media efficiently and comprehensively. We performed quantitative correlations between propagator metrics and underlying phosphorylated tau and astroglial pathology in a cross-sectional study of 10 ex vivo human tissue specimens with 'high chronic traumatic encephalopathy' at 0.25 mm isotropic voxels. Linear mixed effects analysis of regions of interest showed significant relationships of phosphorylated tau with propagator-estimated non-Gaussianity in cortical grey matter (P = 0.002) and of astrogliosis with propagator anisotropy in superficial cortical white matter (P = 0.0009). The positive correlation between phosphorylated tau and non-Gaussianity was found to be modest but significant (R2 = 0.44, P = 6.0 × 10-5) using linear regression. We developed an unsupervised clustering algorithm with non-Gaussianity and propagator anisotropy as inputs, which was able to identify voxels in superficial cortical white matter that corresponded to astrocytes that were accumulated at the grey-white matter interface. Our results suggest that mean apparent propagator MRI at high spatial resolution provides a means to not only identify phosphorylated tau pathology but also detect regions with astrocytic pathology and may therefore prove diagnostically valuable in the evaluation of concomitant pathology in cortical tissue with complex microstructure.

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

Microstructural Cortical Gray Matter Changes Preceding Accelerated Volume Changes in Individuals at Clinical High Risk for Psychosis

Recent studies show that accelerated cortical gray matter (GM) volume reduction seen in anatomical MRI can help distinguish between individuals at clinical high risk (CHR) for psychosis who will develop psychosis and those who will not. This reduction is thought to result from an accumulation of microstructural changes, such as decreased spine density and dendritic arborization. Detecting the microstructural sources of these changes before they accumulate is crucial, as volume reduction likely indicates an underlying neurodegenerative process. Our study aimed to detect these microstructural GM alterations using diffusion MRI (dMRI). We tested for baseline and longitudinal group differences in anatomical and dMRI data from 160 individuals at CHR and 96 healthy controls (HC) acquired in a single imaging site. Eight cortical lobes were examined for GM volume and GM microstructure. A novel dMRI measure, interstitial free water (iFW), was used to quantify GM microstructure by eliminating cerebrospinal fluid contribution. Additionally, we assessed whether these measures differentiated the 33 individuals at CHR who developed psychosis (CHR-P) from the 127 individuals at CHR who did not (CHR-NP). At baseline the CHR group had significantly higher iFW than HC in the prefrontal, temporal, parietal, and occipital lobes, while volume was reduced only in the temporal lobe. Neither iFW nor volume differentiated between the CHR-P and CHR-NP groups at baseline. However, in most brain areas, the CHR-P group demonstrated significantly accelerated iFW increase and volume reduction with time than the CHR-NP group. Our results demonstrate that microstructural GM changes in individuals at CHR have a wider extent than volumetric changes and they predate the acceleration of brain changes that occur around psychosis onset. Microstructural GM changes are thus an early pathology at the prodromal stage of psychosis that may be useful for early detection and a better mechanistic understanding of psychosis development.

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.

Adult lifespan maturation and degeneration patterns in gray and white matter: A mean apparent propagator (MAP) MRI study

The relationship between brain microstructure and aging has been the subject of intense study, with diffusion MRI perhaps the most effective modality for elucidating these associations. Here, we used the mean apparent propagator (MAP)-MRI framework, which is suitable to characterize complex microstructure, to investigate age-related cerebral differences in a cohort of cognitively unimpaired participants and compared the results to those derived using diffusion tensor imaging. We studied MAP-MRI metrics, among them the non-Gaussianity (NG) and propagator anisotropy (PA), and established an opposing pattern in white matter of higher NG alongside lower PA among older adults, likely indicative of axonal degradation. In gray matter, however, these two indices were consistent with one another, and exhibited regional pattern heterogeneity compared to other microstructural parameters, which could indicate fewer neuronal projections across cortical layers along with an increased glial concentration. In addition, we report regional variations in the magnitude of age-related microstructural differences consistent with the posterior-anterior shift in aging paradigm. These results encourage further investigations in cognitive impairments and neurodegeneration.

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.

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.

Comparison of experimental times in T1-D and D-T2 correlation experiments in single-sided NMR

Diffusion-relaxation correlation experiments in nuclear magnetic resonance are a powerful technique for the characterization of fluid dynamics in confined geometries or soft matter, in which relaxation may be either spin-spin (T₂) or spin-lattice (T₁). The general approach is to acquire a set of bidimensional data in which diffusion is codified by the evolution of the magnetization with either Hahn or stimulated echoes (STE) in the presence of a constant magnetic field gradient. While T₂ is codified by a Carr-Purcell-Meiboom-Gil (CPMG) sequence, T₁ is either encoded by saturation or inversion-recovery methods. In this work, we analyse the measurement time of diffusion-relaxation times in single-sided NMR and show that T₁-D acquisition is always shorter than D-T₂. Depending on the hardware characteristics, this time reduction can be up to an order of magnitude. We present analytical calculations and examples in model porous media saturated with water and in a dairy product.

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.

Ultrafast methods for relaxation and diffusion

Relaxation and diffusion NMR measurements offer an approach to studying rotational and translational motion of molecules non-invasively, and they also provide chemical resolution complementary to NMR spectra. Multidimensional experiments enable the correlation of relaxation and diffusion parameters as well as the observation of molecular exchange phenomena through relaxation or diffusion contrast. This review describes how to accelerate multidimensional relaxation and diffusion measurements significantly through spatial encoding. This so-called ultrafast Laplace NMR approach shortens the experiment time to a fraction and makes even single-scan experiments possible. Single-scan experiments, in turn, significantly facilitate the use of nuclear spin hyperpolarization methods to boost sensitivity. The ultrafast Laplace NMR method is also applicable with low-field, mobile NMR instruments, and it can be exploited in many disciplines. For example, it has been used in studies of the dynamics of fluids in porous materials, identification of intra- and extracellular metabolites in cancer cells, and elucidation of aggregation phenomena in atmospheric surfactant solutions.

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.

Low-field and variable-field NMR relaxation studies of H2O and D2O molecular dynamics in articular cartilage

Osteoarthritis (OA) as the main degenerative disease of articular cartilage in joints is accompanied by structural and compositional changes in the tissue. Degeneration is a consequence of a reduction of the amount of macromolecules, the so-called proteoglycans, and of a corresponding increase in water content, both leading to structural weakening of cartilage. NMR investigations of cartilage generally address only the relaxation properties of water. In this study, two-dimensional (T1-T2) measurements of bovine articular cartilage samples were carried out for different stages of hydration, complemented by molecular exchange with D2O and treatment by trypsin which simulates degeneration by OA. Two signal components were identified in all measurements, characterized by very different T2 which suggests liquid-like and solid-like dynamics. These measurements allow the quantification of separate hydrogen components and their assignment to defined physical pools which had been discussed repeatedly in the literature, i.e. bulk-like water and a combination of protein hydrogens and strongly bound water. The first determination of 2H relaxation dispersion in comparison to 1H dispersion suggests intramolecular interactions as the dominating source for the pronounced magnetic field dependence of the longitudinal relaxation time T1.

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.

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.

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.