Mapping heterogenous anisotropic tissue mechanical properties with transverse isotropic nonlinear inversion MR elastography

The Networking Brain: How Extracellular Matrix, Cellular Networks, and Vasculature Shape the In Vivo Mechanical Properties of the Brain

Mechanically, the brain is characterized by both solid and fluid properties. The resulting unique material behavior fosters proliferation, differentiation, and repair of cellular and vascular networks, and optimally protects them from damaging shear forces. Magnetic resonance elastography (MRE) is a noninvasive imaging technique that maps the mechanical properties of the brain in vivo. MRE studies have shown that abnormal processes such as neuronal degeneration, demyelination, inflammation, and vascular leakage lead to tissue softening. In contrast, neuronal proliferation, cellular network formation, and higher vascular pressure result in brain stiffening. In addition, brain viscosity has been reported to change with normal blood perfusion variability and brain maturation as well as disease conditions such as tumor invasion. In this article, the contributions of the neuronal, glial, extracellular, and vascular networks are discussed to the coarse-grained parameters determined by MRE. This reductionist multi-network model of brain mechanics helps to explain many MRE observations in terms of microanatomical changes and suggests that cerebral viscoelasticity is a suitable imaging marker for brain disease.

Post-mortem changes of anisotropic mechanical properties in the porcine brain assessed by MR elastography

Knowledge of the mechanical properties of brain tissue in vivo is essential to understanding the mechanisms underlying traumatic brain injury (TBI) and to creating accurate computational models of TBI and neurosurgical simulation. Brain white matter, which is composed of aligned, myelinated, axonal fibers, is structurally anisotropic. White matter in vivo also exhibits mechanical anisotropy, as measured by magnetic resonance elastography (MRE), but measurements of anisotropy obtained by mechanical testing of white matter ex vivo have been inconsistent. The minipig has a gyrencephalic brain with similar white matter and gray matter proportions to humans and therefore provides a relevant model for human brain mechanics. In this study, we compare estimates of anisotropic mechanical properties of the minipig brain obtained by identical, non-invasive methods in the live (in vivo) and dead animals (in situ). To do so, we combine wave displacement fields from MRE and fiber directions derived from diffusion tensor imaging (DTI) with a finite element-based, transversely-isotropic nonlinear inversion (TI-NLI) algorithm. Maps of anisotropic mechanical properties in the minipig brain were generated for each animal alive and at specific times post-mortem. These maps show that white matter is stiffer, more dissipative, and more anisotropic than gray matter when the minipig is alive, but that these differences largely disappear post-mortem, with the exception of tensile anisotropy. Overall, brain tissue becomes stiffer, less dissipative, and less mechanically anisotropic post-mortem. These findings emphasize the importance of testing brain tissue properties in vivo.

Statement of Significance

In this study, MRE and DTI in the minipig were combined to estimate, for the first time, anisotropic mechanical properties in the living brain and in the same brain after death. Significant differences were observed in the anisotropic behavior of brain tissue post-mortem. These results demonstrate the importance of measuring brain tissue properties in vivo as well as ex vivo, and provide new quantitative data for the development of computational models of brain biomechanics.

Mechanical Properties of White Matter Tracts in Aging Assessed via Anisotropic MR Elastography

Magnetic resonance elastography (MRE) is a promising neuroimaging technique to probe tissue microstructure, which has revealed widespread softening with loss of structural integrity in the aging brain. Traditional MRE approaches assume mechanical isotropy. However, white matter is known to be anisotropic from aligned, myelinated axonal bundles, which can lead to uncertainty in mechanical property estimates in these areas when using isotropic MRE. Recent advances in anisotropic MRE now allow for estimation of shear and tensile anisotropy, along with substrate shear modulus, in white matter tracts. The objective of this study was to investigate age-related differences in anisotropic mechanical properties in human brain white matter tracts for the first time. Anisotropic mechanical properties in all tracts were found to be significantly lower in older adults compared to young adults, with average property differences ranging between 0.028-0.107 for shear anisotropy and between 0.139-0.347 for tensile anisotropy. Stiffness perpendicular to the axonal fiber direction was also significantly lower in older age, but only in certain tracts. When compared with fractional anisotropy measures from diffusion tensor imaging, we found that anisotropic MRE measures provided additional, complementary information in describing differences between the white matter integrity of young and older populations. Anisotropic MRE provides a new tool for studying white matter structural integrity in aging and neurodegeneration.

Full-field noise-correlation elastography for in-plane mechanical anisotropy imaging

Elastography contrast imaging has great potential for the detection and characterization of abnormalities in soft biological tissues to help physicians in diagnosis. Transient shear-waves elastography has notably shown promising results for a range of clinical applications. In biological soft tissues such as muscle, high mechanical anisotropy implies different stiffness estimations depending on the direction of the measurement. In this study, we propose the evolution of a noise-correlation elastography approach for in-plane anisotropy mapping. This method is shown to retrieve anisotropy from simulation images before being validated on agarose anisotropic tissue-mimicking phantoms, and the first results on in-vivo biological fibrous tissues are presented.

Unravelling anisotropic nonlinear shear elasticity in muscles: Towards a non-invasive assessment of stress in living organisms

Acoustoelasticity theory describes propagation of shear waves in uniaxially stressed medium and allows the retrieval of nonlinear elastic coefficients of tissues. In transverse isotropic medium such as muscles the theory leads to 9 different configurations of propagating shear waves (stress axis vs. fibers axis vs. shear wave polarization axis vs. shear wave propagation axis). In this work we propose to use 4 configurations to quantify these nonlinear parameters ex vivo and in vivo. Ex vivo experiments combining ultrasound shear wave elastography and mechanical testing were conducted on iliopsoas pig muscles to quantify three third-order nonlinear coefficients A, H and K that are possibly linked to the architectural structure of muscles. In vivo experiments were performed with human volunteers on biceps brachii during a stretching exercise on an ergometer. A combination of the third order nonlinear elastic parameters was assessed. The knowledge of this nonlinear elastic parameters paves the way to quantify in vivo the local forces produced by muscle during exercise, contraction or movements.

Magnetic Resonance Elastography in the Study of Neurodegenerative Diseases

Neurodegenerative diseases such as Alzheimer's disease (AD) and Parkinson's disease (PD) present a major health burden to society. Changes in brain structure and cognition are generally only observed at the late stage of the disease. Although advanced magnetic resonance imaging (MRI) techniques such as diffusion imaging may allow identification of biomarkers at earlier stages of neurodegeneration, early diagnosis is still challenging. Magnetic resonance elastography (MRE) is a noninvasive MRI technique for studying the mechanical properties of tissues by measuring the wave propagation induced in the tissues using a purpose-built actuator. Here, we present a systematic review of preclinical and clinical studies in which MRE has been applied to study neurodegenerative diseases. Actuator systems for data acquisition, inversion algorithms for data analysis, and sample demographics are described and tissue stiffness measures obtained for the whole brain and internal structures are summarized. A total of six animal studies and eight human studies have been published. The animal studies refer to 123 experimental animals (68 AD and 55 PD) and 121 wild-type animals, while the human studies refer to 142 patients with neurodegenerative disease (including 56 AD and 17 PD) and 166 controls. The animal studies are consistent in the reporting of decreased stiffness of the hippocampal region in AD mice. However, in terms of disease progression, although consistent decreases in either storage modulus or shear modulus magnitude are reported for whole brain, there is variation in the results reported for the hippocampal region. The clinical studies are consistent in reports of a significant decrease in either whole brain storage modulus or shear modulus magnitude, in both AD and PD and with different brain structures affected in different neurodegenerative diseases. MRE studies of neurodegenerative diseases are still in their infancy, and in future it will be interesting to investigate potential relationships between brain mechanical properties and clinical measures, which may help elucidate the mechanisms underlying onset and progression of neurodegenerative diseases. EVIDENCE LEVEL: 1. TECHNICAL EFFICACY: Stage 2.

Individual Muscle Force Estimation in the Human Forearm Using Multi-Muscle MR Elastography (MM-MRE)

OBJECTIVE

To establish the sensitivity of magnetic resonance elastography (MRE) to active muscle contraction in multiple muscles of the forearm.

METHODS

We combined MRE of forearm muscles with an MRI-compatible device, the MREbot, to simultaneously measure the mechanical properties of tissues in the forearm and the torque applied by the wrist joint during isometric tasks. We measured shear wave speed of thirteen forearm muscles via MRE in a series of contractile states and wrist postures and fit these outputs to a force estimation algorithm based on a musculoskeletal model.

RESULTS

Shear wave speed changed significantly upon several factors, including whether the muscle was recruited as an agonist or antagonist (p = 0.0019), torque amplitude (p = <0.0001), and wrist posture (p = 0.0002). Shear wave speed increased significantly during both agonist (p = <0.0001) and antagonist (p = 0.0448) contraction. Additionally, there was a greater increase in shear wave speed at greater levels of loading. The variations due to these factors indicate the sensitivity to functional loading of muscle. Under the assumption of a quadratic relationship between shear wave speed and muscle force, MRE measurements accounted for an average of 70% of the variance in the measured joint torque.

CONCLUSION

This study shows the ability of MM-MRE to capture variations in individual muscle shear wave speed due to muscle activation and presents a method to estimate individual muscle force through MM-MRE derived measurements of shear wave speed.

SIGNIFICANCE

MM-MRE could be used to establish normal and abnormal muscle co-contraction patterns in muscles of the forearm controlling hand and wrist function.

MR Elastography With Optimization-Based Phase Unwrapping and Traveling Wave Expansion-Based Neural Network (TWENN)

Magnetic Resonance Elastography (MRE) can characterize biomechanical properties of soft tissue for disease diagnosis and treatment planning. However, complicated wavefields acquired from MRE coupled with noise pose challenges for accurate displacement extraction and modulus estimation. Using optimization-based displacement extraction and Traveling Wave Expansion-based Neural Network (TWENN) modulus estimation, we propose a new pipeline for processing MRE images. An objective function with Dual Data Consistency (Dual-DC) has been used to ensure accurate phase unwrapping and displacement extraction. For the estimation of complex wavenumbers, a complex-valued neural network with displacement covariance as an input has been developed. A model of traveling wave expansion is used to generate training datasets for the network with varying levels of noise. The complex shear modulus map is obtained through fusion of multifrequency and multidirectional data. Validation using brain and liver simulation images demonstrates the practical value of the proposed pipeline, which can estimate the biomechanical properties with minimal root-mean-square errors when compared to state-of-the-art methods. Applications of the proposed method for processing MRE images of phantom, brain, and liver reveal clear anatomical features, robustness to noise, and good generalizability of the pipeline.

Mechanical stiffness and anisotropy measured by MRE during brain development in the minipig

The relationship between brain development and mechanical properties of brain tissue is important, but remains incompletely understood, in part due to the challenges in measuring these properties longitudinally over time. In addition, white matter, which is composed of aligned, myelinated, axonal fibers, may be mechanically anisotropic. Here we use data from magnetic resonance elastography (MRE) and diffusion tensor imaging (DTI) to estimate anisotropic mechanical properties in six female Yucatan minipigs at ages from 3 to 6 months. Fiber direction was estimated from the principal axis of the diffusion tensor in each voxel. Harmonic shear waves in the brain were excited by three different configurations of a jaw actuator and measured using a motion-sensitive MR imaging sequence. Anisotropic mechanical properties are estimated from displacement field and fiber direction data with a finite element- based, transversely-isotropic nonlinear inversion (TI-NLI) algorithm. TI-NLI finds spatially resolved TI material properties that minimize the error between measured and simulated displacement fields. Maps of anisotropic mechanical properties in the minipig brain were generated for each animal at all four ages. These maps show that white matter is more dissipative and anisotropic than gray matter, and reveal significant effects of brain development on brain stiffness and structural anisotropy. Changes in brain mechanical properties may be a fundamental biophysical signature of brain development.

Imaging the subcellular viscoelastic properties of mouse oocytes

In recent years, cellular biomechanical properties have been investigated as an alternative to morphological assessments for oocyte selection in reproductive science. Despite the high relevance of cell viscoelasticity characterization, the reconstruction of spatially distributed viscoelastic parameter images in such materials remains a major challenge. Here, a framework for mapping viscoelasticity at the subcellular scale is proposed and applied to live mouse oocytes. The strategy relies on the principles of optical microelastography for imaging in combination with the overlapping subzone nonlinear inversion technique for complex-valued shear modulus reconstruction. The three-dimensional nature of the viscoelasticity equations was accommodated by applying an oocyte geometry-based 3D mechanical motion model to the measured wave field. Five domains-nucleolus, nucleus, cytoplasm, perivitelline space, and zona pellucida-could be visually differentiated in both oocyte storage and loss modulus maps, and statistically significant differences were observed between most of these domains in either property reconstruction. The method proposed herein presents excellent potential for biomechanical-based monitoring of oocyte health and complex transformations across lifespan. It also shows appreciable latitude for generalization to cells of arbitrary shape using conventional microscopy equipment.

Transversely-isotropic brain in vivo MR elastography with anisotropic damping

Measuring tissue parameters from increasingly sophisticated mechanical property models may uncover new contrast mechanisms with clinical utility. Building on previous work on in vivo brain MR elastography (MRE) with a transversely-isotropic with isotropic damping (TI-ID) model, we explore a new transversely-isotropic with anisotropic damping (TI-AD) model that involves six independent parameters describing direction-dependent behavior for both stiffness and damping. The direction of mechanical anisotropy is determined by diffusion tensor imaging and we fit three complex-valued moduli distributions across the full brain volume to minimize differences between measured and modeled displacements. We demonstrate spatially accurate property reconstruction in an idealized shell phantom simulation, as well as an ensemble of 20 realistic, randomly-generated simulated brains. We characterize the simulated precisions of all six parameters across major white matter tracts to be high, suggesting that they can be measured independently with acceptable accuracy from MRE data. Finally, we present in vivo anisotropic damping MRE reconstruction data. We perform t-tests on eight repeated MRE brain exams on a single-subject, and find that the three damping parameters are statistically distinct for most tracts, lobes and the whole brain. We also show that population variations in a 17-subject cohort exceed single-subject measurement repeatability for most tracts, lobes and whole brain, for all six parameters. These results suggest that the TI-AD model offers new information that may support differential diagnosis of brain diseases.

In vivoestimation of anisotropic mechanical properties of the gastrocnemius during functional loading with MR elastography

Objective.In vivoimaging assessments of skeletal muscle structure and function allow for longitudinal quantification of tissue health. Magnetic resonance elastography (MRE) non-invasively quantifies tissue mechanical properties, allowing for evaluation of skeletal muscle biomechanics in response to loading, creating a better understanding of muscle functional health.Approach. In this study, we analyze the anisotropic mechanical response of calf muscles using MRE with a transversely isotropic, nonlinear inversion algorithm (TI-NLI) to investigate the role of muscle fiber stiffening under load. We estimate anisotropic material parameters including fiber shear stiffness (μ1), substrate shear stiffness (μ2), shear anisotropy (ϕ), and tensile anisotropy (ζ) of the gastrocnemius muscle in response to both passive and active tension.Main results. In passive tension, we found a significant increase inμ1,ϕ,andζwith increasing muscle length. While in active tension, we observed increasingμ2and decreasingϕandζduring active dorsiflexion and plantarflexion-indicating less anisotropy-with greater effects when the muscles act as agonist.Significance. The study demonstrates the ability of this anisotropic MRE method to capture the multifaceted mechanical response of skeletal muscle to tissue loading from muscle lengthening and contraction.