Brain Magnetic Resonance Elastography Research Articles

Indirect Shear Wave Excitation for Brain Magnetic Resonance Elastography with Minimal Cerebral Blood Flow Alteration.

This work demonstrates the initial prototype of the indirect excitation method, which can generate shear waves in the brain with minimal changes in CBF. A simple system was designed to produce stable vibrations underneath the neck. Instead of directly stimulating the skull, shear waves were indirectly transmitted to the brain through the spine and brainstem. Phantom results showed that the proposed actuator did not interfere with the routine imaging sequence and successfully generated multifrequency shear waves. When compared with the conventional direct head stimulation method, brain MRE results from the proposed actuator showed no significant differences in terms of intraclass correlation coefficients (ICC) and coefficients of variation (CV). Moreover, the octahedral shear strain (OSS) generated by the indirect excitation in the frontal and parietal lobes decreased by 25.96% and 16.73% respectively. Evaluation of CBF in healthy volunteers revealed no significant changes for the indirect excitation method, whereas significant decreases in CBF were observed in four subregions when employing direct excitation. The proposed actuator offers a more accurate and comfortable approach to MRE measurements while causing minimal CBF alterations. This work presents the first demonstration of an indirect excitation brain MRE system that minimizes CBF changes, thus holding potential for future applications of brain MRE.

Full-field, frequency-domain comparison of simulated and measured human brain deformation.

We propose a robust framework for quantitatively comparing model-predicted and experimentally measured strain fields in the human brain during harmonic skull motion. Traumatic brain injuries (TBIs) are typically caused by skull impact or acceleration, but how skull motion leads to brain deformation and consequent neural injury remains unclear and comparison of model predictions to experimental data remains limited. Magnetic resonance elastography (MRE) provides high-resolution, full-field measurements of dynamic brain deformation induced by harmonic skull motion. In the proposed framework, full-field strain measurements from human brain MRE in vivo are compared to simulated strain fields from models with similar harmonic loading. To enable comparison, the model geometry and subject anatomy, and subsequently, the predicted and measured strain fields are nonlinearly registered to the same standard brain atlas. Strain field correlations , both global (over the brain volume) and local (over smaller sub-volumes), are then computed from the inner product of the complex-valued strain tensors from model and experiment at each voxel. To demonstrate our approach, we compare strain fields from MRE in six human subjects to predictions from two previously developed models. Notably, global values are higher when comparing strain fields from different subjects ( ) than when comparing strain fields from either of the two models to strain fields in any subject. The proposed framework provides a quantitative method to assess similarity (and to identify discrepancies) between model predictions and experimental measurements of brain deformation, and thus can aid in the development and evaluation of improved models of brain biomechanics.

Association Between Anatomic and Clinical Indicators of Injury Severity After Moderate-Severe Traumatic Brain Injury: A Pilot Study Using Multiparametric Magnetic Resonance Imaging.

This study sought to identify whether an anatomical indicator of injury severity as measured by multiparametric magnetic resonance imaging (MRI) including magnetic resonance elastography (MRE), is predictive of a clinical measure of injury severity after moderate-severe traumatic brain injury (TBI). Nine individuals who were admitted to acute inpatient rehabilitation after moderate-to-severe TBI completed a comprehensive MRI protocol prior to discharge from rehabilitation, which included conventional MRI with diffusion tensor imaging (DTI). Of those, five of nine also underwent brain MRE to measure the brain parenchyma stiffness. Clinical severity of injury was measured by the length of post-traumatic amnesia (PTA). MRI-assessed non-hemorrhage contusion score and hemorrhage score, DTI-measured white matter fractional anisotropy, and MRE-measured lesion stiffness were all assessed. A higher hemorrhagic score was significantly associated with a longer length of PTA (p = 0.026). Participants with a longer PTA tended to have a higher non-hemorrhage contusion score and softer contusion lesions than the contralateral control side, although the small sample size did not allow for assessment of a significant association. To our knowledge, this is the first report applying MRI/MRE imaging protocol to quantitate altered brain anatomy after moderate-severe TBI and its association with PTA, a known clinical predictor of post-acute outcome. Future larger studies could lead to the development of prediction models that integrate clinical data with anatomical (MRI), structural (DTI), and mechanical (MRE) changes caused by TBI, to inform prognosis and care planning.

Effect of Direction and Frequency of Skull Motion on Mechanical Vulnerability of the Human Brain.

Strain energy and kinetic energy in the human brain were estimated by magnetic resonance elastography (MRE) during harmonic excitation of the head, and compared to characterize the effect of loading direction and frequency on brain deformation. In brain MRE, shear waves are induced by external vibration of the skull and imaged by a modified MR imaging sequence; the resulting harmonic displacement fields are typically "inverted" to estimate mechanical properties, like stiffness or damping. However, measurements of tissue motion from MRE also illuminate key features of the response of the brain to skull loading. In this study, harmonic excitation was applied in two different directions and at five different frequencies from 20 to 90 Hz. Lateral loading induced primarily left-right head motion and rotation in the axial plane; occipital loading induced anterior-posterior head motion and rotation in the sagittal plane. The ratio of strain energy to kinetic energy (SE/KE) depended strongly on both direction and frequency. The ratio of SE/KE was approximately four times larger for lateral excitation than for occipital excitation and was largest at the lowest excitation frequencies studied. These results are consistent with clinical observations that suggest lateral impacts are more likely to cause injury than occipital or frontal impacts, and also with observations that the brain has low-frequency (∼10 Hz) natural modes of oscillation. The SE/KE ratio from brain MRE is potentially a simple and powerful dimensionless metric of brain vulnerability to deformation and injury.

Viscoelastic Property of the Brain Assessed With Magnetic Resonance Elastography and Its Association With Glymphatic System in Neurologically Normal Individuals.

To investigate the feasibility of assessing the viscoelastic properties of the brain using magnetic resonance elastography (MRE) and a novel MRE transducer to determine the relationship between the viscoelastic properties and glymphatic function in neurologically normal individuals. This prospective study included 47 neurologically normal individuals aged 23-74 years (male-to-female ratio, 21:26). The MRE was acquired using a gravitational transducer based on a rotational eccentric mass as the driving system. The magnitude of the complex shear modulus |G*| and the phase angle ϕ were measured in the centrum semiovale area. To evaluate glymphatic function, the Diffusion Tensor Image Analysis Along the Perivascular Space (DTI-ALPS) method was utilized and the ALPS index was calculated. Univariable and multivariable (variables with P < 0.2 from the univariable analysis) linear regression analyses were performed for |G*| and ϕ and included sex, age, normalized white matter hyperintensity (WMH) volume, brain parenchymal volume, and ALPS index as covariates. In the univariable analysis for |G*|, age (P = 0.005), brain parenchymal volume (P = 0.152), normalized WMH volume (P = 0.011), and ALPS index (P = 0.005) were identified as candidates with P < 0.2. In the multivariable analysis, only the ALPS index was independently associated with |G*|, showing a positive relationship (β = 0.300, P = 0.029). For ϕ, normalized WMH volume (P = 0.128) and ALPS index (P = 0.015) were identified as candidates for multivariable analysis, and only the ALPS index was independently associated with ϕ (β = 0.057, P = 0.039). Brain MRE using a gravitational transducer is feasible in neurologically normal individuals over a wide age range. The significant correlation between the viscoelastic properties of the brain and glymphatic function suggests that a more organized or preserved microenvironment of the brain parenchyma is associated with a more unimpeded glymphatic fluid flow.

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

Purpose: 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 (BMI). Methods: Multifrequency MRE with 2D- and 3D-tomoelastography postprocessing was applied to the brains of 31 volunteers (age range: 22-61years) for analyzing the coefficient of variation (CV) and effects of biological factors. Eleven volunteers were rescanned after 1day and 1year to determine intraclass correlation coefficient (ICC) and identify possible long-term changes. Results: 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 1day 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 1year, 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 BMI. Conclusion: 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.

Cover Image

The cover image is based on the Research Article An electromagnetic actuator for brain magnetic resonance elastography with high frequency accuracy by Suhao Qiu et al., https://doi.org/10.1002/nbm.4592.

Correlated noise in brain magnetic resonance elastography.

Magnetic resonance elastography (MRE) uses phase-contrast MRI to generate mechanical property maps of the in vivo brain through imaging of tissue deformation from induced mechanical vibration. The mechanical property estimation process in MRE can be susceptible to noise from physiological and mechanical sources encoded in the phase, which is expected to be highly correlated. This correlated noise has yet to be characterized in brain MRE, and its effects on mechanical property estimates computed using inversion algorithms are undetermined. To characterize the effects of signal noise in MRE, we conducted 3 experiments quantifying (1) physiomechanical sources of signal noise, (2) physiological noise because of cardiac-induced movement, and (3) impact of correlated noise on mechanical property estimates. We use a correlation length metric to estimate the extent that correlated signal persists in MRE images and demonstrate the effect of correlated noise on property estimates through simulations. We found that both physiological noise and vibration noise were greater than image noise and were spatially correlated across all subjects. Added physiological and vibration noise to simulated data resulted in property maps with higher error than equivalent levels of Gaussian noise. Our work provides the foundation to understand contributors to brain MRE data quality and provides recommendations for future work to correct for signal noise in MRE.

An electromagnetic actuator for brain magnetic resonance elastography with high frequency accuracy.

Our goal is to design, test and verify an electromagnetic actuator for brain magnetic resonance elastography (MRE). We proposed a grappler-shaped design that can transmit stable vibrations into the brain. To validate its performance, simulations were carried out to ensure the electromagnetic field generated by the actuator did not interfere with the B0 field. The actuation vibration spectrum was analyzed to verify the actuation accuracy. Phantom and volunteer experiments were carried out to evaluate the performance of the actuator. Simulation of the magnetic field showed that the proposed actuator has a fringe field of less than 3G in the imaging region. The phantom experiments showed that the proposed actuator did not interfere with the routine imaging sequences. The measured vibration spectra demonstrated that the frequency offset was about one third that of a pneumatic device and the transmission efficiency was three times higher. The shear moduli estimated from brain MRE were consistent with those from the literature. The actuation frequency of the proposed actuator has less frequency offset and off-center frequency components compared with the pneumatic counterpart. The whole actuator weighted only 980 g. The actuator can carry out multifrequency MRE on the brain with high accuracy. It is easy to use, comfortable for the patient and portable.

Harnessing brain waves: a review of brain magnetic resonance elastography for clinicians and scientists entering the field.

Brain magnetic resonance elastography (MRE) is an imaging technique capable of accurately and non-invasively measuring the mechanical properties of the living human brain. Recent studies have shown that MRE has potential to provide clinically useful information in patients with intracranial tumors, demyelinating disease, neurodegenerative disease, elevated intracranial pressure, and altered functional states. The objectives of this review are: (1) to give a general overview of the types of measurements that have been obtained with brain MRE in patient populations, (2) to survey the tools currently being used to make these measurements possible, and (3) to highlight brain MRE-based quantitative biomarkers that have the highest potential of being adopted into clinical use within the next 5 to 10 years. The specifics of MRE methodology strategies are described, from wave generation to material parameter estimations. The potential clinical role of MRE for characterizing and planning surgical resection of intracranial tumors and assessing diffuse changes in brain stiffness resulting from diffuse neurological diseases and altered intracranial pressure are described. In addition, the emerging technique of functional MRE, the role of artificial intelligence in MRE, and promising applications of MRE in general neuroscience research are presented.

Sensitivity analysis of effective transverse shear viscoelastic and diffusional properties of myelinated white matter

Motivated by the need to interpret the results from a combined use of in vivo brain Magnetic Resonance Elastography (MRE) and Diffusion Tensor Imaging (DTI), we developed a computational framework to study the sensitivity of single-frequency MRE and DTI metrics to white matter microstructure and cell-level mechanical and diffusional properties. White matter was modeled as a triphasic unidirectional composite, consisting of parallel cylindrical inclusions (axons) surrounded by sheaths (myelin), and embedded in a matrix (glial cells plus extracellular matrix). Only 2D mechanics and diffusion in the transverse plane (perpendicular to the axon direction) was considered, and homogenized (effective) properties were derived for a periodic domain containing a single axon. The numerical solutions of the MRE problem were performed with ABAQUS and by employing a sophisticated boundary-conforming grid generation scheme. Based on the linear viscoelastic response to harmonic shear excitation and steady-state diffusion in the transverse plane, a systematic sensitivity analysis of MRE metrics (effective transverse shear storage and loss moduli) and DTI metric (effective radial diffusivity) was performed for a wide range of microstructural and intrinsic (phase-based) physical properties. The microstructural properties considered were fiber volume fraction, and the myelin sheath/axon diameter ratio. The MRE and DTI metrics are very sensitive to the fiber volume fraction, and the intrinsic viscoelastic moduli of the glial phase. The MRE metrics are nonlinear functions of the fiber volume fraction, but the effective diffusion coefficient varies linearly with it. Finally, the transverse metrics of both MRE and DTI are insensitive to the axon diameter in steady state. Our results are consistent with the limited anisotropic MRE and co-registered DTI measurements, mainly in the corpus callosum, available in the literature. We conclude that isotropic MRE and DTI constitutive models are good approximations for myelinated white matter in the transverse plane. The unidirectional composite model presented here is used for the first time to model harmonic shear stress under MRE-relevant frequency on the cell level. This model can be extended to 3D in order to inform the solution of the inverse problem in MRE, establish the biological basis of MRE metrics, and integrate MRE/DTI with other modalities towards increasing the specificity of neuroimaging.

Robustness of MR Elastography in the Healthy Brain: Repeatability, Reliability, and Effect of Different Reconstruction Methods.

Changes in brain stiffness can be an important biomarker for neurological disease. Magnetic resonance elastography (MRE) quantifies tissue stiffness, but the results vary between acquisition and reconstruction methods. To measure MRE repeatability and estimate the effect of different reconstruction methods and varying data quality on estimated brain stiffness. Prospective. Fifteen healthy subjects. 3T MRI, gradient-echo elastography sequence with a 50 Hz vibration frequency. Imaging was performed twice in each subject. Images were reconstructed using a curl-based and a finite-element-model (FEM)-based method. Stiffness was measured in the whole brain, in white matter, and in four cortical and four deep gray matter regions. Repeatability coefficients (RC), intraclass correlation coefficients (ICC), and coefficients of variation (CV) were calculated. MRE data quality was quantified by the ratio between shear waves and compressional waves. Median values with range are presented. Reconstruction methods were compared using paired Wilcoxon signed-rank tests, and Spearman's rank correlation was calculated between MRE data quality and stiffness. Holm-Bonferroni corrections were employed to adjust for multiple comparisons. In the whole brain, CV was 4.3% and 3.8% for the curl and the FEM reconstruction, respectively, with 4.0-12.8% for subregions. Whole-brain ICC was 0.60-0.74, ranging from 0.20 to 0.89 in different regions. RC for the whole brain was 0.14 kPa and 0.17 kPa for the curl and FEM methods, respectively. FEM reconstruction resulted in 39% higher stiffness than the curl reconstruction (P < 0.05). MRE data quality, defined as shear-compression wave ratio, was higher in peripheral regions than in central regions of the brain (P < 0.05). No significant correlations were observed between MRE data quality and stiffness estimates. MRE of the human brain is a robust technique in terms of repeatability. Caution is warranted when comparing stiffness values obtained with different techniques. 1 TECHNICAL EFFICACY STAGE: 1.

Standard-space atlas of the viscoelastic properties of the human brain.

Standard anatomical atlases are common in neuroimaging because they facilitate data analyses and comparisons across subjects and studies. The purpose of this study was to develop a standardized human brain atlas based on the physical mechanical properties (i.e., tissue viscoelasticity) of brain tissue using magnetic resonance elastography (MRE). MRE is a phase contrast‐based MRI method that quantifies tissue viscoelasticity noninvasively and in vivo thus providing a macroscopic representation of the microstructural constituents of soft biological tissue. The development of standardized brain MRE atlases are therefore beneficial for comparing neural tissue integrity across populations. Data from a large number of healthy, young adults from multiple studies collected using common MRE acquisition and analysis protocols were assembled (N = 134; 78F/ 56 M; 18–35 years). Nonlinear image registration methods were applied to normalize viscoelastic property maps (shear stiffness, μ, and damping ratio, ξ) to the MNI152 standard structural template within the spatial coordinates of the ICBM‐152. We find that average MRE brain templates contain emerging and symmetrized anatomical detail. Leveraging the substantial amount of data assembled, we illustrate that subcortical gray matter structures, white matter tracts, and regions of the cerebral cortex exhibit differing mechanical characteristics. Moreover, we report sex differences in viscoelasticity for specific neuroanatomical structures, which has implications for understanding patterns of individual differences in health and disease. These atlases provide reference values for clinical investigations as well as novel biophysical signatures of neuroanatomy. The templates are made openly available (github.com/mechneurolab/mre134) to foster collaboration across research institutions and to support robust cross‐center comparisons.

Paediatric brain tissue properties measured with magnetic resonance elastography.

The aim of this study is to characterise the stiffness of white and grey matter in paediatric subjects using magnetic resonance elastography (MRE) and to determine whether these properties change throughout normal development. MRE was performed using a clinical 3T MRI scanner at three frequencies (30, 40 and 60Hz) on 36 healthy paediatric subjects aged between 7 and 18years (19 F) and 11 adults aged 23-44years (6 F). Anatomical and diffusion tensor imaging was also collected. The stiffness quantified as the magnitude of the complex shear modulus (G*), fractional anisotropy (FA), mean diffusivity (MD) and volume of white and grey matter were calculated. One-way analysis of variance and Tukey's multiple comparison tests were used to compare data in age groups separated into children (7-12years), adolescents (13-18years) and adults (18+ years), and Spearman's correlations were performed for paediatric data. White and grey matter stiffness for each frequency and their frequency dependence was found to be very similar in paediatric and adult subjects (p > 0.05 all variables). No significant correlations were found when comparing G* with age, FA, MD or volume. Adult G*, FA, MD and volume values were within range of others reported in the literature. Paediatric white and grey matter stiffness values are similar to those of adults. We conclude that clinically, adult values can be used as a baseline measure in paediatric brain MRE.

Reliable preparation of agarose phantoms for use in quantitative magnetic resonance elastography

Agarose phantoms are one type of phantom commonly used in developing in vivo brain magnetic resonance elastography (MRE) sequences because they are inexpensive and easy to work with, store, and dispose of; however, protocols for creating agarose phantoms are non-standardized and often result in inconsistent phantoms with significant variability in mechanical properties. Many magnetic resonance imaging (MRI) and ultrasound studies use phantoms, but often these phantoms are not tailored for desired mechanical properties and as such are too stiff or not mechanically consistent enough to be used in MRE. In this work, we conducted a systematic study of agarose phantom creation parameters to identify those factors that are most conducive to producing mechanically consistent agarose phantoms for MRE research. We found that cooling rate and liquid temperature affected phantom homogeneity. Phantom stiffness is affected by agar concentration (quadratically), by final liquid temperature and salt content in phantoms, and by the interaction of these two metrics each with stir rate. We captured and quantified the implied relationships with a regression model that can be used to estimate stiffness of resulting phantoms. Additionally, we characterized repeatability, stability over time, impact on MR signal parameters, and differences in agar gel microstructure. This protocol and regression model should prove beneficial in future MRE development studies that use phantoms to determine stiffness measurement accuracy.

Shear viscoelasticity in the human brain by MR-elastography at low and high frequencies

Brain Magnetic Resonance Elastography is usually performed in humans with 3 mm isotropic voxels at low mechanical excitation frequency, between 40 Hz and 60 Hz, to ease wave propagation deep inside...

Magnetic resonance elastography of brain: Comparison between anisotropic and isotropic stiffness and its correlation to age.

Noninvasive measurement of mechanical properties of brain tissue using magnetic resonance elastography (MRE) has been a promising method for investigating neurologic disorders such as multiple sclerosis, hydrocephalus, and Alzheimer's. However, because of the regional and directional dependency of brain stiffness, estimating anisotropic stiffness is important. This study investigates isotropic and anisotropic stiffness as a function of age as well as the correlation between isotropic and anisotropic stiffness. MRE and diffusion tensor imaging (DTI) were performed on 28 healthy subjects with age ranges between 18-62 y. Isotropic and anisotropic stiffness was measured and compared with age for different regions of interest such as the thalamus, corpus callosum, gray matter, white matter, and whole brain. Isotropic stiffness in gray matter (rs = -0.57; P = 0.001) showed a significant decrease with age. Anisotropic stiffness in gray matter showed a significant decrease with age in C11 through C66 and in the thalamus, only in C33 . Between anisotropic and isotropic stiffness, gray matter showed a significant positive correlation in C11 through C66 , C22 and C66 showed a significant negative correlation in the thalamus and whole brain, and C44 showed a negative correlation in the corpus callosum. No significant difference between genders was observed in any measurements. This study demonstrated a change in isotropic and anisotropic stiffness with age in different regions of the brain along with a correlation of anisotropic stiffness to isotropic stiffness.

Magnetic resonance elastography of the brain: A study of feasibility and reproducibility using an ergonomic pillow-like passive driver

Magnetic resonance elastography (MRE) can be used to noninvasively resolve the displacement pattern of induced mechanical waves propagating in tissue. The goal of this study is to establish an ergonomically flexible passive-driver design for brain MRE, to evaluate the reproducibility of MRE tissue-stiffness measurements, and to investigate the relationship between tissue-stiffness measurements and driver frequencies.An ergonomically flexible passive pillow-like driver was designed to induce mechanical waves in the brain. Two-dimensional finite-element simulation was used to evaluate mechanical wave propagation patterns in brain tissues. MRE scans were performed on 10 healthy volunteers at mechanical frequencies of 60, 50, and 40 Hz. An axial mid-brain slice was acquired using an echo-planar imaging sequence to map the displacement pattern with the motion-encoding gradient along the through-plane (z) direction. All subjects were scanned and rescanned within 1 h. The Wilcoxon signed-rank test was used to test for differences between white matter and gray matter shear-stiffness values. One-way analysis of variance (ANOVA) was used to test for differences between shear-stiffness measurements made at different frequencies. Scan-rescan reproducibility was evaluated by calculating the within-subject coefficient of variation (CV) for each subject.The finite-element simulation showed that a pillow-like passive driver is capable of efficient shear-wave propagation through brain tissue. No subjects complained about discomfort during MRE acquisitions using the ergonomically designed driver. The white-matter elastic modulus (mean ± standard deviation) across all subjects was 3.85 ± 0.12 kPa, 3.78 ± 0.15 kPa, and 3.36 ± 0.11 kPa at frequencies of 60, 50, and 40 Hz, respectively. The gray-matter elastic modulus across all subjects was 3.33 ± 0.14 kPa, 2.82 ± 0.16 kPa, and 2.24 ± 0.14 kPa at frequencies of 60, 50, and 40 Hz, respectively. The Wilcoxon signed-rank test confirmed that the shear stiffness was significantly higher in white matter than gray matter at all three frequencies. The ranges of within-subject coefficients of variation for white matter, gray matter, and whole-brain shear-stiffness measurements for the three frequencies were 1.8–3.5% (60 Hz), 4.7–6.0% (50 Hz), and 3.7–4.1% (40 Hz).An ergonomic pneumatic pillow-like driver is feasible for highly reproducible in vivo evaluation of brain-tissue shear stiffness. Brain-tissue shear-stiffness values were frequency-dependent, thus emphasizing the importance of standardizing MRE acquisition protocols in multi-center studies.

Magnetic resonance elastography of the human brain using a multiphase DENSE acquisition.

PurposeIn magnetic resonance elastography (MRE), a series of time‐shifted images is acquired at specific phase offsets in relation to an induced mechanical excitation. To efficiently gather the set of phase offset images and to overcome limitations due to prolonged TEs and related susceptibility artifacts at low‐frequency MRE, we developed an improved displacement encoding with a stimulated echoes (DENSE) method.MethodsThe proposed multiphase DENSE‐MRE acquisition scheme allows full sampling of the wave propagation in 1 encoding direction during each TR using multiple readouts at specific phase offsets. With this approach, all phase offsets can be imaged in 1 TR without the need for whole sequence repetitions at time‐shifted offsets relative to the excitation motion. We tested this technique in phantom experiments with 60 Hz and in the brain of 4 volunteers using 20‐Hz harmonic excitation.ResultsThree‐dimensional wave propagation could be acquired in 7 minutes 30 seconds. Following background phase elimination, clear wave images were obtained, showing the propagation of the waves over time. Calculated shear modulus maps of the phantom matched well to the maps obtained by conventional gradient‐echo MRE. In the brain, low‐frequency DENSE‐MRE images were free of susceptibility‐induced artifacts and the calculated maps showed a median global complex shear modulus magnitude of 0.72 kPa and phase angle of 1.03 rad across volunteers.ConclusionThe proposed multiphase DENSE approach allows efficient low‐frequency MRE with short TEs and is well‐suited for low‐frequency MRE of the human brain.

Fast Robust Dejitter and Interslice Discontinuity Removal in MRI Phase Acquisitions: Application to Magnetic Resonance Elastography.

MRI phase contrast imaging methods that assemble slice-wise acquisitions into volumes can contain interslice phase discontinuities (IPDs) over the course of the scan from sources, including unavoidable physiological activity. In magnetic resonance elastography (MRE), this can alter wavelength and tissue stiffness estimates, invalidating the analysis. We first model this behavior as jitter along the z-axis of the phase of 3D complex-valued wave volumes. A two-step image processing pipeline is then proposed that removes IPDs. First, constant slicewise phase shift is removed with a novel, non-convex dejittering algorithm. Then, regional physiological noise artifacts are removed with novel filtering of 3D wavelet coefficients. Calibration of two pipeline coefficients, the dejitter parameter α and the wavelet band high-pass coefficient ωc , was first performed on a finite-element method brain phantom. A comparative investigation was then performed, on a cohort of 48 brain acquisitions, of four approaches to IPDs: 1) the proposed method; 2) a "control" condition of neglect of IPDs; 3) an anisotropic wavelet-based method; and 4) a method of in-plane (2D) processing. The present method showed medians of [Formula: see text] Pa for a multifrequency wave inversion centered at 40 Hz which was within 6% of methods 3) and 4), while neglect produced [Formula: see text] estimates a mean of 17% lower. The proposed method reduced the value range of the cohort against methods 3) and 4) by 29% and 31%, respectively. Such reduction in variance enhances the ability of brain MRE to predict subtler physiological changes. Our theoretical approach further enables more powerful applications of fundamental findings in noise and denoising to MRE.