Abstract
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.
1 Introduction
The brain is a pulsatile organ that expands by approximately half a milliliter due to the arterial pulse wave (). This implies that cerebral vessel walls compliantly expand during the passage of the arterial pulse wave and transmit vascular stroke energy to the surrounding tissue (Wagshul et al., 2011). The transfer of mechanical energy from cerebral vessels to brain tissue decreases the propagation velocity of the pressure wave throughout the vascular tree and contributes to a more continuous perfusion profile within cerebral capillaries. Conversely, decreased absorption of pulsation energy due to aging or pathologic changes in the vessel walls can produce long-term damage to vessel walls and impair the transport of oxygen and nutrients. Probably for these reasons, high-pulse wave velocity and decreased vascular compliance have been associated with a greater risk of dementia (Wolters et al., 2017; Rouch et al., 2018).
Direct measurement of vessel compliance is challenging because it requires knowledge about the mechanical stress associated with a measured change in the blood flow (Zamir et al., 2018). Transcranial Doppler ultrasound or tailored MRI techniques (Yan et al., 2016; Wolters et al., 2017; Zwanenburg and van Osch, 2017) can quantify cerebral blood flow (); however, mechanical stresses exerted by the pulse wave onto vessel walls and surrounding tissue are difficult to detect.
Here, we propose an alternative approach to study brain pulsatility and vascular compliance based on elastography. Elastography uses conventional imaging modalities to measure macroscopic tissue stiffness, for example, by sensing externally induced shear waves in soft tissues non-invasively (). In vascularized tissue, such as the brain, the acquired stiffness maps reflect the mechanical properties of the effective medium, which change with perfusion, microfluid flow, and pressure (). Our hypothesis is that cerebral arterial pulsation (CAP) induces measurable stiffness changes in the properties of the effective medium of the in vivo human brain.
Elastography of the brain has been mainly based on MRI (; Yin et al., 2018; ; ; Sack, 2022), while a few studies employed ultrasound (; Tzschatzsch et al., 2018) or optical methods (). Unlike ultrasound elastography, magnetic resonance elastography (MRE) does not encounter acoustic shading through the bones of the skull and can produce high-resolution maps of the stiffness of human brain tissue (). On the downside, MRE is expensive, limited in availability, and requires acquisition times on the order of minutes, which are too long to detect CAP-induced stiffness changes without cardiac gating (Schrank et al., 2019). Ultrasound-based elastography is cost-effective, widely available, and has short acquisition times, but has limited ability to generate and encode shear waves in the brain (; Shiina et al., 2015; ).
For the aforementioned reason, cerebral time-harmonic elastography (THE) has been developed based on sonography through the temples and acquisition of externally induced harmonic shear waves from 27 to 56 Hz, similar to the frequency range of MRE (Tzschatzsch et al., 2018). THE was shown to be sensitive to physiological alterations in cerebral blood flow due to the Valsalva maneuver (Tzschatzsch et al., 2018), hypercapnia (), or dehydration (), indicating the high sensitivity of brain stiffness to vascular compliance and intracranial pressure. In patients with intracranial hypertension, cerebral THE could detect stiffness changes before and after lumbar puncture (). However, owing to insufficient time resolution, none of these studies explored CAP-induced stiffness changes as a potential biomarker of cerebrovascular compliance.
Therefore, here, we developed time-resolved cerebral THE with a frame rate of 100 Hz to continuously capture the effects of CAP on brain stiffness over several heartbeats. In healthy volunteers, we performed two experiments. In one experiment, we used a high-end, plane-wave ultrasound scanner, which allowed us to acquire color flow images and stiffness maps in an interleaved fashion. This experiment was designed to ensure that established flow measurements and novel elastography captured the same physiological event on a shared time axis. In a second experiment, we used conventional line-by-line imaging to acquire only time-resolved stiffness changes without flow measurements, which facilitates clinical translation. Altogether, we aim to 1) report, for the first time, pulsation-induced changes in brain stiffness in the temporal lobe, basal cistern, and brain stem directly compared with cerebral blood flow and 2) demonstrate the implementation of this method on a clinical and cost-efficient ultrasound system based on conventional line-by-line imaging.
2 Materials and methods
2.1 Volunteers
The study was approved by the institutional review board of Charité—Universitätsmedizin Berlin (EA1/004/19). Written informed consent was obtained from all volunteers before participating in this study. Two groups of 10 healthy volunteers each (group 1: 3 women, age range: 21–52 years; group 2: 4 women, age range: 25–52 years) were prospectively included. Descriptive data of all volunteers (sex, age, BMI, heart rate, and mean arterial pressure) are summarized in Tables 1, 2. Inclusion criterion was no history of psychiatric disorders or head injury.
TABLE 1
| Volunteers | CBF | SWS | Delay time | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Sex | Age | BMI | HR | MAP | Mean | Systole | Diastole | Rel. | Mean | Systole | Diastole | Rel. | τR | τΔ |
| M/F | years | kg/m2 | bpm | mmHg | cm/s | cm/s | cm/s | % | m/s | m/s | m/s | % | ms | ms |
| M | 37 | 21.4 | 56 | 92 | 61 | 87 | 38 | 126 | 1.56 | 1.61 | 1.51 | 6.6 | 246 | 107 |
| M | 31 | 22.6 | 70 | 103 | 41 | 66 | 18 | 262 | 1.89 | 1.93 | 1.83 | 5.5 | 170 | −9 |
| M | 27 | 21.0 | 73 | 84 | 87 | 101 | 71 | 42 | 1.93 | 2.01 | 1.86 | 8.1 | 215 | 9 |
| M | 21 | 22.9 | 86 | 85 | 79 | 101 | 45 | 126 | 1.92 | 2.00 | 1.85 | 8.1 | 162 | −21 |
| M | 52 | 21.7 | 56 | 92 | 58 | 78 | 39 | 101 | 1.67 | 1.76 | 1.58 | 11.4 | 172 | 75 |
| F | 28 | 23.8 | 89 | 82 | 52 | 76 | 31 | 148 | 1.69 | 1.73 | 1.66 | 4.2 | 122 | 0 |
| M | 37 | 27.8 | 67 | 93 | 40 | 61 | 22 | 182 | 1.85 | 1.91 | 1.80 | 6.1 | 160 | −9 |
| M | 28 | 23.5 | 70 | 90 | 40 | 57 | 21 | 172 | 1.86 | 1.95 | 1.78 | 9.6 | 189 | −9 |
| F | 27 | 23.6 | 66 | 88 | 55 | 77 | 30 | 159 | 2.11 | 2.18 | 2.03 | 7.4 | 216 | 18 |
| F | 42 | 20.8 | 64 | 83 | 22 | 50 | 17 | 187 | 1.98 | 2.09 | 1.93 | 8.3 | 171 | 0 |
| Mean | 33 | 22.9 | 70 | 89 | 53 | 75 | 33 | 150 | 1.85 | 1.92 | 1.78 | 7.5 | 182 | 16 |
| SD | 9 | 1.9 | 10 | 6 | 18 | 17 | 15 | 55 | 0.15 | 0.16 | 0.15 | 2.0 | 34 | 39 |
Characteristics of the volunteers and individual CBF, global SWS, and peak delay time values detected by CBF-THE.
TABLE 2
| Volunteers | SWS | ||||||||
|---|---|---|---|---|---|---|---|---|---|
| Sex | Age | BMI | HR | MAP | Mean | Systole | Diastole | Rel. | τΔ |
| M/F | years | kg/m2 | bmp | mmHg | m/s | m/s | m/s | % | ms |
| F | 30 | 21 | 78 | 86 | 1.80 | 1.89 | 1.74 | 8.5 | 170 |
| M | 31 | 30 | 52 | 88 | 1.79 | 1.90 | 1.70 | 12.0 | 270 |
| M | 25 | 23 | 59 | 74 | 1.70 | 1.77 | 1.65 | 7.5 | 230 |
| M | 31 | 23 | 64 | 105 | 1.77 | 1.86 | 1.71 | 8.9 | 210 |
| M | 52 | 22 | 55 | 92 | 1.92 | 2.05 | 1.80 | 13.9 | 155 |
| F | 42 | 21 | 78 | 87 | 1.88 | 1.96 | 1.82 | 7.6 | 190 |
| F | 30 | 32 | 72 | 99 | 1.67 | 1.73 | 1.62 | 7.1 | 255 |
| M | 29 | 20 | 48 | 96 | 1.78 | 1.89 | 1.72 | 9.6 | 235 |
| M | 29 | 24 | 59 | 82 | 1.68 | 1.77 | 1.63 | 8.2 | 195 |
| F | 27 | 24 | 70 | 83 | 1.77 | 1.94 | 1.68 | 15.5 | 200 |
| Mean | 32 | 24 | 62 | 89 | 1.78 | 1.88 | 1.71 | 9.9 | 211 |
| SD | 8 | 4 | 11 | 8 | 0.08 | 0.10 | 0.07 | 2.9 | 36 |
Characteristics and SWS values of volunteers acquired by c-THE.
2.2 Study protocol
In the first experiment, high-end, plane-wave imaging (Vantage 64LE, Verasonics, Washington, United States) for interleaved acquisitions of cerebral blood flow (CBF) using color flow imaging and shear wave speed (SWS) by THE was applied to group 1 (CBF-THE). In the second experiment, cost-efficient, line-by-line imaging (GAMPT mbH, Merseburg, Germany) was applied to the second group of volunteers for stiffness mapping (c-THE). In both CBF-THE and c-THE experiments, ECG data were continuously recorded for assigning CBF and SWS to cardiac phases. Each measurement took 8 s and was repeated five times to increase the statistical power.
2.3 Data acquisition
Volunteers were investigated in a supine position on a custom-designed vibration bed (GAMPT mbH, Merseburg, Germany) with the head placed above the vibration unit for efficient transfer of harmonic vibration energy into the brain (Figure 1A) (; ; ). To avoid limitations in time resolution due to the decomposition of multifrequency wave setups, we used only single-frequency harmonic stimulations at 60 Hz, in contrast to previous THE studies. A high-end Vantage 64LE scanner (Verasonics, Washington, United States) equipped with a phased array transducer (P4-2v) with a center frequency of 3 MHz was used for CBF-THE. CBF (in cm/s) and SWS (in m/s) were acquired in an interleaved fashion (Figure 1B) with an effective frame rate of 100 Hz and up to a depth of 10 cm. The transducer was positioned at the temples, and an image plane in the temporal lobe of the brain was selected using color flow imaging of the middle cerebral artery overlaid to the anatomical B-mode (Figure 1C). For the elastography measurement, plane-wave compounding was used (seven plane-wave angles, −5° to 5°, with a pulse repetition frequency of 4 kHz). For color flow imaging, a burst of 40 frames was used (pulse repetition frequency of 5.5 kHz). For c-THE measurement, elastography data were acquired using a line-by-line ultrasound scanner (GAMPT mbH, Merseburg, Germany), as shown in Figure 1A, which was equipped with a phased-array transducer (P5-1S15-A6) having a center frequency of 3 MHz. The selected field-of-view was in the same region as the field-of-view in CBF-THE, with the same depth of up to 10 cm but a higher frame rate of 200 Hz and 26 lines of sight (Figure 1B), resulting in a slightly decreased insonation sector (Figure 1C). For this study, the entire data processing was newly developed to allow time-resolved SWS analysis.
FIGURE 1
2.4 Data processing
THE processing relied on the phase-difference extraction of the local axial displacement between adjacent in-phase and quadrature (IQ) frames. The temporal Fourier transform was then applied to the computed displacement data, and time-resolved complex wavefields (Figure 1C) were obtained using the filtered Hilbert transform with a bandpass filter centered at 60 Hz (Figure 1C). Due to undersampling, the vibration signal of 60 Hz acquired at a frame rate of 100 Hz is wrapped to the 40-Hz signal in the spectrum. The use of this wrapped signal instead of the original frequency is referred to as “controlled aliasing,” as proposed in our previous work (Tzschatzsch et al., 2016). The basic temporal resolution of the complex wavefields is related to both the frame rate of the full image acquisition (100 Hz) and the bandwidth of the Hilbert transformation. The bandpass filter employed in the Hilbert transform was Gaussian shaped with a width of 10 Hz and, thus, represents an empirically found tradeoff between noise robustness and high temporal resolution.
As the final step, time-resolved SWS maps were generated using wave-number based k-MDEV inversion (; Tzschatzsch et al., 2016). The inversion assumes plane shear waves with complex wave number (), angular frequency ω, and amplitude u0 as solutions of the following wave equation:
k-MDEV first decomposed the complex-valued wavefields obtained from the Hilbert transformation into eight complex wavefields for each propagation direction (d). For this, Gaussian-based spatiotemporal filters with a low-pass threshold of 555 rad/m for noise suppression and first-order derivative for suppression of compression waves were applied. The wave numbers were deduced from the phase gradient of :
The final maps of SWS were obtained by averaging over directions:
Processing took approximately 15 min for the CBF-THE data and 30 s for the c-THE data. Regions of interest (ROIs) for the temporal lobe, basal cistern, and brain stem were manually drawn based on the anatomical B-mode image. SWS variations over time were analyzed after averaging within these regions for each time step. For color flow imaging processing, a region covering the middle cerebral artery was automatically selected through thresholding. A threshold of 5 cm/s was empirically set based on the amplitude of the background noise. A temporal moving boxcar filter of 3-pixel width was applied for noise suppression. The envelope of the flow spectrum was computed to obtain the time variation in CBF within the middle cerebral artery. R–R intervals were extracted from the acquired ECG signals with wavelet-based peak detection () using the Haar wavelet (Sunkaria et al., 2010). For each volunteer, a normalized time axis covering one R–R interval was generated onto which the average of all THE and CBF signals acquired over 8 s and five repetitions was mapped. The resulting time curve of each volunteer and ROI were used to deduce the peak and minimum CBF and SWS values corresponding to cerebral peak systole and end-diastole, respectively. To quantify the pulsation, the relative difference between peak-systolic and end-diastolic SWS and CBF referenced to the end-diastolic value was computed. For comparing with the literature, the mean CBF and SWS values were calculated. Furthermore, the time delay of the SWS peak relative to the R-wave (τR) and relative to the CBF peak (τΔ) was quantified and tabulated.
2.5 Statistics
Relative changes in SWS and CBF over the R–R-interval were tested for each subregion using the paired Wilcoxon signed rank test. Delay times τR and τΔ were tested against zero using the Wilcoxon signed rank test. Pearson’s correlation coefficient was used to analyze possible correlations between the relative change in SWS and delay times τR and τΔ in the global ROI, with demographic parameters, such as age, BMI, heart rate, and mean arterial pressure, and the relative change in CBF. The relative change in SWS, mean SWS, and delay time τR obtained from CBF-THE was compared to c-THE using the Wilcoxon signed rank test. Nonparametric tests were chosen due to the small number of volunteers. All p values <0.05 were considered significant. The group size was estimated based on a power analysis assuming an effect size of 10%–20% at a 5% significance level and a power of 0.8 for the variability of stiffness parameters encountered in healthy adult subjects.
3 Results
In all volunteers, time-resolved SWS and, in the CBF-THE group, time-resolved CBF were obtained. Figure 2 shows the representative CBF and SWS data obtained by CBF-THE. The regions of interest were demarcated based on the anatomical landmarks in the B-mode image. Additionally, the time variations in SWS and CBF averaged within the global ROI are shown over a measure time of 8 s. The periodic change in both parameters in synchrony to the ECG suggests a relationship between brain stiffness and arterial blood flow in the temporal lobe.
FIGURE 2
Group SWS and CBF values on a time axis normalized to one R–R-interval are presented for all four regions in Figure 3. As in Figure 2, both parameters vary synchronously over the cardiac phase, with higher relative changes and less group variation according to confidence intervals in CBF than those in SWS. The mean SWS was 1.85 ± 0.15 m/s in the global ROI, 1.85 ± 0.16 m/s in the temporal parenchyma, 2.05 ± 0.39 m/s in the basal cisterns, and 1.48 ± 0.29 m/s in the brain stem, whereas the mean CBF (54 ± 19 cm/s) was in the range of literature values (), and SWS was slightly higher than that in previous cerebral THE studies due to the higher excitation frequency ().
FIGURE 3
In all regions, a significant relative decrease from peak systolic to end diastolic SWS values was observed (Figure 4A), with 7.5% ± 2.0% in the global ROI, 4.8% ± 1.8% in the temporal parenchyma, 11% ± 5% in the basal cisterns, and 13% ± 9% in the brain stem. For CBF, a corresponding relative decrease of 139% ± 59% was observed, as shown in Table 1. Time delays between the R peak and peak systolic SWS (τR) were 180 ± 30 ms in the global ROI, 210 ± 50 ms in the temporal parenchyma, 180 ± 60 ms in the basal cisterns, and 160 ± 40 ms in the brain stem. In all regions except the temporal lobe, there was no significant time delay between the peak systolic SWS and peak systolic CBF (τΔ). In the temporal lobe, τΔ of 50 ± 50 ms was observed, suggesting a slight delay in the SWS wave after the CBF wave (Figure 4B). No significant correlations were observed between the mean SWS and physiological parameters such as age, sex, and BMI. Moreover, no correlation was found between the relative SWS changes and blood pressure and delay times and heart rate. All SWS and CBF values and delay times acquired by CBF-THE are summarized in Table 1.
FIGURE 4
Group SWS within the global ROI measured by c-THE was similar to CBF-THE, with a mean SWS of 1.78 ± 0.08 m/s, relative SWS changes of 9.9% ± 2.9%, and a delay in the SWS peak to the R peak of 211 ± 36 ms. There was no significant difference in the values obtained by CBF THE (mean SWS: p = 0.21, relative SWS changes: p = 0.06, τR: p = 0.13). Data acquired by c-THE are summarized in Table 2.
4 Discussion
To the best of our knowledge, this is the first study in which changes in brain stiffness due to cerebral arterial pulsation have been directly observed without cardiac gating and time averaging, that is, quasi-real time. We observed a change in SWS in synchrony to CBF pulsation, from which we concluded that CAP significantly affected brain stiffness. This effect was present in all regions, which underscores the relevance of the relative SWS change as a potential biomarker of cerebrovascular compliance. It should be noted that THE was less constrained by accurately locating anatomic landmarks than color flow imaging, which is limited to the middle cerebral artery. As a result, the overall time required, including probe positioning, was shorter for THE than for color flow imaging. Because there was no significant temporal shift between the SWS peak and the CBF peak of arterial pulsation, we can conclude that SWS measurement provides a robust and cost-effective imaging marker, that is, complementary to Doppler imaging.
Given the extensive literature on elastography in highly perfused organs, we speculate that the observed variation in brain stiffness is related to cerebrovascular compliance (; ; ; ; ; ; ; ; Solamen et al., 2021). Specifically, brain stiffness has been shown by MRE and THE to increase with cerebral artery dilation due to hypercapnia (; ) or the Valsalva maneuver (Tzschatzsch et al., 2018; ). As shown by MRE, cerebral stiffness of deep gray matter correlates with perfusion pressure rather than cerebral blood flow alone (). This means that the vessel walls dilate during the passage of the pulse wave, causing a global stiffening of the tissue in which the vessels are embedded. In simple terms, when vascular pressure increases in encapsulated (non-expandible) organs such as the brain or liver, whether by increasing the total blood volume due to water uptake (; ) or by transient passage of the arterial pulse wave as observed in this study, the tissue stiffens. Conversely, tissue softening has been reported in the spleen, which may expand with increased blood volume, likely due to the increased volume of fluid in the tissue with no change in pressure (; ).
An exception here is our own previous work in which cardiac-gated MRE was used to measure CAP-induced stiffness changes (Schrank et al., 2019). In this study, a series of stiffness maps were generated, each composed of signals from 81 heartbeats, which, in combination with finger pulse triggering, led to uncertainties in the timing of stiffness peaks relative to CAP. The partially different results of this previous study compared with those of the current study emphasize the relevance of the simultaneous measurement of blood flow and stiffness—as currently only possible by THE in a plane-wave setup. Nevertheless, the question remains unresolved as to why THE is so sensitive to intracranial pressure, such as in patients before and after lumbar puncture (), whereas mixed results have been reported based on MRE ().
Further studies investigating the effects of different states of perfusion and pressure on the temporal modulation of cerebral stiffness are required to better understand the observed stiffness variation. In this context, the effects of dehydration and water uptake on brain stiffness have been studied based on cerebral THE (). However, these studies were limited to time-averaged measurements, as opposed to the time-resolved measurements made possible through this study. Additionally, it would be interesting to see if an altered cerebral blood flow in patients with vascular male formations, vascular dementia, or intracranial hypertension changes the observed pulsation patterns (; ; ).
The main limitation to our study is the lack of clinical data. It would be interesting to see how stiffness pulsation changes upon the relief of intracranial hypertension in patients undergoing lumbar puncture. A possible change in stiffness pulsation over a wider age range and in conditions such as mild cognitive impairment and dementia is also warranted for future studies. In addition, the relatively small number of subjects precludes a conclusive analysis of physiological effects on changes in brain stiffness, as well as the assessment of possible timing differences between CBF and SWS waves. The insignificant correlation of τR with the heart rate may be due to the fact that we only investigated healthy volunteers with a relatively narrow range of heart rates. In future studies including patients with tachy- or bradycardia, we expect τR to be sensitive to the heart rate. However, this study focused on method development and first-time demonstration of CAP-induced stiffness alterations to provide a reference for larger cohort studies in a clinical setting.
5 Conclusion
In summary, our study demonstrated the feasibility of time-resolved cerebral THE for the assessment of cerebrovascular compliance based on pulsations of tissue stiffness. Plane-wave imaging-based THE showed that brain stiffness in healthy volunteers changes from minimum values during cerebral diastole to maximum values during cerebral systole in the order of 10% with variation across anatomical regions. This direct relationship between the stiffness wave and arterial blood pulse wave could make THE a sensitive imaging marker for brain compliance changes due to cerebrovascular pulsation, which could be used for diagnostics and treatment monitoring. Furthermore, relative stiffness changes measured by line-by-line-based THE was shown to be a cost-effective alternative to plane-wave-based THE that is ready to use in the clinic for the noninvasive measurement of brain stiffness associated with intracranial pressure and cerebrovascular compliance.
Statements
Data availability statement
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving human participants were reviewed and approved by the Ethics Committee of Charité—Universitätsmedizin Berlin (EA1/004/19). The patients/participants provided their written informed consent to participate in this study.
Author contributions
TM, BK, JBe, MA, BW, JBr, MD, HT, and IS contributed to the conception and design of the study. TM, EA, and BK contributed to the technical implementation. BK and JBe performed the data acquisition. TM and BK performed the data evaluation and statistical analysis. TM and BK wrote the first draft of the manuscript. TM, BK, HT, and IS wrote sections of the manuscript. All authors contributed to the article and approved the submitted version.
Funding
The authors declare that this study received funding from Pfizer. The funder was not involved in the study design, collection, analysis, interpretation of data, the writing of this article, or the decision to submit it for publication. Support of the German Research Foundation (DFG, project number 460333672—CRC 1540 Exploring Brain Mechanics (subproject X01), project number TZ 98/1-1, CRC 1340 Matrix in Vision, GRK 2260 BIOQIC) and Pfizer (WP2487656) is gratefully acknowledged.
Conflict of interest
The authors declare that this study received funding from Pfizer. The authors declare that the research was conducted in the absence of any other commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
AdamsA. L.KuijfH. J.ViergeverM. A.LuijtenP. R.ZwanenburgJ. J. M. (2019). Quantifying cardiac-induced brain tissue expansion using DENSE. NMR Biomed.32 (2), e4050. 10.1002/nbm.4050
2
AraniA.MinH. K.FattahiN.WetjenN. M.TrzaskoJ. D.ManducaA.et al (2018). Acute pressure changes in the brain are correlated with MR elastography stiffness measurements: initial feasibility in an in vivo large animal model. Magn. Reson Med.79 (2), 1043–1051. 10.1002/mrm.26738
3
CosgroveD.BarrR.BojungaJ.CantisaniV.ChammasM. C.DigheM.et al (2017). WFUMB guidelines and recommendations on the clinical use of ultrasound elastography: part 4. Thyroid. Ultrasound Med. Biol.43 (1), 4–26. 10.1016/j.ultrasmedbio.2016.06.022
4
CuiweiL.ChongxunZ.ChangfengT. (1995). Detection of ECG characteristic points using wavelet transforms. IEEE Trans. Biomed. Eng.42 (1), 21–28. 10.1109/10.362922
5
ErtlM.RaaschN.HammelG.HarterK.LangC. (2018). Transtemporal investigation of brain parenchyma elasticity using 2-D shear wave elastography: definition of age-matched normal values. Ultrasound Med. Biol.44 (1), 78–84. 10.1016/j.ultrasmedbio.2017.08.1885
6
FerraioliG.FiliceC.CasteraL.ChoiB. I.SporeaI.WilsonS. R.et al (2015). WFUMB guidelines and recommendations for clinical use of ultrasound elastography: part 3: liver. Ultrasound Med. Biol.41 (5), 1161–1179. 10.1016/j.ultrasmedbio.2015.03.007
7
ForouhandehpourR.BernierM.GilbertG.ButlerR.WhittingstallK.Van HoutenE. (2021). Cerebral stiffness changes during visual stimulation: differential physiological mechanisms characterized by opposing mechanical effects. Neuroimage Rep.1 (2), 100014. 10.1016/j.ynirp.2021.100014
8
GeG. R.SongW.NedergaardM.RollandJ. P.ParkerK. J. (2022). Theory of sleep/wake cycles affecting brain elastography. Phys. Med. Biol.67 (22), 225013. 10.1088/1361-6560/ac9e40
9
HerthumH.HetzerS.KreftB.TzschätzschH.ShahryariM.MeyerT.et al (2022). Cerebral tomoelastography based on multifrequency MR elastography in two and three dimensions. Front. Bioeng. Biotechnol.10, 1056131. 10.3389/fbioe.2022.1056131
10
HerthumH.ShahryariM.TzschatzschH.SchrankF.WarmuthC.GornerS.et al (2021). Real-time multifrequency MR elastography of the human brain reveals rapid changes in viscoelasticity in response to the Valsalva maneuver. Front. Bioeng. Biotechnol.9 (335), 666456. 10.3389/fbioe.2021.666456
11
HetzerS.BirrP.FehlnerA.HirschS.DittmannF.BarnhillE.et al (2018). Perfusion alters stiffness of deep gray matter. J. Cereb. Blood Flow. Metab.38 (1), 116–125. 10.1177/0271678x17691530
12
HetzerS.DittmannF.BormannK.HirschS.LippA.WangD. J.et al (2019). Hypercapnia increases brain viscoelasticity. J. Cereb. Blood Flow. Metab.39 (12), 2445–2455. 10.1177/0271678x18799241
13
HirschS.GuoJ.ReiterR.PapazoglouS.KroenckeT.BraunJ.et al (2014a). MR elastography of the liver and the spleen using a piezoelectric driver, single-shot wave-field acquisition, and multifrequency dual parameter reconstruction. Magn. Reson Med.71 (1), 267–277. 10.1002/mrm.24674
14
HirschS.GuoJ.ReiterR.SchottE.BuningC.SomasundaramR.et al (2014b). Towards compression-sensitive magnetic resonance elastography of the liver: sensitivity of harmonic volumetric strain to portal hypertension. J. Magn. Reson Imaging39 (2), 298–306. 10.1002/jmri.24165
15
HirschS.KlattD.FreimannF.ScheelM.BraunJ.SackI. (2012). In vivo measurement of volumetric strain in the human brain induced by arterial pulsation and harmonic waves. Magn. Reson Med.70 (3), 671–683. 10.1002/mrm.24499
16
HiscoxL. V.JohnsonC. L.BarnhillE.McGarryM. D.HustonJ.van BeekE. J.et al (2016). Magnetic resonance elastography (MRE) of the human brain: technique, findings and clinical applications. Phys. Med. Biol.61 (24), R401–R437. 10.1088/0031-9155/61/24/r401
17
Ipek-UgayS.TzschatzschH.HudertC.Marticorena GarciaS. R.FischerT.BraunJ.et al (2016). Time harmonic elastography reveals sensitivity of liver stiffness to water ingestion. Ultrasound Med. Biol.42 (6), 1289–1294. 10.1016/j.ultrasmedbio.2015.12.026
18
KolipakaA.WassenaarP. A.ChaS.MarashdehW. M.MoX.KalraP.et al (2018). Magnetic resonance elastography to estimate brain stiffness: measurement reproducibility and its estimate in pseudotumor cerebri patients. Clin. Imaging51, 114–122. 10.1016/j.clinimag.2018.02.005
19
KrakauskaiteS.ThibeaultC.LaVangieJ.ScheidtM.MartinezL.Seth-HunterD.et al (2018). Normative ranges of transcranial Doppler metrics. Acta Neurochir. Suppl.126, 269–273. 10.1007/978-3-319-65798-1_53
20
KreftB.BergsJ.ShahryariM.DanyelL. A.HetzerS.BraunJ.et al (2020b). Cerebral ultrasound time-harmonic elastography reveals softening of the human brain due to dehydration. Front. Physiol.11, 616984. 10.3389/fphys.2020.616984
21
KreftB.TzschatzschH.SchrankF.BergsJ.StreitbergerK. J.WaldchenS.et al (2020a). Time-resolved response of cerebral stiffness to hypercapnia in humans. Ultrasound Med. Biol.46 (4), 936–943. 10.1016/j.ultrasmedbio.2019.12.019
22
KreftB.TzschatzschH.ShahryariM.HaffnerP.BraunJ.SackI.et al (2022). Noninvasive detection of intracranial hypertension by novel ultrasound time-harmonic elastography. Invest. Radiol.57 (2), 77–84. 10.1097/rli.0000000000000817
23
LauriolaM.MangiacottiA.D'OnofrioG.CascavillaL.ParisF.ParoniG.et al (2018). Neurocognitive disorders and dehydration in older patients: clinical experience supports the hydromolecular hypothesis of dementia. Nutrients10 (5), 562. 10.3390/nu10050562
24
LiuH.WangD.LengX.ZhengD.ChenF.WongL. K. S.et al (2020). State-of-the-Art computational models of circle of willis with physiological applications: A review. IEEE Access8, 156261–156273. 10.1109/access.2020.3007737
25
ManducaA.BaylyP. J.EhmanR. L.KolipakaA.RoystonT. J.SackI.et al (2021). MR elastography: principles, guidelines, and terminology. Magn. Reson Med.85 (5), 2377–2390. 10.1002/mrm.28627
26
Marticorena GarciaS. R.GrossmannM.LangS. T.TzschatzschH.DittmannF.HammB.et al (2018). Tomoelastography of the native kidney: regional variation and physiological effects on in vivo renal stiffness. Magn. Reson Med.79 (4), 2126–2134. 10.1002/mrm.26892
27
MeyerT.TzschatzschH.WellgeB.SackI.KronckeT.MartlA. (2022). Valsalva maneuver decreases liver and spleen stiffness measured by time-harmonic ultrasound elastography. Front. Bioeng. Biotechnol.10, 886363. 10.3389/fbioe.2022.886363
28
NanjappaM.KolipakaA. (2021). Magnetic resonance elastography of the brain. Magnetic Reson. Imaging Clin. N. Am.29 (4), 617–630. 10.1016/j.mric.2021.06.011
29
OrmacheaJ.ParkerK. J. (2020). Elastography imaging: the 30 year perspective. Phys. Med. Biol.65 (24), 1361–6560. 10.1088/1361-6560/abca00
30
ParkerK. J. (2014). A microchannel flow model for soft tissue elasticity. Phys. Med. Biol.59 (15), 4443–4457. 10.1088/0031-9155/59/15/4443
31
ParkerK. J. (2015). Experimental evaluations of the microchannel flow model. Phys. Med. Biol.60 (11), 4227–4242. 10.1088/0031-9155/60/11/4227
32
PoulS. S.OrmacheaJ.HollenbachS. J.ParkerK. J. (2020). Validations of the microchannel flow model for characterizing vascularized tissues. Fluids5 (4), 228. 10.3390/fluids5040228
33
RobbaC.CardimD.SekhonM.BudohoskiK.CzosnykaM. (2018). Transcranial Doppler: a stethoscope for the brain-neurocritical care use. J. Neurosci. Res.96 (4), 720–730. 10.1002/jnr.24148
34
RouchL.CestacP.SallerinB.AndrieuS.BaillyH.BeunardeauM.et al (2018). Pulse wave velocity is associated with greater risk of dementia in mild cognitive impairment patients. Hypertension72 (5), 1109–1116. 10.1161/hypertensionaha.118.11443
35
SackI. (2022). Magnetic resonance elastography from fundamental soft-tissue mechanics to diagnostic imaging. Nat. Rev. Phys.5, 25–42. 10.1038/s42254-022-00543-2
36
SchrankF.WarmuthC.TzschatzschH.KreftB.HirschS.BraunJ.et al (2019). Cardiac-gated steady-state multifrequency magnetic resonance elastography of the brain: effect of cerebral arterial pulsation on brain viscoelasticity. J. Cereb. Blood Flow. Metab.40, 991–1001. 10.1177/0271678x19850936
37
ShiinaT.NightingaleK. R.PalmeriM. L.HallT. J.BamberJ. C.BarrR. G.et al (2015). WFUMB guidelines and recommendations for clinical use of ultrasound elastography: part 1: basic principles and terminology. Ultrasound Med. Biol.41 (5), 1126–1147. 10.1016/j.ultrasmedbio.2015.03.009
38
SolamenL. M.McGarryM. D. J.FriedJ.WeaverJ. B.LollisS. S.PaulsenK. D. (2021). Poroelastic mechanical properties of the brain tissue of normal pressure hydrocephalus patients during lumbar drain treatment using intrinsic actuation MR elastography. Acad. Radiol.28 (4), 457–466. 10.1016/j.acra.2020.03.009
39
SunkariaR. K.SaxenaS. C.KumarV.SinghalA. M. (2010). Wavelet based R-peak detection for heart rate variability studies. J. Med. Eng. Technol.34 (2), 108–115. 10.3109/03091900903281215
40
TzschatzschH.KreftB.SchrankF.BergsJ.BraunJ.SackI. (2018). In vivo time-harmonic ultrasound elastography of the human brain detects acute cerebral stiffness changes induced by intracranial pressure variations. Sci. Rep.8 (1), 17888. 10.1038/s41598-018-36191-9
41
TzschatzschH.Nguyen TrongM.ScheuermannT.Ipek-UgayS.FischerT.SchultzM.et al (2016). Two-dimensional time-harmonic elastography of the human liver and spleen. Ultrasound Med. Biol.42 (11), 2562–2571. 10.1016/j.ultrasmedbio.2016.07.004
42
WagshulM.EideP.MadsenJ. (2011). The pulsating brain: A review of experimental and clinical studies of intracranial pulsatility. Fluids Barriers CNS8, 5. 10.1186/2045-8118-8-5
43
WoltersF. J.ZonneveldH. I.HofmanA.van der LugtA.KoudstaalP. J.VernooijM. W.et al (2017). Cerebral perfusion and the risk of dementia A population-based study. Circulation136 (8), 719–728. 10.1161/circulationaha.117.027448
44
YanL.LiuC. Y.SmithR. X.JogM.LanghamM.KrasilevaK.et al (2016). Assessing intracranial vascular compliance using dynamic arterial spin labeling. Neuroimage124, 433–441. 10.1016/j.neuroimage.2015.09.008
45
YinZ.RomanoA. J.ManducaA.EhmanR. L.HustonJ.3rd (2018). Stiffness and beyond: what MR elastography can tell us about brain structure and function under physiologic and pathologic conditions. Top. Magn. Reson Imaging27 (5), 305–318. 10.1097/rmr.0000000000000178
46
ZamirM.MoirM. E.KlassenS. A.BalestriniC. S.ShoemakerJ. K. (2018). Cerebrovascular compliance within the rigid confines of the skull. Front. Physiology9, 940. 10.3389/fphys.2018.00940
47
ZwanenburgJ. J. M.van OschM. J. P. (2017). Targeting cerebral small vessel disease with MRI. Stroke48 (11), 3175–3182. 10.1161/strokeaha.117.016996
Summary
Keywords
ultrasound time-harmonic elastography, in vivo brain stiffness, brain pulsation, cerebrovascular compliance, intracranial pressure
Citation
Meyer T, Kreft B, Bergs J, Antes E, Anders MS, Wellge B, Braun J, Doyley M, Tzschätzsch H and Sack I (2023) Stiffness pulsation of the human brain detected by non-invasive time-harmonic elastography. Front. Bioeng. Biotechnol. 11:1140734. doi: 10.3389/fbioe.2023.1140734
Received
09 January 2023
Accepted
19 July 2023
Published
15 August 2023
Volume
11 - 2023
Edited by
Salavat Aglyamov, University of Houston, United States
Reviewed by
Haipeng Liu, Coventry University, United Kingdom
Manmohan Singh, University of Houston, United States
Updates
Copyright
© 2023 Meyer, Kreft, Bergs, Antes, Anders, Wellge, Braun, Doyley, Tzschätzsch and Sack.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ingolf Sack, Ingolf.sack@charite.de
† These authors have contributed equally to this work and share first authorship
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.