ORIGINAL RESEARCH article

Front. Neurosci., 24 September 2020

Sec. Autonomic Neuroscience

Volume 14 - 2020 | https://doi.org/10.3389/fnins.2020.547433

Instantaneous Cardiac Baroreflex Sensitivity: xBRS Method Quantifies Heart Rate Blood Pressure Variability Ratio at Rest and During Slow Breathing

  • 1. Department of Physics, Humboldt-Universität zu Berlin, Berlin, Germany

  • 2. Institute of Biomedical Engineering, Technische Universität Dresden, Dresden, Germany

  • 3. Division of Obstetrics, Universitätsklinikum Leipzig, Leipzig, Germany

  • 4. Potsdam Institute for Climate Impact Research, Potsdam, Germany

  • 5. Department of Human and Animal Physiology, Saratov State University, Saratov, Russia

Abstract

Spontaneous baroreflex sensitivity (BRS) is a widely used tool for the quantification of the cardiovascular regulation. Numerous groups use the xBRS method, which calculates the cross-correlation between the systolic beat-to-beat blood pressure and the R-R interval (resampled at 1 Hz) in a 10 s sliding window, with 0–5 s delays for the interval. The delay with the highest correlation is selected and, if significant, the quotient of the standard deviations of the R-R intervals and the systolic blood pressures is recorded as the corresponding xBRS value. In this paper we test the hypothesis that the xBRS method quantifies the causal interactions of spontaneous BRS from non-invasive measurements at rest. We use the term spontaneous BRS in the sense of the sensitivity curve is calculated from non-interventional, i.e., spontaneous, baroreceptor activity. This study includes retrospective analysis of 1828 measurements containing ECG as well as continues blood pressure under resting conditions. Our results show a high correlation between the heart rate – systolic blood pressure variability (HRV/BPV) quotient and the xBRS (r = 0.94, p < 0.001). For a deeper understanding we conducted two surrogate analyses by substituting the systolic blood pressure by its reversed time series. These showed that the xBRS method was not able to quantify causal relationships between the two signals. It was not possible to distinguish between random and baroreflex controlled sequences. It appears xBRS rather determines the HRV/BPV quotient. We conclude that the xBRS method has a potentially large bias in characterizing the capacity of the arterial baroreflex under resting conditions. During slow breathing, estimates for xBRS are significantly increased, which clearly shows that measurements at rest only involve limited baroreflex activity, but does neither challenge, nor show the full range of the arterial baroreflex regulatory capacity. We show that xBRS is exclusively dominated by the heart rate to systolic blood pressure ratio (r = 0.965, p < 0.001). Further investigations should focus on additional autonomous testing procedures such as slow breathing or orthostatic testing to provide a basis for a non-invasive evaluation of baroreflex sensitivity.

Introduction

The baroreflex is an important component of cardiovascular regulation to maintain homeostasis. The idea of spontaneous baroreflex sensitivity is to estimate the concomitant effect of respiration on heart period and blood pressure based on non-invasive and non-pharmacological driven measurements. Originally, BRS has been assessed by the Oxford method, based on analysis of heart rate response to drug-induced blood pressure variations (). This method is still the gold standard for assessing baroreflex control. This method has, however, not found wide application in clinical practice due to its laboriousness. It is invasive and requires the administration of vasoactive substances, which is potentially unsafe and costly. Furthermore, undesirable effects of medications on the state of the ANS cannot be excluded. To overcome these drawbacks, numerous methods for non-invasive assessment of BRS have been developed, based on the analysis of spontaneous fluctuations in systolic blood pressure (SAP) and the RR interval (RR). Despite the fact that the idea of using spontaneous heart rate and pressure variations to assess baroreflex may seems desirable, several problems are emerging: In addition to arterial baroreflex itself there are many other sources for pressure and heart rate variations and it is almost impossible to discern baroreflex-driven variations from this mixture. There is no known possibility to isolate specific sets of stimuli and their corresponding reactions. As blood pressure fluctuations during rest in equilibrium are tiny, the effects contributed to the baroreflex seem to be extremely challenging () and drastically reducing the signal-to-noise ratio. The synchronous fluctuations of heart rate with respiration, known as respiratory sinus arrhythmia (RSA), are a consequence of the rapid fluctuations of parasympathetic nerve activity toward the sinus node (). The origin of RSA is known to have various mechanisms (), including arterial baroreflex and cardiopulmonary baroreceptor responses due to fluctuations of cardiac stroke volume, a direct influence of medullary respiratory neurons on the vagal motor nucleus, and pulmonary stretch receptor response to lung inflation.

Fluctuations in blood pressure and heart period can be of clinical importance as risk markers for cardiovascular morbidity and mortality (; ; ). However, show that spontaneous baroreflex indices do not clearly reflect arterial baroreflex gain. They mainly quantify vagal-mediated heart period oscillations induced by cardiac output fluctuations, and do not reflect barosensory vessel distensibility. Without a clear and consistent relationship the baroreflex gain itself, one can only conclude that spontaneous baroreflex sensitivity cannot be used as proxy for baroreflex gain. Nevertheless, the quantification of this reflex is of great relevance for understanding the cardiovascular system and for risk stratification (; ; ). Recently () we were able to show that the spontaneous sensitivity of the arterial baroreflex (BRS) under resting conditions cannot be estimated by the sequence method (SME), which only quantifies the quotient of heart rate and systolic blood pressure variability. In this paper we test whether the xBRS method (; ) is suitable to quantify the baroreflex sensitivity from non-invasive, non-interventional measurements under resting conditions. Therefore, two surrogate analyses were performed in which, due to the design, no causal relationships between blood pressure and heart rate signal can be present. Furthermore, SME and xBRS were calculated from data collected not only under resting conditions but also under controlled breathing.

Data

To allow comparison to the results in , we reanalyzed the same data from 5 different studies in obstetrics, genetics, cardiology and heart surgery (; ; ; ; ) in a similar manner. Demographic data of all sub-studies were given in : “All patients gave written, informed consent, and all studies were approved by the respective local ethics committees. From obstetrics () 915 measurements of 304 pregnant women were included (mean age 28.4 ± 5.4 years). The data contain 398 recordings of healthy women, 120 from patients with chronic hypertension, 38 from gestational hypertension, 152 from women who later developed pre-eclampsia, 88 from pre-existing hypertension with pre-eclampsia, 12 with other hypertensive disease and 78 from women with intrauterine growth restriction. From genetics () we considered measurements from 367 subjects with an age of 10 to 88 years (45.0 ± 16.3 years), 157 were male (43%). From cardiology () we used the measurements from 75 patients with chronic cardiac diseases referred for primary preventive implantable cardioverter-defibrillator implantation following Multicenter Automatic Defibrillator Implantation Trial study criteria, mean age 70.9 ± 10.1 years, body mass index 27.0 ± 3.5. From 302 measurements from patients before and after aortic (AV) or mitral valve (MV) surgery were included for analysis. The mean age of the AV patients and MV patients was 62 ± 13 years and 59 ± 2 years, respectively. From 169 measurements from 58 consecutive patients undergoing either trans-catheter aortic valve implantation (TAVI) or surgical aortic valve replacement (SAVR) with the heart-lung machine and being in stable sinus rhythm were enrolled. Thirty four of them underwent SAVR and 24 of them TAVI, 28 males, mean age 64.6 ± 13.8 in the SAVR group and 80.5 ± 7.3 in TAVI.

All measurements of the considered studies were performed under supine resting position for 30 min using the Task Force Monitor (CNSystems, Graz) or the PortaPres device (Finapres Medical Systems, Enschede). In total we gathered 1,828 time series containing the beat-to-beat values of heart rate (HR) as well as systolic blood pressure (SBP). Exclusion criteria were atrial fibrillation, pacemaker activity, technical artifacts, as well as ectopy time greater than 10%, reducing the number of time series to 1,576 – careful visual inspection for further technical and physiological artifacts reduced the subjects for reanalysis to 1,439.

In addition to this comprehensive data set of rest measurements, recordings with certain autonomous testing procedures were analyzed in this paper. We used 245 measurements of 44 healthy pregnant women, mean age 30 ± 5.4 years, from the Fetal Autonomic Cardiovascular rEgulation (FACE) study which is currently in progress at the University of Leipzig Medical Center in cooperation with the TU Dresden and the Humboldt-Universität zu Berlin. The study was approved by the committee of ethics of the University of Leipzig Medical Center (357/17-ek). One aim of this study was to characterize the reaction of fetal autonomic regulation to maternal paced breathing based on a context dependent biosignal analysis. The measurements were performed under supine resting position using the PortaPres device (Finapres Medical Systems, Enschede), the fetal ECG signal was recorded from the abdomen of the pregnant woman. Our measurement protocol included 10 min of measurement at rest in supine position, 5 min of paced slow respiration (period 7.5 s – 8 respiration cycles per min), and 5 min of fast respiration (period 3 s – 20 cycles per min). Between both paced respiration conditions was one break of 5 min at rest. Exclusion criteria were rhythm disturbances (many of ventricular or supraventricular ectopic beats), technical artifacts, and incomplete study protocols. The data underwent a careful visual inspection for further artifacts by experts which reduced the analyzed data set to 184 records.

Materials and Methods

Originally, BRS was measured by injecting vasoconstrictive agents to raise blood pressure, i.e., quantifying the reflex-like increased beat-to-beat intervals in the electrocardiogram (ECG). Later, attempts have been made to determine baroreflex sensitivity non-invasively, most often using spontaneous heart rate variability (HRV) and blood pressure variability (BPV) obtained from continuous finger pressure measurement (; ). The underlying hypothesis is that there is always some spontaneous variability in blood pressure that should allow an estimate of BRS. In this paper, the BRS estimation method under consideration was the xBRS method (; ).

xBRS is a time domain method designed for estimation of BRS from non-invasively obtained SBP and beat-to-beat-interval (BBI) data. The original time series is resampled to obtain evenly sampled data at 1 Hz. Instead of sequences, xBRS uses sliding windows with a fixed length of 10 seconds. For every SBP window, there should exist at least one corresponding BBI window (in a lag range of 0–5) with positive and significant cross correlation (p < 0.05, two sided test for zero correlation). The BBI window with highest cross correlation is then selected. xBRS has two advantages: (a) it gives more valid “sequences” for the same time series compared to SME and (b) it allows to observe slow regulation circuits that are presumably arising from sympathetic control. Usually, the resulting number of valid “sequences” in the time series is large enough to enable the calculation of an “instant” measure for xBRS. The xBRS method was developed by and initially it set strict constrains on p-value of the cross correlation (p < 0.01). In the last revision (), the threshold for the p-value was increased to 0.05, it doubles the percentage of necessary valid windows. Once sequences are selected, the ratio between the standard deviations of BBI and of SBP is used as a measure of baroreflex sensitivity. The reasoning for this decision was that the measure based on ratios of standard deviations does not differ significantly from the slopes of the regressions but is computationally more efficient. In our article we implemented the last revision of the xBRS method; the geometric mean of the instant xBRS values was used as BRS estimate per study period reflecting the approximately log-normal distribution of xBRS values ().

Furthermore, different parameters from the time domain were calculated to quantify short term HRV and BPV in this large data set (, cf. Table 1). The RMSSDRATIO defined as RMSSDHRV divided by RMSSDSBP showed the highest correlation to the sequence method for BRS estimation there ().

TABLE 1

Mean ± SDR to xBRS
xBRS [ms/mmHg]8.4 ± 5.31
xBRSS1 [ms/mmHg]8.4 ± 5.20.99
xBRSS2 [ms/mmHg]7.9 ± 4.70.97
meanNN [ms]761 ± 1460.54
meanBP [mmHg]129 ± 22.4−0.2
sdNN [ms] 42.5 ± 18.60.67
sdBP [mmHg]8.1 ± 2.9−0.25
RMSSD [ms] 25.6 ± 14.70.84
RMSSDSBP [mmHg]2.9 ± 1−0.21
RMSSDRATIO [ms/mmHg]9.5 ± 6.10.94

Basic characteristics of considered parameters in the final data set (Mean ± SD: mean value ± standard deviation) and the correlation coefficient r to xBRS (R to xBRS, p < 0.001 for all coefficients).

xBRSS1, surrogate 1 of xBRS; xBRSS2, surrogate 2 of xBRS; meanNN, mean beat-to-beat-interval of HR; meanBP, mean blood pressure of SBP; sdNN, standard deviation of HR; sdBP, standard deviation of SBP; RMSSD, root mean square of successive differences of HR; RMSSDSBP, root mean square of successive differences of SBP; RMSSDRATIO, Ratio of rmssd and rmssd_sbp.

Analogous to () we performed two surrogate analyses in order to ascertain whether xBRS can quantify causal relationships between heart rate and blood pressure and thusly spontaneous BRS:

  • (a)

    Systolic blood pressure time-series were analyzed in reversed order, i.e., the first blood pressure value is now the last, the second now the second last etc. This results in surrogate data with the same distribution as the original data since the values of each point are the same, just in a different time position. However, any causal relationship between heart rate and blood pressure has been removed by this procedure (cf. Figure 1).

    FIGURE 1

    .

  • (b)

    Beat-to-beat-intervals were shuffled using the IAFFT approach (). By applying this procedure causal relationship between heart rate and blood pressure has, again, been removed (xBRSS2).

Both surrogate tests are used to test the following hypothesis: “xBRS does quantify causal relationships between heart rate and blood pressure.”

Statistical analysis was performed using IBM SPSS Statistics version 24. To quantify significant relations between the parameters used in this study we applied Pearson’s correlation as a measure of the linear relationship between two continuous random variables. This measure does not assume normality, but assumes finite variance as well as covariance - properties which were assumed for our data sets.

Results

Table 1 shows the results of the correlation analysis between xBRS and further HRV and BPV parameters in the final data set (n = 1,439). The highest correlation coefficient r was found between xBRS and short-term variability parameters RMSSD (r = 0.84, p < 0.001). A more precise estimate yields the highest correlation coefficient for RMSSDRATIO (r = 0.94, p < 0.001, cf. Figure 2).

FIGURE 2

Surrogate analysis (a) between xBRSS1 and RMSSDHRV showed high correlation coefficient (r = 0.84, p < 0.001), i.e., reversing one time series does not affect the results of xBRS. Moreover, the correlation coefficient between xBRSS1 and RMSSDRATIO was equally high with r = 0.94, p < 0.001. Finally, the correlation coefficient between original xBRS and surrogate xBRSS1 was extremely high with r = 0.99, p < 0.001 (cf. Table 1 and Figure 3).

FIGURE 3

In the surrogate analysis (b) the correlation coefficient between xBRSS2 and RMSSDHRV was r = 0.84, p < 0.001. Moreover, the correlation coefficient between xBRSS2 and RMSSDRATIO was equally high with r = 0.93, p < 0.001. Finally, the correlation coefficient between the original xBRS and surrogate xBRSS2 again was extremely high with r = 0.97, p < 0.001 (cf. Table 1).

To overcome concerns for bias due to the presence of multiple measurements per subjects, we repeated our procedures on a reduced data set with only the first measurement per subject (n = 733) and got similar correlations (xBRS vs. RMSSD: 0.83, xBRS vs. RMSSDRATIO: 0.93, xBRSS1 vs. RMSSD: 0.83, xBRSS1 vs. RMSSDRATIO: 0.93, xBRS vs. xBRSS1: 0.98). We conclude that the comparison of time series with different time bases has no influence on the results of xBRS. xBRS as an estimate for the spontaneous BRS shows a potentially large methodological bias. This contradicts the hypothesis that xBRS at rest quantifies causal relationships between heart rate and blood pressure.

In order to investigate whether any causal relationship could be quantified by xBRS under slow breathing conditions, we analyzed the data from the FACE study. Figure 4 shows the xBRS and the SME values under rest conditions, under rapid breathing as well as under slow breathing. There is a, clearly significant, increase in xBRS and SME with slow breathing (p < 0.001), showing that measurements at rest only involves certain range baroreflex activation, but not the full capacity of the arterial baroreflex. However, these estimates are, again, exclusively dominated by the heart rate - systolic blood pressure ratio (r = 0.965, p < 0.001). Moreover, after performing the surrogate analyzes described above, the xBRS estimates do not change between original and surrogate data: The correlation coefficient at rest: xBRS to xBRSS1r = 0.98, to xBRSS2r = 0.96; 3 s (fast) breathing: xBRS to xBRSS1r = 0.98, to xBRSS2r = 0.94; slow respiration: xBRS to xBRSS1r = 0.99, to xBRSS2r = 0.97.

FIGURE 4

In summary, we found that, even under controlled breathing conditions, no causal relationship between beat-to-beat intervals and the systolic blood pressures could be quantified by xBRS.

Nevertheless, Figure 5 shows that significantly more, as well as more consistent, results are obtained under controlled slow breathing compared to normal resting measurements. While under resting conditions only 31% of all calculated correlations are significant for this example and are included in the calculation of the xBRS, during slow breathing, the percentage of valid xBRS windows increases to 98%. Furthermore, the mean correlation between heart rate and systolic blood pressure windows is significantly higher during slow breathing. Taken together, these values support the hypothesis that the xBRS values at slow breathing are more consistent with the BRS. A higher BRS implies a higher regulatory capacity, meaning large blood pressure fluctuations can be balanced well by BBI changes which in turn supports maintenance of homeostasis. However, a low baroreflex sensitivity can lead to deviations from homeostasis and thus to events such as a hypertensive crises or fainting spells due to low blood pressure.

FIGURE 5

Discussion

In this paper we test the hypothesis that the xBRS method quantifies the causal interactions of spontaneous BRS from non-invasive, non-interventional measurements at rest. For the surrogate analysis we substituted the systolic blood pressure by its reversed time series and thusly removed their causal relationship. Our analysis showed that xBRS remains unchanged. Therefore, we conclude that xBRS is not able to estimate causal relationship between heart rate and systolic blood pressure. Our results of the surrogate analysis show, that xBRS can be mostly explained by the short-term HRV, quantified by the RMSSD. Minor deviations of this univariate model are adequately explained by the simple bivariate model RMSSDRATIO which is also based on non-causal interactions.

In contrast, was able to show that low values of pharmacologically determined baroreflex sensitivity (pBRS < 3 ms per mmHg) carry a significant risk of cardiac mortality after myocardial infarction. This indicates a difference between the xBRS estimation at rest and an invasive pBRS, which not only refers to respiration-induced fluctuations, but also includes carotid extensibility (). We suspect that the estimation of baroreflex sensitivity is unreliable in cases of relatively shallow breathing. In these cases, there is also only a small respiratory induced blood pressure variation, which leads to only small baroreceptor activation and thus to a low HRV. Even then, in cases where any BRS estimation would result in a spurious low result, the baroreflex could still be fully functional and its sensitivity in the normal range. The BRS could just be impossible to estimate using the currently dominant protocol, i.e., relaxed respiration in resting supine position (). and show that baroreflex sensitivity varies greatly from patient to patient and that different mechanical (neck chamber) and pharmacological techniques for measuring baroreflex sensitivity are likely to measure different aspects of baroreflex function. This contradicts the idea of spontaneous baroreflex sensitivity in the sense that the sensitivity should not be affected by the way it is being measured. Recently developed sophisticated methods for BRS estimation disentangle the effects of respiration from heart period and blood pressure (; ; ), however, all methods fail to be reliable estimates of BRS. An increase in SME values during slow breathing lead to systematically review the baroreflex function. His hypothesis was that controlled slow breathing, which causes higher blood pressure fluctuations, increases cardiovagal baroreflex gain in young healthy subjects. Baroreflex enhancement was investigated using both the classical pBRS and the non-invasive SME method. Compared to breathing at rest, slow breathing was associated with a significant increase in the SME index, while the pBRS remained unchanged. The SME values for slow breathing are higher than the pBRS values in the study, which could be a result of overestimation or a systematic error in the pBRS determination. However, showed that pBRS values can be predicted from non-pharmacological indices acquired during slow breathing. From both studies we derive the opinion that autonomic testing should allow a reliable, non-invasive, non-pharmacological driven quantification of the baroreflex gain. This will require large scale medical studies, where the BRS is measured invasively according to the state of the art [modified Oxford method ()] and additional runs of autonomous tests are performed.

To investigate whether a causal relationship between heart rate and blood pressure exist under controlled breathing conditions, we analyzed the data from the FACE study. We found a clearly significant increase in xBRS and SME under slow breathing conditions compared to rest or fast breathing, showing that the latter ones involve only limited baroreflex activity. However, for all conditions these estimates are exclusively dominated by the heart rate – systolic blood pressure ratio. Moreover, after performing the surrogate analyzes, the xBRS does not vary between original and surrogate data. Thus, even under controlled breathing conditions, no causal relationship between beat-to-beat intervals and the systolic blood pressures could be found.

In conclusion, we demonstrated for all short measurements, under resting conditions and controlled breathing, that RMSSDRATIO carries similar vagally mediated information as xBRS. However, we found, under controlled breathing, a potentially large methodological bias in xBRS and SME as estimates for the baroreflex sensitivity. During slow breathing estimates for SME and xBRS are significantly increasing, which clearly shows that measurements at rest are only accompanied by limited baroreflex activity, but not to the full capacity of the arterial baroreflex. Further investigations should focus on additional autonomic testing procedures (e.g., orthostatic test, carotid occlusion, neck suction) to provide a better empirical foundation of non-invasive assessment of the baroreflex sensitivity.

Statements

Data availability statement

The data that support the finding of this study are fully available to the authors but cannot be shared publicly due to risk of violating privacy. Requests to access results datasets should be directed to NW.

Ethics statement

The studies involving human participants were reviewed and approved by the University of Leipzig Medical Center (357/17-ek). The patients/participants provided their written informed consent to participate in this study.

Author contributions

NW, AG, KB, JKr, and JKu contributed conception and design of the study. NW, HM, and HS organized the database. AG performed the statistical analysis. NW wrote the first draft of the manuscript. All authors contributed to manuscript revision, read and approved the submitted version.

Funding

NW, AG, and KB were supported by the DFG grant WE 2834/11-1. HM and JW were supported by the DFG grant MA 2294/8-1. HS was supported by the DFG grant STE 952/12-1. JKu was supported by the project RF Government Grant 075-15-2019-1885. JKr and KB were supported by the BMBF grant 13GW0278E. We acknowledge further support by the DFG for the Open Access Publication Fund of Humboldt-Universität zu Berlin.

Acknowledgments

We thank Hendrik Bonnemeiner, Robert Bauernschmitt, Andreas Voss, Alexander Schirdewan, and Ludmila Sidorenko for long standing cooperation and providing the data.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

Summary

Keywords

barorecepter reflex sensitivity, heart rate variability, blood pressure variability, slow breathing, rest

Citation

Wessel N, Gapelyuk A, Weiß J, Schmidt M, Kraemer JF, Berg K, Malberg H, Stepan H and Kurths J (2020) Instantaneous Cardiac Baroreflex Sensitivity: xBRS Method Quantifies Heart Rate Blood Pressure Variability Ratio at Rest and During Slow Breathing. Front. Neurosci. 14:547433. doi: 10.3389/fnins.2020.547433

Received

31 March 2020

Accepted

04 September 2020

Published

24 September 2020

Volume

14 - 2020

Edited by

Andreas Voss, Institut für Innovative Gesundheitstechnologien (IGHT), Germany

Reviewed by

Lauro C. Vianna, University of Brasilia, Brazil; Claudia Lerma, Instituto Nacional de Cardiologia Ignacio Chavez, Mexico

Updates

Copyright

*Correspondence: Niels Wessel,

This article was submitted to Autonomic Neuroscience, a section of the journal Frontiers in Neuroscience

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics