1
|
Coudert T, Delphin A, Barrier A, Legris L, Warnking JM, Lamalle L, Doneva M, Lemasson B, Barbier EL, Christen T. Relaxometry and contrast-free cerebral microvascular quantification using balanced steady-state free precession MR fingerprinting. Magn Reson Med 2025; 94:302-316. [PMID: 39825561 DOI: 10.1002/mrm.30434] [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/07/2024] [Revised: 12/11/2024] [Accepted: 12/31/2024] [Indexed: 01/20/2025]
Abstract
PURPOSE This study proposes a novel, contrast-free Magnetic Resonance Fingerprinting (MRF) method using balanced Steady-State Free Precession (bSSFP) sequences for the quantification of cerebral blood volume (CBV), vessel radius (R), and relaxometry parameters (T 1 $$ {}_1 $$ , T 2 $$ {}_2 $$ , T 2 $$ {}_2 $$ *) in the brain. METHODS The technique leverages the sensitivity of bSSFP sequences to intra-voxel frequency distributions in both transient and steady-state regimes. A dictionary-matching process is employed, using simulations of realistic mouse microvascular networks to generate the MRF dictionary. The method is validated through in silico and in vivo experiments on six healthy subjects, comparing results with standard MRF methods and literature values. RESULTS The proposed method shows strong correlation and agreement with standard MRF methods for T 1 $$ {}_1 $$ and T 2 $$ {}_2 $$ values. High-resolution maps provide detailed visualizations of CBV and microvascular structures, highlighting differences in white matter (WM) and gray matter (GM) regions. The measured GM/WM ratio for CBV is 1.91, consistent with literature values. CONCLUSION This contrast-free bSSFP-based MRF method offers an new approach for quantifying CBV, vessel radius, and relaxometry parameters. Further validation against DSC imaging and clinical studies in pathological conditions is warranted to confirm its clinical utility.
Collapse
Affiliation(s)
- Thomas Coudert
- Université Grenoble Alpes, INSERM, U1216, Grenoble Institute Neurosciences, GIN, Grenoble, France
| | - Aurélien Delphin
- Université Grenoble Alpes, INSERM, US17, CNRS, UAR3552, CHU Grenoble Alpes, IRMaGe, Grenoble, France
| | - Antoine Barrier
- Université Grenoble Alpes, INSERM, U1216, Grenoble Institute Neurosciences, GIN, Grenoble, France
| | - Loïc Legris
- Université Grenoble Alpes, INSERM, U1216, Grenoble Institute Neurosciences, GIN, Grenoble, France
- Université Grenoble Alpes, Stroke Unit, Department of Neurology, CHU Grenoble Alpes, Grenoble, France
| | - Jan M Warnking
- Université Grenoble Alpes, INSERM, U1216, Grenoble Institute Neurosciences, GIN, Grenoble, France
- Université Grenoble Alpes, INSERM, US17, CNRS, UAR3552, CHU Grenoble Alpes, IRMaGe, Grenoble, France
| | - Laurent Lamalle
- Université Grenoble Alpes, INSERM, US17, CNRS, UAR3552, CHU Grenoble Alpes, IRMaGe, Grenoble, France
| | | | - Benjamin Lemasson
- Université Grenoble Alpes, INSERM, U1216, Grenoble Institute Neurosciences, GIN, Grenoble, France
| | - Emmanuel L Barbier
- Université Grenoble Alpes, INSERM, U1216, Grenoble Institute Neurosciences, GIN, Grenoble, France
| | - Thomas Christen
- Université Grenoble Alpes, INSERM, U1216, Grenoble Institute Neurosciences, GIN, Grenoble, France
| |
Collapse
|
2
|
Oliveira R, Raynaud Q, Jelescu I, Kiselev V, Kirilina E, Lutti A. In Vivo Characterization of Magnetic Inclusions in the Subcortex From Nonexponential Transverse Relaxation Decay. NMR IN BIOMEDICINE 2025; 38:e70051. [PMID: 40325973 PMCID: PMC12053162 DOI: 10.1002/nbm.70051] [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: 10/28/2024] [Revised: 04/16/2025] [Accepted: 04/17/2025] [Indexed: 05/07/2025]
Abstract
According to theoretical studies, MRI signal decay due to transverse relaxation in brain tissue with magnetic inclusions (e.g., blood vessels and iron-rich cells) is expected to follow a transition from Gaussian behaviour at short echo times to exponential behaviour at longer times. The decay parameters carry information about the inclusions (e.g., size and volume fraction) and provide unique insights into brain microstructure. However, gradient-echo decays typically only capture the long-time exponential behaviour. We provide experimental evidence of nonexponential transverse relaxation decay in human subcortical grey matter from in vivo MRI data acquired at 3 T, allowing the subsequent characterization of the magnetic inclusions. Gradient-echo data were collected with short interecho spacings, minimal echo time (1.25 ms) and novel acquisition strategies to mitigate motion and cardiac-induced effects. The data were fitted using exponential and nonexponential models that describe the impact of magnetic inclusions on the MRI signal. Nonexponential models provided superior fits. The strongest deviations from exponential were detected in the substantia nigra and globus pallidus. Numerical simulations of the signal decay from histological maps of iron concentration in the substantia nigra replicated the experimental data, highlighting that non-haem iron can be at the source of the nonexponential decay. To investigate the potential of nonexponential decays to characterize brain microstructure, we estimated the properties of the underlying inclusions using two analytical models. Under the static dephasing regime, the magnetic susceptibility and volume fractions of the inclusions ranged between 1.8-4 and 0.02-0.04 ppm, respectively. Alternatively, under the diffusion narrowing regime, the typical inclusion size was ~2.4 μm. Both simulations and experimental data suggest an intermediate regime with a non-negligible effect of water diffusion. Nonexponential transverse relaxation decay allows to characterize the spatial distribution of magnetic material within subcortical tissue with increased specificity, with potential applications for Parkinson's disease and other pathologies.
Collapse
Affiliation(s)
- Rita Oliveira
- Laboratory for Research in Neuroimaging, Department of Clinical NeuroscienceLausanne University Hospital and University of LausanneLausanneSwitzerland
| | - Quentin Raynaud
- Laboratory for Research in Neuroimaging, Department of Clinical NeuroscienceLausanne University Hospital and University of LausanneLausanneSwitzerland
| | - Ileana Jelescu
- Department of RadiologyLausanne University Hospital (CHUV) and University of LausanneLausanneSwitzerland
| | - Valerij G. Kiselev
- Division of Medical Physics, Department of Radiology, Faculty of MedicineUniversity of FreiburgFreiburgGermany
| | - Evgeniya Kirilina
- Department of NeurophysicsMax Planck Institute for Human Cognitive and Brain SciencesLeipzigGermany
| | - Antoine Lutti
- Laboratory for Research in Neuroimaging, Department of Clinical NeuroscienceLausanne University Hospital and University of LausanneLausanneSwitzerland
| |
Collapse
|
3
|
Şişman M, Nguyen TD, Roberts AG, Romano DJ, Dimov AV, Kovanlikaya I, Spincemaille P, Wang Y. Microstructure-Informed Myelin Mapping (MIMM) from routine multi-echo gradient echo data using multiscale physics modeling of iron and myelin effects and QSM. Magn Reson Med 2025; 93:1499-1515. [PMID: 39552224 PMCID: PMC11910495 DOI: 10.1002/mrm.30369] [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: 06/27/2024] [Revised: 10/08/2024] [Accepted: 10/21/2024] [Indexed: 11/19/2024]
Abstract
PURPOSE Myelin quantification is used in the study of demyelination in neurodegenerative diseases. A novel noninvasive MRI method, Microstructure-Informed Myelin Mapping (MIMM), is proposed to quantify the myelin volume fraction (MVF) from a routine multi-gradient echo sequence (mGRE) using a multiscale biophysical signal model of the effects of microstructural myelin and iron. THEORY AND METHODS In MIMM, the effects of myelin are modeled based on the Hollow Cylinder Fiber Model accounting for anisotropy, while iron is considered as an isotropic paramagnetic point source. This model is used to create a dictionary of mGRE magnitude signal evolution and total voxel susceptibility using finite elements of size 0.2 μm. Next, voxel-by-voxel stochastic matching pursuit between acquired mGRE data (magnitude+QSM) and the pre-computed dictionary generates quantitative MVF and iron susceptibility maps. Dictionary matching was evaluated under three conditions: (1) without fiber orientation (basic), (2) with fiber orientation obtained using DTI, and (3) with fiber orientation obtained using an atlas (atlas). MIMM was compared with the three-pool complex fitting (3PCF) using T2-relaxometry myelin water fraction (MWF) map as reference. RESULTS The DTI MIMM and atlas MIMM approaches were equally effective in reducing the overestimation of MVF in certain white matter tracts observed in the basic MIMM approach, and they both showed good agreement with T2-relaxometry MWF. MIMM MVF reduced myelin overestimation of globus pallidus observed in 3PCF MWF. CONCLUSION MIMM processing of mGRE data can provide MVF maps from routine clinical scans without requiring special sequences.
Collapse
Affiliation(s)
- Mert Şişman
- Department of Electrical and Computer Engineering, Cornell University, Ithaca, New York
- Department of Radiology, Weill Cornel Medicine, New York, New York
| | - Thanh D. Nguyen
- Department of Radiology, Weill Cornel Medicine, New York, New York
| | - Alexandra G. Roberts
- Department of Electrical and Computer Engineering, Cornell University, Ithaca, New York
- Department of Radiology, Weill Cornel Medicine, New York, New York
| | - Dominick J. Romano
- Department of Radiology, Weill Cornel Medicine, New York, New York
- Meinig School of Biomedical Engineering, Cornell University, Ithaca, New York
| | - Alexey V. Dimov
- Department of Radiology, Weill Cornel Medicine, New York, New York
| | | | | | - Yi Wang
- Department of Radiology, Weill Cornel Medicine, New York, New York
- Meinig School of Biomedical Engineering, Cornell University, Ithaca, New York
| |
Collapse
|
4
|
Sandgaard A, Jespersen S. Predicting Mesoscopic Larmor Frequency Shifts in White Matter With Diffusion MRI-A Monte Carlo Study in Axonal Phantoms. NMR IN BIOMEDICINE 2025; 38:e70004. [PMID: 39933490 PMCID: PMC11813543 DOI: 10.1002/nbm.70004] [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: 07/11/2024] [Revised: 12/18/2024] [Accepted: 01/14/2025] [Indexed: 02/13/2025]
Abstract
Magnetic susceptibility MRI offers potential insights into the chemical composition and microstructural organization of tissue. However, estimating magnetic susceptibility in white matter is challenging due to anisotropic subvoxel Larmor frequency shifts caused by axonal microstructure relative to the B0 field orientation. Recent biophysical models have analytically described how axonal microstructure influences the Larmor frequency shifts, relating these shifts to a mesoscopically averaged magnetic field that depends on the axons' fiber orientation distribution function (fODF), typically estimated using diffusion MRI. This study is aimed at validating the use of MRI to estimate mesoscopic magnetic fields and determining whether diffusion MRI can faithfully estimate the orientation dependence of the Larmor frequency shift in realistic axonal microstructure. To achieve this, we developed a framework for performing Monte Carlo simulations of MRI signals in mesoscopically sized white matter axon substrates segmented with electron microscopy. Our simulations demonstrated that with careful experimental design, it is feasible to estimate mesoscopic magnetic fields. Additionally, the fODF estimated by the standard model of diffusion in white matter could predict the orientation dependence of the mesoscopic Larmor frequency shift. We also found that incorporating the intra-axonal axial kurtosis into the standard model could explain a significant amount of signal variance, thereby improving the estimation of the Larmor frequency shift. This factor should not be neglected when fitting the standard model.
Collapse
Affiliation(s)
- Anders Dyhr Sandgaard
- Center of Functionally Integrative Neuroscience, Department of Clinical MedicineAarhus UniversityAarhusDenmark
| | - Sune Nørhøj Jespersen
- Center of Functionally Integrative Neuroscience, Department of Clinical MedicineAarhus UniversityAarhusDenmark
- Department of Physics and AstronomyAarhus UniversityAarhusDenmark
| |
Collapse
|
5
|
Winther S, Lundell H, Rafael-Patiño J, Andersson M, Thiran JP, Dyrby TB. Susceptibility-induced internal gradients reveal axon morphology and cause anisotropic effects in the diffusion-weighted MRI signal. Sci Rep 2024; 14:29636. [PMID: 39609481 PMCID: PMC11605075 DOI: 10.1038/s41598-024-79043-5] [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: 07/11/2023] [Accepted: 11/05/2024] [Indexed: 11/30/2024] Open
Abstract
Diffusion-weighted MRI is our most promising method for estimating microscopic tissue morphology in vivo. The signal acquisition is based on scanner-generated external magnetic gradients. However, it will also be affected by susceptibility-induced internal magnetic gradients caused by interactions between the tissue and the static magnetic field of the scanner. With 3D in silico experiments, we show how internal gradients cause morphology-, compartment-, and orientation-dependence of spin-echo and pulsed-gradient spin-echo experiments in myelinated axons. These effects surpass those observed with previous 2D modelling corresponding to straight cylinders. For an ex vivo monkey brain, we observe the orientation-dependence generated only when including non-circular cross-sections in the in silico morphological configurations, and find orientation-dependent deviation of up to 17% for diffusion tensor metrics. Interestingly, we find that the orientation-dependence not only biases the signal across different brain regions, but also carries a sensitivity to the morphology of axonal cross-sections which is not attainable by the idealised theoretical diffusion-weighted MRI signal.
Collapse
Affiliation(s)
- S Winther
- Department of Applied Mathematics and Computer Science, Technical University of Denmark, Kongens Lyngby, Denmark.
- Danish Research Centre for Magnetic Resonance, Department of Radiology and Nuclear Medicine, Copenhagen University Hospital - Amager and Hvidovre, 2650, Copenhagen, Denmark.
| | - H Lundell
- Danish Research Centre for Magnetic Resonance, Department of Radiology and Nuclear Medicine, Copenhagen University Hospital - Amager and Hvidovre, 2650, Copenhagen, Denmark
- Department of Health Technology, Technical University of Denmark, Kongens Lyngby, Denmark
| | - J Rafael-Patiño
- Signal Processing Laboratory (LTS5), École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
- Radiology Department, Centre Hospitalier Universitaire Vaudois and University of Lausanne, Lausanne, Switzerland
| | - M Andersson
- Danish Research Centre for Magnetic Resonance, Department of Radiology and Nuclear Medicine, Copenhagen University Hospital - Amager and Hvidovre, 2650, Copenhagen, Denmark
| | - J-P Thiran
- Signal Processing Laboratory (LTS5), École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
- CIBM Center for Biomedical Imaging, Lausanne, Switzerland
- Radiology Department, Centre Hospitalier Universitaire Vaudois and University of Lausanne, Lausanne, Switzerland
| | - T B Dyrby
- Department of Applied Mathematics and Computer Science, Technical University of Denmark, Kongens Lyngby, Denmark.
- Danish Research Centre for Magnetic Resonance, Department of Radiology and Nuclear Medicine, Copenhagen University Hospital - Amager and Hvidovre, 2650, Copenhagen, Denmark.
| |
Collapse
|
6
|
Chen L, Shin HG, van Zijl PC, Li X. Exploiting gradient-echo frequency evolution: Probing white matter microstructure and extracting bulk susceptibility-induced frequency for quantitative susceptibility mapping. Magn Reson Med 2024; 91:1676-1693. [PMID: 38102838 PMCID: PMC10880384 DOI: 10.1002/mrm.29958] [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: 05/14/2023] [Revised: 10/08/2023] [Accepted: 11/17/2023] [Indexed: 12/17/2023]
Abstract
PURPOSE This work is to investigate the microstructure-induced frequency shift in white matter (WM) with crossing fibers and to separate the microstructure-related frequency shift from the bulk susceptibility-induced frequency shift by model fitting the gradient-echo (GRE) frequency evolution for potentially more accurate quantitative susceptibility mapping (QSM). METHODS A hollow-cylinder fiber model (HCFM) with two fiber populations was developed to investigate GRE frequency evolutions in WM voxels with microstructural orientation dispersion. The simulated and experimentally measured TE-dependent local frequency shift was then fitted to a simplified frequency evolution model to obtain a microstructure-related frequency difference parameter (∆ f $$ \Delta f $$ ) and a TE-independent bulk susceptibility-induced frequency shift (C f $$ {C}_f $$ ). The obtainedC f $$ {C}_f $$ was then used for QSM reconstruction. Reconstruction performances were evaluated using a numerical head phantom and in vivo data and then compared to other multi-echo combination methods. RESULTS GRE frequency evolutions and∆ f $$ \Delta f $$ -based tissue parameters in both parallel and crossing fibers determined from our simulations were comparable to those observed in vivo. The TE-dependent frequency fitting method outperformed other multi-echo combination methods in estimatingC f $$ {C}_f $$ in simulations. The fitted∆ f $$ \Delta f $$ ,C f $$ {C}_f $$ , and QSM could be improved further by navigator-based B0 fluctuation correction. CONCLUSION A HCFM with two fiber populations can be used to characterize microstructure-induced frequency shifts in WM regions with crossing fibers. HCFM-based TE-dependent frequency fitting provides tissue contrast related to microstructure (∆ f $$ \Delta f $$ ) and in addition may help improve the quantification accuracy ofC f $$ {C}_f $$ and the corresponding QSM.
Collapse
Affiliation(s)
- Lin Chen
- F.M. Kirby Research Center for Functional Brain Imaging, Kennedy Krieger Institute, Baltimore, Maryland, United States
- Department of Radiology and Radiological Sciences, Johns Hopkins University, Baltimore, Maryland, United States
| | - Hyeong-Geol Shin
- F.M. Kirby Research Center for Functional Brain Imaging, Kennedy Krieger Institute, Baltimore, Maryland, United States
- Department of Radiology and Radiological Sciences, Johns Hopkins University, Baltimore, Maryland, United States
| | - Peter C.M. van Zijl
- F.M. Kirby Research Center for Functional Brain Imaging, Kennedy Krieger Institute, Baltimore, Maryland, United States
- Department of Radiology and Radiological Sciences, Johns Hopkins University, Baltimore, Maryland, United States
| | - Xu Li
- F.M. Kirby Research Center for Functional Brain Imaging, Kennedy Krieger Institute, Baltimore, Maryland, United States
- Department of Radiology and Radiological Sciences, Johns Hopkins University, Baltimore, Maryland, United States
| |
Collapse
|
7
|
Biondetti E, Karsa A, Grussu F, Battiston M, Yiannakas MC, Thomas DL, Shmueli K. Multi-echo quantitative susceptibility mapping: how to combine echoes for accuracy and precision at 3 Tesla. Magn Reson Med 2022; 88:2101-2116. [PMID: 35766450 PMCID: PMC9545116 DOI: 10.1002/mrm.29365] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2022] [Revised: 04/29/2022] [Accepted: 05/27/2022] [Indexed: 12/04/2022]
Abstract
Purpose To compare different multi‐echo combination methods for MRI QSM. Given the current lack of consensus, we aimed to elucidate how to optimally combine multi‐echo gradient‐recalled echo signal phase information, either before or after applying Laplacian‐base methods (LBMs) for phase unwrapping or background field removal. Methods Multi‐echo gradient‐recalled echo data were simulated in a numerical head phantom, and multi‐echo gradient‐recalled echo images were acquired at 3 Tesla in 10 healthy volunteers. To enable image‐based estimation of gradient‐recalled echo signal noise, 5 volunteers were scanned twice in the same session without repositioning. Five QSM processing pipelines were designed: 1 applied nonlinear phase fitting over TEs before LBMs; 2 applied LBMs to the TE‐dependent phase and then combined multiple TEs via either TE‐weighted or SNR‐weighted averaging; and 2 calculated TE‐dependent susceptibility maps via either multi‐step or single‐step QSM and then combined multiple TEs via magnitude‐weighted averaging. Results from different pipelines were compared using visual inspection; summary statistics of susceptibility in deep gray matter, white matter, and venous regions; phase noise maps (error propagation theory); and, in the healthy volunteers, regional fixed bias analysis (Bland–Altman) and regional differences between the means (nonparametric tests). Results Nonlinearly fitting the multi‐echo phase over TEs before applying LBMs provided the highest regional accuracy of χ and the lowest phase noise propagation compared to averaging the LBM‐processed TE‐dependent phase. This result was especially pertinent in high‐susceptibility venous regions. Conclusion For multi‐echo QSM, we recommend combining the signal phase by nonlinear fitting before applying LBMs. Click here for author‐reader discussions
Collapse
Affiliation(s)
- Emma Biondetti
- Institute for Advanced Biomedical Technologies, Department of Neuroscience, Imaging and Clinical Sciences, "D'Annunzio University" of Chieti-Pescara, Chieti, Italy.,Department of Medical Physics and Biomedical Engineering, University College London, London, United Kingdom
| | - Anita Karsa
- Department of Medical Physics and Biomedical Engineering, University College London, London, United Kingdom
| | - Francesco Grussu
- NMR Research Unit, Queen Square MS Centre, Department of Neuroinflammation, UCL Queen Square Institute of Neurology, University College London, London, United Kingdom.,Radiomics Group, Vall d'Hebron Institute of Oncology, Barcelona, Spain
| | - Marco Battiston
- NMR Research Unit, Queen Square MS Centre, Department of Neuroinflammation, UCL Queen Square Institute of Neurology, University College London, London, United Kingdom
| | - Marios C Yiannakas
- NMR Research Unit, Queen Square MS Centre, Department of Neuroinflammation, UCL Queen Square Institute of Neurology, University College London, London, United Kingdom
| | - David L Thomas
- Dementia Research Centre, UCL Queen Square Institute of Neurology, University College London, London, United Kingdom.,Wellcome Centre for Human Neuroimaging, UCL Queen Square Institute of Neurology, University College London, London, United Kingdom
| | - Karin Shmueli
- Department of Medical Physics and Biomedical Engineering, University College London, London, United Kingdom
| |
Collapse
|
8
|
Champagne AA, Coverdale NS, Allen MD, Tremblay JC, MacPherson REK, Pyke KE, Olver TD, Cook DJ. The physiological basis underlying functional connectivity differences in older adults: A multi-modal analysis of resting-state fMRI. Brain Imaging Behav 2022; 16:1575-1591. [PMID: 35092574 DOI: 10.1007/s11682-021-00570-0] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2020] [Accepted: 09/27/2021] [Indexed: 11/02/2022]
Abstract
The purpose of this study was to determine if differences in functional connectivity strength (FCS) with age were confounded by vascular parameters including resting cerebral blood flow (CBF0), cerebrovascular reactivity (CVR), and BOLD-CBF coupling. Neuroimaging data were collected from 13 younger adults (24 ± 2 years) and 14 older adults (71 ± 4 years). A dual-echo resting state pseudo-continuous arterial spin labeling sequence was performed, as well as a BOLD breath-hold protocol. A group independent component analysis was used to identify networks, which were amalgamated into a region of interest (ROI). Within the ROI, FC strength (FCS) was computed for all voxels and compared across the groups. CBF0, CVR and BOLD-CBF coupling were examined within voxels where FCS was different between young and older adults. FCS was greater in old compared to young (P = 0.001). When the effect of CBF0, CVR and BOLD-CBF coupling on FCS was examined, BOLD-CBF coupling had a significant effect (P = 0.003) and group differences in FCS were not present once all vascular parameters were considered in the statistical model (P = 0.07). These findings indicate that future studies of FCS should consider vascular physiological markers in order to improve our understanding of aging processes on brain connectivity.
Collapse
Affiliation(s)
- Allen A Champagne
- Centre for Neuroscience Studies, Queen's University, Room 260, Kingston, ON, K7L 3N6, Canada
| | - Nicole S Coverdale
- Centre for Neuroscience Studies, Queen's University, Room 260, Kingston, ON, K7L 3N6, Canada
| | - Matti D Allen
- Department of Physical Medicine and Rehabilitation, Queen's University, Kingston, ON, Canada.,School of Kinesiology and Health Studies, Cardiovascular Stress Response Laboratory, Queen's University, Kingston, ON, K7L 3N6, Canada.,Department of Physical Medicine and Rehabilitation, Providence Care Hospital, 752 King St., Ontario, West Kingston, Canada
| | - Joshua C Tremblay
- School of Kinesiology and Health Studies, Cardiovascular Stress Response Laboratory, Queen's University, Kingston, ON, K7L 3N6, Canada
| | - Rebecca E K MacPherson
- Department of Health Sciences, Faculty of Applied Health Sciences, Brock University, 1812 Sir Isaac Brock Way, St Catharines, ON, L2S 3A1, Canada
| | - Kyra E Pyke
- School of Kinesiology and Health Studies, Cardiovascular Stress Response Laboratory, Queen's University, Kingston, ON, K7L 3N6, Canada
| | - T Dylan Olver
- Biomedical Sciences, Western College of Veterinarian Medicine, University of Saskatchewan, 52 Campus Drive, Saskatoon, Saskatchewan, S7N 5B4, Canada
| | - Douglas J Cook
- Centre for Neuroscience Studies, Queen's University, Room 260, Kingston, ON, K7L 3N6, Canada. .,Department of Surgery, Queen's University, Room 232, 18 Stuart St, Kingston, ON, K7L 3N6, Canada.
| |
Collapse
|
9
|
Schäper J, Bauman G, Ganter C, Bieri O. Pure balanced steady-state free precession imaging (pure bSSFP). Magn Reson Med 2021; 87:1886-1893. [PMID: 34775622 PMCID: PMC9299476 DOI: 10.1002/mrm.29086] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2021] [Revised: 10/02/2021] [Accepted: 10/28/2021] [Indexed: 12/26/2022]
Abstract
Purpose To show that for tissues the conspicuous asymmetries in the frequency response function of bSSFP can be mitigated by using a short enough TR. Theory and Methods Configuration theory indicates that bSSFP becomes apparently “pure” (i.e., exhibiting a symmetric profile) in the limit of TR →0. To this end, the frequency profile of bSSFP was measured as a function of the TR using a manganese‐doped aqueous probe, as well as brain tissue that was shown to exhibit a pronounced asymmetry due to its microstructure. The frequency response function was sampled using N=72 (phantom) and N=36 (in vivo) equally distributed linear RF phase increments in the interval [0,2π). Imaging was performed with 2.0 mm isotropic resolution over a TR range of 1.5–8 ms at 3 and 1.5 T. Results As expected, pure substances showed a symmetric TR‐independent frequency profile, whereas brain tissue revealed a pronounced asymmetry. The observed asymmetry for the tissue, however, decreases with decreasing TR and gives strong evidence that the frequency response function of bSSFP becomes symmetric in the limit of TR →0, in agreement with theory. The limit of apparently pure bSSFP imaging can thus be achieved for a TR ∼ 1.5 ms at 1.5 T, whereas at 3 T, tissues still show some residual asymmetry. Conclusion In the limit of short enough TR, tissues become apparently pure for bSSFP. This limit can be reached for brain tissue at 1.5 T with TR ∼ 1–2 ms at clinically relevant resolutions.
Collapse
Affiliation(s)
- Jessica Schäper
- Department of Biomedical Engineering, University of Basel, Basel, Switzerland.,Division of Radiological Physics, Department of Radiology, University Hospital Basel, Basel, Switzerland
| | - Grzegorz Bauman
- Department of Biomedical Engineering, University of Basel, Basel, Switzerland.,Division of Radiological Physics, Department of Radiology, University Hospital Basel, Basel, Switzerland
| | - Carl Ganter
- Department of Diagnostic Radiology, Klinikum rechts der Isar, Technical University of Munich, Munich, Germany
| | - Oliver Bieri
- Department of Biomedical Engineering, University of Basel, Basel, Switzerland.,Division of Radiological Physics, Department of Radiology, University Hospital Basel, Basel, Switzerland
| |
Collapse
|
10
|
Cottaar M, Wu W, Tendler BC, Nagy Z, Miller K, Jbabdi S. Quantifying myelin in crossing fibers using diffusion-prepared phase imaging: Theory and simulations. Magn Reson Med 2021; 86:2618-2634. [PMID: 34254349 PMCID: PMC8581995 DOI: 10.1002/mrm.28907] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2020] [Revised: 06/03/2021] [Accepted: 06/08/2021] [Indexed: 11/24/2022]
Abstract
PURPOSE Myelin has long been the target of neuroimaging research. However, most available techniques can only provide a voxel-averaged estimate of myelin content. In the human brain, white matter fiber pathways connecting different brain areas and carrying different functions often cross each other in the same voxel. A measure that can differentiate the degree of myelination of crossing fibers would provide a more specific marker of myelination. THEORY AND METHODS One MRI signal property that is sensitive to myelin is the phase accumulation. This sensitivity is used by measuring the phase accumulation of the signal remaining after diffusion-weighting, which is called diffusion-prepared phase imaging (DIPPI). Including diffusion-weighting before estimating the phase accumulation has two distinct advantages for estimating the degree of myelination: (1) It increases the relative contribution of intra-axonal water, whose phase is related linearly to the thickness of the surrounding myelin (in particular the log g-ratio); and (2) it gives directional information, which can be used to distinguish between crossing fibers. Here the DIPPI sequence is described, an approach is proposed to estimate the log g-ratio, and simulations are used and DIPPI data acquired in an isotropic phantom to quantify other sources of phase accumulation. RESULTS The expected bias is estimated in the log g-ratio for reasonable in vivo acquisition parameters caused by eddy currents (~4%-10%), remaining extra-axonal signal (~15%), and gradients in the bulk off-resonance field (<10% for most of the brain). CONCLUSION This new sequence may provide a g-ratio estimate per fiber population crossing within a voxel.
Collapse
Affiliation(s)
- Michiel Cottaar
- Wellcome Centre for Integrative Neuroimaging—Centre for Functional Magnetic Resonance Imaging of the BrainJohn Radcliffe HospitalUniversity of OxfordOxfordUnited Kingdom
| | - Wenchuan Wu
- Wellcome Centre for Integrative Neuroimaging—Centre for Functional Magnetic Resonance Imaging of the BrainJohn Radcliffe HospitalUniversity of OxfordOxfordUnited Kingdom
| | - Benjamin C. Tendler
- Wellcome Centre for Integrative Neuroimaging—Centre for Functional Magnetic Resonance Imaging of the BrainJohn Radcliffe HospitalUniversity of OxfordOxfordUnited Kingdom
| | - Zoltan Nagy
- Laboratory for Social and Neural Systems ResearchUniversity of ZurichZurichSwitzerland
| | - Karla Miller
- Wellcome Centre for Integrative Neuroimaging—Centre for Functional Magnetic Resonance Imaging of the BrainJohn Radcliffe HospitalUniversity of OxfordOxfordUnited Kingdom
| | - Saad Jbabdi
- Wellcome Centre for Integrative Neuroimaging—Centre for Functional Magnetic Resonance Imaging of the BrainJohn Radcliffe HospitalUniversity of OxfordOxfordUnited Kingdom
| |
Collapse
|
11
|
Yeo SJ, Lee SH, Lee SK. Rapid calculation of static magnetic field perturbation generated by magnetized objects in arbitrary orientations. Magn Reson Med 2021; 87:1015-1027. [PMID: 34617634 DOI: 10.1002/mrm.29037] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/09/2021] [Revised: 08/17/2021] [Accepted: 09/18/2021] [Indexed: 11/08/2022]
Abstract
PURPOSE Most previous work on the calculation of susceptibility-induced static magnetic field (B0 ) inhomogeneity has considered strictly unidirectional magnetic fields. Here, we present the theory and implementation of a computational method to rapidly calculate static magnetic field vectors produced by an arbitrary distribution of voxelated magnetization vectors. THEORY AND METHODS Two existing B0 calculation methods were systematically extended to include arbitrary orientations of the magnetization and the magnetic field; they are (1) Fourier-domain convolution with k-space-discretized (KD) dipolar field, and (2) generalized susceptibility voxel convolution (gSVC). The methods were tested on an analytical ellipsoid model and a tilted human head model, as well as against experimentally measured B0 fields induced by a stainless-steel implant located in an inhomogeneous region of a clinical 3T MRI magnet. RESULTS Both methods were capable of correctly calculating B0 fields inside a magnetized ellipsoid in all tested orientations. The KD method generally required a larger grid and longer computation time to achieve accuracy comparable to gSVC. Measured B0 fields due to the implant showed a good match with the gSVC-calculated fields that accounted for the spatial variation of the applied magnetic field including the radial components. CONCLUSION Our method can provide a reliable and efficient computational tool to calculate B0 perturbation by magnetized objects under a variety of circumstances, including those with inhomogeneous magnetizing fields, anisotropic susceptibility, and a rotated coordinate system.
Collapse
Affiliation(s)
- Seok-Jin Yeo
- Department of Biomedical Engineering, Sungkyunkwan University, Suwon, South Korea
| | - So-Hee Lee
- Department of Biomedical Engineering, Sungkyunkwan University, Suwon, South Korea.,Department of Intelligent Precision Healthcare Convergence, Sungkyunkwan University, Suwon, South Korea
| | - Seung-Kyun Lee
- Department of Biomedical Engineering, Sungkyunkwan University, Suwon, South Korea.,Department of Intelligent Precision Healthcare Convergence, Sungkyunkwan University, Suwon, South Korea.,Center for Neuroscience Imaging Research, Institute for Basic Science (IBS), Suwon, South Korea.,Department of Physics, Sungkyunkwan University, Suwon, South Korea
| |
Collapse
|
12
|
In vivo assessment of anisotropy of apparent magnetic susceptibility in white matter from a single orientation acquisition. Neuroimage 2021; 241:118442. [PMID: 34339831 DOI: 10.1016/j.neuroimage.2021.118442] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2020] [Revised: 06/02/2021] [Accepted: 07/29/2021] [Indexed: 11/20/2022] Open
Abstract
Multiple studies have reported a significant dependence of the effective transverse relaxation rate constant (R2*) and the phase of gradient-echo based (GRE) signal on the orientation of white matter fibres in the human brain. It has also been hypothesized that magnetic susceptibility, as obtained by single-orientation quantitative susceptibility mapping (QSM), exhibits such a dependence. In this study, we investigated this hypothesized relationship in a cohort of healthy volunteers. We show that R2* follows the predicted orientation dependence consistently across white matter regions, whereas the apparent magnetic susceptibility is related differently to fibre orientation across the brain and often in a complex non-monotonic manner. In addition, we explored the effect of fractional anisotropy measured by diffusion-weighted MRI on the strength of the orientation dependence and observed only a limited influence in many regions. However, with careful consideration of such an impact and the limitations imposed by the ill-posed nature of the dipole inversion process, it is possible to study magnetic susceptibility anisotropy in specific brain regions with a single orientation acquisition.
Collapse
|
13
|
Kaden E, Gyori NG, Rudrapatna SU, Barskaya IY, Dragonu I, Does MD, Jones DK, Clark CA, Alexander DC. Microscopic susceptibility anisotropy imaging. Magn Reson Med 2020; 84:2739-2753. [PMID: 32378746 PMCID: PMC7402021 DOI: 10.1002/mrm.28303] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2020] [Revised: 04/07/2020] [Accepted: 04/08/2020] [Indexed: 11/25/2022]
Abstract
PURPOSE The gradient-echo MR signal in brain white matter depends on the orientation of the fibers with respect to the external magnetic field. To map microstructure-specific magnetic susceptibility in orientationally heterogeneous material, it is thus imperative to regress out unwanted orientation effects. METHODS This work introduces a novel framework, referred to as microscopic susceptibility anisotropy imaging, that disentangles the 2 principal effects conflated in gradient-echo measurements, (a) the susceptibility properties of tissue microenvironments, especially the myelin microstructure, and (b) the axon orientation distribution relative to the magnetic field. Specifically, we utilize information about the orientational tissue structure inferred from diffusion MRI data to factor out the B 0 -direction dependence of the frequency difference signal. RESULTS A human pilot study at 3 T demonstrates proxy maps of microscopic susceptibility anisotropy unconfounded by fiber crossings and orientation dispersion as well as magnetic field direction. The developed technique requires only a dual-echo gradient-echo scan acquired at 1 or 2 head orientations with respect to the magnetic field and a 2-shell diffusion protocol achievable on standard scanners within practical scan times. CONCLUSIONS The quantitative recovery of microscopic susceptibility features in the presence of orientational heterogeneity potentially improves the assessment of microstructural tissue integrity.
Collapse
Affiliation(s)
- Enrico Kaden
- Centre for Medical Image ComputingUniversity College LondonLondonUK
- Great Ormond Street Institute of Child HealthUniversity College LondonLondonUK
| | - Noemi G. Gyori
- Centre for Medical Image ComputingUniversity College LondonLondonUK
- Great Ormond Street Institute of Child HealthUniversity College LondonLondonUK
| | | | | | | | - Mark D. Does
- Institute of Imaging ScienceVanderbilt UniversityNashvilleTNUSA
| | - Derek K. Jones
- Cardiff University Brain Research Imaging CentreCardiff UniversityCardiffUK
- School of PsychologyAustralian Catholic UniversityMelbourneVICAustralia
| | - Chris A. Clark
- Great Ormond Street Institute of Child HealthUniversity College LondonLondonUK
| | | |
Collapse
|
14
|
Foxley S, Wildenberg G, Sampathkumar V, Karczmar GS, Brugarolas P, Kasthuri N. Sensitivity to myelin using model-free analysis of the water resonance line-shape in postmortem mouse brain. Magn Reson Med 2020; 85:667-677. [PMID: 32783262 DOI: 10.1002/mrm.28440] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2020] [Revised: 06/17/2020] [Accepted: 06/30/2020] [Indexed: 12/18/2022]
Abstract
PURPOSE Dysmyelinating diseases are characterized by abnormal myelin formation and function. Such microstructural abnormalities in myelin have been demonstrated to produce measurable effects on the MR signal. This work examines these effects on measurements of voxel-wise, high-resolution water spectra acquired using a 3D echo-planar spectroscopic imaging (EPSI) pulse sequence from both postmortem fixed control mouse brains and a dysmyelination mouse brain model. METHODS Perfusion fixed, resected control (n = 5) and shiverer (n = 4) mouse brains were imaged using 3D-EPSI with 100 µm isotropic resolution. The free induction decay (FID) was sampled every 2.74 ms over 192 echoes, for a total sampling duration of 526.08 ms. Voxel-wise FIDs were Fourier transformed to produce water spectra with 1.9 Hz resolution. Spectral asymmetry was computed and compared between the two tissue types. RESULTS The water resonance is more asymmetrically broadened in the white matter of control mouse brain compared with dysmyelinated white matter. In control brain, this is modulated by and consistent with previously reported orientationally dependent effects of white matter relative to B0 . Similar sensitivity to orientation is observed in dysmyelinated white matter as well; however, the magnitude of the resonance asymmetry is much lower across all directions. CONCLUSION Results demonstrate that components of the spectra are specifically differentially affected by myelin concentration. This suggests that water proton spectra may be sensitive to the presence of myelin, and as such, could serve as a MRI-based biomarker of dysmyelinating disease, free of mathematical models.
Collapse
Affiliation(s)
- Sean Foxley
- Department of Radiology, University of Chicago, Chicago, Illinois, USA
| | - Gregg Wildenberg
- Department of Neurobiology, University of Chicago, Chicago, Illinois, USA
| | | | | | - Pedro Brugarolas
- Department of Radiology, Harvard Medical School, Boston, Maryland, USA.,Gordon Center for Medical Imaging, Massachusetts General Hospital, Boston, Maryland, USA
| | - Narayanan Kasthuri
- Department of Neurobiology, University of Chicago, Chicago, Illinois, USA
| |
Collapse
|
15
|
Biondetti E, Karsa A, Thomas DL, Shmueli K. Investigating the accuracy and precision of TE-dependent versus multi-echo QSM using Laplacian-based methods at 3 T. Magn Reson Med 2020; 84:3040-3053. [PMID: 32491224 DOI: 10.1002/mrm.28331] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2019] [Revised: 03/19/2020] [Accepted: 04/29/2020] [Indexed: 12/14/2022]
Abstract
PURPOSE Multi-echo gradient-recalled echo acquisitions for QSM enable optimizing the SNR for several tissue types through multi-echo (TE) combination or investigating temporal variations in the susceptibility (potentially reflecting tissue microstructure) by calculating one QSM image at each TE (TE-dependent QSM). In contrast with multi-echo QSM, applying Laplacian-based methods (LBMs) for phase unwrapping and background field removal to single TEs could introduce nonlinear temporal variations (independent of tissue microstructure) into the measured susceptibility. Here, we aimed to compare the effect of LBMs on the QSM susceptibilities in TE-dependent versus multi-echo QSM. METHODS TE-dependent recalled echo data simulated in a numerical head phantom and gradient-recalled echo images acquired at 3 T in 10 healthy volunteers. Several QSM pipelines were tested, including four distinct LBMs: sophisticated harmonic artifact reduction for phase data (SHARP), variable-radius sophisticated harmonic artifact reduction for phase data (V-SHARP), Laplacian boundary value background field removal (LBV), and one-step total generalized variation (TGV). Results from distinct pipelines were compared using visual inspection, summary statistics of susceptibility in deep gray matter/white matter/venous regions of interest, and, in the healthy volunteers, regional susceptibility bias analysis and nonparametric tests. RESULTS Multi-echo versus TE-dependent QSM had higher regional accuracy, especially in high-susceptibility regions and at shorter TEs. Everywhere except in the veins, a processing pipeline incorporating TGV provided the most temporally stable TE-dependent QSM results with an accuracy similar to multi-echo QSM. CONCLUSIONS For TE-dependent QSM, carefully choosing LBMs can minimize the introduction of LBM-related nonlinear temporal susceptibility variations.
Collapse
Affiliation(s)
- Emma Biondetti
- Centre de NeuroImagerie de Recherche (CENIR), Team "Movement Investigations and Therapeutics", Institut du Cerveau (ICM), Inserm U 1127, CNRS UMR 7225, Sorbonne Université, Paris, France.,Department of Medical Physics and Biomedical Engineering, University College London, London, United Kingdom
| | - Anita Karsa
- Department of Medical Physics and Biomedical Engineering, University College London, London, United Kingdom
| | - David L Thomas
- Dementia Research Centre, UCL Queen Square Institute of Neurology, University College London, London, United Kingdom.,Wellcome Centre for Human Neuroimaging, UCL Queen Square Institute of Neurology, University College London, London, United Kingdom
| | - Karin Shmueli
- Department of Medical Physics and Biomedical Engineering, University College London, London, United Kingdom
| |
Collapse
|
16
|
Sood S, Reutens DC, Kadamangudi S, Barth M, Vegh V. Field strength influences on gradient recalled echo MRI signal compartment frequency shifts. Magn Reson Imaging 2020; 70:98-107. [PMID: 32360972 DOI: 10.1016/j.mri.2020.04.018] [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: 06/25/2019] [Revised: 04/09/2020] [Accepted: 04/27/2020] [Indexed: 11/20/2022]
Abstract
Different trends of echo time dependent gradient recalled echo MRI signals in different brain regions have been attributed to signal compartments in image voxels. It remains unclear how variations in gradient recalled echo MRI signals change as a function of MRI field strength, and how data processing may impact signal compartment parameters. We used two popular quantitative susceptibility mapping methods of processing raw phase images (Laplacian and path-based unwrapping with V-SHARP) and expressed values in the form of induced frequency shifts (in Hz) in six specific brain regions at 3T and 7T. We found the frequency shift curves to vary with echo time, and a good overlap between 3T and 7T mean frequency shift curves was present. However, the amount of variation across participants was greater at 3T, and we were able to obtain better compartment model fits of the signal at 7T. We also found the temporal trends in the signal and compartment frequency shifts to change with the method used to process images. The inter-participant averaged trends were consistent between 3T and 7T for each quantitative susceptibility pipeline. However, signal compartment frequency shifts generated using different pipelines may not be comparable.
Collapse
Affiliation(s)
- Surabhi Sood
- Centre for Advanced Imaging, University of Queensland, Brisbane, Queensland, Australia
| | - David C Reutens
- Centre for Advanced Imaging, University of Queensland, Brisbane, Queensland, Australia
| | - Shrinath Kadamangudi
- Centre for Advanced Imaging, University of Queensland, Brisbane, Queensland, Australia
| | - Markus Barth
- Centre for Advanced Imaging, University of Queensland, Brisbane, Queensland, Australia
| | - Viktor Vegh
- Centre for Advanced Imaging, University of Queensland, Brisbane, Queensland, Australia.
| |
Collapse
|
17
|
Thapaliya K, Vegh V, Bollmann S, Barth M. Influence of 7T GRE-MRI Signal Compartment Model Choice on Tissue Parameters. Front Neurosci 2020; 14:271. [PMID: 32457565 PMCID: PMC7206227 DOI: 10.3389/fnins.2020.00271] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2019] [Accepted: 03/10/2020] [Indexed: 12/16/2022] Open
Abstract
Quantitative assessment of tissue microstructure is important in studying human brain diseases and disorders. Ultra-high field magnetic resonance imaging (MRI) data obtained using a multi-echo gradient echo sequence have been shown to contain information on myelin, axonal, and extracellular compartments in tissue. Quantitative assessment of water fraction, relaxation time (T2*), and frequency shift using multi-compartment models has been shown to be useful in studying white matter properties via specific tissue parameters. It remains unclear how tissue parameters vary with model selection based on 7T multiple echo time gradient-recalled echo (GRE) MRI data. We applied existing signal compartment models to the corpus callosum and investigated whether a three-compartment model can be reduced to two compartments and still resolve white matter parameters [i.e., myelin water fraction (MWF) and g-ratio]. We show that MWF should be computed using a three-compartment model in the corpus callosum, and the g-ratios obtained using three compartment models are consistent with previous reports. We provide results for other parameters, such as signal compartment frequency shifts.
Collapse
Affiliation(s)
- Kiran Thapaliya
- Centre for Advanced Imaging, The University of Queensland, Brisbane, QLD, Australia
| | - Viktor Vegh
- Centre for Advanced Imaging, The University of Queensland, Brisbane, QLD, Australia.,ARC Centre for Innovation in Biomedical Imaging Technology, Brisbane, QLD, Australia
| | - Steffen Bollmann
- Centre for Advanced Imaging, The University of Queensland, Brisbane, QLD, Australia.,ARC Centre for Innovation in Biomedical Imaging Technology, Brisbane, QLD, Australia
| | - Markus Barth
- Centre for Advanced Imaging, The University of Queensland, Brisbane, QLD, Australia.,ARC Centre for Innovation in Biomedical Imaging Technology, Brisbane, QLD, Australia.,School of Information Technology and Electrical Engineering, The University of Queensland, Brisbane, QLD, Australia
| |
Collapse
|
18
|
Kaczmarz S, Göttler J, Zimmer C, Hyder F, Preibisch C. Characterizing white matter fiber orientation effects on multi-parametric quantitative BOLD assessment of oxygen extraction fraction. J Cereb Blood Flow Metab 2020; 40:760-774. [PMID: 30952200 PMCID: PMC7168796 DOI: 10.1177/0271678x19839502] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 05/07/2018] [Revised: 01/23/2019] [Accepted: 02/22/2019] [Indexed: 12/19/2022]
Abstract
Relative oxygen extraction fraction (rOEF) is a fundamental indicator of cerebral metabolic function. An easily applicable method for magnetic resonance imaging (MRI) based rOEF mapping is the multi-parametric quantitative blood oxygenation level dependent (mq-BOLD) approach with separate acquisitions of transverse relaxation times T 2 * and T2 and dynamic susceptibility contrast (DSC) based relative cerebral blood volume (rCBV). Given that transverse relaxation and rCBV in white matter (WM) strongly depend on nerve fiber orientation, mq-BOLD derived rOEF is expected to be affected as well. To investigate fiber orientation related rOEF artefacts, we present a methodological study characterizing anisotropy effects of WM as measured by diffusion tensor imaging (DTI) on mq-BOLD in 30 healthy volunteers. Using a 3T clinical MRI-scanner, we performed a comprehensive correlation of all parameters ( T 2 * , T2, R 2 ' , rCBV, rOEF, where R 2 ' =1/ T 2 * -1/T2) with DTI-derived fiber orientation towards the main magnetic field (B0). Our results confirm strong dependencies of transverse relaxation and rCBV on the nerve fiber orientation towards B0, with anisotropy-driven variations up to 37%. Comparably weak orientation-dependent variations of mq-BOLD derived rOEF (3.8%) demonstrate partially counteracting influences of R 2 ' and rCBV effects, possibly suggesting applicability of rOEF as an oxygenation sensitive biomarker. However, unresolved issues warrant caution when applying mq-BOLD to WM.
Collapse
Affiliation(s)
- Stephan Kaczmarz
- Department of Diagnostic and Interventional Neuroradiology, Klinikum rechts der Isar, Technical University of Munich, Munich, Germany
- TUM Neuroimaging Center (TUM-NIC), Klinikum rechts der Isar, Technical University of Munich, Munich, Germany
- Departments of Radiology & Biomedical Imaging and of Biomedical Engineering, Magnetic Resonance Research Center, Yale University, New Haven, CT, USA
| | - Jens Göttler
- Department of Diagnostic and Interventional Neuroradiology, Klinikum rechts der Isar, Technical University of Munich, Munich, Germany
- TUM Neuroimaging Center (TUM-NIC), Klinikum rechts der Isar, Technical University of Munich, Munich, Germany
- Departments of Radiology & Biomedical Imaging and of Biomedical Engineering, Magnetic Resonance Research Center, Yale University, New Haven, CT, USA
- Department of Diagnostic and Interventional Radiology, Klinikum rechts der Isar, Technical University of Munich, Munich, Germany
| | - Claus Zimmer
- Department of Diagnostic and Interventional Neuroradiology, Klinikum rechts der Isar, Technical University of Munich, Munich, Germany
| | - Fahmeed Hyder
- Departments of Radiology & Biomedical Imaging and of Biomedical Engineering, Magnetic Resonance Research Center, Yale University, New Haven, CT, USA
| | - Christine Preibisch
- Department of Diagnostic and Interventional Neuroradiology, Klinikum rechts der Isar, Technical University of Munich, Munich, Germany
- TUM Neuroimaging Center (TUM-NIC), Klinikum rechts der Isar, Technical University of Munich, Munich, Germany
- Clinic for Neurology, Klinikum rechts der Isar, Technical University of Munich, Munich, Germany
| |
Collapse
|
19
|
Weber AM, Zhang Y, Kames C, Rauscher A. Myelin water imaging and R 2* mapping in neonates: Investigating R 2* dependence on myelin and fibre orientation in whole brain white matter. NMR IN BIOMEDICINE 2020; 33:e4222. [PMID: 31846134 DOI: 10.1002/nbm.4222] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/08/2019] [Revised: 09/27/2019] [Accepted: 10/20/2019] [Indexed: 06/10/2023]
Abstract
R2* relaxation provides a semiquantitative method of detecting myelin, iron and white matter fibre orientation angles. Compared with standard histogram-based analyses, angle-resolved analysis of R2* has previously been shown to substantially improve the detection of subtle differences in the brain between healthy siblings of subjects with multiple sclerosis and unrelated healthy controls. Neonates, who are born with very little myelin and iron, and an underdeveloped connectome, provide researchers with an opportunity to investigate whether R2* is intimately linked with fibre-angle or myelin content as it is in adults, which may in future studies be explored as a potential white matter developmental biomarker. Five healthy adult volunteers (mean age [±SD] = 31.2 [±8.3] years; three males) were recruited from Vancouver, Canada. Eight term neonates (mean age = 38.6 ± 1.2 weeks; five males) were recruited from the Children's Hospital of Chongqing Medical University neonatal ward. All subjects were scanned on identical 3 T Philips Achieva scanners equipped with an eight-channel SENSE head coil and underwent a multiecho gradient echo scan, a 32-direction DTI scan and a myelin water imaging scan. For both neonates and adults, bin-averaged R2* variation across the brain's white matter was found to be best explained by fibre orientation. For adults, this represented a difference in R2* values of 3.5 Hz from parallel to perpendicular fibres with respect to the main magnetic field. In neonates, the fibre orientation dependency displayed a cosine wave shape, with a small R2* range of 0.4 Hz. This minor relationship in neonates provides further evidence for the key role myelin probably plays in creating this fibre orientation dependence later in life, but suggests limited clinical application in newborn populations. Future studies should investigate fibre-orientation dependency in infants in the first 5 years, when substantial myelin development occurs.
Collapse
Affiliation(s)
- Alexander Mark Weber
- Division of Neurology, Department of Pediatrics, University of British Columbia, Vancouver, BC, Canada
- UBC MRI Research Centre, University of British Columbia, Vancouver, BC, Canada
| | - Yuting Zhang
- Department of Radiology, Children's Hospital of Chongqing Medical University, Chongqing, China
- Ministry of Education Key Laboratory of Child Development and Disorders, Chongqing Medical University, Chongqing, China
- Key Laboratory of Pediatrics in Chongqing, Chongqing Medical University, Chongqing, China
- Chongqing International Science and Technology Cooperation Center for Child Development and Disorders, Chongqing, China
| | - Christian Kames
- UBC MRI Research Centre, University of British Columbia, Vancouver, BC, Canada
- Department of Physics & Astronomy, University of British Columbia, Vancouver, BC, Canada
| | - Alexander Rauscher
- Division of Neurology, Department of Pediatrics, University of British Columbia, Vancouver, BC, Canada
- UBC MRI Research Centre, University of British Columbia, Vancouver, BC, Canada
- Department of Physics & Astronomy, University of British Columbia, Vancouver, BC, Canada
- Department of Radiology, University of British Columbia, Vancouver, BC, Canada
| |
Collapse
|
20
|
Ruh A, Kiselev VG. Larmor frequency dependence on structural anisotropy of magnetically heterogeneous media. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2019; 307:106584. [PMID: 31476632 DOI: 10.1016/j.jmr.2019.106584] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/08/2019] [Revised: 08/21/2019] [Accepted: 08/22/2019] [Indexed: 06/10/2023]
Abstract
The effect of anisotropic magnetic microstructure on the measurable Larmor frequency offset is investigated in media with heterogeneous magnetic susceptibility using Monte Carlo simulations. The focus is on the transition between the regimes of fast and slow diffusion of NMR-reporting molecules. Simulations demonstrate a perfect agreement with the previously developed analytic theory for fast diffusion. Beyond this regime, the frequency offset shows a pronounced dependence on the medium microarchitecture and the diffusivity of NMR-reporting spins in relation to the magnitude of the susceptibility-induced magnetic field.
Collapse
Affiliation(s)
- Alexander Ruh
- Medical Physics, Dept. of Radiology, Medical Center - University of Freiburg, Faculty of Medicine, University of Freiburg, Germany; Dept. of Radiology, Feinberg School of Medicine, Northwestern University, Chicago, IL, United States
| | - Valerij G Kiselev
- Medical Physics, Dept. of Radiology, Medical Center - University of Freiburg, Faculty of Medicine, University of Freiburg, Germany.
| |
Collapse
|
21
|
Abstract
The study of brain plasticity has tended to focus on the synapse, where well-described activity-dependent mechanisms are known to play a key role in learning and memory. However, it is becoming increasingly clear that plasticity occurs beyond the synapse. This review focuses on the emerging concept of white matter plasticity. For example, there is growing evidence, both from animal studies and from human neuroimaging, that activity-dependent regulation of myelin may play a role in learning. This previously overlooked phenomenon may provide a complementary but powerful route through which experience shapes the brain.
Collapse
|
22
|
Santiago I, Santinha J, Ianus A, Galzerano A, Theias R, Maia J, Barata MJ, Loução N, Costa-Silva B, Beltran A, Matos C, Shemesh N. Susceptibility Perturbation MRI Maps Tumor Infiltration into Mesorectal Lymph Nodes. Cancer Res 2019; 79:2435-2444. [PMID: 30894376 DOI: 10.1158/0008-5472.can-18-3682] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/27/2018] [Revised: 02/13/2019] [Accepted: 03/15/2019] [Indexed: 11/16/2022]
Abstract
Noninvasive characterization of lymph node involvement in cancer is an enduring onerous challenge. In rectal cancer, pathologic lymph node status constitutes the most important determinant of local recurrence and overall survival, and patients with involved lymph nodes may benefit from preoperative chemo and/or radiotherapy. However, knowledge of lymph node status before surgery is currently hampered by limited imaging accuracy. Here, we introduce Susceptibility-Perturbation MRI (SPI) as a novel source of contrast to map malignant infiltration into mesorectal lymph nodes. SPI involves multigradient echo (MGE) signal decays presenting a nonmonoexponential nature, which we show is sensitive to the underlying microstructure via susceptibility perturbations. Using numerical simulations, we predicted that the large cell morphology and the high cellularity of tumor within affected mesorectal lymph nodes would induce signature SPI decays. We validated this prediction in mesorectal lymph nodes excised from total mesorectal excision specimens of patients with rectal cancer using ultrahigh field (16.4 T) MRI. SPI signals distinguished benign from malignant nodal tissue, both qualitatively and quantitatively, and our histologic analyses confirmed cellularity and cell size were the likely underlying sources for the differences observed. SPI was then adapted to a clinical 1.5 T scanner, added to patients' staging protocol, and compared with conventional assessment by two expert radiologists. Nonmonoexponential decays, similar to those observed in the ex vivo study, were demonstrated, and SPI classified lymph nodes more accurately than standard high-resolution T2-weighted imaging assessment. These findings suggest this simple, yet highly informative, method can improve rectal cancer patient selection for neoadjuvant therapy. SIGNIFICANCE: These findings introduce an MRI methodology tailored to detect magnetic susceptibility perturbations induced by subtle alterations in tissue microstructure.
Collapse
Affiliation(s)
- Inês Santiago
- Neuroplasticity and Neural Activity Lab, Champalimaud Research, Champalimaud Centre for the Unknown, Lisbon, Portugal.,Radiology Department, Champalimaud Centre for the Unknown, Lisbon, Portugal.,Nova Medical School, Campo Mártires da Pátria, Lisbon, Portugal
| | - João Santinha
- Computational Clinical Imaging Group, Champalimaud Centre for the Unknown, Lisbon, Portugal
| | - Andrada Ianus
- Neuroplasticity and Neural Activity Lab, Champalimaud Research, Champalimaud Centre for the Unknown, Lisbon, Portugal.,Centre for Medical Imaging Computing, Department of Computer Science, University College London, London, United Kingdom
| | - Antonio Galzerano
- Pathology Department, Champalimaud Centre for the Unknown, Lisbon, Portugal
| | - Rita Theias
- Pathology Department, Hospital Fernando Fonseca, Amadora, Portugal
| | - Joana Maia
- Systems Oncology Lab, Champalimaud Research, Champalimaud Centre for the Unknown, Lisbon, Portugal
| | - Maria J Barata
- Radiology Department, Champalimaud Centre for the Unknown, Lisbon, Portugal
| | | | - Bruno Costa-Silva
- Systems Oncology Lab, Champalimaud Research, Champalimaud Centre for the Unknown, Lisbon, Portugal
| | - Antonio Beltran
- Pathology Department, Champalimaud Centre for the Unknown, Lisbon, Portugal
| | - Celso Matos
- Radiology Department, Champalimaud Centre for the Unknown, Lisbon, Portugal
| | - Noam Shemesh
- Neuroplasticity and Neural Activity Lab, Champalimaud Research, Champalimaud Centre for the Unknown, Lisbon, Portugal.
| |
Collapse
|
23
|
Kiselev VG. Larmor frequency in heterogeneous media. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2019; 299:168-175. [PMID: 30639748 DOI: 10.1016/j.jmr.2018.12.008] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/14/2018] [Revised: 11/28/2018] [Accepted: 12/13/2018] [Indexed: 06/09/2023]
Abstract
High-precision signal phase measurements ignited an ongoing discussion of the microstructural correlates of the Larmor frequency shift of water in biological tissues. In a broader context, this is the question about the averaged precession frequency in magnetically heterogenous, in particular, porous media. In this study, the Larmor frequency shift is found analytically for water filling connected pore space between NMR-invisible magnetized inclusions with a constant magnetic susceptibility tensor. The magnetic microstructure that encompasses the inclusions' shape and their spatial arrangement is arbitrary as well as the inclusions' magnetic susceptibility tensor. The result is limited to the case of effectively fast diffusion of water molecules. In this limit, the effect of magnetic microstructure on the Larmor frequency shift is represented by only five relevant parameters in the general case and by a single parameter in the case of axially symmetric microstructure. This single parameter enters the dependence of the Larmor frequency on the sample orientation relative to the main magnetic field in the previously performed experiments. The result can help interpreting known experimental data and developing realistic models of biological tissues.
Collapse
Affiliation(s)
- Valerij G Kiselev
- Dept. of Radiology, Medical Physics, Medical Center, Faculty of Medicine, University of Freiburg, Freiburg, Germany.
| |
Collapse
|
24
|
Kurz FT, Buschle LR, Hahn A, Jende JME, Bendszus M, Heiland S, Ziener CH. Diffusion effects in myelin sheath free induction decay. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2018; 297:61-75. [PMID: 30366221 DOI: 10.1016/j.jmr.2018.10.001] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/14/2018] [Revised: 09/27/2018] [Accepted: 10/01/2018] [Indexed: 06/08/2023]
Abstract
Myelin sheath microstructure and composition produce MR signal decay characteristics that can be used to evaluate status and outcome of demyelinating disease. We extend a recently proposed model of neuronal magnetic susceptibility, that accounts for both the structural and inherent anisotropy of the myelin sheath, by including the whole dynamic range of diffusion effects. The respective Bloch-Torrey equation for local spin dephasing is solved with a uniformly convergent perturbation expansion method, and the resulting magnetization decay is validated with a numerical solution based on a finite difference method. We show that a variation of diffusion strengths can lead to substantially different MR signal decay curves. Our results may be used to adjust or control simulations for water diffusion in neuronal structures.
Collapse
Affiliation(s)
- F T Kurz
- Heidelberg University Hospital, INF 400, D-69120 Heidelberg, Germany; German Cancer Research Center, INF 280, D-69120 Heidelberg, Germany.
| | - L R Buschle
- Heidelberg University Hospital, INF 400, D-69120 Heidelberg, Germany; German Cancer Research Center, INF 280, D-69120 Heidelberg, Germany; Heidelberg University, Faculty of Physics and Astronomy, INF 227, D-69120 Heidelberg, Germany
| | - A Hahn
- Heidelberg University Hospital, INF 400, D-69120 Heidelberg, Germany
| | - J M E Jende
- Heidelberg University Hospital, INF 400, D-69120 Heidelberg, Germany
| | - M Bendszus
- Heidelberg University Hospital, INF 400, D-69120 Heidelberg, Germany
| | - S Heiland
- Heidelberg University Hospital, INF 400, D-69120 Heidelberg, Germany
| | - C H Ziener
- Heidelberg University Hospital, INF 400, D-69120 Heidelberg, Germany; German Cancer Research Center, INF 280, D-69120 Heidelberg, Germany.
| |
Collapse
|
25
|
Microstructural imaging of human neocortex in vivo. Neuroimage 2018; 182:184-206. [DOI: 10.1016/j.neuroimage.2018.02.055] [Citation(s) in RCA: 75] [Impact Index Per Article: 10.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2017] [Revised: 02/13/2018] [Accepted: 02/26/2018] [Indexed: 12/12/2022] Open
|
26
|
Alonso-Ortiz E, Levesque IR, Pike GB. Impact of magnetic susceptibility anisotropy at 3 T and 7 T on T2*-based myelin water fraction imaging. Neuroimage 2018; 182:370-378. [DOI: 10.1016/j.neuroimage.2017.09.040] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2017] [Revised: 08/14/2017] [Accepted: 09/20/2017] [Indexed: 10/18/2022] Open
|
27
|
Lee H, Nam Y, Kim D. Echo time‐range effects on gradient‐echo based myelin water fraction mapping at 3T. Magn Reson Med 2018; 81:2799-2807. [DOI: 10.1002/mrm.27564] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2018] [Revised: 09/15/2018] [Accepted: 09/19/2018] [Indexed: 12/19/2022]
Affiliation(s)
- Hongpyo Lee
- Department of Electrical and Electronic Engineering Yonsei University Seoul Korea
| | - Yoonho Nam
- Department of Radiology Seoul St. Mary’s Hospital, College of Medicine, the Catholic University of Korea Seoul Korea
| | - Dong‐Hyun Kim
- Department of Electrical and Electronic Engineering Yonsei University Seoul Korea
| |
Collapse
|
28
|
Fouladivanda M, Kazemi K, Helfroush MS, Shakibafard A. Morphological active contour driven by local and global intensity fitting for spinal cord segmentation from MR images. J Neurosci Methods 2018; 308:116-128. [DOI: 10.1016/j.jneumeth.2018.07.015] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2017] [Revised: 07/18/2018] [Accepted: 07/18/2018] [Indexed: 10/28/2022]
|
29
|
Foxley S, Karczmar GS, Takahashi K. The effects of variations in tissue microstructure from postmortem rat brain on the asymmetry of the water proton resonance. Magn Reson Med 2018; 81:79-89. [PMID: 30014543 DOI: 10.1002/mrm.27338] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2017] [Revised: 03/23/2018] [Accepted: 04/10/2018] [Indexed: 12/13/2022]
Abstract
PURPOSE This work was performed to investigate the effects of tissue microstructure from postmortem rat brain on the shape of the water proton spectrum. METHODS Perfusion-fixed, resected rat brains (N = 4) were imaged at 9.4T. 3D DTI and 3D echo-planar spectroscopic imaging (EPSI) data were acquired with 150 μm isotropic resolution. DTI data were acquired over 60 directions with b = 3000 s/mm2 . Water spectra were produced from EPSI data acquired over 128 echoes, with 2.9 Hz spectral resolution. A voxel-wise metric reflecting spectral asymmetry about the peak of the resonance was computed and compared with orientation estimates from DTI data by fitting data with the susceptibility anisotropy model. RESULTS Asymmetric broadening of the water resonance was computed for mixed populations of grey and/or white matter as determined by thresholding the fractional anisotropy. Asymmetry was shown to be differentially affected by tract orientation relative to B0 in high FA voxels, whereas low FA voxels exhibited little sensitivity. Anatomic structures in the hippocampus were also found to produce distinct changes in the water resonance. CONCLUSION Present results demonstrate that structural variations in tissue architecture cause characteristic, reproducible changes in the water resonance shape. This suggests that water spectra are sensitive to cytoarchitectural variations in brain tissue.
Collapse
Affiliation(s)
- Sean Foxley
- MRI Research Center, Department of Radiology, University of Chicago, Chicago, Illinois
| | - Gregory S Karczmar
- MRI Research Center, Department of Radiology, University of Chicago, Chicago, Illinois
| | - Kazutaka Takahashi
- Department of Organismal Biology and Anatomy, University of Chicago, Chicago, Illinois
| |
Collapse
|
30
|
Kadamangudi S, Reutens D, Sood S, Vegh V. Signal compartments in ultra-high field multi-echo gradient echo MRI reflect underlying tissue microstructure in the brain. Neuroimage 2018; 178:403-413. [PMID: 29852284 DOI: 10.1016/j.neuroimage.2018.05.061] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2018] [Revised: 05/24/2018] [Accepted: 05/25/2018] [Indexed: 10/14/2022] Open
Abstract
Gradient recalled echo magnetic resonance imaging (GRE-MRI) at ultra-high field holds great promise for new contrast mechanisms and delineation of putative tissue compartments that contribute to the multi-echo GRE-MRI signal may aid structural characterization. Several studies have adopted the three water-pool compartment model to study white matter brain regions, associating individual compartments with myelin, axonal and extracellular water. However, the number and identifiability of GRE-MRI signal compartments has not been fully explored. We undertook this task for human brain imaging data. Multiple echo time GRE-MRI data were acquired in five healthy participants, specific anatomical structures were segmented in each dataset (substantia nigra, caudate, insula, putamen, thalamus, fornix, internal capsule, corpus callosum and cerebrospinal fluid), and the signal fitted with models comprising one to six signal compartments using a complex-valued plane wave formulation. Information criteria and cluster analysis methods were used to ascertain the number of distinct compartments within the signal from each structure and to determine their respective frequency shifts. We identified five principal signal compartments with different relative contributions to each structure's signal. Voxel-based maps of the volume fraction of each of these compartments were generated and demonstrated spatial correlation with brain anatomy.
Collapse
Affiliation(s)
- Shrinath Kadamangudi
- Centre for Advanced Imaging, University of Queensland, Brisbane, Queensland, Australia; Queensland Brain Institute, University of Queensland, Brisbane, Queensland, Australia
| | - David Reutens
- Centre for Advanced Imaging, University of Queensland, Brisbane, Queensland, Australia
| | - Surabhi Sood
- Centre for Advanced Imaging, University of Queensland, Brisbane, Queensland, Australia
| | - Viktor Vegh
- Centre for Advanced Imaging, University of Queensland, Brisbane, Queensland, Australia.
| |
Collapse
|
31
|
MRI gradient-echo phase contrast of the brain at ultra-short TE with off-resonance saturation. Neuroimage 2018; 175:1-11. [PMID: 29604452 DOI: 10.1016/j.neuroimage.2018.03.066] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/27/2017] [Revised: 03/07/2018] [Accepted: 03/28/2018] [Indexed: 02/06/2023] Open
Abstract
Larmor-frequency shift or image phase measured by gradient-echo sequences has provided a new source of MRI contrast. This contrast is being used to study both the structure and function of the brain. So far, phase images of the brain have been largely obtained at long echo times as maximum phase signal-to-noise ratio (SNR) is achieved at TE = T2* (∼40 ms at 3T). The structures of the brain, however, are compartmentalized and complex with a wide range of signal relaxation times. At such long TE, the short-T2 components are largely attenuated and contribute minimally to phase contrast. The purpose of this study was to determine whether proton gradient-echo images of the brain exhibit phase contrast at ultra-short TE (UTE). Our data showed that UTE images acquired at 7 T without off-resonance saturation do not contain significant phase contrast between gray and white matter. However, UTE images of the brain can attain strong phase contrast even at a nominal TE of 106 μs by using off-resonance RF saturation pulses, which provide direct saturation of ultra-short-T2 components and indirect saturation of longer-T2 components via magnetization transfer. In addition, phase contrast between gray and white matter acquired at UTE with off-resonance saturation is reversed compared to that of the long-T2 signals acquired at long TEs. This finding opens up a potential new way to manipulate image phase contrast of the brain. By accessing short and ultra-short-T2 species, MRI phase images may further improve the characterization of tissue microstructure in the brain.
Collapse
|
32
|
Ziser L, Meyer-Schell N, Kurniawan ND, Sullivan R, Reutens D, Chen M, Vegh V. Utility of gradient recalled echo magnetic resonance imaging for the study of myelination in cuprizone mice treated with fingolimod. NMR IN BIOMEDICINE 2018; 31:e3877. [PMID: 29266540 DOI: 10.1002/nbm.3877] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/02/2017] [Revised: 11/10/2017] [Accepted: 11/10/2017] [Indexed: 05/26/2023]
Abstract
The availability of high-field-strength magnetic resonance imaging (MRI) systems has brought about the development of techniques that aim to map myelination via the exploitation of various contrast mechanisms. Myelin mapping techniques have the potential to provide tools for the diagnosis and treatment of diseases, such as multiple sclerosis. In this study, we evaluated the sensitivity of T2 *, frequency shift and susceptibility measures to myelin levels in a cuprizone mouse model of demyelination. The model was supplemented with two different dosages of fingolimod, a drug known to positively affect demyelination. A decrease in grey-white matter contrast with the cuprizone diet was observed for T2 *, frequency shift and susceptibility measures, together with myelin basic protein antibody findings. These results indicate that T2 *, frequency shift and susceptibility measures have the potential to act as biomarkers for myelination. Susceptibility was found to be the most sensitive measure to changes in grey-white matter contrast. In addition, fingolimod treatment was found to reduce the level of demyelination, with a larger dosage exhibiting a greater reduction in demyelination for the in vivo MRI results. Overall, susceptibility mapping appears to be a more promising tool than T2 * or frequency shift mapping for the early diagnosis and treatment of diseases in which myelination is implicated.
Collapse
Affiliation(s)
- Laura Ziser
- Centre for Advanced Imaging, University of Queensland, Brisbane, Qld, Australia
| | - Naja Meyer-Schell
- Centre for Advanced Imaging, University of Queensland, Brisbane, Qld, Australia
| | - Nyoman D Kurniawan
- Centre for Advanced Imaging, University of Queensland, Brisbane, Qld, Australia
| | - Robert Sullivan
- Queensland Brain Institute, University of Queensland, Brisbane, Qld, Australia
| | - David Reutens
- Centre for Advanced Imaging, University of Queensland, Brisbane, Qld, Australia
| | - Min Chen
- Centre for Advanced Imaging, University of Queensland, Brisbane, Qld, Australia
| | - Viktor Vegh
- Centre for Advanced Imaging, University of Queensland, Brisbane, Qld, Australia
| |
Collapse
|
33
|
Schyboll F, Jaekel U, Weber B, Neeb H. The impact of fibre orientation on T1-relaxation and apparent tissue water content in white matter. MAGNETIC RESONANCE MATERIALS IN PHYSICS BIOLOGY AND MEDICINE 2018; 31:501-510. [DOI: 10.1007/s10334-018-0678-8] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/26/2017] [Revised: 02/07/2018] [Accepted: 02/07/2018] [Indexed: 11/29/2022]
|
34
|
Thapaliya K, Vegh V, Bollmann S, Barth M. Assessment of microstructural signal compartments across the corpus callosum using multi-echo gradient recalled echo at 7 T. Neuroimage 2017; 182:407-416. [PMID: 29183776 DOI: 10.1016/j.neuroimage.2017.11.029] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2017] [Revised: 11/06/2017] [Accepted: 11/15/2017] [Indexed: 12/18/2022] Open
Abstract
Quantitative assessment of tissue microstructure is important in studying human brain diseases and disorders in which white matter is implicated, as it has been linked to demyelination, re-myelination, and axonal damage in clinical conditions. Ultra-high field magnetic resonance imaging data obtained using a multi-echo gradient echo sequence has been shown to contain information on myelin, axonal and extracellular compartments in white matter. In this study, we aimed to assess the sensitivity of a three-compartment model to estimate the variation of corresponding compartment parameters (water fraction, relaxation time and frequency shift) of the corpus callosum sub-regions, which are known to have different tissue structure. Additionally, we computed the g-ratio using myelin and axonal water fractions and performed a voxel-by-voxel analysis in the corpus callosum. Based on data acquired for ten participants, we show that the myelin compartment water fraction and T2∗ is consistent across the corpus callosum sub-regions, whilst myelin frequency shift varies. The results show that the variation in water fraction, T2∗ and frequency shift for the myelin signal compartment across the corpus callosum is smaller than for the axonal and extracellular signal compartments. The computed g-ratio was comparable to previously published studies in the corpus callosum. Our study suggests that a multi-echo GRE approach in vivo combined with a complex three-compartment model is sensitive to microstructural parameter variations across the human corpus callosum.
Collapse
Affiliation(s)
- Kiran Thapaliya
- Centre for Advanced Imaging, The University of Queensland, Queensland, Australia
| | - Viktor Vegh
- Centre for Advanced Imaging, The University of Queensland, Queensland, Australia
| | - Steffen Bollmann
- Centre for Advanced Imaging, The University of Queensland, Queensland, Australia
| | - Markus Barth
- Centre for Advanced Imaging, The University of Queensland, Queensland, Australia.
| |
Collapse
|
35
|
An R2* model of white matter for fiber orientation and myelin concentration. Neuroimage 2017; 162:269-275. [DOI: 10.1016/j.neuroimage.2017.08.050] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/09/2017] [Revised: 07/10/2017] [Accepted: 08/16/2017] [Indexed: 01/31/2023] Open
|
36
|
Kaller MS, Lazari A, Blanco-Duque C, Sampaio-Baptista C, Johansen-Berg H. Myelin plasticity and behaviour-connecting the dots. Curr Opin Neurobiol 2017; 47:86-92. [PMID: 29054040 PMCID: PMC5844949 DOI: 10.1016/j.conb.2017.09.014] [Citation(s) in RCA: 69] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2017] [Revised: 08/29/2017] [Accepted: 09/21/2017] [Indexed: 12/31/2022]
Abstract
Changes in white matter and myelin are associated with learning during adulthood across species. The causal link between myelin plasticity and behaviour remains elusive. Preventing the differentiation of new OLs can impair learning within the first few hours. Myelin remodelling may occur through many different routes and mechanism. The functional arrangement of myelination along axons can be complex and diverse.
Myelin sheaths in the vertebrate nervous system enable faster impulse propagation, while myelinating glia provide vital support to axons. Once considered a static insulator, converging evidence now suggests that myelin in the central nervous system can be dynamically regulated by neuronal activity and continues to participate in nervous system plasticity beyond development. While the link between experience and myelination gains increased recognition, it is still unclear what role such adaptive myelination plays in facilitating and shaping behaviour. Additionally, fundamental mechanisms and principles underlying myelin remodelling remain poorly understood. In this review, we will discuss new insights into the link between myelin plasticity and behaviour, as well as mechanistic aspects of myelin remodelling that may help to elucidate this intriguing process.
Collapse
Affiliation(s)
- Malte Sebastian Kaller
- Nuffield Department of Clinical Neurosciences, John Radcliffe Hospital, University of Oxford, Oxford OX1 2JD, United Kingdom.
| | - Alberto Lazari
- Nuffield Department of Clinical Neurosciences, John Radcliffe Hospital, University of Oxford, Oxford OX1 2JD, United Kingdom
| | - Cristina Blanco-Duque
- Nuffield Department of Clinical Neurosciences, John Radcliffe Hospital, University of Oxford, Oxford OX1 2JD, United Kingdom
| | - Cassandra Sampaio-Baptista
- Nuffield Department of Clinical Neurosciences, John Radcliffe Hospital, University of Oxford, Oxford OX1 2JD, United Kingdom
| | - Heidi Johansen-Berg
- Nuffield Department of Clinical Neurosciences, John Radcliffe Hospital, University of Oxford, Oxford OX1 2JD, United Kingdom
| |
Collapse
|
37
|
Álvarez GA, Shemesh N, Frydman L. Internal gradient distributions: A susceptibility-derived tensor delivering morphologies by magnetic resonance. Sci Rep 2017; 7:3311. [PMID: 28607445 PMCID: PMC5468317 DOI: 10.1038/s41598-017-03277-9] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2016] [Accepted: 04/26/2017] [Indexed: 11/23/2022] Open
Abstract
Nuclear magnetic resonance is a powerful tool for probing the structures of chemical and biological systems. Combined with field gradients it leads to NMR imaging (MRI), a widespread tool in non-invasive examinations. Sensitivity usually limits MRI’s spatial resolution to tens of micrometers, but other sources of information like those delivered by constrained diffusion processes, enable one extract morphological information down to micron and sub-micron scales. We report here on a new method that also exploits diffusion – isotropic or anisotropic– to sense morphological parameters in the nm-mm range, based on distributions of susceptibility-induced magnetic field gradients. A theoretical framework is developed to define this source of information, leading to the proposition of internal gradient-distribution tensors. Gradient-based spin-echo sequences are designed to measure these new observables. These methods can be used to map orientations even when dealing with unconstrained diffusion, as is here demonstrated with studies of structured systems, including tissues.
Collapse
Affiliation(s)
- Gonzalo A Álvarez
- Department of Chemical Physics, Weizmann Institute of Science, Rehovot, 76100, Israel.,Centro Atómico Bariloche, CONICET, CNEA, 8400, S. C. de Bariloche, Argentina
| | - Noam Shemesh
- Department of Chemical Physics, Weizmann Institute of Science, Rehovot, 76100, Israel.,Champalimaud Neuroscience Programme, Champalimaud Centre for the Unknown, Lisbon, 1400-138, Portugal
| | - Lucio Frydman
- Department of Chemical Physics, Weizmann Institute of Science, Rehovot, 76100, Israel.
| |
Collapse
|
38
|
Farrher E, Lindemeyer J, Grinberg F, Oros-Peusquens AM, Shah NJ. Concerning the matching of magnetic susceptibility differences for the compensation of background gradients in anisotropic diffusion fibre phantoms. PLoS One 2017; 12:e0176192. [PMID: 28467458 PMCID: PMC5415224 DOI: 10.1371/journal.pone.0176192] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2017] [Accepted: 04/06/2017] [Indexed: 11/20/2022] Open
Abstract
Artificial, anisotropic fibre phantoms are nowadays increasingly used in the field of diffusion-weighted MRI. Such phantoms represent useful tools for, among others, the calibration of pulse sequences and validation of diffusion models since they can mimic well-known structural features of brain tissue on the one hand, but exhibit a reduced complexity, on the other. Among all materials, polyethylene fibres have been widely used due to their excellent properties regarding the restriction of water diffusion and surface relaxation properties. Yet the magnetic susceptibility of polyethylene can be distinctly lower than that of distilled water. This difference produces strong microscopic, background field gradients in the vicinity of fibre bundles which are not parallel to the static magnetic field. This, in turn, modulates the MRI signal behaviour. In the present work we investigate an approach to reduce the susceptibility-induced background gradients via reducing the heterogeneity in the internal magnetic susceptibility. An aqueous solution of magnesium chloride hexahydrate (MgCl2·6H2O) is used for this purpose. Its performance is demonstrated in dedicated anisotropic fibre phantoms with different geometrical configurations.
Collapse
Affiliation(s)
- Ezequiel Farrher
- Institute of Neuroscience and Medicine – 4, Forschungszentrum Jülich GmbH, Jülich, Germany
- * E-mail:
| | - Johannes Lindemeyer
- Institute of Neuroscience and Medicine – 4, Forschungszentrum Jülich GmbH, Jülich, Germany
| | - Farida Grinberg
- Institute of Neuroscience and Medicine – 4, Forschungszentrum Jülich GmbH, Jülich, Germany
- Department of Neurology, Faculty of Medicine, RWTH Aachen University, Aachen, Germany
| | | | - N. Jon Shah
- Institute of Neuroscience and Medicine – 4, Forschungszentrum Jülich GmbH, Jülich, Germany
- Department of Neurology, Faculty of Medicine, RWTH Aachen University, Aachen, Germany
- JARA – BRAIN – Translational Medicine, RWTH Aachen University, Aachen, Germany
- Institute of Neuroscience and Medicine – 11, Forschungszentrum Jülich GmbH, Jülich, Germany
- Department of Electrical and Computer Systems Engineering, and Monash Biomedical Imaging, School of Psychological Sciences, Monash University, Melbourne, Victoria, Australia
| |
Collapse
|
39
|
Xu T, Foxley S, Kleinnijenhuis M, Chen WC, Miller KL. The effect of realistic geometries on the susceptibility-weighted MR signal in white matter. Magn Reson Med 2017; 79:489-500. [PMID: 28394030 PMCID: PMC6585669 DOI: 10.1002/mrm.26689] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2016] [Revised: 02/05/2017] [Accepted: 03/04/2017] [Indexed: 12/15/2022]
Abstract
PURPOSE To investigate the effect of realistic microstructural geometry on the susceptibility-weighted MR signal in white matter (WM), with application to demyelination. METHODS Previous work has modeled susceptibility-weighted signals under the assumption that axons are cylindrical. In this study, we explored the implications of this assumption by considering the effect of more realistic geometries. A three-compartment WM model incorporating relevant properties based on the literature was used to predict the MR signal. Myelinated axons were modeled with several cross-sectional geometries of increasing realism: nested circles, warped/elliptical circles, and measured axonal geometries from electron micrographs. Signal simulations from the different microstructural geometries were compared with measured signals from a cuprizone mouse model with varying degrees of demyelination. RESULTS Simulation results suggest that axonal geometry affects the MR signal. Predictions with realistic models were significantly different compared with circular models under the same microstructural tissue properties, for simulations with and without diffusion. CONCLUSION The geometry of axons affects the MR signal significantly. Literature estimates of myelin susceptibility, which are based on fitting biophysical models to the MR signal, are likely to be biased by the assumed geometry, as will any derived microstructural properties. Magn Reson Med 79:489-500, 2018. © 2017 International Society for Magnetic Resonance in Medicine.
Collapse
Affiliation(s)
- Tianyou Xu
- Oxford Centre for Functional MRI of the Brain, University of Oxford, Oxford, United Kingdom
| | - Sean Foxley
- Oxford Centre for Functional MRI of the Brain, University of Oxford, Oxford, United Kingdom
| | - Michiel Kleinnijenhuis
- Oxford Centre for Functional MRI of the Brain, University of Oxford, Oxford, United Kingdom
| | - Way Cherng Chen
- Oxford Centre for Functional MRI of the Brain, University of Oxford, Oxford, United Kingdom.,Singapore Bioimaging Consortium, A*STAR, Singapore
| | - Karla L Miller
- Oxford Centre for Functional MRI of the Brain, University of Oxford, Oxford, United Kingdom
| |
Collapse
|
40
|
Lee J, Nam Y, Choi JY, Kim EY, Oh SH, Kim DH. Mechanisms of T 2 * anisotropy and gradient echo myelin water imaging. NMR IN BIOMEDICINE 2017; 30:e3513. [PMID: 27060968 DOI: 10.1002/nbm.3513] [Citation(s) in RCA: 35] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/02/2015] [Revised: 01/26/2016] [Accepted: 02/17/2016] [Indexed: 06/05/2023]
Abstract
In MRI, structurally aligned molecular or micro-organization (e.g. axonal fibers) can be a source of substantial signal variations that depend on the structural orientation and the applied magnetic field. This signal anisotropy gives us a unique opportunity to explore information that exists at a resolution several orders of magnitude smaller than that of typical MRI. In this review, one of the signal anisotropies, T2 * anisotropy in white matter, and a related imaging method, gradient echo myelin water imaging (GRE-MWI), are explored. The T2 * anisotropy has been attributed to isotropic and anisotropic magnetic susceptibility of myelin and compartmentalized microstructure of white matter fibers (i.e. axonal, myelin, and extracellular space). The susceptibility and microstructure create magnetic frequency shifts that change with the relative orientation of the fiber and the main magnetic field, generating the T2 * anisotropy. The resulting multi-component magnitude decay and nonlinear phase evolution have been utilized for GRE-MWI, assisting in resolving the signal fraction of the multiple compartments in white matter. The GRE-MWI method has been further improved by signal compensation techniques including physiological noise compensation schemes. The T2 * anisotropy and GRE-MWI provide microstructural information on a voxel (e.g. fiber orientation and tissue composition), and may serve as sensitive biomarkers for microstructural changes in the brain. Copyright © 2016 John Wiley & Sons, Ltd.
Collapse
Affiliation(s)
- Jongho Lee
- Laboratory for Imaging Science and Technology, Department of Electrical and Computer Engineering, Seoul National University, Seoul, Korea
| | - Yoonho Nam
- Department of Radiology, Seoul St. Mary's Hospital, College of Medicine, The Catholic University of Korea, Seoul, Korea
| | - Joon Yul Choi
- Laboratory for Imaging Science and Technology, Department of Electrical and Computer Engineering, Seoul National University, Seoul, Korea
| | - Eung Yeop Kim
- Department of Radiology, Gachon University Gil Medical Center, Incheon, Korea
| | - Se-Hong Oh
- Imaging Institute, Cleveland Clinic, Cleveland, OH, USA
| | - Dong-Hyun Kim
- Department of Electrical and Electronic Engineering, Yonsei University, Seoul, Korea
| |
Collapse
|
41
|
Nunes D, Cruz TL, Jespersen SN, Shemesh N. Mapping axonal density and average diameter using non-monotonic time-dependent gradient-echo MRI. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2017; 277:117-130. [PMID: 28282586 DOI: 10.1016/j.jmr.2017.02.017] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/03/2017] [Revised: 02/23/2017] [Accepted: 02/25/2017] [Indexed: 06/06/2023]
Abstract
White Matter (WM) microstructures, such as axonal density and average diameter, are crucial to the normal function of the Central Nervous System (CNS) as they are closely related with axonal conduction velocities. Conversely, disruptions of these microstructural features may result in severe neurological deficits, suggesting that their noninvasive mapping could be an important step towards diagnosing and following pathophysiology. Whereas diffusion based MRI methods have been proposed to map these features, they typically entail the application of powerful gradients, which are rarely available in the clinic, or extremely long acquisition schemes to extract information from parameter-intensive models. In this study, we suggest that simple and time-efficient multi-gradient-echo (MGE) MRI can be used to extract the axon density from susceptibility-driven non-monotonic decay in the time-dependent signal. We show, both theoretically and with simulations, that a non-monotonic signal decay will occur for multi-compartmental microstructures - such as axons and extra-axonal spaces, which were here used as a simple model for the microstructure - and that, for axons parallel to the main magnetic field, the axonal density can be extracted. We then experimentally demonstrate in ex-vivo rat spinal cords that its different tracts - characterized by different microstructures - can be clearly contrasted using the MGE-derived maps. When the quantitative results are compared against ground-truth histology, they reflect the axonal fraction (though with a bias, as evident from Bland-Altman analysis). As well, the extra-axonal fraction can be estimated. The results suggest that our model is oversimplified, yet at the same time evidencing a potential and usefulness of the approach to map underlying microstructures using a simple and time-efficient MRI sequence. We further show that a simple general-linear-model can predict the average axonal diameters from the four model parameters, and map these average axonal diameters in the spinal cords. While clearly further modelling and theoretical developments are necessary, we conclude that salient WM microstructural features can be extracted from simple, SNR-efficient multi-gradient echo MRI, and that this paves the way towards easier estimation of WM microstructure in vivo.
Collapse
Affiliation(s)
- Daniel Nunes
- Champalimaud Neuroscience Programme, Champalimaud Centre for the Unknown, Av. Brasilia 1400-038, Lisbon, Portugal
| | - Tomás L Cruz
- Champalimaud Neuroscience Programme, Champalimaud Centre for the Unknown, Av. Brasilia 1400-038, Lisbon, Portugal
| | - Sune N Jespersen
- Center of Functionally Integrative Neuroscience (CFIN) and MINDLab, Clinical Institute, Aarhus University, Aarhus, Denmark; Department of Physics and Astronomy, Aarhus University, Aarhus, Denmark
| | - Noam Shemesh
- Champalimaud Neuroscience Programme, Champalimaud Centre for the Unknown, Av. Brasilia 1400-038, Lisbon, Portugal.
| |
Collapse
|
42
|
He Q, Ma Y, Fan S, Shao H, Sheth V, Bydder GM, Du J. Direct magnitude and phase imaging of myelin using ultrashort echo time (UTE) pulse sequences: A feasibility study. Magn Reson Imaging 2017; 39:194-199. [PMID: 28219648 DOI: 10.1016/j.mri.2017.02.009] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/31/2016] [Revised: 02/01/2017] [Accepted: 02/16/2017] [Indexed: 12/19/2022]
Abstract
In this paper, we aimed to investigate the feasibility of direct visualization of myelin, including myelin lipid and myelin basic protein (MBP), using two-dimensional ultrashort echo time (2D UTE) sequences and utilize phase information as a contrast mechanism in phantoms and in volunteers. The standard UTE sequence was used to detect both myelin and long T2 signal. An adiabatic inversion recovery UTE (IR-UTE) sequence was used to selectively detect myelin by suppressing signal from long T2 water protons. Magnitude and phase imaging and T2* were investigated on myelin lipid and MBP in the forms of lyophilized powders as well as paste-like phantoms with the powder mixed with D2O, and rubber phantoms as well as healthy volunteers. Contrast to noise ratio (CNR) between white and gray matter was measured. Both magnitude and phase images were generated for myelin and rubber phantoms as well white matter in vivo using the IR-UTE sequence. T2* values of ~300μs were comparable for myelin paste phantoms and the short T2* component in white matter of the brain in vivo. Mean CNR between white and gray matter in IR-UTE imaging was increased from -7.3 for the magnitude images to 57.4 for the phase images. The preliminary results suggest that the IR-UTE sequence allows simultaneous magnitude and phase imaging of myelin in vitro and in vivo.
Collapse
Affiliation(s)
- Qun He
- Department of Radiology, University of California, San Diego, United States; Ningbo Jansen NMR Technology Co., Ltd., Cixi, Zhejiang Province, China
| | - Yajun Ma
- Department of Radiology, University of California, San Diego, United States
| | - Shujuan Fan
- Department of Radiology, University of California, San Diego, United States
| | - Hongda Shao
- Department of Radiology, University of California, San Diego, United States
| | - Vipul Sheth
- Department of Radiology, University of California, San Diego, United States
| | - Graeme M Bydder
- Department of Radiology, University of California, San Diego, United States
| | - Jiang Du
- Department of Radiology, University of California, San Diego, United States.
| |
Collapse
|
43
|
Cronin MJ, Wang N, Decker KS, Wei H, Zhu WZ, Liu C. Exploring the origins of echo-time-dependent quantitative susceptibility mapping (QSM) measurements in healthy tissue and cerebral microbleeds. Neuroimage 2017; 149:98-113. [PMID: 28126551 DOI: 10.1016/j.neuroimage.2017.01.053] [Citation(s) in RCA: 60] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2016] [Revised: 01/19/2017] [Accepted: 01/22/2017] [Indexed: 12/20/2022] Open
Abstract
Quantitative susceptibility mapping (QSM) is increasingly used to measure variation in tissue composition both in the brain and in other areas of the body in a range of disease pathologies. Although QSM measurements were originally believed to be independent of the echo time (TE) used in the gradient-recalled echo (GRE) acquisition from which they are derived; recent literature (Sood et al., 2016) has shown that these measurements can be highly TE-dependent in a number of brain regions. In this work we systematically investigate possible causes of this effect through analysis of apparent frequency and QSM measurements derived from data acquired at multiple TEs in vivo in healthy brain regions and in cerebral microbleeds (CMBs); QSM data acquired in a gadolinium-doped phantom; and in QSM data derived from idealized simulated phase data. Apparent frequency measurements in the optic radiations (OR) and central corpus callosum (CC) were compared to those predicted by a 3-pool white matter model, however the model failed to fully explain contrasting frequency profiles measured in the OR and CC. Our results show that TE-dependent QSM measurements can be caused by a failure of phase unwrapping algorithms in and around strong susceptibility sources such as CMBs; however, in healthy brain regions this behavior appears to result from intrinsic non-linear phase evolution in the MR signal. From these results we conclude that care must be taken when deriving frequency and QSM measurements in strong susceptibility sources due to the inherent limitations in phase unwrapping; and that while signal compartmentalization due to tissue microstructure and content is a plausible cause of TE-dependent frequency and QSM measurements in healthy brain regions, better sampling of the MR signal and more complex models of tissue are needed to fully exploit this relationship.
Collapse
Affiliation(s)
- Matthew J Cronin
- Department of Electrical Engineering and Computer Sciences, University of California, Berkeley, CA 94720, USA; Brain Imaging and Analysis Center, Duke University, Durham, NC 27710, USA
| | - Nian Wang
- Brain Imaging and Analysis Center, Duke University, Durham, NC 27710, USA
| | - Kyle S Decker
- Brain Imaging and Analysis Center, Duke University, Durham, NC 27710, USA
| | - Hongjiang Wei
- Department of Electrical Engineering and Computer Sciences, University of California, Berkeley, CA 94720, USA; Brain Imaging and Analysis Center, Duke University, Durham, NC 27710, USA
| | - Wen-Zhen Zhu
- Department of Radiology, Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
| | - Chunlei Liu
- Department of Electrical Engineering and Computer Sciences, University of California, Berkeley, CA 94720, USA; Brain Imaging and Analysis Center, Duke University, Durham, NC 27710, USA.
| |
Collapse
|
44
|
Caporale A, Palombo M, Macaluso E, Guerreri M, Bozzali M, Capuani S. The γ-parameter of anomalous diffusion quantified in human brain by MRI depends on local magnetic susceptibility differences. Neuroimage 2016; 147:619-631. [PMID: 28011255 DOI: 10.1016/j.neuroimage.2016.12.051] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2016] [Revised: 11/22/2016] [Accepted: 12/19/2016] [Indexed: 12/15/2022] Open
Abstract
Motivated by previous results obtained in vitro, we investigated the dependence of the anomalous diffusion (AD) MRI technique on local magnetic susceptibility differences (Δχ) driven by magnetic field inhomogeneity in human brains. The AD-imaging contrast investigated here is quantified by the stretched-exponential parameter γ, extracted from diffusion weighted (DW) data collected by varying diffusion gradient strengths. We performed T2* and DW experiments in eight healthy subjects at 3.0T. T2*-weighted images at different TEs=(10,20,35,55)ms and DW-EPI images with fourteen b-values from 0 to 5000s/mm2 were acquired. AD-metrics and Diffusion Tensor Imaging (DTI) parameters were compared and correlated to R2* and to Δχ values taken from literature for the gray (GM) and the white (WM) matter. Pearson's correlation test and Analysis of Variance with Bonferroni post-hoc test were used. Significant strong linear correlations were found between AD γ-metrics and R2* in both GM and WM of the human brain, but not between DTI-metrics and R2*. Depending on Δχ driven magnetic field inhomogeneity, the new contrast provided by AD-γ imaging reflects Δχ due to differences in myelin orientation and iron content within selected regions in the WM and GM, respectively. This feature of the AD-γ imaging due to the fact that γ is quantified by using MRI, may be an alternative strategy to investigate, at high magnetic fields, microstructural changes in myelin, and alterations due to iron accumulation. Possible clinical applications might be in the field of neurodegenerative diseases.
Collapse
Affiliation(s)
- A Caporale
- Morpho-functional Sciences, Department of Anatomical, Histological, Forensic and of the Locomotor System Science, Sapienza University of Rome, Italy; CNR ISC UOS Roma Sapienza, Physics Department Sapienza University of Rome, Rome, Italy.
| | - M Palombo
- CNR ISC UOS Roma Sapienza, Physics Department Sapienza University of Rome, Rome, Italy; MIRCen, CEA/DSV/I(2)BM, Fontenay-aux-Roses, France
| | - E Macaluso
- ImpAct Team, Lyon Neuroscience Research Center, Lyon, France
| | - M Guerreri
- CNR ISC UOS Roma Sapienza, Physics Department Sapienza University of Rome, Rome, Italy; Morphogenesis & Tissue Engineering, Department of Anatomical, Histological, Forensic and of the Locomotor System Science, Sapienza University of Rome, Italy
| | - M Bozzali
- Neuroimaging Laboratory Santa Lucia Foundation, Rome, Italy
| | - S Capuani
- CNR ISC UOS Roma Sapienza, Physics Department Sapienza University of Rome, Rome, Italy
| |
Collapse
|
45
|
van der Zwaag W, Schäfer A, Marques JP, Turner R, Trampel R. Recent applications of UHF-MRI in the study of human brain function and structure: a review. NMR IN BIOMEDICINE 2016; 29:1274-1288. [PMID: 25762497 DOI: 10.1002/nbm.3275] [Citation(s) in RCA: 63] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/01/2014] [Revised: 12/19/2014] [Accepted: 01/22/2015] [Indexed: 06/04/2023]
Abstract
The increased availability of ultra-high-field (UHF) MRI has led to its application in a wide range of neuroimaging studies, which are showing promise in transforming fundamental approaches to human neuroscience. This review presents recent work on structural and functional brain imaging, at 7 T and higher field strengths. After a short outline of the effects of high field strength on MR images, the rapidly expanding literature on UHF applications of blood-oxygenation-level-dependent-based functional MRI is reviewed. Structural imaging is then discussed, divided into sections on imaging weighted by relaxation time, including quantitative relaxation time mapping, phase imaging and quantitative susceptibility mapping, angiography, diffusion-weighted imaging, and finally magnetization-transfer imaging. The final section discusses studies using the high spatial resolution available at UHF to identify explicit links between structure and function. Copyright © 2015 John Wiley & Sons, Ltd.
Collapse
Affiliation(s)
- Wietske van der Zwaag
- Centre d'Imagerie Biomédicale, Ecole Polytechnique Fédérale de Lausanne, Switzerland
| | - Andreas Schäfer
- Max Planck Institute for Human Cognitive and Brain Sciences, Leipzig, Germany
| | - José P Marques
- Centre d'Imagerie Biomédicale, Ecole Polytechnique Fédérale de Lausanne, Switzerland
- Donders Institute for Brain, Cognition and Behaviour, Radboud University, Nijmegen, The Netherlands
| | - Robert Turner
- Max Planck Institute for Human Cognitive and Brain Sciences, Leipzig, Germany
- Spinoza Centre, University of Amsterdam, The Netherlands
- SPMMRC, School of Physics and Astronomy, University of Nottingham, UK
| | - Robert Trampel
- Max Planck Institute for Human Cognitive and Brain Sciences, Leipzig, Germany
| |
Collapse
|
46
|
Barbara TM. White matter shifts in MRI: Rehabilitating the Lorentz sphere in magnetic resonance. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2016; 270:40-46. [PMID: 27393892 DOI: 10.1016/j.jmr.2016.06.016] [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: 03/16/2016] [Revised: 06/22/2016] [Accepted: 06/25/2016] [Indexed: 06/06/2023]
Abstract
A thorough exposition and analysis of the role of the Lorentz sphere in magnetic resonance is presented from the fundamental standpoint of macroscopic magnetostatics. The analysis will be useful to those interested in understanding susceptibility and chemical shift contributions to frequency shifts in magnetic resonance. Though the topic is mature, recent research on white matter shifts in the brain promotes the notion of replacing the Lorentz sphere with a generalized Lorentzian cylinder, and has put into question the long standing spherical approach when elongated structures are present. The cavity shape issue can be resolved by applying Helmholtz's theorem, which can be expressed in a differential and an integral formulation. The general validity of the Lorentz sphere for any situation is confirmed. Furthermore, a clear exposition of the "generalized approach" is offered, using the language of Lorentz's theory. With the rehabilitation of the Lorentz sphere settled, one must consider alternative contributions to white matter shifts and a likely candidate is the effect of molecular environment on chemical shifts.
Collapse
Affiliation(s)
- Thomas M Barbara
- Advanced Imaging Research Center, Oregon Health and Sciences University, Portland, OR 97239, United States.
| |
Collapse
|
47
|
Sood S, Urriola J, Reutens D, O'Brien K, Bollmann S, Barth M, Vegh V. Echo time-dependent quantitative susceptibility mapping contains information on tissue properties. Magn Reson Med 2016; 77:1946-1958. [PMID: 27221590 DOI: 10.1002/mrm.26281] [Citation(s) in RCA: 49] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2016] [Revised: 04/29/2016] [Accepted: 04/29/2016] [Indexed: 12/12/2022]
Abstract
PURPOSE Magnetic susceptibility is a physical property of matter that varies depending on chemical composition and abundance of different molecular species. Interest is growing in mapping of magnetic susceptibility in the human brain using magnetic resonance imaging techniques, but the influences affecting the mapped values are not fully understood. METHODS We performed quantitative susceptibility mapping on 7 Tesla (T) multiple echo time gradient recalled echo data and evaluated the trend in 10 regions of the human brain. Temporal plots of susceptibility were performed in the caudate, pallidum, putamen, thalamus, insula, red nucleus, substantia nigra, internal capsule, corpus callosum, and fornix. We implemented an existing three compartment signal model and used optimization to fit the experimental result to assess the influences that could be responsible for our findings. RESULTS The temporal trend in susceptibility is different for different brain regions, and subsegmentation of specific regions suggests that differences are likely to be attributable to variations in tissue structure and composition. Using a signal model, we verified that a nonlinear temporal behavior in experimentally computed susceptibility within imaging voxels may be the result of the heterogeneous composition of tissue properties. CONCLUSIONS Decomposition of voxel constituents into meaningful parameters may lead to informative measures that reflect changes in tissue microstructure. Magn Reson Med 77:1946-1958, 2017. © 2016 International Society for Magnetic Resonance in Medicine.
Collapse
Affiliation(s)
- Surabhi Sood
- Centre for Advanced Imaging, the University of Queensland, Brisbane, Australia
| | - Javier Urriola
- Centre for Advanced Imaging, the University of Queensland, Brisbane, Australia
| | - David Reutens
- Centre for Advanced Imaging, the University of Queensland, Brisbane, Australia
| | | | - Steffen Bollmann
- Centre for Advanced Imaging, the University of Queensland, Brisbane, Australia
| | - Markus Barth
- Centre for Advanced Imaging, the University of Queensland, Brisbane, Australia
| | - Viktor Vegh
- Centre for Advanced Imaging, the University of Queensland, Brisbane, Australia
| |
Collapse
|
48
|
Abstract
PURPOSE OF REVIEW Current computational neuroanatomy based on MRI focuses on morphological measures of the brain. We present recent methodological developments in quantitative MRI (qMRI) that provide standardized measures of the brain, which go beyond morphology. We show how biophysical modelling of qMRI data can provide quantitative histological measures of brain tissue, leading to the emerging field of in-vivo histology using MRI (hMRI). RECENT FINDINGS qMRI has greatly improved the sensitivity and specificity of computational neuroanatomy studies. qMRI metrics can also be used as direct indicators of the mechanisms driving observed morphological findings. For hMRI, biophysical models of the MRI signal are being developed to directly access histological information such as cortical myelination, axonal diameters or axonal g-ratio in white matter. Emerging results indicate promising prospects for the combined study of brain microstructure and function. SUMMARY Non-invasive brain tissue characterization using qMRI or hMRI has significant implications for both research and clinics. Both approaches improve comparability across sites and time points, facilitating multicentre/longitudinal studies and standardized diagnostics. hMRI is expected to shed new light on the relationship between brain microstructure, function and behaviour, both in health and disease, and become an indispensable addition to computational neuroanatomy.
Collapse
|
49
|
Aggarwal M, Kageyama Y, Li X, van Zijl PC. B0 -orientation dependent magnetic susceptibility-induced white matter contrast in the human brainstem at 11.7T. Magn Reson Med 2016; 75:2455-63. [PMID: 27018784 DOI: 10.1002/mrm.26208] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2015] [Revised: 02/12/2016] [Accepted: 02/18/2016] [Indexed: 01/17/2023]
Abstract
PURPOSE To investigate B0 -field-orientation dependent white matter contrast in the human brainstem based on R2 * and frequency difference (Δf) mapping from gradient echo (GRE) imaging at 11.7T. METHODS Multi-echo GRE data were acquired from two fixed human brainstem specimens at multiple orientations with respect to the static B0 field. The B0 -orientation dependent modulation curves of R2 * and Δf measurements between short and long echo time regimes were used to reconstruct maps of three-dimensional (3D) white matter orientation vectors. The results were compared with maps from diffusion MRI, susceptibility tensor imaging, and histological staining of the same specimens. RESULTS R2 * and Δf maps demonstrated distinct and significant contrast modulation between the corticospinal tract (CST) and transverse pontine fibers (TPF) dependent on B0 orientation. Interleaved fiber orientations of the CST and TPF could be sensitively resolved based on field-orientation-dependent fitting of the R2 * and Δf measurements. The fitted 3D orientation vector maps and peak-to-peak amplitude of R2 * and Δf modulation exhibited close correspondence to primary eigenvector and anisotropy maps derived from diffusion MRI. The amplitude of B0 -orientation dependent R2 * modulation was significantly (P < 0.005) higher in the CST compared with TPF, while fractional anisotropies were comparable. CONCLUSION The findings of this study demonstrate the potential of B0 -orientation dependent susceptibility-induced R2 * and Δf contrasts to probe tract-specific orientation and microstructure in white matter. Magn Reson Med 75:2455-2463, 2016. © 2016 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Manisha Aggarwal
- Russell H. Morgan Department of Radiology and Radiological Science, Division of MR Research, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA
| | - Yusuke Kageyama
- Department of Pathology, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA
| | - Xu Li
- F. M. Kirby Research Center for Functional Brain Imaging, Kennedy Krieger Institute, Baltimore, Maryland, USA
| | - Peter C van Zijl
- F. M. Kirby Research Center for Functional Brain Imaging, Kennedy Krieger Institute, Baltimore, Maryland, USA
| |
Collapse
|
50
|
Wisnieff C, Liu T, Wang Y, Spincemaille P. The influence of molecular order and microstructure on the R2* and the magnetic susceptibility tensor. Magn Reson Imaging 2015; 34:682-9. [PMID: 26692502 DOI: 10.1016/j.mri.2015.12.003] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2015] [Revised: 12/03/2015] [Accepted: 12/07/2015] [Indexed: 01/07/2023]
Abstract
In this work, we demonstrate that in the presence of ordered sub-voxel structure such as tubular organization, biomaterials with molecular isotropy exhibits only apparent R2* anisotropy, while biomaterials with molecular anisotropy exhibit both apparent R2* and susceptibility anisotropy by means of susceptibility tensor imaging (STI). To this end, R2* and STI from gradient echo magnitude and phase data were examined in phantoms made from carbon fiber and Gadolinium (Gd) solutions with and without intrinsic molecular order and sub-voxel structure as well as in the in vivo brain. Confidence in the tensor reconstructions was evaluated with a wild bootstrap analysis. Carbon fiber showed both apparent anisotropy in R2* and anisotropy in STI, while the Gd filled capillary tubes only showed apparent anisotropy on R2*. Similarly, white matter showed anisotropic R2* and magnetic susceptibility with higher confidence, while the cerebral veins displayed only strong apparent R2* tensor anisotropy. Ordered sub-voxel tissue microstructure leads to apparent R2* anisotropy, which can be found in both white matter tracts and cerebral veins. However, additional molecular anisotropy is required for magnetic susceptibility anisotropy, which can be found in white matter tracts but not in cerebral veins.
Collapse
Affiliation(s)
- Cynthia Wisnieff
- Department of Biomedical Engineering, Cornell University, Ithaca, NY, USA
| | - Tian Liu
- Department of Radiology, Weill Cornell Medical College, New York, NY, USA; Medimagemetric, LLC, New York, NY, USA
| | - Yi Wang
- Department of Biomedical Engineering, Cornell University, Ithaca, NY, USA; Department of Radiology, Weill Cornell Medical College, New York, NY, USA
| | | |
Collapse
|