Characterization of a hyper-viscoelastic phantom mimicking biological soft tissue using an abdominal pneumatic driver with magnetic resonance elastography (MRE)

Evaluation of the Reproducibility of MR Elastography Measurements of the Lumbar Back Muscles

BACKGROUND

MR elastography (MRE) may provide quantitative imaging biomarkers of lumbar back muscles (LBMs), complementing MRI in spinal diseases by assessing muscle mechanical properties. However, reproducibility analyses for MRE of LBM are lacking.

PURPOSE

To assess technical failure, within-day and inter-day reproducibility, robustness with the excitation source positioning, and inter-observer agreement of MRE of muscles.

STUDY TYPE

Prospective.

SUBJECTS

Seventeen healthy subjects (mean age 28 ± 4 years; 11 females).

FIELD STRENGTH/SEQUENCE

1.5 T, gradient-echo MRE, T1-weighted turbo spin echo.

ASSESSMENT

The pneumatic driver was centered at L3 level. Four MRE were performed during two visits, 2-4 weeks apart, each consisting of two MRE with less than 10 minutes inter-scan interval. At Visit 1, after the first MRE, the coil and driver were removed, then reinstalled. The MRE was repeated. At Visit 2, following the first MRE, only the driver was moved down 5 cm. The MRE was repeated. Two radiologists segmented the multifidus and erector spinae muscles.

STATISTICAL TESTS

Paired t-test, analysis of variance, intraclass correlation coefficients (ICCs). P-values <0.05 were considered statistically significant.

RESULTS

Mean stiffness of LBM ranged from 1.44 to 1.60 kPa. Mean technical failure rate was 2.5%. Inter-observer agreement was excellent (ICC ranging from 0.82 [0.64-0.96] to 0.99 [0.98-0.99] in the multifidus, and from 0.85 [0.69-0.92] to 0.99 [0.97-0.99] in the erector spinae muscles). Within-day reproducibility was fair in the multifidus (ICC: 0.53 [0.47-0.77]) and good in the erector spinae muscles (ICC: 0.74 [0.48-0.88]). Reproducibility after moving the driver was excellent in both multifidus (ICC: 0.85 [0.69-0.93]) and erector spinae muscles (ICC: 0.84 [0.67-0.92]). Inter-day reproducibility was excellent in the multifidus (ICC: 0.76 [0.48-0.89]) and poor in the erector spinae muscles (ICC: 0.23 [-0.61 to 0.63]).

DATA CONCLUSION

MRE of LBM provides measurements of stiffness with fair to excellent reproducibility and excellent inter-observer agreement. However, inter-day reproducibility in the multifidus muscles indicated that the herein used MRE protocol may not be optimal for this muscle.

EVIDENCE LEVEL

2 TECHNICAL EFFICACY: Stage 1.

Impact of axisymmetric deformation on MR elastography of a nonlinear tissue-mimicking material and implications in peri-tumour stiffness quantification

Solid tumour growth is often associated with the accumulation of mechanical stresses acting on the surrounding host tissue. Due to tissue nonlinearity, the shear modulus of the peri-tumoural region inherits a signature from the tumour expansion which depends on multiple factors, including the soft tissue constitutive behaviour and its stress/strain state. Shear waves used in MR-elastography (MRE) sense the apparent change in shear modulus along their propagation direction, thereby probing the anisotropic stiffness field around the tumour. We developed an analytical framework for a heterogeneous shear modulus distribution using a thick-shelled sphere approximation of the tumour and soft tissue ensemble. A hyperelastic material (plastisol) was identified to validate the proposed theory in a phantom setting. A balloon-catheter connected to a pressure sensor was used to replicate the stress generated from tumour pressure and growth while MRE data were acquired. The shear modulus anisotropy retrieved from the reconstructed elastography data confirmed the analytically predicted patterns at various levels of inflation. An alternative measure, combining the generated deformation and the local wave direction and independent of the reconstruction strategy, was also proposed to correlate the analytical findings with the stretch probed by the waves. Overall, this work demonstrates that MRE in combination with non-linear mechanics, is able to identify the apparent shear modulus variation arising from the strain generated by a growth within tissue, such as an idealised model of tumour. Investigation in real tissue represents the next step to further investigate the implications of endogenous forces in tissue characterisation through MRE.

Quantification of iris elasticity using acoustic radiation force optical coherence elastography

Careful quantification of the changes in biomechanical properties of the iris can offer insight into the pathophysiology of some ocular diseases. However, to date there has not been much information available regarding this subject because clinical detection for iris elasticity remains challenging. To overcome this limitation, we explore, for the first time to our knowledge, the potential of measuring iris elasticity using acoustic radiation force optical coherence elastography (ARF-OCE). The resulting images and shear wave propagation, as well as the corresponding shear modulus and Young's modulus from ex vivo and in vivo rabbit models confirmed the feasibility of this method. With features of noninvasive imaging, micrometer-scale resolution, high acquisition speed and real-time processing, ARF-OCE is a promising method for reconstruction of iris elasticity and may have great potential to be applied in clinical ophthalmology with further refinement.

Ex-vivo quantification of ovine pia arachnoid complex biomechanical properties under uniaxial tension

Background

The pia arachnoid complex (PAC) is a cerebrospinal fluid-filled tissue conglomerate that surrounds the brain and spinal cord. Pia mater adheres directly to the surface of the brain while the arachnoid mater adheres to the deep surface of the dura mater. Collagen fibers, known as subarachnoid trabeculae (SAT) fibers, and microvascular structure lie intermediately to the pia and arachnoid meninges. Due to its structural role, alterations to the biomechanical properties of the PAC may change surface stress loading in traumatic brain injury (TBI) caused by sub-concussive hits. The aim of this study was to quantify the mechanical and morphological properties of ovine PAC.

Methods

Ovine brain samples (n = 10) were removed from the skull and tissue was harvested within 30 min post-mortem. To access the PAC, ovine skulls were split medially from the occipital region down the nasal bone on the superior and inferior aspects of the skull. A template was used to remove arachnoid samples from the left and right sides of the frontal and occipital regions of the brain. 10 ex-vivo samples were tested with uniaxial tension at 2 mm s⁻¹, average strain rate of 0.59 s⁻¹, until failure at < 5 h post extraction. The force and displacement data were acquired at 100 Hz. PAC tissue collagen fiber microstructure was characterized using second-harmonic generation (SHG) imaging on a subset of n = 4 stained tissue samples. To differentiate transverse blood vessels from SAT by visualization of cell nuclei and endothelial cells, samples were stained with DAPI and anti-von Willebrand Factor, respectively. The Mooney-Rivlin model for average stress–strain curve fit was used to model PAC material properties.

Results

The elastic modulus, ultimate stress, and ultimate strain were found to be 7.7 ± 3.0, 2.7 ± 0.76 MPa, and 0.60 ± 0.13, respectively. No statistical significance was found across brain dissection locations in terms of biomechanical properties. SHG images were post-processed to obtain average SAT fiber intersection density, concentration, porosity, tortuosity, segment length, orientation, radial counts, and diameter as 0.23, 26.14, 73.86%, 1.07 ± 0.28, 17.33 ± 15.25 µm, 84.66 ± 49.18°, 8.15%, 3.46 ± 1.62 µm, respectively.

Conclusion

For the sizes, strain, and strain rates tested, our results suggest that ovine PAC mechanical behavior is isotropic, and that the Mooney-Rivlin model is an appropriate curve-fitting constitutive equation for obtaining material parameters of PAC tissues.

Reliable preparation of agarose phantoms for use in quantitative magnetic resonance elastography

Agarose phantoms are one type of phantom commonly used in developing in vivo brain magnetic resonance elastography (MRE) sequences because they are inexpensive and easy to work with, store, and dispose of; however, protocols for creating agarose phantoms are non-standardized and often result in inconsistent phantoms with significant variability in mechanical properties. Many magnetic resonance imaging (MRI) and ultrasound studies use phantoms, but often these phantoms are not tailored for desired mechanical properties and as such are too stiff or not mechanically consistent enough to be used in MRE. In this work, we conducted a systematic study of agarose phantom creation parameters to identify those factors that are most conducive to producing mechanically consistent agarose phantoms for MRE research. We found that cooling rate and liquid temperature affected phantom homogeneity. Phantom stiffness is affected by agar concentration (quadratically), by final liquid temperature and salt content in phantoms, and by the interaction of these two metrics each with stir rate. We captured and quantified the implied relationships with a regression model that can be used to estimate stiffness of resulting phantoms. Additionally, we characterized repeatability, stability over time, impact on MR signal parameters, and differences in agar gel microstructure. This protocol and regression model should prove beneficial in future MRE development studies that use phantoms to determine stiffness measurement accuracy.

3D multicellular model of shock wave-cell interaction

Understanding the interaction between shock waves and tissue is critical for advancing the use of shock waves for medical applications, such as cancer therapy. This work aims to study shock wave-cell interaction in a more realistic environment, relevant to in vitro and in vivo studies, by using 3D computational models of healthy and cancerous cells. The results indicate that for a single cell embedded in an extracellular environment, the cellular geometry does not influence significantly the membrane strain but does influence the von Mises stress. On the contrary, the presence of neighbouring cells has a strong effect on the cell response, by increasing fourfold both quantities. The membrane strain response of a cell converges with more than three neighbouring cell layers, indicating that a cluster of four layers of cells is sufficient to model the membrane strain in a large domain of tissue. However, a full 3D tissue model is needed if the stress evaluation is of main interest. A tumour mimicking multicellular spheroid model is also proposed to study mutual interaction between healthy and cancer cells and shows that cancer cells can be specifically targeted in an early stage tumour-mimicking environment.

STATEMENT OF SIGNIFICANCE

This work presents 3D computational models of shock-wave/cell interaction in a biophysically realistic environment using real cell morphology in tissue-mimicking phantoms and multicellular spheroids. Results show that cell morphology does not strongly influence the membrane strain but influences the von Mises stress. While the presence of neighbouring cells significantly increases the cell response, four cell layers are enough to capture the membrane strain change in tissue. However, a full tissue model is necessary if accurate stress analysis is needed. The work also shows that cancer cells can be specifically targeted in early stage tumour mimicking environment. This work is a step towards realistic modelling of shock-wave/cell interactions in tissues and provides insight on the use of 3D models for different scenarios.

Phantom evaluations of nonlinear inversion MR elastography

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.

Comparison of viscoelastic property characterization of plastisol phantoms with magnetic resonance elastography and high-frequency rheometry

This study aims at evaluating Magnetic Resonance Elastography (MRE) as a reliable technique for the characterization of viscoelastic properties of soft tissues. Three phantoms with different concentrations of plastisol and softener were prepared in order to mechanically mimic a broad panel of healthy and pathological soft tissues. Once placed in a MRI device, each sample was excited by a homemade external driver, inducing shear waves within the medium. The storage (G') and loss (G") moduli of each phantom were then reconstructed from MRE acquisitions over a frequency range from 300 to 1,000 Hz, by applying a 2D Helmholtz inversion algorithm. At the same time, mechanical tests were performed on four samples of each phantom with a High-Frequency piezo-Rheometer (HFR) over an overlapping frequency range (from 160 to 630 Hz) with the same test conditions (temperature, ageing). The comparison between both techniques shows a good agreement in the measurement of the storage and loss moduli, underlying the capability of MRE to noninvasively assess the complex shear modulus G* of a medium and its interest for investigating the viscoelastic properties of living tissues. Moreover, the phantoms with varying concentrations of plastisol used in this study show interesting rheological properties, which make them good candidates to simulate the broad variety of viscoelastic behaviors of healthy and pathological soft tissues.

Development of a phantom mimicking the functional and structural behaviors of the thigh muscles characterized with magnetic resonance elastography technique

Magnetic Resonance Elastography (MRE) is a non invasive technique based on the propagation of shear waves in soft tissues providing the quantification of the mechanical properties [1]. MRE was successfully applied to healthy and pathological muscles. However, the MRE muscle methods must be further improved to characterize the deep muscles. A way will be to develop phantom mimicking the muscle behavior in order to set up new MRE protocol. Thus, the purpose of this study is to create a phantom composed of a similar skeletal muscle architecture (fiber, aponorosis) and equivalent elastic properties as a function of the muscle state (passive or active). Two homogeneous phantoms were manufactured with different concentrations of plastisol to simulate the elastic properties in relaxed (50% of plastisol) and contracted (70% of plastisol) muscle conditions. Moreover, teflon tubing pipes (D = 0.9 mm) were thread in the upper part of the phantom (50%) to represent the muscle fibers and a plastic sheet (8 × 15 cm) was also included in the middle of the phantom to mimic the aponeurosis structure. Subsequently, MRE tests were performed with two different pneumatic drivers, tube and round, (f = 90Hz) to analyze the effect of the type of driver on the wave propagation. Then, the wavelength was measured from the phase images to obtain the elastic properties (shear modulus). Both phantoms revealed elastic properties which were in the same range as in vivo muscle in passive (μ(50%) = 2.40 ± 0.18 kPa ) and active (6.24 ± 0.21 kPa) states. The impact of the type of driver showed higher values (about 1.2kPa) with the tube. The analysis of the wave behavior revealed a sliding along the plastic sheet as it was observed for in vivo muscle study. The wave was also sensitive to the presence of the fibers where gaps were identified. The present study demonstrates the ability of the phantom to mimic the structural and functional properties of the muscle.

Analysis of shear wave propagation derived from MR elastography in 3D thigh skeletal muscle using subject specific finite element model

The purpose of this study was to develop a subject specific finite element model derived from MRI images to numerically analyze the MRE (magnetic resonance elastography) shear wave propagation within skeletal thigh muscles. A sagittal T2 CUBE MRI sequence was performed on the 20-cm thigh segment of a healthy male subject. Skin, adipose tissue, femoral bone and 11 muscles were manually segmented in order to have 3D smoothed solid and meshed models. These tissues were modeled with different constitutive laws. A transient modal dynamics analysis was applied to simulate the shear wave propagation within the thigh tissues. The effects of MRE experimental parameters (frequency, force) and the muscle material properties (shear modulus: C10) were analyzed through the simulated shear wave displacement within the vastus medialis muscle. The results showed a plausible range of frequencies (from 90Hz to 120 Hz), which could be used for MRE muscle protocol. The wave amplitude increased with the level of the force, revealing the importance of the boundary condition. Moreover, different shear displacement patterns were obtained as a function of the muscle mechanical properties. The present study is the first to analyze the shear wave propagation in skeletal muscles using a 3D subject specific finite element model. This study could be of great value to assist the experimenters in the set-up of MRE protocols.

Parameter identification in a generalized time-harmonic Rayleigh damping model for elastography

The identifiability of the two damping components of a Generalized Rayleigh Damping model is investigated through analysis of the continuum equilibrium equations as well as a simple spring-mass system. Generalized Rayleigh Damping provides a more diversified attenuation model than pure Viscoelasticity, with two parameters to describe attenuation effects and account for the complex damping behavior found in biological tissue. For heterogeneous Rayleigh Damped materials, there is no equivalent Viscoelastic system to describe the observed motions. For homogeneous systems, the inverse problem to determine the two Rayleigh Damping components is seen to be uniquely posed, in the sense that the inverse matrix for parameter identification is full rank, with certain conditions: when either multi-frequency data is available or when both shear and dilatational wave propagation is taken into account. For the multi-frequency case, the frequency dependency of the elastic parameters adds a level of complexity to the reconstruction problem that must be addressed for reasonable solutions. For the dilatational wave case, the accuracy of compressional wave measurement in fluid saturated soft tissues becomes an issue for qualitative parameter identification. These issues can be addressed with reasonable assumptions on the negligible damping levels of dilatational waves in soft tissue. In general, the parameters of a Generalized Rayleigh Damping model are identifiable for the elastography inverse problem, although with more complex conditions than the simpler Viscoelastic damping model. The value of this approach is the additional structural information provided by the Generalized Rayleigh Damping model, which can be linked to tissue composition as well as rheological interpretations.

Identification process based on shear wave propagation within a phantom using finite element modelling and magnetic resonance elastography

Magnetic resonance elastography (MRE), based on shear wave propagation generated by a specific driver, is a non-invasive exam performed in clinical practice to improve the liver diagnosis. The purpose was to develop a finite element (FE) identification method for the mechanical characterisation of phantom mimicking soft tissues investigated with MRE technique. Thus, a 3D FE phantom model, composed of the realistic MRE liver boundary conditions, was developed to simulate the shear wave propagation with the software ABAQUS. The assumptions of homogeneity and elasticity were applied to the FE phantom model. Different ranges of mesh size, density and Poisson's ratio were tested in order to develop the most representative FE phantom model. The simulated wave displacement was visualised with a dynamic implicit analysis. Subsequently, an identification process was performed with a cost function and an optimisation loop provided the optimal elastic properties of the phantom. The present identification process was validated on a phantom model, and the perspective will be to apply this method on abdominal tissues for the set-up of new clinical MRE protocols that could be applied for the follow-up of the effects of treatments.

Cutoff values for alcoholic liver fibrosis using magnetic resonance elastography technique

BACKGROUND

Due to the lack of cutoff values validated for specific liver diseases, the purpose of this study was to set up specific magnetic resonance elastography (MRE) cutoff values for asymptomatic liver fibrosis in alcoholic patients.

METHODS

Ninety patients underwent 3 clinical exams. The liver stiffness was measured locally with the Fibroscan, and globally through cartographies of shear modulus generated with MRE. The Fibroscan method was chosen as the gold standard to classify the fibrosis. The liver score was also obtained with the Fibrometer A, and the diagnostic performance of the methods was analyzed with receiver-operating characteristic (ROC) curves and cutoff values were calculated.

RESULTS

Spearman correlation and area under the ROC curve revealed that MRE is a better diagnostic method than the Fibrometer A, to identify various levels of fibrosis. The results showed that the Fibrometer A was adapted for severe fibrosis. The MRE cutoff values are F1:2.20 kPa, F2:2.57 kPa, F3:3.31 kPa, and F4:4 kPa and were not influenced by the glutamic oxaloacetic transaminase level. By using the ultrasound cutoff values attributed for alcoholism, 66% of patients had a similar liver fibrosis diagnosis as the MRE cutoffs. However, both imaging techniques did not provide the same distribution for minor fibrosis.

CONCLUSIONS

None of the imaging techniques (Fibroscan, MRE) could replace the gold standard of the biopsy. However, due to the risk and the unnecessary procedure for the present recruited alcoholic patients, the Fibroscan method was chosen as the reference. Since MRE is currently being used as a clinical exam, the present MRE cutoffs could aid clinicians with their diagnosis of liver fibrosis for alcoholism disease.