Microstructure Imaging of Crossing (MIX) White Matter Fibers from diffusion MRI

Assessment of training-associated changes of the lumbar back muscle using a multiparametric MRI protocol

Adaptations in muscle physiology due to long-term physical training have been monitored using various methods: ranging from invasive techniques, such as biopsy, to less invasive approaches, such as electromyography (EMG), to various quantitative magnetic resonance imaging (qMRI) parameters. Typically, these latter parameters are assessed immediately after exercise. In contrast, this work assesses such adaptations in a set of qMRI parameters obtained at rest in the lumbar spine muscles of volunteers. To this end, we developed a multiparametric measurement protocol to extract quantitative values of (water) T2, fat fraction, T1, and Intra Voxel Incoherent Motion (IVIM) diffusion parameters in the lumbar back muscle. The protocol was applied to 31 healthy subjects divided into three differently trained cohorts: two groups of athletes (endurance athletes and powerlifters) and a control group with a sedentary lifestyle. Significant differences in muscle water T2, fat fraction, and pseudo-diffusion coefficient linked to microcirculatory blood flow in muscle tissue were found between the trained and untrained cohorts. At the same time, diffusion coefficients (resolved along different directions) provided additional differentiation between the two groups of athletes. Specifically, the strength-trained athletes showed lower axial and higher radial diffusion components compared to the endurance-trained cohort, which may indicate muscle hypertrophy. In conclusion, utilizing multiparametric information revealed new insights into the potential of quantitative MR parameters to detect and quantify long-term effects associated with training in differently trained cohorts, even at rest.

MR Imaging in Ataxias: Consensus Recommendations by the Ataxia Global Initiative Working Group on MRI Biomarkers

With many viable strategies in the therapeutic pipeline, upcoming clinical trials in hereditary and sporadic degenerative ataxias will benefit from non-invasive MRI biomarkers for patient stratification and the evaluation of therapies. The MRI Biomarkers Working Group of the Ataxia Global Initiative therefore devised guidelines to facilitate harmonized MRI data acquisition in clinical research and trials in ataxias. Recommendations are provided for a basic structural MRI protocol that can be used for clinical care and for an advanced multi-modal MRI protocol relevant for research and trial settings. The advanced protocol consists of modalities with demonstrated utility for tracking brain changes in degenerative ataxias and includes structural MRI, magnetic resonance spectroscopy, diffusion MRI, quantitative susceptibility mapping, and resting-state functional MRI. Acceptable ranges of acquisition parameters are provided to accommodate diverse scanner hardware in research and clinical contexts while maintaining a minimum standard of data quality. Important technical considerations in setting up an advanced multi-modal protocol are outlined, including the order of pulse sequences, and example software packages commonly used for data analysis are provided. Outcome measures most relevant for ataxias are highlighted with use cases from recent ataxia literature. Finally, to facilitate access to the recommendations by the ataxia clinical and research community, examples of datasets collected with the recommended parameters are provided and platform-specific protocols are shared via the Open Science Framework.

V-Mail: 3D-Enabled Correspondence About Spatial Data on (Almost) All Your Devices

We present V-Mail, a framework of cross-platform applications, interactive techniques, and communication protocols for improved multi-person correspondence about spatial 3D datasets. Inspired by the daily use of e-mail, V-Mail seeks to enable a similar style of rapid, multi-person communication accessible on any device; however, it aims to do this in the new context of spatial 3D communication, where limited access to 3D graphics hardware typically prevents such communication. The approach integrates visual data storytelling with data exploration, spatial annotations, and animated transitions. V-Mail "data stories" are exported in a standard video file format to establish a common baseline level of access on (almost) any device. The V-Mail framework also includes a series of complementary client applications and plugins that enable different degrees of story co-authoring and data exploration, adjusted automatically to match the capabilities of various devices. A lightweight, phone-based V-Mail app makes it possible to annotate data by adding captions to the video. These spatial annotations are then immediately accessible to team members running high-end 3D graphics visualization systems that also include a V-Mail client, implemented as a plugin. Results and evaluation from applying V-Mail to assist communication within an interdisciplinary science team studying Antarctic ice sheets confirm the utility of the asynchronous, cross-platform collaborative framework while also highlighting some current limitations and opportunities for future work.

Longitudinal Assessment of Intravoxel Incoherent Motion Diffusion-Weighted MRI Metrics in Cognitive Decline

BACKGROUND

Advanced diffusion-based MRI biomarkers may provide insight into microstructural and perfusion changes associated with neurodegeneration and cognitive decline.

PURPOSE

To assess longitudinal microstructural and perfusion changes using apparent diffusion coefficient (ADC) and intravoxel incoherent motion diffusion-weighted imaging (IVIM-DWI) parameters in cognitively impaired (CI) and healthy control (HC) groups.

STUDY TYPE

Prospective/longitudinal.

POPULATION

Twelve CI patients (75% female) and 13 HC subjects (69% female).

FIELD STRENGTH/SEQUENCE

3 T; Spin-Echo-IVIM-DWI.

ASSESSMENT

Two MRI scans were performed with a 12-month interval. ADC and IVIM-DWI metrics (diffusion coefficient [D] and perfusion fraction [f]) were generated from monoexponential and biexponential fits, respectively. Additionally, voxel-based correlations were evaluated between change in Montreal Cognitive Assessment (ΔMoCA) and baseline imaging parameters.

STATISTICAL TESTS

Analysis of covariance with sex and age as covariates was performed for main effects of group and time (false discovery rate [FDR] corrected) with post hoc comparisons using Bonferroni correction. Partial-η² and Hedges' g were used for effect-size analysis. Spearman's correlations (FDR corrected) were used for the relationship between ΔMoCA score and imaging. P < 0.05 was considered statistically significant.

RESULTS

Significant differences were found for the main effects of group (HC vs. CI) and time. For group effects, higher ADC, IVIM-D, and IVIM-f were observed in the CI group compared to HC (ADC: 1.23 ± 0.08^. 10^-3 vs. 1.09 ± 0.07^. 10^-3  mm² /sec; IVIM-D: 0.82 ± 0.01^. 10^-3 vs. 0.73 ± 0.01^. 10^-3  mm² /sec; and IVIM-f: 0.317 ± 0.008 vs. 0.253 ± 0.009). Significantly higher ADC, IVIM-D, and IVIM-f values were observed in the CI group after 12 months (ADC: 1.45 ± 0.05^. 10^-3 vs. 1.50 ± 0.07^. 10^-3  mm² /sec; IVIM-D: 0.87 ± 0.01^. 10^-3 vs. 0.94 ± 0.02^. 10^-3  mm² /sec; and IVIM-f: 0.303 ± 0.007 vs. 0.332 ± 0.008), but not in the HC group at large effect size. ADC, IVIM-D, and IVIM-f negatively correlated with ΔMoCA score (ρ = -0.49, -0.51, and -0.50, respectively).

DATA CONCLUSION

These findings demonstrate that longitudinal differences between CI and HC cohorts can be measured using IVIM-based metrics.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY STAGE: 2.

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.

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.

Bifurcated Topological Optimization for IVIM

In this work, we shed light on the issue of estimating Intravoxel Incoherent Motion (IVIM) for diffusion and perfusion estimation by characterizing the objective function using simplicial homology tools. We provide a robust solution via topological optimization of this model so that the estimates are more reliable and accurate. Estimating the tissue microstructure from diffusion MRI is in itself an ill-posed and a non-linear inverse problem. Using variable projection functional (VarPro) to fit the standard bi-exponential IVIM model we perform the optimization using simplicial homology based global optimization to better understand the topology of objective function surface. We theoretically show how the proposed methodology can recover the model parameters more accurately and consistently by casting it in a reduced subspace given by VarPro. Additionally we demonstrate that the IVIM model parameters cannot be accurately reconstructed using conventional numerical optimization methods due to the presence of infinite solutions in subspaces. The proposed method helps uncover multiple global minima by analyzing the local geometry of the model enabling the generation of reliable estimates of model parameters.

Bundle-Specific Axon Diameter Index as a New Contrast to Differentiate White Matter Tracts

In the central nervous system of primates, several pathways are characterized by different spectra of axon diameters. In vivo methods, based on diffusion-weighted magnetic resonance imaging, can provide axon diameter index estimates non-invasively. However, such methods report voxel-wise estimates, which vary from voxel-to-voxel for the same white matter bundle due to partial volume contributions from other pathways having different microstructure properties. Here, we propose a novel microstructure-informed tractography approach, COMMITAxSize, to resolve axon diameter index estimates at the streamline level, thus making the estimates invariant along trajectories. Compared to previously proposed voxel-wise methods, our formulation allows the estimation of a distinct axon diameter index value for each streamline, directly, furnishing a complementary measure to the existing calculation of the mean value along the bundle. We demonstrate the favourable performance of our approach comparing our estimates with existing histologically-derived measurements performed in the corpus callosum and the posterior limb of the internal capsule. Overall, our method provides a more robust estimation of the axon diameter index of pathways by jointly estimating the microstructure properties of the tissue and the macroscopic organisation of the white matter connectivity.

Decoding the microstructural properties of white matter using realistic models

Multi-echo gradient echo (ME-GRE) magnetic resonance signal evolution in white matter has a strong dependence on the orientation of myelinated axons with respect to the main static field. Although analytical solutions have been able to predict some of the white matter (WM) signal behaviour of the hollow cylinder model, it has been shown that realistic models of WM offer a better description of the signal behaviour observed. In this work, we present a pipeline to (i) generate realistic 2D WM models with their microstructure based on real axon morphology with adjustable fiber volume fraction (FVF) and g-ratio. We (ii) simulate their interaction with the static magnetic field to be able to simulate their MR signal. For the first time, we (iii) demonstrate that realistic 2D WM models can be used to simulate a MR signal that provides a good approximation of the signal obtained from a real 3D WM model derived from electron microscopy. We then (iv) demonstrate in silico that 2D WM models can be used to predict microstructural parameters in a robust way if ME-GRE multi-orientation data is available and the main fiber orientation in each pixel is known using DTI. A deep learning network was trained and characterized in its ability to recover the desired microstructural parameters such as FVF, g-ratio, free and bound water transverse relaxation and magnetic susceptibility. Finally, the network was trained to recover these micro-structural parameters from an ex vivo dataset acquired in 9 orientations with respect to the magnetic field and 12 echo times. We demonstrate that this is an overdetermined problem and that as few as 3 orientations can already provide comparable results for some of the decoded metrics.

Efficient whole-brain tract-specific T1 mapping at 3T with slice-shuffled inversion-recovery diffusion-weighted imaging

PURPOSE

Most voxels in white matter contain multiple fiber populations with different orientations and levels of myelination. Conventional T₁ mapping measures 1 T₁ value per voxel, representing a weighted average of the multiple tract T₁ times. Inversion-recovery diffusion-weighted imaging (IR-DWI) allows the T₁ times of multiple tracts in a voxel to be disentangled, but the scan time is prohibitively long. Recently, slice-shuffled IR-DWI implementations have been proposed to significantly reduce scan time. In this work, we demonstrate that we can measure tract-specific T₁ values in the whole brain using simultaneous multi-slice slice-shuffled IR-DWI at 3T.

METHODS

We perform simulations to evaluate the accuracy and precision of our crossing fiber IR-DWI signal model for various fiber parameters. The proposed sequence and signal model are tested in a phantom consisting of crossing asparagus pieces doped with gadolinium to vary T₁ , and in 2 human subjects.

RESULTS

Our simulations show that tract-specific T₁ times can be estimated within 5% of the nominal fiber T₁ values. Tract-specific T₁ values were resolved in subvoxel 2 fiber crossings in the asparagus phantom. Tract-specific T₁ times were resolved in 2 different tract crossings in the human brain where myelination differences have previously been reported; the crossing of the cingulum and genu of the corpus callosum and the crossing of the corticospinal tract and pontine fibers.

CONCLUSION

Whole-brain tract-specific T₁ mapping is feasible using slice-shuffled IR-DWI at 3T. This technique has the potential to improve the microstructural characterization of specific tracts implicated in neurodevelopment, aging, and demyelinating disorders.

Beyond the diffusion standard model in fixed rat spinal cord with combined linear and planar encoding

Information about tissue on the microscopic and mesoscopic scales can be accessed by modelling diffusion MRI signals, with the aim of extracting microstructure-specific biomarkers. The standard model (SM) of diffusion, currently the most broadly adopted microstructural model, describes diffusion in white matter (WM) tissues by two Gaussian components, one of which has zero radial diffusivity, to represent diffusion in intra- and extra-axonal water, respectively. Here, we reappraise these SM assumptions by collecting comprehensive double diffusion encoded (DDE) MRI data with both linear and planar encodings, which was recently shown to substantially enhance the ability to estimate SM parameters. We find however, that the SM is unable to account for data recorded in fixed rat spinal cord at an ultrahigh field of 16.4 T, suggesting that its underlying assumptions are violated in our experimental data. We offer three model extensions to mitigate this problem: first, we generalize the SM to accommodate finite radii (axons) by releasing the constraint of zero radial diffusivity in the intra-axonal compartment. Second, we include intracompartmental kurtosis to account for non-Gaussian behaviour. Third, we introduce an additional (third) compartment. The ability of these models to account for our experimental data are compared based on parameter feasibility and Bayesian information criterion. Our analysis identifies the three-compartment description as the optimal model. The third compartment exhibits slow diffusion with a minor but non-negligible signal fraction (∼12%). We demonstrate how failure to take the presence of such a compartment into account severely misguides inferences about WM microstructure. Our findings bear significance for microstructural modelling at large and can impact the interpretation of biomarkers extracted from the standard model of diffusion.

The sensitivity of diffusion MRI to microstructural properties and experimental factors

Diffusion MRI is a non-invasive technique to study brain microstructure. Differences in the microstructural properties of tissue, including size and anisotropy, can be represented in the signal if the appropriate method of acquisition is used. However, to depict the underlying properties, special care must be taken when designing the acquisition protocol as any changes in the procedure might impact on quantitative measurements. This work reviews state-of-the-art methods for studying brain microstructure using diffusion MRI and their sensitivity to microstructural differences and various experimental factors. Microstructural properties of the tissue at a micrometer scale can be linked to the diffusion signal at a millimeter-scale using modeling. In this paper, we first give an introduction to diffusion MRI and different encoding schemes. Then, signal representation-based methods and multi-compartment models are explained briefly. The sensitivity of the diffusion MRI signal to the microstructural components and the effects of curvedness of axonal trajectories on the diffusion signal are reviewed. Factors that impact on the quality (accuracy and precision) of derived metrics are then reviewed, including the impact of random noise, and variations in the acquisition parameters (i.e., number of sampled signals, b-value and number of acquisition shells). Finally, yet importantly, typical approaches to deal with experimental factors are depicted, including unbiased measures and harmonization. We conclude the review with some future directions and recommendations on this topic.

HARDI-ZOOMit protocol improves specificity to microstructural changes in presymptomatic myelopathy

Diffusion magnetic resonance imaging (dMRI) proved promising in patients with non-myelopathic degenerative cervical cord compression (NMDCCC), i.e., without clinically manifested myelopathy. Aim of the study is to present a fast multi-shell HARDI-ZOOMit dMRI protocol and validate its usability to detect microstructural myelopathy in NMDCCC patients. In 7 young healthy volunteers, 13 age-comparable healthy controls, 18 patients with mild NMDCCC and 15 patients with severe NMDCCC, the protocol provided higher signal-to-noise ratio, enhanced visualization of white/gray matter structures in microstructural maps, improved dMRI metric reproducibility, preserved sensitivity (SE = 87.88%) and increased specificity (SP = 92.31%) of control-patient group differences when compared to DTI-RESOLVE protocol (SE = 87.88%, SP = 76.92%). Of the 56 tested microstructural parameters, HARDI-ZOOMit yielded significant patient-control differences in 19 parameters, whereas in DTI-RESOLVE data, differences were observed in 10 parameters, with mostly lower robustness. Novel marker the white-gray matter diffusivity gradient demonstrated the highest separation. HARDI-ZOOMit protocol detected larger number of crossing fibers (5–15% of voxels) with physiologically plausible orientations than DTI-RESOLVE protocol (0–8% of voxels). Crossings were detected in areas of dorsal horns and anterior white commissure. HARDI-ZOOMit protocol proved to be a sensitive and practical tool for clinical quantitative spinal cord imaging.

Challenges for biophysical modeling of microstructure

The biophysical modeling efforts in diffusion MRI have grown considerably over the past 25 years. In this review, we dwell on the various challenges along the journey of bringing a biophysical model from initial design to clinical implementation, identifying both hurdles that have been already overcome and outstanding issues. First, we describe the critical initial task of selecting which features of tissue microstructure can be estimated using a model and which acquisition protocol needs to be implemented to make the estimation possible. The model performance should necessarily be tested in realistic numerical simulations and in experimental data - adapting the fitting strategy accordingly, and parameter estimates should be validated against complementary techniques, when/if available. Secondly, the model performance and validity should be explored in pathological conditions, and, if appropriate, dedicated models for pathology should be developed. We build on examples from tumors, ischemia and demyelinating diseases. We then discuss the challenges associated with clinical translation and added value. Finally, we single out four major unresolved challenges that are related to: the availability of a microstructural ground truth, the validation of model parameters which cannot be accessed with complementary techniques, the development of a generalized standard model for any brain region and pathology, and the seamless communication between different parties involved in the development and application of biophysical models of diffusion.

Close Reading for Visualization Evaluation

Visualizations produced by collaborations between artists, scientists, and visualization experts lay claim to being not only more effective in delivering information but also more effective in their abilities to elicit qualities like human connection. However, as prior work in the visualization community has demonstrated, it is difficult to evaluate these claims because characteristics associated with human connection are not easily measured quantitatively. In this Visualization Viewpoints piece, we address this problem in the context of our work to develop methods of evaluating visualizations created by Sculpting Visualization, a multidisciplinary project that incorporates art and design theory and practice into the process of scientific visualization. We present the design and results of a study in which we used close reading, a formal methodology used by humanities scholars, as a way to test reactions and analyses from evaluation participants related to an image created using Sculpting Visualization. In addition to specific suggestions about how to improve future iterations of the visualization, we discuss key findings of the evaluation related to contextual information, visual perspective, and associations that individual viewers brought to bear on their experience with the visualization.

Diffusion Imaging in the Post HCP Era

Diffusion imaging is a critical component in the pursuit of developing a better understanding of the human brain. Recent technical advances promise enabling the advancement in the quality of data that can be obtained. In this review the context for different approaches relative to the Human Connectome Project are compared. Significant new gains are anticipated from the use of high-performance head gradients. These gains can be particularly large when the high-performance gradients are employed together with ultrahigh magnetic fields. Transmit array designs are critical in realizing high accelerations in diffusion-weighted (d)MRI acquisitions, while maintaining large field of view (FOV) coverage, and several techniques for optimal signal-encoding are now available. Reconstruction and processing pipelines that precisely disentangle the acquired neuroanatomical information are established and provide the foundation for the application of deep learning in the advancement of dMRI for complex tissues. Level of Evidence: 3 Technical Efficacy Stage: Stage 3.

Leveraging multi-shell diffusion for studies of brain development in youth and young adulthood

•

Multi-shell imaging sequences may improve sensitivity to developmental effects.

•

Models that leverage multi-shell information are often less sensitive to the confounding effects of motion.

•

Multi-shell sequences and models that leverage this data may be of particular utility for studying the developing brain.

Diffusion weighted imaging (DWI) has advanced our understanding of brain microstructure evolution over development. Recently, the use of multi-shell diffusion imaging sequences has coincided with advances in modeling the diffusion signal, such as Neurite Orientation Dispersion and Density Imaging (NODDI) and Laplacian-regularized Mean Apparent Propagator MRI (MAPL). However, the relative utility of recently-developed diffusion models for understanding brain maturation remains sparsely investigated. Additionally, despite evidence that motion artifact is a major confound for studies of development, the vulnerability of metrics derived from contemporary models to in-scanner motion has not been described. Accordingly, in a sample of 120 youth and young adults (ages 12–30) we evaluated metrics derived from diffusion tensor imaging (DTI), NODDI, and MAPL for associations with age and in-scanner head motion at multiple scales. Specifically, we examined mean white matter values, white matter tracts, white matter voxels, and connections in structural brain networks. Our results revealed that multi-shell diffusion imaging data can be leveraged to robustly characterize neurodevelopment, and demonstrate stronger age effects than equivalent single-shell data. Additionally, MAPL-derived metrics were less sensitive to the confounding effects of head motion. Our findings suggest that multi-shell imaging data and contemporary modeling techniques confer important advantages for studies of neurodevelopment.

Improved fibre dispersion estimation using b-tensor encoding

Measuring fibre dispersion in white matter with diffusion magnetic resonance imaging (MRI) is limited by an inherent degeneracy between fibre dispersion and microscopic diffusion anisotropy (i.e., the diffusion anisotropy expected for a single fibre orientation). This means that estimates of fibre dispersion rely on strong assumptions, such as constant microscopic anisotropy throughout the white matter or specific biophysical models. Here we present a simple approach for resolving this degeneracy using measurements that combine linear (conventional) and spherical tensor diffusion encoding. To test the accuracy of the fibre dispersion when our microstructural model is only an approximation of the true tissue structure, we simulate multi-compartment data and fit this with a single-compartment model. For such overly simplistic tissue assumptions, we show that the bias in fibre dispersion is greatly reduced (~5x) for single-shell linear and spherical tensor encoding data compared with single-shell or multi-shell conventional data. In in-vivo data we find a consistent estimate of fibre dispersion as we reduce the b-value from 3 to 1.5 ms/μm², increase the repetition time, increase the echo time, or increase the diffusion time. We conclude that the addition of spherical tensor encoded data to conventional linear tensor encoding data greatly reduces the sensitivity of the estimated fibre dispersion to the model assumptions of the tissue microstructure.

Fast and accurate initialization of the free-water imaging model parameters from multi-shell diffusion MRI

Cerebrospinal fluid partial volume effect is a known bias in the estimation of Diffusion Tensor Imaging (DTI) parameters from diffusion MRI data. The Free-Water Imaging model for diffusion MRI data adds a second compartment to the DTI model, which explicitly accounts for the signal contribution of extracellular free-water, such as cerebrospinal fluid. As a result the DTI parameters obtained through the free-water model are corrected for partial volume effects, and thus better represent tissue microstructure. In addition, the model estimates the fractional volume of free-water, and can be used to monitor changes in the extracellular space. Under certain assumptions, the model can be estimated from single-shell diffusion MRI data. However, by using data from multi-shell diffusion acquisitions, these assumptions can be relaxed, and the fit becomes more robust. Nevertheless, fitting the model to multi-shell data requires high computational cost, with a non-linear iterative minimization, which has to be initialized close enough to the global minimum to avoid local minima and to robustly estimate the model parameters. Here we investigate the properties of the main initialization approaches that are currently being used, and suggest new fast approaches to improve the initial estimates of the model parameters. We show that our proposed approaches provide a fast and accurate initial approximation of the model parameters, which is very close to the final solution. We demonstrate that the proposed initializations improve the final outcome of non-linear model fitting.

Artifact-Based Rendering: Harnessing Natural and Traditional Visual Media for More Expressive and Engaging 3D Visualizations

We introduce Artifact-Based Rendering (ABR), a framework of tools, algorithms, and processes that makes it possible to produce real, data-driven 3D scientific visualizations with a visual language derived entirely from colors, lines, textures, and forms created using traditional physical media or found in nature. A theory and process for ABR is presented to address three current needs: (i) designing better visualizations by making it possible for non-programmers to rapidly design and critique many alternative data-to-visual mappings; (ii) expanding the visual vocabulary used in scientific visualizations to depict increasingly complex multivariate data; (iii) bringing a more engaging, natural, and human-relatable handcrafted aesthetic to data visualization. New tools and algorithms to support ABR include front-end applets for constructing artifact-based colormaps, optimizing 3D scanned meshes for use in data visualization, and synthesizing textures from artifacts. These are complemented by an interactive rendering engine with custom algorithms and interfaces that demonstrate multiple new visual styles for depicting point, line, surface, and volume data. A within-the-research-team design study provides early evidence of the shift in visualization design processes that ABR is believed to enable when compared to traditional scientific visualization systems. Qualitative user feedback on applications to climate science and brain imaging support the utility of ABR for scientific discovery and public communication.

The Dmipy Toolbox: Diffusion MRI Multi-Compartment Modeling and Microstructure Recovery Made Easy

Non-invasive estimation of brain microstructure features using diffusion MRI (dMRI)—known as Microstructure Imaging—has become an increasingly diverse and complicated field over the last decades. Multi-compartment (MC)-models, representing the measured diffusion signal as a linear combination of signal models of distinct tissue types, have been developed in many forms to estimate these features. However, a generalized implementation of MC-modeling as a whole, providing deeper insights in its capabilities, remains missing. To address this fact, we present Diffusion Microstructure Imaging in Python (Dmipy), an open-source toolbox implementing PGSE-based MC-modeling in its most general form. Dmipy allows on-the-fly implementation, signal modeling, and optimization of any user-defined MC-model, for any PGSE acquisition scheme. Dmipy follows a “building block”-based philosophy to Microstructure Imaging, meaning MC-models are modularly constructed to include any number and type of tissue models, allowing simultaneous representation of a tissue's diffusivity, orientation, volume fractions, axon orientation dispersion, and axon diameter distribution. In particular, Dmipy is geared toward facilitating reproducible, reliable MC-modeling pipelines, often allowing the whole process from model construction to parameter map recovery in fewer than 10 lines of code. To demonstrate Dmipy's ease of use and potential, we implement a wide range of well-known MC-models, including IVIM, AxCaliber, NODDI(x), Bingham-NODDI, the spherical mean-based SMT and MC-MDI, and spherical convolution-based single- and multi-tissue CSD. By allowing parameter cascading between MC-models, Dmipy also facilitates implementation of advanced approaches like CSD with voxel-varying kernels and single-shell 3-tissue CSD. By providing a well-tested, user-friendly toolbox that simplifies the interaction with the otherwise complicated field of dMRI-based Microstructure Imaging, Dmipy contributes to more reproducible, high-quality research.

Optimizing the intrinsic parallel diffusivity in NODDI: An extensive empirical evaluation

Purpose

NODDI is widely used in parameterizing microstructural brain properties. The model includes three signal compartments: intracellular, extracellular, and free water. The neurite compartment intrinsic parallel diffusivity (d_∥) is set to 1.7 μm²⋅ms⁻¹, though the effects of this assumption have not been extensively explored. This work investigates the optimality of d_∥ = 1.7 μm²⋅ms⁻¹ under varying imaging protocol, age groups, sex, and tissue type in comparison to other biologically plausible values of d_∥.

Methods

Model residuals were used as the optimality criterion. The model residuals were evaluated in function of d_∥ over the range from 0.5 to 3.0 μm²⋅ms⁻¹. This was done with respect to tissue type (i.e., white matter versus gray matter), sex, age (infancy to late adulthood), and diffusion-weighting protocol (maximum b-value). Variation in the estimated parameters with respect to d_∥ was also explored.

Results

Results show d_∥ = 1.7 μm²⋅ms⁻¹ is appropriate for adult brain white matter but it is suboptimal for gray matter with optimal values being significantly lower. d_∥ = 1.7 μm²⋅ms⁻¹ was also suboptimal in the infant brain for both white and gray matter with optimal values being significantly lower. Minor optimum d_∥ differences were observed versus diffusion protocol. No significant sex effects were observed. Additionally, changes in d_∥ resulted in significant changes to the estimated NODDI parameters.

Conclusion

The default (d_∥) of 1.7 μm²⋅ms⁻¹ is suboptimal in gray matter and infant brains.

A deep network for tissue microstructure estimation using modified LSTM units

Diffusion magnetic resonance imaging (dMRI) offers a unique tool for noninvasively assessing tissue microstructure. However, accurate estimation of tissue microstructure described by complicated signal models can be challenging when a reduced number of diffusion gradients are used. Deep learning based microstructure estimation has recently been developed and achieved promising results. In particular, optimization-based learning, where deep network structures are constructed by unfolding the iterative processes performed for solving optimization problems, has demonstrated great potential in accurate microstructure estimation with a reduced number of diffusion gradients. In this work, using the optimization-based learning strategy, we propose a deep network structure that is motivated by the use of historical information in iterative optimization for tissue microstructure estimation, and such incorporation of historical information has not been previously explored in the design of deep networks for microstructure estimation. We assume that (1) diffusion signals can be sparsely represented by a dictionary and its coefficients jointly in the spatial and angular domain, and (2) tissue microstructure can be computed from the sparse representation. Following these assumptions, our network comprises two cascaded stages. The first stage takes image patches as input and computes the spatial-angular sparse representation of the input with learned weights. Specifically, the network structure in the first stage is constructed by unfolding an iterative process for solving sparse reconstruction problems, where historical information is incorporated. The components in this network can be shown to correspond to modified long short-term memory (LSTM) units. In the second stage, fully connected layers are added to compute the mapping from the sparse representation to tissue microstructure. The weights in the two stages are learned jointly by minimizing the mean squared error of microstructure estimation. Experiments were performed on dMRI scans with a reduced number of diffusion gradients. For demonstration, we evaluated the estimation of tissue microstructure described by three signal models: the neurite orientation dispersion and density imaging (NODDI) model, the spherical mean technique (SMT) model, and the ensemble average propagator (EAP) model. The results indicate that the proposed approach outperforms competing methods.

Imaging brain microstructure with diffusion MRI: practicality and applications

This article gives an overview of microstructure imaging of the brain with diffusion MRI and reviews the state of the art. The microstructure-imaging paradigm aims to estimate and map microscopic properties of tissue using a model that links these properties to the voxel scale MR signal. Imaging techniques of this type are just starting to make the transition from the technical research domain to wide application in biomedical studies. We focus here on the practicalities of both implementing such techniques and using them in applications. Specifically, the article summarizes the relevant aspects of brain microanatomy and the range of diffusion-weighted MR measurements that provide sensitivity to them. It then reviews the evolution of mathematical and computational models that relate the diffusion MR signal to brain tissue microstructure, as well as the expanding areas of application. Next we focus on practicalities of designing a working microstructure imaging technique: model selection, experiment design, parameter estimation, validation, and the pipeline of development of this class of technique. The article concludes with some future perspectives on opportunities in this topic and expectations on how the field will evolve in the short-to-medium term.

Towards microstructure fingerprinting: Estimation of tissue properties from a dictionary of Monte Carlo diffusion MRI simulations

Many closed-form analytical models have been proposed to relate the diffusion-weighted magnetic resonance imaging (DW-MRI) signal to microstructural features of white matter tissues. These models generally make assumptions about the tissue and the diffusion processes which often depart from the biophysical reality, limiting their reliability and interpretability in practice. Monte Carlo simulations of the random walk of water molecules are widely recognized to provide near groundtruth for DW-MRI signals. However, they have mostly been limited to the validation of simpler models rather than used for the estimation of microstructural properties. This work proposes a general framework which leverages Monte Carlo simulations for the estimation of physically interpretable microstructural parameters, both in single and in crossing fascicles of axons. Monte Carlo simulations of DW-MRI signals, or fingerprints, are pre-computed for a large collection of microstructural configurations. At every voxel, the microstructural parameters are estimated by optimizing a sparse combination of these fingerprints. Extensive synthetic experiments showed that our approach achieves accurate and robust estimates in the presence of noise and uncertainties over fixed or input parameters. In an in vivo rat model of spinal cord injury, our approach provided microstructural parameters that showed better correspondence with histology than five closed-form models of the diffusion signal: MMWMD, NODDI, DIAMOND, WMTI and MAPL. On whole-brain in vivo data from the human connectome project (HCP), our method exhibited spatial distributions of apparent axonal radius and axonal density indices in keeping with ex vivo studies. This work paves the way for microstructure fingerprinting with Monte Carlo simulations used directly at the modeling stage and not only as a validation tool.

Robust and Fast Markov Chain Monte Carlo Sampling of Diffusion MRI Microstructure Models

In diffusion MRI analysis, advances in biophysical multi-compartment modeling have gained popularity over the conventional Diffusion Tensor Imaging (DTI), because they can obtain a greater specificity in relating the dMRI signal to underlying cellular microstructure. Biophysical multi-compartment models require a parameter estimation, typically performed using either the Maximum Likelihood Estimation (MLE) or the Markov Chain Monte Carlo (MCMC) sampling. Whereas, the MLE provides only a point estimate of the fitted model parameters, the MCMC recovers the entire posterior distribution of the model parameters given in the data, providing additional information such as parameter uncertainty and correlations. MCMC sampling is currently not routinely applied in dMRI microstructure modeling, as it requires adjustment and tuning, specific to each model, particularly in the choice of proposal distributions, burn-in length, thinning, and the number of samples to store. In addition, sampling often takes at least an order of magnitude, more time than non-linear optimization. Here we investigate the performance of the MCMC algorithm variations over multiple popular diffusion microstructure models, to examine whether a single, well performing variation could be applied efficiently and robustly to many models. Using an efficient GPU-based implementation, we showed that run times can be removed as a prohibitive constraint for the sampling of diffusion multi-compartment models. Using this implementation, we investigated the effectiveness of different adaptive MCMC algorithms, burn-in, initialization, and thinning. Finally we applied the theory of the Effective Sample Size, to the diffusion multi-compartment models, as a way of determining a relatively general target for the number of samples needed to characterize parameter distributions for different models and data sets. We conclude that adaptive Metropolis methods increase MCMC performance and select the Adaptive Metropolis-Within-Gibbs (AMWG) algorithm as the primary method. We furthermore advise to initialize the sampling with an MLE point estimate, in which case 100 to 200 samples are sufficient as a burn-in. Finally, we advise against thinning in most use-cases and as a relatively general target for the number of samples, we recommend a multivariate Effective Sample Size of 2,200.

Cortical fibers orientation mapping using in-vivo whole brain 7 T diffusion MRI

Diffusion MRI of the cortical gray matter is challenging because the micro-environment probed by water molecules is much more complex than within the white matter. High spatial and angular resolutions are therefore necessary to uncover anisotropic diffusion patterns and laminar structures, which provide complementary (e.g. to anatomical and functional MRI) microstructural information about the cortex architectonic. Several ex-vivo and in-vivo MRI studies have recently addressed this question, however predominantly with an emphasis on specific cortical areas. There is currently no whole brain in-vivo data leveraging multi-shell diffusion MRI acquisition at high spatial resolution, and depth dependent analysis, to characterize the complex organization of cortical fibers. Here, we present unique in-vivo human 7T diffusion MRI data, and a dedicated cortical depth dependent analysis pipeline. We leverage the high spatial (1.05 mm isotropic) and angular (198 diffusion gradient directions) resolution of this whole brain dataset to improve cortical fiber orientations mapping, and study neurites (axons and/or dendrites) trajectories across cortical depths. Tangential fibers in superficial cortical depths and crossing fiber configurations in deep cortical depths are identified. Fibers gradually inserting into the gyral walls are visualized, which contributes to mitigating the gyral bias effect. Quantitative radiality maps and histograms in individual subjects and cortex-based aligned datasets further support our results.

Multidimensional encoding of brain connectomes

The ability to map brain networks in living individuals is fundamental in efforts to chart the relation between human behavior, health and disease. Advances in network neuroscience may benefit from developing new frameworks for mapping brain connectomes. We present a framework to encode structural brain connectomes and diffusion-weighted magnetic resonance (dMRI) data using multidimensional arrays. The framework integrates the relation between connectome nodes, edges, white matter fascicles and diffusion data. We demonstrate the utility of the framework for in vivo white matter mapping and anatomical computing by evaluating 1,490 connectomes, thirteen tractography methods, and three data sets. The framework dramatically reduces storage requirements for connectome evaluation methods, with up to 40x compression factors. Evaluation of multiple, diverse datasets demonstrates the importance of spatial resolution in dMRI. We measured large increases in connectome resolution as function of data spatial resolution (up to 52%). Moreover, we demonstrate that the framework allows performing anatomical manipulations on white matter tracts for statistical inference and to study the white matter geometrical organization. Finally, we provide open-source software implementing the method and data to reproduce the results.

Tissue microstructure estimation using a deep network inspired by a dictionary-based framework

Diffusion magnetic resonance imaging (dMRI) captures the anisotropic pattern of water displacement in the neuronal tissue and allows noninvasive investigation of the complex tissue microstructure. A number of biophysical models have been proposed to relate the tissue organization with the observed diffusion signals, so that the tissue microstructure can be inferred. The Neurite Orientation Dispersion and Density Imaging (NODDI) model has been a popular choice and has been widely used for many neuroscientific studies. It models the diffusion signal with three compartments that are characterized by distinct diffusion properties, and the parameters in the model describe tissue microstructure. In NODDI, these parameters are estimated in a maximum likelihood framework, where the nonlinear model fitting is computationally intensive. Therefore, efforts have been made to develop efficient and accurate algorithms for NODDI microstructure estimation, which is still an open problem. In this work, we propose a deep network based approach that performs end-to-end estimation of NODDI microstructure, which is named Microstructure Estimation using a Deep Network (MEDN). MEDN comprises two cascaded stages and is motivated by the AMICO algorithm, where the NODDI microstructure estimation is formulated in a dictionary-based framework. The first stage computes the coefficients of the dictionary. It resembles the solution to a sparse reconstruction problem, where the iterative process in conventional estimation approaches is unfolded and truncated, and the weights are learned instead of predetermined by the dictionary. In the second stage, microstructure properties are computed from the output of the first stage, which resembles the weighted sum of normalized dictionary coefficients in AMICO, and the weights are also learned. Because spatial consistency of diffusion signals can be used to reduce the effect of noise, we also propose MEDN+, which is an extended version of MEDN. MEDN+ allows incorporation of neighborhood information by inserting a stage with learned weights before the MEDN structure, where the diffusion signals in the neighborhood of a voxel are processed. The weights in MEDN or MEDN+ are jointly learned from training samples that are acquired with diffusion gradients densely sampling the q-space. We performed MEDN and MEDN+ on brain dMRI scans, where two shells each with 30 gradient directions were used, and measured their accuracy with respect to the gold standard. Results demonstrate that the proposed networks outperform the competing methods.

Evaluating g-ratio weighted changes in the corpus callosum as a function of age and sex

Recent years have seen a growing interest in relating MRI measurements to the structural-biophysical properties of white matter fibers. The fiber g-ratio, defined as the ratio between the inner and outer radii of the axon myelin sheath, is an important structural property of white matter, affecting signal conduction. Recently proposed modeling methods that use a combination of quantitative-MRI signals, enable a measurement of the fiber g-ratio in vivo. Here we use an MRI-based g-ratio estimation to observe the variance of the g-ratio within the corpus callosum, and evaluate sex and age related differences. To estimate the g-ratio we used a model (Stikov et al., 2011; Duval et al., 2017) based on two different WM microstructure parameters: the relative amounts of myelin (myelin volume fraction, MVF) and fibers (fiber volume fraction, FVF) in a voxel. We derived the FVF from the fractional anisotropy (FA), and estimated the MVF by using the lipid and macromolecular tissue volume (MTV), calculated from the proton density (Mezer et al., 2013). In comparison to other methods of estimating the MVF, MTV represents a stable parameter with a straightforward route of acquisition. To establish our model, we first compared histological MVF measurements (West et al., 2016) with the MRI derived MTV. We then implemented our model on a large database of 92 subjects (44 males), aged 7 to 81, in order to evaluate age and sex related changes within the corpus callosum. Our results show that the MTV provides a good estimation of MVF for calculating g-ratio, and produced values from the corpus callosum that correspond to those found in animals ex vivo and are close to the theoretical optimum, as well as to published in vivo data. Our results demonstrate that the MTV derived g-ratio provides a simple and reliable in vivo g-ratio-weighted (GR*) measurement in humans. In agreement with theoretical predictions, and unlike other tissue parameters measured with MRI, the g-ratio estimations were found to be relatively stable with age, and we found no support for a significant sexual dimorphism with age.