The link between diffusion MRI and tumor heterogeneity: Mapping cell eccentricity and density by diffusional variance decomposition (DIVIDE)

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.

Comparison of methods for intravoxel incoherent motion parameter estimation in the brain from flow-compensated and non-flow-compensated diffusion-encoded data

PURPOSE

Joint analysis of flow-compensated (FC) and non-flow-compensated (NC) diffusion MRI (dMRI) data has been suggested for increased robustness of intravoxel incoherent motion (IVIM) parameter estimation. For this purpose, a set of methods commonly used or previously found useful for IVIM analysis of dMRI data obtained with conventional diffusion encoding were evaluated in healthy human brain.

METHODS

Five methods for joint IVIM analysis of FC and NC dMRI data were compared: (1) direct non-linear least squares fitting, (2) a segmented fitting algorithm with estimation of the diffusion coefficient from higher b-values of NC data, (3) a Bayesian algorithm with uniform prior distributions, (4) a Bayesian algorithm with spatial prior distributions, and (5) a deep learning-based algorithm. Methods were evaluated on brain dMRI data from healthy subjects and simulated data at multiple noise levels. Bipolar diffusion encoding gradients were used with b-values 0-200 s/mm2 and corresponding flow weighting factors 0-2.35 s/mm for NC data and by design 0 for FC data. Data were acquired twice for repeatability analysis.

RESULTS

Measurement repeatability as well as estimation bias and variability were at similar levels or better with the Bayesian algorithm with spatial prior distributions and the deep learning-based algorithm for IVIM parametersD $$ D $$ andf $$ f $$ , and for the Bayesian algorithm only forv d $$ {v}_d $$ , relative to the other methods.

CONCLUSION

A Bayesian algorithm with spatial prior distributions is preferable for joint IVIM analysis of FC and NC dMRI data in the healthy human brain, but deep learning-based algorithms appear promising.

Tensor-valued diffusion MRI of human acute stroke

PURPOSE

Tensor-valued diffusion encoding can disentangle orientation dispersion and subvoxel anisotropy, potentially offering insight into microstructural changes after cerebral ischemia. The purpose was to evaluate tensor-valued diffusion MRI in human acute ischemic stroke, assess potential confounders from diffusion time dependencies, and compare to Monte Carlo diffusion simulations of axon beading.

METHODS

Linear (LTE) and spherical (STE) b-tensor encoding with inherently different effective diffusion times were acquired in 21 acute ischemic stroke patients between 3 and 57 h post-onset at 3 T in 2.5 min. In an additional 10 patients, STE with 2 LTE yielding different effective diffusion times were acquired for comparison. Diffusional variance decomposition (DIVIDE) was used to estimate microscopic anisotropy (μFA), as well as anisotropic, isotropic, and total diffusional variance (MKA , MKI , MKT ). DIVIDE parameters, and diffusion tensor imaging (DTI)-derived mean diffusivity and fractional anisotropy (FA) were compared in lesion versus contralateral white matter. Monte Carlo diffusion simulations of various cylindrical geometries for all b-tensor protocols were used to interpret parameter measurements.

RESULTS

MD was ˜40% lower in lesions for all LTE/STE protocols. The DIVIDE parameters varied with effective diffusion time: higher μFA and MKA in lesion versus contralateral white matter for STE with longer effective diffusion time LTE, whereas the shorter effective diffusion time LTE protocol yielded lower μFA and MKA in lesions. Both protocols, regardless of diffusion time, were consistent with simulations of greater beading amplitude and intracellular volume fraction.

CONCLUSION

DIVIDE parameters depend on diffusion time in acute stroke but consistently indicate neurite beading and larger intracellular volume fraction.

The diagnostic efficiency of integration of 2HG MRS and IVIM versus individual parameters for predicting IDH mutation status in gliomas in clinical scenarios: A retrospective study

PURPOSE

Currently, there remains a scarcity of established preoperative tests to accurately predict the isocitrate dehydrogenase (IDH) mutation status in clinical scenarios, with limited research has explored the potential synergistic diagnostic performance among metabolite, perfusion, and diffusion parameters. To address this issue, we aimed to develop an imaging protocol that integrated 2-hydroxyglutarate (2HG) magnetic resonance spectroscopy (MRS) and intravoxel incoherent motion (IVIM) by comprehensively assessing metabolic, cellular, and angiogenic changes caused by IDH mutations, and explored the diagnostic efficiency of this imaging protocol for predicting IDH mutation status in clinical scenarios.

METHODS

Patients who met the inclusion criteria were categorized into two groups: IDH-wild type (IDH-WT) group and IDH-mutant (IDH-MT) group. Subsequently, we quantified the 2HG concentration, the relative apparent diffusion coefficient (rADC), the relative true diffusion coefficient value (rD), the relative pseudo-diffusion coefficient (rD*) and the relative perfusion fraction value (rf). Intergroup differences were estimated using t-test and Mann-Whitney U test. Finally, we performed receiver operating characteristic (ROC) curve and DeLong's test to evaluate and compare the diagnostic performance of individual parameters and their combinations.

RESULTS

64 patients (female, 21; male, 43; age, 47.0 ± 13.7 years) were enrolled. Compared with IDH-WT gliomas, IDH-MT gliomas had higher 2HG concentration, rADC and rD (P < 0.001), and lower rD* (P = 0.013). The ROC curve demonstrated that 2HG + rD + rD* exhibited the highest areas under curve (AUC) value (0.967, 95%CI 0.889-0.996) for discriminating IDH mutation status. Compared with each individual parameter, the predictive efficiency of 2HG + rADC + rD* and 2HG + rD + rD* shows a statistically significant enhancement (DeLong's test: P < 0.05).

CONCLUSIONS

The integration of 2HG MRS and IVIM significantly improves the diagnostic efficiency for predicting IDH mutation status in clinical scenarios.

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.

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.

MRI with generalized diffusion encoding reveals damaged white matter in patients previously hospitalized for COVID-19 and with persisting symptoms at follow-up

There is mounting evidence of the long-term effects of COVID-19 on the central nervous system, with patients experiencing diverse symptoms, often suggesting brain involvement. Conventional brain MRI of these patients shows unspecific patterns, with no clear connection of the symptomatology to brain tissue abnormalities, whereas diffusion tensor studies and volumetric analyses detect measurable changes in the brain after COVID-19. Diffusion MRI exploits the random motion of water molecules to achieve unique sensitivity to structures at the microscopic level, and new sequences employing generalized diffusion encoding provide structural information which are sensitive to intravoxel features. In this observational study, a total of 32 persons were investigated: 16 patients previously hospitalized for COVID-19 with persisting symptoms of post-COVID condition (mean age 60 years: range 41-79, all male) at 7-month follow-up and 16 matched controls, not previously hospitalized for COVID-19, with no post-COVID symptoms (mean age 58 years, range 46-69, 11 males). Standard MRI and generalized diffusion encoding MRI were employed to examine the brain white matter of the subjects. To detect possible group differences, several tissue microstructure descriptors obtainable with the employed diffusion sequence, the fractional anisotropy, mean diffusivity, axial diffusivity, radial diffusivity, microscopic anisotropy, orientational coherence (Cc) and variance in compartment's size (CMD) were analysed using the tract-based spatial statistics framework. The tract-based spatial statistics analysis showed widespread statistically significant differences (P < 0.05, corrected for multiple comparisons using the familywise error rate) in all the considered metrics in the white matter of the patients compared to the controls. Fractional anisotropy, microscopic anisotropy and Cc were lower in the patient group, while axial diffusivity, radial diffusivity, mean diffusivity and CMD were higher. Significant changes in fractional anisotropy, microscopic anisotropy and CMD affected approximately half of the analysed white matter voxels located across all brain lobes, while changes in Cc were mainly found in the occipital parts of the brain. Given the predominant alteration in microscopic anisotropy compared to Cc, the observed changes in diffusion anisotropy are mostly due to loss of local anisotropy, possibly connected to axonal damage, rather than white matter fibre coherence disruption. The increase in radial diffusivity is indicative of demyelination, while the changes in mean diffusivity and CMD are compatible with vasogenic oedema. In summary, these widespread alterations of white matter microstructure are indicative of vasogenic oedema, demyelination and axonal damage. These changes might be a contributing factor to the diversity of central nervous system symptoms that many patients experience after COVID-19.

Extra-axonal contribution to double diffusion encoding-based pore size estimates in the corticospinal tract

OBJECTIVE

To study the origin of compartment size overestimation in double diffusion encoding MRI (DDE) in vivo experiments in the human corticospinal tract. Here, the extracellular space is hypothesized to be the origin of the DDE signal. By exploiting the DDE sensitivity to pore shape, it could be possible to identify the origin of the measured signal. The signal difference between parallel and perpendicular diffusion gradient orientation can indicate if a compartment is regular or eccentric in shape. As extracellular space can be considered an eccentric compartment, a positive difference would mean a high contribution to the compartment size estimates.

MATERIALS AND METHODS

Computer simulations using MISST and in vivo experiments in eight healthy volunteers were performed. DDE experiments using a double spin-echo preparation with eight perpendicular directions were measured in vivo. The difference between parallel and perpendicular gradient orientations was analyzed using a Wilcoxon signed-rank test and a Mann-Whitney U test.

RESULTS

Simulations and MR experiments showed a statistically significant difference between parallel and perpendicular diffusion gradient orientation signals ([Formula: see text]).

CONCLUSION

The results suggest that the DDE-based size estimate may be considerably influenced by the extra-axonal compartment. However, the experimental results are also consistent with purely intra-axonal contributions in combination with a large fiber orientation dispersion.

Neuroinflammation and Mitochondrial Dysfunction in Parkinson's Disease: Connecting Neuroimaging with Pathophysiology

There is a pressing need for disease-modifying therapies in patients suffering from neurodegenerative diseases, including Parkinson's disease (PD). However, these disorders face unique challenges in clinical trial designs to assess the neuroprotective properties of potential drug candidates. One of these challenges relates to the often unknown individual disease mechanisms that would, however, be relevant for targeted treatment strategies. Neuroinflammation and mitochondrial dysfunction are two proposed pathophysiological hallmarks and are considered to be highly interconnected in PD. Innovative neuroimaging methods can potentially help to gain deeper insights into one's predominant disease mechanisms, can facilitate patient stratification in clinical trials, and could potentially map treatment responses. This review aims to highlight the role of neuroinflammation and mitochondrial dysfunction in patients with PD (PwPD). We will specifically introduce different neuroimaging modalities, their respective technical hurdles and challenges, and their implementation into clinical practice. We will gather preliminary evidence for their potential use in PD research and discuss opportunities for future clinical trials.

Cardiac q-space trajectory imaging by motion-compensated tensor-valued diffusion encoding in human heart in vivo

PURPOSE

Tensor-valued diffusion encoding can probe more specific features of tissue microstructure than what is available by conventional diffusion weighting. In this work, we investigate the technical feasibility of tensor-valued diffusion encoding at high b-values with q-space trajectory imaging (QTI) analysis, in the human heart in vivo.

METHODS

Ten healthy volunteers were scanned on a 3T scanner. We designed time-optimal gradient waveforms for tensor-valued diffusion encoding (linear and planar) with second-order motion compensation. Data were analyzed with QTI. Normal values and repeatability were investigated for the mean diffusivity (MD), fractional anisotropy (FA), microscopic FA (μFA), isotropic, anisotropic and total mean kurtosis (MKi, MKa, and MKt), and orientation coherence (Cc ). A phantom, consisting of two fiber blocks at adjustable angles, was used to evaluate sensitivity of parameters to orientation dispersion and diffusion time.

RESULTS

QTI data in the left ventricular myocardium were MD = 1.62 ± 0.07 μm2 /ms, FA = 0.31 ± 0.03, μFA = 0.43 ± 0.07, MKa = 0.20 ± 0.07, MKi = 0.13 ± 0.03, MKt = 0.33 ± 0.09, and Cc  = 0.56 ± 0.22 (mean ± SD across subjects). Phantom experiments showed that FA depends on orientation dispersion, whereas μFA was insensitive to this effect.

CONCLUSION

We demonstrated the first tensor-valued diffusion encoding and QTI analysis in the heart in vivo, along with first measurements of myocardial μFA, MKi, MKa, and Cc . The methodology is technically feasible and provides promising novel biomarkers for myocardial tissue characterization.

Higher order diffusion imaging as a putative index of human sleep-related microstructural changes and glymphatic clearance

The brain has a unique macroscopic waste clearance system, termed the glymphatic system which utilises perivascular tunnels surrounded by astroglia to promote cerebrospinal-interstitial fluid exchange. Rodent studies have demonstrated a marked increase in glymphatic clearance during sleep which has been linked to a sleep-induced expansion of the extracellular space and concomitant reduction in intracellular volume. However, despite being implicated in the pathophysiology of multiple human neurodegenerative disorders, non-invasive techniques for imaging glymphatic clearance in humans are currently limited. Here we acquired multi-shell diffusion weighted MRI (dwMRI) in twenty-one healthy young participants (6 female, 22.3 ± 3.2 years) each scanned twice, once during wakefulness and once during sleep induced by a combination of one night of sleep deprivation and 10 mg of the hypnotic zolpidem 30 min before scanning. To capture hypothesised sleep-associated changes in intra/extracellular space, dwMRI were analysed using higher order diffusion modelling with the prediction that sleep-associated increases in interstitial (extracellular) fluid volume would result in a decrease in diffusion kurtosis, particularly in areas associated with slow wave generation at the onset of sleep. In line with our hypothesis, we observed a global reduction in diffusion kurtosis (t15=2.82, p = 0.006) during sleep as well as regional reductions in brain areas associated with slow wave generation during early sleep and default mode network areas that are highly metabolically active during wakefulness. Analysis with a higher-order representation of diffusion (MAP-MRI) further indicated that changes within the intra/extracellular domain rather than membrane permeability likely underpin the observed sleep-associated decrease in kurtosis. These findings identify higher-order modelling of dwMRI as a potential new non-invasive method for imaging glymphatic clearance and extend rodent findings to suggest that sleep is also associated with an increase in interstitial fluid volume in humans.

Coregistered histology sections with diffusion tensor imaging data at 200 µm resolution in meningioma tumors

A significant problem in diffusion MRI (dMRI) is the lack of understanding regarding which microstructural features account for the variability in the diffusion tensor imaging (DTI) parameters observed in meningioma tumors. A common assumption is that mean diffusivity (MD) and fractional anisotropy (FA) from DTI are inversely proportional to cell density and proportional to tissue anisotropy, respectively. Although these associations have been established across a wide range of tumors, they have been challenged for interpreting within-tumor variations where several additional microstructural features have been suggested as contributing to MD and FA. To facilitate the investigation of the biological underpinnings of DTI parameters, we performed ex-vivo DTI at 200 µm isotropic resolution on sixteen excised meningioma tumor samples. The samples exhibit a variety of microstructural features because the dataset includes meningiomas of six different meningioma types and two different grades. Diffusion-weighted signal (DWI) maps, DWI maps averaged over all directions for given b-value, signal intensities without diffusion encoding (S₀) as well as DTI parameters: MD, FA, in-plane FA (FA_IP), axial diffusivity (AD) and radial diffusivity (RD), were coregistered to Hematoxylin & Eosin- (H&E) and Elastica van Gieson-stained (EVG) histological sections by a non-linear landmark-based approach. Here, we provide DWI signal and DTI maps coregistered to histology sections and describe the pipeline for processing the raw DTI data and the coregistration. The raw, processed, and coregistered data are hosted by Analytic Imaging Diagnostics Arena (AIDA) data hub registry, and software tools for processing are provided via GitHub. We hope that data can be used in research and education concerning the link between the meningioma microstructure and parameters obtained by DTI.

Comprehensive Brain Tumour Characterisation with VERDICT-MRI: Evaluation of Cellular and Vascular Measures Validated by Histology

The aim of this work was to extend the VERDICT-MRI framework for modelling brain tumours, enabling comprehensive characterisation of both intra- and peritumoural areas with a particular focus on cellular and vascular features. Diffusion MRI data were acquired with multiple b-values (ranging from 50 to 3500 s/mm2), diffusion times, and echo times in 21 patients with brain tumours of different types and with a wide range of cellular and vascular features. We fitted a selection of diffusion models that resulted from the combination of different types of intracellular, extracellular, and vascular compartments to the signal. We compared the models using criteria for parsimony while aiming at good characterisation of all of the key histological brain tumour components. Finally, we evaluated the parameters of the best-performing model in the differentiation of tumour histotypes, using ADC (Apparent Diffusion Coefficient) as a clinical standard reference, and compared them to histopathology and relevant perfusion MRI metrics. The best-performing model for VERDICT in brain tumours was a three-compartment model accounting for anisotropically hindered and isotropically restricted diffusion and isotropic pseudo-diffusion. VERDICT metrics were compatible with the histological appearance of low-grade gliomas and metastases and reflected differences found by histopathology between multiple biopsy samples within tumours. The comparison between histotypes showed that both the intracellular and vascular fractions tended to be higher in tumours with high cellularity (glioblastoma and metastasis), and quantitative analysis showed a trend toward higher values of the intracellular fraction (fic) within the tumour core with increasing glioma grade. We also observed a trend towards a higher free water fraction in vasogenic oedemas around metastases compared to infiltrative oedemas around glioblastomas and WHO 3 gliomas as well as the periphery of low-grade gliomas. In conclusion, we developed and evaluated a multi-compartment diffusion MRI model for brain tumours based on the VERDICT framework, which showed agreement between non-invasive microstructural estimates and histology and encouraging trends for the differentiation of tumour types and sub-regions.

Optimal experimental design and estimation for q-space trajectory imaging

Tensor-valued diffusion encoding facilitates data analysis by q-space trajectory imaging. By modeling the diffusion signal of heterogeneous tissues with a diffusion tensor distribution (DTD) and modulating the encoding tensor shape, this novel approach allows disentangling variations in diffusivity from microscopic anisotropy, orientation dispersion, and mixtures of multiple isotropic diffusivities. To facilitate the estimation of the DTD parameters, a parsimonious acquisition scheme coupled with an accurate and precise estimation of the DTD is needed. In this work, we create two precision-optimized acquisition schemes: one that maximizes the precision of the raw DTD parameters, and another that maximizes the precision of the scalar measures derived from the DTD. The improved precision of these schemes compared to a naïve sampling scheme is demonstrated in both simulations and real data. Furthermore, we show that the weighted linear least squares (WLLS) estimator that uses the squared reciprocal of the noisy signal as weights can be biased, whereas the iteratively WLLS estimator with the squared reciprocal of the predicted signal as weights outperforms the conventional unweighted linear LS and nonlinear LS estimators in terms of accuracy and precision. Finally, we show that the use of appropriate constraints can considerably increase the precision of the estimator with only a limited decrease in accuracy.

Deep learning characterization of brain tumours with diffusion weighted imaging

Glioblastoma multiforme (GBM) is one of the most deadly forms of cancer. Methods of characterizing these tumours are valuable for improving predictions of their progression and response to treatment. A mathematical model called the proliferation-invasion (PI) model has been used extensively in the literature to model the growth of these tumours, though it relies on known values of two key parameters: the tumour cell diffusivity and proliferation rate. Unfortunately, these parameters are difficult to estimate in a patient-specific manner, making personalized tumour forecasting challenging. In this paper, we develop and apply a deep learning model capable of making accurate estimates of these key GBM-characterizing parameters while simultaneously producing a full prediction of the tumour progression curve. Our method uses two sets of multi sequence MRI in order to produce estimations and relies on a preprocessing pipeline which includes brain tumour segmentation and conversion to tumour cellularity. We first apply our deep learning model to synthetic tumours to showcase the model's capabilities and identify situations where prediction errors are likely to occur. We then apply our model to a clinical dataset consisting of five patients diagnosed with GBM. For all patients, we derive evidence-based estimates for each of the PI model parameters and predictions for the future progression of the tumour, along with estimates of the parameter uncertainties. Our work provides a new, easily generalizable method for the estimation of patient-specific tumour parameters, which can be built upon to aid physicians in designing personalized treatments.

Combining Multi-Shell Diffusion with Conventional MRI Improves Molecular Diagnosis of Diffuse Gliomas with Deep Learning

The WHO classification since 2016 confirms the importance of integrating molecular diagnosis for prognosis and treatment decisions of adult-type diffuse gliomas. This motivates the development of non-invasive diagnostic methods, in particular MRI, to predict molecular subtypes of gliomas before surgery. At present, this development has been focused on deep-learning (DL)-based predictive models, mainly with conventional MRI (cMRI), despite recent studies suggesting multi-shell diffusion MRI (dMRI) offers complementary information to cMRI for molecular subtyping. The aim of this work is to evaluate the potential benefit of combining cMRI and multi-shell dMRI in DL-based models. A model implemented with deep residual neural networks was chosen as an illustrative example. Using a dataset of 146 patients with gliomas (from grade 2 to 4), the model was trained and evaluated, with nested cross-validation, on pre-operative cMRI, multi-shell dMRI, and a combination of the two for the following classification tasks: (i) IDH-mutation; (ii) 1p/19q-codeletion; and (iii) three molecular subtypes according to WHO 2021. The results from a subset of 100 patients with lower grades gliomas (2 and 3 according to WHO 2016) demonstrated that combining cMRI and multi-shell dMRI enabled the best performance in predicting IDH mutation and 1p/19q codeletion, achieving an accuracy of 75 ± 9% in predicting the IDH-mutation status, higher than using cMRI and multi-shell dMRI separately (both 70 ± 7%). Similar findings were observed for predicting the 1p/19q-codeletion status, with the accuracy from combining cMRI and multi-shell dMRI (72 ± 4%) higher than from each modality used alone (cMRI: 65 ± 6%; multi-shell dMRI: 66 ± 9%). These findings remain when we considered all 146 patients for predicting the IDH status (combined: 81 ± 5% accuracy; cMRI: 74 ± 5%; multi-shell dMRI: 73 ± 6%) and for the diagnosis of the three molecular subtypes according to WHO 2021 (combined: 60 ± 5%; cMRI: 57 ± 8%; multi-shell dMRI: 56 ± 7%). Together, these findings suggest that combining cMRI and multi-shell dMRI can offer higher accuracy than using each modality alone for predicting the IDH and 1p/19q status and in diagnosing the three molecular subtypes with DL-based models.

Meningioma microstructure assessed by diffusion MRI: An investigation of the source of mean diffusivity and fractional anisotropy by quantitative histology

BACKGROUND

Mean diffusivity (MD) and fractional anisotropy (FA) from diffusion MRI (dMRI) have been associated with cell density and tissue anisotropy across tumors, but it is unknown whether these associations persist at the microscopic level.

PURPOSE

To quantify the degree to which cell density and anisotropy, as determined from histology, account for the intra-tumor variability of MD and FA in meningioma tumors. Furthermore, to clarify whether other histological features account for additional intra-tumor variability of dMRI parameters.

MATERIALS AND METHODS

We performed ex-vivo dMRI at 200 μm isotropic resolution and histological imaging of 16 excised meningioma tumor samples. Diffusion tensor imaging (DTI) was used to map MD and FA, as well as the in-plane FA (FA_IP). Histology images were analyzed in terms of cell nuclei density (CD) and structure anisotropy (SA; obtained from structure tensor analysis) and were used separately in a regression analysis to predict MD and FA_IP, respectively. A convolutional neural network (CNN) was also trained to predict the dMRI parameters from histology patches. The association between MRI and histology was analyzed in terms of out-of-sample (R²_OS) on the intra-tumor level and within-sample R² across tumors. Regions where the dMRI parameters were poorly predicted from histology were analyzed to identify features apart from CD and SA that could influence MD and FA_IP, respectively.

RESULTS

Cell density assessed by histology poorly explained intra-tumor variability of MD at the mesoscopic level (200 μm), as median R²_OS = 0.04 (interquartile range 0.01-0.26). Structure anisotropy explained more of the variation in FA_IP (median R²_OS = 0.31, 0.20-0.42). Samples with low R²_OS for FA_IP exhibited low variations throughout the samples and thus low explainable variability, however, this was not the case for MD. Across tumors, CD and SA were clearly associated with MD (R² = 0.60) and FA_IP (R² = 0.81), respectively. In 37% of the samples (6 out of 16), cell density did not explain intra-tumor variability of MD when compared to the degree explained by the CNN. Tumor vascularization, psammoma bodies, microcysts, and tissue cohesivity were associated with bias in MD prediction based solely on CD. Our results support that FA_IP is high in the presence of elongated and aligned cell structures, but low otherwise.

CONCLUSION

Cell density and structure anisotropy account for variability in MD and FA_IP across tumors but cell density does not explain MD variations within the tumor, which means that low or high values of MD locally may not always reflect high or low tumor cell density. Features beyond cell density need to be considered when interpreting MD.

Comparison of two types of microscopic diffusion anisotropy in mouse brain

Two distinct types of microscopic diffusion anisotropy (MA) are compared in brain for both normal control and transgenic (3xTg-AD) mice, which develop Alzheimer's disease pathology. The first type of MA is the commonly used microscopic fractional anisotropy (μFA), and the second is a new MA measure referred to as μFA'. These two MA parameters have different symmetry properties that are central to their physical interpretations. Specifically, μFA is invariant with respect to local rotations of compartmental diffusion tensors while μFA' is invariant with respect to global diffusion tensor deformations. A key distinction between μFA and μFA' is that μFA is affected by the same type of orientationally coherent diffusion anisotropy as the conventional fractional anisotropy (FA) while μFA' is not. Furthermore, μFA can be viewed as having independent contributions from FA and μFA', as is quantified by an equation relating all three anisotropies. The normal control and transgenic mice are studied at ages ranging from 2 to 15 months, with double diffusion encoding MRI being used to estimate μFA and μFA'. μFA and μFA' are nearly identical in low FA brain regions, but they show notable differences when FA is large. In particular, μFA and FA are found to be strongly correlated in the fimbria, but μFA' and FA are not. In addition, both μFA and μFA' are seen to increase with age in the corpus callosum and external capsule, and modest differences between normal control and transgenic mice are observed for μFA and μFA' in the corpus callosum and for μFA in the fimbria. The triad of FA, μFA, and μFA' is proposed as a useful combination of parameters for assessing diffusion anisotropy in brain.

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.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

The relationship between radiomics and pathomics in Glioblastoma patients: Preliminary results from a cross-scale association study

Glioblastoma multiforme (GBM) typically exhibits substantial intratumoral heterogeneity at both microscopic and radiological resolution scales. Diffusion Weighted Imaging (DWI) and dynamic contrast-enhanced (DCE) magnetic resonance imaging (MRI) are two functional MRI techniques that are commonly employed in clinic for the assessment of GBM tumor characteristics. This work presents initial results aiming at determining if radiomics features extracted from preoperative ADC maps and post-contrast T1 (T1C) images are associated with pathomic features arising from H&E digitized pathology images. 48 patients from the public available CPTAC-GBM database, for which both radiology and pathology images were available, were involved in the study. 91 radiomics features were extracted from ADC maps and post-contrast T1 images using PyRadiomics. 65 pathomic features were extracted from cell detection measurements from H&E images. Moreover, 91 features were extracted from cell density maps of H&E images at four different resolutions. Radiopathomic associations were evaluated by means of Spearman's correlation (ρ) and factor analysis. p values were adjusted for multiple correlations by using a false discovery rate adjustment. Significant cross-scale associations were identified between pathomics and ADC, both considering features (n = 186, 0.45 < ρ < 0.74 in absolute value) and factors (n = 5, 0.48 < ρ < 0.54 in absolute value). Significant but fewer ρ values were found concerning the association between pathomics and radiomics features (n = 53, 0.5 < ρ < 0.65 in absolute value) and factors (n = 2, ρ = 0.63 and ρ = 0.53 in absolute value). The results of this study suggest that cross-scale associations may exist between digital pathology and ADC and T1C imaging. This can be useful not only to improve the knowledge concerning GBM intratumoral heterogeneity, but also to strengthen the role of radiomics approach and its validation in clinical practice as "virtual biopsy", introducing new insights for omics integration toward a personalized medicine approach.

Reproducibility of the Standard Model of diffusion in white matter on clinical MRI systems

Estimating intra- and extra-axonal microstructure parameters, such as volume fractions and diffusivities, has been one of the major efforts in brain microstructure imaging with MRI. The Standard Model (SM) of diffusion in white matter has unified various modeling approaches based on impermeable narrow cylinders embedded in locally anisotropic extra-axonal space. However, estimating the SM parameters from a set of conventional diffusion MRI (dMRI) measurements is ill-conditioned. Multidimensional dMRI helps resolve the estimation degeneracies, but there remains a need for clinically feasible acquisitions that yield robust parameter maps. Here we find optimal multidimensional protocols by minimizing the mean-squared error of machine learning-based SM parameter estimates for two 3T scanners with corresponding gradient strengths of 40and80mT/m. We assess intra-scanner and inter-scanner repeatability for 15-minute optimal protocols by scanning 20 healthy volunteers twice on both scanners. The coefficients of variation all SM parameters except free water fraction are ≲10% voxelwise and 1-4% for their region-averaged values. As the achieved SM reproducibility outcomes are similar to those of conventional diffusion tensor imaging, our results enable robust in vivo mapping of white matter microstructure in neuroscience research and in the clinic.

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

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

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

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

Efficacy of MR diffusion kurtosis imaging for differentiating low-grade from high-grade glioma before surgery: A systematic review and meta-analysis

BACKGROUND

Accurate discrimination and diagnosis of low-grade glioma (LGG) and high-grade glioma (HGG) before surgery is clinically important because it affects the patient's outcome and guides the clinicians to select appropriate management. The aim of this study was to evaluate the diagnostic performance of diffusion kurtosis imaging (DKI) for differentiating LGG from HGG.

METHODS

A literature search of the PubMed, Web of Science, Cochrane Library and EMBASE databases was conducted up to December 15, 2020. Studies that evaluated the diagnostic performance of DKI for differentiating LGG from HGG were selected. Retrieved hits were evaluated by the Quality Assessment of Diagnostic Accuracy Studies-2 tool. Summary sensitivity and specificity were determined, and the data analysis was performed using Stata 14.0 and Review Manager 5.3.

RESULTS

Thirteen studies with 705 patients were included. The individual sensitivity and specificity of the 13 studies varied from 71% to 100% for sensitivity and 73-100% for specificity. The pooled sensitivity of DKI was 88% (95% confidence interval [CI], 83-91%), and the pooled specificity was 91% (95% CI, 86-95%). The area under the summary receiver operating characteristic curve was 0.93 (95% CI, 0.90-0.95). The pooled diagnostic odds ratio of DKI was 64.85 (95% CI 38.52-109.19). The levels of heterogeneity for sensitivity and specificity across the included studies were high (I² =66%) and mild (I² =47.04%), respectively. The multiple subgroup analyses were driven by DKI technique and study region.

CONCLUSIONS

DKI demonstrated a high diagnostic performance for differentiation of LGG from HGG.

A review on investigation of the basic contrast mechanism underlying multidimensional diffusion MRI in assessment of neurological disorders

INTRODUCTION

Multidimensional diffusion MRI (MDD MRI) is a novel diffusion technique that uses advanced gradient waveforms for microstructural tissue characterization to provide information about average rate, anisotropy and orientation of the diffusion and to disentangle the signal fraction from specific cell types i.e., elongated cells, isotropic cells and free water.

AIM

To review the diagnostic potential of MDD MRI in the clinical setting for microstructural tissue characterization in patients with neurological disorders to aid in patient care and treatment.

METHOD

A scoping review on the clinical applications of MDD MRI was conducted from original articles published in PubMed and Scopus from 2015 to 2021 using the keywords "Multidimensional diffusion MRI" OR "diffusion tensor distribution" OR "Tensor-Valued Diffusion" OR "b-tensor encoding" OR "microscopic diffusion anisotropy" OR "microscopic anisotropy" OR "microscopic fractional anisotropy" OR "double diffusion encoding" OR "triple diffusion encoding" OR "double pulsed field gradients" OR "double wave vector" OR "correlation tensor imaging" AND "brain" OR "axons".

RESULTS

Initially 145 articles were screened and after applying inclusion and exclusion criteria, nine articles were included in the final analysis. In most of these studies, microscopic diffusion anisotropy within the lesion showed deviation from the normal-appearing tissue.

CONCLUSION

Multidimensional diffusion MRI can provide better quantification and visualization of tissue microstructure than conventional diffusion MRI and can be used in the clinical setting for diagnosis of neurological disorders.

Separating Glioma Hyperintensities From White Matter by Diffusion-Weighted Imaging With Spherical Tensor Encoding

Background

Tumor-related hyperintensities in high b-value diffusion-weighted imaging (DWI) are radiologically important in the workup of gliomas. However, the white matter may also appear as hyperintense, which may conflate interpretation.

Purpose

To investigate whether DWI with spherical b-tensor encoding (STE) can be used to suppress white matter and enhance the conspicuity of glioma hyperintensities unrelated to white matter.

Materials and Methods

Twenty-five patients with a glioma tumor and at least one pathology-related hyperintensity on DWI underwent conventional MRI at 3 T. The DWI was performed both with linear and spherical tensor encoding (LTE-DWI and STE-DWI). The LTE-DWI here refers to the DWI obtained with conventional diffusion encoding and averaged across diffusion-encoding directions. Retrospectively, the differences in contrast between LTE-DWI and STE-DWI, obtained at a b-value of 2,000 s/mm², were evaluated by comparing hyperintensities and contralateral normal-appearing white matter (NAWM) both visually and quantitatively in terms of the signal intensity ratio (SIR) and contrast-to-noise ratio efficiency (CNR_eff).

Results

The spherical tensor encoding DWI was more effective than LTE-DWI at suppressing signals from white matter and improved conspicuity of pathology-related hyperintensities. The median SIR improved in all cases and on average by 28%. The median (interquartile range) SIR was 1.9 (1.6 – 2.1) for STE and 1.4 (1.3 – 1.7) for LTE, with a significant difference of 0.4 (0.3 –0.5) (p < 10^–4, paired U-test). In 40% of the patients, the SIR was above 2 for STE-DWI, but with LTE-DWI, the SIR was below 2 for all patients. The CNR_eff of STE-DWI was significantly higher than of LTE-DWI: 2.5 (2 – 3.5) vs. 2.3 (1.7 – 3.1), with a significant difference of 0.4 (−0.1 –0.6) (p < 10^–3, paired U-test). The STE improved CNR_eff in 70% of the cases. We illustrate the benefits of STE-DWI in three patients, where STE-DWI may facilitate an improved radiological description of tumor-related hyperintensity, including one case that could have been missed out if only LTE-DWI was inspected.

Conclusion

The contrast mechanism of high b-value STE-DWI results in a stronger suppression of white matter than conventional LTE-DWI, and may, therefore, be more sensitive and specific for assessment of glioma tumors and DWI-hyperintensities.

Multi-parametric PET/MRI for enhanced tumor characterization of patients with cervical cancer

Aim

The concept of personalized medicine has brought increased awareness to the importance of inter- and intra-tumor heterogeneity for cancer treatment. The aim of this study was to explore simultaneous multi-parametric PET/MRI prior to chemoradiotherapy for cervical cancer for characterization of tumors and tumor heterogeneity.

Methods

Ten patients with histologically proven primary cervical cancer were examined with multi-parametric 68Ga-NODAGA-E[c(RGDyK)]2-PET/MRI for radiation treatment planning after diagnostic 18F-FDG-PET/CT. Standardized uptake values (SUV) of RGD and FDG, diffusion weighted MRI and the derived apparent diffusion coefficient (ADC), and pharmacokinetic maps obtained from dynamic contrast-enhanced MRI with the Tofts model (iAUC60, Ktrans, ve, and kep) were included in the analysis. The spatial relation between functional imaging parameters in tumors was examined by a correlation analysis and joint histograms at the voxel level. The ability of multi-parametric imaging to identify tumor tissue classes was explored using an unsupervised 3D Gaussian mixture model-based cluster analysis.

Results

Functional MRI and PET of cervical cancers appeared heterogeneous both between patients and spatially within the tumors, and the relations between parameters varied strongly within the patient cohort. The strongest spatial correlation was observed between FDG uptake and ADC (median r =  − 0.7). There was moderate voxel-wise correlation between RGD and FDG uptake, and weak correlations between all other modalities. Distinct relations between the ADC and RGD uptake as well as the ADC and FDG uptake were apparent in joint histograms. A cluster analysis using the combination of ADC, FDG and RGD uptake suggested tissue classes which could potentially relate to tumor sub-volumes.

Conclusion

A multi-parametric PET/MRI examination of patients with cervical cancer integrated with treatment planning and including estimation of angiogenesis and glucose metabolism as well as MRI diffusion and perfusion parameters is feasible. A combined analysis of functional imaging parameters indicates a potential of multi-parametric PET/MRI to contribute to a better characterization of tumor heterogeneity than the modalities alone. However, the study is based on small patient numbers and further studies are needed prior to the future design of individually adapted treatment approaches based on multi-parametric functional imaging.

In vivo Correlation Tensor MRI reveals microscopic kurtosis in the human brain on a clinical 3T scanner

Diffusion MRI (dMRI) has become one of the most important imaging modalities for noninvasively probing tissue microstructure. Diffusional Kurtosis MRI (DKI) quantifies the degree of non-gaussian diffusion, which in turn has been shown to increase sensitivity towards, e.g., disease and orientation mapping in neural tissue. However, the specificity of DKI is limited as different sources can contribute to the total intravoxel diffusional kurtosis, including: variance in diffusion tensor magnitudes (K_iso), variance due to diffusion anisotropy (K_aniso), and microscopic kurtosis (μK) related to restricted diffusion, microstructural disorder, and/or exchange. Interestingly, μK is typically ignored in diffusion MRI signal modeling as it is assumed to be negligible in neural tissues. However, recently, Correlation Tensor MRI (CTI) based on Double-Diffusion-Encoding (DDE) was introduced for kurtosis source separation, revealing non negligible μK in preclinical imaging. Here, we implemented CTI for the first time on a clinical 3T scanner and investigated the sources of total kurtosis in healthy subjects. A robust framework for kurtosis source separation in humans is introduced, followed by estimation of μK (and the other kurtosis sources) in the healthy brain. Using this clinical CTI approach, we find that μK significantly contributes to total diffusional kurtosis both in gray and white matter tissue but, as expected, not in the ventricles. The first μK maps of the human brain are presented, revealing that the spatial distribution of μK provides a unique source of contrast, appearing different from isotropic and anisotropic kurtosis counterparts. Moreover, group average templates of these kurtosis sources have been generated for the first time, which corroborated our findings at the underlying individual-level maps. We further show that the common practice of ignoring μK and assuming the multiple gaussian component approximation for kurtosis source estimation introduces significant bias in the estimation of other kurtosis sources and, perhaps even worse, compromises their interpretation. Finally, a twofold acceleration of CTI is discussed in the context of potential future clinical applications. We conclude that CTI has much potential for future in vivo microstructural characterizations in healthy and pathological tissue.

Histogram analysis of tensor-valued diffusion MRI in meningiomas: Relation to consistency, histological grade and type

•

Tensor-valued dMRI facilitates prediction of meningioma consistency, grade and type.

•

Tensor-valued dMRI corroborates findings of diffusion tensor and kurtosis imaging.

•

MK and MK_A is associated with firm and MD with variable meningioma consistency.

•

Variability of MK_I in the vicinity of the tumor is associated with meningioma grade.

•

MK_{A 50} and MK_{I 50} separates psammomatous meningiomas from other meningioma types.

Background

Preoperative radiological assessment of meningioma characteristics is of value for pre- and post-operative patient management, counselling, and surgical approach.

Purpose

To investigate whether tensor-valued diffusion MRI can add to the preoperative prediction of meningioma consistency, grade and type.

Materials and methods

30 patients with intracranial meningiomas (22 WHO grade I, 8 WHO grade II) underwent MRI prior to surgery. Diffusion MRI was performed with linear and spherical b-tensors with b-values up to 2000 s/mm². The data were used to estimate mean diffusivity (MD), fractional anisotropy (FA), mean kurtosis (MK) and its components—the anisotropic and isotropic kurtoses (MK_A and MK_I). Meningioma consistency was estimated for 16 patients during resection based on ultrasonic aspiration intensity, ease of resection with instrumentation or suction. Grade and type were determined by histopathological analysis. The relation between consistency, grade and type and dMRI parameters was analyzed inside the tumor (“whole-tumor”) and within brain tissue in the immediate periphery outside the tumor (“rim”) by histogram analysis.

Results

Lower 10th percentiles of MK and MK_A in the whole-tumor were associated with firm consistency compared with pooled soft and variable consistency (n = 7 vs 9; U test, p = 0.02 for MK_{A 10} and p = 0.04 for MK₁₀) and lower 10th percentile of MD with variable against soft and firm (n = 5 vs 11; U test, p = 0.02). Higher standard deviation of MK_I in the rim was associated with lower grade (n = 22 vs 8; U test, p = 0.04) and in the MK_I maps we observed elevated rim-like structure that could be associated with grade. Higher median MK_A and lower median MK_I distinguished psammomatous type from other pooled meningioma types (n = 5 vs 25; U test; p = 0.03 for MK_{A 50} and p = 0.03 and p = 0.04 for MK_{I 50}).

Conclusion

Parameters from tensor-valued dMRI can facilitate prediction of consistency, grade and type.

Clinical experience of tensor-valued diffusion encoding for microstructure imaging by diffusional variance decomposition in patients with breast cancer

BACKGROUND

Diffusion-weighted imaging plays a key role in magnetic resonance imaging (MRI) of breast tumors. However, it remains unclear how to interpret single diffusion encoding with respect to its link with tissue microstructure. The purpose of this retrospective cross-sectional study was to use tensor-valued diffusion encoding to investigate the underlying microstructure of invasive ductal carcinoma (IDC) and evaluate its potential value in a clinical setting.

METHODS

We retrospectively reviewed biopsy-proven breast cancer patients who underwent preoperative breast MRI examination from July 2020 to March 2021. We reviewed the MRI of 29 patients with 30 IDCs, including analysis by diffusional variance decomposition enabled by tensor-valued diffusion encoding. The diffusion parameters of mean diffusivity (MD), total mean kurtosis (MKT), anisotropic mean kurtosis (MKA), isotropic mean kurtosis (MKI), macroscopic fractional anisotropy (FA), and microscopic fractional anisotropy (µFA) were estimated. The parameter differences were compared between IDC and normal fibroglandular breast tissue (FGBT), as well as the association between the diffusion parameters and histopathologic items.

RESULTS

The mean value of MD in IDCs was significantly lower than that of normal FGBT (1.07±0.27 vs. 1.34±0.29, P<0.001); however, MKT, MKA, MKI, FA, and µFA were significantly higher (P<0.005). Among all the diffusion parameters, MKI was positively correlated with the tumor size on both MRI and pathological specimen (rs=0.38, P<0.05 vs. rs=0.54, P<0.01), whereas MKT had a positive correlation with the tumor size in the pathological specimen only (rs=0.47, P<0.02). In addition, the lymph node (LN) metastasis group had significantly higher MKT, MKA, and µFA compared to the metastasis negative group (P<0.05).

CONCLUSIONS

Tensor-valued diffusion encoding enables a useful non-invasive method for characterizing breast cancers with information on tissue microstructures. Particularly, µFA could be a potential imaging biomarker for evaluating breast cancers prior to surgery or chemotherapy.

Mapping the Human Connectome using Diffusion MRI at 300 mT/m Gradient Strength: Methodological Advances and Scientific Impact

Tremendous efforts have been made in the last decade to advance cutting-edge MRI technology in pursuit of mapping structural connectivity in the living human brain with unprecedented sensitivity and speed. The first Connectom 3T MRI scanner equipped with a 300 mT/m whole-body gradient system was installed at the Massachusetts General Hospital in 2011 and was specifically constructed as part of the Human Connectome Project. Since that time, numerous technological advances have been made to enable the broader use of the Connectom high gradient system for diffusion tractography and tissue microstructure studies and leverage its unique advantages and sensitivity to resolving macroscopic and microscopic structural information in neural tissue for clinical and neuroscientific studies. The goal of this review article is to summarize the technical developments that have emerged in the last decade to support and promote large-scale and scientific studies of the human brain using the Connectom scanner. We provide a brief historical perspective on the development of Connectom gradient technology and the efforts that led to the installation of three other Connectom 3T MRI scanners worldwide - one in the United Kingdom in Cardiff, Wales, another in Continental Europe in Leipzig, Germany, and the latest in Asia in Shanghai, China. We summarize the key developments in gradient hardware and image acquisition technology that have formed the backbone of Connectom-related research efforts, including the rich array of high-sensitivity receiver coils, pulse sequences, image artifact correction strategies and data preprocessing methods needed to optimize the quality of high-gradient strength dMRI data for subsequent analyses. Finally, we review the scientific impact of the Connectom MRI scanner, including advances in diffusion tractography, tissue microstructural imaging, ex vivo validation, and clinical investigations that have been enabled by Connectom technology. We conclude with brief insights into the unique value of strong gradients for dMRI and where the field is headed in the coming years.

Diffusion MRI signal cumulants and hepatocyte microstructure at fixed diffusion time: Insights from simulations, 9.4T imaging, and histology

Purpose

Relationships between diffusion‐weighted MRI signals and hepatocyte microstructure were investigated to inform liver diffusion MRI modeling, focusing on the following question: Can cell size and diffusivity be estimated at fixed diffusion time, realistic SNR, and negligible contribution from extracellular/extravascular water and exchange?

Methods

Monte Carlo simulations were performed within synthetic hepatocytes for varying cell size/diffusivity L/D0, and clinical protocols (single diffusion encoding; maximum b‐value: {1000, 1500, 2000} s/mm²; 5 unique gradient duration/separation pairs; SNR = {∞, 100, 80, 40, 20}), accounting for heterogeneity in (D0,L) and perfusion contamination. Diffusion (D) and kurtosis (K) coefficients were calculated, and relationships between (D0,L) and (D,K) were visualized. Functions mapping (D,K) to (D0,L) were computed to predict unseen (D0,L) values, tested for their ability to classify discrete cell‐size contrasts, and deployed on 9.4T ex vivo MRI‐histology data of fixed mouse livers

Results

Relationships between (D,K) and (D0,L) are complex and depend on the diffusion encoding. Functions mapping D,K to (D0,L) captures salient characteristics of D0(D,K) and L(D,K) dependencies. Mappings are not always accurate, but they enable just under 70% accuracy in a three‐class cell‐size classification task (for SNR = 20, bmax = 1500 s/mm², δ = 20 ms, and Δ = 75 ms). MRI detects cell‐size contrasts in the mouse livers that are confirmed by histology, but overestimates the largest cell sizes.

Conclusion

Salient information about liver cell size and diffusivity may be retrieved from minimal diffusion encodings at fixed diffusion time, in experimental conditions and pathological scenarios for which extracellular, extravascular water and exchange are negligible.

Potential Pitfalls of Using Fractional Anisotropy, Axial Diffusivity, and Radial Diffusivity as Biomarkers of Cerebral White Matter Microstructure

Fractional anisotropy (FA), axial diffusivity (AD), and radial diffusivity (RD) are commonly used as MRI biomarkers of white matter microstructure in diffusion MRI studies of neurodevelopment, brain aging, and neurologic injury/disease. Some of the more frequent practices include performing voxel-wise or region-based analyses of these measures to cross-sectionally compare individuals or groups, longitudinally assess individuals or groups, and/or correlate with demographic, behavioral or clinical variables. However, it is now widely recognized that the majority of cerebral white matter voxels contain multiple fiber populations with different trajectories, which renders these metrics highly sensitive to the relative volume fractions of the various fiber populations, the microstructural integrity of each constituent fiber population, and the interaction between these factors. Many diffusion imaging experts are aware of these limitations and now generally avoid using FA, AD or RD (at least in isolation) to draw strong reverse inferences about white matter microstructure, but based on the continued application and interpretation of these metrics in the broader biomedical/neuroscience literature, it appears that this has perhaps not yet become common knowledge among diffusion imaging end-users. Therefore, this paper will briefly discuss the complex biophysical underpinnings of these measures in the context of crossing fibers, provide some intuitive “thought experiments” to highlight how conventional interpretations can lead to incorrect conclusions, and suggest that future studies refrain from using (over-interpreting) FA, AD, and RD values as standalone biomarkers of cerebral white matter microstructure.

Correlation Tensor MRI deciphers underlying kurtosis sources in stroke

Noninvasively detecting and characterizing modulations in cellular scale micro-architecture remains a desideratum for contemporary neuroimaging. Diffusion MRI (dMRI) has become the mainstay methodology for probing microstructure, and, in ischemia, its contrasts have revolutionized stroke management. Diffusion kurtosis imaging (DKI) has been shown to significantly enhance the sensitivity of stroke detection compared to its diffusion tensor imaging (DTI) counterparts. However, the interpretation of DKI remains ambiguous as its contrast may arise from competing kurtosis sources related to the anisotropy of tissue components, diffusivity variance across components, and microscopic kurtosis (e.g., arising from cross-sectional variance, structural disorder, and restriction). Resolving these sources may be fundamental for developing more specific imaging techniques for stroke management, prognosis, and understanding its pathophysiology. In this study, we apply Correlation Tensor MRI (CTI) - a double diffusion encoding (DDE) methodology recently introduced for deciphering kurtosis sources based on the unique information captured in DDE's diffusion correlation tensors - to investigate the underpinnings of kurtosis measurements in acute ischemic lesions. Simulations for the different kurtosis sources revealed specific signatures for cross-sectional variance (representing neurite beading), edema, and cell swelling. Ex vivo CTI experiments at 16.4 T were then performed in an experimental photothrombotic stroke model 3 h post-stroke (N = 10), and successfully separated anisotropic, isotropic, and microscopic non-Gaussian diffusion sources in the ischemic lesions. Each of these kurtosis sources provided unique contrasts in the stroked area. Particularly, microscopic kurtosis was shown to be a primary "driver" of total kurtosis upon ischemia; its large increases, coupled with decreases in anisotropic kurtosis, are consistent with the expected elevation in cross-sectional variance, likely linked to beading effects in small objects such as neurites. In vivo experiments at 9.4 T at the same time point (3 h post ischemia, N = 5) demonstrated the stability and relevance of the findings and showed that fixation is not a dominant confounder in our findings. In future studies, the different CTI contrasts may be useful to address current limitations of stroke imaging, e.g., penumbra characterization, distinguishing lesion progression form tissue recovery, and elucidating pathophysiological correlates.

Multi-tissue spherical deconvolution of tensor-valued diffusion MRI

Multi-tissue constrained spherical deconvolution (MT-CSD) leverages the characteristic b-value dependency of each tissue type to estimate both the apparent tissue densities and the white matter fiber orientation distribution function from diffusion MRI data. In this work, we generalize MT-CSD to tensor-valued diffusion encoding with arbitrary b-tensor shapes. This enables the use of data encoded with mixed b-tensors, rather than being limited to the subset of linear (conventional) b-tensors. Using the complete set of data, including all b-tensor shapes, provides a categorical improvement in the estimation of apparent tissue densities, fiber ODF, and resulting tractography. Furthermore, we demonstrate that including multiple b-tensor shapes in the analysis provides improved contrast between tissue types, in particular between gray matter and white matter. We also show that our approach provides high-quality apparent tissue density maps and high-quality fiber tracking from data, even with sparse sampling across b-tensors that yield whole-brain coverage at 2 mm isotropic resolution in approximately 5:15 min.

A Pilot Study of Multidimensional Diffusion MRI for Assessment of Tissue Heterogeneity in Prostate Cancer

OBJECTIVES

The objectives of this exploratory study were to investigate the feasibility of multidimensional diffusion magnetic resonance imaging (MddMRI) in assessing diffusion heterogeneity at both a macroscopic and microscopic level in prostate cancer (PCa).

MATERIALS AND METHODS

Informed consent was obtained from 46 subjects who underwent 3.0-T prostate multiparametric MRI, complemented with a prototype spin echo-based MddMRI sequence in this institutional review board-approved study. Prostate cancer tumors and comparative normal tissue from each patient were contoured on both apparent diffusion coefficient and MddMRI-derived mean diffusivity (MD) maps (from which microscopic diffusion heterogeneity [MKi] and microscopic diffusion anisotropy were derived) using 3D Slicer. The discriminative ability of MddMRI-derived parameters to differentiate PCa from normal tissue was determined using the Friedman test. To determine if tumor diffusion heterogeneity is similar on macroscopic and microscopic scales, the linear association between SD of MD and mean MKi was estimated using robust regression (bisquare weighting). Hypothesis testing was 2 tailed; P values less than 0.05 were considered statistically significant.

RESULTS

All MddMRI-derived parameters could distinguish tumor from normal tissue in the fixed-effects analysis (P < 0.0001). Tumor MKi was higher (P < 0.05) compared with normal tissue (median, 0.40; interquartile range, 0.29-0.52 vs 0.20-0.18; 0.25), as was tumor microscopic diffusion anisotropy (0.55; 0.36-0.81 vs 0.20-0.15; 0.28). The MKi could not be predicted (no significant association) by SD of MD. There was a significant correlation between tumor volume and SD of MD (R2 = 0.50, slope = 0.008 μm2/ms per millimeter, P < 0.001) but not between tumor volume and MKi.

CONCLUSIONS

This explorative study demonstrates that MddMRI provides novel information on MKi and microscopic anisotropy, which differ from measures at the macroscopic level. MddMRI has the potential to characterize tumor tissue heterogeneity at different spatial scales.

Current emerging MRI tools for radionecrosis and pseudoprogression diagnosis

PURPOSE OF REVIEW

This review aims to cover current MRI techniques for assessing treatment response in brain tumors, with a focus on radio-induced lesions.

RECENT FINDINGS

Pseudoprogression and radionecrosis are common radiological entities after brain tumor irradiation and are difficult to distinguish from real progression, with major consequences on daily patient care. To date, shortcomings of conventional MRI have been largely recognized but morphological sequences are still used in official response assessment criteria. Several complementary advanced techniques have been proposed but none of them have been validated, hampering their clinical use. Among advanced MRI, brain perfusion measures increase diagnostic accuracy, especially when added with spectroscopy and susceptibility-weighted imaging. However, lack of reproducibility, because of several hard-to-control variables, is still a major limitation for their standardization in routine protocols. Amide Proton Transfer is an emerging molecular imaging technique that promises to offer new metrics by indirectly quantifying intracellular mobile proteins and peptide concentration. Preliminary studies suggest that this noncontrast sequence may add key biomarkers in tumor evaluation, especially in posttherapeutic settings.

SUMMARY

Benefits and pitfalls of conventional and advanced imaging on posttreatment assessment are discussed and the potential added value of APT in this clinicoradiological evolving scenario is introduced.

Diffusion-Weighted Imaging: Recent Advances and Applications

Quantitative diffusion imaging techniques enable the characterization of tissue microstructural properties of the human brain "in vivo", and are widely used in neuroscientific and clinical contexts. In this review, we present the basic physical principles behind diffusion imaging and provide an overview of the current diffusion techniques, including standard and advanced techniques as well as their main clinical applications. Standard diffusion tensor imaging (DTI) offers sensitivity to changes in microstructure due to diseases and enables the characterization of single fiber distributions within a voxel as well as diffusion anisotropy. Nonetheless, its inability to represent complex intravoxel fiber topologies and the limited biological specificity of its metrics motivated the development of several advanced diffusion MRI techniques. For example, high-angular resolution diffusion imaging (HARDI) techniques enabled the characterization of fiber crossing areas and other complex fiber topologies in a single voxel and supported the development of higher-order signal representations aiming to decompose the diffusion MRI signal into distinct microstructure compartments. Biophysical models, often known by their acronym (e.g., CHARMED, WMTI, NODDI, DBSI, DIAMOND) contributed to capture the diffusion properties from each of such tissue compartments, enabling the computation of voxel-wise maps of axonal density and/or morphology that hold promise as clinically viable biomarkers in several neurological and neuroscientific applications; for example, to quantify tissue alterations due to disease or healthy processes. Current challenges and limitations of state-of-the-art models are discussed, including validation efforts. Finally, novel diffusion encoding approaches (e.g., b-tensor or double diffusion encoding) may increase the biological specificity of diffusion metrics towards intra-voxel diffusion heterogeneity in clinical settings, holding promise in neurological applications.

Comparative analysis of signal models for microscopic fractional anisotropy estimation using q-space trajectory encoding

Microscopic diffusion anisotropy imaging using diffusion-weighted MRI and multidimensional diffusion encoding is a promising method for quantifying clinically and scientifically relevant microstructural properties of neural tissue. Several methods for estimating microscopic fractional anisotropy (µFA), a normalized measure of microscopic diffusion anisotropy, have been introduced but the differences between the methods have received little attention thus far. In this study, the accuracy and precision of µFA estimation using q-space trajectory encoding and different signal models were assessed using imaging experiments and simulations. Three healthy volunteers and a microfibre phantom were imaged with five non-zero b-values and gradient waveforms encoding linear and spherical b-tensors. Since the ground-truth µFA was unknown in the imaging experiments, Monte Carlo random walk simulations were performed using axon-mimicking fibres for which the ground truth was known. Furthermore, parameter bias due to time-dependent diffusion was quantified by repeating the simulations with tuned waveforms, which have similar power spectra, and with triple diffusion encoding, which, unlike q-space trajectory encoding, is not based on the assumption of time-independent diffusion. The truncated cumulant expansion of the powder-averaged signal, gamma-distributed diffusivities assumption, and q-space trajectory imaging, a generalization of the truncated cumulant expansion to individual signals, were used to estimate µFA. The gamma-distributed diffusivities assumption consistently resulted in greater µFA values than the second order cumulant expansion, 0.1 greater when averaged over the whole brain. In the simulations, the generalized cumulant expansion provided the most accurate estimates. Importantly, although time-dependent diffusion caused significant overestimation of µFA using all the studied methods, the simulations suggest that the resulting bias in µFA is less than 0.1 in human white matter.

Evidence for microscopic kurtosis in neural tissue revealed by correlation tensor MRI

PURPOSE

The impact of microscopic diffusional kurtosis (µK), arising from restricted diffusion and/or structural disorder, remains a controversial issue in contemporary diffusion MRI (dMRI). Recently, correlation tensor imaging (CTI) was introduced to disentangle the sources contributing to diffusional kurtosis, without relying on a-priori multi-gaussian component (MGC) or other microstructural assumptions. Here, we investigated µK in in vivo rat brains and assessed its impact on state-of-the-art methods ignoring µK.

THEORY AND METHODS

CTI harnesses double diffusion encoding (DDE) experiments, which were here improved for speed and minimal bias using four different sets of acquisition parameters. The robustness of the improved CTI protocol was assessed via simulations. In vivo CTI acquisitions were performed in healthy rat brains using a 9.4T pre-clinical scanner equipped with a cryogenic coil, and targeted the estimation of µK, anisotropic kurtosis, and isotropic kurtosis.

RESULTS

The improved CTI acquisition scheme substantially reduces scan time and importantly, also minimizes higher-order-term biases, thus enabling robust µK estimation, alongside K_aniso and K_iso metrics. Our CTI experiments revealed positive µK both in white and gray matter of the rat brain in vivo; µK is the dominant kurtosis source in healthy gray matter tissue. The non-negligible µK substantially were found to bias prior MGC analyses of K_iso and K_aniso .

CONCLUSIONS

Correlation Tensor MRI offers a more accurate and robust characterization of kurtosis sources than its predecessors. µK is non-negligible in vivo in healthy white and gray matter tissues and could be an important biomarker for future studies. Our findings thus have both theoretical and practical implications for future dMRI research.

Mapping complex cell morphology in the grey matter with double diffusion encoding MR: A simulation study

This paper investigates the impact of cell body (namely soma) size and branching of cellular projections on diffusion MR imaging (dMRI) and spectroscopy (dMRS) signals for both standard single diffusion encoding (SDE) and more advanced double diffusion encoding (DDE) measurements using numerical simulations. The aim is to investigate the ability of dMRI/dMRS to characterize the complex morphology of brain cells focusing on these two distinctive features of brain grey matter. To this end, we employ a recently developed computational framework to create three dimensional meshes of neuron-like structures for Monte Carlo simulations, using diffusion coefficients typical of water and brain metabolites. Modelling the cellular structure as realistically connected spherical soma and cylindrical cellular projections, we cover a wide range of combinations of sphere radii and branching order of cellular projections, characteristic of various grey matter cells. We assess the impact of spherical soma size and branching order on the b-value dependence of the SDE signal as well as the time dependence of the mean diffusivity (MD) and mean kurtosis (MK). Moreover, we also assess the impact of spherical soma size and branching order on the angular modulation of DDE signal at different mixing times, together with the mixing time dependence of the apparent microscopic anisotropy (μA), a promising contrast derived from DDE measurements. The SDE results show that spherical soma size has a measurable impact on both the b-value dependence of the SDE signal and the MD and MK diffusion time dependence for both water and metabolites. On the other hand, we show that branching order has little impact on either, especially for water. In contrast, the DDE results show that spherical soma size has a measurable impact on the DDE signal's angular modulation at short mixing times and the branching order of cellular projections significantly impacts the mixing time dependence of the DDE signal's angular modulation as well as of the derived μA, for both water and metabolites. Our results confirm that SDE based techniques may be sensitive to spherical soma size, and most importantly, show for the first time that DDE measurements may be more sensitive to the dendritic tree complexity (as parametrized by the branching order of cellular projections), paving the way for new ways of characterizing grey matter morphology, non-invasively using dMRS and potentially dMRI.

Diffusional Kurtosis Imaging in the Diffusion Imaging in Python Project

Diffusion-weighted magnetic resonance imaging (dMRI) measurements and models provide information about brain connectivity and are sensitive to the physical properties of tissue microstructure. Diffusional Kurtosis Imaging (DKI) quantifies the degree of non-Gaussian diffusion in biological tissue from dMRI. These estimates are of interest because they were shown to be more sensitive to microstructural alterations in health and diseases than measures based on the total anisotropy of diffusion which are highly confounded by tissue dispersion and fiber crossings. In this work, we implemented DKI in the Diffusion in Python (DIPY) project-a large collaborative open-source project which aims to provide well-tested, well-documented and comprehensive implementation of different dMRI techniques. We demonstrate the functionality of our methods in numerical simulations with known ground truth parameters and in openly available datasets. A particular strength of our DKI implementations is that it pursues several extensions of the model that connect it explicitly with microstructural models and the reconstruction of 3D white matter fiber bundles (tractography). For instance, our implementations include DKI-based microstructural models that allow the estimation of biophysical parameters, such as axonal water fraction. Moreover, we illustrate how DKI provides more general characterization of non-Gaussian diffusion compatible with complex white matter fiber architectures and gray matter, and we include a novel mean kurtosis index that is invariant to the confounding effects due to tissue dispersion. In summary, DKI in DIPY provides a well-tested, well-documented and comprehensive reference implementation for DKI. It provides a platform for wider use of DKI in research on brain disorders and in cognitive neuroscience.

Survival Prediction Analysis in Glioblastoma With Diffusion Kurtosis Imaging

SIMPLE SUMMARY

Glioblastoma (GBM) is the most common and aggressive primary brain tumor. Diffusion kurtosis imaging (DKI) has characterized non-Gaussian diffusion behaviors in brain normal tissue and gliomas, but there are very limited efforts in investigating treatment responses of kurtosis in GBM. This study aimed to investigate whether any parameter derived from the DKI is a significant predictor of overall survival (OS). We found that the large mean, 80 and 90 percentile kurtosis values in the contrast enhanced gross tumor volume (Gd-GTV) on post-Gd T1-weighted images pre-RT were significantly associated with reduced OS. In the multivariate Cox model, the mean kurtosis Gd-GTV pre-RT after considering effects of age, extent of surgery, and methylation were significant predictors of OS. In addition, the 80 and 90 percentile kurtosis values in Gd-GTV post RT were significantly associated with progression free survival (PFS). The DKI model demonstrates the potential to predict outcomes in the patients with GBM.

PURPOSE

Non-Gaussian diffusion behaviors in gliomas have been characterized by diffusion kurtosis imaging (DKI). But there are very limited efforts in investigating the kurtosis in glioblastoma (GBM) and its prognostic and predictive values. This study aimed to investigate whether any of the diffusion kurtosis parameters derived from DKI is a significant predictor of overall survival.

METHODS AND MATERIALS

Thirty-three patients with GBM had pre-radiation therapy (RT) and mid-RT diffusion weighted (DW) images. Kurtosis and diffusion coefficient (DC) values in the contrast enhanced gross tumor volume (Gd-GTV) on post-Gd T1 weighted images pre-RT and mid-RT were calculated. Univariate and multivariate Cox models were used to evaluate the DKI parameters and clinical factors for prediction of OS and PFS.

RESULTS

The large mean kurtosis values in the Gd-GTV pre-RT were significantly associated with reduced OS (p = 0.02), but the values at mid-RT were not (p > 0.8). In the multivariate Cox model, the mean kurtosis in the Gd-GTV pre-RT (p = 0.009) was still a significant predictor of OS after adjusting effects of age, O6-Methylguanine-DNA Methyl transferase (MGMT) methylation and extent of resection. In Gd-GTV post-RT, 80 and 90 percentile kurtosis values were significant predictors (p ≤ 0.05) for progression free survival (PFS).

CONCLUSION

The DKI model demonstrates the potential to predict OS and PFS in the patients with GBM. Further development and histopathological validation of the DKI model will warrant its role in clinical management of GBM.

On the potential for mapping apparent neural soma density via a clinically viable diffusion MRI protocol

•

B-tensor encoding enables estimation of spherical cellular structures in the brain.

•

Spherical compartments may provide markers for apparent neural soma density.

•

Model parameters can be estimated in a fast and robust way using deep learning.

•

Practical acquisition times are achievable on widely available clinical scanners.

Diffusion MRI is a valuable tool for probing tissue microstructure in the brain noninvasively. Today, model-based techniques are widely available and used for white matter characterisation where their development is relatively mature. Conversely, tissue modelling in grey matter is more challenging, and no generally accepted models exist. With advances in measurement technology and modelling efforts, a clinically viable technique that reveals salient features of grey matter microstructure, such as the density of quasi-spherical cell bodies and quasi-cylindrical cell projections, is an exciting prospect. As a step towards capturing the microscopic architecture of grey matter in clinically feasible settings, this work uses a biophysical model that is designed to disentangle the diffusion signatures of spherical and cylindrical structures in the presence of orientation heterogeneity, and takes advantage of B-tensor encoding measurements, which provide additional sensitivity compared to standard single diffusion encoding sequences. For the fast and robust estimation of microstructural parameters, we leverage recent advances in machine learning and replace conventional fitting techniques with an artificial neural network that fits complex biophysical models within seconds. Our results demonstrate apparent markers of spherical and cylindrical geometries in healthy human subjects, and in particular an increased volume fraction of spherical compartments in grey matter compared to white matter. We evaluate the extent to which spherical and cylindrical geometries may be interpreted as correlates of neural soma and neural projections, respectively, and quantify parameter estimation errors in the presence of various departures from the modelling assumptions. While further work is necessary to translate the ideas presented in this work to the clinic, we suggest that biomarkers focussing on quasi-spherical cellular geometries may be valuable for the enhanced assessment of neurodevelopmental disorders and neurodegenerative diseases.

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.

Mapping prostatic microscopic anisotropy using linear and spherical b-tensor encoding: A preliminary study

Purpose:

Tensor-valued diffusion encoding provides more specific information than conventional diffusion-weighted imaging (DWI), but has mainly been applied in neuroimaging studies. This study aimed to assess its potential for the imaging of prostate cancer (PCa).

Methods:

Seventeen patients with histologically proven PCa were enrolled. DWI of the prostate was performed with linear and spherical tensor encoding using a maximal b-value of 1.5 ms/μm² and a voxel size of 3 × 3 × 4 mm³. The gamma-distribution model was used to estimate the mean diffusivity (MD), the isotropic kurtosis (MK_I), and the anisotropic kurtosis (MK_A). Regions of interest were placed in MR-defined cancerous tissues, as well as in apparently healthy tissues in the peripheral and transitional zones (PZs and TZs).

Results:

DWI with linear and spherical encoding yielded different image contrasts at high b-values, which enabled the estimation of MK_A and MK_I. Compared with healthy tissue (PZs and TZs combined) the cancers displayed a significantly lower MD (P < .05), higher MK_I (P < 10⁻⁵), and lower MK_A (P < .05). Compared with the TZ, tissue in the PZ showed lower MD (P < 10⁻³) and higher MK_A (P < 10⁻³). No significant differences were found between cancers of different Gleason scores, possibly because of the limited sample size.

Conclusion:

Tensor-valued diffusion encoding enabled mapping of MK_A and MK_I in the prostate. The elevated MK_I in PCa compared with normal tissues suggests an elevated heterogeneity in the cancers. Increased in-plane resolution could improve tumor delineation in future studies.

Q-space trajectory imaging with positivity constraints (QTI+)

Q-space trajectory imaging (QTI) enables the estimation of useful scalar measures indicative of the local tissue structure. This is accomplished by employing generalized gradient waveforms for diffusion sensitization alongside a diffusion tensor distribution (DTD) model. The first two moments of the underlying DTD are made available by acquisitions at low diffusion sensitivity (b-values). Here, we show that three independent conditions have to be fulfilled by the mean and covariance tensors associated with distributions of symmetric positive semidefinite tensors. We introduce an estimation framework utilizing semi-definite programming (SDP) to guarantee that these conditions are met. Applying the framework on simulated signal profiles for diffusion tensors distributed according to non-central Wishart distributions demonstrates the improved noise resilience of QTI+ over the commonly employed estimation methods. Our findings on a human brain data set also reveal pronounced improvements, especially so for acquisition protocols featuring few number of volumes. Our method’s robustness to noise is expected to not only improve the accuracy of the estimates, but also enable a meaningful interpretation of contrast in the derived scalar maps. The technique’s performance on shorter acquisitions could make it feasible in routine clinical practice.