Matthew D. J. McGarry

Local mechanical properties of white matter structures in the human brain

The noninvasive measurement of the mechanical properties of brain tissue using magnetic resonance elastography (MRE) has emerged as a promising method for investigating neurological disorders. To date, brain MRE investigations have been limited to reporting global mechanical properties, though quantification of the stiffness of specific structures in the white matter architecture may be valuable in assessing the localized effects of disease. This paper reports the mechanical properties of the corpus callosum and corona radiata measured in healthy volunteers using MRE and atlas-based segmentation. Both structures were found to be significantly stiffer than overall white matter, with the corpus callosum exhibiting greater stiffness and less viscous damping than the corona radiata. Reliability of both local and global measures was assessed through repeated experiments, and the coefficient of variation for each measure was less than 10%. Mechanical properties within the corpus callosum and corona radiata demonstrated correlations with measures from diffusion tensor imaging pertaining to axonal microstructure.

An octahedral shear strain-based measure of SNR for 3D MR elastography

A signal-to-noise ratio (SNR) measure based on the octahedral shear strain (the maximum shear strain in any plane for a 3D state of strain) is presented for magnetic resonance elastography (MRE), where motion-based SNR measures are commonly used. The shear strain, γ, is directly related to the shear modulus, μ, through the definition of shear stress, τ = μγ. Therefore, noise in the strain is the important factor in determining the quality of motion data, rather than the noise in the motion. Motion and strain SNR measures were found to be correlated for MRE of gelatin phantoms and the human breast. Analysis of the stiffness distributions of phantoms reconstructed from the measured motion data revealed a threshold for both strain and motion SNR where MRE stiffness estimates match independent mechanical testing. MRE of the feline brain showed significantly less correlation between the two SNR measures. The strain SNR measure had a threshold above which the reconstructed stiffness values were consistent between cases, whereas the motion SNR measure did not provide a useful threshold, primarily due to rigid body motion effects.

Magnetic Resonance Elastography of the Brain Using Multishot Spiral Readouts with Self-Navigated Motion Correction

Magnetic resonance elastography (MRE) has been introduced in clinical practice as a possible surrogate for mechanical palpation, but its application to study the human brain in vivo has been limited by low spatial resolution and the complexity of the inverse problem associated with biomechanical property estimation. Here, we report significant improvements in brain MRE data acquisition by reporting images with high spatial resolution and signal‐to‐noise ratio as quantified by octahedral shear strain metrics. Specifically, we have developed a sequence for brain MRE based on multishot, variable‐density spiral imaging, and three‐dimensional displacement acquisition and implemented a correction scheme for any resulting phase errors. A Rayleigh damped model of brain tissue mechanics was adopted to represent the parenchyma and was integrated via a finite element‐based iterative inversion algorithm. A multiresolution phantom study demonstrates the need for obtaining high‐resolution MRE data when estimating focal mechanical properties. Measurements on three healthy volunteers demonstrate satisfactory resolution of gray and white matter, and mechanical heterogeneities correspond well with white matter histoarchitecture. Together, these advances enable MRE scans that result in high‐fidelity, spatially resolved estimates of in vivo brain tissue mechanical properties, improving upon lower resolution MRE brain studies that only report volume averaged stiffness values. Magn Reson Med 70:404–412, 2013.

Subzone based magnetic resonance elastography using a Rayleigh damped material model.

PURPOSE Recently, the attenuating behavior of soft tissue has been addressed in magnetic resonance elastography by the inclusion of a damping mechanism in the methods used to reconstruct the resulting mechanical property image. To date, this mechanism has been based on a viscoelastic model for material behavior. Rayleigh, or proportional, damping provides a more generalized model for elastic energy attenuation that uses two parameters to characterize contributions proportional to elastic and inertial forces. In the case of time-harmonic vibration, these two parameters lead to both the elastic modulus and the density being complex valued (as opposed to the case of pure viscoelasticity, where only the elastic modulus is complex valued). METHODS This article presents a description of Rayleigh damping in the time-harmonic case, discussing the differences between this model and the viscoelastic damping models. In addition, the results from a subzone based Rayleigh damped elastography study of gelatin and tofu phantoms are discussed, along with preliminary results from in vivo breast data. RESULTS Both the phantom and the tissue studies presented here indicate a change in the Rayleigh damping structure, described as Rayleigh composition, between different material types, with tofu and healthy tissue showing lower Rayleigh composition values than gelatin or cancerous tissue. CONCLUSIONS It is possible that Rayleigh damping elastography and the concomitant Rayleigh composition images provide a mechanism for differentiating tissue structure in addition to measuring elastic stiffness and attenuation. Such information could be valuable in the use of Rayleigh damped magnetic resonance elastography as a diagnostic imaging tool.

3D multislab, multishot acquisition for fast, whole-brain MR elastography with high signal-to-noise efficiency

To develop an acquisition scheme for generating MR elastography (MRE) displacement data with whole‐brain coverage, high spatial resolution, and adequate signal‐to‐noise ratio (SNR) in a short scan time.

Use of a Rayleigh damping model in elastography

A Rayleigh damping model applied to magnetic resonance elastography incorporates attenuation behavior proportionally related to both elastic and inertial forces, and allows two damping parameters to be extracted from an MRI motion dataset. Under time-harmonic conditions, the model can be implemented by the use of complex shear modulus and density, whereas viscoelastic damping models commonly used in elastography consist of only a complex shear modulus, and model only a single damping effect. Simulation studies reveal that the differences between damped elastic behavior resulting from a purely complex shear modulus (CSM damping) and from a purely complex density (CD damping) become larger as the overall level of damping present (indicated by the damping ratio) increases. A plot of results generated from the finite element (FE) model indicate the relative motion differences estimated for a range of damping ratios and CSM/CD damping combinations increase with damping ratio, and can be up to 15% at a damping ratio of 50% and therefore using the correct model for a Rayleigh damped material becomes increasingly important as damping levels increase. Resonance-related effects cause values from this plot to vary by as much as 3% as parameters such as wave speed, frequency, and problem size are altered. These motion differences can be compared to expected noise levels to estimate the parameter resolution achievable by a reconstruction algorithm. An optimization-based global property reconstruction algorithm was developed, and used for testing Rayleigh damping parameter reconstructions with gaussian noise added to the simulated motion input data. The coherent motion errors resulting from altering the combination of the two damping parameters are large enough to allow accurate determination of both of the Rayleigh damping parameters with incoherent noise levels comparable to MR measurements. The accuracy achieved by the global reconstructions was significantly better than would be predicted by examining the motion differences for differing CSM/CD damping combinations, which is likely to be due to the low ratio between number of reconstructed parameters and number of noisy measurements.

Time-harmonic magnetic resonance elastography of the normal feline brain

Imaging of the mechanical properties of in vivo brain tissue could eventually lead to non-invasive diagnosis of hydrocephalus, Alzheimers disease and other pathologies known to alter the intracranial environment. The purpose of this work is to (1) use time-harmonic magnetic resonance elastography (MRE) to estimate the mechanical property distribution of cerebral tissue in the normal feline brain and (2) compare the recovered properties of grey and white matter. Various in vivo and ex vivo brain tissue property measurement strategies have led to the highly variable results that have been reported in the literature. MR elastography is an imaging technique that can estimate mechanical properties of tissue non-invasively and in vivo. Data was acquired in 14 felines and elastic parameters were estimated using a globo-regional nonlinear image reconstruction algorithm. Results fell within the range of values reported in the literature and showed a mean shear modulus across the subject group of 7-8 kPa with all but one animal falling within 5-15 kPa. White matter was statistically stiffer (p<0.01) than grey matter by about 1 kPa on a per subject basis. To the best of our knowledge, the results reported represent the most extensive set of estimates in the in vivo brain which have been based on MRE acquisition of the three-dimensional displacement field coupled to volumetric shear modulus image reconstruction achieved through nonlinear parameter estimation. However, the inter-subject variation in mean shear modulus indicates the need for further study, including the possibility of applying more advanced models to estimate the relevant tissue mechanical properties from the data.

Medial temporal lobe viscoelasticity and relational memory performance.

Structural and functional imaging studies have been among converging lines of evidence demonstrating the importance of the hippocampus in successful memory performance. The advent of a novel neuroimaging technique - magnetic resonance elastography (MRE) - now makes it possible for us to investigate the relationship between the microstructural integrity of hippocampal tissue and successful memory processing. Mechanical properties of brain tissue estimated with MRE provide a measure of the integrity of the underlying tissue microstructure and have proven to be sensitive measures of tissue health in neurodegeneration. However, until recently, MRE methods lacked sufficient resolution necessary to accurately examine specific neuroanatomical structures in the brain, and thus could not contribute to examination of specific structure-function relationships. In this study, we took advantage of recent developments in MRE spatial resolution and mechanical inversion techniques to measure the viscoelastic properties of the human hippocampus in vivo, and investigated how these properties reflect hippocampal function. Our data reveal a strong relationship between relative elastic/viscous behavior of the hippocampus and relational memory performance (N=20). This is the first report linking the mechanical properties of brain tissue with functional performance.

Including Spatial Information in Nonlinear Inversion MR Elastography Using Soft Prior Regularization

Tissue displacements required for mechanical property reconstruction in magnetic resonance elastography (MRE) are acquired in a magnetic resonance imaging (MRI) scanner, therefore, anatomical information is available from other imaging sequences. Despite its availability, few attempts to incorporate prior spatial information in the MRE reconstruction process have been reported. This paper implements and evaluates soft prior regularization (SPR), through which homogeneity in predefined spatial regions is enforced by a penalty term in a nonlinear inversion strategy. Phantom experiments and simulations show that when predefined regions are spatially accurate, recovered property values are stable for SPR weighting factors spanning several orders of magnitude, whereas inaccurate segmentation results in bias in the reconstructed properties that can be mitigated through proper choice of regularization weighting. The method was evaluated in vivo by estimating viscoelastic mechanical properties of frontal lobe gray and white matter for five repeated scans of a healthy volunteer. Segmentations of each tissue type were generated using automated software, and statistically significant differences between frontal lobe gray and white matter were found for both the storage modulus and loss modulus . Provided homogeneous property assumptions are reasonable, SPR produces accurate quantitative property estimates for tissue structures which are finer than the resolution currently achievable with fully distributed MRE.

Observation of direction-dependent mechanical properties in the human brain with multi-excitation MR elastography

Magnetic resonance elastography (MRE) has shown promise in noninvasively capturing changes in mechanical properties of the human brain caused by neurodegenerative conditions. MRE involves vibrating the brain to generate shear waves, imaging those waves with MRI, and solving an inverse problem to determine mechanical properties. Despite the known anisotropic nature of brain tissue, the inverse problem in brain MRE is based on an isotropic mechanical model. In this study, distinct wave patterns are generated in the brain through the use of multiple excitation directions in order to characterize the potential impact of anisotropic tissue mechanics on isotropic inversion methods. Isotropic inversions of two unique displacement fields result in mechanical property maps that vary locally in areas of highly aligned white matter. Investigation of the corpus callosum, corona radiata, and superior longitudinal fasciculus, three highly ordered white matter tracts, revealed differences in estimated properties between excitations of up to 33%. Using diffusion tensor imaging to identify dominant fiber orientation of bundles, relationships between estimated isotropic properties and shear asymmetry are revealed. This study has implications for future isotropic and anisotropic MRE studies of white matter tracts in the human brain.