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.

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.

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.

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.

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.

Biomarkers of liver fibrosis: prospective comparison of multimodal magnetic resonance, serum algorithms and transient elastography.

BACKGROUND AND AIMS: Accurate biomarkers for quantifying liver fibrosis are important for clinical practice and trial end-points. We compared the diagnostic performance of magnetic resonance imaging (MRI), including gadoxetate-enhanced MRI and 31P-MR spectroscopy, with fibrosis stage and serum fibrosis algorithms in a clinical setting. Also, in a subset of patients, MR- and transient elastography (MRE and TE) was evaluated when available. METHODS: Patients were recruited prospectively if they were scheduled to undergo liver biopsy on a clinical indication due to elevated liver enzyme levels without decompensated cirrhosis. Within a month of the clinical work-up, an MR-examination and liver needle biopsy were performed on the same day. Based on late-phase gadoxetate-enhanced MRI, a mathematical model calculated hepatobiliary function (relating to OATP1 and MRP2). The hepatocyte gadoxetate uptake rate (KHep) and the normalised liver-to-spleen contrast ratio (LSC_N10) were also calculated. Nine serum fibrosis algorithms were investigated (GUCI, King's Score, APRI, FIB-4, Lok-Index, NIKEI, NASH-CRN regression score, Forns' score, and NAFLD-fibrosis score). RESULTS: The diagnostic performance (AUROC) for identification of significant fibrosis (F2-4) was 0.78, 0.80, 0.69, and 0.78 for MRE, TE, LSC_N10, and GUCI, respectively. For the identification of advanced fibrosis (F3-4), the AUROCs were 0.93, 0.84, 0.81, and 0.82 respectively. CONCLUSION: MRE and TE were superior for non-invasive identification of significant fibrosis. Serum fibrosis algorithms developed for specific liver diseases are applicable in this cohort of diverse liver diseases aetiologies. Gadoxetate-MRI was sufficiently sensitive to detect the low function losses associated with fibrosis. None was able to efficiently distinguish between stages within the low fibrosis stages.Lay summaryExcessive accumulation of scar tissue, fibrosis, in the liver is an important aspect in chronic liver disease. To replace the invasive needle biopsy, we have explored non-invasive methods to assess liver fibrosis. In our study we found that elastographic methods, which assess the mechanical properties of the liver, are superior in assessing fibrosis in a clinical setting. Of interest from a clinical trial point-of-view, none of the tested methods was sufficiently accurate to distinguish between adjacent moderate fibrosis stages.

Magnetic Resonance Elastography reveals effects of anti-angiogenic glioblastoma treatment on tumor stiffness and captures progression in an orthotopic mouse model.

BACKGROUND: Anti-angiogenic treatment of glioblastoma (GBM) complicates radiologic monitoring. We evaluated magnetic resonance elastography (MRE) as an imaging tool for monitoring the efficacy of anti-VEGF treatment of GBM. METHODS: Longitudinal studies were performed in an orthotopic GBM xenograft mouse model. Animals treated with B20 anti-VEGF antibody were compared to untreated controls regarding survival (n = 13), classical MRI-contrasts and biomechanics as quantified via MRE (n = 15). Imaging was performed on a 7 T small animal horizontal bore MRI scanner. MRI and MRE parameters were compared to histopathology. RESULTS: Anti-VEGF-treated animals survived longer than untreated controls (p = 0.0011) with progressively increased tumor volume in controls (p = 0.0001). MRE parameters viscoelasticity |G*| and phase angle Y significantly decreased in controls (p = 0.02 for |G*| and p = 0.0071 for Y). This indicates that untreated tumors became softer and more elastic than viscous with progression. Tumor volume in treated animals increased more slowly than in controls, indicating efficacy of the therapy, reaching significance only at the last time point (p = 0.02). Viscoelasticity and phase angle Y tended to decrease throughout therapy, similar as for control animals. However, in treated animals, the decrease in phase angle Y was significantly attenuated and reached statistical significance at the last time point (p = 0.04). Histopathologically, control tumors were larger and more heterogeneous than treated tumors. Vasculature was normalized in treated tumors compared with controls, which showed abnormal vasculature and necrosis. In treated tumors, a higher amount of myelin was observed within the tumor area (p = 0.03), likely due to increased tumor invasion. Stiffness of the contralateral hemisphere was influenced by tumor mass effect and edema. CONCLUSIONS: Anti-angiogenic GBM treatment prolonged animal survival, slowed tumor growth and softening, but did not prevent progression. MRE detected treatment effects on tumor stiffness; the decrease of viscoelasticity and phase angle in GBM was attenuated in treated animals, which might be explained by normalized vasculature and greater myelin preservation within treated tumors. Thus, further investigation of MRE is warranted to understand the potential for MRE in monitoring treatment in GBM patients by complementing existing MRI techniques.

Towards noninvasive estimation of tumour pressure by utilising MR elastography and nonlinear biomechanical models: a simulation and phantom study.

The solid and fluid pressures of tumours are often elevated relative to surrounding tissue. This increased pressure is known to correlate with decreased treatment efficacy and potentially with tumour aggressiveness and therefore, accurate noninvasive estimates of tumour pressure would be of great value. We present a proof-of-concept method to infer the total tumour pressure, that is the sum of the fluid and solid parts, by examining stiffness in the peritumoural tissue with MR elastography and utilising nonlinear biomechanical models. The pressure from the tumour deforms the surrounding tissue leading to changes in stiffness. Understanding and accounting for these biases in stiffness has the potential to enable estimation of total tumour pressure. Simulations are used to validate the method with varying pressure levels, tumour shape, tumour size, and noise levels. Results show excellent matching in low noise cases and still correlate well with higher noise. Percent error remains near or below 10% for higher pressures in all noise level cases. Reconstructed pressures were also calculated from experiments with a catheter balloon embedded in a plastisol phantom at multiple inflation levels. Here the reconstructed pressures generally match the increases in pressure measured during the experiments. Percent errors between average reconstructed and measured pressures at four inflation states are 17.9%, 52%, 23.2%, and 0.9%. Future work will apply this method to in vivo data, potentially providing an important biomarker for cancer diagnosis and treatment.

Nonlinear viscoelastic constitutive model for bovine liver tissue.

Soft tissue mechanical characterisation is important in many areas of medical research. Examples span from surgery training, device design and testing, sudden injury and disease diagnosis. The liver is of particular interest, as it is the most commonly injured organ in frontal and side motor vehicle crashes, and also assessed for inflammation and fibrosis in chronic liver diseases. Hence, an extensive rheological characterisation of liver tissue would contribute to advancements in these areas, which are dependent upon underlying biomechanical models. The aim of this paper is to define a liver constitutive equation that is able to characterise the nonlinear viscoelastic behaviour of liver tissue under a range of deformations and frequencies. The tissue response to large amplitude oscillatory shear (1-50%) under varying preloads (1-20%) and frequencies (0.5-2 Hz) is modelled using viscoelastic-adapted forms of the Mooney-Rivlin, Ogden and exponential models. These models are fit to the data using classical or modified objective norms. The results show that all three models are suitable for capturing the initial nonlinear regime, with the latter two being capable of capturing, simultaneously, the whole deformation range tested. The work presented here provides a comprehensive analysis across several material models and norms, leading to an identifiable constitutive equation that describes the nonlinear viscoelastic behaviour of the liver.

REVIEW: MR elastography of brain tumors.

MR elastography allows non-invasive quantification of the shear modulus of tissue, i.e. tissue stiffness and viscosity, information that offers the potential to guide presurgical planning for brain tumor resection. Here, we review brain tumor MRE studies with particular attention to clinical applications. Studies that investigated MRE in patients with intracranial tumors, both malignant and benign as well as primary and metastatic, were queried from the Pubmed/Medline database in August 2018. Reported tumor and normal appearing white matter stiffness values were extracted and compared as a function of tumor histopathological diagnosis and MRE vibration frequencies. Because different studies used different elastography hardware, pulse sequences, reconstruction inversion algorithms, and different symmetry assumptions about the mechanical properties of tissue, effort was directed to ensure that similar quantities were used when making inter-study comparisons. In addition, because different methodologies and processing pipelines will necessarily bias the results, when pooling data from different studies, whenever possible, tumor values were compared with the same subject's contralateral normal appearing white matter to minimize any study-dependent bias. The literature search yielded 10 studies with a total of 184 primary and metastatic brain tumor patients. The group mean tumor stiffness, as measured with MRE, correlated with intra-operatively assessed stiffness of meningiomas and pituitary adenomas. Pooled data analysis showed significant overlap between shear modulus values across brain tumor types. When adjusting for the same patient normal appearing white matter shear modulus values, meningiomas were the stiffest tumor-type. MRE is increasingly being examined for potential in brain tumor imaging and might have value for surgical planning. However, significant overlap of shear modulus values between a number of different tumor types limits applicability of MRE for diagnostic purposes. Thus, further rigorous studies are needed to determine specific clinical applications of MRE for surgical planning, disease monitoring and molecular stratification of brain tumors.

Investigating the Contribution of Collagen to the Tumor Biomechanical Phenotype with Noninvasive Magnetic Resonance Elastography.

Increased stiffness in the extracellular matrix (ECM) contributes to tumor progression and metastasis. Therefore, stromal modulating therapies and accompanying biomarkers are being developed to target ECM stiffness. Magnetic resonance (MR) elastography can noninvasively and quantitatively map the viscoelastic properties of tumors in vivo and thus has clear clinical applications. Herein, we used MR elastography, coupled with computational histopathology, to interrogate the contribution of collagen to the tumor biomechanical phenotype and to evaluate its sensitivity to collagenase-induced stromal modulation. Elasticity (G d) and viscosity (G l) were significantly greater for orthotopic BT-474 (G d = 5.9 ± 0.2 kPa, G l = 4.7 ± 0.2 kPa, n = 7) and luc-MDA-MB-231-LM2-4 (G d = 7.9 ± 0.4 kPa, G l = 6.0 ± 0.2 kPa, n = 6) breast cancer xenografts, and luc-PANC1 (G d = 6.9 ± 0.3 kPa, G l = 6.2 ± 0.2 kPa, n = 7) pancreatic cancer xenografts, compared with tumors associated with the nervous system, including GTML/Trp53KI/KI medulloblastoma (G d = 3.5 ± 0.2 kPa, G l = 2.3 ± 0.2 kPa, n = 7), orthotopic luc-D-212-MG (G d = 3.5 ± 0.2 kPa, G l = 2.3 ± 0.2 kPa, n = 7), luc-RG2 (G d = 3.5 ± 0.2 kPa, G l = 2.3 ± 0.2 kPa, n = 5), and luc-U-87-MG (G d = 3.5 ± 0.2 kPa, G l = 2.3 ± 0.2 kPa, n = 8) glioblastoma xenografts, intracranially propagated luc-MDA-MB-231-LM2-4 (G d = 3.7 ± 0.2 kPa, G l = 2.2 ± 0.1 kPa, n = 7) breast cancer xenografts, and Th-MYCN neuroblastomas (G d = 3.5 ± 0.2 kPa, G l = 2.3 ± 0.2 kPa, n = 5). Positive correlations between both elasticity (r = 0.72, P < 0.0001) and viscosity (r = 0.78, P < 0.0001) were determined with collagen fraction, but not with cellular or vascular density. Treatment with collagenase significantly reduced G d (P = 0.002) and G l (P = 0.0006) in orthotopic breast tumors. Texture analysis of extracted images of picrosirius red staining revealed significant negative correlations of entropy with G d (r = -0.69, P < 0.0001) and G l (r = -0.76, P < 0.0001), and positive correlations of fractal dimension with G d (r = 0.75, P < 0.0001) and G l (r = 0.78, P < 0.0001). MR elastography can thus provide sensitive imaging biomarkers of tumor collagen deposition and its therapeutic modulation. SIGNIFICANCE: MR elastography enables noninvasive detection of tumor stiffness and will aid in the development of ECM-targeting therapies.

Assessing Tumor Mechanics by MR Elastography at Different Strain Levels.

BACKGROUND: Malignant tumors are associated with increased tissue rigidity, which can be an indicator of tumor progression. MR elastography (MRE) has the potential to study the variations of tumor mechanical properties. ex vivo studies have shown the ability of MRE to assess increase of mechanical properties; nevertheless, it has not yet been observed in vivo. PURPOSE: To propose a method to assess the increase in mechanical properties of tumors in vivo under static external compression using MRE. STUDY TYPE: Prospective, experimental study. ANIMAL MODEL: Forty-six SCID mice with subcutaneous tumor implantation (patient-derived hepatocellular carcinoma xenografts, Model 1, n = 13, and Model 2, n = 33). FIELD STRENGTH/SEQUENCE: 7.0T; a spin echo sequence was used for anatomical images and a modified spin echo sequence for elastography acquisitions with a vibration frequency of 600 Hz. ASSESSMENT: An inflatable balloon was placed on the abdomen to apply a load to the tumor. MRE acquisitions were performed at the basal state and at increasing compression levels. Anatomical images were used to calculate the octahedral shear strain between the tumor at the basal strain state and each strain level. For six mice (Model 2), each static preloading scan was acquired twice consecutively without moving the mouse to evaluate repeatability. Statistical Tests: The Bland-Altman method was used to assess repeatability. Correlations between tumor stiffness and deformation were evaluated with Pearson correlation coefficients. RESULTS: For stiffness (G*), a good repeatability was obtained between the acquisitions; the limits of agreement of the Bland-Altman test were [-10.17%; 11.49%] with an absolute bias of 0.66%. A significant correlation between tumor stiffness and deformation was observed for both models (Model 1: r = 0.57, P < 0.0001 and Model 2: r = 0.31, P < 0.0001). DATA CONCLUSION: We establish that tumor mechanical properties can increase under mechanical compression. This increase can effectively be monitored using a proposed MRE setup. LEVEL OF EVIDENCE: 2 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2019;50:1982-1989.

Rheological determinants for simultaneous staging of hepatic fibrosis and inflammation in patients with chronic liver disease.

The purpose of this study is to investigate the use of fundamental rheological parameters as quantified by MR elastography (MRE) to measure liver fibrosis and inflammation simultaneously in humans. MRE was performed on 45 patients at 3 T using a vibration frequency of 56 Hz. Fibrosis and inflammation scores were obtained from liver biopsies. Biomechanical properties were quantified in terms of complex shear modulus G* as well as shear wave phase velocity c and shear wave attenuation α. A rheological fractional derivative order model was used to investigate the linear dependence of the free model parameters (dispersion slope y, intrinsic speed c0 , and intrinsic relaxation time τ) on histopathology. Leave-one-out cross-validation was then utilized to demonstrate the effectiveness of the model. The intrinsic speed c0 increases with hepatic fibrosis, while an increased relaxation time τ is reflective of more inflammation of the liver parenchyma. The dispersion slope y does not depend either on fibrosis or on inflammation. The proposed rheological model, given this specific parameterization, establishes the functional dependences of biomechanical parameters on histological fibrosis and inflammation. The leave-one-out cross-validation demonstrates that the model allows identification, from the MRE measurements, of the histology scores when grouped into low-/high-grade fibrosis and low-/high-grade inflammation with significance levels of P = 0.0004 (fibrosis) and P = 0.035 (inflammation). The functional dependences of intrinsic speed and relaxation time on fibrosis and inflammation, respectively, shed new light onto the impact hepatic pathological changes on liver tissue biomechanics in humans. The dispersion slope y appears to represent a structural parameter of liver parenchyma not impacted by the severity of fibrosis/inflammation present in this patient cohort. This specific parametrization of the well-established rheological fractional order model is valuable for the clinical assessment of both fibrosis and inflammation scores, going beyond the capability of the plain shear modulus measurement commonly used for MRE.

Stiffness reconstruction methods for MR elastography.

Assessment of tissue stiffness is desirable for clinicians and researchers, as it is well established that pathophysiological mechanisms often alter the structural properties of tissue. Magnetic resonance elastography (MRE) provides an avenue for measuring tissue stiffness and has a long history of clinical application, including staging liver fibrosis and stratifying breast cancer malignancy. A vital component of MRE consists of the reconstruction algorithms used to derive stiffness from wave-motion images by solving inverse problems. A large range of reconstruction methods have been presented in the literature, with differing computational expense, required user input, underlying physical assumptions, and techniques for numerical evaluation. These differences, in turn, have led to varying accuracy, robustness, and ease of use. While most reconstruction techniques have been validated against in silico or in vitro phantoms, performance with real data is often more challenging, stressing the robustness and assumptions of these algorithms. This article reviews many current MRE reconstruction methods and discusses the aforementioned differences. The material assumptions underlying the methods are developed and various approaches for noise reduction, regularization, and numerical discretization are discussed. Reconstruction methods are categorized by inversion type, underlying assumptions, and their use in human and animal studies. Future directions, such as alternative material assumptions, are also discussed.

Characterization of Small Renal Tumors With Magnetic Resonance Elastography: A Feasibility Study.

OBJECTIVES: The aim of this study was to explore the feasibility of magnetic resonance elastography (MRE) for characterizing indeterminate small renal tumors (SRTs) as part of a multiparametric magnetic resonance (MR) imaging protocol. MATERIALS AND METHODS: After institutional review board approval and informed consent were obtained, 21 prospective adults (15 men; median age, 55 years; age range, 25-72 years) with SRT were enrolled. Tumors (2-5 cm Ø) were imaged using 3-directional, gradient echo MRE. Viscoelastic parametric maps (shear wave velocity [c] and attenuation [α]) were analyzed by 2 independent radiologists. Interobserver agreement (Bland-Altman statistics and intraclass correlation coefficients) was assessed. Anatomical T2-weighted, dynamic contrast-enhanced (DCE) and diffusion sequences completed the acquisition protocol. Imaging parameters were compared between groups (Mann-Whitney U test). RESULTS: Quality of MRE was good in 18 cases (mean nonlinearity <50%), including 1 papillary renal cell carcinoma and 1 metanephric adenoma. A cohort of 5 oncocytomas and 11 clear-cell renal cell carcinomas (ccRCCs) was analyzed for statistical differences. The MRE viscoelastic parameters were the strongest imaging discriminators: oncocytomas displayed significantly lower shear velocity c (median, 0.77 m/s; interquartile range [IQR], 0.76-0.79) (P = 0.007) and higher shear attenuation α (median, 0.087 mm; IQR, 0.082-0.087) (P = 0.008) than ccRCC (medians, 0.92 m/s and 0.066 mm; IQR, 0.84-0.97 and 0.054-0.074, respectively). T2 signal intensity ratio (tumor/renal cortex) was lower in oncocytomas (P = 0.02). The DCE and diffusion MR parameters overlapped substantially (P ≥ 0.1). Oncocytomas displayed a consistent MRE viscoelastic profile, corresponding to data point clustering in a bidimensional scatter plot. Values for MRE intraclass correlation coefficient were 0.982 for c and 0.984 for α, indicating excellent interobserver agreement. CONCLUSIONS: Magnetic resonance elastography is feasible for SRT characterization; MRE viscoelastic parameters were stronger discriminators between oncocytoma and ccRCC than anatomical, DCE and diffusion MR imaging parameters.

Thymidine Metabolism as a Confounding Factor for 3'-Deoxy-3'-18F-Fluorothymidine Uptake After Therapy in a Colorectal Cancer Model.

Noninvasive monitoring of tumor therapy response helps in developing personalized treatment strategies. Here, we performed sequential PET and diffusion-weighted MRI to evaluate changes induced by a FOLFOX-like combination chemotherapy in colorectal cancer xenografts, to identify the cellular and molecular determinants of these imaging biomarkers. Methods: Tumor-bearing CD1 nude mice, engrafted with FOLFOX-sensitive Colo205 colorectal cancer xenografts, were treated with FOLFOX (5-fluorouracil, leucovorin, and oxaliplatin) weekly. On days 1, 2, 6, 9, and 13 of therapy, tumors were assessed by in vivo imaging and ex vivo analyses. In addition, HCT116 xenografts, which did not respond to the FOLFOX treatment, were imaged on day 1 of therapy. Results: In Colo205 xenografts, FOLFOX induced a profound increase in uptake of the proliferation PET tracer 3'-deoxy-3'-18F-fluorothymidine (18F-FLT) accompanied by increases in markers for proliferation (Ki-67, thymidine kinase 1) and for activated DNA damage response (γH2AX), whereas the effect on cell death was minimal. Because tracer uptake was unaltered in the HCT116 model, these changes appear to be specific for tumor response. Conclusion: We demonstrated that 18F-FLT PET can noninvasively monitor cancer treatment-induced molecular alterations, including thymidine metabolism and DNA damage response. The cellular or imaging changes may not, however, be directly related to therapy response as assessed by volumetric measurements.