Multidimensional diffusion MRI

Tuned exchange imaging: Can the filter exchange imaging pulse sequence be adapted for applications with thin slices and restricted diffusion?

Filter exchange imaging (FEXI) is a double diffusion-encoding (DDE) sequence that is specifically sensitive to exchange between sites with different apparent diffusivities. FEXI uses a diffusion-encoding filtering block followed by a detection block at varying mixing times to map the exchange rate. Long mixing times enhance the sensitivity to exchange, but they pose challenges for imaging applications that require a stimulated echo sequence with crusher gradients. Thin imaging slices require strong crushers, which can introduce significant diffusion weighting and bias exchange rate estimates. Here, we treat the crushers as an additional encoding block and consider FEXI as a triple diffusion-encoding sequence. This allows the bias to be corrected in the case of multi-Gaussian diffusion, but not easily in the presence of restricted diffusion. Our approach addresses challenges in the presence of restricted diffusion and relies on the ability to independently gauge sensitivities to exchange and restricted diffusion for arbitrary gradient waveforms. It follows two principles: (i) the effects of crushers are included in the forward model using signal cumulant expansion; and (ii) timing parameters of diffusion gradients in filter and detection blocks are adjusted to maintain the same level of restriction encoding regardless of the mixing time. This results in the tuned exchange imaging (TEXI) protocol. The accuracy of exchange mapping with TEXI was assessed through Monte Carlo simulations in spheres of identical sizes and gamma-distributed sizes, and in parallel hexagonally packed cylinders. The simulations demonstrate that TEXI provides consistent exchange rates regardless of slice thickness and restriction size, even with strong crushers. However, the accuracy depends on b-values, mixing times, and restriction geometry. The constraints and limitations of TEXI are discussed, including suggestions for protocol adaptations. Further studies are needed to optimize the precision of TEXI and assess the approach experimentally in realistic, heterogeneous substrates.

Direct segmentation of cortical cytoarchitectonic domains using ultra-high-resolution whole-brain diffusion MRI

We assess the potential of detecting cortical laminar patterns and areal borders by directly clustering voxel values of microstructural parameters derived from high-resolution mean apparent propagator (MAP) magnetic resonance imaging (MRI), as an alternative to conventional template-warping-based cortical parcellation methods. We acquired MAP-MRI data with 200μm resolution in a fixed macaque monkey brain. To improve the sensitivity to cortical layers, we processed the data with a local anisotropic Gaussian filter determined voxel-wise by the plane tangent to the cortical surface. We directly clustered all cortical voxels using only the MAP-derived microstructural imaging biomarkers, with no information regarding their relative spatial location or dominant diffusion orientations. MAP-based 3D cytoarchitectonic segmentation revealed laminar patterns similar to those observed in the corresponding histological images. Moreover, transition regions between these laminar patterns agreed more accurately with histology than the borders between cortical areas estimated using conventional atlas/template-warping cortical parcellation. By cross-tabulating all cortical labels in the atlas- and MAP-based segmentations, we automatically matched the corresponding MAP-derived clusters (i.e., cytoarchitectonic domains) across the left and right hemispheres. Our results demonstrate that high-resolution MAP-MRI biomarkers can effectively delineate three-dimensional cortical cytoarchitectonic domains in single individuals. Their intrinsic tissue microstructural contrasts enable the construction of whole-brain mesoscopic cortical atlases.

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.

Menstrual cycle-driven hormone concentrations co-fluctuate with white and gray matter architecture changes across the whole brain

Cyclic fluctuations in hypothalamic-pituitary-gonadal axis (HPG-axis) hormones exert powerful behavioral, structural, and functional effects through actions on the mammalian central nervous system. Yet, very little is known about how these fluctuations alter the structural nodes and information highways of the human brain. In a study of 30 naturally cycling women, we employed multidimensional diffusion and T1-weighted imaging during three estimated menstrual cycle phases (menses, ovulation, and mid-luteal) to investigate whether HPG-axis hormone concentrations co-fluctuate with alterations in white matter (WM) microstructure, cortical thickness (CT), and brain volume. Across the whole brain, 17β-estradiol and luteinizing hormone (LH) concentrations were directly proportional to diffusion anisotropy (μFA; 17β-estradiol: β1 = 0.145, highest density interval (HDI) = [0.211, 0.4]; LH: β1 = 0.111, HDI = [0.157, 0.364]), while follicle-stimulating hormone (FSH) was directly proportional to CT (β1 = 0 .162, HDI = [0.115, 0.678]). Within several individual regions, FSH and progesterone demonstrated opposing relationships with mean diffusivity (Diso) and CT. These regions mainly reside within the temporal and occipital lobes, with functional implications for the limbic and visual systems. Finally, progesterone was associated with increased tissue (β1 = 0.66, HDI = [0.607, 15.845]) and decreased cerebrospinal fluid (CSF; β1 = -0.749, HDI = [-11.604, -0.903]) volumes, with total brain volume remaining unchanged. These results are the first to report simultaneous brain-wide changes in human WM microstructure and CT coinciding with menstrual cycle-driven hormone rhythms. Effects were observed in both classically known HPG-axis receptor-dense regions (medial temporal lobe, prefrontal cortex) and in other regions located across frontal, occipital, temporal, and parietal lobes. Our results suggest that HPG-axis hormone fluctuations may have significant structural impacts across the entire brain.

Signal drift in diffusion MRI of the brain: effects on intravoxel incoherent motion parameter estimates

OBJECTIVES

Signal drift has been put forward as one of the fundamental confounding factors in diffusion MRI (dMRI) of the brain. This study characterizes signal drift in dMRI of the brain, evaluates correction methods, and exemplifies its impact on parameter estimation for three intravoxel incoherent motion (IVIM) protocols.

MATERIALS AND METHODS

dMRI of the brain was acquired in ten healthy subjects using protocols designed to enable retrospective characterization and correction of signal drift. All scans were acquired twice for repeatability analysis. Three temporal polynomial correction methods were evaluated: (1) global, (2) voxelwise, and (3) spatiotemporal. Effects of acquisition order were simulated using estimated drift fields.

RESULTS

Signal drift was around 2% per 5 min in the brain as a whole, but reached above 5% per 5 min in the frontal regions. Only correction methods taking spatially varying signal drift into account could achieve effective corrections. Altered acquisition order introduced both systematic changes and differences in repeatability in the presence of signal drift.

DISCUSSION

Signal drift in dMRI of the brain was found to be spatially varying, calling for correction methods taking this into account. Without proper corrections, choice of protocol can affect dMRI parameter estimates and their repeatability.

Neuronal microstructural changes in the human brain are associated with neurocognitive aging

Gray matter (GM) alterations play a role in aging-related disorders like Alzheimer's disease and related dementias, yet MRI studies mainly focus on macroscopic changes. Although reliable indicators of atrophy, morphological metrics like cortical thickness lack the sensitivity to detect early changes preceding visible atrophy. Our study aimed at exploring the potential of diffusion MRI in unveiling sensitive markers of cortical and subcortical age-related microstructural changes and assessing their associations with cognitive and behavioral deficits. We leveraged the Human Connectome Project-Aging cohort that included 707 participants (394 female; median age = 58, range = 36-90 years) and applied the powerful mean apparent diffusion propagator model to measure microstructural parameters, along with comprehensive behavioral and cognitive test scores. Both macro- and microstructural GM characteristics were strongly associated with age, with widespread significant microstructural correlations reflective of cellular morphological changes, reduced cellular density, increased extracellular volume, and increased membrane permeability. Importantly, when correlating MRI and cognitive test scores, our findings revealed no link between macrostructural volumetric changes and neurobehavioral performance. However, we found that cellular and extracellular alterations in cortical and subcortical GM regions were associated with neurobehavioral performance. Based on these findings, it is hypothesized that increased microstructural heterogeneity and decreased neurite orientation dispersion precede macrostructural changes, and that they play an important role in subsequent cognitive decline. These alterations are suggested to be early markers of neurocognitive performance that may distinctly aid in identifying the mechanisms underlying phenotypic aging and subsequent age-related functional decline.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Multiparametric MRI for characterization of the tumour microenvironment

Our understanding of tumour biology has evolved over the past decades and cancer is now viewed as a complex ecosystem with interactions between various cellular and non-cellular components within the tumour microenvironment (TME) at multiple scales. However, morphological imaging remains the mainstay of tumour staging and assessment of response to therapy, and the characterization of the TME with non-invasive imaging has not yet entered routine clinical practice. By combining multiple MRI sequences, each providing different but complementary information about the TME, multiparametric MRI (mpMRI) enables non-invasive assessment of molecular and cellular features within the TME, including their spatial and temporal heterogeneity. With an increasing number of advanced MRI techniques bridging the gap between preclinical and clinical applications, mpMRI could ultimately guide the selection of treatment approaches, precisely tailored to each individual patient, tumour and therapeutic modality. In this Review, we describe the evolving role of mpMRI in the non-invasive characterization of the TME, outline its applications for cancer detection, staging and assessment of response to therapy, and discuss considerations and challenges for its use in future medical applications, including personalized integrated diagnostics.

A genetic optimisation and iterative reconstruction framework for sparse multi-dimensional diffusion-relaxation correlation MRI

Multi-dimensional diffusion-relaxation correlation (DRC) magnetic resonance imaging (MRI) techniques have recently been developed to investigate tissue microstructures. Sub-voxel tissue heterogeneity is resolved from the local correlation distributions of relaxation time and molecular diffusivity. However, the implementation of these techniques considerably increases the total acquisition time, and simply reducing the scan time may be at the expense of detailed structural resolution. To overcome these limitations, an optimised framework was proposed for acquiring microstructural maps of the human brain on a clinically feasible timescale. First, the acquisition parameters of the multi-dimensional DRC MRI method were sparsely optimised using a genetic algorithm with a fitness function according to the spectral resolution of the correlation map, hardware requirements, and total scan time. Next, the acquired DRC MRI data were processed using a proposed numerical algorithm based on the dynamic inverse Laplace transform (ILT). Prior knowledge from one-dimensional data was then utilised in the iterative procedure to improve the spectral resolution. Finally, the proposed framework was validated using Monte Carlo simulations and experimental data acquired from healthy participants on an MRI scanner. The results demonstrated that the suggested approach is feasible for offering high-resolution DRC maps that correspond to distinct microstructures with a limited amount of optimised acquisition data from two-dimensional DRC sampling space. By significantly reducing scan time while retaining structural resolution, this approach may enable multi-dimensional DRC MRI to be more widely used for quantitative evaluation in biological and medical settings.

Influence of Magnetic Field Strength on Intravoxel Incoherent Motion Parameters in Diffusion MRI of the Calf

Background: The purpose of this study was to investigate the dependence of Intravoxel Incoherent Motion (IVIM) parameters measured in the human calf on B0. Methods: Diffusion-weighted image data of eight healthy volunteers were acquired using five b-values (0-600 s/mm2) at rest and after muscle activation at 0.55 and 7 T. The musculus gastrocnemius mediale (GM, activated) was assessed. The perfusion fraction f and diffusion coefficient D were determined using segmented fits. The dependence on field strength was assessed using Student's t-test for paired samples and the Wilcoxon signed-rank test. A biophysical model built on the three non-exchanging compartments of muscle, venous blood, and arterial blood was used to interpret the data using literature relaxation times. Results: The measured perfusion fraction of the GM was significantly lower at 7 T, both for the baseline measurement and after muscle activation. For 0.55 and 7 T, the mean f values were 7.59% and 3.63% at rest, and 14.03% and 6.92% after activation, respectively. The biophysical model estimations for the mean proton-density-weighted perfusion fraction were 3.37% and 6.50% for the non-activated and activated states, respectively. Conclusions: B0 may have a significant effect on the measured IVIM parameters. The blood relaxation times suggest that 7 T IVIM may be arterial-weighted whereas 0.55 T IVIM may exhibit an approximately equal weighting of arterial and venous blood.

An open relaxation-diffusion MRI dataset in neurosurgical studies

Diffusion MRI (dMRI) is a safe and noninvasive technique that provides insight into the microarchitecture of brain tissue. Relaxation-diffusion MRI (rdMRI) is an extension of traditional dMRI that captures diffusion imaging data at multiple TEs to detect tissue heterogeneity between relaxation and diffusivity. rdMRI has great potential in neurosurgical research including brain tumor grading and treatment response evaluation. However, the lack of available data has limited the exploration of rdMRI in clinical settings. To address this, we are sharing a high-quality rdMRI dataset from 18 neurosurgical patients with different types of lesions, as well as two healthy individuals as controls. The rdMRI data was acquired using 7 TEs, where at each TE multi-shell dMRI with high spatial and angular resolutions is obtained at each TE. Each rdMRI scan underwent thorough artifact and distortion corrections using a specially designed processing pipeline. The dataset's quality was assessed using standard practices, including quality control and assurance. This resource is a valuable addition to neurosurgical studies, and all data are openly accessible.

Neuronal microstructural changes in the human brain are associated with neurocognitive aging

Gray matter (GM) alterations play a role in aging-related disorders like Alzheimer's disease and related dementias, yet MRI studies mainly focus on macroscopic changes. Although reliable indicators of atrophy, morphological metrics like cortical thickness lack the sensitivity to detect early changes preceding visible atrophy. Our study aimed at exploring the potential of diffusion MRI in unveiling sensitive markers of cortical and subcortical age-related microstructural changes and assessing their associations with cognitive and behavioral deficits. We leveraged the Human Connectome Project-Aging cohort that included 707 unimpaired participants (394 female; median age = 58, range = 36-90 years) and applied the powerful mean apparent diffusion propagator model to measure microstructural parameters, along with comprehensive behavioral and cognitive test scores. Both macro- and microstructural GM characteristics were strongly associated with age, with widespread significant microstructural correlations reflective of cellular morphological changes, reduced cellular density, increased extracellular volume, and increased membrane permeability. Importantly, when correlating MRI and cognitive test scores, our findings revealed no link between macrostructural volumetric changes and neurobehavioral performance. However, we found that cellular and extracellular alterations in cortical and subcortical GM regions were associated with neurobehavioral performance. Based on these findings, it is hypothesized that increased microstructural heterogeneity and decreased neurite orientation dispersion precede macrostructural changes, and that they play an important role in subsequent cognitive decline. These alterations are suggested to be early markers of neurocognitive performance that may distinctly aid in identifying the mechanisms underlying phenotypic aging and subsequent age-related functional decline.

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

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

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

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

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.

Quantitative MRI protocol and decision model for a 'one stop shop' early-stage Parkinsonism diagnosis: Study design

Differentiating among early-stage parkinsonisms is a challenge in clinical practice. Quantitative MRI can aid the diagnostic process, but studies with singular MRI techniques have had limited success thus far. Our objective is to develop a multi-modal MRI method for this purpose. In this review we describe existing methods and present a dedicated quantitative MRI protocol, a decision model and a study design to validate our approach ahead of a pilot study. We present example imaging data from patients and a healthy control, which resemble related literature.

Differentiation of white matter histopathology using b-tensor encoding and machine learning

Diffusion-weighted magnetic resonance imaging (DW-MRI) is a non-invasive technique that is sensitive to microstructural geometry in neural tissue and is useful for the detection of neuropathology in research and clinical settings. Tensor-valued diffusion encoding schemes (b-tensor) have been developed to enrich the microstructural data that can be obtained through DW-MRI. These advanced methods have proven to be more specific to microstructural properties than conventional DW-MRI acquisitions. Additionally, machine learning methods are particularly useful for the study of multidimensional data sets. In this work, we have tested the reach of b-tensor encoding data analyses with machine learning in different histopathological scenarios. We achieved this in three steps: 1) We induced different levels of white matter damage in rodent optic nerves. 2) We obtained ex vivo DW-MRI data with b-tensor encoding schemes and calculated quantitative metrics using Q-space trajectory imaging. 3) We used a machine learning model to identify the main contributing features and built a voxel-wise probabilistic classification map of histological damage. Our results show that this model is sensitive to characteristics of microstructural damage. In conclusion, b-tensor encoded DW-MRI data analyzed with machine learning methods, have the potential to be further developed for the detection of histopathology and neurodegeneration.

A novel framework for in-vivo diffusion tensor distribution MRI of the human brain

Neural tissue microstructure plays an important role in developmental, physiological and pathophysiological processes. Diffusion tensor distribution (DTD) MRI helps probe subvoxel heterogeneity by describing water diffusion within a voxel using an ensemble of non-exchanging compartments characterized by a probability density function of diffusion tensors. In this study, we provide a new framework for acquiring multiple diffusion encoding (MDE) images and estimating DTD from them in the human brain in vivo. We interfused pulsed field gradients (iPFG) in a single spin echo to generate arbitrary b-tensors of rank one, two, or three without introducing concomitant gradient artifacts. Employing well-defined diffusion encoding parameters we show that iPFG retains salient features of a traditional multiple-PFG (mPFG/MDE) sequence while reducing the echo time and coherence pathway artifacts thereby extending its applications beyond DTD MRI. Our DTD is a maximum entropy tensor-variate normal distribution whose tensor random variables are constrained to be positive definite to ensure their physicality. In each voxel, the second-order mean and fourth-order covariance tensors of the DTD are estimated using a Monte Carlo method that synthesizes micro-diffusion tensors with corresponding size, shape, and orientation distributions to best fit the measured MDE images. From these tensors we obtain the spectrum of diffusion tensor ellipsoid sizes and shapes, and the microscopic orientation distribution function (μODF) and microscopic fractional anisotropy (μFA) that disentangle the underlying heterogeneity within a voxel. Using the DTD-derived μODF, we introduce a new method to perform fiber tractography capable of resolving complex fiber configurations. The results revealed microscopic anisotropy in various gray and white matter regions and skewed MD distributions in cerebellar gray matter not observed previously. DTD MRI tractography captured complex white matter fiber organization consistent with known anatomy. DTD MRI also resolved some degeneracies associated with diffusion tensor imaging (DTI) and elucidated the source of diffusion heterogeneity which may help improve the diagnosis of various neurological diseases and disorders.

Using mesoscopic tract-tracing data to guide the estimation of fiber orientation distributions in the mouse brain from diffusion MRI

Diffusion MRI (dMRI) tractography is the only tool for non-invasive mapping of macroscopic structural connectivity over the entire brain. Although it has been successfully used to reconstruct large white matter tracts in the human and animal brains, the sensitivity and specificity of dMRI tractography remained limited. In particular, the fiber orientation distributions (FODs) estimated from dMRI signals, key to tractography, may deviate from histologically measured fiber orientation in crossing fibers and gray matter regions. In this study, we demonstrated that a deep learning network, trained using mesoscopic tract-tracing data from the Allen Mouse Brain Connectivity Atlas, was able to improve the estimation of FODs from mouse brain dMRI data. Tractography results based on the network generated FODs showed improved specificity while maintaining sensitivity comparable to results based on FOD estimated using a conventional spherical deconvolution method. Our result is a proof-of-concept of how mesoscale tract-tracing data can guide dMRI tractography and enhance our ability to characterize brain connectivity.

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

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

Multidimensional encoding of restricted and anisotropic diffusion by double rotation of the q vector

Diffusion NMR and MRI methods building on the classic pulsed gradient spin-echo sequence are sensitive to many aspects of translational motion, including time and frequency dependence ("restriction"), anisotropy, and flow, leading to ambiguities when interpreting experimental data from complex heterogeneous materials such as living biological tissues. While the oscillating gradient technique specifically targets frequency dependence and permits control of the sensitivity to flow, tensor-valued encoding enables investigations of anisotropy in orientationally disordered materials. Here, we propose a simple scheme derived from the "double-rotation" technique in solid-state NMR to generate a family of modulated gradient waveforms allowing for comprehensive exploration of the 2D frequency-anisotropy space and convenient investigation of both restricted and anisotropic diffusion with a single multidimensional acquisition protocol, thereby combining the desirable characteristics of the oscillating gradient and tensor-valued encoding techniques. The method is demonstrated by measuring multicomponent isotropic Gaussian diffusion in simple liquids, anisotropic Gaussian diffusion in a polydomain lyotropic liquid crystal, and restricted diffusion in a yeast cell sediment.

Profiling specific cell populations within the inflammatory tumor microenvironment by oscillating-gradient diffusion-weighted MRI

BACKGROUND

The inflammatory tumor microenvironment (TME) is formed by various immune cells, being closely associated with tumorigenesis. Especially, the interaction between tumor-infiltrating T-cells and macrophages has a crucial impact on tumor progression and metastatic spread. The purpose of this study was to investigate whether oscillating-gradient diffusion-weighted MRI (OGSE-DWI) enables a cell size-based discrimination between different cell populations of the TME.

METHODS

Sine-shaped OGSE-DWI was combined with the Imaging Microstructural Parameters Using Limited Spectrally Edited Diffusion (IMPULSED) approach to measure microscale diffusion distances, here relating to cell sizes. The accuracy of IMPULSED-derived cell radii was evaluated using in vitro spheroid models, consisting of either pure cancer cells, macrophages, or T-cells. Subsequently, in vivo experiments aimed to assess changes within the TME and its specific immune cell composition in syngeneic murine breast cancer models with divergent degrees of malignancy (4T1, 67NR) during tumor progression, clodronate liposome-mediated depletion of macrophages, and immune checkpoint inhibitor (ICI) treatment. Ex vivo analysis of IMPULSED-derived cell radii was conducted by immunohistochemical wheat germ agglutinin staining of cell membranes, while intratumoral immune cell composition was analyzed by CD3 and F4/80 co-staining.

RESULTS

OGSE-DWI detected mean cell radii of 8.8±1.3 µm for 4T1, 8.2±1.4 µm for 67NR, 13.0±1.7 for macrophage, and 3.8±1.8 µm for T-cell spheroids. While T-cell infiltration during progression of 4T1 tumors was observed by decreasing mean cell radii from 9.7±1.0 to 5.0±1.5 µm, increasing amount of intratumoral macrophages during progression of 67NR tumors resulted in increasing mean cell radii from 8.9±1.2 to 12.5±1.1 µm. After macrophage depletion, mean cell radii decreased from 6.3±1.7 to 4.4±0.5 µm. T-cell infiltration after ICI treatment was captured by decreasing mean cell radii in both tumor models, with more pronounced effects in the 67NR tumor model.

CONCLUSIONS

OGSE-DWI provides a versatile tool for non-invasive profiling of the inflammatory TME by assessing the dominating cell type T-cells or macrophages.

HYDI-DSI revisited: Constrained non-parametric EAP imaging without q-space re-gridding

Hybrid Diffusion Imaging (HYDI) was one of the first attempts to use multi-shell samplings of the q-space to infer diffusion properties beyond Diffusion Tensor Imaging (DTI) or High Angular Resolution Diffusion Imaging (HARDI). HYDI was intended as a flexible protocol embedding both DTI (for lower b-values) and HARDI (for higher b-values) processing, as well as Diffusion Spectrum Imaging (DSI) when the entire data set was exploited. In the latter case, the spherical sampling of the q-space is re-gridded by interpolation to a Cartesian lattice whose extent covers the range of acquired b-values, hence being acquisition-dependent. The Discrete Fourier Transform (DFT) is afterwards used to compute the corresponding Cartesian sampling of the Ensemble Average Propagator (EAP) in an entirely non-parametric way. From this lattice, diffusion markers such as the Return To Origin Probability (RTOP) or the Mean Squared Displacement (MSD) can be numerically estimated. We aim at re-formulating this scheme by means of a Fourier Transform encoding matrix that eliminates the need for q-space re-gridding at the same time it preserves the non-parametric nature of HYDI-DSI. The encoding matrix is adaptively designed at each voxel according to the underlying DTI approximation, so that an optimal sampling of the EAP can be pursued without being conditioned by the particular acquisition protocol. The estimation of the EAP is afterwards carried out as a regularized Quadratic Programming (QP) problem, which allows to impose positivity constraints that cannot be trivially embedded within the conventional HYDI-DSI. We demonstrate that the definition of the encoding matrix in the adaptive space allows to analytically (as opposed to numerically) compute several popular descriptors of diffusion with the unique source of error being the cropping of high frequency harmonics in the Fourier analysis of the attenuation signal. They include not only RTOP and MSD, but also Return to Axis/Plane Probabilities (RTAP/RTPP), which are defined in terms of specific spatial directions and are not available with the former HYDI-DSI. We report extensive experiments that suggest the benefits of our proposal in terms of accuracy, robustness and computational efficiency, especially when only standard, non-dedicated q-space samplings are available.

Metabolic activity diffusion imaging (MADI): II. Noninvasive, high-resolution human brain mapping of sodium pump flux and cell metrics

We introduce a new ¹ H₂ O magnetic resonance approach: metabolic activity diffusion imaging (MADI). Numerical diffusion-weighted imaging decay simulations characterized by the mean cellular water efflux (unidirectional) rate constant (k_io ), mean cell volume (V), and cell number density (ρ) are produced from Monte Carlo random walks in virtual stochastically sized/shaped cell ensembles. Because of active steady-state trans-membrane water cycling (AWC), k_io reflects the cytolemmal Na⁺ , K⁺ ATPase (NKA) homeostatic cellular metabolic rate (^c MR_NKA ). A digital 3D "library" contains thousands of simulated single diffusion-encoded (SDE) decays. Library entries match well with disparate, animal, and human experimental SDE decays. The V and ρ values are consistent with estimates from pertinent in vitro cytometric and ex vivo histopathological literature: in vivo V and ρ values were previously unavailable. The library allows noniterative pixel-by-pixel experimental SDE decay library matchings that can be used to advantage. They yield proof-of-concept MADI parametric mappings of the awake, resting human brain. These reflect the tissue morphology seen in conventional MRI. While V is larger in gray matter (GM) than in white matter (WM), the reverse is true for ρ. Many brain structures have k_io values too large for current, invasive methods. For example, the median WM k_io is 22s^-1 ; likely reflecting mostly exchange within myelin. The k_io •V product map displays brain tissue ^c MR_NKA variation. The GM activity correlates, quantitatively and qualitatively, with the analogous resting-state brain ¹⁸ FDG-PET tissue glucose consumption rate (^t MR_glucose ) map; but noninvasively, with higher spatial resolution, and no pharmacokinetic requirement. The cortex, thalamus, putamen, and caudate exhibit elevated metabolic activity. MADI accuracy and precision are assessed. The results are contextualized with literature overall homeostatic brain glucose consumption and ATP production/consumption measures. The MADI/PET results suggest different GM and WM metabolic pathways. Preliminary human prostate results are also presented.

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.

Metabolic activity diffusion imaging (MADI): I. Metabolic, cytometric modeling and simulations

Evidence mounts that the steady-state cellular water efflux (unidirectional) first-order rate constant (k_io [s^-1 ]) magnitude reflects the ongoing, cellular metabolic rate of the cytolemmal Na⁺ , K⁺ -ATPase (NKA), ^c MR_NKA (pmol [ATP consumed by NKA]/s/cell), perhaps biology's most vital enzyme. Optimal ¹ H₂ O MR k_io determinations require paramagnetic contrast agents (CAs) in model systems. However, results suggest that the homeostatic metabolic k_io biomarker magnitude in vivo is often too large to be reached with allowable or possible CA living tissue distributions. Thus, we seek a noninvasive (CA-free) method to determine k_io in vivo. Because membrane water permeability has long been considered important in tissue water diffusion, we turn to the well-known diffusion-weighted MRI (DWI) modality. To analyze the diffusion tensor magnitude, we use a parsimoniously primitive model featuring Monte Carlo simulations of water diffusion in virtual ensembles comprising water-filled and -immersed randomly sized/shaped contracted Voronoi cells. We find this requires two additional, cytometric properties: the mean cell volume (V [pL]) and the cell number density (ρ [cells/μL]), important biomarkers in their own right. We call this approach metabolic activity diffusion imaging (MADI). We simulate water molecule displacements and transverse MR signal decays covering the entirety of b-space from pure water (ρ = V = 0; k_io undefined; diffusion coefficient, D₀ ) to zero diffusion. The MADI model confirms that, in compartmented spaces with semipermeable boundaries, diffusion cannot be described as Gaussian: the nanoscopic D (D_n ) is diffusion time-dependent, a manifestation of the "diffusion dispersion". When the "well-mixed" (steady-state) condition is reached, diffusion becomes limited, mainly by the probabilities of (1) encountering (ρ, V), and (2) permeating (k_io ) cytoplasmic membranes, and less so by D_n magnitudes. Importantly, for spaces with large area/volume (A/V; claustrophobia) ratios, this can happen in less than a millisecond. The model matches literature experimental data well, with implications for DWI interpretations.

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

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

Spatiotemporal encoding MRI using subspace-constrained sampling and locally-low-rank regularization: Applications to diffusion weighted and diffusion kurtosis imaging of human brain and prostate

The benefits of performing locally low-rank (LLR) reconstructions on subsampled diffusion weighted and diffusion kurtosis imaging data employing spatiotemporal encoding (SPEN) methods, is investigated. SPEN allows for self-referenced correction of motion-induced phase errors in case of interleaved diffusion-oriented acquisitions, and allows one to overcome distortions otherwise observed along EPI's phase-encoded dimension. In combination with LLR-based reconstructions of the pooled imaging data and with a joint subsampling of b-weighted and interleaved images, additional improvements in terms of sensitivity as well as shortened acquisition times are demonstrated, without noticeable penalties. Details on how the LLR-regularized, subspace-constrained image reconstructions were adapted to SPEN are given; the improvements introduced by adopting these reconstruction frameworks for the accelerated acquisition of diffusivity and of kurtosis imaging data in both relatively homogeneous regions like the human brain and in more challenging regions like the human prostate, are presented and discussed within the context of similar efforts in the field.

High-resolution cortical MAP-MRI reveals areal borders and laminar substructures observed with histological staining

The variations in cellular composition and tissue architecture measured with histology provide the biological basis for partitioning the brain into distinct cytoarchitectonic areas and for characterizing neuropathological tissue alterations. Clearly, there is an urgent need to develop whole-brain neuroradiological methods that can assess cortical cyto- and myeloarchitectonic features non-invasively. Mean apparent propagator (MAP) MRI is a clinically feasible diffusion MRI method that quantifies efficiently and comprehensively the net microscopic displacements of water molecules diffusing in tissues. We investigate the sensitivity of high-resolution MAP-MRI to detecting areal and laminar variations in cortical cytoarchitecture and compare our results with observations from corresponding histological sections in the entire brain of a rhesus macaque monkey. High-resolution images of MAP-derived parameters, in particular the propagator anisotropy (PA), non-gaussianity (NG), and the return-to-axis probability (RTAP) reveal cortical area-specific lamination patterns in good agreement with the corresponding histological stained sections. In a few regions, the MAP parameters provide superior contrast to the five histological stains used in this study, delineating more clearly boundaries and transition regions between cortical areas and laminar substructures. Throughout the cortex, various MAP parameters can be used to delineate transition regions between specific cortical areas observed with histology and to refine areal boundaries estimated using atlas registration-based cortical parcellation. Using surface-based analysis of MAP parameters we quantify the cortical depth dependence of diffusion propagators in multiple regions-of-interest in a consistent and rigorous manner that is largely independent of the cortical folding geometry. The ability to assess cortical cytoarchitectonic features efficiently and non-invasively, its clinical feasibility, and translatability make high-resolution MAP-MRI a promising 3D imaging tool for studying whole-brain cortical organization, characterizing abnormal cortical development, improving early diagnosis of neurodegenerative diseases, identifying targets for biopsies, and complementing neuropathological investigations.

Advanced spinal cord MRI in multiple sclerosis: Current techniques and future directions

Spinal cord magnetic resonance imaging (MRI) has a central role in multiple sclerosis (MS) clinical practice for diagnosis and disease monitoring. Advanced MRI sequences capable of visualizing and quantifying tissue macro- and microstructure and reflecting different pathological disease processes have been used in MS research; however, the spinal cord remains under-explored, partly due to technical obstacles inherent to imaging this structure. We propose that the study of the spinal cord merits equal ambition in overcoming technical challenges, and that there is much information to be exploited to make valuable contributions to our understanding of MS. We present a narrative review on the latest progress in advanced spinal cord MRI in MS, covering in the first part structural, functional, metabolic and vascular imaging methods. We focus on recent studies of MS and those making significant technical steps, noting the challenges that remain to be addressed and what stands to be gained from such advances. Throughout we also refer to other works that presend more in-depth review on specific themes. In the second part, we present several topics that, in our view, hold particular potential. The need for better imaging of gray matter is discussed. We stress the importance of developing imaging beyond the cervical spinal cord, and explore the use of ultra-high field MRI. Finally, some recommendations are given for future research, from study design to newer developments in analysis, and the need for harmonization of sequences and methods within the field. This review is aimed at researchers and clinicians with an interest in gaining an overview of the current state of advanced MRI research in this field and what is primed to be the future of spinal cord imaging in MS research.

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

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

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.

Histopathological modeling of status epilepticus-induced brain damage based on in vivo diffusion tensor imaging in rats

Non-invasive magnetic resonance imaging (MRI) methods have proved useful in the diagnosis and prognosis of neurodegenerative diseases. However, the interpretation of imaging outcomes in terms of tissue pathology is still challenging. This study goes beyond the current interpretation of in vivo diffusion tensor imaging (DTI) by constructing multivariate models of quantitative tissue microstructure in status epilepticus (SE)-induced brain damage. We performed in vivo DTI and histology in rats at 79 days after SE and control animals. The analyses focused on the corpus callosum, hippocampal subfield CA3b, and layers V and VI of the parietal cortex. Comparison between control and SE rats indicated that a combination of microstructural tissue changes occurring after SE, such as cellularity, organization of myelinated axons, and/or morphology of astrocytes, affect DTI parameters. Subsequently, we constructed a multivariate regression model for explaining and predicting histological parameters based on DTI. The model revealed that DTI predicted well the organization of myelinated axons (cross-validated R = 0.876) and astrocyte processes (cross-validated R = 0.909) and possessed a predictive value for cell density (CD) (cross-validated R = 0.489). However, the morphology of astrocytes (cross-validated R > 0.05) was not well predicted. The inclusion of parameters from CA3b was necessary for modeling histopathology. Moreover, the multivariate DTI model explained better histological parameters than any univariate model. In conclusion, we demonstrate that combining several analytical and statistical tools can help interpret imaging outcomes to microstructural tissue changes, opening new avenues to improve the non-invasive diagnosis and prognosis of brain tissue damage.

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.

Diffusion time dependence, power-law scaling, and exchange in gray matter

Characterizing neural tissue microstructure is a critical goal for future neuroimaging. Diffusion MRI (dMRI) provides contrasts that reflect diffusing spins' interactions with myriad microstructural features of biological systems. However, the specificity of dMRI remains limited due to the ambiguity of its signals vis-à-vis the underlying microstructure. To improve specificity, biophysical models of white matter (WM) typically express dMRI signals according to the Standard Model (SM) and have more recently in gray matter (GM) taken spherical compartments into account (the SANDI model) in attempts to represent cell soma. The validity of the assumptions underlying these models, however, remains largely undetermined, especially in GM. To validate these assumptions experimentally, observing their unique, functional properties, such as the b-1/2 power-law associated with one-dimensional diffusion, has emerged as a fruitful strategy. The absence of this signature in GM, in turn, has been explained by neurite water exchange, non-linear morphology, and/or by obscuring soma signal contributions. Here, we present diffusion simulations in realistic neurons demonstrating that curvature and branching does not destroy the stick power-law behavior in impermeable neurites, but also that their signal is drowned by the soma signal under typical experimental conditions. Nevertheless, by studying the GM dMRI signal's behavior as a function of diffusion weighting as well as time, we identify an attainable experimental regime in which the neurite signal dominates. Furthermore, we find that exchange-driven time dependence produces a signal behavior opposite to that which would be expected from restricted diffusion, thereby providing a functional signature that disambiguates the two effects. We present data from dMRI experiments in ex vivo rat brain at ultrahigh field of 16.4T and observe a time dependence that is consistent with substantial exchange but also with a GM stick power-law. The first finding suggests significant water exchange between neurites and the extracellular space while the second suggests a small sub-population of impermeable neurites. To quantify these observations, we harness the Kärger exchange model and incorporate the corresponding signal time dependence in the SM and SANDI models.

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

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

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.

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.

Advanced Diffusion MR Imaging for Multiple Sclerosis in the Brain and Spinal Cord

Diffusion tensor imaging (DTI) has been established its usefulness in evaluating normal-appearing white matter (NAWM) and other lesions that are difficult to evaluate with routine clinical MRI in the evaluation of the brain and spinal cord lesions in multiple sclerosis (MS), a demyelinating disease. With the recent advances in the software and hardware of MRI systems, increasingly complex and sophisticated MRI and analysis methods, such as q-space imaging, diffusional kurtosis imaging, neurite orientation dispersion and density imaging, white matter tract integrity, and multiple diffusion encoding, referred to as advanced diffusion MRI, have been proposed. These are capable of capturing in vivo microstructural changes in the brain and spinal cord in normal and pathological states in greater detail than DTI.

This paper reviews the current status of recent advanced diffusion MRI for assessing MS in vivo as part of an issue celebrating two decades of magnetic resonance in medical sciences (MRMS), an official journal of the Japanese Society of Magnetic Resonance in Medicine.

3D Relaxation-Assisted Separation of Wideline Solid-State NMR Patterns for Achieving Site Resolution

There are currently no methods for the acquisition of ultra-wideline (UW) solid-state NMR spectra under static conditions that enable reliable separation and resolution of overlapping powder patterns arising from magnetically...

Optogenetic activation of striatal D1R and D2R cells differentially engages downstream connected areas beyond the basal ganglia

The basal ganglia (BG) are a group of subcortical nuclei responsible for motor and executive function. Central to BG function are striatal cells expressing D1 (D1R) and D2 (D2R) dopamine receptors. D1R and D2R cells are considered functional antagonists that facilitate voluntary movements and inhibit competing motor patterns, respectively. However, whether they maintain a uniform function across the striatum and what influence they exert outside the BG is unclear. Here, we address these questions by combining optogenetic activation of D1R and D2R cells in the mouse ventrolateral caudoputamen with fMRI. Striatal D1R/D2R stimulation evokes distinct activity within the BG-thalamocortical network and differentially engages cerebellar and prefrontal regions. Computational modeling of effective connectivity confirms that changes in D1R/D2R output drive functional relationships between these regions. Our results suggest a complex functional organization of striatal D1R/D2R cells and hint toward an interconnected fronto-BG-cerebellar network modulated by striatal D1R and D2R cells.

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.

Data-driven separation of MRI signal components for tissue characterization

PURPOSE

MRI can be utilized for quantitative characterization of tissue. To assess e.g. water fractions or diffusion coefficients for compartments in the brain, a decomposition of the signal is necessary. Imposing standard models carries the risk of estimating biased parameters if model assumptions are violated. This work introduces a data-driven multicomponent analysis, the monotonous slope non-negative matrix factorization (msNMF), tailored to extract data features expected in MR signals.

METHODS

The msNMF was implemented by extending the standard NMF with monotonicity constraints on the signal profiles and their first derivatives. The method was validated using simulated data, and subsequently applied to both ex vivo DWI data and in vivo relaxometry data. Reproducibility of the method was tested using the latter.

RESULTS

The msNMF recovered the multi-exponential signals in the simulated data and showed superiority to standard NMF (based on the explained variance, area under the ROC curve, and coefficient of variation). Diffusion components extracted from the DWI data reflected the cell density of the underlying tissue. The relaxometry analysis resulted in estimates of edema water fractions (EWF) highly correlated with published results, and demonstrated acceptable reproducibility.

CONCLUSION

The msNMF can robustly separate MR signals into components with relation to the underlying tissue composition, and may potentially be useful for e.g. tumor tissue characterization.

Neural networks for parameter estimation in microstructural MRI: Application to a diffusion-relaxation model of white matter

Specific features of white matter microstructure can be investigated by using biophysical models to interpret relaxation-diffusion MRI brain data. Although more intricate models have the potential to reveal more details of the tissue, they also incur time-consuming parameter estimation that may converge to inaccurate solutions due to a prevalence of local minima in a degenerate fitting landscape. Machine-learning fitting algorithms have been proposed to accelerate the parameter estimation and increase the robustness of the attained estimates. So far, learning-based fitting approaches have been restricted to microstructural models with a reduced number of independent model parameters where dense sets of training data are easy to generate. Moreover, the degree to which machine learning can alleviate the degeneracy problem is poorly understood. For conventional least-squares solvers, it has been shown that degeneracy can be avoided by acquisition with optimized relaxation-diffusion-correlation protocols that include tensor-valued diffusion encoding. Whether machine-learning techniques can offset these acquisition requirements remains to be tested. In this work, we employ artificial neural networks to vastly accelerate the parameter estimation for a recently introduced relaxation-diffusion model of white matter microstructure. We also develop strategies for assessing the accuracy and sensitivity of function fitting networks and use those strategies to explore the impact of the acquisition protocol. The developed learning-based fitting pipelines were tested on relaxation-diffusion data acquired with optimal and sub-optimal acquisition protocols. Networks trained with an optimized protocol were observed to provide accurate parameter estimates within short computational times. Comparing neural networks and least-squares solvers, we found the performance of the former to be less affected by sub-optimal protocols; however, model fitting networks were still susceptible to degeneracy issues and their use could not fully replace a careful design of the acquisition protocol.

High-resolution microscopic diffusion anisotropy imaging in the human hippocampus at 3T

PURPOSE

Several neurological conditions are associated with microstructural changes in the hippocampus that can be observed using DWI. Imaging studies often use protocols with whole-brain coverage, imposing limits on image resolution and worsening partial-volume effects. Also, conventional single-diffusion-encoding methods confound microscopic diffusion anisotropy with size variance of microscopic diffusion environments. This study addresses these issues by implementing a multidimensional diffusion-encoding protocol for microstructural imaging of the hippocampus at high resolution.

METHODS

The hippocampus of 8 healthy volunteers was imaged at 1.5-mm isotropic resolution with a multidimensional diffusion-encoding sequence developed in house. Microscopic fractional anisotropy (µFA) and normalized size variance (C_MD ) were estimated using q-space trajectory imaging, and their values were compared with DTI metrics. The overall scan time was 1 hour. The reproducibility of the protocol was confirmed with scan-rescan experiments, and a shorter protocol (14 minutes) was defined for situations with time constraints.

RESULTS

Mean µFA (0.47) was greater than mean FA (0.20), indicating orientation dispersion in hippocampal tissue microstructure. Mean C_MD was 0.17. The reproducibility of q-space trajectory imaging metrics was comparable to DTI, and microstructural metrics in the healthy hippocampus are reported.

CONCLUSION

This work shows the feasibility of high-resolution microscopic anisotropy imaging in the human hippocampus at 3 T and provides reference values for microstructural metrics in a healthy hippocampus.

Connectome 2.0: Developing the next-generation ultra-high gradient strength human MRI scanner for bridging studies of the micro-, meso- and macro-connectome

The first phase of the Human Connectome Project pioneered advances in MRI technology for mapping the macroscopic structural connections of the living human brain through the engineering of a whole-body human MRI scanner equipped with maximum gradient strength of 300 mT/m, the highest ever achieved for human imaging. While this instrument has made important contributions to the understanding of macroscale connectional topology, it has also demonstrated the potential of dedicated high-gradient performance scanners to provide unparalleled in vivo assessment of neural tissue microstructure. Building on the initial groundwork laid by the original Connectome scanner, we have now embarked on an international, multi-site effort to build the next-generation human 3T Connectome scanner (Connectome 2.0) optimized for the study of neural tissue microstructure and connectional anatomy across multiple length scales. In order to maximize the resolution of this in vivo microscope for studies of the living human brain, we will push the diffusion resolution limit to unprecedented levels by (1) nearly doubling the current maximum gradient strength from 300 mT/m to 500 mT/m and tripling the maximum slew rate from 200 T/m/s to 600 T/m/s through the design of a one-of-a-kind head gradient coil optimized to minimize peripheral nerve stimulation; (2) developing high-sensitivity multi-channel radiofrequency receive coils for in vivo and ex vivo human brain imaging; (3) incorporating dynamic field monitoring to minimize image distortions and artifacts; (4) developing new pulse sequences to integrate the strongest diffusion encoding and highest spatial resolution ever achieved in the living human brain; and (5) calibrating the measurements obtained from this next-generation instrument through systematic validation of diffusion microstructural metrics in high-fidelity phantoms and ex vivo brain tissue at progressively finer scales with accompanying diffusion simulations in histology-based micro-geometries. We envision creating the ultimate diffusion MRI instrument capable of capturing the complex multi-scale organization of the living human brain - from the microscopic scale needed to probe cellular geometry, heterogeneity and plasticity, to the mesoscopic scale for quantifying the distinctions in cortical structure and connectivity that define cyto- and myeloarchitectonic boundaries, to improvements in estimates of macroscopic connectivity.

Multidimensional MRI for Characterization of Subtle Axonal Injury Accelerated Using an Adaptive Nonlocal Multispectral Filter

Multidimensional MRI is an emerging approach that simultaneously encodes water relaxation (T₁ and T₂) and mobility (diffusion) and replaces voxel-averaged values with subvoxel distributions of those MR properties. While conventional (i.e., voxel-averaged) MRI methods cannot adequately quantify the microscopic heterogeneity of biological tissue, using subvoxel information allows to selectively map a specific T₁-T₂-diffusion spectral range that corresponds to a group of tissue elements. The major obstacle to the adoption of rich, multidimensional MRI protocols for diagnostic or monitoring purposes is the prolonged scan time. Our main goal in the present study is to evaluate the performance of a nonlocal estimation of multispectral magnitudes (NESMA) filter on reduced datasets to limit the total acquisition time required for reliable multidimensional MRI characterization of the brain. Here we focused and reprocessed results from a recent study that identified potential imaging biomarkers of axonal injury pathology from the joint analysis of multidimensional MRI, in particular voxelwise T₁-T₂ and diffusion-T₂ spectra in human Corpus Callosum, and histopathological data. We tested the performance of NESMA and its effect on the accuracy of the injury biomarker maps, relative to the co-registered histological reference. Noise reduction improved the accuracy of the resulting injury biomarker maps, while permitting data reduction of 35.7 and 59.6% from the full dataset for T₁-T₂ and diffusion-T₂ cases, respectively. As successful clinical proof-of-concept applications of multidimensional MRI are continuously being introduced, reliable and robust noise removal and consequent acquisition acceleration would advance the field towards clinically-feasible diagnostic multidimensional MRI protocols.

SPHERIOUSLY? The challenges of estimating sphere radius non-invasively in the human brain from diffusion MRI

The Soma and Neurite Density Imaging (SANDI) three-compartment model was recently proposed to disentangle cylindrical and spherical geometries, attributed to neurite and soma compartments, respectively, in brain tissue. There are some recent advances in diffusion-weighted MRI signal encoding and analysis (including the use of multiple so-called 'b-tensor' encodings and analysing the signal in the frequency-domain) that have not yet been applied in the context of SANDI. In this work, using: (i) ultra-strong gradients; (ii) a combination of linear, planar, and spherical b-tensor encodings; and (iii) analysing the signal in the frequency domain, three main challenges to robust estimation of sphere size were identified: First, the Rician noise floor in magnitude-reconstructed data biases estimates of sphere properties in a non-uniform fashion. It may cause overestimation or underestimation of the spherical compartment size and density. This can be partly ameliorated by accounting for the noise floor in the estimation routine. Second, even when using the strongest diffusion-encoding gradient strengths available for human MRI, there is an empirical lower bound on the spherical signal fraction and radius that can be detected and estimated robustly. For the experimental setup used here, the lower bound on the sphere signal fraction was approximately 10%. We employed two different ways of establishing the lower bound for spherical radius estimates in white matter. The first, examining power-law relationships between the DW-signal and diffusion weighting in empirical data, yielded a lower bound of 7μm, while the second, pure Monte Carlo simulations, yielded a lower limit of 3μm and in this low radii domain, there is little differentiation in signal attenuation. Third, if there is sensitivity to the transverse intra-cellular diffusivity in cylindrical structures, e.g., axons and cellular projections, then trying to disentangle two diffusion-time-dependencies using one experimental parameter (i.e., change in frequency-content of the encoding waveform) makes spherical radii estimates particularly challenging. We conclude that due to the aforementioned challenges spherical radii estimates may be biased when the corresponding sphere signal fraction is low, which must be considered.

Diffuse axonal injury has a characteristic multidimensional MRI signature in the human brain

Axonal injury is a major contributor to the clinical symptomatology in patients with traumatic brain injury. Conventional neuroradiological tools, such as CT and MRI, are insensitive to diffuse axonal injury (DAI) caused by trauma. Diffusion tensor MRI parameters may change in DAI lesions; however, the nature of these changes is inconsistent. Multidimensional MRI is an emerging approach that combines T1, T2, and diffusion, and replaces voxel-averaged values with distributions, which allows selective isolation of specific potential abnormal components. By performing a combined post-mortem multidimensional MRI and histopathology study, we aimed to investigate T1-T2-diffusion changes linked to DAI and to define their histopathological correlates. Corpora callosa derived from eight subjects who had sustained traumatic brain injury, and three control brain donors underwent post-mortem ex vivo MRI at 7 T. Multidimensional, diffusion tensor, and quantitative T1 and T2 MRI data were acquired and processed. Following MRI acquisition, slices from the same tissue were tested for amyloid precursor protein (APP) immunoreactivity to define DAI severity. A robust image co-registration method was applied to accurately match MRI-derived parameters and histopathology, after which 12 regions of interest per tissue block were selected based on APP density, but blind to MRI. We identified abnormal multidimensional T1-T2, diffusion-T2, and diffusion-T1 components that are strongly associated with DAI and used them to generate axonal injury images. We found that compared to control white matter, mild and severe DAI lesions contained significantly larger abnormal T1-T2 component (P = 0.005 and P < 0.001, respectively), and significantly larger abnormal diffusion-T2 component (P = 0.005 and P < 0.001, respectively). Furthermore, within patients with traumatic brain injury the multidimensional MRI biomarkers differentiated normal-appearing white matter from mild and severe DAI lesions, with significantly larger abnormal T1-T2 and diffusion-T2 components (P = 0.003 and P < 0.001, respectively, for T1-T2; P = 0.022 and P < 0.001, respectively, for diffusion-T2). Conversely, none of the conventional quantitative MRI parameters were able to differentiate lesions and normal-appearing white matter. Lastly, we found that the abnormal T1-T2, diffusion-T1, and diffusion-T2 components and their axonal damage images were strongly correlated with quantitative APP staining (r = 0.876, P < 0.001; r = 0.727, P < 0.001; and r = 0.743, P < 0.001, respectively), while producing negligible intensities in grey matter and in normal-appearing white matter. These results suggest that multidimensional MRI may provide non-invasive biomarkers for detection of DAI, which is the pathological substrate for neurological disorders ranging from concussion to severe traumatic brain injury.