Quantification of microscopic diffusion anisotropy disentangles effects of orientation dispersion from microstructure: applications in healthy volunteers and in brain tumors

Linear rotationally invariant kurtosis measures from double diffusion encoding MRI

PURPOSE

To characterize the complete set of linear rotationally invariant kurtosis measures provided by double diffusion encoding (DDE) MRI, show their utility in distinguishing different types of multiple Gaussian compartment (MGC) models, and demonstrate simplified acquisition and analysis schemes for their estimation.

THEORY AND METHODS

The lowest order novel information obtainable with DDE MRI can be encapsulated in a six-dimensional kurtosis tensor. The most basic DDE MRI kurtosis measures are rotational invariants that are linear in this tensor while depending on no other physical quantities. We identify four such invariants and show that any others must be linear combinations of these. The invariants are applied to classify MGC models according to whether they include microscopic anisotropy or intercompartmental water exchange. In addition, they are used to investigate the effect of exchange on estimates of the microscopic fractional anisotropy (μFA). Simplified acquisition and analysis schemes for the invariants are proposed and demonstrated with human brain data obtained at 3 T.

RESULTS

For the considered brain regions, the kurtosis invariants are found to be largely consistent with MGC models having microscopic anisotropy. They also indicate that water exchange in gray matter may affect estimates of μFA.

CONCLUSION

The kurtosis measures can classify MGC models according to whether they have microscopic anisotropy or water exchange, and they can be estimated with simple acquisition and analysis schemes. Measurements of the invariants in brain support the validity of MGC models with microscopic anisotropy and the importance of water exchange for modeling diffusion in gray matter.

Basis function compression for field probe monitoring

PURPOSE

Field monitoring using field probes allows for accurate measurement of magnetic field perturbations, such as from eddy currents, during MRI scanning. However, errors may result when the spatial variation of the fields is not well-described by the conventionally used spherical harmonics model that has the maximum order constrained by the number of probes. The objective of this work was to develop and validate a field monitoring approach that compresses higher order spherical harmonics into a smaller set of new basis functions that can be characterized using fewer probes.

METHODS

Field monitoring of acquisitions was repeated with probes in different locations. High-order field dynamics were computed from this "calibration" data assembled from provided scans, from which compression matrices could be devised using principal component analysis. Compression matrices were then used to fit field dynamics using "compressed" basis functions with data from 16 probes, which were then used in image reconstruction. Performance was evaluated by assessing the accuracy of computed field dynamics as well as in vivo image quality. Technique generalizability was also assessed by using various acquisition and diffusion encoding strategies in the calibration.

RESULTS

Qualitative and quantitative improvements in accuracy were observed when using the proposed fitting method compared to the conventional approach. However, compression effectiveness was influenced by the probe quantity and arrangement, and the specific acquisition data included in the calibration.

CONCLUSION

The ability to tailor basis functions to more compactly describe the spatial variation of field perturbations enables improved characterization of fields with rapid spatial variations.

Pancreatic Intraepithelial Neoplasia Revealed by Diffusion-Tensor MRI

OBJECTIVES

Detecting premalignant lesions for pancreatic ductal adenocarcinoma, mainly pancreatic intraepithelial neoplasia (PanIN), is critical for early diagnosis and for understanding PanIN biology. Based on PanIN's histology, we hypothesized that diffusion tensor imaging (DTI) and T2* could detect PanIN.

MATERIALS AND METHODS

DTI was explored for the detection and characterization of PanIN in genetically engineered mice (KC, KPC). Following in vivo DTI, ex vivo ultrahigh-field (16.4 T) MR microscopy using DTI, T2* was performed with histological validation. Sources of MR contrasts and histological features were investigated, including histological scoring for disease burden (lesion span) and severity (adjusted score). To test if findings in mice can be translated to humans, human pancreas specimens were imaged.

RESULTS

DTI detected PanIN and pancreatic ductal adenocarcinoma in vivo (6 KPC, 4 KC, 6 controls) with high discriminative ability: fractional anisotropy (FA) and radial diffusivity with area under the curve = 0.983 (95% confidence interval: 0.932-1.000); mean diffusivity and axial diffusivity (AD) with area under the curve = 1 (95% confidence interval: 1.000-1.000). MR microscopy with histological correlation (20 KC/KPC; 5 controls) revealed that sources of MR contrasts likely arise from microarchitectural signatures: high FA, AD in fibrotic areas surrounding lesions, high diffusivities within cysts, and high T2* within lesions' stroma. The strongest histological correlations for lesion span and adjusted score were obtained with AD ( R = 0.708, P < 0.001; R = 0.789, P < 0.001, respectively). Ex vivo observations in 5 human pancreases matched our findings in mice, revealing substantial contrast between PanIN and normal pancreas.

CONCLUSIONS

DTI and T2* are useful for detecting and characterizing PanIN in genetically engineered mice and in the human pancreas, especially with AD and FA. These are encouraging findings for future clinical applications of pancreatic imaging.

Diffusion MRI with double diffusion encoding and variable mixing times disentangles water exchange from transient kurtosis

Double diffusion encoding (DDE) makes diffusion MRI sensitive to a wide range of microstructural features, and the acquired data can be analysed using different approaches. Correlation tensor imaging (CTI) uses DDE to resolve three components of the diffusional kurtosis: isotropic, anisotropic, and microscopic kurtosis. The microscopic kurtosis is estimated from the contrast between single diffusion encoding (SDE) and parallel DDE signals at the same b-value. Another approach is multi-Gaussian exchange (MGE), which employs DDE to measure exchange. Sensitivity to exchange is obtained by contrasting SDE and DDE signals at the same b-value. CTI and MGE exploit the same signal contrast to quantify microscopic kurtosis and exchange, and this study investigates the interplay between these two quantities. We perform Monte Carlo simulations in different geometries with varying levels of exchange and study the behaviour of the parameters from CTI and MGE. We conclude that microscopic kurtosis from CTI is sensitive to the exchange rate and that intercompartmental exchange and the transient kurtosis of individual compartments are distinct sources of microscopic kurtosis. In an attempt to disentangle these two sources, we propose a heuristic signal representation referred to as tMGE (MGE incorporating transient kurtosis) that accounts for both effects by exploiting the distinct signatures of exchange and transient kurtosis with varying mixing time: exchange causes a slow dependence of the signal on mixing time while transient kurtosis arguably has a much faster dependence. We find that applying tMGE to data acquired with multiple mixing times for both parallel and orthogonal DDE may enable estimation of the exchange rate as well as isotropic, anisotropic, and transient kurtosis.

Investigating microstructural changes between in vivo and perfused ex vivo marmoset brains using oscillating gradient and b-tensor encoded diffusion MRI at 9.4 T

PURPOSE

To investigate microstructural alterations induced by perfusion fixation in brain tissues using advanced diffusion MRI techniques and estimate their potential impact on the application of ex vivo models to in vivo microstructure.

METHODS

We used oscillating gradient spin echo (OGSE) and b-tensor encoding diffusion MRI to examine in vivo and ex vivo microstructural differences in the marmoset brain. OGSE was used to shorten effective diffusion times, whereas b-tensor encoding allowed for the differentiation of isotropic and anisotropic kurtosis. Additionally, we performed Monte Carlo simulations to estimate the potential microstructural changes in the tissues.

RESULTS

We report large changes (˜50%-60%) in kurtosis frequency dispersion (OGSE) and in both anisotropic and isotropic kurtosis (b-tensor encoding) after perfusion fixation. Structural MRI showed an average volume reduction of about 10%. Monte Carlo simulations indicated that these alterations could likely be attributed to extracellular fluid loss possibly combined with axon beading and increased dot compartment signal fraction. Little evidence was observed for reductions in axonal caliber.

CONCLUSION

Our findings shed light on advanced MRI parameter changes that are induced by perfusion fixation and potential microstructural sources for these changes. This work also suggests that caution should be exercised when applying ex vivo models to infer in vivo tissue microstructure, as significant differences may arise.

Tutorial: a guide to diffusion MRI and structural connectomics

Diffusion magnetic resonance imaging (dMRI) is a versatile imaging technique that has gained popularity thanks to its sensitive ability to measure displacement of water molecules within a living tissue on a micrometer scale. Although dMRI has been around since the early 1990s, its applications are constantly evolving, primarily regarding the inference of structural connectomics from nerve fiber trajectories. However, these applications require expertise in image processing and statistics, and it can be difficult for a newcomer to choose an appropriate pipeline to fit their research needs, not least because dMRI is such a flexible methodology that dozens of acquisition and analysis pipelines have been developed over the years. This introductory guide is designed for graduate students and researchers in the neuroscience community who are interested in integrating this new methodology regardless of their background in neuroimaging and computational tools. The guide provides a brief overview of the basic dMRI methodologies but focuses on its applications in neuroplasticity and connectomics. The guide starts with dMRI experimental designs and a complete step-by-step pipeline for structural connectomics. The following section covers the basics of dMRI, including parameters and clinical applications (apparent diffusion coefficient, mean diffusivity, fractional anisotropy and microscopic fractional anisotropy), as well as different approaches and models. The final section focuses on structural connectomics, covering subjects from fiber tracking (techniques, evaluation and limitations) to structural networks (constructing, analyzing and visualizing a network).

IDH Status in Brain Gliomas Can Be Predicted by the Spherical Mean MRI Technique

BACKGROUND AND PURPOSE

Diffuse gliomas, a heterogeneous group of primary brain tumors, have traditionally been stratified by histology, but recent insights into their molecular features, especially the IDH mutation status, have fundamentally changed their classification and prognosis. Current diagnostic methods, still predominantly relying on invasive biopsy, necessitate the exploration of noninvasive imaging alternatives for glioma characterization.

MATERIALS AND METHODS

In this prospective study, we investigated the utility of the spherical mean technique (SMT) in predicting the IDH status and histologic grade of adult-type diffuse gliomas. Patients with histologically confirmed adult-type diffuse glioma underwent a multiparametric MRI examination using a 3T system, which included a multishell diffusion sequence. Advanced diffusion parameters were obtained using SMT, diffusional kurtosis imaging, and ADC modeling. The diagnostic performance of studied parameters was evaluated by plotting receiver operating characteristic curves with associated area under curve, specificity, and sensitivity values.

RESULTS

A total of 80 patients with a mean age of 48 (SD, 16) years were included in the study. SMT metrics, particularly microscopic fractional anisotropy (μFA), intraneurite voxel fraction, and μFA to the third power (μFA3), demonstrated strong diagnostic performance (all AUC = 0.905, 95% CI, 0.835-0.976; P < .001) in determining IDH status and compared favorably with diffusional kurtosis imaging and ADC models. These parameters also showed a strong predictive capability for tumor grade, with intraneurite voxel fraction and μFA achieving the highest diagnostic accuracy (AUC = 0.937, 95% CI, 0.880-0.993; P < .001). Control analyses on normal-appearing brain tissue confirmed the specificity of these metrics for tumor tissue.

CONCLUSIONS

Our study highlights the potential of SMT for noninvasive characterization of adult-type diffuse gliomas, with a potential to predict IDH status and tumor grade more accurately than traditional ADC metrics. SMT offers a promising addition to the current diagnostic toolkit, enabling more precise preoperative assessments and contributing to personalized treatment planning.

In vivo rat brain mapping of multiple gray matter water populations using nonparametric D(ω)-R1-R2 distributions MRI

Massively multidimensional diffusion magnetic resonance imaging combines tensor-valued encoding, oscillating gradients, and diffusion-relaxation correlation to provide multicomponent subvoxel parameters depicting some tissue microstructural features. This method was successfully implemented ex vivo in microimaging systems and clinical conditions with tensor-valued gradient waveform of variable duration giving access to a narrow diffusion frequency (ω) range. We demonstrate here its preclinical in vivo implementation with a protocol of 389 contrast images probing a wide diffusion frequency range of 18 to 92 Hz at b-values up to 2.1 ms/μm2 enabled by the use of modulated gradient waveforms and combined with multislice high-resolution and low-distortion echo planar imaging acquisition with segmented and full reversed phase-encode acquisition. This framework allows the identification of diffusion ω-dependence in the rat cerebellum and olfactory bulb gray matter (GM), and the parameter distributions are shown to resolve two water pools in the cerebellum GM with different diffusion coefficients, shapes, ω-dependence, relaxation rates, and spatial repartition whose attribution to specific microstructure could modify the current understanding of the origin of restriction in GM.

Tensor-valued diffusion MRI detects brain microstructural abnormalities in HIV infected individuals with cognitive impairment

Despite advancements, the prevalence of HIV-associated neurocognitive impairment remains at approximately 40%, attributed to factors like pre-cART (combination antiretroviral therapy) irreversible brain injury. People with HIV (PWH) treated with cART do not show significant neurocognitive changes over relatively short follow-up periods. However, quantitative neuroimaging may be able to detect ongoing subtle microstructural changes. In this study, we hypothesized that tensor-valued diffusion encoding metrics would provide greater sensitivity than conventional diffusion tensor imaging (DTI) metrics in detecting HIV-associated brain microstructural injury. We further hypothesized that tensor-valued metrics would exhibit stronger associations with blood markers of neuronal and glial injury, such as neurofilament light chain (NFL) and glial fibrillary acidic protein (GFAP), as well as with cognitive performance. Using MRI at 3T, 24 PWH and 31 healthy controls underwent cross-sectional examination. The results revealed significant variations in tensor-valued diffusion encoding metrics across white matter regions, with associations observed between these metrics, cognitive performance, NFL and GFAP. Moreover, a significant interaction between HIV status and imaging metrics in gray and white matter was observed, particularly impacting total cognitive scores. Of interest, DTI metrics were less likely to be associated with HIV status than tensor-valued diffusion metrics. These findings suggest that tensor-valued diffusion encoding metrics offer heightened sensitivity in detecting subtle changes associated with axonal injury in HIV infection. Longitudinal studies are needed to further evaluate responsiveness of tensor-valued diffusion b-tensor encoding metrics in the contest HIV-associate mild chronic neuroinflammation.

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.

Nonparametric distributions of tensor-valued Lorentzian diffusion spectra for model-free data inversion in multidimensional diffusion MRI

Magnetic resonance imaging (MRI) is the method of choice for noninvasive studies of micrometer-scale structures in biological tissues via their effects on the time- and frequency-dependent (restricted) and anisotropic self-diffusion of water. While new designs of time-dependent magnetic field gradient waveforms have enabled disambiguation between different aspects of translational motion that are convolved in traditional MRI methods relying on single pairs of field gradient pulses, data analysis for complex heterogeneous materials remains a challenge. Here, we propose and demonstrate nonparametric distributions of tensor-valued Lorentzian diffusion spectra, or "D(ω) distributions," as a general representation with sufficient flexibility to describe the MRI signal response from a wide range of model systems and biological tissues investigated with modulated gradient waveforms separating and correlating the effects of restricted and anisotropic diffusion.

Distinct Longitudinal Brain White Matter Microstructure Changes and Associated Polygenic Risk of Common Psychiatric Disorders and Alzheimer's Disease in the UK Biobank

Background

During the course of adulthood and aging, white matter (WM) structure and organization are characterized by slow degradation processes such as demyelination and shrinkage. An acceleration of such aging processes has been linked to the development of a range of diseases. Thus, an accurate description of healthy brain maturation, particularly in terms of WM features, is fundamental to the understanding of aging.

Methods

We used longitudinal diffusion magnetic resonance imaging to provide an overview of WM changes at different spatial and temporal scales in the UK Biobank (UKB) (n = 2678; agescan 1 = 62.38 ± 7.23 years; agescan 2 = 64.81 ± 7.1 years). To examine the genetic overlap between WM structure and common clinical conditions, we tested the associations between WM structure and polygenic risk scores for the most common neurodegenerative disorder, Alzheimer's disease, and common psychiatric disorders (unipolar and bipolar depression, anxiety, obsessive-compulsive disorder, autism, schizophrenia, attention-deficit/hyperactivity disorder) in longitudinal (n = 2329) and cross-sectional (n = 31,056) UKB validation data.

Results

Our findings indicate spatially distributed WM changes across the brain, as well as distributed associations of polygenic risk scores with WM. Importantly, brain longitudinal changes reflected genetic risk for disorder development better than the utilized cross-sectional measures, with regional differences giving more specific insights into gene-brain change associations than global averages.

Conclusions

We extend recent findings by providing a detailed overview of WM microstructure degeneration on different spatial levels, helping to understand fundamental brain aging processes. Further longitudinal research is warranted to examine aging-related gene-brain associations.

Tensor-valued diffusion MRI detects brain microstructure changes in HIV infected individuals with cognitive impairment

Despite advancements, the prevalence of HIV-associated neurocognitive impairment remains at approximately 40%, attributed to factors like pre-cART (combination antiretroviral therapy) irreversible brain injury. People with HIV (PWH) treated with cART do not show significant neurocognitive changes over relatively short follow-up periods. However, quantitative neuroimaging may be able to detect ongoing subtle microstructural changes. This study aimed to investigate the sensitivity of tensor-valued diffusion encoding in detecting such changes in brain microstructural integrity in cART-treated PWH. Additionally, it explored relationships between these metrics, neurocognitive scores, and plasma levels of neurofilament light (NFL) chain and glial fibrillary acidic protein (GFAP). Using MRI at 3T, 24 PWH and 31 healthy controls underwent cross-sectional examination. The results revealed significant variations in b-tensor encoding metrics across white matter regions, with associations observed between these metrics, cognitive performance, and blood markers of neuronal and glial injury (NFL and GFAP). Moreover, a significant interaction between HIV status and imaging metrics was observed, particularly impacting total cognitive scores in both gray and white matter. These findings suggest that b-tensor encoding metrics offer heightened sensitivity in detecting subtle changes associated with axonal injury in HIV infection, underscoring their potential clinical relevance in understanding neurocognitive impairment in PWH.

FlexDTI: flexible diffusion gradient encoding scheme-based highly efficient diffusion tensor imaging using deep learning

Objective. Most deep neural network-based diffusion tensor imaging methods require the diffusion gradients' number and directions in the data to be reconstructed to match those in the training data. This work aims to develop and evaluate a novel dynamic-convolution-based method called FlexDTI for highly efficient diffusion tensor reconstruction with flexible diffusion encoding gradient scheme.Approach. FlexDTI was developed to achieve high-quality DTI parametric mapping with flexible number and directions of diffusion encoding gradients. The method used dynamic convolution kernels to embed diffusion gradient direction information into feature maps of the corresponding diffusion signal. Furthermore, it realized the generalization of a flexible number of diffusion gradient directions by setting the maximum number of input channels of the network. The network was trained and tested using datasets from the Human Connectome Project and local hospitals. Results from FlexDTI and other advanced tensor parameter estimation methods were compared.Main results. Compared to other methods, FlexDTI successfully achieves high-quality diffusion tensor-derived parameters even if the number and directions of diffusion encoding gradients change. It reduces normalized root mean squared error by about 50% on fractional anisotropy and 15% on mean diffusivity, compared with the state-of-the-art deep learning method with flexible diffusion encoding gradient scheme.Significance. FlexDTI can well learn diffusion gradient direction information to achieve generalized DTI reconstruction with flexible diffusion gradient scheme. Both flexibility and reconstruction quality can be taken into account in this network.

Deep learning with diffusion MRI as in vivo microscope reveals sex-related differences in human white matter microstructure

Biological sex is a crucial variable in neuroscience studies where sex differences have been documented across cognitive functions and neuropsychiatric disorders. While gross statistical differences have been previously documented in macroscopic brain structure such as cortical thickness or region size, less is understood about sex-related cellular-level microstructural differences which could provide insight into brain health and disease. Studying these microstructural differences between men and women paves the way for understanding brain disorders and diseases that manifest differently in different sexes. Diffusion MRI is an important in vivo, non-invasive methodology that provides a window into brain tissue microstructure. Our study develops multiple end-to-end classification models that accurately estimates the sex of a subject using volumetric diffusion MRI data and uses these models to identify white matter regions that differ the most between men and women. 471 male and 560 female healthy subjects (age range, 22-37 years) from the Human Connectome Project are included. Fractional anisotropy, mean diffusivity and mean kurtosis are used to capture brain tissue microstructure characteristics. Diffusion parametric maps are registered to a standard template to reduce bias that can arise from macroscopic anatomical differences like brain size and contour. This study employ three major model architectures: 2D convolutional neural networks, 3D convolutional neural networks and Vision Transformer (with self-supervised pretraining). Our results show that all 3 models achieve high sex classification performance (test AUC 0.92-0.98) across all diffusion metrics indicating definitive differences in white matter tissue microstructure between males and females. We further use complementary model architectures to inform about the pattern of detected microstructural differences and the influence of short-range versus long-range interactions. Occlusion analysis together with Wilcoxon signed-rank test is used to determine which white matter regions contribute most to sex classification. The results indicate that sex-related differences manifest in both local features as well as global features / longer-distance interactions of tissue microstructure. Our highly consistent findings across models provides new insight supporting differences between male and female brain cellular-level tissue organization particularly in the central white matter.

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.

Advanced magnetic resonance neuroimaging techniques: feasibility and applications in long or post-COVID-19 syndrome - a review

Long-term or post-COVID-19 syndrome (PCS) is a condition that affects people infected with SARS‑CoV‑2, the virus that causes COVID-19. PCS is characterized by a wide range of persistent or new symptoms that last months after the initial infection, such as fatigue, shortness of breath, cognitive dysfunction, and pain. Advanced magnetic resonance (MR) neuroimaging techniques can provide valuable information on the structural and functional changes in the brain associated with PCS as well as potential biomarkers for diagnosis and prognosis. In this review, we discuss the feasibility and applications of various advanced MR neuroimaging techniques in PCS, including perfusion-weighted imaging (PWI), diffusion-weighted imaging (DWI), susceptibility-weighted imaging (SWI), functional MR imaging (fMRI), diffusion tensor imaging (DTI), and tractography. We summarize the current evidence on neuroimaging findings in PCS, the challenges and limitations of these techniques, and the future directions for research and clinical practice. Although still uncertain, advanced MRI techniques show promise for gaining insight into the pathophysiology and guiding the management of COVID-19 syndrome, pending larger validation studies.

Advanced diffusion imaging reveals microstructural characteristics of primary CNS lymphoma, allowing differentiation from glioblastoma

Background

Primary CNS lymphoma (PCNSL) and glioblastoma (GBM) both represent frequent intracranial malignancies with differing clinical management. However, distinguishing PCNSL from GBM with conventional MRI can be challenging when atypical imaging features are present. We employed advanced dMRI for noninvasive characterization of the microstructure of PCNSL and differentiation from GBM as the most frequent primary brain malignancy.

Methods

Multiple dMRI metrics including Diffusion Tensor Imaging, Neurite Orientation Dispersion and Density Imaging, and Diffusion Microstructure Imaging were extracted from the contrast-enhancing tumor component in 10 PCNSL and 10 age-matched GBM on 3T MRI. Imaging findings were correlated with cell density and axonal markers obtained from histopathology.

Results

We found significantly increased intra-axonal volume fractions (V-intra and intracellular volume fraction) and microFA in PCNSL compared to GBM (all P < .001). In contrast, mean diffusivity (MD), axial diffusivity (aD), and microADC (all P < .001), and also free water fractions (V-CSF and V-ISO) were significantly lower in PCNSL (all P < .01). Receiver-operating characteristic analysis revealed high predictive values regarding the presence of a PCNSL for MD, aD, microADC, V-intra, ICVF, microFA, V-CSF, and V-ISO (area under the curve [AUC] in all >0.840, highest for MD and ICVF with an AUC of 0.960). Comparative histopathology between PCNSL and GBM revealed a significantly increased cell density in PCNSL and the presence of axonal remnants in a higher proportion of samples.

Conclusions

Advanced diffusion imaging enables the characterization of the microstructure of PCNSL and reliably distinguishes PCNSL from GBM. Both imaging and histopathology revealed a relatively increased cell density and a preserved axonal microstructure in PCNSL.

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.

Mesoscopic Assessment of Microstructure in Glioblastomas and Metastases by Merging Advanced Diffusion Imaging with Immunohistopathology

BACKGROUND AND PURPOSE

Glioblastomas and metastases are the most common malignant intra-axial brain tumors in adults and can be difficult to distinguish on conventional MR imaging due to similar imaging features. We used advanced diffusion techniques and structural histopathology to distinguish these tumor entities on the basis of microstructural axonal and fibrillar signatures in the contrast-enhancing tumor component.

MATERIALS AND METHODS

Contrast-enhancing tumor components were analyzed in 22 glioblastomas and 21 brain metastases on 3T MR imaging using DTI-fractional anisotropy, neurite orientation dispersion and density imaging-orientation dispersion, and diffusion microstructural imaging-micro-fractional anisotropy. Available histopathologic specimens (10 glioblastomas and 9 metastases) were assessed for the presence of axonal structures and scored using 4-level scales for Bielschowsky staining (0: no axonal structures, 1: minimal axonal fragments preserved, 2: decreased axonal density, 3: no axonal loss) and glial fibrillary acid protein expression (0: no glial fibrillary acid protein positivity, 1: limited expression, 2: equivalent to surrounding parenchyma, 3: increased expression).

RESULTS

When we compared glioblastomas and metastases, fractional anisotropy was significantly increased and orientation dispersion was decreased in glioblastomas (each P < .001), with a significant shift toward increased glial fibrillary acid protein and Bielschowsky scores. Positive associations of fractional anisotropy and negative associations of orientation dispersion with glial fibrillary acid protein and Bielschowsky scores were revealed, whereas no association between micro-fractional anisotropy with glial fibrillary acid protein and Bielschowsky scores was detected. Receiver operating characteristic curves revealed high predictive values of both fractional anisotropy (area under the curve = 0.8463) and orientation dispersion (area under the curve = 0.8398) regarding the presence of a glioblastoma.

CONCLUSIONS

Diffusion imaging fractional anisotropy and orientation dispersion metrics correlated with histopathologic markers of directionality and may serve as imaging biomarkers in contrast-enhancing tumor components.

Mapping the individual human cortex using multidimensional MRI and unsupervised learning

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

An open resource combining multi-contrast MRI and microscopy in the macaque brain

Understanding brain structure and function often requires combining data across different modalities and scales to link microscale cellular structures to macroscale features of whole brain organisation. Here we introduce the BigMac dataset, a resource combining in vivo MRI, extensive postmortem MRI and multi-contrast microscopy for multimodal characterisation of a single whole macaque brain. The data spans modalities (MRI and microscopy), tissue states (in vivo and postmortem), and four orders of spatial magnitude, from microscopy images with micrometre or sub-micrometre resolution, to MRI signals on the order of millimetres. Crucially, the MRI and microscopy images are carefully co-registered together to facilitate quantitative multimodal analyses. Here we detail the acquisition, curation, and first release of the data, that together make BigMac a unique, openly-disseminated resource available to researchers worldwide. Further, we demonstrate example analyses and opportunities afforded by the data, including improvement of connectivity estimates from ultra-high angular resolution diffusion MRI, neuroanatomical insight provided by polarised light imaging and myelin-stained histology, and the joint analysis of MRI and microscopy data for reconstruction of the microscopy-inspired connectome. All data and code are made openly available.

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.

CrOssing fiber Modeling in the Peritumoral Area using dMRI (COMPARI)

Visualization of fiber tracts around the tumor is critical for neurosurgical planning and preservation of crucial structural connectivity during tumor resection. Biophysical modeling approaches estimate fiber tract orientations from differential water diffusivity information of diffusion MRI. However, the presence of edema and tumor infiltration presents a challenge to visualize crossing fiber tracts in the peritumoral region. Previous approaches proposed free water modeling to compensate for the effect of water diffusivity in edema, but those methods were limited in estimating complex crossing fiber tracts. We propose a new cascaded multi-compartment model to estimate tissue microstructure in the presence of edema and pathological contaminants in the area surrounding brain tumors. In our model (COMPARI), the isotropic components of diffusion signal, including free water and hindered water, were eliminated, and the fiber orientation distribution (FOD) of the remaining signal was estimated. In simulated data, COMPARI accurately recovered fiber orientations in the presence of extracellular water. In a dataset of 23 patients with highly edematous brain tumors, the amplitudes of FOD and anisotropic index distribution within the peritumoral region were higher with COMPARI than with a recently proposed multi-compartment constrained deconvolution model. In a selected patient with metastatic brain tumor, we demonstrated COMPARI's ability to effectively model and eliminate water from the peritumoral region. The white matter bundles reconstructed with our model were qualitatively improved compared to those of other models, and allowed the identification of crossing fibers. In conclusion, the removal of isotropic components as proposed with COMPARI improved the bio-physical modeling of dMRI in edema, thus providing information on crossing fibers, thereby enabling improved tractography in a highly edematous brain tumor. This model may improve surgical planning tools to help achieve maximal safe resection of brain tumors.

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.

A longitudinal microstructural MRI dataset in healthy C57Bl/6 mice at 9.4 Tesla

Multimodal microstructural MRI has shown increased sensitivity and specificity to changes in various brain disease and injury models in the preclinical setting. Here, we present an in vivo longitudinal dataset, including a subset of ex vivo data, acquired as control data and to investigate microstructural changes in the healthy mouse brain. The dataset consists of structural T2-weighted imaging, magnetization transfer ratio and saturation imaging, and advanced quantitative diffusion MRI (dMRI) methods. The dMRI methods include oscillating gradient spin echo (OGSE) dMRI and microscopic anisotropy (μA) dMRI, which provide additional insight by increasing sensitivity to smaller spatial scales and disentangling fiber orientation dispersion from true microstructural changes, respectively. The technical skills required to analyze microstructural MRI data are complex and include MRI sequence development, acquisition, and computational neuroimaging expertise. Here, we share unprocessed and preprocessed data, and scalar maps of quantitative MRI metrics. We envision utility of this dataset in the microstructural MRI field to develop and test biophysical models, methods that model temporal brain dynamics, and registration and preprocessing pipelines.

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.

Estimation of free water-corrected microscopic fractional anisotropy

Water diffusion anisotropy MRI is sensitive to microstructural changes in the brain that are hallmarks of various neurological conditions. However, conventional metrics like fractional anisotropy are confounded by neuron fiber orientation dispersion, and the relatively low resolution of diffusion-weighted MRI gives rise to significant free water partial volume effects in many brain regions that are adjacent to cerebrospinal fluid. Microscopic fractional anisotropy is a recent metric that can report water diffusion anisotropy independent of neuron fiber orientation dispersion but is still susceptible to free water contamination. In this paper, we present a free water elimination (FWE) technique to estimate microscopic fractional anisotropy and other related diffusion indices by implementing a signal representation in which the MRI signal within a voxel is assumed to come from two distinct sources: a tissue compartment and a free water compartment. A two-part algorithm is proposed to rapidly fit a set of diffusion-weighted MRI volumes containing both linear- and spherical-tensor encoding acquisitions to the representation. Simulations and in vivo acquisitions with four healthy volunteers indicated that the FWE method may be a feasible technique for measuring microscopic fractional anisotropy and other indices with greater specificity to neural tissue characteristics than conventional methods.

Diffusion MRI with pulsed and free gradient waveforms: Effects of restricted diffusion and exchange

Monitoring time dependence with diffusion MRI yields observables sensitive to compartment sizes (restricted diffusion) and membrane permeability (water exchange). However, restricted diffusion and exchange have opposite effects on the diffusion-weighted signal, which can lead to errors in parameter estimates. In this work, we propose a signal representation that incorporates the effects of both restricted diffusion and exchange up to second order in b-value and is compatible with gradient waveforms of arbitrary shape. The representation features mappings from a gradient waveform to two scalars that separately control the sensitivity to restriction and exchange. We demonstrate that these scalars span a two-dimensional space that can be used to choose waveforms that selectively probe restricted diffusion or exchange, eliminating the correlation between the two phenomena. We found that waveforms with specific but unconventional shapes provide an advantage over conventional pulsed and oscillating gradient acquisitions. We also show that parametrization of waveforms into a two-dimensional space can be used to understand protocols from other approaches that probe restricted diffusion and exchange. For example, we found that the variation of mixing time in filter-exchange imaging corresponds to variation of our exchange-weighting scalar at a fixed value of the restriction-weighting scalar. The proposed signal representation was evaluated using Monte Carlo simulations in identical parallel cylinders with hexagonal and random packing as well as parallel cylinders with gamma-distributed radii. Results showed that the approach is sensitive to sizes in the interval 4-12 μm and exchange rates in the simulated range of 0 to 20 s - 1 , but also that there is a sensitivity to the extracellular geometry. The presented theory constitutes a simple and intuitive description of how restricted diffusion and exchange influence the signal as well as a guide to protocol design capable of separating the two effects.

Estimating axial diffusivity in the NODDI model

To estimate microstructure-related parameters from diffusion MRI data, biophysical models make strong, simplifying assumptions about the underlying tissue. The extent to which many of these assumptions are valid remains an open research question. This study was inspired by the disparity between the estimated intra-axonal axial diffusivity from literature and that typically assumed by the Neurite Orientation Dispersion and Density Imaging (NODDI) model (d_∥=1.7μm²/ms). We first demonstrate how changing the assumed axial diffusivity results in considerably different NODDI parameter estimates. Second, we illustrate the ability to estimate axial diffusivity as a free parameter of the model using high b-value data and an adapted NODDI framework. Using both simulated and in vivo data we investigate the impact of fitting to either real-valued or magnitude data, with Gaussian and Rician noise characteristics respectively, and what happens if we get the noise assumptions wrong in this high b-value and thus low SNR regime. Our results from real-valued human data estimate intra-axonal axial diffusivities of ∼2-2.5μm²/ms, in line with current literature. Crucially, our results demonstrate the importance of accounting for both a rectified noise floor and/or a signal offset to avoid biased parameter estimates when dealing with low SNR data.

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.

Multi-compartment diffusion magnetic resonance imaging models link tract-related characteristics with working memory performance in healthy older adults

Multi-compartment diffusion MRI metrics [such as metrics from free water elimination diffusion tensor imaging (FWE-DTI) and neurite orientation dispersion and density imaging (NODDI)] may reflect more specific underlying white-matter tract characteristics than traditional, single-compartment metrics [i.e., metrics from Diffusion Tensor Imaging (DTI)]. However, it remains unclear if multi-compartment metrics are more closely associated with age and/or cognitive performance than single-compartment metrics. Here we compared the associations of single-compartment [Fractional Anisotropy (FA)] and multi-compartment diffusion MRI metrics [FWE-DTI metrics: Free Water Eliminated Fractional Anisotropy (FWE-FA) and Free Water (FW); NODDI metrics: Intracellular Volume Fraction (ICVF), Orientation Dispersion Index (ODI), and CSF-Fraction] with both age and working memory performance. A functional magnetic resonance imaging (fMRI) guided, white matter tractography approach was employed to compute diffusion metrics within a network of tracts connecting functional regions involved in working memory. Ninety-nine healthy older adults (aged 60-85) performed an in-scanner working memory task while fMRI was performed and also underwent multi-shell diffusion acquisition. The network of white matter tracts connecting functionally-activated regions was identified using probabilistic tractography. Diffusion metrics were extracted from skeletonized white matter tracts connecting fMRI activation peaks. Diffusion metrics derived from both single and multi-compartment models were associated with age (p s ≤ 0.011 for FA, FWE-FA, ICVF and ODI). However, only multi-compartment metrics, specifically FWE-FA (p = 0.045) and ICVF (p = 0.020), were associated with working memory performance. Our results suggest that while most current diffusion metrics are sensitive to age, several multi-compartment metrics (i.e., FWE-FA and ICVF) appear more sensitive to cognitive performance in healthy older adults.

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

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

Performance of orientation distribution function-fingerprinting with a biophysical multicompartment diffusion model

PURPOSE

Orientation Distribution Function (ODF) peak finding methods typically fail to reconstruct fibers crossing at shallow angles below 40°, leading to errors in tractography. ODF-Fingerprinting (ODF-FP) with the biophysical multicompartment diffusion model allows for breaking this barrier.

METHODS

A randomized mechanism to generate a multidimensional ODF-dictionary that covers biologically plausible ranges of intra- and extra-axonal diffusivities and fraction volumes is introduced. This enables ODF-FP to address the high variability of brain tissue. The performance of the proposed approach is evaluated on both numerical simulations and a reconstruction of major fascicles from high- and low-resolution in vivo diffusion images.

RESULTS

ODF-FP with the suggested modifications correctly identifies fibers crossing at angles as shallow as 10 degrees in the simulated data. In vivo, our approach reaches 56% of true positives in determining fiber directions, resulting in visibly more accurate reconstruction of pyramidal tracts, arcuate fasciculus, and optic radiations than the state-of-the-art techniques. Moreover, the estimated diffusivity values and fraction volumes in corpus callosum conform with the values reported in the literature.

CONCLUSION

The modified ODF-FP outperforms commonly used fiber reconstruction methods at shallow angles, which improves deterministic tractography outcomes of major fascicles. In addition, the proposed approach allows for linearization of the microstructure parameters fitting problem.

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 diffusion MRI 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 diffusion MRI and where the field is headed in the coming years.

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.

Sensitivity of Diffusion MRI to White Matter Pathology: Influence of Diffusion Protocol, Magnetic Field Strength, and Processing Pipeline in Systemic Lupus Erythematosus

There are many ways to acquire and process diffusion MRI (dMRI) data for group studies, but it is unknown which maximizes the sensitivity to white matter (WM) pathology. Inspired by this question, we analyzed data acquired for diffusion tensor imaging (DTI) and diffusion kurtosis imaging (DKI) at 3T (3T-DTI and 3T-DKI) and DTI at 7T in patients with systemic lupus erythematosus (SLE) and healthy controls (HC). Parameter estimates in 72 WM tracts were obtained using TractSeg. The impact on the sensitivity to WM pathology was evaluated for the diffusion protocol, the magnetic field strength, and the processing pipeline. Sensitivity was quantified in terms of Cohen's d for group comparison. Results showed that the choice of diffusion protocol had the largest impact on the effect size. The effect size in fractional anisotropy (FA) across all WM tracts was 0.26 higher when derived by DTI than by DKI and 0.20 higher in 3T compared with 7T. The difference due to the diffusion protocol was larger than the difference due to magnetic field strength for the majority of diffusion parameters. In contrast, the difference between including or excluding different processing steps was near negligible, except for the correction of distortions from eddy currents and motion which had a clearly positive impact. For example, effect sizes increased on average by 0.07 by including motion and eddy correction for FA derived from 3T-DTI. Effect sizes were slightly reduced by the incorporation of denoising and Gibbs-ringing removal (on average by 0.011 and 0.005, respectively). Smoothing prior to diffusion model fitting generally reduced effect sizes. In summary, 3T-DTI in combination with eddy current and motion correction yielded the highest sensitivity to WM pathology in patients with SLE. However, our results also indicated that the 3T-DKI and 7T-DTI protocols used here may be adjusted to increase effect sizes.

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.

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

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

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.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

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.

Microscopic Fractional Anisotropy Detects Cognitive Training-Induced Microstructural Brain Changes

Cognitive training-induced neuroplastic brain changes have been reported. This prospective study evaluated whether microscopic fractional anisotropy (μFA) derived from double diffusion encoding (DDE) MRI could detect brain changes following a 4 week cognitive training. Twenty-nine healthy volunteers were recruited and randomly assigned into the training (n = 21) and control (n = 8) groups. Both groups underwent brain MRI including DDE MRI and 3D-T1-weighted imaging twice at an interval of 4–6 weeks, during which the former underwent the training. The training consisted of hour-long dual N-back and attention network tasks conducted five days per week. Training and time-related changes of DDE MRI indices (μFA, fractional anisotropy (FA), and mean diffusivity (MD)) and the gray and white matter volume were evaluated using mixed-design analysis of variance. In addition, any significant imaging indices were tested for correlation with cognitive training-induced task performance changes, using partial correlation analyses. μFA in the left middle frontal gyrus decreased upon the training (53 voxels, uncorrected p < 0.001), which correlated moderately with response time changes in the orienting component of attention (r = −0.521, uncorrected p = 0.032). No significant training and time-related changes were observed for other imaging indices. Thus, μFA can become a sensitive index to detect cognitive training-induced neuroplastic changes.

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.

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

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

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.