Yuxuan Chen
Julia Tang2
Han Yuan- 1School of Electrical and Computer Engineering, University of Oklahoma, Norman, OK, United States
- 2Stephenson School of Biomedical Engineering, University of Oklahoma, Norman, OK, United States
- 3Fran and Earl Ziegler College of Nursing, University of Oklahoma Health Sciences Center, Oklahoma City, OK, United States
- 4Institute for Biomedical Engineering, Science, and Technology, University of Oklahoma, Norman, OK, United States
Recently, functional near-infrared spectroscopy (fNIRS) has been utilized to image the hemodynamic activities and connectivity in the human brain. With the advantage of economic efficiency, portability, and fewer physical constraints, fNIRS enables studying of the human brain at versatile environment and various body positions, including at bed side and during exercise, which complements the use of functional magnetic resonance imaging (fMRI). However, like fMRI, fNIRS imaging can be influenced by the presence of a strong global component. Yet, the nature of the global signal in fNIRS has not been established. In this study, we investigated the relationship between fNIRS global signal and electroencephalogram (EEG) vigilance using simultaneous recordings in resting healthy subjects in high-density and whole-head montage. In Experiment 1, data were acquired at supine, sitting, and standing positions. Results found that the factor of body positions significantly affected the amplitude of the resting-state fNIRS global signal, prominently in the frequency range of 0.05–0.1 Hz but not in the very low frequency range of less than 0.05 Hz. As a control, the task-induced fNIRS or EEG responses to auditory stimuli did not differ across body positions. However, EEG vigilance plays a modulatory role in the fNIRS signals in the frequency range of less than 0.05 Hz: resting-state sessions of low EEG vigilance measures are associated with high amplitudes of fNIRS global signals. Moreover, in Experiment 2, we further examined the epoch-to-epoch fluctuations in concurrent fNIRS and EEG data acquired from a separate group of subjects and found a negative temporal correlation between EEG vigilance measures and fNIRS global signal amplitudes. Our study for the first time revealed that vigilance as a neurophysiological factor modulates the resting-state dynamics of fNIRS, which have important implications for understanding and processing the noises in fNIRS signals.
Introduction
Functional near-infrared spectroscopy (fNIRS) is a non-invasive functional neuroimaging technique that can monitor concentration changes in oxygenated and deoxygenated hemoglobin (HbO and HbR) in the cerebral cortex. fNIRS measurement is based on the absorption of light in near-infrared spectrum from 700 to 1000 nm by biological tissues. Different chromophores, such as hemoglobin, myoglobin, and cytochrome aa3, have different absorptivity (Sood et al., 2015). With the advantage of low-cost, portability, and ease to co-register with other neural recording modalities, such as an EEG and fNIRS has become an attractive means for imaging and monitoring hemodynamic signals in the human brain, which complements the use of fMRI in versatile environment. fNIRS has been widely applied in functional neuroimaging (Torricelli et al., 2014; Chen et al., 2020), cerebral monitoring in neonates (Sood et al., 2015; Hu et al., 2020) and brain-computer interface (Naseer and Hong, 2015; Ahn and Jun, 2017; Shin et al., 2017). Unlike fMRI constraining subjects to lying down on a scanner bed, fNIRS poses fewer physical constraints on the participants, thereby permitting them to be studied at flexible body positions during recordings.
Particularly, imaging of resting-state functional connectivity (RSFC) in the human brain has been a recent focus for neuroimaging studies, including fNIRS (Mohammadi-Nejad et al., 2018; Pinti et al., 2018). The activity of the resting brain exhibits spontaneous and large-amplitude fluctuations, which have been observed in a number of imaging modalities such as fMRI (Biswal et al., 1995), positron emission tomography (Raichle et al., 2001; Watabe and Hatazawa, 2019), and direct measures of neuronal activity with electro- or magneto-encephalography (EEG or MEG) (Goldman et al., 2002; Mantini et al., 2007; Brookes et al., 2011; Yuan et al., 2012, 2016). The measures of resting-state cerebral hemodynamics, mostly using fMRI based on the blood-oxygenation-level dependent (BOLD) contrast, show fluctuations predominantly at a low frequency band of <0.1 Hz (Cordes et al., 2001). The temporal synchrony across brain regions have been revealed (Beckmann et al., 2005; Damoiseaux et al., 2006), and demonstrated to be important biomarkers for the brain at diseased conditions (Zhang and Raichle, 2010). Prior studies of RSFC in both healthy and diseased conditions can be influenced by the presence of a strong global component, which is usually observed throughout sampled voxels or sensors, thereby dominating the RSFC (Greicius et al., 2003; Fox et al., 2005, 2009). However, the approach of removing global signal has recently been shown to induce systematic biases and the anti-correlation enhanced by global signal regression (GSR) becomes the main concern (Fox et al., 2009; Murphy et al., 2009). Furthermore, evidences show that a neural component (Scholvinck et al., 2010; Wong et al., 2013, 2016) and even diagnostic information (Hahamy et al., 2014; Murphy and Fox, 2017; Yang et al., 2017) exist in the global signal, which challenges the assumption of removing it in the first place.
Like fMRI signals, fNIRS also offers the potential to examine the human brain at resting state by measuring concentration changes of HbO and HbR in the vasculature of the cortical tissues below sensing channels (Obrig and Villringer, 2003; Scholkmann et al., 2014). fNIRS has been effectively employed to characterize the resting-state brain in adults (Obrig et al., 2000; White et al., 2009; Lu et al., 2010; Mesquita et al., 2010; Zhang et al., 2010; Sasai et al., 2011; de Souza Rodrigues et al., 2019), infants or children (Homae et al., 2010; White et al., 2012; Molavi et al., 2013; Watanabe et al., 2017; Bulgarelli et al., 2019, 2020; Wang et al., 2020), and to assess differences between experimental groups (Keehn et al., 2013; Ieong et al., 2019; Arun et al., 2020). The most common RSFC analysis of fNIRS data involves evaluating the temporal relationship between time series of the preprocessed data from recording units, for example, through the Pearson’s correlation. A global component has been observed in fNIRS measurements and commonly removed for the purpose of attenuating systematic noises at the resting state (White et al., 2009; Gregg et al., 2010; Mesquita et al., 2010; Eggebrecht et al., 2014; Tachtsidis and Scholkmann, 2016; Duan et al., 2018; Sherafati et al., 2020; Wyser et al., 2020). Whereas removing superficial contributions from short-distanced channels to fNIRS is increasingly employed to attenuate the systematic noises (Saager and Berger, 2005; Gagnon et al., 2011), data from both long-distanced and short-distanced channels commonly suggest a global component exist in fNIRS measurements and distribute across wide regions (Zhang et al., 2005, 2007, 2009; Kohno et al., 2007; Gregg et al., 2010; Tong and Frederick, 2010; Novi et al., 2016; Sato et al., 2016). However, the physiological nature of the fNIRS global signal has not been fully established, since the neurophysiological components in the resting-state global fNIRS signal have not been systematically investigated. Therefore, whether or not to remove the global signal in fNIRS-based RSFC analysis remains not clear.
In this study, we aimed to investigate the physiological underpinning of resting-state fNIRS global signal by concurrently acquiring fNIRS and EEG in whole-head, high-density montage. Previous studies using simultaneous EEG and fMRI recordings have revealed a negative association between the amplitude of resting state fMRI global signal and EEG vigilance level (Wong et al., 2013; Chang et al., 2016; Falahpour et al., 2018). These studies have shown that subjects at lower vigilance states are characteristic of higher global signal amplitudes, indicating that neurophysiological covariates exist in the global signal. Moreover, the temporal fluctuations of vigilance levels has been linked to the spontaneous activities in regions constituting the default mode network (DMN) (Olbrich et al., 2009) and also linked to the fluctuations of fMRI global signal (Chang et al., 2016; Falahpour et al., 2018), suggesting that regressing out the resting-state global signal could potentially impact the dynamic connectivity in resting state networks. Based on the prior studies using BOLD fMRI, in the current study we hypothesize that the fNIRS global signal has a neurological component and is related to the EEG vigilance. Two experiments were included: the 1st is a 10-min resting study from which we calculated the stationary metrics; the 2nd is a 45-min resting study from which we examined the epoch-to-epoch dynamics during the wakeful epochs. Furthermore, considering that fNIRS is a promising technology for imaging the human brain at versatile body positions, Experiment 1 also examined whether and how different body positions affect the fNIRS global signal at resting state conditions and, as a control, the impact of body positions on evoked activities to auditory stimuli was studied. To our knowledge, our study is the first of its kind to examine the relationship between fNIRS global signal and EEG vigilance with an advanced simultaneous EEG and fNIRS system in high-density and whole-head montage.
Materials and Methods
Protocol
All study procedures were completed according to the Declaration of Helsinki guidelines and approved by the Institutional Review Board at the University of Oklahoma Health Sciences Center.
Experiment 1: 10-min Resting at Different Body Positions
Twenty-four healthy subjects were recruited after giving informed consent. All subjects were right-handed. Two subjects were excluded due to excessive movements during the recording. Another three subjects were excluded due to high electrode impedance in recording sessions. Thus, nineteen subjects’ data were included in the analysis (11 males and 8 females, aged 19 to 55 years old, average age ± STD = 30.8 ± 12.2 years). Each subject participated in two separate sessions, eyes-open (EO) and eyes-closed (EC), the order of which was randomized. For each subject, two sessions occurred on different days that were within a 4-week period (mean interval ± STD = 18.2 ± 29.5 days). Each session contained three recording blocks at different body positions: standing, sitting, and supine. The order of these blocks was randomized among all subjects but was kept the same for each subject at their consecutive visits. Each block contained a 10-min resting-state part and an auditory task part that lasted 6 min and 30 s. In the resting-state part, subjects were instructed to keep as still as possible and not fall asleep. Specifically, in EO resting condition, subjects were instructed to focus on a black cross on a white background. In the auditory task part, subjects were instructed to keep still and listen to the auditory stimuli from a pair of earbuds. In terms of the presentation, periods of 30-s stimuli on and 30-s stimuli off were interleaved. Within a period of stimuli on, subjects listened to a sequence of 15 brief one-kilohertz tones. One tone lased for 100 ms duration and sampled at 44.1 kHz. Tones within a stimuli-on period were separated by an inter-stimulus interval of 2 s. Six datasets (two eye conditions by three body positions) were obtained for each subject, yielding a total of 114 datasets in the current study, which included concurrent EEG and fNIRS data of both resting state and task conditions.
Experiment 2: 45-min Resting at Supine Position
Because Experiment 1 used a stationary design of 10-min resting study to investigate the fNIRS resting state signals, we further included Experiment 2 in the design of a 45-min resting study to examine the temporal dynamics in fNIRS global signal. Specifically, subjects were instructed to rest still and allowed to fall asleep during a 45-min recording, while subjects laid supine in an adjustable recliner. A total of 20 healthy subjects (sex: 12 females and 8 males; aged 28 to 63 years old; average age ± STD: 42.8 ± 11.7 years) were studied in Experiment 2 and no subjects overlapped in Experiments 1 and 2. The recording began and ended with bio-calibration, which were used to identify artifacts in the EEG recordings. The bio-calibration procedure was performed in a standard order of instructing subjects to (1) open and close their eyes, (2) blink, (3) perform lateral eye movements, (4) take deep breaths, (5) clench their teeth, and to (6) speak.
fNIRS Data Acquisition
An identical configuration of acquisition was used in Experiments 1 and 2. The fNIRS measurements were acquired with a NIRScout system (NIRx Medical Technologies, LLC, New York, United States). Thirty-two source probes and 32 detector probes were plugged into holders and arranged into a cap based on the international 10–5 system (Jasper, 1958). A total of 105 channels (i.e., 105 pairs of sources and detectors) were defined, covering the areas from the forehead to the occipital lobe. The inter-optode distance varied between 25, 27, and 30 mm, corresponding to three different sizes of caps (54, 58, and 60 cm). The intersection between the left and right tragus and the Nasion and Inion was the center of the cap, which was denoted by the Cz position. A dark black over-cap covered the cap to block external light luminance. The absorption of near-infrared light of two different wavelengths (760 and 850 nm) was measured with a sampling rate of 1.95 Hz.
EEG Data Acquisition
A 64-channel, fNIRS-compatible EEG system (BrainProducts, München, Germany) was utilized to record the EEG data. In order to couple the EEG signal with the fNIRS hemodynamic signal, the montage of the EEG electrodes was designed to match the fNIRS montage. Every EEG channel was crossed by an adjacent pair of light source and detector. Sixty-four EEG electrodes were also mounted onto corresponding holders. The electrode at FCz position was selected as the reference point. Two 32-channel amplifiers, which were powered by a rechargeable battery, were included in our EEG system. Electrically conductive gel was added to decrease the impedance between scalp and electrodes. The impedances of EEG electrodes were kept under 20 kΩ throughout the recordings. All the EEG datasets were digitized with a wide band of 0.1–250 Hz at a 500 Hz sampling rate.
fNIRS and EEG Data Pre-processing
Figure 1 shows the analysis flowchart of EEG and fNIRS data. EEGLAB (Delorme and Makeig, 2004) was used for pre-processing of EEG data. After loading the raw datasets, the data was re-referenced to the common average reference. A basic FIR bandpass filter from 0.1 Hz to 70 Hz was used to filter the data in addition to a notch filter of 60 Hz. Additional ocular and muscular artifacts were removed by the independent component analysis implemented in EEGLAB. The ocular components, muscle movement components, and other artifacts were manually inspected and removed (Chaumon et al., 2015). Preprocessed EEG data were down-sampled to 250 Hz.
Figure 1. Flowchart of the data processing flowchart. (Left) fNIRS signal processing procedures, and (Right) EEG signal processing procedures. Dashed line circles the pre-processing procedures.
fNIRS data was pre-processed in HOMER2 (Huppert et al., 2009). Channels consisted of a source electrode and adjacent detector electrodes. Montages were created according to the setup of sources and detectors. Preprocessing of fNIRS data included converting raw light intensity to optical density, principle component analysis (PCA) removal, and motion artifact detection and correction. The PCA algorithm we performed here is to filter out the first principal component (Novi et al., 2016). Discontinuities and spikes existing in recordings were replaced by an average of its adjacent data segment. All channels were bandpass filtered from 0.01 to 0.2 Hz. The resulted time series were subject to hemodynamic computation via the modified Beer-Lambert law (Kocsis et al., 2006), yielding relative changes in concentrations of Oxy-Hemoglobin (HbO) and Deoxy-Hemoglobin (HbR) (Gratzer and Allison, 1960; Putze et al., 2014).
Calculation of Resting-State fNIRS Global Signal
After pre-processing, the fNIRS data became a measure of the relative concentration changes of HbO and HbR in units of μM. Then the preprocessed fNIRS were separated into two frequency bands: the lower range of <0.05 Hz and the upper range of >0.05 Hz that contains the Mayer wave (Julien, 2006), guided by inspection of power spectrum and analysis of variance (ANOVA) tests. To calculate the global signal, the time series of relative changes in HbO or HbR were averaged across all channels covering the whole brain. Then, the amplitude of the global signal was defined by the standard deviation of the global signal time series.
Quantification of EEG Vigilance Level
For EEG data, after removing artifacts, a spectrum was calculated by using Welch’s power spectral density estimate with segments of 10 s and 50% overlap for each EEG channel. Then the spectrum was divided by its overall root mean square (RMS) amplitude of all frequency bins, resulting in the relative amplitude spectrum. The relative amplitude spectrums were then averaged across all channels. Three frequency bands (delta: 1 – 4 Hz, theta: 4 – 7 Hz, alpha: 7 – 13 Hz) were delineated, and the RMS amplitudes were calculated separately for each band. A measure of EEG vigilance was defined as the RMS amplitude in the alpha band divided by the sum of RMS amplitudes in the delta and theta bands (Horovitz et al., 2008; Wong et al., 2013), which is equivalent to the alpha slow-wave index (ASI) (Jobert et al., 1994).
Auditory fNIRS and EEG Data
Whereas a key investigation of the Experiment 1 was to examine the impact of body position on resting-state neural recordings, we also included investigation of the body position on task-induced responses in EEG and fNIRS, in order to control systematic and environmental nuisances. A mixed block and event-related design were used for the concurrent EEG and fNIRS recordings. One session contained six task blocks and each block contained 15 auditory stimuli. The auditory stimuli were controlled by E-Prime software (Psychology Software Tools, Pennsylvania, United States). Stimuli was sent to earbuds by the stimulation computer. The trigger pulse corresponding to the sound then marked EEG and fNIRS synchronously via a parallel control box. There was a total of six task conditions (standing/sitting/supine body positions X EO/EC conditions).
For fNIRS analysis, block average was obtained after preprocessing in the same way described in resting data. The first marker of a task block was kept as the start of each task block. Based on all available auditory markers (as time 0 s), the time series were demeaned with reference to the time window from −5 to 0 s and averaged, resulting in the auditory response waveform. Segments containing detected motion artifacts were excluded from the average. To visualize the time courses of hemodynamic responses, the fNIRS auditory response was plotted from −10 to 50 s.
For EEG analysis, auditory evoked potentials (AEP) were obtained. Specifically, the recordings were band-pass filtered between 0.1 and 30 Hz, down sampled to 250 Hz, referenced to common average reference, and segmented into epochs from −100 to 500 ms. For every single segment, the t = 0 s denotes the onset of auditory stimuli. The mean of the baseline (averaged from −100 to 0 ms) was subtracted from the time series. Ocular and muscular artifacts were removed by the independent component analysis implemented in EEGLAB. Visual inspection further excluded the trials containing motion artifacts. Remaining trials of EEG epochs were averaged for each auditory task condition, resulting the AEP waveforms.
Sleep Stage Scoring
In Experiment 2, the 45-min recording was reviewed and manually scored into sleep stages by a certified expert (BWC), using standard scoring criteria by the American Academy of Sleep Medicine (Berry et al., 2017). Briefly, EEG data were first segmented into epochs of 30-s length. Based on the frequency and amplitude of the signal, each segment was assigned as awake, non-rapid-eye-movement sleep, Slow Wave Sleep, or rapid-eye-movement sleep. Only epochs of awake stages before first sleep onsite were included in the analysis.
Statistical Analysis
In Experiment 1, in order to explore the effect of body position on neural recordings (i.e., standing, sitting and supine), ANOVA was applied on the EEG or fNIRS quantities, separately for the EO and EC conditions. We performed the statistical test on each frequency bin along a continuous spectrum (Figures 3, 4). Then, based on the delineation of frequencies, we segregated the quantities of the fNIRS global signal in two bands: f < 0.05 and f > 0.05 Hz.
Next, two-way repeated measures ANOVA (standing/sitting/supine body positions X EO/EC) was applied to assess if any main effect of body position or eye condition, or interaction between the body position and the eye condition, separately in the frequency range of <0.05 and >0.05 Hz and separately for HbO and HbR. Likewise, two-way repeated measures ANOVA (supine/sitting/standing body positions X EO/EC) was tested on the EEG vigilance scores. Furthermore, post hoc analysis assessed the difference between conditions using a paired, two-sided t-test. Bonferroni correction was used to correct the multiple comparison problem.
After delineating the position and eye effects, we assessed whether fNIRS global signal amplitudes were associated with the EEG vigilances. Particularly, motivated by a negative relationship between fMRI global signal and EEG vigilance reported in the literature, we tested whether higher fNIRS global signals are associated with lower vigilance. The analysis has excluded the frequency band of greater than 0.05 Hz that contains the Mayer wave. Also, the analysis only considered the eyes-open condition to exclude the body position factor on EEG or fNIRS. Then, per each individual, the EEG vigilance at three body positions was sorted into highest, medium and lowest levels; and the fNIRS global signal associated with the highest vigilance levels were compared to fNIRS global signals at the medium or lowest level using a paired, two-sample and one-sided t-test. Furthermore, for the purpose of determining whether vigilance variations underlie the fluctuations in fNIRS global signal, the co-variation was assessed in one eye condition at one body position across all subjects, because practically resting state experiments are conducted in a single experiment condition rather than combined. Partial correlations between global signal amplitude and vigilance measures were calculated, by controlling age and gender as confounding factors.
While Experiment 1 focused on the stationary properties of fNIRS and EEG, we further evaluated whether the epoch-by-epoch fluctuations in fNIRS global signal and EEG vigilance are associated in Experiment 2. For each subject, we calculated the fNIRS global signal amplitudes and EEG vigilance measures in 30-s epochs. Then the temporal correlation between the fNIRS global signal and EEG vigilance was calculated using Pearson’s correlation coefficient across all epochs per each subject. To assess the temporal correlations at a group level, the correlation coefficients were converted to z scores using the Fisher’s transform. Afterward, a one-sample t-test against 0 was performed on all individuals’ z scores to evaluate the significance of temporal correlation at a group level.
Results
The aim of the study was to investigate the neurological basis of the fNIRS resting-state global signal, if any, and the impact of body positions on the resting-state signals. The results are organized as such: In Experiment 1, the frequency-dependent impact of body positions on fNIRS and EEG signals was explored, then the factors of body positions and eye conditions were assessed in fNIRS global signal in two delineated frequency bands, and finally, the co-variation in the amplitude of fNIRS global signal and EEG vigilance was analyzed. As control results, the fNIRS and EEG task responses to auditory stimuli were included. In Experiment 2, the epoch-to-epoch fluctuations of fNIRS global signal amplitudes and EEG vigilance measures across 30-s epochs were examined and their temporal correlation was reported.
Experiment 1: 10-Min Resting at Different Body Positions
Firstly, spontaneous fluctuations were observed in the fNIRS global signal when subjects rested with their eyes open and closed, without any external stimuli. Representative single-session traces of fNIRS global signals are shown in Figure 2, at an EC resting condition. Notably, the global signal at all three positions exhibit fluctuations with a peak frequency of ∼0.02 Hz. Meanwhile, the data acquired from these body positions exhibited different patterns of fluctuations in the time domain, i.e., slower fluctuations are observed in the supine position and faster fluctuations in the sitting and standing positions. In terms of the amplitude, we noted that the power spectrum at the supine position showed a lowest amplitude in the frequency range of 0.05 – 0.1 Hz than those at sitting and standing positions, in the representative subject.
Figure 2. Representative traces of fNIRS global signal derived from HbO signal when the subject closed the eyes and rested in (A) supine position (in gray color), (B) sitting position (in orange color), and (C) standing position (in blue color), exhibited in the time domain (left panel) and frequency domain (right panel).
Figure 3. The grand average of the power spectrum of fNIRS global signals at resting state. (A) and (B) show HbO signal at EO and EC conditions. (C) and (D) show HbR signal at EO and EC conditions, respectively. The gray, orange, and blue curves represent supine, sitting, and standing position, respectively, in all panels.
Figure 4. The grand average of the power spectrum of EEG resting-state signals at (A) EO condition and (B) EC condition. The gray, orange, and blue curves represent supine, sitting, and standing position, respectively, in both panels.
Furthermore, the position-dependent profile of the resting-state fNIRS global signal is also prominent at the group level. Figure 3 shows the grand average of the power spectrum of fNIRS global signal at various resting state conditions. Notably, in the frequency range from 0.05 Hz to 0.1 Hz, the spectrums at three different body positions show largely different amplitudes. The spectrum at the standing position appears to be of highest amplitudes in 0.05 – 0.1 Hz, at both EO and EC conditions (Figures 3A,B) in blue curves, whereas spectrum at supine are of lowest amplitudes (Figures 3A,B) in gray curves. In order to delineate the frequency-dependent effect of body position, we performed one-way ANOVA on the amplitude of fNIRS global signal separately in each frequency bin. At the EO condition, between 0.05 and 0.09 Hz, the effect of body position was significant on HbO (p < 0.05, uncorrected). Similarly, at the EC condition, the effect of body position was significant in the range from 0.07 to 0.09 Hz on HbO (p < 0.05, uncorrected). Since the fNIRS signal in the frequency range of 0.05 – 0.1 Hz has been related to a physiological noise known as the Mayer wave (Julien, 2006), our later analysis of the fNIRS global signal then focused on two distinct frequency bands, i.e., f < 0.05 and f > 0.05 Hz, in order to distinguish a position-dependent impact that may be attributed to physiological noises. In HbR (Figures 3C,D), we used the same frequency bands as with HbO. Noteworthy, none of HbO or HbR data showed any significant effect of body positions in a frequency bin less than 0.05 Hz (p > 0.05, uncorrected).
Likewise, in the resting-state EEG, our analysis explored whether a position-dependent profile exists on the spectrum. Figures 4A,B show the grand average of the power spectrum at EC and EO conditions, respectively. ANOVA revealed that the body position was not significant in any of the frequency bins at either EO or EC conditions (p > 0.05, uncorrected). Notably, although the grand average at the EC conditions appears with different amplitudes for three different conditions, it did not reach a significance level (p = 0.067 at f = 10.8 Hz, uncorrected).
Next, we aggregated the fNIRS and EEG quantities as the amplitude of global signal and vigilance scores, respectively. In particular, we averaged the amplitude of the fNIRS global signal (as root-mean-square) in the range of f < 0.05 Hz, which excludes the Mayer wave, and then separately in the range of f > 0.05 Hz. Meanwhile, EEG vigilance scores were calculated based on the power spectrum of resting state EEG as the ratio of alpha-band RMS divided by the sum of delta- and theta-band RMS. Two-way Repeated Measures ANOVA (body positions × eye conditions) revealed that the effect of body positions was not significant on fNIRS global signal amplitude in the very low frequency range of f < 0.05 Hz (q = 0.10). Meanwhile, the effect of body position was significant on the fNIRS global signal in the range of f > 0.05 Hz (q < 0.001). Noteworthy, the interaction of body positions and eye condition was not significant in fNIRS global signal amplitude in either frequency range.
Post hoc comparison on HbO in the range of f > 0.05 Hz was then conducted to assess the difference between pairs of body positions (i.e., standing vs. supine, sitting vs. supine, and standing vs. sitting) (Figure 5B). Analysis showed that the amplitude of fNIRS global signal for the supine position was significantly lower than the sitting position (q < 0.01) and standing position (q < 0.001), after multiple comparison correction. But amplitude of fNIRS global signal for sitting position did not differ from the standing position. Noteworthy, neither the eye factor or the eye-position interaction was significant in fNIRS HbO or HbR data.
Figure 5. Effects of body positions (supine, sitting, and standing) and eye conditions (EC and EO) on (A–D) fNIRS global signal amplitude, and (E) EEG vigilance measurement. The effect of body position was significant in HbO global signal amplitude, in the range of f > 0.05 Hz that contains Mayer wave (B). EEG vigilance measurement showed significance in the effect of body position, the effect of eye conditions and the position-eye interaction (E). Stars indicate post hoc significance after multiple comparison correction (∗ indicates q < 0.05, ∗∗ indicates q < 0.01, ∗∗∗ indicates q < 0.001). Error bars indicate standard error.
In terms of EEG vigilance scores (Figure 5E), the two-way repeated measure ANOVA found that the effect of body position, the effect of eye condition and the eye-position interaction was all significant (q < 0.001). Post hoc comparisons on the EEG vigilance scores were then conducted to assess the differences. Informed by the significant interaction factor, we performed separate ANOVA analysis on the effect of body positions at separate eye condition and also performed separate t-test on pairs of body positions and eye conditions. Only under EC condition, the supine position had significant smaller EEG vigilance than sitting (q < 0.05) and standing position (q < 0.001). Furthermore, regarding the eye factor (EO vs. EC), the EEG vigilance showed significance at both sitting (q < 0.01) and standing positions (q < 0.001), but not in the supine position. However, under EO condition, there was no significant effect of body positions.
As a next step, we examined the relationship between the fNIRS global signal and the EEG at resting state. Firstly, we tested whether higher fNIRS global signals are associated with lower vigilance, which was motivated by a negative relationship between fMRI global signal and EEG vigilance reported in the literature. Particularly, we compared the amplitude of fNIRS global signal in the frequency range of f < 0.05 Hz against the EEG vigilance scores, only at EO state when either fNIRS or EEG quantities were not impacted by the factor of body positions. Results in Figure 6 shows a reversed association was identified between EEG vigilance and fNIRS global signals. After sorting the vigilance measures within each individual, the resting sessions of lowest vigilance were associated with significantly higher HbO global signals [t(18) = 2.02, p < 0.05] and also higher HbR global signals [t(18) = 2.98, p < 0.01] than those of highest vigilance. Importantly, note that neither the vigilance nor the global signal differed between body positions at the eyes-open condition; nonetheless, a reversed association was still found between EEG and fNIRS.
Figure 6. A reversed association between EEG vigilance and fNIRS global signals. Resting state sessions of lowest vigilance measures were associated with significantly higher amplitudes of fNIRS global signals in the (A) HbO (p < 0.05) and (B) HbR (p < 0.01). ∗p < 0.05, ∗∗q < 0.01, and error bars indicate standard error.
In addition, we examined the co-variation between fNIRS global signal and EEG per each body position across individuals. Results showed a consistent negative trend such that higher global signals are associated with lower vigilance states. In particular, both HbO and HbR at the standing position significantly co-varied with the EEG vigilance after controlling age and gender as confounding factor (HbO: r = −0.51, p < 0.05; HbR: r = −0.57, p < 0.05), as shown in Figure 7. However, at other positions, the covariation did not reach significance after multiple comparison correction, although a negative trend in the association was consistently noted. HbR at the supine position showed a significance-approaching covariation with EEG vigilance (r = −0.36, p = 0.1) and HbO at the sitting position also approached significance (r = –0.34, p = 0.1). HbO at supine position (r = −0.06) and HbR at sitting position (r = −0.20) did not reach a significant covariation with EEG vigilance.
Figure 7. Amplitude of fNIRS global signal in the range of f < 0.05 Hz co-varied with EEG-based measure of vigilance, when subjects rested in the standing position with their eyes open. fNIRS global signal amplitudes derived from HbO signal (A) and from HbR signal (B) both significantly co-varied with EEG vigilance after controlling age and gender as confounding factors (HbO: r = –0.51, p < 0.05, HbR: r = –0.57, p < 0.05).
In order to control systematic and environmental nuisances, task-induced responses in EEG and fNIRS were investigated. Figure 8 shows the grand average of EEG AEP curves at three different body positions, from the FCz electrode. In either EO and EC conditions, the AEP curves at three body positions followed a very similar profile: negative activities at the ∼100 ms (N1) and positive activities at the ∼200 ms (P2). The factor of body positions has no significant effect on N1 or P2 (q > 0.1). Furthermore, the eye-position interaction is not significant, either (q > 0.1).
Figure 8. Grand average of EEG auditory evoked potentials at (A) eyes-open (EO) condition, and at (B) eye-closed (EC) condition. The gray, orange, and blue curves represent supine, sitting, and standing positions in both panels.
Meanwhile, task-related fNIRS responses were averaged across the blocks after subtracting the activities between −5 s and 0 s, with time 0-s as the beginning of the block. Figure 9 shows the grand average of fNIRS response to auditory stimuli. Representative time courses from the channels located over the left and right auditory cortex regions are shown in Figure 9B. When the auditory stimuli were on (shaded gray area in Figure 9B), relative changes of the HbO increased while the relative changes of the HbR decreased. In order to visualize the topography of auditory task responses, the block-average time series of fNIRS was selected from 10 to 30 s, averaged, and shown in Figure 9A. Positive activations of HbO are shown in left and right auditory cortex in the topography of HbO (Figure 9B). Results revealed no significant effect of body positions on the fNIRS auditory response.
Figure 9. Grand average of fNIRS auditory response derived from HbO signals. (A) shows the topography of averaged HbO responses between 10 to 30 s after the stimulus onsite). (B) plots the time courses of fNIRS auditory response obtained from representative channels over left and right auditory cortex (orange: left auditory cortex HbO, blue: right auditory cortex HbO, gray: left auditory cortex HbR, yellow: right auditory cortex HbR). Shaded gray area indicates the time window of auditory stimulus. Time zero is the onsite of stimulus.
Experiment 2: 45-Min Resting at Supine Position
In Experiment 2, subjects were instructed to rest for a total duration of 45 min while allowed to fall asleep. Sleep scoring found that all 19 were able to fall asleep. Only data at wakeful resting before any sleep were used in the current analysis (Mean ± SD = 711 ± 602 s, ranging from 150 to 2430 s). After quantifying vigilance and fNIRS global signal amplitude in 30-s epochs, temporal fluctuations were observed. The time courses of HbO, HbR, and EEG vigilance measurement in a representative subject are displayed in Figure 10A. Every dot represents the global signal amplitude calculated from a windowed 30-s fNIRS signal and vigilance measurement calculated from 30-s EEG in the same aligned time window. As the subject gradually fell into the sleep, the vigilance exhibited a slowly decreasing trend, in reversed synchrony with surges of increases in fNIRS global signal and temporal. In the same time course, scattered moments of rebounds in vigilance are also aligned with drops of global signals, especially toward the later duration before the subject fell into sleep. In the same subject (shown in Figures 10B,C), a negative temporal correlation was observed between fNIRS global signal amplitude and EEG vigilance measurement (HbO vs. vigilance: r = −0.46, p = 0.004; HbR vs. vigilance: r = −0.61, p < 0.001).
Figure 10. Temporal correlation between epoch-to-epoch fluctuations of vigilance and fNIRS global signal in a representative subject at wakeful rest. (A) plots the time courses of EEG vigilance measure (black curve) and global signal amplitude of HbO (orange curve) and HbR (blue curve) calculated in 30-s epochs over 18 min. (B) shows the negative correlation between vigilance and HbO global signal amplitude is significant (r = –0.46, p = 0.004). (C) shows the negative correlation between vigilance and HbR global signal amplitude is significant (r = –0.62, p < 0.001).
At the group level, all subjects exhibited temporal variations in EEG vigilance. Specifically, the standard derivation of EEG vigilance over the period of wakeful rest ranged from 0.07 to 0.46 across all subjects (Mean ± SD = 0.25 ± 0.13). In terms of association, analysis showed that there existed significant correlation between HbO global signal amplitude and vigilance [t(19) = −2.57, p = 0.02), as well as significant for HbR and vigilance [t(19) = −2.80, p = 0.01]. Despite that subjects had different durations of wakefulness, the window length of wakefulness was not associated with the temporal correlation between fNIRS global and vigilance (p > 0.1 for HbO and HbR). When restricting the wakeful epochs to be within the 10 min before falling into sleep, the temporal correlation was still significant between HbO and vigilance [t(19) = −2.32, p = 0.03], and also significant between HbR and vigilance [t(19) = −2.83, p = 0.01]. Furthermore, considering that individuals exhibited different levels of vigilance fluctuations (i.e., standard derivation of EEG vigilance ranged from 0.07 to 0.46), we examined whether the standard derivation of EEG vigilance modulated the association between EEG vigilance measures and fNIRS global signal; yet the analysis showed that the scale of vigilance fluctuation levels was not relevant (p > 0.1).
Discussion
Our study has investigated the neurophysiological nature of the global signal of fNIRS measured at resting state. The results for the first time have demonstrated that the amplitude of the fNIRS global signal, particularly in the frequency range of 0.01 < f < 0.05 Hz, is reversely associated with EEG vigilance measures. The discovery of a neurological origin for fNIRS global signal has important implications for the processing of fNIRS signal acquired at resting state.
One of the most fundamental and critical issues in analyzing neuroimaging data is how to handle the global signal, which is defined as the time series of intensity averaged across imaging units in PET (Fox et al., 1988; Friston et al., 1990), fMRI (Desjardins et al., 2001; Macey et al., 2004), and more recently in fNIRS (Franceschini et al., 2006; Zeff et al., 2007; White et al., 2009; Mesquita et al., 2010). A strong presence of global signal in fMRI may lead to a massive and diffused activation pattern in task-based studies if the time series of the global signal is of similar profile with the task modulation (Kay et al., 2013; Power et al., 2015). Likewise, fNIRS studies of various tasks commonly removed a global component derived from the measurements to reveal focal activations, via linear regression or spatial filtering based on PCA/ICA decomposition (Zhang et al., 2005, 2007, 2016; Franceschini et al., 2006; Kohno et al., 2007; Zeff et al., 2007; Eggebrecht et al., 2012, 2014; Sato et al., 2016). Nonetheless, the impact of global signal is more problematic in task-free, resting state studies, as the global signal may lead to a perfusive connectivity pattern that is attributed to the global signal, no matter whichever seed region of interest is selected. Because region-specific connectivity is more desirable and because non-neuronal sources can dominantly contribute to the global signal (Glover et al., 2000; Wise et al., 2004; Birn et al., 2006; Yuan et al., 2013), the analysis of resting state fMRI data has commonly included steps to attenuate the impact of a global signal. For example, GSR removes an averaged signal of all recording units from the time series through linear regression. This procedure was originally developed for and applied to task-based fMRI data (Zarahn et al., 1997; Aguirre et al., 1998; Macey et al., 2004). Later, most resting-state fMRI studies have adopted GSR as a pre-processing approach: the global signal component is regressed out of preprocessed BOLD signals prior to computation of connectivity measures and therefore regionally focused connectivity patterns are reported (Fox et al., 2009). Similarly, in recent fNIRS studies of resting state brain, a global component has been recognized in the measurements from regularly distanced optodes (White et al., 2009; Mesquita et al., 2010; Tong and Frederick, 2010; Eggebrecht et al., 2014; Tachtsidis and Scholkmann, 2016; Duan et al., 2018; Wyser et al., 2020) and from short-distanced optodes (White et al., 2009; Gregg et al., 2010; Mesquita et al., 2010; Eggebrecht et al., 2014; Tachtsidis and Scholkmann, 2016; Duan et al., 2018; Sherafati et al., 2020; Wyser et al., 2020). To date, there is no well-established pre-processing routine in resting state fNIRS studies although multiple efforts are being made (Huppert et al., 2009; Ye et al., 2009; Xu et al., 2014; Santosa et al., 2018). Approaches such as GSR and spatial filtering via PCA and ICA decomposition that were used in task-based fNIRS studies are also commonly adapted in resting state fNIRS studies to remove the global component, yielding regionally focused connectivity pattern (Mesquita et al., 2010; Zhang et al., 2010, 2011; Sakakibara et al., 2016).
However, the removal of global signal in neuroimaging data has encountered controversial critiques, particularly in the studies of resting state functional connectivity. Because a global neurophysiological component may be present in direct neural recordings (Scholvinck et al., 2010; Wong et al., 2013, 2016), removing the global signal is shown to cause loss of such neural components, thereby confounding the resulted pattern of resting state functional connectivity. For example, Chen et al. (2012) found that the global signal is highly correlated with DMN component. Further evidences indicated that the global signal resembles the resting-state fMRI time courses of the largest cluster when the level of global noise is low (Chen et al., 2012). Under such circumstances, GSR could mathematically mandate the presence of anti-correlation network in fMRI studies (Murphy et al., 2009). Other studies have further linked the fluctuations of global signals to the varying levels of vigilance or arousal (Chang et al., 2016; Falahpour et al., 2018), which suggests that removing the global signal in those situations could remove an underlying behavioral factor. Therefore, the GSR should be very carefully applied when studying resting-state MRI (Murphy et al., 2009; Saad et al., 2012; Murphy and Fox, 2017). Until now, the nature of the fNIRS global signal has not been fully established since the neurophysiological components in the resting-state global fNIRS signal have not been systematically investigated. Our current study is the first of its kind to investigate the neuronal and non-neuronal sources in the fNIRS global signal by using concurrent fNIRS and EEG in whole-brain and high-density setup. Because both fNIRS and BOLD fMRI measure the cerebral hemodynamics, they carry similar substrates for neuronal activities while they also share common caveats due to non-neuronal sources, including respiration, cardiac pulsations, motion, etc. Like in the case of fMRI, removal of fNIRS global signal may lead to spurious results in the functional connectivity pattern, depending on whether or not there exists any neural component in the global signal of fNIRS and the amplitude level of global signal.
In this study, we have shown that fNIRS global signals acquired from the resting human brain are periodical oscillations. As shown in Figures 2, 3 at respective individual and group level, the resting-state fNIRS global signal resides in three ranges: dominantly less than 0.05 Hz with a peak component at ∼0.02 Hz, a second peak between 0.05 and 0.1 Hz (also known as the Mayer wave) and greater than 0.1 Hz. Furthermore, our study extended investigations of the fNIRS global signal at standing, sitting and supine positions. Indeed, periodic fluctuations were observed in the global signal at all body positions. The presence of a fluctuating fNIRS global signal with dominate activities of <0.1 Hz suggests that the RSFC pattern may be affected by the global signal. Comparing with intracranial neural recordings (Leopold et al., 2003; He et al., 2008; Shmuel and Leopold, 2008), fNIRS global signal and spontaneous neural activities overlap their peaks in the range of <0.1 Hz. Meanwhile, in comparison with fMRI, the fNIRS global signal shows a very similar spectral profile with those from BOLD fMRI. Especially, the spectrum of fNIRS at the supine position (Figures 3A,B) for both EO and EC conditions are almost identical to those reported in fMRI (e.g., Figure 1 in Biswal et al., 1995). Since in our study the whole-head fNIRS montage were sampled at 1.95 Hz, which is a higher frequency than BOLD fMRI (usually 0.5 Hz), the spectrum of fNIRS global signal revealed a more accurate spectrum.
Importantly, for the first time our study reported a negative association between the amplitude of fNIRS global signal in the range of <0.05 Hz and the EEG vigilance based on the simultaneous recording in the Experiment 1 (Figures 6,7). Within individuals, the resting state sessions with lowest EEG vigilance measures were associated with significantly higher fNIRS global (HbO and HbR in Figure 6), which was observed at eyes open condition and neither EEG nor fNIRS was affected by body positions. Furthermore, in a single body position at eyes-open condition, a negative covariation between fNIRS global signal amplitude and EEG vigilance was also confirmed across individuals (HbO and HbR in Figure 7). The selection of frequencies f < 0.05 Hz for fNIRS is critical: it is within the range of resting state fMRI data but distinctly narrower to exclude the Mayer wave. Previous fMRI study has demonstrated that the functional connectivity in auditory, visual and sensorimotor cortices is characterized 90% by the low-frequency band from 0 to 0.1 Hz (Cordes et al., 2001). Meanwhile, the fractional amplitude of low-frequency fluctuation (fALFF) is defined as the ratio of power spectrum of 0.01 – 0.08 Hz to that of the whole frequency band (Zou et al., 2008). Noteworthy, one of the most studied networks – DMN – has significantly higher fALFF than other brain regions, which indicates DMN has higher intensity of regional spontaneous brain activity in the range of 0.01 – 0.08 Hz (Zou et al., 2008). More importantly, our fNIRS signal was further narrowed to the range of <0.05 Hz, in order to avoid the Mayor wave which is shown to depend on body positions. Because of a high sampling frequency, fNIRS was effective in preventing aliasing of high-frequencies related to pulse and respiration into the range of <0.05 Hz.
In addition, our results revealed that the power spectrum of HbO global signal depends on body positions in the range between 0.05 – 0.1 Hz, regardless eyes were opened and closed (shown in Figures 2, 3 at respective individual and group level). Data at the standing position show the largest amplitude than the others, while the supine position is associated with lowest amplitude. These findings are consistent with previous reports by Tachtsidis et al. (2004), who compared three different positions’ effect on cerebral blood pressure with fNIRS. Their results showed that standing position has the highest mean blood pressure (MBP) and supine has the lowest MBP. They followed the Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology (1994) to separate the frequency spectrum into three standard frequency bands: very low frequency (VLF: 0.02–0.04 Hz), low frequency (LF: 0.04–0.15 Hz) and high frequency (HF: 0.15–0.4 Hz). Although VLF did not reveal any significant impact of body position, their results reported that the magnitude of low frequency oscillation in HbO in the resting brain shows a significant difference between different postures in LF. Coincidentally, Mayer wave, i.e., the cyclic changes in arterial blood pressure, fall into this LF range (Muller et al., 2003; Julien, 2006). Mayer wave appears to have a close relationship with fNIRS global signal. It is observed as oscillations of arterial pressure at ∼0.1 Hz in conscious humans (Julien, 2006). Besides, it is positively related with the strength of the corresponding sympathetic nervous activity and the mean level of sympathetic nerve activity (Furlan et al., 2000). More importantly, prone, supine, and sitting have significantly different effect on autonomic regulation of cardiovascular function (Watanabe et al., 2007). One rational speculation is that different body positions, especially the up-tilt positions, significantly affect autonomic regulation includes SNA which set the level of sympathetic vasoconstrictor tone, hence contributing to sustain arterial pressure (Julien, 2006; Scholkmann et al., 2014; Mohammadi-Nejad et al., 2018). Therefore, we regarded position-dependent effect in the Mayer wave range to be of physiological origin and discarded them for comparison against EEG. Aside from the Mayer wave range, our analysis further eliminated the factor of body positions and revealed a negative association between the EEG vigilance and fNIRS global signal in the frequency range of <0.05 Hz (HbO and HbR in Figures 6, 7). Such EEG-fNIRS association for the first time revealed a neurophysiological contribution to the fluctuations of fNIRS global signal (due to EEG vigilance), rather than a physiological factor (due to body positions). As control data in the Experiment 1, we conducted qualitative analysis and statistical analysis on auditory EEG and fNIRS responses. Our results did not observe the different body positions’ effect on AEP of EEG data or task-related average of fNIRS data, at both EO and EC conditions. This excludes the concerns of environmental and systematic biases, such as the quality of data recording when subjects were positioned differently.
Our findings of a negative association between fNIRS global signal and EEG vigilance measures have important implications for the analysis and interpretation of fNIRS-based resting state functional connectivity. The reversed association between EEG vigilance and fNIRS global signal observed within individuals (Figure 6) indicates that removal of fNIRS global signal will also remove the neurological effect of vigilance in the signals. Therefore, in resting state functional connectivity studies using repeated fNIRS measures within an individual, especially if the subjects’ conditions are related to the vigilance levels, global signal removal should not be performed. In addition, a negative covariation across individuals at the standing position (Figure 7) indicates that removal of global signal will also remove the neurological effect of vigilance in the signals. Therefore, in resting state studies using fNIRS (such as biomarkers of disease across individuals), especially when vigilance is an individual-level trait relevant to disease symptoms (Yang et al., 2017) or behaviors (Li et al., 2019), global signal removal could become problematic.
Furthermore, our study in the Experiment 2 for the first time reported a negative temporal correlation between the epoch-by-epoch fluctuations of fNIRS global signal and EEG vigilance, which further confirmed the negative association observed across subjects in Experiment 1. In wakeful rest periods that were verified by sleep scoring criteria, subjects exhibited momentary upsurges and drops of vigilance, as they stayed awake but were falling into sleep. Such vigilance fluctuations were then shown to be in reversed synchrony with the HbO and HbR global signal: epochs of decreased vigilance were associated with surges in fNIRS global signal, and vice versa. Interestingly, such findings of negative association are consistent with other studies that have examined the BOLD fMRI signals and EEG or behavior measures of vigilance. For example, Olbrich et al. (2009) of simultaneous EEG and fMRI study have reported that decreased of EEG vigilance measures are associated with a perfusive increase of BOLD signals in widespread cortices (except the thalamus). Similarly, Liu and colleges have observed a negative temporal correlation between EEG vigilance and fMRI global signal calculated from whole-brain average, which was reported in 23 sessions out of the 25 sessions in total and ranged between 0 and −0.5 (Falahpour et al., 2018, Figure 4), which is similar with our observation. In another study performed in unanesthetized monkeys, Chang et al. (2016) investigated the behavior measure of vigilance, indicated as opening and closure of eyes, and reported again that the fluctuations of vigilance level have negative temporal correlation with BOLD signals in widespread cortices, in a similar spatial extent and consistent temporal manner with those observed in human studies (Olbrich et al., 2009; Falahpour et al., 2018). Our study, however, reported for the first time the negative temporal correlation exists in concurrent and whole-head fNIRS and EEG recordings in human. Our findings of the epoch-by-epoch association are important for the interpretation of dynamic resting state functional connectivity, which commonly used a windowed approach of 30-s to 120-s duration. While the resting state brain connectivity is increasingly recognized to possess rich information of dynamics (Hutchison et al., 2013), some studies removed the global signal (Allen et al., 2014) whereas other studies did not (Chang et al., 2013). By showing that the fNIRS global signal amplitude is negatively correlated with EEG vigilance, our findings suggest that the removal the global signal should not be performed in the investigation of dynamic connectivity using fNIRS, especially in conditions affected by vigilance fluctuations. Beyond that, DMN has been reported to be correlated with EEG vigilance scores (Olbrich et al., 2009). Removing fNIRS global signal therefore may attenuate activities of DMN that are related with vigilance fluctuations. Evidence has shown that the working memory plays a critical role in both visual rehearsal and vigilance performance (Baddeley et al., 1999). And age-related alterations and disease-related decrements (such as Alzheimer’s disease) in DMN have significantly impacted working memory performance (Baddeley et al., 1999; Sambataro et al., 2010). Therefore, the fNIRS global signal should not be treated as non-neural confound, and its removal should be carefully considered via a frequency delineation.
Noteworthy, the calculation of the fNIRS global signal amplitude in our study is a reasonable adaption from the definition of global signal amplitude in previous fMRI study (Wong et al., 2013). Considering that the fNIRS optical density is converted to relative changes of HbO/HbR concentration in the stage of hemodynamic computation, the normalization in fNIRS equates the normalization in fMRI analysis (i.e., divided by the mean of fMRI time course), the calculation of fNIRS global signal in our study followed exactly the same definition in Wong et al. (2013). Our findings are consistent with previous findings on the relationship between fMRI global signal and EEG vigilance (Wong et al., 2013, 2016). Such discovery of a neurological component in fNIRS global in our study is novel. Importantly, our investigation adds findings from a unique perspective by showing a covariation relationship in a carefully constrained frequency range that has excluded the possible physiological noise of blood pressure regulation. Our studies of two experiment datasets have demonstrated the reversed association exiting in both static and dynamic manner.
Additionally, it is worthy to note that the quantification of EEG vigilance has limitations, due to interindividual variance in EEG activity that is commonly observed in many EEG studies (Jobert et al., 1994; de Munck et al., 2007; Olbrich et al., 2009). For example, certain subjects may exhibit almost no EEG alpha peak in the power spectrum of resting EEG at eyes open and sometimes, even at eyes closed states. In the meanwhile, there are subjects that show strong alpha peak power at both eyes open and eyes closed states. Here in our analysis we have taken multiple steps to mitigate the factor of interindividual variance in EEG activity. Firstly, we have calculated a normalized spectrum of EEG accounting all frequency bins; and the vigilance measures was calculated as the ration between the alpha-band amplitude divided by the sum of amplitudes in the delta and theta bands based on the normalized spectrum. Then, in Experiment 1, we accounted the interindividual variance by contrasting between the lowest vigilance condition and higher vigilance conditions within individuals, at which subjects’ eyes were all open. Moreover, in Experiment 2, we examined the temporal fluctuations to determine the association between EEG vigilance and fNIRS, while we concluded that the levels of vigilance fluctuations did not affect the temporal association. Nonetheless, because we also examined the across-individual covariation (Figure 7), the observation of the association between EEG vigilance and fNIRS global signal could be attributed to interindividual variance in EEG activities rather than the neurological factor of vigilance, although removal the global signal under such situation could still introduce confounds to the fNIRS resting state functional connectivity analysis.
Conclusion
With the advantage of economic efficiency and portability, fNIRS has been proposed as a complementary option to fMRI, especially to be used in populations with contraindications to MRI scanner and in challenged environment (such as brain monitoring at bed-side or during surgery). The current study for the first time revealed a negative relationship between fNIRS global signal amplitudes and EEG vigilance in human participants, based on concurrent EEG and fNIRS recordings at high-density and whole-head montage. Our results stressed the significant effect of body positions on the fNIRS resting-state global signal, primarily in the frequency range of greater than 0.05 Hz yet not in the range of less than 0.05 Hz. However, EEG vigilance plays a modulatory role in the fNIRS signals in the frequency range of less than 0.05 Hz: resting-state sessions of low EEG vigilance measures are associated with high amplitudes of fNIRS global signals. Moreover, the epoch-by-epoch fluctuations of EEG vigilance and fNIRS global signals are significantly correlated in a negative manner at a wakeful resting period. The findings of a neural component, i.e., EEG vigilance, in fNIRS global signal suggests that such global signal should not be removed as non-neural physiological signal, especially in studies and conditions where vigilance and related brain networks are of interest.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics Statement
The studies involving human participants were reviewed and approved by the Institutional Review Board at the University of Oklahoma Health Sciences Center. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
YC contributed to data collection, the data analysis, and manuscript writing. JT, YC, and MC contributed to data collection. JF and BC contributed to data analysis. HY contributed to study design, the data analysis, and the manuscript writing. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by NSF RII Track-2 FEC 1539068, Oklahoma Center for the Advancement of Science and Technology HR16-057, and Institute for Biomedical Engineering, Science, and Technology at University of Oklahoma.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The authors would like to thank Dr. Lei Ding, Dr. Guofa Shou, Johnny O’Keeffe, Josiah Rippetoe, and Junwei Ma for their assistance in the data collection.
References
Aguirre, G. K., Zarahn, E., and D’Esposito, M. (1998). The variability of human. BOLD hemodynamic responses. Neuroimage 8, 360–369. doi: 10.1006/nimg.1998.0369
Ahn, S., and Jun, S. C. (2017). Multi-Modal Integration of EEG-fNIRS for brain-computer interfaces - current limitations and future directions. Front. Hum. Neurosci. 11:503. doi: 10.3389/fnhum.2017.00503
Allen, E. A., Damaraju, E., Plis, S. M., Erhardt, E. B., Eichele, T., and Calhoun, V. D. (2014). Tracking whole-brain connectivity dynamics in the resting state. Cereb. Cortex 24, 663–676. doi: 10.1093/cercor/bhs352
Arun, K. M., Smitha, K. A., Sylaja, P. N., and Kesavadas, C. (2020). Identifying resting-state functional connectivity changes in the motor cortex using fNIRS during recovery from stroke. Brain Topogr. 33, 710–719. doi: 10.1007/s10548-020-00785-782
Baddeley, A., Cocchini, G., Della Sala, S., Logie, R. H., and Spinnler, H. (1999). Working memory and vigilance: evidence from normal aging and Alzheimer’s disease. Brain Cogn. 41, 87–108. doi: 10.1006/brcg.1999.1097
Beckmann, C. F., DeLuca, M., Devlin, J. T., and Smith, S. M. (2005). Investigations into resting-state connectivity using independent component analysis. Philos. Trans. R. Soc. Lond. B Biol. Sci. 360, 1001–1013. doi: 10.1098/rstb.2005.1634
Berry, R. B., Brooks, R., Gamaldo, C., Harding, S. M., Lloyd, R. M., Quan, S. F., et al. (2017). AASM scoring manual updates for 2017 (version 2.4). J. Clin. Sleep Med. 13, 665–666. doi: 10.5664/jcsm.6576
Birn, R. M., Diamond, J. B., Smith, M. A., and Bandettini, P. A. (2006). Separating respiratory-variation-related fluctuations from neuronal-activity-related fluctuations in fMRI. Neuroimage 31, 1536–1548. doi: 10.1016/j.neuroimage.2006.02.048
Biswal, B., Yetkin, F. Z., Haughton, V. M., and Hyde, J. S. (1995). Functional connectivity in the motor cortex of resting human brain using echo-planar MRI. Magn. Reson. Med. 34, 537–541. doi: 10.1002/mrm.1910340409
Brookes, M. J., Woolrich, M., Luckhoo, H., Price, D., Hale, J. R., Stephenson, M. C., et al. (2011). Investigating the electrophysiological basis of resting state networks using magnetoencephalography. Proc. Natl. Acad. Sci. U.S.A. 108, 16783–16788. doi: 10.1073/pnas.1112685108
Bulgarelli, C., Blasi, A., de Klerk, C., Richards, J. E., Hamilton, A., and Southgate, V. (2019). Fronto-temporoparietal connectivity and self-awareness in 18-month-olds: a resting state fNIRS study. Dev. Cogn. Neurosci. 38:100676. doi: 10.1016/j.dcn.2019.100676
Bulgarelli, C., de Klerk, C., Richards, J. E., Southgate, V., Hamilton, A., and Blasi, A. (2020). The developmental trajectory of fronto-temporoparietal connectivity as a proxy of the default mode network: a longitudinal fNIRS investigation. Hum. Brain Mapp. 41, 2717–2740. doi: 10.1002/hbm.24974
Chang, C., Leopold, D. A., Scholvinck, M. L., Mandelkow, H., Picchioni, D., Liu, X., et al. (2016). Tracking brain arousal fluctuations with fMRI. Proc. Natl. Acad. Sci. U.S.A. 113, 4518–4523. doi: 10.1073/pnas.1520613113
Chang, C., Liu, Z., Chen, M. C., Liu, X., and Duyn, J. H. (2013). EEG correlates of time-varying BOLD functional connectivity. NeuroImage 72, 227–236. doi: 10.1016/j.neuroimage.2013.01.049
Chaumon, M., Bishop, D. V., and Busch, N. A. (2015). A practical guide to the selection of independent components of the electroencephalogram for artifact correction. J. Neurosci. Methods 250, 47–63. doi: 10.1016/j.jneumeth.2015.02.025
Chen, G., Chen, G., Xie, C., Ward, B. D., Li, W., Antuono, P., et al. (2012). A method to determine the necessity for global signal regression in resting-state fMRI studies. Magn. Reson. Med. 68, 1828–1835. doi: 10.1002/mrm.24201
Chen, W. L., Wagner, J., Heugel, N., Sugar, J., Lee, Y. W., Conant, L., et al. (2020). Functional near-infrared spectroscopy and its clinical application in the field of neuroscience: advances and future directions. Front. Neurosci. 14:724. doi: 10.3389/fnins.2020.00724
Cordes, D., Haughton, V. M., Arfanakis, K., Carew, J. D., Turski, P. A., Moritz, C. H., et al. (2001). Frequencies contributing to functional connectivity in the cerebral cortex in “resting-state” data. AJNR Am. J. Neuroradiol. 22, 1326–1333.
Damoiseaux, J. S., Rombouts, S. A., Barkhof, F., Scheltens, P., Stam, C. J., Smith, S. M., et al. (2006). Consistent resting-state networks across healthy subjects. Proc. Natl. Acad. Sci. U.S.A. 103, 13848–13853. doi: 10.1073/pnas.0601417103
de Munck, J. C., Gonçalves, S. I, Huijboom, L., Kuijer, J. P. A., Pouwels, P. J. W., Heethaar, R. M., et al. (2007). The hemodynamic response of the alpha rhythm: an EEG/fMRI study. NeuroImage 35, 1142–1151. doi: 10.1016/j.neuroimage.2007.01.022
de Souza Rodrigues, J., Ribeiro, F. L., Sato, J. R., Mesquita, R. C., and Júnior, C. E. B. (2019). Identifying individuals using fNIRS-based cortical connectomes. Biomed. Opt. Express 10, 2889–2897. doi: 10.1364/boe.10.002889
Delorme, A., and Makeig, S. (2004). EEGLAB: an open source toolbox for analysis of single-trial EEG dynamics including independent component analysis. J. Neurosci. Methods 134, 9–21. doi: 10.1016/j.jneumeth.2003.10.009
Desjardins, A. E., Kiehl, K. A., and Liddle, P. F. (2001). Removal of confounding effects of global signal in functional MRI analyses. Neuroimage 13, 751–758. doi: 10.1006/nimg.2000.0719
Duan, L., Zhao, Z., Lin, Y., Wu, X., Luo, Y., and Xu, P. (2018). Wavelet-based method for removing global physiological noise in functional near-infrared spectroscopy. Biomed. Opt. Express 9, 3805–3820. doi: 10.1364/BOE.9.003805
Eggebrecht, A. T., Ferradal, S. L., Robichaux-Viehoever, A., Hassanpour, M. S., Dehghani, H., Snyder, A. Z., et al. (2014). Mapping distributed brain function and networks with diffuse optical tomography. Nat. Photonics 8, 448–454. doi: 10.1038/nphoton.2014.107
Eggebrecht, A. T., White, B. R., Ferradal, S. L., Chen, C., Zhan, Y., Snyder, A. Z., et al. (2012). A quantitative spatial comparison of high-density diffuse optical tomography and fMRI cortical mapping. Neuroimage 61, 1120–1128. doi: 10.1016/j.neuroimage.2012.01.124
Falahpour, M., Chang, C., Wong, C. W., and Liu, T. T. (2018). Template-based prediction of vigilance fluctuations in resting-state fMRI. Neuroimage 174, 317–327. doi: 10.1016/j.neuroimage.2018.03.012
Fox, M. D., Snyder, A. Z., Vincent, J. L., Corbetta, M., Van Essen, D. C., and Raichle, M. E. (2005). The human brain is intrinsically organized into dynamic, anticorrelated functional networks. Proc. Natl. Acad. Sci. U.S.A. 102, 9673–9678. doi: 10.1073/pnas.0504136102
Fox, M. D., Zhang, D., Snyder, A. Z., and Raichle, M. E. (2009). The global signal and observed anticorrelated resting state brain networks. J. Neurophysiol. 101, 3270–3283. doi: 10.1152/jn.90777.2008
Fox, P. T., Mintun, M. A., Reiman, E. M., and Raichle, M. E. (1988). Enhanced detection of focal brain responses using intersubject averaging and change-distribution analysis of subtracted PET images. J. Cereb. Blood Flow Metab. 8, 642–653. doi: 10.1038/jcbfm.1988.111
Franceschini, M. A., Joseph, D. K., Huppert, T. J., Diamond, S. G., and Boas, D. A. (2006). Diffuse optical imaging of the whole head. J. Biomed. Opt. 11:054007. doi: 10.1117/1.2363365
Friston, K. J., Frith, C. D., Liddle, P. F., Dolan, R. J., Lammertsma, A. A., and Frackowiak, R. S. (1990). The relationship between global and local changes in PET scans. J. Cereb. Blood Flow Metab. 10, 458–466. doi: 10.1038/jcbfm.1990.88
Furlan, R., Porta, A., Costa, F., Tank, J., Baker, L., Schiavi, R., et al. (2000). Oscillatory patterns in sympathetic neural discharge and cardiovascular variables during orthostatic stimulus. Circulation 101, 886–892. doi: 10.1161/01.cir.101.8.886
Gagnon, L., Perdue, K., Greve, D. N., Goldenholz, D., Kaskhedikar, G., and Boas, D. A. (2011). Improved recovery of the hemodynamic response in diffuse optical imaging using short optode separations and state-space modeling. Neuroimage 56, 1362–1371. doi: 10.1016/j.neuroimage.2011.03.001
Glover, G. H., Li, T. Q., and Ress, D. (2000). Image-based method for retrospective correction of physiological motion effects in fMRI: RETROICOR. Magn. Reson. Med. 44, 162–167. doi: 10.1002/1522-2594(200007)44:1<162::aid-mrm23>3.0.co;2-e
Goldman, R. I., Stern, J. M., Engel, J. Jr., and Cohen, M. S. (2002). Simultaneous EEG and fMRI of the alpha rhythm. Neuroreport 13, 2487–2492. doi: 10.1097/01.wnr.0000047685.08940.d0
Gratzer, W. B., and Allison, A. C. (1960). Multiple haemoglobins. Biol. Rev. 35, 459–503. doi: 10.1111/j.1469-185X.1960.tb01523.x
Gregg, N. M., White, B. R., Zeff, B. W., Berger, A. J., and Culver, J. P. (2010). Brain specificity of diffuse optical imaging: improvements from superficial signal regression and tomography. Front. Neuroenerget. 2:14. doi: 10.3389/fnene.2010.00014
Greicius, M. D., Krasnow, B., Reiss, A. L., and Menon, V. (2003). Functional connectivity in the resting brain: a network analysis of the default mode hypothesis. Proc. Natl. Acad. Sci. U.S.A. 100, 253–258. doi: 10.1073/pnas.0135058100
Hahamy, A., Calhoun, V., Pearlson, G., Harel, M., Stern, N., Attar, F., et al. (2014). Save the global: global signal connectivity as a tool for studying clinical populations with functional magnetic resonance imaging. Brain Connect. 4, 395–403. doi: 10.1089/brain.2014.0244
He, B. J., Snyder, A. Z., Zempel, J. M., Smyth, M. D., and Raichle, M. E. (2008). Electrophysiological correlates of the brain’s intrinsic large-scale functional architecture. Proc. Natl. Acad. Sci. U.S.A. 105, 16039–16044. doi: 10.1073/pnas.0807010105
Homae, F., Watanabe, H., Otobe, T., Nakano, T., Go, T., Konishi, Y., et al. (2010). Development of global cortical networks in early infancy. J. Neurosci. 30, 4877–4882. doi: 10.1523/JNEUROSCI.5618-09.2010
Horovitz, S. G., Fukunaga, M., de Zwart, J. A., van Gelderen, P., Fulton, S. C., Balkin, T. J., et al. (2008). Low frequency BOLD fluctuations during resting wakefulness and light sleep: a simultaneous EEG-fMRI study. Hum. Brain Mapp. 29, 671–682. doi: 10.1002/hbm.20428
Hu, Z., Liu, G., Dong, Q., and Niu, H. (2020). Applications of resting-state fNIRS in the developing brain: a review from the connectome perspective. Front. Neurosci. 14:476. doi: 10.3389/fnins.2020.00476
Huppert, T. J., Diamond, S. G., Franceschini, M. A., and Boas, D. A. (2009). HomER: a review of time-series analysis methods for near-infrared spectroscopy of the brain. Appl. Opt. 48, D280–D298.
Hutchison, R. M., Womelsdorf, T., Allen, E. A., Bandettini, P. A., Calhoun, V. D., Corbetta, M., et al. (2013). Dynamic functional connectivity: promise, issues, and interpretations. NeuroImage 80, 360–378. doi: 10.1016/j.neuroimage.2013.05.079
Ieong, H. F., Gao, F., and Yuan, Z. (2019). Machine learning: assessing neurovascular signals in the prefrontal cortex with non-invasive bimodal electro-optical neuroimaging in opiate addiction. Sci. Rep. 9:18262. doi: 10.1038/s41598-019-54316-54316
Jasper, H. H. (1958). The ten-twenty electrode system of the international federation. Electroencephalogr. Clin. Neurophysiol. 10, 371–375.
Jobert, M., Schulz, H., Jahnig, P., Tismer, C., Bes, F., and Escola, H. (1994). A computerized method for detecting episodes of wakefulness during sleep based on the alpha slow-wave index (ASI). Sleep 17, 37–46.
Julien, C. (2006). The enigma of Mayer waves: facts and models. Cardiovasc. Res. 70, 12–21. doi: 10.1016/j.cardiores.2005.11.008
Kay, K. N., Rokem, A., Winawer, J., Dougherty, R. F., and Wandell, B. A. (2013). GLMdenoise: a fast, automated technique for denoising task-based fMRI data. Front. Neurosci. 7:247. doi: 10.3389/fnins.2013.00247
Keehn, B., Wagner, J. B., Tager-Flusberg, H., and Nelson, C. A. (2013). Functional connectivity in the first year of life in infants at-risk for autism: a preliminary near-infrared spectroscopy study. Front. Hum. Neurosci. 7:444. doi: 10.3389/fnhum.2013.00444
Kocsis, L., Herman, P., and Eke, A. (2006). The modified Beer-Lambert law revisited. Phys. Med. Biol. 51, N91–N98. doi: 10.1088/0031-9155/51/5/N02
Kohno, S., Miyai, I., Seiyama, A., Oda, I., Ishikawa, A., Tsuneishi, S., et al. (2007). Removal of the skin blood flow artifact in functional near-infrared spectroscopic imaging data through independent component analysis. J. Biomed. Opt. 12:062111. doi: 10.1117/1.2814249
Leopold, D. A., Murayama, Y., and Logothetis, N. K. (2003). Very slow activity fluctuations in monkey visual cortex: implications for functional brain imaging. Cereb. Cortex 13, 422–433. doi: 10.1093/cercor/13.4.422
Li, J., Bolt, T., Bzdok, D., Nomi, J. S., Yeo, B. T. T., Spreng, R. N., et al. (2019). Topography and behavioral relevance of the global signal in the human brain. Sci. Rep. 9:14286. doi: 10.1038/s41598-019-50750-50758
Lu, C.-M., Zhang, Y.-J., Biswal, B. B., Zang, Y.-F., Peng, D.-L., and Zhu, C.-Z. (2010). Use of fNIRS to assess resting state functional connectivity. J. Neurosci. Methods 186, 242–249. doi: 10.1016/j.jneumeth.2009.11.010
Macey, P. M., Macey, K. E., Kumar, R., and Harper, R. M. (2004). A method for removal of global effects from fMRI time series. Neuroimage 22, 360–366. doi: 10.1016/j.neuroimage.2003.12.042
Mantini, D., Perrucci, M. G., Del Gratta, C., Romani, G. L., and Corbetta, M. (2007). Electrophysiological signatures of resting state networks in the human brain. Proc. Natl. Acad. Sci. U.S.A. 104, 13170–13175. doi: 10.1073/pnas.0700668104
Mesquita, R. C., Franceschini, M. A., and Boas, D. A. (2010). Resting state functional connectivity of the whole head with near-infrared spectroscopy. Biomed. Opt. Express 1, 324–336. doi: 10.1364/BOE.1.000324
Mohammadi-Nejad, A. R., Mahmoudzadeh, M., Hassanpour, M. S., Wallois, F., Muzik, O., Papadelis, C., et al. (2018). Neonatal brain resting-state functional connectivity imaging modalities. Photoacoustics 10, 1–19. doi: 10.1016/j.pacs.2018.01.003
Molavi, B., May, L., Gervain, J., Carreiras, M., Werker, J. F., and Dumont, G. A. (2013). Analyzing the resting state functional connectivity in the human language system using near infrared spectroscopy. Front. Hum. Neurosci. 7:921. doi: 10.3389/fnhum.2013.00921
Muller, T., Reinhard, M., Oehm, E., Hetzel, A., and Timmer, J. (2003). Detection of very low-frequency oscillations of cerebral haemodynamics is influenced by data detrending. Med. Biol. Eng. Comput. 41, 69–74. doi: 10.1007/bf02343541
Murphy, K., Birn, R. M., Handwerker, D. A., Jones, T. B., and Bandettini, P. A. (2009). The impact of global signal regression on resting state correlations: are anti-correlated networks introduced? Neuroimage 44, 893–905. doi: 10.1016/j.neuroimage.2008.09.036
Murphy, K., and Fox, M. D. (2017). Towards a consensus regarding global signal regression for resting state functional connectivity MRI. Neuroimage 154, 169–173. doi: 10.1016/j.neuroimage.2016.11.052
Naseer, N., and Hong, K. S. (2015). fNIRS-based brain-computer interfaces: a review. Front. Hum. Neurosci. 9:3. doi: 10.3389/fnhum.2015.00003
Novi, S. L., Rodrigues, R. B., and Mesquita, R. C. (2016). Resting state connectivity patterns with near-infrared spectroscopy data of the whole head. Biomed. Opt. Express 7, 2524–2537. doi: 10.1364/BOE.7.002524
Obrig, H., Neufang, M., Wenzel, R., Kohl, M., Steinbrink, J., Einhaupl, K., et al. (2000). Spontaneous low frequency oscillations of cerebral hemodynamics and metabolism in human adults. Neuroimage 12, 623–639. doi: 10.1006/nimg.2000.0657
Obrig, H., and Villringer, A. (2003). Beyond the visible–imaging the human brain with light. J. Cereb. Blood Flow Metab. 23, 1–18. doi: 10.1097/01.WCB.0000043472.45775.29
Olbrich, S., Mulert, C., Karch, S., Trenner, M., Leicht, G., Pogarell, O., et al. (2009). EEG-vigilance and BOLD effect during simultaneous EEG/fMRI measurement. Neuroimage 45, 319–332. doi: 10.1016/j.neuroimage.2008.11.014
Pinti, P., Tachtsidis, I., Hamilton, A., Hirsch, J., Aichelburg, C., Gilbert, S., et al. (2018). The present and future use of functional near-infrared spectroscopy (fNIRS) for cognitive neuroscience. Ann. N. Y. Acad. Sci. 1464, 5–29. doi: 10.1111/nyas.13948
Power, J. D., Schlaggar, B. L., and Petersen, S. E. (2015). Recent progress and outstanding issues in motion correction in resting state fMRI. Neuroimage 105, 536–551. doi: 10.1016/j.neuroimage.2014.10.044
Putze, F., Hesslinger, S., Tse, C.-Y., Huang, Y., Herff, C., Guan, C., et al. (2014). Hybrid fNIRS-EEG based classification of auditory and visual perception processes. Front. Neurosci. 8:373–373. doi: 10.3389/fnins.2014.00373
Raichle, M. E., MacLeod, A. M., Snyder, A. Z., Powers, W. J., Gusnard, D. A., and Shulman, G. L. (2001). A default mode of brain function. Proc. Natl. Acad. Sci. U.S.A. 98, 676–682. doi: 10.1073/pnas.98.2.676
Saad, Z. S., Gotts, S. J., Murphy, K., Chen, G., Jo, H. J., Martin, A., et al. (2012). Trouble at rest: how correlation patterns and group differences become distorted after global signal regression. Brain Connect. 2, 25–32. doi: 10.1089/brain.2012.0080
Saager, R. B., and Berger, A. J. (2005). Direct characterization and removal of interfering absorption trends in two-layer turbid media. J. Opt. Soc. Am. A Opt. Image Sci. Vis. 22, 1874–1882. doi: 10.1364/josaa.22.001874
Sakakibara, E., Homae, F., Kawasaki, S., Nishimura, Y., Takizawa, R., Koike, S., et al. (2016). Detection of resting state functional connectivity using partial correlation analysis: a study using multi-distance and whole-head probe near-infrared spectroscopy. Neuroimage 142, 590–601. doi: 10.1016/j.neuroimage.2016.08.011
Sambataro, F., Murty, V. P., Callicott, J. H., Tan, H. Y., Das, S., Weinberger, D. R., et al. (2010). Age-related alterations in default mode network: impact on working memory performance. Neurobiol. Aging 31, 839–852. doi: 10.1016/j.neurobiolaging.2008.05.022
Santosa, H., Zhai, X., Fishburn, F., and Huppert, T. (2018). The NIRS Brain AnalyzIR Toolbox. Algorithms 11:73. doi: 10.3390/a11050073
Sasai, S., Homae, F., Watanabe, H., and Taga, G. (2011). Frequency-specific functional connectivity in the brain during resting state revealed by NIRS. Neuroimage 56, 252–257. doi: 10.1016/j.neuroimage.2010.12.075
Sato, T., Nambu, I., Takeda, K., Aihara, T., Yamashita, O., Isogaya, Y., et al. (2016). Reduction of global interference of scalp-hemodynamics in functional near-infrared spectroscopy using short distance probes. Neuroimage 141, 120–132. doi: 10.1016/j.neuroimage.2016.06.054
Scholkmann, F., Kleiser, S., Metz, A. J., Zimmermann, R., Mata Pavia, J., Wolf, U., et al. (2014). A review on continuous wave functional near-infrared spectroscopy and imaging instrumentation and methodology. Neuroimage 85(Pt 1), 6–27. doi: 10.1016/j.neuroimage.2013.05.004
Scholvinck, M. L., Maier, A., Ye, F. Q., Duyn, J. H., and Leopold, D. A. (2010). Neural basis of global resting-state fMRI activity. Proc. Natl. Acad. Sci. U.S.A. 107, 10238–10243. doi: 10.1073/pnas.0913110107
Sherafati, A., Snyder, A. Z., Eggebrecht, A. T., Bergonzi, K. M., Burns-Yocum, T. M., Lugar, H. M., et al. (2020). Global motion detection and censoring in high-density diffuse optical tomography. Hum. Brain Mapp. doi: 10.1002/hbm.25111 [Online ahead of print]
Shin, J., Kwon, J., Choi, J., and Im, C. H. (2017). Performance enhancement of a brain-computer interface using high-density multi-distance NIRS. Sci. Rep. 7:16545. doi: 10.1038/s41598-017-16639-16630
Shmuel, A., and Leopold, D. A. (2008). Neuronal correlates of spontaneous fluctuations in fMRI signals in monkey visual cortex: Implications for functional connectivity at rest. Hum. Brain Mapp. 29, 751–761. doi: 10.1002/hbm.20580
Sood, B. G., McLaughlin, K., and Cortez, J. (2015). Near-infrared spectroscopy: applications in neonates. Semin. Fetal Neonatal Med. 20, 164–172. doi: 10.1016/j.siny.2015.03.008
Tachtsidis, I., Elwell, C. E., Leung, T. S., Lee, C. W., Smith, M., and Delpy, D. T. (2004). Investigation of cerebral haemodynamics by near-infrared spectroscopy in young healthy volunteers reveals posture-dependent spontaneous oscillations. Physiol. Meas. 25, 437–445. doi: 10.1088/0967-3334/25/2/003
Tachtsidis, I., and Scholkmann, F. (2016). False positives and false negatives in functional near-infrared spectroscopy: issues, challenges, and the way forward. Neurophotonics 3:031405. doi: 10.1117/1.NPh.3.3.031405
Tong, Y., and Frederick, B. D. (2010). Time lag dependent multimodal processing of concurrent fMRI and near-infrared spectroscopy (NIRS) data suggests a global circulatory origin for low-frequency oscillation signals in human brain. Neuroimage 53, 553–564. doi: 10.1016/j.neuroimage.2010.06.049
Torricelli, A., Contini, D., Pifferi, A., Caffini, M., Re, R., Zucchelli, L., et al. (2014). Time domain functional NIRS imaging for human brain mapping. Neuroimage 85(Pt 1), 28–50. doi: 10.1016/j.neuroimage.2013.05.106
Wang, M., Hu, Z., Liu, L., Li, H., Qian, Q., and Niu, H. (2020). Disrupted functional brain connectivity networks in children with attention-deficit/hyperactivity disorder: evidence from resting-state functional near-infrared spectroscopy. Neurophotonics 7:015012. doi: 10.1117/1.NPh.7.1.015012
Watabe, T., and Hatazawa, J. (2019). Evaluation of functional connectivity in the brain using positron emission tomography: a mini-review. Front. Neurosci. 13:775. doi: 10.3389/fnins.2019.00775
Watanabe, H., Shitara, Y., Aoki, Y., Inoue, T., Tsuchida, S., Takahashi, N., et al. (2017). Hemoglobin phase of oxygenation and deoxygenation in early brain development measured using fNIRS. Proc. Natl. Acad. Sci. U.S.A. 114, E1737–E1744. doi: 10.1073/pnas.1616866114
Watanabe, N., Reece, J., and Polus, B. I. (2007). Effects of body position on autonomic regulation of cardiovascular function in young, healthy adults. Chiropr. Osteopat. 15:19. doi: 10.1186/1746-1340-15-19
White, B. R., Liao, S. M., Ferradal, S. L., Inder, T. E., and Culver, J. P. (2012). Bedside optical imaging of occipital resting-state functional connectivity in neonates. Neuroimage 59, 2529–2538. doi: 10.1016/j.neuroimage.2011.08.094
White, B. R., Snyder, A. Z., Cohen, A. L., Petersen, S. E., Raichle, M. E., Schlaggar, B. L., et al. (2009). Resting-state functional connectivity in the human brain revealed with diffuse optical tomography. Neuroimage 47, 148–156. doi: 10.1016/j.neuroimage.2009.03.058
Wise, R. G., Ide, K., Poulin, M. J., and Tracey, I. (2004). Resting fluctuations in arterial carbon dioxide induce significant low frequency variations in BOLD signal. Neuroimage 21, 1652–1664. doi: 10.1016/j.neuroimage.2003.11.025
Wong, C. W., DeYoung, P. N., and Liu, T. T. (2016). Differences in the resting-state fMRI global signal amplitude between the eyes open and eyes closed states are related to changes in EEG vigilance. Neuroimage 124(Pt A), 24–31. doi: 10.1016/j.neuroimage.2015.08.053
Wong, C. W., Olafsson, V., Tal, O., and Liu, T. T. (2013). The amplitude of the resting-state fMRI global signal is related to EEG vigilance measures. Neuroimage 83, 983–990. doi: 10.1016/j.neuroimage.2013.07.057
Wyser, D., Mattille, M., Wolf, M., Lambercy, O., Scholkmann, F., and Gassert, R. (2020). Short-channel regression in functional near-infrared spectroscopy is more effective when considering heterogeneous scalp hemodynamics. Neurophotonics 7:035011. doi: 10.1117/1.NPh.7.3.035011
Xu, Y., Graber, H. L., and Barbour, R. L. (2014). “nirsLAB: a Computing(Environment)for fNIRS Neuroimaging Data Analysis,” in Proceedings of the conference Biomedical Optics 2014 (Washington, DC: Optical Society of America).
Yang, G. J., Murray, J. D., Glasser, M., Pearlson, G. D., Krystal, J. H., Schleifer, C., et al. (2017). Altered Global Signal Topography in Schizophrenia. Cereb. Cortex 27, 5156–5169. doi: 10.1093/cercor/bhw297
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
Yuan, H., Ding, L., Zhu, M., Zotev, V., Phillips, R., and Bodurka, J. (2016). Reconstructing large-scale brain resting-state networks from high-resolution EEG: spatial and temporal comparisons with fMRI. Brain Connect. 6, 122–135. doi: 10.1089/brain.2014.0336
Yuan, H., Zotev, V., Phillips, R., and Bodurka, J. (2013). Correlated slow fluctuations in respiration, EEG, and BOLD fMRI. Neuroimage 79, 81–93. doi: 10.1016/j.neuroimage.2013.04.068
Yuan, H., Zotev, V., Phillips, R., Drevets, W. C., and Bodurka, J. (2012). Spatiotemporal dynamics of the brain at rest–exploring EEG microstates as electrophysiological signatures of BOLD resting state networks. Neuroimage 60, 2062–2072. doi: 10.1016/j.neuroimage.2012.02.031
Zarahn, E., Aguirre, G., and D’Esposito, M. (1997). A trial-based experimental design for fMRI. Neuroimage 6, 122–138. doi: 10.1006/nimg.1997.0279
Zeff, B. W., White, B. R., Dehghani, H., Schlaggar, B. L., and Culver, J. P. (2007). Retinotopic mapping of adult human visual cortex with high-density diffuse optical tomography. Proc. Natl. Acad. Sci. U.S.A. 104, 12169–12174. doi: 10.1073/pnas.0611266104
Zhang, D., and Raichle, M. E. (2010). Disease and the brain’s dark energy. Nat. Rev. Neurol. 6, 15–28. doi: 10.1038/nrneurol.2009.198
Zhang, H., Duan, L., Zhang, Y. J., Lu, C. M., Liu, H., and Zhu, C. Z. (2011). Test-retest assessment of independent component analysis-derived resting-state functional connectivity based on functional near-infrared spectroscopy. Neuroimage 55, 607–615. doi: 10.1016/j.neuroimage.2010.12.007
Zhang, H., Zhang, Y.-J., Lu, C.-M., Ma, S.-Y., Zang, Y.-F., and Zhu, C.-Z. (2010). Functional connectivity as revealed by independent component analysis of resting-state fNIRS measurements. Neuroimage 51, 1150–1161. doi: 10.1016/j.neuroimage.2010.02.080
Zhang, Q., Brown, E. N., and Strangman, G. E. (2007). Adaptive filtering to reduce global interference in evoked brain activity detection: a human subject case study. J. Biomed. Opt. 12:064009. doi: 10.1117/1.2804706
Zhang, Q., Strangman, G. E., and Ganis, G. (2009). Adaptive filtering to reduce global interference in non-invasive NIRS measures of brain activation: how well and when does it work? Neuroimage 45, 788–794. doi: 10.1016/j.neuroimage.2008.12.048
Zhang, X., Noah, J. A., and Hirsch, J. (2016). Separation of the global and local components in functional near-infrared spectroscopy signals using principal component spatial filtering. Neurophotonics 3:015004. doi: 10.1117/1.NPh.3.1.015004
Zhang, Y., Brooks, D. H., Franceschini, M. A., and Boas, D. A. (2005). Eigenvector-based spatial filtering for reduction of physiological interference in diffuse optical imaging. J. Biomed. Opt. 10:11014. doi: 10.1117/1.1852552
Keywords: electroencephalography, functional near-infrared spectroscopy, resting state, vigilance, global signal
Citation: Chen Y, Tang J, Chen Y, Farrand J, Craft MA, Carlson BW and Yuan H (2020) Amplitude of fNIRS Resting-State Global Signal Is Related to EEG Vigilance Measures: A Simultaneous fNIRS and EEG Study. Front. Neurosci. 14:560878. doi: 10.3389/fnins.2020.560878
Received: 10 May 2020; Accepted: 11 November 2020;
Published: 03 December 2020.
Edited by:
Waldemar Karwowski, University of Central Florida, United StatesReviewed by:
Vassiliy Tsytsarev, University of Maryland, College Park, United StatesAleksandra Dagmara Kawala-Sterniuk, Opole University of Technology, Poland
Copyright © 2020 Chen, Tang, Chen, Farrand, Craft, Carlson and Yuan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Han Yuan, hanyuan@ou.edu