Estimating the viscoelastic properties of the human brain at 7 T MRI using intrinsic MRE and nonlinear inversion

Intrinsic actuation magnetic resonance elastography (MRE) is a phase-contrast MRI technique that allows for in vivo quantification of mechanical properties of the brain by exploiting brain motion that arise naturally due to the cardiac pulse. The mechanical properties of the brain reflect its tissue microstructure, making it a potentially valuable parameter in studying brain disease. The main purpose of this study was to assess the feasibility of reconstructing the viscoelastic properties of the brain using high-quality 7 T MRI displacement measurements, obtained using displacement encoding with stimulated echoes (DENSE) and intrinsic actuation. The repeatability and sensitivity of the method for detecting normal regional variation in brain tissue properties was assessed as secondary goal. The displacement measurements used in this analysis were previously acquired for a separate study, where eight healthy subjects (27 +/- 7 years) were imaged with repeated scans (spatial resolution approx. 2 mm isotropic, temporal resolution 75 ms, motion sensitivity 0.35 mm/2 pi for displacements in anterior-posterior and left-right directions, and 0.7 mm/2 pi for feet-head displacements). The viscoelastic properties of the brain were estimated using a subzone based non-linear inversion scheme. The results show comparable consistency to that of extrinsic MRE between the viscoelastic property maps obtained from repeated displacement measurements. The shear stiffness maps showed fairly consistent spatial patterns. The whole-brain repeatability coefficient (RC) for shear stiffness was (mean +/- standard deviation) 8 +/- 8% relative to the mean whole-brain stiffness, and the damping ratio RC was 28 +/- 17% relative to the whole-brain damping ratio. The shear stiffness maps showed similar statistically significant regional trends as demonstrated in a publicly available atlas of viscoelastic properties obtained with extrinsic actuation MRE at 50 Hz. The damping ratio maps showed less consistency, likely due to data-model mismatch of describing the brain as a viscoelastic material under low frequencies. While artifacts induced by fluid flow within the brain remain a limitation of the technique in its current state, intrinsic actuation based MRE allow for consistent and repeatable estimation of the mechanical properties of the brain. The method provides enough sensitivity to investigate regional variation in such properties in the normal brain, which is likely sufficient to also investigate pathological changes.