Impact Factor 3.707 | CiteScore 5.1
More on impact ›

Original Research ARTICLE

Front. Neurosci., 18 April 2019 |

Baroreflex Sensitivity Measured by Pulse Photoplethysmography

  • 1Department of Biomedical Engineering, University of Connecticut, Storrs, CT, United States
  • 2Biomedical Signal Interpretation and Computational Simulation (BSICoS) Group, Aragón Institute of Engineering Research (I3A), IIS Aragón, University of Zaragoza, Zaragoza, Spain
  • 3Centro de Investigación Biomédica en Red en Bioingeniería, Biomateriales y Nanomedicina (CIBER-BBN), Madrid, Spain
  • 4Department of Mechanical Engineering, University College London, London, United Kingdom

Novel methods for assessing baroreflex sensitivity (BRS) using only pulse photoplethysmography (PPG) signals are presented. Proposed methods were evaluated with a data set containing electrocardiogram (ECG), blood pressure (BP), and PPG signals from 17 healthy subjects during a tilt table test. The methods are based on a surrogate of α index, which is defined as the power ratio of RR interval variability (RRV) and that of systolic arterial pressure series variability (SAPV). The proposed α index surrogates use pulse-to-pulse interval series variability (PPV) as a surrogate of RRV, and different morphological features of the PPG pulse which have been hypothesized to be related to BP, as series surrogates of SAPV. A time-frequency technique was used to assess BRS, taking into account the non-stationarity of the protocol. This technique identifies two time-varying frequency bands where RRV and SAPV (or their surrogates) are expected to be coupled: the low frequency (LF, inside 0.04–0.15 Hz range), and the high frequency (HF, inside 0.15–0.4 Hz range) bands. Furthermore, time-frequency coherence is used to identify the time intervals when the RRV and SAPV (or their surrogates) are coupled. Conventional α index based on RRV and SAPV was used as Gold Standard. Spearman correlation coefficients between conventional α index and its PPG-based surrogates were computed and the paired Wilcoxon statistical test was applied in order to assess whether the indices can find significant differences (p < 0.05) between different stages of the protocol. The highest correlations with the conventional α index were obtained by the α-index-surrogate based on PPV and pulse up-slope (PUS), with 0.74 for LF band, and 0.81 for HF band. Furthermore, this index found significant differences between rest stages and tilt stage in both LF and HF bands according to the paired Wilcoxon test, as the conventional α index also did. These results suggest that BRS changes induced by the tilt test can be assessed with high correlation by only a PPG signal using PPV as RRV surrogate, and PPG morphological features as SAPV surrogates, being PUS the most convenient SAPV surrogate among the studied ones.

1. Introduction

The baroreflex system plays an important role in regulating short-term fluctuations of arterial blood pressure (BP) (La Rovere et al., 2008; Robertson et al., 2012). Arterial baroreceptors (placed in the wall of the carotid sinuses and aortic arch) sense changes in BP and modulate efferent autonomic neural activity to the central nervous system accordingly. A rise in sensed BP leads to an increase of vagal neurons discharge and a decrease in the discharge of sympathetic neurons, resulting in decreased heart rate (HR), cardiac contractility and peripheral vascular resistance. On the contrary, decreased BP enhances sympathetic and inhibits vagal activity, leading to increased HR, cardiac contractility and peripheral vascular resistance.

Cardiovascular diseases are frequently associated to an impairment of baroreflex mechanisms, resulting in chronic adrenergic activation. Reduced baroreflex control of HR has been reported in coronary artery disease, heart failure, hypertension and myocardial infarction (La Rovere et al., 2008; Pinna et al., 2017). Assessment of baroreflex in humans is usually approached measuring the changes in HR in response to changes in BP, the so-called baroreflex sensitivity (BRS). Alternatively, spontaneous beat-to-beat fluctuations of systolic arterial pressure and RR interval can be analyzed, allowing BRS assessment during daily-life. A wide spectrum of techniques has been used for spontaneous beat-to-beat BRS assessment. Traditional approaches, such as the sequence technique and those based on the spectral analysis of systolic arterial pressure and RR interval series (α index), were reviewed in La Rovere et al. (2008).

In order to deal with the nonstationary nature of cardiovascular variability, methods based on wavelet transform (Nowak et al., 2008; Keissar et al., 2010) and quadratic time-frequency representations (Orini et al., 2011) have been proposed. In Orini et al. (2012) a framework for nonstationary BRS assessment, based on a time-frequency distribution, was presented, taking into account the strength and prevalent direction of local coupling between RR variability (RRV) and systolic arterial pressure variability (SAPV) series. Alternatively, in Chen et al. (2011) dynamic assessment of BRS is accomplished based on a closed loop model within a point process framework. A critical review of clinical studies using spontaneous BRS was reported in Pinna et al. (2017). Despite some limitations, such as the lack of standards and the poor measurability in some patient populations, published studies support spontaneous BRS as a powerful tool for prognostic prediction in diseases such as hypertension, myocardial infarction, chronic heart failure and diabetes (La Rovere et al., 2008; Di Rienzo et al., 2009; de Moura-Tonello et al., 2016).

Spontaneous BRS assessment and monitoring during daily life is limited by the requirement of continuous BP recording, which is usually accomplished by the volume-clamp method or tonometry method, neither of them being suitable for ubiquitous monitoring (Mukkamala et al., 2015). This limitation may be overcome by using a surrogate of systolic arterial pressure which does not require the BP recording. Many works have attempted BP estimation based on pulse transit time (PTT), which is the time delay for the pressure wave to travel between two arterial sites. Most of these approaches, reviewed in Mukkamala et al. (2015), are based on models of arterial wall mechanics and wave propagation in the artery. Due to ease of measurement, pulse arrival time (PAT), which is the time delay between the electrocardiogram (ECG) waveform and a distal arterial waveform, has been widely used instead of PTT for BP estimation. PAT is the sum of PTT and the pre-ejection period (PEP), which varies beat-to-beat depending on ventricular and arterial pressures, short-term physiologic control and medication. Although the effect of PEP modulation makes PAT more inconvenient than PTT for BP estimation, half of the studies reviewed in Mukkamala et al. (2015) used PAT as a surrogate of PTT. Some of these methods have been used for BRS assessment. For instance, in Abe et al. (2015) it was proposed to evaluate baroreflex function using the maximum normalized cross-correlation between the LF components of HRV and PAT, derived from ECG and pulse photoplethysmographic (PPG) signals.

In Liu et al. (2011) it was suggested that PAT can track BP variations in HF range, but was inadequate to follow the LF variations. To overcome this limitation (Ding et al., 2016) proposed to estimate BP combining PAT with a new index, the photoplethysmogram intensity ratio (PIR), which can reflect changes in arterial diameter due to arterial vasomotion. In order to avoid PEP influence in BP estimation, PTT has been derived from impedance plethysmography recorded at the wrist and PPG at the finger (Huynh et al., 2018), or from a ballistocardiogram and PPG at the foot (Martin et al., 2016). Alternatively, PTT was estimated from two PPG signals recorded at ear and toe in Chen et al. (2009) and at forearm and wrist in Wang et al. (2018). Some works have investigated the correlation between PAT and PTT estimated from PPG signals at finger and forehead at rest (Liu et al., 2015) and during a tilt-test (Lázaro et al., 2016). In Li et al. (2014) different PPG indices were investigated for BP estimation. The time ratio of systole to diastole, time span of PPG cycle, diastolic time duration and area ratio of systole to diastole are at least as good as PTT for BP estimation, and can be derived from just one PPG signal.

The PPG signal can be acquired with a sensor placed in many places of the body. Furthermore, its recording is very simple, economical, and comfortable for the subject. Thus, PPG signal is a very interesting signal for ambulatory scenarios and wearable devices, and assessing BRS from PPG signal may have significant impact in such applications. Moreover, several studies have compared pulse rate variability (PRV), derived from the PPG to HRV derived from the ECG, reporting good agreement even in non-stationary situations and during abrupt autonomic nervous system changes (Gil et al., 2010; Wong et al., 2012; Posada-Quintero et al., 2013; Schfer and Vagedes, 2013). In this work we investigate the feasibility of assessing BRS solely from one PPG signal. The proposed approach is based on using PPG-based surrogates of RRV and SAPV series. On one hand, pulse-to-pulse variability (PPV) series was used as surrogate of RRV series. On the other hand, different PPG morphological features which are hypothesized to be related to the BP were used for generating series that were used as surrogates of SAPV. The ability of the proposed methods to capture changes in autonomic nervous system control was evaluated in a tilt-test database.

2. Materials and Methods

2.1. Data and Preprocessing

A data set containing ECG, BP, and PPG recordings from 17 healthy subjects (11 men), aged 28.5 ± 2.5 years, during a tilt table test was used for method evaluation. The protocol started with 4 min in supine position (Rest1), followed with 5 min in 70°-tilt-up position (Tilt), and ended with 4 min back to supine position (Rest2). The table took 18 s for automatic transitions between stages.

ECG lead V4 was recorded by Biopac ECG100C with a sampling rate of 1,000 Hz, BP signal (xBP(n)) was recorded by Finometer system with a sampling rate of 250 Hz, and PPG signal was recorded from the index finger by BIOPAC OXY100C with a sampling rate of Fs = 250 Hz. A low-pass filter with a cut-off frequency of 35 Hz was applied to the PPG in order to attenuate noise. This preprocessed PPG signal is denoted xPPG(n) in this paper. Several points were measured over the PPG pulses. Some of them were measured directly over the pulse as those described in section 2.1.1, and others over the waves extracted from the pulse by the pulse decomposition analysis (PDA) technique described in section 2.1.2.

2.1.1. Pulse Delineation

Several points of the ith PPG pulse were detected in order to take different morphological measurements. All these points are illustrated in Figure 1. First, PPG pulses were detected by an algorithm based on a low-pass derivative and a time-varying threshold (Lázaro et al., 2014). This algorithm detects the maximum up-slope point (nUi), and later it is used for detecting the pulse apex point (nAi) and the pulse basal point (nBi) as:

nAi=arg maxn{xPPG(n)},    n[nUi,nUi+0.3Fs]    (1)
nBi=arg minn{xPPG(n)},    n[nUi-0.3Fs,nUi].    (2)

Subsequently, nAi and nBi are used to compute the medium-amplitude point (nMi). This point is considered as a robust measure of PPG pulse location because it is located during the interval of the steepest slope of the PPG pulse, and it is set as:

nMi=arg minn{|xPPG(n)-xPPG(nAi)+xPPG(nBi)2|},              n[nBi,nAi].    (3)

Pulse onset nOi and end nEi points were detected based on the first derivative (Lázaro et al., 2013). In addition, pulse up-slope end nSEi was detected in a similar way. Let xPPG(n) be the first derivative of xPPG(n) computed by successive differences, after a 5-Hz-low-pass filter. Then, nSEi is set as:

nSEi=arg minn{|xPPG(n)-ηxPPG(nUi)|},    n[nUi,nAi],    (4)

where η was set to 0.05 similarly to the case of nOi and nEi (Lázaro et al., 2013).


Figure 1. Example of PPG pulse with its automatically detected points, and morphological measures taken from them: pulse amplitude (PA), pulse width (PW), and pulse slope transit time (PSTT).

2.1.2. Pulse Decomposition Analysis

PDA is a field in PPG signal proccessing that consits of modeling the PPG pulse as a main wave superposed with several reflected waves, increasing the robustness of some morphological measurements and even allowing others that would not be possible directly over the pulse. Several models can be found in the literature, based on different shapes including Gaussians (Baruch et al., 2011), LogNormal (Huotari et al., 2011), and Rayleight (Goswami et al., 2010). In this work, a modification to the PDA technique presented in Lázaro et al. (2018) is proposed. The main difference of this technique with respect to other PDA techniques in the literature is that the waves are extracted one-by-one, instead of fitting a several-waves-model at once. The modification proposed in this paper consists of not assuming a specific shape for the superposed waves, although it is assumed that they are symmetrical.

First, the baseline of the PPG signal was estimated by cubic-spline-interpolation of xPPG(nBi), and subsequently subtracted from xPPG(n). This baseline-removed version of PPG signal is denoted xPPGb(n) in this manuscript. Then, the beginning and the end of the ith PPG pulse were considered to be nBi and nBi+1, respectively. Note that this criterion ensures that each PPG pulse begins and ends with zero amplitude, as subtracted baseline was estimated at those nBi. Later, the algorithm extracts recursively the jth inner wave of the pulse by the following steps:

1. Set the beginning of the up-slope of the jth wave (nSOj,ib) as the previous to the first non-zero-amplitude sample. Note that in case of j = 1 (the main wave), this corresponds to nBi.

2. Set the end of the up-slope of the jth wave (nSEj,ib) as the first relative maximum.

3. Estimate the jth wave yj,ib(n) by concatenating the up-slope with itself horizontally flipped, assuming that it is symmetric:

xj,ib(n)={xPPGb(n),n[nSOj,ib,nSEj,ib]0,otherwise    (5)
yj,ib(n)=xj,ib(n)+xj,ib(-n+2nSEj,ib+1),    (6)

4. Substract yj,ib(n) to xPPGb(n) and go back to step 1 to continue extracting the (j + 1)th wave.

Once the desired number of waves have been extracted, they can be modeled in order to measure morphological features. In this work, three waves were extracted per PPG pulse. Subsequently, these yj,ib(n) were normalized to the unit in amplitude and to 1,000 samples by spline interpolation, and then they were modeled as Gaussian waves, each one defined by an amplitude, a mean, and a standard deviation (SD). Once these values are estimated, they were re-converted to the original scales of amplitude and time. An illustration of the steps of this algorithm can be observed in Figure 2.


Figure 2. Example of the steps for the pulse decompostion analysis that lead to extraction of main wave (first row), first reflected wave (second row), and second reflected wave (third row). In addition, the subsequent modeling of the extracted waves can be observed in the third column in magenta, red, and blue for the main, first reflected, and second reflected waves, respectively.

2.2. PPG-Based Surrogates of Systolic Arterial Pressure Variability for BRS Estimation

2.2.1. Systolic Arterial Pressure Variability Surrogates Based on Pulse Signal

Four pulse morphological features that have been related to the BP and/or to the arterial stiffness in the literature were measured from each PPG pulse: amplitude (PA), width (PW), up-slope (PUS), and slope transit time (PSTT). Pulse amplitude and width were measured as in Lázaro et al. (2013). The pulse amplitude corresponds to that amplitude reached by nAi with respect to nBi, and the pulse width was measured as the time interval between nOi and nEi. Pulse up-slope was measured as the first derivative value at nUi, and PSTT was measured as the time interval between nOi and nSEi. Later, PA, PW, PUS, and PSTT series were computed as:

dPAu(n)=i[xPPG(nAi)-xPPG(nBi)]δ(n-nMi)    (7)
dPWu(n)=i[nEi-nOi]δ(n-nMi)    (8)
dPUSu(n)=i[xPPG(nUi)]δ(n-nMi)    (9)
dPSTTu(n)=i[nSEi-nOi]δ(n-nMi),    (10)

where δ(·) denotes the Kronecker delta function, and the superscript “u” denotes that the signals are unevenly sampled, as the PPG pulses occur unevenly in time. A median-absolute-deviation outlier-rejection (Bailón et al., 2006) rule was applied to each one of these series, rejecting those points of the series that are outside the boundaries defined as the median ± 5 times the SD of the previous 50 points. Subsequently, a 4-Hz-evenly sampled version of each one of them was obtained by linear interpolation. The resulting signals are denoted using the same nomenclature, this time without the superscript “u” [e.g., dPA(n)].

2.2.2. Systolic Arterial Pressure Variability Surrogates Based on Pulse Decomposition Analysis

Seven morphological features were extracted from each PDA-based modeled PPG pulse. Specifically, the amplitude, mean, and twice the SD of the Gaussian-function fitted to the main wave were studied (mA1, mB1, and mC1, respectively). Moreover, the feature related to twice the SD of the first reflected wave were also studied (mC2), as well as the time delay between the main wave occurrence mB1 and those of reflected ones mB2 and mB3 (mT12 and mT13, respectively). Furthermore, the percentage of amplitude that it is lost in the first reflection was also estimated as:

mA12=mA1-mA2mA1.    (11)

Figure 3 illustrate these measures. These features extracted from the PDA are also hypothesized to be related to the BP and/or to the arterial stiffness since they are related to amplitude, relative position between the waves, and waves dispersion by SD. Their associated series were computed as:

d{A1, B1, C1, C2, T12, T13, A12}u(n)=im{A1, B1, C1, C2, T12, T13, A12}δ(n-nBi).    (12)

The outliers of these series were rejected by the same median-absolute-deviation-based rule applied in the case of the features which were measured over the pulse (see section 2.2.1), and similarly, they were linearly interpolated obtaining a 4 Hz evenly sampled version of each one of them denoted without the superscript “u”.


Figure 3. Measures over the extracted waves of an example PPG pulse. Gaussians fitted to these waves are shown in color.

2.3. Baroreflex Sensitivity Indices

The BRS indices were computed based on the α index, which is computed from a spectral analysis of RRV and SAPV. Several α-index surrogates based on PPG signal were computed, using PPV as RRV surrogate, and the SAPV surrogates described above.

The PPV was estimated using nMi as fiducial point:

dPPVu(n)=i1Fs[nMi-nMi-1]δ(n-nMi).    (13)

These series were also outlier-rejected and linearly interpolated to an even sampling rate of 4 Hz. Then, a power spectrum was computed from dPPV(n), obtaining S̄PPV(f), and from the kth SAPV surrogates, obtaining Sk(f), for each one of the stages of the protocol, where k can be PA, PW, PUS, PSTT, A1, B1, C1, C2, T12, T13, and A12. These power spectra were obtained by the Welch periodogram, using a 2 min Hamming window and 50% of overlap. Then, the PPG-based surrogates of the α index were extracted from these spectra, within LF ([0.04, 0.15] Hz) and HF ([0.15, 0.4] Hz) bands:

αk{LF,HF}=Ω{LF,HF}SPPV(f)df/Ω{LF,HF}Sk(f)df,    (14)

where ΩLF and ΩHF denote the LF and HF bands, respectively.

In addition, in order to take into account the non-stationarity of the protocol, the BRS indices were computed using a time-frequency technique for instantaneous measurement of α index, described in Orini et al. (2012). A time-frequency distribution was applied to dPPV(n) obtaining SPPV(n, f), and to each one of the PPG-morphology series used as SAPV surrogates obtaining Sk(n, f). In addition, a cross time-frequency spectrum SPPV,k(n, f) was also computed as in Orini et al. (2012). The instantaneous frequencies of the main components of SPPV,k(n, f) within [ 0.04, 0.15] Hz [for LF band, fLF(n)] and [0.15, 0.4] Hz [for HF band fHF(n)] were computed as the frequencies where SPPV,k(n, f) is maximum within those bands. Then, ΩLF(n) and ΩHF(n) were defined as the frequency bands centered at fLF(n) and fHF(n), respectively, with a bandwidth equal to the frequency resolution of the used time-frequency distribution. Then, the PPG-based surrogate of α index was computed for each Sk(n, f) as the square root of the ratio between the powers of dPPV(n) (as a surrogate of RRV) and dk(n), for each one of the defined bands:

αk{LF,HF}(n)=Ω{LF,HF}SPPV(n,f)df/Ω{LF,HF}Sk(n,f)df.    (15)

Figure 4 shows inter-subject median and interquartile range (IQR) of αk{LF,HF}(n) during the protocol. For BRS assessment, it is convenient to measure these indices only when PPV and k series are coupled. In order to detect these time courses, a time-frequency coherence (γPPV,k(n, f)) was computed, and PPV and k series were considered to be coupled in those areas where γPPV,k(n, f) is over a significance level. The indices αkLF(n) and αkHF(n) measured only when γPPV,k(n, f) is significant within ΩLF and ΩHF, respectively, are denoted αkLFγ(n) and αkHFγ(n), respectively, in this paper. Further details are given in Orini et al. (2012).


Figure 4. Median (black) and IQR (shaded area) of αLF(n), αHF(n), and their PPG-based surrogates, during the whole protocol. First 3 rows are related to the LF band and last 3 rows are related to the HF band.

For validation purposes, the conventional α index was also computed from the RRV and the SAPV, denoted with no subindex [α{LF,HF,LFγ,HFγ}(n)], and taken as reference. The RRV was computed by the interval function using the R points (nRi) determined from the ECG by Martínez et al. (2004):

dRRVu(n)=i1Fs[nRi-nRi-1]δ(n-nRi).    (16)

The SAPV was computed from the maximum of BP pulses (nAi), which were detected similarly to the case of the PPG pulses (see section 2.1.1):

dSAPVu(n)=ixBP(nAi)δ(n-nAi).    (17)

2.4. Performance Metrics

A unique value per subject and stage of the protocol (Rest1, Tilt, and Rest2) was obtained for each one of the three studied α-index estimation methods:

1. Welch-periodogram approach (α{LF, HF}): As it is based on a non-time-frequency technique, a unique value per subject and stage is available.

2. Time-frequency approach (ᾱ{LF,HF}): The median of α{LF, HF}(n) within each stage and each subject was taken as the unique value per subject and stage.

3. Time-frequency-coherence approach (ᾱ{LFγ,HFγ}): The median of α{LFγ,HFγ}(n) within each stage and each subject was taken as the unique value per subject and stage.

Then, correlation between the indices (α and αk) obtained from the 17 subjects and the 3 stages of the protocol (Rest1, Tilt, and Rest2) were computed. The distributions of these indices were found to be not normal by the Kolmogorov-Smirnov test. Thus, the Spearman's correlation coefficient was used. Furthermore, the Wilcoxon signed-rank test was applied to see if the indices can find significant (p < 0.05) differences between the different stages of the protocol.

As the SAPV surrogates, the α-index surrogates have different units and magnitude than classical α index. Thus, these surrogates cannot be directly compared to the classical α index, but their evolution can be compared. In order to do this, the relative variation of the α-index between consecutive stages was computed for each subject:

Δα=αS2-αS1αS1,    (18)

where αS1 and αS2 represent the studied index within stages S1 and S2, respectively.

A linear regression of the α-index surrogates which obtained best results in terms of correlation (those based on PUS, as it can be observed in section 3) was performed, obtaining similar units than those of the conventional α index (ms/mmHg). This linear regression was performed in order to compare those indices in a Bland-Altman plot. In addition, a multiple linear regression was performed using all the studied α-index surrogates in order to study whether their information is complementary or redundant. The combined α-index surrogates are denoted α^{LF,HF} (Welch-periodogram approach), ᾱ^{LF,HF} (time-frequency approach), and ᾱ^{LFγ,HFγ} (time-frequency-coherence approach).

3. Results

Table 1 shows inter-subject Spearman's correlation coefficients between α{LF, HF} and αk{LF,HF}, between ᾱ{LF,HF} and ᾱk{LF,HF}, and between ᾱ{LFγ,HFγ} and ᾱk{LFγ,HFγ}. The highest correlations were obtained for the α-index surrogates based on PUS. A scatterplot of these indices is shown in Figure 5. In addition, a Bland-Altmant plot of these indices and their associated conventional α indices is shown in Figure 6, after a linear regression in order to obtain similar units and magnitudes. The obtained limits of agreement were 0.94 ± 21.90 ms/mmHg (mean of the two values ± 1.96 × SD), –8.99E–15 ± 60.90, 1.29 ± 20.13, 4.56E–15 ± 40.49, 1.40 ± 18.43, 5.92E–15 ± 46.89 ms/mmHg, for αPUSLF, αPUSHF, ᾱPUSLF, ᾱPUSHF, ᾱPUSLFγ, and ᾱPUSHFγ, respectively.


Table 1. Inter-subject Spearman correlations of α{LF, HF} and αk{LF,HF} obtained in the different stages of the protocol.


Figure 5. Scatterplots of α vs. αPUS indices (first column), of ᾱ vs. ᾱPUS indices (second column), and of ᾱγ vs. ᾱPUSγ indices (third column).


Figure 6. Bland-Altman plots of α vs. αPUS indices (Welch-periodogram approach, first column), of ᾱ vs. ᾱPUS indices (time-frequency approach, second column), and of ᾱγ vs. ᾱPUSγ indices (time-frequency coherence approach, third column), after a linear regression to convert all units to ms/mmHg. Note that scales are not the same for LF band (first row) than for HF band (second row).

The Bland-Altmant plot obtained from the multiple linear regression using all the studied indices is shown in Figure 7. The obtained limits of agreement were 1.31 ± 20.38, 2.61E–15 ± 25.48, 1.12 ± 19.68, –6.30E–15 ± 20.00, 0.89 ± 10.25, –6.85E–15 ± 15.80 ms/mmHg, for αPUSLF, αPUSHF, ᾱPUSLF, ᾱPUSHF, ᾱPUSLFγ, and ᾱPUSHFγ, respectively.


Figure 7. Bland-Altman plots of α vs. its multiple-linear-regression-based combination of surrogates α^ (Welch-periodogram approach, first column), of ᾱ vs. its multiple-linear-regression-based combination of surrogates ᾱ^ (time-frequency approach, second column), and of ᾱγ vs. its multiple-linear-regression-based combination of surrogates ᾱ^γ (time-frequency coherence approach, third column). Note that scales are not the same for LF band (first row) than for HF band (second row).

Table 2 shows the inter-subject median and interquartile ranges of α{LF, HF} and αk{LF,HF}, for those indices which showed at least moderate Spearman's correlation coefficients (>0.50), and their relative increments (Δ). Significant differences (p < 0.05) of these indices between Tilt and rest stages are denoted with*. Similarly, the inter-subject median and interquartile ranges of medians of ᾱ{LF,HF} and ᾱk{LF,HF} are shown in Table 3, and interquartile ranges of medians of ᾱ{LFγ,HFγ} and ᾱk{LFγ,HFγ} are shown in Table 4.


Table 2. Inter-subject median and interquartile ranges of α{LF, HF} and αk{LF,HF}.


Table 3. Inter-subject median and interquartile ranges of intra-subject median of ᾱ{LF,HF} and ᾱk{LF,HF}.


Table 4. Inter-subject median and interquartile ranges of intra-subject median of ᾱ{LFγ,HFγ} and ᾱk{LFγ,HFγ}.

4. Discussion

Novel methods for measuring BRS using a PPG signal have been presented. They are based on surrogates of the α index, defined as the ratio of the power of RRV series and the power of SAPV series. In this work, PPV is used as a surrogate of RRV, and the SAPV is surrogated by different morphological features of the PPG pulse which have been related to BP in the literature. Some of these features are based on a novel PDA technique that has been presented in this paper. Many modeling functions have been applied to fit the PPG pulses in the literature. The novelty of the proposed PDA technique is that the used modeling function does not affect to the decomposition, as it is applied individually to the already extracted waves. It is worthy to note that the goal in this paper is not to obtain a very accurate measure of the studied morphological features, but in deriving a measure which is proportional to those features (as only their variability is needed). Keeping this in mind, a Gaussian function was used because it satisfies de goal while being a simple function that makes sense from the physiological point of view.

Three approaches were studied for estimating the α index from RRV and SAPV (or their surrogates): one based on Welch periodogram (α{LF, HF}), and two based on a time-frequency distribution which takes into account the non-stationarity of the protocol (Orini et al., 2012). This method redefines both LF and HF bands, making them time-varying following the dominant frequencies in such bands (ᾱ{LF,HF}). Alternatively, this method computes a time-frequency coherence between the RRV and the SAPV (or their surrogates), and estimates the α index in restricted areas where the obtained coherence is statistically significant, i.e., in those areas evidencing that RRV and SAPV (or their surrogates) are coupled (ᾱ{LFγ,HFγ}).

The correlation analysis shows how the PPG-based surrogates of α index track the changes of the conventional (ECG-and-BP-based) α index. Five out of the eleven SAPV surrogates leaded to α-index surrogates which obtained at least moderate correlation (>0.5). Those SAPV surrogates are, in order of cases getting the highest correlation: PUS, A1, PA, A12, and T12. Specifically, those α-index surrogates based on PUS obtained high correlation (>0.7) in all the cases. Those α-index surrogates based on A1 also obtained high correlation in four out of the six cases (αA1LF, αA1HF, ᾱA1HF and ᾱA1HFγ), while the remaining two (ᾱA1LF and ᾱA1LFγ) were very close to obtain it (correlation was 0.69 in both cases).

Results regarding the BRS assessment are shown in Table 2 for the Welch-periodogram approach (α{LF, HF}), Table 3 for the time-frequency approach (ᾱ{LF,HF}), and Table 4 for the time-frequency-coherence (ᾱ{LFγ,HFγ}). The conventional α indices showed significant differences between both rest stages and Tilt within both LF and HF, and for the 3 approaches. The highest difference was observed between Rest1 and Tilt within HF band (-220.44% for αHF, -178.17% for ᾱHF, and -189.00% for ᾱHFγ), while the smallest difference was observed between Rest2 and Tilt within LF band (37.87% for αLF, 31.21% for ᾱLF, and 35.37% for ᾱLFγ).

The only SAPV surrogate which led to α-index surrogates showing the same behavior than the conventional α-index was PUS. None of the other SAPV surrogates led to α-index surrogates finding significant differences between Rest2 and Tilt within LF (the smallest observed change) with the exception of PA when using the Welch-periodogram approach (αPALF). However, PUS, A1, PA, and A12-based α-index surrogates found significant Rest1 and Tilt within both the LF and HF band, and between Rest2 and Tilt within the HF band, for the 3 approaches. The T12-based α-index surrogates also found these differences for both the time-frequency and the time-frequency-coherence approaches while they found significant differences only within HF band for the Welch-periodogram approach. In general, those PPG-based-α-index surrogates exploiting the pulse amplitude (PA, A1, and A12) obtained better results than those exploiting the pulse dispersion (PSTT, B1, C1, C2, T12, and T13) for BRS assessment, with the exception of T12. However, the best results were obtained for the index derived from PUS, which exploits both PPG amplitude and pulse dispersion. Another possible reason of the better results obtained by PUS is that it is measured at the beginning of the pulse, which would be the part related to a unique wave (main wave, before superposition of reflections) containing the BP information better expressed than the reflected waves.

The Bland-Altman plots (Figure 6) for PUS-based α-index surrogates (after converting units to ms/mmHg by a linear regression) are wider for HF (±60.90 ms/mmHg, ±40.49 ms/mmHg, and ±46.89 ms/mmHg, for Welch-periodogram, time-frequency, and time-frequency-coherence approaches, respectively) than for LF (±21.90, ±20.13, and ±18.43 ms/mmHg, for Welch-periogram-, time-frequency-, and time-frequency-coherence approaches, respectively). When combining all the PPG-based α-index surrogates by a multiple-linear regression, these limits of agreement are narrower, specially for within HF (±25.48, ±20.00, and ±15.80 mm/mmHg, for Welch-periogram-, time-frequency-, and time-frequency-coherence approaches, respectively). These results suggest that there is complementary information among the SAPV surrogates and thus, they could be combined for improving the α-index surrogate. However, this combination may require a calibration process which may be subject-specific in a final application. Further studies including data from same subjects during different days must be elaborated in order to explore techniques to combine the information of the different α-index surrogates.

Comparing the correlations obtained by the PUS-based α-index surrogates among the three α-index estimation approaches, the highest correlation within LF was obtained when using the Welch-periodogram approach (0.81), while the highest correlation within HF was obtained when using the time-frequency-coherence approach (0.81). However, given the intrinsic non-stationarity of the cardiovascular system, our recommendation is to use the time-frequency-coherence approach (Orini et al., 2012) because it takes into account the time-varying dominant frequencies and the strength of the coupling between RRV and SAPV (or their surrogates) and thus, its estimates are more related to the BRS than the estimates from the other two approaches.

Based on these results, our recommendation for PPG-based BRS assessment is ᾱPUS{LFγ,HFγ}. First, ᾱPUSLFγ presented a significant decrease of more than 100% in median in tilt with respect to supine, which is in concordance with the decrease in reference ᾱLFγ. Second, ᾱPUSHFγ also presented a significant decrease in tilt with respect to supine, in this case around 2 times lower with respect to Rest1 (26.90%) than to Rest2 (61.31%), and these results are also in accordance to the reference ᾱHFγ (with 35.37% and 57.92%, respectively). It is worthy to note that the best α surrogate may be not derived from the best SAPV surrogate, because PPV was used as RRV surrogate while it is the sum of RRV and PAT variability (PATV) (Gil et al., 2010). Thus, for obtaining a exact surrogate for the ratio RRV/SAPV using PPV as numerator of the ratio, the best denominator is not exactly SAPV, but SAPV × (1+PATV/RRV).

These results support the potential value of the proposed index as a surrogate of BRS to monitor baroreflex impairment in certain applications. For example, in de Moura-Tonello et al. (2016) the square root of the RR and systolic BP series power (α index) at rest was significantly reduced (around 50%) in type 2 diabetes mellitus patients without cardiovascular autonomic neuropathy with respect to healthy controls of similar age and antropometric characteristics. In Ranucci et al. (2017) preoperative BRS was evaluated in 150 patients undergoing coronary surgery and related to postoperative complications such as atrial fibrillation, renal function impairment and low cardiac output syndrome. The α index was significantly lower (around 30% in median) in patients experiencing postoperative acute kidney dysfunction, as well as in patients with low cardiac output state (around 50% in median). However, clinical studies have to be elaborated in order to evaluate the proposed indices in different applications. To the best of our knowledge, this is the first time that these indices are studied for BRS assessment, so healthy volunteers with presumably efficient baroreflex were used in order to observe actual changes along the protocol. Different results may be obtained with patients of different diseases, specially taking into account that coherence is reduced in heart disease patients.

Results reported in this work suggest that BRS can be assessed with high correlation by only a PPG signal based on PPV (as RRV surrogate), and PPG-amplitude-based and/or PPG-dispersion-based features (as SAPV surrogates), being PUS the most convenient SAPV surrogate for BRS assessment. The PPG signal recording is simple, economical, and comfortable for the subject. Moreover, PPG signal can be acquired in many places of the body. Thus, these results are very interesting for ambulatory scenarios and for wearable devices. Future studies may include an surrogate of the α index using a combination of different PPG-based SAPV surrogates, specially amplitude- and dispersion-based features.

Ethics Statement

This study was carried out in accordance with the recommendations of Comité Ètico de Investigación Clínica de Aragón (CEICA) with written informed consent from all subjects. All subjects gave written informed consent in accordance with the Declaration of Helsinki. This study was exempt from approval of the ethical committee because no new data were registered for its development.

Author Contributions

JL and RB wrote the first draft in the manuscript preparation. All authors participated in the design of the proposed methods, as well of the statistical analysis for their evaluation, reviewed, and critiqued the manuscript preparation.


This project has received funding from the European Union's Framework Programme for Research and Innovation Horizon 2020 (2014-2020) under the Marie Skłodowska-Curie Grant Agreement No. 745755. This work was supported also by projects DPI2016-75458-R and RTI2018-097723-B-I00 funded by MINECO and FEDER, by Gobierno de Aragón (Reference Group BSICoS T39-17R) cofunded by FEDER 2014-2020 Building Europe from Aragon, and by CIBER in Bioengineering, Biomaterials and Nanomedicine (CIBER-BBN) through Instituto de Salud Carlos III. The computation was performed by the ICTS NANBIOSIS, specifically by the High Performance Computing Unit of CIBER-BBN at University of Zaragoza.

Conflict of Interest Statement

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.


BP, Blood pressure; BRS, Baroreflex sensitivity; ECG, Electrocardiograml; HF, High frequency; HR, Heart rate; HRV, Heart rate variability; LF, Low frequency; PA, Pulse amplitude; PAT, Pulse arrival time; PATV, Pulse arrival time variability; PDA, Pulse decomposition analysis; PEP, Pre-ejection period; PPG, Pulse photoplethysmography; PPV, Pulse-to-pulse variability; PSTT, Pulse slope transit time; PTT, Pulse transit time; PUS, Pulse up-slope; PW, Pulse width; RRV, RR interval variability; SAPV, Systolic arterial pressure variability; SD, Standard deviation.


Abe, M., Yoshizawa, M., Obara, K., Sugita, N., Homma, N., and Yambe, T. (2015). Evaluation of baroreflex function using green light photoplethysmogram in consideration of resistance to artifacts. Adv. Biomed. Eng. 4, 1–6. doi: 10.14326/abe.4.1

CrossRef Full Text | Google Scholar

Bailón, R., Sörnmo, L., and Laguna, P. (2006). A robust method for ecg-based estimation of the respiratory frequency during stress testing. IEEE Trans. Biomed. Eng. 53, 1273–1285. doi: 10.1109/TBME.2006.871888

PubMed Abstract | CrossRef Full Text | Google Scholar

Baruch, M. C., Warburton, D. E. R., Bredin, S. S. D., Cote, A., Gerdt, D. W., and Adkins, C. M. (2011). Pulse decomposition analysis of the digital arterial pulse during hemorrhage simulation. Nonlinear Biomed. Phys. 5, 1–15. doi: 10.1186/1753-4631-5-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y., Wen, C., Tao, G., Bi, M., and Li, G. (2009). Continuous and noninvasive blood pressure measurement: a novel modeling methodology of the relationship between blood pressure and pulse wave velocity. Ann. Biomed. Eng. 37, 2222–2233. doi: 10.1007/s10439-009-9759-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Z., Purdon, P. L., Harrell, G., Pierce, E. T., Walsh, J., Brown, E. N., et al. (2011). Dynamic assessment of baroreflex control of heart rate during induction of propofol anesthesia using a point process method. Ann. Biomed. Eng. 39, 260–276. doi: 10.1007/s10439-010-0179-z

PubMed Abstract | CrossRef Full Text | Google Scholar

de Moura-Tonello, S., Porta, A., Marchi, A., de Almeida Fagundes, A., Francisco, C., Rehder-Santos, P., et al. (2016). Analysis and baroreflex estimation in patients with type 2 diabetes in absence of any manifest neuropathy. PLoS ONE 11:e0148903. doi: 10.1371/journal.pone.0148903

PubMed Abstract | CrossRef Full Text | Google Scholar

Di Rienzo, M., Parati, G., Radaelli, A., and Castiglioni, P. (2009). Baroreflex contribution to blood pressure and heart rate oscillations: time scales, time-variant characteristics and nonlinearities. Philos. Trans. A Math. Phys. Eng. Sci., 367:1301–18. doi: 10.1098/rsta.2008.0274

PubMed Abstract | CrossRef Full Text | Google Scholar

Ding, X., Zhang, Y., Liu, J., Dai, W., and Tsang, H. K. (2016). Continuous cuffless blood pressure estimation using pulse transit time and photoplethysmogram intensity ratio. IEEE Trans. Biomed. Eng. 63, 964–972. doi: 10.1109/TBME.2015.2480679

PubMed Abstract | CrossRef Full Text | Google Scholar

Gil, E., Orini, M., Bailn, R., Vergara, J. M., Mainardi, L., and Laguna, P. (2010). Photoplethysmography pulse rate variability as a surrogate measurement of heart rate variability during non-stationary conditions. Physiol. Measure. 31, 1271–1290. doi: 10.1088/0967-3334/31/9/015

PubMed Abstract | CrossRef Full Text | Google Scholar

Goswami, D., Chaudhuri, K., and Mukherjee, J. (2010). A new two-pulse synthesis model for digital volume pulse signal analysis. Cardiovasc. Eng. 10, 109–117. doi: 10.1007/s10558-010-9098-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Huotari, M., Vehkaoja, A., Mtt, K., and Kostamovaara, J. (2011). Photoplethysmography and its detailed pulse waveform analysis for arterial stiffness. J. Mechan. Mater. Struct. 44, 345–362. Available online at:

Google Scholar

Huynh, T. H., Jafari, R., and Chung, W. (2018). Noninvasive cuffless blood pressure estimation using pulse transit time and impedance plethysmography. IEEE Trans. Biomed. Eng. 66, 967-976. doi: 10.1109/TBME.2018.2865751

PubMed Abstract | CrossRef Full Text | Google Scholar

Keissar, K., Maestri, R., Pinna, G., T La Rovere, M., and Gilad, O. (2010). Non-invasive baroreflex sensitivity assessment using wavelet transfer function-based time-frequency analysis. Physiol. Measure. 31, 1021–1036. doi: 10.1088/0967-3334/31/7/011

PubMed Abstract | CrossRef Full Text | Google Scholar

La Rovere, M. T., Pinna, G. D., and Raczak, G. (2008). Baroreflex sensitivity: measurement and clinical implications. Ann. Noninvas. Electrocard. 13, 191–207. doi: 10.1111/j.1542-474X.2008.00219.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Lázaro, J., Gil, E., Bailón, R., Mincholé, A., and Laguna, P. (2013). Deriving respiration from photoplethysmographic pulse width. Med. Biol. Eng. Comput. 51, 233–242. doi: 10.1007/s11517-012-0954-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Lázaro, J., Gil, E., Vergara, J. M., and Laguna, P. (2014). Pulse rate variability analysis for discrimination of sleep-apnea-related decreases in the amplitude fluctuations of pulse photoplethysmographic signal in children. IEEE J. Biomed. Health Inform. 18, 240–246. doi: 10.1109/JBHI.2013.2267096

PubMed Abstract | CrossRef Full Text | Google Scholar

Lázaro, J., Kontaxis, S., Bailón, R., Laguna, P., and Gil, E. (2018). “Respiratory rate derived from pulse photoplethysmographic signal by pulse decomposition analysis,” in 2018 40th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC) (Honolulu, HI), 5282–5285.

Google Scholar

Lázaro, J., Bailón, R., Laguna, P., Marozas, V., Rapalis, A., and Gil, E. (2016). “Difference in pulse arrival time at forehead and at finger as a surrogate of pulse transit time,” in 2016 Computing in Cardiology Conference (CinC) (Vancouver, BC), 269–272.

Google Scholar

Li, Y., Wang, Z., Zhang, L., Yang, X., and Song, J. (2014). Characters available in photoplethysmogram for blood pressure estimation: beyond the pulse transit time. Australas. Phys. Eng. Sci. Med. 37, 367–376. doi: 10.1007/s13246-014-0269-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, H., Ivanov, K., Wang, Y., and Wang, L. (2015). Toward a smartphone application for estimation of pulse transit time. Sensors 15, 27303–27321. doi: 10.3390/s151027303

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Q., Poon, C., and Zhang, Y. (2011). Time frequency analysis of variabilities of heart rate, systolic blood pressure and pulse transit time before and after exercise using the recursive autoregressive model. Biomed. Signal Proc. Control 6, 364–369.

Google Scholar

Martin, S. L.-O., Carek, A. M., Kim, C.-S., Ashouri, H., Inan, O. T., and Hahn Jin-Oh amd Mukkamala, R. (2016). Weighing scale-based pulse transit time is a superior marker of blood pressure than conventional pulse arrival time. Sci. Rep. 6:39273. doi: 10.1038/srep39273

PubMed Abstract | CrossRef Full Text | Google Scholar

Martínez, J., Almeida, R., Olmos, S., and Rocha, A. P. L. P (2004). A wavelet-based ecg delineator: evaluation on standard databases. IEEE Trans. Biomed. Eng. 51, 570–581. doi: 10.1109/TBME.2003.821031

PubMed Abstract | CrossRef Full Text | Google Scholar

Mukkamala, R., Hahn, J., Inan, O. T., Mestha, L. K., Kim, C., Treyin, H., and Kyal, S. (2015). Toward ubiquitous blood pressure monitoring via pulse transit time: Theory and practice. IEEE Trans. Biomed. Eng. 62, 1879–1901. doi: 10.1109/TBME.2015.2441951

PubMed Abstract | CrossRef Full Text | Google Scholar

Nowak, J., Ocon, A., Taneja, I., Medow, M., and Stewart, J. (2008). Multi-resolution wavelet analysis of time-dependent physiological responses in syncopal youths. Am. J. Phys. Heart Circ. Phys. 296, H171–H179. doi: 10.1152/ajpheart.00963.2008

CrossRef Full Text | Google Scholar

Orini, M., Bail, R., Mainardi, L., Laguna, P., and Flandrin, P. (2011). Characterization of dynamic interactions between cardiovascular signals by time-frequency coherence. IEEE Trans. Biomed. Eng. 59, 663–673. doi: 10.1109/TBME.2011.2171959

PubMed Abstract | CrossRef Full Text | Google Scholar

Orini, M., Laguna, P., Mainardi, L. T., and Bailón, R. (2012). Assessment of the dynamic interactions between heart rate and arterial pressure by the cross time frequency analysis. Physiol. Measure. 33, 315–331. doi: 10.1088/0967-3334/33/3/315

PubMed Abstract | CrossRef Full Text | Google Scholar

Pinna, G. D., Porta, A., Maestri, R., De Maria, B., Dalla Vecchia, L. A., and La Rovere, M. T. (2017). Different estimation methods of spontaneous baroreflex sensitivity have different predictive value in heart failure patients. J. Hyperten. 35, 1666–1675. doi: 10.1097/HJH.0000000000001377

PubMed Abstract | CrossRef Full Text | Google Scholar

Posada-Quintero, H. F., Delisle-Rodriguez, D., Cuadra-Sanz, M. B., and de la Vara-Prieto, R. R. F. (2013). Evaluation of pulse rate variability obtained by the pulse onsets of the photoplethysmographic signal. Physiol. Meas. 34, 179–187. doi: 10.1088/0967-3334/34/2/179

PubMed Abstract | CrossRef Full Text | Google Scholar

Ranucci, M., Porta, A., Bari, V., Pistuddi, V., and La Rovere, M. T. (2017). Baroreflex sensitivity and outcomes following coronary surgery. PLoS ONE 12:e175008. doi: 10.1371/journal.pone.0175008

PubMed Abstract | CrossRef Full Text | Google Scholar

Robertson, D., Diedrich, A., and Chapleau, M. (2012). Afferent baroreflex failure. Auton. Neurosci. 172, 1–3. doi: 10.1016/j.autneu.2012.10.010

CrossRef Full Text | Google Scholar

Schfer, A., and Vagedes, J. (2013). How accurate is pulse rate variability as an estimate of heart rate variability?: A review on studies comparing photoplethysmographic technology with an electrocardiogram. Int. J. Cardiol. 166, 15–29. doi: 10.1016/j.ijcard.2012.03.119

CrossRef Full Text | Google Scholar

Wang, Y., Liu, Z., and Ma, S. (2018). Cuff-less blood pressure measurement from dual-channel photoplethysmographic signals via peripheral pulse transit time with singular spectrum analysis. Physiol. Measure. 39:025010. doi: 10.1088/1361-6579/aa996d

PubMed Abstract | CrossRef Full Text | Google Scholar

Wong, J. S., Lu, W. A., Wu, K. T., Liu, M., Chen, G. Y., and Kuo, C. D. (2012). A comparative study of pulse rate variability and heart rate variability in healthy subjects. J. Clin. Monitor. Comput. 26, 107–114. doi: 10.1007/s10877-012-9340-6.

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: baroreflex, photoplethysmography, alpha index, autonomic nervous system, blood pressure, cardiovascular assessment

Citation: Lázaro J, Gil E, Orini M, Laguna P and Bailón R (2019) Baroreflex Sensitivity Measured by Pulse Photoplethysmography. Front. Neurosci. 13:339. doi: 10.3389/fnins.2019.00339

Received: 20 January 2019; Accepted: 22 March 2019;
Published: 18 April 2019.

Edited by:

Alberto Porta, University of Milan, Italy

Reviewed by:

Roberto Maestri, IRCCS Scientific Clinical Institutes Maugeri (ICS Maugeri), Italy
Luiz Eduardo Virgilio Silva, University of São Paulo, Brazil
Mathias Baumert, University of Adelaide, Australia

Copyright © 2019 Lázaro, Gil, Orini, Laguna and Bailón. 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: Jesús Lázaro,