Evaluation of cerebral cortex viscoelastic property estimation with nonlinear inversion magnetic resonance elastography

Viscoelastic polyacrylamide MR elastography phantoms with tunable damping ratio independent of shear stiffness

Physiologically modeled test samples with known properties and characteristics, or phantoms, are essential for developing sensitive, repeatable, and accurate quantitative MRI techniques. Magnetic resonance elastography (MRE) is one such technique used to estimate tissue mechanical properties, and it is advantageous to use phantoms with independently tunable mechanical properties to benchmark the accuracy of MRE methods. Phantoms with tunable shear stiffness are commonly used for MRE, but tuning the viscosity or damping ratio has proven to be difficult. A promising candidate for MRE phantoms with tunable damping ratio is polyacrylamide (PAA). While pure PAA has very low attenuation, viscoelastic hydrogels have been made by entrapping linear polyacrylamide strands (LPAA) within the PAA network. In this study, we evaluate the use of LPAA/PAA gels as physiologically accurate phantoms with tunable damping ratio, independent of shear stiffness, via MRE. Phantoms were made with 15.3 wt% PAA while the LPAA concentration ranged from 4.5 wt% to 8.0 wt%. MRE was performed at 9.4 T with 400 Hz vibration on all phantoms revealing a strong, positive correlation between damping ratio and LPAA content (p < 0.001). There was no significant correlation between shear stiffness and LPAA content, confirming a constant PAA concentration yielded constant shear stiffness. Rheometry at 10 Hz was performed to verify the damping ratio of the phantoms. Nearly identical slopes for damping ratio versus LPAA content were found from both MRE and rheometry (0.0073 and 0.0075 respectively). Ultimately, this study validates the adaptation of polyacrylamide gels into physiologically-relevant MRE phantoms to enable testing of MRE estimates of damping ratio.

Differential effect of dementia etiology on cortical stiffness as assessed by MR elastography

BACKGROUND

Aging and dementia involve the disruption of brain molecular pathways leading to the alterations in tissue composition and gross morphology of the brain. Phenotypic and biomarker overlap between various etiologies of dementia supports a need for new modes of information to more accurately distinguish these disorders. Brain mechanical properties, which can be measured noninvasively by MR elastography, represent one understudied feature that are sensitive to neurodegenerative processes. In this study, we used two stiffness estimation schemes to test the hypothesis that different etiologies of dementia are associated with unique patterns of mechanical alterations across the cerebral cortex.

METHODS

MR elastography data were acquired for six clinical groups including amyloid-negative cognitively unimpaired (CU), amyloid-positive cognitively unimpaired (A + CU), amyloid-positive participants with mild cognitive impairment (A + MCI), amyloid-positive participants with Alzheimer's clinical syndrome (A + ACS), dementia with Lewy bodies (DLB), and frontotemporal dementia (FTD). Stiffness maps were computed using two neural network inversions with the objective to at least partially separate the parenchyma-specific and morphological effects of neurodegeneration on mechanical property estimates. A tissue-confined inversion algorithm was designed to obtain the best estimate of stiffness in the brain parenchyma itself, while a regionally-aware inversion algorithm was used to measure the tissue stiffness along with the surroundings. Mean stiffness of 15 bilateral gray matter cortical regions were considered for statistical analysis. First, we tested the hypothesis that cortical stiffness changes in the aging brain. Next, we tested the overall study hypothesis by first comparing stiffness in each clinical group to the CU group, and then comparing the clinical groups against one another. Finally, we assessed the spatial and statistical overlap between atrophy and stiffness changes for both inversions.

RESULTS

Cortical brain regions become softer with age for both inversions with larger effects observed using regionally-aware stiffness. Stiffness decreases in the range 0.010-0.027 kPa per year were observed. Pairwise comparisons of each clinical group with cognitively unimpaired participants demonstrated 5 statistically significant differences in stiffness for tissue-confined measurements and 19 statistically different stiffness changes for the regionally-aware stiffness measurements. Pairwise comparisons between clinical groups further demonstrated unique patterns of stiffness differences. Analysis of the atrophy-versus-stiffness relationship showed that regionally-aware stiffness measurements exhibit higher sensitivity to neurodegeneration with findings that are not fully explained by partial volume effects or atrophy.

CONCLUSIONS

Both tissue-confined and regionally-aware stiffness estimates exhibited unique and complementary stiffness differences in various etiologies of dementia. Our results suggest that mechanical alterations measured by MRE reflect both tissue-specific differences as well as environmental effects. Multi-inversion schemes in MRE may provide new insights into the relationships between neuropathology and brain biomechanics.

Mapping brain mechanical property maturation from childhood to adulthood

Magnetic resonance elastography (MRE) is a phase contrast MRI technique which uses external palpation to create maps of brain mechanical properties noninvasively and in vivo. These mechanical properties are sensitive to tissue microstructure and reflect tissue integrity. MRE has been used extensively to study aging and neurodegeneration, and to assess individual cognitive differences in adults, but little is known about mechanical properties of the pediatric brain. Here we use high-resolution MRE imaging in participants of ages ranging from childhood to adulthood to understand brain mechanical properties across brain maturation. We find that brain mechanical properties differ considerably between childhood and adulthood, and that neuroanatomical subregions have differing maturational trajectories. Overall, we observe lower brain stiffness and greater brain damping ratio with increasing age from 5 to 35 years. Gray and white matter change differently during maturation, with larger changes occurring in gray matter for both stiffness and damping ratio. We also found that subregions of cortical and subcortical gray matter change differently, with the caudate and thalamus changing the most with age in both stiffness and damping ratio, while cortical subregions have different relationships with age, even between neighboring regions. Understanding how brain mechanical properties mature using high-resolution MRE will allow for a deeper understanding of the neural substrates supporting brain function at this age and can inform future studies of atypical maturation.

OSCILLATE: A low-rank approach for accelerated magnetic resonance elastography

PURPOSE

MR elastography (MRE) is a technique to characterize brain mechanical properties in vivo. Due to the need to capture tissue deformation in multiple directions over time, MRE is an inherently long acquisition, which limits achievable resolution and use in challenging populations. The purpose of this work is to develop a method for accelerating MRE acquisition by using low-rank image reconstruction to exploit inherent spatiotemporal correlations in MRE data.

THEORY AND METHODS

The proposed MRE sampling and reconstruction method, OSCILLATE (Observing Spatiotemporal Correlations for Imaging with Low-rank Leveraged Acceleration in Turbo Elastography), involves alternating which k-space points are sampled between each repetition by a reduction factor, ROSC. Using a predetermined temporal basis from a low-resolution navigator in a joint low-rank image reconstruction, all images can be accurately reconstructed from a reduced amount of k-space data.

RESULTS

Decomposition of MRE displacement data demonstrated that, on average, 96.1% of all energy from an MRE dataset is captured at rank L = 12 (reduced from a full rank of 24). Retrospectively undersampling data with ROSC  = 2 and reconstructing at low-rank (L = 12) yields highly accurate stiffness maps with voxel-wise error of 5.8% ± 0.7%. Prospectively undersampled data at ROSC  = 2 were successfully reconstructed without loss of material property map fidelity, with average global stiffness error of 1.0% ± 0.7% compared to fully sampled data.

CONCLUSIONS

OSCILLATE produces whole-brain MRE data at 2 mm isotropic resolution in 1 min 48 s.