1
|
Winther S, Peulicke O, Andersson M, Kjer HM, Bærentzen JA, Dyrby TB. Exploring white matter dynamics and morphology through interactive numerical phantoms: the White Matter Generator. Front Neuroinform 2024; 18:1354708. [PMID: 39144684 PMCID: PMC11322502 DOI: 10.3389/fninf.2024.1354708] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2023] [Accepted: 06/25/2024] [Indexed: 08/16/2024] Open
Abstract
Brain white matter is a dynamic environment that continuously adapts and reorganizes in response to stimuli and pathological changes. Glial cells, especially, play a key role in tissue repair, inflammation modulation, and neural recovery. The movements of glial cells and changes in their concentrations can influence the surrounding axon morphology. We introduce the White Matter Generator (WMG) tool to enable the study of how axon morphology is influenced through such dynamical processes, and how this, in turn, influences the diffusion-weighted MRI signal. This is made possible by allowing interactive changes to the configuration of the phantom generation throughout the optimization process. The phantoms can consist of myelinated axons, unmyelinated axons, and cell clusters, separated by extra-cellular space. Due to morphological flexibility and computational advantages during the optimization, the tool uses ellipsoids as building blocks for all structures; chains of ellipsoids for axons, and individual ellipsoids for cell clusters. After optimization, the ellipsoid representation can be converted to a mesh representation which can be employed in Monte-Carlo diffusion simulations. This offers an effective method for evaluating tissue microstructure models for diffusion-weighted MRI in controlled bio-mimicking white matter environments. Hence, the WMG offers valuable insights into white matter's adaptive nature and implications for diffusion-weighted MRI microstructure models, and thereby holds the potential to advance clinical diagnosis, treatment, and rehabilitation strategies for various neurological disorders and injuries.
Collapse
Affiliation(s)
- Sidsel Winther
- Department of Applied Mathematics and Computer Science, Technical University of Denmark, Kongens Lyngby, Denmark
- Danish Research Centre for Magnetic Resonance, Centre for Functional and Diagnostic Imaging and Research, Copenhagen University Hospital—Amager and Hvidovre, Hvidovre, Denmark
| | - Oscar Peulicke
- Department of Applied Mathematics and Computer Science, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Mariam Andersson
- Danish Research Centre for Magnetic Resonance, Centre for Functional and Diagnostic Imaging and Research, Copenhagen University Hospital—Amager and Hvidovre, Hvidovre, Denmark
| | - Hans M. Kjer
- Department of Applied Mathematics and Computer Science, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Jakob A. Bærentzen
- Department of Applied Mathematics and Computer Science, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Tim B. Dyrby
- Department of Applied Mathematics and Computer Science, Technical University of Denmark, Kongens Lyngby, Denmark
- Danish Research Centre for Magnetic Resonance, Centre for Functional and Diagnostic Imaging and Research, Copenhagen University Hospital—Amager and Hvidovre, Hvidovre, Denmark
| |
Collapse
|
2
|
Zhou M, Stobbe R, Szczepankiewicz F, Budde M, Buck B, Kate M, Lloret M, Fairall P, Butcher K, Shuaib A, Emery D, Nilsson M, Westin CF, Beaulieu C. Tensor-valued diffusion MRI of human acute stroke. Magn Reson Med 2024; 91:2126-2141. [PMID: 38156813 DOI: 10.1002/mrm.29975] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2023] [Revised: 11/18/2023] [Accepted: 11/27/2023] [Indexed: 01/03/2024]
Abstract
PURPOSE Tensor-valued diffusion encoding can disentangle orientation dispersion and subvoxel anisotropy, potentially offering insight into microstructural changes after cerebral ischemia. The purpose was to evaluate tensor-valued diffusion MRI in human acute ischemic stroke, assess potential confounders from diffusion time dependencies, and compare to Monte Carlo diffusion simulations of axon beading. METHODS Linear (LTE) and spherical (STE) b-tensor encoding with inherently different effective diffusion times were acquired in 21 acute ischemic stroke patients between 3 and 57 h post-onset at 3 T in 2.5 min. In an additional 10 patients, STE with 2 LTE yielding different effective diffusion times were acquired for comparison. Diffusional variance decomposition (DIVIDE) was used to estimate microscopic anisotropy (μFA), as well as anisotropic, isotropic, and total diffusional variance (MKA , MKI , MKT ). DIVIDE parameters, and diffusion tensor imaging (DTI)-derived mean diffusivity and fractional anisotropy (FA) were compared in lesion versus contralateral white matter. Monte Carlo diffusion simulations of various cylindrical geometries for all b-tensor protocols were used to interpret parameter measurements. RESULTS MD was ˜40% lower in lesions for all LTE/STE protocols. The DIVIDE parameters varied with effective diffusion time: higher μFA and MKA in lesion versus contralateral white matter for STE with longer effective diffusion time LTE, whereas the shorter effective diffusion time LTE protocol yielded lower μFA and MKA in lesions. Both protocols, regardless of diffusion time, were consistent with simulations of greater beading amplitude and intracellular volume fraction. CONCLUSION DIVIDE parameters depend on diffusion time in acute stroke but consistently indicate neurite beading and larger intracellular volume fraction.
Collapse
Affiliation(s)
- Mi Zhou
- Biomedical Engineering, University of Alberta, Edmonton, Alberta, Canada
| | - Robert Stobbe
- Biomedical Engineering, University of Alberta, Edmonton, Alberta, Canada
- Radiology and Diagnostic Imaging, University of Alberta, Edmonton, Alberta, Canada
| | | | - Matthew Budde
- Neurosurgery, Medical College of Wisconsin, Milwaukee, Wisconsin, USA
| | - Brian Buck
- Neurology, University of Alberta, Edmonton, Alberta, Canada
| | - Mahesh Kate
- Neurology, University of Alberta, Edmonton, Alberta, Canada
| | - Mar Lloret
- Neurology, University of Alberta, Edmonton, Alberta, Canada
| | - Paige Fairall
- Neurology, University of Alberta, Edmonton, Alberta, Canada
| | - Ken Butcher
- School of Clinical Medicine, University of New South Wales, Sydney, New South Wales, Australia
| | - Ashfaq Shuaib
- Neurology, University of Alberta, Edmonton, Alberta, Canada
| | - Derek Emery
- Radiology and Diagnostic Imaging, University of Alberta, Edmonton, Alberta, Canada
| | - Markus Nilsson
- Clinical Sciences Lund, Lund University, Lund, Scania, Sweden
| | - Carl-Fredrik Westin
- Radiology, Brigham and Women's Hospital, Harvard Medical School, Boston, MA, USA
| | - Christian Beaulieu
- Biomedical Engineering, University of Alberta, Edmonton, Alberta, Canada
- Radiology and Diagnostic Imaging, University of Alberta, Edmonton, Alberta, Canada
| |
Collapse
|
3
|
Fang C, Yang Z, Wassermann D, Li JR. A simulation-driven supervised learning framework to estimate brain microstructure using diffusion MRI. Med Image Anal 2023; 90:102979. [PMID: 37827109 DOI: 10.1016/j.media.2023.102979] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2023] [Revised: 09/13/2023] [Accepted: 09/22/2023] [Indexed: 10/14/2023]
Abstract
We propose a framework to train supervised learning models on synthetic data to estimate brain microstructure parameters using diffusion magnetic resonance imaging (dMRI). Although further validation is necessary, the proposed framework aims to seamlessly incorporate realistic simulations into dMRI microstructure estimation. Synthetic data were generated from over 1,000 neuron meshes converted from digital neuronal reconstructions and linked to their neuroanatomical parameters (such as soma volume and neurite length) using an optimized diffusion MRI simulator that produces intracellular dMRI signals from the solution of the Bloch-Torrey partial differential equation. By combining random subsets of simulated neuron signals with a free diffusion compartment signal, we constructed a synthetic dataset containing dMRI signals and 40 tissue microstructure parameters of 1.45 million artificial brain voxels. To implement supervised learning models we chose multilayer perceptrons (MLPs) and trained them on a subset of the synthetic dataset to estimate some microstructure parameters, namely, the volume fractions of soma, neurites, and the free diffusion compartment, as well as the area fractions of soma and neurites. The trained MLPs perform satisfactorily on the synthetic test sets and give promising in-vivo parameter maps on the MGH Connectome Diffusion Microstructure Dataset (CDMD). Most importantly, the estimated volume fractions showed low dependence on the diffusion time, the diffusion time independence of the estimated parameters being a desired property of quantitative microstructure imaging. The synthetic dataset we generated will be valuable for the validation of models that map between the dMRI signals and microstructure parameters. The surface meshes and microstructures parameters of the aforementioned neurons have been made publicly available.
Collapse
Affiliation(s)
- Chengran Fang
- INRIA Saclay, Equipe IDEFIX, UMA, ENSTA Paris, 828, Boulevard des Maréchaux, 91762 Palaiseau, France; INRIA Saclay, Equipe MIND, 1 Rue Honoré d'Estienne d'Orves, 91120 Palaiseau, France
| | - Zheyi Yang
- INRIA Saclay, Equipe IDEFIX, UMA, ENSTA Paris, 828, Boulevard des Maréchaux, 91762 Palaiseau, France
| | - Demian Wassermann
- INRIA Saclay, Equipe MIND, 1 Rue Honoré d'Estienne d'Orves, 91120 Palaiseau, France
| | - Jing-Rebecca Li
- INRIA Saclay, Equipe IDEFIX, UMA, ENSTA Paris, 828, Boulevard des Maréchaux, 91762 Palaiseau, France.
| |
Collapse
|
4
|
Villarreal-Haro JL, Gardier R, Canales-Rodríguez EJ, Fischi-Gomez E, Girard G, Thiran JP, Rafael-Patiño J. CACTUS: a computational framework for generating realistic white matter microstructure substrates. Front Neuroinform 2023; 17:1208073. [PMID: 37603781 PMCID: PMC10434236 DOI: 10.3389/fninf.2023.1208073] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2023] [Accepted: 07/13/2023] [Indexed: 08/23/2023] Open
Abstract
Monte-Carlo diffusion simulations are a powerful tool for validating tissue microstructure models by generating synthetic diffusion-weighted magnetic resonance images (DW-MRI) in controlled environments. This is fundamental for understanding the link between micrometre-scale tissue properties and DW-MRI signals measured at the millimetre-scale, optimizing acquisition protocols to target microstructure properties of interest, and exploring the robustness and accuracy of estimation methods. However, accurate simulations require substrates that reflect the main microstructural features of the studied tissue. To address this challenge, we introduce a novel computational workflow, CACTUS (Computational Axonal Configurator for Tailored and Ultradense Substrates), for generating synthetic white matter substrates. Our approach allows constructing substrates with higher packing density than existing methods, up to 95% intra-axonal volume fraction, and larger voxel sizes of up to 500μm3 with rich fibre complexity. CACTUS generates bundles with angular dispersion, bundle crossings, and variations along the fibres of their inner and outer radii and g-ratio. We achieve this by introducing a novel global cost function and a fibre radial growth approach that allows substrates to match predefined targeted characteristics and mirror those reported in histological studies. CACTUS improves the development of complex synthetic substrates, paving the way for future applications in microstructure imaging.
Collapse
Affiliation(s)
- Juan Luis Villarreal-Haro
- Signal Processing Laboratory (LTS5), École Polytechnique Frale de Lausanne (EPFL), Lausanne, Switzerland
| | - Remy Gardier
- Signal Processing Laboratory (LTS5), École Polytechnique Frale de Lausanne (EPFL), Lausanne, Switzerland
| | - Erick J. Canales-Rodríguez
- Signal Processing Laboratory (LTS5), École Polytechnique Frale de Lausanne (EPFL), Lausanne, Switzerland
| | - Elda Fischi-Gomez
- Signal Processing Laboratory (LTS5), École Polytechnique Frale de Lausanne (EPFL), Lausanne, Switzerland
- CIBM Center for Biomedical Imaging, Lausanne, Switzerland
- Radiology Department, Centre Hospitalier Universitaire Vaudois, University of Lausanne, Lausanne, Switzerland
| | - Gabriel Girard
- Signal Processing Laboratory (LTS5), École Polytechnique Frale de Lausanne (EPFL), Lausanne, Switzerland
- CIBM Center for Biomedical Imaging, Lausanne, Switzerland
- Radiology Department, Centre Hospitalier Universitaire Vaudois, University of Lausanne, Lausanne, Switzerland
- Department of Computer Science, University of Sherbrooke, Sherbrooke, QC, Canada
| | - Jean-Philippe Thiran
- Signal Processing Laboratory (LTS5), École Polytechnique Frale de Lausanne (EPFL), Lausanne, Switzerland
- CIBM Center for Biomedical Imaging, Lausanne, Switzerland
- Radiology Department, Centre Hospitalier Universitaire Vaudois, University of Lausanne, Lausanne, Switzerland
| | - Jonathan Rafael-Patiño
- Signal Processing Laboratory (LTS5), École Polytechnique Frale de Lausanne (EPFL), Lausanne, Switzerland
- Radiology Department, Centre Hospitalier Universitaire Vaudois, University of Lausanne, Lausanne, Switzerland
| |
Collapse
|
5
|
Using ex vivo arterial injection and dissection to assess white matter vascularization. Sci Rep 2023; 13:809. [PMID: 36646713 PMCID: PMC9842749 DOI: 10.1038/s41598-022-26227-6] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2022] [Accepted: 12/12/2022] [Indexed: 01/18/2023] Open
Abstract
Advances in the techniques for assessing human cerebral white matter have recently contributed to greater attention to structural connectivity. Yet, little is known about the vascularization of most white matter fasciculi and the fascicular composition of the vascular territories. This paper presents an original method to label the arterial supply of macroscopic white matter fasciculi based on a standardized protocol for post-mortem injection of colored material into main cerebral arteries combined with a novel fiber dissection technique. Twelve whole human cerebral hemispheres obtained post-mortem were included. A detailed description of every step, from obtaining the specimen to image acquisition of its dissection, is provided. Injection and dissection were reproducible and manageable without any sophisticated equipment. They successfully showed the arterial supply of the dissected fasciculi. In addition, we discuss the challenges we faced and overcame during the development of the presented method, highlight its originality. Henceforth, this innovative method serves as a tool to provide a precise anatomical description of the vascularization of the main white matter tracts.
Collapse
|
6
|
Jelescu IO, de Skowronski A, Geffroy F, Palombo M, Novikov DS. Neurite Exchange Imaging (NEXI): A minimal model of diffusion in gray matter with inter-compartment water exchange. Neuroimage 2022; 256:119277. [PMID: 35523369 PMCID: PMC10363376 DOI: 10.1016/j.neuroimage.2022.119277] [Citation(s) in RCA: 46] [Impact Index Per Article: 23.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2022] [Revised: 04/26/2022] [Accepted: 05/01/2022] [Indexed: 01/18/2023] Open
Abstract
Biophysical models of diffusion in white matter have been center-stage over the past two decades and are essentially based on what is now commonly referred to as the "Standard Model" (SM) of non-exchanging anisotropic compartments with Gaussian diffusion. In this work, we focus on diffusion MRI in gray matter, which requires rethinking basic microstructure modeling blocks. In particular, at least three contributions beyond the SM need to be considered for gray matter: water exchange across the cell membrane - between neurites and the extracellular space; non-Gaussian diffusion along neuronal and glial processes - resulting from structural disorder; and signal contribution from soma. For the first contribution, we propose Neurite Exchange Imaging (NEXI) as an extension of the SM of diffusion, which builds on the anisotropic Kärger model of two exchanging compartments. Using datasets acquired at multiple diffusion weightings (b) and diffusion times (t) in the rat brain in vivo, we investigate the suitability of NEXI to describe the diffusion signal in the gray matter, compared to the other two possible contributions. Our results for the diffusion time window 20-45 ms show minimal diffusivity time-dependence and more pronounced kurtosis decay with time, which is well fit by the exchange model. Moreover, we observe lower signal for longer diffusion times at high b. In light of these observations, we identify exchange as the mechanism that best explains these signal signatures in both low-b and high-b regime, and thereby propose NEXI as the minimal model for gray matter microstructure mapping. We finally highlight multi-b multi-t acquisition protocols as being best suited to estimate NEXI model parameters reliably. Using this approach, we estimate the inter-compartment water exchange time to be 15 - 60 ms in the rat cortex and hippocampus in vivo, which is of the same order or shorter than the diffusion time in typical diffusion MRI acquisitions. This suggests water exchange as an essential component for interpreting diffusion MRI measurements in gray matter.
Collapse
Affiliation(s)
- Ileana O Jelescu
- CIBM Center for Biomedical Imaging, Animal Imaging and Technology, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland; Department of Radiology, Lausanne University Hospital (CHUV), Lausanne, Switzerland; School of Biology and Medicine, University of Lausanne (UNIL), Lausanne, Switzerland.
| | - Alexandre de Skowronski
- CIBM Center for Biomedical Imaging, Animal Imaging and Technology, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
| | | | - Marco Palombo
- School of Psychology, Cardiff University Brain Research Imaging Centre (CUBRIC), Cardiff University, Cardiff, UK; School of Computer Science and Informatics, Cardiff University, Cardiff, UK; Department of Computer Science, Centre for Medical Image Computing, University College London, London, UK
| | - Dmitry S Novikov
- Department of Radiology, Center for Biomedical Imaging, New York University School of Medicine, New York, NY, USA
| |
Collapse
|
7
|
Lu Q, Liu W, Zhuo Z, Li Y, Duan Y, Yu P, Qu L, Ye C, Liu Y. A Transfer Learning Approach to Few-shot Segmentation of Novel White Matter Tracts. Med Image Anal 2022; 79:102454. [DOI: 10.1016/j.media.2022.102454] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2021] [Revised: 03/19/2022] [Accepted: 04/08/2022] [Indexed: 12/20/2022]
|
8
|
Li Z, Pang Z, Cheng J, Hsu YC, Sun Y, Özarslan E, Bai R. The direction-dependence of apparent water exchange rate in human white matter. Neuroimage 2021; 247:118831. [PMID: 34923129 DOI: 10.1016/j.neuroimage.2021.118831] [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: 08/22/2021] [Revised: 12/10/2021] [Accepted: 12/15/2021] [Indexed: 11/29/2022] Open
Abstract
Transmembrane water exchange is a potential biomarker in the diagnosis and understanding of cancers, brain disorders, and other diseases. Filter-exchange imaging (FEXI), a special case of diffusion exchange spectroscopy adapted for clinical applications, has the potential to reveal different physiological water exchange processes. However, it is still controversial whether modulating the diffusion encoding gradient direction can affect the apparent exchange rate (AXR) measurements of FEXI in white matter (WM) where water diffusion shows strong anisotropy. In this study, we explored the diffusion-encoding direction dependence of FEXI in human brain white matter by performing FEXI with 20 diffusion-encoding directions on a clinical 3T scanner in-vivo. The results show that the AXR values measured when the gradients are perpendicular to the fiber orientation (0.77 ± 0.13 s - 1, mean ± standard deviation of all the subjects) are significantly larger than the AXR estimates when the gradients are parallel to the fiber orientation (0.33 ± 0.14 s - 1, p < 0.001) in WM voxels with coherently-orientated fibers. In addition, no significant correlation is found between AXRs measured along these two directions, indicating that they are measuring different water exchange processes. What's more, only the perpendicular AXR rather than the parallel AXR shows dependence on axonal diameter, indicating that the perpendicular AXR might reflect transmembrane water exchange between intra-axonal and extra-cellular spaces. Further finite difference (FD) simulations having three water compartments (intra-axonal, intra-glial, and extra-cellular spaces) to mimic WM micro-environments also suggest that the perpendicular AXR is more sensitive to the axonal water transmembrane exchange than parallel AXR. Taken together, our results show that AXR measured along different directions could be utilized to probe different water exchange processes in WM.
Collapse
Affiliation(s)
- Zhaoqing Li
- Department of Physical Medicine and Rehabilitation of the Affiliated Sir Run Shaw Hospital AND Interdisciplinary Institute of Neuroscience and Technology, School of Medicine, Zhejiang University, Hangzhou, China; Key Laboratory of Biomedical Engineering of Education Ministry, College of Biomedical Engineering and Instrument Science, Zhejiang University, Hangzhou, China
| | - Zhenfeng Pang
- Department of Chemistry, Zhejiang University, Hangzhou, China
| | - Juange Cheng
- Department of Physical Medicine and Rehabilitation of the Affiliated Sir Run Shaw Hospital AND Interdisciplinary Institute of Neuroscience and Technology, School of Medicine, Zhejiang University, Hangzhou, China
| | - Yi-Cheng Hsu
- MR Collaboration, Siemens Healthcare, Shanghai, China
| | - Yi Sun
- MR Collaboration, Siemens Healthcare, Shanghai, China
| | - Evren Özarslan
- Department of Biomedical Engineering, Linköping University, Linköping, Sweden; Center for Medical Image Science and Visualization, Linköping University, Linköping, Sweden
| | - Ruiliang Bai
- Department of Physical Medicine and Rehabilitation of the Affiliated Sir Run Shaw Hospital AND Interdisciplinary Institute of Neuroscience and Technology, School of Medicine, Zhejiang University, Hangzhou, China; Key Laboratory of Biomedical Engineering of Education Ministry, College of Biomedical Engineering and Instrument Science, Zhejiang University, Hangzhou, China.
| |
Collapse
|
9
|
Drobnjak I, Neher P, Poupon C, Sarwar T. Physical and digital phantoms for validating tractography and assessing artifacts. Neuroimage 2021; 245:118704. [PMID: 34748954 DOI: 10.1016/j.neuroimage.2021.118704] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2021] [Revised: 10/01/2021] [Accepted: 11/01/2021] [Indexed: 11/17/2022] Open
Abstract
Fiber tractography is widely used to non-invasively map white-matter bundles in vivo using diffusion-weighted magnetic resonance imaging (dMRI). As it is the case for all scientific methods, proper validation is a key prerequisite for the successful application of fiber tractography, be it in the area of basic neuroscience or in a clinical setting. It is well-known that the indirect estimation of the fiber tracts from the local diffusion signal is highly ambiguous and extremely challenging. Furthermore, the validation of fiber tractography methods is hampered by the lack of a real ground truth, which is caused by the extremely complex brain microstructure that is not directly observable non-invasively and that is the basis of the huge network of long-range fiber connections in the brain that are the actual target of fiber tractography methods. As a substitute for in vivo data with a real ground truth that could be used for validation, a widely and successfully employed approach is the use of synthetic phantoms. In this work, we are providing an overview of the state-of-the-art in the area of physical and digital phantoms, answering the following guiding questions: "What are dMRI phantoms and what are they good for?", "What would the ideal phantom for validation fiber tractography look like?" and "What phantoms, phantom datasets and tools used for their creation are available to the research community?". We will further discuss the limitations and opportunities that come with the use of dMRI phantoms, and what future direction this field of research might take.
Collapse
Affiliation(s)
- Ivana Drobnjak
- Center for Medical Image Computing, Department of Computer Science, University College London, UK.
| | - Peter Neher
- Division of Medical Image Computing, German Cancer Research Center (DKFZ), Heidelberg, Germany
| | - Cyril Poupon
- BAOBAB, NeuroSpin, Commissariat à l'Energie Atomique, Institut des Sciences du Vivant Frédéric Joliot, Gif-sur-Yvette, France
| | - Tabinda Sarwar
- School of Computing Technologies, RMIT University, Australia
| |
Collapse
|
10
|
Smirnov M, Destrieux C, Maldonado IL. Cerebral white matter vasculature: still uncharted? Brain 2021; 144:3561-3575. [PMID: 34718425 DOI: 10.1093/brain/awab273] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2021] [Revised: 06/17/2021] [Accepted: 07/11/2021] [Indexed: 11/14/2022] Open
Abstract
White matter vasculature plays a major role in the pathophysiology of permanent neurological deficits following a stroke or progressive cognitive alteration related to small vessel disease. Thus, knowledge of the complex vascularization and functional aspects of the deep white matter territories is paramount to comprehend clinical manifestations of brain ischemia. This review provides a structured presentation of the existing knowledge of the vascularization of the human cerebral white matter from seminal historical studies to the current literature. First, we revisit the highlights of prenatal development of the endoparenchymal telencephalic vascular system that are crucial for the understanding of vessel organization in the adult. Second, we reveal the tangled history of debates on the existence, clinical significance, and physiological role of leptomeningeal anastomoses. Then, we present how conceptions on white matter vascularization transitioned from the mixed ventriculopetal/ventriculofugal theory, in which a low-flow area was interposed in between concurrent arterial flows, to the purely ventriculopetal theory. The latter model explains variable white matter sensitivity to ischemia by various organizations of ventriculopetal vessel terminals having different origin/length properties and interconnection patterns. Next, arteries supplying primarily the white matter are described according to their length and overall structure. Furthermore, the known distribution territories, to date, are studied in relation to primary anatomical structures of the human cerebral white matter, emphasizing the sparsity of the "ground-truth" data available in the literature. Finally, the implications for both large vessel occlusion and chronic small vessel disease are discussed, as well as the insights from neuroimaging. All things considered, we identify the need for further research on deep white matter vascularization, especially regarding the arterial supply of white matter fiber tracts.
Collapse
Affiliation(s)
- Mykyta Smirnov
- UMR 1253, iBrain, Université de Tours, Inserm, Tours, France
| | - Christophe Destrieux
- UMR 1253, iBrain, Université de Tours, Inserm, Tours, France.,CHRU de Tours, Tours, France
| | | |
Collapse
|
11
|
Ianus A, Alexander DC, Zhang H, Palombo M. Mapping complex cell morphology in the grey matter with double diffusion encoding MR: A simulation study. Neuroimage 2021; 241:118424. [PMID: 34311067 PMCID: PMC8961003 DOI: 10.1016/j.neuroimage.2021.118424] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2021] [Revised: 07/13/2021] [Accepted: 07/21/2021] [Indexed: 01/18/2023] Open
Abstract
This paper investigates the impact of cell body (namely soma) size and branching of cellular projections on diffusion MR imaging (dMRI) and spectroscopy (dMRS) signals for both standard single diffusion encoding (SDE) and more advanced double diffusion encoding (DDE) measurements using numerical simulations. The aim is to investigate the ability of dMRI/dMRS to characterize the complex morphology of brain cells focusing on these two distinctive features of brain grey matter. To this end, we employ a recently developed computational framework to create three dimensional meshes of neuron-like structures for Monte Carlo simulations, using diffusion coefficients typical of water and brain metabolites. Modelling the cellular structure as realistically connected spherical soma and cylindrical cellular projections, we cover a wide range of combinations of sphere radii and branching order of cellular projections, characteristic of various grey matter cells. We assess the impact of spherical soma size and branching order on the b-value dependence of the SDE signal as well as the time dependence of the mean diffusivity (MD) and mean kurtosis (MK). Moreover, we also assess the impact of spherical soma size and branching order on the angular modulation of DDE signal at different mixing times, together with the mixing time dependence of the apparent microscopic anisotropy (μA), a promising contrast derived from DDE measurements. The SDE results show that spherical soma size has a measurable impact on both the b-value dependence of the SDE signal and the MD and MK diffusion time dependence for both water and metabolites. On the other hand, we show that branching order has little impact on either, especially for water. In contrast, the DDE results show that spherical soma size has a measurable impact on the DDE signal's angular modulation at short mixing times and the branching order of cellular projections significantly impacts the mixing time dependence of the DDE signal's angular modulation as well as of the derived μA, for both water and metabolites. Our results confirm that SDE based techniques may be sensitive to spherical soma size, and most importantly, show for the first time that DDE measurements may be more sensitive to the dendritic tree complexity (as parametrized by the branching order of cellular projections), paving the way for new ways of characterizing grey matter morphology, non-invasively using dMRS and potentially dMRI.
Collapse
Affiliation(s)
- A Ianus
- Centre for Medical Image Computing and Department of Computer Science, University College London, London, United Kingdom; Champalimaud Research, Champalimaud Centre for the Unknown, Lisbon, Portugal
| | - D C Alexander
- Centre for Medical Image Computing and Department of Computer Science, University College London, London, United Kingdom
| | - H Zhang
- Centre for Medical Image Computing and Department of Computer Science, University College London, London, United Kingdom
| | - M Palombo
- Centre for Medical Image Computing and Department of Computer Science, University College London, London, United Kingdom.
| |
Collapse
|
12
|
Yang JYM, Yeh CH, Poupon C, Calamante F. Diffusion MRI tractography for neurosurgery: the basics, current state, technical reliability and challenges. Phys Med Biol 2021; 66. [PMID: 34157706 DOI: 10.1088/1361-6560/ac0d90] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2021] [Accepted: 06/22/2021] [Indexed: 01/20/2023]
Abstract
Diffusion magnetic resonance imaging (dMRI) tractography is currently the only imaging technique that allows for non-invasive delineation and visualisation of white matter (WM) tractsin vivo,prompting rapid advances in related fields of brain MRI research in recent years. One of its major clinical applications is for pre-surgical planning and intraoperative image guidance in neurosurgery, where knowledge about the location of WM tracts nearby the surgical target can be helpful to guide surgical resection and optimise post-surgical outcomes. Surgical injuries to these WM tracts can lead to permanent neurological and functional deficits, making the accuracy of tractography reconstructions paramount. The quality of dMRI tractography is influenced by many modifiable factors, ranging from MRI data acquisition through to the post-processing of tractography output, with the potential of error propagation based on decisions made at each and subsequent processing steps. Research over the last 25 years has significantly improved the anatomical accuracy of tractography. An updated review about tractography methodology in the context of neurosurgery is now timely given the thriving research activities in dMRI, to ensure more appropriate applications in the clinical neurosurgical realm. This article aims to review the dMRI physics, and tractography methodologies, highlighting recent advances to provide the key concepts of tractography-informed neurosurgery, with a focus on the general considerations, the current state of practice, technical challenges, potential advances, and future demands to this field.
Collapse
Affiliation(s)
- Joseph Yuan-Mou Yang
- Department of Neurosurgery, The Royal Children's Hospital, Melbourne, Australia.,Neuroscience Research, Murdoch Children's Research Institute, Melbourne, Australia.,Developmental Imaging, Murdoch Children's Research Institute, Melbourne, Australia.,Department of Paediatrics, The University of Melbourne, Melbourne, Australia
| | - Chun-Hung Yeh
- Institute for Radiological Research, Chang Gung University and Chang Gung Memorial Hospital, Taoyuan, Taiwan.,Department of Child and Adolescent Psychiatry, Chang Gung Memorial Hospital, Linkou Medical Center, Taoyuan, Taiwan
| | - Cyril Poupon
- NeuroSpin, Frédéric Joliot Life Sciences Institute, CEA, CNRS, Paris-Saclay University, Gif-sur-Yvette, France
| | - Fernando Calamante
- The University of Sydney, Sydney Imaging, Sydney, Australia.,The University of Sydney, School of Biomedical Engineering, Sydney, Australia
| |
Collapse
|
13
|
A simulation study investigating potential diffusion-based MRI signatures of microstrokes. Sci Rep 2021; 11:14229. [PMID: 34244549 PMCID: PMC8271016 DOI: 10.1038/s41598-021-93503-2] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2020] [Accepted: 06/22/2021] [Indexed: 02/06/2023] Open
Abstract
Recent studies suggested that cerebrovascular micro-occlusions, i.e. microstokes, could lead to ischemic tissue infarctions and cognitive deficits. Due to their small size, identifying measurable biomarkers of these microvascular lesions remains a major challenge. This work aims to simulate potential MRI signatures combining arterial spin labeling (ASL) and multi-directional diffusion-weighted imaging (DWI). Driving our hypothesis are recent observations demonstrating a radial reorientation of microvasculature around the micro-infarction locus during recovery in mice. Synthetic capillary beds, randomly- and radially-oriented, and optical coherence tomography (OCT) angiograms, acquired in the barrel cortex of mice (n = 5) before and after inducing targeted photothrombosis, were analyzed. Computational vascular graphs combined with a 3D Monte-Carlo simulator were used to characterize the magnetic resonance (MR) response, encompassing the effects of magnetic field perturbations caused by deoxyhemoglobin, and the advection and diffusion of the nuclear spins. We quantified the minimal intravoxel signal loss ratio when applying multiple gradient directions, at varying sequence parameters with and without ASL. With ASL, our results demonstrate a significant difference (p < 0.05) between the signal-ratios computed at baseline and 3 weeks after photothrombosis. The statistical power further increased (p < 0.005) using angiograms measured at week 4. Without ASL, no reliable signal change was found. We found that higher ratios, and accordingly improved significance, were achieved at lower magnetic field strengths (e.g., B0 = 3T) and shorter echo time TE (< 16 ms). Our simulations suggest that microstrokes might be characterized through ASL-DWI sequence, providing necessary insights for posterior experimental validations, and ultimately, future translational trials.
Collapse
|
14
|
Abdollahzadeh A, Belevich I, Jokitalo E, Sierra A, Tohka J. DeepACSON automated segmentation of white matter in 3D electron microscopy. Commun Biol 2021; 4:179. [PMID: 33568775 PMCID: PMC7876004 DOI: 10.1038/s42003-021-01699-w] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/16/2020] [Accepted: 01/12/2021] [Indexed: 01/30/2023] Open
Abstract
Tracing the entirety of ultrastructures in large three-dimensional electron microscopy (3D-EM) images of the brain tissue requires automated segmentation techniques. Current segmentation techniques use deep convolutional neural networks (DCNNs) and rely on high-contrast cellular membranes and high-resolution EM volumes. On the other hand, segmenting low-resolution, large EM volumes requires methods to account for severe membrane discontinuities inescapable. Therefore, we developed DeepACSON, which performs DCNN-based semantic segmentation and shape-decomposition-based instance segmentation. DeepACSON instance segmentation uses the tubularity of myelinated axons and decomposes under-segmented myelinated axons into their constituent axons. We applied DeepACSON to ten EM volumes of rats after sham-operation or traumatic brain injury, segmenting hundreds of thousands of long-span myelinated axons, thousands of cell nuclei, and millions of mitochondria with excellent evaluation scores. DeepACSON quantified the morphology and spatial aspects of white matter ultrastructures, capturing nanoscopic morphological alterations five months after the injury.
Collapse
Affiliation(s)
- Ali Abdollahzadeh
- A. I. Virtanen Institute for Molecular Sciences, University of Eastern Finland, Kuopio, Finland
| | - Ilya Belevich
- Electron Microscopy Unit, Institute of Biotechnology, University of Helsinki, Helsinki, Finland
| | - Eija Jokitalo
- Electron Microscopy Unit, Institute of Biotechnology, University of Helsinki, Helsinki, Finland
| | - Alejandra Sierra
- A. I. Virtanen Institute for Molecular Sciences, University of Eastern Finland, Kuopio, Finland.
| | - Jussi Tohka
- A. I. Virtanen Institute for Molecular Sciences, University of Eastern Finland, Kuopio, Finland
| |
Collapse
|
15
|
Andersson M, Kjer HM, Rafael-Patino J, Pacureanu A, Pakkenberg B, Thiran JP, Ptito M, Bech M, Bjorholm Dahl A, Andersen Dahl V, Dyrby TB. Axon morphology is modulated by the local environment and impacts the noninvasive investigation of its structure-function relationship. Proc Natl Acad Sci U S A 2020; 117:33649-33659. [PMID: 33376224 PMCID: PMC7777205 DOI: 10.1073/pnas.2012533117] [Citation(s) in RCA: 37] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Abstract
Axonal conduction velocity, which ensures efficient function of the brain network, is related to axon diameter. Noninvasive, in vivo axon diameter estimates can be made with diffusion magnetic resonance imaging, but the technique requires three-dimensional (3D) validation. Here, high-resolution, 3D synchrotron X-ray nano-holotomography images of white matter samples from the corpus callosum of a monkey brain reveal that blood vessels, cells, and vacuoles affect axonal diameter and trajectory. Within single axons, we find that the variation in diameter and conduction velocity correlates with the mean diameter, contesting the value of precise diameter determination in larger axons. These complex 3D axon morphologies drive previously reported 2D trends in axon diameter and g-ratio. Furthermore, we find that these morphologies bias the estimates of axon diameter with diffusion magnetic resonance imaging and, ultimately, impact the investigation and formulation of the axon structure-function relationship.
Collapse
Affiliation(s)
- Mariam Andersson
- Danish Research Centre for Magnetic Resonance, Centre for Functional and Diagnostic Imaging and Research, Copenhagen University Hospital Hvidovre, 2650 Hvidovre, Denmark;
- Department of Applied Mathematics and Computer Science, Technical University of Denmark, 2800 Kongens Lyngby, Denmark
| | - Hans Martin Kjer
- Danish Research Centre for Magnetic Resonance, Centre for Functional and Diagnostic Imaging and Research, Copenhagen University Hospital Hvidovre, 2650 Hvidovre, Denmark
- Department of Applied Mathematics and Computer Science, Technical University of Denmark, 2800 Kongens Lyngby, Denmark
| | - Jonathan Rafael-Patino
- Signal Processing Laboratory (LTS5), École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | | | - Bente Pakkenberg
- Research Laboratory for Stereology and Neuroscience, Copenhagen University Hospital, Bispebjerg, 2400 Copenhagen, Denmark
| | - Jean-Philippe Thiran
- Signal Processing Laboratory (LTS5), École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
- Radiology Department, Centre Hospitalier Universitaire Vaudois and University of Lausanne, 1011 Lausanne, Switzerland
- Center for Biomedical Imaging, 1015 Lausanne, Switzerland
| | - Maurice Ptito
- School of Optometry, University of Montreal, Montreal, QC H3T 1P1, Canada
- Department of Neuroscience, Faculty of Health Science, University of Copenhagen, 2200 Copenhagen, Denmark
| | - Martin Bech
- Division of Medical Radiation Physics, Department of Clinical Sciences, Lund University, 221 85 Lund, Sweden
| | - Anders Bjorholm Dahl
- Department of Applied Mathematics and Computer Science, Technical University of Denmark, 2800 Kongens Lyngby, Denmark
| | - Vedrana Andersen Dahl
- Department of Applied Mathematics and Computer Science, Technical University of Denmark, 2800 Kongens Lyngby, Denmark
| | - Tim B Dyrby
- Danish Research Centre for Magnetic Resonance, Centre for Functional and Diagnostic Imaging and Research, Copenhagen University Hospital Hvidovre, 2650 Hvidovre, Denmark;
- Department of Applied Mathematics and Computer Science, Technical University of Denmark, 2800 Kongens Lyngby, Denmark
| |
Collapse
|
16
|
Harkins KD, Beaulieu C, Xu J, Gore JC, Does MD. A simple estimate of axon size with diffusion MRI. Neuroimage 2020; 227:117619. [PMID: 33301942 PMCID: PMC7949481 DOI: 10.1016/j.neuroimage.2020.117619] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2020] [Revised: 11/06/2020] [Accepted: 11/29/2020] [Indexed: 12/18/2022] Open
Abstract
Noninvasive estimation of mean axon diameter presents a new opportunity to explore white matter plasticity, development, and pathology. Several diffusion-weighted MRI (DW-MRI) methods have been proposed to measure the average axon diameter in white matter, but they typically require many diffusion encoding measurements and complicated mathematical models to fit the signal to multiple tissue compartments, including intra- and extra-axonal spaces. Here, Monte Carlo simulations uncovered a straightforward DW-MRI metric of axon diameter: the change in radial apparent diffusion coefficient estimated at different effective diffusion times, ΔD⊥. Simulations indicated that this metric increases monotonically within a relevant range of effective mean axon diameter while being insensitive to changes in extra-axonal volume fraction, axon diameter distribution, g-ratio, and influence of myelin water. Also, a monotonic relationship was found to exist for signals coming from both intra- and extra-axonal compartments. The slope in ΔD⊥ with effective axon diameter increased with the difference in diffusion time of both oscillating and pulsed gradient diffusion sequences.
Collapse
Affiliation(s)
- Kevin D Harkins
- Biomedical Engineering, Vanderbilt University, United States; Institute of Imaging Science, Vanderbilt University, United States.
| | | | - Junzhong Xu
- Institute of Imaging Science, Vanderbilt University, United States; Radiology and Radiological Sciences, Vanderbilt University Medical Center, United States
| | - John C Gore
- Biomedical Engineering, Vanderbilt University, United States; Institute of Imaging Science, Vanderbilt University, United States; Radiology and Radiological Sciences, Vanderbilt University Medical Center, United States
| | - Mark D Does
- Biomedical Engineering, Vanderbilt University, United States; Institute of Imaging Science, Vanderbilt University, United States
| |
Collapse
|
17
|
Lee HH, Fieremans E, Novikov DS. Realistic Microstructure Simulator (RMS): Monte Carlo simulations of diffusion in three-dimensional cell segmentations of microscopy images. J Neurosci Methods 2020; 350:109018. [PMID: 33279478 DOI: 10.1016/j.jneumeth.2020.109018] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2020] [Revised: 11/16/2020] [Accepted: 11/29/2020] [Indexed: 12/12/2022]
Abstract
BACKGROUND Monte Carlo simulations of diffusion are commonly used as a model validation tool as they are especially suitable for generating the diffusion MRI signal in complicated tissue microgeometries. NEW METHOD Here we describe the details of implementing Monte Carlo simulations in three-dimensional (3d) voxelized segmentations of cells in microscopy images. Using the concept of the corner reflector, we largely reduce the computational load of simulating diffusion within and exchange between multiple cells. Precision is further achieved by GPU-based parallel computations. RESULTS Our simulation of diffusion in white matter axons segmented from a mouse brain demonstrates its value in validating biophysical models. Furthermore, we provide the theoretical background for implementing a discretized diffusion process, and consider the finite-step effects of the particle-membrane reflection and permeation events, needed for efficient simulation of interactions with irregular boundaries, spatially variable diffusion coefficient, and exchange. COMPARISON WITH EXISTING METHODS To our knowledge, this is the first Monte Carlo pipeline for MR signal simulations in a substrate composed of numerous realistic cells, accounting for their permeable and irregularly-shaped membranes. CONCLUSIONS The proposed RMS pipeline makes it possible to achieve fast and accurate simulations of diffusion in realistic tissue microgeometry, as well as the interplay with other MR contrasts. Presently, RMS focuses on simulations of diffusion, exchange, and T1 and T2 NMR relaxation in static tissues, with a possibility to straightforwardly account for susceptibility-induced T2* effects and flow.
Collapse
Affiliation(s)
- Hong-Hsi Lee
- Center for Biomedical Imaging, Department of Radiology, New York University School of Medicine, New York, NY, USA; Center for Advanced Imaging Innovation and Research (CAI2R), New York University School of Medicine, New York, NY, USA.
| | - Els Fieremans
- Center for Biomedical Imaging, Department of Radiology, New York University School of Medicine, New York, NY, USA; Center for Advanced Imaging Innovation and Research (CAI2R), New York University School of Medicine, New York, NY, USA
| | - Dmitry S Novikov
- Center for Biomedical Imaging, Department of Radiology, New York University School of Medicine, New York, NY, USA; Center for Advanced Imaging Innovation and Research (CAI2R), New York University School of Medicine, New York, NY, USA
| |
Collapse
|
18
|
Fang C, Nguyen VD, Wassermann D, Li JR. Diffusion MRI simulation of realistic neurons with SpinDoctor and the Neuron Module. Neuroimage 2020; 222:117198. [PMID: 32730957 DOI: 10.1016/j.neuroimage.2020.117198] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2019] [Revised: 06/30/2020] [Accepted: 07/22/2020] [Indexed: 02/08/2023] Open
Abstract
The diffusion MRI signal arising from neurons can be numerically simulated by solving the Bloch-Torrey partial differential equation. In this paper we present the Neuron Module that we implemented within the Matlab-based diffusion MRI simulation toolbox SpinDoctor. SpinDoctor uses finite element discretization and adaptive time integration to solve the Bloch-Torrey partial differential equation for general diffusion-encoding sequences, at multiple b-values and in multiple diffusion directions. In order to facilitate the diffusion MRI simulation of realistic neurons by the research community, we constructed finite element meshes for a group of 36 pyramidal neurons and a group of 29 spindle neurons whose morphological descriptions were found in the publicly available neuron repository NeuroMorpho.Org. These finite elements meshes range from having 15,163 nodes to 622,553 nodes. We also broke the neurons into the soma and dendrite branches and created finite elements meshes for these cell components. Through the Neuron Module, these neuron and cell components finite element meshes can be seamlessly coupled with the functionalities of SpinDoctor to provide the diffusion MRI signal attributable to spins inside neurons. We make these meshes and the source code of the Neuron Module available to the public as an open-source package. To illustrate some potential uses of the Neuron Module, we show numerical examples of the simulated diffusion MRI signals in multiple diffusion directions from whole neurons as well as from the soma and dendrite branches, and include a comparison of the high b-value behavior between dendrite branches and whole neurons. In addition, we demonstrate that the neuron meshes can be used to perform Monte-Carlo diffusion MRI simulations as well. We show that at equivalent accuracy, if only one gradient direction needs to be simulated, SpinDoctor is faster than a GPU implementation of Monte-Carlo, but if many gradient directions need to be simulated, there is a break-even point when the GPU implementation of Monte-Carlo becomes faster than SpinDoctor. Furthermore, we numerically compute the eigenfunctions and the eigenvalues of the Bloch-Torrey and the Laplace operators on the neuron geometries using a finite elements discretization, in order to give guidance in the choice of the space and time discretization parameters for both finite elements and Monte-Carlo approaches. Finally, we perform a statistical study on the set of 65 neurons to test some candidate biomakers that can potentially indicate the soma size. This preliminary study exemplifies the possible research that can be conducted using the Neuron Module.
Collapse
Affiliation(s)
- Chengran Fang
- INRIA Saclay, Equipe DEFI, CMAP, Ecole Polytechnique, 91128 Palaiseau Cedex, France; INRIA Saclay, Equipe Parietal, 1 Rue Honoré d'Estienne d'Orves, 91120 Palaiseau, France
| | - Van-Dang Nguyen
- Department of Computational Science and Technology, KTH Royal Institute of Technology, Sweden
| | - Demian Wassermann
- INRIA Saclay, Equipe Parietal, 1 Rue Honoré d'Estienne d'Orves, 91120 Palaiseau, France
| | - Jing-Rebecca Li
- INRIA Saclay, Equipe DEFI, CMAP, Ecole Polytechnique, 91128 Palaiseau Cedex, France.
| |
Collapse
|
19
|
Henriques RN, Palombo M, Jespersen SN, Shemesh N, Lundell H, Ianuş A. Double diffusion encoding and applications for biomedical imaging. J Neurosci Methods 2020; 348:108989. [PMID: 33144100 DOI: 10.1016/j.jneumeth.2020.108989] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/09/2020] [Revised: 09/25/2020] [Accepted: 10/20/2020] [Indexed: 12/11/2022]
Abstract
Diffusion Magnetic Resonance Imaging (dMRI) is one of the most important contemporary non-invasive modalities for probing tissue structure at the microscopic scale. The majority of dMRI techniques employ standard single diffusion encoding (SDE) measurements, covering different sequence parameter ranges depending on the complexity of the method. Although many signal representations and biophysical models have been proposed for SDE data, they are intrinsically limited by a lack of specificity. Advanced dMRI methods have been proposed to provide additional microstructural information beyond what can be inferred from SDE. These enhanced contrasts can play important roles in characterizing biological tissues, for instance upon diseases (e.g. neurodegenerative, cancer, stroke), aging, learning, and development. In this review we focus on double diffusion encoding (DDE), which stands out among other advanced acquisitions for its versatility, ability to probe more specific diffusion correlations, and feasibility for preclinical and clinical applications. Various DDE methodologies have been employed to probe compartment sizes (Section 3), decouple the effects of microscopic diffusion anisotropy from orientation dispersion (Section 4), probe displacement correlations, study exchange, or suppress fast diffusing compartments (Section 6). DDE measurements can also be used to improve the robustness of biophysical models (Section 5) and study intra-cellular diffusion via magnetic resonance spectroscopy of metabolites (Section 7). This review discusses all these topics as well as important practical aspects related to the implementation and contrast in preclinical and clinical settings (Section 9) and aims to provide the readers a guide for deciding on the right DDE acquisition for their specific application.
Collapse
Affiliation(s)
- Rafael N Henriques
- Champalimaud Research, Champalimaud Centre for the Unknown, Lisbon, Portugal
| | - Marco Palombo
- Centre for Medical Image Computing and Dept. of Computer Science, University College London, London, UK
| | - Sune N Jespersen
- Center of Functionally Integrative Neuroscience (CFIN) and MINDLab, Department of Clinical Medicine, Aarhus University, Aarhus, Denmark; Department of Physics and Astronomy, Aarhus University, Aarhus, Denmark
| | - Noam Shemesh
- Champalimaud Research, Champalimaud Centre for the Unknown, Lisbon, Portugal
| | - Henrik Lundell
- Danish Research Centre for Magnetic Resonance, Centre for Functional and Diagnostic Imaging and Research, Copenhagen University Hospital Hvidovre, Denmark
| | - Andrada Ianuş
- Champalimaud Research, Champalimaud Centre for the Unknown, Lisbon, Portugal.
| |
Collapse
|
20
|
Hill I, Palombo M, Santin M, Branzoli F, Philippe AC, Wassermann D, Aigrot MS, Stankoff B, Baron-Van Evercooren A, Felfli M, Langui D, Zhang H, Lehericy S, Petiet A, Alexander DC, Ciccarelli O, Drobnjak I. Machine learning based white matter models with permeability: An experimental study in cuprizone treated in-vivo mouse model of axonal demyelination. Neuroimage 2020; 224:117425. [PMID: 33035669 DOI: 10.1016/j.neuroimage.2020.117425] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2019] [Revised: 09/29/2020] [Accepted: 09/30/2020] [Indexed: 01/14/2023] Open
Abstract
The intra-axonal water exchange time (τi), a parameter associated with axonal permeability, could be an important biomarker for understanding and treating demyelinating pathologies such as Multiple Sclerosis. Diffusion-Weighted MRI (DW-MRI) is sensitive to changes in permeability; however, the parameter has so far remained elusive due to the lack of general biophysical models that incorporate it. Machine learning based computational models can potentially be used to estimate such parameters. Recently, for the first time, a theoretical framework using a random forest (RF) regressor suggests that this is a promising new approach for permeability estimation. In this study, we adopt such an approach and for the first time experimentally investigate it for demyelinating pathologies through direct comparison with histology. We construct a computational model using Monte Carlo simulations and an RF regressor in order to learn a mapping between features derived from DW-MRI signals and ground truth microstructure parameters. We test our model in simulations, and find strong correlations between the predicted and ground truth parameters (intra-axonal volume fraction f: R2 =0.99, τi: R2 =0.84, intrinsic diffusivity d: R2 =0.99). We then apply the model in-vivo, on a controlled cuprizone (CPZ) mouse model of demyelination, comparing the results from two cohorts of mice, CPZ (N=8) and healthy age-matched wild-type (WT, N=8). We find that the RF model estimates sensible microstructure parameters for both groups, matching values found in literature. Furthermore, we perform histology for both groups using electron microscopy (EM), measuring the thickness of the myelin sheath as a surrogate for exchange time. Histology results show that our RF model estimates are very strongly correlated with the EM measurements (ρ = 0.98 for f, ρ = 0.82 for τi). Finally, we find a statistically significant decrease in τi in all three regions of the corpus callosum (splenium/genu/body) of the CPZ cohort (<τi>=310ms/330ms/350ms) compared to the WT group (<τi>=370ms/370ms/380ms). This is in line with our expectations that τi is lower in regions where the myelin sheath is damaged, as axonal membranes become more permeable. Overall, these results demonstrate, for the first time experimentally and in vivo, that a computational model learned from simulations can reliably estimate microstructure parameters, including the axonal permeability .
Collapse
Affiliation(s)
- Ioana Hill
- Centre for Medical Image Computing and Dept of Computer Science, University College London, London, UK
| | - Marco Palombo
- Centre for Medical Image Computing and Dept of Computer Science, University College London, London, UK.
| | - Mathieu Santin
- Institut du Cerveau et de la Moelle épinière, ICM, Sorbonne Université, Inserm 1127, CNRS UMR 7225, F-75013, Paris, France; Institut du Cerveau et de la Moelle épinière, ICM, Centre de NeuroImagerie de Recherche, CENIR, Paris, France
| | - Francesca Branzoli
- Institut du Cerveau et de la Moelle épinière, ICM, Sorbonne Université, Inserm 1127, CNRS UMR 7225, F-75013, Paris, France; Institut du Cerveau et de la Moelle épinière, ICM, Centre de NeuroImagerie de Recherche, CENIR, Paris, France
| | - Anne-Charlotte Philippe
- Institut du Cerveau et de la Moelle épinière, ICM, Sorbonne Université, Inserm 1127, CNRS UMR 7225, F-75013, Paris, France
| | - Demian Wassermann
- Université Côte d'Azur, Inria, Sophia-Antipolis, France; Parietal, CEA, Inria, Saclay, Île-de-France
| | - Marie-Stephane Aigrot
- Institut du Cerveau et de la Moelle épinière, ICM, Sorbonne Université, Inserm 1127, CNRS UMR 7225, F-75013, Paris, France
| | - Bruno Stankoff
- Institut du Cerveau et de la Moelle épinière, ICM, Sorbonne Université, Inserm 1127, CNRS UMR 7225, F-75013, Paris, France; AP-HP, Hôpital Saint-Antoine, Paris, France
| | - Anne Baron-Van Evercooren
- Institut du Cerveau et de la Moelle épinière, ICM, Sorbonne Université, Inserm 1127, CNRS UMR 7225, F-75013, Paris, France
| | - Mehdi Felfli
- Institut du Cerveau et de la Moelle épinière, ICM, Sorbonne Université, Inserm 1127, CNRS UMR 7225, F-75013, Paris, France
| | - Dominique Langui
- Institut du Cerveau et de la Moelle épinière, ICM, Sorbonne Université, Inserm 1127, CNRS UMR 7225, F-75013, Paris, France
| | - Hui Zhang
- Centre for Medical Image Computing and Dept of Computer Science, University College London, London, UK
| | - Stephane Lehericy
- Institut du Cerveau et de la Moelle épinière, ICM, Sorbonne Université, Inserm 1127, CNRS UMR 7225, F-75013, Paris, France; Institut du Cerveau et de la Moelle épinière, ICM, Centre de NeuroImagerie de Recherche, CENIR, Paris, France
| | - Alexandra Petiet
- Institut du Cerveau et de la Moelle épinière, ICM, Sorbonne Université, Inserm 1127, CNRS UMR 7225, F-75013, Paris, France; Institut du Cerveau et de la Moelle épinière, ICM, Centre de NeuroImagerie de Recherche, CENIR, Paris, France
| | - Daniel C Alexander
- Centre for Medical Image Computing and Dept of Computer Science, University College London, London, UK
| | - Olga Ciccarelli
- Dept. of Neuroinflammation, University College London, Queen Square Institute of Neurology, University College London, London, UK
| | - Ivana Drobnjak
- Centre for Medical Image Computing and Dept of Computer Science, University College London, London, UK
| |
Collapse
|
21
|
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: 67] [Impact Index Per Article: 16.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
|
22
|
ConFiG: Contextual Fibre Growth to generate realistic axonal packing for diffusion MRI simulation. Neuroimage 2020; 220:117107. [PMID: 32622984 PMCID: PMC7903162 DOI: 10.1016/j.neuroimage.2020.117107] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2020] [Revised: 06/17/2020] [Accepted: 06/25/2020] [Indexed: 11/27/2022] Open
Abstract
This paper presents Contextual Fibre Growth (ConFiG), an approach to generate white matter numerical phantoms by mimicking natural fibre genesis. ConFiG grows fibres one-by-one, following simple rules motivated by real axonal guidance mechanisms. These simple rules enable ConFiG to generate phantoms with tuneable microstructural features by growing fibres while attempting to meet morphological targets such as user-specified density and orientation distribution. We compare ConFiG to the state-of-the-art approach based on packing fibres together by generating phantoms in a range of fibre configurations including crossing fibre bundles and orientation dispersion. Results demonstrate that ConFiG produces phantoms with up to 20% higher densities than the state-of-the-art, particularly in complex configurations with crossing fibres. We additionally show that the microstructural morphology of ConFiG phantoms is comparable to real tissue, producing diameter and orientation distributions close to electron microscopy estimates from real tissue as well as capturing complex fibre cross sections. Signals simulated from ConFiG phantoms match real diffusion MRI data well, showing that ConFiG phantoms can be used to generate realistic diffusion MRI data. This demonstrates the feasibility of ConFiG to generate realistic synthetic diffusion MRI data for developing and validating microstructure modelling approaches. We present ConFiG, a biologically motivated numerical phantom generator for white matter. ConFiG produces phantoms with state-of-the-art density and realistic microstructure. Diffusion MRI simulations in ConFiG phantoms are comparable to real dMRI signals.
Collapse
|
23
|
Li JR, Nguyen VD, Tran TN, Valdman J, Trang CB, Nguyen KV, Vu DTS, Tran HA, Tran HTA, Nguyen TMP. SpinDoctor: A MATLAB toolbox for diffusion MRI simulation. Neuroimage 2019; 202:116120. [DOI: 10.1016/j.neuroimage.2019.116120] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2019] [Revised: 08/21/2019] [Accepted: 08/22/2019] [Indexed: 12/15/2022] Open
|