Biomechanical Assessment of Liver Integrity: Prospective Evaluation of Mechanical Versus Acoustic MR Elastography.

BACKGROUND: Magnetic resonance elastography (MRE) can quantify tissue biomechanics noninvasively, including pathological hepatic states like metabolic dysfunction-associated steatohepatitis. PURPOSE: To compare the performance of 2D/3D-MRE using the gravitational (GT) transducer concept with the current commercial acoustic (AC) solution utilizing a 2D-MRE approach. Additionally, quality index markers (QIs) were proposed to identify image pixels with sufficient quality for reliably estimating tissue biomechanics. STUDY TYPE: Prospective. POPULATION: One hundred seventy participants with suspected or confirmed liver disease (median age, 57 years [interquartile range (IQR), 46-65]; 66 females), and 11 healthy volunteers (median age, 31 years [IQR, 27-34]; 5 females). FIELD STRENGTH/SEQUENCE: Participants were scanned twice at 1.5 T and 60 Hz vibration frequency: first, using AC-MRE (2D-MRE, spin-echo EPI sequence, 11 seconds breath-hold), and second, using GT-MRE (2D- and 3D-MRE, gradient-echo sequence, 14 seconds breath-hold). ASSESSMENT: Image analysis was performed by four independent radiologists and one biomedical engineer. Additionally, superimposed analytic plane shear waves of known wavelength and attenuation at fixed shear modulus were used to propose pertinent QIs. STATISTICAL TESTS: Spearman's correlation coefficient (r) was applied to assess the correlation between modalities. Interreader reproducibility was evaluated using Bland-Altman bias and reproducibility coefficients. P-values <0.05 were considered statistically significant. RESULTS: Liver stiffness quantified via GT-2D/3D correlated well with AC-2D (r ≥ 0.89 [95% CI: 0.85-0.92]) and histopathological grading (r ≥ 0.84 [95% CI: 0.72-0.91]), demonstrating excellent agreement in Bland-Altman plots and between readers (κ ≥ 0.86 [95% CI: 0.81-0.91]). However, GT-2D showed a bias in overestimating stiffness compared to GT-3D. Proposed QIs enabled the identification of pixels deviating beyond 10% from true stiffness based on a combination of total wave amplitude, temporal sinusoidal nonlinearity, and wave signal-to-noise ratio for GT-3D. CONCLUSION: GT-MRE represents an alternative to AC-MRE for noninvasive liver tissue characterization. Both GT-2D and 3D approaches correlated strongly with the established commercial approach, offering advanced capabilities in abdominal imaging compared to AC-MRE. EVIDENCE LEVEL: 1 TECHNICAL EFFICACY: Stage 2.

Assessment of Tumor Cell Invasion and Radiotherapy Response in Experimental Glioma by Magnetic Resonance Elastography.

BACKGROUND: Gliomas are highly invasive brain neoplasms. MRI is the most important tool to diagnose and monitor glioma but has shortcomings. In particular, the assessment of tumor cell invasion is insufficient. This is a clinical dilemma, as recurrence can arise from MRI-occult glioma cell invasion. HYPOTHESIS: Tumor cell invasion, tumor growth and radiotherapy alter the brain parenchymal microstructure and thus are assessable by diffusion tensor imaging (DTI) and MR elastography (MRE). STUDY TYPE: Experimental, animal model. ANIMAL MODEL: Twenty-three male NMRI nude mice orthotopically implanted with S24 patient-derived glioma cells (experimental mice) and 9 NMRI nude mice stereotactically injected with 1 µL PBS (sham-injected mice). FIELD STRENGTH/SEQUENCE: 2D and 3D T2-weighted rapid acquisition with refocused echoes (RARE), 2D echo planar imaging (EPI) DTI, 2D multi-slice multi-echo (MSME) T2 relaxometry, 3D MSME MRE at 900 Hz acquired at 9.4 T (675 mT/m gradient strength). ASSESSMENT: Longitudinal 4-weekly imaging was performed for up to 4 months. Tumor volume was assessed in experimental mice (n = 10 treatment-control, n = 13 radiotherapy). The radiotherapy subgroup and 5 sham-injected mice underwent irradiation (3 × 6 Gy) 9 weeks post-implantation/sham injection. MRI-/MRE-parameters were assessed in the corpus callosum and tumor core/injection tract. Imaging data were correlated to light sheet microscopy (LSM) and histology. STATISTICAL TESTS: Paired and unpaired t-tests, a P-value ≤0.05 was considered significant. RESULTS: From week 4 to 8, a significant callosal stiffening (4.44 ± 0.22 vs. 5.31 ± 0.29 kPa) was detected correlating with LSM-proven tumor cell invasion. This was occult to all other imaging metrics. Histologically proven tissue destruction in the tumor core led to an increased T2 relaxation time (41.65 ± 0.34 vs. 44.83 ± 0.66 msec) and ADC (610.2 ± 12.27 vs. 711.2 ± 13.42 × 10-6 mm2/s) and a softening (5.51 ± 0.30 vs. 4.24 ± 0.29 kPa) from week 8 to 12. Radiotherapy slowed tumor progression. DATA CONCLUSION: MRE is promising for the assessment of key glioma characteristics. EVIDENCE LEVEL: NA TECHNICAL EFFICACY: Stage 2.

Tumor biomechanics as a novel imaging biomarker to assess response to immunotherapy in a murine glioma model.

Glioblastoma is the most common and aggressive primary malignant brain tumor with poor prognosis. Novel immunotherapeutic approaches are currently under investigation. Even though magnetic resonance imaging (MRI) is the most important imaging tool for treatment monitoring, response assessment is often hampered by therapy-related tissue changes. As tumor and therapy-associated tissue reactions differ structurally, we hypothesize that biomechanics could be a pertinent imaging proxy for differentiation. Longitudinal MRI and magnetic resonance elastography (MRE) were performed to monitor response to immunotherapy with a toll-like receptor 7/8 agonist in orthotopic syngeneic experimental glioma. Imaging results were correlated to histology and light sheet microscopy data. Here, we identify MRE as a promising non-invasive imaging method for immunotherapy-monitoring by quantifying changes in response-related tumor mechanics. Specifically, we show that a relative softening of treated compared to untreated tumors is linked to the inflammatory processes following therapy-induced re-education of tumor-associated myeloid cells. Mechanistically, combined effects of myeloid influx and inflammation including extracellular matrix degradation following immunotherapy form the basis of treated tumors being softer than untreated glioma. This is a very early indicator of therapy response outperforming established imaging metrics such as tumor volume. The overall anti-tumor inflammatory processes likely have similar effects on human brain tissue biomechanics, making MRE a promising tool for gauging response to immunotherapy in glioma patients early, thereby strongly impacting patient pathway.

Making sense of scattering: Seeing microstructure through shear waves.

The physics of shear waves traveling through matter carries fundamental insights into its structure, for instance, quantifying stiffness for disease characterization. However, the origin of shear wave attenuation in tissue is currently not properly understood. Attenuation is caused by two phenomena: absorption due to energy dissipation and scattering on structures such as vessels fundamentally tied to the material's microstructure. Here, we present a scattering theory in conjunction with magnetic resonance imaging, which enables the unraveling of a material's innate constitutive and scattering characteristics. By overcoming a three-order-of-magnitude scale difference between wavelength and average intervessel distance, we provide noninvasively a macroscopic measure of vascular architecture. The validity of the theory is demonstrated through simulations, phantoms, in vivo mice, and human experiments and compared against histology as gold standard. Our approach expands the field of imaging by using the dispersion properties of shear waves as macroscopic observable proxies for deciphering the underlying ultrastructures.

Magnetic resonance elastography resolving all gross anatomical segments of the kidney during controlled hydration.

Introduction: Magnetic resonance elastography (MRE) is a non-invasive method to quantify biomechanical properties of human tissues. It has potential in diagnosis and monitoring of kidney disease, if established in clinical practice. The interplay of flow and volume changes in renal vessels, tubule, urinary collection system and interstitium is complex, but physiological ranges of in vivo viscoelastic properties during fasting and hydration have never been investigated in all gross anatomical segments simultaneously. Method: Ten healthy volunteers underwent two imaging sessions, one following a 12-hour fasting period and the second after a drinking challenge of >10 mL per kg body weight (60-75 min before the second examination). High-resolution renal MRE was performed using a novel driver with rotating eccentric mass placed at the posterior-lateral wall to couple waves (50 Hz) to the kidney. The biomechanical parameters, shear wave speed (cs in m/s), storage modulus (Gd in kPa), loss modulus (Gl in kPa), phase angle (Υ=2πatanGlGd) and attenuation (α in 1/mm) were derived. Accurate separation of gross anatomical segments was applied in post-processing (whole kidney, cortex, medulla, sinus, vessel). Results: High-quality shear waves coupled into all gross anatomical segments of the kidney (mean shear wave displacement: 163 ± 47 µm, mean contamination of second upper harmonics <23%, curl/divergence: 4.3 ± 0.8). Regardless of the hydration state, median Gd of the cortex and medulla (0.68 ± 0.11 kPa) was significantly higher than that of the sinus and vessels (0.48 ± 0.06 kPa), and consistently, significant differences were found in cs, Υ, and Gl (all p < 0.001). The viscoelastic parameters of cortex and medulla were not significantly different. After hydration sinus exhibited a small but significant reduction in median Gd by -0.02 ± 0.04 kPa (p = 0.01), and, consequently, the cortico-sinusoidal-difference in Gd increased by 0.04 ± 0.07 kPa (p = 0.05). Only upon hydration, the attenuation in vessels became lower (0.084 ± 0.013 1/mm) and differed significantly from the whole kidney (0.095 ± 0.007 1/mm, p = 0.01). Conclusion: High-resolution renal MRE with an innovative driver and well-defined 3D segmentation can resolve all renal segments, especially when including the sinus in the analysis. Even after a prolonged hydration period the approach is sensitive to small hydration-related changes in the sinus and in the cortico-sinusoidal-difference.

Single Breath-Hold MR Elastography for Fast Biomechanical Probing of Pancreatic Stiffness.

BACKGROUND: Pancreatic ductal adenocarcinoma (PDAC) stromal disposition is thought to influence chemotherapy efficacy and increase tissue stiffness, which could be quantified noninvasively via MR elastography (MRE). Current methods cause position-based errors in pancreas location over time, hampering accuracy. It would be beneficial to have a single breath-hold acquisition. PURPOSE: To develop and test a single breath-hold three-dimensional MRE technique utilizing prospective undersampling and a compressed sensing reconstruction (CS-MRE). STUDY TYPE: Prospective. POPULATION: A total of 30 healthy volunteers (HV) (31 ± 9 years; 33% male) and five patients with PDAC (69 ± 5 years; 80% male). FIELD STRENGTH/SEQUENCE: 3-T, GRE Ristretto MRE. ASSESSMENT: First, optimization of multi breath-hold MRE was done in 10 HV using four combinations of vibration frequency, number of measured wave-phase offsets, and TE and looking at MRE quality measures in the pancreas head. Second, viscoelastic parameters delineated in the pancreas head or tumor of CS-MRE were compared against (I) 2D and (II) 3D four breath-hold acquisitions in HV (N = 20) and PDAC patients. Intrasession repeatability was assessed for CS-MRE in a subgroup of healthy volunteers (N = 15). STATISTICAL TESTS: Tests include repeated measures analysis of variance (ANOVA), Bland-Altman analysis, and coefficients of variation (CoVs). A P-value <.05 was considered statistically significant. RESULTS: Optimization of the four breath-hold acquisitions resulted in 40 Hz vibration frequency, five wave-phases, and echo time (TE) = 6.9 msec as the preferred method (4BH-MRE). CS-MRE quantitative results did not differ from 4BH-MRE. Shear wave speed (SWS) and phase angle differed significantly between HV and PDAC patients using 4BH-MRE or CS-MRE. The limits of agreement for SWS were [-0.09, 0.10] m/second and the within-subject CoV was 4.8% for CS-MRE. DATA CONCLUSION: CS-MRE might allow a single breath-hold MRE acquisition with comparable SWS and phase angle as 4BH-MRE, and it may still enable to differentiate between HV and PDAC. LEVEL OF EVIDENCE: 2 Technical Efficacy Stage: 2.

Viscoelastic Property of the Brain Assessed With Magnetic Resonance Elastography and Its Association With Glymphatic System in Neurologically Normal Individuals.

OBJECTIVE: To investigate the feasibility of assessing the viscoelastic properties of the brain using magnetic resonance elastography (MRE) and a novel MRE transducer to determine the relationship between the viscoelastic properties and glymphatic function in neurologically normal individuals. MATERIALS AND METHODS: This prospective study included 47 neurologically normal individuals aged 23-74 years (male-to-female ratio, 21:26). The MRE was acquired using a gravitational transducer based on a rotational eccentric mass as the driving system. The magnitude of the complex shear modulus |G*| and the phase angle ϕ were measured in the centrum semiovale area. To evaluate glymphatic function, the Diffusion Tensor Image Analysis Along the Perivascular Space (DTI-ALPS) method was utilized and the ALPS index was calculated. Univariable and multivariable (variables with P < 0.2 from the univariable analysis) linear regression analyses were performed for |G*| and ϕ and included sex, age, normalized white matter hyperintensity (WMH) volume, brain parenchymal volume, and ALPS index as covariates. RESULTS: In the univariable analysis for |G*|, age (P = 0.005), brain parenchymal volume (P = 0.152), normalized WMH volume (P = 0.011), and ALPS index (P = 0.005) were identified as candidates with P < 0.2. In the multivariable analysis, only the ALPS index was independently associated with |G*|, showing a positive relationship (ß = 0.300, P = 0.029). For ϕ, normalized WMH volume (P = 0.128) and ALPS index (P = 0.015) were identified as candidates for multivariable analysis, and only the ALPS index was independently associated with ϕ (ß = 0.057, P = 0.039). CONCLUSION: Brain MRE using a gravitational transducer is feasible in neurologically normal individuals over a wide age range. The significant correlation between the viscoelastic properties of the brain and glymphatic function suggests that a more organized or preserved microenvironment of the brain parenchyma is associated with a more unimpeded glymphatic fluid flow.

MR elastography identifies regions of extracellular matrix reorganization associated with shorter survival in glioblastoma patients.

Background: Biomechanical tissue properties of glioblastoma tumors are heterogeneous, but the molecular mechanisms involved and the biological implications are poorly understood. Here, we combine magnetic resonance elastography (MRE) measurement of tissue stiffness with RNA sequencing of tissue biopsies to explore the molecular characteristics of the stiffness signal. Methods: MRE was performed preoperatively in 13 patients with glioblastoma. Navigated biopsies were harvested during surgery and classified as "stiff" or "soft" according to MRE stiffness measurements (|G*|norm). Twenty-two biopsies from eight patients were analyzed by RNA sequencing. Results: The mean whole-tumor stiffness was lower than normal-appearing white matter. The surgeon's stiffness evaluation did not correlate with the MRE measurements, which suggests that these measures assess different physiological properties. Pathway analysis of the differentially expressed genes between "stiff" and "soft" biopsies showed that genes involved in extracellular matrix reorganization and cellular adhesion were overexpressed in "stiff" biopsies. Supervised dimensionality reduction identified a gene expression signal separating "stiff" and "soft" biopsies. Using the NIH Genomic Data Portal, 265 glioblastoma patients were divided into those with (n = 63) and without (n = 202) this gene expression signal. The median survival time of patients with tumors expressing the gene signal associated with "stiff" biopsies was 100 days shorter than that of patients not expressing it (360 versus 460 days, hazard ratio: 1.45, P < .05). Conclusion: MRE imaging of glioblastoma can provide noninvasive information on intratumoral heterogeneity. Regions of increased stiffness were associated with extracellular matrix reorganization. An expression signal associated with "stiff" biopsies correlated with shorter survival of glioblastoma patients.

Investigating the contribution of hyaluronan to the breast tumour microenvironment using multiparametric MRI and MR elastography.

Hyaluronan (HA) is a key component of the dense extracellular matrix in breast cancer, and its accumulation is associated with poor prognosis and metastasis. Pegvorhyaluronidase alfa (PEGPH20) enzymatically degrades HA and can enhance drug delivery and treatment response in preclinical tumour models. Clinical development of stromal-targeted therapies would be accelerated by imaging biomarkers that inform on therapeutic efficacy in vivo. Here, PEGPH20 response was assessed by multiparametric magnetic resonance imaging (MRI) in three orthotopic breast tumour models. Treatment of 4T1/HAS3 tumours, the model with the highest HA accumulation, reduced T1 and T2 relaxation times and the apparent diffusion coefficient (ADC), and increased the magnetisation transfer ratio, consistent with lower tissue water content and collapse of the extracellular space. The transverse relaxation rate R2 * increased, consistent with greater erythrocyte accessibility following vascular decompression. Treatment of MDA-MB-231 LM2-4 tumours reduced ADC and dramatically increased tumour viscoelasticity measured by MR elastography. Correlation matrix analyses of data from all models identified ADC as having the strongest correlation with HA accumulation, suggesting that ADC is the most sensitive imaging biomarker of tumour response to PEGPH20.

Diffusion phase-imaging in anisotropic media using non-linear gradients for diffusion encoding.

Diffusion MRI classically uses gradient fields that vary linearly in space to encode the diffusion of water molecules in the signal magnitude by tempering its intensity. In spin ensembles, a presumably equal number of particles move in positive and negative direction, resulting in approximately zero change in net phase. Hence, in classical diffusion weighted MRI with a linear gradient field, the phase does not carry any information as the incoherent motion of the spins only impacts the magnitude of the signal. Conversely, when the linear gradient field is replaced with one that varies quadratically over space, the diffusion of water molecules in anisotropic media does give rise to a change in net phase and preserves large portion of the signal around the saddle point of the gradient field. In this work, the phase evolution of anisotropic fibre phantoms in the presence of quadratic gradient fields was studied in Monte Carlo simulations and diffusion MRI experiments. The simulations confirm the dependence of the phase change on the degree of anisotropy of the media and the diffusion weighting, as predicted by the derived analytic model. First MR experiments show a phase change depending on the diffusion time in an anisotropic synthetic fibre phantom, and approximately zero phase change for the experiment repeated in an isotropic agar phantom. As predicted by the analytic model, an increase of the diffusion time by approximately a factor of two leads to an increase of approximately a factor of two in the signal phase.

Single Breath-Hold 3-Dimensional Magnetic Resonance Elastography Depicts Liver Fibrosis and Inflammation in Obese Patients.

OBJECTIVES: Three-dimensional (3D) magnetic resonance elastography (MRE) measures liver fibrosis and inflammation but requires several breath-holds that hamper clinical acceptance. The aim of this study was to evaluate the technical and clinical feasibility of a single breath-hold 3D MRE sequence as a means of measuring liver fibrosis and inflammation in obese patients. METHODS: From November 2020 to December 2021, subjects were prospectively enrolled and divided into 2 groups. Group 1 included healthy volunteers (n = 10) who served as controls to compare the single breath-hold 3D MRE sequence with a multiple-breath-hold 3D MRE sequence. Group 2 included liver patients (n = 10) who served as participants to evaluate the clinical feasibility of the single breath-hold 3D MRE sequence in measuring liver fibrosis and inflammation. Controls and participants were scanned at 60 Hz mechanical excitation with the single breath-hold 3D MRE sequence to retrieve the magnitude of the complex-valued shear modulus (|G*| [kPa]), the shear wave speed (Cs [m/s]), and the loss modulus (G" [kPa]). The controls were also scanned with a multiple-breath-hold 3D MRE sequence for comparison, and the participants had histopathology (Ishak scores) for correlation with Cs and G". RESULTS: For the 10 controls, 5 were female, and the mean age and body mass index were 33.1 ± 9.5 years and 23.0 ± 2.1 kg/m 2 , respectively. For the 10 participants, 8 were female, and the mean age and body mass index were 45.1 ± 16.5 years and 33.1 ± 4.0 kg/m 2 (obese range), respectively. All participants were suspected of having nonalcoholic fatty liver disease. Bland-Altman analysis of the comparison in controls shows there are nonsignificant differences in |G*|, Cs, and G" below 6.5%, suggesting good consensus between the 2 sequences. For the participants, Cs and G" correlated significantly with Ishak fibrosis and inflammation grades, respectively ( ρ = 0.95, P < 0.001, and ρ = 0.84, P = 0.002). CONCLUSION: The single breath-hold 3D MRE sequence may be effective in measuring liver fibrosis and inflammation in obese patients.

Magnetic resonance elastography to study the effect of amyloid plaque accumulation in a mouse model.

BACKGROUND AND PURPOSE: Biomechanical changes in the brain have not been fully elucidated in Alzheimer's disease (AD). We aimed to investigate the effect of ß-amyloid accumulation on mouse brain viscoelasticity. METHODS: Magnetic resonance elastography was used to calculate magnitude of the viscoelastic modulus (|G*|), elasticity (Gd ), and viscosity (Gl ) in the whole brain parenchyma (WB) and bilateral hippocampi of 9 transgenic J20 (AD) mice (5 males/4 females) and 10 wild-type (WT) C57BL/6 mice (5 males/5 females) at 11 and 14 months of age. RESULTS: Cross-sectional analyses showed no significant difference between AD and WT mice at either timepoints. No sex-specific differences were observed at 11 months of age, but AD females showed significantly higher hippocampal |G*| and Gl and WB |G*|, Gd , and Gl compared to both AD and WT males at 14 months of age. Similar trending differences were found between female AD and female WT animals but did not reach significance. Longitudinal analyses showed significant increases in hippocampal |G*|, Gd , and Gl , and significant decreases in WB |G*|, Gd , and Gl between 11 and 14 months in both AD and WT mice. Each subgroup showed significant increases in all hippocampal and significant decreases in all WB measures, with the exception of AD females, which showed no significant changes in WB |G*|, Gd , or Gl . CONCLUSION: Aging had region-specific effects on cerebral viscoelasticity, namely, WB softening and hippocampal stiffening. Amyloid plaque deposition may have sex-specific effects, which require further scrutiny.

Decreased tissue stiffness in glioblastoma by MR elastography is associated with increased cerebral blood flow.

PURPOSE: Understanding how mechanical properties relate to functional changes in glioblastomas may help explain different treatment response between patients. The aim of this study was to map differences in biomechanical and functional properties between tumor and healthy tissue, to assess any relationship between them and to study their spatial distribution. METHODS: Ten patients with glioblastoma and 17 healthy subjects were scanned using MR Elastography, perfusion and diffusion MRI. Stiffness and viscosity measurements G' and G'', cerebral blood flow (CBF), apparent diffusion coefficient (ADC) and fractional anisotropy (FA) were measured in patients' contrast-enhancing tumor, necrosis, edema, and gray and white matter, and in gray and white matter for healthy subjects. A regression analysis was used to predict CBF as a function of ADC, FA, G' and G''. RESULTS: Median G' and G'' in contrast-enhancing tumor were 13% and 37% lower than in normal-appearing white matter (P < 0.01), and 8% and 6% lower in necrosis than in contrast-enhancing tumor, respectively (P < 0.05). Tumors showed both inter-patient and intra-patient heterogeneity. Measurements approached values in normal-appearing tissue when moving outward from the tumor core, but abnormal tissue properties were still present in regions of normal-appearing tissue. Using both a linear and a random-forest model, prediction of CBF was improved by adding MRE measurements to the model (P < 0.01). CONCLUSIONS: The inclusion of MRE measurements in statistical models helped predict perfusion, with stiffer tissue associated with lower perfusion values.

Magnetic resonance elastography with guided pressure waves.

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

Magnetic Resonance Elastography Reconstruction for Anisotropic Tissues.

Elastography has become widely used clinically for characterising changes in soft tissue mechanics that are associated with altered tissue structure and composition. However, some soft tissues, such as muscle, are not isotropic as is assumed in clinical elastography implementations. This limits the ability of these methods to capture changes in anisotropic tissues associated with disease. The objective of this study was to develop and validate a novel elastography reconstruction technique suitable for estimating the linear viscoelastic mechanical properties of transversely isotropic soft tissues. We derived a divergence-free formulation of the governing equations for acoustic wave propagation through a linearly transversely isotropic viscoelastic material, and transformed this into a weak form. This was then implemented into a finite element framework, enabling the analysis of wave input data and tissue structural fibre orientations, in this case based on diffusion tensor imaging. To validate the material constants obtained with this method, numerous in silico phantom experiments were run which encompassed a range of variations in wave input directions, material properties, fibre structure and noise. The method was also tested on ex vivo muscle and in vivo human volunteer calf muscles, and compared with a previous curl-based inversion method. The new method robustly extracted the transversely isotropic shear moduli (G⊥', G∥', G″) from the in silico phantom tests with minimal bias, including in the presence of experimentally realistic levels of noise in either fibre orientation or wave data. This new method performed better than the previous method in the presence of noise. Anisotropy estimates from the ex vivo muscle phantom agreed well with rheological tests. In vivo experiments on human calf muscles were able to detect increases in muscle shear moduli with passive muscle stretch. This new reconstruction method can be applied to quantify tissue mechanical properties of anisotropic soft tissues, such as muscle, in health and disease.

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.

A viscoelastic model for human myocardium.

Understanding the biomechanics of the heart in health and disease plays an important role in the diagnosis and treatment of heart failure. The use of computational biomechanical models for therapy assessment is paving the way for personalized treatment, and relies on accurate constitutive equations mapping strain to stress. Current state-of-the art constitutive equations account for the nonlinear anisotropic stress-strain response of cardiac muscle using hyperelasticity theory. While providing a solid foundation for understanding the biomechanics of heart tissue, most current laws neglect viscoelastic phenomena observed experimentally. Utilizing experimental data from human myocardium and knowledge of the hierarchical structure of heart muscle, we present a fractional nonlinear anisotropic viscoelastic constitutive model. The model is shown to replicate biaxial stretch, triaxial cyclic shear and triaxial stress relaxation experiments (mean error ∼7.68%), showing improvements compared to its hyperelastic (mean error ∼24%) counterparts. Model sensitivity, fidelity and parameter uniqueness are demonstrated. The model is also compared to rate-dependent biaxial stretch as well as different modes of biaxial stretch, illustrating extensibility of the model to a range of loading phenomena. STATEMENT OF SIGNIFICANCE: The viscoelastic response of human heart tissues has yet to be integrated into common constitutive models describing cardiac mechanics. In this work, a fractional viscoelastic modeling approach is introduced based on the hierarchical structure of heart tissue. From these foundations, the current state-of-the-art biomechanical models of the heart muscle are transformed using fractional viscoelasticity, replicating passive muscle function across multiple experimental tests. Comparisons are drawn with current models to highlight the improvements of this approach and predictive responses show strong qualitative agreement with experimental data. The proposed model presents the first constitutive model aimed at capturing viscoelastic nonlinear response across multiple testing regimes, providing a platform for better understanding the biomechanics of myocardial tissue in health and disease.

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.

Assessment of Imaging Modalities Against Liver Biopsy in Nonalcoholic Fatty Liver Disease: The Amsterdam NAFLD-NASH Cohort.

BACKGROUND: Noninvasive diagnostic methods are urgently required in disease stratification and monitoring in nonalcoholic fatty liver disease (NAFLD). Multiparametric magnetic resonance imaging (MRI) is a promising technique to assess hepatic steatosis, inflammation, and fibrosis, potentially enabling noninvasive identification of individuals with active and advanced stages of NAFLD. PURPOSE: To examine the diagnostic performance of multiparametric MRI for the assessment of disease severity along the NAFLD disease spectrum with comparison to histological scores. STUDY TYPE: Prospective, cohort. POPULATION: Thirty-seven patients with NAFLD. FIELD STRENGTH/SEQUENCE: Multiparametric MRI at 3.0 T consisted of magnetic resonance (MR) spectroscopy (MRS) with multi-echo stimulated-echo acquisition mode, magnitude-based and three-point Dixon using a two-dimensional multi-echo gradient echo, MR elastography (MRE) using a generalized multishot gradient-recalled echo sequence and intravoxel incoherent motion (IVIM) using a multislice diffusion weighted single-shot echo-planar sequence. ASSESSMENT: Histological steatosis grades were compared to proton density fat fraction measured by MRS (PDFFMRS ), magnitude-based MRI (PDFFMRI-M ), and three-point Dixon (PDFFDixon ), as well as FibroScan® controlled attenuation parameter (CAP). Fibrosis and disease activity were compared to IVIM and MRE. FibroScan® liver stiffness measurements were compared to fibrosis levels. Diagnostic performance of all imaging parameters was determined for distinction between simple steatosis and nonalcoholic steatohepatitis (NASH). STATISTICAL TESTS: Spearman's rank test, Kruskal-Wallis test, Dunn's post-hoc test with Holm-Bonferroni P-value adjustment, receiver operating characteristic curve analysis. A P-value <0.05 was considered statistically significant. RESULTS: Histological steatosis grade correlated significantly with PDFFMRS (rs  = 0.66, P < 0.001), PDFFMRI-M (rs  = 0.68, P < 0.001), and PDFFDixon (rs  = 0.67, P < 0.001), whereas no correlation was found with CAP. MRE and IVIM diffusion and perfusion significantly correlated with disease activity (rs  = 0.55, P < 0.001, rs  = -0.40, P = 0.016, rs  = -0.37, P = 0.027, respectively) and fibrosis (rs  = 0.55, P < 0.001, rs  = -0.46, P = 0.0051; rs  = -0.53, P < 0.001, respectively). MRE and IVIM diffusion had the highest area-under-the-curve for distinction between simple steatosis and NASH (0.79 and 0.73, respectively). DATA CONCLUSION: Multiparametric MRI is a promising method for noninvasive, accurate, and sensitive distinction between simple hepatic steatosis and NASH, as well as for the assessment of steatosis and fibrosis severity. LEVEL OF EVIDENCE: 2 TECHNICAL EFFICACY: 2.