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

Vascular determinants of hippocampal viscoelastic properties in healthy adults across the lifespan

Arterial stiffness and cerebrovascular pulsatility are non-traditional risk factors of Alzheimer's disease. However, there is a gap in understanding the earliest mechanisms that link these vascular determinants to brain aging. Changes to mechanical tissue properties of the hippocampus (HC), a brain structure essential for memory encoding, may reflect the impact of vascular dysfunction on brain aging. We tested the hypothesis that arterial stiffness and cerebrovascular pulsatility are related to HC tissue properties in healthy adults across the lifespan. Twenty-five adults underwent measurements of brachial blood pressure (BP), large elastic artery stiffness, middle cerebral artery pulsatility index (MCAv PI), and magnetic resonance elastography (MRE), a sensitive measure of HC viscoelasticity. Individuals with higher carotid pulse pressure (PP) exhibited lower HC stiffness (β = -0.39, r = -0.41, p = 0.05), independent of age and sex. Collectively, carotid PP and MCAv PI significantly explained a large portion of the total variance in HC stiffness (adjusted R2 = 0.41, p = 0.005) in the absence of associations with HC volumes. These cross-sectional findings suggest that the earliest reductions in HC tissue properties are associated with alterations in vascular function.

Stiffness pulsation of the human brain detected by non-invasive time-harmonic elastography

Introduction: Cerebral pulsation is a vital aspect of cerebral hemodynamics. Changes in arterial pressure in response to cardiac pulsation cause cerebral pulsation, which is related to cerebrovascular compliance and cerebral blood perfusion. Cerebrovascular compliance and blood perfusion influence the mechanical properties of the brain, causing pulsation-induced changes in cerebral stiffness. However, there is currently no imaging technique available that can directly quantify the pulsation of brain stiffness in real time. Methods: Therefore, we developed non-invasive ultrasound time-harmonic elastography (THE) technique for the real-time detection of brain stiffness pulsation. We used state-of-the-art plane-wave imaging for interleaved acquisitions of shear waves at a frequency of 60 Hz to measure stiffness and color flow imaging to measure cerebral blood flow within the middle cerebral artery. In the second experiment, we used cost-effective lineby-line B-mode imaging to measure the same mechanical parameters without flow imaging to facilitate future translation to the clinic. Results: In 10 healthy volunteers, stiffness increased during the passage of the arterial pulse wave from 4.8% ± 1.8% in the temporal parenchyma to 11% ± 5% in the basal cisterns and 13% ± 9% in the brain stem. Brain stiffness peaked in synchrony with cerebral blood flow at approximately 180 ± 30 ms after the cardiac R-wave. Line-by-line THE provided the same stiffness values with similar time resolution as high-end plane-wave THE, demonstrating the robustness of brain stiffness pulsation as an imaging marker. Discussion: Overall, this study sets the background and provides reference values for time-resolved THE in the human brain as a cost-efficient and easy-touse mechanical biomarker associated with cerebrovascular compliance.

Measuring brain beats: Cardiac-aligned fast functional magnetic resonance imaging signals

Blood and cerebrospinal fluid (CSF) pulse and flow throughout the brain, driven by the cardiac cycle. These fluid dynamics, which are essential to healthy brain function, are characterized by several noninvasive magnetic resonance imaging (MRI) methods. Recent developments in fast MRI, specifically simultaneous multislice acquisition methods, provide a new opportunity to rapidly and broadly assess cardiac-driven flow, including CSF spaces, surface vessels and parenchymal vessels. We use these techniques to assess blood and CSF flow dynamics in brief (3.5 min) scans on a conventional 3 T MRI scanner in five subjects. Cardiac pulses are measured with a photoplethysmography (PPG) on the index finger, along with functional MRI (fMRI) signals in the brain. We, retrospectively, align the fMRI signals to the heartbeat. Highly reliable cardiac-gated fMRI temporal signals are observed in CSF and blood on the timescale of one heartbeat (test-retest reliability within subjects R²  > 50%). In blood vessels, a local minimum is observed following systole. In CSF spaces, the ventricles and subarachnoid spaces have a local maximum following systole instead. Slower resting-state scans with slice timing, retrospectively, aligned to the cardiac pulse, reveal similar cardiac-gated responses. The cardiac-gated measurements estimate the amplitude and phase of fMRI pulsations in the CSF relative to those in the arteries, an estimate of the local intracranial impedance. Cardiac aligned fMRI signals can provide new insights about fluid dynamics or diagnostics for diseases where these dynamics are important.

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.

Noninvasive Detection of Intracranial Hypertension by Novel Ultrasound Time-Harmonic Elastography

OBJECTIVE

A method for measuring intracranial pressure (ICP) noninvasively has long been sought after in neurology and neurosurgery. Treatment failure in individuals presenting with unspecific symptoms such as headache, gait disturbance, or visual impairment occurring in response to increased ICP can lead to irreversible brain injury, progressive disability, and death. Guidelines for diagnostic ICP measurement recommend intracranial placement of pressure tip catheters or lumbar puncture (LP) despite their invasiveness and possible complications. As ICP fluctuations are closely associated with changes in brain stiffness, ultrasound elastography could be a valid method to detect ICP noninvasively and with short examination times.

MATERIALS AND METHODS

In this pilot study, we have investigated the use of time-harmonic shear waves, introduced into the brain by an external shaker, and measured in real-time by transtemporal ultrasound, for deducing a noninvasive imaging marker sensitive to elevated ICP. To this end, we developed cerebral ultrasound time-harmonic elastography for the noninvasive quantification of shear wave speed (SWS) as a surrogate marker of cerebral stiffness in a short examination time of a few minutes.

RESULTS

We found that SWS in patients enrolled for LP with confirmed intracranial hypertension was 1.81 ± 0.10 m/s, distinguishing them from healthy volunteers with excellent diagnostic accuracy (1.55 ± 0.08 m/s; P < 0.001; area under the curve, 0.99). Interestingly, values in symptomatic patients decreased to normal stiffness immediately after LP (1.56 ± 0.06 m/s, P < 0.001). Moreover, invasively measured opening pressure correlated with SWS measured before LP and liquid volume drained through the spinal tap with the SWS difference between the 2 measurements.

CONCLUSIONS

Collectively, our results suggest a tight link between cerebral stiffness and ICP and demonstrate that intracranial hypertension can be detected noninvasively within short examination times, opening avenues for diagnostic applications of cerebral ultrasound time-harmonic elastography in neurology and emergency medicine.

Fast magnetic resonance elastography with multiphase radial encoding and harmonic motion sparsity based reconstruction

Objective. To achieve fast magnetic resonance elastography (MRE) at a low frequency for better shear modulus estimation of the brain. Approach. We proposed a multiphase radial DENSE MRE (MRD-MRE) sequence and an improved GRASP algorithm utilizing the sparsity of the harmonic motion (SH-GRASP) for fast MRE at 20 Hz. For the MRD-MRE sequence, the initial position encoded by spatial modulation of magnetization (SPAMM) was decoded by an arbitrary number of readout blocks without increasing the number of phase offsets. Based on the harmonic motion, a modified total variation and temporal Fourier transform were introduced to utilize the sparsity in the temporal domain. Both phantom and brain experiments were carried out and compared with that from multiphase Cartesian DENSE-MRE (MCD-MRE), and conventional gradient echo sequence (GRE-MRE). Reconstruction performance was also compared with GRASP and compressed sensing. Main results. Results showed the scanning time of a fully sampled image with four phase offsets for MRD-MRE was only 1/5 of that from GRE-MRE. The wave patterns and estimated stiffness maps were similar to those from MCD-MRE and GRE-MRE. With SH-GRASP, the total scan time could be shortened by additional 4 folds, achieving a total acceleration factor of 20. Better metric values were also obtained using SH-GRASP for reconstruction compared with other algorithms. Significance. The MRD-MRE sequence and SH-GRASP algorithm can be used either in combination or independently to accelerate MRE, showing the potentials for imaging the brain as well as other organs.

In vivo stiffness of multiple sclerosis lesions is similar to that of normal-appearing white matter

In 1868, French neurologist Jean-Martin Charcot coined the term multiple sclerosis (MS) after his observation that numerous white matter (WM) glial scars felt like sclerotic tissue. Nowadays, magnetic resonance elastography (MRE) can generate images with contrast of stiffness (CS) in soft in vivo tissues and may therefore be sensitive to MS lesions, provided that sclerosis is indeed a mechanical signature of this disease. We analyzed CS in a total of 147 lesions in patients with relapsing-remitting MS, compared with control regions in contralateral brain regions, and phantom data as well as performed numerical simulations to determine the delineation limits of multifrequency MRE (20 - 40 Hz) in MS. MRE analysis of simulated waves revealed a delineation limit of approximately 10% CS for detecting 9-mm lesions (mean size in our patient population). Due to inversion bias, this limit is reached when true CS is -11% for soft and 35% for stiff lesions. In vivo MRE identified 35 stiffer lesions and 17 softer lesions compared with surrounding WM (mean stiffness: 934±82 Pa). However, a similar pattern was found in the contralateral brain, suggesting that the range of stiffness changes in WM lesions due to MS is within the normal range of WM variability and normal heterogeneity-related CS. Consequently, Charcot's original intuition that MS is a focal sclerotic disease can neither be dismissed nor confirmed by in vivo MRE. However, the observation that MS lesions do not markedly differ in stiffness from surrounding brain tissue suggests that marked tissue sclerosis is not a mechanical signature of MS. STATEMENT OF SIGNIFICANCE: Multiple sclerosis (MS) was named by J.M. Charcot after the sclerotic changes in brain tissue he found in post-mortem autopsies. Since then, nothing has been revealed about the actual stiffening of MS lesions in vivo. Studying the viscoelastic properties of plaques in their natural environment is a major challenge that can only be overcome by MR elastography (MRE). Therefore, we used multifrequency MRE to answer the question whether MS lesions in patients with a relapsing-remitting disease course are mechanically different than surrounding tissue. Our findings suggest that the range of stiffness changes in white matter lesions due to MS is within the normal range of white matter variability and in vivo tissue sclerosis might not be a mechanical signature of MS.

Preliminary Findings on the Potential Use of Magnetic Resonance Elastography to Diagnose Lacunar Infarction

PURPOSE

Lacunar infarction is usually diagnosed by conventional technologies, such as CT and diffusion weighted imaging (DWI). To improve the accuracy of diagnosis, neurocognitive screening is still needed. Therefore, additional imaging methods that can assist and provide more accurate and rapid diagnostics are urgently needed. As an initial step towards potentially using MR elastography (MRE) for such diagnostic purposes, we tested the hypothesis that the mechanical properties of tissue in the vicinity of cerebral vasculature change following lacunar infarction in a way that can be quantified using MRE.

PATIENTS AND METHODS

MRE and MR angiography (MRA) images from 51 patients diagnosed with lacunar infarction and 54 healthy volunteers were acquired on a 3T scanner. All diagnoses were confirmed by matching neurocognitive test results to locations of flow obstruction in MRA. ROIs of the cerebral vessels segmented on the MRA images were mapped to the MRE images. Interpolation-based inversion was applied to estimate the regional biomechanical properties of ROIs that included cerebral vessels. The effects of lacunar infarction, sex, and age were analyzed using analysis of covariance (ANOCOVA).

RESULTS

Shear moduli over vessel ROIs were significantly lower for the lacunar infarction group than those of the healthy control group. A positive correlation between modulus over vessel ROIs and age was observed. However, no significant correlation was found between sex and the regional biomechanical properties of the vessel ROIs.

CONCLUSION

Results supported the hypothesis and suggest that biomechanical properties may be of utility in diagnosis of lacunar infarction.

Multifrequency magnetic resonance elastography for elasticity quantitation and optimal tissue discrimination: A two-platform liver fibrosis mimicking phantom study

In the framework of algebraic inversion, magnetic resonance elastography (MRE) repeatability, reproducibility and robustness were evaluated on extracted shear velocities (or elastic moduli). The same excitation system was implemented at two sites equipped with clinical MR scanners of 1.5 and 3 T. A set of four elastic, isotropic, homogeneous calibrated phantoms of distinct elasticity representing the spectrum of liver fibrosis severity was mechanically characterized. The repeatability of the measurements and the reproducibility between the two platforms were found to be excellent with mean coefficients of variations of 1.62% for the shear velocity mean values and 1.95% for the associated standard deviations. MRE velocities were robust to the amplitude and pattern variations of the displacement field with virtually no difference between outcomes from both magnets at identical excitation frequencies, even when the displacement field amplitude was six times smaller. However, MRE outcomes were very sensitive to the number of voxels per wavelength, s, of the recorded displacement field, with relative biases reaching 62% and precision loss by a factor of up to 23.5. For both magnetic field strengths, MRE accuracy and precision were largely degraded outside of established conditions of validity (6 ≲ s ≲ 9), resulting in estimated shear velocity values not significantly different between phantoms of increasing elasticity. When fulfilling the spatial sampling conditions, either prospectively in the acquisition or retrospectively before the reconstruction, MRE produced quantitative measurements that allowed to unambiguously discriminate, with infinitesimal p values, between the phantoms mimicking increasing severity of liver fibrosis.

Inversion-recovery MR elastography of the human brain for improved stiffness quantification near fluid-solid boundaries

PURPOSE

In vivo MR elastography (MRE) holds promise as a neuroimaging marker. In cerebral MRE, shear waves are introduced into the brain, which also stimulate vibrations in adjacent CSF, resulting in blurring and biased stiffness values near brain surfaces. We here propose inversion-recovery MRE (IR-MRE) to suppress CSF signal and improve stiffness quantification in brain surface areas.

METHODS

Inversion-recovery MRE was demonstrated in agar-based phantoms with solid-fluid interfaces and 11 healthy volunteers using 31.25-Hz harmonic vibrations. It was performed by standard single-shot, spin-echo EPI MRE following 2800-ms IR preparation. Wave fields were acquired in 10 axial slices and analyzed for shear wave speed (SWS) as a surrogate marker of tissue stiffness by wavenumber-based multicomponent inversion.

RESULTS

Phantom SWS values near fluid interfaces were 7.5 ± 3.0% higher in IR-MRE than MRE (P = .01). In the brain, IR-MRE SNR was 17% lower than in MRE, without influencing parenchymal SWS (MRE: 1.38 ± 0.02 m/s; IR-MRE: 1.39 ± 0.03 m/s; P = .18). The IR-MRE tissue-CSF interfaces appeared sharper, showing 10% higher SWS near brain surfaces (MRE: 1.01 ± 0.03 m/s; IR-MRE: 1.11 ± 0.01 m/s; P < .001) and 39% smaller ventricle sizes than MRE (P < .001).

CONCLUSIONS

Our results show that brain MRE is affected by fluid oscillations that can be suppressed by IR-MRE, which improves the depiction of anatomy in stiffness maps and the quantification of stiffness values in brain surface areas. Moreover, we measured similar stiffness values in brain parenchyma with and without fluid suppression, which indicates that shear wavelengths in solid and fluid compartments are identical, consistent with the theory of biphasic poroelastic media.

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.

TURBINE-MRE: A 3D hybrid radial-Cartesian EPI acquisition for MR elastography

PURPOSE

To develop a novel magnetic resonance elastography (MRE) acquisition using a hybrid radial EPI readout scheme (TURBINE), and to demonstrate its feasibility to obtain wave images and stiffness maps in a phantom and in vivo brain.

METHOD

The proposed 3D TURBINE-MRE is based on a spoiled gradient-echo MRE sequence with the EPI readout radially rotating about the phase-encoding axis to sample a full 3D k-space. A polyvinyl chloride phantom and 6 volunteers were scanned on a compact 3T GE scanner with a 32-channel head coil at 80 Hz and 60 Hz external vibration, respectively. For comparison, a standard 2D, multislice, spin-echo (SE) EPI-MRE acquisition was also performed with the same motion encoding and resolution. The TURBINE-MRE images were off-line reconstructed with iterative SENSE algorithm. The regional ROI analysis was performed on the 6 volunteers, and the median stiffness values were compared between SE-EPI-MRE and TURBINE-MRE.

RESULTS

The 3D wave-field images and the generated stiffness maps were comparable between TURBINE-MRE and standard SE-EPI-MRE for the phantom and the volunteers. The Bland-Altman plot showed no significant difference in the median regional stiffness values between the two methods. The stiffness measured with the 2 methods had a strong linear relationship with a Pearson correlation coefficient of 0.943.

CONCLUSION

We demonstrated the feasibility of the new TURBINE-MRE sequence for acquiring the desired 3D wave-field data and stiffness maps in a phantom and in-vivo brains. This pilot study encourages further exploration of TURBINE-MRE for functional MRE, free-breathing abdominal MRE, and cardiac MRE applications.