1
|
Lee HH, Tian Q, Sheft M, Coronado-Leija R, Ramos-Llorden G, Abdollahzadeh A, Fieremans E, Novikov DS, Huang SY. The effects of axonal beading and undulation on axonal diameter estimation from diffusion MRI: Insights from simulations in human axons segmented from three-dimensional electron microscopy. NMR IN BIOMEDICINE 2024; 37:e5087. [PMID: 38168082 PMCID: PMC10942763 DOI: 10.1002/nbm.5087] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/16/2023] [Revised: 11/19/2023] [Accepted: 11/21/2023] [Indexed: 01/05/2024]
Abstract
The increasing availability of high-performance gradient systems in human MRI scanners has generated great interest in diffusion microstructural imaging applications such as axonal diameter mapping. Practically, sensitivity to axon diameter in diffusion MRI is attained at strong diffusion weightings b , where the deviation from the expected 1 / b scaling in white matter yields a finite transverse diffusivity, which is then translated into an axon diameter estimate. While axons are usually modeled as perfectly straight, impermeable cylinders, local variations in diameter (caliber variation or beading) and direction (undulation) are known to influence axonal diameter estimates and have been observed in microscopy data of human axons. In this study, we performed Monte Carlo simulations of diffusion in axons reconstructed from three-dimensional electron microscopy of a human temporal lobe specimen using simulated sequence parameters matched to the maximal gradient strength of the next-generation Connectome 2.0 human MRI scanner ( ≲ 500 mT/m). We show that axon diameter estimation is accurate for nonbeaded, nonundulating fibers; however, in fibers with caliber variations and undulations, the axon diameter is heavily underestimated due to caliber variations, and this effect overshadows the known overestimation of the axon diameter due to undulations. This unexpected underestimation may originate from variations in the coarse-grained axial diffusivity due to caliber variations. Given that increased axonal beading and undulations have been observed in pathological tissues, such as traumatic brain injury and ischemia, the interpretation of axon diameter alterations in pathology may be significantly confounded.
Collapse
Affiliation(s)
- Hong-Hsi Lee
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, Massachusetts, USA
- Harvard Medical School, Boston, Massachusetts, USA
| | - Qiyuan Tian
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, Massachusetts, USA
- Harvard Medical School, Boston, Massachusetts, USA
| | - Maxina Sheft
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, Massachusetts, USA
- Harvard–MIT Health Sciences and Technology, Cambridge, Massachusetts, USA
| | - Ricardo Coronado-Leija
- Center for Biomedical Imaging, Department of Radiology, New York University School of Medicine, New York, New York, USA
- Center for Advanced Imaging Innovation and Research (CAI2R), New York University School of Medicine, New York, New York, USA
| | - Gabriel Ramos-Llorden
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, Massachusetts, USA
- Harvard Medical School, Boston, Massachusetts, USA
| | - Ali Abdollahzadeh
- Center for Biomedical Imaging, Department of Radiology, New York University School of Medicine, New York, New York, USA
- Center for Advanced Imaging Innovation and Research (CAI2R), New York University School of Medicine, New York, New York, USA
| | - Els Fieremans
- Center for Biomedical Imaging, Department of Radiology, New York University School of Medicine, New York, New York, USA
- Center for Advanced Imaging Innovation and Research (CAI2R), New York University School of Medicine, New York, New York, USA
| | - Dmitry S. Novikov
- Center for Biomedical Imaging, Department of Radiology, New York University School of Medicine, New York, New York, USA
- Center for Advanced Imaging Innovation and Research (CAI2R), New York University School of Medicine, New York, New York, USA
| | - Susie Y. Huang
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, Massachusetts, USA
- Harvard Medical School, Boston, Massachusetts, USA
| |
Collapse
|
2
|
Chowdhury R, Wan J, Gardier R, Rafael-Patino J, Thiran JP, Gibou F, Mukherjee A. Molecular Imaging with Aquaporin-Based Reporter Genes: Quantitative Considerations from Monte Carlo Diffusion Simulations. ACS Synth Biol 2023; 12:3041-3049. [PMID: 37793076 DOI: 10.1021/acssynbio.3c00372] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/06/2023]
Abstract
Aquaporins provide a unique approach for imaging genetic activity in deep tissues by increasing the rate of cellular water diffusion, which generates a magnetic resonance contrast. However, distinguishing aquaporin signals from the tissue background is challenging because water diffusion is influenced by structural factors, such as cell size and packing density. Here, we developed a Monte Carlo model to analyze how cell radius and intracellular volume fraction quantitatively affect aquaporin signals. We demonstrated that a differential imaging approach based on subtracting signals at two diffusion times can improve specificity by unambiguously isolating aquaporin signals from the tissue background. We further used Monte Carlo simulations to analyze the connection between diffusivity and the percentage of cells engineered to express aquaporin and established a mapping that accurately determined the volume fraction of aquaporin-expressing cells in mixed populations. The quantitative framework developed in this study will enable a broad range of applications in biomedical synthetic biology, requiring the use of aquaporins to noninvasively monitor the location and function of genetically engineered devices in live animals.
Collapse
Affiliation(s)
| | | | - Remy Gardier
- Signal Processing Laboratory (LTS5), Ecole Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland
| | - Jonathan Rafael-Patino
- Signal Processing Laboratory (LTS5), Ecole Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland
- Radiology Department, Centre Hospitalier Universitaire Vaudois (CHUV) and University of Lausanne (UNIL), 1005 Lausanne, Switzerland
| | - Jean-Philippe Thiran
- Signal Processing Laboratory (LTS5), Ecole Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland
- Radiology Department, Centre Hospitalier Universitaire Vaudois (CHUV) and University of Lausanne (UNIL), 1005 Lausanne, Switzerland
| | | | | |
Collapse
|
3
|
Yang Z, Fang C, Li JR. Incorporating interface permeability into the diffusion MRI signal representation while using impermeable Laplace eigenfunctions. Phys Med Biol 2023; 68:175036. [PMID: 37579758 DOI: 10.1088/1361-6560/acf022] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2023] [Accepted: 08/14/2023] [Indexed: 08/16/2023]
Abstract
Objective. The complex-valued transverse magnetization due to diffusion-encoding magnetic field gradients acting on a permeable medium can be modeled by the Bloch-Torrey partial differential equation. The diffusion magnetic resonance imaging (MRI) signal has a representation in the basis of the Laplace eigenfunctions of the medium. However, in order to estimate the permeability coefficient from diffusion MRI data, it is desirable that the forward solution can be calculated efficiently for many values of permeability.Approach. In this paper we propose a new formulation of the permeable diffusion MRI signal representation in the basis of the Laplace eigenfunctions of the same medium where the interfaces are made impermeable.Main results.We proved the theoretical equivalence between our new formulation and the original formulation in the case that the full eigendecomposition is used. We validated our method numerically and showed promising numerical results when a partial eigendecomposition is used. Two diffusion MRI sequences were used to illustrate the numerical validity of our new method.Significance.Our approach means that the same basis (the impermeable set) can be used for all permeability values, which reduces the computational time significantly, enabling the study of the effects of the permeability coefficient on the diffusion MRI signal in the future.
Collapse
Affiliation(s)
- Zheyi Yang
- Equipe IDEFIX, INRIA Saclay, UMA, ENSTA PARIS, Palaiseau, France
| | - Chengran Fang
- Equipe IDEFIX, INRIA Saclay, UMA, ENSTA PARIS, Palaiseau, France
| | - Jing-Rebecca Li
- Equipe IDEFIX, INRIA Saclay, UMA, ENSTA PARIS, Palaiseau, France
| |
Collapse
|
4
|
Girard G, Rafael-Patiño J, Truffet R, Aydogan DB, Adluru N, Nair VA, Prabhakaran V, Bendlin BB, Alexander AL, Bosticardo S, Gabusi I, Ocampo-Pineda M, Battocchio M, Piskorova Z, Bontempi P, Schiavi S, Daducci A, Stafiej A, Ciupek D, Bogusz F, Pieciak T, Frigo M, Sedlar S, Deslauriers-Gauthier S, Kojčić I, Zucchelli M, Laghrissi H, Ji Y, Deriche R, Schilling KG, Landman BA, Cacciola A, Basile GA, Bertino S, Newlin N, Kanakaraj P, Rheault F, Filipiak P, Shepherd TM, Lin YC, Placantonakis DG, Boada FE, Baete SH, Hernández-Gutiérrez E, Ramírez-Manzanares A, Coronado-Leija R, Stack-Sánchez P, Concha L, Descoteaux M, Mansour L S, Seguin C, Zalesky A, Marshall K, Canales-Rodríguez EJ, Wu Y, Ahmad S, Yap PT, Théberge A, Gagnon F, Massi F, Fischi-Gomez E, Gardier R, Haro JLV, Pizzolato M, Caruyer E, Thiran JP. Tractography passes the test: Results from the diffusion-simulated connectivity (disco) challenge. Neuroimage 2023; 277:120231. [PMID: 37330025 PMCID: PMC10771037 DOI: 10.1016/j.neuroimage.2023.120231] [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: 03/02/2023] [Revised: 05/12/2023] [Accepted: 06/14/2023] [Indexed: 06/19/2023] Open
Abstract
Estimating structural connectivity from diffusion-weighted magnetic resonance imaging is a challenging task, partly due to the presence of false-positive connections and the misestimation of connection weights. Building on previous efforts, the MICCAI-CDMRI Diffusion-Simulated Connectivity (DiSCo) challenge was carried out to evaluate state-of-the-art connectivity methods using novel large-scale numerical phantoms. The diffusion signal for the phantoms was obtained from Monte Carlo simulations. The results of the challenge suggest that methods selected by the 14 teams participating in the challenge can provide high correlations between estimated and ground-truth connectivity weights, in complex numerical environments. Additionally, the methods used by the participating teams were able to accurately identify the binary connectivity of the numerical dataset. However, specific false positive and false negative connections were consistently estimated across all methods. Although the challenge dataset doesn't capture the complexity of a real brain, it provided unique data with known macrostructure and microstructure ground-truth properties to facilitate the development of connectivity estimation methods.
Collapse
Affiliation(s)
- Gabriel Girard
- CIBM Center for Biomedical Imaging, Switzerland; Radiology Department, Centre Hospitalier Universitaire Vaudois and University of Lausanne, Lausanne, Switzerland; Signal Processing Laboratory (LTS5), École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland.
| | - Jonathan Rafael-Patiño
- Radiology Department, Centre Hospitalier Universitaire Vaudois and University of Lausanne, Lausanne, Switzerland; Signal Processing Laboratory (LTS5), École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland.
| | - Raphaël Truffet
- Univ Rennes, Inria, CNRS, Inserm, IRISA UMR 6074, Empenn ERL U-1228, Rennes, France
| | - Dogu Baran Aydogan
- A.I. Virtanen Institute for Molecular Sciences, University of Eastern Finland, Kuopio, Finland; Department of Neuroscience and Biomedical Engineering, Aalto University, Espoo, Finland; Department of Psychiatry, Helsinki University Hospital, Helsinki, Finland
| | - Nagesh Adluru
- Waisman Center, University of Wisconsin-Madison, Madison, WI, United States; Department of Radiology, University of Wisconsin-Madison, Madison, WI, United States
| | - Veena A Nair
- Department of Radiology, University of Wisconsin-Madison, Madison, WI, United States
| | - Vivek Prabhakaran
- Department of Radiology, University of Wisconsin-Madison, Madison, WI, United States
| | - Barbara B Bendlin
- Department of Medicine, University of Wisconsin-Madison, Madison, WI, United States
| | - Andrew L Alexander
- Waisman Center, University of Wisconsin-Madison, Madison, WI, United States; Department of Medical Physics, University of Wisconsin-Madison, Madison, WI, United States; Department of Psychiatry, University of Wisconsin-Madison, Madison, WI, United States
| | - Sara Bosticardo
- Diffusion Imaging and Connectivity Estimation (DICE) Lab, Department of Computer Science, University of Verona, Verona, Italy; Translational Imaging in Neurology (ThINk), Department of Biomedical Engineering, University Hospital Basel and University of Basel, Basel, Switzerland
| | - Ilaria Gabusi
- Diffusion Imaging and Connectivity Estimation (DICE) Lab, Department of Computer Science, University of Verona, Verona, Italy; Department of Advanced Biomedical Sciences, University of Naples Federico II, Naples, Italy
| | - Mario Ocampo-Pineda
- Diffusion Imaging and Connectivity Estimation (DICE) Lab, Department of Computer Science, University of Verona, Verona, Italy
| | - Matteo Battocchio
- Diffusion Imaging and Connectivity Estimation (DICE) Lab, Department of Computer Science, University of Verona, Verona, Italy; Sherbrooke Connectivity Imaging Laboratory (SCIL), Department of Computer Science, University of Sherbrooke, Sherbrooke, QC, Canada
| | - Zuzana Piskorova
- Diffusion Imaging and Connectivity Estimation (DICE) Lab, Department of Computer Science, University of Verona, Verona, Italy; Brno Faculty of Electrical Engineering and Communication, Department of mathematics, University of Technology, Brno, Czech Republic
| | - Pietro Bontempi
- Diffusion Imaging and Connectivity Estimation (DICE) Lab, Department of Computer Science, University of Verona, Verona, Italy
| | - Simona Schiavi
- Department of Neuroscience, Rehabilitation, Ophthalmology, Genetics, Maternal and Child Health (DINOGMI), University of Genoa, Genoa, Italy
| | - Alessandro Daducci
- Diffusion Imaging and Connectivity Estimation (DICE) Lab, Department of Computer Science, University of Verona, Verona, Italy
| | | | - Dominika Ciupek
- Sano Centre for Computational Personalised Medicine, Kraków, Poland
| | - Fabian Bogusz
- AGH University of Science and Technology, Kraków, Poland
| | - Tomasz Pieciak
- AGH University of Science and Technology, Kraków, Poland; Laboratorio de Procesado de Imagen (LPI), ETSI Telecomunicación, Universidad de Valladolid, Valladolid, Spain
| | - Matteo Frigo
- Athena Project Team, Centre Inria d'Université Côte d'Azur, France
| | - Sara Sedlar
- Athena Project Team, Centre Inria d'Université Côte d'Azur, France
| | | | - Ivana Kojčić
- Athena Project Team, Centre Inria d'Université Côte d'Azur, France
| | - Mauro Zucchelli
- Athena Project Team, Centre Inria d'Université Côte d'Azur, France
| | - Hiba Laghrissi
- Athena Project Team, Centre Inria d'Université Côte d'Azur, France; Institut de Biologie de Valrose, Université Côte d'Azur, Nice, France
| | - Yang Ji
- Athena Project Team, Centre Inria d'Université Côte d'Azur, France
| | - Rachid Deriche
- Athena Project Team, Centre Inria d'Université Côte d'Azur, France
| | - Kurt G Schilling
- Department of Radiology and Radiological Sciences, Vanderbilt University Medical Center, Nashville, TN, United States
| | - Bennett A Landman
- Department of Radiology and Radiological Sciences, Vanderbilt University Medical Center, Nashville, TN, United States; Department of Electrical and Computer Engineering, Vanderbilt University, Nashville, TN, United States
| | - Alberto Cacciola
- Brain Mapping Lab, Department of Biomedical, Dental Sciences and Morphological and Functional Images, University of Messina, Messina, Italy; Center for Complex Network Intelligence (CCNI), Tsinghua Laboratory of Brain and Intelligence (THBI), Tsinghua University, Beijing, China; Department of Biomedical Engineering, Tsinghua University, Beijing, China
| | - Gianpaolo Antonio Basile
- Brain Mapping Lab, Department of Biomedical, Dental Sciences and Morphological and Functional Images, University of Messina, Messina, Italy
| | - Salvatore Bertino
- Brain Mapping Lab, Department of Biomedical, Dental Sciences and Morphological and Functional Images, University of Messina, Messina, Italy
| | - Nancy Newlin
- Department of Electrical and Computer Engineering, Vanderbilt University, Nashville, TN, United States
| | - Praitayini Kanakaraj
- Department of Electrical and Computer Engineering, Vanderbilt University, Nashville, TN, United States
| | - Francois Rheault
- Department of Electrical and Computer Engineering, Vanderbilt University, Nashville, TN, United States
| | - Patryk Filipiak
- Center for Advanced Imaging Innovation and Research (CAI2R), Department of Radiology, NYU Langone Health, New York, NY, United States
| | - Timothy M Shepherd
- Center for Advanced Imaging Innovation and Research (CAI2R), Department of Radiology, NYU Langone Health, New York, NY, United States
| | - Ying-Chia Lin
- Center for Advanced Imaging Innovation and Research (CAI2R), Department of Radiology, NYU Langone Health, New York, NY, United States
| | - Dimitris G Placantonakis
- Department of Neurosurgery, Perlmutter Cancer Center, Neuroscience Institute, Kimmel Center for Stem Cell Biology, NYU Langone Health, New York, NY, United States
| | - Fernando E Boada
- Department of Radiology, Stanford University, Stanford, CA, United States
| | - Steven H Baete
- Center for Advanced Imaging Innovation and Research (CAI2R), Department of Radiology, NYU Langone Health, New York, NY, United States
| | - Erick Hernández-Gutiérrez
- Sherbrooke Connectivity Imaging Laboratory (SCIL), Department of Computer Science, University of Sherbrooke, Sherbrooke, QC, Canada
| | | | - Ricardo Coronado-Leija
- Center for Advanced Imaging Innovation and Research (CAI2R), Department of Radiology, NYU Langone Health, New York, NY, United States
| | - Pablo Stack-Sánchez
- Computer Science Department, Centro de Investigación en Matemáticas A.C, Guanajuato, México
| | - Luis Concha
- Instituto de Neurobiología, Universidad Nacional Autónoma de México, Juriquilla, Querétaro, México
| | - Maxime Descoteaux
- Sherbrooke Connectivity Imaging Laboratory (SCIL), Department of Computer Science, University of Sherbrooke, Sherbrooke, QC, Canada
| | - Sina Mansour L
- Department of Biomedical Engineering, The University of Melbourne, Parkville, Victoria, Australia; Melbourne Neuropsychiatry Centre, Department of Psychiatry, The University of Melbourne and Melbourne Health, Parkville, Victoria, Australia
| | - Caio Seguin
- Melbourne Neuropsychiatry Centre, Department of Psychiatry, The University of Melbourne and Melbourne Health, Parkville, Victoria, Australia; School of Biomedical Engineering, The University of Sydney, Sydney, Australia; Department of Psychological and Brain Sciences, Indiana University, Bloomington, IN, United States
| | - Andrew Zalesky
- Department of Biomedical Engineering, The University of Melbourne, Parkville, Victoria, Australia; Melbourne Neuropsychiatry Centre, Department of Psychiatry, The University of Melbourne and Melbourne Health, Parkville, Victoria, Australia
| | - Kenji Marshall
- Signal Processing Laboratory (LTS5), École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland; McGill University, Montréal, QC, Canada
| | - Erick J Canales-Rodríguez
- Signal Processing Laboratory (LTS5), École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
| | - Ye Wu
- Department of Radiology and Biomedical Research Imaging Center (BRIC), The University of North Carolina at Chapel Hill, Chapel Hill, NC, United States; School of Computer Science and Engineering, Nanjing University of Science and Technology, Nanjing, China
| | - Sahar Ahmad
- Department of Radiology and Biomedical Research Imaging Center (BRIC), The University of North Carolina at Chapel Hill, Chapel Hill, NC, United States
| | - Pew-Thian Yap
- Department of Radiology and Biomedical Research Imaging Center (BRIC), The University of North Carolina at Chapel Hill, Chapel Hill, NC, United States
| | - Antoine Théberge
- Sherbrooke Connectivity Imaging Laboratory (SCIL), Department of Computer Science, University of Sherbrooke, Sherbrooke, QC, Canada
| | - Florence Gagnon
- Sherbrooke Connectivity Imaging Laboratory (SCIL), Department of Computer Science, University of Sherbrooke, Sherbrooke, QC, Canada
| | - Frédéric Massi
- Sherbrooke Connectivity Imaging Laboratory (SCIL), Department of Computer Science, University of Sherbrooke, Sherbrooke, QC, Canada
| | - Elda Fischi-Gomez
- CIBM Center for Biomedical Imaging, Switzerland; Radiology Department, Centre Hospitalier Universitaire Vaudois and University of Lausanne, Lausanne, Switzerland; Signal Processing Laboratory (LTS5), École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
| | - Rémy Gardier
- Signal Processing Laboratory (LTS5), École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
| | - Juan Luis Villarreal Haro
- Signal Processing Laboratory (LTS5), École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
| | - Marco Pizzolato
- Signal Processing Laboratory (LTS5), École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland; Department of Applied Mathematics and Computer Science, Technical University of Denmark, Kgs. Lyngby, Denmark
| | - Emmanuel Caruyer
- Univ Rennes, Inria, CNRS, Inserm, IRISA UMR 6074, Empenn ERL U-1228, Rennes, France
| | - Jean-Philippe Thiran
- CIBM Center for Biomedical Imaging, Switzerland; Radiology Department, Centre Hospitalier Universitaire Vaudois and University of Lausanne, Lausanne, Switzerland; Signal Processing Laboratory (LTS5), École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
| |
Collapse
|
5
|
Wu D, Turnbill V, Lee HH, Wang X, Ba R, Walczak P, Martin LJ, Fieremans E, Novikov DS, Northington FJ, Zhang J. In vivo Mapping of Cellular Resolution Neuropathology in Brain Ischemia by Diffusion MRI. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.08.08.552374. [PMID: 37609182 PMCID: PMC10441332 DOI: 10.1101/2023.08.08.552374] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 08/24/2023]
Abstract
Non-invasive mapping of cellular pathology can provide critical diagnostic and prognostic information. Recent developments in diffusion MRI have produced new tools for examining tissue microstructure at a level well below the imaging resolution. Here, we report the use of diffusion time ( t )-dependent diffusion kurtosis imaging ( t DKI) to simultaneously assess the morphology and transmembrane permeability of cells and their processes in the context of pathological changes in hypoxic-ischemic brain (HI) injury. Through Monte Carlo simulations and cell culture organoid imaging, we demonstrate feasibility in measuring effective size and permeability changes based on the peak and tail of t DKI curves. In a mouse model of HI, in vivo imaging at 11.7T detects a marked shift of the t DKI peak to longer t in brain edema, suggesting swelling and beading associated with the astrocytic processes and neuronal neurites. Furthermore, we observed a faster decrease of the t DKI tail in injured brain regions, reflecting increased membrane permeability that was associated with upregulated water exchange upon astrocyte activation at acute stage as well as necrosis with disrupted membrane integrity at subacute stage. Such information, unavailable with conventional diffusion MRI at a single t, can predict salvageable tissues. For a proof-of-concept, t DKI at 3T on an ischemic stroke patient suggested increased membrane permeability in the stroke region. This work therefore demonstrates the potential of t DKI for in vivo detection of the pathological changes in microstructural morphology and transmembrane permeability after ischemic injury using a clinically translatable protocol.
Collapse
|
6
|
Ulloa P, Methot V, Wottschel V, Koch MA. Extra-axonal contribution to double diffusion encoding-based pore size estimates in the corticospinal tract. MAGMA (NEW YORK, N.Y.) 2023; 36:589-612. [PMID: 36745290 PMCID: PMC10468962 DOI: 10.1007/s10334-022-01058-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/25/2022] [Revised: 12/14/2022] [Accepted: 12/19/2022] [Indexed: 02/07/2023]
Abstract
OBJECTIVE To study the origin of compartment size overestimation in double diffusion encoding MRI (DDE) in vivo experiments in the human corticospinal tract. Here, the extracellular space is hypothesized to be the origin of the DDE signal. By exploiting the DDE sensitivity to pore shape, it could be possible to identify the origin of the measured signal. The signal difference between parallel and perpendicular diffusion gradient orientation can indicate if a compartment is regular or eccentric in shape. As extracellular space can be considered an eccentric compartment, a positive difference would mean a high contribution to the compartment size estimates. MATERIALS AND METHODS Computer simulations using MISST and in vivo experiments in eight healthy volunteers were performed. DDE experiments using a double spin-echo preparation with eight perpendicular directions were measured in vivo. The difference between parallel and perpendicular gradient orientations was analyzed using a Wilcoxon signed-rank test and a Mann-Whitney U test. RESULTS Simulations and MR experiments showed a statistically significant difference between parallel and perpendicular diffusion gradient orientation signals ([Formula: see text]). CONCLUSION The results suggest that the DDE-based size estimate may be considerably influenced by the extra-axonal compartment. However, the experimental results are also consistent with purely intra-axonal contributions in combination with a large fiber orientation dispersion.
Collapse
Affiliation(s)
- Patricia Ulloa
- Institute of Medical Engineering, University of Luebeck, Ratzeburger Allee 160, 23562 Luebeck, Germany
| | - Vincent Methot
- Institute of Medical Engineering, University of Luebeck, Ratzeburger Allee 160, 23562 Luebeck, Germany
| | - Viktor Wottschel
- Department of Radiology and Nuclear Medicine, Amsterdam University Medical Centers, De Boelelaan 1117, 1081, Amsterdam, The Netherlands
| | - Martin A. Koch
- Institute of Medical Engineering, University of Luebeck, Ratzeburger Allee 160, 23562 Luebeck, Germany
| |
Collapse
|
7
|
Chowdhury R, Wan J, Gardier R, Rafael-Patino J, Thiran JP, Gibou F, Mukherjee A. Molecular imaging with aquaporin-based reporter genes: quantitative considerations from Monte Carlo diffusion simulations. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.06.09.544324. [PMID: 37333205 PMCID: PMC10274877 DOI: 10.1101/2023.06.09.544324] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/20/2023]
Abstract
Aquaporins provide a new class of genetic tools for imaging molecular activity in deep tissues by increasing the rate of cellular water diffusion, which generates magnetic resonance contrast. However, distinguishing aquaporin contrast from the tissue background is challenging because water diffusion is also influenced by structural factors such as cell size and packing density. Here, we developed and experimentally validated a Monte Carlo model to analyze how cell radius and intracellular volume fraction quantitatively affect aquaporin signals. We demonstrated that a differential imaging approach based on time-dependent changes in diffusivity can improve specificity by unambiguously isolating aquaporin-driven contrast from the tissue background. Finally, we used Monte Carlo simulations to analyze the connection between diffusivity and the percentage of cells engineered to express aquaporin, and established a simple mapping that accurately determined the volume fraction of aquaporin-expressing cells in mixed populations. This study creates a framework for broad applications of aquaporins, particularly in biomedicine and in vivo synthetic biology, where quantitative methods to measure the location and performance of genetic devices in whole vertebrates are necessary.
Collapse
Affiliation(s)
- Rochishnu Chowdhury
- Department of Mechanical Engineering, University of California, Santa Barbara, CA 93106, USA
| | - Jinyang Wan
- Department of Chemistry, University of California, Santa Barbara, CA 93106, USA
| | - Remy Gardier
- Signal Processing Laboratory (LTS5), Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
| | - Jonathan Rafael-Patino
- Signal Processing Laboratory (LTS5), Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
- Radiology Department, Centre Hospitalier Universitaire Vaudois (CHUV) and University of Lausanne (UNIL), Lausanne, Switzerland
| | - Jean-Philippe Thiran
- Signal Processing Laboratory (LTS5), Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
- Radiology Department, Centre Hospitalier Universitaire Vaudois (CHUV) and University of Lausanne (UNIL), Lausanne, Switzerland
| | - Frederic Gibou
- Department of Mechanical Engineering, University of California, Santa Barbara, CA 93106, USA
- Department of Computer Science, University of California, Santa Barbara, CA 93106, USA
| | - Arnab Mukherjee
- Department of Chemistry, University of California, Santa Barbara, CA 93106, USA
- Department of Chemical Engineering, University of California, Santa Barbara, CA 93106, USA
- Biomolecular Science and Engineering, University of California, Santa Barbara, CA 93106, USA
- Biological Engineering, University of California, Santa Barbara, CA 93106, USA
- Neuroscience Research Institute, University of California, Santa Barbara, CA 93106, USA
| |
Collapse
|
8
|
Lee HH, Tian Q, Sheft M, Coronado-Leija R, Ramos-Llorden G, Abdollahzadeh A, Fieremans E, Novikov DS, Huang SY. The influence of axonal beading and undulation on axonal diameter mapping. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.04.19.537494. [PMID: 37131702 PMCID: PMC10153226 DOI: 10.1101/2023.04.19.537494] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/04/2023]
Abstract
We consider the effect of non-cylindrical axonal shape on axonal diameter mapping with diffusion MRI. Practical sensitivity to axon diameter is attained at strong diffusion weightings b , where the deviation from the 1 / b scaling yields the finite transverse diffusivity, which is then translated into axon diameter. While axons are usually modeled as perfectly straight, impermeable cylinders, the local variations in diameter (caliber variation or beading) and direction (undulation) have been observed in microscopy data of human axons. Here we quantify the influence of cellular-level features such as caliber variation and undulation on axon diameter estimation. For that, we simulate the diffusion MRI signal in realistic axons segmented from 3-dimensional electron microscopy of a human brain sample. We then create artificial fibers with the same features and tune the amplitude of their caliber variations and undulations. Numerical simulations of diffusion in fibers with such tunable features show that caliber variations and undulations result in under- and over-estimation of axon diameters, correspondingly; this bias can be as large as 100%. Given that increased axonal beading and undulations have been observed in pathological tissues, such as traumatic brain injury and ischemia, the interpretation of axon diameter alterations in pathology may be significantly confounded.
Collapse
Affiliation(s)
- Hong-Hsi Lee
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA 02129,USA
- Harvard Medical School, Boston, MA 02115, USA
| | - Qiyuan Tian
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA 02129,USA
- Harvard Medical School, Boston, MA 02115, USA
| | - Maxina Sheft
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA 02129,USA
- Harvard-MIT Health Sciences and Technology, Cambridge, MA 02139, USA
| | - Ricardo Coronado-Leija
- Center for Biomedical Imaging, Department of Radiology, New York University School of Medicine, New York, NY 10016, USA
- Center for Advanced Imaging Innovation and Research (CAI2R), New York University School of Medicine, New York, NY 10016, USA
| | - Gabriel Ramos-Llorden
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA 02129,USA
- Harvard Medical School, Boston, MA 02115, USA
| | - Ali Abdollahzadeh
- Center for Biomedical Imaging, Department of Radiology, New York University School of Medicine, New York, NY 10016, USA
- Center for Advanced Imaging Innovation and Research (CAI2R), New York University School of Medicine, New York, NY 10016, USA
| | - Els Fieremans
- Center for Biomedical Imaging, Department of Radiology, New York University School of Medicine, New York, NY 10016, USA
- Center for Advanced Imaging Innovation and Research (CAI2R), New York University School of Medicine, New York, NY 10016, USA
| | - Dmitry S. Novikov
- Center for Biomedical Imaging, Department of Radiology, New York University School of Medicine, New York, NY 10016, USA
- Center for Advanced Imaging Innovation and Research (CAI2R), New York University School of Medicine, New York, NY 10016, USA
| | - Susie Y. Huang
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA 02129,USA
- Harvard Medical School, Boston, MA 02115, USA
| |
Collapse
|
9
|
Chakwizira A, Westin C, Brabec J, Lasič S, Knutsson L, Szczepankiewicz F, Nilsson M. Diffusion MRI with pulsed and free gradient waveforms: Effects of restricted diffusion and exchange. NMR IN BIOMEDICINE 2023; 36:e4827. [PMID: 36075110 PMCID: PMC10078514 DOI: 10.1002/nbm.4827] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/16/2021] [Revised: 08/27/2022] [Accepted: 09/06/2022] [Indexed: 05/06/2023]
Abstract
Monitoring time dependence with diffusion MRI yields observables sensitive to compartment sizes (restricted diffusion) and membrane permeability (water exchange). However, restricted diffusion and exchange have opposite effects on the diffusion-weighted signal, which can lead to errors in parameter estimates. In this work, we propose a signal representation that incorporates the effects of both restricted diffusion and exchange up to second order in b-value and is compatible with gradient waveforms of arbitrary shape. The representation features mappings from a gradient waveform to two scalars that separately control the sensitivity to restriction and exchange. We demonstrate that these scalars span a two-dimensional space that can be used to choose waveforms that selectively probe restricted diffusion or exchange, eliminating the correlation between the two phenomena. We found that waveforms with specific but unconventional shapes provide an advantage over conventional pulsed and oscillating gradient acquisitions. We also show that parametrization of waveforms into a two-dimensional space can be used to understand protocols from other approaches that probe restricted diffusion and exchange. For example, we found that the variation of mixing time in filter-exchange imaging corresponds to variation of our exchange-weighting scalar at a fixed value of the restriction-weighting scalar. The proposed signal representation was evaluated using Monte Carlo simulations in identical parallel cylinders with hexagonal and random packing as well as parallel cylinders with gamma-distributed radii. Results showed that the approach is sensitive to sizes in the interval 4-12 μm and exchange rates in the simulated range of 0 to 20 s - 1 , but also that there is a sensitivity to the extracellular geometry. The presented theory constitutes a simple and intuitive description of how restricted diffusion and exchange influence the signal as well as a guide to protocol design capable of separating the two effects.
Collapse
Affiliation(s)
- Arthur Chakwizira
- Department of Medical Radiation Physics, LundLund UniversityLundSweden
| | - Carl‐Fredrik Westin
- Department of Radiology, Brigham and Women's HospitalHarvard Medical SchoolBostonMassachusettsUSA
| | - Jan Brabec
- Department of Medical Radiation Physics, LundLund UniversityLundSweden
| | - Samo Lasič
- Danish Research Centre for Magnetic Resonance, Centre for Functional and Diagnostic Imaging and ResearchCopenhagen University Hospital ‐ Amager and HvidovreCopenhagenDenmark
- Random Walk Imaging ABLundSweden
| | - Linda Knutsson
- Department of Medical Radiation Physics, LundLund UniversityLundSweden
- Russell H. Morgan Department of Radiology and Radiological ScienceJohns Hopkins University School of MedicineBaltimoreMarylandUSA
- F. M. Kirby Research Center for Functional Brain ImagingKennedy Krieger InstituteBaltimoreMarylandUSA
| | | | - Markus Nilsson
- Department of Clinical Sciences Lund, RadiologyLund UniversityLundSweden
| |
Collapse
|
10
|
Agdestein SD, Tran TN, Li JR. Practical computation of the diffusion MRI signal based on Laplace eigenfunctions: permeable interfaces. NMR IN BIOMEDICINE 2022; 35:e4646. [PMID: 34796990 DOI: 10.1002/nbm.4646] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/29/2021] [Revised: 09/29/2021] [Accepted: 10/13/2021] [Indexed: 06/13/2023]
Abstract
The complex transverse water proton magnetization subject to diffusion-encoding magnetic field gradient pulses in a heterogeneous medium such as brain tissue can be modeled by the Bloch-Torrey partial differential equation. The spatial integral of the solution of this equation in realistic geometry provides a gold-standard reference model for the diffusion MRI signal arising from different tissue micro-structures of interest. A closed form representation of this reference diffusion MRI signal, called matrix formalism, which makes explicit the link between the Laplace eigenvalues and eigenfunctions of the tissue geometry and its diffusion MRI signal, was derived 20 years ago. In addition, once the Laplace eigendecomposition has been computed and saved, the diffusion MRI signal can be calculated for arbitrary diffusion-encoding sequences and b-values at negligible additional cost. In a previous publication, we presented a simulation framework that we implemented inside the MATLAB-based diffusion MRI simulator SpinDoctor that efficiently computes the matrix formalism representation for biological cells subject to impermeable membrane boundary conditions. In this work, we extend our simulation framework to include geometries that contain permeable cell membranes. We describe the new computational techniques that allowed this generalization and we analyze the effects of the magnitude of the permeability coefficient on the eigendecomposition of the diffusion and Bloch-Torrey operators. This work is another step in bringing advanced mathematical tools and numerical method development to the simulation and modeling of diffusion MRI.
Collapse
Affiliation(s)
| | | | - Jing-Rebecca Li
- INRIA Saclay-Equipe DEFI, CMAP, Ecole Polytechnique, Palaiseau, France
| |
Collapse
|
11
|
Connecting macroscopic diffusion metrics of cardiac diffusion tensor imaging and microscopic myocardial structures based on simulation. Med Image Anal 2022; 77:102325. [DOI: 10.1016/j.media.2021.102325] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2019] [Revised: 11/26/2021] [Accepted: 11/29/2021] [Indexed: 11/20/2022]
|
12
|
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
|
13
|
Cheng H, Calhoun V. Exploring microstructure with diffusion-weighted imaging: From acquisition to modeling. J Neurosci Methods 2021; 363:109335. [PMID: 34450220 DOI: 10.1016/j.jneumeth.2021.109335] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- Hu Cheng
- Imaging Research Facility, Department of Psychological and Brain Sciences, Indiana University, Bloomington, IN 47405, USA.
| | - Vince Calhoun
- Tri-institutional center for Translational Research in Neuroimaging and Data Science (TReNDS) Center, Georgia Institute of Technology, Emory University, Georgia State University, Atlanta 30303, GA, USA.
| |
Collapse
|
14
|
Kerkelä L, Nery F, Callaghan R, Zhou F, Gyori NG, Szczepankiewicz F, Palombo M, Parker GJM, Zhang H, Hall MG, Clark CA. Comparative analysis of signal models for microscopic fractional anisotropy estimation using q-space trajectory encoding. Neuroimage 2021; 242:118445. [PMID: 34375753 DOI: 10.1016/j.neuroimage.2021.118445] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2021] [Revised: 07/06/2021] [Accepted: 08/02/2021] [Indexed: 12/12/2022] Open
Abstract
Microscopic diffusion anisotropy imaging using diffusion-weighted MRI and multidimensional diffusion encoding is a promising method for quantifying clinically and scientifically relevant microstructural properties of neural tissue. Several methods for estimating microscopic fractional anisotropy (µFA), a normalized measure of microscopic diffusion anisotropy, have been introduced but the differences between the methods have received little attention thus far. In this study, the accuracy and precision of µFA estimation using q-space trajectory encoding and different signal models were assessed using imaging experiments and simulations. Three healthy volunteers and a microfibre phantom were imaged with five non-zero b-values and gradient waveforms encoding linear and spherical b-tensors. Since the ground-truth µFA was unknown in the imaging experiments, Monte Carlo random walk simulations were performed using axon-mimicking fibres for which the ground truth was known. Furthermore, parameter bias due to time-dependent diffusion was quantified by repeating the simulations with tuned waveforms, which have similar power spectra, and with triple diffusion encoding, which, unlike q-space trajectory encoding, is not based on the assumption of time-independent diffusion. The truncated cumulant expansion of the powder-averaged signal, gamma-distributed diffusivities assumption, and q-space trajectory imaging, a generalization of the truncated cumulant expansion to individual signals, were used to estimate µFA. The gamma-distributed diffusivities assumption consistently resulted in greater µFA values than the second order cumulant expansion, 0.1 greater when averaged over the whole brain. In the simulations, the generalized cumulant expansion provided the most accurate estimates. Importantly, although time-dependent diffusion caused significant overestimation of µFA using all the studied methods, the simulations suggest that the resulting bias in µFA is less than 0.1 in human white matter.
Collapse
Affiliation(s)
- Leevi Kerkelä
- UCL Great Ormond Street Institute of Child Health, University College London, London, UK.
| | - Fabio Nery
- UCL Great Ormond Street Institute of Child Health, University College London, London, UK
| | - Ross Callaghan
- UCL Centre for Medical Image Computing, University College London, London, UK
| | - Fenglei Zhou
- UCL Centre for Medical Image Computing, University College London, London, UK; UCL School of Pharmacy, University College London, London, UK
| | - Noemi G Gyori
- UCL Centre for Medical Image Computing, University College London, London, UK; UCL Great Ormond Street Institute of Child Health, University College London, London, UK
| | - Filip Szczepankiewicz
- Department of Radiology, Brigham and Women's Hospital, Boston, Massachusetts, US; Harvard Medical School, Boston, Massachusetts, US; Clinical Sciences Lund, Lund University, Lund, Sweden
| | - Marco Palombo
- UCL Centre for Medical Image Computing, University College London, London, UK
| | - Geoff J M Parker
- UCL Centre for Medical Image Computing, University College London, London, UK; Bioxydyn Limited, Manchester, UK; UCL Queen Square Institute of Neurology, University College London, London, UK
| | - Hui Zhang
- UCL Centre for Medical Image Computing, University College London, London, UK
| | - Matt G Hall
- UCL Great Ormond Street Institute of Child Health, University College London, London, UK; National Physical Laboratory, Teddington, UK
| | - Chris A Clark
- UCL Great Ormond Street Institute of Child Health, University College London, London, UK
| |
Collapse
|
15
|
Karimi D, Vasung L, Jaimes C, Machado-Rivas F, Khan S, Warfield SK, Gholipour A. A machine learning-based method for estimating the number and orientations of major fascicles in diffusion-weighted magnetic resonance imaging. Med Image Anal 2021; 72:102129. [PMID: 34182203 PMCID: PMC8320341 DOI: 10.1016/j.media.2021.102129] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2020] [Revised: 05/26/2021] [Accepted: 05/28/2021] [Indexed: 12/29/2022]
Abstract
Accurate modeling of diffusion-weighted magnetic resonance imaging measurements is necessary for accurate brain connectivity analysis. Existing methods for estimating the number and orientations of fascicles in an imaging voxel either depend on non-convex optimization techniques that are sensitive to initialization and measurement noise, or are prone to predicting spurious fascicles. In this paper, we propose a machine learning-based technique that can accurately estimate the number and orientations of fascicles in a voxel. Our method can be trained with either simulated or real diffusion-weighted imaging data. Our method estimates the angle to the closest fascicle for each direction in a set of discrete directions uniformly spread on the unit sphere. This information is then processed to extract the number and orientations of fascicles in a voxel. On realistic simulated phantom data with known ground truth, our method predicts the number and orientations of crossing fascicles more accurately than several classical and machine learning methods. It also leads to more accurate tractography. On real data, our method is better than or compares favorably with other methods in terms of robustness to measurement down-sampling and also in terms of expert quality assessment of tractography results.
Collapse
Affiliation(s)
- Davood Karimi
- Department of Radiology, Boston Children's Hospital, Harvard Medical School, Boston, MA, USA.
| | - Lana Vasung
- Department of Pediatrics at Boston Children's Hospital, and Harvard Medical School, Boston, Massachusetts, USA
| | - Camilo Jaimes
- Department of Radiology, Boston Children's Hospital, Harvard Medical School, Boston, MA, USA
| | - Fedel Machado-Rivas
- Department of Radiology, Boston Children's Hospital, Harvard Medical School, Boston, MA, USA
| | - Shadab Khan
- Department of Radiology, Boston Children's Hospital, Harvard Medical School, Boston, MA, USA
| | - Simon K Warfield
- Department of Radiology, Boston Children's Hospital, Harvard Medical School, Boston, MA, USA
| | - Ali Gholipour
- Department of Radiology, Boston Children's Hospital, Harvard Medical School, Boston, MA, USA
| |
Collapse
|
16
|
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: 19] [Impact Index Per Article: 6.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
|
17
|
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
|