Skip to main content

ORIGINAL RESEARCH article

Front. Bioeng. Biotechnol., 15 August 2023
Sec. Biomechanics
Volume 11 - 2023 | https://doi.org/10.3389/fbioe.2023.1140734

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

www.frontiersin.orgTom Meyer1 www.frontiersin.orgBernhard Kreft1 www.frontiersin.orgJudith Bergs1 www.frontiersin.orgErik Antes1 www.frontiersin.orgMatthias S. Anders1 www.frontiersin.orgBrunhilde Wellge1 www.frontiersin.orgJürgen Braun2 www.frontiersin.orgMarvin Doyley3 www.frontiersin.orgHeiko Tzschätzsch1 www.frontiersin.orgIngolf Sack1*
  • 1Department of Radiology, Charité—University Medicine Berlin, Berlin, Germany
  • 2Institute of Medical Informatics, Charité—University Medicine Berlin, Berlin, Germany
  • 3Hajim School of Engineering and Applied Sciences, University of Rochester, Rochester, NY, United States

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 (Adams et al., 2019). 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 (Robba et al., 2018); 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 (Ormachea and Parker, 2020). 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 (Poul et al., 2020). 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 (Hiscox et al., 2016; Yin et al., 2018; Nanjappa and Kolipaka, 2021; Herthum et al., 2022; Sack, 2022), while a few studies employed ultrasound (Ertl et al., 2018; Tzschatzsch et al., 2018) or optical methods (Ge et al., 2022). 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 (Manduca et al., 2021). 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 (Ferraioli et al., 2015; Shiina et al., 2015; Cosgrove et al., 2017).

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 (Kreft et al., 2020a), or dehydration (Kreft et al., 2020b), 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 (Kreft et al., 2022). 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
www.frontiersin.org

TABLE 1. Characteristics of the volunteers and individual CBF, global SWS, and peak delay time values detected by CBF-THE.

TABLE 2
www.frontiersin.org

TABLE 2. 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) (Kreft et al., 2020a; Kreft et al., 2020b; Kreft et al., 2022). 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
www.frontiersin.org

FIGURE 1. (A) Experimental setup for time-harmonic elastography (THE) with a vibration bed (GAMPT mbH, Merseburg, Germany), high-end, plane-wave ultrasound scanner (Vantage 64LE, Verasonics, Washington, United States) and cost-efficient, line-by-line THE device (GAMPT mbH, Merseburg, Germany). Volunteers were positioned in a supine position with the vibration unit underneath the back of their head. Measurements were performed through the temporal bone window, similar to clinical transcranial Doppler ultrasound. (B) Sequence diagram in the top shows the interleaved acquisition of THE and color flow imaging and in the bottom shows the line-by-line acquisition scheme. One image for plane-wave THE was composed of seven plane-wave angles, while a color flow image was composed of 40 frames. c-THE required 26 lines of sight (LOSs). (C) Representative wave images for both CBF-THE and c-THE, as well as the representative B-mode with CFI for CBF-THE.

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 (Cuiwei et al., 1995; Tzschatzsch et al., 2016). The inversion assumes plane shear waves ur,t with complex wave number k* (k*=k+ik), angular frequency ω, and amplitude u0 as solutions of the following wave equation:

ur,t=u0eik*rωt.(1)

k-MDEV first decomposed the complex-valued wavefields ur,t obtained from the Hilbert transformation into eight complex wavefields ud 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 ud:

kd=argud,(2a)
kd=udud.(2b)

The final maps of SWS were obtained by averaging over directions:

SWS=ωdwdkdw,withw=ud2.(3a)

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 (Cuiwei et al., 1995) 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
www.frontiersin.org

FIGURE 2. Representative data of B-mode, color flow imaging, and THE. (A) B-mode image overlaid with color flow of the middle cerebral artery (blue-to-red color scale). (B) SWS map. In (A) and (B), the regions of interest are demarcated by white lines: temporal lobe parenchyma (TL), basal cistern (BC), and brain stem (BS) global region (dashed line). (C) Time curves of SWS (blue) and CBF (gray spectrum and red envelope) together with ECG (green).

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 (Krakauskaite et al., 2018), and SWS was slightly higher than that in previous cerebral THE studies due to the higher excitation frequency (Kreft et al., 2020b).

FIGURE 3
www.frontiersin.org

FIGURE 3. Group-averaged CBF and SWS values over normalized cardiac cycles. Shear wave speed (SWS, blue) with 95% confidence interval (blue area), cerebral blood flow (CBF, red) with the 95% confidence interval (red area), as well as ECG (green) for (A) global brain, (B) temporal lobe, (C) basal cisterns, and (D) brain stem regions. Cranial arterial systole and diastole are demarcated in (A). In (B), the peak systolic SWS and end diastolic SWS are demarcated by vertical dashed lines. τR indicates the time delay between ECG and SWS, while τΔ is the time delay between the SWS and Doppler peak. Horizontal dashed lines refer to the mean SWS values.

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
www.frontiersin.org

FIGURE 4. Group statistical plots of SWS and delay times. (A) SWS in peak systole and end diastole for all four regions of interest. In all regions, SWS decreased from systole to diastole. (B) Time delays between systolic ECG R-wave and peak SWS (τR), as well as those between the peak SWS and peak CBF (τΔ) for all examined regions of interest. **p < 0.01 and *p < 0.05.

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 (Hirsch et al., 2012; Hirsch et al., 2014a; Hirsch et al., 2014b; Parker, 2014; Parker, 2015; Arani et al., 2018; Marticorena Garcia et al., 2018; Forouhandehpour et al., 2021; Solamen et al., 2021). Specifically, brain stiffness has been shown by MRE and THE to increase with cerebral artery dilation due to hypercapnia (Hetzer et al., 2019; Kreft et al., 2020a) or the Valsalva maneuver (Tzschatzsch et al., 2018; Herthum et al., 2021). As shown by MRE, cerebral stiffness of deep gray matter correlates with perfusion pressure rather than cerebral blood flow alone (Hetzer et al., 2018). 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 (Ipek-Ugay et al., 2016; Kreft et al., 2020b) 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 (Hirsch et al., 2014a; Meyer et al., 2022).

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 (Kreft et al., 2022), whereas mixed results have been reported based on MRE (Kolipaka et al., 2018).

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 (Kreft et al., 2020b). 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 (Lauriola et al., 2018; Liu et al., 2020; Kreft et al., 2022).

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.

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

Adams, A. L., Kuijf, H. J., Viergever, M. A., Luijten, P. R., and Zwanenburg, J. J. M. (2019). Quantifying cardiac-induced brain tissue expansion using DENSE. NMR Biomed. 32 (2), e4050. doi:10.1002/nbm.4050

PubMed Abstract | CrossRef Full Text | Google Scholar

Arani, A., Min, H. K., Fattahi, N., Wetjen, N. M., Trzasko, J. D., Manduca, A., 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. doi:10.1002/mrm.26738

PubMed Abstract | CrossRef Full Text | Google Scholar

Cosgrove, D., Barr, R., Bojunga, J., Cantisani, V., Chammas, M. C., Dighe, M., et al. (2017). WFUMB guidelines and recommendations on the clinical use of ultrasound elastography: part 4. Thyroid. Ultrasound Med. Biol. 43 (1), 4–26. doi:10.1016/j.ultrasmedbio.2016.06.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Cuiwei, L., Chongxun, Z., and Changfeng, T. (1995). Detection of ECG characteristic points using wavelet transforms. IEEE Trans. Biomed. Eng. 42 (1), 21–28. doi:10.1109/10.362922

PubMed Abstract | CrossRef Full Text | Google Scholar

Ertl, M., Raasch, N., Hammel, G., Harter, K., and Lang, C. (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. doi:10.1016/j.ultrasmedbio.2017.08.1885

PubMed Abstract | CrossRef Full Text | Google Scholar

Ferraioli, G., Filice, C., Castera, L., Choi, B. I., Sporea, I., Wilson, S. R., et al. (2015). WFUMB guidelines and recommendations for clinical use of ultrasound elastography: part 3: liver. Ultrasound Med. Biol. 41 (5), 1161–1179. doi:10.1016/j.ultrasmedbio.2015.03.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Forouhandehpour, R., Bernier, M., Gilbert, G., Butler, R., Whittingstall, K., and Van Houten, E. (2021). Cerebral stiffness changes during visual stimulation: differential physiological mechanisms characterized by opposing mechanical effects. Neuroimage Rep. 1 (2), 100014. doi:10.1016/j.ynirp.2021.100014

CrossRef Full Text | Google Scholar

Ge, G. R., Song, W., Nedergaard, M., Rolland, J. P., and Parker, K. J. (2022). Theory of sleep/wake cycles affecting brain elastography. Phys. Med. Biol. 67 (22), 225013. doi:10.1088/1361-6560/ac9e40

CrossRef Full Text | Google Scholar

Herthum, H., Hetzer, S., Kreft, B., Tzschätzsch, H., Shahryari, M., Meyer, T., et al. (2022). Cerebral tomoelastography based on multifrequency MR elastography in two and three dimensions. Front. Bioeng. Biotechnol. 10, 1056131. doi:10.3389/fbioe.2022.1056131

PubMed Abstract | CrossRef Full Text | Google Scholar

Herthum, H., Shahryari, M., Tzschatzsch, H., Schrank, F., Warmuth, C., Gorner, S., 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. doi:10.3389/fbioe.2021.666456

PubMed Abstract | CrossRef Full Text | Google Scholar

Hetzer, S., Birr, P., Fehlner, A., Hirsch, S., Dittmann, F., Barnhill, E., et al. (2018). Perfusion alters stiffness of deep gray matter. J. Cereb. Blood Flow. Metab. 38 (1), 116–125. doi:10.1177/0271678x17691530

PubMed Abstract | CrossRef Full Text | Google Scholar

Hetzer, S., Dittmann, F., Bormann, K., Hirsch, S., Lipp, A., Wang, D. J., et al. (2019). Hypercapnia increases brain viscoelasticity. J. Cereb. Blood Flow. Metab. 39 (12), 2445–2455. doi:10.1177/0271678x18799241

PubMed Abstract | CrossRef Full Text | Google Scholar

Hirsch, S., Guo, J., Reiter, R., Papazoglou, S., Kroencke, T., Braun, J., 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. doi:10.1002/mrm.24674

PubMed Abstract | CrossRef Full Text | Google Scholar

Hirsch, S., Guo, J., Reiter, R., Schott, E., Buning, C., Somasundaram, R., et al. (2014b). Towards compression-sensitive magnetic resonance elastography of the liver: sensitivity of harmonic volumetric strain to portal hypertension. J. Magn. Reson Imaging 39 (2), 298–306. doi:10.1002/jmri.24165

PubMed Abstract | CrossRef Full Text | Google Scholar

Hirsch, S., Klatt, D., Freimann, F., Scheel, M., Braun, J., and Sack, I. (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. doi:10.1002/mrm.24499

PubMed Abstract | CrossRef Full Text | Google Scholar

Hiscox, L. V., Johnson, C. L., Barnhill, E., McGarry, M. D., Huston, J., van Beek, E. J., et al. (2016). Magnetic resonance elastography (MRE) of the human brain: technique, findings and clinical applications. Phys. Med. Biol. 61 (24), R401–R437. doi:10.1088/0031-9155/61/24/r401

PubMed Abstract | CrossRef Full Text | Google Scholar

Ipek-Ugay, S., Tzschatzsch, H., Hudert, C., Marticorena Garcia, S. R., Fischer, T., Braun, J., et al. (2016). Time harmonic elastography reveals sensitivity of liver stiffness to water ingestion. Ultrasound Med. Biol. 42 (6), 1289–1294. doi:10.1016/j.ultrasmedbio.2015.12.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Kolipaka, A., Wassenaar, P. A., Cha, S., Marashdeh, W. M., Mo, X., Kalra, P., et al. (2018). Magnetic resonance elastography to estimate brain stiffness: measurement reproducibility and its estimate in pseudotumor cerebri patients. Clin. Imaging 51, 114–122. doi:10.1016/j.clinimag.2018.02.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Krakauskaite, S., Thibeault, C., LaVangie, J., Scheidt, M., Martinez, L., Seth-Hunter, D., et al. (2018). Normative ranges of transcranial Doppler metrics. Acta Neurochir. Suppl. 126, 269–273. doi:10.1007/978-3-319-65798-1_53

PubMed Abstract | CrossRef Full Text | Google Scholar

Kreft, B., Bergs, J., Shahryari, M., Danyel, L. A., Hetzer, S., Braun, J., et al. (2020b). Cerebral ultrasound time-harmonic elastography reveals softening of the human brain due to dehydration. Front. Physiol. 11, 616984. doi:10.3389/fphys.2020.616984

PubMed Abstract | CrossRef Full Text | Google Scholar

Kreft, B., Tzschatzsch, H., Schrank, F., Bergs, J., Streitberger, K. J., Waldchen, S., et al. (2020a). Time-resolved response of cerebral stiffness to hypercapnia in humans. Ultrasound Med. Biol. 46 (4), 936–943. doi:10.1016/j.ultrasmedbio.2019.12.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Kreft, B., Tzschatzsch, H., Shahryari, M., Haffner, P., Braun, J., Sack, I., et al. (2022). Noninvasive detection of intracranial hypertension by novel ultrasound time-harmonic elastography. Invest. Radiol. 57 (2), 77–84. doi:10.1097/rli.0000000000000817

PubMed Abstract | CrossRef Full Text | Google Scholar

Lauriola, M., Mangiacotti, A., D'Onofrio, G., Cascavilla, L., Paris, F., Paroni, G., et al. (2018). Neurocognitive disorders and dehydration in older patients: clinical experience supports the hydromolecular hypothesis of dementia. Nutrients 10 (5), 562. doi:10.3390/nu10050562

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, H., Wang, D., Leng, X., Zheng, D., Chen, F., Wong, L. K. S., et al. (2020). State-of-the-Art computational models of circle of willis with physiological applications: A review. IEEE Access 8, 156261–156273. doi:10.1109/access.2020.3007737

CrossRef Full Text | Google Scholar

Manduca, A., Bayly, P. J., Ehman, R. L., Kolipaka, A., Royston, T. J., Sack, I., et al. (2021). MR elastography: principles, guidelines, and terminology. Magn. Reson Med. 85 (5), 2377–2390. doi:10.1002/mrm.28627

PubMed Abstract | CrossRef Full Text | Google Scholar

Marticorena Garcia, S. R., Grossmann, M., Lang, S. T., Tzschatzsch, H., Dittmann, F., Hamm, B., 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. doi:10.1002/mrm.26892

PubMed Abstract | CrossRef Full Text | Google Scholar

Meyer, T., Tzschatzsch, H., Wellge, B., Sack, I., Kroncke, T., and Martl, A. (2022). Valsalva maneuver decreases liver and spleen stiffness measured by time-harmonic ultrasound elastography. Front. Bioeng. Biotechnol. 10, 886363. doi:10.3389/fbioe.2022.886363

PubMed Abstract | CrossRef Full Text | Google Scholar

Nanjappa, M., and Kolipaka, A. (2021). Magnetic resonance elastography of the brain. Magnetic Reson. Imaging Clin. N. Am. 29 (4), 617–630. doi:10.1016/j.mric.2021.06.011

CrossRef Full Text | Google Scholar

Ormachea, J., and Parker, K. J. (2020). Elastography imaging: the 30 year perspective. Phys. Med. Biol. 65 (24), 1361–6560. doi:10.1088/1361-6560/abca00

PubMed Abstract | CrossRef Full Text | Google Scholar

Parker, K. J. (2014). A microchannel flow model for soft tissue elasticity. Phys. Med. Biol. 59 (15), 4443–4457. doi:10.1088/0031-9155/59/15/4443

PubMed Abstract | CrossRef Full Text | Google Scholar

Parker, K. J. (2015). Experimental evaluations of the microchannel flow model. Phys. Med. Biol. 60 (11), 4227–4242. doi:10.1088/0031-9155/60/11/4227

PubMed Abstract | CrossRef Full Text | Google Scholar

Poul, S. S., Ormachea, J., Hollenbach, S. J., and Parker, K. J. (2020). Validations of the microchannel flow model for characterizing vascularized tissues. Fluids 5 (4), 228. doi:10.3390/fluids5040228

PubMed Abstract | CrossRef Full Text | Google Scholar

Robba, C., Cardim, D., Sekhon, M., Budohoski, K., and Czosnyka, M. (2018). Transcranial Doppler: a stethoscope for the brain-neurocritical care use. J. Neurosci. Res. 96 (4), 720–730. doi:10.1002/jnr.24148

PubMed Abstract | CrossRef Full Text | Google Scholar

Rouch, L., Cestac, P., Sallerin, B., Andrieu, S., Bailly, H., Beunardeau, M., et al. (2018). Pulse wave velocity is associated with greater risk of dementia in mild cognitive impairment patients. Hypertension 72 (5), 1109–1116. doi:10.1161/hypertensionaha.118.11443

PubMed Abstract | CrossRef Full Text | Google Scholar

Sack, I. (2022). Magnetic resonance elastography from fundamental soft-tissue mechanics to diagnostic imaging. Nat. Rev. Phys. 5, 25–42. doi:10.1038/s42254-022-00543-2

CrossRef Full Text | Google Scholar

Schrank, F., Warmuth, C., Tzschatzsch, H., Kreft, B., Hirsch, S., Braun, J., 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. doi:10.1177/0271678x19850936

PubMed Abstract | CrossRef Full Text | Google Scholar

Shiina, T., Nightingale, K. R., Palmeri, M. L., Hall, T. J., Bamber, J. C., Barr, R. 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. doi:10.1016/j.ultrasmedbio.2015.03.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Solamen, L. M., McGarry, M. D. J., Fried, J., Weaver, J. B., Lollis, S. S., and Paulsen, K. 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. doi:10.1016/j.acra.2020.03.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Sunkaria, R. K., Saxena, S. C., Kumar, V., and Singhal, A. M. (2010). Wavelet based R-peak detection for heart rate variability studies. J. Med. Eng. Technol. 34 (2), 108–115. doi:10.3109/03091900903281215

PubMed Abstract | CrossRef Full Text | Google Scholar

Tzschatzsch, H., Kreft, B., Schrank, F., Bergs, J., Braun, J., and Sack, I. (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. doi:10.1038/s41598-018-36191-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Tzschatzsch, H., Nguyen Trong, M., Scheuermann, T., Ipek-Ugay, S., Fischer, T., Schultz, M., et al. (2016). Two-dimensional time-harmonic elastography of the human liver and spleen. Ultrasound Med. Biol. 42 (11), 2562–2571. doi:10.1016/j.ultrasmedbio.2016.07.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Wagshul, M., Eide, P., and Madsen, J. (2011). The pulsating brain: A review of experimental and clinical studies of intracranial pulsatility. Fluids Barriers CNS 8, 5. doi:10.1186/2045-8118-8-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Wolters, F. J., Zonneveld, H. I., Hofman, A., van der Lugt, A., Koudstaal, P. J., Vernooij, M. W., et al. (2017). Cerebral perfusion and the risk of dementia A population-based study. Circulation 136 (8), 719–728. doi:10.1161/circulationaha.117.027448

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, L., Liu, C. Y., Smith, R. X., Jog, M., Langham, M., Krasileva, K., et al. (2016). Assessing intracranial vascular compliance using dynamic arterial spin labeling. Neuroimage 124, 433–441. doi:10.1016/j.neuroimage.2015.09.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Yin, Z., Romano, A. J., Manduca, A., Ehman, R. L., and Huston, J. (2018). Stiffness and beyond: what MR elastography can tell us about brain structure and function under physiologic and pathologic conditions. Top. Magn. Reson Imaging 27 (5), 305–318. doi:10.1097/rmr.0000000000000178

PubMed Abstract | CrossRef Full Text | Google Scholar

Zamir, M., Moir, M. E., Klassen, S. A., Balestrini, C. S., and Shoemaker, J. K. (2018). Cerebrovascular compliance within the rigid confines of the skull. Front. Physiology 9, 940. doi:10.3389/fphys.2018.00940

PubMed Abstract | CrossRef Full Text | Google Scholar

Zwanenburg, J. J. M., and van Osch, M. J. P. (2017). Targeting cerebral small vessel disease with MRI. Stroke 48 (11), 3175–3182. doi:10.1161/strokeaha.117.016996

PubMed Abstract | CrossRef Full Text | Google Scholar

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.

Edited by:

Salavat Aglyamov, University of Houston, United States

Reviewed by:

Haipeng Liu, Coventry University, United Kingdom
Manmohan Singh, University of Houston, United States

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

Download