Multifrequency magnetic resonance elastography-based tomoelastography of the parotid glands–feasibility and reference values

Objectives: Accurate radiological differentiation of parotid tumors remains challenging despite recent technical advances in quantitative medical imaging. Multifrequency magnetic resonance elastography (MRE) could provide additional information on viscoelastic properties of normal and abnormal biological tissues. This study investigates the feasibility of MRE of the parotid glands in healthy participants and provides first reference values. Methods: 20 healthy participants underwent multifrequency MRE of both parotid glands at 3 Tesla. Shear waves at frequencies of 25, 30, 40, and 50 Hz were introduced into the participants' heads through the occiput using pressurized-air actuators. Shear wave speed (SWS) and loss angle of the shear modulus (φ) were reconstructed by tomoelastography post-processing as surrogate parameters for tissue stiffness and viscosity or fluidity. 10 participants underwent repeated MRE to determine test-retest reliability based on intraclass correlation coefficients. Results: All MRE datasets acquired could be included in the analysis. Mean SWS was 0.97 ± 0.13 m/s, and mean φ was 0.59 ± 0.05 rad, each for both sides combined and without notable lateral difference (p = 0.88/0.87). Test-retest reliability was good for SWS (ICC = 0.84 for both sides/ICC = 0.77 for the right side/ICC = 0.79 for the left side) and good to excellent for φ(ICC = 0.94/0.86/0.90). Conclusions: Multifrequency MRE of the parotid glands is feasible and reliable. This technique, therefore, is a promising method for investigating the viscoelastic properties of salivary gland tumors in future studies

Spatial heterogeneity of hepatic fibrosis in primary sclerosing cholangitis vs. viral hepatitis assessed by MR elastography

Spatial heterogeneity of hepatic fibrosis in primary sclerosing cholangitis (PSC) in comparison to viral hepatitis was assessed as a potential new biomarker using MR elastography (MRE). In this proof-of-concept study, we hypothesized a rather increased heterogeneity in PSC and a rather homogeneous distribution in viral hepatitis. Forty-six consecutive subjects (PSC: n=20, viral hepatitis: n=26) were prospectively enrolled between July 2014 and April 2017. Subjects underwent multifrequency MRE (1.5 T) using drive frequencies of 35-60 Hz and generating shear-wave speed (SWS in m/s) maps as a surrogate of stiffness. The coefficient of variation (CV in %) was determined to quantify fibrosis heterogeneity. Mean SWS and CV were 1.70 m/s and 21% for PSC, and 1.84 m/s and 18% for viral hepatitis. Fibrosis heterogeneity was significantly increased for PSC (P=0.04) while no difference was found for SWS of PSC and viral hepatitis (P=0.17). Global hepatic stiffness was similar in PSC and viral hepatitis groups, but spatial heterogeneity may reveal spatial patterns of stiffness changes towards enhanced biophysics-based diagnosis by MRI

Cerebral Ultrasound Time-Harmonic Elastography Reveals Softening of the Human Brain Due to Dehydration

Hydration influences blood volume, blood viscosity, and water content in soft tissues - variables that determine the biophysical properties of biological tissues including their stiffness. In the brain, the relationship between hydration and stiffness is largely unknown despite the increasing importance of stiffness as a quantitative imaging marker. In this study, we investigated cerebral stiffness (CS) in 12 healthy volunteers using ultrasound time-harmonic elastography (THE) in different hydration states: (i) during normal hydration, (ii) after overnight fasting, and (iii) within 1 h of drinking 12 ml of water per kg body weight. In addition, we correlated shear wave speed (SWS) with urine osmolality and hematocrit. SWS at normal hydration was 1.64 ± 0.02 m/s and decreased to 1.57 ± 0.04 m/s (p < 0.001) after overnight fasting. SWS increased again to 1.63 ± 0.01 m/s within 30 min of water drinking, returning to values measured during normal hydration (p = 0.85). Urine osmolality at normal hydration (324 ± 148 mOsm/kg) increased to 784 ± 107 mOsm/kg (p < 0.001) after fasting and returned to normal (288 ± 128 mOsm/kg, p = 0.83) after water drinking. SWS and urine osmolality correlated linearly (r = -0.68, p < 0.001), while SWS and hematocrit did not correlate (p = 0.31). Our results suggest that mild dehydration in the range of diurnal fluctuations is associated with significant softening of brain tissue, possibly due to reduced cerebral perfusion. To ensure consistency of results, it is important that cerebral elastography with a standardized protocol is performed during normal hydration

Fully automated quantification of in vivo viscoelasticity of prostate zones using magnetic resonance elastography with Dense U-net segmentation

Magnetic resonance elastography (MRE) for measuring viscoelasticity heavily depends on proper tissue segmentation, especially in heterogeneous organs such as the prostate. Using trained network-based image segmentation, we investigated if MRE data suffice to extract anatomical and viscoelastic information for automatic tabulation of zonal mechanical properties of the prostate. Overall, 40 patients with benign prostatic hyperplasia (BPH) or prostate cancer (PCa) were examined with three magnetic resonance imaging (MRI) sequences: T2-weighted MRI (T2w), diffusion-weighted imaging (DWI), and MRE-based tomoelastography yielding six independent sets of imaging data per patient (T2w, DWI, apparent diffusion coefficient (ADC), MRE magnitude, shear wave speed, and loss angle maps). Combinations of these data were used to train Dense U-nets with manually segmented masks of the entire prostate gland (PG), central zone (CZ), and peripheral zone (PZ) in 30 patients and to validate them in 10 patients. Dice score (DS), sensitivity, specificity, and Hausdorff distance were determined. We found that segmentation based on MRE magnitude maps alone (DS, PG: 0.93$\pm$0.04, CZ: 0.95$\pm$0.03, PZ: 0.77$\pm$0.05) was more accurate than magnitude maps combined with T2w and DWI_b (DS, PG: 0.91$\pm$0.04, CZ: 0.91$\pm$0.06, PZ: 0.63$\pm$0.16) or T2w alone (DS, PG: 0.92$\pm$0.03, CZ: 0.91$\pm$0.04, PZ: 0.65$\pm$0.08). Automatically tabulated MRE values were not different from ground-truth values (P>0.05). In conclusion: MRE combined with Dense U-net segmentation allows tabulation of quantitative imaging markers without manual analysis and independent of other MRI sequences and can thus contribute to PCa detection and classification

Microscopic multifrequency MR elastography for mapping viscoelasticity in zebrafish

Purpose: The zebrafish (Danio rerio) has become an important animal model in a wide range of biomedical research disciplines. Growing awareness of the role of biomechanical properties in tumor progression and neuronal development has led to an increasing interest in the noninvasive mapping of the viscoelastic properties of zebrafish by elastography methods applicable to bulky and nontranslucent tissues. Methods: Microscopic multifrequency MR elastography is introduced for mapping shear wave speed (SWS) and loss angle (φ) as markers of stiffness and viscosity of muscle, brain, and neuroblastoma tumors in postmortem zebrafish with 60 µm in-plane resolution. Experiments were performed in a 7 Tesla MR scanner at 1, 1.2, and 1.4 kHz driving frequencies. Results: Detailed zebrafish viscoelasticity maps revealed that the midbrain region (SWS = 3.1 ± 0.7 m/s, φ = 1.2 ± 0.3 radian [rad]) was stiffer and less viscous than telencephalon (SWS = 2.6 ± 0. 5 m/s, φ = 1.4 ± 0.2 rad) and optic tectum (SWS = 2.6 ± 0.5 m/s, φ = 1.3 ± 0.4 rad), whereas the cerebellum (SWS = 2.9 ± 0.6 m/s, φ = 0.9 ± 0.4 rad) was stiffer but less viscous than both (all p < .05). Overall, brain tissue (SWS = 2.9 ± 0.4 m/s, φ = 1.2 ± 0.2 rad) had similar stiffness but lower viscosity values than muscle tissue (SWS = 2.9 ± 0.5 m/s, φ = 1.4 ± 0.2 rad), whereas neuroblastoma (SWS = 2.4 ± 0.3 m/s, φ = 0.7 ± 0.1 rad, all p < .05) was the softest and least viscous tissue. Conclusion: Microscopic multifrequency MR elastography-generated maps of zebrafish show many details of viscoelasticity and resolve tissue regions, of great interest in neuromechanical and oncological research and for which our study provides first reference values

Mechanical properties of murine hippocampal subregions investigated by atomic force microscopy and in vivo magnetic resonance elastography

The hippocampus is a very heterogeneous brain structure with different mechanical properties reflecting its functional variety. In particular, adult neurogenesis in rodent hippocampus has been associated with specific viscoelastic properties in vivo and ex vivo. Here, we study the microscopic mechanical properties of hippocampal subregions using ex vivo atomic force microscopy (AFM) in correlation with the expression of GFP in presence of the nestin promoter, providing a marker of neurogenic activity. We further use magnetic resonance elastography (MRE) to investigate whether in vivo mechanical properties reveal similar spatial patterns, however, on a much coarser scale. AFM showed that tissue stiffness increases with increasing distance from the subgranular zone (p = 0.0069), and that stiffness is 39% lower in GFP than non-GFP regions (p = 0.0004). Consistently, MRE showed that dentate gyrus is, on average, softer than Ammon's horn (shear wave speed = 3.2 ± 0.2 m/s versus 4.4 ± 0.3 m/s, p = 0.01) with another 3.4% decrease towards the subgranular zone (p = 0.0001). The marked reduction in stiffness measured by AFM in areas of high neurogenic activity is consistent with softer MRE values, indicating the sensitivity of macroscopic mechanical properties in vivo to micromechanical structures as formed by the neurogenic niche of the hippocampus

Cerebral tomoelastography based on multifrequency MR elastography in two and three dimensions

Magnetic resonance elastography (MRE) generates quantitative maps of the mechanical properties of biological soft tissues. However, published values obtained by brain MRE vary largely and lack detail resolution, due to either true biological effects or technical challenges. We here introduce cerebral tomoelastography in two and three dimensions for improved data consistency and detail resolution while considering aging, brain parenchymal fraction (BPF), systolic blood pressure, and body-mass-index. Multifrequency MRE with 2D- and 3D-tomoelastography postprocessing was applied to the brains of 31 volunteers (age range: 22-61 years) for analyzing the coefficient of variation (CV) and effects of biological factors. Eleven volunteers were rescanned after one day and one year to determine intraclass correlation coefficient (ICC) and identify possible long-term changes. White matter shear-wave-speed (SWS) was slightly higher in 2D-MRE (1.28±0.02m/s) than 3D-MRE (1.22±0.05m/s, p<0.0001), with less variation after one day in 2D (0.33±0.32%) than in 3D (0.96±0.66%, p=0.004), which was also reflected in a slightly lower CV and higher ICC in 2D (1.84%, 0.97 [0.88-0.99]) than in 3D (3.89%, 0.95 [0.76-0.99]). Remarkably, 3D-MRE was sensitive to a decrease in white matter SWS within only one year, whereas no change in white matter volume was observed during this follow-up period. Across volunteers, stiffness correlated with age and BPF, but not with blood pressure and body-mass-index. Cerebral tomoelastography provides high-resolution viscoelasticity maps with excellent consistency. Brain MRE in 2D shows less variation across volunteers in shorter scan times than 3D-MRE, while 3D-MRE appears to be more sensitive to subtle biological effects such as aging

Adipose cells and tissues soften with lipid accumulation while in diabetes adipose tissue stiffens

Adipose tissue expansion involves both differentiation of new precursors and size increase of mature adipocytes. While the two processes are well balanced in healthy tissues, obesity and diabetes type II are associated with abnormally enlarged adipocytes and excess lipid accumulation. Previous studies suggested a link between cell stiffness, volume and stem cell differentiation, although in the context of preadipocytes, there have been contradictory results regarding stiffness changes with differentiation. Thus, we set out to quantitatively monitor adipocyte shape and size changes with differentiation and lipid accumulation. We quantified by optical diffraction tomography that differentiating preadipocytes increased their volumes drastically. Atomic force microscopy (AFM)-indentation and -microrheology revealed that during the early phase of differentiation, human preadipocytes became more compliant and more fluid-like, concomitant with ROCK-mediated F-actin remodelling. Adipocytes that had accumulated large lipid droplets were more compliant, and further promoting lipid accumulation led to an even more compliant phenotype. In line with that, high fat diet-induced obesity was associated with more compliant adipose tissue compared to lean animals, both for drosophila fat bodies and murine gonadal adipose tissue. In contrast, adipose tissue of diabetic mice became significantly stiffer as shown not only by AFM but also magnetic resonance elastography. Altogether, we dissect relative contributions of the cytoskeleton and lipid droplets to cell and tissue mechanical changes across different functional states, such as differentiation, nutritional state and disease. Our work therefore sets the basis for future explorations on how tissue mechanical changes influence the behaviour of mechanosensitive tissue-resident cells in metabolic disorders

Multiple motion encoding in Phase-Contrast MRI:A general theory and application to elastography imaging

While MRI allows to encode the motion of tissue in the magnetization’s phase, it remains yet a challenge to obtain high fidelity motion images due to wraps in the phase for high encoding efficiencies. Therefore, we propose an optimal multiple motion encoding method (OMME) and exemplify it in Magnetic Resonance Elastography (MRE) data. OMME is formulated as a non-convex least-squares problem for the motion using an arbitrary number of phase-contrast measurements with different motion encoding gradients (MEGs). The mathematical properties of OMME are proved in terms of standard deviation and dynamic range of the motion’s estimate for arbitrary MEGs combination which are confirmed using synthetically generated data. OMME’s performance is assessed on MRE data from in vivo human brain experiments and compared to dual encoding strategies. The unwrapped images are further used to reconstruct stiffness maps and compared to the ones obtained using conventional unwrapping methods. OMME allowed to successfully combine several MRE phase images with different MEGs, outperforming dual encoding strategies in either motion-to-noise ratio (MNR) or number of successfully reconstructed voxels with good noise stability. This lead to stiffness maps with greater resolution of details than obtained with conventional unwrapping methods. The proposed OMME method allows for a flexible and noise robust increase in the dynamic range and thus provides wrap-free phase images with high MNR. In MRE, the method may be especially suitable when high resolution images with high MNR are needed

Optical time-harmonic elastography for multiscale stiffness mapping across the phylogenetic tree

Rapid mapping of the mechanical properties of soft biological tissues from light microscopy to macroscopic imaging could transform fundamental biophysical research by providing clinical biomarkers to complement in vivo elastography. We here introduce superfast optical time-harmonic elastography (OTHE) to remotely encode surface and subsurface shear wave fields for generating maps of tissue stiffness with unprecedented detail resolution. OTHE rigorously exploits the space-time propagation characteristics of time-harmonic waves to address current limitations of biomechanical imaging and elastography. Key solutions are presented for stimulation, encoding, and stiffness reconstruction of time-harmonic, multifrequency shear waves, all tuned to provide consistent stiffness values across resolutions from microns to millimeters. OTHE's versatility is demonstrated in Bacillus subtilis biofilms, zebrafish embryos, adult zebrafish, and human skeletal muscle, reflecting the diversity of the phylogenetic tree from a mechanics perspective. By zooming in on stiffness details from coarse to finer scales, OTHE advances developmental biology and offers a way to perform biomechanics-based tissue histology that consistently matches in vivo time-harmonic elastography in patients