Imaging localized neuronal activity at fast time scales through biomechanics

Making sense of scattering: Seeing microstructure through shear waves

The physics of shear waves traveling through matter carries fundamental insights into its structure, for instance, quantifying stiffness for disease characterization. However, the origin of shear wave attenuation in tissue is currently not properly understood. Attenuation is caused by two phenomena: absorption due to energy dissipation and scattering on structures such as vessels fundamentally tied to the material's microstructure. Here, we present a scattering theory in conjunction with magnetic resonance imaging, which enables the unraveling of a material's innate constitutive and scattering characteristics. By overcoming a three-order-of-magnitude scale difference between wavelength and average intervessel distance, we provide noninvasively a macroscopic measure of vascular architecture. The validity of the theory is demonstrated through simulations, phantoms, in vivo mice, and human experiments and compared against histology as gold standard. Our approach expands the field of imaging by using the dispersion properties of shear waves as macroscopic observable proxies for deciphering the underlying ultrastructures.

Displacement and functional ultrasound (fUS) imaging of displacement-guided focused ultrasound (FUS) neuromodulation in mice

Focused ultrasound (FUS) stimulation is a promising neuromodulation technique with the merits of non-invasiveness, high spatial resolution, and deep penetration depth. However, simultaneous imaging of FUS-induced brain tissue displacement and the subsequent effect of FUS stimulation on brain hemodynamics has proven challenging thus far. In addition, earlier studies lack in situ confirmation of targeting except for the magnetic resonance imaging-guided FUS system-based studies. The purpose of this study is 1) to introduce a fully ultrasonic approach to in situ target, modulate neuronal activity, and monitor the resultant neuromodulation effect by respectively leveraging displacement imaging, FUS, and functional ultrasound (fUS) imaging, and 2) to investigate FUS-evoked cerebral blood volume (CBV) response and the relationship between CBV and displacement. We performed displacement imaging on craniotomized mice to confirm the in situ targeting for neuromodulation site. We recorded hemodynamic responses evoked by FUS while fUS imaging revealed an ipsilateral CBV increase that peaks at 4 s post-FUS. We report a stronger hemodynamic activation in the subcortical region than cortical, showing good agreement with a brain elasticity map that can also be obtained using a similar methodology. We observed dose-dependent CBV responses with peak CBV, activated area, and correlation coefficient increasing with the ultrasonic dose. Furthermore, by mapping displacement and hemodynamic activation, we found that displacement colocalized and linearly correlated with CBV increase. The findings presented herein demonstrated that FUS evokes ipsilateral hemodynamic activation in cortical and subcortical depths while the evoked hemodynamic responses colocalize and correlate with FUS-induced displacement. We anticipate that our findings will help consolidate accurate targeting as well as shedding light on one of the mechanisms behind FUS modulation, i.e., how FUS mechanically displaces brain tissue affecting cerebral hemodynamics and thereby its associated connectivity.

The Networking Brain: How Extracellular Matrix, Cellular Networks, and Vasculature Shape the In Vivo Mechanical Properties of the Brain

Mechanically, the brain is characterized by both solid and fluid properties. The resulting unique material behavior fosters proliferation, differentiation, and repair of cellular and vascular networks, and optimally protects them from damaging shear forces. Magnetic resonance elastography (MRE) is a noninvasive imaging technique that maps the mechanical properties of the brain in vivo. MRE studies have shown that abnormal processes such as neuronal degeneration, demyelination, inflammation, and vascular leakage lead to tissue softening. In contrast, neuronal proliferation, cellular network formation, and higher vascular pressure result in brain stiffening. In addition, brain viscosity has been reported to change with normal blood perfusion variability and brain maturation as well as disease conditions such as tumor invasion. In this article, the contributions of the neuronal, glial, extracellular, and vascular networks are discussed to the coarse-grained parameters determined by MRE. This reductionist multi-network model of brain mechanics helps to explain many MRE observations in terms of microanatomical changes and suggests that cerebral viscoelasticity is a suitable imaging marker for brain disease.

Functional ultrasound (fUS) imaging of displacement-guided focused ultrasound (FUS) neuromodulation in mice

Focused ultrasound (FUS) stimulation is a promising neuromodulation technique with the merits of non-invasiveness, high spatial resolution, and deep penetration depth. However, simultaneous imaging of FUS-induced brain tissue displacement and the subsequent effect of FUS stimulation on brain hemodynamics has proven challenging thus far. In addition, earlier studies lack in situ confirmation of targeting except for the magnetic resonance imaging-guided FUS system-based studies. The purpose of this study is 1) to introduce a fully ultrasonic approach to in situ target, modulate neuronal activity, and monitor the resultant neuromodulation effect by respectively leveraging displacement imaging, FUS, and functional ultrasound (fUS) imaging, and 2) to investigate FUS-evoked cerebral blood volume (CBV) response and the relationship between CBV and displacement. We performed displacement imaging on craniotomized mice to confirm the in targeting for neuromodulation site. We recorded hemodynamic responses evoked by FUS and fUS revealed an ipsilateral CBV increase that peaks at 4 s post-FUS. We saw a stronger hemodynamic activation in the subcortical region than cortical, showing good agreement with the brain elasticity map that can also be obtained using a similar methodology. We observed dose-dependent CBV response with peak CBV, activated area, and correlation coefficient increasing with ultrasonic dose. Furthermore, by mapping displacement and hemodynamic activation, we found that displacement colocalizes and linearly correlates with CBV increase. The findings presented herein demonstrated that FUS evokes ipsilateral hemodynamic activation in cortical and subcortical depths and the evoked hemodynamic responses colocalized and correlate with FUS-induced displacement. We anticipate that our findings will help consolidate accurate targeting as well as an understanding of how FUS displaces brain tissue and affects cerebral hemodynamics.

Mechanical behavior and microstructure of porcine brain tissues under pulsed electric fields

Pulsed electric fields are extensively utilized in clinical treatments, such as subthalamic deep brain stimulation, where electric field loading is in direct contact with brain tissue. However, the alterations in brain tissue's mechanical properties and microstructure due to changes in electric field parameters have not received adequate attention. In this study, the mechanical properties and microstructure of the brain tissue under pulsed electric fields were focused on. Herein, a custom indentation device was equipped with a module for electric field loading. Parameters such as pulse amplitude and frequency were adjusted. The results demonstrated that following an indentation process lasting 5 s and reaching a depth of 1000 μm, and a relaxation process of 175 s, the average shear modulus of brain tissue was reduced, and viscosity decreased. At the same amplitude, high-frequency pulsed electric fields had a smaller effect on brain tissue than low-frequency ones. Furthermore, pulsed electric fields induced cell polarization and reduced the proteoglycan concentration in brain tissue. As pulse frequency increased, cell polarization diminished, and proteoglycan concentration decreased significantly. High-frequency pulsed electric fields applied to brain tissue were found to reduce impedance fluctuation amplitude. This study revealed the effect of pulsed electric fields on the mechanical properties and microstructure of ex vivo brain tissue, providing essential information to promote the advancement of brain tissue electrotherapy in clinical settings.

Performance reserves in brain-imaging-based phenotype prediction

This study examines the impact of sample size on predicting cognitive and mental health phenotypes from brain imaging via machine learning. Our analysis shows a 3- to 9-fold improvement in prediction performance when sample size increases from 1,000 to 1 M participants. However, despite this increase, the data suggest that prediction accuracy remains worryingly low and far from fully exploiting the predictive potential of brain imaging data. Additionally, we find that integrating multiple imaging modalities boosts prediction accuracy, often equivalent to doubling the sample size. Interestingly, the most informative imaging modality often varied with increasing sample size, emphasizing the need to consider multiple modalities. Despite significant performance reserves for phenotype prediction, achieving substantial improvements may necessitate prohibitively large sample sizes, thus casting doubt on the practical or clinical utility of machine learning in some areas of neuroimaging.

Cerebrospinal fluid solute transport associated with sensorimotor brain activity in rodents

Cerebrospinal fluid (CSF) is crucial for maintaining neuronal homeostasis, providing nutrition, and removing metabolic waste from the brain. However, the relationship between neuronal activity and CSF solute transport remains poorly understood. To investigate the effect of regional neuronal activity on CSF solute transport, Sprague-Dawley rats (all male, n = 30) under anesthesia received an intracisternal injection of a fluorescent tracer (Texas Red ovalbumin) and were subjected to unilateral electrical stimulation of a forelimb. Two groups (n = 10 each) underwent two different types of stimulation protocols for 90 min, one including intermittent 7.5-s resting periods and the other without rest. The control group was not stimulated. Compared to the control, the stimulation without resting periods led to increased transport across most of the cortical areas, including the ventricles. The group that received intermittent stimulation showed an elevated level of solute uptake in limited areas, i.e., near/within the ventricles and on the ventral brain surface. Interhemispheric differences in CSF solute transport were also found in the cortical regions that overlap with the forelimb sensorimotor area. These findings suggest that neuronal activity may trigger local and brain-wide increases in CSF solute transport, contributing to waste clearance.

Ultra-wideband optical coherence elastography from acoustic to ultrasonic frequencies

Visualizing viscoelastic waves in materials and tissues through noninvasive imaging is valuable for analyzing their mechanical properties and detecting internal anomalies. However, traditional elastography techniques have been limited by a maximum wave frequency below 1-10 kHz, which hampers temporal and spatial resolution. Here, we introduce an optical coherence elastography technique that overcomes the limitation by extending the frequency range to MHz. Our system can measure the stiffness of hard materials including bones and extract viscoelastic shear moduli for polymers and hydrogels in conventionally inaccessible ranges between 100 Hz and 1 MHz. The dispersion of Rayleigh surface waves across the ultrawide band allowed us to profile depth-dependent shear modulus in cartilages ex vivo and human skin in vivo with sub-mm anatomical resolution. This technique holds immense potential as a noninvasive measurement tool for material sciences, tissue engineering, and medical diagnostics.

Transcranial focused ultrasound stimulation of cortical and thalamic somatosensory areas in human

The effects of transcranial focused ultrasound (FUS) stimulation of the primary somatosensory cortex and its thalamic projection (i.e., ventral posterolateral nucleus) on the generation of electroencephalographic (EEG) responses were evaluated in healthy human volunteers. Stimulation of the unilateral somatosensory circuits corresponding to the non-dominant hand generated EEG evoked potentials across all participants; however, not all perceived stimulation-mediated tactile sensations of the hand. These FUS-evoked EEG potentials (FEP) were observed from both brain hemispheres and shared similarities with somatosensory evoked potentials (SSEP) from median nerve stimulation. Use of a 0.5 ms pulse duration (PD) sonication given at 70% duty cycle, compared to the use of 1 and 2 ms PD, elicited more distinctive FEP peak features from the hemisphere ipsilateral to sonication. Although several participants reported hearing tones associated with FUS stimulation, the observed FEP were not likely to be confounded by the auditory sensation based on a separate measurement of auditory evoked potentials (AEP) to tonal stimulation (mimicking the same repetition frequency as the FUS stimulation). Off-line changes in resting-state functional connectivity (FC) associated with thalamic stimulation revealed that the FUS stimulation enhanced connectivity in a network of sensorimotor and sensory integration areas, which lasted for at least more than an hour. Clinical neurological evaluations, EEG, and neuroanatomical MRI did not reveal any adverse or unintended effects of sonication, attesting its safety. These results suggest that FUS stimulation may induce long-term neuroplasticity in humans, indicating its neurotherapeutic potential for various neurological and neuropsychiatric conditions.

Mechanical behavior of the hippocampus and corpus callosum: An attempt to reconcile ex vivo with in vivo and micro with macro properties

Mechanical properties of brain tissue are very complex and vary with the species, region, method, and dynamic range, and between in vivo and ex vivo measurements. To reconcile this variability, we investigated in vivo and ex vivo stiffness properties of two distinct regions in the human and mouse brain - the hippocampus (HP) and the corpus callosum (CC) - using different methods. Under quasi-static conditions, we examined ex vivo murine HP and CC by atomic force microscopy (AFM). Between 16 and 40Hz, we investigated the in vivo brains of healthy volunteers by magnetic resonance elastography (MRE) in a 3-T clinical scanner. At high-frequency stimulation between 1000 and 1400Hz, we investigated the murine HP and CC ex vivo and in vivo with MRE in a 7-T preclinical system. HP and CC showed pronounced stiffness dispersion, as reflected by a factor of 32-36 increase in shear modulus from AFM to low-frequency human MRE and a 25-fold higher shear wave velocity in murine MRE than in human MRE. At low frequencies, HP was softer than CC, in both ex vivo mouse specimens (p < 0.05) and in vivo human brains (p < 0.01) while, at high frequencies, CC was softer than HP under in vivo (p < 0.01) and ex vivo (p < 0.05) conditions. The standard linear solid model comprising three elements reproduced the observed HP and CC stiffness dispersions, while other two- and three-element models failed. Our results indicate a remarkable consistency of brain stiffness across species, ex vivo and in vivo states, and different measurement techniques when marked viscoelastic dispersion properties combining equilibrium and non-equilibrium mechanical elements are considered.

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-61 years) for analyzing the coefficient of variation (CV) and effects of biological factors. Eleven volunteers were rescanned after 1 day and 1 year 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.02 m/s) than 3D-MRE (1.22 ± 0.05 m/s, p < 0.0001), with less variation after 1 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 1 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 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.

Harmonic viscoelastic response of 3D histology-informed white matter model

White matter (WM) consists of bundles of long axons embedded in a glial matrix, which lead to anisotropic mechanical properties of brain tissue, and this complicates direct numerical simulations of WM viscoelastic response. The detailed axonal geometry contains scales that range from axonal diameter (microscale) to many diameters (mesoscale) imposing an additional challenge to numerical simulations. Here we describe the development of a 3D homogenization model for the central nervous system (CNS) that accounts for the anisotropy introduced by the axon/neuroglia composite, the axonal trace curvature, and the tissue dynamic response in the frequency domain. Homogenized models that allow the incorporation of all the above factors are important for accurately simulating the tissue's mechanical behavior, and this in turn is essential in interpreting non-invasive elastography measurements. Geometric and material parameters affect the material properties and thus the response of the brain tissue. More complex, orthotropic, or anisotropic material properties are to be considered as necessitated by the 3D tissue structure. An assembly of micro-scale 3D representative elemental volumes (REVs) is constructed, leading to an integrated mesoscale WM finite element model. Assemblies of microscopic REVs, with orientations based on geometrical reconstructions driven by confocal microscopy data are employed to form the elements of the WM model. Each REV carries local material properties based on a finite element model of biphasic (axon-glial matrix) unidirectional composite. The viscoelastic response of the microscopic REVs is extracted based on geometric information and fiber volume fractions calculated from the relative distance between the local elements and global axonal trace. The response of the WM tissue model is homogenized by averaging the shear moduli over the total volume (thus deriving effective properties) under realistic external loading conditions. Under harmonic shear loading, it is proven that that the effective transverse shear moduli are higher than the axial moduli when the axon moduli are higher than the glial. Methodologically, the process of using micro-scale 3D REVs to describe more complex axon geometries avoids the partition process in traditional composite finite element methods (based on partition of finite element grids) and constitutes a robust algorithm to automatically build a WM model based on available axonal trace information. Analytically, the model provides unmatched simulation flexibility and computational power as the position, orientation, and the magnitude of each tissue building block is calculated using real tissue data, as are the training and testing processes at each level of the multiscale WM tissue.

Theory of sleep/wake cycles affecting brain elastography

As elastography of the brain finds increasing clinical applications, fundamental questions remain about baseline viscoelastic properties of the brainin vivo. Furthermore, the underlying mechanisms of how and why elastographic measures can change over time are still not well understood. To study these issues, reverberant shear wave elastography using an optical coherence tomography scanner is implemented on a mouse model, both under awake conditions and in a sleep state where there are known changes in the glymphatic fluid flow system in the brain. We find that shear wave speed, a measure of stiffness, changes by approximately 12% between the two states, sleep versus awake, in the entire cortical brain imaging volume. Our microchannel flow model of biphasic (fluid plus solid) tissue provides a plausible rheological model based on the fractal branching vascular and perivascular system, plus a second parallel system representing the finer scale glymphatic fluid microchannels. By adjusting the glymphatic system fluid volume proportional to the known sleep/wake changes, we are able to approximately predict the measured shear wave speeds and their change with the state of the glymphatic system. The advantages of this model are that its main parameters are derived from anatomical measures and are linked to other major derivations of branching fluid structures including Murray's Law. The implications for clinical studies are that elastography of the brain is strongly influenced by the regulation or dysregulation of the vascular, perivascular, and glymphatic systems.

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.

A Minireview on Brain Models Simulating Geometrical, Physical, and Biochemical Properties of the Human Brain

There is a growing body of evidences that brain surrogates will be of great interest for researchers and physicians in the medical field. They are currently mainly used for education and training purposes or to verify the appropriate functionality of medical devices. Depending on the purpose, a variety of materials have been used with specific and accurate mechanical and biophysical properties, More recently they have been used to assess the biocompatibility of implantable devices, but they are still not validated to study the migration of leaching components from devices. This minireview shows the large diversity of approaches and uses of brain phantoms, which converge punctually. All these phantoms are complementary to numeric models, which benefit, reciprocally, of their respective advances. It also suggests avenues of research for the analysis of leaching components from implantable devices.

Simulating Local Deformations in the Human Cortex Due to Blood Flow-Induced Changes in Mechanical Tissue Properties: Impact on Functional Magnetic Resonance Imaging

Investigating human brain tissue is challenging due to the complexity and the manifold interactions between structures across different scales. Increasing evidence suggests that brain function and microstructural features including biomechanical features are related. More importantly, the relationship between tissue mechanics and its influence on brain imaging results remains poorly understood. As an important example, the study of the brain tissue response to blood flow could have important theoretical and experimental consequences for functional magnetic resonance imaging (fMRI) at high spatial resolutions. Computational simulations, using realistic mechanical models can predict and characterize the brain tissue behavior and give us insights into the consequent potential biases or limitations of in vivo, high-resolution fMRI. In this manuscript, we used a two dimensional biomechanical simulation of an exemplary human gyrus to investigate the relationship between mechanical tissue properties and the respective changes induced by focal blood flow changes. The model is based on the changes in the brain’s stiffness and volume due to the vasodilation evoked by neural activity. Modeling an exemplary gyrus from a brain atlas we assessed the influence of different potential mechanisms: (i) a local increase in tissue stiffness (at the level of a single anatomical layer), (ii) an increase in local volume, and (iii) a combination of both effects. Our simulation results showed considerable tissue displacement because of these temporary changes in mechanical properties. We found that the local volume increase causes more deformation and consequently higher displacement of the gyrus. These displacements introduced considerable artifacts in our simulated fMRI measurements. Our results underline the necessity to consider and characterize the tissue displacement which could be responsible for fMRI artifacts.

Aging brain mechanics: Progress and promise of magnetic resonance elastography

Neuroimaging techniques that can sensitivity characterize healthy brain aging and detect subtle neuropathologies have enormous potential to assist in the early detection of neurodegenerative conditions such as Alzheimer's disease. Magnetic resonance elastography (MRE) has recently emerged as a reliable, high-resolution, and especially sensitive technique that can noninvasively characterize tissue biomechanical properties (i.e., viscoelasticity) in vivo in the living human brain. Brain tissue viscoelasticity provides a unique biophysical signature of neuroanatomy that are representative of the composition and organization of the complex tissue microstructure. In this article, we detail how progress in brain MRE technology has provided unique insights into healthy brain aging, neurodegeneration, and structure-function relationships. We further discuss additional promising technical innovations that will enhance the specificity and sensitivity for brain MRE to reveal considerably more about brain aging as well as its potentially valuable role as an imaging biomarker of neurodegeneration. MRE sensitivity may be particularly useful for assessing the efficacy of rehabilitation strategies, assisting in differentiating between dementia subtypes, and in understanding the causal mechanisms of disease which may lead to eventual pharmacotherapeutic development.

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.

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.

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.

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

BACKGROUND

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.

PURPOSE

To measure MRE repeatability and estimate the effect of different reconstruction methods and varying data quality on estimated brain stiffness.

STUDY TYPE

Prospective.

SUBJECTS

Fifteen healthy subjects.

FIELD STRENGTH/SEQUENCE

3T MRI, gradient-echo elastography sequence with a 50 Hz vibration frequency.

ASSESSMENT

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.

STATISTICAL TESTS

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.

RESULTS

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.

DATA CONCLUSION

MRE of the human brain is a robust technique in terms of repeatability. Caution is warranted when comparing stiffness values obtained with different techniques.

LEVEL OF EVIDENCE

1 TECHNICAL EFFICACY STAGE: 1.

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.

Imaging brain function with simultaneous BOLD and viscoelasticity contrast: fMRI/fMRE

Magnetic resonance elastography (MRE) is emerging as a new tool for studying viscoelastic changes in the brain resulting from functional processes. Here, we demonstrate a novel time series method to generate robust functional magnetic resonance elastography (fMRE) activation maps in response to a visual task with a flashing checkerboard stimulus. Using a single-shot spin-echo (SS-SE) pulse sequence, the underlying raw images inherently contain blood-oxygen-level dependent (BOLD) contrast, allowing simultaneous generation of functional magnetic resonance imaging (fMRI) activation maps from the magnitude and functional magnetic resonance elastography (fMRE) maps from the phase. This allows an accurate comparison of the spatially localized stiffness (fMRE) and BOLD (fMRI) changes within a single scan, eliminating confounds inherent in separately acquired scans. Results indicate that tissue stiffness within the visual cortex increases 6-11% with visual stimuli, whereas the BOLD signal change was 1-2%. Furthermore, the fMRE and fMRI activation maps have strong spatial overlap within the visual cortex, providing convincing evidence that fMRE is possible in the brain. However, the fMRE temporal SNR (tSNR_fMRE) maps are heterogeneous across the brain. Using a dictionary matching approach to characterize the time series, the viscoelastic changes are consistent with a viscoelastic response function (VRF) time constant of 12.1 ​s ± 3.0 ​s for a first-order exponential decay, or a shape parameter of 8.1 ​s ± 1.4 ​s for a gamma-variate.

Investigating the Contribution of Collagen to the Tumor Biomechanical Phenotype with Noninvasive Magnetic Resonance Elastography

Increased stiffness in the extracellular matrix (ECM) contributes to tumor progression and metastasis. Therefore, stromal modulating therapies and accompanying biomarkers are being developed to target ECM stiffness. Magnetic resonance (MR) elastography can noninvasively and quantitatively map the viscoelastic properties of tumors in vivo and thus has clear clinical applications. Herein, we used MR elastography, coupled with computational histopathology, to interrogate the contribution of collagen to the tumor biomechanical phenotype and to evaluate its sensitivity to collagenase-induced stromal modulation. Elasticity (G d) and viscosity (G l) were significantly greater for orthotopic BT-474 (G d = 5.9 ± 0.2 kPa, G l = 4.7 ± 0.2 kPa, n = 7) and luc-MDA-MB-231-LM2-4 (G d = 7.9 ± 0.4 kPa, G l = 6.0 ± 0.2 kPa, n = 6) breast cancer xenografts, and luc-PANC1 (G d = 6.9 ± 0.3 kPa, G l = 6.2 ± 0.2 kPa, n = 7) pancreatic cancer xenografts, compared with tumors associated with the nervous system, including GTML/Trp53KI/KI medulloblastoma (G d = 3.5 ± 0.2 kPa, G l = 2.3 ± 0.2 kPa, n = 7), orthotopic luc-D-212-MG (G d = 3.5 ± 0.2 kPa, G l = 2.3 ± 0.2 kPa, n = 7), luc-RG2 (G d = 3.5 ± 0.2 kPa, G l = 2.3 ± 0.2 kPa, n = 5), and luc-U-87-MG (G d = 3.5 ± 0.2 kPa, G l = 2.3 ± 0.2 kPa, n = 8) glioblastoma xenografts, intracranially propagated luc-MDA-MB-231-LM2-4 (G d = 3.7 ± 0.2 kPa, G l = 2.2 ± 0.1 kPa, n = 7) breast cancer xenografts, and Th-MYCN neuroblastomas (G d = 3.5 ± 0.2 kPa, G l = 2.3 ± 0.2 kPa, n = 5). Positive correlations between both elasticity (r = 0.72, P < 0.0001) and viscosity (r = 0.78, P < 0.0001) were determined with collagen fraction, but not with cellular or vascular density. Treatment with collagenase significantly reduced G d (P = 0.002) and G l (P = 0.0006) in orthotopic breast tumors. Texture analysis of extracted images of picrosirius red staining revealed significant negative correlations of entropy with G d (r = -0.69, P < 0.0001) and G l (r = -0.76, P < 0.0001), and positive correlations of fractal dimension with G d (r = 0.75, P < 0.0001) and G l (r = 0.78, P < 0.0001). MR elastography can thus provide sensitive imaging biomarkers of tumor collagen deposition and its therapeutic modulation. SIGNIFICANCE: MR elastography enables noninvasive detection of tumor stiffness and will aid in the development of ECM-targeting therapies.

Viscoelasticity of striatal brain areas reflects variations in body mass index of lean to overweight male adults

Although a variety of MRI studies investigated the link between body mass index (BMI) and parameters of neural gray matter (GM), the technique applied in most of these studies, voxel-based morphometry (VBM), focusses on the regional GM volume, a macroscopic tissue property. Thus, the studies were not able to exploit the BMI-related information contained in the GM microstructure although PET studies suggest that these factors are important. Here, we used cerebral MR Elastography (MRE) to characterize features of tissue microstructure by evaluating the propagation of shear waves applied to the skull and to assess local tissue viscoelasticity to test the link between this parameter and BMI in 22 lean to overweight males. Unlike the majority of existing MRE studies investigating neural viscoelasticity signals averaged across large brain regions, we used the viscoelasticity of individual voxels for our experiment. Our technique revealed a negative link between BMI and viscoelasticity of two areas of the striatal reward system, i.e., right putamen (t = -8.2; p_{FWE-corrected} = 0.005) and left globus pallidus (t = -7.1; p_FWE = 0.037) which was independent of GM volume at these coordinates. Finally, comparison of BMI models based on individual voxels vs. on signals averaged across brain atlas regions demonstrates that voxel-based models explain a significantly higher proportion of variance. Consequently, our findings show that cerebral MRE is suitable to identify medically relevant microstructural tissue properties. Using a voxel-wise analysis approach, we were able to utilize the high spatial resolution of MRE for mapping BMI-related information in the brain.

A Systematic Review on the Possible Relationship Between Bilingualism, Cognitive Decline, and the Onset of Dementia

A systematic review was conducted to investigate whether bilingualism has a protective effect against cognitive decline in aging and can protect against dementia. We searched the Medline, ScienceDirect, Scopus, and ERIC databases with a cut-off date of 31 March, 2019, thereby following the guidelines of the Preferred Reporting Items for Systematic Reviews and Meta-analysis (PRISMA) protocol. Our search resulted in 34 eligible studies. Mixed results were found with respect to the protective effect of bilingualism against cognitive decline. Several studies showed a protective effect whereas other studies failed to find it. Moreover, evidence for a delay of the onset of dementia of between 4 and 5.5 years in bilingual individuals compared to monolinguals was found in several studies, but not in all. Methodological differences in the set-up of the studies seem to explain these mixed results. Lifelong bilingualism is a complex individual process, and many factors seem to influence this and need to be further investigated. This can be best achieved through large longitudinal studies with objective behavioral and neuroimaging measurements. In conclusion, although some evidence was found for a cognitive reserve-enhancing effect of lifelong bilingualism and protection against dementia, to date, no firm conclusions can be drawn.