Needle shear wave driver for magnetic resonance elastography

MR Elastography of the Abdomen: Experimental Protocols

Application of MRE for noninvasive evaluation of renal fibrosis has great potential for noninvasive assessment in patients with chronic kidney disease (CKD). CKD leads to severe complications, which require dialysis or kidney transplant and could even result in death. CKD in native kidneys and interstitial fibrosis in allograft kidneys are the two major kidney fibrotic pathologies where MRE may be clinically useful. Both these conditions can lead to extensive morbidity, mortality, and high health care costs. Currently, biopsy is the standard method for renal fibrosis staging. This method of diagnosis is painful, invasive, limited by sampling bias, exhibits inter- and intraobserver variability, requires prolonged hospitalization, poses risk of complications and significant bleeding, and could even lead to death. MRE based methods can potentially be useful to noninvasively detect, stage, and monitor renal fibrosis, reducing the need for renal biopsy. In this chapter, we describe experimental procedure and step by step instructions to run MRE along with some illustrative applications. We also includes sections on how to perform data quality check and analysis methods.This publication is based upon work from the COST Action PARENCHIMA, a community-driven network funded by the European Cooperation in Science and Technology (COST) program of the European Union, which aims to improve the reproducibility and standardization of renal MRI biomarkers.

Analytical solution for diverging elliptic shear wave in bounded and unbounded transverse isotropic viscoelastic material with nonhomogeneous inner boundary

A theoretical approach was recently introduced by Guidetti and Royston [J. Acoust. Soc. Am. 144, 2312-2323 (2018)] for the radially converging elliptic shear wave pattern in transverse isotropic materials subjected to axisymmetric excitation normal to the fiber axis at the outer boundary of the material. This approach is enabled via a transformation to an elliptic coordinate system with isotropic properties. The approach is extended to the case of diverging shear waves radiating from a cylindrical rod that is axially oscillating perpendicular to the axis of isotropy and parallel to the plane of isotropy.

Frequency-dependent shear properties of annulus fibrosus and nucleus pulposus by magnetic resonance elastography

Aging and degeneration are associated with changes in mechanical properties in the intervertebral disc, generating interest in the establishment of mechanical properties as early biomarkers for the degenerative cascade. Magnetic resonance elastography (MRE) of the intervertebral disc is usually limited to the nucleus pulposus, as the annulus fibrosus is stiffer and less hydrated. The objective of this work was to adapt high-frequency needle MRE to the characterization of the shear modulus of both the nucleus pulposus and annulus fibrosus. Bovine intervertebral discs were removed from fresh oxtails and characterized by needle MRE. The needle was inserted in the center of the disc and vibrations were generated by an amplified piezoelectric actuator. MRE acquisitions were performed on a 4.7-T small-animal MR scanner using a spin echo sequence with sinusoidal motion encoding gradients. Acquisitions were repeated over a frequency range of 1000-1800 Hz. The local frequency estimation inversion algorithm was used to compute the shear modulus. Stiffness maps allowed the visualization of the soft nucleus pulposus surrounded by the stiffer annulus fibrosus surrounded by the homogeneous gel. A significant difference in shear modulus between the nucleus pulposus and annulus fibrosus, and an increase in the shear modulus with excitation frequency, were observed, in agreement with the literature. This study demonstrates that global characterization of both the nucleus pulposus and annulus fibrosus of the intervertebral disc is possible with needle MRE using a preclinical magnetic resonance imaging (MRI) scanner. MRE can be a powerful method for the mapping of the complex properties of the intervertebral disc. The developed method could be adapted for in situ use by preserving adjacent vertebrae and puncturing the side of the intervertebral disc, thereby allowing an assessment of the contribution of osmotic pressure to the mechanical behavior of the intervertebral disc.

A 2D finite element model for shear wave propagation in biological soft tissues: Application to magnetic resonance elastography

Dynamic elastography is a virtual palpation tool that aims at investigating the mechanical response of biological soft tissues in vivo. The objective of this study is to develop a finite element model (FEM) with low computational cost for reproducing realistically wave propagation for magnetic resonance elastography in heterogeneous soft tissues. Based on the first-order shear deformation theory for moderately thick structures, this model is developed and validated through comparison with analytical formulations of wave propagating in heterogeneous, viscoelastic infinite medium. This 2D-FEM is then compared to experimental data and a 3D-FEM using a commercial software. Our FEM is a powerful promising tool for investigations of magnetic resonance elastography.

Ultrasound viscoelasticity assessment using an adaptive torsional shear wave propagation method

PURPOSE

Different approaches have been used in dynamic elastography to assess mechanical properties of biological tissues. Most techniques are based on a simple inversion based on the measurement of the shear wave speed to assess elasticity, whereas some recent strategies use more elaborated analytical or finite element method (FEM) models. In this study, a new method is proposed for the quantification of both shear storage and loss moduli of confined lesions, in the context of breast imaging, using adaptive torsional shear waves (ATSWs) generated remotely with radiation pressure.

METHODS

A FEM model was developed to solve the inverse wave propagation problem and obtain viscoelastic properties of interrogated media. The inverse problem was formulated and solved in the frequency domain and its robustness to noise and geometric constraints was evaluated. The proposed model was validated in vitro with two independent rheology methods on several homogeneous and heterogeneous breast tissue-mimicking phantoms over a broad range of frequencies (up to 400 Hz).

RESULTS

Viscoelastic properties matched benchmark rheology methods with discrepancies of 8%-38% for the shear modulus G' and 9%-67% for the loss modulus G″. The robustness study indicated good estimations of storage and loss moduli (maximum mean errors of 19% on G' and 32% on G″) for signal-to-noise ratios between 19.5 and 8.5 dB. Larger errors were noticed in the case of biases in lesion dimension and position.

CONCLUSIONS

The ATSW method revealed that it is possible to estimate the viscoelasticity of biological tissues with torsional shear waves when small biases in lesion geometry exist.

An automatic differentiation-based gradient method for inversion of the shear wave equation in magnetic resonance elastography: specific application in fibrous soft tissues

Quantitative and accurate measurement of in vivo mechanical properties using dynamic elastography has been the scope of many research efforts over the past two decades. Most of the shear-wave-based inverse approaches for magnetic resonance elastography (MRE) make the assumption of isotropic viscoelasticity. In this paper, we propose a quantitative gradient method for inversion of the shear wave equation in anisotropic media derived from a full waveform description using analytical viscoelastic Green formalism and automatic differentiation. The abilities and performances of the proposed identification method are first evaluated on numerical phantoms calculated in a transversely isotropic medium, and subsequently on experimental MRE data measured on an isotropic hydrogel phantom, on an anisotropic cryogel phantom and on an ex vivo fibrous muscle. The experiments are carried out by coupling circular shear wave profiles generated by acoustic radiation force and MRE acquisition of the wave front. Shear modulus values obtained by our MRE method are compared to those obtained by rheometry in the isotropic hydrogel phantom, and are found to be in good agreement despite non-overlapping frequency ranges. Both the cryogel and the ex vivo muscle are found to be anisotropic. Stiffness values in the longitudinal direction are found to be 1.8 times and 1.9 times higher than those in the transverse direction for the cryogel and the muscle, respectively. The proposed method shows great perspectives and substantial benefits for the in vivo quantitative investigation of complex mechanical properties in fibrous soft tissues.

Application of Elastography for the Noninvasive Assessment of Biomechanics in Engineered Biomaterials and Tissues

The elastic properties of engineered biomaterials and tissues impact their post-implantation repair potential and structural integrity, and are critical to help regulate cell fate and gene expression. The measurement of properties (e.g., stiffness or shear modulus) can be attained using elastography, which exploits noninvasive imaging modalities to provide functional information of a material indicative of the regeneration state. In this review, we outline the current leading elastography methodologies available to characterize the properties of biomaterials and tissues suitable for repair and mechanobiology research. We describe methods utilizing magnetic resonance, ultrasound, and optical coherent elastography, highlighting their potential for longitudinal monitoring of implanted materials in vivo, in addition to spatiotemporal limits of each method for probing changes in cell-laden constructs. Micro-elastography methods now allow acquisitions at length scales approaching 5-100 μm in two and three dimensions. Many of the methods introduced in this review are therefore capable of longitudinal monitoring in biomaterials and tissues approaching the cellular scale. However, critical factors such as anisotropy, heterogeneity and viscoelasity-inherent in many soft tissues-are often not fully described and therefore require further advancements and future developments.

Mechanical contrast in spectroscopic magnetomotive optical coherence elastography

The viscoelastic properties of tissues are altered during pathogenesis of numerous diseases and can therefore be a useful indicator of disease status and progression. Several elastography studies have utilized the mechanical frequency response and the resonance frequencies of tissue samples to characterize their mechanical properties. However, using the resonance frequency as a source of mechanical contrast in heterogeneous samples is complicated because it not only depends on the viscoelastic properties but also on the geometry and boundary conditions. In an elastography technique called magnetomotive optical coherence elastography (MM-OCE), the controlled movement of magnetic nanoparticles (MNPs) within the sample is used to obtain the mechanical properties. Previous demonstrations of MM-OCE have typically used point measurements in elastically homogeneous samples assuming a uniform concentration of MNPs. In this study, we evaluate the feasibility of generating MM-OCE elastograms in heterogeneous samples based on a spectroscopic approach which involves measuring the magnetomotive response at different excitation frequencies. Biological tissues and tissue-mimicking phantoms with two elastically distinct regions placed in side-by-side and bilayer configurations were used for the experiments, and finite element method simulations were used to validate the experimental results.

Magnetomotive optical coherence elastography using magnetic particles to induce mechanical waves

Magnetic particles are versatile imaging agents that have found wide spread applicability in diagnostic, therapeutic, and rheology applications. In this study, we demonstrate that mechanical waves generated by a localized inclusion of magnetic nanoparticles can be used for assessment of the tissue viscoelastic properties using magnetomotive optical coherence elastography. We show these capabilities in tissue mimicking elastic and viscoelastic phantoms and in biological tissues by measuring the shear wave speed under magnetomotive excitation. Furthermore, we demonstrate the extraction of the complex shear modulus by measuring the shear wave speed at different frequencies and fitting to a Kelvin-Voigt model.

Shear wave induced resonance elastography of venous thrombi: a proof-of-concept

Shear wave induced resonance elastography (SWIRE) is proposed for deep venous thrombosis (DVT) elasticity assessment. This new imaging technique takes advantage of properly polarized shear waves to induce resonance of a confined mechanical heterogeneity. Realistic phantoms (n = 9) of DVT total and partial clot occlusions with elasticities from 406 to 3561 Pa were built for in vitro experiments. An ex vivo study was also performed to evaluate the elasticity of two fresh porcine venous thrombi in a pig model. Transient shear waves at 45-205 Hz were generated by the vibration of a rigid plate (plane wavefront) or by a needle to simulate a radiation pressure on a line segment (cylindrical wavefront). Induced propagation of shear waves was imaged with an ultrafast ultrasound scanner and a finite element method was developed to simulate tested experimental conditions. An inverse problem was then formulated considering the first resonance frequency of the DVT inclusion. Elasticity agreements between SWIRE and a reference spectroscopy instrument (RheoSpectris) were found in vitro for total clots either in plane (r(2) = 0.989) or cylindrical (r(2) = 0.986) wavefront configurations. For total and partial clots, elasticity estimation errors were 9.0 ±4.6% and 9.3 ±11.3%, respectively. Ex vivo, the blood clot elasticity was 498 ±58 Pa within the inferior vena cava and 436 ±45 Pa in the right common iliac vein (p = 0.22). To conclude, the SWIRE technique seems feasible to quantitatively assess blood clot elasticity in the context of DVT ultrasound imaging.

Quantitative imaging methods for the development and validation of brain biomechanics models

Rapid deformation of brain tissue in response to head impact or acceleration can lead to numerous pathological changes, both immediate and delayed. Modeling and simulation hold promise for illuminating the mechanisms of traumatic brain injury (TBI) and for developing preventive devices and strategies. However, mathematical models have predictive value only if they satisfy two conditions. First, they must capture the biomechanics of the brain as both a material and a structure, including the mechanics of brain tissue and its interactions with the skull. Second, they must be validated by direct comparison with experimental data. Emerging imaging technologies and recent imaging studies provide important data for these purposes. This review describes these techniques and data, with an emphasis on magnetic resonance imaging approaches. In combination, these imaging tools promise to extend our understanding of brain biomechanics and improve our ability to study TBI in silico.

Viscoelastic properties of soft gels: comparison of magnetic resonance elastography and dynamic shear testing in the shear wave regime

Magnetic resonance elastography (MRE) is used to quantify the viscoelastic shear modulus, G*, of human and animal tissues. Previously, values of G* determined by MRE have been compared to values from mechanical tests performed at lower frequencies. In this study, a novel dynamic shear test (DST) was used to measure G* of a tissue-mimicking material at higher frequencies for direct comparison to MRE. A closed-form solution, including inertial effects, was used to extract G* values from DST data obtained between 20 and 200 Hz. MRE was performed using cylindrical 'phantoms' of the same material in an overlapping frequency range of 100-400 Hz. Axial vibrations of a central rod caused radially propagating shear waves in the phantom. Displacement fields were fit to a viscoelastic form of Navier's equation using a total least-squares approach to obtain local estimates of G*. DST estimates of the storage G' (Re[G*]) and loss modulus G″ (Im[G*]) for the tissue-mimicking material increased with frequency from 0.86 to 0.97 kPa (20-200 Hz, n = 16), while MRE estimates of G' increased from 1.06 to 1.15 kPa (100-400 Hz, n = 6). The loss factor (Im[G*]/Re[G*]) also increased with frequency for both test methods: 0.06-0.14 (20-200 Hz, DST) and 0.11-0.23 (100-400 Hz, MRE). Close agreement between MRE and DST results at overlapping frequencies indicates that G* can be locally estimated with MRE over a wide frequency range. Low signal-to-noise ratio, long shear wavelengths and boundary effects were found to increase residual fitting error, reinforcing the use of an error metric to assess confidence in local parameter estimates obtained by MRE.

Physiology and cell biology of acupuncture observed in calcium signaling activated by acoustic shear wave

This article presents a novel model of acupuncture physiology based on cellular calcium activation by an acoustic shear wave (ASW) generated by the mechanical movement of the needle. An acupuncture needle was driven by a piezoelectric transducer at 100 Hz or below, and the ASW in human calf was imaged by magnetic resonance elastography. At the cell level, the ASW activated intracellular Ca²⁺ transients and oscillations in fibroblasts and endothelial, ventricular myocytes and neuronal PC-12 cells along with frequency–amplitude tuning and memory capabilities. Monitoring in vivo mammalian experiments with ASW, enhancement of endorphin in blood plasma and blocking by Gd³⁺ were observed; and increased Ca²⁺ fluorescence in mouse hind leg muscle was imaged by two-photon microscopy. In contrast with traditional acupuncture models, the signal source is derived from the total acoustic energy. ASW signaling makes use of the anisotropy of elasticity of tissues as its waveguides for transmission and that cell activation is not based on the nervous system.

Piezoelectric actuator design for MR elastography: implementation and vibration issues

BACKGROUND

MR elastography (MRE) is an emerging technique for tumor diagnosis. MRE actuation devices require precise mechanical design and radiofrequency engineering to achieve the required mechanical vibration performance and MR compatibility.

METHOD

A method of designing a general-purpose, compact and inexpensive MRE actuator is presented. It comprises piezoelectric bimorphs arranged in a resonant structure designed to operate at its resonant frequency for maximum vibration amplitude. An analytical model was established to understand the device vibration characteristics.

RESULTS

The model-predicted performance was validated in experiments, showing its accuracy in predicting the actuator resonant frequency with an error < 4%. The device MRI compatibility was shown to cause minimal interference to a 1.5 tesla MRI scanner, with maximum signal-to-noise ratio reduction of 7.8% and generated artefact of 7.9 mm in MR images.

CONCLUSIONS

A piezoelectric MRE actuator is proposed, and its implementation, vibration issues and future work are discussed.

High-frequency mode conversion technique for stiff lesion detection with magnetic resonance elastography (MRE)

A novel imaging technique is described in which the mode conversion of longitudinal waves is used for the qualitative detection of stiff lesions within soft tissue using magnetic resonance elastography (MRE) methods. Due to the viscoelastic nature of tissue, high-frequency shear waves attenuate rapidly in soft tissues but much less in stiff tissues. By introducing minimally-attenuating longitudinal waves at a significantly high frequency into tissue, shear waves produced at interfaces by mode conversion will be detectable in stiff regions, but will be significantly attenuated and thus not detectable in the surrounding soft tissue. This contrast can be used to detect the presence of stiff tissue. The proposed technique is shown to readily depict hard regions (mimicking tumors) present in tissue-simulating phantoms and ex vivo breast tissue. In vivo feasibility is demonstrated on a patient with liver metastases in whom the tumors are readily distinguished. Preliminary evidence also suggests that quantitative stiffness measurements of stiff regions obtained with this technique are more accurate than those from conventional MRE because of the short shear wavelengths. This rapid, qualitative technique may lend itself to applications in which the localization of stiff, suspicious neoplasms is coupled with more sensitive techniques for thorough characterization.

In vivo MR elastography of the prostate gland using a transurethral actuator

Conventional approaches for MR elastography (MRE) using surface drivers have difficulty achieving sufficient shear wave propagation in the prostate gland due to attenuation. In this study we evaluate the feasibility of generating shear wave propagation in the prostate gland using a transurethral device. A novel transurethral actuator design is proposed, and the performance of this device was evaluated in gelatin phantoms and in a canine prostate gland. All MRI was performed on a 1.5T MR imager using a conventional gradient-echo MRE sequence. A piezoceramic actuator was used to vibrate the transurethral device along its length. Shear wave propagation was measured transverse and parallel to the rod at frequencies between 100 and 250 Hz in phantoms and in the prostate gland. The shear wave propagation was cylindrical, and uniform along the entire length of the rod in the gel experiments. The feasibility of transurethral MRE was demonstrated in vivo in a canine model, and shear wave propagation was observed in the prostate gland as well as along the rod. These experiments demonstrate the technical feasibility of transurethral MRE in vivo. Further development of this technique is warranted.

Magnetic resonance elastography hardware design: A survey

Magnetic resonance elastography (MRE) is an emerging technique capable of measuring the shear modulus of tissue. A suspected tumour can be identified by comparing its properties with those of tissues surrounding it; this can be achieved even in deep-lying areas as long as mechanical excitation is possible. This would allow non-invasive methods for cancer-related diagnosis in areas not accessible with conventional palpation. An actuating mechanism is required to generate the necessary tissue displacements directly on the patient in the scanner and three different approaches, in terms of actuator action and position, exist to derive stiffness measurements. However, the magnetic resonance (MR) environment places considerable constraints on the design of such devices, such as the possibility of mutual interference between electrical components, the scanner field, and radio frequency pulses, and the physical space restrictions of the scanner bore. This paper presents a review of the current solutions that have been developed for MRE devices giving particular consideration to the design criteria including the required vibration frequency and amplitude in different applications, the issue of MR compatibility, actuation principles, design complexity, and scanner synchronization issues. The future challenges in this field are also described.

In vivo tumor detection on rabbit with biopsy needle as MRE driver

Magnetic Resonance Elastography (MRE) is a phase contrast imaging technique to quantitatively measure the elasticity of tissues. Typically, the oscillating driver is placed on the surface of objects to generate shear waves. When it is applied to detect tumors in deep location, the depth penetration of the wave is limited by attenuation and the biopsy procedure has to be performed separately. In this study, we describe a method using biopsy needle as the MRE driver to produce shear waves in tissues. We made comparison between the MRE acquisitions obtained with biopsy needle and surface drivers. Because the well-defined propagation wave pattern reduces the error in wavelength calculation, the acquisitions of biopsy needle driver shows better homogeneity in stiffness map. We also performed the experiment with the biopsy needle for in vivo tumor detection in rabbits. This study demonstrates that the biopsy needle driver is more effective than the surface driver for accurately measuring the stiffness and location of the tumor.

Biopsy needle as MRE driver for tumor detection

Magnetic Resonance Elastography (MRE) is a phase contrast imaging technique to quantitatively measure the elasticity of tissues. Typically, an oscillating driver is placed on the surface to generate the shear waves. The depth penetration of the wave is limited by attenuation and the biopsy procedure has to be done separately. In this study, we use a biopsy needle as the driver to detect the 15% porcine gel inclusion in a 10% porcine gel phantom which simulates a tumor in tissues. We also perform the experiment with the biopsy needle for in-vivo tumor detection in rabbits. It is shown that the biopsy needle driver can accurately measure the stiffness and location of the tumor.

Quantitative assessment of hepatic fibrosis in an animal model with magnetic resonance elastography

Chronic liver disease is a world-wide problem that causes progressive hepatic fibrosis as a hallmark of progressive injury. At present, the gold standard for diagnosing hepatic fibrosis is liver biopsy, which is an invasive method with many limitations, including questionable accuracy and risks of complications. MR elastography (MRE), a phase-contrast MRI technique for quantitatively assessing the mechanical properties of soft tissues, is a potential noninvasive diagnostic method to assess hepatic fibrosis. In this work, MRE was evaluated as a quantitative method to assess the in vivo mechanical properties of the liver tissues in a knockout animal model of liver fibrosis. This work demonstrates that the shear stiffness of liver tissue increases systematically with the extent of hepatic fibrosis, as measured by histology. A linear correlation between liver stiffness and fibrosis extent was well-defined in this animal model. An additional finding of the study was that fat infiltration, commonly present in chronic liver disease, does not significantly correlate with liver stiffness at each fibrosis stage and thus does not appear to interfere with the ability of MRE to assess fibrosis extent. In conclusion, MRE has the potential not only for assessing liver stiffness, but also for monitoring potential therapies for hepatic fibrosis.