1
|
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
|
2
|
Schilling KG, Combes AJE, Ramadass K, Rheault F, Sweeney G, Prock L, Sriram S, Cohen-Adad J, Gore JC, Landman BA, Smith SA, O'Grady KP. Influence of preprocessing, distortion correction and cardiac triggering on the quality of diffusion MR images of spinal cord. Magn Reson Imaging 2024; 108:11-21. [PMID: 38309376 PMCID: PMC11218893 DOI: 10.1016/j.mri.2024.01.008] [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: 09/25/2023] [Revised: 01/04/2024] [Accepted: 01/14/2024] [Indexed: 02/05/2024]
Abstract
Diffusion MRI of the spinal cord (SC) is susceptible to geometric distortion caused by field inhomogeneities, and prone to misalignment across time series and signal dropout caused by biological motion. Several modifications of image acquisition and image processing techniques have been introduced to overcome these artifacts, but their specific benefits are largely unproven and warrant further investigations. We aim to evaluate two specific aspects of image acquisition and processing that address image quality in diffusion studies of the spinal cord: susceptibility corrections to reduce geometric distortions, and cardiac triggering to minimize motion artifacts. First, we evaluate 4 distortion preprocessing strategies on 7 datasets of the cervical and lumbar SC and find that while distortion correction techniques increase geometric similarity to structural images, they are largely driven by the high-contrast cerebrospinal fluid, and do not consistently improve the geometry within the cord nor improve white-to-gray matter contrast. We recommend at a minimum to perform bulk-motion correction in preprocessing and posit that improvements/adaptations are needed for spinal cord distortion preprocessing algorithms, which are currently optimized and designed for brain imaging. Second, we design experiments to evaluate the impact of removing cardiac triggering. We show that when triggering is foregone, images are qualitatively similar to triggered sequences, do not have increased prevalence of artifacts, and result in similar diffusion tensor indices with similar reproducibility to triggered acquisitions. When triggering is removed, much shorter acquisitions are possible, which are also qualitatively and quantitatively similar to triggered sequences. We suggest that removing cardiac triggering for cervical SC diffusion can be a reasonable option to save time with minimal sacrifice to image quality.
Collapse
Affiliation(s)
- Kurt G Schilling
- Department of Radiology and Radiological Sciences, Vanderbilt University Medical Center, Nashville, TN, USA; Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, TN, USA.
| | - Anna J E Combes
- Department of Radiology and Radiological Sciences, Vanderbilt University Medical Center, Nashville, TN, USA; Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, TN, USA
| | - Karthik Ramadass
- Department of Electrical and Computer Engineering, Vanderbilt University, Nashville, TN, USA; Department of Computer Science, Vanderbilt University, Nashville, TN, USA
| | - Francois Rheault
- Medical Imaging and Neuroinformatic (MINi) Lab, Department of Computer Science, University of Sherbrooke, Canada
| | - Grace Sweeney
- Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, TN, USA
| | - Logan Prock
- Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, TN, USA
| | - Subramaniam Sriram
- Department of Neurology, Vanderbilt University Medical Center, Nashville, TN, USA
| | - Julien Cohen-Adad
- NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique Montreal, Montreal, QC, Canada; Functional Neuroimaging Unit, CRIUGM, University of Montreal, Montreal, QC, Canada; Mila - Quebec AI Institute, Montreal, QC, Canada; Centre de recherche du CHU Sainte-Justine, Université de Montréal, Montreal, QC, Canada
| | - John C Gore
- Department of Radiology and Radiological Sciences, Vanderbilt University Medical Center, Nashville, TN, USA; Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, TN, USA; Department of Biomedical Engineering, Vanderbilt University, Nashville, TN, USA
| | - Bennett A Landman
- Department of Radiology and Radiological Sciences, Vanderbilt University Medical Center, Nashville, TN, USA; Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, TN, USA; Department of Electrical and Computer Engineering, Vanderbilt University, Nashville, TN, USA; Department of Computer Science, Vanderbilt University, Nashville, TN, USA
| | - Seth A Smith
- Department of Radiology and Radiological Sciences, Vanderbilt University Medical Center, Nashville, TN, USA; Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, TN, USA; Department of Biomedical Engineering, Vanderbilt University, Nashville, TN, USA
| | - Kristin P O'Grady
- Department of Radiology and Radiological Sciences, Vanderbilt University Medical Center, Nashville, TN, USA; Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, TN, USA; Department of Biomedical Engineering, Vanderbilt University, Nashville, TN, USA
| |
Collapse
|
3
|
Barakovic M, Pizzolato M, Tax CMW, Rudrapatna U, Magon S, Dyrby TB, Granziera C, Thiran JP, Jones DK, Canales-Rodríguez EJ. Estimating axon radius using diffusion-relaxation MRI: calibrating a surface-based relaxation model with histology. Front Neurosci 2023; 17:1209521. [PMID: 37638307 PMCID: PMC10457121 DOI: 10.3389/fnins.2023.1209521] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2023] [Accepted: 07/24/2023] [Indexed: 08/29/2023] Open
Abstract
Axon radius is a potential biomarker for brain diseases and a crucial tissue microstructure parameter that determines the speed of action potentials. Diffusion MRI (dMRI) allows non-invasive estimation of axon radius, but accurately estimating the radius of axons in the human brain is challenging. Most axons in the brain have a radius below one micrometer, which falls below the sensitivity limit of dMRI signals even when using the most advanced human MRI scanners. Therefore, new MRI methods that are sensitive to small axon radii are needed. In this proof-of-concept investigation, we examine whether a surface-based axonal relaxation process could mediate a relationship between intra-axonal T2 and T1 times and inner axon radius, as measured using postmortem histology. A unique in vivo human diffusion-T1-T2 relaxation dataset was acquired on a 3T MRI scanner with ultra-strong diffusion gradients, using a strong diffusion-weighting (i.e., b = 6,000 s/mm2) and multiple inversion and echo times. A second reduced diffusion-T2 dataset was collected at various echo times to evaluate the model further. The intra-axonal relaxation times were estimated by fitting a diffusion-relaxation model to the orientation-averaged spherical mean signals. Our analysis revealed that the proposed surface-based relaxation model effectively explains the relationship between the estimated relaxation times and the histological axon radius measured in various corpus callosum regions. Using these histological values, we developed a novel calibration approach to predict axon radius in other areas of the corpus callosum. Notably, the predicted radii and those determined from histological measurements were in close agreement.
Collapse
Affiliation(s)
- Muhamed Barakovic
- Translational Imaging in Neurology (ThINk) Basel, Department of Biomedical Engineering, University Hospital Basel and University of Basel, Basel, Switzerland
- Department of Neurology, University Hospital Basel, Basel, Switzerland
- Cardiff University Brain Research Imaging Centre, Cardiff University, Cardiff, Wales, United Kingdom
- Signal Processing Laboratory 5 (LTS5), Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
- Roche Pharma Research and Early Development, Neuroscience and Rare Diseases, Roche Innovation Center, Basel, Switzerland
| | - Marco Pizzolato
- Department of Applied Mathematics and Computer Science, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Chantal M. W. Tax
- Cardiff University Brain Research Imaging Centre, Cardiff University, Cardiff, Wales, United Kingdom
- Image Sciences Institute, University Medical Center Utrecht, Utrecht, Netherlands
| | - Umesh Rudrapatna
- Cardiff University Brain Research Imaging Centre, Cardiff University, Cardiff, Wales, United Kingdom
| | - Stefano Magon
- Roche Pharma Research and Early Development, Neuroscience and Rare Diseases, Roche Innovation Center, Basel, Switzerland
| | - Tim B. Dyrby
- Department of Applied Mathematics and Computer Science, Technical University of Denmark, Kongens Lyngby, Denmark
- Danish Research Centre for Magnetic Resonance (DRCMR), Centre for Functional and Diagnostic Imaging and Research, Copenhagen University Hospital Amager and Hvidovre, Copenhagen, Denmark
| | - Cristina Granziera
- Translational Imaging in Neurology (ThINk) Basel, Department of Biomedical Engineering, University Hospital Basel and University of Basel, Basel, Switzerland
- Department of Neurology, University Hospital Basel, Basel, Switzerland
- Research Center for Clinical Neuroimmunology and Neuroscience Basel (RC2NB), 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, Cardiff University, Cardiff, Wales, United Kingdom
| | - Erick J. Canales-Rodríguez
- Signal Processing Laboratory 5 (LTS5), Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
| |
Collapse
|
4
|
Alghamdi AJ. The Value of Various Post-Processing Modalities of Diffusion Weighted Imaging in the Detection of Multiple Sclerosis. Brain Sci 2023; 13:brainsci13040622. [PMID: 37190587 DOI: 10.3390/brainsci13040622] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2023] [Revised: 03/31/2023] [Accepted: 04/04/2023] [Indexed: 04/08/2023] Open
Abstract
Diffusion tensor imaging (DTI) showed its adequacy in evaluating the normal-appearing white matter (NAWM) and lesions in the brain that are difficult to evaluate with routine clinical magnetic resonance imaging (MRI) in multiple sclerosis (MS). Recently, MRI systems have been developed with regard to software and hardware, leading to different proposed diffusion analysis methods such as diffusion tensor imaging, q-space imaging, diffusional kurtosis imaging, neurite orientation dispersion and density imaging, and axonal diameter measurement. These methods have the ability to better detect in vivo microstructural changes in the brain than DTI. These different analysis modalities could provide supplementary inputs for MS disease characterization and help in monitoring the disease’s progression as well as treatment efficacy. This paper reviews some of the recent diffusion MRI methods used for the assessment of MS in vivo.
Collapse
Affiliation(s)
- Ahmad Joman Alghamdi
- Radiological Sciences Department, College of Applied Medical Sciences, Taif University, Taif 21944, Saudi Arabia
| |
Collapse
|
5
|
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
|
6
|
Schilling KG, Fadnavis S, Batson J, Visagie M, Combes AJE, By S, McKnight CD, Bagnato F, Garyfallidis E, Landman BA, Smith SA, O'Grady KP. Denoising of diffusion MRI in the cervical spinal cord - effects of denoising strategy and acquisition on intra-cord contrast, signal modeling, and feature conspicuity. Neuroimage 2023; 266:119826. [PMID: 36543265 PMCID: PMC9843739 DOI: 10.1016/j.neuroimage.2022.119826] [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: 04/07/2022] [Revised: 11/02/2022] [Accepted: 12/17/2022] [Indexed: 12/23/2022] Open
Abstract
Quantitative diffusion MRI (dMRI) is a promising technique for evaluating the spinal cord in health and disease. However, low signal-to-noise ratio (SNR) can impede interpretation and quantification of these images. The purpose of this study is to evaluate several dMRI denoising approaches on their ability to improve the quality, reliability, and accuracy of quantitative diffusion MRI of the spinal cord. We evaluate three denoising approaches (Non-Local Means, Marchenko-Pastur PCA, and a newly proposed Patch2Self algorithm) and conduct five experiments to validate the denoising performance on clinical-quality and commonly-acquired dMRI acquisitions: 1) a phantom experiment to assess denoising error and bias; 2) a multi-vendor, multi-acquisition open experiment for both qualitative and quantitative evaluation of noise residuals; 3) a bootstrapping experiment to estimate uncertainty of parametric maps; 4) an assessment of spinal cord lesion conspicuity in a multiple sclerosis group; and 5) an evaluation of denoising for advanced parametric multi-compartment modeling. We find that all methods improve signal-to-noise ratio and conspicuity of MS lesions in individual diffusion weighted images (DWIs), but MPPCA and Patch2Self excel at improving the quality and intra-cord contrast of diffusion weighted images - removing signal fluctuations due to thermal noise while improving precision of estimation of diffusion parameters even with very few DWIs (i.e., 16-32) typical of clinical acquisitions. These denoising approaches hold promise for facilitating reliable diffusion observations and measurements in the spinal cord to investigate biological and pathological processes.
Collapse
Affiliation(s)
- Kurt G Schilling
- Department of Radiology and Radiological Sciences, Vanderbilt University Medical Center, Nashville, TN, United States; Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, TN, United States.
| | - Shreyas Fadnavis
- Department of Intelligent Systems Engineering, Indiana University Bloomington, Bloomington, IN, United States
| | - Joshua Batson
- The Public Health Company, California, United States
| | - Mereze Visagie
- Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, TN, United States
| | - Anna J E Combes
- Department of Radiology and Radiological Sciences, Vanderbilt University Medical Center, Nashville, TN, United States; Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, TN, United States
| | - Samantha By
- Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, TN, United States; Department of Biomedical Engineering, Vanderbilt University, Nashville, TN, United States
| | - Colin D McKnight
- Department of Radiology and Radiological Sciences, Vanderbilt University Medical Center, Nashville, TN, United States
| | - Francesca Bagnato
- Neuroimaging Unit, Neuroimmunology Division, Department of Neurology, Vanderbilt University Medical Center, Nashville, TN, United States; Department of Neurology, VA Hospital, TN Valley Healthcare System, Nashville, TN, United States
| | - Eleftherios Garyfallidis
- Department of Intelligent Systems Engineering, Indiana University Bloomington, Bloomington, IN, United States
| | - Bennett A Landman
- Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, TN, United States; Department of Biomedical Engineering, Vanderbilt University, Nashville, TN, United States; Department of Electrical and Computer Engineering, Vanderbilt University, Nashville, TN, United States
| | - Seth A Smith
- Department of Radiology and Radiological Sciences, Vanderbilt University Medical Center, Nashville, TN, United States; Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, TN, United States; Department of Biomedical Engineering, Vanderbilt University, Nashville, TN, United States
| | - Kristin P O'Grady
- Department of Radiology and Radiological Sciences, Vanderbilt University Medical Center, Nashville, TN, United States; Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, TN, United States; Department of Biomedical Engineering, Vanderbilt University, Nashville, TN, United States.
| |
Collapse
|
7
|
Resolution and b value dependent Structural Connectome in ex vivo Mouse Brain. Neuroimage 2022; 255:119199. [PMID: 35417754 PMCID: PMC9195912 DOI: 10.1016/j.neuroimage.2022.119199] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/19/2021] [Revised: 04/05/2022] [Accepted: 04/07/2022] [Indexed: 12/24/2022] Open
Abstract
Diffusion magnetic resonance imaging has been widely used in both clinical and preclinical studies to characterize tissue microstructure and structural connectivity. The diffusion MRI protocol for the Human Connectome Project (HCP) has been developed and optimized to obtain high-quality, high-resolution diffusion MRI (dMRI) datasets. However, such efforts have not been fully explored in preclinical studies, especially for rodents. In this study, high quality dMRI datasets of mouse brains were acquired at 9.4T system from two vendors. In particular, we acquired a high-spatial resolution dMRI dataset (25 μm isotropic with 126 diffusion encoding directions), which we believe to be the highest spatial resolution yet obtained; and a high-angular resolution dMRI dataset (50 μm isotropic with 384 diffusion encoding directions), which we believe to be the highest angular resolution compared to the dMRI datasets at the microscopic resolution. We systematically investigated the effects of three important parameters that affect the final outcome of the connectome: b value (1000s/mm2 to 8000 s/mm2), angular resolution (10 to 126), and spatial resolution (25 μm to 200 μm). The stability of tractography and connectome increase with the angular resolution, where more than 50 angles is necessary to achieve consistent results. The connectome and quantitative parameters derived from graph theory exhibit a linear relationship to the b value (R2 > 0.99); a single-shell acquisition with b value of 3000 s/mm2 shows comparable results to the multi-shell high angular resolution dataset. The dice coefficient decreases and both false positive rate and false negative rate gradually increase with coarser spatial resolution. Our study provides guidelines and foundations for exploration of tradeoffs among acquisition parameters for the structural connectome in ex vivo mouse brain.
Collapse
|
8
|
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 DOI: 10.1016/j.neuroimage.2022.118958] [Citation(s) in RCA: 22] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Key Words] [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 dMRI 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 dMRI 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, United States; 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, Netherlands
| | - Mathias Davids
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, United States; 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, United States; Harvard Medical School, Boston, MA, USA; Harvard-MIT Division of Health Sciences and Technology, Massachusetts Institute of Technology, Cambridge, MA, United States
| | - 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, United States; Harvard Medical School, Boston, MA, USA
| | - Qiyuan Tian
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, United States; Harvard Medical School, Boston, MA, USA
| | - Chiara Maffei
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, United States; 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, United States; Harvard Medical School, Boston, MA, USA
| | - Aapo Nummenmaa
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, United States; Harvard Medical School, Boston, MA, USA
| | | | - Anastasia Yendiki
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, United States; Harvard Medical School, Boston, MA, USA
| | - Yi-Qiao Song
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, United States; 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, United States; Harvard Medical School, Boston, MA, USA; Harvard-MIT Division of Health Sciences and Technology, Massachusetts Institute of Technology, Cambridge, MA, United States
| | - Lawrence L Wald
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, United States; Harvard Medical School, Boston, MA, USA; Harvard-MIT Division of Health Sciences and Technology, Massachusetts Institute of Technology, Cambridge, MA, United States
| | - Susie Y Huang
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, United States; Harvard Medical School, Boston, MA, USA; Harvard-MIT Division of Health Sciences and Technology, Massachusetts Institute of Technology, Cambridge, MA, United States.
| |
Collapse
|
9
|
Oliva V, Hartley-Davies R, Moran R, Pickering AE, Brooks JC. Simultaneous brain, brainstem and spinal cord pharmacological-fMRI reveals involvement of an endogenous opioid network in attentional analgesia. eLife 2022; 11:71877. [PMID: 35080494 PMCID: PMC8843089 DOI: 10.7554/elife.71877] [Citation(s) in RCA: 20] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2021] [Accepted: 01/25/2022] [Indexed: 11/13/2022] Open
Abstract
Pain perception is decreased by shifting attentional focus away from a threatening event. This attentional analgesia engages parallel descending control pathways from anterior cingulate (ACC) to locus coeruleus, and ACC to periaqueductal grey (PAG) – rostral ventromedial medulla (RVM), indicating possible roles for noradrenergic or opioidergic neuromodulators. To determine which pathway modulates nociceptive activity in humans, we used simultaneous whole brain-spinal cord pharmacological-fMRI (N = 39) across three sessions. Noxious thermal forearm stimulation generated somatotopic-activation of dorsal horn (DH) whose activity correlated with pain report and mirrored attentional pain modulation. Activity in an adjacent cluster reported the interaction between task and noxious stimulus. Effective connectivity analysis revealed that ACC interacts with PAG and RVM to modulate spinal cord activity. Blocking endogenous opioids with Naltrexone impairs attentional analgesia and disrupts RVM-spinal and ACC-PAG connectivity. Noradrenergic augmentation with Reboxetine did not alter attentional analgesia. Cognitive pain modulation involves opioidergic ACC-PAG-RVM descending control which suppresses spinal nociceptive activity.
Collapse
Affiliation(s)
- Valeria Oliva
- Department of Anesthesiology, University of California, San Diego, La Jolla, United States
| | - Ron Hartley-Davies
- School of Psychological Science, University of Bristol, Bristol, United Kingdom
| | - Rosalyn Moran
- Department of Neuroimaging, King's College London, London, United Kingdom
| | - Anthony E Pickering
- School of Physiology, Pharmacology and Neuroscience, University of Bristol, Bristol, United Kingdom
| | | |
Collapse
|
10
|
Tian Q, Fan Q, Witzel T, Polackal MN, Ohringer NA, Ngamsombat C, Russo AW, Machado N, Brewer K, Wang F, Setsompop K, Polimeni JR, Keil B, Wald LL, Rosen BR, Klawiter EC, Nummenmaa A, Huang SY. Comprehensive diffusion MRI dataset for in vivo human brain microstructure mapping using 300 mT/m gradients. Sci Data 2022; 9:7. [PMID: 35042861 PMCID: PMC8766594 DOI: 10.1038/s41597-021-01092-6] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2021] [Accepted: 10/25/2021] [Indexed: 12/27/2022] Open
Abstract
Strong gradient systems can improve the signal-to-noise ratio of diffusion MRI measurements and enable a wider range of acquisition parameters that are beneficial for microstructural imaging. We present a comprehensive diffusion MRI dataset of 26 healthy participants acquired on the MGH-USC 3 T Connectome scanner equipped with 300 mT/m maximum gradient strength and a custom-built 64-channel head coil. For each participant, the one-hour long acquisition systematically sampled the accessible diffusion measurement space, including two diffusion times (19 and 49 ms), eight gradient strengths linearly spaced between 30 mT/m and 290 mT/m for each diffusion time, and 32 or 64 uniformly distributed directions. The diffusion MRI data were preprocessed to correct for gradient nonlinearity, eddy currents, and susceptibility induced distortions. In addition, scan/rescan data from a subset of seven individuals were also acquired and provided. The MGH Connectome Diffusion Microstructure Dataset (CDMD) may serve as a test bed for the development of new data analysis methods, such as fiber orientation estimation, tractography and microstructural modelling.
Collapse
Affiliation(s)
- 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
| | - 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
| | - Thomas Witzel
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, Massachusetts, United States
- Harvard Medical School, Boston, Massachusetts, United States
| | - Maya N Polackal
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, Massachusetts, United States
| | - Ned A Ohringer
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, Massachusetts, United States
| | - Chanon Ngamsombat
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, Massachusetts, United States
| | - Andrew W Russo
- Department of Neurology, Massachusetts General Hospital, Boston, Massachusetts, United States
| | - Natalya Machado
- Department of Neurology, Massachusetts General Hospital, Boston, Massachusetts, United States
| | - Kristina Brewer
- Department of Neurology, Massachusetts General Hospital, Boston, Massachusetts, United States
| | - Fuyixue Wang
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, Massachusetts, United States
- Harvard-MIT Division of Health Sciences and Technology, Massachusetts Institute of Technology, Cambridge, Massachusetts, United States
| | - Kawin Setsompop
- 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
| | - Jonathan R Polimeni
- 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
| | - Boris Keil
- Department of Life Science Engineering, Institute of Medical Physics and Radiation Protection, Giessen, Germany
| | - Lawrence L Wald
- 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
| | - Bruce R Rosen
- 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
| | - Eric C Klawiter
- 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 Neurology, Massachusetts General Hospital, Boston, Massachusetts, United States
| | - Aapo Nummenmaa
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, Massachusetts, United States
- Harvard Medical School, Boston, Massachusetts, United States
| | - 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.
| |
Collapse
|
11
|
High-fidelity diffusion tensor imaging of the cervical spinal cord using point-spread-function encoded EPI. Neuroimage 2021; 236:118043. [PMID: 33857617 DOI: 10.1016/j.neuroimage.2021.118043] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2021] [Revised: 03/16/2021] [Accepted: 03/31/2021] [Indexed: 01/21/2023] Open
Abstract
Diffusion tensor imaging (DTI) of the spinal cord is technically challenging due to the size of its structure and susceptibility-induced field inhomogeneity, which impedes clinical applications. This study aimed to achieve high-fidelity spinal cord DTI with reasonable SNR and practical acquisition efficiency. Particularly, a distortion-free multi-shot EPI technique, namely point-spread-function encoded EPI (PSF-EPI), was adopted for diffusion imaging of the cervical spinal cord (CSC). The shot number can be reduced to six for sagittal scans through titled-CAIPI acceleration and partial Fourier undersampling, consequently rendering this technique beneficial in clinics. Fifteen healthy volunteers and seven patients with metallic implants underwent sagittal scans using tilted-CAIPI PSF-EPI at 3T. Unsuppressed fat signals were further removed by retrospective water/fat separation using the intrinsic chemical-shift encoded signals. Compared with multi-shot interleaved EPI method, highly accelerated PSF-EPI method provided evidently improved distortion reduction and higher consistency with anatomical references even with metallic implants. Additionally, axial DTI scans using PSF-EPI were also evaluated quantitatively, and the measured DTI metrics are similar to those obtained from the zonal oblique multi-slice EPI (ZOOM-EPI) method and reported values. The high anatomical consistency, practical scan time and quantitative reliability indicate PSF-EPI's clinical potential for CSC diffusion imaging.
Collapse
|
12
|
Shahrampour S, De Leener B, Alizadeh M, Middleton D, Krisa L, Flanders AE, Faro SH, Cohen-Adad J, Mohamed FB. Atlas-Based Quantification of DTI Measures in a Typically Developing Pediatric Spinal Cord. AJNR Am J Neuroradiol 2021; 42:1727-1734. [PMID: 34326104 DOI: 10.3174/ajnr.a7221] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2021] [Accepted: 04/19/2021] [Indexed: 11/07/2022]
Abstract
BACKGROUND AND PURPOSE Multi-parametric MRI, provides a variety of biomarkers sensitive to white matter integrity, However, spinal cord MRI data in pediatrics is rare compared to adults. The purpose of this work was 3-fold: 1) to develop a processing pipeline for atlas-based generation of the typically developing pediatric spinal cord WM tracts, 2) to derive atlas-based normative values of the DTI indices for various WM pathways, and 3) to investigate age-related changes in the obtained normative DTI indices along the extracted tracts. MATERIALS AND METHODS DTI scans of 30 typically developing subjects (age range, 6-16 years) were acquired on a 3T MR imaging scanner. The data were registered to the PAM50 template in the Spinal Cord Toolbox. Next, the DTI indices for various WM regions were extracted at a single section centered at the C3 vertebral body in all the 30 subjects. Finally, an ANOVA test was performed to examine the effects of the following: 1) laterality, 2) functionality, and 3) age, with DTI-derived indices in 34 extracted WM regions. RESULTS A postprocessing pipeline was developed and validated to delineate pediatric spinal cord WM tracts. The results of ANOVA on fractional anisotropy values showed no effect for laterality (P = .72) but an effect for functionality (P < .001) when comparing the 30 primary WM labels. There was a significant (P < .05) effect of age and maturity of the left spinothalamic tract on mean diffusivity, radial diffusivity, and axial diffusivity values. CONCLUSIONS The proposed automated pipeline in this study incorporates unique postprocessing steps followed by template registration and quantification of DTI metrics using atlas-based regions. This method eliminates the need for manual ROI analysis of WM tracts and, therefore, increases the accuracy and speed of the measurements.
Collapse
Affiliation(s)
- S Shahrampour
- From the Departments of Radiology (S.S., M.A., D.M., F.B.M.)
| | - B De Leener
- Department of Computer Engineering and Software Engineering (B.D.L.)
| | - M Alizadeh
- From the Departments of Radiology (S.S., M.A., D.M., F.B.M.)
| | - D Middleton
- From the Departments of Radiology (S.S., M.A., D.M., F.B.M.)
| | | | - A E Flanders
- Radiology (A.E.F., S.H.F.), Thomas Jefferson University, Philadelphia, Pennsylvania
| | - S H Faro
- Radiology (A.E.F., S.H.F.), Thomas Jefferson University, Philadelphia, Pennsylvania
| | - J Cohen-Adad
- NeuroPoly Lab (J.C.-A.), Institute of Biomedical Engineering, Polytechnique Montreal, Montreal, Quebec, Canada.,Functional Neuroimaging Unit (J.C.-A.), Centre de Recherche de l'Institut Universitaire de Gériatrie de Montréal, Université de Montréal, Montreal, Quebec, Canada
| | - F B Mohamed
- From the Departments of Radiology (S.S., M.A., D.M., F.B.M.)
| |
Collapse
|
13
|
Structural and resting state functional connectivity beyond the cortex. Neuroimage 2021; 240:118379. [PMID: 34252527 DOI: 10.1016/j.neuroimage.2021.118379] [Citation(s) in RCA: 23] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/08/2021] [Revised: 05/21/2021] [Accepted: 07/07/2021] [Indexed: 12/14/2022] Open
Abstract
Mapping the structural and functional connectivity of the central nervous system has become a key area within neuroimaging research. While detailed network structures across the entire brain have been probed using animal models, non-invasive neuroimaging in humans has thus far been dominated by cortical investigations. Beyond the cortex, subcortical nuclei have traditionally been less accessible due to their smaller size and greater distance from radio frequency coils. However, major neuroimaging developments now provide improved signal and the resolution required to study these structures. Here, we present an overview of the connectivity between the amygdala, brainstem, cerebellum, spinal cord and the rest of the brain. While limitations to their imaging and analyses remain, we also provide some recommendations and considerations for mapping brain connectivity beyond the cortex.
Collapse
|
14
|
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
|
15
|
Tinnermann A, Büchel C, Cohen-Adad J. Cortico-spinal imaging to study pain. Neuroimage 2020; 224:117439. [PMID: 33039624 DOI: 10.1016/j.neuroimage.2020.117439] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2020] [Revised: 09/21/2020] [Accepted: 10/01/2020] [Indexed: 12/15/2022] Open
Abstract
Functional magnetic resonance imaging of the brain has helped to reveal mechanisms of pain perception in health and disease. Recently, imaging approaches have been developed that allow recording neural activity simultaneously in the brain and in the spinal cord. These approaches offer the possibility to examine pain perception in the entire central pain system and in addition, to investigate cortico-spinal interactions during pain processing. Although cortico-spinal imaging is a promising technique, it bears challenges concerning data acquisition and data analysis strategies. In this review, we discuss studies that applied simultaneous imaging of the brain and spinal cord to explore central pain processing. Furthermore, we describe different MR-related acquisition techniques, summarize advantages and disadvantages of approaches that have been implemented so far and present software that has been specifically developed for the analysis of spinal fMRI data to address challenges of spinal data analysis.
Collapse
Affiliation(s)
- Alexandra Tinnermann
- Department for Systems Neuroscience, University Medical Center Hamburg-Eppendorf, Hamburg, Germany; Max Planck School of Cognition, Leipzig, Germany.
| | - Christian Büchel
- Department for Systems Neuroscience, University Medical Center Hamburg-Eppendorf, Hamburg, Germany; Max Planck School of Cognition, Leipzig, Germany
| | - Julien Cohen-Adad
- NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique Montreal, Montreal, Quebec, Canada; Functional Neuroimaging Unit, CRIUGM, Université de Montréal, Montreal, Quebec, Canada.
| |
Collapse
|
16
|
Jelescu IO, Palombo M, Bagnato F, Schilling KG. Challenges for biophysical modeling of microstructure. J Neurosci Methods 2020; 344:108861. [PMID: 32692999 PMCID: PMC10163379 DOI: 10.1016/j.jneumeth.2020.108861] [Citation(s) in RCA: 63] [Impact Index Per Article: 15.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2020] [Revised: 07/10/2020] [Accepted: 07/14/2020] [Indexed: 02/07/2023]
Abstract
The biophysical modeling efforts in diffusion MRI have grown considerably over the past 25 years. In this review, we dwell on the various challenges along the journey of bringing a biophysical model from initial design to clinical implementation, identifying both hurdles that have been already overcome and outstanding issues. First, we describe the critical initial task of selecting which features of tissue microstructure can be estimated using a model and which acquisition protocol needs to be implemented to make the estimation possible. The model performance should necessarily be tested in realistic numerical simulations and in experimental data - adapting the fitting strategy accordingly, and parameter estimates should be validated against complementary techniques, when/if available. Secondly, the model performance and validity should be explored in pathological conditions, and, if appropriate, dedicated models for pathology should be developed. We build on examples from tumors, ischemia and demyelinating diseases. We then discuss the challenges associated with clinical translation and added value. Finally, we single out four major unresolved challenges that are related to: the availability of a microstructural ground truth, the validation of model parameters which cannot be accessed with complementary techniques, the development of a generalized standard model for any brain region and pathology, and the seamless communication between different parties involved in the development and application of biophysical models of diffusion.
Collapse
|
17
|
Grussu F, Battiston M, Veraart J, Schneider T, Cohen-Adad J, Shepherd TM, Alexander DC, Fieremans E, Novikov DS, Gandini Wheeler-Kingshott CAM. Multi-parametric quantitative in vivo spinal cord MRI with unified signal readout and image denoising. Neuroimage 2020; 217:116884. [PMID: 32360689 PMCID: PMC7378937 DOI: 10.1016/j.neuroimage.2020.116884] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2019] [Revised: 03/18/2020] [Accepted: 04/23/2020] [Indexed: 12/11/2022] Open
Abstract
Multi-parametric quantitative MRI (qMRI) of the spinal cord is a promising non-invasive tool to probe early microstructural damage in neurological disorders. It is usually performed in vivo by combining acquisitions with multiple signal readouts, which exhibit different thermal noise levels, geometrical distortions and susceptibility to physiological noise. This ultimately hinders joint multi-contrast modelling and makes the geometric correspondence of parametric maps challenging. We propose an approach to overcome these limitations, by implementing state-of-the-art microstructural MRI of the spinal cord with a unified signal readout in vivo (i.e. with matched spatial encoding parameters across a range of imaging contrasts). We base our acquisition on single-shot echo planar imaging with reduced field-of-view, and obtain data from two different vendors (vendor 1: Philips Achieva; vendor 2: Siemens Prisma). Importantly, the unified acquisition allows us to compare signal and noise across contrasts, thus enabling overall quality enhancement via multi-contrast image denoising methods. As a proof-of-concept, here we provide a demonstration with one such method, known as Marchenko-Pastur (MP) Principal Component Analysis (PCA) denoising. MP-PCA is a singular value (SV) decomposition truncation approach that relies on redundant acquisitions, i.e. such that the number of measurements is large compared to the number of components that are maintained in the truncated SV decomposition. Here we used in vivo and synthetic data to test whether a unified readout enables more efficient MP-PCA denoising of less redundant acquisitions, since these can be denoised jointly with more redundant ones. We demonstrate that a unified readout provides robust multi-parametric maps, including diffusion and kurtosis tensors from diffusion MRI, myelin metrics from two-pool magnetisation transfer, and T1 and T2 from relaxometry. Moreover, we show that MP-PCA improves the quality of our multi-contrast acquisitions, since it reduces the coefficient of variation (i.e. variability) by up to 17% for mean kurtosis, 8% for bound pool fraction (myelin-sensitive), and 13% for T1, while enabling more efficient denoising of modalities limited in redundancy (e.g. relaxometry). In conclusion, multi-parametric spinal cord qMRI with unified readout is feasible and provides robust microstructural metrics with matched resolution and distortions, whose quality benefits from multi-contrast denoising methods such as MP-PCA.
Collapse
Affiliation(s)
- Francesco Grussu
- Queen Square MS Centre, UCL Queen Square Institute of Neurology, Faculty of Brain Sciences, University College London, London, UK; Centre for Medical Image Computing, Department of Computer Science, University College London, London, UK.
| | - Marco Battiston
- Queen Square MS Centre, UCL Queen Square Institute of Neurology, Faculty of Brain Sciences, University College London, London, UK
| | - Jelle Veraart
- Center for Biomedical Imaging, Department of Radiology, New York University School of Medicine, New York, USA
| | | | - Julien Cohen-Adad
- NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique Montreal, Montreal, Canada; Functional Neuroimaging Unit, CRIUGM, Université de Montréal, Montreal, Canada
| | - Timothy M Shepherd
- Center for Biomedical Imaging, Department of Radiology, New York University School of Medicine, New York, USA
| | - Daniel C Alexander
- Centre for Medical Image Computing, Department of Computer Science, University College London, London, UK
| | - Els Fieremans
- Center for Biomedical Imaging, Department of Radiology, New York University School of Medicine, New York, USA
| | - Dmitry S Novikov
- Center for Biomedical Imaging, Department of Radiology, New York University School of Medicine, New York, USA
| | - Claudia A M Gandini Wheeler-Kingshott
- Queen Square MS Centre, UCL Queen Square Institute of Neurology, Faculty of Brain Sciences, University College London, London, UK; Brain MRI 3T Research Centre, IRCCS Mondino Foundation, Pavia, Italy; Department of Brain and Behavioural Sciences, University of Pavia, Pavia, Italy
| |
Collapse
|
18
|
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
|
19
|
Lakhani DA, Schilling KG, Xu J, Bagnato F. Advanced Multicompartment Diffusion MRI Models and Their Application in Multiple Sclerosis. AJNR Am J Neuroradiol 2020; 41:751-757. [PMID: 32354707 DOI: 10.3174/ajnr.a6484] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2019] [Accepted: 02/03/2020] [Indexed: 01/22/2023]
Abstract
Conventional MR imaging techniques are sensitive to pathologic changes of the brain and spinal cord seen in MS, but they lack specificity for underlying axonal and myelin integrity. By isolating the signal contribution from different tissue compartments, newly developed advanced multicompartment diffusion MR imaging models have the potential to detect specific tissue subtypes and associated injuries with increased pathologic specificity. These models include neurite orientation dispersion and density imaging, diffusion basis spectrum imaging, multicompartment microscopic diffusion MR imaging with the spherical mean technique, and models enabled through high-gradient diffusion MR imaging. In this review, we provide an appraisal of the current literature on the physics principles, histopathologic validation, and clinical applications of each of these techniques in both brains and spinal cords of patients with MS. We discuss limitations of each of the methods and directions that future research could take to provide additional validation of their roles as biomarkers of axonal and myelin injury in MS.
Collapse
Affiliation(s)
- D A Lakhani
- From the Neuroimaging Unit (D.A.L., F.B.), Neuroimmunology Division, Department of Neurology
- Division of Internal Medicine (D.A.L.)
- Department of Radiology (D.A.L.), West Virginia University, Morgantown, West Virginia
| | - K G Schilling
- Department of Radiology and Radiological Sciences (K.G.S., J.X.), Vanderbilt University Institute of Imaging Sciences, Vanderbilt University Medical Center, Nashville, Tennessee
| | - J Xu
- Department of Radiology and Radiological Sciences (K.G.S., J.X.), Vanderbilt University Institute of Imaging Sciences, Vanderbilt University Medical Center, Nashville, Tennessee
| | - F Bagnato
- From the Neuroimaging Unit (D.A.L., F.B.), Neuroimmunology Division, Department of Neurology
- Department of Neurology (F.B.), VA Tennessee Valley Healthcare System, Nashville, Tennessee
| |
Collapse
|
20
|
Anterior fissure, central canal, posterior septum and more: New insights into the cervical spinal cord gray and white matter regional organization using T1 mapping at 7T. Neuroimage 2020; 205:116275. [DOI: 10.1016/j.neuroimage.2019.116275] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2019] [Revised: 09/16/2019] [Accepted: 10/10/2019] [Indexed: 12/12/2022] Open
|
21
|
Anaby D, Morozov D, Seroussi I, Hametner S, Sochen N, Cohen Y. Single- and double-Diffusion encoding MRI for studying ex vivo apparent axon diameter distribution in spinal cord white matter. NMR IN BIOMEDICINE 2019; 32:e4170. [PMID: 31573745 DOI: 10.1002/nbm.4170] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/24/2019] [Revised: 07/28/2019] [Accepted: 07/31/2019] [Indexed: 06/10/2023]
Abstract
Mapping average axon diameter (AAD) and axon diameter distribution (ADD) in neuronal tissues non-invasively is a challenging task that may have a tremendous effect on our understanding of the normal and diseased central nervous system (CNS). Water diffusion is used to probe microstructure in neuronal tissues, however, the different water populations and barriers that are present in these tissues turn this into a complex task. Therefore, it is not surprising that recently we have witnessed a burst in the development of new approaches and models that attempt to obtain, non-invasively, detailed microstructural information in the CNS. In this work, we aim at challenging and comparing the microstructural information obtained from single diffusion encoding (SDE) with double diffusion encoding (DDE) MRI. We first applied SDE and DDE MR spectroscopy (MRS) on microcapillary phantoms and then applied SDE and DDE MRI on an ex vivo porcine spinal cord (SC), using similar experimental conditions. The obtained diffusion MRI data were fitted by the same theoretical model, assuming that the signal in every voxel can be approximated as the superposition of a Gaussian-diffusing component and a series of restricted components having infinite cylindrical geometries. The diffusion MRI results were then compared with histological findings. We found a good agreement between the fittings and the experimental data in white matter (WM) voxels of the SC in both diffusion MRI methods. The microstructural information and apparent AADs extracted from SDE MRI were found to be similar or somewhat larger than those extracted from DDE MRI especially when the diffusion time was set to 40 ms. The apparent ADDs extracted from SDE and DDE MRI show reasonable agreement but somewhat weaker correspondence was observed between the diffusion MRI results and histology. The apparent subtle differences between the microstructural information obtained from SDE and DDE MRI are briefly discussed.
Collapse
Affiliation(s)
- Debbie Anaby
- School of Chemistry, The Raymond and Beverly Sackler Faculty of Exact Sciences, Tel Aviv University, Tel Aviv, Israel
- Department of Diagnostic Imaging, Sheba Medical Center, Tel HaShomer, Israel
| | - Darya Morozov
- School of Chemistry, The Raymond and Beverly Sackler Faculty of Exact Sciences, Tel Aviv University, Tel Aviv, Israel
| | - Inbar Seroussi
- School of Mathematical Sciences, The Raymond and Beverly Sackler Faculty of Exact Sciences, Tel Aviv University, Tel Aviv, Israel
| | - Simon Hametner
- Neuroimmunology Department, Center of Brain Research, Medical University of Vienna, Vienna, Austria
| | - Nir Sochen
- School of Mathematical Sciences, The Raymond and Beverly Sackler Faculty of Exact Sciences, Tel Aviv University, Tel Aviv, Israel
- The Sagol School of Neuroscience, Tel Aviv University, Tel Aviv, Israel
| | - Yoram Cohen
- School of Chemistry, The Raymond and Beverly Sackler Faculty of Exact Sciences, Tel Aviv University, Tel Aviv, Israel
- The Sagol School of Neuroscience, Tel Aviv University, Tel Aviv, Israel
| |
Collapse
|
22
|
Saliani A, Zaimi A, Nami H, Duval T, Stikov N, Cohen-Adad J. Construction of a rat spinal cord atlas of axon morphometry. Neuroimage 2019; 202:116156. [PMID: 31491525 DOI: 10.1016/j.neuroimage.2019.116156] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2019] [Revised: 08/26/2019] [Accepted: 09/02/2019] [Indexed: 12/27/2022] Open
Abstract
Atlases of the central nervous system are essential for understanding the pathophysiology of neurological diseases, which remains one of the greatest challenges in neuroscience research today. These atlases provide insight into the underlying white matter microstructure and have been created from a variety of animal models, including rats. Although existing atlases of the rat spinal cord provide some details of axon microstructure, there is currently no histological dataset that quantifies axon morphometry exhaustively in the entire spinal cord. In this study, we created the first comprehensive rat spinal cord atlas of the white matter microstructure with quantifiable axon and myelin morphometrics. Using full-slice scanning electron microscopy images and state-of-the-art segmentation algorithms, we generated an atlas of microstructural metrics such as axon diameter, axonal density and g-ratio. After registering the Watson spinal cord white matter atlas to our template, we computed statistics across metrics, spinal levels and tracts. We notably found that g-ratio is relatively constant, whereas axon diameter showed the greatest variation. The atlas, data and full analysis code are freely available at: https://github.com/neuropoly/atlas-rat.
Collapse
Affiliation(s)
- Ariane Saliani
- NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique Montreal, Montreal, QC, Canada.
| | - Aldo Zaimi
- NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique Montreal, Montreal, QC, Canada
| | - Harris Nami
- NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique Montreal, Montreal, QC, Canada
| | - Tanguy Duval
- NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique Montreal, Montreal, QC, Canada
| | - Nikola Stikov
- NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique Montreal, Montreal, QC, Canada; Montreal Heart Institute, Montreal, QC, Canada
| | - Julien Cohen-Adad
- NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique Montreal, Montreal, QC, Canada; Functional Neuroimaging Unit, CRIUGM, Université de Montréal, Montréal, QC, Canada.
| |
Collapse
|
23
|
Cadotte DW, Akbar MA, Fehlings MG, Stroman PW, Cohen-Adad J. What Has Been Learned from Magnetic Resonance Imaging Examination of the Injured Human Spinal Cord: A Canadian Perspective. J Neurotrauma 2019; 35:1942-1957. [PMID: 30074873 DOI: 10.1089/neu.2018.5903] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Abstract
Magnetic resonance imaging (MRI) has transformed the way surgeons and researchers study and treat spinal cord injury. In this narrative review, we explore the historical context of imaging the human spinal cord and describe how MRI has evolved from providing the first visualization of the human spinal cord in the 1980s to a remarkable set of imaging tools today. The article focuses in particular on the role of Canadian researchers to this field. We begin by outlining the clinical context of traumatic injury to the human spinal cord and describe why current MRI standards fall short when it comes to treating this disabling condition. Parts 2 and 3 of this work explore an exciting and dramatic shift in the use of MRI technology to aid in our understanding and treatment of traumatic injury to the spinal cord. We explore the use of functional imaging (part 2) and structural imaging (part 3) and explore how these techniques have evolved, how they are used, and the challenges that we face for continued refinement and application to patients who live with the neurological and functional deficits caused by injury to the delicate spinal cord.
Collapse
Affiliation(s)
- David W Cadotte
- 1 University of Calgary Spine Program, Division of Neurosurgery, Department of Clinical Neurosciences, University of Calgary , Foothills Medical Centre, Calgary, Alberta, Canada
| | - M Ali Akbar
- 2 Department of Surgery, Division of Neurosurgery and Spinal Program, Toronto Western Hospital, University of Toronto , Toronto, Ontario, Canada
| | - Michael G Fehlings
- 2 Department of Surgery, Division of Neurosurgery and Spinal Program, Toronto Western Hospital, University of Toronto , Toronto, Ontario, Canada
| | - Patrick W Stroman
- 3 Centre for Neuroscience Studies, Queens University , Kingston, Ontario, Canada
| | - Julien Cohen-Adad
- 4 NeuroPoly Lab, Institute of Biomedical Engineering , Polytechnique Montreal, Montreal, Quebéc, Canada .,5 Functional Neuroimaging Unit, CRIUGM, Université de Montréal , Montreal, Quebéc, Canada
| |
Collapse
|
24
|
Schilling KG, By S, Feiler HR, Box BA, O'Grady KP, Witt A, Landman BA, Smith SA. Diffusion MRI microstructural models in the cervical spinal cord - Application, normative values, and correlations with histological analysis. Neuroimage 2019; 201:116026. [PMID: 31326569 DOI: 10.1016/j.neuroimage.2019.116026] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2019] [Revised: 07/12/2019] [Accepted: 07/16/2019] [Indexed: 12/14/2022] Open
Abstract
Multi-compartment tissue modeling using diffusion magnetic resonance imaging has proven valuable in the brain, offering novel indices sensitive to the tissue microstructural environment in vivo on clinical MRI scanners. However, application, characterization, and validation of these models in the spinal cord remain relatively under-studied. In this study, we apply a diffusion "signal" model (diffusion tensor imaging, DTI) and two commonly implemented "microstructural" models (neurite orientation dispersion and density imaging, NODDI; spherical mean technique, SMT) in the human cervical spinal cord of twenty-one healthy controls. We first provide normative values of DTI, SMT, and NODDI indices in a number of white matter ascending and descending pathways, as well as various gray matter regions. We then aim to validate the sensitivity and specificity of these diffusion-derived contrasts by relating these measures to indices of the tissue microenvironment provided by a histological template. We find that DTI indices are sensitive to a number of microstructural features, but lack specificity. The microstructural models also show sensitivity to a number of microstructure features; however, they do not capture the specific microstructural features explicitly modelled. Although often regarded as a simple extension of the brain in the central nervous system, it may be necessary to re-envision, or specifically adapt, diffusion microstructural models for application to the human spinal cord with clinically feasible acquisitions - specifically, adjusting, adapting, and re-validating the modeling as it relates to both theory (i.e. relevant biology, assumptions, and signal regimes) and parameter estimation (for example challenges of acquisition, artifacts, and processing).
Collapse
Affiliation(s)
- Kurt G Schilling
- Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, TN, USA.
| | - Samantha By
- Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, TN, USA
| | - Haley R Feiler
- Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, TN, USA
| | - Bailey A Box
- Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, TN, USA
| | - Kristin P O'Grady
- Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, TN, USA; Department of Radiology and Radiological Sciences, Vanderbilt University Medical Center, Nashville, TN, USA
| | - Atlee Witt
- Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, TN, USA
| | - Bennett A Landman
- Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, TN, USA; Department of Biomedical Engineering, Vanderbilt University, Nashville, TN, USA; Department of Electrical Engineering, Vanderbilt University, Nashville, TN, USA
| | - Seth A Smith
- Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, TN, USA; Department of Biomedical Engineering, Vanderbilt University, Nashville, TN, USA; Department of Radiology and Radiological Sciences, Vanderbilt University Medical Center, Nashville, TN, USA
| |
Collapse
|
25
|
Lee HH, Yaros K, Veraart J, Pathan JL, Liang FX, Kim SG, Novikov DS, Fieremans E. Along-axon diameter variation and axonal orientation dispersion revealed with 3D electron microscopy: implications for quantifying brain white matter microstructure with histology and diffusion MRI. Brain Struct Funct 2019; 224:1469-1488. [PMID: 30790073 PMCID: PMC6510616 DOI: 10.1007/s00429-019-01844-6] [Citation(s) in RCA: 57] [Impact Index Per Article: 11.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2018] [Accepted: 02/01/2019] [Indexed: 10/27/2022]
Abstract
Tissue microstructure modeling of diffusion MRI signal is an active research area striving to bridge the gap between macroscopic MRI resolution and cellular-level tissue architecture. Such modeling in neuronal tissue relies on a number of assumptions about the microstructural features of axonal fiber bundles, such as the axonal shape (e.g., perfect cylinders) and the fiber orientation dispersion. However, these assumptions have not yet been validated by sufficiently high-resolution 3-dimensional histology. Here, we reconstructed sequential scanning electron microscopy images in mouse brain corpus callosum, and introduced a random-walker (RaW)-based algorithm to rapidly segment individual intra-axonal spaces and myelin sheaths of myelinated axons. Confirmed by a segmentation based on human annotations initiated with conventional machine-learning-based carving, our semi-automatic algorithm is reliable and less time-consuming. Based on the segmentation, we calculated MRI-relevant estimates of size-related parameters (inner axonal diameter, its distribution, along-axon variation, and myelin g-ratio), and orientation-related parameters (fiber orientation distribution and its rotational invariants; dispersion angle). The reported dispersion angle is consistent with previous 2-dimensional histology studies and diffusion MRI measurements, while the reported diameter exceeds those in other mouse brain studies. Furthermore, we calculated how these quantities would evolve in actual diffusion MRI experiments as a function of diffusion time, thereby providing a coarse-graining window on the microstructure, and showed that the orientation-related metrics have negligible diffusion time-dependence over clinical and pre-clinical diffusion time ranges. However, the MRI-measured inner axonal diameters, dominated by the widest cross sections, effectively decrease with diffusion time by ~ 17% due to the coarse-graining over axonal caliber variations. Furthermore, our 3d measurement showed that there is significant variation of the diameter along the axon. Hence, fiber orientation dispersion estimated from MRI should be relatively stable, while the "apparent" inner axonal diameters are sensitive to experimental settings, and cannot be modeled by perfectly cylindrical axons.
Collapse
Affiliation(s)
- Hong-Hsi Lee
- Department of Radiology, Center for Biomedical Imaging, New York University School of Medicine, 660 First Avenue, New York, NY, 10016, USA.
- Center for Advanced Imaging Innovation and Research (CAI2R), New York University School of Medicine, 660 First Avenue, New York, NY, 10016, USA.
| | - Katarina Yaros
- Department of Radiology, Center for Biomedical Imaging, New York University School of Medicine, 660 First Avenue, New York, NY, 10016, USA
| | - Jelle Veraart
- Department of Radiology, Center for Biomedical Imaging, New York University School of Medicine, 660 First Avenue, New York, NY, 10016, USA
| | - Jasmine L Pathan
- Department of Radiology, Center for Biomedical Imaging, New York University School of Medicine, 660 First Avenue, New York, NY, 10016, USA
| | - Feng-Xia Liang
- Department of Cell Biology and Microscopy Core, New York University School of Medicine, 540 First Avenue, New York, NY, 10016, USA
| | - Sungheon G Kim
- Department of Radiology, Center for Biomedical Imaging, New York University School of Medicine, 660 First Avenue, New York, NY, 10016, USA
- Center for Advanced Imaging Innovation and Research (CAI2R), New York University School of Medicine, 660 First Avenue, New York, NY, 10016, USA
| | - Dmitry S Novikov
- Department of Radiology, Center for Biomedical Imaging, New York University School of Medicine, 660 First Avenue, New York, NY, 10016, USA
- Center for Advanced Imaging Innovation and Research (CAI2R), New York University School of Medicine, 660 First Avenue, New York, NY, 10016, USA
| | - Els Fieremans
- Department of Radiology, Center for Biomedical Imaging, New York University School of Medicine, 660 First Avenue, New York, NY, 10016, USA
- Center for Advanced Imaging Innovation and Research (CAI2R), New York University School of Medicine, 660 First Avenue, New York, NY, 10016, USA
| |
Collapse
|
26
|
Novikov DS, Fieremans E, Jespersen SN, Kiselev VG. Quantifying brain microstructure with diffusion MRI: Theory and parameter estimation. NMR IN BIOMEDICINE 2019; 32:e3998. [PMID: 30321478 PMCID: PMC6481929 DOI: 10.1002/nbm.3998] [Citation(s) in RCA: 248] [Impact Index Per Article: 49.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/13/2016] [Revised: 06/11/2018] [Accepted: 06/28/2018] [Indexed: 05/18/2023]
Abstract
We review, systematize and discuss models of diffusion in neuronal tissue, by putting them into an overarching physical context of coarse-graining over an increasing diffusion length scale. From this perspective, we view research on quantifying brain microstructure as occurring along three major avenues. The first avenue focusses on transient, or time-dependent, effects in diffusion. These effects signify the gradual coarse-graining of tissue structure, which occurs qualitatively differently in different brain tissue compartments. We show that transient effects contain information about the relevant length scales for neuronal tissue, such as the packing correlation length for neuronal fibers, as well as the degree of structural disorder along the neurites. The second avenue corresponds to the long-time limit, when the observed signal can be approximated as a sum of multiple nonexchanging anisotropic Gaussian components. Here, the challenge lies in parameter estimation and in resolving its hidden degeneracies. The third avenue employs multiple diffusion encoding techniques, able to access information not contained in the conventional diffusion propagator. We conclude with our outlook on future directions that could open exciting possibilities for designing quantitative markers of tissue physiology and pathology, based on methods of studying mesoscopic transport in disordered systems.
Collapse
Affiliation(s)
- Dmitry S. Novikov
- Center for Biomedical Imaging, Department of Radiology, NYU School of Medicine, New York, NY, USA
| | - Els Fieremans
- Center for Biomedical Imaging, Department of Radiology, NYU School of Medicine, New York, NY, USA
| | - Sune N. Jespersen
- CFIN/MINDLab, Department of Clinical Medicine and Department of Physics and Astronomy, Aarhus University, Aarhus, Denmark
| | - Valerij G. Kiselev
- Medical Physics, Deptartment of Radiology, Faculty of Medicine, University of Freiburg, Germany
| |
Collapse
|
27
|
Alexander DC, Dyrby TB, Nilsson M, Zhang H. Imaging brain microstructure with diffusion MRI: practicality and applications. NMR IN BIOMEDICINE 2019; 32:e3841. [PMID: 29193413 DOI: 10.1002/nbm.3841] [Citation(s) in RCA: 204] [Impact Index Per Article: 40.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/28/2016] [Revised: 07/09/2017] [Accepted: 09/11/2017] [Indexed: 05/22/2023]
Abstract
This article gives an overview of microstructure imaging of the brain with diffusion MRI and reviews the state of the art. The microstructure-imaging paradigm aims to estimate and map microscopic properties of tissue using a model that links these properties to the voxel scale MR signal. Imaging techniques of this type are just starting to make the transition from the technical research domain to wide application in biomedical studies. We focus here on the practicalities of both implementing such techniques and using them in applications. Specifically, the article summarizes the relevant aspects of brain microanatomy and the range of diffusion-weighted MR measurements that provide sensitivity to them. It then reviews the evolution of mathematical and computational models that relate the diffusion MR signal to brain tissue microstructure, as well as the expanding areas of application. Next we focus on practicalities of designing a working microstructure imaging technique: model selection, experiment design, parameter estimation, validation, and the pipeline of development of this class of technique. The article concludes with some future perspectives on opportunities in this topic and expectations on how the field will evolve in the short-to-medium term.
Collapse
Affiliation(s)
- Daniel C Alexander
- Centre for Medical Image Computing (CMIC), Department of Computer Science, UCL (University College London), Gower Street, London, UK
| | - Tim B Dyrby
- Danish Research Centre for Magnetic Resonance, Center for Functional and Diagnostic Imaging and Research, Copenhagen University Hospital Hvidovre, Hvidovre, Denmark
- Department of Applied Mathematics and Computer Science, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Markus Nilsson
- Clinical Sciences Lund, Department of Radiology, Lund University, Lund, Sweden
| | - Hui Zhang
- Centre for Medical Image Computing (CMIC), Department of Computer Science, UCL (University College London), Gower Street, London, UK
| |
Collapse
|
28
|
Berman S, Filo S, Mezer AA. Modeling conduction delays in the corpus callosum using MRI-measured g-ratio. Neuroimage 2019; 195:128-139. [PMID: 30910729 DOI: 10.1016/j.neuroimage.2019.03.025] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/26/2018] [Revised: 03/12/2019] [Accepted: 03/12/2019] [Indexed: 11/26/2022] Open
Abstract
Conduction of action potentials along myelinated axons is affected by their structural features, such as the axonal g-ratio, the ratio between the inner and outer diameters of the myelin sheath surrounding the axon. The effect of g-ratio variance on conduction properties has been quantitatively evaluated using single-axon models. It has recently become possible to estimate a g-ratio weighted measurement in vivo using quantitative MRI. Nevertheless, it is still unclear whether the variance in the g-ratio in the healthy human brain leads to significant differences in conduction velocity. In this work we tested whether the g-ratio MRI measurement can be used to predict conduction delays in the corpus callosum. We present a comprehensive framework in which the structural properties of fibers (i.e. length and g-ratio, measured using MRI), are incorporated in a biophysical model of axon conduction, to model conduction delays of long-range white matter fibers. We applied this framework to the corpus callosum, and found conduction delay estimates that are compatible with previously estimated values of conduction delays. We account for the variance in the velocity given the axon diameter distribution in the splenium, mid-body and genu, to further compare the fibers within the corpus callosum. Conduction delays have been suggested to increase with age. Therefore, we investigated whether there are differences in the g-ratio and the fiber length between young and old adults, and whether this leads to a difference in conduction speed and delays. We found very small differences between the predicted delays of the two groups in the motor fibers of the corpus callosum. We also found that the motor fibers of the corpus callosum have the fastest conduction estimates. Using the axon diameter distributions, we found that the occipital fibers have the slowest estimations, while the frontal and motor fiber tracts have similar estimates. Our study provides a framework for predicting conduction latencies in vivo. The framework could have major implications for future studies of white matter diseases and large range network computations. Our results highlight the need for improving additional in vivo measurements of white matter microstructure.
Collapse
Affiliation(s)
- S Berman
- Edmond & Lily Safra Center for Brain Sciences, The Hebrew University of Jerusalem, Jerusalem, Israel.
| | - S Filo
- Edmond & Lily Safra Center for Brain Sciences, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - A A Mezer
- Edmond & Lily Safra Center for Brain Sciences, The Hebrew University of Jerusalem, Jerusalem, Israel
| |
Collapse
|
29
|
Lee HH, Fieremans E, Novikov DS. What dominates the time dependence of diffusion transverse to axons: Intra- or extra-axonal water? Neuroimage 2018; 182:500-510. [PMID: 29253652 PMCID: PMC6004237 DOI: 10.1016/j.neuroimage.2017.12.038] [Citation(s) in RCA: 54] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2017] [Revised: 11/18/2017] [Accepted: 12/13/2017] [Indexed: 12/31/2022] Open
Abstract
Brownian motion of water molecules provides an essential length scale, the diffusion length, commensurate with cell dimensions in biological tissues. Measuring the diffusion coefficient as a function of diffusion time makes in vivo diffusion MRI uniquely sensitive to the cellular features about three orders of magnitude below imaging resolution. However, there is a longstanding debate, regarding which contribution - intra- or extra-cellular - is more relevant in the overall time-dependence of the MRI-derived diffusion metrics. Here we resolve this debate in the human brain white matter. By varying not just the diffusion time, but also the gradient pulse duration of a standard diffusion MRI sequence, we identify a functional form of the measured time-dependent diffusion coefficient transverse to white matter tracts in 10 healthy volunteers. This specific functional form is shown to originate from the extra-axonal space, and provides estimates of the fiber packing correlation length for axons in a bundle. Our results offer a metric for the outer axonal diameter, a promising candidate marker for demyelination in neurodegenerative diseases. From the methodological perspective, our analysis demonstrates how competing models, which describe different physics yet interpolate standard measurements equally well, can be distinguished based on their prediction for an independent "orthogonal" measurement.
Collapse
Affiliation(s)
- Hong-Hsi Lee
- Center for Biomedical Imaging, Department of Radiology, NYU School of Medicine, New York, NY 10016, United States.
| | - Els Fieremans
- Center for Biomedical Imaging, Department of Radiology, NYU School of Medicine, New York, NY 10016, United States
| | - Dmitry S Novikov
- Center for Biomedical Imaging, Department of Radiology, NYU School of Medicine, New York, NY 10016, United States
| |
Collapse
|
30
|
Duval T, Saliani A, Nami H, Nanci A, Stikov N, Leblond H, Cohen-Adad J. Axons morphometry in the human spinal cord. Neuroimage 2018; 185:119-128. [PMID: 30326296 DOI: 10.1016/j.neuroimage.2018.10.033] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2018] [Revised: 10/05/2018] [Accepted: 10/10/2018] [Indexed: 12/15/2022] Open
Abstract
Due to the technical challenges of large-scale microscopy and analysis, to date only limited knowledge has been made available about axon morphometry (diameter, shape, myelin thickness, volume fraction), thereby limiting our understanding of neuronal microstructure and slowing down research on neurodegenerative pathologies. This study addresses this knowledge gap by establishing a state-of-the-art acquisition and analysis framework for mapping axon morphometry, and providing the first comprehensive mapping of axon morphometry in the human spinal cord. We dissected, fixed and stained a human spinal cord with osmium tetroxide, and used a scanning electron microscope to image the entirety of 23 axial slices, covering C1 to L5 spinal levels. An automatic method based on deep learning was then used to segment each axon and myelin sheath to produce maps of axon morphometry. These maps were then registered to a standard spinal cord magnetic resonance imaging (MRI) template. Between 500,000 (lumbar) and 1 million (cervical) myelinated axons were segmented at each level of this human spinal cord. Morphometric features show a large disparity between tracts, but high right-left symmetry. Our results suggest a modality-based organization of the dorsal column in the human, as it has been observed in the rat. The generated axon morphometry template is publicly available at https://osf.io/8k7jr/ and could be used as a reference for quantitative MRI studies. The proposed framework for axon morphometry mapping could be extended to other parts of the central or peripheral nervous system that exhibit coherently-oriented axons.
Collapse
Affiliation(s)
- Tanguy Duval
- NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique Montreal, Montreal, QC, Canada
| | - Ariane Saliani
- NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique Montreal, Montreal, QC, Canada
| | - Harris Nami
- NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique Montreal, Montreal, QC, Canada
| | - Antonio Nanci
- Laboratory for the Study of Calcified Tissues and Biomaterials, Faculty of Dental Medicine, University of Montreal, Montreal, QC, Canada
| | - Nikola Stikov
- NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique Montreal, Montreal, QC, Canada; Montreal Heart Institute, Montreal, QC, Canada
| | - Hugues Leblond
- Anatomy department, Université du Québec à Trois-Rivières, Trois-Rivières, QC, Canada
| | - Julien Cohen-Adad
- NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique Montreal, Montreal, QC, Canada; Functional Neuroimaging Unit, CRIUGM, Université de Montréal, Montréal, QC, Canada.
| |
Collapse
|
31
|
Grussu F, Ianuş A, Tur C, Prados F, Schneider T, Kaden E, Ourselin S, Drobnjak I, Zhang H, Alexander DC, Gandini Wheeler-Kingshott CAM. Relevance of time-dependence for clinically viable diffusion imaging of the spinal cord. Magn Reson Med 2018; 81:1247-1264. [PMID: 30229564 PMCID: PMC6586052 DOI: 10.1002/mrm.27463] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/09/2018] [Revised: 06/28/2018] [Accepted: 06/29/2018] [Indexed: 12/17/2022]
Abstract
Purpose Time‐dependence is a key feature of the diffusion‐weighted (DW) signal, knowledge of which informs biophysical modelling. Here, we study time‐dependence in the human spinal cord, as its axonal structure is specific and different from the brain. Methods We run Monte Carlo simulations using a synthetic model of spinal cord white matter (WM) (large axons), and of brain WM (smaller axons). Furthermore, we study clinically feasible multi‐shell DW scans of the cervical spinal cord (b = 0; b = 711 s mm−2; b = 2855 s mm−2), obtained using three diffusion times (Δ of 29, 52 and 76 ms) from three volunteers. Results Both intra‐/extra‐axonal perpendicular diffusivities and kurtosis excess show time‐dependence in our synthetic spinal cord model. This time‐dependence is reflected mostly in the intra‐axonal perpendicular DW signal, which also exhibits strong decay, unlike our brain model. Time‐dependence of the total DW signal appears detectable in the presence of noise in our synthetic spinal cord model, but not in the brain. In WM in vivo, we observe time‐dependent macroscopic and microscopic diffusivities and diffusion kurtosis, NODDI and two‐compartment SMT metrics. Accounting for large axon calibers improves fitting of multi‐compartment models to a minor extent. Conclusions Time‐dependence of clinically viable DW MRI metrics can be detected in vivo in spinal cord WM, thus providing new opportunities for the non‐invasive estimation of microstructural properties. The time‐dependence of the perpendicular DW signal may feature strong intra‐axonal contributions due to large spinal axon caliber. Hence, a popular model known as “stick” (zero‐radius cylinder) may be sub‐optimal to describe signals from the largest spinal axons.
Collapse
Affiliation(s)
- Francesco Grussu
- Queen Square MS Centre, UCL Institute of Neurology, Faculty of Brain Sciences, University College London, London, United Kingdom.,Centre for Medical Image Computing, Department of Computer Science, University College London, London, United Kingdom
| | - Andrada Ianuş
- Centre for Medical Image Computing, Department of Computer Science, University College London, London, United Kingdom.,Champalimaud Centre for the Unknown, Champalimaud Foundation, Lisbon, Portugal
| | - Carmen Tur
- Queen Square MS Centre, UCL Institute of Neurology, Faculty of Brain Sciences, University College London, London, United Kingdom
| | - Ferran Prados
- Queen Square MS Centre, UCL Institute of Neurology, Faculty of Brain Sciences, University College London, London, United Kingdom.,Centre for Medical Image Computing, Department of Medical Physics and Biomedical Engineering, University College London, London, United Kingdom
| | | | - Enrico Kaden
- Centre for Medical Image Computing, Department of Computer Science, University College London, London, United Kingdom
| | - Sébastien Ourselin
- Centre for Medical Image Computing, Department of Medical Physics and Biomedical Engineering, University College London, London, United Kingdom
| | - Ivana Drobnjak
- Centre for Medical Image Computing, Department of Computer Science, University College London, London, United Kingdom
| | - Hui Zhang
- Centre for Medical Image Computing, Department of Computer Science, University College London, London, United Kingdom
| | - Daniel C Alexander
- Centre for Medical Image Computing, Department of Computer Science, University College London, London, United Kingdom.,Clinical Imaging Research Centre, National University of Singapore, Singapore, Singapore
| | - Claudia A M Gandini Wheeler-Kingshott
- Queen Square MS Centre, UCL Institute of Neurology, Faculty of Brain Sciences, University College London, London, United Kingdom.,Brain MRI 3T Research Centre, C. Mondino National Neurological Institute, Pavia, Italy.,Department of Brain and Behavioural Sciences, University of Pavia, Pavia, Italy
| |
Collapse
|
32
|
Frostell A, Mattsson P, Svensson M. Guiding Device for Precision Grafting of Peripheral Nerves in Complete Thoracic Spinal Cord Injury: Design and Sizing for Clinical Trial. Front Neurol 2018; 9:356. [PMID: 29872421 PMCID: PMC5972322 DOI: 10.3389/fneur.2018.00356] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2017] [Accepted: 05/01/2018] [Indexed: 01/08/2023] Open
Abstract
Background In an effort to translate preclinical success in achieving spinal cord regeneration through peripheral nerve grafts, this study details the design and sizing of a guiding device for precision grafting of peripheral nerves for use in a clinical trial in complete (AIS-A) thoracic spinal cord injury (SCI). The device’s design and sizing are compared to a simulation of human spinal cord sizes based on the best available data. Methods Spinal cord segmental sizes were generated by computer simulation based on data from a meta-analysis recently published by our group. Thoracic segments T2–T12 were plotted, and seven elliptical shapes were positioned across the center of the distribution of sizes. Geometrical measures of error-of-fit were calculated. CAD modeling was used to create cranial and caudal interfaces for the human spinal cord, aiming to guide descending white matter tracts to gray matter at the caudal end of the device and ascending white matter tracts to gray matter at the cranial end of the device. The interfaces were compared qualitatively to the simulated spinal cord sizes and gray-to-white matter delineations. Results The mean error-of-fit comparing simulated spinal cord segments T2–T12 to the best elliptical shape was 0.41 and 0.36 mm, and the 95th percentile was found at 1.3 and 0.98 mm for transverse and anteroposterior diameter, respectively. A guiding device design was reached for capturing the majority of corticospinal axons at the cranial end of the device and guiding them obliquely to gray matter at the caudal end of the device. Based on qualitative comparison, the vast majority of spinal cord sizes generated indicate an excellent fit to the device’s interfaces. Conclusion A set of SCI guiding devices of seven sizes can cover the variability of human thoracic spinal cord segments T2–T12 with an acceptable error-of-fit for the elliptical shape as well as guiding channels. The computational framework developed can be used with other medical technologies involving the human spinal cord where exact sizes and positioning are of importance.
Collapse
Affiliation(s)
- Arvid Frostell
- Department of Clinical Neuroscience, Karolinska Institutet, Stockholm, Sweden
| | - Per Mattsson
- Department of Clinical Neuroscience, Karolinska Institutet, Stockholm, Sweden.,Department of Breast, Endocrine and Sarcoma Tumors, Karolinska University Hospital, Stockholm, Sweden
| | - Mikael Svensson
- Department of Clinical Neuroscience, Karolinska Institutet, Stockholm, Sweden.,Department of Neurosurgery, Karolinska University Hospital, Stockholm, Sweden
| |
Collapse
|
33
|
Duval T, Smith V, Stikov N, Klawiter EC, Cohen-Adad J. Scan-rescan of axcaliber, macromolecular tissue volume, and g-ratio in the spinal cord. Magn Reson Med 2018; 79:2759-2765. [PMID: 28994487 PMCID: PMC5821542 DOI: 10.1002/mrm.26945] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2017] [Revised: 09/01/2017] [Accepted: 09/03/2017] [Indexed: 01/25/2023]
Abstract
PURPOSE Recent MRI techniques have been introduced that can extract microstructural information in the white matter, such as the density or macromolecular content. Translating quantitative MRI to the clinic raises many challenges in terms of acquisition strategy, modeling of the MRI signal, artifact corrections, and metric extraction (template registration and partial volume effects). In this work, we investigated the scan-rescan repeatability of several quantitative MRI techniques in the human spinal cord. METHODS AxCaliber metrics, macromolecular tissue volume, and the fiber g-ratio were estimated in the spinal cord of eight healthy subjects, scanned and rescanned the same day in two different sessions. RESULTS Scan-rescan repeatability deviation was 3% for all metrics, in average in the white matter of all subjects. Intraclass correlation coefficient was up to 0.9. A three-way analysis of variance showed significant effects of white matter pathway, laterality, and subject. CONCLUSION The present study suggests that quantitative MRI gives stable measurements of white matter microstructure in the spinal cord of healthy subjects. Our findings remain to be evaluated in diseased populations. Magn Reson Med 79:2759-2765, 2018. © 2017 International Society for Magnetic Resonance in Medicine.
Collapse
Affiliation(s)
- Tanguy Duval
- NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique Montreal, Montreal, QC, Canada
| | - Victoria Smith
- Department of Neurology, Massachusetts General Hospital, Harvard Medical School, Boston, MA, United States
| | - Nikola Stikov
- NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique Montreal, Montreal, QC, Canada
- Montreal Heart Institute, Montreal, QC, Canada
| | - Eric C. Klawiter
- Department of Neurology, Massachusetts General Hospital, Harvard Medical School, Boston, MA, United States
| | - Julien Cohen-Adad
- NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique Montreal, Montreal, QC, Canada
- Functional Neuroimaging Unit, CRIUGM, Université de Montréal, Montréal, QC, Canada
| |
Collapse
|
34
|
Cohen-Adad J. Microstructural imaging in the spinal cord and validation strategies. Neuroimage 2018; 182:169-183. [PMID: 29635029 DOI: 10.1016/j.neuroimage.2018.04.009] [Citation(s) in RCA: 31] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2017] [Revised: 03/02/2018] [Accepted: 04/06/2018] [Indexed: 12/13/2022] Open
Abstract
In vivo histology using magnetic resonance imaging (MRI) is a newly emerging research field that aims to non-invasively characterize tissue microstructure. The implications of in vivo histology are many, from discovering novel biomarkers to studying human development, to providing tools for disease diagnosis and monitoring the effects of novel treatments on tissue. This review focuses on quantitative MRI (qMRI) techniques that are used to map spinal cord microstructure. Opening with a rationale for non-invasive imaging of the spinal cord, this article continues with a brief overview of the existing MRI techniques for axon and myelin imaging, followed by the specific challenges and potential solutions for acquiring and processing such data. The final part of this review focuses on histological validation, with suggested tissue preparation, acquisition and processing protocols for large-scale microscopy.
Collapse
Affiliation(s)
- J Cohen-Adad
- NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique Montreal, Montreal, QC, Canada; Functional Neuroimaging Unit, CRIUGM, Université de Montréal, Montreal, QC, Canada.
| |
Collapse
|
35
|
By S, Xu J, Box BA, Bagnato FR, Smith SA. Multi-compartmental diffusion characterization of the human cervical spinal cord in vivo using the spherical mean technique. NMR IN BIOMEDICINE 2018; 31:e3894. [PMID: 29388719 PMCID: PMC5854548 DOI: 10.1002/nbm.3894] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/28/2017] [Revised: 12/12/2017] [Accepted: 12/16/2017] [Indexed: 05/02/2023]
Abstract
The purpose of this work was to evaluate the feasibility and reproducibility of the spherical mean technique (SMT), a multi-compartmental diffusion model, in the spinal cord of healthy controls, and to assess its ability to improve spinal cord characterization in multiple sclerosis (MS) patients at 3 T. SMT was applied in the cervical spinal cord of eight controls and six relapsing-remitting MS patients. SMT provides an elegant framework to model the apparent axonal volume fraction vax , intrinsic diffusivity Dax , and extra-axonal transverse diffusivity Dex_perp (which is estimated as a function of vax and Dax ) without confounds related to complex fiber orientation distribution that reside in diffusion MRI modeling. SMT's reproducibility was assessed with two different scans within a month, and SMT-derived indices in healthy and MS cohorts were compared. The influence of acquisition scheme on SMT was also evaluated. SMT's vax , Dax , and Dex_perp measurements all showed high reproducibility. A decrease in vax was observed at the site of lesions and normal appearing white matter (p < 0.05), and trends towards a decreased Dax and increased Dex_perp were seen. Importantly, a twofold reduction in acquisition yielded similarly high accuracy with SMT. SMT provides a fast, reproducible, and accurate method to improve characterization of the cervical spinal cord, and may have clinical potential for MS patients.
Collapse
Affiliation(s)
- Samantha By
- Department of Biomedical Engineering, Vanderbilt University, Nashville, TN, USA
- Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, TN, USA
- Philips Healthcare, Gainesville, FL, USA
| | - Junzhong Xu
- Department of Biomedical Engineering, Vanderbilt University, Nashville, TN, USA
- Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, TN, USA
- Department of Radiology and Radiological Sciences, Vanderbilt University Medical Center, Nashville, TN, USA
| | - Bailey A. Box
- Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, TN, USA
| | - Francesca R. Bagnato
- Department of Neurology, Vanderbilt University Medical Center, Nashville, TN, USA
| | - Seth A. Smith
- Department of Biomedical Engineering, Vanderbilt University, Nashville, TN, USA
- Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, TN, USA
- Department of Radiology and Radiological Sciences, Vanderbilt University Medical Center, Nashville, TN, USA
| |
Collapse
|
36
|
Duyn JH. Studying brain microstructure with magnetic susceptibility contrast at high-field. Neuroimage 2018; 168:152-161. [PMID: 28242317 PMCID: PMC5569005 DOI: 10.1016/j.neuroimage.2017.02.046] [Citation(s) in RCA: 31] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2016] [Revised: 02/03/2017] [Accepted: 02/16/2017] [Indexed: 12/14/2022] Open
Abstract
A rapidly developing application of high field MRI is the study of brain anatomy and function with contrast based on the magnetic susceptibility of tissues. To study the subtle variations in susceptibility contrast between tissues and with changes in brain activity, dedicated scan techniques such as susceptibility-weighted MRI and blood-oxygen level dependent functional MRI have been developed. Particularly strong susceptibility contrast has been observed with systems that operate at 7T and above, and their recent widespread use has led to an improved understanding of contributing sources and mechanisms. To interpret magnetic susceptibility contrast, analysis approaches have been developed with the goal of extracting measures that report on local tissue magnetic susceptibility, a physical quantity that, under certain conditions, allows estimation of blood oxygenation, local tissue iron content, and quantification of its changes with disease. Interestingly, high field studies have also brought to light that not only the makeup of tissues affects MRI susceptibility contrast, but that also a tissue's sub-voxel structure at scales all the way down to the molecular level plays an important role as well. In this review, various ways will be discussed by which sub-voxel structure can affect the MRI signal in general, and magnetic susceptibility in particular, sometimes in a complex fashion. In the light of this complexity, it appears likely that accurate, brain-wide quantification of iron will require the combination of multiple contrasts that may include diffusion and magnetization transfer information with susceptibility-weighted contrast. On the other hand, this complexity also offers opportunities to use magnetic susceptibility contrast to inform about specific microstructural aspects of brain tissue. Details and several examples will be presented in this review.
Collapse
Affiliation(s)
- Jeff H Duyn
- Advanced MRI Section, Laboratory of Functional and Molecular Imaging, National Institutes of Neurological Disorders and Stroke, National Institutes of Health, Bethesda, MD 20892, USA
| |
Collapse
|
37
|
Kafali SG, Çukur T, Saritas EU. Phase-correcting non-local means filtering for diffusion-weighted imaging of the spinal cord. Magn Reson Med 2018; 80:1020-1035. [PMID: 29427379 DOI: 10.1002/mrm.27105] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2017] [Revised: 12/26/2017] [Accepted: 01/02/2018] [Indexed: 11/09/2022]
Affiliation(s)
- Sevgi Gokce Kafali
- Department of Electrical and Electronics Engineering, Bilkent University, Ankara, Turkey.,National Magnetic Resonance Research Center (UMRAM), Bilkent University, Ankara, Turkey
| | - Tolga Çukur
- Department of Electrical and Electronics Engineering, Bilkent University, Ankara, Turkey.,National Magnetic Resonance Research Center (UMRAM), Bilkent University, Ankara, Turkey.,Neuroscience Program, Sabuncu Brain Research Center, Bilkent University, Ankara, Turkey
| | - Emine Ulku Saritas
- Department of Electrical and Electronics Engineering, Bilkent University, Ankara, Turkey.,National Magnetic Resonance Research Center (UMRAM), Bilkent University, Ankara, Turkey.,Neuroscience Program, Sabuncu Brain Research Center, Bilkent University, Ankara, Turkey
| |
Collapse
|
38
|
Fan Q, Nummenmaa A, Wichtmann B, Witzel T, Mekkaoui C, Schneider W, Wald LL, Huang SY. Validation of diffusion MRI estimates of compartment size and volume fraction in a biomimetic brain phantom using a human MRI scanner with 300 mT/m maximum gradient strength. Neuroimage 2018; 182:469-478. [PMID: 29337276 DOI: 10.1016/j.neuroimage.2018.01.004] [Citation(s) in RCA: 26] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2017] [Revised: 12/08/2017] [Accepted: 01/03/2018] [Indexed: 10/18/2022] Open
Abstract
Diffusion microstructural imaging techniques have attracted great interest in the last decade due to their ability to quantify axon diameter and volume fraction in healthy and diseased human white matter. The estimates of compartment size and volume fraction continue to be debated, in part due to the lack of a gold standard for validation and quality control. In this work, we validate diffusion MRI estimates of compartment size and volume fraction using a novel textile axon ("taxon") phantom constructed from hollow polypropylene yarns with distinct intra- and extra-taxonal compartments to mimic white matter in the brain. We acquired a comprehensive set of diffusion MRI measurements in the phantom using multiple gradient directions, diffusion times and gradient strengths on a human MRI scanner equipped with maximum gradient strength (Gmax) of 300 mT/m. We obtained estimates of compartment size and restricted volume fraction through a straightforward extension of the AxCaliber/ActiveAx frameworks that enables estimation of mean compartment size in fiber bundles of arbitrary orientation. The voxel-wise taxon diameter estimates of 12.2 ± 0.9 μm were close to the manufactured inner diameter of 11.8 ± 1.2 μm with Gmax = 300 mT/m. The estimated restricted volume fraction demonstrated an expected decrease along the length of the fiber bundles in accordance with the known construction of the phantom. When Gmax was restricted to 80 mT/m, the taxon diameter was overestimated, and the estimates for taxon diameter and packing density showed greater uncertainty compared to data with Gmax = 300 mT/m. In conclusion, the compartment size and volume fraction estimates resulting from diffusion measurements on a human scanner were validated against ground truth in a phantom mimicking human white matter, providing confidence that this method can yield accurate estimates of parameters in simplified but realistic microstructural environments. Our work also demonstrates the importance of a biologically analogous phantom that can be applied to validate a variety of diffusion microstructural imaging methods in human scanners and be used for standardization of diffusion MRI protocols for neuroimaging research.
Collapse
Affiliation(s)
- Qiuyun Fan
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Harvard Medical School, Charlestown, MA, United States.
| | - Aapo Nummenmaa
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Harvard Medical School, Charlestown, MA, United States
| | - Barbara Wichtmann
- Computer Assisted Clinical Medicine, Medical Faculty Mannheim, Heidelberg University, Mannheim, Germany; Institute of Clinical Radiology and Nuclear Medicine, University Medical Center Mannheim, Medical Faculty Mannheim, Heidelberg University, Mannheim, Germany
| | - Thomas Witzel
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Harvard Medical School, Charlestown, MA, United States
| | - Choukri Mekkaoui
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Harvard Medical School, Charlestown, MA, United States
| | - Walter Schneider
- Learning Research and Development Center, University of Pittsburgh, Pittsburgh, PA, United States
| | - Lawrence L Wald
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Harvard Medical School, Charlestown, MA, United States; Harvard-MIT Division of Health Sciences and Technology, Massachusetts Institute of Technology, Cambridge, MA, United States
| | - Susie Y Huang
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Harvard Medical School, Charlestown, MA, United States; Harvard-MIT Division of Health Sciences and Technology, Massachusetts Institute of Technology, Cambridge, MA, United States
| |
Collapse
|
39
|
Lévy S, Guertin MC, Khatibi A, Mezer A, Martinu K, Chen JI, Stikov N, Rainville P, Cohen-Adad J. Test-retest reliability of myelin imaging in the human spinal cord: Measurement errors versus region- and aging-induced variations. PLoS One 2018; 13:e0189944. [PMID: 29293550 PMCID: PMC5749716 DOI: 10.1371/journal.pone.0189944] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2017] [Accepted: 12/05/2017] [Indexed: 01/06/2023] Open
Abstract
PURPOSE To implement a statistical framework for assessing the precision of several quantitative MRI metrics sensitive to myelin in the human spinal cord: T1, Magnetization Transfer Ratio (MTR), saturation imposed by an off-resonance pulse (MTsat) and Macromolecular Tissue Volume (MTV). METHODS Thirty-three healthy subjects within two age groups (young, elderly) were scanned at 3T. Among them, 16 underwent the protocol twice to assess repeatability. Statistical reliability indexes such as the Minimal Detectable Change (MDC) were compared across metrics quantified within different cervical levels and white matter (WM) sub-regions. The differences between pathways and age groups were quantified and interpreted in context of the test-retest repeatability of the measurements. RESULTS The MDC was respectively 105.7ms, 2.77%, 0.37% and 4.08% for T1, MTR, MTsat and MTV when quantified over all WM, while the standard-deviation across subjects was 70.5ms, 1.34%, 0.20% and 2.44%. Even though particular WM regions did exhibit significant differences, these differences were on the same order as test-retest errors. No significant difference was found between age groups for all metrics. CONCLUSION While T1-based metrics (T1 and MTV) exhibited better reliability than MT-based measurements (MTR and MTsat), the observed differences between subjects or WM regions were comparable to (and often smaller than) the MDC. This makes it difficult to determine if observed changes are due to variations in myelin content, or simply due to measurement error. Measurement error remains a challenge in spinal cord myelin imaging, but this study provides statistical guidelines to standardize the field and make it possible to conduct large-scale multi-center studies.
Collapse
Affiliation(s)
- Simon Lévy
- NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique Montreal, Montreal, QC, Canada
- Centre de Recherche de l'Institut Universitaire de Gériatrie de Montréal (CRIUGM), Montréal, QC, Canada
| | - Marie-Claude Guertin
- Montreal Health Innovations Coordinating Center (MHICC), Montreal Heart Institute, Montreal, QC, Canada
| | - Ali Khatibi
- Centre de Recherche de l'Institut Universitaire de Gériatrie de Montréal (CRIUGM), Montréal, QC, Canada
- Psychology Department, Bilkent University, Ankara, Turkey
- Interdisciplinary program in Neuroscience, Bilkent University, Ankara, Turkey
- National Magnetic Resonance Research Center (UMRAM), Bilkent University, Ankara, Turkey
| | - Aviv Mezer
- The Edmond and Lily Safra Center for Brain Sciences (ELSC), The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Kristina Martinu
- Centre de Recherche de l'Institut Universitaire de Gériatrie de Montréal (CRIUGM), Montréal, QC, Canada
| | - Jen-I Chen
- Centre de Recherche de l'Institut Universitaire de Gériatrie de Montréal (CRIUGM), Montréal, QC, Canada
- Department of Stomatology, Faculty of Dentistry, Université de Montréal, Montreal, QC, Canada
| | - Nikola Stikov
- NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique Montreal, Montreal, QC, Canada
- Montreal Heart Institute, Montreal, QC, Canada
| | - Pierre Rainville
- Centre de Recherche de l'Institut Universitaire de Gériatrie de Montréal (CRIUGM), Montréal, QC, Canada
- Department of Stomatology, Faculty of Dentistry, Université de Montréal, Montreal, QC, Canada
| | - Julien Cohen-Adad
- NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique Montreal, Montreal, QC, Canada
- Functional Neuroimaging Unit, CRIUGM, Université de Montréal, Montreal, QC, Canada
| |
Collapse
|
40
|
PAM50: Unbiased multimodal template of the brainstem and spinal cord aligned with the ICBM152 space. Neuroimage 2018; 165:170-179. [DOI: 10.1016/j.neuroimage.2017.10.041] [Citation(s) in RCA: 96] [Impact Index Per Article: 16.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2017] [Revised: 09/16/2017] [Accepted: 10/20/2017] [Indexed: 11/17/2022] Open
|
41
|
Martin AR, Tadokoro N, Tetreault L, Arocho-Quinones EV, Budde MD, Kurpad SN, Fehlings MG. Imaging Evaluation of Degenerative Cervical Myelopathy. Neurosurg Clin N Am 2018; 29:33-45. [DOI: 10.1016/j.nec.2017.09.003] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
|
42
|
Saliani A, Perraud B, Duval T, Stikov N, Rossignol S, Cohen-Adad J. Axon and Myelin Morphology in Animal and Human Spinal Cord. Front Neuroanat 2017; 11:129. [PMID: 29311857 PMCID: PMC5743665 DOI: 10.3389/fnana.2017.00129] [Citation(s) in RCA: 51] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2017] [Accepted: 12/13/2017] [Indexed: 12/11/2022] Open
Abstract
Characterizing precisely the microstructure of axons, their density, size and myelination is of interest for the neuroscientific community, for example to help maximize the outcome of studies on white matter (WM) pathologies of the spinal cord (SC). The existence of a comprehensive and structured database of axonal measurements in healthy and disease models could help the validation of results obtained by different researchers. The purpose of this article is to provide such a database of healthy SC WM, to discuss the potential sources of variability and to suggest avenues for robust and accurate quantification of axon morphometry based on novel acquisition and processing techniques. The article is organized in three sections. The first section reviews morphometric results across species according to range of densities and counts of myelinated axons, axon diameter and myelin thickness, and characteristics of unmyelinated axons in different regions. The second section discusses the sources of variability across studies, such as age, sex, spinal pathways, spinal levels, statistical power and terminology in regard to tracts and protocols. The third section presents new techniques and perspectives that could benefit histology studies. For example, coherent anti-stokes Raman spectroscopy (CARS) imaging can provide sub-micrometric resolution without the need for fixation and staining, while slide scanners and stitching algorithms can provide full cross-sectional area of SC. In combination with these acquisition techniques, automatic segmentation algorithms for delineating axons and myelin sheath can help provide large-scale statistics on axon morphometry.
Collapse
Affiliation(s)
- Ariane Saliani
- NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique Montreal, Montreal, QC, Canada
| | - Blanche Perraud
- NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique Montreal, Montreal, QC, Canada
| | - Tanguy Duval
- NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique Montreal, Montreal, QC, Canada
| | - Nikola Stikov
- NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique Montreal, Montreal, QC, Canada
- Montreal Heart Institute, Montreal, QC, Canada
| | - Serge Rossignol
- Groupe de Recherche sur le Système Nerveux Central, Department of Neuroscience, Faculty of Medicine, Université de Montréal, Montreal, QC, Canada
| | - Julien Cohen-Adad
- NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique Montreal, Montreal, QC, Canada
- Functionnal Neuroimaging Unit, Centre de Recherche de l'Institut Universitaire de Gériatrie de Montréal, Université de Montréal, Montreal, QC, Canada
| |
Collapse
|
43
|
Yang G, Tian Q, Leuze C, Wintermark M, McNab JA. Double diffusion encoding MRI for the clinic. Magn Reson Med 2017; 80:507-520. [PMID: 29266375 DOI: 10.1002/mrm.27043] [Citation(s) in RCA: 54] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2017] [Revised: 11/17/2017] [Accepted: 11/18/2017] [Indexed: 12/13/2022]
Abstract
PURPOSE The purpose of this study is to develop double diffusion encoding (DDE) MRI methods for clinical use. Microscopic diffusion anisotropy measurements from DDE promise greater specificity to changes in tissue microstructure compared with conventional diffusion tensor imaging, but implementation of DDE sequences on whole-body MRI scanners is challenging because of the limited gradient strengths and lengthy acquisition times. METHODS A custom single-refocused DDE sequence was implemented on a 3T whole-body scanner. The DDE gradient orientation scheme and sequence parameters were optimized based on a Gaussian diffusion assumption. Using an optimized 5-min DDE acquisition, microscopic fractional anisotropy (μFA) maps were acquired for the first time in multiple sclerosis patients. RESULTS Based on simulations and in vivo human measurements, six parallel and six orthogonal diffusion gradient pairs were found to be the minimum number of diffusion gradient pairs necessary to produce a rotationally invariant measurement of μFA. Simulations showed that optimal precision and accuracy of μFA measurements were obtained using b-values between 1500 and 3000 s/mm2 . The μFA maps showed improved delineation of multiple sclerosis lesions compared with conventional fractional anisotropy and distinct contrast from T2 -weighted fluid attenuated inversion recovery and T1 -weighted imaging. CONCLUSION The μFA maps can be measured using DDE in a clinical setting and may provide new opportunities for characterizing multiple sclerosis lesions and other types of tissue degeneration. Magn Reson Med 80:507-520, 2018. © 2017 International Society for Magnetic Resonance in Medicine.
Collapse
Affiliation(s)
- Grant Yang
- Department of Electrical Engineering, Stanford University, Stanford, California, USA.,Department of Radiology, Stanford University, Stanford, California, USA
| | - Qiyuan Tian
- Department of Electrical Engineering, Stanford University, Stanford, California, USA.,Department of Radiology, Stanford University, Stanford, California, USA
| | - Christoph Leuze
- Department of Radiology, Stanford University, Stanford, California, USA
| | - Max Wintermark
- Department of Radiology, Stanford University, Stanford, California, USA
| | - Jennifer A McNab
- Department of Radiology, Stanford University, Stanford, California, USA
| |
Collapse
|
44
|
Berman S, West KL, Does MD, Yeatman JD, Mezer AA. Evaluating g-ratio weighted changes in the corpus callosum as a function of age and sex. Neuroimage 2017; 182:304-313. [PMID: 28673882 DOI: 10.1016/j.neuroimage.2017.06.076] [Citation(s) in RCA: 60] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2017] [Revised: 06/26/2017] [Accepted: 06/27/2017] [Indexed: 11/16/2022] Open
Abstract
Recent years have seen a growing interest in relating MRI measurements to the structural-biophysical properties of white matter fibers. The fiber g-ratio, defined as the ratio between the inner and outer radii of the axon myelin sheath, is an important structural property of white matter, affecting signal conduction. Recently proposed modeling methods that use a combination of quantitative-MRI signals, enable a measurement of the fiber g-ratio in vivo. Here we use an MRI-based g-ratio estimation to observe the variance of the g-ratio within the corpus callosum, and evaluate sex and age related differences. To estimate the g-ratio we used a model (Stikov et al., 2011; Duval et al., 2017) based on two different WM microstructure parameters: the relative amounts of myelin (myelin volume fraction, MVF) and fibers (fiber volume fraction, FVF) in a voxel. We derived the FVF from the fractional anisotropy (FA), and estimated the MVF by using the lipid and macromolecular tissue volume (MTV), calculated from the proton density (Mezer et al., 2013). In comparison to other methods of estimating the MVF, MTV represents a stable parameter with a straightforward route of acquisition. To establish our model, we first compared histological MVF measurements (West et al., 2016) with the MRI derived MTV. We then implemented our model on a large database of 92 subjects (44 males), aged 7 to 81, in order to evaluate age and sex related changes within the corpus callosum. Our results show that the MTV provides a good estimation of MVF for calculating g-ratio, and produced values from the corpus callosum that correspond to those found in animals ex vivo and are close to the theoretical optimum, as well as to published in vivo data. Our results demonstrate that the MTV derived g-ratio provides a simple and reliable in vivo g-ratio-weighted (GR*) measurement in humans. In agreement with theoretical predictions, and unlike other tissue parameters measured with MRI, the g-ratio estimations were found to be relatively stable with age, and we found no support for a significant sexual dimorphism with age.
Collapse
Affiliation(s)
- Shai Berman
- The Edmond and Lily Safra Center for Brain Sciences, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Kathryn L West
- Department of Biomedical Engineering, Vanderbilt University, Nashville, TN, USA; Vanderbilt University Institute of Imaging Science, Vanderbilt University, Nashville, TN, USA
| | - Mark D Does
- Department of Biomedical Engineering, Vanderbilt University, Nashville, TN, USA; Vanderbilt University Institute of Imaging Science, Vanderbilt University, Nashville, TN, USA
| | - Jason D Yeatman
- Institute for Learning & Brain Sciences and Department of Speech & Hearing Sciences, University of Washington, Seattle, WA, USA
| | - Aviv A Mezer
- The Edmond and Lily Safra Center for Brain Sciences, The Hebrew University of Jerusalem, Jerusalem, Israel.
| |
Collapse
|
45
|
By S, Xu J, Box BA, Bagnato FR, Smith SA. Application and evaluation of NODDI in the cervical spinal cord of multiple sclerosis patients. Neuroimage Clin 2017; 15:333-342. [PMID: 28560158 PMCID: PMC5443965 DOI: 10.1016/j.nicl.2017.05.010] [Citation(s) in RCA: 73] [Impact Index Per Article: 10.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2017] [Revised: 05/04/2017] [Accepted: 05/17/2017] [Indexed: 11/02/2022]
Abstract
INTRODUCTION There is a need to develop imaging methods sensitive to axonal injury in multiple sclerosis (MS), given the prominent impact of axonal pathology on disability and outcome. Advanced multi-compartmental diffusion models offer novel indices sensitive to white matter microstructure. One such model, neurite orientation dispersion and density imaging (NODDI), is sensitive to neurite morphology, providing indices of apparent volume fractions of axons (vin), isotropic water (viso) and the dispersion of fibers about a central axis (orientation dispersion index, ODI). NODDI has yet to be studied for its sensitivity to spinal cord pathology. Here, we investigate the feasibility and utility of NODDI in the cervical spinal cord of MS patients. METHODS NODDI was applied in the cervical spinal cord in a cohort of 8 controls and 6 MS patients. Statistical analyses were performed to test the sensitivity of NODDI-derived indices to pathology in MS (both lesion and normal appearing white matter NAWM). Diffusion kurtosis imaging (DKI) and diffusion tensor imaging (DTI) analysis were also performed to compare with NODDI. RESULTS A decrease in NODDI-derived vin was observed at the site of the lesion (p < 0.01), whereas a global increase in ODI was seen throughout white matter (p < 0.001). DKI-derived mean kurtosis (MK) and radial kurtosis (RK) and DTI-derived fractional anisotropy (FA) and radial diffusivity (RD) were all significantly different in MS patients (p < 0.02), however NODDI provided higher contrast between NAWM and lesion in all MS patients. CONCLUSION NODDI provides unique contrast that is not available with DKI or DTI, enabling improved characterization of the spinal cord in MS.
Collapse
Affiliation(s)
- Samantha By
- Department of Biomedical Engineering, Vanderbilt University, Nashville, TN, USA; Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, TN, USA
| | - Junzhong Xu
- Department of Biomedical Engineering, Vanderbilt University, Nashville, TN, USA; Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, TN, USA; Department of Radiology and Radiological Sciences, Vanderbilt University Medical Center, Nashville, TN, USA
| | - Bailey A Box
- Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, TN, USA
| | - Francesca R Bagnato
- Department of Neurology, Vanderbilt University Medical Center, Nashville, TN, USA
| | - Seth A Smith
- Department of Biomedical Engineering, Vanderbilt University, Nashville, TN, USA; Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, TN, USA; Department of Radiology and Radiological Sciences, Vanderbilt University Medical Center, Nashville, TN, USA.
| |
Collapse
|
46
|
Nunes D, Cruz TL, Jespersen SN, Shemesh N. Mapping axonal density and average diameter using non-monotonic time-dependent gradient-echo MRI. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2017; 277:117-130. [PMID: 28282586 DOI: 10.1016/j.jmr.2017.02.017] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/03/2017] [Revised: 02/23/2017] [Accepted: 02/25/2017] [Indexed: 06/06/2023]
Abstract
White Matter (WM) microstructures, such as axonal density and average diameter, are crucial to the normal function of the Central Nervous System (CNS) as they are closely related with axonal conduction velocities. Conversely, disruptions of these microstructural features may result in severe neurological deficits, suggesting that their noninvasive mapping could be an important step towards diagnosing and following pathophysiology. Whereas diffusion based MRI methods have been proposed to map these features, they typically entail the application of powerful gradients, which are rarely available in the clinic, or extremely long acquisition schemes to extract information from parameter-intensive models. In this study, we suggest that simple and time-efficient multi-gradient-echo (MGE) MRI can be used to extract the axon density from susceptibility-driven non-monotonic decay in the time-dependent signal. We show, both theoretically and with simulations, that a non-monotonic signal decay will occur for multi-compartmental microstructures - such as axons and extra-axonal spaces, which were here used as a simple model for the microstructure - and that, for axons parallel to the main magnetic field, the axonal density can be extracted. We then experimentally demonstrate in ex-vivo rat spinal cords that its different tracts - characterized by different microstructures - can be clearly contrasted using the MGE-derived maps. When the quantitative results are compared against ground-truth histology, they reflect the axonal fraction (though with a bias, as evident from Bland-Altman analysis). As well, the extra-axonal fraction can be estimated. The results suggest that our model is oversimplified, yet at the same time evidencing a potential and usefulness of the approach to map underlying microstructures using a simple and time-efficient MRI sequence. We further show that a simple general-linear-model can predict the average axonal diameters from the four model parameters, and map these average axonal diameters in the spinal cords. While clearly further modelling and theoretical developments are necessary, we conclude that salient WM microstructural features can be extracted from simple, SNR-efficient multi-gradient echo MRI, and that this paves the way towards easier estimation of WM microstructure in vivo.
Collapse
Affiliation(s)
- Daniel Nunes
- Champalimaud Neuroscience Programme, Champalimaud Centre for the Unknown, Av. Brasilia 1400-038, Lisbon, Portugal
| | - Tomás L Cruz
- Champalimaud Neuroscience Programme, Champalimaud Centre for the Unknown, Av. Brasilia 1400-038, Lisbon, Portugal
| | - 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 Neuroscience Programme, Champalimaud Centre for the Unknown, Av. Brasilia 1400-038, Lisbon, Portugal.
| |
Collapse
|
47
|
Abstract
Quantitative magnetic resonance imaging can be combined with advanced biophysical models to measure microstructural features of white matter. Non-invasive microstructural imaging has the potential to revolutionize neuroscience, and acquiring these measures in clinically feasible times would greatly improve patient monitoring and clinical studies of drug efficacy. However, a good understanding of microstructural imaging techniques is essential to set realistic expectations and to prevent over-interpretation of results. This review explains the methodology behind microstructural modeling and imaging, and gives an overview of the breakthroughs and challenges associated with it.
Collapse
|
48
|
Cohen Y, Anaby D, Morozov D. Diffusion MRI of the spinal cord: from structural studies to pathology. NMR IN BIOMEDICINE 2017; 30:e3592. [PMID: 27598689 DOI: 10.1002/nbm.3592] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/30/2015] [Revised: 06/01/2016] [Accepted: 07/05/2016] [Indexed: 05/27/2023]
Abstract
Diffusion MRI is extensively used to study brain microarchitecture and pathologies, and water diffusion appears highly anisotropic in the white matter (WM) of the spinal cord (SC). Despite these facts, the use of diffusion MRI to study the SC, which has increased in recent years, is much less common than that in the brain. In the present review, after a brief outline of early studies of diffusion MRI (DWI) and diffusion tensor MRI (DTI) of the SC, we provide a short survey on DTI and on diffusion MRI methods beyond the tensor that have been used to study SC microstructure and pathologies. After introducing the porous view of WM and describing the q-space approach and q-space diffusion MRI (QSI), we describe other methodologies that can be applied to study the SC. Selected applications of the use of DTI, QSI, and other more advanced diffusion MRI methods to study SC microstructure and pathologies are presented, with some emphasis on the use of less conventional diffusion methodologies. Because of length constraints, we concentrate on structural studies and on a few selected pathologies. Examples of the use of diffusion MRI to study dysmyelination, demyelination as in experimental autoimmune encephalomyelitis and multiple sclerosis, amyotrophic lateral sclerosis, and traumatic SC injury are presented. We conclude with a brief summary and a discussion of challenges and future directions for diffusion MRI of the SC. Copyright © 2016 John Wiley & Sons, Ltd.
Collapse
Affiliation(s)
- Yoram Cohen
- The Sackler School of Chemistry, The Raymond and Beverly Sackler Faculty of Exact Sciences, Tel Aviv University, Tel Aviv, Israel
- The Sagol School of Neuroscience, Tel Aviv University, Tel Aviv, Israel
| | - Debbie Anaby
- The Sackler School of Chemistry, The Raymond and Beverly Sackler Faculty of Exact Sciences, Tel Aviv University, Tel Aviv, Israel
| | - Darya Morozov
- The Sackler School of Chemistry, The Raymond and Beverly Sackler Faculty of Exact Sciences, Tel Aviv University, Tel Aviv, Israel
| |
Collapse
|
49
|
Mingasson T, Duval T, Stikov N, Cohen-Adad J. AxonPacking: An Open-Source Software to Simulate Arrangements of Axons in White Matter. Front Neuroinform 2017; 11:5. [PMID: 28197091 PMCID: PMC5281605 DOI: 10.3389/fninf.2017.00005] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2016] [Accepted: 01/13/2017] [Indexed: 11/24/2022] Open
Abstract
HIGHLIGHTSAxonPacking: Open-source software for simulating white matter microstructure. Validation on a theoretical disk packing problem. Reproducible and stable for various densities and diameter distributions. Can be used to study interplay between myelin/fiber density and restricted fraction.
Quantitative Magnetic Resonance Imaging (MRI) can provide parameters that describe white matter microstructure, such as the fiber volume fraction (FVF), the myelin volume fraction (MVF) or the axon volume fraction (AVF) via the fraction of restricted water (fr). While already being used for clinical application, the complex interplay between these parameters requires thorough validation via simulations. These simulations required a realistic, controlled and adaptable model of the white matter axons with the surrounding myelin sheath. While there already exist useful algorithms to perform this task, none of them combine optimisation of axon packing, presence of myelin sheath and availability as free and open source software. Here, we introduce a novel disk packing algorithm that addresses these issues. The performance of the algorithm is tested in term of reproducibility over 50 runs, resulting density, and stability over iterations. This tool was then used to derive multiple values of FVF and to study the impact of this parameter on fr and MVF in light of the known microstructure based on histology sample. The standard deviation of the axon density over runs was lower than 10−3 and the expected hexagonal packing for monodisperse disks was obtained with a density close to the optimal density (obtained: 0.892, theoretical: 0.907). Using an FVF ranging within [0.58, 0.82] and a mean inter-axon gap ranging within [0.1, 1.1] μm, MVF ranged within [0.32, 0.44] and fr ranged within [0.39, 0.71], which is consistent with the histology. The proposed algorithm is implemented in the open-source software AxonPacking (https://github.com/neuropoly/axonpacking) and can be useful for validating diffusion models as well as for enabling researchers to study the interplay between microstructure parameters when evaluating qMRI methods.
Collapse
Affiliation(s)
- Tom Mingasson
- NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique MontrealMontreal, QC, Canada; Signal Processing Department, École Centrale de NantesNantes, France
| | - Tanguy Duval
- NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique Montreal Montreal, QC, Canada
| | - Nikola Stikov
- NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique MontrealMontreal, QC, Canada; Department of Biomedical Engineering, Montreal Heart Institute, University of MontrealMontreal, QC, Canada
| | - Julien Cohen-Adad
- NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique MontrealMontreal, QC, Canada; Functional Neuroimaging Unit, Centre de Recherche de l'Institut Universitaire de Gériatrie de Montréal (CRIUGM), Université de MontréalMontreal, QC, Canada
| |
Collapse
|
50
|
Duval T, Le Vy S, Stikov N, Campbell J, Mezer A, Witzel T, Keil B, Smith V, Wald LL, Klawiter E, Cohen-Adad J. g-Ratio weighted imaging of the human spinal cord in vivo. Neuroimage 2017; 145:11-23. [PMID: 27664830 PMCID: PMC5179300 DOI: 10.1016/j.neuroimage.2016.09.018] [Citation(s) in RCA: 63] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2016] [Revised: 08/22/2016] [Accepted: 09/08/2016] [Indexed: 12/13/2022] Open
Abstract
The fiber g-ratio is defined as the ratio of the inner to the outer diameter of the myelin sheath. This ratio provides a measure of the myelin thickness that complements axon morphology (diameter and density) for assessment of demyelination in diseases such as multiple sclerosis. Previous work has shown that an aggregate g-ratio map can be computed using a formula that combines axon and myelin density measured with quantitative MRI. In this work, we computed g-ratio weighted maps in the cervical spinal cord of nine healthy subjects. We utilized the 300mT/m gradients from the CONNECTOM scanner to estimate the fraction of restricted water (fr) with high accuracy, using the CHARMED model. Myelin density was estimated using the lipid and macromolecular tissue volume (MTV) method, derived from normalized proton density (PD) mapping. The variability across spinal level, laterality and subject were assessed using a three-way ANOVA. The average g-ratio value obtained in the white matter was 0.76+/-0.03, consistent with previous histology work. Coefficients of variation of fr and MTV were respectively 4.3% and 13.7%. fr and myelin density were significantly different across spinal tracts (p=3×10-7 and 0.004 respectively) and were positively correlated in the white matter (r=0.42), suggesting shared microstructural information. The aggregate g-ratio did not show significant differences across tracts (p=0.6). This study suggests that fr and myelin density can be measured in vivo with high precision and that they can be combined to produce a g-ratio-weighted map robust to free water pool contamination from cerebrospinal fluid or veins. Potential applications include the study of early demyelination in multiple sclerosis, and the quantitative assessment of remyelination drugs.
Collapse
Affiliation(s)
- T Duval
- NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique Montreal, Montreal, QC, Canada
| | - S Le Vy
- NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique Montreal, Montreal, QC, Canada; Functional Neuroimaging Unit, CRIUGM, Université de Montréal, Montréal, QC, Canada
| | - N Stikov
- NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique Montreal, Montreal, QC, Canada; Montreal Heart Institute, Montreal, QC, Canada
| | - J Campbell
- Montreal Neurological Institute, McGill University, Montreal, QC, Canada
| | - A Mezer
- Edmond and Lily Safra Center for Brain Sciences (ELSC), The Hebrew University of Jerusalem, Jerusalem, Israel
| | - T Witzel
- A.A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Harvard Medical School, Charlestown, MA, United States
| | - B Keil
- A.A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Harvard Medical School, Charlestown, MA, United States
| | - V Smith
- A.A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Harvard Medical School, Charlestown, MA, United States
| | - L L Wald
- A.A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Harvard Medical School, Charlestown, MA, United States
| | - E Klawiter
- A.A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Harvard Medical School, Charlestown, MA, United States
| | - J Cohen-Adad
- NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique Montreal, Montreal, QC, Canada; Functional Neuroimaging Unit, CRIUGM, Université de Montréal, Montréal, QC, Canada.
| |
Collapse
|