Search

Rationale and Objectives: Hydrocephalus (HC) is caused by accumulating cerebrospinal fluid resulting in enlarged ventricles and neurological symptoms. HC can be treated via a shunt in a subset of patients; identifying which individuals will respond through noninvasive imaging would avoid complications from unsuccessful treatments. This preliminary work is a longitudinal study applying MR Elastography (MRE) to HC patients with a focus on normal pressure hydrocephalus (NPH)., Materials and Methods: Twenty-two ventriculomegaly patients were imaged and subsequently received a lumbar drain placement for cerebrospinal fluid (CSF) drainage. NPH lumbar drain responders and NPH syndrome nonresponders were categorized by clinical presentation. Displacement images were acquired using intrinsic activation (IA) MRE and poroelastic inversion recovered shear stiffness and hydraulic conductivity values. A stable IA-MRE inversion protocol was developed to produce unique solutions for both recovered properties, independent of initial estimates., Results: Property images showed significantly increased shear modulus (p = 0.003 in periventricular region, p = 0.005 in remaining cerebral tissue) and hydraulic conductivity (p = 0.04 in periventricular region) in ventriculomegaly patients compared to healthy volunteers. Baseline MRE imaging did not detect significant differences between NPH lumbar drain responders and NPH syndrome nonresponders; however, MRE time series analysis demonstrated consistent trends in average poroelastic shear modulus values over the course of the lumbar drain process in responders (initial increase, followed by a later decrease) which did not occur in nonresponders., Conclusion: These findings are indicative of acute mechanical changes in the brain resulting from CSF drainage in NPH patients., (Copyright © 2020 The Association of University Radiologists. Published by Elsevier Inc. All rights reserved.)

Magnetic resonance elastography (MRE) has been developed to noninvasively reconstruct mechanical properties for tissue and tissue-like materials over a frequency range of 10 ~200 Hz. In this work, low frequency (1~1.5 Hz) MRE activations were employed to estimate mechanical property distributions of simulated data and experimental phantoms. Nonlinear inversion (NLI) MRE algorithms based on viscoelastic and poroelastic material models were used to solve the inverse problems and recover images of the shear modulus and hydraulic conductivity. Data from a simulated phantom containing an inclusion with property contrast was carried out to study the feasibility of our low frequency actuated approach. To verify the stability of NLI algorithms for low frequency actuation, different levels of synthetic noise were added to the displacement data. Spatial distributions and property values were recovered well for noise level less than 5%. For the presented experimental phantom reconstructions with regularizations, the computed storage moduli from viscoelastic and poroelastic MRE gave similar results. Contrast was detected between inclusions and background in recovered hydraulic conductivity images. Results and findings confirm the feasibility of future in vivo neuroimaging examinations using natural cerebrovascular pulsations at cardiac frequencies, which can eliminate specialized equipment for high frequency actuation.

Intrinsic activation MR elastography (IA-MRE) is a novel technique which seeks to estimate brain mechanical properties non-invasively and without external mechanical drivers. The method eliminates actuation hardware and patient discomfort while capitalizing on the brain's intrinsic low frequency motion. This study explores low frequency actuation (1 Hz) MR elastography in phantoms and analyzes performance of non-linear inversion (NLI) of viscoelastic and poroelastic mechanical models as a framework for assessing clinical results from IA-MRE. We present results from four gelatin phantoms and report stiffness resolution of 6 mm (two measurement voxels) with a stiffness contrast ratio of 4.21 relative to the background and 9 mm (three measurement voxels) with a lower stiffness contrast ratio of near 1.77. Stiffness edge resolution was also evaluated using edge spread and line spread functions and yielded a stiffness edge response distance of 9 mm. The intraclass correlation coefficient was high (0.93) between mechanical testing and poroelastic estimates, although quantitative agreement was affected by model-data mismatch. Viscoelastic MRE at low frequencies has issues with non-uniqueness due to small inertial forces, and performed worse than poroelastic MRE in terms of inclusion detection and consistency with mechanical testing. These results present the first evaluation of MR elastography using displacement measurements from an actuation frequency less than 5 Hz and support the validity of brain IA-MRE to recover spatially resolved stiffness changes. They provide a baseline of performance in terms of standard metrics for future animal and human brain stiffness studies and analyses based on intrinsic motion.

Magnetic Resonance Elastography (MRE) detects induced periodic motions in biological tissues allowing maps of tissue mechanical properties to be derived. In-vivo MRE is commonly performed at frequencies of 30-100 Hz using external actuation, however, using cerebro-vascular pulsation at 1 Hz as a form of intrinsic actuation (IA-MRE) eliminates the need for external motion sources and simplifies data acquisition. In this study a hydraulic actuation system was developed to drive 1 Hz motions in gelatin as a tool for investigating the performance limits of IA-MRE image reconstruction under controlled conditions. Quantitative flow (QFLOW) MR techniques were used to phase encode 1 Hz motions as a function of gradient direction using 3D or 4D acquisition; 4D acquisition was twice as fast and yielded comparable motion field and concomitant image reconstruction results provided the motion signal was sufficiently strong. Per voxel motion noise floor corresponded to a displacement amplitude of about 20-30 μm. Signal to noise ratio (SNR) was 94 ± 17 for 3D and dropped to 69 ± 10 for the faster 4D acquisition, but yielded octahedral shear stress and shear modulus maps of high quality that differed by only about 20% on average. QFLOW measurements in gel phantoms were improved significantly by adding Mn(II) to mimic relaxation rates found in brain. Overall, the hydraulic 1 Hz actuation system when coupled with 4D sequence acquisition produced a fast reliable approach for future IA-MRE phantom evaluation and contrast detail studies needed to benchmark imaging performance., (Copyright © 2018 Elsevier Inc. All rights reserved.)

This study evaluated non-linear inversion MRE (NLI-MRE) based on viscoelastic governing equations to determine its sensitivity to small, low contrast inclusions and interface changes in shear storage modulus and damping ratio. Reconstruction parameters identical to those used in recent in vivo MRE studies of mechanical property variations in small brain structures were applied. NLI-MRE was evaluated on four phantoms with contrast in stiffness and damping ratio. Image contrast to noise ratio was assessed as a function of inclusion diameter and property contrast, and edge and line spread functions were calculated as measures of imaging resolution. Phantoms were constructed from silicone, agar, and tofu materials. Reconstructed property estimates were compared with independent mechanical testing using dynamic mechanical analysis (DMA). The NLI-MRE technique detected inclusions as small as 8 mm with a stiffness contrast as low as 14%. Storage modulus images also showed an interface edge response distance of 11 mm. Damping ratio images distinguished inclusions with a diameter as small as 8 mm, and yielded an interface edge response distance of 10 mm. Property differences relative to DMA tests were in the 15%-20% range in most cases. In this study, NLI-MRE storage modulus estimates resolved the smallest inclusion with the lowest stiffness contrast, and spatial resolution of attenuation parameter images was quantified for the first time. These experiments and image quality metrics establish quantitative guidelines for the accuracy expected in vivo for MRE images of small brain structures, and provide a baseline for evaluating future improvements to the NLI-MRE pipeline.

The noninvasive imaging technique of magnetic resonance elastography (MRE) was used to estimate the power law behavior of the viscoelastic properties of the human brain in vivo. The mechanical properties for four volunteers are investigated using shear waves induced over a frequency range of 10-50Hz to produce a displacement field measured by magnetic resonance motion-encoding gradients. The average storage modulus (μ R ) reconstructed with non-linear inversion (NLI) increased from approximately 0.95 to 2.58kPa over the 10-50Hz span; the average loss modulus (μ I ) also increased from 0.29 to 1.25kPa over the range. These increases were modeled by independent power law (PL) relations for μ R and μ I returning whole brain, group mean exponent values of 0.88 and 1.07 respectively. Investigation of these exponents also showed distinctly different behavior in the region of the cerebral falx compared to other brain structures., (Copyright © 2017 Elsevier Ltd. All rights reserved.)

A numerical framework for interstitial fluid pressure imaging (IFPI) in biphasic materials is investigated based on three-dimensional nonlinear finite element poroelastic inversion. The objective is to reconstruct the time-harmonic pore-pressure field from tissue excitation in addition to the elastic parameters commonly associated with magnetic resonance elastography (MRE). The unknown pressure boundary conditions (PBCs) are estimated using the available full-volume displacement data from MRE. A subzone-based nonlinear inversion (NLI) technique is then used to update mechanical and hydrodynamical properties, given the appropriate subzone PBCs, by solving a pressure forward problem (PFP). The algorithm was evaluated on a single-inclusion phantom in which the elastic property and hydraulic conductivity images were recovered. Pressure field and material property estimates had spatial distributions reflecting their true counterparts in the phantom geometry with RMS errors around 20% for cases with 5% noise, but degraded significantly in both spatial distribution and property values for noise levels > 10%. When both shear moduli and hydraulic conductivity were estimated along with the pressure field, property value error rates were as high as 58%, 85% and 32% for the three quantities, respectively, and their spatial distributions were more distorted. Opportunities for improving the algorithm are discussed.

Magnetic resonance elastography aims to non-invasively and remotely characterize the mechanical properties of living tissues. To quantitatively and regionally map the shear viscoelastic moduli in vivo, the technique must achieve proper mechanical excitation throughout the targeted tissues. Although it is straightforward, ante manibus, in close organs such as the liver or the breast, which practitioners clinically palpate already, it is somewhat fortunately highly challenging to trick the natural protective barriers of remote organs such as the brain. So far, mechanical waves have been induced in the latter by shaking the surrounding cranial bones. Here, the skull was circumvented by guiding pressure waves inside the subject's buccal cavity so mechanical waves could propagate from within through the brainstem up to the brain. Repeatable, reproducible and robust displacement fields were recorded in phantoms and in vivo by magnetic resonance elastography with guided pressure waves such that quantitative mechanical outcomes were extracted in the human brain., (© 2022 John Wiley & Sons, Ltd.)

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., (© 2023 International Society for Magnetic Resonance in Medicine.)

Background: Magnetic resonance elastography (MRE) of the brain allows quantitative measurement of tissue mechanics. Multiple studies are exploring possible applications in normal pressure hydrocephalus (NPH) in clinical and paraclinical contexts. This is of great interest in neurological surgery due to challenges related to diagnosis and prediction of treatment effects. In this scoping review, we present a topical overview and discuss the current literature, with particular attention to clinical implications and current challenges., Methods: The protocol was based on the PRISMA extension for scoping reviews. After a systematic database search (PubMed, Embase, and Web of Science), the articles were screened for relevance. Thirty articles were subject to detailed screening, and key technical and clinical data items were extracted. The inclusion criteria included the use of MRE on human subjects with NPH., Results: Seven articles were included in the final study. These studies had various objectives including the role of MRE in the assessment of regional elastic changes in NPH, shunt effect, and evaluation of NPH symptoms. MRE revealed patterns of mechanical changes in NPH that differed from other dementias. Regional MRE changes were associated with specific NPH signs and symptoms. Neurosurgical shunting caused partial normalization in tissue scaffold parameters. The studies were highly heterogeneous in technical aspects and design., Conclusion: MRE studies in NPH are still limited by few participants, variable cohorts, inconsistent methodologies, and technical challenges, but the approach shows great potential for future clinical application., (© 2021. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

We describe an efficient gradient computation for solving inverse problems arising in magnetic resonance elastography (MRE). The algorithm can be considered as a generalized 'adjoint method' based on a Lagrangian formulation. One requirement for the classic adjoint method is assurance of the self-adjoint property of the stiffness matrix in the elasticity problem. In this paper, we show this property is no longer a necessary condition in our algorithm, but the computational performance can be as efficient as the classic method, which involves only two forward solutions and is independent of the number of parameters to be estimated. The algorithm is developed and implemented in material property reconstructions using poroelastic and viscoelastic modeling. Various gradient- and Hessian-based optimization techniques have been tested on simulation, phantom and in vivo brain data. The numerical results show the feasibility and the efficiency of the proposed scheme for gradient calculation.

Purpose: In vivo MR elastography (MRE) holds promise as a neuroimaging marker. In cerebral MRE, shear waves are introduced into the brain, which also stimulate vibrations in adjacent CSF, resulting in blurring and biased stiffness values near brain surfaces. We here propose inversion-recovery MRE (IR-MRE) to suppress CSF signal and improve stiffness quantification in brain surface areas., Methods: Inversion-recovery MRE was demonstrated in agar-based phantoms with solid-fluid interfaces and 11 healthy volunteers using 31.25-Hz harmonic vibrations. It was performed by standard single-shot, spin-echo EPI MRE following 2800-ms IR preparation. Wave fields were acquired in 10 axial slices and analyzed for shear wave speed (SWS) as a surrogate marker of tissue stiffness by wavenumber-based multicomponent inversion., Results: Phantom SWS values near fluid interfaces were 7.5 ± 3.0% higher in IR-MRE than MRE (P = .01). In the brain, IR-MRE SNR was 17% lower than in MRE, without influencing parenchymal SWS (MRE: 1.38 ± 0.02 m/s; IR-MRE: 1.39 ± 0.03 m/s; P = .18). The IR-MRE tissue-CSF interfaces appeared sharper, showing 10% higher SWS near brain surfaces (MRE: 1.01 ± 0.03 m/s; IR-MRE: 1.11 ± 0.01 m/s; P < .001) and 39% smaller ventricle sizes than MRE (P < .001)., Conclusions: Our results show that brain MRE is affected by fluid oscillations that can be suppressed by IR-MRE, which improves the depiction of anatomy in stiffness maps and the quantification of stiffness values in brain surface areas. Moreover, we measured similar stiffness values in brain parenchyma with and without fluid suppression, which indicates that shear wavelengths in solid and fluid compartments are identical, consistent with the theory of biphasic poroelastic media., (© 2021 The Authors. Magnetic Resonance in Medicine published by Wiley Periodicals LLC on behalf of International Society for Magnetic Resonance in Medicine.)

Purpose: Descriptions of the structure of brain tissue as a porous cellular matrix support application of a poroelastic (PE) mechanical model which includes both solid and fluid phases. However, the majority of brain magnetic resonance elastography (MRE) studies use a single phase viscoelastic (VE) model to describe brain tissue behavior, in part due to availability of relatively simple direct inversion strategies for mechanical property estimation. A notable exception is low frequency intrinsic actuation MRE, where PE mechanical properties are imaged with a nonlinear inversion algorithm., Methods: This paper investigates the effect of model choice at each end of the spectrum of in vivo human brain actuation frequencies. Repeat MRE examinations of the brains of healthy volunteers were used to compare image quality and repeatability for each inversion model for both 50 Hz externally produced motion and ≈1 Hz intrinsic motions. Additionally, realistic simulated MRE data were generated with both VE and PE finite element solvers to investigate the effect of inappropriate model choice for ideal VE and PE materials., Results: In vivo, MRE data revealed that VE inversions appear more representative of anatomical structure and quantitatively repeatable for 50 Hz induced motions, whereas PE inversion produces better results at 1 Hz. Reasonable VE approximations of PE materials can be derived by equating the equivalent wave velocities for the two models, provided that the timescale of fluid equilibration is not similar to the period of actuation. An approximation of the equilibration time for human brain reveals that this condition is violated at 1 Hz but not at 50 Hz. Additionally, simulation experiments when using the "wrong" model for the inversion demonstrated reasonable shear modulus reconstructions at 50 Hz, whereas cross-model inversions at 1 Hz were poor quality. Attenuation parameters were sensitive to changes in the forward model at both frequencies, however, no spatial information was recovered because the mechanisms of VE and PE attenuation are different., Conclusions: VE inversions are simpler with fewer unknown properties and may be sufficient to capture the mechanical behavior of PE brain tissue at higher actuation frequencies. However, accurate modeling of the fluid phase is required to produce useful mechanical property images at the lower frequencies of intrinsic brain motions.

Identification of the mechanical properties of a viscoelastic material depends on characteristics of the observed motion field within the object in question. For certain physical and experimental configurations and certain resolutions and variance within the measurement data, the viscoelastic properties of an object may become non-identifiable. Elastographic imaging methods seek to provide maps of these viscoelastic properties based on displacement data measured by traditional imaging techniques, such as magnetic resonance and ultrasound. Here, 1D analytic solutions of the viscoelastic wave equation are used to generate displacement fields over wave conditions representative of diverse time-harmonic elastography applications. These solutions are tested through the minimization of a least squares objective function suitable for framing the elastography inverse calculation. Analysis shows that the damping ratio and the ratio of the viscoelastic wavelength to the size of the domain play critical roles in the form of this least squares objective function. In addition, it can be shown analytically that this objective function will contain local minima, which hinder discovery of the global minima via gradient descent methods., (© 2023 The Authors. International Journal for Numerical Methods in Biomedical Engineering published by John Wiley & Sons Ltd.)

Poroelastic magnetic resonance elastography is an imaging technique that could recover mechanical and hydrodynamical material properties of in vivo tissue. To date, mechanical properties have been estimated while hydrodynamical parameters have been assumed homogeneous with literature-based values. Estimating spatially-varying hydraulic conductivity would likely improve model accuracy and provide new image information related to a tissue's interstitial fluid compartment. A poroelastic model was reformulated to recover hydraulic conductivity with more appropriate fluid-flow boundary conditions. Simulated and physical experiments were conducted to evaluate the accuracy and stability of the inversion algorithm. Simulations were accurate (property errors were < 2%) even in the presence of Gaussian measurement noise up to 3%. The reformulated model significantly decreased variation in the shear modulus estimate (p << 0.001) and eliminated the homogeneity assumption and the need to assign hydraulic conductivity values from literature. Material property contrast was recovered experimentally in three different tofu phantoms and the accuracy was improved through soft-prior regularization. A frequency-dependence in hydraulic conductivity contrast was observed suggesting that fluid-solid interactions may be more prominent at low frequency. In vivo recovery of both structural and hydrodynamical characteristics of tissue could improve detection and diagnosis of neurological disorders such as hydrocephalus and brain tumors.

Purpose: 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., Theory and Methods: A 3D multislab, multishot acquisition for whole-brain MRE with 2.0 mm isotropic spatial resolution is proposed. The multislab approach allowed for the use of short repetition time to achieve very high SNR efficiency. High SNR efficiency allowed for a reduced acquisition time of only 6 min while the minimum SNR needed for inversion was maintained., Results: The mechanical property maps estimated from whole-brain displacement data with nonlinear inversion (NLI) demonstrated excellent agreement with neuroanatomical features, including the cerebellum and brainstem. A comparison with an equivalent 2D acquisition illustrated the improvement in SNR efficiency of the 3D multislab acquisition. The flexibility afforded by the high SNR efficiency allowed for higher resolution with a 1.6 mm isotropic voxel size, which generated higher estimates of brainstem stiffness compared with the 2.0 mm isotropic acquisition., Conclusion: The acquisition presented allows for the capture of whole-brain MRE displacement data in a short scan time, and may be used to generate local mechanical property estimates of neuroanatomical features throughout the brain., (Copyright © 2013 Wiley Periodicals, Inc.)

Purpose: Biological soft tissues often have a porous architecture comprising fluid and solid compartments. Upon displacement through physiological or externally induced motion, the relative motion of these compartments depends on poroelastic parameters, such as coupling density ( ρ 12 ) and tissue porosity. This study introduces inversion recovery MR elastography (IR-MRE) (1) to quantify porosity defined as fluid volume over total volume, (2) to separate externally induced shear strain fields of fluid and solid compartments, and (3) to quantify coupling density assuming a biphasic behavior of in vivo brain tissue., Theory and Methods: Porosity was measured in eight tofu phantoms and gray matter (GM) and white matter (WM) of 21 healthy volunteers. Porosity of tofu was compared to values obtained by fluid draining and microscopy. Solid and fluid shear-strain amplitudes and ρ 12 were estimated both in phantoms and in in vivo brain., Results: T 1 -based measurement of tofu porosity agreed well with reference values (R = 0.99, P < .01). Brain tissue porosity was 0.14 ± 0.02 in GM and 0.05 ± 0.01 in WM (P < .001). Fluid shear strain was found to be phase-locked with solid shear strain but had lower amplitudes in both tofu phantoms and brain tissue (P < .05). In accordance with theory, tofu and brain ρ 12 were negative., Conclusion: IR-MRE allowed for the first time separation of shear strain fields of solid and fluid compartments for measuring coupling density according to the biphasic theory of poroelasticity. Thus, IR-MRE opens horizons for poroelastography-derived imaging markers that can be used in basic research and diagnostic applications., (© 2020 The Authors. Magnetic Resonance in Medicine published by Wiley Periodicals LLC on behalf of International Society for Magnetic Resonance in Medicine.)

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., (Copyright © 2013 Elsevier Inc. All rights reserved.)

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.

In tomoelastography, to achieve a final wave speed map by combining reconstructions obtained from all spatial directions and excitation frequencies, the use of weights is inevitable. Here, a new weighting scheme, which maximizes the signal-to-noise ratio (SNR) of the final wave speed map, has been proposed. To maximize the SNR of the final wave speed map, the use of squares of estimated SNR values of reconstructed individual maps has been proposed. Therefore, derivations of the SNR of the reconstructed wave speed maps have become necessary. Considering the noise on the complex MRI signal, the SNR of the reconstructed wave speed map was formulated by an analytical approach assuming a high SNR, and the results were verified using Monte Carlo simulations (MCSs). It has been assumed that the noise remains approximately Gaussian when the image SNR is high enough, despite the nonlinear operations in tomoelastography inversion. Hence, the SNR threshold was determined by comparing the SNR computed by MCSs and analytical approximations. The weighting scheme was evaluated for accuracy, spatial resolution and SNR performances on simulated phantoms. MR elastography (MRE) experiments on two different phantoms were conducted. Wave speed maps were generated for simulated 3D human abdomen MRE data and experimental human abdomen MRE data. The simulation results demonstrated that the SNR-weighted inversion improved the SNR performance of the wave speed map by a factor of two compared to the performance of the original (i.e., amplitude-weighted) reconstruction. In the case of a low SNR, no bias occurred in the wave speed map when SNR weighting was used, whereas 10% bias occurred when the original weighting (i.e., amplitude weighting) was used. Thus, while not altering the accuracy or spatial resolution of the wave speed map with the proposed weighting method, the SNR of the wave speed map has been significantly improved., (© 2020 John Wiley & Sons, Ltd.)

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., (© 2012 Wiley Periodicals, Inc.)

Purpose: Magnetic resonance elastography (MRE) measures stiffness of soft tissues by analyzing their spatial harmonic response to externally induced shear vibrations. Many MRE methods use inversion-based reconstruction approaches, which invoke first- or second-order derivatives by finite difference operators (first- and second-FDOs) and thus give rise to a biased frequency dispersion of stiffness estimates., Methods: We here demonstrate analytically, numerically, and experimentally that FDO-based stiffness estimates are affected by (1) noise-related underestimation of values in the range of high spatial wave support, that is, at lower vibration frequencies, and (2) overestimation of values due to wave discretization at low spatial support, that is, at higher vibration frequencies., Results: Our results further demonstrate that second-FDOs are more susceptible to noise than first-FDOs and that FDO dispersion depends both on signal-to-noise ratio (SNR) and on a lumped parameter A, which is defined as wavelength over pixel size and over a number of pixels per stencil of the FDO. Analytical FDO dispersion functions are derived for optimizing A parameters at a given SNR. As a simple rule of thumb, we show that FDO artifacts are minimized when A/2 is in the range of the square root of 2SNR for the first-FDO or cubic root of 5SNR for the second-FDO., Conclusions: Taken together, the results of our study provide an analytical solution to a long-standing, well-recognized, yet unsolved problem in MRE postprocessing and might thus contribute to the ongoing quest for minimizing inversion artifacts in MRE., (© 2020 The Authors. Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine.)

Many pathologies alter the mechanical properties of tissue. Magnetic resonance elastography (MRE) has been developed to noninvasively characterize these quantities in vivo. Typically, small vibrations are induced in the tissue of interest with an external mechanical actuator. The resulting displacements are measured with phase contrast sequences and are then used to estimate the underlying mechanical property distribution. Several MRE studies have quantified brain tissue properties. However, the cranium and meninges, especially the dura, are very effective at damping externally applied vibrations from penetrating deeply into the brain. Here, we report a method, termed 'intrinsic activation', that eliminates the requirement for external vibrations by measuring the motion generated by natural blood vessel pulsation. A retrospectively gated phase contrast MR angiography sequence was used to record the tissue velocity at eight phases of the cardiac cycle. The velocities were numerically integrated via the Fourier transform to produce the harmonic displacements at each position within the brain. The displacements were then reconstructed into images of the shear modulus based on both linear elastic and poroelastic models. The mechanical properties produced fall within the range of brain tissue estimates reported in the literature and, equally important, the technique yielded highly reproducible results. The mean shear modulus was 8.1 kPa for linear elastic reconstructions and 2.4 kPa for poroelastic reconstructions where fluid pressure carries a portion of the stress. Gross structures of the brain were visualized, particularly in the poroelastic reconstructions. Intra-subject variability was significantly less than the inter-subject variability in a study of six asymptomatic individuals. Further, larger changes in mechanical properties were observed in individuals when examined over time than when the MRE procedures were repeated on the same day. Cardiac pulsation, termed intrinsic activation, produces sufficient motion to allow mechanical properties to be recovered. The poroelastic model is more consistent with the measured data from brain at low frequencies than the linear elastic model. Intrinsic activation allows MRE to be performed without a device shaking the head so the patient notices no differences between it and the other sequences in an MR examination.

Purpose: Nonlinear inversion (NLI) in MR elastography requires discretization of the displacement field for a finite element (FE) solution of the "forward problem", and discretization of the unknown mechanical property field for the iterative solution of the "inverse problem". The resolution requirements for these two discretizations are different: the forward problem requires sufficient resolution of the displacement FE mesh to ensure convergence, whereas lowering the mechanical property resolution in the inverse problem stabilizes the mechanical property estimates in the presence of measurement noise. Previous NLI implementations use the same FE mesh to support the displacement and property fields, requiring a trade-off between the competing resolution requirements., Methods: This work implements and evaluates multiresolution FE meshes for NLI elastography, allowing independent discretizations of the displacements and each mechanical property parameter to be estimated. The displacement resolution can then be selected to ensure mesh convergence, and the resolution of the property meshes can be independently manipulated to control the stability of the inversion., Results: Phantom experiments indicate that eight nodes per wavelength (NPW) are sufficient for accurate mechanical property recovery, whereas mechanical property estimation from 50 Hz in vivo brain data stabilizes once the displacement resolution reaches 1.7 mm (approximately 19 NPW). Viscoelastic mechanical property estimates of in vivo brain tissue show that subsampling the loss modulus while holding the storage modulus resolution constant does not substantially alter the storage modulus images. Controlling the ratio of the number of measurements to unknown mechanical properties by subsampling the mechanical property distributions (relative to the data resolution) improves the repeatability of the property estimates, at a cost of modestly decreased spatial resolution., Conclusions: Multiresolution NLI elastography provides a more flexible framework for mechanical property estimation compared to previous single mesh implementations.

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 (MRE) uses accumulated phases that are acquired at multiple, uniformly spaced relative phase offsets, to estimate harmonic motion information. Heavily wrapped phase occurs when the motion is large and unwrapping procedures are necessary to estimate the displacements required by MRE. Two unwrapping methods were developed and compared in this paper. The first method is a sequentially applied approach. The three-dimensional MRE phase image block for each slice was processed by two-dimensional unwrapping followed by a one-dimensional phase unwrapping approach along the phase-offset direction. This unwrapping approach generally works well for low noise data. However, there are still cases where the two-dimensional unwrapping method fails when noise is high. In this case, the baseline of the corrupted regions within an unwrapped image will not be consistent. Instead of separating the two-dimensional and one-dimensional unwrapping in a sequential approach, an interleaved three-dimensional quality-guided unwrapping method was developed to combine both the two-dimensional phase image continuity and one-dimensional harmonic motion information. The quality of one-dimensional harmonic motion unwrapping was used to guide the three-dimensional unwrapping procedures and it resulted in stronger guidance than in the sequential method. In this work, in vivo results generated by the two methods were compared.

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.

The mechanical model commonly used in magnetic resonance elastography (MRE) is linear elasticity. However, soft tissue may exhibit frequency- and direction-dependent (FDD) shear moduli in response to an induced excitation causing a purely linear elastic model to provide an inaccurate image reconstruction of its mechanical properties. The goal of this study was to characterize the effects of reconstructing FDD data using a linear elastic inversion (LEI) algorithm. Linear and FDD phantoms were manufactured and LEI images were obtained from time-harmonic MRE acquisitions with variations in frequency and driving signal amplitude. LEI responses to artificially imposed uniform phase shifts in the displacement data from both purely linear elastic and FDD phantoms were also evaluated. Of the variety of FDD phantoms considered, LEI appeared to tolerate viscoelastic data-model mismatch better than deviations caused by poroelastic and anisotropic mechanical properties in terms of visual image contrast. However, the estimated shear modulus values were substantially incorrect relative to independent mechanical measurements even in the successful viscoelastic cases and the variations in mean values with changes in experimental conditions associated with uniform phase shifts, driving signal frequency and amplitude were unpredictable. Overall, use of LEI to reconstruct data acquired in phantoms with FDD material properties provided biased results under the best conditions and significant artifacts in the worst cases. These findings suggest that the success with which LEI is applied to MRE data in tissue will depend on the underlying mechanical characteristics of the tissues and/or organs systems of clinical interest.

Imaging of the mechanical properties of in vivo brain tissue could eventually lead to non-invasive diagnosis of hydrocephalus, Alzheimer's 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., (Copyright © 2010 Elsevier Ltd. All rights reserved.)

Purpose: The clinical efficacy of breast elastography may be limited when the authors employ the assumption that soft tissues exhibit linear, frequency-independent isotropic mechanical properties during the recovery of shear modulus. Thus, the purpose of this research was to evaluate the degradation in performance incurred when linear-elastic MR reconstruction methods are applied to phantoms that are fabricated using viscoelastic materials., Methods: To develop phantoms with frequency-dependent mechanical properties, the authors measured the complex modulus of two groups of cylindrical-shaped gelatin samples over a wide frequency range (up to 1 kHz) with the established principles of time-temperature superposition (TTS). In one group of samples, the authors added varying amounts of agar (1%-4%); in the other group, the authors added varying amounts of sucrose (2.5%-20%). To study how viscosity affected the performance of the linear-elastic reconstruction method, the authors constructed an elastically heterogeneous MR phantom to simulate the case where small viscoelastic lesions were surrounded by relatively nonviscous breast tissue. The breast phantom contained four linear, viscoelastic spherical inclusions (10 mm diameter) that were embedded in normal gelatin. The authors imaged the breast phantom with a clinical prototype of a MRE system and recovered the shear-modulus distribution using the overlapping-subzone-linear-elastic image-reconstruction method. The authors compared the recovered shear modulus to that measured using the TTS method., Results: The authors demonstrated that viscoelastic phantoms could be fabricated by including sucrose in the gelation process and that small viscoelastic inclusions were visible in MR elastograms recovered using a linear-elastic MR reconstruction process; however, artifacts that degraded contrast and spatial resolution were more prominent in highly viscoelastic inclusions. The authors also established that the accuracy of the MR elastograms depended on the degree of viscosity that the inclusion exhibited., Conclusions: The results demonstrated that reconstructing shear modulus from other constitutive laws, such as viscosity, should improve both the accuracy and quality of MR elastograms of the breast.

Purpose: Recent interest in the poroelastic behavior of tissues has led to the development of magnetic resonance poroelastography (MRPE) as an alternative to single-phase MR elastographic image reconstruction. In addition to the elastic parameters (i.e., Lamé's constants) commonly associated with magnetic resonance elastography (MRE), MRPE enables estimation of the time-harmonic pore-pressure field induced by external mechanical vibration., Methods: This study presents numerical simulations that demonstrate the sensitivity of the computed displacement and pore-pressure fields to a priori estimates of the experimentally derived model parameters. In addition, experimental data collected in three poroelastic phantoms are used to assess the quantitative accuracy of MR poroelastographic imaging through comparisons with both quasistatic and dynamic mechanical tests., Results: The results indicate hydraulic conductivity to be the dominant parameter influencing the deformation behavior of poroelastic media under conditions applied during MRE. MRPE estimation of the matrix shear modulus was bracketed by the values determined from independent quasistatic and dynamic mechanical measurements as expected, whereas the contrast ratios for embedded inclusions were quantitatively similar (10%-15% difference between the reconstructed images and the mechanical tests)., Conclusions: The findings suggest that the addition of hydraulic conductivity and a viscoelastic solid component as parameters in the reconstruction may be warranted.

Magnetic resonance poroelastography (MRPE) is introduced as an alternative to single-phase model-based elastographic reconstruction methods. A 3-D finite element poroelastic inversion algorithm was developed to recover the mechanical properties of fluid-saturated tissues. The performance of this algorithm was assessed through a variety of numerical experiments, using synthetic data to probe its stability and sensitivity to the relevant model parameters. Preliminary results suggest the algorithm is robust in the presence of noise and capable of producing accurate assessments of the underlying mechanical properties in simulated phantoms. Furthermore, a 3-D time-harmonic motion field was recorded for a poroelastic phantom containing a single cylindrical inclusion and used to assess the feasibility of MRPE image reconstruction from experimental data. The elastograms obtained from the proposed poroelastic algorithm demonstrate significant improvement over linearly elastic MRE images generated using the same data. In addition, MRPE offers the opportunity to estimate the time-harmonic pressure field resulting from tissue excitation, highlighting the potential for its application in the diagnosis and monitoring of disease processes associated with changes in interstitial pressure.

Elastography is an emerging imaging technique that focuses on assessing the resistance to deformation of soft biological tissues in vivo. Magnetic resonance elastography (MRE) uses measured displacement fields resulting from low-amplitude, low-frequency (10 Hz-1 kHz) time-harmonic vibration to recover images of the elastic property distribution of tissues including breast, liver, muscle, prostate, and brain. While many soft tissues display complex time-dependent behavior not described by linear elasticity, the models most commonly employed in MRE parameter reconstructions are based on elastic assumptions. Further, elasticity models fail to include the interstitial fluid phase present in vivo. Alternative continuum models, such as consolidation theory, are able to represent tissue and other materials comprising two distinct phases, generally consisting of a porous elastic solid and penetrating fluid. MRE reconstructions of simulated elastic and poroelastic phantoms were performed to investigate the limitations of current-elasticity-based methods in producing accurate elastic parameter estimates in poroelastic media. The results indicate that linearly elastic reconstructions of fluid-saturated porous media at amplitudes and frequencies relevant to steady-state MRE can yield misleading effective property distributions resulting from the complex interaction between their solid and fluid phases.

Motion estimation is an essential step common to all magnetic resonance elastography (MRE) methods. For dynamic techniques, the motion is obtained from a sinusoidal fit of the image phase at multiple, uniformly spaced relative phase offsets, phi, between the motion and the motion encoding gradients (MEGs). Generally, eight values of phi sampled at the Nyquist interval pi/4 over [0, 2pi). We introduce a method, termed reduced motion encoding (RME), that reduces the number of phi required, thereby reducing the imaging time for an MRE acquisition. A frequency-domain algorithm was implemented using the discrete Fourier transform (DFT) to derive the general least-squares solution for the motion amplitude and phase given an arbitrary number of phi. A closed form representation of the condition number of the transformation matrix which is used for estimating motion was introduced to determine the sensitivity to noise for different sampling patterns of phi. Simulation results confirmed the minimum error sampling patterns suggested from the condition number maps. The minimum noise in the motion estimate is obtained when the sampled phi are essentially evenly distributed over the range [0, pi) with an interval pi/n, where n is the number of phi sampled, or alternatively with an interval 2pi/n over the range [0, 2pi) which represents the Nyquist interval. Simulations also show that the noise level decreases as n increases as expected. The decrease in noise is the largest when n is small and it becomes less significant as n increases. The algorithm also makes it possible to estimate the motion from only two values of phi, which cannot be accomplished with traditional methods because sampling at the Nyquist interval is indeterminate. Finally, noise levels in motion estimated from phantom studies and in vivo results taken with different n agreed with that predicted by simulation and condition number calculations.

Background: Magnetic resonance elastography (MRE) is an imaging technique that can noninvasively assess the shear properties of the intervertebral disc (IVD). Unlike the standard gradient recalled echo (GRE) MRE technique, a spin-echo echo-planar imaging (SE-EPI) sequence has the potential to improve imaging efficiency and patient compliance., Purpose: To validate the use of an SE-EPI sequence for MRE of the IVD compared against the standard GRE sequence., Study Type: Cross-over., Subjects: Twenty-eight healthy volunteers (15 males and 13 females, age range: 19-55)., Field Strength/sequence: 3 T; GRE, SE-EPI with breath holds (SE-EPI-BH) and SE-EPI with free breathing (SE-EPI-FB) MRE sequences., Assessment: MRE-derived shear stiffnesses were calculated via principal frequency analysis. SE-EPI derived shear stiffness and octahedral shear strain signal-to-noise ratios (OSS-SNR) were compared against those derived using the GRE sequence. The reproducibility and repeatability of SE-EPI stiffness measurements were determined. Shear stiffness was evaluated in the nucleus pulposus (NP) and annulus fibrosus (AF) regions of the disc. Scan times between sequences were compared., Statistical Tests: Linear mixed models, Bland-Altman plots, and Lin's concordance correlation coefficients (CCCs) were used with P < 0.05 considered statistically significant., Results: Good correlation was observed between shear stiffnesses derived from the SE-EPI sequences with those derived from the GRE sequence with CCC values greater than 0.73 and 0.78 for the NP and AF regions, respectively. OSS-SNR was not significantly different between GRE and SE-EPI sequences (P > 0.05). SE-EPI sequences generated highly reproducible and repeatable stiffness measurements with CCC values greater than 0.97 in the NP and AF regions and reduced scan time by at least 51% compared to GRE. SE-EPI-BH and SE-EPI-FB stiffness measurements were similar with CCC values greater than 0.98 for both regions., Data Conclusion: SE-EPI-based MRE-derived stiffnesses were highly reproducible and repeatable and correlated with current standard GRE MRE-derived stiffness estimates while reducing scan times., Level of Evidence: 1 TECHNICAL EFFICACY STAGE: 1., (© 2022 International Society for Magnetic Resonance in Medicine.)

Purpose: Magnetic resonance elastography (MRE) maps the viscoelastic properties of soft tissues for diagnostic purposes. However, different MRE inversion methods yield different results, which hinder comparison of values, standardization, and establishment of quantitative MRE markers. Here, we introduce an expandable, open-access, webserver-based platform that offers multiple inversion techniques for multifrequency, 3D MRE data., Methods: The platform comprises a data repository and standard MRE inversion methods including local frequency estimation (LFE), direct-inversion based multifrequency dual elasto-visco (MDEV) inversion, and wavenumber-based (k-) MDEV. The use of the platform is demonstrated in phantom data and in vivo multifrequency MRE data of the kidneys and brains of healthy volunteers., Results: Detailed maps of stiffness were generated by all inversion methods showing similar detail of anatomy. Specifically, the inner renal cortex had higher shear wave speed (SWS) than renal medulla and outer cortex without lateral differences. k-MDEV yielded higher SWS values than MDEV or LFE (full kidney/brain k-MDEV: 2.71 ± 0.19/1.45 ± 0.14 m/s, MDEV: 2.14 ± 0.16/0.99 ± 0.11 m/s, LFE: 2.12 ± 0.15/0.89 ± 0.06 m/s)., Conclusion: The freely accessible platform supports the comparison of MRE results obtained with different inversion methods, filter thresholds, or excitation frequencies, promoting reproducibility in MRE across community-developed methods., (© 2022 The Authors. Magnetic Resonance in Medicine published by Wiley Periodicals LLC on behalf of International Society for Magnetic Resonance in Medicine.)

Objectives: To evaluate the potential of MR elastography (MRE)-based shear strain mapping to noninvasively predict the presence of microvascular invasion (MVI) in hepatocellular carcinoma (HCC)., Methods: Fifty-nine histopathology-proven HCC patients with conventional 60-Hz MRE examinations (+/-MVI, n = 34/25) were enrolled retrospectively between December 2016 and October 2019, with one subgroup comprising 29/59 patients (+/-MVI, n = 16/13) who also underwent 40- and 30-Hz MRE examinations. Octahedral shear strain (OSS) maps were calculated, and the percentage of peritumoral interface length with low shear strain (i.e., a low-shear-strain length, pLSL, %) was recorded. For OSS-pLSL, differences between the MVI (+) and MVI (-) groups and diagnostic performance at different MRE frequencies were analyzed using the Mann-Whitney test and area under the receiver operating characteristic curve (AUC), respectively., Results: The peritumor OSS-pLSL was significantly higher in the MVI (+) group than in the MVI (-) group at the three frequencies (all p < 0.01). The AUC of peritumor OSS-pLSL for predicting MVI was good/excellent in all frequency groups (60-Hz: 0.73 (n = 59)/0.80 (n = 29); 40-Hz: 0.84; 30-Hz: 0.90). On further analysis of the 29 cases with all frequencies, the AUCs were not significantly different. As the frequency decreased from 60-Hz, the specificity of OSS increased at 40-Hz (53.8-61.5%) and further increased at 30-Hz (53.8-76.9%), and the sensitivity remained high at lower frequencies (100.0-93.8%) (all p > 0.05)., Conclusions: MRE-based shear strain mapping is a promising technique for noninvasively predicting the presence of MVI in patients with HCC, and the most recommended frequency for OSS is 30-Hz., Key Points: • MR elastography (MRE)-based shear strain mapping has the potential to predict the presence of microvascular invasion (MVI) in hepatocellular carcinoma preoperatively. • The low interface shear strain identified at tumor-liver boundaries was highly correlated with the presence of MVI., (© 2022. The Author(s), under exclusive licence to European Society of Radiology.)

Our goal is to design, test and verify an electromagnetic actuator for brain magnetic resonance elastography (MRE). We proposed a grappler-shaped design that can transmit stable vibrations into the brain. To validate its performance, simulations were carried out to ensure the electromagnetic field generated by the actuator did not interfere with the B 0 field. The actuation vibration spectrum was analyzed to verify the actuation accuracy. Phantom and volunteer experiments were carried out to evaluate the performance of the actuator. Simulation of the magnetic field showed that the proposed actuator has a fringe field of less than 3 G in the imaging region. The phantom experiments showed that the proposed actuator did not interfere with the routine imaging sequences. The measured vibration spectra demonstrated that the frequency offset was about one third that of a pneumatic device and the transmission efficiency was three times higher. The shear moduli estimated from brain MRE were consistent with those from the literature. The actuation frequency of the proposed actuator has less frequency offset and off-center frequency components compared with the pneumatic counterpart. The whole actuator weighted only 980 g. The actuator can carry out multifrequency MRE on the brain with high accuracy. It is easy to use, comfortable for the patient and portable., (© 2021 John Wiley & Sons, Ltd.)

Objectives: To develop an objective quantitative method to characterize and visualize meningioma-brain adhesion using MR elastography (MRE)-based slip interface imaging (SII)., Methods: This retrospective study included 47 meningiomas (training dataset: n = 35; testing dataset: n = 12) with MRE/SII examinations. Normalized octahedral shear strain (NOSS) values were calculated from the acquired MRE displacement data. The change in NOSS at the tumor boundary (ΔNOSS bdy ) was computed, from which a 3D ΔNOSS bdy map of the tumor surface was created and the probability distribution of ΔNOSS bdy over the entire tumor surface was calculated. Statistical features were calculated from the probability histogram. After eliminating highly correlated features, the capability of the remaining feature for tumor adhesion classification was assessed using a one-way ANOVA and ROC analysis., Results: The magnitude and location of the tumor adhesion can be visualized by the reconstructed 3D ΔNOSS bdy surface map. The entropy of the ΔNOSS bdy histogram was significantly different between adherent tumors and partially/completely non-adherent tumors in both the training (AUC: 0.971) and testing datasets (AUC: 0.900). Based on the cutoff values obtained from the training set, the ΔNOSS bdy entropy in the testing dataset yielded an accuracy of 0.83 for distinguishing adherent versus partially/non-adherent tumors, and 0.67 for distinguishing non-adherent versus completely/partially adherent tumors., Conclusions: SII-derived ΔNOSS bdy values are useful for quantification and classification of meningioma-brain adhesion. The reconstructed 3D ΔNOSS bdy surface map presents the state and location of tumor adhesion in a "clinician-friendly" manner, and can identify meningiomas with a high risk of adhesion to adjacent brain parenchyma., Key Points: • MR elastography (MRE)-based slip interface imaging shows promise as an objective tool to preoperatively discriminate meningiomas with a high risk of intraoperative adhesion. • Measurement of the change of shear strain at meningioma boundaries can provide quantitative metrics depicting the state of adhesion at the tumor-brain interface. • The surface map of tumor adhesion shows promise in assisting precise adhesion localization, using a comprehensible, "clinician-friendly" 3D visualization., (© 2021. European Society of Radiology.)

Background: Changes in brain stiffness can be an important biomarker for neurological disease. Magnetic resonance elastography (MRE) quantifies tissue stiffness, but the results vary between acquisition and reconstruction methods., Purpose: To measure MRE repeatability and estimate the effect of different reconstruction methods and varying data quality on estimated brain stiffness., Study Type: Prospective., Subjects: Fifteen healthy subjects., Field Strength/sequence: 3T MRI, gradient-echo elastography sequence with a 50 Hz vibration frequency., Assessment: Imaging was performed twice in each subject. Images were reconstructed using a curl-based and a finite-element-model (FEM)-based method. Stiffness was measured in the whole brain, in white matter, and in four cortical and four deep gray matter regions. Repeatability coefficients (RC), intraclass correlation coefficients (ICC), and coefficients of variation (CV) were calculated. MRE data quality was quantified by the ratio between shear waves and compressional waves., Statistical Tests: Median values with range are presented. Reconstruction methods were compared using paired Wilcoxon signed-rank tests, and Spearman's rank correlation was calculated between MRE data quality and stiffness. Holm-Bonferroni corrections were employed to adjust for multiple comparisons., Results: In the whole brain, CV was 4.3% and 3.8% for the curl and the FEM reconstruction, respectively, with 4.0-12.8% for subregions. Whole-brain ICC was 0.60-0.74, ranging from 0.20 to 0.89 in different regions. RC for the whole brain was 0.14 kPa and 0.17 kPa for the curl and FEM methods, respectively. FEM reconstruction resulted in 39% higher stiffness than the curl reconstruction (P < 0.05). MRE data quality, defined as shear-compression wave ratio, was higher in peripheral regions than in central regions of the brain (P < 0.05). No significant correlations were observed between MRE data quality and stiffness estimates., Data Conclusion: MRE of the human brain is a robust technique in terms of repeatability. Caution is warranted when comparing stiffness values obtained with different techniques., Level of Evidence: 1 TECHNICAL EFFICACY STAGE: 1., (© 2021 The Authors. Journal of Magnetic Resonance Imaging published by Wiley Periodicals LLC. on behalf of International Society for Magnetic Resonance in Medicine.)

Objectives: To demonstrate the feasibility of 2-dimensional brain ultrasound shear wave elastography (SWE) and to define the average elasticity values of the gray and white matter in term neonates., Methods: This work was a prospective observational single-center study including 55 healthy term neonates consecutively recruited in the maternity ward between the second and third postnatal days. All were successfully evaluated with a cerebral SWE examination performed with a multifrequency 4-9-MHz transducer. Bilateral sagittal planes of the thalamus and corona radiata were used to measure stiffness using a quantitative SWE method. Several elastograms with 5 to 15 nonoverlapping areas were obtained from the 2 different anatomic locations. The 5 most central measurements were averaged as representative values., Results: The 55 neonates ranged from 37 to 40 weeks' gestation. The estimated mean velocity values of the thalamus (1.17 m/s; 95% confidence interval, 1.13, 1.22 m/s) and corona radiata (1.60 m/s; 95% confidence interval, 1.57, 1.64 m/s) were statistically different (P < .001). There was no significant influence of laterality, gestational age, cephalic perimeter, sex, length, or type of delivery on the stiffness measurements., Conclusions: Brain ultrasound SWE is feasible and allows measurements of neonatal brain elasticity. The elasticity of the thalamus and corona radiata at the frontal white matter in healthy term neonates is different. The knowledge of normal SWE ranges in term neonates allows comparative studies under pathologic conditions., (© 2020 American Institute of Ultrasound in Medicine.)

Objectives: This study exploited finite-element modeling (FEM) to simulate breast tissue multicompression during ultrasound elastography to classify breast tumors based on their nonlinear biomechanical properties., Methods: Numeric simulations were first calculated by using 3-dimensional (3D) virtual models with an assumed tumor's geometric dimensions but with actual material properties to test and validate the FEM. Further numeric simulations were used to construct 3D models based on in vivo experimental data to verify our models. The models were designed for each individual in vivo case, emphasizing the geometry, position, and biomechanical properties of the breast tissue. At different compression levels, tissue strains were analyzed between the tumors and the background normal tissues to explore their nonlinearity and classify the tumor type. Tumor classification parameters were deduced by using a power-law relationship between the applied compressive forces and strain differences., Results: Classification parameters were compared between benign and malignant tumors, for which they were found to be statistically significant in classifying the tumor types (P < .05) by both the validation and verification of FEM. We compared the classification parameters between the in vivo and FEM classifications, for which they were found to be strongly correlated (R = 0.875; P < .001), with no statistical differences between their outcomes (P = .909)., Conclusions: Good agreement between the model outcomes and the in vivo diagnostics was reported. The implemented models were validated and verified. The introduced 3D modeling method may augment elastographic methods to preliminary classify breast tumors at an early stage., (© 2020 American Institute of Ultrasound in Medicine.)

Background: Malignant tumors are associated with increased tissue rigidity, which can be an indicator of tumor progression. MR elastography (MRE) has the potential to study the variations of tumor mechanical properties. ex vivo studies have shown the ability of MRE to assess increase of mechanical properties; nevertheless, it has not yet been observed in vivo., Purpose: To propose a method to assess the increase in mechanical properties of tumors in vivo under static external compression using MRE., Study Type: Prospective, experimental study., Animal Model: Forty-six SCID mice with subcutaneous tumor implantation (patient-derived hepatocellular carcinoma xenografts, Model 1, n = 13, and Model 2, n = 33)., Field Strength/sequence: 7.0T; a spin echo sequence was used for anatomical images and a modified spin echo sequence for elastography acquisitions with a vibration frequency of 600 Hz., Assessment: An inflatable balloon was placed on the abdomen to apply a load to the tumor. MRE acquisitions were performed at the basal state and at increasing compression levels. Anatomical images were used to calculate the octahedral shear strain between the tumor at the basal strain state and each strain level. For six mice (Model 2), each static preloading scan was acquired twice consecutively without moving the mouse to evaluate repeatability. Statistical Tests: The Bland-Altman method was used to assess repeatability. Correlations between tumor stiffness and deformation were evaluated with Pearson correlation coefficients., Results: For stiffness (G*), a good repeatability was obtained between the acquisitions; the limits of agreement of the Bland-Altman test were [-10.17%; 11.49%] with an absolute bias of 0.66%. A significant correlation between tumor stiffness and deformation was observed for both models (Model 1: r = 0.57, P < 0.0001 and Model 2: r = 0.31, P < 0.0001)., Data Conclusion: We establish that tumor mechanical properties can increase under mechanical compression. This increase can effectively be monitored using a proposed MRE setup., Level of Evidence: 2 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2019;50:1982-1989., (© 2019 International Society for Magnetic Resonance in Medicine.)

Magnetic resonance elastography is a relatively new, rapidly evolving quantitative magnetic resonance imaging technique which can be used for mapping the viscoelastic mechanical properties of soft tissues. MR elastography measurements are akin to manual palpation but with the advantages of both being quantitative and being useful for regions which are not available for palpation, such as the human brain. MR elastography is noninvasive, well tolerated, and complements standard radiological and histopathological studies by providing in vivo measurements that reflect tissue microstructural integrity. While brain MR elastography studies in adults are becoming frequent, published studies on the utility of MR elastography in children are sparse. In this review, we have summarized the major scientific principles and recent clinical applications of brain MR elastography in diagnostic neuroscience and discuss avenues for impact in assessing the pediatric brain. [ABSTRACT FROM AUTHOR]

Objective: Magnetic resonance elastography (MRE) requires substantial data processing based on phase image reconstruction, wave enhancement, and inverse problem solving. The objective of this study is to propose a new, fast MRE method based on MR raw data processing, particularly adapted to applications requiring fast MRE measurement or high elastogram update rate. Materials and methods: The proposed method allows measuring tissue elasticity directly from raw data without prior phase image reconstruction and without phase unwrapping. Experimental feasibility is assessed both in a gelatin phantom and in the liver of a porcine model in vivo. Elastograms are reconstructed with the raw MRE method and compared to those obtained using conventional MRE. In a third experiment, changes in elasticity are monitored in real-time in a gelatin phantom during its solidification by using both conventional MRE and raw MRE. Results: The raw MRE method shows promising results by providing similar elasticity values to the ones obtained with conventional MRE methods while decreasing the number of processing steps and circumventing the delicate step of phase unwrapping. Limitations of the proposed method are the influence of the magnitude on the elastogram and the requirement for a minimum number of phase offsets. Conclusion: This study demonstrates the feasibility of directly reconstructing elastograms from raw data. [ABSTRACT FROM AUTHOR]

Tissue-mimicking phantoms with well-defined properties can help in identifying the potential weaknesses in medical imaging systems. Among the imaging systems, magnetic resonance elastography is a new noninvasive technique used to quantify the shear modulus of biological tissues, and therefore has shown promise in studying liver and brain pathologies. Polyvinyl alcohol (PVA) cryogel prepared by the freeze-thaw technique is a potential candidate for mimicking the mechanical properties of soft tissues and has been extensively used as a phantom material. However, large PVA cryogels suffer from variations in properties, partly due to the low thermal conductivity of PVA solution. The loss of homogeneity in cryogel phantoms is also attributed to inhomogeneous thawing rates during the freeze-thaw cycle. We have used a modified freeze-thaw process that imposes multiple isotherms so as to enhance the homogeneity of the produced cryogels. In addition, we have developed a finite-element modeling tool (a virtual controller) to optimize the temperature profile during the freeze-thaw cycle. Our experimental validations demonstrated the potential of the virtual controller in predicting the optimal tempera-ture profile for the freeze-thaw process (phantom diameters: 60 and 100 mm). A robust simulation framework can fill the gap in the scientific literature with regard to phantom design for medical imaging and will help to reduce phantom development time and cost. [ABSTRACT FROM AUTHOR]