1
|
Wu D, Lee HH, Ba R, Turnbill V, Wang X, Luo Y, Walczak P, Fieremans E, Novikov DS, Martin LJ, Northington FJ, Zhang J. In vivo mapping of cellular resolution neuropathology in brain ischemia with diffusion MRI. SCIENCE ADVANCES 2024; 10:eadk1817. [PMID: 39018390 PMCID: PMC466947 DOI: 10.1126/sciadv.adk1817] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/08/2023] [Accepted: 06/11/2024] [Indexed: 07/19/2024]
Abstract
Noninvasive mapping of cellular pathology can provide critical diagnostic and prognostic information. Recent advances in diffusion magnetic resonance imaging enabled in vivo examination of tissue microstructures well beyond the imaging resolution. Here, we proposed to use diffusion time-dependent diffusion kurtosis imaging (tDKI) to simultaneously assess cellular morphology and transmembrane permeability in hypoxic-ischemic (HI) brain injury. Through numerical simulations and organoid imaging, we demonstrated the feasibility of capturing effective size and permeability changes using tDKI. In vivo MRI of HI-injured mouse brains detected a shift of the tDKI peak to longer diffusion times, suggesting swelling of the cellular processes. Furthermore, we observed a faster decrease of the tDKI tail, reflecting increased transmembrane permeability associated with up-regulated water exchange or necrosis. Such information, unavailable from a single diffusion time, can predict salvageable tissues. Preliminary applications of tDKI in patients with ischemic stroke suggested increased transmembrane permeability in stroke regions, illustrating tDKI's potential for detecting pathological changes in the clinics.
Collapse
Affiliation(s)
- Dan Wu
- Key Laboratory for Biomedical Engineering of Ministry of Education, Department of Biomedical Engineering, College of Biomedical Engineering and Instrument Science, Zhejiang University, Hangzhou, Zhejiang, China
- Department of Radiology, Children’s Hospital, Zhejiang University School of Medicine, National Clinical Research Center for Child Health, Hangzhou, Zhejiang, China
- Binjiang Institute of Zhejiang University, Hangzhou, Zhejiang, China
| | - Hong-Hsi Lee
- Department of Radiology, Massachusetts General Hospital, Boston, MA, USA
| | - Ruicheng Ba
- Key Laboratory for Biomedical Engineering of Ministry of Education, Department of Biomedical Engineering, College of Biomedical Engineering and Instrument Science, Zhejiang University, Hangzhou, Zhejiang, China
| | - Victoria Turnbill
- Department of Pediatrics, Johns Hopkins University School of Medicine, Baltimore, MD, USA
| | - Xiaoli Wang
- School of Medical Imaging, Weifang Medical School, Weifang, Shandong, China
| | - Yu Luo
- Department of Radiology, Shanghai Fourth People’s Hospital, School of Medicine, Tongji University, Shanghai, China
| | - Piotr Walczak
- Department of Radiology, University of Maryland, Baltimore, MD, USA
| | - Els Fieremans
- Center for Biomedical Imaging, Department of Radiology, New York University Grossman School of Medicine, New York, NY, USA
| | - Dmitry S. Novikov
- Center for Biomedical Imaging, Department of Radiology, New York University Grossman School of Medicine, New York, NY, USA
| | - Lee J. Martin
- Department of Pathology, Johns Hopkins University School of Medicine, Baltimore, MD, USA
| | - Frances J. Northington
- Department of Pediatrics, Johns Hopkins University School of Medicine, Baltimore, MD, USA
| | - Jiangyang Zhang
- Center for Biomedical Imaging, Department of Radiology, New York University Grossman School of Medicine, New York, NY, USA
| |
Collapse
|
2
|
Canales-Rodríguez EJ, Pizzolato M, Zhou FL, Barakovic M, Thiran JP, Jones DK, Parker GJM, Dyrby TB. Pore size estimation in axon-mimicking microfibers with diffusion-relaxation MRI. Magn Reson Med 2024; 91:2579-2596. [PMID: 38192108 DOI: 10.1002/mrm.29991] [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/17/2023] [Revised: 12/04/2023] [Accepted: 12/12/2023] [Indexed: 01/10/2024]
Abstract
PURPOSE This study aims to evaluate two distinct approaches for fiber radius estimation using diffusion-relaxation MRI data acquired in biomimetic microfiber phantoms that mimic hollow axons. The methods considered are the spherical mean power-law approach and a T2-based pore size estimation technique. THEORY AND METHODS A general diffusion-relaxation theoretical model for the spherical mean signal from water molecules within a distribution of cylinders with varying radii was introduced, encompassing the evaluated models as particular cases. Additionally, a new numerical approach was presented for estimating effective radii (i.e., MRI-visible mean radii) from the ground truth radii distributions, not reliant on previous theoretical approximations and adaptable to various acquisition sequences. The ground truth radii were obtained from scanning electron microscope images. RESULTS Both methods show a linear relationship between effective radii estimated from MRI data and ground-truth radii distributions, although some discrepancies were observed. The spherical mean power-law method overestimated fiber radii. Conversely, the T2-based method exhibited higher sensitivity to smaller fiber radii, but faced limitations in accurately estimating the radius in one particular phantom, possibly because of material-specific relaxation changes. CONCLUSION The study demonstrates the feasibility of both techniques to predict pore sizes of hollow microfibers. The T2-based technique, unlike the spherical mean power-law method, does not demand ultra-high diffusion gradients, but requires calibration with known radius distributions. This research contributes to the ongoing development and evaluation of neuroimaging techniques for fiber radius estimation, highlights the advantages and limitations of both methods, and provides datasets for reproducible research.
Collapse
Affiliation(s)
- Erick J Canales-Rodríguez
- Signal Processing Laboratory 5 (LTS5), Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
- Danish Research Centre for Magnetic Resonance (DRCMR), Centre for Functional and Diagnostic Imaging and Research, Copenhagen University Hospital Amager and Hvidovre, Copenhagen, Denmark
| | - Marco Pizzolato
- Danish Research Centre for Magnetic Resonance (DRCMR), Centre for Functional and Diagnostic Imaging and Research, Copenhagen University Hospital Amager and Hvidovre, Copenhagen, Denmark
- Department of Applied Mathematics and Computer Science, Technical University of Denmark (DTU), Kongens Lyngby, Denmark
| | - Feng-Lei Zhou
- Centre for Medical Image Computing, Department of Medical Physics and Biomedical Engineering, University College London (UCL), London, UK
- MicroPhantoms Limited, Cambridge, UK
| | - Muhamed Barakovic
- Translational Imaging in Neurology (ThINk) Basel, Department of Biomedical Engineering, University Hospital Basel and University of Basel, Basel, Switzerland
| | - Jean-Philippe Thiran
- Signal Processing Laboratory 5 (LTS5), Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
- Radiology Department, Centre Hospitalier Universitaire Vaudois and University of Lausanne, Lausanne, Switzerland
- Centre d'Imagerie Biomédicale (CIBM), EPFL, Lausanne, Switzerland
| | - Derek K Jones
- Cardiff University Brain Research Imaging Centre (CUBRIC), Cardiff University, Cardiff, UK
| | - Geoffrey J M Parker
- Centre for Medical Image Computing, Department of Medical Physics and Biomedical Engineering, University College London (UCL), London, UK
- Department of Neuroinflammation, Queen Square Institute of Neurology, University College London (UCL), London, UK
- Bioxydyn Limited, Manchester, UK
| | - Tim B Dyrby
- Danish Research Centre for Magnetic Resonance (DRCMR), Centre for Functional and Diagnostic Imaging and Research, Copenhagen University Hospital Amager and Hvidovre, Copenhagen, Denmark
- Department of Applied Mathematics and Computer Science, Technical University of Denmark (DTU), Kongens Lyngby, Denmark
| |
Collapse
|
3
|
Shen Y, Zhao X, Wang K, Sun Y, Zhang X, Wang C, Yang Z, Feng Z, Zhang X. Exploring White Matter Abnormalities in Young Children with Autism Spectrum Disorder: Integrating Multi-shell Diffusion Data and Machine Learning Analysis. Acad Radiol 2024; 31:2074-2084. [PMID: 38185571 DOI: 10.1016/j.acra.2023.12.023] [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: 11/18/2023] [Revised: 12/09/2023] [Accepted: 12/14/2023] [Indexed: 01/09/2024]
Abstract
RATIONALE AND OBJECTIVES This study employed tract-based spatial statistics (TBSS) to investigate abnormalities in the white matter microstructure among children with autism spectrum disorder (ASD). Additionally, an eXtreme Gradient Boosting (XGBoost) model was developed to effectively classify individuals with ASD and typical developing children (TDC). METHODS AND MATERIALS Multi-shell diffusion weighted images were acquired from 62 children with ASD and 44 TDC. Using the Pydesigner procedure, diffusion tensor (DT), diffusion kurtosis (DK), and white matter tract integrity (WMTI) metrics were computed. Subsequently, TBSS analysis was applied to discern differences in these diffusion parameters between ASD and TDC groups. The XGBoost model was then trained using metrics showing significant differences, and Shapley Additive explanations (SHAP) values were computed to assess the feature importance in the model's predictions. RESULTS TBSS analysis revealed a significant reduction in axonal diffusivity (AD) in the left posterior corona radiata and the right superior corona radiata. Among the DK indicators, mean kurtosis, axial kurtosis, and kurtosis fractional anisotropy were notably increased in children with ASD, with no significant difference in radial kurtosis. WMTI metrics such as axonal water fraction, axonal diffusivity of the extra-axonal space (EAS_AD), tortuosity of the extra-axonal space (EAS_TORT), and diffusivity of intra-axonal space (IAS_Da) were significantly increased, primarily in the corpus callosum and fornix. Notably, there was no significant difference in radial diffusivity of the extra-axial space (EAS_RD). The XGBoost model demonstrated excellent classification ability, and the SHAP analysis identified EAS_TORT as the feature with the highest importance in the model's predictions. CONCLUSION This study utilized TBSS analyses with multi-shell diffusion data to examine white matter abnormalities in pediatric autism. Additionally, the developed XGBoost model showed outstanding performance in classifying ASD and TDC. The ranking of SHAP values based on the XGBoost model underscored the significance of features in influencing model predictions.
Collapse
Affiliation(s)
- Yanyong Shen
- Department of Radiology, The Third Affiliated Hospital of Zhengzhou University, Zhengzhou, 450052, China (Y.S., X.Z., X.Z., C.W., Z.Y., Z.F., X.Z.); Henan International Joint Laboratory of Neuroimaging, Zhengzhou, 450052, China (Y.S., X.Z., X.Z., C.W., Z.Y., Z.F., X.Z.)
| | - Xin Zhao
- Department of Radiology, The Third Affiliated Hospital of Zhengzhou University, Zhengzhou, 450052, China (Y.S., X.Z., X.Z., C.W., Z.Y., Z.F., X.Z.); Henan International Joint Laboratory of Neuroimaging, Zhengzhou, 450052, China (Y.S., X.Z., X.Z., C.W., Z.Y., Z.F., X.Z.)
| | - Kaiyu Wang
- MR Research China, GE Healthcare, Beijing, 100000, PR China (K.W.)
| | - Yongbing Sun
- Department of Radiology, Henan Provincial People's Hospital, Zhengzhou, 450000, China (Y.S.)
| | - Xiaoxue Zhang
- Department of Radiology, The Third Affiliated Hospital of Zhengzhou University, Zhengzhou, 450052, China (Y.S., X.Z., X.Z., C.W., Z.Y., Z.F., X.Z.); Henan International Joint Laboratory of Neuroimaging, Zhengzhou, 450052, China (Y.S., X.Z., X.Z., C.W., Z.Y., Z.F., X.Z.)
| | - Changhao Wang
- Department of Radiology, The Third Affiliated Hospital of Zhengzhou University, Zhengzhou, 450052, China (Y.S., X.Z., X.Z., C.W., Z.Y., Z.F., X.Z.); Henan International Joint Laboratory of Neuroimaging, Zhengzhou, 450052, China (Y.S., X.Z., X.Z., C.W., Z.Y., Z.F., X.Z.)
| | - Zhexuan Yang
- Department of Radiology, The Third Affiliated Hospital of Zhengzhou University, Zhengzhou, 450052, China (Y.S., X.Z., X.Z., C.W., Z.Y., Z.F., X.Z.); Henan International Joint Laboratory of Neuroimaging, Zhengzhou, 450052, China (Y.S., X.Z., X.Z., C.W., Z.Y., Z.F., X.Z.)
| | - Zhanqi Feng
- Department of Radiology, The Third Affiliated Hospital of Zhengzhou University, Zhengzhou, 450052, China (Y.S., X.Z., X.Z., C.W., Z.Y., Z.F., X.Z.); Henan International Joint Laboratory of Neuroimaging, Zhengzhou, 450052, China (Y.S., X.Z., X.Z., C.W., Z.Y., Z.F., X.Z.)
| | - Xiaoan Zhang
- Department of Radiology, The Third Affiliated Hospital of Zhengzhou University, Zhengzhou, 450052, China (Y.S., X.Z., X.Z., C.W., Z.Y., Z.F., X.Z.); Henan International Joint Laboratory of Neuroimaging, Zhengzhou, 450052, China (Y.S., X.Z., X.Z., C.W., Z.Y., Z.F., X.Z.).
| |
Collapse
|
4
|
Zhou M, Stobbe R, Szczepankiewicz F, Budde M, Buck B, Kate M, Lloret M, Fairall P, Butcher K, Shuaib A, Emery D, Nilsson M, Westin CF, Beaulieu C. Tensor-valued diffusion MRI of human acute stroke. Magn Reson Med 2024; 91:2126-2141. [PMID: 38156813 DOI: 10.1002/mrm.29975] [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: 08/08/2023] [Revised: 11/18/2023] [Accepted: 11/27/2023] [Indexed: 01/03/2024]
Abstract
PURPOSE Tensor-valued diffusion encoding can disentangle orientation dispersion and subvoxel anisotropy, potentially offering insight into microstructural changes after cerebral ischemia. The purpose was to evaluate tensor-valued diffusion MRI in human acute ischemic stroke, assess potential confounders from diffusion time dependencies, and compare to Monte Carlo diffusion simulations of axon beading. METHODS Linear (LTE) and spherical (STE) b-tensor encoding with inherently different effective diffusion times were acquired in 21 acute ischemic stroke patients between 3 and 57 h post-onset at 3 T in 2.5 min. In an additional 10 patients, STE with 2 LTE yielding different effective diffusion times were acquired for comparison. Diffusional variance decomposition (DIVIDE) was used to estimate microscopic anisotropy (μFA), as well as anisotropic, isotropic, and total diffusional variance (MKA , MKI , MKT ). DIVIDE parameters, and diffusion tensor imaging (DTI)-derived mean diffusivity and fractional anisotropy (FA) were compared in lesion versus contralateral white matter. Monte Carlo diffusion simulations of various cylindrical geometries for all b-tensor protocols were used to interpret parameter measurements. RESULTS MD was ˜40% lower in lesions for all LTE/STE protocols. The DIVIDE parameters varied with effective diffusion time: higher μFA and MKA in lesion versus contralateral white matter for STE with longer effective diffusion time LTE, whereas the shorter effective diffusion time LTE protocol yielded lower μFA and MKA in lesions. Both protocols, regardless of diffusion time, were consistent with simulations of greater beading amplitude and intracellular volume fraction. CONCLUSION DIVIDE parameters depend on diffusion time in acute stroke but consistently indicate neurite beading and larger intracellular volume fraction.
Collapse
Affiliation(s)
- Mi Zhou
- Biomedical Engineering, University of Alberta, Edmonton, Alberta, Canada
| | - Robert Stobbe
- Biomedical Engineering, University of Alberta, Edmonton, Alberta, Canada
- Radiology and Diagnostic Imaging, University of Alberta, Edmonton, Alberta, Canada
| | | | - Matthew Budde
- Neurosurgery, Medical College of Wisconsin, Milwaukee, Wisconsin, USA
| | - Brian Buck
- Neurology, University of Alberta, Edmonton, Alberta, Canada
| | - Mahesh Kate
- Neurology, University of Alberta, Edmonton, Alberta, Canada
| | - Mar Lloret
- Neurology, University of Alberta, Edmonton, Alberta, Canada
| | - Paige Fairall
- Neurology, University of Alberta, Edmonton, Alberta, Canada
| | - Ken Butcher
- School of Clinical Medicine, University of New South Wales, Sydney, New South Wales, Australia
| | - Ashfaq Shuaib
- Neurology, University of Alberta, Edmonton, Alberta, Canada
| | - Derek Emery
- Radiology and Diagnostic Imaging, University of Alberta, Edmonton, Alberta, Canada
| | - Markus Nilsson
- Clinical Sciences Lund, Lund University, Lund, Scania, Sweden
| | - Carl-Fredrik Westin
- Radiology, Brigham and Women's Hospital, Harvard Medical School, Boston, MA, USA
| | - Christian Beaulieu
- Biomedical Engineering, University of Alberta, Edmonton, Alberta, Canada
- Radiology and Diagnostic Imaging, University of Alberta, Edmonton, Alberta, Canada
| |
Collapse
|
5
|
Dong Z, Reese TG, Lee HH, Huang SY, Polimeni JR, Wald LL, Wang F. Romer-EPTI: rotating-view motion-robust super-resolution EPTI for SNR-efficient distortion-free in-vivo mesoscale dMRI and microstructure imaging. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.01.26.577343. [PMID: 38352481 PMCID: PMC10862730 DOI: 10.1101/2024.01.26.577343] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 02/19/2024]
Abstract
Purpose To overcome the major challenges in dMRI acquisition, including low SNR, distortion/blurring, and motion vulnerability. Methods A novel Romer-EPTI technique is developed to provide distortion-free dMRI with significant SNR gain, high motion-robustness, sharp spatial resolution, and simultaneous multi-TE imaging. It introduces a ROtating-view Motion-robust supEr-Resolution technique (Romer) combined with a distortion/blurring-free EPTI encoding. Romer enhances SNR by a simultaneous multi-thick-slice acquisition with rotating-view encoding, while providing high motion-robustness through a motion-aware super-resolution reconstruction, which also incorporates slice-profile and real-value diffusion, to resolve high-isotropic-resolution volumes. The in-plane encoding is performed using distortion/blurring-free EPTI, which further improves effective spatial resolution and motion robustness by preventing not only T2/T2*-blurring but also additional blurring resulting from combining encoded volumes with inconsistent geometries caused by dynamic distortions. Self-navigation was incorporated to enable efficient phase correction. Additional developments include strategies to address slab-boundary artifacts, achieve minimal TE for SNR gain at 7T, and achieve high robustness to strong phase variations at high b-values. Results Using Romer-EPTI, we demonstrate distortion-free whole-brain mesoscale in-vivo dMRI at both 3T (500-μm-iso) and 7T (485-μm-iso) for the first time, with high SNR efficiency (e.g., 25 × ), and high image quality free from distortion and slab-boundary artifacts with minimal blurring. Motion experiments demonstrate Romer-EPTI's high motion-robustness and ability to recover sharp images in the presence of motion. Romer-EPTI also demonstrates significant SNR gain and robustness in high b-value (b=5000s/mm2) and time-dependent dMRI. Conclusion Romer-EPTI significantly improves SNR, motion-robustness, and image quality, providing a highly efficient acquisition for high-resolution dMRI and microstructure imaging.
Collapse
Affiliation(s)
- Zijing Dong
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Charlestown, Massachusetts, USA
- Department of Radiology, Harvard Medical School, Boston, Massachusetts, USA
| | - Timothy G. Reese
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Charlestown, Massachusetts, USA
- Department of Radiology, Harvard Medical School, Boston, Massachusetts, USA
| | - Hong-Hsi Lee
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Charlestown, Massachusetts, USA
- Department of Radiology, Harvard Medical School, Boston, Massachusetts, USA
| | - Susie Y. Huang
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Charlestown, Massachusetts, USA
- Department of Radiology, Harvard Medical School, Boston, Massachusetts, USA
- Harvard-MIT Health Sciences and Technology, MIT, Cambridge, Massachusetts, USA
| | - Jonathan R. Polimeni
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Charlestown, Massachusetts, USA
- Department of Radiology, Harvard Medical School, Boston, Massachusetts, USA
- Harvard-MIT Health Sciences and Technology, MIT, Cambridge, Massachusetts, USA
| | - Lawrence L. Wald
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Charlestown, Massachusetts, USA
- Department of Radiology, Harvard Medical School, Boston, Massachusetts, USA
- Harvard-MIT Health Sciences and Technology, MIT, Cambridge, Massachusetts, USA
| | - Fuyixue Wang
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Charlestown, Massachusetts, USA
- Department of Radiology, Harvard Medical School, Boston, Massachusetts, USA
| |
Collapse
|
6
|
Zhu A, Shih R, Huang RY, DeMarco JK, Bhushan C, Morris HD, Kohls G, Yeo DTB, Marinelli L, Mitra J, Hood M, Ho VB, Foo TKF. Revealing tumor microstructure with oscillating diffusion encoding MRI in pre-surgical and post-treatment glioma patients. Magn Reson Med 2023; 90:1789-1801. [PMID: 37335831 DOI: 10.1002/mrm.29758] [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: 03/06/2023] [Revised: 05/09/2023] [Accepted: 05/24/2023] [Indexed: 06/21/2023]
Abstract
PURPOSE We hypothesized that the time-dependent diffusivity at short diffusion times, as measured by oscillating gradient spin echo (OGSE) diffusion MRI, can characterize tissue microstructures in glioma patients. THEORY AND METHODS Five adult patients with known diffuse glioma, including two pre-surgical and three with new enhancing lesions after treatment for high-grade glioma, were scanned in an ultra-high-performance gradient 3.0T MRI system. OGSE diffusion MRI at 30-100 Hz and pulsed gradient spin echo diffusion imaging (approximated as 0 Hz) were obtained. The ADC and trace-diffusion-weighted image at each acquired frequency were calculated, that is, ADC (f) and TraceDWI (f). RESULTS In pre-surgical patients, biopsy-confirmed solid enhancing tumor in a high-grade glioblastoma showed higherADC ( f ) ADC ( 0 Hz ) $$ \frac{\mathrm{ADC}\ (f)}{\mathrm{ADC}\ \left(0\ \mathrm{Hz}\right)} $$ and lowerTraceDWI ( f ) TraceDWI ( 0 Hz ) $$ \frac{\mathrm{TraceDWI}\ (f)}{\mathrm{TraceDWI}\ \left(0\ \mathrm{Hz}\right)} $$ , compared to that at same OGSE frequency in a low-grade astrocytoma. In post-treatment patients, the enhancing lesions of two patients who were diagnosed with tumor progression contained more voxels with highADC ( f ) ADC ( 0 Hz ) $$ \frac{\mathrm{ADC}\ (f)}{\mathrm{ADC}\ \left(0\ \mathrm{Hz}\right)} $$ and lowTraceDWI ( f ) TraceDWI ( 0 Hz ) $$ \frac{\mathrm{TraceDWI}\left(\mathrm{f}\right)}{\mathrm{TraceDWI}\left(0\ \mathrm{Hz}\right)} $$ , compared to the enhancing lesions of a patient who was diagnosed with treatment effect. Non-enhancing T2 signal abnormality lesions in both the pre-surgical high-grade glioblastoma and post-treatment tumor progressions showed regions with highADC ( f ) ADC ( 0 Hz ) $$ \frac{\mathrm{ADC}\ (f)}{\mathrm{ADC}\ \left(0\ \mathrm{Hz}\right)} $$ and lowTraceDWI ( f ) TraceDWI ( 0 Hz ) $$ \frac{\mathrm{TraceDWI}\ \left(\mathrm{f}\right)}{\mathrm{TraceDWI}\ \left(0\ \mathrm{Hz}\right)} $$ , consistent with infiltrative tumor. The solid tumor of the glioblastoma, the enhancing lesions of post-treatment tumor progressions, and the suspected infiltrative tumors showed high diffusion time-dependency from 30 to 100 Hz, consistent with high intra-tumoral volume fraction (cellular density). CONCLUSION Different characteristics of OGSE-based time-dependent diffusivity can reveal heterogenous tissue microstructures that indicate cellular density in glioma patients.
Collapse
Affiliation(s)
- Ante Zhu
- GE Research, Niskayuna, New York, USA
| | - Robert Shih
- Uniformed Services University, Bethesda, Maryland, USA
- Walter Reed National Military Medical Center, Bethesda, Maryland, USA
| | - Raymond Y Huang
- Department of Radiology, Brigham and Women's Hospital, Boston, Massachusetts, USA
| | - J Kevin DeMarco
- Uniformed Services University, Bethesda, Maryland, USA
- Walter Reed National Military Medical Center, Bethesda, Maryland, USA
| | | | - H Douglas Morris
- Uniformed Services University, Bethesda, Maryland, USA
- Walter Reed National Military Medical Center, Bethesda, Maryland, USA
| | - Gail Kohls
- Uniformed Services University, Bethesda, Maryland, USA
- Walter Reed National Military Medical Center, Bethesda, Maryland, USA
| | | | | | | | - Maureen Hood
- Uniformed Services University, Bethesda, Maryland, USA
- Walter Reed National Military Medical Center, Bethesda, Maryland, USA
| | - Vincent B Ho
- Uniformed Services University, Bethesda, Maryland, USA
- Walter Reed National Military Medical Center, Bethesda, Maryland, USA
| | - Thomas K F Foo
- GE Research, Niskayuna, New York, USA
- Uniformed Services University, Bethesda, Maryland, USA
| |
Collapse
|
7
|
Boito D, Eklund A, Tisell A, Levi R, Özarslan E, Blystad I. MRI with generalized diffusion encoding reveals damaged white matter in patients previously hospitalized for COVID-19 and with persisting symptoms at follow-up. Brain Commun 2023; 5:fcad284. [PMID: 37953843 PMCID: PMC10638510 DOI: 10.1093/braincomms/fcad284] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2023] [Revised: 08/25/2023] [Accepted: 10/26/2023] [Indexed: 11/14/2023] Open
Abstract
There is mounting evidence of the long-term effects of COVID-19 on the central nervous system, with patients experiencing diverse symptoms, often suggesting brain involvement. Conventional brain MRI of these patients shows unspecific patterns, with no clear connection of the symptomatology to brain tissue abnormalities, whereas diffusion tensor studies and volumetric analyses detect measurable changes in the brain after COVID-19. Diffusion MRI exploits the random motion of water molecules to achieve unique sensitivity to structures at the microscopic level, and new sequences employing generalized diffusion encoding provide structural information which are sensitive to intravoxel features. In this observational study, a total of 32 persons were investigated: 16 patients previously hospitalized for COVID-19 with persisting symptoms of post-COVID condition (mean age 60 years: range 41-79, all male) at 7-month follow-up and 16 matched controls, not previously hospitalized for COVID-19, with no post-COVID symptoms (mean age 58 years, range 46-69, 11 males). Standard MRI and generalized diffusion encoding MRI were employed to examine the brain white matter of the subjects. To detect possible group differences, several tissue microstructure descriptors obtainable with the employed diffusion sequence, the fractional anisotropy, mean diffusivity, axial diffusivity, radial diffusivity, microscopic anisotropy, orientational coherence (Cc) and variance in compartment's size (CMD) were analysed using the tract-based spatial statistics framework. The tract-based spatial statistics analysis showed widespread statistically significant differences (P < 0.05, corrected for multiple comparisons using the familywise error rate) in all the considered metrics in the white matter of the patients compared to the controls. Fractional anisotropy, microscopic anisotropy and Cc were lower in the patient group, while axial diffusivity, radial diffusivity, mean diffusivity and CMD were higher. Significant changes in fractional anisotropy, microscopic anisotropy and CMD affected approximately half of the analysed white matter voxels located across all brain lobes, while changes in Cc were mainly found in the occipital parts of the brain. Given the predominant alteration in microscopic anisotropy compared to Cc, the observed changes in diffusion anisotropy are mostly due to loss of local anisotropy, possibly connected to axonal damage, rather than white matter fibre coherence disruption. The increase in radial diffusivity is indicative of demyelination, while the changes in mean diffusivity and CMD are compatible with vasogenic oedema. In summary, these widespread alterations of white matter microstructure are indicative of vasogenic oedema, demyelination and axonal damage. These changes might be a contributing factor to the diversity of central nervous system symptoms that many patients experience after COVID-19.
Collapse
Affiliation(s)
- Deneb Boito
- Department of Biomedical Engineering, Linköping University, S-58183 Linköping, Sweden
- Centre for Medical Image Science and Visualization (CMIV), Linköping University, S-58183 Linköping, Sweden
| | - Anders Eklund
- Department of Biomedical Engineering, Linköping University, S-58183 Linköping, Sweden
- Centre for Medical Image Science and Visualization (CMIV), Linköping University, S-58183 Linköping, Sweden
- Division of Statistics and Machine learning, Department of Computer and Information Science, Linköping University, S-58183 Linköping, Sweden
| | - Anders Tisell
- Centre for Medical Image Science and Visualization (CMIV), Linköping University, S-58183 Linköping, Sweden
- Department of Radiation Physics, Linköping University, S-58185 Linköping, Sweden
- Department of Health, Medicine and Caring Sciences, Linköping University, S58183 Linköping, Sweden
| | - Richard Levi
- Centre for Medical Image Science and Visualization (CMIV), Linköping University, S-58183 Linköping, Sweden
- Department of Health, Medicine and Caring Sciences, Linköping University, S58183 Linköping, Sweden
- Department of Rehabilitation Medicine in Linköping, Linköping University, S-58185 Linköping, Sweden
| | - Evren Özarslan
- Department of Biomedical Engineering, Linköping University, S-58183 Linköping, Sweden
- Centre for Medical Image Science and Visualization (CMIV), Linköping University, S-58183 Linköping, Sweden
| | - Ida Blystad
- Centre for Medical Image Science and Visualization (CMIV), Linköping University, S-58183 Linköping, Sweden
- Department of Health, Medicine and Caring Sciences, Linköping University, S58183 Linköping, Sweden
- Department of Radiology in Linköping, Linköping University, S-58185 Linköping, Sweden
| |
Collapse
|
8
|
Dan G, Sun K, Luo Q, Zhou XJ. Time-dependent diffusion MRI using multiple stimulated echoes. Magn Reson Med 2023; 90:910-921. [PMID: 37103885 PMCID: PMC10330017 DOI: 10.1002/mrm.29677] [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/05/2022] [Revised: 03/31/2023] [Accepted: 04/03/2023] [Indexed: 04/28/2023]
Abstract
PURPOSE To develop a time-efficient pulse sequence that acquires multiple diffusion-weighted images with distinct diffusion times in a single shot by using multiple stimulated echoes (mSTE) with variable flip angles (VFA). METHODS The proposed diffusion-weighted mSTE with VFA (DW-mSTE-VFA) sequence begins with two 90° RF pulses that straddle a diffusion gradient lobe (GD ) to excite and restore one half of the magnetization into the longitudinal axis. The restored longitudinal magnetization was successively re-excited by a series of RF pulses with VFA, each followed by another GD , to generate a set of stimulated echoes. Each of the multiple stimulated echoes was acquired with an EPI echo train. As such, the train of multiple stimulated echoes produced a set of diffusion-weighted images with varying diffusion times in a single shot. This technique was experimentally demonstrated on a diffusion phantom, a fruit, and healthy human brain and prostate at 3 T. RESULTS In the phantom experiment, the mean ADC measured at different diffusion times using DW-mSTE-VFA were highly consistent (r = 0.999) with those from a commercial spin-echo diffusion-weighted EPI sequence. In the fruit and brain experiments, DW-mSTE-VFA exhibited similar diffusion-time dependence to a standard diffusion-weighted stimulated echo sequence. The ADC showed significant time dependence in the human brain (p = 0.003 in both white matter and gray matter) and prostate tissues (p = 0.003 in both peripheral zone and central gland). CONCLUSION DW-mSTE-VFA offers a time-efficient tool for investigating the diffusion-time dependency in diffusion MRI studies.
Collapse
Affiliation(s)
- Guangyu Dan
- Center for Magnetic Resonance Research, University of Illinois College of Medicine at Chicago, Chicago, IL, United States
- Department of Biomedical Engineering, University of Illinois at Chicago, Chicago, IL, United States
| | - Kaibao Sun
- Center for Magnetic Resonance Research, University of Illinois College of Medicine at Chicago, Chicago, IL, United States
| | - Qingfei Luo
- Center for Magnetic Resonance Research, University of Illinois College of Medicine at Chicago, Chicago, IL, United States
| | - Xiaohong Joe Zhou
- Center for Magnetic Resonance Research, University of Illinois College of Medicine at Chicago, Chicago, IL, United States
- Department of Biomedical Engineering, University of Illinois at Chicago, Chicago, IL, United States
- Departments of Radiology and Neurosurgery, University of Illinois College of Medicine at Chicago, Chicago, IL, United States
| |
Collapse
|
9
|
Warner W, Palombo M, Cruz R, Callaghan R, Shemesh N, Jones DK, Dell'Acqua F, Ianus A, Drobnjak I. Temporal Diffusion Ratio (TDR) for imaging restricted diffusion: Optimisation and pre-clinical demonstration. Neuroimage 2023; 269:119930. [PMID: 36750150 PMCID: PMC7615244 DOI: 10.1016/j.neuroimage.2023.119930] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2022] [Revised: 01/12/2023] [Accepted: 02/02/2023] [Indexed: 02/07/2023] Open
Abstract
Temporal Diffusion Ratio (TDR) is a recently proposed dMRI technique (Dell'Acqua et al., proc. ISMRM 2019) which provides contrast between areas with restricted diffusion and areas either without restricted diffusion or with length scales too small for characterisation. Hence, it has a potential for informing on pore sizes, in particular the presence of large axon diameters or other cellular structures. TDR employs the signal from two dMRI acquisitions obtained with the same, large, b-value but with different diffusion gradient waveforms. TDR is advantageous as it employs standard acquisition sequences, does not make any assumptions on the underlying tissue structure and does not require any model fitting, avoiding issues related to model degeneracy. This work for the first time introduces and optimises the TDR method in simulation for a range of different tissues and scanner constraints and validates it in a pre-clinical demonstration. We consider both substrates containing cylinders and spherical structures, representing cell soma in tissue. Our results show that contrasting an acquisition with short gradient duration, short diffusion time and high gradient strength with an acquisition with long gradient duration, long diffusion time and low gradient strength, maximises the TDR contrast for a wide range of pore configurations. Additionally, in the presence of Rician noise, computing TDR from a subset (50% or fewer) of the acquired diffusion gradients rather than the entire shell as proposed originally further improves the contrast. In the last part of the work the results are demonstrated experimentally on rat spinal cord. In line with simulations, the experimental data shows that optimised TDR improves the contrast compared to non-optimised TDR. Furthermore, we find a strong correlation between TDR and histology measurements of axon diameter. In conclusion, we find that TDR has great potential and is a very promising alternative (or potentially complement) to model-based approaches for informing on pore sizes and restricted diffusion in general.
Collapse
Affiliation(s)
- William Warner
- Centre for Medical Image Computing (CMIC), Computer Science Department, University College London, United Kingdom
| | - Marco Palombo
- Cardiff University Brain Research Imaging Centre, School of Psychology, Cardiff University, Cardiff, United Kingdom; School of Computer Science and Informatics, Cardiff University, Cardiff, United Kingdom
| | - Renata Cruz
- Champalimaud Research, Champalimaud Foundation, Lisbon, Portugal
| | | | - Noam Shemesh
- Champalimaud Research, Champalimaud Foundation, Lisbon, Portugal
| | - Derek K Jones
- Cardiff University Brain Research Imaging Centre, School of Psychology, Cardiff University, Cardiff, United Kingdom
| | - Flavio Dell'Acqua
- NatBrainLab, Institute of Psychiatry, Psychology and Neuroscience, King's College London, London, United Kingdom
| | - Andrada Ianus
- Champalimaud Research, Champalimaud Foundation, Lisbon, Portugal.
| | - Ivana Drobnjak
- Centre for Medical Image Computing (CMIC), Computer Science Department, University College London, United Kingdom.
| |
Collapse
|
10
|
Pizzolato M, Canales-Rodríguez EJ, Andersson M, Dyrby TB. Axial and radial axonal diffusivities and radii from single encoding strongly diffusion-weighted MRI. Med Image Anal 2023; 86:102767. [PMID: 36867913 DOI: 10.1016/j.media.2023.102767] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2022] [Revised: 12/13/2022] [Accepted: 02/08/2023] [Indexed: 02/18/2023]
Abstract
We enable the estimation of the per-axon axial diffusivity from single encoding, strongly diffusion-weighted, pulsed gradient spin echo data. Additionally, we improve the estimation of the per-axon radial diffusivity compared to estimates based on spherical averaging. The use of strong diffusion weightings in magnetic resonance imaging (MRI) allows to approximate the signal in white matter as the sum of the contributions from only axons. At the same time, spherical averaging leads to a major simplification of the modeling by removing the need to explicitly account for the unknown distribution of axonal orientations. However, the spherically averaged signal acquired at strong diffusion weightings is not sensitive to the axial diffusivity, which cannot therefore be estimated although needed for modeling axons - especially in the context of multi-compartmental modeling. We introduce a new general method for the estimation of both the axial and radial axonal diffusivities at strong diffusion weightings based on kernel zonal modeling. The method could lead to estimates that are free from partial volume bias with gray matter or other isotropic compartments. The method is tested on publicly available data from the MGH Adult Diffusion Human Connectome project. We report reference values of axonal diffusivities based on 34 subjects, and derive estimates of axonal radii from only two shells. The estimation problem is also addressed from the angle of the required data preprocessing, the presence of biases related to modeling assumptions, current limitations, and future possibilities.
Collapse
Affiliation(s)
- Marco Pizzolato
- Department of Applied Mathematics and Computer Science, Technical University of Denmark, Kgs. Lyngby, Denmark; Danish Research Centre for Magnetic Resonance, Centre for Functional and Diagnostic Imaging and Research, Copenhagen University Hospital Amager and Hvidovre, Copenhagen, Denmark.
| | | | - Mariam Andersson
- Danish Research Centre for Magnetic Resonance, Centre for Functional and Diagnostic Imaging and Research, Copenhagen University Hospital Amager and Hvidovre, Copenhagen, Denmark
| | - Tim B Dyrby
- Department of Applied Mathematics and Computer Science, Technical University of Denmark, Kgs. Lyngby, Denmark; Danish Research Centre for Magnetic Resonance, Centre for Functional and Diagnostic Imaging and Research, Copenhagen University Hospital Amager and Hvidovre, Copenhagen, Denmark
| |
Collapse
|
11
|
Howard AF, Cottaar M, Drakesmith M, Fan Q, Huang SY, Jones DK, Lange FJ, Mollink J, Rudrapatna SU, Tian Q, Miller KL, Jbabdi S. Estimating axial diffusivity in the NODDI model. Neuroimage 2022; 262:119535. [PMID: 35931306 PMCID: PMC9802007 DOI: 10.1016/j.neuroimage.2022.119535] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2022] [Revised: 07/20/2022] [Accepted: 08/01/2022] [Indexed: 01/03/2023] Open
Abstract
To estimate microstructure-related parameters from diffusion MRI data, biophysical models make strong, simplifying assumptions about the underlying tissue. The extent to which many of these assumptions are valid remains an open research question. This study was inspired by the disparity between the estimated intra-axonal axial diffusivity from literature and that typically assumed by the Neurite Orientation Dispersion and Density Imaging (NODDI) model (d∥=1.7μm2/ms). We first demonstrate how changing the assumed axial diffusivity results in considerably different NODDI parameter estimates. Second, we illustrate the ability to estimate axial diffusivity as a free parameter of the model using high b-value data and an adapted NODDI framework. Using both simulated and in vivo data we investigate the impact of fitting to either real-valued or magnitude data, with Gaussian and Rician noise characteristics respectively, and what happens if we get the noise assumptions wrong in this high b-value and thus low SNR regime. Our results from real-valued human data estimate intra-axonal axial diffusivities of ∼2-2.5μm2/ms, in line with current literature. Crucially, our results demonstrate the importance of accounting for both a rectified noise floor and/or a signal offset to avoid biased parameter estimates when dealing with low SNR data.
Collapse
Affiliation(s)
- Amy Fd Howard
- FMRIB Centre, Wellcome Centre for Integrative Neuroimaging, Nuffield Department of Clinical Neurosciences, University of Oxford, Oxford, United Kingdom.
| | - Michiel Cottaar
- FMRIB Centre, Wellcome Centre for Integrative Neuroimaging, Nuffield Department of Clinical Neurosciences, University of Oxford, Oxford, United Kingdom
| | - Mark Drakesmith
- Cardiff University Brain Research Imaging Centre, Cardiff University, Cardiff, United Kingdom; Neuroscience and Mental Health Research Institute, Cardiff University, Cardiff, United Kingdom
| | - Qiuyun Fan
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, Massachusetts, United States; Harvard Medical School, Boston, Massachusetts, United States; Department of Biomedical Engineering, College of Precision Instruments and Optoelectronics Engineering, Tianjin University, Tianjin, China; Academy of Medical Engineering and Translational Medicine, Medical College, Tianjin University, Tianjin, China
| | - Susie Y Huang
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, Massachusetts, United States; Harvard Medical School, Boston, Massachusetts, United States; Harvard-MIT Division of Health Sciences and Technology, Massachusetts Institute of Technology, Cambridge, Massachusetts, United States
| | - Derek K Jones
- Cardiff University Brain Research Imaging Centre, Cardiff University, Cardiff, United Kingdom; Neuroscience and Mental Health Research Institute, Cardiff University, Cardiff, United Kingdom
| | - Frederik J Lange
- FMRIB Centre, Wellcome Centre for Integrative Neuroimaging, Nuffield Department of Clinical Neurosciences, University of Oxford, Oxford, United Kingdom
| | - Jeroen Mollink
- FMRIB Centre, Wellcome Centre for Integrative Neuroimaging, Nuffield Department of Clinical Neurosciences, University of Oxford, Oxford, United Kingdom
| | - Suryanarayana Umesh Rudrapatna
- Cardiff University Brain Research Imaging Centre, Cardiff University, Cardiff, United Kingdom; Philips Innovation Campus, Bangalore, India
| | - Qiyuan Tian
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, Massachusetts, United States; Harvard Medical School, Boston, Massachusetts, United States
| | - Karla L Miller
- FMRIB Centre, Wellcome Centre for Integrative Neuroimaging, Nuffield Department of Clinical Neurosciences, University of Oxford, Oxford, United Kingdom
| | - Saad Jbabdi
- FMRIB Centre, Wellcome Centre for Integrative Neuroimaging, Nuffield Department of Clinical Neurosciences, University of Oxford, Oxford, United Kingdom
| |
Collapse
|
12
|
Jensen JH. Impact of intra-axonal kurtosis on fiber orientation density functions estimated with fiber ball imaging. Magn Reson Med 2022; 88:1347-1354. [PMID: 35436362 PMCID: PMC9246967 DOI: 10.1002/mrm.29270] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2021] [Revised: 02/27/2022] [Accepted: 03/27/2022] [Indexed: 11/09/2022]
Abstract
PURPOSE To determine the impact of an intra-axonal kurtosis on estimates of the fiber orientation density function (fODF) obtained with fiber ball imaging (FBI). THEORY AND METHODS Standard FBI assumes Gaussian diffusion within individual axons and estimates the fODF by applying an inverse generalized Funk transform to diffusion MRI data for b-values of 4000 s/mm2 or higher. However, recent work based on numeric simulations shows that diffusion inside axons is non-Gaussian with an intra-axonal kurtosis of ∼ 0.4. Here, the theory underlying FBI is extended to incorporate an intra-axonal kurtosis. This is done to first order in the intra-axonal kurtosis without making assumptions about the details of the diffusion dynamics and to all orders for a specific model based on a gamma distribution of diffusivities. The first order approximation is used to assess the effect of an intra-axonal kurtosis on FBI estimates for the fODF and axonal water fraction. The gamma distribution model is used to test the validity of the approximation. RESULTS The first order approximation indicates the estimated fODF is altered by a few percent for an intra-axonal kurtosis of 0.4 in comparison to predictions of standard FBI. If one neglects the intra-axonal kurtosis, the angular resolution of the point spread function for the fODF is changed by <1°, whereas the axonal water fraction is overestimated by ∼ 5%. The gamma distribution model shows that the first order approximation is accurate to within a few percent. CONCLUSION The intra-axonal kurtosis has a small impact on fODFs estimated with FBI.
Collapse
Affiliation(s)
- Jens H. Jensen
- Center for Biomedical Imaging, Medical University of South Carolina, Charleston, South Carolina
- Department of Neuroscience, Medical University of South Carolina, Charleston, South Carolina
- Department of Radiology and Radiological Science, Medical University of South Carolina, Charleston, South Carolina
| |
Collapse
|
13
|
Jelescu IO, de Skowronski A, Geffroy F, Palombo M, Novikov DS. Neurite Exchange Imaging (NEXI): A minimal model of diffusion in gray matter with inter-compartment water exchange. Neuroimage 2022; 256:119277. [PMID: 35523369 PMCID: PMC10363376 DOI: 10.1016/j.neuroimage.2022.119277] [Citation(s) in RCA: 46] [Impact Index Per Article: 23.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2022] [Revised: 04/26/2022] [Accepted: 05/01/2022] [Indexed: 01/18/2023] Open
Abstract
Biophysical models of diffusion in white matter have been center-stage over the past two decades and are essentially based on what is now commonly referred to as the "Standard Model" (SM) of non-exchanging anisotropic compartments with Gaussian diffusion. In this work, we focus on diffusion MRI in gray matter, which requires rethinking basic microstructure modeling blocks. In particular, at least three contributions beyond the SM need to be considered for gray matter: water exchange across the cell membrane - between neurites and the extracellular space; non-Gaussian diffusion along neuronal and glial processes - resulting from structural disorder; and signal contribution from soma. For the first contribution, we propose Neurite Exchange Imaging (NEXI) as an extension of the SM of diffusion, which builds on the anisotropic Kärger model of two exchanging compartments. Using datasets acquired at multiple diffusion weightings (b) and diffusion times (t) in the rat brain in vivo, we investigate the suitability of NEXI to describe the diffusion signal in the gray matter, compared to the other two possible contributions. Our results for the diffusion time window 20-45 ms show minimal diffusivity time-dependence and more pronounced kurtosis decay with time, which is well fit by the exchange model. Moreover, we observe lower signal for longer diffusion times at high b. In light of these observations, we identify exchange as the mechanism that best explains these signal signatures in both low-b and high-b regime, and thereby propose NEXI as the minimal model for gray matter microstructure mapping. We finally highlight multi-b multi-t acquisition protocols as being best suited to estimate NEXI model parameters reliably. Using this approach, we estimate the inter-compartment water exchange time to be 15 - 60 ms in the rat cortex and hippocampus in vivo, which is of the same order or shorter than the diffusion time in typical diffusion MRI acquisitions. This suggests water exchange as an essential component for interpreting diffusion MRI measurements in gray matter.
Collapse
Affiliation(s)
- Ileana O Jelescu
- CIBM Center for Biomedical Imaging, Animal Imaging and Technology, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland; Department of Radiology, Lausanne University Hospital (CHUV), Lausanne, Switzerland; School of Biology and Medicine, University of Lausanne (UNIL), Lausanne, Switzerland.
| | - Alexandre de Skowronski
- CIBM Center for Biomedical Imaging, Animal Imaging and Technology, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
| | | | - Marco Palombo
- School of Psychology, Cardiff University Brain Research Imaging Centre (CUBRIC), Cardiff University, Cardiff, UK; School of Computer Science and Informatics, Cardiff University, Cardiff, UK; Department of Computer Science, Centre for Medical Image Computing, University College London, London, UK
| | - Dmitry S Novikov
- Department of Radiology, Center for Biomedical Imaging, New York University School of Medicine, New York, NY, USA
| |
Collapse
|
14
|
Yendiki A, Aggarwal M, Axer M, Howard AF, van Cappellen van Walsum AM, Haber SN. Post mortem mapping of connectional anatomy for the validation of diffusion MRI. Neuroimage 2022; 256:119146. [PMID: 35346838 PMCID: PMC9832921 DOI: 10.1016/j.neuroimage.2022.119146] [Citation(s) in RCA: 43] [Impact Index Per Article: 21.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2021] [Revised: 03/02/2022] [Accepted: 03/23/2022] [Indexed: 01/13/2023] Open
Abstract
Diffusion MRI (dMRI) is a unique tool for the study of brain circuitry, as it allows us to image both the macroscopic trajectories and the microstructural properties of axon bundles in vivo. The Human Connectome Project ushered in an era of impressive advances in dMRI acquisition and analysis. As a result of these efforts, the quality of dMRI data that could be acquired in vivo improved substantially, and large collections of such data became widely available. Despite this progress, the main limitation of dMRI remains: it does not image axons directly, but only provides indirect measurements based on the diffusion of water molecules. Thus, it must be validated by methods that allow direct visualization of axons but that can only be performed in post mortem brain tissue. In this review, we discuss methods for validating the various features of connectional anatomy that are extracted from dMRI, both at the macro-scale (trajectories of axon bundles), and at micro-scale (axonal orientations and other microstructural properties). We present a range of validation tools, including anatomic tracer studies, Klingler's dissection, myelin stains, label-free optical imaging techniques, and others. We provide an overview of the basic principles of each technique, its limitations, and what it has taught us so far about the accuracy of different dMRI acquisition and analysis approaches.
Collapse
Affiliation(s)
- Anastasia Yendiki
- Department of Radiology, Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital and Harvard Medical School, Charlestown, MA, United States,Corresponding author (A. Yendiki)
| | - Manisha Aggarwal
- Department of Radiology and Radiological Science, Johns Hopkins University School of Medicine, Baltimore, MD, United States
| | - Markus Axer
- Forschungszentrum Jülich, Institute of Neuroscience and Medicine, Jülich, Germany,Department of Physics, University of Wuppertal Germany
| | - Amy F.D. Howard
- Wellcome Centre for Integrative Neuroimaging, Nuffield Department of Clinical Neurosciences, University of Oxford, Oxford, United Kingdom
| | - Anne-Marie van Cappellen van Walsum
- Department of Medical Imaging, Anatomy, Radboud University Medical Center, Nijmegen, the Netherland,Cognition and Behaviour, Donders Institute for Brain, Nijmegen, the Netherland
| | - Suzanne N. Haber
- Department of Pharmacology and Physiology, University of Rochester, Rochester, NY, United States,McLean Hospital, Belmont, MA, United States
| |
Collapse
|
15
|
Fan Q, Eichner C, Afzali M, Mueller L, Tax CMW, Davids M, Mahmutovic M, Keil B, Bilgic B, Setsompop K, Lee HH, Tian Q, Maffei C, Ramos-Llordén G, Nummenmaa A, Witzel T, Yendiki A, Song YQ, Huang CC, Lin CP, Weiskopf N, Anwander A, Jones DK, Rosen BR, Wald LL, Huang SY. Mapping the human connectome using diffusion MRI at 300 mT/m gradient strength: Methodological advances and scientific impact. Neuroimage 2022; 254:118958. [PMID: 35217204 PMCID: PMC9121330 DOI: 10.1016/j.neuroimage.2022.118958] [Citation(s) in RCA: 22] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2021] [Revised: 01/27/2022] [Accepted: 01/31/2022] [Indexed: 12/20/2022] Open
Abstract
Tremendous efforts have been made in the last decade to advance cutting-edge MRI technology in pursuit of mapping structural connectivity in the living human brain with unprecedented sensitivity and speed. The first Connectom 3T MRI scanner equipped with a 300 mT/m whole-body gradient system was installed at the Massachusetts General Hospital in 2011 and was specifically constructed as part of the Human Connectome Project. Since that time, numerous technological advances have been made to enable the broader use of the Connectom high gradient system for diffusion tractography and tissue microstructure studies and leverage its unique advantages and sensitivity to resolving macroscopic and microscopic structural information in neural tissue for clinical and neuroscientific studies. The goal of this review article is to summarize the technical developments that have emerged in the last decade to support and promote large-scale and scientific studies of the human brain using the Connectom scanner. We provide a brief historical perspective on the development of Connectom gradient technology and the efforts that led to the installation of three other Connectom 3T MRI scanners worldwide - one in the United Kingdom in Cardiff, Wales, another in continental Europe in Leipzig, Germany, and the latest in Asia in Shanghai, China. We summarize the key developments in gradient hardware and image acquisition technology that have formed the backbone of Connectom-related research efforts, including the rich array of high-sensitivity receiver coils, pulse sequences, image artifact correction strategies and data preprocessing methods needed to optimize the quality of high-gradient strength diffusion MRI data for subsequent analyses. Finally, we review the scientific impact of the Connectom MRI scanner, including advances in diffusion tractography, tissue microstructural imaging, ex vivo validation, and clinical investigations that have been enabled by Connectom technology. We conclude with brief insights into the unique value of strong gradients for diffusion MRI and where the field is headed in the coming years.
Collapse
Affiliation(s)
- Qiuyun Fan
- Department of Biomedical Engineering, College of Precision Instruments and Optoelectronics Engineering, Tianjin University, Tianjin, China; Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, USA; Harvard Medical School, Boston, MA, USA
| | - Cornelius Eichner
- Max Planck Institute for Human Cognitive and Brain Sciences, Department of Neuropsychology, Leipzig, Germany
| | - Maryam Afzali
- Cardiff University Brain Research Imaging Centre (CUBRIC), Cardiff University, Cardiff, Wales, UK; Leeds Institute of Cardiovascular and Metabolic Medicine, University of Leeds, Leeds LS2 9JT, UK
| | - Lars Mueller
- Leeds Institute of Cardiovascular and Metabolic Medicine, University of Leeds, Leeds LS2 9JT, UK
| | - Chantal M W Tax
- Cardiff University Brain Research Imaging Centre (CUBRIC), Cardiff University, Cardiff, Wales, UK; Image Sciences Institute, University Medical Center (UMC) Utrecht, Utrecht, the Netherlands
| | - Mathias Davids
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, USA; Harvard Medical School, Boston, MA, USA; Computer Assisted Clinical Medicine, Medical Faculty Mannheim, Heidelberg University, Mannheim, Germany
| | - Mirsad Mahmutovic
- Institute of Medical Physics and Radiation Protection (IMPS), TH-Mittelhessen University of Applied Sciences (THM), Giessen, Germany
| | - Boris Keil
- Institute of Medical Physics and Radiation Protection (IMPS), TH-Mittelhessen University of Applied Sciences (THM), Giessen, Germany
| | - Berkin Bilgic
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, USA; Harvard Medical School, Boston, MA, USA; Harvard-MIT Division of Health Sciences and Technology, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Kawin Setsompop
- Department of Radiology, Stanford University, Stanford, CA, USA; Department of Electrical Engineering, Stanford University, Stanford, CA, USA
| | - Hong-Hsi Lee
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, USA; Harvard Medical School, Boston, MA, USA
| | - Qiyuan Tian
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, USA; Harvard Medical School, Boston, MA, USA
| | - Chiara Maffei
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, USA; Harvard Medical School, Boston, MA, USA
| | - Gabriel Ramos-Llordén
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, USA; Harvard Medical School, Boston, MA, USA
| | - Aapo Nummenmaa
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, USA; Harvard Medical School, Boston, MA, USA
| | | | - Anastasia Yendiki
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, USA; Harvard Medical School, Boston, MA, USA
| | - Yi-Qiao Song
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, USA; John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA USA
| | - Chu-Chung Huang
- Key Laboratory of Brain Functional Genomics (MOE & STCSM), Affiliated Mental Health Center (ECNU), School of Psychology and Cognitive Science, East China Normal University, Shanghai, China; Shanghai Changning Mental Health Center, Shanghai, China
| | - Ching-Po Lin
- Institute of Neuroscience, National Yang Ming Chiao Tung University, Taipei, Taiwan; Institute of Science and Technology for Brain-Inspired Intelligence, Fudan University, Shanghai, China
| | - Nikolaus Weiskopf
- Department of Neurophysics, Max Planck Institute for Human Cognitive and Brain Sciences, Leipzig, Germany; Felix Bloch Institute for Solid State Physics, Faculty of Physics and Earth Sciences, Leipzig University, Leipzig, Germany
| | - Alfred Anwander
- Max Planck Institute for Human Cognitive and Brain Sciences, Department of Neuropsychology, Leipzig, Germany
| | - Derek K Jones
- Cardiff University Brain Research Imaging Centre (CUBRIC), Cardiff University, Cardiff, Wales, UK
| | - Bruce R Rosen
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, USA; Harvard Medical School, Boston, MA, USA; Harvard-MIT Division of Health Sciences and Technology, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Lawrence L Wald
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, USA; Harvard Medical School, Boston, MA, USA; Harvard-MIT Division of Health Sciences and Technology, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Susie Y Huang
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, USA; Harvard Medical School, Boston, MA, USA; Harvard-MIT Division of Health Sciences and Technology, Massachusetts Institute of Technology, Cambridge, MA, USA.
| |
Collapse
|
16
|
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] [Key Words] [MESH Headings] [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
Affiliation(s)
- Jonas L Olesen
- 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
| | - Leif Østergaard
- Center of Functionally Integrative Neuroscience (CFIN) and MINDLab, Department of Clinical Medicine, Aarhus University, Aarhus, Denmark
| | - Noam Shemesh
- Champalimaud Research, Champalimaud Centre for the Unknown, Lisbon, Portugal
| | - 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.
| |
Collapse
|
17
|
Spees WM, Sukstanskii AL, Bretthorst GL, Neil JJ, Ackerman JJH. Rat Brain Global Ischemia-Induced Diffusion Changes Revisited: Biophysical Modeling of the Water and NAA MR "Diffusion Signal". Magn Reson Med 2022; 88:1333-1346. [PMID: 35452137 DOI: 10.1002/mrm.29262] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2021] [Revised: 03/18/2022] [Accepted: 03/20/2022] [Indexed: 11/06/2022]
Abstract
PURPOSE To assess changes in intracellular diffusion as a mechanism for the reduction in water ADC that accompanies brain injury. Using NAA as a marker of neuronal cytoplasmic diffusion, NAA diffusion was measured before and after global ischemia (immediately postmortem) in the female Sprague-Dawley rat. METHODS Diffusion-weighted PRESS spectra, with diffusion encoding in a single direction, were acquired from large voxels of rat brain gray matter in vivo and postischemia employing either pairs of pulsed half-sine-shaped gradients (in vivo and postischemia, bmax = 19 ms/μm2 ) or sinusoidal oscillating gradients (in vivo only) with frequencies of 99.2-250 Hz. A 2D randomly oriented cylinder (neurite) model gave estimates of longitudinal and transverse diffusivities (DL and DT , respectively). In this model, DL represents the "free" diffusivity of NAA, whereas DT reflects highly restricted diffusion. Using oscillating gradients, the frequency dependence of DT [DT (ω)] gave estimates of the cylinder (axon/dendrite) radius. RESULTS A 10% decrease in DL,NAA followed global ischemia, dropping from 0.391 ± 0.012 μm2 /ms to 0.350 ± 0.009 μm2 /ms. Modeling DT,NAA (ω) provided an estimate of the neurite radius of 1.0 ± 0.6 μm. CONCLUSION Whereas the increase in apparent intraneuronal viscosity suggested by changes in DL,NAA may contribute to the overall reduction in water ADC associated with brain injury, it is not sufficient to be the sole explanation. Estimates of neurite radius based on DT (ω) were consistent with literature values.
Collapse
Affiliation(s)
- William M Spees
- Biomedical MR Laboratory, Mallinckrodt Institute of Radiology, Washington University School of Medicine, St. Louis, Missouri
| | - Alex L Sukstanskii
- Biomedical MR Laboratory, Mallinckrodt Institute of Radiology, Washington University School of Medicine, St. Louis, Missouri
| | - G Larry Bretthorst
- Biomedical MR Laboratory, Mallinckrodt Institute of Radiology, Washington University School of Medicine, St. Louis, Missouri
| | - Jeffrey J Neil
- Biomedical MR Laboratory, Mallinckrodt Institute of Radiology, Washington University School of Medicine, St. Louis, Missouri.,Department of Neurology, Washington University School of Medicine, St. Louis, Missouri.,Department of Pediatrics, Washington University School of Medicine, St. Louis, Missouri
| | - Joseph J H Ackerman
- Biomedical MR Laboratory, Mallinckrodt Institute of Radiology, Washington University School of Medicine, St. Louis, Missouri.,Department of Chemistry, Washington University, St. Louis, Missouri.,Department of Internal Medicine, Washington University School of Medicine, St. Louis, Missouri.,Alvin J Siteman Cancer Center, Washington University School of Medicine, St. Louis, Missouri
| |
Collapse
|
18
|
Yang Q, Reutens DC, Vegh V. Generalisation of continuous time random walk to anomalous diffusion MRI models with an age-related evaluation of human corpus callosum. Neuroimage 2022; 250:118903. [PMID: 35033674 DOI: 10.1016/j.neuroimage.2022.118903] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2021] [Revised: 12/07/2021] [Accepted: 01/10/2022] [Indexed: 12/22/2022] Open
Abstract
Diffusion MRI measures of the human brain provide key insight into microstructural variations across individuals and into the impact of central nervous system diseases and disorders. One approach to extract information from diffusion signals has been to use biologically relevant analytical models to link millimetre scale diffusion MRI measures with microscale influences. The other approach has been to represent diffusion as an anomalous transport process and infer microstructural information from the different anomalous diffusion equation parameters. In this study, we investigated how parameters of various anomalous diffusion models vary with age in the human brain white matter, particularly focusing on the corpus callosum. We first unified several established anomalous diffusion models (the super-diffusion, sub-diffusion, quasi-diffusion and fractional Bloch-Torrey models) under the continuous time random walk modelling framework. This unification allows a consistent parameter fitting strategy to be applied from which meaningful model parameter comparisons can be made. We then provided a novel way to derive the diffusional kurtosis imaging (DKI) model, which is shown to be a degree two approximation of the sub-diffusion model. This link between the DKI and sub-diffusion models led to a new robust technique for generating maps of kurtosis and diffusivity using the sub-diffusion parameters βSUB and DSUB. Superior tissue contrast is achieved in kurtosis maps based on the sub-diffusion model. 7T diffusion weighted MRI data for 65 healthy participants in the age range 19-78 years was used in this study. Results revealed that anomalous diffusion model parameters α and β have shown consistent positive correlation with age in the corpus callosum, indicating α and β are sensitive to tissue microstructural changes in ageing.
Collapse
Affiliation(s)
- Qianqian Yang
- School of Mathematical Sciences, Faculty of Science, Queensland University of Technology, Brisbane 4000, Australia.
| | - David C Reutens
- Centre for Advanced Imaging, University of Queensland, Brisbane 4072, Australia; ARC Training Centre for Innovation in Biomedical Imaging Technology, Brisbane 4072, Australia
| | - Viktor Vegh
- Centre for Advanced Imaging, University of Queensland, Brisbane 4072, Australia; ARC Training Centre for Innovation in Biomedical Imaging Technology, Brisbane 4072, Australia
| |
Collapse
|
19
|
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
|
20
|
Connecting macroscopic diffusion metrics of cardiac diffusion tensor imaging and microscopic myocardial structures based on simulation. Med Image Anal 2022; 77:102325. [DOI: 10.1016/j.media.2021.102325] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2019] [Revised: 11/26/2021] [Accepted: 11/29/2021] [Indexed: 11/20/2022]
|
21
|
Li X, Li M, Wang M, Wu F, Liu H, Sun Q, Zhang Y, Liu C, Jin C, Yang J. Mapping white matter maturational processes and degrees on neonates by diffusion kurtosis imaging with multiparametric analysis. Hum Brain Mapp 2022; 43:799-815. [PMID: 34708903 PMCID: PMC8720196 DOI: 10.1002/hbm.25689] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2021] [Revised: 10/03/2021] [Accepted: 10/07/2021] [Indexed: 11/10/2022] Open
Abstract
White matter maturation has been characterized by diffusion tensor (DT) metrics. However, maturational processes and degrees are not fully investigated due to limitations of univariate approaches and limited specificity/sensitivity. Diffusion kurtosis imaging (DKI) provides kurtosis tensor (KT) and white matter tract integrity (WMTI) metrics, besides DT metrics. Therefore, we tried to investigate performances of DKI with the multiparametric analysis in characterizing white matter maturation. Developmental changes in metrics were investigated by using tract-based spatial statistics and the region of interest analysis on 50 neonates with postmenstrual age (PMA) from 37.43 to 43.57 weeks. Changes in metrics were combined into various patterns to reveal different maturational processes. Mahalanobis distance based on DT metrics (DM,DT ) and that combing DT and KT metrics (DM,DT-KT ) were computed, separately. Performances of DM,DT-KT and DM,DT were compared in revealing correlations with PMA and the neurobehavioral score. Compared with DT metrics, WMTI metrics demonstrated additional changing patterns. Furthermore, variations of DM,DT-KT across regions were in agreement with the maturational sequence. Additionally, DM,DT-KT demonstrated stronger negative correlations with PMA and the neurobehavioral score in more regions than DM,DT . Results suggest that DKI with the multiparametric analysis benefits the understanding of white matter maturational processes and degrees on neonates.
Collapse
Affiliation(s)
- Xianjun Li
- Department of Radiology, The First Affiliated Hospital of Xi'an Jiaotong University, Xi'an, China
| | - Mengxuan Li
- Department of Radiology, The First Affiliated Hospital of Xi'an Jiaotong University, Xi'an, China
| | - Miaomiao Wang
- Department of Radiology, The First Affiliated Hospital of Xi'an Jiaotong University, Xi'an, China
| | - Fan Wu
- Department of Radiology, The First Affiliated Hospital of Xi'an Jiaotong University, Xi'an, China
| | - Heng Liu
- Department of Radiology, The First Affiliated Hospital of Xi'an Jiaotong University, Xi'an, China.,Department of Biomedical Engineering, The Key Laboratory of Biomedical Information Engineering of the Ministry of Education, School of Life Science and Technology, Xi'an Jiaotong University, Xi'an, China
| | - Qinli Sun
- Department of Radiology, The First Affiliated Hospital of Xi'an Jiaotong University, Xi'an, China.,Department of Biomedical Engineering, The Key Laboratory of Biomedical Information Engineering of the Ministry of Education, School of Life Science and Technology, Xi'an Jiaotong University, Xi'an, China
| | - Yuli Zhang
- Department of Radiology, The First Affiliated Hospital of Xi'an Jiaotong University, Xi'an, China
| | - Congcong Liu
- Department of Radiology, The First Affiliated Hospital of Xi'an Jiaotong University, Xi'an, China
| | - Chao Jin
- Department of Radiology, The First Affiliated Hospital of Xi'an Jiaotong University, Xi'an, China
| | - Jian Yang
- Department of Radiology, The First Affiliated Hospital of Xi'an Jiaotong University, Xi'an, China.,Department of Biomedical Engineering, The Key Laboratory of Biomedical Information Engineering of the Ministry of Education, School of Life Science and Technology, Xi'an Jiaotong University, Xi'an, China
| |
Collapse
|
22
|
Wang W, Metzler R, Cherstvy AG. Anomalous diffusion, aging, and nonergodicity of scaled Brownian motion with fractional Gaussian noise: overview of related experimental observations and models. Phys Chem Chem Phys 2022; 24:18482-18504. [DOI: 10.1039/d2cp01741e] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
How does a systematic time-dependence of the diffusion coefficient $D (t)$ affect the ergodic and statistical characteristics of fractional Brownian motion (FBM)? Here, we examine how the behavior of the...
Collapse
|
23
|
Correlation Tensor MRI deciphers underlying kurtosis sources in stroke. Neuroimage 2021; 247:118833. [PMID: 34929382 DOI: 10.1016/j.neuroimage.2021.118833] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2021] [Revised: 12/15/2021] [Accepted: 12/16/2021] [Indexed: 02/06/2023] Open
Abstract
Noninvasively detecting and characterizing modulations in cellular scale micro-architecture remains a desideratum for contemporary neuroimaging. Diffusion MRI (dMRI) has become the mainstay methodology for probing microstructure, and, in ischemia, its contrasts have revolutionized stroke management. Diffusion kurtosis imaging (DKI) has been shown to significantly enhance the sensitivity of stroke detection compared to its diffusion tensor imaging (DTI) counterparts. However, the interpretation of DKI remains ambiguous as its contrast may arise from competing kurtosis sources related to the anisotropy of tissue components, diffusivity variance across components, and microscopic kurtosis (e.g., arising from cross-sectional variance, structural disorder, and restriction). Resolving these sources may be fundamental for developing more specific imaging techniques for stroke management, prognosis, and understanding its pathophysiology. In this study, we apply Correlation Tensor MRI (CTI) - a double diffusion encoding (DDE) methodology recently introduced for deciphering kurtosis sources based on the unique information captured in DDE's diffusion correlation tensors - to investigate the underpinnings of kurtosis measurements in acute ischemic lesions. Simulations for the different kurtosis sources revealed specific signatures for cross-sectional variance (representing neurite beading), edema, and cell swelling. Ex vivo CTI experiments at 16.4 T were then performed in an experimental photothrombotic stroke model 3 h post-stroke (N = 10), and successfully separated anisotropic, isotropic, and microscopic non-Gaussian diffusion sources in the ischemic lesions. Each of these kurtosis sources provided unique contrasts in the stroked area. Particularly, microscopic kurtosis was shown to be a primary "driver" of total kurtosis upon ischemia; its large increases, coupled with decreases in anisotropic kurtosis, are consistent with the expected elevation in cross-sectional variance, likely linked to beading effects in small objects such as neurites. In vivo experiments at 9.4 T at the same time point (3 h post ischemia, N = 5) demonstrated the stability and relevance of the findings and showed that fixation is not a dominant confounder in our findings. In future studies, the different CTI contrasts may be useful to address current limitations of stroke imaging, e.g., penumbra characterization, distinguishing lesion progression form tissue recovery, and elucidating pathophysiological correlates.
Collapse
|
24
|
Andersson M, Pizzolato M, Kjer HM, Skodborg KF, Lundell H, Dyrby TB. Does powder averaging remove dispersion bias in diffusion MRI diameter estimates within real 3D axonal architectures? Neuroimage 2021; 248:118718. [PMID: 34767939 DOI: 10.1016/j.neuroimage.2021.118718] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2021] [Revised: 10/26/2021] [Accepted: 11/08/2021] [Indexed: 11/26/2022] Open
Abstract
Noninvasive estimation of axon diameter with diffusion MRI holds the potential to investigate the dynamic properties of the brain network and pathology of neurodegenerative diseases. Recent studies use powder averaging to account for complex white matter architectures, but these have not been validated for real axonal geometries from regions that contain fibre crossings. Here, we present 120-304μm long segmented axons from X-ray nano-holotomography volumes of a splenium and crossing fibre region of a vervet monkey brain. We show that the axons in the complex crossing fibre region, which contains callosal, association, and corticospinal connections, are larger and exhibit a wider distribution than those of the splenium region. To accurately estimate the axon diameter in these regions, therefore, sensitivity to a wide range of diameters is required. We demonstrate how the q-value, b-value, signal-to-noise ratio and the assumed intra-axonal parallel diffusivity influence the range of measurable diameters with powder average approaches. Furthermore, we show how Gaussian distributed noise results in a wider range of measurable diameter at high b-values than Rician distributed noise, even at high signal-to-noise ratios of 100. The number of gradient directions is also shown to impose a lower bound on measurable diameter. Our results indicate that axon diameter estimation can be performed with only few b-shells, and that additional shells do not improve the accuracy of the estimate. For strong gradients available on human Connectom and preclinical scanners, Monte Carlo simulations of diffusion confirm that powder averaging techniques succeed in providing accurate estimates of axon diameter across a range of sequence parameters and diffusion times, even in complex white matter architectures. At relatively low b-values, the diameter estimate becomes sensitive to axonal microdispersion and the intra-axonal parallel diffusivity shows time dependency at both in vivo and ex vivo intrinsic diffusivities.
Collapse
Affiliation(s)
- Mariam Andersson
- Danish Research Centre for Magnetic Resonance, Center for Functional and Diagnostic Imaging and Research, Copenhagen University Hospital Amager and Hvidovre, Hvidovre 2650, Denmark; Department of Applied Mathematics and Computer Science, Technical University of Denmark, Kongens Lyngby 2800, Denmark
| | - Marco Pizzolato
- Danish Research Centre for Magnetic Resonance, Center for Functional and Diagnostic Imaging and Research, Copenhagen University Hospital Amager and Hvidovre, Hvidovre 2650, Denmark; Department of Applied Mathematics and Computer Science, Technical University of Denmark, Kongens Lyngby 2800, Denmark; Signal Processing Laboratory (LTS5), École Polytechnique Fédérale de Lausanne, Lausanne 1015, Switzerland
| | - Hans Martin Kjer
- Danish Research Centre for Magnetic Resonance, Center for Functional and Diagnostic Imaging and Research, Copenhagen University Hospital Amager and Hvidovre, Hvidovre 2650, Denmark; Department of Applied Mathematics and Computer Science, Technical University of Denmark, Kongens Lyngby 2800, Denmark
| | - Katrine Forum Skodborg
- Danish Research Centre for Magnetic Resonance, Center for Functional and Diagnostic Imaging and Research, Copenhagen University Hospital Amager and Hvidovre, Hvidovre 2650, Denmark; Department of Applied Mathematics and Computer Science, Technical University of Denmark, Kongens Lyngby 2800, Denmark
| | - Henrik Lundell
- Danish Research Centre for Magnetic Resonance, Center for Functional and Diagnostic Imaging and Research, Copenhagen University Hospital Amager and Hvidovre, Hvidovre 2650, Denmark; Department of Applied Mathematics and Computer Science, Technical University of Denmark, Kongens Lyngby 2800, Denmark
| | - Tim B Dyrby
- Danish Research Centre for Magnetic Resonance, Center for Functional and Diagnostic Imaging and Research, Copenhagen University Hospital Amager and Hvidovre, Hvidovre 2650, Denmark; Department of Applied Mathematics and Computer Science, Technical University of Denmark, Kongens Lyngby 2800, Denmark
| |
Collapse
|
25
|
Huang SY, Witzel T, Keil B, Scholz A, Davids M, Dietz P, Rummert E, Ramb R, Kirsch JE, Yendiki A, Fan Q, Tian Q, Ramos-Llordén G, Lee HH, Nummenmaa A, Bilgic B, Setsompop K, Wang F, Avram AV, Komlosh M, Benjamini D, Magdoom KN, Pathak S, Schneider W, Novikov DS, Fieremans E, Tounekti S, Mekkaoui C, Augustinack J, Berger D, Shapson-Coe A, Lichtman J, Basser PJ, Wald LL, Rosen BR. Connectome 2.0: Developing the next-generation ultra-high gradient strength human MRI scanner for bridging studies of the micro-, meso- and macro-connectome. Neuroimage 2021; 243:118530. [PMID: 34464739 PMCID: PMC8863543 DOI: 10.1016/j.neuroimage.2021.118530] [Citation(s) in RCA: 55] [Impact Index Per Article: 18.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2021] [Revised: 08/10/2021] [Accepted: 08/27/2021] [Indexed: 11/26/2022] Open
Abstract
The first phase of the Human Connectome Project pioneered advances in MRI technology for mapping the macroscopic structural connections of the living human brain through the engineering of a whole-body human MRI scanner equipped with maximum gradient strength of 300 mT/m, the highest ever achieved for human imaging. While this instrument has made important contributions to the understanding of macroscale connectional topology, it has also demonstrated the potential of dedicated high-gradient performance scanners to provide unparalleled in vivo assessment of neural tissue microstructure. Building on the initial groundwork laid by the original Connectome scanner, we have now embarked on an international, multi-site effort to build the next-generation human 3T Connectome scanner (Connectome 2.0) optimized for the study of neural tissue microstructure and connectional anatomy across multiple length scales. In order to maximize the resolution of this in vivo microscope for studies of the living human brain, we will push the diffusion resolution limit to unprecedented levels by (1) nearly doubling the current maximum gradient strength from 300 mT/m to 500 mT/m and tripling the maximum slew rate from 200 T/m/s to 600 T/m/s through the design of a one-of-a-kind head gradient coil optimized to minimize peripheral nerve stimulation; (2) developing high-sensitivity multi-channel radiofrequency receive coils for in vivo and ex vivo human brain imaging; (3) incorporating dynamic field monitoring to minimize image distortions and artifacts; (4) developing new pulse sequences to integrate the strongest diffusion encoding and highest spatial resolution ever achieved in the living human brain; and (5) calibrating the measurements obtained from this next-generation instrument through systematic validation of diffusion microstructural metrics in high-fidelity phantoms and ex vivo brain tissue at progressively finer scales with accompanying diffusion simulations in histology-based micro-geometries. We envision creating the ultimate diffusion MRI instrument capable of capturing the complex multi-scale organization of the living human brain - from the microscopic scale needed to probe cellular geometry, heterogeneity and plasticity, to the mesoscopic scale for quantifying the distinctions in cortical structure and connectivity that define cyto- and myeloarchitectonic boundaries, to improvements in estimates of macroscopic connectivity.
Collapse
Affiliation(s)
- Susie Y Huang
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Harvard Medical School, Boston, MA, USA.
| | | | - Boris Keil
- Institute of Medical Physics and Radiation Protection (IMPS), TH-Mittelhessen University of Applied Sciences (THM), Giessen, Germany
| | - Alina Scholz
- Institute of Medical Physics and Radiation Protection (IMPS), TH-Mittelhessen University of Applied Sciences (THM), Giessen, Germany
| | - Mathias Davids
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Harvard Medical School, Boston, MA, USA
| | | | | | | | - John E Kirsch
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Harvard Medical School, Boston, MA, USA
| | - Anastasia Yendiki
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Harvard Medical School, Boston, MA, USA
| | - Qiuyun Fan
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Harvard Medical School, Boston, MA, USA
| | - Qiyuan Tian
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Harvard Medical School, Boston, MA, USA
| | - Gabriel Ramos-Llordén
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Harvard Medical School, Boston, MA, USA
| | - Hong-Hsi Lee
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Harvard Medical School, Boston, MA, USA
| | - Aapo Nummenmaa
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Harvard Medical School, Boston, MA, USA
| | - Berkin Bilgic
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Harvard Medical School, Boston, MA, USA
| | - Kawin Setsompop
- Radiological Sciences Laboratory, Department of Radiology, Stanford University, Stanford, CA, USA
| | - Fuyixue Wang
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Harvard Medical School, Boston, MA, USA
| | - Alexandru V Avram
- Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health, Bethesda, MD, USA
| | - Michal Komlosh
- 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
| | - Kulam Najmudeen Magdoom
- Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health, Bethesda, MD, USA
| | - Sudhir Pathak
- Learning Research and Development Center, University of Pittsburgh, Pittsburgh, PA, USA
| | - Walter Schneider
- Learning Research and Development Center, University of Pittsburgh, Pittsburgh, PA, 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
| | - Slimane Tounekti
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Harvard Medical School, Boston, MA, USA
| | - Choukri Mekkaoui
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Harvard Medical School, Boston, MA, USA
| | - Jean Augustinack
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Harvard Medical School, Boston, MA, USA
| | - Daniel Berger
- Department of Molecular and Cell Biology and Center for Brain Science, Harvard University, Cambridge, MA, USA
| | - Alexander Shapson-Coe
- Department of Molecular and Cell Biology and Center for Brain Science, Harvard University, Cambridge, MA, USA
| | - Jeff Lichtman
- Department of Molecular and Cell Biology and Center for Brain Science, Harvard University, Cambridge, MA, USA
| | - Peter J Basser
- Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health, Bethesda, MD, USA
| | - Lawrence L Wald
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Harvard Medical School, Boston, MA, USA
| | - Bruce R Rosen
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Harvard Medical School, Boston, MA, USA
| |
Collapse
|
26
|
Huang CC, Hsu CCH, Zhou FL, Kusmia S, Drakesmith M, Parker GJM, Lin CP, Jones DK. Validating pore size estimates in a complex microfiber environment on a human MRI system. Magn Reson Med 2021; 86:1514-1530. [PMID: 33960501 PMCID: PMC7613441 DOI: 10.1002/mrm.28810] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2020] [Revised: 03/18/2021] [Accepted: 03/26/2021] [Indexed: 12/19/2022]
Abstract
PURPOSE Recent advances in diffusion-weighted MRI provide "restricted diffusion signal fraction" and restricting pore size estimates. Materials based on co-electrospun oriented hollow cylinders have been introduced to provide validation for such methods. This study extends this work, exploring accuracy and repeatability using an extended acquisition on a 300 mT/m gradient human MRI scanner, in substrates closely mimicking tissue, that is, non-circular cross-sections, intra-voxel fiber crossing, intra-voxel distributions of pore-sizes, and smaller pore-sizes overall. METHODS In a single-blind experiment, diffusion-weighted data were collected from a biomimetic phantom on a 3T Connectom system using multiple gradient directions/diffusion times. Repeated scans established short-term and long-term repeatability. The total scan time (54 min) matched similar protocols used in human studies. The number of distinct fiber populations was estimated using spherical deconvolution, and median pore size estimated through the combination of CHARMED and AxCaliber3D framework. Diffusion-based estimates were compared with measurements derived from scanning electron microscopy. RESULTS The phantom contained substrates with different orientations, fiber configurations, and pore size distributions. Irrespective of one or two populations within the voxel, the pore-size estimates (~5 μm) and orientation-estimates showed excellent agreement with the median values of pore-size derived from scanning electron microscope and phantom configuration. Measurement repeatability depended on substrate complexity, with lower values seen in samples containing crossing-fibers. Sample-level repeatability was found to be good. CONCLUSION While no phantom mimics tissue completely, this study takes a step closer to validating diffusion microstructure measurements for use in vivo by demonstrating the ability to quantify microgeometry in relatively complex configurations.
Collapse
Affiliation(s)
- Chu-Chung Huang
- Key Laboratory of Brain Functional Genomics (MOE & STCSM), Affiliated Mental Health Center (ECNU), Institute of Cognitive Neuroscience, School of Psychology and Cognitive Science, East China Normal University, Shanghai, China
- Shanghai Changning Mental Health Center, Shanghai, China
| | - Chih-Chin Heather Hsu
- Center for Geriatrics and Gerontology, Taipei Veterans General Hospital, Taipei, Taiwan
- Institute of Neuroscience, National Yang Ming Chiao Tung University, Taipei, Taiwan
| | - Feng-Lei Zhou
- Centre for Medical Image Computing, Department of Computer Science, University College London, London, United Kingdom
- School of Pharmacy, University College London, London, United Kingdom
| | - Slawomir Kusmia
- Centre for Medical Image Computing, Department of Computer Science, University College London, London, United Kingdom
- Cardiff University Brain Research Imaging Centre (CUBRIC), Cardiff University, Cardiff, United Kingdom
- Epilepsy Society MRI Unit, Chalfont St Peter, United Kingdom
| | - Mark Drakesmith
- Cardiff University Brain Research Imaging Centre (CUBRIC), Cardiff University, Cardiff, United Kingdom
| | - Geoff J. M. Parker
- Centre for Medical Image Computing, Department of Computer Science, University College London, London, United Kingdom
- Department of Neuroinflammation, Queen Square Institute of Neurology, University College London, London, United Kingdom
- Bioxydyn Limited, Manchester, United Kingdom
| | - Ching-Po Lin
- Institute of Neuroscience, National Yang Ming Chiao Tung University, Taipei, Taiwan
| | - Derek K. Jones
- Cardiff University Brain Research Imaging Centre (CUBRIC), Cardiff University, Cardiff, United Kingdom
| |
Collapse
|
27
|
Afzali M, Nilsson M, Palombo M, Jones DK. SPHERIOUSLY? The challenges of estimating sphere radius non-invasively in the human brain from diffusion MRI. Neuroimage 2021; 237:118183. [PMID: 34020013 PMCID: PMC8285594 DOI: 10.1016/j.neuroimage.2021.118183] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/26/2020] [Revised: 04/25/2021] [Accepted: 05/16/2021] [Indexed: 11/16/2022] Open
Abstract
The Soma and Neurite Density Imaging (SANDI) three-compartment model was recently proposed to disentangle cylindrical and spherical geometries, attributed to neurite and soma compartments, respectively, in brain tissue. There are some recent advances in diffusion-weighted MRI signal encoding and analysis (including the use of multiple so-called 'b-tensor' encodings and analysing the signal in the frequency-domain) that have not yet been applied in the context of SANDI. In this work, using: (i) ultra-strong gradients; (ii) a combination of linear, planar, and spherical b-tensor encodings; and (iii) analysing the signal in the frequency domain, three main challenges to robust estimation of sphere size were identified: First, the Rician noise floor in magnitude-reconstructed data biases estimates of sphere properties in a non-uniform fashion. It may cause overestimation or underestimation of the spherical compartment size and density. This can be partly ameliorated by accounting for the noise floor in the estimation routine. Second, even when using the strongest diffusion-encoding gradient strengths available for human MRI, there is an empirical lower bound on the spherical signal fraction and radius that can be detected and estimated robustly. For the experimental setup used here, the lower bound on the sphere signal fraction was approximately 10%. We employed two different ways of establishing the lower bound for spherical radius estimates in white matter. The first, examining power-law relationships between the DW-signal and diffusion weighting in empirical data, yielded a lower bound of 7μm, while the second, pure Monte Carlo simulations, yielded a lower limit of 3μm and in this low radii domain, there is little differentiation in signal attenuation. Third, if there is sensitivity to the transverse intra-cellular diffusivity in cylindrical structures, e.g., axons and cellular projections, then trying to disentangle two diffusion-time-dependencies using one experimental parameter (i.e., change in frequency-content of the encoding waveform) makes spherical radii estimates particularly challenging. We conclude that due to the aforementioned challenges spherical radii estimates may be biased when the corresponding sphere signal fraction is low, which must be considered.
Collapse
Affiliation(s)
- Maryam Afzali
- Cardiff University Brain Research Imaging Centre (CUBRIC), School of Psychology, Cardiff University, Cardiff, United Kingdom.
| | - Markus Nilsson
- Clinical Sciences Lund, Radiology, Lund University, Lund, Sweden.
| | - Marco Palombo
- Centre for Medical Image Computing, Department of Computer Science, University College London, London, United Kingdom.
| | - Derek K Jones
- Cardiff University Brain Research Imaging Centre (CUBRIC), School of Psychology, Cardiff University, Cardiff, United Kingdom.
| |
Collapse
|
28
|
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
|
29
|
Ianus A, Alexander DC, Zhang H, Palombo M. Mapping complex cell morphology in the grey matter with double diffusion encoding MR: A simulation study. Neuroimage 2021; 241:118424. [PMID: 34311067 PMCID: PMC8961003 DOI: 10.1016/j.neuroimage.2021.118424] [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: 01/08/2021] [Revised: 07/13/2021] [Accepted: 07/21/2021] [Indexed: 01/18/2023] Open
Abstract
This paper investigates the impact of cell body (namely soma) size and branching of cellular projections on diffusion MR imaging (dMRI) and spectroscopy (dMRS) signals for both standard single diffusion encoding (SDE) and more advanced double diffusion encoding (DDE) measurements using numerical simulations. The aim is to investigate the ability of dMRI/dMRS to characterize the complex morphology of brain cells focusing on these two distinctive features of brain grey matter. To this end, we employ a recently developed computational framework to create three dimensional meshes of neuron-like structures for Monte Carlo simulations, using diffusion coefficients typical of water and brain metabolites. Modelling the cellular structure as realistically connected spherical soma and cylindrical cellular projections, we cover a wide range of combinations of sphere radii and branching order of cellular projections, characteristic of various grey matter cells. We assess the impact of spherical soma size and branching order on the b-value dependence of the SDE signal as well as the time dependence of the mean diffusivity (MD) and mean kurtosis (MK). Moreover, we also assess the impact of spherical soma size and branching order on the angular modulation of DDE signal at different mixing times, together with the mixing time dependence of the apparent microscopic anisotropy (μA), a promising contrast derived from DDE measurements. The SDE results show that spherical soma size has a measurable impact on both the b-value dependence of the SDE signal and the MD and MK diffusion time dependence for both water and metabolites. On the other hand, we show that branching order has little impact on either, especially for water. In contrast, the DDE results show that spherical soma size has a measurable impact on the DDE signal's angular modulation at short mixing times and the branching order of cellular projections significantly impacts the mixing time dependence of the DDE signal's angular modulation as well as of the derived μA, for both water and metabolites. Our results confirm that SDE based techniques may be sensitive to spherical soma size, and most importantly, show for the first time that DDE measurements may be more sensitive to the dendritic tree complexity (as parametrized by the branching order of cellular projections), paving the way for new ways of characterizing grey matter morphology, non-invasively using dMRS and potentially dMRI.
Collapse
Affiliation(s)
- A Ianus
- Centre for Medical Image Computing and Department of Computer Science, University College London, London, United Kingdom; Champalimaud Research, Champalimaud Centre for the Unknown, Lisbon, Portugal
| | - D C Alexander
- Centre for Medical Image Computing and Department of Computer Science, University College London, London, United Kingdom
| | - H Zhang
- Centre for Medical Image Computing and Department of Computer Science, University College London, London, United Kingdom
| | - M Palombo
- Centre for Medical Image Computing and Department of Computer Science, University College London, London, United Kingdom.
| |
Collapse
|
30
|
Falangola MF, Nie X, Ward R, Dhiman S, Voltin J, Nietert PJ, Jensen JH. Diffusion MRI detects basal forebrain cholinergic abnormalities in the 3xTg-AD mouse model of Alzheimer's disease. Magn Reson Imaging 2021; 83:1-13. [PMID: 34229088 DOI: 10.1016/j.mri.2021.06.022] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2021] [Revised: 06/24/2021] [Accepted: 06/29/2021] [Indexed: 12/17/2022]
Abstract
Degeneration of the basal forebrain (BF) is detected early in the course of Alzheimer's disease (AD). Reduction in the number of BF cholinergic (ChAT) neurons associated with age-related hippocampal cholinergic neuritic dystrophy is described in the 3xTg-AD mouse model; however, no prior diffusion MRI (dMRI) study has explored the presence of BF alterations in this model. Here we investigated the ability of diffusion MRI (dMRI) to detect abnormalities in BF microstructure for the 3xTg-AD mouse model, along with related pathology in the hippocampus (HP) and white matter (WM) tracks comprising the septo-hippocampal pathway. 3xTg-AD and normal control (NC) mice were imaged in vivo using the specific dMRI technique known as diffusional kurtosis imaging (DKI) at 2, 8, and 15 months of age, and 8 dMRI parameters were measured at each time point. Our results revealed significant lower dMRI values in the BF of 2 months-old 3xTg-AD mice compared with NC mice, most likely related to the increased number of ChAT neurons seen in this AD mouse model at this age. They also showed significant age-related dMRI changes in the BF of both groups between 2 and 8 months of age, mainly a decrease in fractional anisotropy and axial diffusivity, and an increase in radial kurtosis. These dMRI changes in the BF may be reflecting the complex aging and pathological microstructural changes described in this region. Group differences and age-related changes were also observed in the HP, fimbria (Fi) and fornix (Fx). In the HP, diffusivity values were significantly higher in the 2 months-old 3xTg-AD mice, and the HP of NC mice showed a significant increase in axial kurtosis after 8 months, reflecting a normal pattern of increased fiber density complexity, which was not seen in the 3xTg-AD mice. In the Fi, mean and radial diffusivity values were significantly higher, and fractional anisotropy, radial kurtosis and kurtosis fractional anisotropy were significantly lower in the 2 months-old 3xTg-AD mice. The age trajectories for both NC and TG mice in the Fi and Fx were similar between 2 and 8 months, but after 8 months there was a significant decrease in diffusivity metrics associated with an increase in kurtosis metrics in the 3xTg-AD mice. These later HP, Fi and Fx dMRI changes probably reflect the growing number of dystrophic neurites and AD pathology progression in the HP, accompanied by WM disruption in the septo-hippocampal pathway. Our results demonstrate that dMRI can detect early cytoarchitectural abnormalities in the BF, as well as related aging and neurodegenerative changes in the HP, Fi and Fx of the 3xTg-AD mice. Since DKI is widely available on clinical scanners, these results also support the potential of the considered dMRI parameters as in vivo biomarkers for AD disease progression.
Collapse
Affiliation(s)
- Maria Fatima Falangola
- Department of Neuroscience, Medical University of South Carolina, Charleston, SC, USA; Center for Biomedical Imaging, Medical University of South Carolina, Charleston, SC, USA.
| | - Xingju Nie
- Department of Neuroscience, Medical University of South Carolina, Charleston, SC, USA; Center for Biomedical Imaging, Medical University of South Carolina, Charleston, SC, USA
| | - Ralph Ward
- Department of Public Health Sciences, Medical University of South Carolina, Charleston, SC, USA
| | - Siddhartha Dhiman
- Department of Neuroscience, Medical University of South Carolina, Charleston, SC, USA
| | - Joshua Voltin
- Department of Neuroscience, Medical University of South Carolina, Charleston, SC, USA
| | - Paul J Nietert
- Department of Public Health Sciences, Medical University of South Carolina, Charleston, SC, USA
| | - Jens H Jensen
- Department of Neuroscience, Medical University of South Carolina, Charleston, SC, USA; Center for Biomedical Imaging, Medical University of South Carolina, Charleston, SC, USA; Department of Radiology and Radiological Science, Medical University of South Carolina, Charleston, SC, USA
| |
Collapse
|
31
|
Scan-rescan repeatability of axonal imaging metrics using high-gradient diffusion MRI and statistical implications for study design. Neuroimage 2021; 240:118323. [PMID: 34216774 PMCID: PMC8646020 DOI: 10.1016/j.neuroimage.2021.118323] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/07/2021] [Revised: 05/12/2021] [Accepted: 06/26/2021] [Indexed: 11/29/2022] Open
Abstract
Axon diameter mapping using diffusion MRI in the living human brain has attracted growing interests with the increasing availability of high gradient strength MRI systems. A systematic assessment of the consistency of axon diameter estimates within and between individuals is needed to gain a comprehensive understanding of how such methods extend to quantifying differences in axon diameter index between groups and facilitate the design of neurobiological studies using such measures. We examined the scan-rescan repeatability of axon diameter index estimation based on the spherical mean technique (SMT) approach using diffusion MRI data acquired with gradient strengths up to 300 mT/m on a 3T Connectom system in 7 healthy volunteers. We performed statistical power analyses using data acquired with the same protocol in a larger cohort consisting of 15 healthy adults to investigate the implications for study design. Results revealed a high degree of repeatability in voxel-wise restricted volume fraction estimates and tract-wise estimates of axon diameter index derived from high-gradient diffusion MRI data. On the region of interest (ROI) level, across white matter tracts in the whole brain, the Pearson’s correlation coefficient of the axon diameter index estimated between scan and rescan experiments was r = 0.72 with an absolute deviation of 0.18 μm. For an anticipated 10% effect size in studies of axon diameter index, most white matter regions required a sample size of less than 15 people to observe a measurable difference between groups using an ROI-based approach. To facilitate the use of high-gradient strength diffusion MRI data for neuroscientific studies of axonal microstructure, the comprehensive multi-gradient strength, multi-diffusion time data used in this work will be made publicly available, in support of open science and increasing the accessibility of such data to the greater scientific community.
Collapse
|
32
|
Mushtaha FN, Kuehn TK, El-Deeb O, Rohani SA, Helpard LW, Moore J, Ladak H, Moehring A, Baron CA, Khan AR. Design and characterization of a 3D-printed axon-mimetic phantom for diffusion MRI. Magn Reson Med 2021; 86:2482-2496. [PMID: 34196049 PMCID: PMC8596689 DOI: 10.1002/mrm.28886] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2020] [Revised: 05/17/2021] [Accepted: 05/18/2021] [Indexed: 01/05/2023]
Abstract
PURPOSE To introduce and characterize inexpensive and easily produced 3D-printed axon-mimetic diffusion MRI phantoms in terms of pore geometry and diffusion kurtosis imaging metrics. METHODS Phantoms were 3D-printed with a composite printing material that, after the dissolution of the polyvinyl alcohol, exhibits microscopic fibrous pores. Confocal microscopy and synchrotron phase-contrast micro-CT imaging were performed to visualize and assess the pore sizes. Diffusion MRI scans of four identical phantoms and phantoms with varying print parameters in water were performed at 9.4 T. Diffusion kurtosis imaging was fit to both data sets and used to assess the reproducibility between phantoms and effects of print parameters on diffusion kurtosis imaging metrics. Identical scans were performed 25 and 76 days later, to test their stability. RESULTS Segmentation of pores in three microscopy images yielded a mean, median, and SD of equivalent pore diameters of 7.57 μm, 3.51 μm, and 12.13 μm, respectively. Phantoms had T1 /T2 = 2 seconds/180 ms, and those with identical parameters showed a low coefficient of variation (~10%) in mean diffusivity (1.38 × 10-3 mm2 /s) and kurtosis (0.52) metrics and radial diffusivity (1.01 × 10-3 mm2 /s) and kurtosis (1.13) metrics. Printing temperature and speed had a small effect on diffusion kurtosis imaging metrics (< 16%), whereas infill density had a larger and more variable effect (> 16%). The stability analysis showed small changes over 2.5 months (< 7%). CONCLUSION Three-dimension-printed axon-mimetic phantoms can mimic the fibrous structure of axon bundles on a microscopic scale, serving as complex, anisotropic diffusion MRI phantoms.
Collapse
Affiliation(s)
- Farah N Mushtaha
- Centre for Functional and Metabolic Mapping, Robarts Research Institute, Western University, London, Canada
| | - Tristan K Kuehn
- Centre for Functional and Metabolic Mapping, Robarts Research Institute, Western University, London, Canada.,School of Biomedical Engineering, Western University, London, Canada
| | - Omar El-Deeb
- Department of Biology, Western University, London, Canada
| | - Seyed A Rohani
- School of Biomedical Engineering, Western University, London, Canada
| | - Luke W Helpard
- School of Biomedical Engineering, Western University, London, Canada
| | - John Moore
- Imaging Research Laboratories, Robarts Research Institute, Western University, London, Canada
| | - Hanif Ladak
- School of Biomedical Engineering, Western University, London, Canada.,Department of Medical Biophysics, Schulich School of Medicine and Dentistry, Western University, London, Canada.,Department of Electrical and Computer Engineering, Western University, London, Canada
| | | | - Corey A Baron
- Centre for Functional and Metabolic Mapping, Robarts Research Institute, Western University, London, Canada.,School of Biomedical Engineering, Western University, London, Canada.,Imaging Research Laboratories, Robarts Research Institute, Western University, London, Canada.,Department of Medical Biophysics, Schulich School of Medicine and Dentistry, Western University, London, Canada.,The Brain and Mind Institute, Western University, London, Canada
| | - Ali R Khan
- Centre for Functional and Metabolic Mapping, Robarts Research Institute, Western University, London, Canada.,Department of Biology, Western University, London, Canada.,Imaging Research Laboratories, Robarts Research Institute, Western University, London, Canada.,Department of Medical Biophysics, Schulich School of Medicine and Dentistry, Western University, London, Canada.,The Brain and Mind Institute, Western University, London, Canada
| |
Collapse
|
33
|
Barakovic M, Girard G, Schiavi S, Romascano D, Descoteaux M, Granziera C, Jones DK, Innocenti GM, Thiran JP, Daducci A. Bundle-Specific Axon Diameter Index as a New Contrast to Differentiate White Matter Tracts. Front Neurosci 2021; 15:646034. [PMID: 34211362 PMCID: PMC8239216 DOI: 10.3389/fnins.2021.646034] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/24/2020] [Accepted: 05/17/2021] [Indexed: 12/30/2022] Open
Abstract
In the central nervous system of primates, several pathways are characterized by different spectra of axon diameters. In vivo methods, based on diffusion-weighted magnetic resonance imaging, can provide axon diameter index estimates non-invasively. However, such methods report voxel-wise estimates, which vary from voxel-to-voxel for the same white matter bundle due to partial volume contributions from other pathways having different microstructure properties. Here, we propose a novel microstructure-informed tractography approach, COMMITAxSize, to resolve axon diameter index estimates at the streamline level, thus making the estimates invariant along trajectories. Compared to previously proposed voxel-wise methods, our formulation allows the estimation of a distinct axon diameter index value for each streamline, directly, furnishing a complementary measure to the existing calculation of the mean value along the bundle. We demonstrate the favourable performance of our approach comparing our estimates with existing histologically-derived measurements performed in the corpus callosum and the posterior limb of the internal capsule. Overall, our method provides a more robust estimation of the axon diameter index of pathways by jointly estimating the microstructure properties of the tissue and the macroscopic organisation of the white matter connectivity.
Collapse
Affiliation(s)
- Muhamed Barakovic
- Signal Processing Lab 5, École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
- Cardiff University Brain Research Imaging Centre, Cardiff University, Cardiff, United Kingdom
- Translational Imaging in Neurology (ThINk) Basel, Department of Biomedical Engineering, University Hospital Basel and University of Basel, Basel, Switzerland
- Neurologic Clinic and Policlinic, Departments of Medicine, Clinical Research and Biomedical Engineering, University Hospital Basel and University of Basel, Basel, Switzerland
| | - Gabriel Girard
- Signal Processing Lab 5, École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
- CIBM Center for BioMedical Imaging, Lausanne, Switzerland
- Radiology Department, Centre Hospitalier Universitaire Vaudois and University of Lausanne, Lausanne, Switzerland
| | - Simona Schiavi
- Signal Processing Lab 5, École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
- Department of Computer Science, University of Verona, Verona, Italy
| | - David Romascano
- Signal Processing Lab 5, École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
| | - Maxime Descoteaux
- Sherbrooke Connectivity Imaging Lab, Université de Sherbrooke, Sherbrooke, QC, Canada
| | - Cristina Granziera
- Translational Imaging in Neurology (ThINk) Basel, Department of Biomedical Engineering, University Hospital Basel and University of Basel, Basel, Switzerland
- Neurologic Clinic and Policlinic, Departments of Medicine, Clinical Research and Biomedical Engineering, University Hospital Basel and University of Basel, Basel, Switzerland
| | - Derek K. Jones
- Cardiff University Brain Research Imaging Centre, Cardiff University, Cardiff, United Kingdom
- Neuroscience and Mental Health Research Institute, Cardiff University, Cardiff, United Kingdom
- Mary MacKillop Institute for Health Research, Australian Catholic University, Melbourne, VIC, Australia
| | - Giorgio M. Innocenti
- Signal Processing Lab 5, École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
- Department of Neuroscience, Karolinska Institutet, Stockholm, Sweden
- Brain and Mind Institute, École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
| | - Jean-Philippe Thiran
- Signal Processing Lab 5, École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
- CIBM Center for BioMedical Imaging, Lausanne, Switzerland
- Radiology Department, Centre Hospitalier Universitaire Vaudois and University of Lausanne, Lausanne, Switzerland
| | | |
Collapse
|
34
|
Veraart J, Raven EP, Edwards LJ, Weiskopf N, Jones DK. The variability of MR axon radii estimates in the human white matter. Hum Brain Mapp 2021; 42:2201-2213. [PMID: 33576105 PMCID: PMC8046139 DOI: 10.1002/hbm.25359] [Citation(s) in RCA: 23] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2020] [Revised: 01/07/2021] [Accepted: 01/21/2021] [Indexed: 12/13/2022] Open
Abstract
The noninvasive quantification of axonal morphology is an exciting avenue for gaining understanding of the function and structure of the central nervous system. Accurate non-invasive mapping of micron-sized axon radii using commonly applied neuroimaging techniques, that is, diffusion-weighted MRI, has been bolstered by recent hardware developments, specifically MR gradient design. Here the whole brain characterization of the effective MR axon radius is presented and the inter- and intra-scanner test-retest repeatability and reproducibility are evaluated to promote the further development of the effective MR axon radius as a neuroimaging biomarker. A coefficient-of-variability of approximately 10% in the voxelwise estimation of the effective MR radius is observed in the test-retest analysis, but it is shown that the performance can be improved fourfold using a customized along-tract analysis.
Collapse
Affiliation(s)
- Jelle Veraart
- Bernard and Irene Schwartz Center for Biomedical Imaging, Department of RadiologyNew York University Grossman School of MedicineNew YorkNew YorkUSA
| | - Erika P. Raven
- Bernard and Irene Schwartz Center for Biomedical Imaging, Department of RadiologyNew York University Grossman School of MedicineNew YorkNew YorkUSA
- CUBRIC, School of PsychologyCardiff UniversityCardiffUK
| | - Luke J. Edwards
- Department of NeurophysicsMax Planck Institute for Human Cognitive and Brain SciencesLeipzigGermany
| | - Nikolaus Weiskopf
- Department of NeurophysicsMax Planck Institute for Human Cognitive and Brain SciencesLeipzigGermany
- Felix Bloch Institute for Solid State Physics, Faculty of Physics and Earth SciencesLeipzig UniversityLeipzigGermany
| | - Derek K. Jones
- CUBRIC, School of PsychologyCardiff UniversityCardiffUK
- Mary MacKillop Institute for Health ResearchAustralian Catholic UniversityMelbourneVictoriaAustralia
| |
Collapse
|
35
|
Abstract
Diffusion-weighted images provide a unique contrast that shows the ability to assess tissue structure and condition on a micrometer scale. Notably, these equations are necessary to understand diffusion MR imaging as a theory but not for real imaging, particularly in clinical practice. The diffusion phenomenon can be observed only through MR measurements. One of the emerging fields of diffusion MRI is to probe the tissue microstructure by altering the diffusion time t, the time interval over which spin displacements are sampled. However, the diffusion time is, in a sense, more important than the b-value for diffusion-weighted images and their quantitative metrics.
Collapse
|
36
|
Ludwig D, Laun FB, Ladd ME, Bachert P, Kuder TA. Apparent exchange rate imaging: On its applicability and the connection to the real exchange rate. Magn Reson Med 2021; 86:677-692. [PMID: 33749019 DOI: 10.1002/mrm.28714] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2020] [Revised: 01/11/2021] [Accepted: 01/15/2021] [Indexed: 12/15/2022]
Abstract
PURPOSE Water exchange between the intracellular and extracellular space can be measured using apparent exchange rate (AXR) imaging. The aim of this study was to investigate the relationship between the measured AXR and the geometry of diffusion restrictions, membrane permeability, and the real exchange rate, as well as to explore the applicability of AXR for typical human measurement settings. METHODS The AXR measurements and the underlying exchange rates were simulated using the Monte Carlo method with different geometries, size distributions, packing densities, and a broad range of membrane permeabilities. Furthermore, the influence of SNR and sequence parameters was analyzed. RESULTS The estimated AXR values correspond to the simulated values and show the expected proportionality to membrane permeability, except for fast exchange (ie, AXR > 20 - 30 s - 1 ) and small packing densities. Moreover, it was found that the duration of the filter gradient must be shorter than 2 · AX R - 1 . In cell size and permeability distributions, AXR depends on the average surface-to-volume ratio, permeability, and the packing density. Finally, AXR can be reliably determined in the presence of orientation dispersion in axon-like structures with sufficient gradient sampling (ie, 30 gradient directions). CONCLUSION Currently used experimental settings for in vivo human measurements are well suited for determining AXR, with the exception of single-voxel analysis, due to limited SNR. The detection of changes in membrane permeability in diseased tissue is nonetheless challenging because of the AXR dependence on further factors, such as packing density and geometry, which cannot be disentangled without further knowledge of the underlying cell structure.
Collapse
Affiliation(s)
- Dominik Ludwig
- Department of Medical Physics in Radiology, German Cancer Research Center (DKFZ), Heidelberg, Germany.,Faculty of Physics and Astronomy, University of Heidelberg, Heidelberg, Germany
| | - Frederik Bernd Laun
- Institute of Radiology, University Hospital Erlangen, Friedrich-Alexander-University Erlangen-Nürnberg, Erlangen, Germany
| | - Mark Edward Ladd
- Department of Medical Physics in Radiology, German Cancer Research Center (DKFZ), Heidelberg, Germany.,Faculty of Physics and Astronomy, University of Heidelberg, Heidelberg, Germany.,Faculty of Medicine, University of Heidelberg, Heidelberg, Germany
| | - Peter Bachert
- Department of Medical Physics in Radiology, German Cancer Research Center (DKFZ), Heidelberg, Germany.,Faculty of Physics and Astronomy, University of Heidelberg, Heidelberg, Germany
| | - Tristan Anselm Kuder
- Department of Medical Physics in Radiology, German Cancer Research Center (DKFZ), Heidelberg, Germany
| |
Collapse
|
37
|
Olesen JL, Østergaard L, Shemesh N, Jespersen SN. Beyond the diffusion standard model in fixed rat spinal cord with combined linear and planar encoding. Neuroimage 2021; 231:117849. [PMID: 33582270 DOI: 10.1016/j.neuroimage.2021.117849] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2020] [Revised: 01/20/2021] [Accepted: 02/04/2021] [Indexed: 10/22/2022] Open
Abstract
Information about tissue on the microscopic and mesoscopic scales can be accessed by modelling diffusion MRI signals, with the aim of extracting microstructure-specific biomarkers. The standard model (SM) of diffusion, currently the most broadly adopted microstructural model, describes diffusion in white matter (WM) tissues by two Gaussian components, one of which has zero radial diffusivity, to represent diffusion in intra- and extra-axonal water, respectively. Here, we reappraise these SM assumptions by collecting comprehensive double diffusion encoded (DDE) MRI data with both linear and planar encodings, which was recently shown to substantially enhance the ability to estimate SM parameters. We find however, that the SM is unable to account for data recorded in fixed rat spinal cord at an ultrahigh field of 16.4 T, suggesting that its underlying assumptions are violated in our experimental data. We offer three model extensions to mitigate this problem: first, we generalize the SM to accommodate finite radii (axons) by releasing the constraint of zero radial diffusivity in the intra-axonal compartment. Second, we include intracompartmental kurtosis to account for non-Gaussian behaviour. Third, we introduce an additional (third) compartment. The ability of these models to account for our experimental data are compared based on parameter feasibility and Bayesian information criterion. Our analysis identifies the three-compartment description as the optimal model. The third compartment exhibits slow diffusion with a minor but non-negligible signal fraction (∼12%). We demonstrate how failure to take the presence of such a compartment into account severely misguides inferences about WM microstructure. Our findings bear significance for microstructural modelling at large and can impact the interpretation of biomarkers extracted from the standard model of diffusion.
Collapse
Affiliation(s)
- Jonas L Olesen
- 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
| | - Leif Østergaard
- Center of Functionally Integrative Neuroscience (CFIN) and MINDLab, Department of Clinical Medicine, Aarhus University, Aarhus, Denmark
| | - Noam Shemesh
- Champalimaud Research, Champalimaud Centre for the Unknown, Lisbon, Portugal
| | - 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.
| |
Collapse
|
38
|
Reymbaut A, Caron AV, Gilbert G, Szczepankiewicz F, Nilsson M, Warfield SK, Descoteaux M, Scherrer B. Magic DIAMOND: Multi-fascicle diffusion compartment imaging with tensor distribution modeling and tensor-valued diffusion encoding. Med Image Anal 2021; 70:101988. [PMID: 33611054 DOI: 10.1016/j.media.2021.101988] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2020] [Revised: 01/25/2021] [Accepted: 01/29/2021] [Indexed: 01/05/2023]
Abstract
Diffusion tensor imaging provides increased sensitivity to microstructural tissue changes compared to conventional anatomical imaging but also presents limited specificity. To tackle this problem, the DIAMOND model subdivides the voxel content into diffusion compartments and draws from diffusion-weighted data to estimate compartmental non-central matrix-variate Gamma distributions of diffusion tensors. It models each sub-voxel fascicle separately, resolving crossing white-matter pathways and allowing for a fascicle-element (fixel) based analysis of microstructural features. Alternatively, specific features of the intra-voxel diffusion tensor distribution can be selectively measured using tensor-valued diffusion-weighted acquisition schemes. However, the impact of such schemes on estimating brain microstructural features has only been studied in a handful of parametric single-fascicle models. In this work, we derive a general Laplace transform for the non-central matrix-variate Gamma distribution, which enables the extension of DIAMOND to tensor-valued encoded data. We then evaluate this "Magic DIAMOND" model in silico and in vivo on various combinations of tensor-valued encoded data. Assessing uncertainty on parameter estimation via stratified bootstrap, we investigate both voxel-based and fixel-based metrics by carrying out multi-peak tractography. We demonstrate using in silico evaluations that tensor-valued diffusion encoding significantly improves Magic DIAMOND's accuracy. Most importantly, we show in vivo that our estimated metrics can be robustly mapped along tracks across regions of fiber crossing, which opens new perspectives for tractometry and microstructure mapping along specific white-matter tracts.
Collapse
Affiliation(s)
| | | | - Guillaume Gilbert
- MR Clinical Science, Philips Healthcare Canada, Markham, ON L6C 2S3, Canada
| | - Filip Szczepankiewicz
- Department of Clinical Sciences, Lund University, 22184, Lund, Sweden; Random Walk Imaging AB, 22224, Lund, Sweden
| | - Markus Nilsson
- Department of Clinical Sciences, Lund University, 22184, Lund, Sweden
| | - Simon K Warfield
- Department of Radiology, Boston Children's Hospital, Boston, MA 02115, United States
| | | | - Benoit Scherrer
- Department of Radiology, Boston Children's Hospital, Boston, MA 02115, United States
| |
Collapse
|
39
|
Harkins KD, Beaulieu C, Xu J, Gore JC, Does MD. A simple estimate of axon size with diffusion MRI. Neuroimage 2020; 227:117619. [PMID: 33301942 PMCID: PMC7949481 DOI: 10.1016/j.neuroimage.2020.117619] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2020] [Revised: 11/06/2020] [Accepted: 11/29/2020] [Indexed: 12/18/2022] Open
Abstract
Noninvasive estimation of mean axon diameter presents a new opportunity to explore white matter plasticity, development, and pathology. Several diffusion-weighted MRI (DW-MRI) methods have been proposed to measure the average axon diameter in white matter, but they typically require many diffusion encoding measurements and complicated mathematical models to fit the signal to multiple tissue compartments, including intra- and extra-axonal spaces. Here, Monte Carlo simulations uncovered a straightforward DW-MRI metric of axon diameter: the change in radial apparent diffusion coefficient estimated at different effective diffusion times, ΔD⊥. Simulations indicated that this metric increases monotonically within a relevant range of effective mean axon diameter while being insensitive to changes in extra-axonal volume fraction, axon diameter distribution, g-ratio, and influence of myelin water. Also, a monotonic relationship was found to exist for signals coming from both intra- and extra-axonal compartments. The slope in ΔD⊥ with effective axon diameter increased with the difference in diffusion time of both oscillating and pulsed gradient diffusion sequences.
Collapse
Affiliation(s)
- Kevin D Harkins
- Biomedical Engineering, Vanderbilt University, United States; Institute of Imaging Science, Vanderbilt University, United States.
| | | | - Junzhong Xu
- Institute of Imaging Science, Vanderbilt University, United States; Radiology and Radiological Sciences, Vanderbilt University Medical Center, United States
| | - John C Gore
- Biomedical Engineering, Vanderbilt University, United States; Institute of Imaging Science, Vanderbilt University, United States; Radiology and Radiological Sciences, Vanderbilt University Medical Center, United States
| | - Mark D Does
- Biomedical Engineering, Vanderbilt University, United States; Institute of Imaging Science, Vanderbilt University, United States
| |
Collapse
|
40
|
Lee HH, Fieremans E, Novikov DS. Realistic Microstructure Simulator (RMS): Monte Carlo simulations of diffusion in three-dimensional cell segmentations of microscopy images. J Neurosci Methods 2020; 350:109018. [PMID: 33279478 DOI: 10.1016/j.jneumeth.2020.109018] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2020] [Revised: 11/16/2020] [Accepted: 11/29/2020] [Indexed: 12/12/2022]
Abstract
BACKGROUND Monte Carlo simulations of diffusion are commonly used as a model validation tool as they are especially suitable for generating the diffusion MRI signal in complicated tissue microgeometries. NEW METHOD Here we describe the details of implementing Monte Carlo simulations in three-dimensional (3d) voxelized segmentations of cells in microscopy images. Using the concept of the corner reflector, we largely reduce the computational load of simulating diffusion within and exchange between multiple cells. Precision is further achieved by GPU-based parallel computations. RESULTS Our simulation of diffusion in white matter axons segmented from a mouse brain demonstrates its value in validating biophysical models. Furthermore, we provide the theoretical background for implementing a discretized diffusion process, and consider the finite-step effects of the particle-membrane reflection and permeation events, needed for efficient simulation of interactions with irregular boundaries, spatially variable diffusion coefficient, and exchange. COMPARISON WITH EXISTING METHODS To our knowledge, this is the first Monte Carlo pipeline for MR signal simulations in a substrate composed of numerous realistic cells, accounting for their permeable and irregularly-shaped membranes. CONCLUSIONS The proposed RMS pipeline makes it possible to achieve fast and accurate simulations of diffusion in realistic tissue microgeometry, as well as the interplay with other MR contrasts. Presently, RMS focuses on simulations of diffusion, exchange, and T1 and T2 NMR relaxation in static tissues, with a possibility to straightforwardly account for susceptibility-induced T2* effects and flow.
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.
| | - 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
| | - 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
| |
Collapse
|
41
|
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: 43] [Impact Index Per Article: 10.8] [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
|
42
|
Novikov DS. The present and the future of microstructure MRI: From a paradigm shift to normal science. J Neurosci Methods 2020; 351:108947. [PMID: 33096152 DOI: 10.1016/j.jneumeth.2020.108947] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2020] [Revised: 08/29/2020] [Accepted: 09/10/2020] [Indexed: 12/29/2022]
Abstract
The aspiration of imaging tissue microstructure with MRI is to uncover micrometer-scale tissue features within millimeter-scale imaging voxels, in vivo. This kind of super-resolution has fueled a paradigm shift within the biomedical imaging community. However, what feels like an ongoing revolution in MRI, has been conceptually experienced in physics decades ago; from this point of view, our current developments can be seen as Thomas Kuhn's "normal science" stage of progress. While the concept of model-based quantification below the nominal imaging resolution is not new, its possibilities in neuroscience and neuroradiology are only beginning to be widely appreciated. This disconnect calls for communicating the progress of tissue microstructure MR imaging to its potential users. Here, a number of recent research developments are outlined in terms of the overarching concept of coarse-graining the tissue structure over an increasing diffusion length. A variety of diffusion models and phenomena are summarized on the phase diagram of diffusion MRI, with the unresolved problems and future directions corresponding to its unexplored domains.
Collapse
Affiliation(s)
- Dmitry S Novikov
- Bernard and Irene Schwartz Center for Biomedical Imaging, Department of Radiology, New York University School of Medicine, New York, NY, USA.
| |
Collapse
|
43
|
Kamiya K, Kamagata K, Ogaki K, Hatano T, Ogawa T, Takeshige-Amano H, Murata S, Andica C, Murata K, Feiweier T, Hori M, Hattori N, Aoki S. Brain White-Matter Degeneration Due to Aging and Parkinson Disease as Revealed by Double Diffusion Encoding. Front Neurosci 2020; 14:584510. [PMID: 33177985 PMCID: PMC7594529 DOI: 10.3389/fnins.2020.584510] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2020] [Accepted: 09/22/2020] [Indexed: 11/16/2022] Open
Abstract
Microstructure imaging by means of multidimensional diffusion encoding is increasingly applied in clinical research, with expectations that it yields a parameter that better correlates with clinical disability than current methods based on single diffusion encoding. Under the assumption that diffusion within a voxel can be well described by a collection of diffusion tensors, several parameters of this diffusion tensor distribution can be derived, including mean size, variance of sizes, orientational dispersion, and microscopic anisotropy. The information provided by multidimensional diffusion encoding also enables us to decompose the sources of the conventional fractional anisotropy and mean kurtosis. In this study, we explored the utility of the diffusion tensor distribution approach for characterizing white-matter degeneration in aging and in Parkinson disease by using double diffusion encoding. Data from 23 healthy older subjects and 27 patients with Parkinson disease were analyzed. Advanced age was associated with greater mean size and size variances, as well as smaller microscopic anisotropy. By analyzing the parameters underlying diffusion kurtosis, we found that the reductions of kurtosis in aging and Parkinson disease reported in the literature are likely driven by the reduction in microscopic anisotropy. Furthermore, microscopic anisotropy correlated with the severity of motor impairment in the patients with Parkinson disease. The present results support the use of multidimensional diffusion encoding in clinical studies and are encouraging for its future clinical implementation.
Collapse
Affiliation(s)
- Kouhei Kamiya
- Department of Radiology, Juntendo University School of Medicine, Tokyo, Japan.,Department of Radiology, Toho University, Tokyo, Japan
| | - Koji Kamagata
- Department of Radiology, Juntendo University School of Medicine, Tokyo, Japan
| | - Kotaro Ogaki
- Department of Neurology, Juntendo University School of Medicine, Tokyo, Japan
| | - Taku Hatano
- Department of Neurology, Juntendo University School of Medicine, Tokyo, Japan
| | - Takashi Ogawa
- Department of Neurology, Juntendo University School of Medicine, Tokyo, Japan
| | | | - Syo Murata
- Department of Radiology, Juntendo University School of Medicine, Tokyo, Japan
| | - Christina Andica
- Department of Radiology, Juntendo University School of Medicine, Tokyo, Japan
| | | | | | - Masaaki Hori
- Department of Radiology, Juntendo University School of Medicine, Tokyo, Japan.,Department of Radiology, Toho University, Tokyo, Japan
| | - Nobutaka Hattori
- Department of Neurology, Juntendo University School of Medicine, Tokyo, Japan
| | - Shigeki Aoki
- Department of Radiology, Juntendo University School of Medicine, Tokyo, Japan
| |
Collapse
|
44
|
Lee HH, Jespersen SN, Fieremans E, Novikov DS. The impact of realistic axonal shape on axon diameter estimation using diffusion MRI. Neuroimage 2020; 223:117228. [PMID: 32798676 PMCID: PMC7806404 DOI: 10.1016/j.neuroimage.2020.117228] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2020] [Accepted: 07/29/2020] [Indexed: 11/24/2022] Open
Abstract
To study axonal microstructure with diffusion MRI, axons are typically modeled as straight impermeable cylinders, whereby the transverse diffusion MRI signal can be made sensitive to the cylinder’s inner diameter. However, the shape of a real axon varies along the axon direction, which couples the longitudinal and transverse diffusion of the overall axon direction. Here we develop a theory of the intra-axonal diffusion MRI signal based on coarse-graining of the axonal shape by 3-dimensional diffusion. We demonstrate how the estimate of the inner diameter is confounded by the diameter variations (beading), and by the local variations in direction (undulations) along the axon. We analytically relate diffusion MRI metrics, such as time-dependent radial diffusivity D⊥(t) and kurtosis K⊥(t), to the axonal shape, and validate our theory using Monte Carlo simulations in synthetic undulating axons with randomly positioned beads, and in realistic axons reconstructed from electron microscopy images of mouse brain white matter. We show that (i) In the narrow pulse limit, the inner diameter from D⊥(t) is overestimated by about twofold due to a combination of axon caliber variations and undulations (each contributing a comparable effect size); (ii) The narrow-pulse kurtosis K⊥∣t→∞ deviates from that in an ideal cylinder due to caliber variations; we also numerically calculate the fourth-order cumulant for an ideal cylinder in the wide pulse limit, which is relevant for inner diameter overestimation; (iii) In the wide pulse limit, the axon diameter overestimation is mainly due to undulations at low diffusion weightings b; and (iv) The effect of undulations can be considerably reduced by directional averaging of high-b signals, with the apparent inner diameter given by a combination of the axon caliber (dominated by the thickest axons), caliber variations, and the residual contribution of undulations.
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.
| | - Sune N Jespersen
- CFIN/MINDLab, Department of Clinical Medicine, Aarhus University, Aarhus, Denmark; Department of Physics and Astronomy, Aarhus University, Aarhus, Denmark
| | - 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
| | - 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
| |
Collapse
|