1
|
Chakwizira A, Zhu A, Foo T, Westin CF, Szczepankiewicz F, Nilsson M. Diffusion MRI with free gradient waveforms on a high-performance gradient system: Probing restriction and exchange in the human brain. Neuroimage 2023; 283:120409. [PMID: 37839729 DOI: 10.1016/j.neuroimage.2023.120409] [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: 04/05/2023] [Revised: 09/29/2023] [Accepted: 10/12/2023] [Indexed: 10/17/2023] Open
Abstract
The dependence of the diffusion MRI signal on the diffusion time carries signatures of restricted diffusion and exchange. Here we seek to highlight these signatures in the human brain by performing experiments using free gradient waveforms designed to be selectively sensitive to the two effects. We examine six healthy volunteers using both strong and ultra-strong gradients (80, 200 and 300 mT/m). In an experiment featuring a large set of 150 gradient waveforms with different sensitivities to restricted diffusion and exchange, our results reveal unique and different time-dependence signatures in grey and white matter. Grey matter was characterised by both restricted diffusion and exchange and white matter predominantly by restricted diffusion. Exchange in grey matter was at least twice as fast as in white matter, across all subjects and all gradient strengths. The cerebellar cortex featured relatively short exchange times (115 ms). Furthermore, we show that gradient waveforms with tailored designs can be used to map exchange in the human brain. We also assessed the feasibility of clinical applications of the method used in this work and found that the exchange-related contrast obtained with a 25-minute protocol at 300 mT/m was preserved in a 4-minute protocol at 300 mT/m and a 10-minute protocol at 80 mT/m. Our work underlines the utility of free waveforms for detecting time dependence signatures due to restricted diffusion and exchange in vivo, which may potentially serve as a tool for studying diseased tissue.
Collapse
Affiliation(s)
- Arthur Chakwizira
- Department of Medical Radiation Physics, Clinical Sciences Lund, Lund University, Lund, Sweden.
| | - Ante Zhu
- GE Research, Niskayuna, New York, United States
| | - Thomas Foo
- GE Research, Niskayuna, New York, United States
| | - Carl-Fredrik Westin
- Department of Radiology, Brigham and Women's Hospital, Harvard Medical School, Boston, MA, United States
| | - Filip Szczepankiewicz
- Department of Medical Radiation Physics, Clinical Sciences Lund, Lund University, Lund, Sweden
| | - Markus Nilsson
- Department of Clinical Sciences Lund, Radiology, Lund University, Lund, Sweden; Department of Radiology, Skåne University Hospital, Lund, Sweden
| |
Collapse
|
2
|
Lampinen B, Szczepankiewicz F, Lätt J, Knutsson L, Mårtensson J, Björkman-Burtscher IM, van Westen D, Sundgren PC, Ståhlberg F, Nilsson M. Probing brain tissue microstructure with MRI: principles, challenges, and the role of multidimensional diffusion-relaxation encoding. Neuroimage 2023; 282:120338. [PMID: 37598814 DOI: 10.1016/j.neuroimage.2023.120338] [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: 02/02/2023] [Revised: 06/30/2023] [Accepted: 08/17/2023] [Indexed: 08/22/2023] Open
Abstract
Diffusion MRI uses the random displacement of water molecules to sensitize the signal to brain microstructure and to properties such as the density and shape of cells. Microstructure modeling techniques aim to estimate these properties from acquired data by separating the signal between virtual tissue 'compartments' such as the intra-neurite and the extra-cellular space. A key challenge is that the diffusion MRI signal is relatively featureless compared with the complexity of brain tissue. Another challenge is that the tissue microstructure is wildly different within the gray and white matter of the brain. In this review, we use results from multidimensional diffusion encoding techniques to discuss these challenges and their tentative solutions. Multidimensional encoding increases the information content of the data by varying not only the b-value and the encoding direction but also additional experimental parameters such as the shape of the b-tensor and the echo time. Three main insights have emerged from such encoding. First, multidimensional data contradict common model assumptions on diffusion and T2 relaxation, and illustrates how the use of these assumptions cause erroneous interpretations in both healthy brain and pathology. Second, many model assumptions can be dispensed with if data are acquired with multidimensional encoding. The necessary data can be easily acquired in vivo using protocols optimized to minimize Cramér-Rao lower bounds. Third, microscopic diffusion anisotropy reflects the presence of axons but not dendrites. This insight stands in contrast to current 'neurite models' of brain tissue, which assume that axons in white matter and dendrites in gray matter feature highly similar diffusion. Nevertheless, as an axon-based contrast, microscopic anisotropy can differentiate gray and white matter when myelin alterations confound conventional MRI contrasts.
Collapse
Affiliation(s)
- Björn Lampinen
- Clinical Sciences Lund, Diagnostic Radiology, Lund University, Lund, Sweden.
| | | | - Jimmy Lätt
- Department of Medical Imaging and Physiology, Skåne University Hospital Lund, Lund, Sweden
| | - Linda Knutsson
- Clinical Sciences Lund, Medical Radiation Physics, Lund University, Lund, Sweden; Russell H. Morgan Department of Radiology and Radiological Science, Johns Hopkins University School of Medicine, Baltimore, MD, United States; F.M. Kirby Research Center for Functional Brain Imaging, Kennedy Krieger Institute, Baltimore, MD, United States
| | - Johan Mårtensson
- Clinical Sciences Lund, Logopedics, Phoniatrics and Audiology, Lund University, Lund, Sweden
| | - Isabella M Björkman-Burtscher
- Department of Radiology, Institute of Clinical Sciences, Sahlgrenska Academy, University of Gothenburg and Sahlgrenska University Hospital, Gothenburg, Sweden
| | - Danielle van Westen
- Clinical Sciences Lund, Diagnostic Radiology, Lund University, Lund, Sweden; Department of Medical Imaging and Physiology, Skåne University Hospital Lund, Lund, Sweden
| | - Pia C Sundgren
- Clinical Sciences Lund, Diagnostic Radiology, Lund University, Lund, Sweden; Department of Medical Imaging and Physiology, Skåne University Hospital Lund, Lund, Sweden; Lund University BioImaging Centre (LBIC), Lund University, Lund, Sweden
| | - Freddy Ståhlberg
- Clinical Sciences Lund, Diagnostic Radiology, Lund University, Lund, Sweden; Clinical Sciences Lund, Medical Radiation Physics, Lund University, Lund, Sweden
| | - Markus Nilsson
- Clinical Sciences Lund, Diagnostic Radiology, Lund University, Lund, Sweden
| |
Collapse
|
3
|
Chakwizira A, Zhu A, Foo T, Westin CF, Szczepankiewicz F, Nilsson M. Diffusion MRI with free gradient waveforms on a high-performance gradient system: Probing restriction and exchange in the human brain. ARXIV 2023:arXiv:2304.02764v1. [PMID: 37064535 PMCID: PMC10104199] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Subscribe] [Scholar Register] [Indexed: 04/18/2023]
Abstract
The dependence of the diffusion MRI signal on the diffusion time carries signatures of restricted diffusion and exchange. Here we seek to highlight these signatures in the human brain by performing experiments using free gradient waveforms that are selectively sensitive to the two effects. We examine six healthy volunteers using both strong and ultra-strong gradients (80, 200 and 300 mT/m). In an experiment featuring a large set of gradient waveforms with different sensitivities to restricted diffusion and exchange (150 samples), our results reveal unique time-dependence signatures in grey and white matter, where the former is characterised by both restricted diffusion and exchange and the latter predominantly exhibits restricted diffusion. Furthermore, we show that gradient waveforms with independently varying sensitivities to restricted diffusion and exchange can be used to map exchange in the human brain. We consistently find that exchange in grey matter is at least twice as fast as in white matter, across all subjects and all gradient strengths. The shortest exchange times observed in this study were in the cerebellar cortex (115 ms). We also assess the feasibility of future clinical applications of the method used in this work, where we find that the grey-white matter exchange contrast obtained with a 25-minute 300 mT/m protocol is preserved by a 4-minute 300 mT/m and a 10-minute 80 mT/m protocol. Our work underlines the utility of free waveforms for detecting time-dependence signatures due to restricted diffusion and exchange in vivo, which may potentially serve as a tool for studying diseased tissue.
Collapse
Affiliation(s)
- Arthur Chakwizira
- Medical Radiation Physics, Clinical Sciences Lund, Lund University, Lund, Sweden
| | - Ante Zhu
- GE Research, Niskayuna, New York, USA
| | | | - Carl-Fredrik Westin
- Department of Radiology, Brigham and Women’s Hospital, Harvard Medical School, Boston, MA, United States
| | | | - Markus Nilsson
- Department of Clinical Sciences Lund, Radiology, Lund University, Lund, Sweden
| |
Collapse
|
4
|
Li C, Fieremans E, Novikov DS, Ge Y, Zhang J. Measuring water exchange on a preclinical MRI system using filter exchange and diffusion time dependent kurtosis imaging. Magn Reson Med 2023; 89:1441-1455. [PMID: 36404493 PMCID: PMC9892228 DOI: 10.1002/mrm.29536] [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: 12/19/2021] [Revised: 11/01/2022] [Accepted: 11/02/2022] [Indexed: 11/22/2022]
Abstract
PURPOSE Filter exchange imaging (FEXI) and diffusion time (t)-dependent diffusion kurtosis imaging (DKI(t)) are both sensitive to water exchange between tissue compartments. The restrictive effects of tissue microstructure, however, introduce bias to the exchange rate obtained by these two methods, as their interpretation conventionally rely on the Kärger model of barrier limited exchange between Gaussian compartments. Here, we investigated whether FEXI and DKI(t) can provide comparable exchange rates in ex vivo mouse brains. THEORY AND METHODS FEXI and DKI(t) data were acquired from ex vivo mouse brains on a preclinical MRI system. Phase cycling and negative slice prewinder gradients were used to minimize the interferences from imaging gradients. RESULTS In the corpus callosum, apparent exchange rate (AXR) from FEXI correlated with the exchange rate (the inverse of exchange time, 1/τex ) from DKI(t) along the radial direction. In comparison, discrepancies between FEXI and DKI(t) were found in the cortex due to low filter efficiency and confounding effects from tissue microstructure. CONCLUSION The results suggest that FEXI and DKI(t) are sensitive to the same exchange processes in white matter when separated from restrictive effects of microstructure. The complex microstructure in gray matter, with potential exchange among multiple compartments and confounding effects of microstructure, still pose a challenge for FEXI and DKI(t).
Collapse
Affiliation(s)
- Chenyang Li
- Department of Radiology, Center for Biomedical Imaging, NYU Grossman School of Medicine, New York, NY, USA
- Vilcek Institute of Graduate Biomedical Sciences, NYU Grossman School of Medicine, New York, NY, USA
| | - Els Fieremans
- Department of Radiology, Center for Biomedical Imaging, NYU Grossman School of Medicine, New York, NY, USA
| | - Dmitry S. Novikov
- Department of Radiology, Center for Biomedical Imaging, NYU Grossman School of Medicine, New York, NY, USA
| | - Yulin Ge
- Department of Radiology, Center for Biomedical Imaging, NYU Grossman School of Medicine, New York, NY, USA
| | - Jiangyang Zhang
- Department of Radiology, Center for Biomedical Imaging, NYU Grossman School of Medicine, New York, NY, USA
| |
Collapse
|
5
|
Chakwizira A, Westin C, Brabec J, Lasič S, Knutsson L, Szczepankiewicz F, Nilsson M. Diffusion MRI with pulsed and free gradient waveforms: Effects of restricted diffusion and exchange. NMR IN BIOMEDICINE 2023; 36:e4827. [PMID: 36075110 PMCID: PMC10078514 DOI: 10.1002/nbm.4827] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/16/2021] [Revised: 08/27/2022] [Accepted: 09/06/2022] [Indexed: 05/06/2023]
Abstract
Monitoring time dependence with diffusion MRI yields observables sensitive to compartment sizes (restricted diffusion) and membrane permeability (water exchange). However, restricted diffusion and exchange have opposite effects on the diffusion-weighted signal, which can lead to errors in parameter estimates. In this work, we propose a signal representation that incorporates the effects of both restricted diffusion and exchange up to second order in b-value and is compatible with gradient waveforms of arbitrary shape. The representation features mappings from a gradient waveform to two scalars that separately control the sensitivity to restriction and exchange. We demonstrate that these scalars span a two-dimensional space that can be used to choose waveforms that selectively probe restricted diffusion or exchange, eliminating the correlation between the two phenomena. We found that waveforms with specific but unconventional shapes provide an advantage over conventional pulsed and oscillating gradient acquisitions. We also show that parametrization of waveforms into a two-dimensional space can be used to understand protocols from other approaches that probe restricted diffusion and exchange. For example, we found that the variation of mixing time in filter-exchange imaging corresponds to variation of our exchange-weighting scalar at a fixed value of the restriction-weighting scalar. The proposed signal representation was evaluated using Monte Carlo simulations in identical parallel cylinders with hexagonal and random packing as well as parallel cylinders with gamma-distributed radii. Results showed that the approach is sensitive to sizes in the interval 4-12 μm and exchange rates in the simulated range of 0 to 20 s - 1 , but also that there is a sensitivity to the extracellular geometry. The presented theory constitutes a simple and intuitive description of how restricted diffusion and exchange influence the signal as well as a guide to protocol design capable of separating the two effects.
Collapse
Affiliation(s)
- Arthur Chakwizira
- Department of Medical Radiation Physics, LundLund UniversityLundSweden
| | - Carl‐Fredrik Westin
- Department of Radiology, Brigham and Women's HospitalHarvard Medical SchoolBostonMassachusettsUSA
| | - Jan Brabec
- Department of Medical Radiation Physics, LundLund UniversityLundSweden
| | - Samo Lasič
- Danish Research Centre for Magnetic Resonance, Centre for Functional and Diagnostic Imaging and ResearchCopenhagen University Hospital ‐ Amager and HvidovreCopenhagenDenmark
- Random Walk Imaging ABLundSweden
| | - Linda Knutsson
- Department of Medical Radiation Physics, LundLund UniversityLundSweden
- Russell H. Morgan Department of Radiology and Radiological ScienceJohns Hopkins University School of MedicineBaltimoreMarylandUSA
- F. M. Kirby Research Center for Functional Brain ImagingKennedy Krieger InstituteBaltimoreMarylandUSA
| | | | - Markus Nilsson
- Department of Clinical Sciences Lund, RadiologyLund UniversityLundSweden
| |
Collapse
|
6
|
Afzali M, Pieciak T, Jones DK, Schneider JE, Özarslan E. Cumulant expansion with localization: A new representation of the diffusion MRI signal. FRONTIERS IN NEUROIMAGING 2022; 1:958680. [PMID: 37555138 PMCID: PMC10406302 DOI: 10.3389/fnimg.2022.958680] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 05/31/2022] [Accepted: 07/19/2022] [Indexed: 08/10/2023]
Abstract
Diffusion MR is sensitive to the microstructural features of a sample. Fine-scale characteristics can be probed by employing strong diffusion gradients while the low b-value regime is determined by the cumulants of the distribution of particle displacements. A signal representation based on the cumulants, however, suffers from a finite convergence radius and cannot represent the 'localization regime' characterized by a stretched exponential decay that emerges at large gradient strengths. Here, we propose a new representation for the diffusion MR signal. Our method provides not only a robust estimate of the first three cumulants but also a meaningful extrapolation of the entire signal decay.
Collapse
Affiliation(s)
- Maryam Afzali
- Leeds Institute of Cardiovascular and Metabolic Medicine, University of Leeds, Leeds, United Kingdom
- Cardiff University Brain Research Imaging Centre (CUBRIC), School of Psychology, Cardiff University, Cardiff, United Kingdom
| | - Tomasz Pieciak
- LPI, ETSI Telecomunicación, Universidad de Valladolid, Valladolid, Spain
| | - Derek K. Jones
- Cardiff University Brain Research Imaging Centre (CUBRIC), School of Psychology, Cardiff University, Cardiff, United Kingdom
| | - Jürgen E. Schneider
- Leeds Institute of Cardiovascular and Metabolic Medicine, University of Leeds, Leeds, United Kingdom
| | - Evren Özarslan
- Department of Biomedical Engineering, Linköping University, Linköping, Sweden
- Center for Medical Image Science and Visualization, Linköping University, Linköping, Sweden
| |
Collapse
|
7
|
Novello L, Henriques RN, Ianuş A, Feiweier T, Shemesh N, Jovicich J. In vivo Correlation Tensor MRI reveals microscopic kurtosis in the human brain on a clinical 3T scanner. Neuroimage 2022; 254:119137. [PMID: 35339682 DOI: 10.1016/j.neuroimage.2022.119137] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2021] [Revised: 02/17/2022] [Accepted: 03/22/2022] [Indexed: 12/15/2022] Open
Abstract
Diffusion MRI (dMRI) has become one of the most important imaging modalities for noninvasively probing tissue microstructure. Diffusional Kurtosis MRI (DKI) quantifies the degree of non-gaussian diffusion, which in turn has been shown to increase sensitivity towards, e.g., disease and orientation mapping in neural tissue. However, the specificity of DKI is limited as different sources can contribute to the total intravoxel diffusional kurtosis, including: variance in diffusion tensor magnitudes (Kiso), variance due to diffusion anisotropy (Kaniso), and microscopic kurtosis (μK) related to restricted diffusion, microstructural disorder, and/or exchange. Interestingly, μK is typically ignored in diffusion MRI signal modeling as it is assumed to be negligible in neural tissues. However, recently, Correlation Tensor MRI (CTI) based on Double-Diffusion-Encoding (DDE) was introduced for kurtosis source separation, revealing non negligible μK in preclinical imaging. Here, we implemented CTI for the first time on a clinical 3T scanner and investigated the sources of total kurtosis in healthy subjects. A robust framework for kurtosis source separation in humans is introduced, followed by estimation of μK (and the other kurtosis sources) in the healthy brain. Using this clinical CTI approach, we find that μK significantly contributes to total diffusional kurtosis both in gray and white matter tissue but, as expected, not in the ventricles. The first μK maps of the human brain are presented, revealing that the spatial distribution of μK provides a unique source of contrast, appearing different from isotropic and anisotropic kurtosis counterparts. Moreover, group average templates of these kurtosis sources have been generated for the first time, which corroborated our findings at the underlying individual-level maps. We further show that the common practice of ignoring μK and assuming the multiple gaussian component approximation for kurtosis source estimation introduces significant bias in the estimation of other kurtosis sources and, perhaps even worse, compromises their interpretation. Finally, a twofold acceleration of CTI is discussed in the context of potential future clinical applications. We conclude that CTI has much potential for future in vivo microstructural characterizations in healthy and pathological tissue.
Collapse
Affiliation(s)
- Lisa Novello
- Center for Mind/Brain Sciences - CIMeC, University of Trento, Rovereto, Italy.
| | | | - Andrada Ianuş
- Champalimaud Research, Champalimaud Foundation, Lisbon, Portugal
| | | | - Noam Shemesh
- Champalimaud Research, Champalimaud Foundation, Lisbon, Portugal
| | - Jorge Jovicich
- Center for Mind/Brain Sciences - CIMeC, University of Trento, Rovereto, Italy
| |
Collapse
|
8
|
Olesen JL, Østergaard L, Shemesh N, Jespersen SN. Diffusion time dependence, power-law scaling, and exchange in gray matter. Neuroimage 2022; 251:118976. [PMID: 35168088 PMCID: PMC8961002 DOI: 10.1016/j.neuroimage.2022.118976] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/25/2021] [Revised: 12/24/2021] [Accepted: 02/04/2022] [Indexed: 12/27/2022] Open
Abstract
Characterizing neural tissue microstructure is a critical goal for future neuroimaging. Diffusion MRI (dMRI) provides contrasts that reflect diffusing spins’ interactions with myriad microstructural features of biological systems. However, the specificity of dMRI remains limited due to the ambiguity of its signals vis-à-vis the underlying microstructure. To improve specificity, biophysical models of white matter (WM) typically express dMRI signals according to the Standard Model (SM) and have more recently in gray matter (GM) taken spherical compartments into account (the SANDI model) in attempts to represent cell soma. The validity of the assumptions underlying these models, however, remains largely undetermined, especially in GM. To validate these assumptions experimentally, observing their unique, functional properties, such as the b−1/2 power-law associated with one-dimensional diffusion, has emerged as a fruitful strategy. The absence of this signature in GM, in turn, has been explained by neurite water exchange, non-linear morphology, and/or by obscuring soma signal contributions. Here, we present diffusion simulations in realistic neurons demonstrating that curvature and branching does not destroy the stick power-law behavior in impermeable neurites, but also that their signal is drowned by the soma signal under typical experimental conditions. Nevertheless, by studying the GM dMRI signal's behavior as a function of diffusion weighting as well as time, we identify an attainable experimental regime in which the neurite signal dominates. Furthermore, we find that exchange-driven time dependence produces a signal behavior opposite to that which would be expected from restricted diffusion, thereby providing a functional signature that disambiguates the two effects. We present data from dMRI experiments in ex vivo rat brain at ultrahigh field of 16.4T and observe a time dependence that is consistent with substantial exchange but also with a GM stick power-law. The first finding suggests significant water exchange between neurites and the extracellular space while the second suggests a small sub-population of impermeable neurites. To quantify these observations, we harness the Kärger exchange model and incorporate the corresponding signal time dependence in the SM and SANDI models.
Collapse
|
9
|
Henriques RN, Jespersen SN, Shemesh N. Evidence for microscopic kurtosis in neural tissue revealed by correlation tensor MRI. Magn Reson Med 2021; 86:3111-3130. [PMID: 34329509 PMCID: PMC9290035 DOI: 10.1002/mrm.28938] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2021] [Revised: 07/01/2021] [Accepted: 07/04/2021] [Indexed: 12/14/2022]
Abstract
PURPOSE The impact of microscopic diffusional kurtosis (µK), arising from restricted diffusion and/or structural disorder, remains a controversial issue in contemporary diffusion MRI (dMRI). Recently, correlation tensor imaging (CTI) was introduced to disentangle the sources contributing to diffusional kurtosis, without relying on a-priori multi-gaussian component (MGC) or other microstructural assumptions. Here, we investigated µK in in vivo rat brains and assessed its impact on state-of-the-art methods ignoring µK. THEORY AND METHODS CTI harnesses double diffusion encoding (DDE) experiments, which were here improved for speed and minimal bias using four different sets of acquisition parameters. The robustness of the improved CTI protocol was assessed via simulations. In vivo CTI acquisitions were performed in healthy rat brains using a 9.4T pre-clinical scanner equipped with a cryogenic coil, and targeted the estimation of µK, anisotropic kurtosis, and isotropic kurtosis. RESULTS The improved CTI acquisition scheme substantially reduces scan time and importantly, also minimizes higher-order-term biases, thus enabling robust µK estimation, alongside Kaniso and Kiso metrics. Our CTI experiments revealed positive µK both in white and gray matter of the rat brain in vivo; µK is the dominant kurtosis source in healthy gray matter tissue. The non-negligible µK substantially were found to bias prior MGC analyses of Kiso and Kaniso . CONCLUSIONS Correlation Tensor MRI offers a more accurate and robust characterization of kurtosis sources than its predecessors. µK is non-negligible in vivo in healthy white and gray matter tissues and could be an important biomarker for future studies. Our findings thus have both theoretical and practical implications for future dMRI research.
Collapse
Affiliation(s)
| | - Sune N Jespersen
- Center of Functionally Integrative Neuroscience (CFIN) and MINDLab, Clinical Institute, Aarhus University, Aarhus, Denmark.,Department of Physics and Astronomy, Aarhus University, Aarhus, Denmark
| | - Noam Shemesh
- Champalimaud Research, Champalimaud Centre for the Unknown, Lisbon, Portugal
| |
Collapse
|
10
|
Lampinen B, Lätt J, Wasselius J, van Westen D, Nilsson M. Time dependence in diffusion MRI predicts tissue outcome in ischemic stroke patients. Magn Reson Med 2021; 86:754-764. [PMID: 33755261 PMCID: PMC8445077 DOI: 10.1002/mrm.28743] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2020] [Revised: 01/29/2021] [Accepted: 01/30/2021] [Indexed: 12/18/2022]
Abstract
Purpose: Reperfusion therapy enables effective treatment of ischemic stroke presenting within 4–6 hours. However, tissue progression from ischemia to infarction is variable, and some patients benefit from treatment up until 24 hours. Improved imaging techniques are needed to identify these patients. Here, it was hypothesized that time dependence in diffusion MRI may predict tissue outcome in ischemic stroke. Methods: Diffusion MRI data were acquired with multiple diffusion times in five non-reperfused patients at 2, 9, and 100 days after stroke onset. Maps of “rate of kurtosis change” (k), mean kurtosis, ADC, and fractional anisotropy were derived. The ADC maps defined lesions, normal-appearing tissue, and the lesion tissue that would either be infarcted or remain viable by day 100. Diffusion parameters were compared (1) between lesions and normal-appearing tissue, and (2) between lesion tissue that would be infarcted or remain viable. Results: Positive values of k were observed within stroke lesions on day 2 (P = .001) and on day 9 (P = .023), indicating diffusional exchange. On day 100, high ADC values indicated infarction of 50 ± 20% of the lesion volumes. Tissue infarction was predicted by high k values both on day 2 (P = .026) and on day 9 (P = .046), by low mean kurtosis values on day 2 (P = .043), and by low fractional anisotropy values on day 9 (P = .029), but not by low ADC values. Conclusions: Diffusion time dependence predicted tissue outcome in ischemic stroke more accurately than the ADC, and may be useful for predicting reperfusion benefit.
Collapse
Affiliation(s)
- Björn Lampinen
- Clinical Sciences Lund, Medical Radiation Physics, Lund University, Lund, Sweden
| | - Jimmy Lätt
- Center for Medical Imaging and Physiology, Skåne University Hospital Lund, Lund, Sweden
| | - Johan Wasselius
- Clinical Sciences Lund, Radiology, Lund University, Lund, Sweden
| | | | - Markus Nilsson
- Clinical Sciences Lund, Radiology, Lund University, Lund, Sweden
| |
Collapse
|
11
|
Szczepankiewicz F, Westin CF, Nilsson M. Gradient waveform design for tensor-valued encoding in diffusion MRI. J Neurosci Methods 2021; 348:109007. [PMID: 33242529 PMCID: PMC8443151 DOI: 10.1016/j.jneumeth.2020.109007] [Citation(s) in RCA: 34] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2020] [Revised: 11/17/2020] [Accepted: 11/19/2020] [Indexed: 12/13/2022]
Abstract
Diffusion encoding along multiple spatial directions per signal acquisition can be described in terms of a b-tensor. The benefit of tensor-valued diffusion encoding is that it unlocks the 'shape of the b-tensor' as a new encoding dimension. By modulating the b-tensor shape, we can control the sensitivity to microscopic diffusion anisotropy which can be used as a contrast mechanism; a feature that is inaccessible by conventional diffusion encoding. Since imaging methods based on tensor-valued diffusion encoding are finding an increasing number of applications we are prompted to highlight the challenge of designing the optimal gradient waveforms for any given application. In this review, we first establish the basic design objectives in creating field gradient waveforms for tensor-valued diffusion MRI. We also survey additional design considerations related to limitations imposed by hardware and physiology, potential confounding effects that cannot be captured by the b-tensor, and artifacts related to the diffusion encoding waveform. Throughout, we discuss the expected compromises and tradeoffs with an aim to establish a more complete understanding of gradient waveform design and its impact on accurate measurements and interpretations of data.
Collapse
Affiliation(s)
- Filip Szczepankiewicz
- Radiology, Brigham and Women's Hospital, Boston, MA, United States; Harvard Medical School, Boston, MA, United States; Clinical Sciences, Lund University, Lund, Sweden.
| | - Carl-Fredrik Westin
- Radiology, Brigham and Women's Hospital, Boston, MA, United States; Harvard Medical School, Boston, MA, United States
| | | |
Collapse
|
12
|
Afzali M, Pieciak T, Newman S, Garyfallidis E, Özarslan E, Cheng H, Jones DK. The sensitivity of diffusion MRI to microstructural properties and experimental factors. J Neurosci Methods 2021; 347:108951. [PMID: 33017644 PMCID: PMC7762827 DOI: 10.1016/j.jneumeth.2020.108951] [Citation(s) in RCA: 38] [Impact Index Per Article: 12.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2020] [Revised: 08/27/2020] [Accepted: 09/15/2020] [Indexed: 12/13/2022]
Abstract
Diffusion MRI is a non-invasive technique to study brain microstructure. Differences in the microstructural properties of tissue, including size and anisotropy, can be represented in the signal if the appropriate method of acquisition is used. However, to depict the underlying properties, special care must be taken when designing the acquisition protocol as any changes in the procedure might impact on quantitative measurements. This work reviews state-of-the-art methods for studying brain microstructure using diffusion MRI and their sensitivity to microstructural differences and various experimental factors. Microstructural properties of the tissue at a micrometer scale can be linked to the diffusion signal at a millimeter-scale using modeling. In this paper, we first give an introduction to diffusion MRI and different encoding schemes. Then, signal representation-based methods and multi-compartment models are explained briefly. The sensitivity of the diffusion MRI signal to the microstructural components and the effects of curvedness of axonal trajectories on the diffusion signal are reviewed. Factors that impact on the quality (accuracy and precision) of derived metrics are then reviewed, including the impact of random noise, and variations in the acquisition parameters (i.e., number of sampled signals, b-value and number of acquisition shells). Finally, yet importantly, typical approaches to deal with experimental factors are depicted, including unbiased measures and harmonization. We conclude the review with some future directions and recommendations on this topic.
Collapse
Affiliation(s)
- Maryam Afzali
- Cardiff University Brain Research Imaging Centre (CUBRIC), School of Psychology, Cardiff University, Cardiff, United Kingdom.
| | - Tomasz Pieciak
- AGH University of Science and Technology, Kraków, Poland; LPI, ETSI Telecomunicación, Universidad de Valladolid, Valladolid, Spain.
| | - Sharlene Newman
- Department of Psychological and Brain Sciences, Indiana University, Bloomington, IN 47405, USA; Program of Neuroscience, Indiana University, Bloomington, IN 47405, USA.
| | - Eleftherios Garyfallidis
- Program of Neuroscience, Indiana University, Bloomington, IN 47405, USA; Department of Intelligent Systems Engineering, Indiana University, Bloomington, IN 47408, USA.
| | - Evren Özarslan
- Department of Biomedical Engineering, Linköping University, Linköping, Sweden; Center for Medical Image Science and Visualization, Linköping University, Linköping, Sweden.
| | - Hu Cheng
- Department of Psychological and Brain Sciences, Indiana University, Bloomington, IN 47405, USA; Program of Neuroscience, Indiana University, Bloomington, IN 47405, USA.
| | - Derek K Jones
- Cardiff University Brain Research Imaging Centre (CUBRIC), School of Psychology, Cardiff University, Cardiff, United Kingdom.
| |
Collapse
|
13
|
Lee HH, Papaioannou A, Novikov DS, Fieremans E. In vivo observation and biophysical interpretation of time-dependent diffusion in human cortical gray matter. Neuroimage 2020; 222:117054. [PMID: 32585341 PMCID: PMC7736473 DOI: 10.1016/j.neuroimage.2020.117054] [Citation(s) in RCA: 38] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2019] [Revised: 06/08/2020] [Accepted: 06/10/2020] [Indexed: 12/25/2022] Open
Abstract
The dependence of the diffusion MRI signal on the diffusion time t is a hallmark of tissue microstructure at the scale of the diffusion length. Here we measure the time-dependence of the mean diffusivity D(t) and mean kurtosis K(t) in cortical gray matter and in 25 gray matter sub-regions, in 10 healthy subjects. Significant diffusivity and kurtosis time-dependence is observed for t=21.2-100 ms, and is characterized by a power-law tail ∼t-ϑ with dynamical exponent ϑ. To interpret our measurements, we systematize the relevant scenarios and mechanisms for diffusion time-dependence in the brain. Using the effective medium theory formalism, we derive an exact relation between the power-law tails in D(t) and K(t). The estimated dynamical exponent ϑ≃1/2 in both D(t) and K(t) is consistent with one-dimensional diffusion in the presence of randomly positioned restrictions along neurites. We analyze the short-range disordered statistics of synapses on axon collaterals in the cortex, and perform one-dimensional Monte Carlo simulations of diffusion restricted by permeable barriers with a similar randomness in their placement, to confirm the ϑ=1/2 exponent. In contrast, the Kärger model of exchange is less consistent with the data since it does not capture the diffusivity time-dependence, and the estimated exchange time from K(t) falls below our measured t-range. Although we cannot exclude exchange as a contributing factor, we argue that structural disorder along neurites is mainly responsible for the observed time-dependence of diffusivity and kurtosis. Our observation and theoretical interpretation of the t-1/2 tail in D(t) and K(t) altogether establish the sensitivity of a macroscopic MRI signal to micrometer-scale structural heterogeneities along neurites in human gray matter in vivo.
Collapse
Affiliation(s)
- Hong-Hsi Lee
- Center for Biomedical Imaging, Department of Radiology, New York University School of Medicine, New York, NY, USA; Center for Advanced Imaging Innovation and Research (CAI2R), New York University School of Medicine, New York, NY, USA.
| | - Antonios Papaioannou
- Center for Biomedical Imaging, Department of Radiology, New York University School of Medicine, New York, NY, USA; Center for Advanced Imaging Innovation and Research (CAI2R), New York University School of Medicine, New York, NY, USA
| | - Dmitry S Novikov
- Center for Biomedical Imaging, Department of Radiology, New York University School of Medicine, New York, NY, USA; Center for Advanced Imaging Innovation and Research (CAI2R), New York University School of Medicine, New York, NY, USA
| | - Els Fieremans
- Center for Biomedical Imaging, Department of Radiology, New York University School of Medicine, New York, NY, USA; Center for Advanced Imaging Innovation and Research (CAI2R), New York University School of Medicine, New York, NY, USA
| |
Collapse
|
14
|
Henriques RN, Palombo M, Jespersen SN, Shemesh N, Lundell H, Ianuş A. Double diffusion encoding and applications for biomedical imaging. J Neurosci Methods 2020; 348:108989. [PMID: 33144100 DOI: 10.1016/j.jneumeth.2020.108989] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/09/2020] [Revised: 09/25/2020] [Accepted: 10/20/2020] [Indexed: 12/11/2022]
Abstract
Diffusion Magnetic Resonance Imaging (dMRI) is one of the most important contemporary non-invasive modalities for probing tissue structure at the microscopic scale. The majority of dMRI techniques employ standard single diffusion encoding (SDE) measurements, covering different sequence parameter ranges depending on the complexity of the method. Although many signal representations and biophysical models have been proposed for SDE data, they are intrinsically limited by a lack of specificity. Advanced dMRI methods have been proposed to provide additional microstructural information beyond what can be inferred from SDE. These enhanced contrasts can play important roles in characterizing biological tissues, for instance upon diseases (e.g. neurodegenerative, cancer, stroke), aging, learning, and development. In this review we focus on double diffusion encoding (DDE), which stands out among other advanced acquisitions for its versatility, ability to probe more specific diffusion correlations, and feasibility for preclinical and clinical applications. Various DDE methodologies have been employed to probe compartment sizes (Section 3), decouple the effects of microscopic diffusion anisotropy from orientation dispersion (Section 4), probe displacement correlations, study exchange, or suppress fast diffusing compartments (Section 6). DDE measurements can also be used to improve the robustness of biophysical models (Section 5) and study intra-cellular diffusion via magnetic resonance spectroscopy of metabolites (Section 7). This review discusses all these topics as well as important practical aspects related to the implementation and contrast in preclinical and clinical settings (Section 9) and aims to provide the readers a guide for deciding on the right DDE acquisition for their specific application.
Collapse
Affiliation(s)
- Rafael N Henriques
- Champalimaud Research, Champalimaud Centre for the Unknown, Lisbon, Portugal
| | - Marco Palombo
- Centre for Medical Image Computing and Dept. of Computer Science, University College London, London, UK
| | - Sune N Jespersen
- Center of Functionally Integrative Neuroscience (CFIN) and MINDLab, Department of Clinical Medicine, Aarhus University, Aarhus, Denmark; Department of Physics and Astronomy, Aarhus University, Aarhus, Denmark
| | - Noam Shemesh
- Champalimaud Research, Champalimaud Centre for the Unknown, Lisbon, Portugal
| | - Henrik Lundell
- Danish Research Centre for Magnetic Resonance, Centre for Functional and Diagnostic Imaging and Research, Copenhagen University Hospital Hvidovre, Denmark
| | - Andrada Ianuş
- Champalimaud Research, Champalimaud Centre for the Unknown, Lisbon, Portugal.
| |
Collapse
|
15
|
Szczepankiewicz F, Sjölund J, Dall'Armellina E, Plein S, Schneider JE, Teh I, Westin CF. Motion-compensated gradient waveforms for tensor-valued diffusion encoding by constrained numerical optimization. Magn Reson Med 2020; 85:2117-2126. [PMID: 33048401 PMCID: PMC7821235 DOI: 10.1002/mrm.28551] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2020] [Revised: 09/18/2020] [Accepted: 09/21/2020] [Indexed: 12/11/2022]
Abstract
PURPOSE Diffusion-weighted MRI is sensitive to incoherent tissue motion, which may confound the measured signal and subsequent analysis. We propose a "motion-compensated" gradient waveform design for tensor-valued diffusion encoding that negates the effects bulk motion and incoherent motion in the ballistic regime. METHODS Motion compensation was achieved by constraining the magnitude of gradient waveform moment vectors. The constraint was incorporated into a numerical optimization framework, along with existing constraints that account for b-tensor shape, hardware restrictions, and concomitant field gradients. We evaluated the efficacy of encoding and motion compensation in simulations, and we demonstrated the approach by linear and planar b-tensor encoding in a healthy heart in vivo. RESULTS The optimization framework produced asymmetric motion-compensated waveforms that yielded b-tensors of arbitrary shape with improved efficiency compared with previous designs for tensor-valued encoding, and equivalent efficiency to previous designs for linear (conventional) encoding. Technical feasibility was demonstrated in the heart in vivo, showing vastly improved data quality when using motion compensation. The optimization framework is available online in open source. CONCLUSION Our gradient waveform design is both more flexible and efficient than previous methods, facilitating tensor-valued diffusion encoding in tissues in which motion would otherwise confound the signal. The proposed design exploits asymmetric encoding times, a single refocusing pulse or multiple refocusing pulses, and integrates compensation for concomitant gradient effects throughout the imaging volume.
Collapse
Affiliation(s)
- Filip Szczepankiewicz
- Harvard Medical School, Boston, Massachusetts, USA.,Radiology, Brigham and Women's Hospital, Boston, Massachusetts, USA.,Diagnostic Radiology, Clinical Sciences Lund, Lund University, Lund, Sweden
| | - Jens Sjölund
- Elekta Instrument AB, Stockholm, Sweden.,Department of Information Technology, Uppsala University, Uppsala, Sweden
| | - Erica Dall'Armellina
- Leeds Institute of Cardiovascular and Metabolic Medicine, University of Leeds, Leeds, United Kingdom
| | - Sven Plein
- Leeds Institute of Cardiovascular and Metabolic Medicine, University of Leeds, Leeds, United Kingdom
| | - Jürgen E Schneider
- Leeds Institute of Cardiovascular and Metabolic Medicine, University of Leeds, Leeds, United Kingdom
| | - Irvin Teh
- Leeds Institute of Cardiovascular and Metabolic Medicine, University of Leeds, Leeds, United Kingdom
| | - Carl-Fredrik Westin
- Harvard Medical School, Boston, Massachusetts, USA.,Radiology, Brigham and Women's Hospital, Boston, Massachusetts, USA
| |
Collapse
|
16
|
Zhou Z, Tong Q, Zhang L, Ding Q, Lu H, Jonkman LE, Yao J, He H, Zhu K, Zhong J. Evaluation of the diffusion MRI white matter tract integrity model using myelin histology and Monte-Carlo simulations. Neuroimage 2020; 223:117313. [PMID: 32882384 DOI: 10.1016/j.neuroimage.2020.117313] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2019] [Revised: 08/21/2020] [Accepted: 08/24/2020] [Indexed: 12/13/2022] Open
Abstract
Quantitative evaluation of brain myelination has drawn considerable attention. Conventional diffusion-based magnetic resonance imaging models, including diffusion tensor imaging and diffusion kurtosis imaging (DKI),1 have been used to infer the microstructure and its changes in neurological diseases. White matter tract integrity (WMTI) was proposed as a biophysical model to relate the DKI-derived metrics to the underlying microstructure. Although the model has been validated on ex vivo animal brains, it was not well evaluated with ex vivo human brains. In this study, histological samples (namely corpus callosum) from postmortem human brains have been investigated based on WMTI analyses on a clinical 3T scanner and comparisons with gold standard myelin staining in proteolipid protein and Luxol fast blue. In addition, Monte Carlo simulations were conducted to link changes from ex vivo to in vivo conditions based on the microscale parameters of water diffusivity and permeability. The results show that WMTI metrics, including axonal water fraction AWF, radial extra-axonal diffusivity De⊥, and intra-axonal diffusivity Dawere needed to characterize myelin content alterations. Thus, WMTI model metrics are shown to be promising candidates as sensitive biomarkers of demyelination.
Collapse
Affiliation(s)
- Zihan Zhou
- Center for Brain Imaging Science and Technology, Key Laboratory for Biomedical Engineering of Ministry of Education, College of Biomedical Engineering and Instrumental Science, Zhejiang University, Zhouyiqing Building, Room 314, Yuquan Campus, Hangzhou 310027, China
| | - Qiqi Tong
- Center for Brain Imaging Science and Technology, Key Laboratory for Biomedical Engineering of Ministry of Education, College of Biomedical Engineering and Instrumental Science, Zhejiang University, Zhouyiqing Building, Room 314, Yuquan Campus, Hangzhou 310027, China
| | - Lei Zhang
- China Brain Bank and Department of Neurology in Second Affiliated Hospital, Key Laboratory of Medical Neurobiology of Zhejiang Province, and Department of Neurobiology, Zhejiang University School of Medicine, Hangzhou 310058, China; Department of Pathology, Zhejiang University School of Medicine, Hangzhou 310058, China
| | - Qiuping Ding
- Center for Brain Imaging Science and Technology, Key Laboratory for Biomedical Engineering of Ministry of Education, College of Biomedical Engineering and Instrumental Science, Zhejiang University, Zhouyiqing Building, Room 314, Yuquan Campus, Hangzhou 310027, China
| | - Hui Lu
- China Brain Bank and Department of Neurology in Second Affiliated Hospital, Key Laboratory of Medical Neurobiology of Zhejiang Province, and Department of Neurobiology, Zhejiang University School of Medicine, Hangzhou 310058, China
| | - Laura E Jonkman
- Department of Anatomy and Neurosciences, Amsterdam Neuroscience, Amsterdam UMC, location VUmc, the Netherlands
| | - Junye Yao
- Center for Brain Imaging Science and Technology, Key Laboratory for Biomedical Engineering of Ministry of Education, College of Biomedical Engineering and Instrumental Science, Zhejiang University, Zhouyiqing Building, Room 314, Yuquan Campus, Hangzhou 310027, China
| | - Hongjian He
- Center for Brain Imaging Science and Technology, Key Laboratory for Biomedical Engineering of Ministry of Education, College of Biomedical Engineering and Instrumental Science, Zhejiang University, Zhouyiqing Building, Room 314, Yuquan Campus, Hangzhou 310027, China.
| | - Keqing Zhu
- China Brain Bank and Department of Neurology in Second Affiliated Hospital, Key Laboratory of Medical Neurobiology of Zhejiang Province, and Department of Neurobiology, Zhejiang University School of Medicine, Hangzhou 310058, China; Department of Pathology, Zhejiang University School of Medicine, Hangzhou 310058, China.
| | - Jianhui Zhong
- Center for Brain Imaging Science and Technology, Key Laboratory for Biomedical Engineering of Ministry of Education, College of Biomedical Engineering and Instrumental Science, Zhejiang University, Zhouyiqing Building, Room 314, Yuquan Campus, Hangzhou 310027, China; Department of Imaging Sciences, University of Rochester, United States
| |
Collapse
|
17
|
Williamson NH, Ravin R, Cai TX, Benjamini D, Falgairolle M, O'Donovan MJ, Basser PJ. Real-time measurement of diffusion exchange rate in biological tissue. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2020; 317:106782. [PMID: 32679514 PMCID: PMC7427561 DOI: 10.1016/j.jmr.2020.106782] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/06/2020] [Revised: 06/22/2020] [Accepted: 06/26/2020] [Indexed: 05/05/2023]
Abstract
Diffusion exchange spectroscopy (DEXSY) provides a means to isolate the signal attenuation associated with exchange from other sources of signal loss. With the total diffusion weighting b1+b2=bs held constant, DEXSY signals acquired with b1=0 or b2=0 have no exchange weighting, while a DEXSY signal acquired with b1=b2 has maximal exchange weighting. The exchange rate can be estimated by fitting a diffusion exchange model to signals acquired with variable mixing times. Conventionally, acquired signals are normalized by a signal with b1=0 and b2=0 to remove the decay due to spin-lattice relaxation. Instead, division by a signal with equal bs but b1=0 or b2=0 reduces spin-lattice relaxation weighting of the apparent exchange rate (AXR). Furthermore, apparent diffusion-weighted R1 relaxation rates can be estimated from non-exchange-weighted DEXSY signals. Estimated R1 values are utilized to remove signal decay due to spin-lattice relaxation from exchange-weighted signals, permitting a more precise estimate of AXR with less data. Data reduction methods are proposed and tested with regards to statistical accuracy and precision of AXR estimates on simulated and experimental data. Simulations show that the methods are capable of accurately measuring the ground-truth exchange rate. The methods remain accurate even when the assumption that DEXSY signals attenuate with b is violated, as occurs for restricted diffusion. Experimental data was collected from fixed neonatal mouse spinal cord samples at 25 and 7°C using the strong static magnetic field gradient produced by a single-sided permanent magnet (i.e., an NMR MOUSE). The most rapid method for exchange measurements requires only five data points (an 80 s experiment as implemented) and achieves a similar level of accuracy and precision to the baseline method using 44 data points. This represents a significant improvement in acquisition speed, overcoming a barrier which has limited the use of DEXSY on living specimen.
Collapse
Affiliation(s)
- Nathan H Williamson
- National Institute of General Medical Sciences, National Institutes of Health, Bethesda, MD, USA; Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health, Bethesda, MD, USA.
| | - Rea Ravin
- Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health, Bethesda, MD, USA; Celoptics, Rockville, MD, USA
| | - Teddy X Cai
- Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health, Bethesda, MD, USA
| | - Dan Benjamini
- Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health, Bethesda, MD, USA; The Center for Neuroscience and Regenerative Medicine, Uniformed Service University of the Health Sciences, Bethesda, MD 20814, USA
| | - Melanie Falgairolle
- National Institute of Neurological Disorders and Stroke, National Institutes of Health, Bethesda, MD, USA
| | - Michael J O'Donovan
- National Institute of Neurological Disorders and Stroke, National Institutes of Health, Bethesda, MD, USA
| | - Peter J Basser
- Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health, Bethesda, MD, USA.
| |
Collapse
|
18
|
Ning L, Gagoski B, Szczepankiewicz F, Westin CF, Rathi Y. Joint RElaxation-Diffusion Imaging Moments to Probe Neurite Microstructure. IEEE TRANSACTIONS ON MEDICAL IMAGING 2020; 39:668-677. [PMID: 31398113 PMCID: PMC7164590 DOI: 10.1109/tmi.2019.2933982] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/23/2023]
Abstract
Joint relaxation-diffusion measurements can provide new insight about the tissue microstructural properties. Most recent methods have focused on inverting the Laplace transform to recover the joint distribution of relaxation-diffusion. However, as is well-known, this problem is notoriously ill-posed and numerically unstable. In this work, we address this issue by directly computing the joint moments of transverse relaxation rate and diffusivity, which can be robustly estimated. To zoom into different parts of the joint distribution, we further enhance our method by applying multiplicative filters to the joint probability density function of relaxation and diffusion and compute the corresponding moments. We propose an approach to use these moments to compute several novel scalar indices to characterize specific properties of the underlying tissue microstructure. Furthermore, for the first time, we propose an algorithm to estimate diffusion signals that are independent of echo time based on the moments of the marginal probability density function of diffusion. We demonstrate its utility in extracting tissue information not contaminated with multiple intra-voxel relaxation rates. We compare the performance of four types of filters that zoom into tissue components with different relaxation and diffusion properties and demonstrate it on an in-vivo human dataset. Experimental results show that these filters are able to characterize heterogeneous tissue microstructure. Moreover, the filtered diffusion signals are also able to distinguish fiber bundles with similar orientations but different relaxation rates. The proposed method thus allows to characterize the neural microstructure information in a robust and unique manner not possible using existing techniques.
Collapse
|
19
|
Nilsson M, Szczepankiewicz F, Brabec J, Taylor M, Westin CF, Golby A, van Westen D, Sundgren PC. Tensor-valued diffusion MRI in under 3 minutes: an initial survey of microscopic anisotropy and tissue heterogeneity in intracranial tumors. Magn Reson Med 2019; 83:608-620. [PMID: 31517401 PMCID: PMC6900060 DOI: 10.1002/mrm.27959] [Citation(s) in RCA: 44] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2019] [Revised: 07/05/2019] [Accepted: 07/30/2019] [Indexed: 12/23/2022]
Abstract
PURPOSE To evaluate the feasibility of a 3-minutes protocol for assessment of the microscopic anisotropy and tissue heterogeneity based on tensor-valued diffusion MRI in a wide range of intracranial tumors. METHODS B-tensor encoding was performed in 42 patients with intracranial tumors (gliomas, meningiomas, adenomas, and metastases). Microscopic anisotropy and tissue heterogeneity were evaluated by estimating the anisotropic kurtosis (MKA ) and isotropic kurtosis (MKI ), respectively. An extensive imaging protocol was compared with a 3-minutes protocol. RESULTS The fast imaging protocol yielded parameters with characteristics in terms of bias and precision similar to the full protocol. Glioblastomas had lower microscopic anisotropy than meningiomas (MKA = 0.29 ± 0.06 vs. 0.45 ± 0.08, P = 0.003). Metastases had higher tissue heterogeneity (MKI = 0.57 ± 0.07) than both the glioblastomas (0.44 ± 0.06, P < 0.001) and meningiomas (0.46 ± 0.06, P = 0.03). CONCLUSION Evaluation of the microscopic anisotropy and tissue heterogeneity in intracranial tumor patients is feasible in clinically relevant times frames.
Collapse
Affiliation(s)
- Markus Nilsson
- Department of Clinical Sciences Lund, Radiology, Lund University, Lund, Sweden
| | | | - Jan Brabec
- Department of Clinical Sciences Lund, Medical Radiation Physics, Lund University, Lund, Sweden
| | - Marie Taylor
- Department of Clinical Sciences Lund, Radiology, Lund University, Lund, Sweden
| | | | - Alexandra Golby
- Brigham and Women's Hospital, Harvard Medical School, Boston, Massachusetts
| | - Danielle van Westen
- Department of Clinical Sciences Lund, Radiology, Lund University, Lund, Sweden
| | - Pia C Sundgren
- Department of Clinical Sciences Lund, Radiology, Lund University, Lund, Sweden.,Lund University Bioimaging Center (LBIC), Lund University, Lund, Sweden
| |
Collapse
|
20
|
Brusini L, Menegaz G, Nilsson M. Monte Carlo Simulations of Water Exchange Through Myelin Wraps: Implications for Diffusion MRI. IEEE TRANSACTIONS ON MEDICAL IMAGING 2019; 38:1438-1445. [PMID: 30835213 DOI: 10.1109/tmi.2019.2894398] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
Diffusion magnetic resonance imaging (dMRI) yields parameters sensitive to brain tissue microstructure. A structurally important aspect of this microstructure is the myelin wrapping around the axons. This paper investigated the forward problem concerning whether water exchange via the spiraling structure of the myelin can meaningfully contribute to the signal in dMRI. Monte Carlo simulations were performed in a system with intra-axonal, myelin, and extra-axonal compartments. Diffusion in the myelin was simulated as a spiral wrapping the axon, with a custom number of wraps. Exchange (or intra-axonal residence) times were analyzed for various number of wraps and axon diameters. Pulsed gradient sequences were employed to simulate the dMRI signal, which was analyzed using different methods. Diffusional kurtosis imaging analysis yielded the radial diffusivity (RD) and radial kurtosis (RK), while the two-compartment Kärger model yielded estimates the intra-axonal volume fraction ( ν ic ) and exchange time ( τ ). Results showed that τ was on the sub-second level for geometries with axon diameters below 1.0 μ m and less than eight wraps. Otherwise, exchange was negligible compared to typical experimental durations, with τ of seconds or longer. In situations where exchange influenced the signal, estimates of RK and ν ic increased with the number of wraps, while RD decreased. τ estimates from simulated signals were in agreement with predicted ones. In conclusion, exchange through spiraling myelin permits sub-second τ for small diameters and low number of wraps. Such conditions may arise in the developing brain or in neurodegenerative disease, and thus the results could aid the interpretation of dMRI studies.
Collapse
|
21
|
Moutal N, Nilsson M, Topgaard D, Grebenkov D. The Kärger vs bi-exponential model: Theoretical insights and experimental validations. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2018; 296:72-78. [PMID: 30223153 DOI: 10.1016/j.jmr.2018.08.015] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/17/2018] [Revised: 08/28/2018] [Accepted: 08/30/2018] [Indexed: 06/08/2023]
Abstract
We revise three common models accounting for water exchange in pulsed-gradient spin-echo measurements: a bi-exponential model with time-dependent water fractions, the Kärger model, and a modified Kärger model designed for restricted diffusion, e.g. inside cells. The three models are compared and applied to experimental data from yeast cell suspensions. The Kärger model and the modified Kärger model yield very close results and accurately fit the data. The bi-exponential model, although less rigorous, has a natural physical interpretation and suggests a new experimental modality to estimate the water exchange time.
Collapse
Affiliation(s)
- Nicolas Moutal
- PMC, CNRS - Ecole Polytechnique, F-91128 Palaiseau, France.
| | - Markus Nilsson
- Physical Chemistry, Lund University, P.O.B. 124, SE-22100 Lund, Sweden
| | - Daniel Topgaard
- Physical Chemistry, Lund University, P.O.B. 124, SE-22100 Lund, Sweden
| | | |
Collapse
|