1
|
Cai TX, Williamson NH, Ravin R, Basser PJ. The Diffusion Exchange Ratio (DEXR): A minimal sampling of diffusion exchange spectroscopy to probe exchange, restriction, and time-dependence. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2024; 366:107745. [PMID: 39126819 DOI: 10.1016/j.jmr.2024.107745] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/05/2024] [Revised: 07/31/2024] [Accepted: 08/02/2024] [Indexed: 08/12/2024]
Abstract
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.
Collapse
Affiliation(s)
- Teddy X Cai
- Section on Quantitative Imaging and Tissue Sciences, Eunice Kennedy Shriver National Institute of Child Health and Human Development, Bethesda, 20892, MD, USA
| | - Nathan H Williamson
- Section on Quantitative Imaging and Tissue Sciences, Eunice Kennedy Shriver National Institute of Child Health and Human Development, Bethesda, 20892, MD, USA
| | - Rea Ravin
- Section on Quantitative Imaging and Tissue Sciences, Eunice Kennedy Shriver National Institute of Child Health and Human Development, Bethesda, 20892, MD, USA; Celoptics, Inc., Rockville, 20850, MD, USA
| | - Peter J Basser
- Section on Quantitative Imaging and Tissue Sciences, Eunice Kennedy Shriver National Institute of Child Health and Human Development, Bethesda, 20892, MD, USA.
| |
Collapse
|
2
|
Storey P, Novikov DS. Signatures of microstructure in gradient-echo and spin-echo signals. Magn Reson Med 2024; 92:269-288. [PMID: 38520259 PMCID: PMC11178261 DOI: 10.1002/mrm.30022] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/26/2023] [Revised: 01/06/2024] [Accepted: 01/08/2024] [Indexed: 03/25/2024]
Abstract
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.
Collapse
Affiliation(s)
- Pippa Storey
- Center for Biomedical Imaging, Department of Radiology, New York University School of Medicine, New York, New York, USA
| | - Dmitry S Novikov
- Center for Biomedical Imaging, Department of Radiology, New York University School of Medicine, New York, New York, USA
| |
Collapse
|
3
|
Fang C, Yang Z, Wassermann D, Li JR. A simulation-driven supervised learning framework to estimate brain microstructure using diffusion MRI. Med Image Anal 2023; 90:102979. [PMID: 37827109 DOI: 10.1016/j.media.2023.102979] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2023] [Revised: 09/13/2023] [Accepted: 09/22/2023] [Indexed: 10/14/2023]
Abstract
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.
Collapse
Affiliation(s)
- Chengran Fang
- INRIA Saclay, Equipe IDEFIX, UMA, ENSTA Paris, 828, Boulevard des Maréchaux, 91762 Palaiseau, France; INRIA Saclay, Equipe MIND, 1 Rue Honoré d'Estienne d'Orves, 91120 Palaiseau, France
| | - Zheyi Yang
- INRIA Saclay, Equipe IDEFIX, UMA, ENSTA Paris, 828, Boulevard des Maréchaux, 91762 Palaiseau, France
| | - Demian Wassermann
- INRIA Saclay, Equipe MIND, 1 Rue Honoré d'Estienne d'Orves, 91120 Palaiseau, France
| | - Jing-Rebecca Li
- INRIA Saclay, Equipe IDEFIX, UMA, ENSTA Paris, 828, Boulevard des Maréchaux, 91762 Palaiseau, France.
| |
Collapse
|
4
|
Diez-Gómez V, Sobrados I, Ruiz-Santaquiteria C, Bucheli W, Jiménez R, Sanz J. Self-diffusion in polycrystalline Li 1+x Ti 2-x Al x (PO 4) 3 (0.2 ≤ x ≤ 0.4) samples followed by 7Li PFG (pulse field gradient) NMR spectroscopy. RSC Adv 2023; 13:13040-13051. [PMID: 37124023 PMCID: PMC10133837 DOI: 10.1039/d3ra02094k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2023] [Accepted: 04/19/2023] [Indexed: 05/02/2023] Open
Abstract
Short and long range lithium motions in powder Li1+x Ti2-x Al x (PO4)3 (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 m2 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.
Collapse
Affiliation(s)
- Virginia Diez-Gómez
- Dpto. Energy, Materials Science Institute of Madrid, ICMM-CSIC Cantoblanco 28049 Madrid Spain
| | - Isabel Sobrados
- Dpto. Energy, Materials Science Institute of Madrid, ICMM-CSIC Cantoblanco 28049 Madrid Spain
| | | | - Wilmer Bucheli
- Dpto. Energy, Materials Science Institute of Madrid, ICMM-CSIC Cantoblanco 28049 Madrid Spain
| | - Ricardo Jiménez
- Dpto. Energy, Materials Science Institute of Madrid, ICMM-CSIC Cantoblanco 28049 Madrid Spain
- Dpto. Materials for the Information Technologies, Materials Science Institute of Madrid, ICMM-CSIC Cantoblanco 28049 Madrid Spain
| | - Jesús Sanz
- Dpto. Energy, Materials Science Institute of Madrid, ICMM-CSIC Cantoblanco 28049 Madrid Spain
| |
Collapse
|
5
|
A new perspective of molecular diffusion by nuclear magnetic resonance. Sci Rep 2023; 13:1703. [PMID: 36717666 PMCID: PMC9887074 DOI: 10.1038/s41598-023-27389-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2022] [Accepted: 01/02/2023] [Indexed: 01/31/2023] Open
Abstract
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.
Collapse
|
6
|
Cai TX, Williamson NH, Ravin R, Basser PJ. Disentangling the effects of restriction and exchange with diffusion exchange spectroscopy. FRONTIERS IN PHYSICS 2022; 10:805793. [PMID: 37063496 PMCID: PMC10104504 DOI: 10.3389/fphy.2022.805793] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/17/2023]
Abstract
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.
Collapse
Affiliation(s)
- Teddy X. Cai
- Section on Quantitative Imaging and Tissue Sciences, Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health, Bethesda, Maryland, USA
- Wellcome Centre for Integrative Neuroimaging, FMRIB, Nuffield Department of Clinical Neurosciences, University of Oxford, Oxford, UK
| | - Nathan H. Williamson
- Section on Quantitative Imaging and Tissue Sciences, Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health, Bethesda, Maryland, USA
- National Institute of General Medical Sciences, National Institutes of Health, Bethesda, Maryland, USA
| | - Rea Ravin
- Section on Quantitative Imaging and Tissue Sciences, Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health, Bethesda, Maryland, USA
- Celoptics, Rockville, Maryland, USA
| | - Peter J. Basser
- Section on Quantitative Imaging and Tissue Sciences, Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health, Bethesda, Maryland, USA
- Correspondence: Peter J. Basser, Section on Quantitative Imaging and Tissue Sciences, Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health, Building 13, Room 3W16, 13 South Drive, Bethesda, Maryland 20892-5772, USA,
| |
Collapse
|