MR Elastography of the Liver and the Spleen Using a Piezoelectric Driver, Single-Shot Wave-Field Acquisition, and Multifrequency Dual Parameter Reconstruction

Perfusion alters stiffness of deep gray matter

Viscoelastic properties of the brain reflect tissue architecture at multiple length scales. However, little is known about the relation between vital tissue functions, such as perfusion, and the macroscopic mechanical properties of cerebral tissue. In this study, arterial spin labelling is paired with magnetic resonance elastography to investigate the relationship between tissue stiffness and cerebral blood flow (CBF) in the in vivo human brain. The viscoelastic modulus, | G*|, and CBF were studied in deep gray matter (DGM) of 14 healthy male volunteers in the following sub-regions: putamen, nucleus accumbens, hippocampus, thalamus, globus pallidus, and amygdala. CBF was further normalized by vessel area data to obtain the flux rate q which is proportional to the perfusion pressure gradient. The striatum (represented by putamen and nucleus accumbens) was distinct from the other DGM regions by displaying markedly higher stiffness and perfusion values. q was a predictive marker for DGM stiffness as analyzed by linear regression | G*| =  q·(4.2 ± 0.6)kPa·s + (0.80 ± 0.06)kPa ( R²= 0.92, P = 0.006). These results suggest a high sensitivity of MRE in DGM to perfusion pressure. The distinct mechano-vascular properties of striatum tissue, as compared to the rest of DGM, may reflect elevated perfusion pressure, which could explain the well-known susceptibility of the putamen to hemorrhages.

Combining viscoelasticity, diffusivity and volume of the hippocampus for the diagnosis of Alzheimer's disease based on magnetic resonance imaging

Dementia due to Alzheimer's Disease (AD) is a neurodegenerative disease for which treatment strategies at an early stage are of great clinical importance. So far, there is still a lack of non-invasive diagnostic tools to sensitively detect AD in early stages and to predict individual disease progression. Magnetic resonance elastography (MRE) of the brain may be a promising novel tool. In this proof-of-concept study, we investigated whether multifrequency-MRE (MMRE) can detect differences in hippocampal stiffness between patients with clinical diagnosis of dementia due to AD and healthy controls (HC). Further, we analyzed if the combination of three MRI-derived parameters, i.e., hippocampal stiffness, hippocampal volume and mean diffusivity (MD), improves diagnostic accuracy.

Diagnostic criteria for probable dementia due to AD were in line with the NINCDS-ADRDA criteria and were verified through history-taking (patient and informant), neuropsychological testing, routine blood results and routine MRI to exclude other medical causes of a cognitive decline.

21 AD patients and 21 HC (median age 75 years) underwent MMRE and structural MRI, from which hippocampal volume and MD were calculated. From the MMRE-images maps of the magnitude |G*| and phase angle φ of the complex shear modulus were reconstructed using multifrequency inversion. Median values of |G*| and φ were extracted within three regions of interest (hippocampus, thalamus and whole brain white matter). To test the predictive value of the main outcome parameters, we performed receiver operating characteristic (ROC) curve analyses.

Hippocampal stiffness (|G*|) and viscosity (φ) were significantly lower in the patient group (both p < 0.001). ROC curve analyses showed an area under the curve (AUC) for | G*| of 0.81 [95%CI 0.68–0.94]; with sensitivity 86%, specificity 67% for cutoff at |G*| = 980 Pa) and for φ an AUC of 0.79 [95%CI 0.66–0.93]. In comparison, the AUC of MD and hippocampal volume were 0.83 [95%CI 0.71–0.95] and 0.86 [95%CI 0.74–0.97], respectively. A combined ROC curve of |G*|, MD and hippocampal volume yielded a significantly improved AUC of 0.90 [95%CI 0.81–0.99].

In conclusion, we demonstrated reduced hippocampal stiffness and reduced hippocampal viscosity, as determined by MMRE, in patients with clinical diagnosis of dementia of the AD type. Diagnostic sensitivity was further improved by the combination with two other MRI-based hippocampal parameters. These findings motivate further investigation whether MMRE can detect decreased brain stiffness already in pre-dementia stages, and whether these changes predict cognitive decline.

•

Non-invasive methods for early detection of AD are lacking.

•

MRE of the brain is a promising new non-invasive diagnostic tool.

•

We demonstrate reduced hippocampal stiffness (|G*|) in AD patients.

•

|G*| distinguishes healthy and demented with 86% sensitivity and 67% specificity.

•

Combining hippocampal stiffness, MD and volume improved diagnostic accuracy.

Artificial neural networks for stiffness estimation in magnetic resonance elastography

PURPOSE

To investigate the feasibility of using artificial neural networks to estimate stiffness from MR elastography (MRE) data.

METHODS

Artificial neural networks were fit using model-based training patterns to estimate stiffness from images of displacement using a patch size of ∼1 cm in each dimension. These neural network inversions (NNIs) were then evaluated in a set of simulation experiments designed to investigate the effects of wave interference and noise on NNI accuracy. NNI was also tested in vivo, comparing NNI results against currently used methods.

RESULTS

In 4 simulation experiments, NNI performed as well or better than direct inversion (DI) for predicting the known stiffness of the data. Summary NNI results were also shown to be significantly correlated with DI results in the liver (R²  = 0.974) and in the brain (R²  = 0.915), and also correlated with established biological effects including fibrosis stage in the liver and age in the brain. Finally, repeatability error was lower in the brain using NNI compared to DI, and voxel-wise modeling using NNI stiffness maps detected larger effects than using DI maps with similar levels of smoothing.

CONCLUSION

Artificial neural networks represent a new approach to inversion of MRE data. Summary results from NNI and DI are highly correlated and both are capable of detecting biologically relevant signals. Preliminary evidence suggests that NNI stiffness estimates may be more resistant to noise than an algebraic DI approach. Taken together, these results merit future investigation into NNIs to improve the estimation of stiffness in small regions. Magn Reson Med 80:351-360, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

MR elastography of the brain and its application in neurological diseases

Magnetic resonance elastography (MRE) is an imaging technique for noninvasively and quantitatively assessing tissue stiffness, akin to palpation. MRE is further able assess the mechanical properties of tissues that cannot be reached by hand including the brain. The technique is a three-step process beginning with the introduction of shear waves into the tissue of interest by applying an external vibration. Next, the resulting motion is imaged using a phase-contrast MR pulse sequence with motion encoding gradients that are synchronized to the vibration. Finally, the measured displacement images are mathematically inverted to compute a map of the estimated stiffness. In the brain, the technique has demonstrated strong test-retest repeatability with typical errors of 1% for measuring global stiffness, 2% for measuring stiffness in the lobes of the brain, and 3-7% for measuring stiffness in subcortical gray matter. In healthy volunteers, multiple studies have confirmed that stiffness decreases with age, while more recent studies have demonstrated a strong relationship between viscoelasticity and behavioral performance. Furthermore, several studies have demonstrated the sensitivity of brain stiffness to neurodegeneration, as stiffness has been shown to decrease in multiple sclerosis and in several forms of dementia. Moreover, the spatial pattern of stiffness changes varies among these different classes of dementia. Finally, MRE is a promising tool for the preoperative assessment of intracranial tumors, as it can measure both tumor consistency and adherence to surrounding tissues. These factors are important predictors of surgical difficulty. In brief, MRE demonstrates potential value in a number of neurological diseases. However, significant opportunity remains to further refine the technique and better understand the underlying physiology.

Magnetic resonance elastography for examining developmental changes in the mechanical properties of the brain

Magnetic resonance elastography (MRE) is a quantitative imaging technique for noninvasively characterizing tissue mechanical properties, and has recently emerged as a valuable tool for neuroimaging. The measured mechanical properties reflect the microstructural composition and organization of neural tissue, and have shown significant effects in many neurological conditions and normal, healthy aging, and evidence has emerged supporting novel relationships between mechanical structure and cognitive function. The sensitivity of MRE to brain structure, function, and health make it an ideal technique for studying the developing brain; however, brain MRE studies on children and adolescents have only just begun. In this article, we review brain MRE and its findings, discuss its potential role in developmental neuroimaging, and provide suggestions for researchers interested in adopting this technique.

Finite element analysis to investigate variability of MR elastography in the human thigh

PURPOSE

To develop finite element analysis (FEA) of magnetic resonance elastography (MRE) in the human thigh and investigate inter-individual variability of measurement of muscle mechanical properties.

METHODS

Segmentation was performed on MRI datasets of the human thigh from 5 individuals and FEA models consisting of 12 muscles and surrounding tissue created. The same material properties were applied to each tissue type and a previously developed transient FEA method of simulating MRE using Abaqus was performed at 4 frequencies. Synthetic noise was applied to the simulated data at various levels before inversion was performed using the Elastography Software Pipeline. Maps of material properties were created and visually assessed to determine key features. The coefficient of variation (CoV) was used to assess the variability of measurements in each individual muscle and in the groups of muscles across the subjects. Mean measurements for the set of muscles were ranked in size order and compared with the expected ranking.

RESULTS

At noise levels of 2% the CoV in measurements of |G^*| ranged from 5.3 to 21.9% and from 7.1 to 36.1% for measurements of ϕ in the individual muscles. A positive correlation (R² value 0.80) was attained when the expected and measured |G^*| ranking were compared, whilst a negative correlation (R² value 0.43) was found for ϕ.

CONCLUSIONS

Created elastograms demonstrated good definition of muscle structure and were robust to noise. Variability of measurements across the 5 subjects was dramatically lower for |G^*| than it was for ϕ. This large variability in ϕ measurements was attributed to artefacts.

MR elastography measurement of the effect of passive warmup prior to eccentric exercise on thigh muscle mechanical properties

Purpose

To investigate the effect of warmup by application of the thermal agent Deep Heat (DH) on muscle mechanical properties using magnetic resonance elastography (MRE) at 3T before and after exercise‐induced muscle damage (EIMD).

Materials and Methods

Twenty male participants performed an individualized protocol designed to induce EIMD in the quadriceps. DH was applied to the thigh in 50% of the participants before exercise. MRE, T₂‐weighted MRI, maximal voluntary contraction (MVC), creatine kinase (CK) concentration, and muscle soreness were measured before and after the protocol to assess EIMD effects. Five participants were excluded: four having not experienced EIMD and one due to incidental findings.

Results

Total workload performed during the EIMD protocol was greater in the DH group than the control group (P < 0.03), despite no significant differences in baseline MVC (P = 0.23). Shear stiffness |G*| increased in the rectus femoris (RF) muscle in both groups (P < 0.03); however, DH was not a significant between‐group factor (P =  0.15). MVC values returned to baseline faster in the DH group (5 days) than the control group (7 days). Participants who displayed hyperintensity on T₂‐weighted images had a greater stiffness increase following damage than those without: RF; 0.61 kPa vs. 0.15 kPa, P < 0.006, vastus intermedius; 0.34 kPa vs. 0.03 kPa, P = 0.06.

Conclusion

EIMD produces increased muscle stiffness as measured by MRE, with the change in |G*| significantly increased when T₂ hyperintensity was present. DH did not affect CK concentration or soreness; however, DH participants produced greater workload during the EIMD protocol and exhibited accelerated MVC recovery.

Level of Evidence: 1

Technical Efficacy: Stage 2

J. Magn. Reson. Imaging 2017;46:1115–1127.

Viscoelasticity of subcortical gray matter structures

Viscoelastic mechanical properties of the brain assessed with magnetic resonance elastography (MRE) are sensitive measures of microstructural tissue health in neurodegenerative conditions. Recent efforts have targeted measurements localized to specific neuroanatomical regions differentially affected in disease. In this work, we present a method for measuring the viscoelasticity in subcortical gray matter (SGM) structures, including the amygdala, hippocampus, caudate, putamen, pallidum, and thalamus. The method is based on incorporating high spatial resolution MRE imaging (1.6 mm isotropic voxels) with a mechanical inversion scheme designed to improve local measures in pre-defined regions (soft prior regularization [SPR]). We find that in 21 healthy, young volunteers SGM structures differ from each other in viscoelasticity, quantified as the shear stiffness and damping ratio, but also differ from the global viscoelasticity of the cerebrum. Through repeated examinations on a single volunteer, we estimate the uncertainty to be between 3 and 7% for each SGM measure. Furthermore, we demonstrate that the use of specific methodological considerations-higher spatial resolution and SPR-both decrease uncertainty and increase sensitivity of the SGM measures. The proposed method allows for reliable MRE measures of SGM viscoelasticity for future studies of neurodegenerative conditions. Hum Brain Mapp 37:4221-4233, 2016. © 2016 Wiley Periodicals, Inc.

Nonlinear multiscale regularisation in MR elastography: Towards fine feature mapping

Fine-featured elastograms may provide additional information of radiological interest in the context of in vivo elastography. Here a new image processing pipeline called ESP (Elastography Software Pipeline) is developed to create Magnetic Resonance Elastography (MRE) maps of viscoelastic parameters (complex modulus magnitude |G^*| and loss angle ϕ) that preserve fine-scale information through nonlinear, multi-scale extensions of typical MRE post-processing techniques.

METHODS

A new MRE image processing pipeline was developed that incorporates wavelet-domain denoising, image-driven noise estimation, and feature detection. ESP was first validated using simulated data, including viscoelastic Finite Element Method (FEM) simulations, at multiple noise levels. ESP images were compared with MDEV pipeline images, both in the FEM models and in three ten-subject cohorts of brain, thigh, and liver acquisitions. ESP and MDEV mean values were compared to 2D local frequency estimation (LFE) mean values for the same cohorts as a benchmark. Finally, the proportion of spectral energy at fine frequencies was quantified using the Reduced Energy Ratio (RER) for both ESP and MDEV.

RESULTS

Blind estimates of added noise (σ) were within 5.3% ± 2.6% of prescribed, and the same technique estimated σ in the in vivo cohorts at 1.7 ± 0.8%. A 5 × 5 × 5 truncated Gabor filter bank effectively detects local spatial frequencies at wavelengths λ ≤ 10px. For FEM inversions, mean |G^*| of hard target, soft target, and background remained within 8% of prescribed up to σ=20%, and mean ϕ results were within 10%, excepting hard target ϕ, which required redrawing around a ring artefact to achieve similar accuracy. Inspection of FEM |G^*| images showed some spatial distortion around hard target boundaries and inspection of ϕ images showed ring artefacts around the same target. For the in vivo cohorts, ESP results showed mean correlation of R=0.83 with MDEV and liver stiffness estimates within 7% of 2D-LFE results. Finally, ESP showed statistically significant increase in fine feature spectral energy as measured with RER for both |G^*| (p<1×10^-9) and ϕ (p<1×10^-3).

CONCLUSION

Information at finer frequencies can be recovered in ESP elastograms in typical experimental conditions, however scatter- and boundary-related artefacts may cause the fine features to have inaccurate values. In in vivo cohorts, ESP delivers an increase in fine feature spectral energy, and better performance with longer wavelengths, than MDEV while showing similar stability and robustness.

Tomoelastography by multifrequency wave number recovery from time-harmonic propagating shear waves

Palpation is one of the most sensitive, effective diagnostic practices, motivating the quantitative and spatially resolved determination of soft tissue elasticity parameters by medical ultrasound or MRI. However, this so-called elastography often suffers from limited anatomical resolution due to noise and insufficient elastic deformation, currently precluding its use as a tomographic modality on its own. We here introduce an efficient way of processing wave images acquired by multifrequency magnetic resonance elastography (MMRE), which relies on wave number reconstruction at different harmonic frequencies followed by their amplitude-weighted averaging prior to inversion. This results in compound maps of wave speed, which reveal variations in tissue elasticity in a tomographic fashion, i.e. an unmasked, slice-wise display of anatomical details at pixel-wise resolution. The method is demonstrated using MMRE data from the literature including abdominal and pelvic organs such as the liver, spleen, uterus body and uterus cervix. Even in small regions with low wave amplitudes, such as nucleus pulposus and spinal cord, elastic parameters consistent with literature values were obtained. Overall, the proposed method provides a simple and noise-robust strategy of in-plane wave analysis of MMRE data, with a pixel-wise resolution producing superior detail to MRE direct inversion methods.

Simultaneous measurement of hepatic and splenic stiffness using MR elastography: preliminary experience

PURPOSE

To compare MR elastography (MRE) using a single and a dual driver excitation for the quantification of hepatic and splenic stiffness (HS and SS), and to investigate the performance of HS and SS measured with single or dual driver excitation for the detection of liver cirrhosis in subjects with liver disease.

PATIENTS AND METHODS

This prospective HIPAA compliant and IRB approved study involved 49 subjects who underwent MRE at 3.0T, comparing three different acquisition methods (single driver on the liver, single driver on the spleen and dual driver acoustic excitation). A Mann-Whitney test was used to assess changes in stiffness values. Bland-Altman analysis was used to compare single and dual driver configurations for each organ. Performance for detection of liver cirrhosis was assessed using ROC analysis. Pearson correlation was used to estimate the dependence of HS and SS on spleen size.

RESULTS

There were 40 noncirrhotic and 9 cirrhotic patients. There was good agreement between stiffness values measured with a single or a dual driver (Bland-Altman limits of agreement -14.3 % to 18.9 % and -18.1 % to 29.7 %, CV 6.4 % and 9.4 %, for HS and SS. respectively). HS and SS were higher in subjects with liver cirrhosis (p < 0.001), with excellent detection performance (AUROC range 0.87-0.93). SS correlated strongly with spleen size (r = 0.69, p < 0.001), while HS showed weak correlation (r = 0.38, p = 0.006).

CONCLUSION

Using a dual acoustic driver configuration, hepatic and splenic stiffness can be simultaneously estimated with good concordance with single driver measurement.

In vivo abdominal magnetic resonance elastography for the assessment of portal hypertension before and after transjugular intrahepatic portosystemic shunt implantation

OBJECTIVE

The objective of this study was to investigate the correlation between hepatic venous pressure gradient (HVPG) and in vivo viscoelasticity of the liver and spleen before and after transjugular intrahepatic portosystemic shunt (TIPS) implantation.

MATERIALS AND METHODS

Ten patients with portal hypertension were examined twice by 3-dimensional multifrequency magnetic resonance elastography as well as prior and subsequent TIPS intervention; HVPG was also measured during TIPS placement. Five harmonic vibrations (25-60 Hz) were transferred to the abdominal region and recorded for the reconstruction of 2 viscoelastic constants, |G*| and φ, corresponding to the magnitude and the phase angle of the complex shear modulus G* of the liver and spleen.

RESULTS

All patients had cirrhosis, yielding high |G*| values in the liver (8.34 ± 2.18 kPa) and spleen (8.44 ± 1.36 kPa). In both organs, a decrease of |G*| after TIPS placement was observed (liver: 8.34 ± 2.18 kPa vs 7.02 ± 1.46 kPa, P = 0.01; spleen: 8.44 ± 1.36 kPa vs 7.06 ± 1.32 kPa, P = 0.01), whereas φ was insensitive to TIPS. Relative changes in |G*| of the spleen were correlated with the relative change of HVPG (R² = 0.659, P = 0.013).

CONCLUSIONS

The observed linear correlation between spleen viscoelasticity HVPG raises the prospect of an image-based noninvasive assessment of portal pressure by magnetic resonance elastography in the follow-up of TIPS placements.

Magnetic resonance elastography (MRE) in cancer: Technique, analysis, and applications

Tissue mechanical properties are significantly altered with the development of cancer. Magnetic resonance elastography (MRE) is a noninvasive technique capable of quantifying tissue mechanical properties in vivo. This review describes the basic principles of MRE and introduces some of the many promising MRE methods that have been developed for the detection and characterization of cancer, evaluation of response to therapy, and investigation of the underlying mechanical mechanisms associated with malignancy.

Quantifying the Elastic Property of Nine Thigh Muscles Using Magnetic Resonance Elastography

BACKGROUND

Pathologies of the muscles can manifest different physiological and functional changes. To adapt treatment, it is necessary to characterize the elastic property (shear modulus) of single muscles. Previous studies have used magnetic resonance elastography (MRE), a technique based on MRI technology, to analyze the mechanical behavior of healthy and pathological muscles. The purpose of this study was to develop protocols using MRE to determine the shear modulus of nine thigh muscles at rest.

METHODS

Twenty-nine healthy volunteers (mean age = 26 ± 3.41 years) with no muscle abnormalities underwent MRE tests (1.5 T MRI). Five MRE protocols were developed to quantify the shear moduli of the nine following thigh muscles at rest: rectus femoris (RF), vastus medialis (VM), vastus intermedius (VI), vastus lateralis (VL), sartorius (Sr), gracilis (Gr), semimembranosus (SM), semitendinosus (ST), and biceps (BC). In addition, the shear modulus of the subcutaneous adipose tissue was analyzed.

RESULTS

The gracilis, sartorius, and semitendinosus muscles revealed a significantly higher shear modulus (μ_Gr = 6.15 ± 0.45 kPa, μ_ Sr = 5.15 ± 0.19 kPa, and μ_ ST = 5.32 ± 0.10 kPa, respectively) compared to other tissues (from μ_ RF = 3.91 ± 0.16 kPa to μ_VI = 4.23 ± 0.25 kPa). Subcutaneous adipose tissue had the lowest value (μ_adipose tissue = 3.04 ± 0.12 kPa) of all the tissues tested.

CONCLUSION

The different elasticities measured between the tissues may be due to variations in the muscles' physiological and architectural compositions. Thus, the present protocol could be applied to injured muscles to identify their behavior of elastic property. Previous studies on muscle pathology found that quantification of the shear modulus could be used as a clinical protocol to identify pathological muscles and to follow-up effects of treatments and therapies. These data could also be used for modelling purposes.

Ex vivo engineered immune organoids for controlled germinal center reactions

Ex vivo engineered three-dimensional organotypic cultures have enabled the real-time study and control of biological functioning of mammalian tissues. Organs of broad interest where its architectural, cellular, and molecular complexity has prevented progress in ex vivo engineering are the secondary immune organs. Ex vivo immune organs can enable mechanistic understanding of the immune system and more importantly, accelerate the translation of immunotherapies as well as a deeper understanding of the mechanisms that lead to their malignant transformation into a variety of B and T cell malignancies. However, till date, no modular ex vivo immune organ has been developed with an ability to control the rate of immune reaction through tunable design parameter. Here we describe a B cell follicle organoid made of nanocomposite biomaterials, which recapitulates the anatomical microenvironment of a lymphoid tissue that provides the basis to induce an accelerated germinal center (GC) reaction by continuously providing extracellular matrix (ECM) and cell-cell signals to naïve B cells. Compared to existing co-cultures, immune organoids provide a control over primary B cell proliferation with ∼100-fold higher and rapid differentiation to the GC phenotype with robust antibody class switching.

Multifrequency time-harmonic elastography for the measurement of liver viscoelasticity in large tissue windows

Elastography of the liver for the non-invasive diagnosis of hepatic fibrosis is an established method. However, investigations of obese patients or patients with ascites are often limited by small and superficial elastographic windows. Therefore, multifrequency time-harmonic elastography (THE) based on time-resolved A-line ultrasound has recently been developed for measuring liver viscoelasticity in wide soft tissue windows and at greater depths. In this study, THE was integrated into a clinical B-mode scanner connected to a dedicated actuator bed driven by superimposed vibrations of 30- to 60-Hz frequencies. The resulting shear waves in the liver were captured along multiple profiles 7 to 14 cm in width and automatically processed for reconstruction of mean efficient shear wave speed and shear wave dispersion slope. This new modality was tested in healthy volunteers and 22 patients with clinically proven cirrhosis. Patients could be separated from controls by higher shear wave speeds (3.11 ± 0.64 m/s, 2.14-4.81 m/s, vs. 1.74 ± 0.10 m/s, 1.60-1.91 m/s) without significant degradation of data by high body mass index or ascites. Furthermore, the wave speed dispersion slope was significantly (p = 0.0025) lower in controls (5.2 ± 1.8 m/s/kHz) than in patients (39.1 ± 32.2 m/s/kHz). In conclusion, THE is useful for the diagnosis of cirrhosis in large tissue windows.

Ultra wideband (0.5-16 kHz) MR elastography for robust shear viscoelasticity model identification

Changes in the viscoelastic parameters of soft biological tissues often correlate with progression of disease, trauma or injury, and response to treatment. Identifying the most appropriate viscoelastic model, then estimating and monitoring the corresponding parameters of that model can improve insight into the underlying tissue structural changes. MR Elastography (MRE) provides a quantitative method of measuring tissue viscoelasticity. In a previous study by the authors (Yasar et al 2013 Magn. Reson. Med. 70 479-89), a silicone-based phantom material was examined over the frequency range of 200 Hz-7.75 kHz using MRE, an unprecedented bandwidth at that time. Six viscoelastic models including four integer order models and two fractional order models, were fit to the wideband viscoelastic data (measured storage and loss moduli as a function of frequency). The 'fractional Voigt' model (spring and springpot in parallel) exhibited the best fit and was even able to fit the entire frequency band well when it was identified based only on a small portion of the band. This paper is an extension of that study with a wider frequency range from 500 Hz to 16 kHz. Furthermore, more fractional order viscoelastic models are added to the comparison pool. It is found that added complexity of the viscoelastic model provides only marginal improvement over the 'fractional Voigt' model. And, again, the fractional order models show significant improvement over integer order viscoelastic models that have as many or more fitting parameters.

In vivo multifrequency magnetic resonance elastography of the human intervertebral disk

PURPOSE

To test in vivo magnetic resonance elastography (MRE) of the human intervertebral disk (IVD).

METHODS

The feasibility of MRE in IVD was demonstrated in ex vivo bovine disks. Sixteen asymptomatic volunteers underwent multifrequency MRE of the lumbar spine (IVD L3/4 and L4/5, n = 32) using a posterior plate transducer connected to a loudspeaker and operated at five frequencies from 50 to 70 Hz. Full wave field data were acquired in 10 transverse slices of 2 × 2 × 2 mm(3) resolution. High-resolution maps of magnitude |G*| and phase angle φ of complex shear modulus G* were generated by multifrequency dual elasto visco (MDEV) inversion. Disk morphology was assessed by the Pfirrmann score (Pf).

RESULTS

Morphological Pf was 1 in 25, 2 in 3, and 3 in 4 disks. |G*| decreased with Pf by a Pearson's linear correlation coefficient of R = -0.592 (P = 0.0004), while φ remained unchanged. Group mean mechanical parameters for Pf = 1 to 3 were |G*| = 6.51 ± 1.27, 5.29 ± 0.95, 4.03 ± 0.99 kPa, and φ = 1.190 ± 0.181, 1.170 ± 0.156, 1.088 ± 0.084 rad, respectively (p[Pf1-Pf3] < 0.001). The variability of mechanical parameters in one volunteer including diurnal changes was approximately 11%.

CONCLUSION

Multifrequency MRE with MDEV inversion allows measurement of in vivo mechanical properties of IVDs and may provide additional information in disc degeneration beyond standard morphological changes.

High-resolution mechanical imaging of glioblastoma by multifrequency magnetic resonance elastography

Objective

To generate high-resolution maps of the viscoelastic properties of human brain parenchyma for presurgical quantitative assessment in glioblastoma (GB).

Methods

Twenty-two GB patients underwent routine presurgical work-up supplemented by additional multifrequency magnetic resonance elastography. Two three-dimensional viscoelastic parameter maps, magnitude |G*|, and phase angle φ of the complex shear modulus were reconstructed by inversion of full wave field data in 2-mm isotropic resolution at seven harmonic drive frequencies ranging from 30 to 60 Hz.

Results

Mechanical brain maps confirmed that GB are composed of stiff and soft compartments, resulting in high intratumor heterogeneity. GB could be easily differentiated from healthy reference tissue by their reduced viscous behavior quantified by φ (0.37±0.08 vs. 0.58±0.07). |G*|, which in solids more relates to the material's stiffness, was significantly reduced in GB with a mean value of 1.32±0.26 kPa compared to 1.54±0.27 kPa in healthy tissue (P = 0.001). However, some GB (5 of 22) showed increased stiffness.

Conclusion

GB are generally less viscous and softer than healthy brain parenchyma. Unrelated to the morphology-based contrast of standard magnetic resonance imaging, elastography provides an entirely new neuroradiological marker and contrast related to the biomechanical properties of tumors.

In vivo high-resolution magnetic resonance elastography of the uterine corpus and cervix

OBJECTIVES

To apply 3D multifrequency MR elastography (3DMMRE) to the uterus and analyse the viscoelasticity of the uterine tissue in healthy volunteers considering individual variations and variations over the menstrual cycle.

METHODS

Sixteen healthy volunteers participated in the study, one of whom was examined 12 times over two menstrual cycles. Pelvic 3DMMRE was performed on a 1.5-T scanner with seven vibration frequencies (30-60 Hz) using a piezoelectric driver. Two mechanical parameter maps were obtained corresponding to the magnitude (|G (*) |) and the phase angle (φ) of the complex shear modulus.

RESULTS

On average, the uterine corpus had higher elasticity, but similar viscosity compared with the cervix, reflected by |G (*) |uterine corpus = 2.58 ± 0.52 kPa vs. |G (*) |cervix = 2.00 ± 0.34 kPa (p < 0.0001) and φ uterine corpus = 0.54 ± 0.08, φ cervix = 0.57 ± 0.12 (p = 0.428). With 2.23 ± 0.26 kPa, |G (*) | of the myometrium was lower in the secretory phase (SP) compared with that of the proliferative phase (PP, |G (*) | = 3.01 ± 0.26 kPa). For the endometrium, the value of |G (*) | in SP was 68% lower than during PP (PP, |G (*) | = 3.34 ± 0.42 kPa; SP, |G (*) | = 1.97 ± 0.34 kPa; p = 0.0061).

CONCLUSION

3DMMRE produces high-resolution mechanical parameter maps of the uterus and cervix and shows sensitivity to structural and functional changes of the endometrium and myometrium during the menstrual cycle.

KEY POINTS

MR elastography provided for the first time spatially resolved viscoelasticity maps of uterus. Uterine corpus had a higher elasticity, but similar viscosity compared with cervix. The stiffness of both endometrium and myometrium decreases during the menstrual cycle.

High Resolution Imaging of Viscoelastic Properties of Intracranial Tumours by Multi-Frequency Magnetic Resonance Elastography

PURPOSE

In recent years Magnetic Resonance Elastography (MRE) emerged into a clinically applicable imaging technique. It has been shown that MRE is capable of measuring global changes of the viscoelastic properties of cerebral tissue. The purpose of our study was to evaluate a spatially resolved three-dimensional multi-frequent MRE (3DMMRE) for assessment of the viscoelastic properties of intracranial tumours.

METHODS

A total of 27 patients (63 ± 13 years) were included. All examinations were performed on a 3.0 T scanner, using a modified phase-contrast echo planar imaging sequence. We used 7 vibration frequencies in the low acoustic range with a temporal resolution of 8 dynamics per wave cycle. Post-processing included multi-frequency dual elasto-visco (MDEV) inversion to generate high-resolution maps of the magnitude |G*| and the phase angle φ of the complex valued shear modulus.

RESULTS

The tumour entities included in this study were: glioblastoma (n = 11), anaplastic astrocytoma (n = 3), meningioma (n = 7), cerebral metastasis (n = 5) and intracerebral abscess formation (n = 1). Primary brain tumours and cerebral metastases were not distinguishable in terms of |G*| and φ. Glioblastoma presented the largest range of |G*| values and a trend was delineable that glioblastoma were slightly softer than WHO grade III tumours. In terms of φ, meningiomas were clearly distinguishable from all other entities.

CONCLUSIONS

In this pilot study, while analysing the viscoelastic constants of various intracranial tumour entities with an improved spatial resolution, it was possible to characterize intracranial tumours by their mechanical properties. We were able to clearly delineate meningiomas from intraaxial tumours, while for the latter group an overlap remains in viscoelastic terms.

Magnetic resonance elastography of liver: clinical applications

Magnetic resonance elastography (MRE) has been successfully implemented in the assessment of diffuse liver diseases. Currently, MRE is the most accurate noninvasive technique for detection and staging of liver fibrosis with a potential to replace liver biopsy. Magnetic resonance elastography is able to differentiate isolated fatty liver disease from steatohepatitis with or without fibrosis. Potential clinical applications include the differentiation of benign and malignant focal liver masses and the assessment of treatment response in diffuse liver diseases.

Sample interval modulation for the simultaneous acquisition of displacement vector data in magnetic resonance elastography: theory and application

SampLe Interval Modulation-magnetic resonance elastography (SLIM-MRE) is introduced for simultaneously encoding all three displacement projections of a monofrequency vibration into the MR signal phase. In SLIM-MRE, the individual displacement components are observed using different sample intervals. In doing so, the components are modulated with different apparent frequencies in the MR signal phase expressed as a harmonic function of the start time of the motion encoding gradients and can thus be decomposed by applying a Fourier transform to the sampled multidirectional MR phases. In this work, the theoretical foundations of SLIM-MRE are presented and the new idea is implemented using a high field (11.7 T) vertical bore magnetic resonance imaging system on an inhomogeneous agarose gel phantom sample. The local frequency estimation-derived stiffness values were the same within the error margins for both the new SLIM-MRE method and for conventional MRE, while the number of temporally-resolved MRE experiments needed for each study was reduced from three to one. In this work, we present for the first time, monofrequency displacement data along three sensitization directions that were acquired simultaneously and stored in the same k-space.

Towards an elastographic atlas of brain anatomy

Cerebral viscoelastic constants can be measured in a noninvasive, image-based way by magnetic resonance elastography (MRE) for the detection of neurological disorders. However, MRE brain maps of viscoelastic constants are still limited by low spatial resolution. Here we introduce three-dimensional multifrequency MRE of the brain combined with a novel reconstruction algorithm based on a model-free multifrequency inversion for calculating spatially resolved viscoelastic parameter maps of the human brain corresponding to the dynamic range of shear oscillations between 30 and 60 Hz. Maps of two viscoelastic parameters, the magnitude and the phase angle of the complex shear modulus, |G*| and φ, were obtained and normalized to group templates of 23 healthy volunteers in the age range of 22 to 72 years. This atlas of the anatomy of brain mechanics reveals a significant contrast in the stiffness parameter |G*| between different anatomical regions such as white matter (WM; 1.252±0.260 kPa), the corpus callosum genu (CCG; 1.104±0.280 kPa), the thalamus (TH; 1.058±0.208 kPa) and the head of the caudate nucleus (HCN; 0.649±0.101 kPa). φ, which is sensitive to the lossy behavior of the tissue, was in the order of CCG (1.011±0.172), TH (1.037±0.173), CN (0.906±0.257) and WM (0.854±0.169). The proposed method provides the first normalized maps of brain viscoelasticity with anatomical details in subcortical regions and provides useful background data for clinical applications of cerebral MRE.