Effective Connectivity Analysis of the Brain Network in Drivers during Actual Driving Using Near-Infrared Spectroscopy

Introduction

Driving a vehicle is a complex activity that requires high-level brain functions, such as planning, decision making, visual attention, motor control and high cognitive activity to make fast cognitive decisions in a complex and rapidly changing environment (Derosière et al., 2014). Some individuals may be able to drive under simple conditions, but may be incapable of driving safely when cognitive workloads become heavier, such as following a proceeding car, talking on the phone and performing other subtasks (Uchiyama et al., 2003; Schweizer et al., 2013). Several studies have reported that some brain regions, such as the prefrontal cortex (PFC), motor-related areas (MA), parietal cortex, vision-related areas (VA), thalamus and cerebellum, engage in high-level cognitive activity during car driving with functional near-infrared spectroscopy (fNIRS; Yoshino et al., 2013a,b; Orino et al., 2015). However, little information is known on how the functional networks are affected by the cognitive workload during actual driving. Therefore, it makes sense to analyze driving data in a manner that evaluates the brain activation among regions, which enables us to study how the brain is functionally connected and how these intrinsic networks are modulated by the cognitive workload (Calhoun and Pearlson, 2012).

The brain activation was identified with functional magnetic resonance imaging (fMRI) and their modulation with speed was investigated during simulated driving (Calhoun et al., 2002). The authors found that the signal in frontoparietal regions decreased exponentially with a rate proportional to the driving speed. They also found that increases in the cerebella and occipital areas, presumably related to the complex visuomotor integration, were activated during driving but not associated with the driving speed (Calhoun et al., 2002).

The brain activation among cortex areas was found to differ with familiar or unfamiliar routes. For example, significant activation was found in the middle temporal and occipital cortex and in the cerebellum for the unfamiliar route. A training period and a familiar, monotonous route may lead to a reduction in attention and perception (Mader et al., 2009). Furthermore, the brain activity assessed in terms of regional response and regional interactions in highly trained racing-car drivers was found to differ from that of subjects with an ordinary driving experience (Bernardi et al., 2014).

The brain activation during actual car-driving on the road was demonstrated to be similar to that of simulated driving and visual perception and visuomotor coordination were the main brain functions while driving (Jeong et al., 2006). However, autonomic and emotional responses such as attention and autonomic arousal should be considered using actual driving (Jeong et al., 2006).

fNIRS has been proven to be a reliable method to represent the cortical activities of the brain network during driving (Tsunashima and Yanagisawa, 2009; Yoshino et al., 2013b; Lin and Lin, 2016), which monitors spontaneous hemodynamic oscillations from cortical regions. Compared with other neuroimaging methods such as fMRI, fNIRS requires relatively few physical constraints and is suitable for measuring moving subjects. Moreover, fNIRS exhibits a higher sampling rate (approximately 10 Hz) than fMRI (approximately 1 Hz), implying a better time resolution.

Several studies show that functional connectivity (FC) can be assessed using fNIRS data in drivers (Wang W. et al., 2016; Xu et al., 2017). FC is generally inferred by the correlation between nodal activities on the basis of blood oxygenation level-dependent fMRI or coherence in electro- or magneto- encephalogram signals acquired during task performance or the resting state (Park and Friston, 2013). However, FC does not provide the mechanisms of neuronal coupling. Effective connectivity (EC) is defined as the influence that a node exerts over another under a network model of causal dynamics and is inferred from a model of neuronal integration, which defines the mechanisms of neuronal coupling (Friston, 2011).

EC is usually computed using several methods, including dynamic causal modeling (DCM; Friston et al., 2003), Granger causality (GC; Granger, 1969), and structural equation modeling (SEM; Anderson and Gerbing, 1988). The basic idea of GC was proposed by Wiener and was then successively formalized by Granger using an autoregressive model (Granger, 1969; Bressler and Seth, 2011). Subsequently, Geweke introduced the conditional GC (CGC) algorithm to extend GC to a conditional case (Geweke, 1982). Recently, GC has been used to assess the directionality of neuronal interactions in both the frequency and the time domain (Im et al., 2010; Sitaram et al., 2014; Anwar et al., 2016; Cai et al., 2016). Compared with GC, SEM becomes unstable when the data is presented as a time series, while DCM always requires the selection of a prior model in advance, which may prove to be a very complex process because of the presence of multiple variables. The EC networks among overlapping core regions recruited by motor execution and motor imagery were explored by using CGC and on fMRI data. It was demonstrated that more circuits of EC among the selected seed regions were activated during right-hand performance than during left-hand performance (Gao et al., 2011).

This study hypothesizes that driving as a complex activity would activate a brain interaction among different brain areas and strengthen the EC of the brain network compared to the resting state. Given that PFC significantly contributes to cognitive behavior (Miller and Cohen, 2001), MA is highly associated with sensation and movement control (Peterka, 2002), and VA plays an important role in processing visual information (Uchiyama et al., 2003). Understanding the cooperation mechanism of these brain areas during driving can help explain how the cognitive workload influences the brain network. In this study, EC was measured using CGC, which is deemed suitable for multivariate time series data and requires no prior knowledge of the directions of connectivity among different brain areas (Gao et al., 2011). This study aimed to achieve the following goals: (1) to investigate the activation status between actual driving and the resting state; and (2) to assess the EC among different brain areas and examine how cognitive workload affects such connectivity.

Materials and Methods

Subjects

Twelve young male subjects were recruited from the university for this study (age: 24.5 ± 2.8 years, BMI: 23.6 ± 2.3). The subjects were not under psychotropic medication (e.g., stimulants, anti-depressants and anxiolytics) and had no history of neurological injury or disease, seizure disorder, or psychiatric diagnosis based on an interview. All subjects were right-handed according to the Edinburgh Handedness Inventory (Oldfield, 1971) and were clear of the following characteristics: hypertension, neurological or psychiatric diseases, smoking or drinking habits and abnormal heart, lung, and kidney functions. All subjects were paid for their participation. The experimental procedure was approved by the Human Ethics Committee of the National Research Center for Rehabilitation Technical Aids and was conducted in accordance with the ethical standards specified by the Helsinki Declaration of 1975 (revised in 2008). Written informed consent was obtained from all subjects before study enrollment.

Experimental Procedure and Measurements

This experiment was performed on an annular route in the Xinglong Mountain Campus of Shandong University, Jinan, China (Figure 1). This experimental route is approximately 1400 m long and covers five corners (A-E), where the A-B-C section is slightly downhill, the C-D-E section is flat, and the E-A section is slightly uphill.

FIGURE 1

Figure 1. Experiment location. The driving route is cycling as A-B-C-D-E-A.

A Toyota Yaris-2011 subcompact car and a Toyota E’Z-2011 FUV were used for this experiment. The former was used as the test car and the latter was used as the leading car. For the test car, a 2000W inverter was used to convert the 12V DC of the vehicle battery to 220V AC for the power supply of the NIRS instrument and laptop. A multi-channel tissue oxygenation monitor (NirSan Danyang Huichuang Medical Equipment Co. Ltd.) was utilized to record the delta [HbO₂] of the subjects. The optodes were positioned over left PFC (l-PFC: FP1 and AF3), right PFC (r-PFC: FP2 and AF4), left MA (l-MA: FC1, FC3, C1 and C3), right MA (r-MA: FC2, FC4, C2 and C4), left VA (l-VA: O1 and PO3), and right VA (r-VA: O2 and PO4) with 16 channels in accordance with the International 10/10 System (Table 1, Figure 2). The sampling rate was set to 10 Hz.

TABLE 1

Table 1. Distribution of the channels in the six brain areas.

FIGURE 2

Figure 2. Optode placement of 16 channels. Red dots denote the source optodes, green dots denote the detector optodes, and blue dashed circles represent the corresponding brain areas. The measurement channels (black circles) in the prefrontal cortex (PFC), motor-related areas (MA) and vision-related areas (VA) are in accordance with the international 10/10 system.

The experiment was divided into three sessions, namely, resting, task_1 (simple-driving) and task_2 (car-following) states. Each session lasted 5 min following a 5-min rest after each session. The subjects were instructed to familiarize themselves with the protocol for performing the tasks prior to the experiment. In the resting state, the subjects were instructed to remain still and relaxed with their eyes open, hands on the steering wheel and feet on the pedals while remaining seated on the driver’s seat. In the task_1 state, the subjects were required to drive in the right lane within a speed of around 25 km/h. In the task_2 state, the driver in the leading car was asked to drive at a constant speed of 25 km/h. The subjects in the test car were asked to drive at a safe distance from the leading car and maintain the distance as far as possible. In the resting state, the engine was turned on as air conditioning was needed to maintain a comfortable temperature in the car. After each session, subjects were asked to complete the NASA Task Load Index (NASA-TLX), whose rating range was from 0 (no cognitive workload) to 100 (heavy cognitive workload), to rate the subjective cognitive workload (Hart and Staveland, 1988). Two persons (the subject and the recorder) were in the car during the entire experiment session.

Data Processing

Pre-Processing

The data pre-processing and wavelet-based coherence analysis methods have been described in our previous studies (Bu et al., 2016; Tan et al., 2016; Wang B. et al., 2016). In the present study, first, a five-order Butterworth band-pass filter was used: the frequency elements below 0.021 Hz, which were mainly physiological noises of endothelial activities, and the elements above 0.145 Hz, which mainly reflected the respiration and cardiac activities, were removed (Shiogai et al., 2010). Then the moving average method was used to remove the outlier terms in the delta [HbO₂] signals.

Wavelet Amplitude

Wavelet transform is a common method for transforming time series data from the time domain to the time-frequency domain with appropriate time and frequency resolution. The Morlet wavelet, which is a Gaussian function modulated by a sine wave with basic frequency, was used as a mother wavelet to detect the frequency content. The wavelet amplitude (WA) was the time-averaged wavelet transform amplitude of delta [HbO₂], which can also indicate the frequency properties over the time domain.

Conditional Granger Causal Analysis

The following procedure was mainly accomplished with the GCCA toolbox (Seth, 2010). The basic idea of GC is that for time series X and Y, if knowing the past information of X helps to predict the future Y, X is assumed to “cause” Y, which is shown as F_X→Y. CGC is an extension of GC in multivariate autoregressive models by “conditioning out” the influences of the other time series (Barnett and Seth, 2014; Wang L. et al., 2016). For time series X and Y chosen as the source and the target, respectively, all the other time series are composed of Z values, shown as F_X→Y|Z. Each step of this method is described in detail as follows.

The data was first detrended and demeaned; then, the Augmented-Dickey-Fuller test and the Kwiatkowski-Phillips-Schmidt-Shin test were performed to check the covariance stationarity (Chatfield and Fuller, 1977; Kwiatkowski et al., 1990). The model order was identified using the Akaike information criterion and Bayesian information criterion (Akaike, 1974; Schwarz, 1978).

Second, a multivariate autoregression model was established to calculate the CGC of the data. The GC values of each connection between two channels were calculated under a Bonferroni-corrected significance threshold of p = 0.01 (Seth, 2010).

Three tests were performed to test the validity of the proposed model: the Durbin-Watson test checked whether the residuals were serially uncorrelated (Durbin and Watson, 1950), the consistency test was applied to assess the portion of data captured by this model (Ding et al., 2000), and a third test was performed as a supplementary test for consistency to analyze the adjusted sum square error (Seth, 2010).

The existence and the strength of a connection must be ensured to construct a connectivity graph. On the one hand, for a single subject, the existence of one connection between two channels was determined by the GC values under a Bonferroni-corrected significance threshold of p = 0.01. At group level, if the connections between two channels exist over 75% of subjects (Park and Friston, 2013), it was considered to have a connection between these two channels. In the connectivity graph, if the connections between two brain areas exist over 75% of connections between any two channels from these two brain areas, a connection was considered to exist a connection between these two brain areas in this study. On the other hand, the strength of each connection is defined as the average GC values between two brain areas. The GC values of the subjects in the three states were then averaged to establish three matrices, whose elements denoted the average GC values from one brain area to another.

Causal Flow

One additional index, causal flow, was applied after establishing an EC graph. The causal flow of a brain area, which is denoted as node i in a connectivity graph, is defined as the difference between its in-degree and its out-degree (Seth, 2009). A brain area with a positive causal flow exerts a strong causal influence on a causal network and can be referred to as a causal source. By contrast, a node with negative causal flow can be referred to as a causal target.

Statistical Analysis

The Kolmogorov-Smirnov test and the Levene test were performed in this study to ensure that the obtained values met the assumption required by the analysis of variance (ANOVA) analysis. One-way repeated ANOVA was performed to assess the main differences in the NASA-TLX scores, WA, matrices of average GC values, and causal flows among the resting, task_1 and task_2 states. Then the post hoc test was performed. The statistical significance was set to p < 0.05 with Bonferroni-correction for pair-wise comparisons.

Results

Cognitive Workload

Figure 3 shows the comparison of the cognitive workload among three states through NASA-TLX. The NASA-TLX scores significantly increased as the states changed from the resting state to the task_2 state. Significant changes were observed among the three states (resting-task_1-task_2: p < 0.05, F_(2,33) = 45.3305; resting-task_1: p < 0.05, F_(1,22) = 55.5834; resting-task_2: p < 0.05, F_(1,22) = 92.2196; and task_1-task_2: p = 0.0428, F_(1,22) = 4.6235).

FIGURE 3

Figure 3. Comparison of the NASA-task load index (TLX) among three states. “*” Indicates a statistically significant change.

Wavelet Amplitude

Figure 4 gives an example of the delta [HbO₂] signals of one subject after pre-processing. As shown in the picture, for most channels, especially the channels in VA, the fluctuations of delta [HbO₂] signals in the task_2 state were larger than in the other two states. Figure 5 shows a comparison of the WA of the subjects in six brain areas among the three states. The WA of right MA and VA was higher than that of the left part of the brain in MA and VA. The WA increased with an increase in the cognitive workload. According to the ANOVA test results, statistically significant differences were observed in bilateral PFC among the three states (resting-task_1-task_2; r-PFC: p = 0.0071, F_(2,33) = 5.7682; l-PFC: p = 0.0298, F_(2,33) = 3.9172), resting-task_1 state (r-PFC: p = 0.0072, F_(1,22) = 8.7867; l-PFC: p = 0.0157, F_(1,22) = 6.8510), and resting-task_2 state (r-PFC: p = 0.0028, F_(1,22) = 11.3614 ; l-PFC: p = 0.0102, F_(1,22) = 7.8985).

FIGURE 4

Figure 4. Delta [HbO₂] signals of one subject. The distribution of the channels is the same as that in Table 1 and Figure 2. For each channel, the horizontal X-axis shows the time and the vertical Y-axis shows the delta [HbO₂] signals. The signals of different colors represent different delta [HbO₂] signals in different states.

FIGURE 5

Figure 5. Comparison of the WA in the six brain areas among the three states. “*” Indicates a statistically significant change. Abbreviations: r-PFC, right prefrontal cortex; l-PFC, left prefrontal cortex; r-MA, right motor-related areas; l-MA, right motor-related areas; r-VA, right vision-related areas; l-VA, left vision-related areas; WA, wavelet amplitude.

EC Graphs

Figures 6A–C presents the GC connectivity graphs of six brain areas in the three states, respectively. Figures 7A–C presents the average GC value matrices of six brain areas in the three states, respectively. Table 2 compares the connection strength among the three states.

FIGURE 6

Figure 6. EC graphs of six brain areas in the (A) (resting), (B) (task_1) and (C) (task_2) states, respectively. The connection with double-headed arrows denotes the connection in EC is bidirectional. The connection with single head arrow denotes the connection in EC is unidirectional. Different colors represent different connection strength. The colder color represents the stronger connection strength. Abbreviations: r-PFC, right prefrontal cortex; l-PFC, left prefrontal cortex; r-MA, right motor-related areas; l-MA, right motor-related areas; r-VA, right vision-related areas; l-VA, left vision-related areas; EC, effective connectivity.

FIGURE 7

Figure 7. Average GC value matrices of six brain areas in the (A) (resting), (B) (task_1) and (C) (task_2) states, respectively. Each element in the matrix denotes the average GC values across all subjects. The causal line is from bottom row to the left column. The darker color of the block represents the higher GC values of the connection. Abbreviations: r-PFC, right prefrontal cortex; l-PFC, left prefrontal cortex; r-MA, right motor-related areas; l-MA, right motor-related areas; r-VA, right vision-related areas; l-VA, left vision-related areas; GC, granger causality.

TABLE 2

Table 2. Comparison in the connection strength of the subjects among three states (A; resting-task_1-task_2), (B; resting-task_1), (C; resting- task_2) and (D; task_1-task_2).

There were more connections among brain areas in the task states than in the resting state. The connections from bilateral PFC to bilateral VA, from bilateral VA to l-MA, and from l-MA to r-VA “appeared” from the resting state to the task_1 state. The connection from r-MA to r-VA was “lost” from the task_1 state to the task_2 state. For the connection strength, among the three states (resting-task_1-task_2), the connection strength of connections from l-PFC to r-PFC (p = 0.0184, F_(2,33) = 4.5173), from r-VA to r-PFC (p = 0.0375, F_(2,33) = 3.6327), and l-PFC (p = 0.0423, F_(2,33) = 3.4863) showed significant difference. From the resting state to the task_1 state, the connection strength of the connection from l-PFC to r-PFC increased significantly (p = 0.0158, F_(1,22) = 6.8365). From the task_1 state to the task_2 state, the connection strength of the connection from l-PFC to r-PFC (p = 0.0137, F_(1,22) = 7.1812), from r-VA to r-PFC (p = 0.0152, F_(1,22) = 6.9311) and l-PFC (p = 0.0109, F_(1,22) = 7.7326) decreased significantly.

Causal Flow

The granger causal flows are shown in Figure 8. In the resting state, the bilateral PFC can be considered as causal targets while bilateral MA and VA can be considered as causal sources. In the task_1 state, the status remained the same. In the task_2 state, the l-MA became a causal target and the in-out degree of l-MA decreased significantly from the task_1 state to the task_2 state (p = 0.0114, F_(1,22) = 7.6337).

FIGURE 8

Figure 8. Comparison of the causal flows in the six brain areas among the three states. “*” Indicates a statistically significant change. Abbreviations: r-PFC, right prefrontal cortex; l-PFC, left prefrontal cortex; r-MA, right motor-related areas; l-MA, right motor-related areas; r-VA, right vision-related areas; l-VA, left vision-related areas.

Discussion

The changes in the EC of the brain network among PFC, MA and VA in the resting, simple-driving and car-following states were assessed in this study using CGC. Several interesting findings were obtained: (1) the WA changes showed the same trend as that of the cognitive workload from the resting state to the task states. The WA of left MA and VA was higher than that of the right brain areas, and the WA of the bilateral PFC changed statistically significantly along with the change in the cognitive workload; (2) the connection strength of EC increased in the task_1 state as compared with that in the resting state, whereas it decreased in the task_2 state as compared with that in the task_1 state; and (3) the bilateral PFC were causal targets, and the bilateral VA and MA were causal sources except for l-MA became a causal target in the task_2 state. These findings suggest that the EC of the brain network can be strengthened by a cognitive workload, and can also be weakened by a superfluous cognitive workload such as driving with a car-following task.

The typical hemodynamic response to brain activation is the basis for NIRS measurement. When a specific brain area is activated, neural metabolism is supported through a localized vascular response that causes an influx of oxygen-rich blood to the active area and the surrounding tissue. This phenomenon leads to an increase in [HbO₂] and a decrease in [HbR] in the active brain area (Matthews et al., 2008). The functional hyperemia mechanism adjusts the distribution of cerebral blood flow on the basis of the functional activities of different brain regions (Iadecola, 2004). Therefore, when the activity level of the cerebral cortex region increases or decreases, the blood flow in this area changes accordingly, and this change can be reflected by the fNIRS signals ([HbO₂]).

It was considered that the cerebral NIRS signals originated from neurovascular coupling and systemic activity components (Holper et al., 2014). The sympathetic nervous system and vascular myogenic responses could play a part in neurovascular coupling (Schroeter et al., 2004; Hamner et al., 2010). The signals originate from the intrinsic myogenic activity of smooth muscle cells in resistance vessels and are regulated by the neurogenic activity in the vessel wall (Shiogai et al., 2010). A change in the WA reflects a change in the brain hemodynamic activity level in the neurovascular systemic activity components among different cognitive driving workloads.

Previous studies have highlighted the important role that PFC plays in executive attentional control, eye movement generation, cognitive control and other abstract reasoning functions (Miller and Cohen, 2001). MA have been correlated with the planning, control, and execution of body movements (Zhang et al., 2010). VA, as the visual processing center, receive visual information stemming from the eyes and elicits visual stimuli as object recognition and visuospatial guidance to the other cortex areas (Uchiyama et al., 2003). The WA of bilateral PFC changed significantly with the change of the cognitive workload indicated by the NASA-TLX scores. The increases in the WA were statistically significant from the resting state to the task states, while the increase between the task_1 and the task_2 states did not exhibit such change. This finding may be attributed to the fact that PFC was more sensitive to the change of cognitive work state. The car-following task needed more visuomotor control than the simple-driving task, which was related more to VA than to PFC and MA (Uchiyama et al., 2008). The WA of VA was clearly higher than that of PFC and MA. The significant right lateralization of the WA in VA was possibly attributed to the fact that the left hemisphere directs attention to the right space, but the right hemisphere directs attention to both spaces (Heilman and Van Den Abell, 1980; Kashiwagi et al., 1990). The right lateralization of the WA in MA was probably due to the fact that all subjects were right-handed, and therefore manipulation with their left hand was considerably harder than that with their right hand. Furthermore, according to Chinese laws, the vehicle is driven in the right lane, thus the visual stimuli come mainly from the left side of the view. This may explain the significant right lateralization of the WA in VA.

The directions and the connection strength of the connections of EC can reflect the influence of one brain region on another (Park and Friston, 2013). Some studies have revealed that movement impulses generated by the supplementary motor and premotor cortexes are induced by two sources: one is the occipital-parietal lobes of the posterior attention system, which correlates with the visual-spatial orientation and integration functions; and the other is PFC from the anterior attention system, which serves a relatively higher-level attention function such as executive attentional control in more complex cognitive tasks associated with problem-solving and decision-making (Posner, 2012). For the visual motion information, it is considered to be generated by VA, and then follows two main streams: one follows the dorsal stream to the parietal lobe, and the other follows the ventral stream to the temporal lobe, then the two streams eventually converge in PFC. PFC generates the stimuli for eye movements to obtain visual information (Schiller and Chou, 1998; Rizzolatti and Matelli, 2003).

In this article, the existence of a group-level effective connection between two areas was ensured by a threshold of 75% of subjects. As shown on the constructed EC graphs in this study, the brain activity information transitions among PFC, MA and VA were well connected in the task states. Bilateral PFC to bilateral VA connections “appeared” as the need for visual information increased from the resting state to the task_1 state. This is the same reason that the connection “appeared” from l-MA to r-VA. The connections from bilateral VA to l-MA “appeared” as driving manipulations were mainly implemented by the right hand under the guidance of visual information. During the task_2 state, the connection from r-MA to r-VA became “lost”, which implies that for the car-following task, the feedback neural information transition from r-MA to r-VA decreased, and the influence exertion from r-VA to l-MA relatively increased.

The connection strength from r-VA to bilateral PFC decreased significantly from the task_1 state to the task_2 state. When performing the car-following task, the subjects would concentrate more on the view in front of the car and less on the view on the left and right sides, resulting in a decrease in the influence strength from r-VA to bilateral PFC. Furthermore, the connection strength from l-PFC to r-PFC increased from the resting state to the task_1 state, whereas it decreased during the task_2 state compared with that in the task_1 state. The above mentioned increase implies that driving as an intense cognitive behavior can help strengthen the connection between the bilateral PFC for attention control, planning, decision making and other functions (Yoshino et al., 2013b). The car-following task imposed a larger cognitive workload on subjects than the simple-driving task. This increase in the cognitive level led to a decrease in the coordination of bilateral PFC for attention maintenance and cognitive control.

The positive in-out degree implies that brain area is a causal source in the EC, and the negative degree implies that the brain area served as a causal target. From the resting state to the task states, the in-out degree of PFC exhibited an increasing trend. This trend implies that with an increase in the cognitive workload, the activity of PFC for cognitive function also increases. Additionally, the increase corresponded with the decreasing connection strength from l-PFC to r-PFC. Furthermore, left dorsolateral PFC is necessary for manipulating information in the working memory and right dorsolateral PFC is critical for manipulating information in a broader range of reasoning contexts (Barbey et al., 2013). The in-out degree asymmetry between r-PFC and l-PFC was consistent with the different levels of need for spatial reasoning and visuomotor control for different cognitive workload behaviors. As VA execute the function of receiving visual information and eliciting visual stimuli to other areas (Uchiyama et al., 2003), they always presented as causal sources, and the in-out degree decrease of r-VA matched the decrease in the connection strength from r-VA to bilateral PFC. From the task_1 state to the task_2 state, the in-out degree of l-MA, which controlled the movement function of the contralateral part of the body, exhibited a statistically significant decrease. In the task_2 state, the right foot performed more manipulations for the throttle and brake than in the task_1 state, and l-MA was influenced considerably by the visual stimuli.

Conclusion

In this study, the EC was calculated and analyzed using CGC for actual driving situations. The hemodynamic activity level of brain areas showed an increasing trend with the increase in the cognitive workload, particularly in bilateral PFC. The connection strength of EC was enhanced from the resting state to the simple-driving state, but it deteriorated during the car-following task. The PFC in EC appeared as the causal targets while MA and VA appeared as the causal sources. However, with a subtask of car-following, l-MA turned into causal targets. These findings indicate that the hemodynamic activity level of brain increases linearly with an increase in the cognitive workload. Furthermore, the EC of the brain network can be strengthened by a cognitive workload, and can also be weakened by a superfluous cognitive workload such as driving with a car-following task. This study proposes a method for evaluating the effects of the cognitive workload of drivers.

Limitations of The Study

One limitation of this study was the interference of Mayer waves (0.08–0.1 Hz). Mayer waves are oscillations of arterial pressure occurring spontaneously in conscious subjects at a frequency lower than that of respiration and are tightly coupled with the synchronized oscillations of the efferent sympathetic nervous activity (Julien, 2006). As physiological noise induced by the global systemic physiological processes in fNIRS measurement, Mayer waves may affect the EC. Although the task-related changes in connectivity mainly reflects the synchronization of neurovascular coupling, the interference of Mayer waves should be taken into account in a future study. In addition, the number of subjects included in the present work was relatively limited in light of the current standards for neuroimaging experiments (Friston, 2012). Third, future research will focus on activation level in the standard space using NIRS-SPM (Ye et al., 2009).

Author Contributions

ZLi and QY conceived and designed the study. ZLiu did the experiment, analyzed the data and drafted the manuscript. GX did the experiment and analyzed the data. CH performed the statistical analysis. QT designed the driving task. MZ contributed to the physiological interpretation of the results and edited the manuscript.

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.

Acknowledgments

This project was supported by the National Natural Science Foundation of China (Grant No. 31771071, 11732015) and the Fundamental Research Funds for Central Public Welfare Research Institutes (118009001000160001).

References

Akaike, H. (1974). A new look at the statistical model identification. IEEE. Trans. Automat. Contr. 19, 716–723. doi: 10.1109/TAC.1974.1100705

CrossRef Full Text | Google Scholar

Anderson, J. C., and Gerbing, D. W. (1988). Structural equation modeling in practice: a review and recommended two-step approach. Psychol. Bull. 103, 411–423. doi: 10.1037//0033-2909.103.3.411

CrossRef Full Text | Google Scholar

Anwar, A. R., Muthalib, M., Perrey, S., Galka, A., Granert, O., Wolff, S., et al. (2016). Effective connectivity of cortical sensorimotor networks during finger movement tasks: a simultaneous fNIRS, fMRI, EEG study. Brain Topogr. 29, 645–660. doi: 10.1007/s10548-016-0507-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Barbey, A. K., Koenigs, M., and Grafman, J. (2013). Dorsolateral prefrontal contributions to human working memory. Cortex 49, 1195–1205. doi: 10.1016/j.cortex.2012.05.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Barnett, L., and Seth, A. K. (2014). The MVGC multivariate Granger causality toolbox: a new approach to Granger-causal inference. J. Neurosci. Methods 223, 50–68. doi: 10.1016/j.jneumeth.2013.10.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Bernardi, G., Cecchetti, L., Handjaras, G., Sani, L., Gaglianese, A., Ceccarelli, R., et al. (2014). It’s not all in your car: functional and structural correlates of exceptional driving skills in professional racers. Front. Hum. Neurosci. 8:888. doi: 10.3389/fnhum.2014.00888

PubMed Abstract | CrossRef Full Text | Google Scholar

Bressler, S. L., and Seth, A. K. (2011). Wiener-Granger causality: a well established methodology. Neuroimage 58, 323–329. doi: 10.1016/j.neuroimage.2010.02.059

PubMed Abstract | CrossRef Full Text | Google Scholar

Bu, L., Li, J., Li, F., Liu, H., and Li, Z. (2016). Wavelet coherence analysis of cerebral oxygenation signals measured by near-infrared spectroscopy in sailors: an exploratory, experimental study. BMJ Open 6:e013357. doi: 10.1136/bmjopen-2016-013357

PubMed Abstract | CrossRef Full Text | Google Scholar

Cai, W., Chen, T., Ryali, S., Kochalka, J., Li, C. S., and Menon, V. (2016). Causal interactions within a frontal-cingulate-parietal network during cognitive control: convergent evidence from a multisite-multitask investigation. Cereb. Cortex 26, 2140–2153. doi: 10.1093/cercor/bhv046

PubMed Abstract | CrossRef Full Text | Google Scholar

Calhoun, V. D., and Pearlson, G. D. (2012). A selective review of simulated driving studies: combining naturalistic and hybrid paradigms, analysis approaches, and future directions. Neuroimage 59, 25–35. doi: 10.1016/j.neuroimage.2011.06.037

PubMed Abstract | CrossRef Full Text | Google Scholar

Calhoun, V. D., Pekar, J. J., McGinty, V. B., Adali, T., Watson, T. D., and Pearlson, G. D. (2002). Different activation dynamics in multiple neural systems during simulated driving. Hum. Brain Mapp. 16, 158–167. doi: 10.1002/hbm.10032

PubMed Abstract | CrossRef Full Text | Google Scholar

Chatfield, C., and Fuller, W. A. (1977). Introduction to statistical time series. J. R. Stat. Soc. 140, 379–380.

Derosière, G., Dalhoumi, S., Perrey, S., Dray, G., and Ward, T. (2014). Towards a near infrared spectroscopy-based estimation of operator attentional state. PLoS One 9:e92045. doi: 10.1371/journal.pone.0092045

PubMed Abstract | CrossRef Full Text | Google Scholar

Ding, M., Bressler, S. L., Yang, W., and Liang, H. (2000). Short-window spectral analysis of cortical event-related potentials by adaptive multivariate autoregressive modeling: data preprocessing, model validation, and variability assessment. Biol. Cybern. 83, 35–45. doi: 10.1007/s004229900137

PubMed Abstract | CrossRef Full Text | Google Scholar

Durbin, J., and Watson, G. S. (1950). Testing for serial correlation in least squares regression: I. Biometrika 37, 409–428. doi: 10.2307/2332391

PubMed Abstract | CrossRef Full Text | Google Scholar

Friston, K. J. (2011). Functional and effective connectivity: a review. Brain Connect. 1, 13–36. doi: 10.1089/brain.2011.0008

PubMed Abstract | CrossRef Full Text | Google Scholar

Friston, K. (2012). Ten ironic rules for non-statistical reviewers. Neuroimage 61, 1300–1310. doi: 10.1016/j.neuroimage.2012.04.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Friston, K. J., Harrison, L., and Penny, W. D. (2003). Dynamic causal modeling. Neuroimage 19, 1273–1302. doi: 10.1016/S1053-8119(03)00202-7

CrossRef Full Text | Google Scholar

Gao, Q., Duan, X. J., and Chen, H. F. (2011). Evaluation of effective connectivity of motor areas during motor imagery and execution using conditional Granger causality. Neuroimage 54, 1280–1288. doi: 10.1016/j.neuroimage.2010.08.071

PubMed Abstract | CrossRef Full Text | Google Scholar

Geweke, J. (1982). Measurement of linear dependence and feedback between multiple time series. J. Am. Stat. Assoc. 77, 315–316. doi: 10.2307/2287240

CrossRef Full Text

Granger, C. W. J. (1969). Investigating causal relations by econometric models and cross-spectral methods. Econometrica 37, 424–438. doi: 10.2307/1912791

CrossRef Full Text | Google Scholar

Hamner, J. W., Tan, C. O., Lee, K., Cohen, M. A., and Taylor, J. A. (2010). Sympathetic control of the cerebral vasculature in humans. Stroke 41, 102–109. doi: 10.1161/strokeaha.109.557132

PubMed Abstract | CrossRef Full Text | Google Scholar

Hart, S. G., and Staveland, L. E. (1988). Development of NASA-TLX (Task Load Index): results of empirical and theoretical research. Adv. Psychol. 52, 139–183. doi: 10.1016/s0166-4115(08)62386-9

CrossRef Full Text | Google Scholar

Heilman, K. M., and Van Den Abell, T. (1980). Right-hemisphere dominance for attention: mechanism underlying hemispheric asymmetries of inattention (neglect). Neurology 30, 327–330. doi: 10.1212/wnl.30.3.327

PubMed Abstract | CrossRef Full Text | Google Scholar

Holper, L., Scholkmann, F., and Wolf, M. (2014). The relationship between sympathetic nervous activity and cerebral hemodynamics and oxygenation: a study using skin conductance measurement and functional near-infrared spectroscopy. Behav. Brain Res. 270, 95–107. doi: 10.1016/j.bbr.2014.04.056

PubMed Abstract | CrossRef Full Text | Google Scholar

Iadecola, C. (2004). Neurovascular regulation in the normal brain and in Alzheimer’s disease. Nat. Rev.Neurosci. 5, 347–360. doi: 10.1038/nrn1387

PubMed Abstract | CrossRef Full Text | Google Scholar

Im, C. H., Jung, Y. J., Lee, S., Koh, D., Kim, D. W., and Kim, B. M. (2010). Estimation of directional coupling between cortical areas using Near-Infrared Spectroscopy (NIRS). Opt. Express 18, 5730–5739. doi: 10.1364/OE.18.005730

PubMed Abstract | CrossRef Full Text | Google Scholar

Jeong, M., Tashiro, M., Singh, L. N., Yamaguchi, K., Horikawa, E., Miyake, M., et al. (2006). Functional brain mapping of actual car-driving using [¹⁸F]FDG-PET. Ann. Nucl. Med. 20, 623–628. doi: 10.1007/bf02984660

PubMed Abstract | CrossRef Full Text | Google Scholar

Julien, C. (2006). The enigma of Mayer waves: facts and models. Cardiovasc. Res. 70, 12–21. doi: 10.1016/j.cardiores.2005.11.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Kashiwagi, A., Kashiwagi, T., Nishikawa, T., Tanabe, H., and Okuda, J. I. (1990). Hemispatial neglect in a patient with callosal infarction. Brain 113, 1005–1023. doi: 10.1093/brain/113.4.1005

PubMed Abstract | CrossRef Full Text | Google Scholar

Kwiatkowski, D., Phillips, P. C. B., Schmidt, P., and Shin, Y. (1990). Testing the null hypothesis of stationarity against the alternative of a unit root: how sure are we that economic time series have a unit root? Papers 54, 159–178. doi: 10.1016/0304-4076(92)90104-y

CrossRef Full Text | Google Scholar

Lin, M. I. B., and Lin, K. H. (2016). Walking while performing working memory tasks changes the prefrontal cortex hemodynamic activations and gait kinematics. Front. Behav. Neurosci. 10:92. doi: 10.3389/fnbeh.2016.00092

PubMed Abstract | CrossRef Full Text | Google Scholar

Mader, M., Bresges, A., Topal, R., Busse, A., Forsting, M., and Gizewski, E. R. (2009). Simulated car driving in fMRI—Cerebral activation patterns driving an unfamiliar and a familiar route. Neurosci. Lett. 464, 222–227. doi: 10.1016/j.neulet.2009.08.056

PubMed Abstract | CrossRef Full Text | Google Scholar

Matthews, F., Pearlmutter, B. A., Wards, T. E., Soraghan, C., and Markham, C. (2008). Hemodynamics for brain-computer interfaces. IEEE Signal Process. Mag. 25, 87–94. doi: 10.1109/MSP.2008.4408445

CrossRef Full Text | Google Scholar

Miller, E. K., and Cohen, J. D. (2001). An integrative theory of prefrontal cortex function. Annu. Rev. Neurosci. 24, 167–202. doi: 10.1146/annurev.neuro.24.1.167

PubMed Abstract | CrossRef Full Text | Google Scholar

Oldfield, R. C. (1971). The assessment and analysis of handedness: the Edinburgh inventory. Neuropsychologia 9, 97–113. doi: 10.1016/0028-3932(71)90067-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Orino, Y., Yoshino, K., Oka, N., Yamamoto, K., Takahashi, H., and Kato, T. (2015). Brain activity involved in vehicle velocity changes in a sag vertical curve on an expressway. Transp. Res. Rec. J. Transp. Res. Board 2518, 18–26. doi: 10.3141/2518-03

CrossRef Full Text | Google Scholar

Park, H. J., and Friston, K. (2013). Structural and functional brain networks: from connections to cognition. Science 342:1238411. doi: 10.1126/science.1238411

PubMed Abstract | CrossRef Full Text | Google Scholar

Peterka, R. (2002). Sensorimotor integration in human postural control. J. Neurophysiol. 88, 1097–1118. doi: 10.1152/jn.00605.2001

PubMed Abstract | CrossRef Full Text | Google Scholar

Posner, M. I. (2012). Imaging attention networks. Neuroimage 61, 450–456. doi: 10.1016/j.neuroimage.2011.12.040

PubMed Abstract | CrossRef Full Text | Google Scholar

Rizzolatti, G., and Matelli, M. (2003). Two different streams form the dorsal visual system: anatomy and functions. Exp. Brain Res. 153, 146–157. doi: 10.1007/s00221-003-1588-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Schiller, P. H., and Chou, I. H. (1998). The effects of frontal eye field and dorsomedial frontal cortex lesions on visually guided eye movements. Nat. Neurosci. 1, 248–253. doi: 10.1038/693

PubMed Abstract | CrossRef Full Text | Google Scholar

Schroeter, M. L., Schmiedel, O., and von Cramon, D. Y. (2004). Spontaneous low-frequency oscillations decline in the aging brain. J. Cereb. Blood Flow Metab. 24, 1183–1191. doi: 10.1097/01.wcb.0000135231.90164.40

PubMed Abstract | CrossRef Full Text | Google Scholar

Schwarz, G. (1978). Estimating the dimension of a model. Ann. Stat. 6, 461–464. doi: 10.1214/aos/1176344136

CrossRef Full Text | Google Scholar

Schweizer, T. A., Kan, K., Hung, Y., Tam, F., Naglie, G., and Graham, S. J. (2013). Brain activity during driving with distraction: an immersive fMRI study. Front. Hum. Neurosci. 7:53. doi: 10.3389/fnhum.2013.00053

PubMed Abstract | CrossRef Full Text | Google Scholar

Seth, A. K. (2009). Causal connectivity of evolved neural networks during behavior. Network 16, 35–54. doi: 10.1080/09548980500238756

PubMed Abstract | CrossRef Full Text | Google Scholar

Seth, A. K. (2010). A MATLAB toolbox for Granger causal connectivity analysis. J. Neurosci. Methods 186, 262–273. doi: 10.1016/j.jneumeth.2009.11.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Shiogai, Y., Stefanovska, A., and McClintock, P. V. E. (2010). Nonlinear dynamics of cardiovascular ageing. Phys. Rep. 488, 51–110. doi: 10.1016/j.physrep.2009.12.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Sitaram, R., Caria, A., Veit, R., Gaber, T., Ruiz, S., and Birbaumer, N. (2014). Volitional control of the anterior insula in criminal psychopaths using real-time fMRI neurofeedback: a pilot study. Front. Behav. Neurosci. 8:344. doi: 10.3389/fnbeh.2014.00344

PubMed Abstract | CrossRef Full Text | Google Scholar

Tan, Q., Zhang, M., Wang, Y., Zhang, M., Wang, B., Xin, Q., et al. (2016). Age-related alterations in phase synchronization of oxyhemoglobin concentration changes in prefrontal tissues as measured by near-infrared spectroscopy signals. Microvasc. Res. 103, 19–25. doi: 10.1016/j.mvr.2015.10.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Tsunashima, H., and Yanagisawa, K. (2009). Measurement of brain function of car driver using functional near-infrared spectroscopy (fNIRS). Comput. Intell. Neurosci. 2009:164958. doi: 10.1155/2009/164958

PubMed Abstract | CrossRef Full Text | Google Scholar

Uchiyama, Y., Ebe, K., Kozato, A., Okada, T., and Sadato, N. (2003). The neural substrates of driving at a safe distance: a functional MRI study. Neurosci. Lett. 352, 199–202. doi: 10.1016/j.neulet.2003.08.072

PubMed Abstract | CrossRef Full Text | Google Scholar

Uchiyama, Y., Toyoda, H., Honda, M., Yoshida, H., Kochiyama, T., Ebe, K., et al. (2008). Functional segregation of the inferior frontal gyrus for syntactic processes: a functional magnetic-resonance imaging study. Neurosci. Res. 61, 309–318. doi: 10.1016/j.neures.2008.03.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, W., Wang, B., Bu, L., Xu, L., Li, Z., and Fan, Y. (2016). Vigilance task-related change in brain functional connectivity as revealed by wavelet phase coherence analysis of near-infrared spectroscopy signals. Front. Hum. Neurosci. 10:400. doi: 10.3389/fnhum.2016.00400

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, B., Zhang, M., Bu, L. G., Xu, L. W., Wang, W., and Li, Z. Y. (2016). Posture-related changes in brain functional connectivity as assessed by wavelet phase coherence of NIRS signals in elderly subjects. Behav. Brain Res. 312, 238–245. doi: 10.1016/j.bbr.2016.06.037

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, L., Zhang, J., Zhang, Y., Yan, R., Liu, H., and Qiu, M. (2016). Conditional granger causality analysis of effective connectivity during motor imagery and motor execution in stroke patients. Biomed. Res. Int. 2016:3870863. doi: 10.1155/2016/3870863

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, L., Wang, B., Xu, G., Wang, W., Liu, Z., and Li, Z. (2017). Functional connectivity analysis using fNIRS in healthy subjects during prolonged simulated driving. Neurosci. Lett. 640, 21–28. doi: 10.1016/j.neulet.2017.01.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Ye, J. C., Tak, S., Jang, K. E., Jung, J., and Jang, J. (2009). NIRS-SPM: statistical parametric mapping for near-infrared spectroscopy. Neuroimage 44, 428–447. doi: 10.1016/j.neuroimage.2008.08.036

PubMed Abstract | CrossRef Full Text | Google Scholar

Yoshino, K., Oka, N., Yamamoto, K., Takahashi, H., and Kato, T. (2013a). Correlation of prefrontal cortical activation with changing vehicle speeds in actual driving: a vector-based functional near-infrared spectroscopy study. Front. Hum. Neurosci. 7:895. doi: 10.3389/fnhum.2013.00895

PubMed Abstract | CrossRef Full Text | Google Scholar

Yoshino, K., Oka, N., Yamamoto, K., Takahashi, H., and Kato, T. (2013b). Functional brain imaging using near-infrared spectroscopy during actual driving on an expressway. Front. Hum. Neurosci. 7:882. doi: 10.3389/fnhum.2013.00882

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, L., Zhong, G., Wu, Y., Vangel, M. G., Jiang, B., and Kong, J. (2010). Using Granger-Geweke causality model to evaluate the effective connectivity of primary motor cortex (M1), supplementary motor area (SMA) and cerebellum. J. Biomed. Sci. Eng. 3, 848–860. doi: 10.4236/jbise.2010.39115

PubMed Abstract | CrossRef Full Text | Google Scholar