The localization regime in a nutshell

The Diffusion Exchange Ratio (DEXR): A minimal sampling of diffusion exchange spectroscopy to probe exchange, restriction, and time-dependence

Water exchange is increasingly recognized as an important biological process that can affect the study of biological tissue using diffusion MR. Methods to measure exchange, however, remain immature as opposed to those used to characterize restriction, with no consensus on the optimal pulse sequence (s) or signal model (s). In general, the trend has been towards data-intensive fitting of highly parameterized models. We take the opposite approach and show that a judicious sub-sample of diffusion exchange spectroscopy (DEXSY) data can be used to robustly quantify exchange, as well as restriction, in a data-efficient manner. This sampling produces a ratio of two points per mixing time: (i) one point with equal diffusion weighting in both encoding periods, which gives maximal exchange contrast, and (ii) one point with the same total diffusion weighting in just the first encoding period, for normalization. We call this quotient the Diffusion EXchange Ratio (DEXR). Furthermore, we show that it can be used to probe time-dependent diffusion by estimating the velocity autocorrelation function (VACF) over intermediate to long times (∼2-500ms). We provide a comprehensive theoretical framework for the design of DEXR experiments in the case of static or constant gradients. Data from Monte Carlo simulations and experiments acquired in fixed and viable ex vivo neonatal mouse spinal cord using a permanent magnet system are presented to test and validate this approach. In viable spinal cord, we report the following apparent parameters from just 6 data points: τk=17±4ms, fNG=0.72±0.01, Reff=1.05±0.01μm, and κeff=0.19±0.04μm/ms, which correspond to the exchange time, restricted or non-Gaussian signal fraction, an effective spherical radius, and permeability, respectively. For the VACF, we report a long-time, power-law scaling with ≈t-2.4, which is approximately consistent with disordered domains in 3-D. Overall, the DEXR method is shown to be highly efficient, capable of providing valuable quantitative diffusion metrics using minimal MR data.

Signatures of microstructure in gradient-echo and spin-echo signals

PURPOSE

To determine whether the spatial scale and magnetic susceptibility of microstructure can be evaluated robustly from the decay of gradient-echo and spin-echo signals.

THEORY AND METHODS

Gradient-echo and spin-echo images were acquired from suspensions of spherical polystyrene microbeads of 10, 20, and 40 μm nominal diameter. The sizes of the beads and their magnetic susceptibility relative to the medium were estimated from the signal decay curves, using a lookup table generated from Monte Carlo simulations and an analytic model based on the Gaussian phase approximation.

RESULTS

Fitting Monte Carlo predictions to spin-echo data yielded acceptable estimates of microstructural parameters for the 20 and 40 μm microbeads. Using gradient-echo data, the Monte Carlo lookup table provided satisfactory parameter estimates for the 20 μm beads but unstable results for the diameter of the largest beads. Neither spin-echo nor gradient-echo data allowed accurate parameter estimation for the smallest beads. The analytic model performed poorly over all bead sizes.

CONCLUSIONS

Microstructural sources of magnetic susceptibility produce distinctive non-exponential signatures in the decay of gradient-echo and spin-echo signals. However, inverting the problem to extract microstructural parameters from the signals is nontrivial and, in certain regimes, ill-conditioned. For microstructure with small characteristic length scales, parameter estimation is hampered by the difficulty of acquiring accurate data at very short echo times. For microstructure with large characteristic lengths, the gradient-echo signal approaches the static-dephasing regime, where it becomes insensitive to size. Applicability of the analytic model was further limited by failure of the Gaussian phase approximation for all but the smallest beads.

A simulation-driven supervised learning framework to estimate brain microstructure using diffusion MRI

We propose a framework to train supervised learning models on synthetic data to estimate brain microstructure parameters using diffusion magnetic resonance imaging (dMRI). Although further validation is necessary, the proposed framework aims to seamlessly incorporate realistic simulations into dMRI microstructure estimation. Synthetic data were generated from over 1,000 neuron meshes converted from digital neuronal reconstructions and linked to their neuroanatomical parameters (such as soma volume and neurite length) using an optimized diffusion MRI simulator that produces intracellular dMRI signals from the solution of the Bloch-Torrey partial differential equation. By combining random subsets of simulated neuron signals with a free diffusion compartment signal, we constructed a synthetic dataset containing dMRI signals and 40 tissue microstructure parameters of 1.45 million artificial brain voxels. To implement supervised learning models we chose multilayer perceptrons (MLPs) and trained them on a subset of the synthetic dataset to estimate some microstructure parameters, namely, the volume fractions of soma, neurites, and the free diffusion compartment, as well as the area fractions of soma and neurites. The trained MLPs perform satisfactorily on the synthetic test sets and give promising in-vivo parameter maps on the MGH Connectome Diffusion Microstructure Dataset (CDMD). Most importantly, the estimated volume fractions showed low dependence on the diffusion time, the diffusion time independence of the estimated parameters being a desired property of quantitative microstructure imaging. The synthetic dataset we generated will be valuable for the validation of models that map between the dMRI signals and microstructure parameters. The surface meshes and microstructures parameters of the aforementioned neurons have been made publicly available.

Self-diffusion in polycrystalline Li1+x Ti2-x Al x (PO4)3 (0.2 ≤ x ≤ 0.4) samples followed by 7Li PFG (pulse field gradient) NMR spectroscopy

Short and long range lithium motions in powder Li_1+x Ti_2-x Al _x (PO₄)₃ (LTAP) NASICON compounds prepared by ceramic (x = 0.2 and 0.4) and sol-gel (x = 0.3 and 0.4) routes are discussed. ND diffraction and MAS-NMR spectroscopy were previously used to investigate structural features of these compounds. In particular, Fourier map differences showed that the amount of Li atoms allocated at M3 increases at the expense of M1 sites when the Li content increases. In this work, PFG-NMR results show that diffusion coefficients rise with the amount of lithium and temperature. The restricted diffusion inside NASICON particles is compared with "free" diffusion processes. At 300 K, diffusion coefficients D _PFG ∼ 5 × 10^-12 m² s^-1 have been deduced in ceramic x = 0.2 and 0.4 samples, decreasing with diffusion time Δ used in PFG experiments. In sol-gel samples, diffusion coefficients are near those of ceramic samples, but decrease faster with diffusion Δ times, as a consequence of the Li confinement inside sub-micrometric crystallites. The NMR spin-echo signal displays minima at specific q(γgδ) values that are related to the crystallite size. From R _dif ∼ q _m ^-1 distances, calculated from the position of minima, and from diffusion coefficients deduced for high Δ values, the mean crystallite size was estimated. Finally, from the temperature dependence of conductivity and diffusion coefficients, the activation energy and charge carriers concentrations were determined.

A new perspective of molecular diffusion by nuclear magnetic resonance

The diffusion-weighted NMR signal acquired using Pulse Field Gradient (PFG) techniques, allows for extrapolating microstructural information from porous materials and biological tissues. In recent years there has been a multiplication of diffusion models expressed by parametric functions to fit the experimental data. However, clear-cut criteria for the model selection are lacking. In this paper, we develop a theoretical framework for the interpretation of NMR attenuation signals in the case of Gaussian systems with stationary increments. The full expression of the Stejskal-Tanner formula for normal diffusing systems is devised, together with its extension to the domain of anomalous diffusion. The range of applicability of the relevant parametric functions to fit the PFG data can be fully determined by means of appropriate checks to ascertain the correctness of the fit. Furthermore, the exact expression for diffusion weighted NMR signals pertaining to Brownian yet non-Gaussian processes is also derived, accompanied by the proper check to establish its contextual relevance. The analysis provided is particularly useful in the context of medical MRI and clinical practise where the hardware limitations do not allow the use of narrow pulse gradients.

Disentangling the effects of restriction and exchange with diffusion exchange spectroscopy

Diffusion exchange spectroscopy (DEXSY) is a multidimensional NMR technique that can reveal how water molecules exchange between compartments within heterogeneous media, such as biological tissue. Data from DEXSY experiments is typically processed using numerical inverse Laplace transforms (ILTs) to produce a diffusion-diffusion spectrum. A tacit assumption of this ILT approach is that the signal behavior is Gaussian - i.e., the spin echo intensity decays exponentially with the degree of diffusion weighting. The assumptions that underlie Gaussian signal behavior may be violated, however, depending on the gradient strength applied and the sample under study. We argue that non-Gaussian signal behavior due to restrictions is to be expected in the study of biological tissue using diffusion NMR. Further, we argue that this signal behavior can produce confounding features in the diffusion-diffusion spectra obtained from numerical ILTs of DEXSY data - entangling the effects of restriction and exchange. Specifically, restricted signal behavior can result in broadening of peaks and in the appearance of illusory exchanging compartments with distributed diffusivities, which pearl into multiple peaks if not highly regularized. We demonstrate these effects on simulated data. That said, we suggest the use of features in the signal acquisition domain that can be used to rapidly probe exchange without employing an ILT. We also propose a means to characterize the non-Gaussian signal behavior due to restrictions within a sample using DEXSY measurements with a near zero mixing time or storage interval. We propose a combined acquisition scheme to independently characterize restriction and exchange with various DEXSY measurements, which we term Restriction and Exchange from Equally-weighted Double and Single Diffusion Encodings (REEDS-DE). We test this method on ex vivo neonatal mouse spinal cord - a sample consisting primarily of gray matter - using a low-field, static gradient NMR system. In sum, we highlight critical shortcomings of prevailing DEXSY analysis methods that conflate the effects of restriction and exchange, and suggest a viable experimental approach to disentangle them.