1
|
Narvaez O, Yon M, Jiang H, Bernin D, Forssell-Aronsson E, Sierra A, Topgaard D. Nonparametric distributions of tensor-valued Lorentzian diffusion spectra for model-free data inversion in multidimensional diffusion MRI. J Chem Phys 2024; 161:084201. [PMID: 39171708 DOI: 10.1063/5.0213252] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2024] [Accepted: 07/09/2024] [Indexed: 08/23/2024] Open
Abstract
Magnetic resonance imaging (MRI) is the method of choice for noninvasive studies of micrometer-scale structures in biological tissues via their effects on the time- and frequency-dependent (restricted) and anisotropic self-diffusion of water. While new designs of time-dependent magnetic field gradient waveforms have enabled disambiguation between different aspects of translational motion that are convolved in traditional MRI methods relying on single pairs of field gradient pulses, data analysis for complex heterogeneous materials remains a challenge. Here, we propose and demonstrate nonparametric distributions of tensor-valued Lorentzian diffusion spectra, or "D(ω) distributions," as a general representation with sufficient flexibility to describe the MRI signal response from a wide range of model systems and biological tissues investigated with modulated gradient waveforms separating and correlating the effects of restricted and anisotropic diffusion.
Collapse
Affiliation(s)
- Omar Narvaez
- A.I.Virtanen Institute for Molecular Sciences, University of Eastern Finland, Kuopio, Finland
| | - Maxime Yon
- Department of Chemistry, Lund University, Lund, Sweden
| | - Hong Jiang
- Department of Chemistry, Lund University, Lund, Sweden
| | - Diana Bernin
- Department of Chemistry and Chemical Engineering, Chalmers University of Technology, Gothenburg, Sweden
| | - Eva Forssell-Aronsson
- Department of Medical Radiation Sciences, University of Gothenburg, Gothenburg, Sweden
- Medical Physics and Biomedical Engineering, Sahlgrenska University Hospital, Gothenburg, Sweden
- Sahlgrenska Center for Cancer Research, Sahlgrenska Academy at the University of Gothenburg, Gothenburg, Sweden
| | - Alejandra Sierra
- A.I.Virtanen Institute for Molecular Sciences, University of Eastern Finland, Kuopio, Finland
| | | |
Collapse
|
2
|
Stuprich CM, Loh M, Nemerth JT, Nagel AM, Uder M, Laun FB. Velocity-compensated intravoxel incoherent motion of the human calf muscle. Magn Reson Med 2024; 92:543-555. [PMID: 38688865 DOI: 10.1002/mrm.30059] [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: 10/19/2023] [Revised: 01/15/2024] [Accepted: 02/03/2024] [Indexed: 05/02/2024]
Abstract
PURPOSE To determine whether intravoxel incoherent motion (IVIM) describes the blood perfusion in muscles better, assuming pseudo diffusion (Bihan Model 1) or ballistic motion (Bihan Model 2). METHODS IVIM parameters were measured in 18 healthy subjects with three different diffusion gradient time profiles (bipolar with two diffusion times and one with velocity compensation) and 17 b-values (0-600 s/mm2) at rest and after muscle activation. The diffusion coefficient, perfusion fraction, and pseudo-diffusion coefficient were estimated with a segmented fit in the gastrocnemius medialis (GM) and tibialis anterior (TA) muscles. RESULTS Velocity-compensated gradients resulted in a decreased perfusion fraction (6.9% ± 1.4% vs. 4.4% ± 1.3% in the GM after activation) and pseudo-diffusion coefficient (0.069 ± 0.046 mm2/s vs. 0.014 ± 0.006 in the GM after activation) compared to the bipolar gradients with the longer diffusion encoding time. Increased diffusion coefficients, perfusion fractions, and pseudo-diffusion coefficients were observed in the GM after activation for all gradient profiles. However, the increase was significantly smaller for the velocity-compensated gradients. A diffusion time dependence was found for the pseudo-diffusion coefficient in the activated muscle. CONCLUSION Velocity-compensated diffusion gradients significantly suppress the IVIM effect in the calf muscle, indicating that the ballistic limit is mostly reached, which is supported by the time dependence of the pseudo-diffusion coefficient.
Collapse
Affiliation(s)
- Christoph M Stuprich
- Institute of Radiology, University Hospital Erlangen, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Erlangen, Germany
| | - Martin Loh
- Institute of Radiology, University Hospital Erlangen, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Erlangen, Germany
| | - Johannes T Nemerth
- Institute of Radiology, University Hospital Erlangen, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Erlangen, Germany
| | - Armin M Nagel
- Institute of Radiology, University Hospital Erlangen, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Erlangen, Germany
| | - Michael Uder
- Institute of Radiology, University Hospital Erlangen, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Erlangen, Germany
| | - Frederik B Laun
- Institute of Radiology, University Hospital Erlangen, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Erlangen, Germany
| |
Collapse
|
3
|
Jiang H, Svenningsson L, Topgaard D. Multidimensional encoding of restricted and anisotropic diffusion by double rotation of the q vector. MAGNETIC RESONANCE (GOTTINGEN, GERMANY) 2023; 4:73-85. [PMID: 37904800 PMCID: PMC10583292 DOI: 10.5194/mr-4-73-2023] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/22/2022] [Accepted: 02/11/2023] [Indexed: 11/01/2023]
Abstract
Diffusion NMR and MRI methods building on the classic pulsed gradient spin-echo sequence are sensitive to many aspects of translational motion, including time and frequency dependence ("restriction"), anisotropy, and flow, leading to ambiguities when interpreting experimental data from complex heterogeneous materials such as living biological tissues. While the oscillating gradient technique specifically targets frequency dependence and permits control of the sensitivity to flow, tensor-valued encoding enables investigations of anisotropy in orientationally disordered materials. Here, we propose a simple scheme derived from the "double-rotation" technique in solid-state NMR to generate a family of modulated gradient waveforms allowing for comprehensive exploration of the 2D frequency-anisotropy space and convenient investigation of both restricted and anisotropic diffusion with a single multidimensional acquisition protocol, thereby combining the desirable characteristics of the oscillating gradient and tensor-valued encoding techniques. The method is demonstrated by measuring multicomponent isotropic Gaussian diffusion in simple liquids, anisotropic Gaussian diffusion in a polydomain lyotropic liquid crystal, and restricted diffusion in a yeast cell sediment.
Collapse
Affiliation(s)
- Hong Jiang
- Physical Chemistry, Lund University, P.O. Box 124, 22100 Lund, Sweden
| | - Leo Svenningsson
- Physical Chemistry, Lund University, P.O. Box 124, 22100 Lund, Sweden
| | - Daniel Topgaard
- Physical Chemistry, Lund University, P.O. Box 124, 22100 Lund, Sweden
| |
Collapse
|
4
|
Rosenberg JT, Grant SC, Topgaard D. Nonparametric 5D D-R 2 distribution imaging with single-shot EPI at 21.1 T: Initial results for in vivo rat brain. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2022; 341:107256. [PMID: 35753184 PMCID: PMC9339475 DOI: 10.1016/j.jmr.2022.107256] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/04/2022] [Revised: 05/27/2022] [Accepted: 06/10/2022] [Indexed: 06/15/2023]
Abstract
In vivo human diffusion MRI is by default performed using single-shot EPI with greater than 50-ms echo times and associated signal loss from transverse relaxation. The individual benefits of the current trends of increasing B0 to boost SNR and employing more advanced signal preparation schemes to improve the specificity for selected microstructural properties eventually may be cancelled by increased relaxation rates at high B0 and echo times with advanced encoding. Here, initial attempts to translate state-of-the-art diffusion-relaxation correlation methods from 3 T to 21.1 T are made to identify hurdles that need to be overcome to fulfill the promises of both high SNR and readily interpretable microstructural information.
Collapse
Affiliation(s)
- Jens T Rosenberg
- National High Magnetic Field Laboratory, Florida State University, Tallahassee FL, United States.
| | - Samuel C Grant
- National High Magnetic Field Laboratory, Florida State University, Tallahassee FL, United States; Chemical and Biomedical Engineering, FAMU-FSU College of Engineering, Florida State University, Tallahassee, FL, United States.
| | | |
Collapse
|
5
|
Nonparametric D-R 1-R 2 distribution MRI of the living human brain. Neuroimage 2021; 245:118753. [PMID: 34852278 DOI: 10.1016/j.neuroimage.2021.118753] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2021] [Revised: 11/17/2021] [Accepted: 11/22/2021] [Indexed: 11/23/2022] Open
Abstract
Diffusion-relaxation correlation NMR can simultaneously characterize both the microstructure and the local chemical composition of complex samples that contain multiple populations of water. Recent developments on tensor-valued diffusion encoding and Monte Carlo inversion algorithms have made it possible to transfer diffusion-relaxation correlation NMR from small-bore scanners to clinical MRI systems. Initial studies on clinical MRI systems employed 5D D-R1 and D-R2 correlation to characterize healthy brain in vivo. However, these methods are subject to an inherent bias that originates from not including R2 or R1 in the analysis, respectively. This drawback can be remedied by extending the concept to 6D D-R1-R2 correlation. In this work, we present a sparse acquisition protocol that records all data necessary for in vivo 6D D-R1-R2 correlation MRI across 633 individual measurements within 25 min-a time frame comparable to previous lower-dimensional acquisition protocols. The data were processed with a Monte Carlo inversion algorithm to obtain nonparametric 6D D-R1-R2 distributions. We validated the reproducibility of the method in repeated measurements of healthy volunteers. For a post-therapy glioblastoma case featuring cysts, edema, and partially necrotic remains of tumor, we present representative single-voxel 6D distributions, parameter maps, and artificial contrasts over a wide range of diffusion-, R1-, and R2-weightings based on the rich information contained in the D-R1-R2 distributions.
Collapse
|
6
|
Reymbaut A, Critchley J, Durighel G, Sprenger T, Sughrue M, Bryskhe K, Topgaard D. Toward nonparametric diffusion- T1 characterization of crossing fibers in the human brain. Magn Reson Med 2021; 85:2815-2827. [PMID: 33301195 PMCID: PMC7898694 DOI: 10.1002/mrm.28604] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2020] [Revised: 10/26/2020] [Accepted: 10/27/2020] [Indexed: 12/24/2022]
Abstract
PURPOSE To estimate T 1 for each distinct fiber population within voxels containing multiple brain tissue types. METHODS A diffusion- T 1 correlation experiment was carried out in an in vivo human brain using tensor-valued diffusion encoding and multiple repetition times. The acquired data were inverted using a Monte Carlo algorithm that retrieves nonparametric distributions P ( D , R 1 ) of diffusion tensors and longitudinal relaxation rates R 1 = 1 / T 1 . Orientation distribution functions (ODFs) of the highly anisotropic components of P ( D , R 1 ) were defined to visualize orientation-specific diffusion-relaxation properties. Finally, Monte Carlo density-peak clustering (MC-DPC) was performed to quantify fiber-specific features and investigate microstructural differences between white matter fiber bundles. RESULTS Parameter maps corresponding to P ( D , R 1 ) 's statistical descriptors were obtained, exhibiting the expected R 1 contrast between brain tissue types. Our ODFs recovered local orientations consistent with the known anatomy and indicated differences in R 1 between major crossing fiber bundles. These differences, confirmed by MC-DPC, were in qualitative agreement with previous model-based works but seem biased by the limitations of our current experimental setup. CONCLUSIONS Our Monte Carlo framework enables the nonparametric estimation of fiber-specific diffusion- T 1 features, thereby showing potential for characterizing developmental or pathological changes in T 1 within a given fiber bundle, and for investigating interbundle T 1 differences.
Collapse
Affiliation(s)
- Alexis Reymbaut
- Department of Physical ChemistryLund UniversityLundSweden
- Random Walk Imaging ABLundSweden
| | | | | | - Tim Sprenger
- Karolinska InstituteStockholmSweden
- GE HealthcareStockholmSweden
| | | | | | - Daniel Topgaard
- Department of Physical ChemistryLund UniversityLundSweden
- Random Walk Imaging ABLundSweden
| |
Collapse
|
7
|
Reymbaut A, Mezzani P, de Almeida Martins JP, Topgaard D. Accuracy and precision of statistical descriptors obtained from multidimensional diffusion signal inversion algorithms. NMR IN BIOMEDICINE 2020; 33:e4267. [PMID: 32067322 DOI: 10.1002/nbm.4267] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/18/2019] [Revised: 01/15/2020] [Accepted: 01/16/2020] [Indexed: 05/22/2023]
Abstract
In biological tissues, typical MRI voxels comprise multiple microscopic environments, the local organization of which can be captured by microscopic diffusion tensors. The measured diffusion MRI signal can, therefore, be written as the multidimensional Laplace transform of an intravoxel diffusion tensor distribution (DTD). Tensor-valued diffusion encoding schemes have been designed to probe specific features of the DTD, and several algorithms have been introduced to invert such data and estimate statistical descriptors of the DTD, such as the mean diffusivity, the variance of isotropic diffusivities, and the mean squared diffusion anisotropy. However, the accuracy and precision of these estimations have not been assessed systematically and compared across methods. In this article, we perform and compare such estimations in silico for a one-dimensional Gamma fit, a generalized two-term cumulant approach, and two-dimensional and four-dimensional Monte-Carlo-based inversion techniques, using a clinically feasible tensor-valued acquisition scheme. In particular, we compare their performance at different signal-to-noise ratios (SNRs) for voxel contents varying in terms of the aforementioned statistical descriptors, orientational order, and fractions of isotropic and anisotropic components. We find that all inversion techniques share similar precision (except for a lower precision of the two-dimensional Monte Carlo inversion) but differ in terms of accuracy. While the Gamma fit exhibits infinite-SNR biases when the signal deviates strongly from monoexponentiality and is unaffected by orientational order, the generalized cumulant approach shows infinite-SNR biases when this deviation originates from the variance in isotropic diffusivities or from the low orientational order of anisotropic diffusion components. The two-dimensional Monte Carlo inversion shows remarkable accuracy in all systems studied, given that the acquisition scheme possesses enough directions to yield a rotationally invariant powder average. The four-dimensional Monte Carlo inversion presents no infinite-SNR bias, but suffers significantly from noise in the data, while preserving good contrast in most systems investigated.
Collapse
Affiliation(s)
- Alexis Reymbaut
- Physical Chemistry Department, Lund University, Lund, Sweden
- Random Walk Imaging AB, Lund, Sweden
| | - Paolo Mezzani
- Physical Chemistry Department, Lund University, Lund, Sweden
- Physics Department, Università degli Studi di Milano, Milan, Italy
| | | | - Daniel Topgaard
- Physical Chemistry Department, Lund University, Lund, Sweden
- Random Walk Imaging AB, Lund, Sweden
| |
Collapse
|
8
|
Yon M, de Almeida Martins JP, Bao Q, Budde MD, Frydman L, Topgaard D. Diffusion tensor distribution imaging of an in vivo mouse brain at ultrahigh magnetic field by spatiotemporal encoding. NMR IN BIOMEDICINE 2020; 33:e4355. [PMID: 32812669 PMCID: PMC7583469 DOI: 10.1002/nbm.4355] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/23/2020] [Revised: 05/27/2020] [Accepted: 05/29/2020] [Indexed: 05/08/2023]
Abstract
Diffusion tensor distribution (DTD) imaging builds on principles from diffusion, solid-state and low-field NMR spectroscopies, to quantify the contents of heterogeneous voxels as nonparametric distributions, with tensor "size", "shape" and orientation having direct relations to corresponding microstructural properties of biological tissues. The approach requires the acquisition of multiple images as a function of the magnitude, shape and direction of the diffusion-encoding gradients, leading to long acquisition times unless fast image read-out techniques like EPI are employed. While in previous in vivo human brain studies performed at 3 T this proved a viable option, porting these measurements to very high magnetic fields and/or to heterogeneous organs induces B0 - and B1 -inhomogeneity artifacts that challenge the limits of EPI. To overcome such challenges, we demonstrate here that high spatial resolution DTD of mouse brain can be carried out at 15.2 T with a surface-cryoprobe, by relying on SPatiotemporal ENcoding (SPEN) imaging sequences. These new acquisition and data-processing protocols are demonstrated with measurements on in vivo mouse brain, and validated with synthetic phantoms designed to mimic the diffusion properties of white matter, gray matter and cerebrospinal fluid. While still in need of full extensions to 3D mappings and of scanning additional animals to extract more general physiological conclusions, this work represents another step towards the model-free, noninvasive in vivo characterization of tissue microstructure and heterogeneity in animal models, at ≈0.1 mm resolutions.
Collapse
Affiliation(s)
- Maxime Yon
- Department of Chemical and Biological PhysicsWeizmann InstituteRehovotIsrael
| | - João P. de Almeida Martins
- Division of Physical Chemistry, Department of ChemistryLund UniversityLundSweden
- Random Walk Imaging ABLundSweden
| | - Qingjia Bao
- Department of Chemical and Biological PhysicsWeizmann InstituteRehovotIsrael
| | | | - Lucio Frydman
- Department of Chemical and Biological PhysicsWeizmann InstituteRehovotIsrael
| | - Daniel Topgaard
- Division of Physical Chemistry, Department of ChemistryLund UniversityLundSweden
- Random Walk Imaging ABLundSweden
| |
Collapse
|