Measurement of in vivo local shear modulus using MR elastography multiple-phase patchwork offsets

Design, Construction, and Implementation of a Magnetic Resonance Elastography Actuator for Research Purposes

Magnetic resonance elastography (MRE) is a technique for determining the mechanical response of soft materials using applied harmonic deformation of the material and a motion-sensitive magnetic resonance imaging sequence. This technique can elucidate significant information about the health and development of human tissue such as liver and brain, and has been used on phantom models (e.g., agar, silicone) to determine their suitability for use as a mechanical surrogate for human tissues in experimental models. The applied harmonic deformation used in MRE is generated by an actuator, transmitted in bursts of a specified duration, and synchronized with the magnetic resonance signal excitation. These actuators are most often a pneumatic design (common for human tissues or phantoms) or a piezoelectric design (common for small animal tissues or phantoms). Here, we describe how to design and assemble both a pneumatic and a piezoelectric MRE actuator for research purposes. For each of these actuator types, we discuss displacement requirements, end-effector options and challenges, electronics and electronic-driving requirements and considerations, and full MRE implementation. We also discuss how to choose the actuator type, size, and power based on the intended material and use. © 2022 Wiley Periodicals LLC. Basic Protocol 1: Design, construction, and implementation of a convertible pneumatic MRE actuator for use with tissues and phantom models Basic Protocol 2: Design, construction, and implementation of a piezoelectric MRE actuator for localized excitation in phantom models.

Quantitative analysis of vibration waves based on Fourier transform in magnetic resonance elastography

We developed a novel magnetic resonance elastography (MRE) analysis method based on Fourier transform to assess the responsive characteristics for different tissue stiffness and degree of transmission of the vibration wave emanating from a passive driver during MRE. A phantom tissue study was conducted with an MRE sequence and vibration wave system using a clinical MR scanner. The phantom tissue consisted of two layers of agar: 0.75 wt% and 1.0 wt%. Phase-unwrapped images derived from acquired MRE phase images were used to generate a phase profile curve, with a line plotted for the phase-unwrapped images. Fourier transform was performed, and the peak value of the power spectrum was derived. The damping rate/ratio was calculated using the Hilbert transform of the phase profile. We found that the mean shear stiffness value of 1.0 wt% agar was higher than that of 0.75 wt% agar. The responsive frequency of the 0.75 wt% agar layer showed a wider range and the damping rate of the signal showed a higher value than the respective values of the 1.0 wt% agar layer. In conclusion, Fourier transform analysis of MRE enabled us to obtain more detailed information of the tissue characteristics and vibration-wave conditions.

Narrowband Shear Wave Generation Using Sinusoidally Modulated Acoustic Radiation Force

Most transient ultrasound elastography methods use high-intensity ultrasound "push" pulses that generate a shear wave with a wide frequency spectrum. However, it is difficult to control how the energy of the wave is distributed within that spectrum. For this reason, the shear-wave group velocity may not match that of harmonic methods like magnetic resonance elastography (MRE). The objective of this study was to introduce a narrowband shear wave generation method produced by "push" pulses with sinusoidally modulated intensity. The method, named harmonic shear wave imaging (HSWI), successively transmits a series of push pulses with a periodic change in duration. The excited shear waves form a continuous shear wave with a known main frequency that can be controlled by the user. Push pulses are interleaved with imaging pulses so only one clinical transducer is used to generate and record the shear waves. The proposed method was compared to MRE and a transient shear wave elastography method using phantoms and in vivo measurements. It was found that HSWI produces narrowband waves with a speed that closely matches that measured by MRE. Measurement of the acoustic output parameters indicated that the acoustic intensities in HSWI are suitable for clinical applications. The ability of HSWI to generate narrowband shear waves using a single linear array transducer makes it amenable for clinical translation. HSWI can potentially use the same thresholds as MRE for diagnosis of diseases affecting the stiffness of soft tissues.

Mechanics of ultrasound elastography

Ultrasound elastography enables in vivo measurement of the mechanical properties of living soft tissues in a non-destructive and non-invasive manner and has attracted considerable interest for clinical use in recent years. Continuum mechanics plays an essential role in understanding and improving ultrasound-based elastography methods and is the main focus of this review. In particular, the mechanics theories involved in both static and dynamic elastography methods are surveyed. They may help understand the challenges in and opportunities for the practical applications of various ultrasound elastography methods to characterize the linear elastic, viscoelastic, anisotropic elastic and hyperelastic properties of both bulk and thin-walled soft materials, especially the in vivo characterization of biological soft tissues.

Model independent MRE data analysis

For the diagnosing modality called MRE (magnetic resonance elastography), the displacement vector of a wave propagating in a human tissue can be measured. The average of the local wavelength from this measured data could be an index for the diagnosing, because the local wave length becomes larger when the tissue is stiffer. By assuming that the local form of the wave is given approximately as multiple complex plane waves, we identify the real part of the complex linear phase of the strongest plane wave of this multiple complex plane waves, by first applying the FBI transform (Fourier-Bros-Iagolnitzer transform) with an appropriate size of Gaussian window and then taking the maximum of the modulus of the transform with respect to the Fourier variable. The real part of the linear phase is nothing but the real inner product of the wave vector and the position vector. Similarly the imaginary part of the linear phase describes the attenuation of the wave and it is given as a real inner product of a real vector and the position vector. This vector can also be recovered by our method. We also apply these methods to design some denoising and filtering for noisy MRE data.

Contrast detection in fluid-saturated media with magnetic resonance poroelastography

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.

Muscle performance in a soft-bodied terrestrial crawler: constitutive modelling of strain-rate dependency

Experimental data on the passive mechanical properties of the ventral interior lateral muscle of the tobacco hornworm caterpillar, Manduca sexta, are reported. The stress-deformation response of the Manduca muscle is shown to be nonlinear pseudo-elastic, capable of large deformations and subject to stress softening during initial loading cycles. The muscle passive mechanical properties also depend on multiple time-dependent processes. In particular, we show new experimental data from cyclic loading tests of an unstimulated muscle with constant maximum stretch and different, constant engineering strain rates. Then, on the basis of these data a constitutive model is derived to reproduce the main characteristics of this behaviour. In formulating the constitutive model, we consider the muscle as a complex macromolecular structure with fibrous components at numerous size scales. The model uses a phenomenological approach to account for different mechanisms by which passive force changes during applied deformation and how the muscle properties recover after unloading.

Optimized motion estimation for MRE data with reduced motion encodes

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.

Magnetic resonance elastography to observe deep areas: comparison of external vibration systems

MRE methods deform the sample using an external vibration system. We have been using a transverse driver, which generates shear waves at the object surface. One of the problems is that shear waves rapidly attenuate at the surface of tissue and do not propagate into the body. In this study, we compared the shear waves generated by transverse and longitudinal drivers. The longitudinal driver was found to induce shear waves deep inside a porcine liver phantom. These results suggest that the longitudinal driver will allow measurement of the shear modulus deep inside the body.

A magnetic resonance elastography (MRE) approach for liver investigation

Due to its key role in carbohydrate metabolisation, blood detoxification, filtering and substance storing, liver is prone to a wide variety of diseases: hepatitis, cirrhosis, cancer. As such, liver investigation is a major field of interest. We describe here a non invasive way to obtain important information on liver elasticity properties. In vitro Magnetic Resonance Elastography (MRE) experiments are realized using a normalized substance (Agar gel in different concentrations) and this technique is also demonstrated in the context of an in vivo investigation of the elastic properties of a normal control liver.

Dynamic MRS and MRI of skeletal muscle function and biomechanics

MR is a powerful technique for studying the biomechanical and functional properties of skeletal muscle in vivo in health and disease. This review focuses on 31P, 1H and 13C MR spectroscopy for assessment of the dynamics of muscle metabolism and on dynamic 1H MRI methods for non-invasive measurement of the biomechanical and functional properties of skeletal muscle. The information thus obtained ranges from the microscopic level of the metabolism of the myocyte to the macroscopic level of the contractile function of muscle complexes. The MR technology presented plays a vital role in achieving a better understanding of many basic aspects of muscle function, including the regulation of mitochondrial activity and the intricate interplay between muscle fiber organization and contractile function. In addition, these tools are increasingly being employed to establish novel diagnostic procedures as well as to monitor the effects of therapeutic and lifestyle interventions for muscle disorders that have an increasing impact in modern society.

Mechanical transient-based magnetic resonance elastography

Magnetic resonance elastography (MRE) is a technique for quantifying material properties by measuring cyclic displacements of propagating shear waves. As an alternative to dynamic harmonic wave MRE or quasi-steady-state methods, the idea of using a transient impulse for mechanical excitation is introduced. Two processing methods to calculate shear stiffness from transient data were developed. The techniques were tested in phantom studies, and the transient results were found to be comparable to the harmonic wave results. Transient wave based analysis was applied to the brains of six healthy volunteers in order to assess the method in areas of complex wave patterns and geometry. The results demonstrated the feasibility of measuring brain stiffness in vivo using a transient mechanical excitation. Transient and harmonic methods both measure white matter (approximately 12 kPa) to be stiffer than gray matter ( approximately 8 kPa). There were some anatomic differences between harmonic and transient MRE, specifically where the transient results better depicted the deeper structures of the brain.