1
|
Plähn NMJ, Poli S, Peper ES, Açikgöz BC, Kreis R, Ganter C, Bastiaansen JAM. Getting the phase consistent: The importance of phase description in balanced steady-state free precession MRI of multi-compartment systems. Magn Reson Med 2024; 92:215-225. [PMID: 38321594 DOI: 10.1002/mrm.30033] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2023] [Revised: 01/09/2024] [Accepted: 01/12/2024] [Indexed: 02/08/2024]
Abstract
PURPOSE Determine the correct mathematical phase description for balanced steady-state free precession (bSSFP) signals in multi-compartment systems. THEORY AND METHODS Based on published bSSFP signal models, different phase descriptions can be formulated: one predicting the presence and the other predicting the absence of destructive interference effects in multi-compartment systems. Numerical simulations of bSSFP signals of water and acetone were performed to evaluate the predictions of these different phase descriptions. For experimental validation, bSSFP profiles were measured at 3T using phase-cycled bSSFP acquisitions performed in a phantom containing mixtures of water and acetone, which replicates a system with two signal components. Localized single voxel MRS was performed at 7T to determine the relative chemical shift of the acetone-water mixtures. RESULTS Based on the choice of phase description, the simulated bSSFP profiles of water-acetone mixtures varied significantly, either displaying or lacking destructive interference effects, as predicted theoretically. In phantom experiments, destructive interference was consistently observed in the measured bSSFP profiles of water-acetone mixtures, supporting the theoretical description that predicts such interference effects. The connection between the choice of phase description and predicted observation enables unambiguous experimental identification of the correct phase description for multi-compartment bSSFP profiles, which is consistent with the Bloch equations. CONCLUSION The study emphasizes that consistent phase descriptions are crucial for accurately describing multi-compartment bSSFP signals, as incorrect phase descriptions result in erroneous predictions.
Collapse
Affiliation(s)
- Nils M J Plähn
- Department of Diagnostic, Interventional and Pediatric Radiology (DIPR), Inselspital, Bern University Hospital, University of Bern, Bern, Switzerland
- Translation Imaging Center (TIC), Swiss Institute for Translational and Entrepreneurial Medicine, Bern, Switzerland
- Graduate School for Cellular and Biomedical Sciences, University of Bern, Bern, Switzerland
| | - Simone Poli
- Translation Imaging Center (TIC), Swiss Institute for Translational and Entrepreneurial Medicine, Bern, Switzerland
- Graduate School for Cellular and Biomedical Sciences, University of Bern, Bern, Switzerland
- MR Methodology, Department for Diagnostic and Interventional Neuroradiology, University of Bern, Bern, Switzerland
| | - Eva S Peper
- Department of Diagnostic, Interventional and Pediatric Radiology (DIPR), Inselspital, Bern University Hospital, University of Bern, Bern, Switzerland
- Translation Imaging Center (TIC), Swiss Institute for Translational and Entrepreneurial Medicine, Bern, Switzerland
| | - Berk C Açikgöz
- Department of Diagnostic, Interventional and Pediatric Radiology (DIPR), Inselspital, Bern University Hospital, University of Bern, Bern, Switzerland
- Translation Imaging Center (TIC), Swiss Institute for Translational and Entrepreneurial Medicine, Bern, Switzerland
- Graduate School for Cellular and Biomedical Sciences, University of Bern, Bern, Switzerland
| | - Roland Kreis
- Translation Imaging Center (TIC), Swiss Institute for Translational and Entrepreneurial Medicine, Bern, Switzerland
- MR Methodology, Department for Diagnostic and Interventional Neuroradiology, University of Bern, Bern, Switzerland
- Department for Biomedical Research, University of Bern, Bern, Switzerland
| | - Carl Ganter
- Department of Diagnostic and Interventional Radiology, TUM School of Medicine and Health, Klinikum rechts der Isar, Technical University of Munich, Munich, Germany
| | - Jessica A M Bastiaansen
- Department of Diagnostic, Interventional and Pediatric Radiology (DIPR), Inselspital, Bern University Hospital, University of Bern, Bern, Switzerland
- Translation Imaging Center (TIC), Swiss Institute for Translational and Entrepreneurial Medicine, Bern, Switzerland
| |
Collapse
|
2
|
Schäper J, Bieri O. Myelin water imaging at 0.55 T using a multigradient-echo sequence. Magn Reson Med 2024; 91:1043-1056. [PMID: 38010053 DOI: 10.1002/mrm.29949] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2023] [Revised: 10/19/2023] [Accepted: 11/12/2023] [Indexed: 11/29/2023]
Abstract
PURPOSE To investigate the prospects of a multigradient-echo (mGRE) acquisition for in vivo myelin water imaging at 0.55 T. METHODS Scans were performed on the brain of four healthy volunteers at 0.55 and 3 T, using a 3D mGRE sequence. The myelin water fraction (MWF) was calculated for both field strengths using a nonnegative least squares (NNLS) algorithm, implemented in the qMRLab suite. The quality of these maps as well as single-voxel fits were compared visually for 0.55 and 3 T. RESULTS The obtained MWF values at 0.55 T are consistent with previously reported ones at higher field strengths. The MWF maps are a considerable improvement over the ones at 3 T. Example fits show that 0.55 T data is better described by an exponential model than 3 T data, making the assumed multi-exponential model of the NNLS algorithm more accurate. CONCLUSION This first assessment shows that mGRE myelin water imaging at 0.55 T is feasible and has the potential to yield better results than at higher fields.
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
| | - Oliver Bieri
- Department of Biomedical Engineering, University of Basel, Basel, Switzerland
- Division of Radiological Physics, Department of Radiology, University Hospital Basel, Basel, Switzerland
| |
Collapse
|
3
|
Rossi GMC, Mackowiak ALC, Açikgöz BC, Pierzchała K, Kober T, Hilbert T, Bastiaansen JAM. SPARCQ: A new approach for fat fraction mapping using asymmetries in the phase-cycled balanced SSFP signal profile. Magn Reson Med 2023; 90:2348-2361. [PMID: 37496187 DOI: 10.1002/mrm.29813] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2023] [Revised: 06/19/2023] [Accepted: 07/12/2023] [Indexed: 07/28/2023]
Abstract
PURPOSE To develop SPARCQ (Signal Profile Asymmetries for Rapid Compartment Quantification), a novel approach to quantify fat fraction (FF) using asymmetries in the phase-cycled balanced SSFP (bSSFP) profile. METHODS SPARCQ uses phase-cycling to obtain bSSFP frequency profiles, which display asymmetries in the presence of fat and water at certain TRs. For each voxel, the measured signal profile is decomposed into a weighted sum of simulated profiles via multi-compartment dictionary matching. Each dictionary entry represents a single-compartment bSSFP profile with a specific off-resonance frequency and relaxation time ratio. Using the results of dictionary matching, the fractions of the different off-resonance components are extracted for each voxel, generating quantitative maps of water and FF and banding-artifact-free images for the entire image volume. SPARCQ was validated using simulations, experiments in a water-fat phantom and in knees of healthy volunteers. Experimental results were compared with reference proton density FFs obtained with 1 H-MRS (phantoms) and with multiecho gradient-echo MRI (phantoms and volunteers). SPARCQ repeatability was evaluated in six scan-rescan experiments. RESULTS Simulations showed that FF quantification is accurate and robust for SNRs greater than 20. Phantom experiments demonstrated good agreement between SPARCQ and gold standard FFs. In volunteers, banding-artifact-free quantitative maps and water-fat-separated images obtained with SPARCQ and ME-GRE demonstrated the expected contrast between fatty and non-fatty tissues. The coefficient of repeatability of SPARCQ FF was 0.0512. CONCLUSION SPARCQ demonstrates potential for fat quantification using asymmetries in bSSFP profiles and may be a promising alternative to conventional FF quantification techniques.
Collapse
Affiliation(s)
- Giulia M C Rossi
- Department of Diagnostic, Interventional and Pediatric Radiology (DIPR), Inselspital, Bern University Hospital, University of Bern, Bern, Switzerland
- Translational Imaging Center, Swiss Institute for Translational and Entrepreneurial Medicine, Bern, Switzerland
- Department of Diagnostic and Interventional Radiology, Lausanne University Hospital and University of Lausanne, Lausanne, Switzerland
| | - Adèle L C Mackowiak
- Department of Diagnostic, Interventional and Pediatric Radiology (DIPR), Inselspital, Bern University Hospital, University of Bern, Bern, Switzerland
- Translational Imaging Center, Swiss Institute for Translational and Entrepreneurial Medicine, Bern, Switzerland
- Department of Diagnostic and Interventional Radiology, Lausanne University Hospital and University of Lausanne, Lausanne, Switzerland
| | - Berk Can Açikgöz
- Department of Diagnostic, Interventional and Pediatric Radiology (DIPR), Inselspital, Bern University Hospital, University of Bern, Bern, Switzerland
- Translational Imaging Center, Swiss Institute for Translational and Entrepreneurial Medicine, Bern, Switzerland
| | - Katarzyna Pierzchała
- CIBM Center for Biomedical Imaging, Lausanne, Switzerland
- Animal Imaging and Technology, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
| | - Tobias Kober
- Department of Diagnostic and Interventional Radiology, Lausanne University Hospital and University of Lausanne, Lausanne, Switzerland
- Advanced Clinical Imaging Technology, Siemens Healthineers International AG, Lausanne, Switzerland
- LTS5, École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
| | - Tom Hilbert
- Department of Diagnostic and Interventional Radiology, Lausanne University Hospital and University of Lausanne, Lausanne, Switzerland
- Advanced Clinical Imaging Technology, Siemens Healthineers International AG, Lausanne, Switzerland
- LTS5, École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
| | - Jessica A M Bastiaansen
- Department of Diagnostic, Interventional and Pediatric Radiology (DIPR), Inselspital, Bern University Hospital, University of Bern, Bern, Switzerland
- Translational Imaging Center, Swiss Institute for Translational and Entrepreneurial Medicine, Bern, Switzerland
| |
Collapse
|
4
|
Birk F, Glang F, Loktyushin A, Birkl C, Ehses P, Scheffler K, Heule R. High-resolution neural network-driven mapping of multiple diffusion metrics leveraging asymmetries in the balanced steady-state free precession frequency profile. NMR IN BIOMEDICINE 2022; 35:e4669. [PMID: 34964998 DOI: 10.1002/nbm.4669] [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: 09/17/2021] [Revised: 12/04/2021] [Accepted: 12/06/2021] [Indexed: 06/14/2023]
Abstract
We propose to utilize the rich information content about microstructural tissue properties entangled in asymmetric balanced steady-state free precession (bSSFP) profiles to estimate multiple diffusion metrics simultaneously by neural network (NN) parameter quantification. A 12-point bSSFP phase-cycling scheme with high-resolution whole-brain coverage is employed at 3 and 9.4 T for NN input. Low-resolution target diffusion data are derived based on diffusion-weighted spin-echo echo-planar-imaging (SE-EPI) scans, that is, mean, axial, and radial diffusivity (MD, AD, and RD), fractional anisotropy (FA), as well as the spherical coordinates (azimuth Φ and inclination ϴ) of the principal diffusion eigenvector. A feedforward NN is trained with incorporated probabilistic uncertainty estimation. The NN predictions yielded highly reliable results in white matter (WM) and gray matter structures for MD. The quantification of FA, AD, and RD was overall in good agreement with the reference but the dependence of these parameters on WM anisotropy was somewhat biased (e.g. in corpus callosum). The inclination ϴ was well predicted for anisotropic WM structures, while the azimuth Φ was overall poorly predicted. The findings were highly consistent across both field strengths. Application of the optimized NN to high-resolution input data provided whole-brain maps with rich structural details. In conclusion, the proposed NN-driven approach showed potential to provide distortion-free high-resolution whole-brain maps of multiple diffusion metrics at high to ultrahigh field strengths in clinically relevant scan times.
Collapse
Affiliation(s)
- Florian Birk
- High Field Magnetic Resonance, Max Planck Institute for Biological Cybernetics, Tübingen, Germany
| | - Felix Glang
- High Field Magnetic Resonance, Max Planck Institute for Biological Cybernetics, Tübingen, Germany
| | - Alexander Loktyushin
- High Field Magnetic Resonance, Max Planck Institute for Biological Cybernetics, Tübingen, Germany
- Empirical Inference, Max Planck Institute for Intelligent Systems, Tübingen, Germany
| | - Christoph Birkl
- Department of Neuroradiology, Medical University of Innsbruck, Innsbruck, Austria
| | - Philipp Ehses
- German Center for Neurodegenerative Diseases (DZNE), Bonn, Germany
| | - Klaus Scheffler
- High Field Magnetic Resonance, Max Planck Institute for Biological Cybernetics, Tübingen, Germany
- Department of Biomedical Magnetic Resonance, University of Tübingen, Tübingen, Germany
| | - Rahel Heule
- High Field Magnetic Resonance, Max Planck Institute for Biological Cybernetics, Tübingen, Germany
- Department of Biomedical Magnetic Resonance, University of Tübingen, Tübingen, Germany
| |
Collapse
|
5
|
Keskin K, Yilmaz U, Cukur T. Constrained Ellipse Fitting for Efficient Parameter Mapping With Phase-Cycled bSSFP MRI. IEEE TRANSACTIONS ON MEDICAL IMAGING 2022; 41:14-26. [PMID: 34351856 DOI: 10.1109/tmi.2021.3102852] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Balanced steady-state free precession (bSSFP) imaging enables high scan efficiency in MRI, but differs from conventional sequences in terms of elevated sensitivity to main field inhomogeneity and nonstandard [Formula: see text]-weighted tissue contrast. To address these limitations, multiple bSSFP images of the same anatomy are commonly acquired with a set of different RF phase-cycling increments. Joint processing of phase-cycled acquisitions serves to mitigate sensitivity to field inhomogeneity. Recently phase-cycled bSSFP acquisitions were also leveraged to estimate relaxation parameters based on explicit signal models. While effective, these model-based methods often involve a large number of acquisitions (N ≈ 10-16), degrading scan efficiency. Here, we propose a new constrained ellipse fitting method (CELF) for parameter estimation with improved efficiency and accuracy in phase-cycled bSSFP MRI. CELF is based on the elliptical signal model framework for complex bSSFP signals; and it introduces geometrical constraints on ellipse properties to improve estimation efficiency, and dictionary-based identification to improve estimation accuracy. CELF generates maps of [Formula: see text], [Formula: see text], off-resonance and on-resonant bSSFP signal by employing a separate [Formula: see text] map to mitigate sensitivity to flip angle variations. Our results indicate that CELF can produce accurate off-resonance and banding-free bSSFP maps with as few as N = 4 acquisitions, while estimation accuracy for relaxation parameters is notably limited by biases from microstructural sensitivity of bSSFP imaging.
Collapse
|
6
|
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: 1.0] [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
|
7
|
From signal-based to comprehensive magnetic resonance imaging. Sci Rep 2021; 11:17216. [PMID: 34446804 PMCID: PMC8390767 DOI: 10.1038/s41598-021-96791-w] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2021] [Accepted: 08/09/2021] [Indexed: 12/03/2022] Open
Abstract
We present and evaluate a new insight into magnetic resonance imaging (MRI). It is based on the algebraic description of the magnetization during the transient response—including intrinsic magnetic resonance parameters such as longitudinal and transverse relaxation times (T1, T2) and proton density (PD) and experimental conditions such as radiofrequency field (B1) and constant/homogeneous magnetic field (B0) from associated scanners. We exploit the correspondence among three different elements: the signal evolution as a result of a repetitive sequence of blocks of radiofrequency excitation pulses and encoding gradients, the continuous Bloch equations and the mathematical description of a sequence as a linear system. This approach simultaneously provides, in a single measurement, all quantitative parameters of interest as well as associated system imperfections. Finally, we demonstrate the in-vivo applicability of the new concept on a clinical MRI scanner.
Collapse
|
8
|
Iyyakkunnel S, Schäper J, Bieri O. Configuration-based electrical properties tomography. Magn Reson Med 2020; 85:1855-1864. [PMID: 33107082 DOI: 10.1002/mrm.28542] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2020] [Revised: 09/09/2020] [Accepted: 09/14/2020] [Indexed: 11/07/2022]
Abstract
PURPOSE To introduce phase-based conductivity mapping from a configuration space analysis. METHODS The frequency response function of balanced SSFP (bSSFP) is used to perform a configuration space analysis. It is shown that the transceive phase for conductivity mapping can be directly obtained by a simple fast Fourier transform of a series of phase-cycled bSSFP scans. For validation, transceive phase and off-resonance mapping with fast Fourier transform is compared with phase estimation using a recently proposed method, termed PLANET. Experiments were performed in phantoms and for in vivo brain imaging at 3 T using a quadrature head coil. RESULTS For fast Fourier transform, aliasing can lead to systematic phase errors. This bias, however, decreases rapidly with increasing sampling points. Interestingly, Monte Carlo simulations revealed a lower uncertainty for the transceive phase and the off-resonance using fast Fourier transform as compared with PLANET. Both methods, however, essentially retrieve the same phase information from a set of phase-cycled bSSFP scans. As a result, configuration-based conductivity mapping was successfully performed using eight phase-cycled bSSFP scans in the phantoms and for brain tissues. Overall, the retrieved values were in good agreement with expectations. Conductivity estimation and mapping of the field inhomogeneities can therefore be performed in conjunction with the estimation of other quantitative parameters, such as relaxation, using configuration theory. CONCLUSIONS Phase-based conductivity mapping can be estimated directly from a simple Fourier analysis, such as in conjunction with relaxometry, using a series of phase-cycled bSSFP scans.
Collapse
Affiliation(s)
- Santhosh Iyyakkunnel
- Division of Radiological Physics, Department of Radiology, University Hospital Basel, Basel, Switzerland.,Department of Biomedical Engineering, University of Basel, Basel, Switzerland
| | - Jessica Schäper
- Division of Radiological Physics, Department of Radiology, University Hospital Basel, Basel, Switzerland.,Department of Biomedical Engineering, University of Basel, Basel, Switzerland
| | - Oliver Bieri
- Division of Radiological Physics, Department of Radiology, University Hospital Basel, Basel, Switzerland.,Department of Biomedical Engineering, University of Basel, Basel, Switzerland
| |
Collapse
|
9
|
Whitaker ST, Nataraj G, Nielsen JF, Fessler JA. Myelin water fraction estimation using small-tip fast recovery MRI. Magn Reson Med 2020; 84:1977-1990. [PMID: 32281179 PMCID: PMC7478173 DOI: 10.1002/mrm.28259] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2019] [Revised: 02/05/2020] [Accepted: 02/26/2020] [Indexed: 11/09/2022]
Abstract
PURPOSE To demonstrate the feasibility of an optimized set of small-tip fast recovery (STFR) MRI scans for rapidly estimating myelin water fraction (MWF) in the brain. METHODS We optimized a set of STFR scans to minimize the Cramér-Rao Lower Bound of MWF estimates. We evaluated the RMSE of MWF estimates from the optimized scans in simulation. We compared STFR-based MWF estimates (both modeling exchange and not modeling exchange) to multi-echo spin echo (MESE)-based estimates. We used the optimized scans to acquire in vivo data from which a MWF map was estimated. We computed the STFR-based MWF estimates using PERK, a recently developed kernel regression technique, and the MESE-based MWF estimates using both regularized non-negative least squares (NNLS) and PERK. RESULTS In simulation, the optimized STFR scans led to estimates of MWF with low RMSE across a range of tissue parameters and across white matter and gray matter. The STFR-based MWF estimates that modeled exchange compared well to MESE-based MWF estimates in simulation. When the optimized scans were tested in vivo, the MWF map that was estimated using a 3-compartment model with exchange was closer to the MESE-based MWF map. CONCLUSIONS The optimized STFR scans appear to be well suited for estimating MWF in simulation and in vivo when we model exchange in training. In this case, the STFR-based MWF estimates are close to the MESE-based estimates.
Collapse
Affiliation(s)
- Steven T. Whitaker
- Department of Electrical Engineering and Computer Science, University of Michigan, Ann Arbor, Michigan, USA
| | - Gopal Nataraj
- Department of Medical Physics, Memorial Sloan Kettering Cancer Center, New York, New York, USA
| | - Jon-Fredrik Nielsen
- Department of Biomedical Engineering, University of Michigan, Ann Arbor, Michigan, USA
| | - Jeffrey A. Fessler
- Department of Electrical Engineering and Computer Science, University of Michigan, Ann Arbor, Michigan, USA
| |
Collapse
|
10
|
Heule R, Bause J, Pusterla O, Scheffler K. Multi-parametric artificial neural network fitting of phase-cycled balanced steady-state free precession data. Magn Reson Med 2020; 84:2981-2993. [PMID: 32479661 DOI: 10.1002/mrm.28325] [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: 01/17/2020] [Revised: 04/22/2020] [Accepted: 04/27/2020] [Indexed: 12/16/2022]
Abstract
PURPOSE Standard relaxation time quantification using phase-cycled balanced steady-state free precession (bSSFP), eg, motion-insensitive rapid configuration relaxometry (MIRACLE), is subject to a considerable underestimation of tissue T1 and T2 due to asymmetric intra-voxel frequency distributions. In this work, an artificial neural network (ANN) fitting approach is proposed to simultaneously extract accurate reference relaxation times (T1 , T2 ) and robust field map estimates ( B 1 + , ΔB0 ) from the bSSFP profile. METHODS Whole-brain bSSFP data acquired at 3T were used for the training of a feedforward ANN with N = 12, 6, and 4 phase-cycles. The magnitude and phase of the Fourier transformed complex bSSFP frequency response served as input and the multi-parametric reference set [T1 , T2 , B 1 + , ∆B0 ] as target. The ANN predicted relaxation times were validated against the target and MIRACLE. RESULTS The ANN prediction of T1 and T2 for trained and untrained data agreed well with the reference, even for only four acquired phase-cycles. In contrast, relaxometry based on 4-point MIRACLE was prone to severe off-resonance-related artifacts. ANN predicted B 1 + and ∆B0 maps showed the expected spatial inhomogeneity patterns in high agreement with the reference measurements for 12-point, 6-point, and 4-point bSSFP phase-cycling schemes. CONCLUSION ANNs show promise to provide accurate brain tissue T1 and T2 values as well as reliable field map estimates. Moreover, the bSSFP acquisition can be accelerated by reducing the number of phase-cycles while still delivering robust T1 , T2 , B 1 + , and ∆B0 estimates.
Collapse
Affiliation(s)
- Rahel Heule
- High Field Magnetic Resonance, Max Planck Institute for Biological Cybernetics, Tübingen, Germany
| | - Jonas Bause
- High Field Magnetic Resonance, Max Planck Institute for Biological Cybernetics, Tübingen, Germany
| | - Orso Pusterla
- Institute for Biomedical Engineering, University and ETH Zurich, Zurich, Switzerland.,Division of Radiological Physics, Department of Radiology, University Hospital Basel, Basel, Switzerland.,Department of Biomedical Engineering, University of Basel, Basel, Switzerland
| | - Klaus Scheffler
- High Field Magnetic Resonance, Max Planck Institute for Biological Cybernetics, Tübingen, Germany.,Department of Biomedical Magnetic Resonance, University of Tübingen, Tübingen, Germany
| |
Collapse
|
11
|
Heule R, Deshmane A, Zaiss M, Herz K, Ehses P, Scheffler K. Structure or Exchange? On the Feasibility of Chemical Exchange Detection with Balanced Steady-State Free Precession in Tissue - An In Vitro Study. NMR IN BIOMEDICINE 2020; 33:e4200. [PMID: 31833130 DOI: 10.1002/nbm.4200] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/02/2019] [Revised: 08/15/2019] [Accepted: 09/16/2019] [Indexed: 06/10/2023]
Abstract
Balanced steady-state free precession imaging has recently been suggested for chemical exchange detection (bSSFPX). The objective of this work is to investigate the contributions of microstructural, chemical shift and chemical exchange effects to the asymmetry of the bSSFP profile at field strengths of 3 T and 9.4 T. To this end, in vitro bSSFPX experiments are performed for a range of repetition times and flip angles in glucose water solutions with different MnCl2 concentrations and tissue homogenates obtained from the brainstem of pig brains. The experimental results are compared to multi-pool Bloch-McConnell simulations. Additionally, the influence of white matter tract geometry is analyzed ex vivo in pig brain hemispheres measured at two different angles with respect to B0 . The detectable bSSFP profile asymmetry in glucose solutions with tissue-like relaxation times and white matter homogenates was consistent with Bloch-McConnell simulations but relatively small. In intact white matter tracts, the asymmetry was dominated by structural effects with a strong dependency on tract orientation relative to B0 . In tracts perpendicular to B0 , the asymmetry was ≈ 3-4 times higher than in the homogenates, thus barely affected by chemical exchange effects. In conclusion, chemical exchange-related bSSFP profile asymmetries are detectable in tissue homogenates, however, the observed asymmetry level is generally low and prone to confounding structural effects rendering in vivo chemical exchange detection with bSSFP challenging in the brain.
Collapse
Affiliation(s)
- Rahel Heule
- High Field Magnetic Resonance, Max Planck Institute for Biological Cybernetics, Tübingen, Germany
| | - Anagha Deshmane
- High Field Magnetic Resonance, Max Planck Institute for Biological Cybernetics, Tübingen, Germany
| | - Moritz Zaiss
- High Field Magnetic Resonance, Max Planck Institute for Biological Cybernetics, Tübingen, Germany
| | - Kai Herz
- High Field Magnetic Resonance, Max Planck Institute for Biological Cybernetics, Tübingen, Germany
| | - Philipp Ehses
- German Center for Neurodegenerative Diseases (DZNE), Bonn, Germany
| | - Klaus Scheffler
- High Field Magnetic Resonance, Max Planck Institute for Biological Cybernetics, Tübingen, Germany
- Department of Biomedical Magnetic Resonance, University of Tübingen, Tübingen, Germany
| |
Collapse
|
12
|
Wood TC, Teixeira RPAG, Malik SJ. Magnetization transfer and frequency distribution effects in the SSFP ellipse. Magn Reson Med 2019; 84:857-865. [PMID: 31872921 PMCID: PMC7216875 DOI: 10.1002/mrm.28149] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2019] [Revised: 11/15/2019] [Accepted: 12/06/2019] [Indexed: 01/08/2023]
Abstract
Purpose To demonstrate that quantitative magnetization transfer (qMT) parameters can be extracted from steady‐state free‐precession (SSFP) data with no external T1 map or banding artifacts. Methods SSFP images with multiple MT weightings were acquired and qMT parameters fitted with a two‐stage elliptical signal model. Results Monte Carlo simulations and data from a 3T scanner indicated that most qMT parameters could be recovered with reasonable accuracy. Systematic deviations from theory were observed in white matter, consistent with previous literature on frequency distribution effects. Conclusions qMT parameters can be extracted from SSFP data alone, in a manner robust to banding artifacts, despite several confounds.
Collapse
Affiliation(s)
- Tobias C Wood
- Department of Neuroimaging, King's College London, London, UK
| | - Rui P A G Teixeira
- School of Imaging Sciences and Biomedical Engineering, King's College London, London, UK
| | - Shaihan J Malik
- School of Imaging Sciences and Biomedical Engineering, King's College London, London, UK
| |
Collapse
|
13
|
Gavazzi S, Shcherbakova Y, Bartels LW, Stalpers LJA, Lagendijk JJW, Crezee H, van den Berg CAT, van Lier ALHMW. Transceive phase mapping using the PLANET method and its application for conductivity mapping in the brain. Magn Reson Med 2019; 83:590-607. [PMID: 31483520 PMCID: PMC6900152 DOI: 10.1002/mrm.27958] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2019] [Revised: 07/25/2019] [Accepted: 07/30/2019] [Indexed: 12/23/2022]
Abstract
Purpose To demonstrate feasibility of transceive phase mapping with the PLANET method and its application for conductivity reconstruction in the brain. Methods Accuracy and precision of transceive phase (ϕ±) estimation with PLANET, an ellipse fitting approach to phase‐cycled balanced steady state free precession (bSSFP) data, were assessed with simulations and measurements and compared to standard bSSFP. Measurements were conducted on a homogeneous phantom and in the brain of healthy volunteers at 3 tesla. Conductivity maps were reconstructed with Helmholtz‐based electrical properties tomography. In measurements, PLANET was also compared to a reference technique for transceive phase mapping, i.e., spin echo. Results Accuracy and precision of ϕ± estimated with PLANET depended on the chosen flip angle and TR. PLANET‐based ϕ± was less sensitive to perturbations induced by off‐resonance effects and partial volume (e.g., white matter + myelin) than bSSFP‐based ϕ±. For flip angle = 25° and TR = 4.6 ms, PLANET showed an accuracy comparable to that of reference spin echo but a higher precision than bSSFP and spin echo (factor of 2 and 3, respectively). The acquisition time for PLANET was ~5 min; 2 min faster than spin echo and 8 times slower than bSSFP. However, PLANET simultaneously reconstructed T1, T2, B0 maps besides mapping ϕ±. In the phantom, PLANET‐based conductivity matched the true value and had the smallest spread of the three methods. In vivo, PLANET‐based conductivity was similar to spin echo‐based conductivity. Conclusion Provided that appropriate sequence parameters are used, PLANET delivers accurate and precise ϕ± maps, which can be used to reconstruct brain tissue conductivity while simultaneously recovering T1, T2, and B0 maps.
Collapse
Affiliation(s)
- Soraya Gavazzi
- Department of Radiotherapy, University Medical Center Utrecht, Utrecht, The Netherlands
| | - Yulia Shcherbakova
- Center for Image Sciences, University Medical Center Utrecht, Utrecht, The Netherlands
| | - Lambertus W Bartels
- Center for Image Sciences, University Medical Center Utrecht, Utrecht, The Netherlands.,Image Sciences Institute, Imaging Division, University Medical Center Utrecht, Utrecht, The Netherlands
| | - Lukas J A Stalpers
- Department of Radiotherapy, Amsterdam UMC, University of Amsterdam, Amsterdam, The Netherlands
| | - Jan J W Lagendijk
- Department of Radiotherapy, University Medical Center Utrecht, Utrecht, The Netherlands
| | - Hans Crezee
- Department of Radiotherapy, Amsterdam UMC, University of Amsterdam, Amsterdam, The Netherlands
| | - Cornelis A T van den Berg
- Department of Radiotherapy, University Medical Center Utrecht, Utrecht, The Netherlands.,Center for Image Sciences, University Medical Center Utrecht, Utrecht, The Netherlands
| | | |
Collapse
|
14
|
Reynaud O, da Silva AR, Gruetter R, Jelescu IO. Multi-slice passband bSSFP for human and rodent fMRI at ultra-high field. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2019; 305:31-40. [PMID: 31195214 DOI: 10.1016/j.jmr.2019.05.010] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/07/2019] [Revised: 05/07/2019] [Accepted: 05/28/2019] [Indexed: 06/09/2023]
Abstract
Balanced steady-state free precession (bSSFP) can be used as an alternative to gradient-echo (GE) EPI for BOLD functional MRI when image distortions and signal drop-outs are severe such as at ultra-high field. However, 3D-bSSFP acquisitions have distinct drawbacks on either human or animal MR systems. On clinical scanners, 3D imaging is suboptimal for localized fMRI applications. It can also display distortions when acceleration methods such as spiral read-outs are used, and, compared to multi-slice acquisitions, suffers from increased sensitivity to motion or physiological noise which further results in blurring. On pre-clinical systems, 3D acquisitions have low temporal resolution due to limited acceleration options, while single slice often results in insufficient coverage. The aim of the present study was to implement a multi-slice bSSFP acquisition with Cartesian read-out to obtain non-distorted BOLD fMRI activation maps in the human and rat brain at ultra-high field. We show that, when using a new pseudo-steady-state, the bSSFP signal characteristics are preserved. In the human brain at 7 T, we demonstrate that both task- and resting-state fMRI can be performed with multi-slice bSSFP, with a temporal SNR that matches that of 3D-bSSFP, resulting in - at least - equal performance. In the rat brain at 14 T, we show that the multi-slice bSSFP protocol has similar sensitivity to gradient-echo EPI for task fMRI, while benefitting from much reduced distortions and drop-outs. The advantages of passband bSSFP at 14 T in comparison with GE-EPI are expected to be even more marked for mouse brain.
Collapse
Affiliation(s)
- Olivier Reynaud
- Centre d'Imagerie Biomédicale, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
| | - Analina R da Silva
- Centre d'Imagerie Biomédicale, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
| | - Rolf Gruetter
- Centre d'Imagerie Biomédicale, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
| | - Ileana O Jelescu
- Centre d'Imagerie Biomédicale, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland.
| |
Collapse
|
15
|
Shcherbakova Y, van den Berg CAT, Moonen CTW, Bartels LW. Investigation of the influence of B 0 drift on the performance of the PLANET method and an algorithm for drift correction. Magn Reson Med 2019; 82:1725-1740. [PMID: 31317584 PMCID: PMC6772029 DOI: 10.1002/mrm.27860] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2019] [Revised: 04/26/2019] [Accepted: 05/23/2019] [Indexed: 01/03/2023]
Abstract
Purpose The PLANET method was designed to simultaneously reconstruct maps of T1 and T2, the off‐resonance, the RF phase, and the banding free signal magnitude. The method requires a stationary B0 field over the course of a phase‐cycled balanced SSFP acquisition. In this work we investigated the influence of B0 drift on the performance of the PLANET method for single‐component and two‐component signal models, and we propose a strategy for drift correction. Methods The complex phase‐cycled balanced SSFP signal was modeled with and without frequency drift. The behavior of the signal influenced by drift was mathematically interpreted as a sum of drift‐dependent displacement of the data points along an ellipse and drift‐dependent rotation around the origin. The influence of drift on parameter estimates was investigated experimentally on a phantom and on the brain of healthy volunteers and was verified by numerical simulations. A drift correction algorithm was proposed and tested on a phantom and in vivo. Results Drift can be assumed to be linear over the typical duration of a PLANET acquisition. In a phantom (a single‐component signal model), drift induced errors of 4% and 8% in the estimated T1 and T2 values. In the brain, where multiple components are present, drift only had a minor effect. For both single‐component and two‐component signal models, drift‐induced errors were successfully corrected by applying the proposed drift correction algorithm. Conclusion We have demonstrated theoretically and experimentally the sensitivity of the PLANET method to B0 drift and have proposed a drift correction method.
Collapse
Affiliation(s)
- Yulia Shcherbakova
- Center for Image Sciences, Imaging Division, University Medical Center Utrecht, Utrecht, Netherlands
| | - Cornelis A T van den Berg
- Department of Radiotherapy, Imaging Division, University Medical Center Utrecht, Utrecht, Netherlands
| | - Chrit T W Moonen
- Center for Image Sciences, Imaging Division, University Medical Center Utrecht, Utrecht, Netherlands
| | - Lambertus W Bartels
- Center for Image Sciences, Imaging Division, University Medical Center Utrecht, Utrecht, Netherlands
| |
Collapse
|
16
|
Wang CY, Coppo S, Mehta BB, Seiberlich N, Yu X, Griswold MA. Magnetic resonance fingerprinting with quadratic RF phase for measurement of T 2 * simultaneously with δ f , T 1 , and T 2. Magn Reson Med 2018; 81:1849-1862. [PMID: 30499221 DOI: 10.1002/mrm.27543] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2018] [Revised: 08/29/2018] [Accepted: 08/30/2018] [Indexed: 11/07/2022]
Abstract
PURPOSE This study explores the possibility of using a gradient moment balanced sequence with a quadratically varied RF excitation phase in the magnetic resonance fingerprinting (MRF) framework to quantify T2 * in addition to δ f , T1 , and T2 tissue properties. METHODS The proposed quadratic RF phase-based MRF method (qRF-MRF) combined a varied RF excitation phase with the existing balanced SSFP (bSSFP)-based MRF method to generate signals that were uniquely sensitive to δ f , T1 , T2 , as well as the distribution width of intravoxel frequency dispersion, Γ . A dictionary, generated through Bloch simulation, containing possible signal evolutions within the physiological range of δ f , T1 , T2 , and Γ , was used to perform parameter estimation. The estimated T2 and Γ were subsequently used to estimate T2 * . The proposed method was evaluated in phantom experiments and healthy volunteers (N = 5). RESULTS The T1 and T2 values from the phantom by qRF-MRF demonstrated good agreement with values obtained by traditional gold standard methods (r2 = 0.995 and 0.997, respectively; concordance correlation coefficient = 0.978 and 0.995, respectively). The T2 * values from the phantom demonstrated good agreement with values obtained through the multi-echo gradient-echo method (r2 = 0.972, concordance correlation coefficient = 0.983). In vivo qRF-MRF-measured T1 , T2 , and T2 * values were compared with measurements by existing methods and literature values. CONCLUSION The proposed qRF-MRF method demonstrated the potential for simultaneous quantification of δ f , T1 , T2 , and T2 * tissue properties.
Collapse
Affiliation(s)
- Charlie Yi Wang
- Department of Biomedical Engineering, Case Western Reserve University, Cleveland, Ohio
| | - Simone Coppo
- Department of Radiology, Case Western Reserve University, Cleveland, Ohio
| | | | - Nicole Seiberlich
- Department of Biomedical Engineering, Case Western Reserve University, Cleveland, Ohio.,Department of Radiology, Case Western Reserve University, Cleveland, Ohio
| | - Xin Yu
- Department of Biomedical Engineering, Case Western Reserve University, Cleveland, Ohio.,Department of Radiology, Case Western Reserve University, Cleveland, Ohio
| | - Mark Alan Griswold
- Department of Biomedical Engineering, Case Western Reserve University, Cleveland, Ohio.,Department of Radiology, Case Western Reserve University, Cleveland, Ohio
| |
Collapse
|
17
|
Shcherbakova Y, van den Berg CAT, Moonen CTW, Bartels LW. On the accuracy and precision of PLANET for multiparametric MRI using phase-cycled bSSFP imaging. Magn Reson Med 2018; 81:1534-1552. [PMID: 30303562 PMCID: PMC6585657 DOI: 10.1002/mrm.27491] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2018] [Revised: 07/17/2018] [Accepted: 07/19/2018] [Indexed: 11/09/2022]
Abstract
Purpose In this work we demonstrate how sequence parameter settings influence the accuracy and precision in T1, T2, and off‐resonance maps obtained with the PLANET method for a single‐component signal model. In addition, the performance of the method for the particular case of a two‐component relaxation model for white matter tissue was assessed. Methods Numerical simulations were performed to investigate the influence of sequence parameter settings on the accuracy and precision in the estimated parameters for a single‐component model, as well as for a two‐component white matter model. Phantom and in vivo experiments were performed for validation. In addition, the effects of Gibbs ringing were investigated. Results By making a proper choice for sequence parameter settings, accurate and precise parameter estimation can be achieved for a single‐component signal model over a wide range of relaxation times at realistic SNR levels. Due to the presence of a second myelin‐related signal component in white matter, an underestimation of approximately 30% in T1 and T2 was observed, predicted by simulations and confirmed by measurements. Gibbs ringing artifacts correction improved the precision and accuracy of the parameter estimates. Conclusion For a single‐component signal model there is a broad “sweet spot” of sequence parameter combinations for which a high accuracy and precision in the parameter estimates is achieved over a wide range of relaxation times. For a multicomponent signal model, the single‐component PLANET reconstruction results in systematic errors in the parameter estimates as expected.
Collapse
Affiliation(s)
- Yulia Shcherbakova
- Center for Image Sciences, Imaging Division, University Medical Center Utrecht, Utrecht, the Netherlands
| | - Cornelis A T van den Berg
- 2Department of Radiotherapy, Imaging Division, University Medical Center Utrecht, Utrecht, the Netherlands
| | - Chrit T W Moonen
- Center for Image Sciences, Imaging Division, University Medical Center Utrecht, Utrecht, the Netherlands
| | - Lambertus W Bartels
- Center for Image Sciences, Imaging Division, University Medical Center Utrecht, Utrecht, the Netherlands
| |
Collapse
|
18
|
Malik SJ, Teixeira RPAG, Hajnal JV. Extended phase graph formalism for systems with magnetization transfer and exchange. Magn Reson Med 2017; 80:767-779. [PMID: 29243295 PMCID: PMC5947218 DOI: 10.1002/mrm.27040] [Citation(s) in RCA: 54] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2017] [Revised: 11/02/2017] [Accepted: 11/19/2017] [Indexed: 01/23/2023]
Abstract
Purpose An extended phase graph framework (EPG‐X) for modeling systems with exchange or magnetization transfer (MT) is proposed. Theory EPG‐X models coupled two‐compartment systems by describing each compartment with separate phase graphs that exchange during evolution periods. There are two variants: EPG‐X(BM) for systems governed by the Bloch‐McConnell equations, and EPG‐X(MT) for the pulsed MT formalism. For the MT case, the “bound” protons have no transverse components, so their phase graph consists of only longitudinal states. Methods The EPG‐X model was validated against steady‐state solutions and isochromat‐based simulation of gradient‐echo sequences. Three additional test cases were investigated: (i) MT effects in multislice turbo spin‐echo; (ii) variable flip angle gradient‐echo imaging of the type used for MR fingerprinting; and (iii) water exchange in multi‐echo spin‐echo T2 relaxometry. Results EPG‐X was validated successfully against isochromat based transient simulations and known steady‐state solutions. EPG‐X(MT) simulations matched in‐vivo measurements of signal attenuation in white matter in multislice turbo spin‐echo images. Magnetic resonance fingerprinting–style experiments with a bovine serum albumin (MT) phantom showed that the data were not consistent with a single‐pool model, but EPG‐X(MT) could be used to fit the data well. The EPG‐X(BM) simulations of multi‐echo spin‐echo T2 relaxometry suggest that exchange could lead to an underestimation of the myelin‐water fraction. Conclusions The EPG‐X framework can be used for modeling both steady‐state and transient signal response of systems exhibiting exchange or MT. This may be particularly beneficial for relaxometry approaches that rely on characterizing transient rather than steady‐state sequences. Magn Reson Med 80:767–779, 2018. © 2017 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.
Collapse
Affiliation(s)
- Shaihan J Malik
- Department of Biomedical Engineering, School of Biomedical Engineering and Imaging Sciences, King's College London, St. Thomas' Hospital, London, SE1 7EH, United Kingdom.,Centre for the Developing Brain, School of Biomedical Engineering and Imaging Sciences, King's College London, St. Thomas' Hospital, London, SE1 7EH, United Kingdom
| | - Rui Pedro A G Teixeira
- Department of Biomedical Engineering, School of Biomedical Engineering and Imaging Sciences, King's College London, St. Thomas' Hospital, London, SE1 7EH, United Kingdom.,Centre for the Developing Brain, School of Biomedical Engineering and Imaging Sciences, King's College London, St. Thomas' Hospital, London, SE1 7EH, United Kingdom
| | - Joseph V Hajnal
- Department of Biomedical Engineering, School of Biomedical Engineering and Imaging Sciences, King's College London, St. Thomas' Hospital, London, SE1 7EH, United Kingdom.,Centre for the Developing Brain, School of Biomedical Engineering and Imaging Sciences, King's College London, St. Thomas' Hospital, London, SE1 7EH, United Kingdom
| |
Collapse
|
19
|
Shcherbakova Y, van den Berg CAT, Moonen CTW, Bartels LW. PLANET: An ellipse fitting approach for simultaneous T 1 and T 2 mapping using phase-cycled balanced steady-state free precession. Magn Reson Med 2017; 79:711-722. [PMID: 28543430 PMCID: PMC5811804 DOI: 10.1002/mrm.26717] [Citation(s) in RCA: 37] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2016] [Revised: 03/23/2017] [Accepted: 03/26/2017] [Indexed: 11/09/2022]
Abstract
Purpose To demonstrate the feasibility of a novel, ellipse fitting approach, named PLANET, for simultaneous estimation of relaxation times T1 and T2 from a single 3D phase‐cycled balanced steady‐state free precession (bSSFP) sequence. Methods A method is presented in which the elliptical signal model is used to describe the phase‐cycled bSSFP steady‐state signal. The fitting of the model to the acquired data is reformulated into a linear convex problem, which is solved directly by a linear least squares method, specific to ellipses. Subsequently, the relaxation times T1 and T2, the banding free magnitude, and the off‐resonance are calculated from the fitting results. Results Maps of T1 and T2, as well as an off‐resonance and a banding free magnitude can be simultaneously, quickly, and robustly estimated from a single 3D phase‐cycled bSSFP sequence. The feasibility of the method was demonstrated in a phantom and in the brain of healthy volunteers on a clinical MR scanner. The results were in good agreement for the phantom, but a systematic underestimation of T1 was observed in the brain. Conclusion The presented method allows for accurate mapping of relaxation times and off‐resonance, and for the reconstruction of banding free magnitude images at realistic signal‐to‐noise ratios. Magn Reson Med 79:711–722, 2018. © 2017 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution‐NonCommercial‐NoDerivs License, which permits use and distribution in any medium, provided the original work is properly cited, the use is non‐commercial and no modifications or adaptations are made.
Collapse
Affiliation(s)
- Yulia Shcherbakova
- Center for Image Sciences, Imaging Division, University Medical Center Utrecht, Utrecht, The Netherlands
| | - Cornelis A T van den Berg
- Department of Radiotherapy, Imaging Division, University Medical Center Utrecht, Utrecht, The Netherlands
| | - Chrit T W Moonen
- Center for Image Sciences, Imaging Division, University Medical Center Utrecht, Utrecht, The Netherlands
| | - Lambertus W Bartels
- Center for Image Sciences, Imaging Division, University Medical Center Utrecht, Utrecht, The Netherlands
| |
Collapse
|
20
|
Nguyen D, Bieri O. Motion-insensitive rapid configuration relaxometry. Magn Reson Med 2016; 78:518-526. [DOI: 10.1002/mrm.26384] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2016] [Revised: 07/19/2016] [Accepted: 07/26/2016] [Indexed: 11/09/2022]
Affiliation(s)
- Damien Nguyen
- Division of Radiological Physics; Department of Radiology, University of Basel Hospital; Basel Switzerland
- Department of Biomedical Engineering; University of Basel; Basel Switzerland
| | - Oliver Bieri
- Division of Radiological Physics; Department of Radiology, University of Basel Hospital; Basel Switzerland
- Department of Biomedical Engineering; University of Basel; Basel Switzerland
| |
Collapse
|
21
|
Oh SH, Bilello M, Schindler M, Markowitz CE, Detre JA, Lee J. Direct visualization of short transverse relaxation time component (ViSTa). Neuroimage 2013; 83:485-92. [PMID: 23796545 DOI: 10.1016/j.neuroimage.2013.06.047] [Citation(s) in RCA: 66] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2013] [Revised: 05/14/2013] [Accepted: 06/13/2013] [Indexed: 10/26/2022] Open
Abstract
White matter of the brain has been demonstrated to have multiple relaxation components. Among them, the short transverse relaxation time component (T2<40 ms; T2⁎<25 ms at 3 T) has been suggested to originate from myelin water whereas long transverse relaxation time components have been associated with axonal and/or interstitial water. In myelin water imaging, T2 or T2⁎ signal decay is measured to estimate myelin water fraction based on T2 or T2⁎ differences among the water components. This method has been demonstrated to be sensitive to demyelination in the brain but suffers from low SNR and image artifacts originating from ill-conditioned multi-exponential fitting. In this study, a novel approach that selectively acquires short transverse relaxation time signal is proposed. The method utilizes a double inversion RF pair to suppress a range of long T1 signal. This suppression leaves short T2⁎ signal, which has been suggested to have short T1, as the primary source of the image. The experimental results confirm that after suppression of long T1 signals, the image is dominated by short T2⁎ in the range of myelin water, allowing us to directly visualize the short transverse relaxation time component in the brain. Compared to conventional myelin water imaging, this new method of direct visualization of short relaxation time component (ViSTa) provides high quality images. When applied to multiple sclerosis patients, chronic lesions show significantly reduced signal intensity in ViSTa images suggesting sensitivity to demyelination.
Collapse
Affiliation(s)
- Se-Hong Oh
- Department of Radiology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | | | | | | | | | | |
Collapse
|
22
|
Sati P, van Gelderen P, Silva AC, Reich DS, Merkle H, de Zwart JA, Duyn JH. Micro-compartment specific T2* relaxation in the brain. Neuroimage 2013; 77:268-78. [PMID: 23528924 DOI: 10.1016/j.neuroimage.2013.03.005] [Citation(s) in RCA: 158] [Impact Index Per Article: 14.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2012] [Revised: 01/29/2013] [Accepted: 03/05/2013] [Indexed: 12/12/2022] Open
Abstract
MRI at high field can be sensitized to the magnetic properties of tissues, which introduces a signal dependence on the orientation of white matter (WM) fiber bundles relative to the magnetic field. In addition, study of the NMR relaxation properties of this signal has indicated contributions from compartmentalized water environments inside and outside the myelin sheath that may be separable. Here we further investigated the effects of water compartmentalization on the MRI signal with the goal of extracting compartment-specific information. By comparing MRI measurements of human and marmoset brain at 7T with magnetic field modeling, we show that: (1) water between the myelin lipid bilayers, in the axonal, and in the interstitial space each experience characteristic magnetic field effects that depend on fiber orientation (2) these field effects result in characteristic relaxation properties and frequency shifts for these compartments; and (3) compartmental contributions may be separated by multi-component fitting of the MRI signal relaxation (i.e. decay) curve. We further show the potential application of these findings to the direct mapping of myelin content and assessment of WM fiber integrity with high field MRI.
Collapse
Affiliation(s)
- Pascal Sati
- Translational Neuroradiology Unit, Neuroimmunology Branch, National Institute of Neurological Disorders and Stroke, National Institutes of Health, Bethesda, MD 20892, USA
| | | | | | | | | | | | | |
Collapse
|
23
|
Chen WC, Foxley S, Miller KL. Detecting microstructural properties of white matter based on compartmentalization of magnetic susceptibility. Neuroimage 2012; 70:1-9. [PMID: 23268785 PMCID: PMC3591254 DOI: 10.1016/j.neuroimage.2012.12.032] [Citation(s) in RCA: 61] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/25/2012] [Revised: 11/18/2012] [Accepted: 12/18/2012] [Indexed: 11/24/2022] Open
Abstract
The microscopic structure of neuronal tissue is crucial to brain function, with axon diameter, axonal density and myelination directly influencing signal conduction in the white matter. There is increasing evidence that these microstructural properties alter signal in magnetic resonance imaging (MRI) driven by magnetic susceptibility of different compartments (e.g., myelin sheaths and iron-laden cells). To explain these observations, we have developed a multi-compartmental geometric model of whitematter microstructure. Using a single set of literature parameters, this forward model predicts experimentally observed orientation dependence and temporal evolution of the MRI signal. Where previous models have aimed to explain only the orientation dependence of signal phase, the proposed approach encapsulates the full repertoire of signal behavior. The frequency distribution underlying signal behavior is predicted to be a rich source of microstructural information with relevance to neuronal pathology.
Collapse
Affiliation(s)
- Way Cherng Chen
- Centre for Functional MRI of the Brain, University of Oxford, Oxford, UK.
| | | | | |
Collapse
|
24
|
Challenges of High-resolution Diffusion Imaging of the Human Medial Temporal Lobe in Alzheimer Disease. Top Magn Reson Imaging 2012; 21:355-65. [PMID: 22158129 DOI: 10.1097/rmr.0b013e31823f6413] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
The human medial temporal lobe performs an essential role in memory formation and retrieval. Diseases involving the hippocampus such as Alzheimer disease present a unique opportunity for advanced imaging techniques to detect abnormalities at an early stage. In particular, it is possible that diffusion imaging will measure abnormal microarchitecture beyond the realm of macroscopic imaging. However, this task is formidable because of the detailed anatomy of the medial temporal lobe, the difficulties in obtaining high-quality diffusion images of adequate resolution, and the challenges in diffusion data processing. Moreover, it is unclear if any differences will be significant for an individual patient or simply groups of patients. Successful endeavors will need to address each of these challenges in an integrated fashion. The rewards of such analysis may be detection of microscopic disease in vivo, which could represent a landmark accomplishment for the field of neuroradiology.
Collapse
|
25
|
van Gelderen P, de Zwart JA, Lee J, Sati P, Reich DS, Duyn JH. Nonexponential T₂ decay in white matter. Magn Reson Med 2011; 67:110-7. [PMID: 21630352 DOI: 10.1002/mrm.22990] [Citation(s) in RCA: 83] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2011] [Revised: 04/08/2011] [Accepted: 04/12/2011] [Indexed: 11/05/2022]
Abstract
Visualizing myelin in human brain may help the study of diseases such as multiple sclerosis. Previous studies based on T(1) and T(2) relaxation contrast have suggested the presence of a distinct water pool that may report directly on local myelin content. Recent work indicates that T(2) contrast may offer particular advantages over T(1) and T(2) contrast, especially at high field. However, the complex mechanism underlying T(2) relaxation may render interpretation difficult. To address this issue, T(2) relaxation behavior in human brain was studied at 3 and 7 T. Multiple gradient echoes covering most of the decay curve were analyzed for deviations from mono-exponential behavior. The data confirm the previous finding of a distinct rapidly relaxing signal component (T(2) ∼ 6 ms), tentatively attributed to myelin water. However, in extension to previous findings, this rapidly relaxing component displayed a substantial resonance frequency shift, reaching 36 Hz in the corpus callosum at 7 T. The component's fractional amplitude and frequency shift appeared to depend on both field strength and fiber orientation, consistent with a mechanism originating from magnetic susceptibility effects. The findings suggest that T(2) contrast at high field may be uniquely sensitive to tissue myelin content and that proper interpretation will require modeling of susceptibility-induced resonance frequency shifts.
Collapse
Affiliation(s)
- Peter van Gelderen
- Advanced MRI Section, Laboratory of Functional and Molecular Imaging, National Institute of Neurological Disorders and Stroke, National Institutes of Health, Bethesda, Maryland, USA.
| | | | | | | | | | | |
Collapse
|
26
|
Miller KL, Tijssen RHN, Stikov N, Okell TW. Steady-state MRI: methods for neuroimaging. ACTA ACUST UNITED AC 2011. [DOI: 10.2217/iim.10.66] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|