Characterizing magnetic resonance signal decay due to Gaussian diffusion: the path integral approach and a convenient computational method

MRI with generalized diffusion encoding reveals damaged white matter in patients previously hospitalized for COVID-19 and with persisting symptoms at follow-up

There is mounting evidence of the long-term effects of COVID-19 on the central nervous system, with patients experiencing diverse symptoms, often suggesting brain involvement. Conventional brain MRI of these patients shows unspecific patterns, with no clear connection of the symptomatology to brain tissue abnormalities, whereas diffusion tensor studies and volumetric analyses detect measurable changes in the brain after COVID-19. Diffusion MRI exploits the random motion of water molecules to achieve unique sensitivity to structures at the microscopic level, and new sequences employing generalized diffusion encoding provide structural information which are sensitive to intravoxel features. In this observational study, a total of 32 persons were investigated: 16 patients previously hospitalized for COVID-19 with persisting symptoms of post-COVID condition (mean age 60 years: range 41-79, all male) at 7-month follow-up and 16 matched controls, not previously hospitalized for COVID-19, with no post-COVID symptoms (mean age 58 years, range 46-69, 11 males). Standard MRI and generalized diffusion encoding MRI were employed to examine the brain white matter of the subjects. To detect possible group differences, several tissue microstructure descriptors obtainable with the employed diffusion sequence, the fractional anisotropy, mean diffusivity, axial diffusivity, radial diffusivity, microscopic anisotropy, orientational coherence (Cc) and variance in compartment's size (CMD) were analysed using the tract-based spatial statistics framework. The tract-based spatial statistics analysis showed widespread statistically significant differences (P < 0.05, corrected for multiple comparisons using the familywise error rate) in all the considered metrics in the white matter of the patients compared to the controls. Fractional anisotropy, microscopic anisotropy and Cc were lower in the patient group, while axial diffusivity, radial diffusivity, mean diffusivity and CMD were higher. Significant changes in fractional anisotropy, microscopic anisotropy and CMD affected approximately half of the analysed white matter voxels located across all brain lobes, while changes in Cc were mainly found in the occipital parts of the brain. Given the predominant alteration in microscopic anisotropy compared to Cc, the observed changes in diffusion anisotropy are mostly due to loss of local anisotropy, possibly connected to axonal damage, rather than white matter fibre coherence disruption. The increase in radial diffusivity is indicative of demyelination, while the changes in mean diffusivity and CMD are compatible with vasogenic oedema. In summary, these widespread alterations of white matter microstructure are indicative of vasogenic oedema, demyelination and axonal damage. These changes might be a contributing factor to the diversity of central nervous system symptoms that many patients experience after COVID-19.

Q-space trajectory imaging with positivity constraints (QTI+)

Q-space trajectory imaging (QTI) enables the estimation of useful scalar measures indicative of the local tissue structure. This is accomplished by employing generalized gradient waveforms for diffusion sensitization alongside a diffusion tensor distribution (DTD) model. The first two moments of the underlying DTD are made available by acquisitions at low diffusion sensitivity (b-values). Here, we show that three independent conditions have to be fulfilled by the mean and covariance tensors associated with distributions of symmetric positive semidefinite tensors. We introduce an estimation framework utilizing semi-definite programming (SDP) to guarantee that these conditions are met. Applying the framework on simulated signal profiles for diffusion tensors distributed according to non-central Wishart distributions demonstrates the improved noise resilience of QTI+ over the commonly employed estimation methods. Our findings on a human brain data set also reveal pronounced improvements, especially so for acquisition protocols featuring few number of volumes. Our method’s robustness to noise is expected to not only improve the accuracy of the estimates, but also enable a meaningful interpretation of contrast in the derived scalar maps. The technique’s performance on shorter acquisitions could make it feasible in routine clinical practice.

Multimodal integration of diffusion MRI for better characterization of tissue biology

The contrast in diffusion-weighted MR images is due to variations of diffusion properties within the examined specimen. Certain microstructural information on the underlying tissues can be inferred through quantitative analyses of the diffusion-sensitized MR signals. In the first part of the paper, we review two types of approach for characterizing diffusion MRI signals: Bloch's equations with diffusion terms, and statistical descriptions. Specifically, we discuss expansions in terms of cumulants and orthogonal basis functions, the confinement tensor formalism and tensor distribution models. Further insights into the tissue properties may be obtained by integrating diffusion MRI with other techniques, which is the subject of the second part of the paper. We review examples involving magnetic susceptibility, structural tensors, internal field gradients, transverse relaxation and functional MRI. Integrating information provided by other imaging modalities (MR based or otherwise) could be a key to improve our understanding of how diffusion MRI relates to physiology and biology.

Dynamics of local magnetization in the eigenbasis of the Bloch-Torrey operator

We consider diffusion within pores with general shapes in the presence of spatially linear magnetic field profiles. The evolution of local magnetization of the spin bearing particles can be described by the Bloch-Torrey equation. We study the diffusive process in the eigenbasis of the non-Hermitian Bloch-Torrey operator. It is possible to find expressions for some special temporal gradient waveforms employed to sensitize the nuclear magnetic resonance (NMR) signal to diffusion. For more general gradient waveforms, we derive an efficient numerical solution by introducing a novel matrix formalism. Compared to previous methods, this new approach requires a fewer number of eigenfunctions to achieve the same accuracy. This shows that these basis functions are better suited to the problem studied. The new framework could provide new important insights into the fundamentals of diffusion sensitization, which could further the development of the field of NMR.

Effective potential for magnetic resonance measurements of restricted diffusion

The signature of diffusive motion on the NMR signal has been exploited to characterize the mesoscopic structure of specimens in numerous applications. For compartmentalized specimens comprising isolated subdomains, a representation of individual pores is necessary for describing restricted diffusion within them. When gradient waveforms with long pulse durations are employed, a quadratic potential profile is identified as an effective energy landscape for restricted diffusion. The dependence of the stochastic effective force on the center-of-mass position is indeed found to be approximately linear (Hookean) for restricted diffusion even when the walls are sticky. We outline the theoretical basis and practical advantages of our picture involving effective potentials.

NMR signal for particles diffusing under potentials: From path integrals and numerical methods to a model of diffusion anisotropy

We study the influence of diffusion on NMR experiments when the molecules undergo random motion under the influence of a force field and place special emphasis on parabolic (Hookean) potentials. To this end, the problem is studied using path integral methods. Explicit relationships are derived for commonly employed gradient waveforms involving pulsed and oscillating gradients. The Bloch-Torrey equation, describing the temporal evolution of magnetization, is modified by incorporating potentials. A general solution to this equation is obtained for the case of parabolic potential by adopting the multiple correlation function (MCF) formalism, which has been used in the past to quantify the effects of restricted diffusion. Both analytical and MCF results were found to be in agreement with random walk simulations. A multidimensional formulation of the problem is introduced that leads to a new characterization of diffusion anisotropy. Unlike the case of traditional methods that employ a diffusion tensor, anisotropy originates from the tensorial force constant, and bulk diffusivity is retained in the formulation. Our findings suggest that some features of the NMR signal that have traditionally been attributed to restricted diffusion are accommodated by the Hookean model. Under certain conditions, the formalism can be envisioned to provide a viable approximation to the mathematically more challenging restricted diffusion problems.