1
|
Herthum H, Hetzer S, Kreft B, Tzschätzsch H, Shahryari M, Meyer T, Görner S, Neubauer H, Guo J, Braun J, Sack I. Cerebral tomoelastography based on multifrequency MR elastography in two and three dimensions. Front Bioeng Biotechnol 2022; 10:1056131. [PMID: 36532573 PMCID: PMC9755504 DOI: 10.3389/fbioe.2022.1056131] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2022] [Accepted: 11/21/2022] [Indexed: 09/01/2023] Open
Abstract
Purpose: Magnetic resonance elastography (MRE) generates quantitative maps of the mechanical properties of biological soft tissues. However, published values obtained by brain MRE vary largely and lack detail resolution, due to either true biological effects or technical challenges. We here introduce cerebral tomoelastography in two and three dimensions for improved data consistency and detail resolution while considering aging, brain parenchymal fraction (BPF), systolic blood pressure, and body mass index (BMI). Methods: Multifrequency MRE with 2D- and 3D-tomoelastography postprocessing was applied to the brains of 31 volunteers (age range: 22-61 years) for analyzing the coefficient of variation (CV) and effects of biological factors. Eleven volunteers were rescanned after 1 day and 1 year to determine intraclass correlation coefficient (ICC) and identify possible long-term changes. Results: White matter shear wave speed (SWS) was slightly higher in 2D-MRE (1.28 ± 0.02 m/s) than 3D-MRE (1.22 ± 0.05 m/s, p < 0.0001), with less variation after 1 day in 2D (0.33 ± 0.32%) than in 3D (0.96 ± 0.66%, p = 0.004), which was also reflected in a slightly lower CV and higher ICC in 2D (1.84%, 0.97 [0.88-0.99]) than in 3D (3.89%, 0.95 [0.76-0.99]). Remarkably, 3D-MRE was sensitive to a decrease in white matter SWS within only 1 year, whereas no change in white matter volume was observed during this follow-up period. Across volunteers, stiffness correlated with age and BPF, but not with blood pressure and BMI. Conclusion: Cerebral tomoelastography provides high-resolution viscoelasticity maps with excellent consistency. Brain MRE in 2D shows less variation across volunteers in shorter scan times than 3D-MRE, while 3D-MRE appears to be more sensitive to subtle biological effects such as aging.
Collapse
Affiliation(s)
- Helge Herthum
- Berlin Center for Advanced Neuroimaging, Charité—Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin, Berlin Institute of Health, Humboldt-Universität zu Berlin, Berlin, Germany
- Bernstein Center for Computational Neuroscience, Charité—Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin, Berlin Institute of Health, Humboldt-Universität zu Berlin, Berlin, Germany
- Institute of Medical Informatics, Charité—Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin, Berlin Institute of Health, Humboldt-Universität zu Berlin, Berlin, Germany
| | - Stefan Hetzer
- Berlin Center for Advanced Neuroimaging, Charité—Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin, Berlin Institute of Health, Humboldt-Universität zu Berlin, Berlin, Germany
- Bernstein Center for Computational Neuroscience, Charité—Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin, Berlin Institute of Health, Humboldt-Universität zu Berlin, Berlin, Germany
| | - Bernhard Kreft
- Department of Radiology, Charité—Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin, Berlin Institute of Health, Humboldt-Universität zu Berlin, Berlin, Germany
| | - Heiko Tzschätzsch
- Department of Radiology, Charité—Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin, Berlin Institute of Health, Humboldt-Universität zu Berlin, Berlin, Germany
| | - Mehrgan Shahryari
- Department of Radiology, Charité—Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin, Berlin Institute of Health, Humboldt-Universität zu Berlin, Berlin, Germany
| | - Tom Meyer
- Department of Radiology, Charité—Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin, Berlin Institute of Health, Humboldt-Universität zu Berlin, Berlin, Germany
| | - Steffen Görner
- Department of Radiology, Charité—Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin, Berlin Institute of Health, Humboldt-Universität zu Berlin, Berlin, Germany
| | - Hennes Neubauer
- Department of Radiology, Charité—Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin, Berlin Institute of Health, Humboldt-Universität zu Berlin, Berlin, Germany
| | - Jing Guo
- Department of Radiology, Charité—Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin, Berlin Institute of Health, Humboldt-Universität zu Berlin, Berlin, Germany
| | - Jürgen Braun
- Institute of Medical Informatics, Charité—Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin, Berlin Institute of Health, Humboldt-Universität zu Berlin, Berlin, Germany
| | - Ingolf Sack
- Department of Radiology, Charité—Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin, Berlin Institute of Health, Humboldt-Universität zu Berlin, Berlin, Germany
| |
Collapse
|
2
|
Du Q, Bel-Brunon A, Lambert SA, Hamila N. Numerical simulation of wave propagation through interfaces using the extended finite element method for magnetic resonance elastography. THE JOURNAL OF THE ACOUSTICAL SOCIETY OF AMERICA 2022; 151:3481. [PMID: 35649898 PMCID: PMC9381142 DOI: 10.1121/10.0011392] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/03/2023]
Abstract
Magnetic resonance elastography (MRE) is an elasticity imaging technique for quantitatively assessing the stiffness of human tissues. In MRE, finite element method (FEM) is widely used for modeling wave propagation and stiffness reconstruction. However, in front of inclusions with complex interfaces, FEM can become burdensome in terms of the model partition and computationally expensive. In this work, we implement a formulation of FEM, known as the eXtended finite element method (XFEM), which is a method used for modeling discontinuity like crack and heterogeneity. Using a level-set method, it makes the interface independent of the mesh, thus relieving the meshing efforts. We investigate this method in two studies: wave propagation across an oblique linear interface and stiffness reconstruction of a random-shape inclusion. In the first study, numerical results by XFEM and FEM models revealing the wave conversion rules at linear interface are presented and successfully compared to the theoretical predictions. The second study, investigated in a pseudo-practical application, demonstrates further the applicability of XFEM in MRE and the convenience, accuracy, and speed of XFEM with respect to FEM. XFEM can be regarded as a promising alternative to FEM for inclusion modeling in MRE.
Collapse
Affiliation(s)
- Quanshangze Du
- Univ Lyon, INSA Lyon, CNRS, LaMCoS, UMR5259, 69621 Villeurbanne, France
| | - Aline Bel-Brunon
- Univ Lyon, INSA Lyon, CNRS, LaMCoS, UMR5259, 69621 Villeurbanne, France
- Electronic mail:
| | - Simon Auguste Lambert
- Université de Lyon, INSA Lyon, Université Claude Bernard Lyon 1, Ecole Centrale de Lyon, CNRS, Ampère UMR5005, Villeurbanne, France
| | | |
Collapse
|
3
|
Hiscox LV, McGarry MDJ, Johnson CL. Evaluation of cerebral cortex viscoelastic property estimation with nonlinear inversion magnetic resonance elastography. Phys Med Biol 2022; 67. [PMID: 35316794 PMCID: PMC9208651 DOI: 10.1088/1361-6560/ac5fde] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2021] [Accepted: 03/22/2022] [Indexed: 12/25/2022]
Abstract
Objective. Magnetic resonance elastography (MRE) of the brain has shown promise as a sensitive neuroimaging biomarker for neurodegenerative disorders; however, the accuracy of performing MRE of the cerebral cortex warrants investigation due to the unique challenges of studying thinner and more complex geometries.Approach. A series of realistic, whole-brain simulation experiments are performed to examine the accuracy of MRE to measure the viscoelasticity (shear stiffness,μ, and damping ratio, ξ) of cortical structures predominantly effected in aging and neurodegeneration. Variations to MRE spatial resolution and the regularization of a nonlinear inversion (NLI) approach are examined.Main results. Higher-resolution MRE displacement data (1.25 mm isotropic resolution) and NLI with a low soft prior regularization weighting provided minimal measurement error compared to other studied protocols. With the optimized protocol, an average error inμand ξ was 3% and 11%, respectively, when compared with the known ground truth. Mid-line structures, as opposed to those on the cortical surface, generally display greater error. Varying model boundary conditions and reducing the thickness of the cortex by up to 0.67 mm (which is a realistic portrayal of neurodegenerative pathology) results in no loss in reconstruction accuracy.Significance. These experiments establish quantitative guidelines for the accuracy expected ofin vivoMRE of the cortex, with the proposed method providing valid MRE measures for future investigations into cortical viscoelasticity and relationships with health, cognition, and behavior.
Collapse
Affiliation(s)
- Lucy V Hiscox
- Department of Biomedical Engineering, University of Delaware, Newark, DE, United States of America.,Department of Psychology, University of Bath, Bath, United Kingdom
| | - Matthew D J McGarry
- Thayer School of Engineering, Dartmouth College, Hanover, NH, United States of America
| | - Curtis L Johnson
- Department of Biomedical Engineering, University of Delaware, Newark, DE, United States of America
| |
Collapse
|
4
|
Herthum H, Carrillo H, Osses A, Uribe S, Sack I, Bertoglio C. Multiple motion encoding in phase-contrast MRI: A general theory and application to elastography imaging. Med Image Anal 2022; 78:102416. [PMID: 35334444 DOI: 10.1016/j.media.2022.102416] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2021] [Revised: 12/23/2021] [Accepted: 03/01/2022] [Indexed: 01/04/2023]
Abstract
While MRI allows to encode the motion of tissue in the magnetization's phase, it remains yet a challenge to obtain high fidelity motion images due to wraps in the phase for high encoding efficiencies. Therefore, we propose an optimal multiple motion encoding method (OMME) and exemplify it in Magnetic Resonance Elastography (MRE) data. OMME is formulated as a non-convex least-squares problem for the motion using an arbitrary number of phase-contrast measurements with different motion encoding gradients (MEGs). The mathematical properties of OMME are proved in terms of standard deviation and dynamic range of the motion's estimate for arbitrary MEGs combination which are confirmed using synthetically generated data. OMME's performance is assessed on MRE data from in vivo human brain experiments and compared to dual encoding strategies. The unwrapped images are further used to reconstruct stiffness maps and compared to the ones obtained using conventional unwrapping methods. OMME allowed to successfully combine several MRE phase images with different MEGs, outperforming dual encoding strategies in either motion-to-noise ratio (MNR) or number of successfully reconstructed voxels with good noise stability. This lead to stiffness maps with greater resolution of details than obtained with conventional unwrapping methods. The proposed OMME method allows for a flexible and noise robust increase in the dynamic range and thus provides wrap-free phase images with high MNR. In MRE, the method may be especially suitable when high resolution images with high MNR are needed.
Collapse
Affiliation(s)
- Helge Herthum
- Institute of Medical Informatics, Charité - Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin, Humboldt-Universitt zu Berlin, and Berlin Institute of Health, Berlin, 10117, Germany
| | - Hugo Carrillo
- Center for Mathematical Modeling, Universidad de Chile, Santiago, 8370456, Chile; Bernoulli Institute, University of Groningen, Groningen, 9747AG, the Netherlands
| | - Axel Osses
- Center for Mathematical Modeling, Universidad de Chile, Santiago, 8370456, Chile; Department of Mathematical Engineering, Universidad de Chile, Santiago, 8370456, Chile; ANID - Millennium Nucleus in Cardiovascular Magnetic Resonance, Santiago, 7820436, Chile; ANID - Millenium Nucleus in Applied Control and Inverse Problems ACIP, Santiago, 7820436, Chile
| | - Sergio Uribe
- ANID - Millennium Nucleus in Cardiovascular Magnetic Resonance, Santiago, 7820436, Chile; Biomedical Imaging Center, Pontificia Universidad Católica de Chile, Santiago, 7820436, Chile
| | - Ingolf Sack
- Institute of Medical Informatics, Charité - Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin, Humboldt-Universitt zu Berlin, and Berlin Institute of Health, Berlin, 10117, Germany
| | - Cristóbal Bertoglio
- Bernoulli Institute, University of Groningen, Groningen, 9747AG, the Netherlands.
| |
Collapse
|