Numerical study of a cylinder model of the diffusion MRI signal for neuronal dendrite trees

Diffusion-weighted MR spectroscopy: Consensus, recommendations, and resources from acquisition to modeling

Brain cell structure and function reflect neurodevelopment, plasticity, and aging; and changes can help flag pathological processes such as neurodegeneration and neuroinflammation. Accurate and quantitative methods to noninvasively disentangle cellular structural features are needed and are a substantial focus of brain research. Diffusion-weighted MRS (dMRS) gives access to diffusion properties of endogenous intracellular brain metabolites that are preferentially located inside specific brain cell populations. Despite its great potential, dMRS remains a challenging technique on all levels: from the data acquisition to the analysis, quantification, modeling, and interpretation of results. These challenges were the motivation behind the organization of the Lorentz Center workshop on "Best Practices & Tools for Diffusion MR Spectroscopy" held in Leiden, the Netherlands, in September 2021. During the workshop, the dMRS community established a set of recommendations to execute robust dMRS studies. This paper provides a description of the steps needed for acquiring, processing, fitting, and modeling dMRS data, and provides links to useful resources.

What does FEXI measure?

Filter-exchange imaging (FEXI) has already been utilized in several biomedical studies for evaluating the permeability of cell membranes. The method relies on suppressing the extracellular signal using strong diffusion weighting (the mobility filter causing a reduction in the overall diffusivity) and monitoring the subsequent diffusivity recovery. Using Monte Carlo simulations, we demonstrate that FEXI is sensitive not uniquely to the transcytolemmal exchange but also to the geometry of involved compartments: complex geometry offers locations where spins remain unaffected by the mobility filter; moving to other locations afterwards, such spins contribute to the diffusivity recovery without actually permeating any membrane. This exchange mechanism is a warning for those who aim to use FEXI in complex media such as brain gray matter and opens wide scope for investigation towards crystallizing the genuine membrane permeation and characterizing the compartment geometry.

Highlights from the Era of Open Source Web-Based Tools

High digital connectivity and a focus on reproducibility are contributing to an open science revolution in neuroscience. Repositories and platforms have emerged across the whole spectrum of subdisciplines, paving the way for a paradigm shift in the way we share, analyze, and reuse vast amounts of data collected across many laboratories. Here, we describe how open access web-based tools are changing the landscape and culture of neuroscience, highlighting six free resources that span subdisciplines from behavior to whole-brain mapping, circuits, neurons, and gene variants.

Practical computation of the diffusion MRI signal of realistic neurons based on Laplace eigenfunctions

The complex transverse water proton magnetization subject to diffusion-encoding magnetic field gradient pulses in a heterogeneous medium such as brain tissue can be modeled by the Bloch-Torrey partial differential equation. The spatial integral of the solution of this equation in realistic geometry provides a gold-standard reference model for the diffusion MRI signal arising from different tissue micro-structures of interest. A closed form representation of this reference diffusion MRI signal called matrix formalism, which makes explicit the link between the Laplace eigenvalues and eigenfunctions of the biological cell and its diffusion MRI signal, was derived 20 years ago. In addition, once the Laplace eigendecomposition has been computed and saved, the diffusion MRI signal can be calculated for arbitrary diffusion-encoding sequences and b-values at negligible additional cost. Up to now, this representation, though mathematically elegant, has not been often used as a practical model of the diffusion MRI signal, due to the difficulties of calculating the Laplace eigendecomposition in complicated geometries. In this paper, we present a simulation framework that we have implemented inside the MATLAB-based diffusion MRI simulator SpinDoctor that efficiently computes the matrix formalism representation for realistic neurons using the finite element method. We show that the matrix formalism representation requires a few hundred eigenmodes to match the reference signal computed by solving the Bloch-Torrey equation when the cell geometry originates from realistic neurons. As expected, the number of eigenmodes required to match the reference signal increases with smaller diffusion time and higher b-values. We also convert the eigenvalues to a length scale and illustrate the link between the length scale and the oscillation frequency of the eigenmode in the cell geometry. We give the transformation that links the Laplace eigenfunctions to the eigenfunctions of the Bloch-Torrey operator and compute the Bloch-Torrey eigenfunctions and eigenvalues. This work is another step in bringing advanced mathematical tools and numerical method development to the simulation and modeling of diffusion MRI.

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.

Diffusion MRI simulation in thin-layer and thin-tube media using a discretization on manifolds

The Bloch-Torrey partial differential equation can be used to describe the evolution of the transverse magnetization of the imaged sample under the influence of diffusion-encoding magnetic field gradients inside the MRI scanner. The integral of the magnetization inside a voxel gives the simulated diffusion MRI signal. This paper proposes a finite element discretization on manifolds in order to efficiently simulate the diffusion MRI signal in domains that have a thin layer or a thin tube geometrical structure. The variable thickness of the three-dimensional domains is included in the weak formulation established on the manifolds. We conducted a numerical study of the proposed approach by simulating the diffusion MRI signals from the extracellular space (a thin layer medium) and from neurons (a thin tube medium), comparing the results with the reference signals obtained using a standard three-dimensional finite element discretization. We show good agreements between the simulated signals using our proposed method and the reference signals for a wide range of diffusion MRI parameters. The approximation becomes better as the diffusion time increases. The method helps to significantly reduce the required simulation time, computational memory, and difficulties associated with mesh generation, thus opening the possibilities to simulating complicated structures at low cost for a better understanding of diffusion MRI in the brain.

Influence of the size and curvedness of neural projections on the orientationally averaged diffusion MR signal

Neuronal and glial projections can be envisioned to be tubes of infinitesimal diameter as far as diffusion magnetic resonance (MR) measurements via clinical scanners are concerned. Recent experimental studies indicate that the decay of the orientationally-averaged signal in white-matter may be characterized by the power-law, Ē(q) ∝ q^-1, where q is the wavenumber determined by the parameters of the pulsed field gradient measurements. One particular study by McKinnon et al. [1] reports a distinctively faster decay in gray-matter. Here, we assess the role of the size and curvature of the neurites and glial arborizations in these experimental findings. To this end, we studied the signal decay for diffusion along general curves at all three temporal regimes of the traditional pulsed field gradient measurements. We show that for curvy projections, employment of longer pulse durations leads to a disappearance of the q^-1 decay, while such decay is robust when narrow gradient pulses are used. Thus, in clinical acquisitions, the lack of such a decay for a fibrous specimen can be seen as indicative of fibers that are curved. We note that the above discussion is valid for an intermediate range of q-values as the true asymptotic behavior of the signal decay is Ē(q) ∝ q^-4 for narrow pulses (through Debye-Porod law) or steeper for longer pulses. This study is expected to provide insights for interpreting the diffusion-weighted images of the central nervous system and aid in the design of acquisition strategies.

Attenuation of the NMR signal in a field gradient due to stochastic dynamics with memory

The attenuation function S(t) for an ensemble of spins in a magnetic-field gradient is calculated by accumulation of the phase shifts in the rotating frame resulting from the displacements of spin-bearing particles. The found S(t), expressed through the particle mean square displacement, is applicable for any kind of stationary stochastic motion of spins, including their non-markovian dynamics with memory. The known expressions valid for normal and anomalous diffusion are obtained as special cases in the long time approximation. The method is also applicable to the NMR pulse sequences based on the refocusing principle. This is demonstrated by describing the Hahn spin echo experiment. The attenuation of the NMR signal is also evaluated providing that the random motion of particle is modeled by the generalized Langevin equation with the memory kernel exponentially decaying in time. The models considered in our paper assume massive particles driven by much smaller particles.

Win-win data sharing in neuroscience

Most neuroscientists have yet to embrace a culture of data sharing. Using our decade-long experience at NeuroMorpho.Org as an example, we discuss how publicly available repositories may benefit data producers and end-users alike. We outline practical recipes for resource developers to maximize the research impact of data sharing platforms for both contributors and users.

New paradigm to assess brain cell morphology by diffusion-weighted MR spectroscopy in vivo

The brain is one of the most complex organs, and tools are lacking to assess its cellular morphology in vivo. Here we combine original diffusion-weighted magnetic resonance (MR) spectroscopy acquisition and novel modeling strategies to explore the possibility of quantifying brain cell morphology noninvasively. First, the diffusion of cell-specific metabolites is measured at ultra-long diffusion times in the rodent and primate brain in vivo to observe how cell long-range morphology constrains metabolite diffusion. Massive simulations of particles diffusing in synthetic cells parameterized by morphometric statistics are then iterated to fit experimental data. This method yields synthetic cells (tentatively neurons and astrocytes) that exhibit striking qualitative and quantitative similarities with histology (e.g., using Sholl analysis). With our approach, we measure major interspecies difference regarding astrocytes, whereas dendritic organization appears better conserved throughout species. This work suggests that the time dependence of metabolite diffusion coefficient allows distinguishing and quantitatively characterizing brain cell morphologies noninvasively.

Sharing Neuron Data: Carrots, Sticks, and Digital Records

Routine data sharing is greatly benefiting several scientific disciplines, such as molecular biology, particle physics, and astronomy. Neuroscience data, in contrast, are still rarely shared, greatly limiting the potential for secondary discovery and the acceleration of research progress. Although the attitude toward data sharing is non-uniform across neuroscience subdomains, widespread adoption of data sharing practice will require a cultural shift in the community. Digital reconstructions of axonal and dendritic morphology constitute a particularly "sharable" kind of data. The popularity of the public repository NeuroMorpho.Org demonstrates that data sharing can benefit both users and contributors. Increased data availability is also catalyzing the grassroots development and spontaneous integration of complementary resources, research tools, and community initiatives. Even in this rare successful subfield, however, more data are still unshared than shared. Our experience as developers and curators of NeuroMorpho.Org suggests that greater transparency regarding the expectations and consequences of sharing (or not sharing) data, combined with public disclosure of which datasets are shared and which are not, may expedite the transition to community-wide data sharing.