Next Article in Journal
Do Rural Second Homes Shape Commensal Microbiota of Urban Dwellers? A Pilot Study among Urban Elderly in Finland
Next Article in Special Issue
Exploring Factors Impacting on the Lane Choice of Riders of Non-Motorized Vehicles at Exit Legs of Signalized At-Grade Intersections
Previous Article in Journal
Impact of Individual, Familial and Parental Factors on Adolescent Smoking in Turkey
Previous Article in Special Issue
Drivers’ Attention Strategies before Eyes-off-Road in Different Traffic Scenarios: Adaptation and Anticipation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Electrophysiological Brain-Cardiac Coupling in Train Drivers during Monotonous Driving

1
Edna Bennett Pierce Prevention Research Center, The Pennsylvania State University, 115 HHD Building, University Park, PA 16802, USA
2
Neuroscience Research Unit, School of Life Sciences, University of Technology Sydney, P.O. Box 123, Broadway, NSW 2007, Australia
3
Compumedics Ltd., 30-40 Flockhart Street, Abbotsford, VIC 3067, Australia
4
Interdisciplinary Sleep Medicine Center, Charité Universitätsmedizin Berlin, 10117 Berlin, Germany
*
Author to whom correspondence should be addressed.
Denotes an equal first author contribution.
Int. J. Environ. Res. Public Health 2021, 18(7), 3741; https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph18073741
Submission received: 10 March 2021 / Revised: 24 March 2021 / Accepted: 31 March 2021 / Published: 2 April 2021
(This article belongs to the Special Issue Driving Behaviors and Road Safety)

Abstract

:
Electrophysiological research has previously investigated monotony and the cardiac health of drivers independently; however, few studies have explored the association between the two. As such the present study aimed to examine the impact of monotonous train driving (indicated by electroencephalogram (EEG) activity) on an individual’s cardiac health as measured by heart rate variability (HRV). Sixty-three train drivers participated in the present study, and were required to complete a monotonous train driver simulator task. During this task, a 32 lead EEG and a three-lead electrocardiogram were recorded from each participant. In the present analysis, the low (LF) and high frequency (HF) HRV parameters were associated with delta (p < 0.05), beta (p = 0.03) and gamma (p < 0.001) frequency EEG variables. Further, total HRV was associated with gamma activity, while sympathovagal balance (i.e., LF:HF ratio) was best associated fronto-temporal delta activity (p = 0.02). HRV and EEG parameters appear to be coupled, with the parameters of the delta and gamma EEG frequency bands potentially being the most important to this coupling. These relationships provide insight into the impact of a monotonous task on the cardiac health of train drivers, and may also be indicative of strategies employed to combat fatigue or engage with the driving task.

1. Introduction

Professional driving is often characterised by lengthy periods of monotony, as drivers are frequently exposed to long stretches of relatively little stimulation and activity. Train driving is no exception, and is often considered more inherently monotonous than other driving industries because of the predetermined nature (with respect to time and destination/s) of the journey [1]. Indeed, the subjective experience of monotony has been associated with psychological distress and absence from work [2], cardiac risk factors [3], as well as reductions in task efficiency and performance [4]. Furthermore, train driving is often characterised by high workload, which results from rigid driving requirements and faced impediments, and has subsequently been linked to the development of driver fatigue. Driver fatigue has been associated with periods of hypovigilance [5,6], and increased likelihood of accidents [7,8,9]. This can have significant impact in train driving due to it inherently being a larger and heavier vehicle, which has a slower braking time and travels faster [10], and in the instance of passenger trains, the natural risk of carrying commuters.
Moreover, beyond the experience of driver fatigue, professional driving has been associated with a number of negative health consequences including stroke [11], coronary disease [12], ischemic heart disease [13], and myocardial infarction [14], as well as low physical activity [14] and other risk factors [15], which in turn have been associated with an increased overall disease burden [16] and risk of cardiovascular and metabolic diseases [16,17,18]. Additionally, it is important to note that health status has been shown to have an impact upon the safety of drivers [19], and consequently may contribute to accident incidence [20,21].
Numerous theories exist that aim at explaining this neural-cardiac link; for example, the psychophysiological coherence model postulates that specific cardiac activity patterns modulate the coherence between the cardiac and neural systems [22]. Whilst, neurovisceral integration theory [23,24] posits that cardiac vagal tone (as measured by heart rate variability; HRV) is associated with neural network functionality in affective-cognitive contexts. In all instances, electrophysiological research bridges the gap between driver fatigue and monotony and the possible impact on cardiac function of repeated/prolonged exposure to such conditions.
Previous research has primarily utilised spectral electroencephalography (EEG) to investigate monotonous driving and predict the onset and impact of driver fatigue [25,26,27,28,29,30]. Collectively, this research indicates that low frequency delta and theta brain activity both increase in fatigue and sleepiness [31,32,33,34], and are accompanied by a decrease in high frequency beta and gamma activity [35,36,37]. Heart rate variability has perhaps been more widely utilized with respect to its in electrophysiological research, having been used to investigate a wide variety of physiological conditions (e.g., diabetes and stroke), and psychological states (e.g., stress and fatigue) [38,39,40,41]. Broadly, HRV has been demonstrated to reflect the activity of the autonomic nervous system and its sympathetic and parasympathetic branches, [42] which control almost all visceral, vascular and metabolic functions [43]. Most importantly with respect to train driver health, reductions in HRV parameters have been shown to function as predictors for general mortality [44] and the development of cardiovascular risk factors including hypertension and obesity [45].
However despite the widespread use of both spectral EEG and HRV in electrophysiological and health research, studies investigating the association between these two approaches is often confined to acutely stressful or otherwise affective situations [46], and the investigation of such pairing in monotonous contexts is limited. Therefore, the present research aimed to investigate the relationship between brain activity and heart rate variability recorded during monotonous driving in train drivers. Based upon the aforementioned literature, e.g., [34,35,37] it was hypothesised that: (1) slow-wave EEG parameters (i.e., delta and theta activity) that are representative of monotony, fatigue and sleepiness would be associated with parasympathetic HRV, and; (2) fast-wave EEG parameters (i.e., alpha and beta activity) that are representative of alert and active states would be associated with sympathetic HRV.

2. Materials and Methods

2.1. Participants

Sixty-three train drivers (54 male, 9 female) aged between 24 and 69 years (average 38.56 ± 9.96 years) were recruited from the train driver population, and via Sydney Trains in Australia to participate in the current study. The Lifestyle Appraisal Questionnaire [47] was used to screen all recruited individuals for pre-existing medical conditions, as well as regular medication and other non-medical substance use that may have affected physiological measures. Additionally, pre-study blood pressure values of greater than 160 mmHg (systolic) and/or 100 mmHg (diastolic) were also used as an exclusion criteria, as these values are associated with an increased risk of cardiovascular events [48]. No recruited individual was excluded from participating in the present study. Finally, the study protocol was approved by the Human Research Ethics Committee, University of Technology Sydney. All participants provided written informed consent prior to participating.

2.2. Experimental Tasks and Procedures

Initially, participants were seated for a 5 min resting period, after which three blood pressure measurements using an automated sphygmomanometer (OMRON AIB) were obtained with a two minute rest between each reading as is recommended [49]; these values were then averaged to provide a pre-study blood pressure value for each participant. Following this, the participant was asked to completed the Lifestyle Appraisal questionnaire [47], which enabled the collection of participant demographic data.
Upon completing the questionnaire, the electrophysiological monitors including a three lead electrocardiogram, and a 32 lead electroencephalogram were connected to the participant. ECG monitoring was completed using a FlexComp Infinity encoder (Thought Technology Ltd., Montreal, QC, Canada) supported by the BioGraph Infinity software package (Thought Technology Ltd., Canada). Prior to fixation, the adhesion sites for the ECG electrodes were cleaned with an alcohol wipe in order to remove surface debris and decrease any noise artefacts associated with an ECG signal. Following this, the ECG electrodes were placed with the active electrodes positioned at the intercostal space between the 4th and 5th ribs, two centimetres laterally from each side of the sternum, and the reference electrode as secured underneath the shoulder.
With respect to the EEG data capture, a 32 lead unipolar EEG montage with a sampling rate of 1000 Hz was recorded using a SynAmp2 amplifier (Compumedics Ltd., Abbotsford, VIC, Australia) supported by the Scan software (Version 4.3, Compumedics Ltd., Abbotsford, VIC, Australia). A 64 channel electrode QuikCap (Compumedics Ltd., Australia) was used; filled with Signa Gel (Parker Laboratories Inc, Fairfield, NJ, USA) at the following locations: Fp1, Fp2, F7, F3, Fz, F4, F8, FT7, FC3, FCz, FC4, FT8, T7, C3, Cz, C4, T8, TP7, CP3, CPz, CP4, TP8, P7, P3, Pz, P4, P8, O1, Oz, and O2. Further, the reference electrode was placed at the vertex, and the ground electrode at position AFz. Additionally, a bipolar vertical electrooculogram (EOG) was set up with one electrode above the left eye, and a second electrode below.
Electrophysiological data collection was conducted during a 10-min resting baseline, as well as an active phase, which involved 10 min of train driving using the Trainz Classics (N3V Games, Helensvale, Australia) driving simulator program. During this active phase, participants engaged in a first-person train driving simulation task, in which they were required to correctly adhere to signals whilst undertaking straight-line track driving between stations and occasionally stopping at stations. After the active phase had been completed, the electrophysiological data collection was concluded. Finally, as an ethics requirement, three post-study blood pressure measures were again obtained and the experimental protocol was concluded.

2.3. Electrophysiological Data Processing

Prior to statistical analysis, both sets of raw electrophysiological signals (ECG and EEG) were pre-processed, so as to obtain the relevant and summarizing variables in the time and frequency domains.

2.3.1. Heart Rate Variability

The recorded ECG data was processed utilising the Kubios HRV Premium software (Version 3.1.0, Kubios Oy, Kuopio, Finland) to derive the activity for the following HRV frequency bands: low frequency (LF; 0.04–0.15 Hz) and high frequency (HF; 0.15–0.4 Hz). The automatic artefact correction process built into the Kubios software [50], as well as the Smoothn Priors method [51] of trend component rejection were utilised to correct any artefacts and ectopic beats present in the raw ECG recordings.
All variables were derived from approximately 10 min of ECG [52], utilising Welch’s periodogram method [53] applied to a 300 s window with a 50% overlap. Further, it should be noted that for each of these variables raw power units (ms2) as well as natural log transformed (ln) and normalised units (n.u.) were utilised. Additionally, the ratio between LF and HF (LF/HF) as well as overall HRV power or total power (TPow) was also computed.

2.3.2. Electroencephalography Data

The raw EEG data sets were filtered using a Butterworth IIR bandpass filter set at 0.5 and 50 Hz, which was followed by a Hanning window. The Aligned-artefact average procedure utilising vertical EOG data [54,55] was then applied to minimise any eye movement artefacts, and any muscle artefacts and movement artefacts were also removed as well. After pre-processing, both the baseline and active recordings were sectioned into approximately 300 two second epochs, and the EEG activity in the delta (0.5–4 Hz), theta (4–8 Hz), alpha (8–13 Hz), beta (13–35 Hz) and gamma (35–50 Hz) frequency bands [56] was calculated via Periodogram power spectral density estimate. The epoch values of each recording were scanned for outliers, which were removed using a modified z-score statistic [57] greater than or equal to 5, which was calculated using the following equations:
Z = X x ~ MAD
where: X = Epoch value; x ~ = Median value; MAD = Median Absolute Deviation.
Median Absolute deviation can be calculated using the following equation:
MAD = x ~ ( | X i x ~ ( X i ) | )
where: x ~ = Median; Xi = Epoch value.
After the removal of epoch outliers, the activity values of each epoch were averaged to derive a single value in each frequency band for each electrode location. Additionally, all activity values in each EEG frequency band were also averaged, providing an overall average EEG activity value. Finally, after all EEG data had been collated, outlying values were again removed using the modified Z-score statistic previously applied to the epoch values; however, the threshold at which values were removed in this instance was set at a z statistic greater than or equal to 10 [57].

2.4. Data Analysis

The present statistical analysis was conducted using STATISTICA (Version 10, 2011, Statsoft Inc., Tulsa, OK, USA) with statistical significance set at p < 0.05. Partial Pearson’s correlation analysis controlling for participant age and body mass index (BMI) was used to determine any associations between EEG variables and HRV variables. If a significant association between three or more independent EEG variables and a single dependent HRV variable were identified via correlation analysis, multiple general linear regression analysis was utilised to determine the strongest individual predictor.

3. Results

Log-transformed low frequency HRV (LFln) was significantly and negatively correlated to frontal pole theta and alpha activity as well as fronto-temporal beta activity (Table 1). Furthermore, positive associations were found between LFln and occipital delta and gamma activity, as well as parietal gamma activity (Table 1).
A forward stepwise general linear regression analysis informed by the significant correlations was performed to examine the relationship between LFln HRV and EEG parameters (Table 2). Of the 6 originally entered variables, the regression analysis for LFln retained 3 (O1 delta activity, FT8 beta activity, and O1 gamma activity), and had an overall significance of p = 0.01. Together, these 3 variables explained 52.4% of the variance in LFln (F = 5.14; df = 3; p = 0.01; R = 0.72; R2 = 0.52; AR2 = 0.42). Furthermore, FT8 beta activity also presented as independently significant predictors of LFln (p = 0.03) in the regression.
When examining normalised LF HRV, the present analysis determined positive correlations for central delta activity, and frontal pole alpha activity (Table 3). Furthermore, parietal beta and gamma activity, and occipital gamma activity were significantly and negatively correlated with the LFn.u. parameter.
To further examine the relationship between LFn.u. HRV and EEG parameters, a forward stepwise general linear regression analysis informed by the significant correlations was performed (Table 4). Of the variables entered into the regression analysis, two were retained (C3 delta activity and O1 gamma activity) with an overall significance of p < 0.001. Together, these 2 variables explained 57% of the variance in LFn.u. (F = 13.11; df = 2; p < 0.001; R = 0.75; R2 = 0.57; AR2 = 0.52). Furthermore, both variables also presented as independently significant predictors of LFn.u. (p < 0.05 and p < 0.001, respectively).
With reference to high frequency HRV, HFln was significantly negatively correlated to fronto-temporal and central delta activity, as well as frontal pole alpha and fronto-temporal beta activity. Furthermore, parietal and occipital gamma activity were positively correlated to the HFln parameter (Table 5).
The associations of log-transformed HF HRV with EEG parameters were further investigated using a forward stepwise general linear regression analysis informed by the significant correlations (Table 6). Of the 7 originally entered variables, the regression analysis for HFln retained 3 (C3 delta activity, FT8 beta activity and O1 gamma activity), and had an overall significance of p < 0.001. Together, these 3 variables explained 62% of the variance in HFln (F = 9.14; df = 3; p < 0.001; R = 0.79; R2 = 0.62; AR2 = 0.55). Furthermore, O1 gamma activity also presented as an independently significant predictor of HFln. (p < 0.001).
Similar results were determined for normalised HF HRV, which was also significantly and negatively correlated to central delta activity and frontal pole alpha activity (Table 7). Moreover, parietal beta and gamma activity, and occipital gamma activity were positively correlated with HFn.u.
Similarly, a forward stepwise general linear regression analysis investigated the associations between normalised HF HRV with EEG parameters (Table 8). Of the 6 originally entered variables, the regression analysis for HFn.u. retained 2 (C3 delta activity, and O1 gamma activity), and had an overall significance of p < 0.001. Together, these variables explained 56% of the variance in HFn.u. (F = 12.97; df = 2; p < 0.001; R = 0.75; R2 = 0.56; AR2 = 0.52). Furthermore, both variables also presented as independently significant predictors of HFn.u. (p =< 0.05 and p < 0.001, respectively).
Total HRV power was negatively correlated with frontal pole delta, theta and alpha activity, as well as fronto-temporal beta activity (Table 9). Furthermore, positive associations between TPow and parietal and occipital gamma activity were also determined.
Log-transformed total heart rate variability power and its association with EEG parameters was also further investigated using forward stepwise general linear regression analysis (Table 10). Of the 6 originally entered variables, the regression analysis for TPowln retained 2 (P4 and O1 gamma activity), and had an overall significance of p < 0.001. Together, these 2 variables explained 65% of the variance in TPowln (F = 15.98; df = 2; p < 0.001; R = 0.81; R2 = 0.65; AR2 = 0.61). Furthermore, O1 gamma activity also presented as an independently significant predictor of TPln. (p < 0.001).
Finally, the present analysis determined that the LF/HF ratio was significantly and positively correlated to fronto-temporal delta and theta activity, as well as fronto-central theta activity, frontal pole alpha activity, centro-parietal beta activity, and fronto-central and centro-parietal gamma activity (Table 11).
The associations between the LF/HF ratio and EEG parameters were further investigated using a forward stepwise general linear regression analysis informed by the significant correlations (Table 12). Of the 8 originally entered variables, the regression analysis retained 2 (FT8 delta activity, and FC4 theta activity), and had an overall significance of p < 0.001. Together, these 2 variables explained 61% of the variance in LF/HF (F = 13.07; df = 2; p < 0.001; R = 0.78; R2 = 0.61; AR2 = 0.56). Furthermore, FT8 delta activity also presented as an independently significant predictor of LF/HF (p = 0.02).

4. Discussion

The present study aimed to investigate the relationship between brain activity and HRV recorded in train drivers during a monotonous train driving task. Associations between independent frequency domain HRV parameters (LF and HF) and both low and high frequency EEG variables were observed. More specifically, increased LF HRV was associated with increased delta activity, reduced beta activity, and a location dependent increase or decrease in gamma activity. Increased HF HRV was associated with reductions in delta and beta activity, as well as increased gamma activity. Interestingly, sympathovagal balance (LF:HF) was only correlated with low frequency EEG parameters, being positively associated with fronto-temporal delta and fronto-central theta activity.

4.1. Low and High Frequency Heart Rate Variability

Within frequency domain HRV research, the LF parameter has typically been considered to be in part representative of sympathetic nervous system activity, particularly once it is transformed or normalised [52]. Presently, increased sympathetic activity (as indexed by LF HRV) was associated with increased central and occipital delta activity, reduced fronto-temporal and parietal beta activity, and a potential gamma oscillation, whereby both an increase and decrease in occipital gamma was found.
Previous electrophysiological investigations of fatigue and sleepiness have demonstrated that delta activity increases [32,33,34], and high frequency beta and gamma activity decreases [35,36,37] during these states. When considered in this context, it stands to reason that the observed correlation and increase in low frequency EEG, i.e., delta activity and reduction in high frequency EEG (i.e., beta and gamma) activity may be indicative of the train drivers actively experiencing monotony-related fatigue [31]. Alternatively or perhaps supplementarily, past research has noted that professional drivers, including train drivers, engage and utilise various tactics to actively combat and compensate for fatigue, most of which are typically associated with increases in the amount of consciously applied effort [58,59]. Activation of the sympathetic nervous system fundamentally represents the fight or flight response and it’s underlying mechanics, whereby an individual actively either engages (i.e., fights) or disengages (i.e., flights) with a distressing stimuli [60]; and so, it could be postulated that these combative efforts may be facilitated by or at the very least associated with increases in sympathetic activity. Thusly, it could be suggested that the neural-cardiac coupling between an increase in LF HRV and an associated increase in slow-wave and decrease in fast-wave EEG activity may electrophysiologically represent attempts to more actively engage with the driving task, or the combative responses that train drivers often employ in an attempt to overcome experienced monotony/fatigue. However, it should be noted that this is speculative, and future studies may wish to actively image both the task engagement and fatigue ‘countermeasures’ used by train drivers to explicitly examine these effects.
Continuing the focus on LF HRV, previous research has associated elevated LF HRV to have been linked with the increased presence of cardiac risk factors and cardiac disease [61,62]. Furthermore, physiologically constant and prolonged sympathetic activity has been show to place potentially damaging pressure on cardiac health [63]. As such, it could be inferred that the experience of monotony and fatigue is increasing the cardiac risk and threatening the cardiac health of these train drivers. A particularly poignant inference, as professional driving as an industry itself has also been previously associated with increased cardiac risk and chance of developing cardiac disease [12,13,14]. Further, it is possible that these bodies of literature are describing the same effect, i.e., the increase in LF HRV and cardiac risk, is the same industry specific increase in cardiac risk. However, it may also be possible that the observed cardiac effects could compound and/or potentially exacerbate one another, and future research would do well to consider separating the two effects (i.e., the cardiac impact of monotony, and the cardiac impact of professional driving) to explicitly test this. Nonetheless, not only do the present associations between LF HRV and EEG provide insight into neural-cardiac coupling during monotony and the potential cardiac risk of such conditions, but also demonstrates the potential to form a basis from which such electrophysiological monitoring could be applied towards the active assessment of driver health in situ, and in monotonous environments.
With respect to HF HRV, it is thought to be the complement to LF HRV and so is considered to be representative of parasympathetic nervous system activity [38]. In the present analysis increased parasympathetic activity (i.e., HF HRV) was associated with reductions in central delta and fronto-temporal beta activity, and increases in occipital gamma activity. Interestingly, the traditional “rest and digest” characteristics that are attributed to parasympathetic activity would suggest that rather than the aforementioned observed associations, one would expect to see an increase in delta activity and a reduction in high frequency EEG activity. Furthermore, these results are close to the inverse (i.e., decreased slow-wave activity and increased fast-wave activity) of the correlations that were observed when examining the neural coupling of LF HRV, potentially providing a look at the push-pull interaction of the arms of the autonomic nervous system, and perhaps an alternate parasympathetic view of monotony characterized by reductions in delta and beta activity combined with potential gamma bursts. As to what these correlations may represent, perhaps given the separate nature of LF and HF with respect to autonomic function, this second set of correlations reflects monotony-related task disengagement, which may be driving the previously purported increase in conscious effort. Alternatively, it is quite possible that these correlations simply represent some unmonitored cognitive component of the simulator task like the driving events (e.g., track signal changes), in particular the increase in gamma activity, which has often been correlated with cognitive events [64]. On that note, future research investigating not only neural-cardiac coupling but monotony and fatigue more broadly should consider eye-tracking or other techniques that enable the monitoring of participant attention to visualise both disengagement and any increases is cognitive effort, as well as the flagging of all events (e.g., track signals, station stops, etc.) that participants experience.

4.2. Total Heart Rate Variability and Sympathovagal Balance

Total heart rate variability was also considered in the present analysis and was most strongly associated with parietal and occipital gamma activity, whereby an increase in total HRV power would be accompanied by increases in these high frequency EEG parameters. Additionally, the present analysis also strongly associated increases in sympathovagal balance (as represented by the LF/HF parameter) with increases in fronto-temporal delta and fronto-central theta activity. Given the observed correlation analyses for both the LF and HF HRV parameters, it seems reasonable to suggest that the series of correlations between TP and LF/HF HRV and EEG parameters are most likely attributable to the influence of LF HRV, and could be similarly indicative of a state of sympathetic dominance, whereby an individual reports LF HRV greater than their HF HRV.
Within this context, it again should be noted that sympathetic dominant HRV has been associated with cardiac risk [62], and so the results observed for total HRV power and sympathovagal balance not only provide further support for neural-cardiac coupling experienced during monotony, but in-turn also provide some insight into the cardiac health of train drivers, as well as potentially that of other workers and operators in other industries that also frequently experience monotonous or otherwise stressful work settings.

4.3. Limitations & Future Directions

As with all research, the present study is not without its limitations. The present protocol utilised an approximately 10-min monotonous train driving task to collect and examine HRV and brain activity. Whilst the task itself is designed to be inherently monotonous, it is relatively short, and a lengthier driving task [9,29] may provide greater insight into monotony and its cardiac implications. This could be particularly important if the electrophysiological parameters capture any potential phasic shifts related to individual awareness of experiencing monotony and fatigue, and/or any combative actions and/or increases in conscious effort made by drivers to counter their experience of fatigue [59].
In relation to the electrophysiological parameters, the present work utilised a relatively small 32-lead monopolar EEG montage and future research may wish to expand on the electrode montage to provide greater nuance or open up the possibility of using source-localisation approaches. Further, expanding the derived variables of interest from channel specific spectra to include synchronicity, coherence [65], interaction parameters, i.e., power spectral ratios, and/or other non-traditional EEG [66], may allow for a better representation of monotony and fatigue [26,67], and hold additional insight into the cardiac impact of this experience [40].
Similarly, the present experiment confined the derived HRV parameters to only the traditional frequency domain variables, and an expansion into examining time-domain and other non-linear HRV parameters, such as entropy values, may also provide additional insight into brain-cardiac coupling in monotony and its implications in train drivers. Additionally, beyond the specific parameters themselves, auto-regressive HRV generation/pre-processing techniques also warrant consideration in future research. Further, with respect to the HRV data, the two datasets (i.e., log-transformed and normalised) differently represent the activity of the autonomic nervous system. More specifically, log-transformed data broadly comments on autonomic data, whilst normalised HRV focuses more specifically on the individual autonomic arms. To that end, while this distinction has no bearing on the present work, future studies may wish to engage with this difference more actively.
Moving forward, with respect to analytical options, future research should consider the value of both repeated measures designs and the application of data normalization techniques (e.g., independent/principle component analysis) in their handling of electrophysiological data as it is possible that such approaches may provide additional insight into the link between the two physiological measures. Moreover, future research investigating brain-cardiac coupling within a monotonous or fatigued context should consider moving towards more realistic testing paradigms/conditions. A reasonable starting point would be to capture data either across or at multiple time points during the day/night cycle, as night driving has been associated with variations in markers of sleepiness and fatigue when compared to daytime driving [68,69,70]. Furthermore, electrophysiological parameters like EEG activity and HRV can be prone to circadian influence, and whilst this variance is accounted for in the present experimental design (by restricting the time point of data collection across all participants), it would be poignant for future research to track the relationship of EEG parameters and HRV over the course of day. Finally, even though simulated tasks cover many of the fundamental components of driving, they are inherently different, e.g., simulators make use of computer controlled or AI peers rather than true confederates. And so, future studies should consider the value of transitioning to the use of both controlled and “real” driving conditions and/or closed circuits for their driving tasks rather than the presently used lab-based simulator task.

5. Conclusions

The present research investigated the relationship between HRV and brain activity recorded during monotonous driving in train drivers, and coupled the two sets of electrophysiological parameters. The analyses found a number of associations between both low and high frequency HRV and EEG activity primarily in the delta and gamma frequency EEG bands. The coupling between these electrophysiological parameters provides some insight, not only into the direct cardiac impact of monotony experienced by train drivers, but also the subsequent longer term effects that train drivers may experience due to their prolonged and repeated exposure to a monotonous work environment. However, it is equally possible that these relationships and the related neural-cardiac coupling reflect combative tactics employed by drivers as they attempt to actively overcome fatigue and/or improve their engagement with the driving task. Finally, looking more broadly the interrelationship demonstrated between neural and cardiac biomarker monitoring could have major implications, subject to further investigations, as it relates to the routine use of fatigue or drowsiness monitoring systems.

Author Contributions

S.L. (Sara Lal), D.B., E.Z., S.L. (Shail Lal) and T.P. contributed to the original study conception and design, as well as set-up, technical, and administrative work; and were pivotal in obtaining funding for the project. S.L. (Sara Lal) supervised hired staff and students on the project. T.C. prepared all experimental materials, and conducted all data collection. T.L. conducted the data analysis, and drafted the manuscript. All authors contributed to and edited multiple versions of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded under the Australian Research council linkage grant scheme (LP0989708 and LP100200842).

Institutional Review Board Statement

All components of the study and study protocol were performed in line with the principle of the Declaration of Helsinki, and were approved by the Human Research Ethics Committee, University of Technology Sydney (Ethics number: 2014000110; Date of Approval: 2 May 2014).

Informed Consent Statement

All participants provided written informed consent prior to participating in the study, and consented to the publication of their anonymised data.

Data Availability Statement

The data relevant to the current study are available from the corresponding author on reasonable request.

Acknowledgments

The authors acknowledge all volunteers for participating, and also wish to acknowledge Budi Jap for software support for EEG data pre-processing.

Conflicts of Interest

The authors declare no conflict of interest. Further, the funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Dorrian, J.; Baulk, S.D.; Dawson, D. Work hours, workload, sleep and fatigue in Australian Rail Industry employees. Appl. Ergon. 2011, 42, 202–209. [Google Scholar] [CrossRef]
  2. Melamed, S.; Ben-Avi, I.; Luz, J.; Green, M.S. Objective and subjective work monotony: Effects on job satisfaction, psychological distress, and absenteeism in blue-collar workers. J. Appl. Psychol. 1995, 80, 29–42. [Google Scholar] [CrossRef]
  3. Melamed, S.; Ben-Avi, I.; Luz, J.; Green, M.S. Repetitive work, work underload and coronary heart disease risk factors among blue-collar workers—The CORDIS Study. J. Psychosom. Res. 1995, 39, 19–29. [Google Scholar] [CrossRef]
  4. Molleman, E.; van den Beukel, A. Worker flexibility and its perceived contribution to performance: The moderating role of task characteristics. Hum. Factors Ergon. Manuf. Serv. Ind. 2007, 17, 117–135. [Google Scholar] [CrossRef]
  5. Thiffault, P.; Bergeron, J. Fatigue and individual differences in monotonous simulated driving. Personal. Individ. Differ. 2003, 34, 159–176. [Google Scholar] [CrossRef]
  6. Larue, G.S.; Rakotonirainy, A.; Pettitt, A.N. Driving performance impairments due to hypovigilance on monotonous roads. Accid. Anal. Prev. 2011, 43, 2037–2046. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Brown, I.D. Driver Fatigue. Hum. Factors 1994, 36, 298–314. [Google Scholar] [CrossRef] [PubMed]
  8. Horne, J.A.; Reyner, L.A. Sleep related vehicle accidents. BMJ 1995, 310, 565. [Google Scholar] [CrossRef] [Green Version]
  9. Lal, S.K.L.; Craig, A. Driver fatigue: Electroencephalography and psychological assessment. Psychophysiology 2002, 39, 313–321. [Google Scholar] [CrossRef]
  10. Ashiya, M.; Sone, S.; Sato, Y.; Kaga, A. Application of pure electric braking system to electric railcars. In Proceedings of the 6th International Workshop on Advanced Motion Control, Noyoga, Japan, 30 March–1 April 2000; pp. 163–168. [Google Scholar]
  11. Tüchsen, F.; Hannerz, H.; Roepstorff, C.; Krause, N. Stroke among male professional drivers in Denmark, 1994–2003. Occup. Environ. Med. 2006, 63, 456–460. [Google Scholar] [CrossRef] [Green Version]
  12. Rosengren, A.; Anderson, K.; Wilhelmsen, L. Risk of coronary heart disease in middle-aged male bus and tram drivers compared to men in other occupations: A prospective study. Int. J. Epidemiol. 1991, 20, 82–87. [Google Scholar] [CrossRef] [PubMed]
  13. Hedberg, G.E.; Jacobsson, K.A.; Janlert, U.; Langendoen, S. Risk indicators of ischemic heart disease among male professional drivers in Sweden. Scand. J. Work Environ. Health 1993, 19, 326–333. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Bigert, C.; Gustavsson, P.; Hallqvist, J.; Hogstedt, C.; Lewné, M.; Plato, N.; Reuterwall, C.; Schéele, P. Myocardial Infarction among Professional Drivers. Epidemiology 2003, 14, 333–339. [Google Scholar] [CrossRef]
  15. Korelitz, J.J.; Fernandez, A.A.; Uyeda, V.J.; Spivey, G.H.; Browdy, B.L.; Schmidt, R.T. Health habits and risk factors among truck drivers visiting a health booth during a trucker trade show. Am. J. Health Promot. AJHP 1993, 8, 117–123. [Google Scholar] [CrossRef]
  16. Lee, I.M.; Shiroma, E.J.; Lobelo, F.; Puska, P.; Blair, S.N.; Katzmarzyk, P.T. Effect of physical inactivity on major non-communicable diseases worldwide: An analysis of burden of disease and life expectancy. Lancet 2012, 380, 219–229. [Google Scholar] [CrossRef] [Green Version]
  17. Blair, S.N.; Brodney, S. Effects of physical inactivity and obesity on morbidity and mortality: Current evidence and research issues. Med. Sci. Sports Exerc. 1999, 31 (Suppl. 11), S646–S662. [Google Scholar] [CrossRef] [PubMed]
  18. Sullivan, P.W.; Morrato, E.H.; Ghushchyan, V.; Wyatt, H.R.; Hill, J.O. Obesity, Inactivity, and the Prevalence of Diabetes and Diabetes-Related Cardiovascular Comorbidities in the U.S. 2000–2002. Diabetes Care 2005, 28, 1599. [Google Scholar] [CrossRef] [Green Version]
  19. Hakkanen, H.; Summala, H. Driver Sleepiness-Related Problems, Health Status, and Prolonged Driving Among Professional Heavy-Vehicle Drivers. Transp. Hum. Factors 2000, 2, 151–171. [Google Scholar] [CrossRef]
  20. Jovanovic, J.; Batanjac, J.; Jovanovic, M.; Bulat, P.; Torbica, N. Occupational profile and cardiac risks: Mechanisms and implications for professional drivers. Int. J. Occup. Med. Environ. Health 1998, 11, 145–152. [Google Scholar]
  21. Laberge-Nadeau, C.; Dionne, G.; Ekoé, J.M.; Hamet, P.; Desjardins, D.; Messier, S.; Maag, U. Impact of diabetes on crash risks of truck-permit holders and commercial drivers. Diabetes Care 2000, 23, 612. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. McCraty, R.; Atkinson, M.; Tomasino, D.; Bradley, R.T. The Coherent Heart: Heart-Brain Interactions, Psychophysiological coherence, and the Emergence of System-Wide Order. Integral Rev. 2009, 5, 10–115. [Google Scholar]
  23. Thayer, J.F.; Lane, R.D. A model of neurovisceral integration in emotion regulation and dysregulation. J. Affect. Disord. 2000, 61, 201–216. [Google Scholar] [CrossRef] [Green Version]
  24. Thayer, J.F.; Hansen, A.L.; Saus-Rose, E.; Johnsen, B.H. Heart Rate Variability, Prefrontal Neural Function, and Cognitive Performance: The Neurovisceral Integration Perspective on Self-regulation, Adaptation, and Health. Ann. Behav. Med. 2009, 37, 141–153. [Google Scholar] [CrossRef]
  25. Idogawa, K. On the brain wave activity of professional drivers during monotonous work. Behaviormetrika 1991, 18, 23–34. [Google Scholar] [CrossRef] [Green Version]
  26. Lal, S.K.L.; Craig, A.; Boord, P.; Kirkup, L.; Nguyen, H. Development of an algorithm for an EEG-based driver fatigue countermeasure. J. Saf. Res. 2003, 34, 321–328. [Google Scholar] [CrossRef]
  27. Schmidt, E.A.; Schrauf, M.; Simon, M.; Fritzsche, M.; Buchner, A.; Kincses, W.E. Drivers’ misjudgement of vigilance state during prolonged monotonous daytime driving. Accid. Anal. Prev. 2009, 41, 1087–1093. [Google Scholar] [CrossRef]
  28. Zilberg, E.; Xu, Z.M.; Burton, D.; Karrar, M.; Lal, S. Statistical validation of physiological indicators for noninvasive and hybrid driver drowsiness detection system. Afr. J. Inf. Commun. Technol. 2009, 5, 75–83. [Google Scholar]
  29. Jap, B.T.; Lal, S.; Fischer, P. Comparing combinations of EEG activity in train drivers during monotonous driving. Expert Syst. Appl. 2011, 38, 996–1003. [Google Scholar] [CrossRef]
  30. Lees, T.; Chalmers, T.; Burton, D.; Zilberg, E.; Penzel, T.; Lal, S.; Lal, S. Electroencephalography as a predictor of self-report fatigue/sleepiness during monotonous driving in train drivers. Physiol. Meas. 2018, 39, 105012. [Google Scholar] [CrossRef]
  31. Santamaria, J.; Chiappa, K.H. The EEG of Drowsiness; DEMOS Publications: New York, NY, USA, 1987. [Google Scholar]
  32. Finelli, L.A.; Baumann, H.; Borbély, A.A.; Achermann, P. Dual electroencephalogram markers of human sleep homeostasis: Correlation between theta activity in waking and slow-wave activity in sleep. Neuroscience 2000, 101, 523–529. [Google Scholar] [CrossRef]
  33. Lal, S.K.; Craig, A. A critical review of the psychophysiology of driver fatigue. Biol. Psychol. 2001, 55, 173–194. [Google Scholar] [CrossRef]
  34. Kaida, K.; Takahashi, M.; Akerstedt, T.; Nakata, A.; Otsuka, Y.; Haratani, T.; Fukasawa, K. Validation of the Karolinska sleepiness scale against performance and EEG variables. Clin. Neurophysiol. Off. J. Int. Fed. Clin. Neurophysiol. 2006, 117, 1574–1581. [Google Scholar] [CrossRef]
  35. Lal, S.K.; Craig, A. Electroencephalography Activity Associated with Driver Fatigue: Implications for a Fatigue Countermeasure Device. J. Psychophysiol. 2001, 15, 183–189. [Google Scholar] [CrossRef]
  36. Phipps-Nelson, J.; Redman, J.R.; Rajaratnam, S.M. Temporal profile of prolonged, night-time driving performance: Breaks from driving temporarily reduce time-on-task fatigue but not sleepiness. J. Sleep Res. 2011, 20, 404–415. [Google Scholar] [CrossRef]
  37. Zhao, C.; Zhao, M.; Liu, J.; Zheng, C. Electroencephalogram and electrocardiograph assessment of mental fatigue in a driving simulator. Accid. Anal. Prev. 2012, 45, 83–90. [Google Scholar] [CrossRef] [PubMed]
  38. Berntson, G.G.; Thomas Bigger, J.; Eckberg, D.L.; Grossman, P.; Kaufmann, P.G.; Malik, M.; Nagaraja, H.N.; Porges, S.W.; Saul, J.P.; Stone, P.H.; et al. Heart rate variability: Origins, methods, and interpretive caveats. Psychophysiology 1997, 34, 623–648. [Google Scholar] [CrossRef] [PubMed]
  39. Greiser, K.H.; Kluttig, A.; Schumann, B.; Kors, J.A.; Swenne, C.A.; Kuss, O.; Werdan, K.; Haerting, J. Cardiovascular disease, risk factors and heart rate variability in the elderly general population: Design and objectives of the CARdiovascular disease, Living and Ageing in Halle (CARLA) Study. BMC Cardiovasc. Disord. 2005, 5, 33. [Google Scholar] [CrossRef] [Green Version]
  40. Patel, M.; Lal, S.K.L.; Kavanagh, D.; Rossiter, P. Applying neural network analysis on heart rate variability data to assess driver fatigue. Expert Syst. Appl. 2011, 38, 7235–7242. [Google Scholar] [CrossRef]
  41. Rothberg, L.J.; Lees, T.; Clifton-Bligh, R.; Lal, S. Association between Heart Rate Variability Measures and Blood Glucose Levels: Implications for Noninvasive Glucose Monitoring for Diabetes. Diabetes Technol. Ther. 2016, 18, 366–376. [Google Scholar] [CrossRef]
  42. Nguyen, L.; Su, S.; Nguyen, H.T. Effects of hyperglycemia on variability of RR, QT and corrected QT intervals in Type 1 diabetic patients. In Proceedings of the Annual International Conference of the IEEE Engineering in Medicine and Biology Society, Osaka, Japan, 3–7 July 2013; Volume 2013, pp. 1819–1822. [Google Scholar]
  43. Vinik, A.I.; Erbas, T.; Casellini, C.M. Diabetic cardiac autonomic neuropathy, inflammation and cardiovascular disease. J. Diabetes Investig. 2013, 4, 4–18. [Google Scholar] [CrossRef]
  44. Bruyne, M.C.D.; Kors, J.A.; Hoes, A.W.; Klootwijk, P.; Dekker, J.M.; Hofman, A.; van Bemmel, J.H.; Grobbee, D.E. Both Decreased and Increased Heart Rate Variability on the Standard 10-Second Electrocardiogram Predict Cardiac Mortality in the ElderlyThe Rotterdam Study. Am. J. Epidemiol. 1999, 150, 1282–1288. [Google Scholar] [CrossRef] [PubMed]
  45. Thayer, J.F.; Yamamoto, S.S.; Brosschot, J.F. The relationship of autonomic imbalance, heart rate variability and cardiovascular disease risk factors. Int. J. Cardiol. 2010, 141, 122–131. [Google Scholar] [CrossRef] [PubMed]
  46. Chalmers, T.; Maharaj, S.; Lees, T.; Lin, C.T.; Newton, P.; Clifton-Bligh, R.; McLachlan, C.S.; Gustin, S.M.; Lal, S. Impact of acute stress on cortical electrical activity and cardiac autonomic coupling. J. Integr. Neurosci. 2020, 19, 239–248. [Google Scholar] [PubMed]
  47. Craig, A.; Hancock, K.; Craig, M. The lifestyle appraisal questionnaire: A comprehensive assessment of health and stress. Psychol. Health 1996, 11, 331–343. [Google Scholar] [CrossRef]
  48. National Heart Foundation of Australia Guideline for the Diagnosis and Management of Hypertension in Adults—2016. Available online: https://www.heartfoundation.org.au/images/uploads/publications/PRO-167_Hypertension-guideline-2016_WEB.pdf (accessed on 17 November 2020).
  49. Shapiro, D.; Jamner, L.D.; Lane, J.D.; Light, K.C.; Myrtek, M.; Sawada, Y.; Steptoe, A. Blood pressure publication guidelines. Psychophysiology 1996, 33, 1–12. [Google Scholar] [CrossRef]
  50. Tarvainen, M.P.; Lipponen, J.; Niskanen, J.; Ranta-aho, P.O. Kubios HRV User’s Guide; Kubios Oy: Kuopio, Finland, 2018; pp. 1–64. [Google Scholar]
  51. Tarvainen, M.P.; Ranta-aho, P.O.; Karjalainen, P.A. An advanced detrending method with application to HRV analysis. IEEE Trans. Biomed. Eng. 2002, 49, 172–175. [Google Scholar] [CrossRef] [PubMed]
  52. Task Force of the European Society of Cardiology the North American Society of Pacing Electrophysiology, Heart Rate Variability: Standards of Measurement, Physiological Interpretation, and Clinical Use. Circulation 1996, 93, 1043–1065. [CrossRef] [Green Version]
  53. Welch, P. The use of fast Fourier transform for the estimation of power spectra: A method based on time averaging over short, modified periodograms. IEEE Trans. Audio Electroacoust. 1967, 15, 70–73. [Google Scholar] [CrossRef] [Green Version]
  54. Croft, R.J.; Barry, R.J. EOG correction: A new aligned-artifact average solution. Electroencephalogr. Clin. Neurophysiol. 1998, 107, 395–401. [Google Scholar] [CrossRef]
  55. Croft, R.J.; Barry, R.J. EOG correction of blinks with saccade coefficients: A test and revision of the aligned-artefact average solution. Clin. Neurophysiol. 2000, 111, 444–451. [Google Scholar] [CrossRef]
  56. Rowan, A.J.; Tolunsky, E. Primer of EEG: With a Mini-Atlas; Butterworth-Heinemann: Oxford, UK, 2003; p. 224. [Google Scholar]
  57. Leys, C.; Ley, C.; Klein, O.; Bernard, P.; Licata, L. Detecting outliers: Do not use standard deviation around the mean, use absolute deviation around the median. J. Exp. Soc. Psychol. 2013, 49, 764–766. [Google Scholar] [CrossRef] [Green Version]
  58. Vanlaar, W.; Simpson, H.; Mayhew, D.; Robertson, R. Fatigued and drowsy driving: A survey of attitudes, opinions and behaviors. J. Saf. Res. 2008, 39, 303–309. [Google Scholar] [CrossRef]
  59. Wascher, E.; Rasch, B.; Sänger, J.; Hoffmann, S.; Schneider, D.; Rinkenauer, G.; Heuer, H.; Gutberlet, I. Frontal theta activity reflects distinct aspects of mental fatigue. Biol. Psychol. 2014, 96 (Suppl. C), 57–65. [Google Scholar] [CrossRef]
  60. Bear, M.F.; Connors, B.W.; Paradiso, M.A. Neuroscience: Exploring the Brain; Lippincott Williams & Wilkins: Philadelphia, PA, USA, 2007. [Google Scholar]
  61. Lombardi, F.; Sandrone, G.; Pernpruner, S.; Sala, R.; Garimoldi, M.; Cerutti, S.; Baselli, G.; Pagani, M.; Malliani, A. Heart rate variability as an index of sympathovagal interaction after acute myocardial infarction. Am. J. Cardiol. 1987, 60, 1239–1245. [Google Scholar] [CrossRef]
  62. Lombardi, F.; Sandrone, G.; Mortara, A.; La Rovere, M.T.; Colombo, E.; Guzzetti, S.; Malliani, A. Circadian variation of spectral indices of heart rate variability after myocardial infarction. Am. Heart J. 1992, 123, 1521–1529. [Google Scholar] [CrossRef]
  63. Triposkiadis, F.; Karayannis, G.; Giamouzis, G.; Skoularigis, J.; Louridas, G.; Butler, J. The Sympathetic Nervous System in Heart Failure: Physiology, Pathophysiology, and Clinical Implications. J. Am. Coll. Cardiol. 2009, 54, 1747–1762. [Google Scholar] [CrossRef] [Green Version]
  64. Tallon-Baudry, C. The roles of gamma-band oscillatory synchrony in human visual cognition. Front. Biosci. (Landmark Ed.) 2009, 14, 321–332. [Google Scholar] [CrossRef]
  65. Jap, B.T.; Lal, S.; Fischer, P. Inter-hemispheric electroencephalography coherence analysis: Assessing brain activity during monotonous driving. Int. J. Psychophysiol. 2010, 76, 169–173. [Google Scholar] [CrossRef] [Green Version]
  66. Putilov, A.A.; Donskaya, O.G.; Verevkin, E.G. Generalizability of Frequency Weighting Curve for Extraction of Spectral Drowsy Component from the EEG Signals Recorded in Eyes-Closed Condition. Clin. EEG Neurosci. 2017, 48, 259–269. [Google Scholar] [CrossRef] [PubMed]
  67. Jap, B.T.; Lal, S.; Fischer, P.; Bekiaris, E. Using EEG spectral components to assess algorithms for detecting fatigue. Expert Syst. Appl. 2009, 36 Pt 1, 2352–2359. [Google Scholar] [CrossRef]
  68. Kecklund, G.; Åkerstedt, T. Sleepiness in long distance truck driving: An ambulatory EEG study of night driving. Ergonomics 1993, 36, 1007–1017. [Google Scholar] [CrossRef] [PubMed]
  69. Gillberg, M.; Kecklund, G.; Åkerstedt, T. Sleepiness and performance of professional drivers in a truck simulator—Comparisons between day and night driving. J. Sleep Res. 1996, 5, 12–15. [Google Scholar] [CrossRef] [PubMed]
  70. Sandberg, D.; Anund, A.; Fors, C.; Kecklund, G.; Karlsson, J.G.; Wahde, M.; Åkerstedt, T. The Characteristics of Sleepiness During Real Driving at Night—A Study of Driving Performance, Physiology and Subjective Experience. Sleep 2011, 34, 1317–1325. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Table 1. Significant associations between electroencephalography variables and log transformed low frequency heart rate variability of train drivers.
Table 1. Significant associations between electroencephalography variables and log transformed low frequency heart rate variability of train drivers.
Dependent VariableIndependent Variablenrp
LFlnO1 − δ340.430.04
Fp2 − θ46−0.400.03
Fp1 − α45−0.430.02
FT8 − β38−0.58<0.01
P4 − γ390.410.03
O1 − γ370.39<0.05
Key: F = Frontal; Fp = Frontal Pole; LFln = log-transformed low frequency heart rate variability; n = Sample size; O = Occipital; P = Parietal; T = Temporal; α = Alpha; β = Beta; δ = Delta; γ = Gamma; θ = Theta.
Table 2. Regression analysis for log-transformed low frequency heart rate variability and electroencephalography variables in train drivers.
Table 2. Regression analysis for log-transformed low frequency heart rate variability and electroencephalography variables in train drivers.
R = 0.72; R2 = 0.52; AR2 = 0.42; df = 3, 14; F = 5.14; p = 0.01 *
VariableβSE of βBSE of Btp
Intercept 6.720.3420.00<0.001 *
O1 − δ−0.340.19<−0.01<0.01−1.830.08
FT8 − β−0.440.18−0.050.02−2.400.03 *
O1 − γ0.390.190.130.062.090.05
Key: df = Degrees of freedom; F = Frontal; O = Occipital; SE = Standard Error; T = Temporal; β = Beta; δ = Delta; γ = Gamma; * = Statistical significance.
Table 3. Significant associations between electroencephalography variables and normalised low frequency heart rate variability of train drivers.
Table 3. Significant associations between electroencephalography variables and normalised low frequency heart rate variability of train drivers.
Dependent VariableIndependent Variablenrp
LFn.u.C3 − δ420.380.04
Fp1 − α450.430.02
P4 − β36−0.430.02
P4 − γ39−0.430.02
P8 − γ40−0.390.04
O1 − γ37−0.460.02
Key: C = Central; Fp = Frontal Pole; LFn.u. = Normalised low frequency heart rate variability; n = Sample size; O = Occipital; P = Parietal; α = Alpha; β = Beta; δ = Delta; γ = Gamma.
Table 4. Regression analysis for normalised low frequency heart rate variability and electroencephalography variables in train drivers.
Table 4. Regression analysis for normalised low frequency heart rate variability and electroencephalography variables in train drivers.
R = 0.75; R2 = 0.57; AR2 = 0.52; df = 2, 20; F = 13.11; p < 0.001 *
VariableβSE of βBSE of Btp
Intercept 73.612.5628.74<0.001 *
C3 − δ0.310.150.120.062.11<0.05 *
O1 − γ−0.650.15−2.050.46−4.42<0.001 *
Key: C = Central; df = Degrees of freedom; O = Occipital; SE = Standard Error; δ = Delta; γ = Gamma; * = Statistical significance.
Table 5. Significant associations between electroencephalography variables and log transformed high frequency heart rate variability of train drivers.
Table 5. Significant associations between electroencephalography variables and log transformed high frequency heart rate variability of train drivers.
Dependent VariableIndependent Variablenrp
HFlnFT7 − δ46−0.390.04
C3 − δ42−0.36<0.05
Fp1 − α45−0.56<0.01
FT8 − β38−0.410.03
P4 − γ390.470.01
P8 − γ400.38<0.05
O1 − γ370.470.01
Key: C = Central; F = Frontal; Fp = Frontal Pole; HFln = log-transformed high frequency heart rate variability; n = Sample size; O = Occipital; P = Parietal; T = Temporal; α = Alpha; β = Beta; δ = Delta; γ = Gamma.
Table 6. Regression analysis for log-transformed high frequency heart rate variability and electroencephalography variables in train drivers.
Table 6. Regression analysis for log-transformed high frequency heart rate variability and electroencephalography variables in train drivers.
R = 0.79; R2 = 0.62; AR2 = 0.55; df = 3, 17; F = 9.14; p < 0.001 *
VariableβSE of βBSE of Btp
Intercept 5.400.4312.62<0.001 *
C3 − δ−0.170.15−0.010.01−1.130.27
FT8 − β−0.320.15−0.060.03−2.090.05
O1 − γ0.670.150.270.064.40<0.001 *
Key: C = Central; df = Degrees of freedom; F = Frontal; O = Occipital; SE = Standard Error; T = Temporal; β = Beta; δ = Delta; γ = Gamma; * = Statistical significance.
Table 7. Significant associations between electroencephalography variables and normalised high frequency heart rate variability of train drivers.
Table 7. Significant associations between electroencephalography variables and normalised high frequency heart rate variability of train drivers.
Dependent VariableIndependent Variablenrp
HFn.u.C3 − δ42−0.380.04
Fp1 − α45−0.430.02
P4 − β360.430.03
P4 − γ390.430.02
P8 − γ400.380.04
O1 − γ370.460.02
Key: C = Central; Fp = Frontal Pole; HFn.u. = Normalised high frequency heart rate variability; n = Sample size; O = Occipital; P = Parietal; α = Alpha; β = Beta; δ = Delta; γ = Gamma.
Table 8. Regression analysis for normalised high frequency heart rate variability and electroencephalography variables in train drivers.
Table 8. Regression analysis for normalised high frequency heart rate variability and electroencephalography variables in train drivers.
R = 0.75; R2 = 0.56; AR2 = 0.52; df = 2, 20; F = 12.97; p < 0.001 *
VariableβSE of βBSE of Btp
Intercept 26.342.5510.32<0.001 *
C3 − δ−0.310.15−0.120.06−2.11<0.05 *
O1 − γ0.650.152.030.464.38<0.001 *
Key: C = Central; df = Degrees of freedom; O = Occipital; SE = Standard Error; δ = Delta; γ = Gamma; * = Statistical significance.
Table 9. Significant associations between electroencephalography variables and log transformed total heart rate variability power of train drivers.
Table 9. Significant associations between electroencephalography variables and log transformed total heart rate variability power of train drivers.
Dependent VariableIndependent Variablenrp
TPowlnFp1 − δ45−0.380.04
Fp2 − θ46−0.390.03
Fp1 − α45−0.470.01
FT8 − β38−0.52<0.01
P4 − γ390.460.01
O1 − γ370.460.02
Key: F = Frontal; Fp = Frontal Pole; n = Sample size; O = Occipital; P = Parietal; T = Temporal; TPowln = log-transformed total heart rate variability power; α = Alpha; β = Beta; δ = Delta; γ = Gamma; θ = Theta.
Table 10. Regression analysis for log-transformed total heart rate variability power and electroencephalography variables in train drivers.
Table 10. Regression analysis for log-transformed total heart rate variability power and electroencephalography variables in train drivers.
R = 0.81; R2 = 0.65; AR2 = 0.61;df = 2,17; F = 15.98; p < 0.001 *
VariableβSE of βBSE of Btp
Intercept 930.22851.691.090.29
P4 − γ-0.200.14−1830.531345.77−1.360.19
O1 − γ0.810.14875.16155.645.62<0.001 *
Key: dfM = Degrees of freedom; O = Occipital; P = Parietal; SE = Standard Error; γ = Gamma; * = Statistical significance.
Table 11. Significant associations between electroencephalography variables and the low/high frequency heart rate variability ratio of train drivers.
Table 11. Significant associations between electroencephalography variables and the low/high frequency heart rate variability ratio of train drivers.
Dependent VariableIndependent Variablenrp
LF/HFFT7 − δ460.53<0.01
FT8 − δ460.48<0.01
FT7 − θ470.390.04
FC4 − θ400.63<0.001
Fp1 − α450.52<0.01
CPz − β400.490.01
FC3 − γ350.54<0.01
CPz − γ400.450.02
Key: C = Central; F = Frontal; Fp = Frontal Pole; HF = High-frequency heart rate variability; LF = Low-frequency heart rate variability; n = Sample size; O = Occipital; P = Parietal; T = Temporal; z = Midline; α = Alpha; β = Beta; δ = Delta; γ = Gamma; θ = Theta.
Table 12. Regression analysis for the low/high frequency heart rate variability ratio and electroencephalography variables in train drivers.
Table 12. Regression analysis for the low/high frequency heart rate variability ratio and electroencephalography variables in train drivers.
R = 0.78; R2 = 0.61; AR2 = 0.56;df= 2, 17; F = 13.07; p < 0.001 *
VariableβSE of βBSE of Btp
Intercept 1.590.831.920.07
FT8 − δ0.590.220.030.012.610.02 *
FC4 − θ0.240.220.130.121.060.30
Key: C = Central; df = Degrees of freedom; F = Frontal; O = Occipital; T = Temporal; δ = Delta; θ = Theta; * = Statistical significance.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lees, T.; Chalmers, T.; Burton, D.; Zilberg, E.; Penzel, T.; Lal, S.; Lal, S. Electrophysiological Brain-Cardiac Coupling in Train Drivers during Monotonous Driving. Int. J. Environ. Res. Public Health 2021, 18, 3741. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph18073741

AMA Style

Lees T, Chalmers T, Burton D, Zilberg E, Penzel T, Lal S, Lal S. Electrophysiological Brain-Cardiac Coupling in Train Drivers during Monotonous Driving. International Journal of Environmental Research and Public Health. 2021; 18(7):3741. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph18073741

Chicago/Turabian Style

Lees, Ty, Taryn Chalmers, David Burton, Eugene Zilberg, Thomas Penzel, Shail Lal, and Sara Lal. 2021. "Electrophysiological Brain-Cardiac Coupling in Train Drivers during Monotonous Driving" International Journal of Environmental Research and Public Health 18, no. 7: 3741. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph18073741

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop