Abstract
Emerging neuroimaging studies emphasize the dynamic organization of spontaneous brain activity in both human and non-human primates, even under anesthesia. In a recent study, we were able to characterize the heterogeneous architecture of intrinsic functional flexibility in the awake, resting human brain using time-resolved analysis and a probabilistic model. However, it is unknown whether this organizational principle is preserved in the anesthetized monkey brain, and how anesthesia affects dynamic and static measurements of spontaneous brain activity. To investigate these issues, we collected resting-state functional magnetic resonance imaging (fMRI) datasets from 178 awake humans and 11 anesthetized monkeys (all healthy). Our recently established method, a complexity measurement (i.e., Shannon entropy) of dynamic functional connectivity patterns of each brain region, was used to map the intrinsic functional flexibility across the cerebral cortex. To further explore the potential effects of anesthesia, we performed time series analysis and correlation analysis between dynamic and static measurements within awake human and anesthetized monkey brains, respectively. We observed a heterogeneous profile of intrinsic functional flexibility in the anesthetized monkey brain, which showed some similarities to that of awake humans (r = 0.30, p = 0.007). However, we found that brain activity in anesthetized monkeys generally shifted toward random fluctuations. Moreover, there is a negative correlation between nodal entropy for the distribution of dynamic functional connectivity patterns and static functional connectivity strength in anesthetized monkeys, but not in awake humans. Our findings indicate that the heterogeneous architecture of intrinsic functional flexibility across cortex probably reflects an evolutionarily conserved aspect of functional brain organization, which persists across levels of cognitive processing (states of consciousness). The coupling between nodal entropy for the distribution of dynamic functional connectivity patterns and static functional connectivity strength may serve as a potential signature of anesthesia. This study not only offers fresh insight into the evolution of brain functional architecture, but also advances our understanding of the dynamics of spontaneous brain activity.
Introduction
A fundamental goal of comparative neuroscience is to determine conservation and evolution-driven changes in functional brain organization between species. Functional magnetic resonance imaging (fMRI), a non-invasive technique, has been utilized to identify functionally homologous or unique areas across primate species based on specific experimental tasks (Nakahara et al., 2002; Vanduffel et al., 2002; Wang et al., 2015). This technique has the advantage of providing a direct cross-species comparison using a common physiological measurement, i.e., blood oxygen level-dependent (BOLD) signal. Owing to ease of implementation and the robustness of findings, resting-state fMRI, simply a period of recording of BOLD signal in the absence of any explicit tasks, has become an attractive tool for studying large-scale brain functional organization (Zhang and Raichle, 2010; Power et al., 2014). Resting-state functional connectivity (FC) is typically used to index the interregional coherence in spontaneous low-frequency fluctuations (e.g., 0.01∼0.1 Hz) of BOLD signals (; ). Through this approach, many functional brain networks such as the sensorimotor and default mode networks have been identified in both human and non-human primates (Vincent et al., 2007; Margulies et al., 2009; ; ; Mantini et al., 2013; Miranda-Dominguez et al., 2014; Neubert et al., 2014). In particular, the patterns of resting-state FC have showed some similarities across species and appear to transcend levels of consciousness, being present under anesthesia, see Raichle (2015) for a review. However, conventional resting-state FC is frequently evaluated in a time-averaged sense, under the potential assumption of stationary functional organization. Moreover, the origins and functional significance of resting-state FC patterns remain to be further understood.
Recently, emerging studies have emphasized the dynamic organization of brain function, suggesting that understanding brain function and dysfunction requires an integrated framework linking brain connectivity and brain dynamics (; ; ; ; ; ). Therefore, dynamic FC analysis, e.g., taking into account the temporal fluctuations of FC in different time windows of BOLD signals, has been proposed to characterize spontaneous brain activity (; Zalesky et al., 2014; ; ; ). Although a number of challenges in techniques and interpretation remain, time-resolved analysis allows researchers to extract more in-depth information about brain function than static FC analysis (; ; ; Preti et al., 2017).
Developing new analytic tools to describe spatiotemporal characteristics of resting-state FC patterns across species and states may provide deeper insight into functional organization of spontaneous brain activity. Based on dynamic FC analysis and clustering method, many discrete, reproducible functional states over the time of scan have been identified in both humans and monkeys (; ), see ; for reviews. In contrast, our recent work (Yin et al., 2016) focused on quantifying the flexibility of the connectivity pattern for each brain region over time, using a complexity measurement (i.e., Shannon entropy) and probabilistic model. Consistent with task-induced functional reconfigurations (; ), we revealed the heterogeneous organization of functional flexibility in the resting human brain. However, it is unknown whether this organizational principle is preserved in the anesthetized monkey brain, and how it corresponds with the human brain during wakeful rest. Although the spatiotemporal dynamics of brain activity have been demonstrated in the anesthetic state (; ; Zhang et al., 2018), few studies pay attention to the cross-species correspondence of brain-wide dynamic organizational structure. Moreover, how anesthesia affects dynamic and static measurements of spontaneous brain activity still needs to be clarified.
To investigate these issues, we collected resting-state fMRI datasets from 178 awake humans and 11 anesthetized monkeys (healthy subjects). Our recently established method (Yin et al., 2016), a complexity measurement (i.e., Shannon entropy) of dynamic FC patterns of each brain region, was used to map the intrinsic functional flexibility across the cerebral cortex. Brain regions with high entropy for the distribution of dynamic FC patterns indicate high functional flexibility, and vice versa. For comparison, we conducted another complexity measurement based on distribution of correlation values, which reflects the functional complexity of a system (Zhao et al., 2010; Zamora-Lopez et al., 2016). A temporal variability analysis as a measure of functional flexibility was also carried out (Mueller et al., 2013; Zhang et al., 2016). To further explore the potential effects of anesthesia on the functional organization of spontaneous brain activity, we performed time series analysis and correlation analysis between dynamic and static measurements within awake human and anesthetized monkey brains, respectively. We hypothesized that the heterogeneous organization of intrinsic functional flexibility in the awake human brain persists in the anesthetized monkey brain.
Materials and Methods
Participants
All experimental procedures for non-human primate research in this study were approved by the Institutional Animal Care and Use Committee at the Institute of Neuroscience and the Biomedical Research Ethics Committee, Shanghai Institutes for Biological Sciences, Chinese Academy of Sciences, and conformed to National Institutes of Health guidelines for the humane care and use of laboratory animals.
We recruited 11 wild-type monkeys (age 4.68 ± 0.46 years; weight 3.97 ± 1.36 kg; 7 female). In addition, 178 healthy human subjects (age 14.4 ± 3.6 years; 44 female; 35 subjects with eyes closed and 143 subjects with eyes open) were collected from the Autism Brain Imaging Data Exchange (ABIDE)1. We screened the human data based on demographic and diagnostic information provided in the ABIDE database (). The inclusion criteria in the present study are briefly described as follows: (1) right-handedness, (2) age between 7 and 22 years, (3) a full-scale IQ score greater than 70; and time resolution of fMRI data equal to 2 s (datasets from four sites fit this criterion, including NYU, YALE, TRINITY, and UM).
Monkey Data Acquisition
Magnetic resonance imaging images of monkeys were acquired at the Institute of Neuroscience on a 3T whole-body scanner (Trio; Siemens Healthcare, Erlangen, Germany) running with an enhanced gradient coil insert (AC88; 80 mT/m maximum gradient strength, 800 mT/m/s maximum slew rate). A custom-built 8-channel phased-array transceiver coil was used for animal imaging sessions. Whole-brain resting-state fMRI data were collected using a gradient-echo echo-planar imaging (EPI) sequence (TR = 2000 ms; TE = 29 ms; flip angle = 77°; slices = 32; matrix = 64 × 64; field of view = 96 mm × 96 mm; 1.5 mm × 1.5 mm in plane resolution; slice thickness = 2.5 mm; GRAPPA factor = 2). For each session, 5–10 runs were acquired and each run consisted of 200 functional volumes. A pair of gradient echo images (echo time: 4.22 and 6.68 ms) with the same orientation and resolution as EPI images were acquired to generate a field map for distortion correction of EPI images. High-resolution T1-weighted anatomical images were acquired using a MPRAGE sequence (TR = 2500 ms; TE = 3.12 ms; inversion time = 1100 ms; flip angle = 9°; acquisition voxel size = 0.5 mm × 0.5 mm × 0.5 mm; 144 sagittal slices). Six whole-brain anatomical volumes were acquired and further averaged for better brain segmentation.
For MRI scanning, animals were prepared and maintained in a stable brain state under light anesthesia. The animal preparation procedure was conducted in a manner similar to our previous work (Wang et al., 2013; ). Induction of anesthesia was achieved by intramuscular injection with ketamine (10 mg/kg, Gutian Pharma Co., Ltd., China) before MRI scanning sessions, supplemented with atropine sulfate (0.05 mg/kg, Shanghai Harvest Pharma Co., Ltd., China) to decrease bronchial and salivary secretions. After intubation, animals were ventilated with a mixture of isoflurane (2–2.5%, Lunan Pharma Co., Ltd., China) and oxygen via either a standard ventilator (CWE, Inc., Ardmore, PA, United States) outside the scanner room or an MRI-compatible ventilator (CWE Inc., Weston, WI, United States) inside the scanner room. Macaques were maintained with intermittent positive-pressure ventilation to ensure a constant respiration rate (25–35 breaths/min). The concentration of isoflurane was adjusted based on continuously monitored vital signs, including blood oxygenation, electrocardiogram (ECG), rectal temperature (Small Animal Instruments, Inc., Stony Brook, NY, United States), respiration rate and end-tidal CO2 (Smiths Medical ASD Inc., Dublin, OH, United States). Oxygen saturation was kept over 95% and body temperature was kept constant using a heated water blanket (Gaymar Industries Inc., Orchard Park, NY, United States). Lactated Ringer’s solution was given with a maximum rate of 10 ml/kg/h during the anesthesia process (). We removed the runs that showed erratic vital signs, image artifacts, as well as burst suppression according to the recordings of MRI-compatible electroencephalograph (Brain Products GmbH, Gilching, Germany) during functional data acquisition. In total, 99 runs were left for the final analyses. We treated each run independently following previous studies (; ).
Human Data Acquisition
Human MRI data were acquired from multiple sites with different parameters of pulse sequences (see text footnote 1). In the present study, one of the inclusion criteria was a time resolution of fMRI data equal to 2 s, to match the monkey data. In addition, we kept the same number of time points (i.e., 150 brain volumes) used across different sites.
Preprocessing of Monkey and Human fMRI Data
Functional brain images of monkey and human were preprocessed using the same steps, including slice timing correction, motion correction, coregistration with individual T1-weighted image, normalization to the corresponding standard space, resampling and spatial smoothing, regression of nuisance signals, removal of linear drift, and temporal filtering.
Specifically, the preprocessing of the monkey data were done using the SPM 8.0 toolbox2 and the FMRIB Software Library toolbox (FSL3). The first 10 volumes were discarded. The field map images of each participant were then applied to compensate for the geometric distortion of EPI images caused by magnetic field inhomogeneity using FSL FUGUE. After slice timing correction and motion correction, the corrected images were normalized to standard space of the monkey F99 atlas4 using an optimum 12-parameter affine transformation and non-linear deformations, and then resampled to 2-mm cubic voxels and spatially smoothed with a 4 mm full-width at half-maximum isotropic Gaussian kernel. Six head motion parameters, ventricle, and white matter signals were removed from the smoothed volumes using linear regression. Linear drift of the volumes was removed and temporal filtering (0.0025–0.05 Hz) (Vincent et al., 2007; ; ) was performed.
The preprocessing of human data was performed by the Preprocessed Connectomes Project (PCP5) using the Data Processing Assistant for Resting-State fMRI (DPARSF) Toolbox (Yan and Zang, 2010). Preprocessing steps included slice timing correction, motion correction, spatial normalization into MNI space, resampled to 3 mm × 3 mm × 3 mm voxels and smoothing with a Gaussian kernel (full-width at half-maximum = 6 mm). Friston-24 parameters of head motion, white matter and ventricle signals were regressed out, followed by linear drift correction and temporal filtering (0.01–0.1 Hz). For more details, readers may refer to the description in the PCP6.
Static and Dynamic FC Analysis
For both preprocessed monkey and human datasets, we first divided the brain into different areas. For a direct cross-species comparison, here we adopted a regional map (RM) parcellation (82 cortical regions) for both monkeys and humans (Table 1), which is based on a combination of microstructural, functional, and topographic features (; ; Reid et al., 2016). This parcellation has the same terminology for the monkey and human brains, but the topographic assignments originate from the monkey cortex. Pearson’s correlation coefficients between the mean time courses of any pair of regions over the whole scan were then calculated to represent static FC, resulting in an 82 × 82 connectivity matrix. Finally, Fisher’s Z-transformation was applied to the connectivity matrix so that their distributions could better satisfy normality.
Table 1
| Label | Abbreviation | Full name |
|---|---|---|
| 1/2 | TCpol | Temporal polar cortex |
| 3/4 | Amyg | Amygdala |
| 5/6 | PHC | Parahippocampal cortex |
| 7/8 | TCi | Inferior temporal cortex |
| 9/10 | TCv | Ventral temporal cortex |
| 11/12 | HC | Hippocampal cortex |
| 13/14 | TCc | Central temporal cortex |
| 15/16 | TCs | Superior temporal cortex |
| 17/18 | VACv | Anterior visual cortex, ventral part |
| 19/20 | V1 | Primary visual cortex |
| 21/22 | PFCoi | Orbital inferior prefrontal cortex |
| 23/24 | V2 | Secondary visual cortex |
| 25/26 | PFCom | Orbitomedial prefrontal cortex |
| 27/28 | Ia | Anterior insula |
| 29/30 | Ip | Posterior insula |
| 31/32 | CCs | Subgenual cingulate cortex |
| 33/34 | PMCvl | Ventrolateral premotor cortex |
| 35/36 | CCp | Posterior cingulate cortex |
| 37/38 | CCr | Retrosplenial cingulate cortex |
| 39/40 | G | Gustatory cortex |
| 41/42 | PFCol | Orbitolateral prefrontal cortex |
| 43/44 | A2 | Secondary auditory cortex |
| 45/46 | PFCvl | Ventrolateral prefrontal cortex |
| 47/48 | A1 | Primary auditory cortex |
| 49/50 | VACd | Anterior visual cortex, dorsal part |
| 51/52 | S2 | Secondary somatosensory cortex |
| 53/54 | PFCpol | Prefrontal pole cortex |
| 55/56 | S1 | Primary somatosensory cortex |
| 57/58 | PFCm | Medial prefrontal cortex |
| 59/60 | PCm | Medial parietal cortex |
| 61/62 | M1 | Primary motor cortex |
| 63/64 | FEF | Frontal eye field |
| 65/66 | CCa | Anterior cingulate cortex |
| 67/68 | PFCcl | Centrolateral prefrontal cortex |
| 69/70 | PCip | Intraparietal cortex |
| 71/72 | PCi | Inferior parietal cortex |
| 73/74 | PCs | Superior parietal cortex |
| 75/76 | PFCdm | Dorsomedial prefrontal cortex |
| 77/78 | PFCdl | Dorsolateral prefrontal cortex |
| 79/80 | PMCdl | Dorsolateral premotor cortex |
| 81/82 | PMCm | Medial premotor cortex |
Brain parcellation of regional map.
Odd numbers denote regions in the left hemisphere, and even numbers denote regions in right hemisphere.
To calculate dynamic FC, we applied a commonly used sliding window approach following our previous study (Yin et al., 2016). Briefly, a tapered window was selected and slid 1 TR, resulting in 168 windows for monkeys and 123 windows for humans. For each time window, Pearson’s correlation coefficients between the mean time courses of any pair of regions were calculated and then a symmetric 82 × 82 connectivity matrix was generated. Thus, dynamic FC matrices were obtained for each participant.
Mapping Intrinsic Functional Flexibility of Brain
Based on the dynamic FC matrices of each participant, we computed the normalized probability distribution Pi (j…n) for a given brain region i as follows:
where n(cij) denotes how many times the connection between i and j emerged across temporal windows, k is a predefined threshold indicating number of the strongest connections reserved for region i at each time window, and w denotes the number of temporal windows. Pi(j) denotes the probability of occurrence for the connection between regions i and j across all temporal windows. The greater the value of Pi(j), the more frequent the interaction between region i and j across the temporal windows, and vice versa.
Regarding the threshold k, we have justified the choice of k for the human dataset in our previous study (Yin et al., 2016) as follows. Taking into account that brain is organized as a sparse and economical functional network, we first considered a wide range of k from 1 to 10. We then calculated the entropy (see definition below) for all brain regions for each k. Subsequently, two parameters were calculated, including contrast (identifying the value of k most sensitive to differences in entropy across the whole brain) and consistency (identify the value of k where the resulting entropy distribution is most representative of the distributions at other thresholds). We finally summed the two metrics, contrast and consistency, at each threshold k, and the peak value of this total was considered as corresponding to the optimal threshold (a peak value emerges at k = 3 for human dataset). Considering that the optimal threshold k may be different for different states or species, we conducted the above analysis for the anesthetized monkeys in this study. We found that the peak value also emerges at k = 3, although the maximum value was at k = 1 for the anesthetized monkeys (Supplementary Figure S1). We therefore used the same threshold k = 3 for both humans and monkeys.
Subsequently, Shannon entropy Ei was applied to the probability distribution of each brain region i:
Here, Ei was used to quantify functional flexibility, which characterizes heterogeneous connectivity between region i and others over time. A higher value of Ei indicates greater functional flexibility, and vice versa. Readers can see our previous study for details regarding methodology (Yin et al., 2016).
Mapping Intrinsic Functional Complexity of Brain
The measure of functional complexity based on the distribution of whole-brain correlation values (for static network) without the need for thresholding has been proposed in earlier work (Zhao et al., 2010; Zamora-Lopez et al., 2016). Inspired by this, we calculated complexity Ci for the distribution of correlation values rij (for dynamical FC) of a node i. Here, we chose to define complexity as the difference between the observed distribution p(rij) and the uniform distribution, which is most robust to variations in the number of bins compared with alternatives such as entropy (Zamora-Lopez et al., 2016). The formula is as follows:
where || means the absolute value, Cm = 2∗(m-1)/m, and m indicates number of bins (here using 50 bins). The Ci reflects functional complexity of a node i.
Mapping Intrinsic Temporal Variability of Brain
A recent study (Zhang et al., 2016) used temporal variability analysis to characterize dynamic functional reconfiguration of each brain region, which can be expressed as follows.
where E[ ] denotes mean value, Fi,j indicates FC profile of node i at time window j and w denotes number of temporal windows. This linear measure does not need a threshold for FC values and is indicative of functional flexibility.
Time Series Analysis of Brain Activity
To evaluate fluctuations of brain activity, we performed a time series analysis. For the time series of each brain region, we first calculated the distribution of the BOLD signal values with m bins (here using 30 bins). Then, we used entropy to quantify the randomness of fluctuations of brain activity. The higher the entropy of time series, the more random the fluctuations of brain activity. To further test the statistical significance of the randomness, we finally compared real entropy and entropies of 1000 random time series with the same number of time points and number of bins, and a Z-score was obtained using the following formula.
where Hreal denotes entropy of the observed time series, Hrand denotes entropy of the random time series, and E[] indicates mean value. If the Z-score approaches zero (theoretically Z-score ≤0), the fluctuations of brain activity tend to be random.
Coupling Between Dynamic and Static Measurements of Spontaneous Brain Activity
To explore the potential effects of anesthesia, we calculated correlations between dynamic measurement (nodal entropy E, complexity C, and temporal variability V) and static measurement (nodal strength, i.e., sum of nodal static FC) within anesthetized monkey and awake human brains, respectively.
Validation Analysis
To validate our results, we considered the effects of different data processing on interspecies comparisons. Following previous studies (Vincent et al., 2007; ; ), we used temporal filtering (0.0025–0.05 Hz) for monkeys, which is different from that commonly used for humans (0.01–0.1 Hz). To test the effect of different temporal filtering, we first compared the brain map of entropy E, complexity C, and temporal variability V between different temporal filters in monkeys. We then performed interspecies comparisons using the same temporal filtering (i.e., 0.01–0.1 Hz).
In our main analyses, we used a different number of time points for monkey (n = 190) and human (n = 145) datasets. To test the effect of a different number of time points, we performed the interspecies comparisons using the same number of time points (n = 145).
In the human subject sample, there are 35 subjects with eyes closed and 143 subjects with eyes open. To consider the effect of eye status, we first compared the brain map of entropy E, complexity C, and temporal variability V between the two human subgroups with different eye status. We then performed the interspecies comparisons for both anesthetized monkeys versus human subjects with eyes closed and anesthetized monkeys versus human subjects with eyes open.
Results
Static FC Patterns in Anesthetized Monkey and Awake Human Brains
We calculated static, time-averaged FC for both human and monkey datasets. We found that the FC in anesthetized monkeys (mean ±SD = 0.25 ± 0.16) was generally weaker than in awake humans (mean ±SD = 0.44 ± 0.13) (effect size: Cohen’s d = -1.3). However, FC between anesthetized monkeys and awake humans was significantly correlated (r = 0.60, p < 0.00001) (Figure 1).
FIGURE 1
Similarity of Intrinsic Functional Flexibility Between Anesthetized Monkey and Awake Human Brains
For the anesthetized monkeys, we found that the brain regions with higher entropy E mainly involved the lateral prefrontal cortex, anterior insula, and medial temporal lobe. The brain regions that showed lower entropy E included primary sensory areas (e.g., auditory and somatosensory regions) and midline default mode regions (e.g., posterior cingulate cortex/retrosplenial cingulate cortex) (Figure 2A).
FIGURE 2
For the awake humans, we observed that brain regions that showed higher entropy E mainly involved the lateral prefrontal, parietal, and temporal cortex, anterior insula, as well as supplementary motor area. The brain regions that showed lower entropy E included primary sensory areas (e.g., auditory, visual, and somatosensory regions) and midline default mode regions (e.g., posterior cingulate cortex/retrosplenial cingulate cortex) (Figure 2B). This result is consistent with our previous study (Yin et al., 2016), despite the use of different brain parcellation and dataset.
Quantitatively, we found a significant correlation of brain-wide entropy E between anesthetized monkeys and awake humans (r = 0.30, p = 0.007) (Figure 3A), although averaged entropy E across the whole brain of anesthetized monkeys (mean ±SD = 5.18 ± 0.18) was higher than that of awake humans (mean ±SD = 4.27 ± 0.41) (effect size: Cohen’s d = 2.9) (Figure 3B). These findings indicate that the heterogeneous flexibility across brain regions is preserved in anesthetized monkeys. However, there are some inconsistencies. For example, the primary visual cortex showed relatively higher entropy E in the anesthetized monkeys, and relatively lower entropy E in the awake humans. In contrast, the inferior parietal cortex exhibited relatively lower entropy E in the anesthetized monkeys, but higher in the awake humans. By comparing brain regions with top 30% and bottom 30% entropy E values in anesthetized monkeys and awake humans, we found that overlapping regions with higher entropy E between species included the left dorsolateral prefrontal cortex, left frontal eye field, left orbitolateral prefrontal cortex, bilateral anterior insula, bilateral orbital inferior prefrontal cortex, bilateral hippocampus, and bilateral parahippocampal cortex (Figure 4A); and that overlapping regions with lower entropy E between species included the bilateral posterior cingulate cortex, bilateral retrosplenial cingulate cortex, left primary auditory cortex, left secondary auditory cortex, right anterior visual area, bilateral secondary somatosensory cortex, bilateral medial parietal cortex, and right posterior insula (Figure 4B).
FIGURE 3
FIGURE 4
Comparison of Functional Complexity Measurement With Our Method
We found that overall complexity C was higher in awake humans compared with anesthetized monkeys, and the correlation of brain-wide complexity C between anesthetized monkeys and awake humans was similar with that obtained using our method (Figure 2C,D, 3C,D). However, few of the well-known flexible cognitive control regions such as lateral prefrontal cortex, anterior insula, and hippocampal cortex showed higher complexity C in anesthetized monkeys (Figure 4C). Instead, we observed unimodal regions such as secondary somatosensory cortex, auditory cortex, and visual cortex showed higher complexity C (Figure 4D). It is possible that this complexity measurement C based on the distribution of correlation values is not suitable for quantifying the heterogeneous functional flexibility of the brain, although it can better assess the functional complexity of the system at different states. For instance, we found the hippocampal cortex exhibited narrower distribution of correlation values (lower complexity) than that of secondary somatosensory cortex in both anesthetized monkeys and awake humans, whereas the time-varying strongest connections of hippocampal cortex was more uniform (higher flexibility) across the whole brain than that of secondary somatosensory cortex (Figure 5).
FIGURE 5
Additionally, we found right inferior parietal cortex showed higher complexity C and primary visual cortex showed lower complexity C in both anesthetized monkeys and awake humans (Figure 4C,D). However, we observed that the strongest connections of primary visual cortex were local and stereotyped in awake humans while distributed and variable across brain in anesthetized monkeys. Moreover, we found inferior parietal cortex more frequently connected with local brain regions in parietal, temporal, and visual cortices in monkeys while more frequently connected with broad brain regions in parietal, temporal, and frontal cortices in humans (Figure 6).
FIGURE 6
Comparison of Time Variability Measurement With Our Method
We found the results obtained by temporal variability measurement V were highly consistent with that using our method, but not complexity measurement C based on distribution of correlation values (Figure 2E,F, 3E,F, 4E,F). This suggests that information of spatial connectivity patterns is more important than distribution of correlation values for describing functional flexibility, even with k (= 3) strongest connections at each time window.
Time Series Analysis of Brain Activity
We found the entropy H of time series to be globally higher in anesthetized monkeys compared with awake humans. In particular, we observed that the primary visual cortex showed the highest entropy H in anesthetized monkeys. However, entropy H of time series for all brain regions were remarkably lower than that of random time series in both anesthetized monkeys and awake humans (Figure 7). This result suggests that brain activity in anesthetized monkeys generally shifts toward the random fluctuations, but still differs from random fluctuations. In addition, the effect of anesthesia on fluctuations of brain activity is probably non-uniform across brain.
FIGURE 7
Distinct Relationships Between Dynamic and Static Measurements Within Awake Human and Anesthetized Monkey Brains
To explore the impacts of anesthesia, we calculated correlations between dynamic measurement and static measurement in awake humans and anesthetized monkeys, respectively. We found a negative correlation (r = -0.47, p = 0.00001) between nodal entropy E and strength in anesthetized monkeys, but not in awake humans (r = -0.19, p = 0.085). In contrast, we found a positive correlation (r = 0.69, p < 0.00001) between nodal complexity C and strength in anesthetized monkeys, but not in awake humans (r = -0.14, p = 0.22). Consistent with our method, we found a negative correlation (r = -0.68, p < 0.00001) between nodal temporal variability V and strength in anesthetized monkeys, but not in awake humans (r = -0.19, p = 0.089) (Figure 8).
FIGURE 8
Validation Analysis
Regarding the effect of temporal filtering, we found significant correlations between two different temporal filters in monkeys for entropy E (r = 0.97, p < 0.00001), complexity C (r = 0.73, p < 0.00001), and variability V (r = 0.96, p < 0.00001) (Supplementary Figure S2). Moreover, we consistently observed significant correlations between species for entropy E (r = 0.26, p = 0.02), complexity C (r = 0.30, p = 0.007), and variability V (r = 0.37, p = 0.0007) with the same time filtering (Supplementary Figure S3). However, the interspecies correlations based on the same temporal filtering were a little weaker than when using different temporal filtering for entropy E and variability V, but the same for complexity C. This result indicates that different temporal filtering for monkeys and humans may give better correspondence between species.
For the effect of different number of time points, we consistently observed significant correlations between species for entropy E (r = 0.28, p = 0.01), complexity C (r = 0.30, p = 0.005), and variability V (r = 0.44, p = 0.00003) with the same number of time points (Supplementary Figure S4). This result suggests that the different number of time points used for monkeys and humans did not significantly affect correlations between species.
For the effect of eye status, we found significant correlations between the two human subgroups for entropy E (r = 0.96, p < 0.00001), complexity C (r = 0.75, p < 0.00001), and variability V (r = 0.91, p < 0.00001) with different eye status (Supplementary Figure S5). For the interspecies comparisons, we found the correlations were similar for the two conditions of human subjects with different eye status (Supplementary Figure S6). This result suggests that our main findings are not significantly affected by eye status of human subjects during resting-state fMRI scanning. Specifically, we observed that the visual cortex showed lowest entropy E in the human subjects for both closed and open eyes. This further implies that eye status does not change the rank of functional flexibility of visual cortex in the brain. In the anesthetized monkeys, the primary visual cortex showed relatively high entropy E. This is probably because anesthesia causes brain activity in primary visual cortex to shift much more toward random fluctuations.
Discussion
Although resting-state connectivity networks have been used for decades to probe functional brain organization (; Raichle, 2015; ), the origins and functional significance of resting-state connectivity patterns require further understanding. In a recent study, we were able to characterize the heterogeneous architecture of intrinsic functional flexibility in the awake, resting human brain using dynamic FC analysis and a probabilistic model (Yin et al., 2016). To further understand functional organization of spontaneous brain activity, in this study, we performed an interspecies comparison of intrinsic functional flexibility between anesthetized monkeys and awake humans.
For reference, we compared conventional, static FC patterns between anesthetized monkeys and awake humans. We found that functional coupling between brain regions in anesthetized monkeys was generally weaker than in awake humans; however, brain-wide connectivity patterns were correlated between species. Previous studies showed that patterns of static FC, such as the default mode network, persist even after loss of consciousness in rodents (; ) and primates (Vincent et al., 2007; ; ). In contrast, other studies suggested a breakdown of both within- and between-network resting-state connectivity in the anesthetized state (; Stamatakis et al., 2010) and deep sleep state (Tagliazucchi et al., 2013). These studies raised two hypotheses about the origins of resting-state FC patterns: reflecting a continuous stream of ongoing cognitive process and random fluctuations constrained by a stable anatomical skeleton (Vincent et al., 2007; ; ; Raichle, 2015). It is possible that persisting static FC patterns during anesthesia can be attributed to anatomical constraints (Vincent et al., 2007; ). On the other hand, consciousness is probably indexed by global integration with strong couplings between long-range brain regions ().
For our dynamic analysis, we found that the brain regions that showed higher entropy E in both anesthetized monkeys and awake humans mainly involved the higher-order association cortex such as the lateral prefrontal cortex, and regions that showed lower entropy E included primary sensory areas and midline default mode regions. This result is consistent with the flexible hub theory: the FC patterns of frontoparietal regions shift more than those of other regions across a variety of task states (). A previous study focusing on the temporal dynamics of resting-state functional networks suggested that wakefulness is characterized by the dynamical exploration of a richer repertoire of functional configurations or states (). Moreover, individuals with brain networks showing greater dynamics perform more favorably in behavioral tasks (). Through comparing anesthetized monkeys with awake humans, Hutchison and his colleagues demonstrated that the temporal dynamics of resting-state FC are also an intrinsic property of brain organization and not simply a consequence of conscious or cognitive processing (). Moreover, there is accumulating evidence that many cognitive processes can occur in the absence of awareness (MacDonald et al., 2015). We have demonstrated that brain regions showing higher entropy E may represent a more flexible exploration of functional configurations, even under anesthesia. Expanding upon previous studies (; Yin et al., 2016; Zhang et al., 2016), this study reveals that heterogeneous functional flexibility across the cortex is evolutionarily conserved and persists across brain states.
Although there is significant correlation of brain-wide entropy E between anesthetized monkeys and awake humans, the whole-brain average entropy E is remarkably higher in the anesthetized monkeys. Barttfeld et al. study indicated that anesthesia may lead to a stable brain state that is more similar to the structure, in which time series of brain activity resemble random fluctuations shaped by fixed anatomical connectivity (). Using dynamical systems modeling, a previous study further suggested that low coupling strength between brain regions can coexist with a single stable spontaneous connectivity pattern (). Because it becomes the only available attractor, the sedated brain cannot depart from it and remains confined to a semirandom exploration of the valley surrounding it, thus simultaneously exhibiting interregional correlations along with fixed anatomical connectivity and a memoryless trajectory (; ). Consistently, through time series analysis, we found that brain activity in anesthetized monkeys generally shifted toward random fluctuations, but it was still different from random fluctuations. We speculate that the general increase in entropy E and decrease in static FC strength in the anesthetized monkey brain is likely attributable to anesthesia-induced random fluctuations of brain activity.
Notably, previous electrophysiological studies have suggested that brain neuronal activity is dominant with slow oscillations under anesthetic states as well as during deep sleep state (; ; Ni Mhuircheartaigh et al., 2013). A remarkable feature of slow oscillation is the synchrony over large cortical areas (; ). have demonstrated that stable functional architecture and dynamic neural activity are concurrent during induction of anesthesia. It is possible that the stable functional architecture, i.e., large-scale functional networks frequently observed even in anesthesia (Vincent et al., 2007; ), is result from synchrony induced by slow oscillation. Regarding dynamics of brain activity, a breakdown of long-range temporal correlations was observed in BOLD signals during both anesthesia and deep sleep states, suggesting that the dynamics of time series is close to white noise (Tagliazucchi et al., 2013; ). In despite of different temporal scales, the findings that BOLD signal appears to be more random can provide a supplement for understanding anesthesia-induced slow oscillations of electrophysiological activity.
Moreover, previous evidence suggests that the dynamic complexity of the system under anesthetized state is reduced (). Accordingly, we found an overall reduction of functional complexity in anesthetized monkeys based on a complexity measure C for distribution of correlation values, but this was not the case for our method and time variability measurement. This suggests that the complexity measure C for distribution of correlation values is better to quantify dynamic complexity of system. In contrast, our method and time variability measurement are more suitable for describing functional flexibility, because they enable the capture of information regarding dynamical spatial connectivity patterns of a node, while this is not the case for the complexity measurement based on distribution of correlation values. For instance, the hippocampal cortex exhibited narrower distribution of correlation values (lower complexity) than the secondary somatosensory cortex in both anesthetized monkeys and awake humans, whereas the time-varying strongest connections of hippocampal cortex were more uniform (higher flexibility) than that of secondary somatosensory cortex across the whole brain. Our findings suggest that combining different methods could provide more complete information for in-depth understanding of functional brain organization.
We further found a negative correlation between nodal entropy E and strength in anesthetized monkeys, but not in awake humans. Previous simulating and empirical data suggest that brain dynamics may change from a single stable state to multi-stable state, as coupling strength between brain areas increases from the sedated to the conscious condition (; ; , ; ). In agreement, conscious processing is supported by global integration with strong coupling between long-distance brain regions as well as a diversity of cognitive states (; ). Loss of consciousness due to anesthesia may lack both strong coupling and a rich repertoire of cognitive states (; ). A possible explanation is that the negative correlation between nodal entropy E and strength (i.e., the weaker connectivity strength, the higher entropy E) is primarily dominated by random fluctuations of brain activity induced by anesthesia. In contrast, the wakeful condition or conscious access with dominance of heterogeneous cognitive states may lead to decoupling between nodal entropy E and strength.
Consistently, we observed no correlation between the complexity measurement C for distribution of correlation values and static connectivity strength for the awake humans, but there was a positive correlation for the anesthetized monkeys. Zamora-Lopez et al. (2016) reported a reverse U-shaped relationship between functional complexity and coupling strength, with an optimal functional organization at the peak complexity. In other words, the functional complexity increases and then decreases during the increase of coupling strength. We therefore speculate that positive correlation between the complexity measurement C and static connectivity strength in anesthetized monkeys is attributed to the anesthesia-induced lower coupling strength (at the left part of reverse U-shape). This result further suggests that the coupling between dynamic and static measures may serve as a potential signature of anesthesia, whereas the direction of correlation is probably dependent on specific metrics.
Although brain-wide correlation exists between species, it should be noted that there were some divergences in intrinsic functional flexibility between the anesthetized monkey and awake human brains. For instance, the primary visual cortex exhibited relatively high entropy E in the anesthetized monkeys, although it was low in the awake humans. In contrast, the inferior parietal cortex showed relatively low entropy E in the anesthetized monkeys, whereas it was high in the awake humans. A previous human study indicated that anesthesia preferentially modulates higher-order connections, but not low-level sensory connections (Martuzzi et al., 2010). One rodent study also showed that anesthesia profoundly impacted the dynamic resting-state FC of neural circuits subserving higher-order functions but had less effect on sensory systems (). On the other hand, a study by reported that the largest reduction of temporal variance of BOLD signals occurred in the visual cortex and parietal cortex in anesthetized rats. Although conflicting conclusions were drawn in previous studies, the converging evidence suggests a non-uniform impact of anesthesia on brain systems. In anesthetized monkeys, we found brain activity of primary visual cortex and sensorimotor cortex to be much more close to random fluctuations. In addition, the strongest connections of primary visual cortex were local and stereotyped in awake humans, but distributed and variable across the brain in anesthetized monkeys. It is possible the much more random fluctuations may contribute to the difference of functional flexibility observed in primary visual cortex between species.
Regarding inferior parietal cortex, we found it to be more frequently connected with local brain regions in parietal, temporal, and visual cortices in monkeys while more frequently connected with broad brain regions in parietal, temporal, and frontal cortices in humans. From an evolutionary perspective, the inferior parietal cortex in the human brain mainly contains Brodmann areas 39 and 40, but monkeys do not have a comparable area (; Raichle, 2015). We speculate that the difference of functional flexibility observed in inferior parietal cortex of monkeys and humans is likely attributed to evolution. Although evolution may indeed result in functional reorganization of specific brain regions, it is hard to separate the contributions of evolution and anesthesia in the current study. A further study with awake monkeys and anesthetized humans may help clarify this question.
In addition, there are some limitations to this study. First, our analysis is based on dynamic FC and is affected by the general limits of this technique, such as temporal resolution of fMRI (; ). Using simultaneous imaging and electrophysiological recording is helpful for interpretation of dynamic FC (). Second, for a direct comparison of brain network between monkeys and humans, we used a regional map template with the same cortical partitions. Previous studies have suggested evolutionary differences in anatomy between monkeys and humans (; Raichle, 2015). Our findings may be potentially affected by the anatomical differences resulting from primate evolution. Finally, the human data that was collected from multiple centers with different acquisition parameters likely contained non-trivial variability across sites and individuals.
Summary
This study combined dynamical complexity measurements and static connectivity strength measurement to understand functional brain organization in anesthetized monkeys and awake humans. Cross-species comparison suggests that the heterogeneous brain map of intrinsic functional flexibility persists during primate evolution and transcends levels of consciousness, remaining present under anesthesia. Moreover, the coupling between dynamic and static measurements can provide a potential signature of loss of consciousness due to anesthesia. However, each method may capture different biological information and have its own limitation. Specifically, our method and temporal variability approach might be more suitable for describing functional flexibility of a node, whereas there is a potential flaw in characterizing the changes of system complexity induced by anesthesia. In contrast, the complexity measurement based on distribution of correlation values is better for evaluating functional complexity of systems with different states or coupling strength, but is probably not suitable for describing functional flexibility of a node due to failure in capturing information of spatial connectivity patterns. Combining different methods could provide more complete information for in-depth understanding of functional brain organization. This study not only offers fresh insight into evolution of functional brain organization, but also advances our understanding of dynamics of spontaneous brain activity.
Statements
Data availability statement
The datasets generated for this study are available on request to the corresponding author.
Author contributions
DY, YW, and ZheW designed the research. DY, ZZ, and ZhiW performed the research. DY, ZhiW, QL, and DC analyzed the data. DY, KZ, and ZheW wrote the paper.
Funding
This work was supported by the National Key Research and Development Program of China (2017YFC1310400 to DY), Youth Innovation Promotion Association of the Chinese Academy of Sciences (to DY), Strategic Priority Research Program of the Chinese Academy of Sciences (XDB32000000 to ZheW), National Natural Science Foundation of China (81571300 and 31771174 to ZheW, 31600869 and 81527901 to DY), and Natural Science Foundation and Major Basic Research Program of Shanghai (16JC1420100 to ZheW).
Acknowledgments
We thank Qinying Jiang and Wenwen Yu for their assistance in monkey data collection and also thank Jinqiang Peng for his assistance in data management. DY gratefully acknowledges the support of SA-SIBS Scholarship Program.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2019.00174/full#supplementary-material
Footnotes
1.^http://fcon_1000.projects.nitrc.org/indi/abide/
2.^http://www.fil.ion.ucl.ac.uk/spm
4.^http://sumsdb.wustl.edu/sums/macaquemore.do
5.^http://preprocessed-connectomes-project.org/abide/index.html
6.^http://preprocessed-connectomes-project.org/abide/dparsf.html
References
1
AchermannP.BorbelyA. A. (1997). Low-frequency (< 1 Hz) oscillations in the human sleep electroencephalogram.Neuroscience81213–222. 10.1016/S0306-4522(97)00186-3
2
AlkireM. T.HudetzA. G.TononiG. (2008). Consciousness and anesthesia.Science322876–880. 10.1126/science.1149213
3
AllenE. A.DamarajuE.PlisS. M.ErhardtE. B.EicheleT.CalhounV. D. (2014). Tracking whole-brain connectivity dynamics in the resting state.Cereb. Cortex24663–676. 10.1093/cercor/bhs352
4
BarttfeldP.UhrigL.SittJ. D.SigmanM.JarrayaB.DehaeneS. (2015). Signature of consciousness in the dynamics of resting-state brain activity.Proc. Natl. Acad. Sci. U.S.A.112887–892. 10.1073/pnas.1418031112
5
BassettD. S.SpornsO. (2017). Network neuroscience.Nat. Neurosci.20353–364. 10.1038/nn.4502
6
BezginG.VakorinV. A.van OpstalA. J.McIntoshA. R.BakkerR. (2012). Hundreds of brain maps in one atlas: registering coordinate-independent primate neuro-anatomical data to a standard brain.Neuroimage6267–76. 10.1016/j.neuroimage.2012.04.013
7
BiswalB.YetkinF. Z.HaughtonV. M.HydeJ. S. (1995). Functional connectivity in the motor cortex of resting human brain using echo-planar MRI.Magn. Reson. Med.34537–541. 10.1002/mrm.1910340409
8
BoverouxP.VanhaudenhuyseA.BrunoM. A.NoirhommeQ.LauwickS.LuxenA.et al (2010). Breakdown of within- and between-network resting state functional magnetic resonance imaging connectivity during propofol-induced loss of consciousness.Anesthesiology1131038–1053. 10.1097/ALN.0b013e3181f697f5
9
BraunU.SchaferA.WalterH.ErkS.Romanczuk-SeiferthN.HaddadL.et al (2015). Dynamic reconfiguration of frontal brain networks during executive cognition in humans.Proc. Natl. Acad. Sci. U.S.A.11211678–11683. 10.1073/pnas.1422487112
10
BreshearsJ. D.RolandJ. L.SharmaM.GaonaC. M.FreudenburgZ. V.TempelhoffR.et al (2010). Stable and dynamic cortical electrophysiology of induction and emergence with propofol anesthesia.Proc. Natl. Acad. Sci. U.S.A.10721170–21175. 10.1073/pnas.1011949107
11
CalhounV. D.MillerR.PearlsonG.AdaliT. (2014). The chronnectome: time-varying connectivity networks as the next frontier in fMRI data discovery.Neuron84262–274. 10.1016/j.neuron.2014.10.015
12
ChangC.GloverG. H. (2010). Time-frequency dynamics of resting-state brain connectivity measured with fMRI.Neuroimage5081–98. 10.1016/j.neuroimage.2009.12.011
13
ChenT.CaiW.RyaliS.SupekarK.MenonV. (2016). Distinct global brain dynamics and spatiotemporal organization of the salience network.PLoS Biol.14:e1002469. 10.1371/journal.pbio.1002469
14
ChristoffK.IrvingZ. C.FoxK. C.SprengR. N.Andrews-HannaJ. R. (2016). Mind-wandering as spontaneous thought: a dynamic framework.Nat. Rev. Neurosci.17718–731. 10.1038/nrn.2016.113
15
ColeM. W.ReynoldsJ. R.PowerJ. D.RepovsG.AnticevicA.BraverT. S. (2013). Multi-task connectivity reveals flexible hubs for adaptive task control.Nat. Neurosci.161348–1355. 10.1038/nn.3470
16
de PasqualeF.Della PennaS.SpornsO.RomaniG. L.CorbettaM. (2015). A Dynamic core network and global efficiency in the resting human brain.Cereb. Cortex264015–4033. 10.1093/cercor/bhv185
17
DecoG.JirsaV. K.McIntoshA. R. (2011). Emerging concepts for the dynamical organization of resting-state activity in the brain.Nat. Rev. Neurosci.1243–56. 10.1038/nrn2961
18
DecoG.JirsaV. K.McIntoshA. R. (2013). Resting brains never rest: computational insights into potential cognitive architectures.Trends Neurosci.36268–274. 10.1016/j.tins.2013.03.001
19
DehaeneS.ChangeuxJ. P. (2005). Ongoing spontaneous activity controls access to consciousness: a neuronal model for inattentional blindness.PLoS Biol.3:e141. 10.1371/journal.pbio.0030141
20
DehaeneS.ChangeuxJ. P. (2011). Experimental and theoretical approaches to conscious processing.Neuron70200–227. 10.1016/j.neuron.2011.03.018
21
DestexheA.ContrerasD.SteriadeM. (1999). Spatiotemporal analysis of local field potentials and unit discharges in cat cerebral cortex during natural wake and sleep states.J. Neurosci.194595–4608. 10.1523/JNEUROSCI.19-11-04595.1999
22
Di MartinoA.YanC. G.LiQ.DenioE.CastellanosF. X.AlaertsK.et al (2014). The autism brain imaging data exchange: towards a large-scale evaluation of the intrinsic brain architecture in autism.Mol. Psychiatry19659–667. 10.1038/mp.2013.78
23
FoxM. D.RaichleM. E. (2007). Spontaneous fluctuations in brain activity observed with functional magnetic resonance imaging.Nat. Rev. Neurosci.8700–711. 10.1038/nrn2201
24
GhoshA.RhoY.McIntoshA. R.KotterR.JirsaV. K. (2008). Noise during rest enables the exploration of the brain’s dynamic repertoire.PLoS Comput. Biol.4:e1000196. 10.1371/journal.pcbi.1000196
25
HindriksR.AdhikariM. H.MurayamaY.GanzettiM.MantiniD.LogothetisN. K.et al (2016). Can sliding-window correlations reveal dynamic functional connectivity in resting-state fMRI?Neuroimage127242–256. 10.1016/j.neuroimage.2015.11.055
26
HoneyC. J.SpornsO.CammounL.GigandetX.ThiranJ. P.MeuliR.et al (2009). Predicting human resting-state functional connectivity from structural connectivity.Proc. Natl. Acad. Sci. U.S.A.1062035–2040. 10.1073/pnas.0811168106
27
HudetzA. G.LiuX.PillayS. (2015). Dynamic repertoire of intrinsic brain states is reduced in propofol-induced unconsciousness.Brain Connect.510–22. 10.1089/brain.2014.0230
28
HutchisonR. M.EverlingS. (2012). Monkey in the middle: why non-human primates are needed to bridge the gap in resting-state investigations.Front. Neuroanat.6:29. 10.3389/fnana.2012.00029
29
HutchisonR. M.LeungL. S.MirsattariS. M.GatiJ. S.MenonR. S.EverlingS. (2011). Resting-state networks in the macaque at 7 T.Neuroimage561546–1555. 10.1016/j.neuroimage.2011.02.063
30
HutchisonR. M.WomelsdorfT.AllenE. A.BandettiniP. A.CalhounV. D.CorbettaM.et al (2013a). Dynamic functional connectivity: promise, issues, and interpretations.Neuroimage80360–378. 10.1016/j.neuroimage.2013.05.079
31
HutchisonR. M.WomelsdorfT.GatiJ. S.EverlingS.MenonR. S. (2013b). Resting-state networks show dynamic functional connectivity in awake humans and anesthetized macaques.Hum. Brain Mapp.342154–2177. 10.1002/hbm.22058
32
IsomuraY.SirotaA.OzenS.MontgomeryS.MizusekiK.HenzeD. A.et al (2006). Integration and segregation of activity in entorhinal-hippocampal subregions by neocortical slow oscillations.Neuron52871–882. 10.1016/j.neuron.2006.10.023
33
JiaH.HuX.DeshpandeG. (2014). Behavioral relevance of the dynamics of the functional brain connectome.Brain Connect.4741–759. 10.1089/brain.2014.0300
34
KarahanogluF. I.Van De VilleD. (2015). Transient brain activity disentangles fMRI resting-state dynamics in terms of spatially and temporally overlapping networks.Nat. Commun.6:7751. 10.1038/ncomms8751
35
KeilholzS.Caballero-GaudesC.BandettiniP.DecoG.CalhounV. (2017). Time-resolved resting-state functional magnetic resonance imaging analysis: current status, challenges, and new directions.Brain Connect.7465–481. 10.1089/brain.2017.0543
36
KeilholzS. D. (2014). The neural basis of time-varying resting-state functional connectivity.Brain Connect.4769–779. 10.1089/brain.2014.0250
37
KopellN. J.GrittonH. J.WhittingtonM. A.KramerM. A. (2014). Beyond the connectome: the dynome.Neuron831319–1328. 10.1016/j.neuron.2014.08.016
38
KotterR.WankeE. (2005). Mapping brains without coordinates.Philos. Trans. R. Soc. Lond. B Biol. Sci.360751–766. 10.1098/rstb.2005.1625
39
LiangZ.KingJ.ZhangN. (2012). Intrinsic organization of the anesthetized brain.J. Neurosci.3210183–10191. 10.1523/JNEUROSCI.1020-12.2012
40
LiangZ.LiuX.ZhangN. (2015). Dynamic resting state functional connectivity in awake and anesthetized rodents.Neuroimage10489–99. 10.1016/j.neuroimage.2014.10.013
41
LiegeoisR.LaumannT. O.SnyderA. Z.ZhouJ.YeoB. T. T. (2017). Interpreting temporal fluctuations in resting-state functional connectivity MRI.Neuroimage163437–455. 10.1016/j.neuroimage.2017.09.012
42
LiuX.DuynJ. H. (2013). Time-varying functional network information extracted from brief instances of spontaneous brain activity.Proc. Natl. Acad. Sci. U.S.A.1104392–4397. 10.1073/pnas.1216856110
43
LogothetisN. K.GuggenbergerH.PeledS.PaulsJ. (1999). Functional imaging of the monkey brain.Nat. Neurosci.2555–562. 10.1038/9210
44
LuH.ZouQ.GuH.RaichleM. E.SteinE. A.YangY. (2012). Rat brains also have a default mode network.Proc. Natl. Acad. Sci. U.S.A.1093979–3984. 10.1073/pnas.1200506109
45
LvQ.YangL.LiG.WangZ.ShenZ.YuW.et al (2016). Large-scale persistent network reconfiguration induced by ketamine in anesthetized monkeys: relevance to mood disorders.Biol. Psychiatry79765–775. 10.1016/j.biopsych.2015.02.028
46
MacDonaldA. A.NaciL.MacDonaldP. A.OwenA. M. (2015). Anesthesia and neuroimaging: investigating the neural correlates of unconsciousness.Trends Cogn. Sci.19100–107. 10.1016/j.tics.2014.12.005
47
MantiniD.CorbettaM.RomaniG. L.OrbanG. A.VanduffelW. (2013). Evolutionarily novel functional networks in the human brain?J. Neurosci.333259–3275. 10.1523/JNEUROSCI.4392-12.2013
48
MarguliesD. S.VincentJ. L.KellyC.LohmannG.UddinL. Q.BiswalB. B.et al (2009). Precuneus shares intrinsic functional architecture in humans and monkeys.Proc. Natl. Acad. Sci. U.S.A.10620069–20074. 10.1073/pnas.0905314106
49
MartuzziR.RamaniR.QiuM.RajeevanN.ConstableR. T. (2010). Functional connectivity and alterations in baseline brain state in humans.Neuroimage49823–834. 10.1016/j.neuroimage.2009.07.028
50
Miranda-DominguezO.MillsB. D.GraysonD.WoodallA.GrantK. A.KroenkeC. D.et al (2014). Bridging the gap between the human and macaque connectome: a quantitative comparison of global interspecies structure-function relationships and network topology.J. Neurosci.345552–5563. 10.1523/JNEUROSCI.4229-13.2014
51
MuellerS.WangD.FoxM. D.YeoB. T.SepulcreJ.SabuncuM. R.et al (2013). Individual variability in functional connectivity architecture of the human brain.Neuron77586–595. 10.1016/j.neuron.2012.12.028
52
NakaharaK.HayashiT.KonishiS.MiyashitaY. (2002). Functional MRI of macaque monkeys performing a cognitive set-shifting task.Science2951532–1536. 10.1126/science.1067653
53
NeubertF. X.MarsR. B.ThomasA. G.SalletJ.RushworthM. F. (2014). Comparison of human ventral frontal cortex areas for cognitive control and language with areas in monkey frontal cortex.Neuron81700–713. 10.1016/j.neuron.2013.11.012
54
Ni MhuircheartaighR.WarnabyC.RogersR.JbabdiS.TraceyI. (2013). Slow-wave activity saturation and thalamocortical isolation during propofol anesthesia in humans.Sci. Transl. Med.5:208ra148. 10.1126/scitranslmed.3006007
55
PowerJ. D.SchlaggarB. L.PetersenS. E. (2014). Studying brain organization via spontaneous fMRI signal.Neuron84681–696. 10.1016/j.neuron.2014.09.007
56
PretiM. G.BoltonT. A.Van De VilleD. (2017). The dynamic functional connectome: state-of-the-art and perspectives.Neuroimage16041–54. 10.1016/j.neuroimage.2016.12.061
57
RaichleM. E. (2015). The brain’s default mode network.Annu. Rev. Neurosci.38433–447. 10.1146/annurev-neuro-071013-014030
58
ReidA. T.LewisJ.BezginG.KhundrakpamB.EickhoffS. B.McIntoshA. R.et al (2016). A cross-modal, cross-species comparison of connectivity measures in the primate brain.Neuroimage125311–331. 10.1016/j.neuroimage.2015.10.057
59
StamatakisE. A.AdapaR. M.AbsalomA. R.MenonD. K. (2010). Changes in resting neural connectivity during propofol sedation.PLoS One5:e14224. 10.1371/journal.pone.0014224
60
TagliazucchiE.von WegnerF.MorzelewskiA.BrodbeckV.JahnkeK.LaufsH. (2013). Breakdown of long-range temporal dependence in default mode and attention networks during deep sleep.Proc. Natl. Acad. Sci. U.S.A.11015419–15424. 10.1073/pnas.1312848110
61
VanduffelW.FizeD.PeuskensH.DenysK.SunaertS.ToddJ. T.et al (2002). Extracting 3D from motion: differences in human and monkey intraparietal cortex.Science298413–415. 10.1126/science.1073574
62
VincentJ. L.PatelG. H.FoxM. D.SnyderA. Z.BakerJ. T.Van EssenD. C.et al (2007). Intrinsic functional architecture in the anaesthetized monkey brain.Nature44783–86. 10.1038/nature05758
63
WangL.UhrigL.JarrayaB.DehaeneS. (2015). Representation of numerical and sequential patterns in macaque and human brains.Curr. Biol.251966–1974. 10.1016/j.cub.2015.06.035
64
WangZ.ChenL. M.NegyessyL.FriedmanR. M.MishraA.GoreJ. C.et al (2013). The relationship of anatomical and functional connectivity to resting-state connectivity in primate somatosensory cortex.Neuron781116–1126. 10.1016/j.neuron.2013.04.023
65
YanC. G.ZangY. F. (2010). DPARSF: a MATLAB Toolbox for ”Pipeline” Data Analysis of Resting-State fMRI.Front. Syst. Neurosci.4:13. 10.3389/fnsys.2010.00013
66
YinD.LiuW.ZeljicK.WangZ.LvQ.FanM.et al (2016). Dissociable changes of frontal and parietal cortices in inherent functional flexibility across the human life span.J. Neurosci.3610060–10074. 10.1523/JNEUROSCI.1476-16.2016
67
ZaleskyA.FornitoA.CocchiL.GolloL. L.BreakspearM. (2014). Time-resolved resting-state brain networks.Proc. Natl. Acad. Sci. U.S.A.11110341–10346. 10.1073/pnas.1400181111
68
Zamora-LopezG.ChenY.DecoG.KringelbachM. L.ZhouC. (2016). Functional complexity emerging from anatomical constraints in the brain: the significance of network modularity and rich-clubs.Sci. Rep.6:38424. 10.1038/srep38424
69
ZhangD.RaichleM. E. (2010). Disease and the brain’s dark energy.Nat. Rev. Neurol.615–28. 10.1038/nrneurol.2009.198
70
ZhangJ.ChengW.LiuZ.ZhangK.LeiX.YaoY.et al (2016). Neural, electrophysiological and anatomical basis of brain-network variability and its characteristic changes in mental disorders.Brain1392307–2321. 10.1093/brain/aww143
71
ZhangJ.HuangZ.ChenY.ZhangJ.GhindaD.NikolovaY.et al (2018). Breakdown in the temporal and spatial organization of spontaneous brain activity during general anesthesia.Hum. Brain Mapp.392035–2046. 10.1002/hbm.23984
72
ZhaoM.ZhouC.ChenY.HuB.WangB. H. (2010). Complexity versus modularity and heterogeneity in oscillatory networks: combining segregation and integration in neural systems.Phys. Rev. E Stat. Nonlin. Soft Matter Phys.82:046225. 10.1103/PhysRevE.82.046225
Summary
Keywords
dynamic brain organization, intrinsic functional flexibility, evolution, monkey, anesthesia, resting-state fMRI
Citation
Yin D, Zhang Z, Wang Z, Zeljic K, Lv Q, Cai D, Wang Y and Wang Z (2019) Brain Map of Intrinsic Functional Flexibility in Anesthetized Monkeys and Awake Humans. Front. Neurosci. 13:174. doi: 10.3389/fnins.2019.00174
Received
26 October 2018
Accepted
14 February 2019
Published
28 February 2019
Volume
13 - 2019
Edited by
Shella Keilholz, Emory University, United States
Reviewed by
Changsong Zhou, Hong Kong Baptist University, Hong Kong; Scott Peltier, University of Michigan, United States
Updates
Copyright
© 2019 Yin, Zhang, Wang, Zeljic, Lv, Cai, Wang and Wang.
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: Zheng Wang, zheng.wang@ion.ac.cn
This article was submitted to Brain Imaging Methods, a section of the journal Frontiers in Neuroscience
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.