Quantifying cardiac-induced brain tissue expansion using DENSE

3D Quantitative-Amplified Magnetic Resonance Imaging (3D q-aMRI)

Amplified MRI (aMRI) is a promising new technique that can visualize pulsatile brain tissue motion by amplifying sub-voxel motion in cine MRI data, but it lacks the ability to quantify the sub-voxel motion field in physical units. Here, we introduce a novel post-processing algorithm called 3D quantitative amplified MRI (3D q-aMRI). This algorithm enables the visualization and quantification of pulsatile brain motion. 3D q-aMRI was validated and optimized on a 3D digital phantom and was applied in vivo on healthy volunteers for its ability to accurately measure brain parenchyma and CSF voxel displacement. Simulation results show that 3D q-aMRI can accurately quantify sub-voxel motions in the order of 0.01 of a voxel size. The algorithm hyperparameters were optimized and tested on in vivo data. The repeatability and reproducibility of 3D q-aMRI were shown on six healthy volunteers. The voxel displacement field extracted by 3D q-aMRI is highly correlated with the displacement measurements estimated by phase contrast (PC) MRI. In addition, the voxel displacement profile through the cerebral aqueduct resembled the CSF flow profile reported in previous literature. Differences in brain motion was observed in patients with dementia compared with age-matched healthy controls. In summary, 3D q-aMRI is a promising new technique that can both visualize and quantify pulsatile brain motion. Its ability to accurately quantify sub-voxel motion in physical units holds potential for the assessment of pulsatile brain motion as well as the indirect assessment of CSF homeostasis. While further research is warranted, 3D q-aMRI may provide important diagnostic information for neurological disorders such as Alzheimer's disease.

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.

Dynamic light-responsive RhoA activity regulates mechanosensitive stem cell fate decision in 3D matrices

The behavior of stem cells is regulated by mechanical cues in their niche that continuously vary due to extracellular matrix (ECM) remodeling, pulsated mechanical stress exerted by blood flow, and/or cell migration. However, it is still unclear how dynamics of mechanical cues influence stem cell lineage commitment, especially in a 3D microenvironment where mechanosensing differs from that in a 2D microenvironment. In the present study, we investigated how temporally varying mechanical signaling regulates expression of the early growth response 1 gene (Egr1), which we recently discovered to be a 3D matrix-specific mediator of mechanosensitive neural stem cell (NSC) lineage commitment. Specifically, we temporally controlled the activity of Ras homolog family member A (RhoA), which is known to have a central role in mechanotransduction, using our previously developed Arabidopsis thaliana cryptochrome-2-based optoactivation system. Interestingly, pulsed RhoA activation induced Egr1 upregulation in stiff 3D gels only, whereas static light stimulation induced an increase in Egr1 expression across a wide range of 3D gel stiffnesses. Actin assembly inhibition limited Egr1 upregulation upon RhoA activation, implying that RhoA signaling requires an actin-involved process to upregulate Egr1. Consistently, static-light RhoA activation rather than pulsed-light activation restricted neurogenesis in soft gels. Our findings indicate that the dynamics of RhoA activation influence Egr1-mediated stem cell fate within 3D matrices in a matrix stiffness-dependent manner.

Accelerated multislice MRI with patterned excitation

PURPOSE

Accelerate multislice 2D MRI by using RF pulses that simultaneously act on different slices to combine contrast preparation and image acquisition.

THEORY AND METHODS

MRI applications often require the use of multiple RF pulses to generate desired contrast and prepare the signal for readout. Examples are the use of inversion prepulses to generate T1 contrast, or the use of spin-echo preparations to generate T2 or diffusion contrast. In multislice MRI, this separation of contrast preparation and readout can render scans time-inefficient and lengthy. We introduce a class of pulse sequences that overcomes this inefficiency by combining contrast preparation and signal readout. This is accomplished by using RF pulses that manipulate the magnetization of multiple slices simultaneously and a gradient crusher scheme that selects a target slice for readout.

RESULTS

Feasibility of the method was demonstrated for spin echo-based measurement of water diffusion and tissue pulsation in human brain at 3 T. Increases in time-efficiency and reductions in scan time were highly dependent on specific implementation and reached as high as 25% and 53%, respectively.

CONCLUSION

A novel approach to multislice MRI is demonstrated that reduces scan time for specific applications.

Estimating the viscoelastic properties of the human brain at 7 T MRI using intrinsic MRE and nonlinear inversion

Intrinsic actuation magnetic resonance elastography (MRE) is a phase-contrast MRI technique that allows for in vivo quantification of mechanical properties of the brain by exploiting brain motion that arise naturally due to the cardiac pulse. The mechanical properties of the brain reflect its tissue microstructure, making it a potentially valuable parameter in studying brain disease. The main purpose of this study was to assess the feasibility of reconstructing the viscoelastic properties of the brain using high-quality 7 T MRI displacement measurements, obtained using displacement encoding with stimulated echoes (DENSE) and intrinsic actuation. The repeatability and sensitivity of the method for detecting normal regional variation in brain tissue properties was assessed as secondary goal. The displacement measurements used in this analysis were previously acquired for a separate study, where eight healthy subjects (27 ± 7 years) were imaged with repeated scans (spatial resolution approx. 2 mm isotropic, temporal resolution 75 ms, motion sensitivity 0.35 mm/2π for displacements in anterior-posterior and left-right directions, and 0.7 mm/2π for feet-head displacements). The viscoelastic properties of the brain were estimated using a subzone based non-linear inversion scheme. The results show comparable consistency to that of extrinsic MRE between the viscoelastic property maps obtained from repeated displacement measurements. The shear stiffness maps showed fairly consistent spatial patterns. The whole-brain repeatability coefficient (RC) for shear stiffness was (mean ± standard deviation) 8 ± 8% relative to the mean whole-brain stiffness, and the damping ratio RC was 28 ± 17% relative to the whole-brain damping ratio. The shear stiffness maps showed similar statistically significant regional trends as demonstrated in a publicly available atlas of viscoelastic properties obtained with extrinsic actuation MRE at 50 Hz. The damping ratio maps showed less consistency, likely due to data-model mismatch of describing the brain as a viscoelastic material under low frequencies. While artifacts induced by fluid flow within the brain remain a limitation of the technique in its current state, intrinsic actuation based MRE allow for consistent and repeatable estimation of the mechanical properties of the brain. The method provides enough sensitivity to investigate regional variation in such properties in the normal brain, which is likely sufficient to also investigate pathological changes.

Association between resistance to cerebrospinal fluid flow and cardiac-induced brain tissue motion for Chiari malformation type I

PURPOSE

Chiari malformation type I (CMI) patients have been independently shown to have both increased resistance to cerebrospinal fluid (CSF) flow in the cervical spinal canal and greater cardiac-induced neural tissue motion compared to healthy controls. The goal of this paper is to determine if a relationship exists between CSF flow resistance and brain tissue motion in CMI subjects.

METHODS

Computational fluid dynamics (CFD) techniques were employed to compute integrated longitudinal impedance (ILI) as a measure of unsteady resistance to CSF flow in the cervical spinal canal in thirty-two CMI subjects and eighteen healthy controls. Neural tissue motion during the cardiac cycle was assessed using displacement encoding with stimulated echoes (DENSE) magnetic resonance imaging (MRI) technique.

RESULTS

The results demonstrate a positive correlation between resistance to CSF flow and the maximum displacement of the cerebellum for CMI subjects (r = 0.75, p = 6.77 × 10-10) but not for healthy controls. No correlation was found between CSF flow resistance and maximum displacement in the brainstem for CMI or healthy subjects. The magnitude of resistance to CSF flow and maximum cardiac-induced brain tissue motion were not statistically different for CMI subjects with and without the presence of five CMI symptoms: imbalance, vertigo, swallowing difficulties, nausea or vomiting, and hoarseness.

CONCLUSION

This study establishes a relationship between CSF flow resistance in the cervical spinal canal and cardiac-induced brain tissue motion in the cerebellum for CMI subjects. Further research is necessary to understand the importance of resistance and brain tissue motion in the symptomatology of CMI.

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.

The Relationship Between Imbalance Symptom and Cardiac Pulsation Induced Mechanical Strain in the Brainstem and Cerebellum for Chiari Malformation Type I

Chiari malformation Type I (CMI) is known to have an altered biomechanical environment for the brainstem and cerebellum; however, it is unclear whether these altered biomechanics play a role in the development of CMI symptoms. We hypothesized that CMI subjects have a higher cardiac-induced strain in specific neurological tracts pertaining to balance, and postural control. We measured displacement over the cardiac cycle using displacement encoding with stimulated echoes magnetic resonance imaging in the cerebellum, brainstem, and spinal cord in 37 CMI subjects and 25 controls. Based on these measurements, we computed strain, translation, and rotation in tracts related to balance. The global strain on all tracts was small (<1%) for CMI subject and controls. Strain was found to be nearly doubled in three tracts for CMI subjects compared to controls (p < 0.03). The maximum translation and rotation were ∼150 μm and ∼1 deg, respectively and 1.5-2 times greater in CMI compared to controls in four tracts (p < 0.005). There was no significant difference between strain, translation, and rotation on the analyzed tracts in CMI subjects with imbalance compared to those without imbalance. A moderate correlation was found between cerebellar tonsillar position and strain on three tracts. The lack of statistically significant difference between strain in CMI subjects with and without imbalance could imply that the magnitude of the observed cardiac-induced strain was too small to cause substantial damage to the tissue (<1%). Activities such as coughing, or Valsalva may produce a greater strain.

Abnormalities in cardiac-induced brain tissue deformations are now detectable with MRI: A case-report of a patient who underwent craniotomy after trauma

BACKGROUND

Heartbeat and respiration induce cyclic brain tissue deformations, which receive increasing attention as potential driving force for brain clearance. These deformations can now be assessed using a novel 3D strain tensor imaging (STI) method at 7 T MRI.

METHODS

An 18-year-old man had suffered a traumatic brain injury and was treated with a craniotomy with a maximal diameter of 12 cm. STI was employed to capture cardiac-induced brain tissue deformations and additional time-resolved 2D flow measurements were acquired to capture cerebrospinal fluid (CSF) flow towards the spinal canal.

RESULTS

The craniotomy caused major changes in all aspects of the brain's mechanical dynamics as compared to healthy volunteer references. Tissue strains increased, particularly around the craniotomy, and directionality of deformations showed large abnormalities, also in the contralateral hemisphere. As the brain tissue could pulsate outward from the skull, physiological pulsatile CSF flow at the foramen magnum was abolished.

CONCLUSIONS

This work illustrates how STI can assess physiological patterns of brain tissue deformation and how craniotomy leads to widespread deformation abnormalities that can be detected at a single patient level. While this case is meant to provide proof of concept, application of STI in other conditions of abnormal brain mechanical dynamics warrants further study.

Magnetic resonance imaging of the pulsing brain: a systematic review

OBJECTIVE

To perform a systematic review of the literature exploring magnetic resonance imaging (MRI) methods for measuring natural brain tissue pulsations (BTPs) in humans.

METHODS

A prospective systematic search of MEDLINE, SCOPUS and OpenGrey databases was conducted by two independent reviewers using a pre-determined strategy. The search focused on identifying reported measurements of naturally occurring BTP motion in humans. Studies involving non-human participants, MRI in combination with other modalities, MRI during invasive procedures and MRI studies involving externally applied tests were excluded. Data from the retrieved records were combined to create Forest plots comparing brain tissue displacement between Chiari-malformation type 1 (CM-I) patients and healthy controls using an independent samples t-test.

RESULTS

The search retrieved 22 eligible articles. Articles described 5 main MRI techniques for visualisation or quantification of intrinsic brain motion. MRI techniques generally agreed that the amplitude of BTPs varies regionally from 0.04 mm to ~ 0.80 mm, with larger tissue displacements occurring closer to the centre and base of the brain compared to peripheral regions. Studies of brain pathology using MRI BTP measurements are currently limited to tumour characterisation, idiopathic intracranial hypertension (IIH), and CM-I. A pooled analysis confirmed that displacement of tissue in the cerebellar tonsillar region of CM-I patients was + 0.31 mm [95% CI 0.23, 0.38, p < 0.0001] higher than in healthy controls.

DISCUSSION

MRI techniques used for measurements of brain motion are at an early stage of development with high heterogeneity across the methods used. Further work is required to provide normative data to support systematic BTPs characterisation in health and disease.

Dynamic brain ADC variations over the cardiac cycle and their relation to tissue strain assessed with DENSE at high-field MRI

Purpose

The ADC of brain tissue slightly varies over the cardiac cycle. This variation could reflect physiology, including mixing of the interstitial fluid, relevant for brain waste clearance. However, it is known from cardiac diffusion imaging that tissue deformation by itself affects the magnitude of the MRI signal, leading to artificial ADC variations as well. This study investigates to what extent tissue deformation causes artificial ADC variations in the brain.

Theory and Methods

We implemented a high‐field MRI sequence with stimulated echo acquisition mode that simultaneously measures brain tissue deformation and ADC. Based on the measured tissue deformation, we simulated the artificial ADC variation by combining established theoretical frameworks and compared the results with the measured ADC variation. We acquired data in 8 healthy volunteers with diffusion weighting b = 300 and b = 1000 s/mm².

Results

Apparent diffusion coefficient variation was largest in the feet‐to‐head direction and showed the largest deviation from the mean ADC at peak systole. Artificial ADC variation estimated from tissue deformation was 1.3 ± 0.37·10⁻⁵ mm²/s in the feet‐to‐head direction for gray matter, and 0.75 ± 0.29·10⁻⁵ mm²/s for white matter. The measured ADC variation in the feet‐to‐head direction was 5.6·10⁻⁵ ± 1.5·10⁻⁵ mm²/s for gray matter and 3.2·10⁻⁵ ± 1.0·10⁻⁵ mm²/s for white matter, which was a factor of 3.5 ± 0.82 and 3.4 ± 0.57 larger than the artificial diffusion variations. The measured diffusion variations in the right‐to‐left/anterior‐to‐posterior direction were a factor of 1.5 ± 1.0/1.7 ± 1.4 and 2.0 ± 0.91/2.5 ± 0.94 larger than the artificial diffusion variations for gray matter and white matter, respectively.

Conclusion

Apparent diffusion coefficient variations in the brain likely largely reflect physiology.

Medullary vein architecture modulates the white matter BOLD cerebrovascular reactivity signal response to CO2: Observations from high-resolution T2* weighted imaging at 7T

Brain stress testing using blood oxygenation level-dependent (BOLD) MRI to evaluate changes in cerebrovascular reactivity (CVR) is of growing interest for evaluating white matter integrity. However, even under healthy conditions, the white matter BOLD-CVR response differs notably from that observed in the gray matter. In addition to actual arterial vascular control, the venous draining topology may influence the WM-CVR response leading to signal delays and dispersions. These types of alterations in hemodynamic parameters are sometimes linked with pathology, but may also arise from differences in normal venous architecture. In this work, high-resolution T2*weighted anatomical images combined with BOLD imaging during a hypercapnic breathing protocol were acquired using a 7 tesla MRI system. Hemodynamic parameters including base CVR, hemodynamic lag, lag-corrected CVR, response onset and signal dispersion, and finally ΔCVR (corrected CVR minus base CVR) were calculated in 8 subjects. Parameter maps were spatially normalized and correlated against an MNI-registered white matter medullary vein atlas. Moderate correlations (Pearson's rho) were observed between medullary vessel frequency (MVF) and ΔCVR (0.52; 0.58 for total WM), MVF and hemodynamic lag (0.42; 0.54 for total WM), MVF and signal dispersion (0.44; 0.53 for total WM), and finally MVF and signal onset (0.43; 0.52 for total WM). Results indicate that, when assessed in the context of the WM venous architecture, changes in the response shape may only be partially reflective of the actual vascular reactivity response occurring further upstream by control vessels. This finding may have implications when attributing diseases mechanisms and/or progression to presumed impaired WM BOLD-CVR.

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.

Cerebellar and Brainstem Displacement Measured with DENSE MRI in Chiari Malformation Following Posterior Fossa Decompression Surgery

Background Posterior fossa decompression (PFD) surgery is a treatment for Chiari malformation type I (CMI). The goals of surgery are to reduce cerebellar tonsillar crowding and restore posterior cerebral spinal fluid flow, but regional tissue biomechanics may also change. MRI-based displacement encoding with stimulated echoes (DENSE) can be used to assess neural tissue displacement. Purpose To assess neural tissue displacement by using DENSE MRI in participants with CMI before and after PFD surgery and examine associations between tissue displacement and symptoms. Materials and Methods In a prospective, HIPAA-compliant study of patients with CMI, midsagittal DENSE MRI was performed before and after PFD surgery between January 2017 and June 2020. Peak tissue displacement over the cardiac cycle was quantified in the cerebellum and brainstem, averaged over each structure, and compared before and after surgery. Paired t tests and nonparametric Wilcoxon signed-rank tests were used to identify surgical changes in displacement, and Spearman correlations were determined between tissue displacement and presurgery symptoms. Results Twenty-three participants were included (mean age ± standard deviation, 37 years ± 10; 19 women). Spatially averaged (mean) peak tissue displacement demonstrated reductions of 46% (79/171 µm) within the cerebellum and 22% (46/210 µm) within the brainstem after surgery (P < .001). Maximum peak displacement, calculated within a circular 30-mm² area, decreased by 64% (274/427 µm) in the cerebellum and 33% (100/300 µm) in the brainstem (P < .001). No significant associations were identified between tissue displacement and CMI symptoms (r < .74 and P > .012 for all; Bonferroni-corrected P = .0002). Conclusion Neural tissue displacement was reduced after posterior fossa decompression surgery, indicating that surgical intervention changes brain tissue biomechanics. For participants with Chiari malformation type I, no relationship was identified between presurgery tissue displacement and presurgical symptoms. © RSNA, 2021 Online supplemental material is available for this article.

Regional Brain Tissue Displacement and Strain is Elevated in Subjects with Chiari Malformation Type I Compared to Healthy Controls: A Study Using DENSE MRI

While the degree of cerebellar tonsillar descent is considered the primary radiologic marker of Chiari malformation type I (CMI), biomechanical forces acting on the brain tissue in CMI subjects are less studied and poorly understood. In this study, regional brain tissue displacement and principal strains in 43 CMI subjects and 25 controls were quantified using a magnetic resonance imaging (MRI) methodology known as displacement encoding with stimulated echoes (DENSE). Measurements from MRI were obtained for seven different brain regions-the brainstem, cerebellum, cingulate gyrus, corpus callosum, frontal lobe, occipital lobe, and parietal lobe. Mean displacements in the cerebellum and brainstem were found to be 106 and 64% higher, respectively, for CMI subjects than controls (p < .001). Mean compression and extension strains in the cerebellum were 52 and 50% higher, respectively, in CMI subjects (p < .001). Brainstem mean extension strain was 41% higher in CMI subjects (p < .001), but no significant difference in compression strain was observed. The other brain structures revealed no significant differences between CMI and controls. These findings demonstrate that brain tissue displacement and strain in the cerebellum and brainstem might represent two new biomarkers to distinguish between CMI subjects and controls.

Cortical thinning is associated with brain pulsatility in older adults: An MRI and NIRS study

Aging is accompanied by global brain atrophy occurring unequally across the brain. Cortical thinning is seen with aging with a larger loss in the frontal and temporal subregions. We explored the link between regional cortical thickness and regional cerebral pulsatility. Sixty healthy individuals were divided into two age groups, young (aged 19-31) and older (aged 65-75) adults. Each participant underwent a near-infrared spectroscopy (NIRS) scan to index regional brain pulsatility from cerebral pulse-transit-time-to-the peak-of-the-pulse (PTTp), an anatomical magnetic resonance imaging (MRI) and a phase-contrast MRI (PC-MRI) scan to measure arterial and cerebrospinal fluid (CSF) pulsatility. In older adults, the greatest association between cerebral pulsatility and cortical thickness was found in superior and middle temporal and superior, middle and inferior frontal areas, which are the regions perfused first by the internal carotid arteries. This association dropped in the postcentral and superior parietal regions. These findings suggest higher brain pulsatility as a potential risk factor contributing to cortical thinning for some brain regions more than others.

3D amplified MRI (aMRI)

Purpose

Amplified MRI (aMRI) has been introduced as a new method of detecting and visualizing pulsatile brain motion in 2D. Here, we improve aMRI by introducing a novel 3D aMRI approach.

Methods

3D aMRI was developed and tested for its ability to amplify sub‐voxel motion in all three directions. In addition, 3D aMRI was qualitatively compared to 2D aMRI on multi‐slice and 3D (volumetric) balanced steady‐state free precession cine data and phase contrast (PC‐MRI) acquired on healthy volunteers at 3T. Optical flow maps and 4D animations were produced from volumetric 3D aMRI data.

Results

3D aMRI exhibits better image quality and fewer motion artifacts compared to 2D aMRI. The tissue motion was seen to match that of PC‐MRI, with the predominant brain tissue displacement occurring in the cranial‐caudal direction. Optical flow maps capture the brain tissue motion and display the physical change in shape of the ventricles by the relative movement of the surrounding tissues. The 4D animations show the complete brain tissue and cerebrospinal fluid (CSF) motion, helping to highlight the “piston‐like” motion of the ventricles.

Conclusions

Here, we introduce a novel 3D aMRI approach that enables one to visualize amplified cardiac‐ and CSF‐induced brain motion in striking detail. 3D aMRI captures brain motion with better image quality than 2D aMRI and supports a larger amplification factor. The optical flow maps and 4D animations of 3D aMRI may open up exciting applications for neurological diseases that affect the biomechanics of the brain and brain fluids.

Strain Tensor Imaging: Cardiac-induced brain tissue deformation in humans quantified with high-field MRI

The cardiac cycle induces blood volume pulsations in the cerebral microvasculature that cause subtle deformation of the surrounding tissue. These tissue deformations are highly relevant as a potential source of information on the brain's microvasculature as well as of tissue condition. Besides, cyclic brain tissue deformations may be a driving force in clearance of brain waste products. We have developed a high-field magnetic resonance imaging (MRI) technique to capture these tissue deformations with full brain coverage and sufficient signal-to-noise to derive the cardiac-induced strain tensor on a voxel by voxel basis, that could not be assessed non-invasively before. We acquired the strain tensor with 3 mm isotropic resolution in 9 subjects with repeated measurements for 8 subjects. The strain tensor yielded both positive and negative eigenvalues (principle strains), reflecting the Poison effect in tissue. The principle strain associated with expansion followed the known funnel shaped brain motion pattern pointing towards the foramen magnum. Furthermore, we evaluate two scalar quantities from the strain tensor: the volumetric strain and octahedral shear strain. These quantities showed consistent patterns between subjects, and yielded repeatable results: the peak systolic volumetric strain (relative to end-diastolic strain) was 4.19⋅10^-4 ± 0.78⋅10^-4 and 3.98⋅10^-4 ± 0.44⋅10^-4 (mean ± standard deviation for first and second measurement, respectively), and the peak octahedral shear strain was 2.16⋅10^-3 ± 0.31⋅10^-3 and 2.31⋅10^-3 ± 0.38⋅10^-3, for the first and second measurement, respectively. The volumetric strain was typically highest in the cortex and lowest in the periventricular white matter, while anisotropy was highest in the subcortical white matter and basal ganglia. This technique thus reveals new, regional information on the brain's cardiac-induced deformation characteristics, and has the potential to advance our understanding of the role of microvascular pulsations in health and disease.

Accuracy of cardiac-induced brain motion measurement using displacement-encoding with stimulated echoes (DENSE) magnetic resonance imaging (MRI): A phantom study

PURPOSE

The goal of this study was to determine the accuracy of displacement-encoding with stimulated echoes (DENSE) MRI in a tissue motion phantom with displacements representative of those observed in human brain tissue.

METHODS

The phantom was comprised of a plastic shaft rotated at a constant speed. The rotational motion was converted to a vertical displacement through a camshaft. The phantom generated repeatable cyclical displacement waveforms with a peak displacement ranging from 92 µm to 1.04 mm at 1-Hz frequency. The surface displacement of the tissue was obtained using a laser Doppler vibrometer (LDV) before and after the DENSE MRI scans to check for repeatability. The accuracy of DENSE MRI displacement was assessed by comparing the laser Doppler vibrometer and DENSE MRI waveforms.

RESULTS

Laser Doppler vibrometer measurements of the tissue motion demonstrated excellent cycle-to-cycle repeatability with a maximum root mean square error of 9 µm between the ensemble-averaged displacement waveform and the individual waveforms over 180 cycles. The maximum difference between DENSE MRI and the laser Doppler vibrometer waveforms ranged from 15 to 50 µm. Additionally, the peak-to-peak difference between the 2 waveforms ranged from 1 to 18 µm.

CONCLUSION

Using a tissue phantom undergoing cyclical motion, we demonstrated the percent accuracy of DENSE MRI to measure displacement similar to that observed for in vivo cardiac-induced brain tissue.

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

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

Cardiac and respiration-induced brain deformations in humans quantified with high-field MRI

Microvascular blood volume pulsations due to the cardiac and respiratory cycles induce brain tissue deformation and, as such, are considered to drive the brain's waste clearance system. We have developed a high-field magnetic resonance imaging (MRI) technique to quantify both cardiac and respiration-induced tissue deformations, which could not be assessed noninvasively before. The technique acquires motion encoded snapshot images in which various forms of motion and confounders are entangled. First, we optimized the motion sensitivity for application in the human brain. Next, we isolated the heartbeat and respiration-related deformations, by introducing a linear model that fits the snapshot series to the recorded physiological information. As a result, we obtained maps of the physiological tissue deformation with 3mm isotropic spatial resolution. Heartbeat and respiration-induced volumetric strain were significantly different from zero in the basal ganglia (median (25-75% interquartile range): 0.85·10^-3 (0.39·10^-3-1.05·10^-3), p ​= ​0.0008 and -0.28·10^-3 (-0.41·10^-3-0.06·10^-3), p ​= ​0.047, respectively. Smaller volumetric strains were observed in the white matter of the centrum semi ovale (0.28·10^-3 (0-0.59·10^-3) and -0.06·10^-3 (-0.17·10^-3-0.20·10^-3)), which was only significant for the heartbeat (p ​= ​0.02 and p ​= ​0.7, respectively). Furthermore, heartbeat-induced volumetric strain was about three times larger than respiration-induced volumetric strain. This technique opens a window on the driving forces of the human brain clearance system.

Validating faster DENSE measurements of cardiac-induced brain tissue expansion as a potential tool for investigating cerebral microvascular pulsations

Displacement Encoding with Stimulated Echoes (DENSE) has recently shown potential for measuring cardiac-induced cerebral volumetric strain in the human brain. As such, it may provide a powerful tool for investigating the cerebral small vessels. However, further development and validation are necessary. This study aims, first, to validate a retrospectively-gated implementation of the DENSE method for assessing brain tissue pulsations as a physiological marker, and second, to use the acquired measurements to explore intracranial volume dynamics. We acquired repeated measurements of cerebral volumetric strain in 8 healthy subjects, and internally validated these measurements by comparing them to spinal CSF stroke volumes obtained in the same scan session. Peak volumetric strain was found to be highly repeatable between scan sessions. First/second measured peak volumetric strains were: (6.4 ​± ​1.7)x10^-4/(6.7 ​± ​1.6)x10^-4 for whole brain, (9.5 ​± ​2.5)x10^-4/(9.6 ​± ​2.4)x10^-4 for grey matter, and (4.4 ​± ​1.7)x10^-4/(4.1 ​± ​0.8)x10^-4 for white matter. Grey matter showed significantly higher peak strain (p ​< ​0.001) and earlier time-to-peak strain (p ​< ​0.02) than white matter. An approximately linear relationship was found between CSF and brain tissue volume pulsations over the cardiac cycle (mean slope and R² of 0.88 ​± ​0.23 and 0.89 ​± ​0.07, respectively). The close similarity between CSF and brain tissue volume pulsations implies limited contributions from large intracranial vessel pulsations, providing further evidence for venous compression as an additional mechanism for maintaining stable intracranial pressures over the cardiac cycle. Cerebral pulsatility showed consistent inter-subject peak values in healthy subjects, and was strongly correlated to CSF stroke volumes. These results strengthen the potential of brain tissue volumetric strain as a means for investigating the intracranial dynamics of the ageing brain in normal or diseased states.

Novel strain analysis informs about injury susceptibility of the corpus callosum to repeated impacts

Increasing evidence for the cumulative effects of head trauma on structural integrity of the brain has emphasized the need to understand the relationship between tissue mechanic properties and injury susceptibility. Here, diffusion tensor imaging, helmet accelerometers and amplified magnetic resonance imaging were combined to gather insight about the region-specific vulnerability of the corpus callosum to microstructural changes in white-matter integrity upon exposure to sub-concussive impacts. A total of 33 male Canadian football players (mean_age = 20.3 ± 1.4 years) were assessed at three time points during a football season (baseline pre-season, mid-season and post-season). The athletes were split into a LOW (N = 16) and HIGH (N = 17) exposure group based on the frequency of sub-concussive impacts sustained on a per-session basis, measured using the helmet-mounted accelerometers. Longitudinal decreases in fractional anisotropy were observed in anterior and posterior regions of the corpus callosum (average cluster size = 40.0 ± 4.4 voxels; P < 0.05, corrected) for athletes from the HIGH exposure group. These results suggest that the white-matter tract may be vulnerable to repetitive sub-concussive collisions sustained over the course of a football season. Using these findings as a basis for further investigation, a novel exploratory analysis of strain derived from sub-voxel motion of brain tissues in response to cardiac impulses was developed using amplified magnetic resonance imaging. This approach revealed specific differences in strain (and thus possibly stiffness) along the white-matter tract (P < 0.0001) suggesting a possible signature relationship between changes in white-matter integrity and tissue mechanical properties. In light of these findings, additional information about the viscoelastic behaviour of white-matter tissues may be imperative in elucidating the mechanisms responsible for region-specific differences in injury susceptibility observed, for instance, through changes in microstructural integrity following exposure to sub-concussive head impacts.