Differential Analysis of Heart Rate Variability in Repeated Continuous Performance Tests Among Healthy Young Men
Article information
Abstract
Objective
Executive function correlates with the parasympathetic nervous system (PNS) based on static heart rate variability (HRV) measurements. Our study advances this understanding by employing dynamic assessments of the PNS to explore and quantify its relationship with inhibitory control (IC).
Methods
We recruited 31 men aged 20–35 years. We monitored their electrocardiogram (ECG) signals during the administration of the Conners’ Continuous Performance Test-II (CCPT-II) on a weekly basis over 2 weeks. HRV analysis was performed on ECG-derived RR intervals using 5-minute windows, each overlapping for the next 4 minutes to establish 1-minute intervals. For each time window, the HRV metrics extracted were: mean RR intervals, standard deviation of NN intervals (SDNN), low-frequency power with logarithm (lnLF), and high-frequency power with logarithm (lnHF). Each value was correlated with detectability and compared to the corresponding baseline value at t0.
Results
Compared with the baseline level, SDNN and lnLF showed marked decreases during CCPT-II. The mean values of HRV showed significant correlation with d’, including mean SDNN (R=0.474, p=0.012), mean lnLF (R=0.390, p=0.045), and mean lnHF (R=0.400, p=0.032). In the 14th time window, the significant correlations included SDNN (R=0.578, p=0.002), lnLF (R=0.493, p=0.012), and lnHF (R=0.432, p=0.031). Significant correlation between d’ and HRV parameters emerged only during the initial CCPT-II.
Conclusion
A significant correlation between PNS and IC was observed in the first session alone. The IC in the repeated CCPT-II needs to consider the broader neural network.
INTRODUCTION
Executive function, a fundamental domain of attention, encompasses top-down control processes that regulate behavior and facilitate adaptation to situational changes [1]. This domain includes critical components such as inhibitory control (IC), decision-making, working memory, and set shifting. Among them, the IC serves as the core component [2]. There are various tests employed to assess IC, with the Continuous Performance Test (CPT) being a notable tool and widely used clinically [3-5]. The CPT requires participants to respond to sequences of signals presented at varying intervals, a process that primarily evaluates IC and related impairments [6-9].
IC is linked to the prefrontal cortex and the autonomic nervous system (ANS), encompassing both the sympathetic (SNS) and parasympathetic nervous systems (PNS). This relationship is elaborated in the neurovisceral integration model, which articulates the interplay between the prefrontal cortex and ANS in top-down regulation [10-12]. Furthermore, both the polyvagal theory and the neurovisceral integration model emphasize collaboration between the prefrontal cortex and PNS across various hierarchical structures [12-14].
Heart rate variability (HRV) serves as a non-invasive indicator of the ANS, and its parameters have been identified as biomarkers of top-down self-regulation. The high-frequency (HF) component of HRV related to the IC has been explored in numerous studies [11,15]. It has been reported that the repeated stimuli over a one-week interval produce physiological variations in heart rate [16]. Moreover, an earlier study involving a 90-minute task of correcting transcription errors repeated five times daily over 3 consecutive days revealed a decrease in the mean RR interval and an increase in the low-frequency (LF) component of HRV, which was inconsistent with the previous studies. However, it remains uncertain whether these outcomes stem from prolonged testing fatigue or direct effects of task execution on different days [17]. It is hard to accurately determine the variations in HRV during the CPT measurements taken on different days for the lack of momentary analysis.
In our study, we analyzed the HRV parameters during relatively short-term executive tasks spaced 1 week apart. To capture moment-to-moment HRV fluctuations, we employed a standard 5-minute HRV time-window approach, with each window overlapping the next by 4 minutes, thereby creating 1-minute intervals between measurements. This methodology enabled us to assess (1) the involvement and correlation of the PNS with IC throughout the executive process and (2) the impact of repeated task execution on HRV dynamics during the task.
METHODS
Participants
This study was approved by the Institutional Review Board of the National Yang Ming Chiao Tung University (IRB: YM 108115EF). The authors assert that all procedures contributing to this work complied with the ethical standards of the relevant national and institutional committees on human experimentation and the Declaration of Helsinki of 1975, as revised in 2008. All participants were university or graduate students. We recruited 31 healthy men with body mass index between 18.5–28 kg/m2. They had never received the Conners’ Continuous Performance Test-II (CCPT-II) before. The study protocols were explained in detail to the participants, and their basic information was obtained after they provided written informed consent to participate. By querying participants, we excluded individuals who were receiving treatment for conditions such as cancer, renal or cardiovascular diseases, diabetes, depression, anxiety, psychotic disorders, or any condition related to substance withdrawal or use disorders. Participants were instructed to take a 10-minute break in a quiet room where the temperature was maintained at 25°C by an air conditioner. Then we attached the electrocardiogram (ECG) patch to obtain continuous ECG data on the day consent was obtained (Day 1). After a 5-minute ECG recording, the participants performed a complete 14-minute test. The participants were asked to sit still for 1 minute after the test. The ECG was recorded throughout the entire process, which was repeated on day 7.
CCPT-II
The CCPT-II is a widely used assessment tool for measuring the IC of the EF and is known for its internal consistency and test-retest reliability [18]. This test, conducted using CCPT-II software, spans 14 minutes and features alphabetic characters displayed for 250 ms each. Participants are seated before a computer in a quiet room, responding to the rapidly changing stimuli by pressing or not pressing the space bar for “nonX” and “X” signals, respectively. The CCPT-II software evaluates omissions, commissions, and detectability (d’) at the conclusion of the test [19]. Among these parameters, d’ serves as a psychophysical indicator of attentive executive performance, derived mathematically from signal detection theory [20-22] and widely applied across various conditions [23-25].
HRV
HRV parameters were primarily analyzed using both time-and frequency-domain methods, derived from 5-minute ECG signal windows. The significance of these parameters has been comprehensively reviewed in recent literature [26,27]. Among the time-domain HRV parameters, we measured the mean RR interval and standard deviation of NN intervals (SDNN) for each time window. The SDNN is posited to be more indicative of PNS activity when assessed over a short-term resting condition [26]. In the frequency-domain analysis of HRV, we examined the HF power (0.15–0.40 Hz), which is associated with PNS modulation and respiratory influences, also referred to as the respiratory band. Additionally, we analyzed the LF power (0.04–0.15 Hz), which was considered to reflect a combination of SNS and PNS activities in the past. However, recent evidence suggests that this is primarily due to baroreflex mechanisms in resting conditions [28-31]. We applied a logarithmic transformation to the HF and LF power data to correct potential skewness in their distributions [32].
HRV measurement
Throughout the test, the ECG was continuously recorded, with the participants maintaining a steady sitting position without significant movement. They were instructed to breathe naturally, avoiding both hyperventilation and paced breathing during the study [33]. An HRV analyzer (WG-103A, wegene) was used to acquire and process the ECG signals, employing an 8-bit analog-to-digital converter at a sampling rate of 256 Hz, enabling real-time ECG signal analysis, the fast Fourier transformation was applied to convert the data from the stationary RR intervals into the power spectrum of the frequency domain subsequently, as previously outlined [32]. To track immediate HRV fluctuations during repeated CCPT-II sessions, we employed a 5-minute time window approach. Each window overlapped the subsequent one by 4 minutes, thereby creating a 1-minute interval between consecutive time windows (as shown in Figure 1) [34]. The HRV value from the initial time window was established as baseline (t0). Additionally, we selected the data from the 14th minute time-window (t14) as the reactive HRV during CCPT-II. All the HRV measurements were included in Supplementary Table 1.

HRV was analyzed through a continuous approach, employing 5-minute ECG time-windows. Each time-window overlapped the subsequent one by 4 minutes, resulting in a 1-minute interval between time-windows. The dot plot displayed at the bottom illustrates the dispersion of RR intervals at the relevant time points. ECG, electrocardiogram; HRV, heart rate variability; CCPT-II, Conners’ Continuous Performance Test-II.
Statistical analysis
Statistical analyses were performed using the SPSS Statistics 20.0 (IBM Corp.). The Shapiro-Wilk test did not show a significant deviation from normality in the distribution of the investigated values. We employed one-way analysis of variance (ANOVA) to assess whether there were any significant differences among the tests for the CCPT-II parameters. Additionally, repeated measures ANOVA with Bonferroni correction was utilized to evaluate if there were significant differences among the tests for the HRV parameters. We evaluated the fluctuation of the HRV parameters (including the RR, SDNN, and both the logarithmically transformed LF and HF) using the 5-minute time-window with 1-minute interval. The difference between the corresponding HRV parameters and baseline (t0 time window) was analyzed using GEE. The correlation between the CCPT-II and HRV parameters in the entire process was obtained using the mean HRV data to represent the general HRV value in the entire CCPT-II process. To explore the association between HRV and executive function, we calculated the correlation coefficients between the mean HRV and HRV at the t0 and t14 time-windows, and the difference between the t0 and t14 time-windows during each CCPT-II. The above correlation coefficients with the corresponding CCPT-II parameters were calculated by Pearson’s correlation and all corrected with the participants’ age. Significance was tested at a 95% confidence interval. p<0.05 indicated statistical significance.
RESULTS
The demographic data and HRV metrics are presented in Table 1. There were no significant differences in CCPT-II performance and the corresponding HRV parameters across the repeated tests. This lack of variance suggests stability or homeostasis in cardiac autonomic modulation throughout repeated CCPT-II sessions.

Comparison of demographic data and variables of the study in different intervention times (N=31, male only)
Figure 2 illustrates the HRV parameters during each CCPT-II, highlighting the temporal patterns of HRV metrics in each time window. The analysis revealed significant differences in HRV parameters compared with the baseline time-window (t0) as follows: SDNN showed a significant decrease at t4 and from t6 to t14 in the first CCPT-II, with similar changes in the second CCPT-II. Low-frequency power with logarithm (lnlf) significantly decreased from t5 to t13 in the first CCPT-II and from t5 to t14 in the second CCPT-II that hinted the baroreflex might involve in the CCPT-II process. Neither the RR interval nor high-frequency power with logarithm (lnHF) at each time-window exhibited significant deviations from their respective baseline values at t0 time window.

Analysis parameters of heart rate variability from the t0 in each CCPT-II by GEE. The data is represented by the value derived from each minute time-window. *Denotes the significant difference at each time-window in comparison with t0 time-window in the first CCPT-II; †Denotes the significant difference at each time-window in comparison with t0 time-window in the second CCPT-II. p<0.05 indicates significance. RR, mean R-R interval of electrocardiogram; SDNN, standard deviation of NN intervals; lnLF, low-frequency power with logarithm; lnHF, high-frequency power with logarithm; CCPT-II, Conners’ Continuous Performance Test-II.
In Table 2, we presented the mean HRV parameters for each CCPT-II session and their correlations, as determined using Pearson’s correlation coefficient corrected by the participants’ age. There was no significant correlation between the HRV and the CCPT-II parameters at the baseline (t0) in both tests. Notably, the mean SDNN, LF, and HF demonstrated a significant positive correlation with detectability (d’) in the first CCPT-II session (p=0.012 for SDNN, p=0.045 for lnLF, and p=0.032 for lnHF).

Correlation between the HRV at different time-window and the CCPT-II parameters of each CCPT-II corrected with age
To elucidate the relationship between momentary HRV and d’ during each CCPT-II session, we computed the correlation coefficients between CCPT-II parameters and HRV at t14 and the difference baseline (t0), as detailed in Table 2. Notable correlations were predominantly observed in the first CCPT-II at different time windows: d’ showed a positive correlation at t14 time window with SDNN (R=0.578, p=0.002), lnLF (R=0.493, p=0.012), and lnHF (R=0.432, p=0.031); d’ showed a positive correlation by the difference between t14 and the baseline with lnLF (R=0.444, p=0.044). In the second CCPT-II, no HRV parameters showed significant correlation with CCPT-II parameters at different time-window.
DISCUSSION
This study is the first to meticulously track momentary HRV fluctuations during repeated CCPT-II sessions using a continuous 5-minute time-window approach with a 1-minute interval. We observed a significant correlation between lnHF and detectability (d’) in the first CCPT-II session. Furthermore, momentary changes in lnLF showed a significant decrease compared to baseline in both CCPT-II sessions, whereas lnHF did not exhibit significant variations. These findings are discussed in the context of their implications for understanding HRV dynamics in CCPT-II.
This study has several noteworthy findings. First, lnHF reflects cardiac vagal control mediated by the central autonomic network. The dynamic modulation of IC in response to external stimuli can be inferred from the observed increase or decrease in vagal control. This relationship illustrates how IC adaptively responds to environmental challenges, as indicated by changes in lnHF [10,35]. This adaptive modulation was further elucidated using a recently proposed “vagal tank” model, which suggested that IC in response to external stimuli induces changes in cardiac vagal control. According to this model, reduced vagal withdrawal during cognitive tasks indicates more effective self-regulation. This concept highlights the nuanced interplay between IC and autonomic responses, particularly in the context of cognitive demand [36-39]. In our study, lnHF did not exhibit significant changes compared to baseline during both CCPT-II sessions (Figure 2), suggesting few phasic fluctuations in vagal control during CCPT-II. We regard this stability in lnHFs as indicative of the PNS-efficient regulatory IC rather than potential fatigue during CCPT-II [40,41].
Second, prior research conducted by Siennicka et al. [42] demonstrated an association between higher resting HRV and attentional maintenance, although not explicitly with attentive performance. Importantly, their study collected HRV data exclusively in a resting state without considering HRV dynamics during task engagement. Indeed, the assessment of attentive executive performance is more accurate based on HRV measurements during task execution than at rest. In our study, we utilized d’ as an indicator of IC during the CCPT-II. A distinctive feature of this study was continuous HRV monitoring during the test. This approach allowed for a deeper understanding of the role of HRV in the assessment of IC, a core component of EF. Our results indicate significant correlations between both mean SDNN and mean lnHF with d’, aligning with the neurovisceral integration model and supporting previous assertions regarding the involvement of the PNS in regulating EF [10,43-45]. In this study, we furthered these assertions to IC in the first CCPT-II.
Third, our study found no significant correlation between lnHF and d’ in the second CCPT-II session, a finding that, at first glance, appears inconsistent with Thayer’s integrated neurovisceral model. This model typically posits a close relationship between vagal activity (as indicated by lnHF) and IC. The lack of correlation in the second session prompts further examination and suggests potential nuances in the model’s applicability across different contexts or repeated testing scenarios of IC [10,46]. Thayer’s model emphasizes that the prefrontal cortex is a key neural component in the regulation of goaldirected IC. However, prefrontal cortex is not critical for the performance under the repetition of likely tasks [47]. Indeed, IC is not only directed by the prefrontal cortex, but also the presupplementary motor area of cortex and the inhibitory level of cortex that attributed to the activation of cortical inhibitory neurons [48,49]. This suggests that involvement of the prefrontal cortex alone may not fully explain the observed correlations (or lack thereof) between repeated IC test and PNS activity. Our study put further evidence for the need to consider a broader neural network beyond the prefrontal cortex to thoroughly understand IC [50-53].
In summary, our study posits the following viewpoints: 1) there is efficient regulation of PNS in CCPT-II process; 2) the significantly PNS correlates with IC is only in the first CCPT-II; and 3) the relationship between IC and the PNS derived from repeated tests must considers the influence of broader neural networks other than the prefrontal cortex alone.
Limitations
This study had several limitations. This study was limited by its small sample size and narrow demographic range and focused on a specific age group and sex. This choice is necessitated by the sensitivity of HRV to age and hormonal fluctuations. Despite these constraints, the use of continuous HRV analysis, which is a more data-intensive and complex method than traditional approaches, underscores the value of our research. HRV was chosen for its clinical practicality and objective nature over other physiological measures such as blood pressure, respiration, or pupillometry. Additionally, the absence of blood pressure measurements necessitates a more cautious interpretation of baroreflex data. It is crucial to recognize that the CCPT-II is designed to assess IC over a brief period. Consequently, further studies are essential to explore the longitudinal relationships between IC and HRV, as well as the interactions between baroreflex and other facets.
Conclusion
Our study revealed that PNS modulation, as indicated by lnHF, efficiently modulated IC. PNS significantly correlated with IC only during the initial CCPT-II session. The lack of a significant correlation between HRV parameters and IC in subsequent sessions suggests that repeated tasks may involve broader neural networks beyond the prefrontal cortex. This study contributed to a deeper understanding of the complex relationship between the PNS and IC, highlighting the need to further explore this intricate interplay.
Supplementary Materials
The Supplement is available with this article at https://doi.org/10.30773/pi.2024.0251.
Comparison of demographic data and variables of the study in different intervention times
Notes
Availability of Data and Material
The datasets generated or analyzed during the study are available from the corresponding author on reasonable request.
Conflicts of Interest
Cheryl C.H. Yang receives a grant (112W32101) from the Brain Research Center, National Yang Ming Chiao Tung University and a grant (MOST 110-2629-B-A49A-501) from the National Science and Technology Council of Taiwan. All remaining authors have declared no conflicts of interest.
Author Contributions
Conceptualization: Chung-Chih Hsu, Tien-Yu Chen. Data curation: Chung-Chih Hsu, Terry B. J. Kuo. Methodology: Cheryl C. H. Yang. Project administration: Cheryl C. H. Yang. Resources: Terry B. J. Kuo. Software: Chung-Chih Hsu. Visualization: Jia-Yi Li. Writing—original draft: Chung-Chih Hsu. Writing—review & editing: Tien-Yu Chen.
Funding Statement
This study was supported by the Tri-Service General Hospital Research Foundation (TSGH-D-113147 and TSGH-D-114144).
Acknowledgments
None