Viscoelastic parameters for quantifying liver fibrosis: three-dimensional multifrequency MR elastography study on thin liver rat slices

Clinical utility of 3D magnetic resonance elastography in patients with biliary obstruction

OBJECTIVES

Three-dimensional magnetic resonance elastography (3D-MRE) allows for multiparametric modeling of both elastic and viscous tissue characteristics. Our aim was to compare 3D-MRE with conventional liver shear stiffness assessment in gauging obstructive jaundice (OJ), predicting the adequacy of biliary decompression after drainage, and discriminating OJ from liver fibrosis.

METHODS

Patients with no histories of liver disease (n = 201) were studied in retrospect, grouped by bilirubin levels as no jaundice (NJ ≤ 2 mg/dL; n = 75), mild OJ (>2 mg/dL and ≤ 4 mg/dL; n = 56), and severe OJ (> 4 mg/dL; n = 70). For comparison, another 75 patients with chronic hepatitis B and C infections and histologically proven liver fibrosis were similarly analyzed. Each patient underwent spin-echo echo-planar-imaging MRE at 60 Hz with 3D wave postprocessing. Logistic regression and ordinary regression models were used to compare the 3D-MRE model with liver shear stiffness.

RESULTS

Liver shear stiffness, loss modulus, and damping ratio were incorporated into a 3D-MRE model, which significantly outperformed shear stiffness in predicting OJ severity (accuracy: 0.801 vs 0.672; p < 0.001). Both the 3D-MRE model and liver shear stiffness performed equally well in predicting the outcome of biliary drainage procedure (C-statistics: 0.852 vs 0.847; p = 0.48). The 3D-MRE model also demonstrated significantly better C-statistics than that of liver shear stiffness in discriminating mild OJ from F1-F2 liver fibrosis (0.765 vs 0.641; p = 0.005) and severe OJ from F3-F4 liver fibrosis (0.750 vs 0.635; p = 0.031).

CONCLUSIONS

3D-MRE is an innovative imaging method for gauging OJ severity, predicting the outcome of biliary drainage procedure, and discriminating OJ from liver fibrosis.

KEY POINTS

• 3D-MR elastography achieved promising results for predicting the severity of obstructive jaundice. • Advanced parameters of 3D-MR elastography demonstrated significantly better performance than that of shear stiffness of 2D-MR elastography in discriminating obstructive jaundice from liver fibrosis caused by chronic hepatitis B/C. • Both 3D-MR elastography and 2D-MR elastography were equivalent in predicting the outcome of biliary drainage procedure.

Microscopic multifrequency MR elastography for mapping viscoelasticity in zebrafish

PURPOSE

The zebrafish (Danio rerio) has become an important animal model in a wide range of biomedical research disciplines. Growing awareness of the role of biomechanical properties in tumor progression and neuronal development has led to an increasing interest in the noninvasive mapping of the viscoelastic properties of zebrafish by elastography methods applicable to bulky and nontranslucent tissues.

METHODS

Microscopic multifrequency MR elastography is introduced for mapping shear wave speed (SWS) and loss angle (φ) as markers of stiffness and viscosity of muscle, brain, and neuroblastoma tumors in postmortem zebrafish with 60 µm in-plane resolution. Experiments were performed in a 7 Tesla MR scanner at 1, 1.2, and 1.4 kHz driving frequencies.

RESULTS

Detailed zebrafish viscoelasticity maps revealed that the midbrain region (SWS = 3.1 ± 0.7 m/s, φ = 1.2 ± 0.3 radian [rad]) was stiffer and less viscous than telencephalon (SWS = 2.6 ± 0. 5 m/s, φ = 1.4 ± 0.2 rad) and optic tectum (SWS = 2.6 ± 0.5 m/s, φ = 1.3 ± 0.4 rad), whereas the cerebellum (SWS = 2.9 ± 0.6 m/s, φ = 0.9 ± 0.4 rad) was stiffer but less viscous than both (all p < .05). Overall, brain tissue (SWS = 2.9 ± 0.4 m/s, φ = 1.2 ± 0.2 rad) had similar stiffness but lower viscosity values than muscle tissue (SWS = 2.9 ± 0.5 m/s, φ = 1.4 ± 0.2 rad), whereas neuroblastoma (SWS = 2.4 ± 0.3 m/s, φ = 0.7 ± 0.1 rad, all p < .05) was the softest and least viscous tissue.

CONCLUSION

Microscopic multifrequency MR elastography-generated maps of zebrafish show many details of viscoelasticity and resolve tissue regions, of great interest in neuromechanical and oncological research and for which our study provides first reference values.

Necro-inflammatory activity grading in chronic viral hepatitis with three-dimensional multifrequency MR elastography

The purpose of this study was to assess the diagnostic value of multifrequency MR elastography for grading necro-inflammation in the liver. Fifty participants with chronic hepatitis B or C were recruited for this institutional review board-approved study. Their liver was examined with multifrequency MR elastography. The storage, shear and loss moduli, and the damping ratio were measured at 56 Hz. The multifrequency wave dispersion coefficient of the shear modulus was calculated. The measurements were compared to reference markers of necro-inflammation and fibrosis with Spearman correlations and multiple regression analysis. Diagnostic accuracy was assessed. At multiple regression analysis, necro-inflammation was the only determinant of the multifrequency dispersion coefficient, whereas fibrosis was the only determinant of the storage, loss and shear moduli. The multifrequency dispersion coefficient had the largest AUC for necro-inflammatory activity A ≥ 2 [0.84 (0.71-0.93) vs. storage modulus AUC: 0.65 (0.50-0.79), p = 0.03], whereas the storage modulus had the largest AUC for fibrosis F ≥ 2 [AUC (95% confidence intervals) 0.91 (0.79-0.98)] and cirrhosis F4 [0.97 (0.88-1.00)]. The measurement of the multifrequency dispersion coefficient at three-dimensional MR elastography has the potential to grade liver necro-inflammation in patients with chronic vial hepatitis.

Application of Magnetic Resonance Imaging in Liver Biomechanics: A Systematic Review

MRI-based biomechanical studies can provide a deep understanding of the mechanisms governing liver function, its mechanical performance but also liver diseases. In addition, comprehensive modeling of the liver can help improve liver disease treatment. Furthermore, such studies demonstrate the beginning of an engineering-level approach to how the liver disease affects material properties and liver function. Aimed at researchers in the field of MRI-based liver simulation, research articles pertinent to MRI-based liver modeling were identified, reviewed, and summarized systematically. Various MRI applications for liver biomechanics are highlighted, and the limitations of different viscoelastic models used in magnetic resonance elastography are addressed. The clinical application of the simulations and the diseases studied are also discussed. Based on the developed questionnaire, the papers' quality was assessed, and of the 46 reviewed papers, 32 papers were determined to be of high-quality. Due to the lack of the suitable material models for different liver diseases studied by magnetic resonance elastography, researchers may consider the effect of liver diseases on constitutive models. In the future, research groups may incorporate various aspects of machine learning (ML) into constitutive models and MRI data extraction to further refine the study methodology. Moreover, researchers should strive for further reproducibility and rigorous model validation and verification.

Viscoelasticity Imaging of Biological Tissues and Single Cells Using Shear Wave Propagation

Changes in biomechanical properties of biological soft tissues are often associated with physiological dysfunctions. Since biological soft tissues are hydrated, viscoelasticity is likely suitable to represent its solid-like behavior using elasticity and fluid-like behavior using viscosity. Shear wave elastography is a non-invasive imaging technology invented for clinical applications that has shown promise to characterize various tissue viscoelasticity. It is based on measuring and analyzing velocities and attenuations of propagated shear waves. In this review, principles and technical developments of shear wave elastography for viscoelasticity characterization from organ to cellular levels are presented, and different imaging modalities used to track shear wave propagation are described. At a macroscopic scale, techniques for inducing shear waves using an external mechanical vibration, an acoustic radiation pressure or a Lorentz force are reviewed along with imaging approaches proposed to track shear wave propagation, namely ultrasound, magnetic resonance, optical, and photoacoustic means. Then, approaches for theoretical modeling and tracking of shear waves are detailed. Following it, some examples of applications to characterize the viscoelasticity of various organs are given. At a microscopic scale, a novel cellular shear wave elastography method using an external vibration and optical microscopy is illustrated. Finally, current limitations and future directions in shear wave elastography are presented.

Influence of fibrosis progression on the viscous properties of in vivo liver tissue elucidated by shear wave dispersion in multifrequency MR elastography

PURPOSE

Many elastography studies have shown that liver stiffness increases with fibrosis and thus can be used as a reliable marker for noninvasively staging fibrosis. However, the sensitivity of viscosity-related mechanical parameters, such as shear wave dispersion, to liver fibrosis is less well understood.

METHODS

In this proof-of-concept study, 15 healthy volunteers and 37 patients with chronic liver disease and biopsy-proven fibrosis were prospectively investigated by MR elastography at six drive frequencies of 35-60 Hz. Maps of shear wave speed (SWS, in m/s) and loss angle (φ, in rad), as a marker of stiffness and viscous properties, respectively, were generated using tomoelastography data processing. The Child-Pugh score was used to assess cirrhosis severity.

RESULTS

While SWS increased with fibrosis (F0: 1.53 ± 0.11 m/s, F1-F3: 1.71 ± 0.17 m/s, F4: 2.50 ± 0.39 m/s; P < 0.001), φ remained unchanged during mild to severe fibrosis (F0: 0.63 ± 0.05 rad, F1-F3: 0.60 ± 0.05 rad, P = 0.21) but increased in cirrhosis (F4: 0.81 ± 0.16 rad; P < 0.001). Correspondingly, the slope of SWS-dispersion within the investigated range of vibration frequencies increased from insignificant (F0-F3: 0.010 ± 0.007 m/s/Hz) to significant (F4: 0.038 ± 0.025 m/s/Hz; P = 0.005). Significant correlation with the Child-Pugh score was found for φ (R = 0.60, P = 0.01) but not for SWS.

CONCLUSION

Although cirrhosis is associated with liver stiffening and, intuitively, transition towards more rigid material properties, the observed increases in φ and slope of SWS-dispersion indicate abnormally high mechanical friction in cirrhotic livers. This biophysical signature might provide a prognostic imaging marker for the detection of pathological processes associated with fibrosis independent of stiffness.

Obesity-induced pancreatopathy in rats is reversible after bariatric surgery

Obesity is a risk factor for pancreatic diseases. Bariatric surgery is one of the most efficient treatments of morbid obesity. The aims were to assess pancreatic endocrine and exocrine lesions in obese rats, to analyze effects of bariatric surgery. Sixty-three male Wistar rats were included in five groups: 2 fed with high fat diet (HFD) or normal diet for 3 months, 2 fed with HFD or normal diet for 6 months; 1 group fed with HFD and undergoing bariatric surgery (n = 30). Quantitative MR imaging was performed in HFD6, ND6 and HFD3-BS. Pancreas specimens were analyzed after sacrifice for adipocyte infiltration, fibrosis, acinar-ductal metaplasia, abnormality of Langerhans islets (HHF: hypertrophy, hypervascularisation, fibrosis), and hemosiderin deposits in acinar or endocrine locations. We found that HFD6 rats had more fibro-inflammatory islets (P = 0.0139) and acinar-ducal metaplasia (P = 0.0843) than HFD3 rats. Rats with HFD3+6 had more fibro-inflammatory islets (P < 0.0001), hemosiderin deposits (p < 0.0001), fat infiltration (P = 0.0008) and acinar-ductal metaplasia lesions (P = 0.0424). Weight increase was associated with glycoregulation abnormalities (r = 0.44, P = 0.08) and adipocyte infiltrations (P = 0.009). After surgery, less fibro-inflammatory islets (P = 0.0004), fat and iron infiltrates (P = 0.005 and P = 0.06), and acino-ductal metaplasia (P = 0.05) were observed compared to HFD6 rats. MR image quantifications revealed increased elasticity, fat fraction, and R2 and a decreased elasticity wave dispersion coefficient in the high fat groups that reversed after surgery. MRI parameters were in strong correlation with respective histological counterparts. In conclusion, obese rats develop pancreatic inflammatory lesions with acinar-ductal metaplasia in acinar location and the endocrine-exocrine interface. These changes can be prevented by bariatric surgery. Quantitative MR imaging is accurate in identifying early pancreatic lesions.

Magnetic resonance elastography of liver and spleen: Methods and applications

The viscoelastic properties of the liver and spleen can be assessed with magnetic resonance elastography (MRE). Several actuators, MRI acquisition sequences and reconstruction algorithms have been proposed for this purpose. Reproducible results are obtained, especially when the examination is performed in standard conditions with the patient fasting. Accurate staging of liver fibrosis can be obtained by measuring liver stiffness or elasticity with MRE. Moreover, emerging evidence shows that assessing the tissue viscous parameters with MRE is useful for characterizing liver inflammation, non-alcoholic steatohepatitis, hepatic congestion, portal hypertension, and hepatic tumors. Further advances such as multifrequency acquisitions and compression-sensitive MRE may provide novel quantitative markers of hepatic and splenic mechanical properties that may improve the diagnosis of hepatic and splenic diseases.

Measurement of shear wave speed dispersion in the placenta by transient elastography: A preliminary ex vivo study

Background

Placental elasticity may be modified in women with placental insufficiency. Shear wave elastography (SWE) can measure this, using acoustic radiation force, but the safety of its use in pregnant women has not yet been demonstrated. Transient elastography (TE) is a safer alternative, but has not yet been applied to the placenta. Moreover, the dispersion of shear wave speed (SWS) as a function of frequency has received relatively little study for placental tissue, although it might improve the accuracy of biomechanical assessment.

Objective

To explore the feasibility and reproducibility of TE for placental analysis, to compare the values of SWS and Young’s modulus (YM) from TE and SWE, and to analyze SWS dispersion as a function of frequency ex vivo in normal placentas.

Materials and methods

Ten normal placentas were analyzed ex vivo by an Aixplorer ultrasound system as shear waves were generated by a vibrating plate and by using an Aixplorer system.

The frequency analysis provided the value of the exponent n from a fractional rheological model applied to the TE method. We calculated intra- and interobserver agreement for SWS and YM with 95% prediction intervals, created Bland-Altman plots with 95% limits of agreement, and estimated the intraclass correlation coefficient (ICC).

Main results

The mean SWS was 1.80 m/s +/- 0.28 (standard deviation) with the TE method at 50 Hz and 1.82 m/s +/-0.13 with SWE (P = 0.912). No differences were observed between the central and peripheral regions of placentas with either TE or SWE. With TE, the intraobserver ICC for SWS was 0.68 (0.50–0.82), and the interobserver ICC for SWS 0.65 (0.37–0.85). The mean parameter n obtained from the fractional rheological model was 1.21 +/- 0.12, with variable values of n for any given SWS.

Conclusions

TE is feasible and reproducible on placentas ex vivo. The frequency analysis of SWS provides additional information about placental elasticity and appears to be able to distinguish differences between placental structures.

Assessing the Non-tumorous Liver: Implications for Patient Management and Surgical Therapy

INTRODUCTION

Hepatic resection is performed for various benign and malignant liver tumors. Over the last several decades, there have been improvements in the surgical technique and postoperative care of patients undergoing liver surgery. Despite this, liver failure following an extended hepatic resection remains a critical potential postoperative complication. Patients with underlying parenchymal liver diseases are at particular risk of liver failure due to impaired liver regeneration with an associated mortality risk as high as 60 to 90%. In addition, live donor liver transplantation requires a thorough presurgical assessment of the donor liver to minimize the risk of postoperative complications.

RESULTS AND CONCLUSION

Recently, cross-sectional imaging assessment of diffuse liver diseases has gained momentum due to its ability to provide both anatomical and functional assessments of normal and abnormal tissues. Various imaging techniques are being employed to assess diffuse liver diseases including magnetic resonance imaging (MRI), computed tomography (CT), and ultrasound (US). MRI has the ability to detect abnormal intracellular and molecular processes and tissue architecture. CT has a high spatial resolution, while US provides real-time imaging, is inexpensive, and readily available. We herein review current state-of-the-art techniques to assess the underlying non-tumorous liver. Specifically, we summarize current approaches to evaluating diffuse liver diseases including fatty liver alcoholic or non-alcoholic (NAFLD, AFLD), hepatic fibrosis (HF), and iron deposition (ID) with a focus on advanced imaging techniques for non-invasive assessment along with their implications for patient management. In addition, the role of and techniques to assess hepatic volume in hepatic surgery are discussed.

Distinguishing between Hepatic Inflammation and Fibrosis with MR Elastography

Purpose To investigate the utility of magnetic resonance (MR) elastography-derived mechanical properties in the discrimination of hepatic inflammation and fibrosis in the early stages of chronic liver diseases. Materials and Methods All studies were approved by the institutional animal care and use committee. A total of 187 animals were studied, including 182 mice and five pigs. These animals represented five different liver diseases with a varying combination and extent of hepatic inflammation, fibrosis, congestion, and portal hypertension. Multifrequency three-dimensional MR elastography was performed, and shear stiffness, storage modulus, shear loss modulus, and damping ratio were calculated for all animals. Necroinflammation, fibrosis, and portal pressure were either histologically scored or biochemically and physically quantified in all animals. Two-sided Welch t tests were used to evaluate mean differences between disease and control groups. Spearman correlation analyses were used to evaluate the relationships between mechanical parameters and quantitative fibrosis extent (hydroxyproline concentration) and portal pressure. Results Liver stiffness and storage modulus increased with progressively developed fibrosis and portal hypertension (mean stiffness at 80 Hz and 48-week feeding, 0.51 kPa ± 0.12 in the steatohepatitis group vs 0.29 kPa ± 0.01 in the control group; P = .02). Damping ratio and shear loss modulus can be used to distinguish inflammation from fibrosis at early stages of disease, even before the development of histologically detectable necroinflammation and fibrosis (mean damping ratio at 80 Hz and 20-week feeding, 0.044 ± 0.012 in the steatohepatitis group vs 0.014 ± 0.008 in the control group; P < .001). Damping ratio and liver stiffness vary differently with respect to cause of portal hypertension (ie, congestion- or cirrhosis-induced hypertension). These differentiation abilities have frequency-dependent variations. Conclusion Liver stiffness and damping ratio measurements can extend hepatic MR elastography to potentially enable assessment of necroinflammatory, congestive, and fibrotic processes of chronic liver diseases. © RSNA, 2017 Online supplemental material is available for this article.

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.

Microvasculature alters the dispersion properties of shear waves--a multi-frequency MR elastography study

Magnetic Resonance Elastography (MRE) uses macroscopic shear wave propagation to quantify mechanical properties of soft tissues. Micro-obstacles are capable of affecting the macroscopic dispersion properties of shear waves. Since disease or therapy can change the mechanical integrity and organization of vascular structures, MRE should be able to sense these changes if blood vessels represent a source for wave scattering. To verify this, MRE was performed to quantify alteration of the shear wave speed cs due to the presence of vascular outgrowths using an aortic ring model. Eighteen fragments of rat aorta included in a Matrigel matrix (n=6 without outgrowths, n=6 with a radial outgrowth extent of ~600 µm and n=6 with ~850 µm) were imaged using a 7 Tesla MR scanner (Bruker, PharmaScan). High resolution anatomical images were acquired in addition to multi-frequency MRE (ν = 100, 115, 125, 135 and 150 Hz). Average cs was measured within a ring of ~900 µm thickness encompassing the aorta and were normalized to cs0 of the corresponding Matrigel. The frequency dependence was fit to the power law model cs ~ν(y). After scanning, optical microscopy was performed to visualize outgrowths. Results demonstrated that in presence of vascular outgrowths (1) normalized cs significantly increased for the three highest frequencies (Kruskal-Wallis test, P = 0.0002 at 125 Hz and P = 0.002 at 135 Hz and P = 0.003 at 150 Hz) but not for the two lowest (Kruskal-Wallis test, P = 0.63 at 100 Hz and P = 0.87 at 115 Hz), and (2) normalized cs followed a power law behavior not seen in absence of vascular outgrowths (ANOVA test, P < 0.0001). These results showed that vascular outgrowths acted as micro-obstacles altering the dispersion relationships of propagating shear waves and that MRE could provide valuable information about microvascular changes.

New imaging techniques for liver diseases

Newly developed or advanced methods of ultrasonography and MR imaging provide combined anatomical and quantitative functional information about diffuse and focal liver diseases. Ultrasound elastography has a central role for staging liver fibrosis and an increasing role in grading portal hypertension; dynamic contrast-enhanced ultrasonography may improve tumor characterization. In clinical practice, MR imaging examinations currently include diffusion-weighted and dynamic MR imaging, enhanced with extracellular or hepatobiliary contrast agents. Moreover, quantitative parameters obtained with diffusion-weighted MR imaging, dynamic contrast-enhanced MR imaging and MR elastography have the potential to characterize further diffuse and focal liver diseases, by adding information about tissue cellularity, perfusion, hepatocyte transport function and visco-elasticity. The multiparametric capability of ultrasonography and more markedly of MR imaging gives the opportunity for high diagnostic performance by combining imaging biomarkers. However, image acquisition and post-processing methods should be further standardized and validated in multicenter trials.

MR elastography and diffusion-weighted imaging of ex vivo prostate cancer: quantitative comparison to histopathology

The purpose of this work was (1) to develop a magnetic resonance elastography (MRE) system for imaging of the ex vivo human prostate and (2) to assess the diagnostic power of mono-frequency and multi-frequency MRE and diffusion weighted imaging (DWI) alone and combined as correlated with histopathology in a patient study. An electromagnetic driver was designed specifically for MRE studies in small-bore MR scanners. Ex vivo prostate specimens (post-fixation) of 14 patients who underwent radical prostatectomy were imaged with MRE at 7 T (nine cases had DWI). In six patients, the MRE examination was performed at three frequencies (600, 800, 1000 Hz) to extract the power-law exponent Gamma. The images were registered to wholemount pathology slides marked with the Gleason score. The areas under the receiver-operator-characteristic curves (AUC) were calculated. The methods were validated in a phantom study and it was demonstrated that (i) the driver does not interfere with the acquisition process and (ii) the driver can generate amplitudes greater than 100 µm for frequencies less than 1 kHz. In the quantitative study, cancerous tissue with Gleason score at least 3 + 3 was distinguished from normal tissue in the peripheral zone (PZ) with an average AUC of 0.75 (Gd ), 0.75 (Gl ), 0.70 (Gamma-Gd ), 0.68 (apparent diffusion coefficient, ADC), and 0.82 (Gd  + Gl  + ADC). The differentiation between PZ and central gland was modest for Gd (p < 0.07), Gl (p < 0.06) but not significant for Gamma (p < 0.2). A correlation of 0.4 kPa/h was found between the fixation time of the prostate specimen and the stiffness of the tissue, which could affect the diagnostic power results. DWI and MRE may provide complementary information; in fact MRE performed better than ADC in distinguishing normal from cancerous tissue in some cases. Multi-frequency (Gamma) analysis did not appear to improve the results. However, in light of the effect of tissue fixation, the clinical implication of our results may be inconclusive and more experiments are needed.