Liquid-Liver Phantom: Mimicking the Viscoelastic Dispersion of Human Liver for Ultrasound- and MRI-Based Elastography.

OBJECTIVES: Tissue stiffness can guide medical diagnoses and is exploited as an imaging contrast in elastography. However, different elastography devices show different liver stiffness values in the same subject, hindering comparison of values and establishment of system-independent thresholds for disease detection. There is a need for standardized phantoms that specifically address the viscosity-related dispersion of stiffness over frequency. To improve standardization of clinical elastography across devices and platforms including ultrasound and magnetic resonance imaging (MRI), a comprehensively characterized phantom is introduced that mimics the dispersion of stiffness of the human liver and can be generated reproducibly. MATERIALS AND METHODS: The phantom was made of linear polymerized polyacrylamide (PAAm) calibrated to the viscoelastic properties of healthy human liver in vivo as reported in the literature. Stiffness dispersion was analyzed using the 2-parameter springpot model fitted to the dispersion of shear wave speed of PAAm, which was measured by shear rheometry, ultrasound-based time-harmonic elastography, clinical magnetic resonance elastography (MRE), and tabletop MRE in the frequency range of 5 to 3000 Hz. Imaging parameters for ultrasound and MRI, reproducibility, aging behavior, and temperature dependency were assessed. In addition, the frequency bandwidth of shear wave speed of clinical elastography methods (Aplio i900, Canon; Acuson Sequoia, Siemens; FibroScan, EchoSense) was characterized. RESULTS: Within the entire frequency range analyzed in this study, the PAAm phantom reproduced well the stiffness dispersion of human liver in vivo despite its fluid properties under static loading (springpot stiffness parameter, 2.14 [95% confidence interval, 2.08-2.19] kPa; springpot powerlaw exponent, 0.367 [95% confidence interval, 0.362-0.373]). Imaging parameters were close to those of liver in vivo with only slight variability in stiffness values of 0.5% (0.4%, 0.6%), 4.1% (3.9%, 4.5%), and -0.63% (-0.67%, -0.58%), respectively, between batches, over a 6-month period, and per °C increase in temperature. CONCLUSIONS: The liquid-liver phantom has useful properties for standardization and development of liver elastography. First, it can be used across clinical and experimental elastography devices in ultrasound and MRI. Second, being a liquid, it can easily be adapted in size and shape to specific technical requirements, and by adding inclusions and scatterers. Finally, because the phantom is based on noncrosslinked linear PAAm constituents, it is easy to produce, indicating potential widespread use among researchers and vendors to standardize liver stiffness measurements.

Real-Time Multifrequency MR Elastography of the Human Brain Reveals Rapid Changes in Viscoelasticity in Response to the Valsalva Maneuver.

Modulation of cerebral blood flow and vascular compliance plays an important role in the regulation of intracranial pressure (ICP) and also influences the viscoelastic properties of brain tissue. Therefore, magnetic resonance elastography (MRE), the gold standard for measuring in vivo viscoelasticity of brain tissue, is potentially sensitive to cerebral autoregulation. In this study, we developed a multifrequency MMRE technique that provides serial maps of viscoelasticity at a frame rate of nearly 6 Hz without gating, i.e., in quasi-real time (rt-MMRE). This novel method was used to monitor rapid changes in the viscoelastic properties of the brains of 17 volunteers performing the Valsalva maneuver (VM). rt-MMRE continuously sampled externally induced vibrations comprising three frequencies of 30.03, 30.91, and 31.8 Hz were over 90 s using a steady-state, spiral-readout gradient-echo sequence. Data were processed by multifrequency dual elasto-visco (MDEV) inversion to generate maps of magnitude shear modulus | G∗| (stiffness) and loss angle φ at a frame rate of 5.4 Hz. As controls, the volunteers were examined to study the effects of breath-hold following deep inspiration and breath-hold following expiration. We observed that | G∗| increased while φ decreased due to VM and, less markedly, due to breath-hold in inspiration. Group mean VM values showed an early overshoot of | G∗| 2.4 ± 1.2 s after the onset of the maneuver with peak values of 6.7 ± 4.1% above baseline, followed by a continuous increase in stiffness during VM. A second overshoot of | G∗| occurred 5.5 ± 2.0 s after the end of VM with peak values of 7.4 ± 2.8% above baseline, followed by 25-s sustained recovery until the end of image acquisition. φ was constantly reduced by approximately 2% during the entire VM without noticeable peak values. This is the first report of viscoelasticity changes in brain tissue induced by physiological maneuvers known to alter ICP and detected by clinically applicable rt-MMRE. Our results show that apnea and VM slightly alter brain properties toward a more rigid-solid behavior. Overshooting stiffening reactions seconds after onset and end of VM reveal rapid autoregulatory processes of brain tissue viscoelasticity.

Superviscous properties of the in vivo brain at large scales.

There is growing awareness that brain mechanical properties are important for neural development and health. However, published values of brain stiffness differ by orders of magnitude between static measurements and in vivo magnetic resonance elastography (MRE), which covers a dynamic range over several frequency decades. We here show that there is no fundamental disparity between static mechanical tests and in vivo MRE when considering large-scale properties, which encompass the entire brain including fluid filled compartments. Using gradient echo real-time MRE, we investigated the viscoelastic dispersion of the human brain in, so far, unexplored dynamic ranges from intrinsic brain pulsations at 1 Hz to ultralow-frequency vibrations at 5, 6.25, 7.8 and 10 Hz to the normal frequency range of MRE of 40 Hz. Surprisingly, we observed variations in brain stiffness over more than two orders of magnitude, suggesting that the in vivo human brain is superviscous on large scales with very low shear modulus of 42±13 Pa and relatively high viscosity of 6.6±0.3 Pa∙s according to the two-parameter solid model. Our data shed light on the crucial role of fluid compartments including blood vessels and cerebrospinal fluid (CSF) for whole brain properties and provide, for the first time, an explanation for the variability of the mechanical brain responses to manual palpation, local indentation, and high-dynamic tissue stimulation as used in elastography.

Steady-State Multifrequency Magnetic Resonance Elastography of the Thoracic and Abdominal Human Aorta-Validation and Reference Values.

OBJECTIVES: The aim of this study was to investigate the potential of stroboscopic-wavefield-sampling-based multifrequency magnetic resonance elastography (sMRE) for quantifying the stiffness of the human thoracic and abdominal aorta in vivo. MATERIALS AND METHODS: The sMRE of the thoracic and abdominal aorta was performed at 1.5 T field strength in 20 healthy volunteers aged 27 to 77 years (3 women; median age, 33 years; interquartile range [IQR], 16 years). Compound maps of shear wave speed (SWS) were reconstructed and evaluated during the diastolic phase in 3 anatomical regions: ascending thoracic aorta (AA), descending thoracic aorta (AD), and abdominal aorta (AAb). The SWS maps were read by 2 readers. Blood pressure and pulse wave velocity were determined noninvasively before sMRE. Data are given as median (IQR) and were compared using the Kruskal-Wallis and Wilcoxon rank sum tests. Intraclass correlation was used to determine interobserver and intraobserver agreement, as well as reproducibility. Multiple linear regression analysis was performed to evaluate effects of age, sex, vessel diameter, blood pressure, pulse wave velocity, and aortic segment on measured SWS. RESULTS: All 20 participants underwent successful sMRE, resulting in a total of 60 aortic segments. The median SWS (IQR) of AA, AD, and AAb was 1.62 (0.16) m/s, 2.40 (0.24) m/s, and 2.48 (0.58) m/s, respectively. The SWS in AA was significantly lower (P < 0.001), and no differences in SWS (P = 0.67) were found between AD and AAb. Interobserver and intraobserver agreement, as well as reproducibility, was excellent, with intraclass correlation coefficients ranging between 0.957 and 0.998. A significant but weak influence of age on measured SWS was found, which increased from AA to AD and AAb (R = 0.229, 0.275, 0.377, respectively; P = 0.001-0.005). CONCLUSIONS: Quantification of aortic stiffness in different segments of the human aorta is possible with sMRE. Our results correlate well with known aortic stiffness differences in different anatomical locations and demonstrate the potential of sMRE for clinical stiffness measurement of the thoracoabdominal aorta, which may allow detection of physiological variation and cardiovascular diseases.

An analytical solution to the dispersion-by-inversion problem in magnetic resonance elastography.

PURPOSE: Magnetic resonance elastography (MRE) measures stiffness of soft tissues by analyzing their spatial harmonic response to externally induced shear vibrations. Many MRE methods use inversion-based reconstruction approaches, which invoke first- or second-order derivatives by finite difference operators (first- and second-FDOs) and thus give rise to a biased frequency dispersion of stiffness estimates. METHODS: We here demonstrate analytically, numerically, and experimentally that FDO-based stiffness estimates are affected by (1) noise-related underestimation of values in the range of high spatial wave support, that is, at lower vibration frequencies, and (2) overestimation of values due to wave discretization at low spatial support, that is, at higher vibration frequencies. RESULTS: Our results further demonstrate that second-FDOs are more susceptible to noise than first-FDOs and that FDO dispersion depends both on signal-to-noise ratio (SNR) and on a lumped parameter A, which is defined as wavelength over pixel size and over a number of pixels per stencil of the FDO. Analytical FDO dispersion functions are derived for optimizing A parameters at a given SNR. As a simple rule of thumb, we show that FDO artifacts are minimized when A/2 is in the range of the square root of 2SNR for the first-FDO or cubic root of 5SNR for the second-FDO. CONCLUSIONS: Taken together, the results of our study provide an analytical solution to a long-standing, well-recognized, yet unsolved problem in MRE postprocessing and might thus contribute to the ongoing quest for minimizing inversion artifacts in MRE.

Time-Resolved Response of Cerebral Stiffness to Hypercapnia in Humans.

Cerebral blood flow, cerebral stiffness (CS) and intracranial pressure are tightly linked variables of cerebrovascular reactivity and cerebral autoregulation. Transtemporal ultrasound time-harmonic elastography was used for rapid measurement of CS changes in 10 volunteers before, during and after administration of a gas mixture of 95% O2 and 5% CO2 (carbogen). Within the first 2.2 ± 2.0 min of carbogen breathing, shear wave speed determined as a surrogate parameter of CS increased from 1.57 ± 0.04 to 1.66 ± 0.05 m/s (p < 0.01) in synchrony with end-tidal CO2 while post-hypercapnic CS recovery was delayed by 2.7 ± 1.4 min in relation to end-tidal CO2. Our results indicate that CS is highly sensitive to changes in CO2 levels of inhaled air. Possible mechanisms underlying the observed CS changes might be associated with cerebrovascular reactivity, cerebral blood flow adaptation and intracranial regulation, all of which are potentially relevant for future diagnostic applications of transtemporal time-harmonic elastography in a wide spectrum of neurologic diseases.

Cardiac-gated steady-state multifrequency magnetic resonance elastography of the brain: Effect of cerebral arterial pulsation on brain viscoelasticity.

In-vivo brain viscoelasticity measured by magnetic resonance elastography (MRE) is a sensitive imaging marker for long-term biophysical changes in brain tissue due to aging and disease; however, it is still unknown whether MRE can reveal short-term periodic alterations of brain viscoelasticity related to cerebral arterial pulsation (CAP). We developed cardiac-gated steady-state MRE (ssMRE) with spiral readout and stroboscopic sampling of continuously induced mechanical vibrations in the brain at 20, 31.25, and 40 Hz frequencies. Maps of magnitude |G*| and phase ϕ of the complex shear modulus were generated by multifrequency dual visco-elasto inversion with a temporal resolution of 40 ms over 4 s. The method was tested in 12 healthy volunteers. During cerebral systole, |G*| decreased by 6.6 ± 1.9% (56 ± 22 Pa, p < 0.001, mean ± SD), whereas ϕ increased by 0.5 ± 0.5% (0.006 ± 0.005 rad, p = 0.002). The effect size of CAP-induced softening slightly decreased with age by 0.10 ± 0.05% per year (p = 0.04), indicating lower cerebral vascular compliance in older individuals. Our data show for the first time that the brain softens and becomes more viscous during systole, possibly due to an effect of CAP-induced arterial expansion and increased blood volume on effective-medium tissue properties. This sensitivity to vascular-solid tissue interactions makes ssMRE potentially useful for detection of cerebral vascular disease.

Real-time MR elastography for viscoelasticity quantification in skeletal muscle during dynamic exercises.

PURPOSE: To develop and test real-time MR elastography for viscoelastic parameter quantification in skeletal muscle during dynamic exercises. METHODS: In 15 healthy participants, 6 groups of lower-leg muscles (tibialis anterior, tibialis posterior, peroneus, extensor digitorum longus, soleus, gastrocnemius) were investigated by real-time MR elastography using a single-shot, steady-state spiral gradient-echo pulse sequence and stroboscopic undersampling of harmonic vibrations at 40 Hz frequency. One hundred and eighty consecutive maps of shear-wave speed and loss angle (φ) covering 30.6 s of total acquisition time at 5.9-Hz frame rate were reconstructed from 360 wave images encoding 2 in-plane wave components in an interleaved manner. The experiment was carried out twice to investigate 2 exercises-isometric plantar flexion and isometric dorsiflexion-each performed over 10 s between 2 resting periods. RESULTS: Activation of lower-extremity muscles was associated with increasing viscoelastic parameters shear-wave speed and φ, both reflecting properties related to the transverse direction relative to fiber orientation. Major viscoelastic changes were observed in soleus muscle during plantar flexion (shear-wave speed: 20.0% ± 3.6%, φ: 41.3% ± 12.0%) and in the tibialis anterior muscle during dorsiflexion (41.8% ± 10.2%, φ: 27.9% ± 2.8%; all P < .0001). Two of the muscles analyzed were significantly activated by plantar flexion and 4 by dorsiflexion based on shear-wave speed, whereas φ changed significantly in 5 muscles during both exercises. CONCLUSION: Real-time MR elastography allows mapping of dynamic, nonperiodic viscoelasticity changes in soft tissues such as voluntary muscle with high spatial and temporal resolution. Real-time MR elastography thus opens new horizons for the in vivo study of physiological processes in soft tissues toward functional elastography.

How tissue fluidity influences brain tumor progression.

Mechanical properties of biological tissues and, above all, their solid or fluid behavior influence the spread of malignant tumors. While it is known that solid tumors tend to have higher mechanical rigidity, allowing them to aggressively invade and spread in solid surrounding healthy tissue, it is unknown how softer tumors can grow within a more rigid environment such as the brain. Here, we use in vivo magnetic resonance elastography (MRE) to elucidate the role of anomalous fluidity for the invasive growth of soft brain tumors, showing that aggressive glioblastomas (GBMs) have higher water content while behaving like solids. Conversely, our data show that benign meningiomas (MENs), which contain less water than brain tissue, are characterized by fluid-like behavior. The fact that the 2 tumor entities do not differ in their soft properties suggests that fluidity plays an important role for a tumor's aggressiveness and infiltrative potential. Using tissue-mimicking phantoms, we show that the anomalous fluidity of neurotumors physically enables GBMs to penetrate surrounding tissue, a phenomenon similar to Saffman-Taylor viscous-fingering instabilities, which occur at moving interfaces between fluids of different viscosity. Thus, targeting tissue fluidity of malignant tumors might open horizons for the diagnosis and treatment of cancer.

Sensitivity of multifrequency magnetic resonance elastography and diffusion-weighted imaging to cellular and stromal integrity of liver tissue.

Microscopic structural alterations of liver tissue induced by freeze-thaw cycles give rise to palpable property changes. However, the underlying damage to tissue architecture is difficult to quantify histologically, and published data on macroscopic changes in biophysical properties are sparse. To better understand the influence of hepatic cells and stroma on global biophysical parameters, we studied rat liver specimens freshly taken (within 30 min after death) and treated by freeze-thaw cycles overnight at either -20 °C or -80 °C using diffusion-weighted imaging (DWI) and multifrequency magnetic resonance elastography (MRE) performed at 0.5 T in a tabletop MRE scanner. Tissue structure was analyzed histologically and rheologic data were analyzed using fractional order derivatives conceptualized by a called spring-pot component that interpolates between pure elastic and viscous responses. Overnight freezing and thawing induced membrane disruptions and cell detachment in the space of Disse, resulting in a markedly lower shear modulus µ and apparent diffusion coefficient (ADC) (µ[-20 °C] = 1.23 ±â€¯0.73 kPa, µ[-80 °C] = 0.66 ±â€¯0.75 kPa; ADC[-20 °C] = 0.649 ±â€¯0.028 µm2/s, ADC[-80 °C] = 0.626 ±â€¯0.025 µm2/s) compared to normal tissue (µâ€¯= 9.92 ±â€¯3.30 kPa, ADC = 0.770 ±â€¯0.023 µm2/s, all p < 0.001). Furthermore, we analyzed the springpot-powerlaw coefficient and observed a reduction in -20 °C specimens (0.22 ±â€¯0.14) compared to native tissue (0.40 ±â€¯0.10, p = 0.033) and -80 °C specimens (0.54 ±â€¯0.22, p = 0.002), that correlated with histological observations of sinusoidal dilation and collagen distortion within the space of Disse. Overall, the results suggest that shear modulus and water diffusion in liver tissue markedly decrease due to cell membrane degradation and cell detachment while viscosity-related properties appear to be more sensitive to distorted stromal and microvascular architecture.

Collagen networks determine viscoelastic properties of connective tissues yet do not hinder diffusion of the aqueous solvent.

Collagen accounts for the major extracellular matrix (ECM) component in many tissues and provides mechanical support for cells. Magnetic Resonance (MR) Imaging, MR based diffusion measurements and MR Elastography (MRE) are considered sensitive to the microstructure of tissues including collagen networks of the ECM. However, little is known whether water diffusion interacts with viscoelastic properties of tissues. This study combines highfield MR based diffusion measurements, novel compact tabletop MRE and confocal microscopy in collagen networks of different cross-linking states (untreated collagen gels versus additional treatment with glutaraldehyde). The consistency of bulk rheology and MRE within a wide dynamic range is demonstrated in heparin gels, a viscoelastic standard for MRE. Additional crosslinking of collagen led to an 8-fold increased storage modulus, a 4-fold increased loss modulus and a significantly decreased power law exponent, describing multi-relaxational behavior, corresponding to a pronounced transition from viscous-soft to elastic-rigid properties. Collagen network changes were not detectable by MR based diffusion measurements and microscopy which are sensitive to the micrometer scale. The MRE-measured shear modulus is sensitive to collagen fiber interactions which take place on the intrafiber level such as fiber stiffness. The insensitivity of MR based diffusion measurements to collagen hydrogels of different cross-linking states alludes that congeneric collagen structures in connective tissues do not hinder extracellular diffusive water transport. Furthermore, the glutaraldehyde induced rigorous changes in viscoelastic properties indicate that intrafibrillar dissipation is the dominant mode of viscous dissipation in collagen-dominated connective tissue.

In vivo time-harmonic ultrasound elastography of the human brain detects acute cerebral stiffness changes induced by intracranial pressure variations.

Cerebral stiffness (CS) reflects the biophysical environment in which neurons grow and function. While long-term CS changes can occur in the course of chronic neurological disorders and aging, little is known about acute variations of CS induced by intracranial pressure variations. Current gold standard methods for CS and intracranial pressure such as magnetic resonance elastography and direct pressure recordings are either expensive and slow or invasive. The study objective was to develop a real-time method for in vivo CS measurement and to demonstrate its sensitivity to physiological aging and intracranial pressure variations induced by the Valsalva maneuver in healthy volunteers. We used trans-temporal ultrasound time-harmonic elastography (THE) with external shear-wave stimulation by continuous and superimposed vibrations in the frequency range from 27 to 56 Hz. Multifrequency wave inversion generated maps of shear wave speed (SWS) as a surrogate maker of CS. On average, cerebral SWS was 1.56 ± 0.08 m/s with a tendency to reduce with age (R = -0.76, p < 0.0001) while Valsalva maneuver induced an immediate stiffening of the brain as reflected by a 10.8 ± 2.5% increase (p < 0.0001) in SWS. Our results suggest that CS is tightly linked to intracranial pressure and might be used in the future as non-invasive surrogate marker for intracranial pressure, which otherwise requires invasive measurements.