Spontaneous Brain Activity Relates to Autonomic Arousal

Although possible sources and functions of the resting-state networks (RSNs) of the brain have been proposed, most evidence relies on circular logic and reverse inference. We propose that autonomic arousal provides an objective index of psychophysiological states during rest that may also function as a driving source of the activity and connectivity of RSNs. Recording blood oxygenation level-dependent (BOLD) signal using functional magnetic resonance imaging and skin conductance simultaneously during rest in human subjects, we found that the spontaneous fluctuations of BOLD signals in key nodes of RSNs are associated with changes in nonspecific skin conductance response, a sensitive psychophysiological index of autonomic arousal. Our findings provide evidence of an important role for the autonomic nervous system to the spontaneous activity of the brain during “rest.”


Introduction
Despite an explosion of interest in the default mode network (DMN) (Raichle et al., 2001) and the anticorrelated task-positive network (TPN) (Fox et al., 2005), which compose the restingstate networks (RSNs) of the brain during "rest" (Deco et al., 2011), hypotheses about the functional role of spontaneous brain activity nearly always rely on the logically flawed practice of reverse inference (Deco et al., 2011;Poldrack, 2011). Among others, one attempt to explore the causal basis of RSNs has been the simultaneous recording of the brain's electrophysiological activity (Mantini et al., 2007). However, the electrical activity of the brain is known to fundamentally relate to blood oxygenation level-dependent (BOLD) activity (Logothetis et al., 2001). Thus, electrophysiological measures may merely provide a more specific explanation of the neural activity relating to spontaneous BOLD fluctuations (Fox and Raichle, 2007), rather than a causal explanation or objective index of mental states.
By definition, rest suggests no specific instructions. Therefore, we cannot use the common approach of task-induced changes in psychological processes to examine brain activity; activity could no longer be labeled "spontaneous." However, it has been recognized that rest is more akin to a task (with relatively unknown psychological correlates) than to a baseline (Deco et al., 2011). One tenable approach to objectively measure mental states during rest is to use psychophysiological indices, much as has been done for decades in emotion research (Schachter and Singer, 1962). In relation to RSNs, the "signal" is primarily thought to comprise a combination of anatomical and functional neural connectivity patterns (Fox and Raichle, 2007;Mantini et al., 2007;Honey et al., 2009), while the "noise" is related to purportedly confounding sources, such as activity related to the autonomic nervous system (ANS; e.g., heart-rate, respiration) (Chang and Glover, 2009;Iacovella and Hasson, 2011). However, evidence suggests that variation in arousal and other ANS activity can affect RSNs (Birn et al., 2008;Iacovella and Hasson, 2011). A significant relationship between ANS activity and RSNs would provide evidence that psychophysiological states during restingstate data collection are a potential source and/or functional explanation of correlational patterns in spontaneous brain activity.
We studied the contribution of ANS to RSN activity by measuring nonspecific (nontask) skin conductance response (SCR) and related brain activity and connectivity during rest. SCR shares common neural regions with the TPN [e.g., anterior insular (AI) and anterior cingulate cortices (ACC)] for autonomic, affective, and cognitive integration (Critchley, 2002;Critchley et al., 2011). While respiratory and heart rate variability (HRV) can serve as indices of ANS activity at the same filter band for restingstate functional connectivity MRI (rs-fcMRI) (Shmueli et al., 2007;Birn et al., 2008), the typical SCR curve and hemodynamic response function exhibit similar waveforms (Boucsein, 1992), sparing the need for convolution (Patterson et al., 2002) and facilitating modeling. Because the brain is part of a dynamic homeostatic system (Thompson and Varela, 2001;Craig, 2002;Deco et al., 2011), we hypothesized that autonomic arousal would be associated with resting-state functional activity and connectivity of the brain.

Participants
Healthy volunteers (N ϭ 15, male only, mean Ϯ SD age, 27.1 Ϯ 8.2 years) participated in this study. The consent procedure was approved by the institutional review board of Mount Sinai School of Medicine, and written informed consent was obtained from each participant.

Data acquisition
Skin conductance recording. GSR100C (BIOPAC Systems), together with the base module MP150 and the AcqKnowledge software (version 3.9.1.6), was used to record skin conductance during the rs-fcMRI session. The GSR100C applies a constant voltage (0.5 V) between the two electrodes attached to the skin to measure skin conductance. It measures the skin conductance level (SCL) and SCR, which vary with sweat gland activity due to stress, arousal, or emotional excitement. Skin conductance, measured in S, was recorded using a 2000 Hz sampling rate (gain ϭ 2 S/V, both high-pass filters ϭ DC, low-pass filter ϭ 10 Hz). Two EL507 disposable EDA (isotonic gel) electrodes were placed on the palmar surface of the distal phalanges of the big and second toes of the left foot after cleaning with alcohol preps. The signal was low-pass filtered (using the MRI-Compatible MRI CBL/FILTER System MECMRI-TRANS) to reduce radio frequency interference from the scanner. Digital event markers were recorded to enable precise time alignment of skin conductance recording with scan onsets. BIOPAC recording was synchronized to the E-Prime program via the parallel port of the computer.
Image data acquisition. MRI acquisitions were obtained on a 3 T Siemens Allegra MRI system at Mount Sinai School of Medicine. All participants underwent one session with rs-fcMRI first and then fMRI. The scan session lasted ϳ1 h in total length. Foam padding was used to reduce participant head motion. All images were acquired along axial planes parallel to the anterior commissure (AC)-posterior commissure (PC) line. A high-resolution T2-weighted anatomical volume of the whole  brain was acquired on an axial plane parallel to the AC-PC line with a turbo spin-echo pulse sequence with the following parameters: 40 axial slices of 4 mm thickness; skip ϭ 0 mm; repetition time (TR) ϭ 4050 ms; echo time (TE) ϭ 99 ms; flip angle ϭ 170°; field of view (FOV) ϭ 240 mm; matrix size ϭ 448 ϫ 512, voxel size ϭ 0.47 ϫ 0.47 ϫ 4 mm. One run of T2*-weighted images was acquired for rs-fcMRI. Slices were obtained corresponding to the T2-weighted images. The rs-fcMRI was performed using a gradient echoplanar imaging (EPI) sequence with the following parameters: 40 axial slices, 4 mm thick; skip ϭ 0 mm; TR ϭ 2500 ms; TE ϭ 27 ms; flip angle ϭ 82°; FOV ϭ 240 mm; and matrix size ϭ 64 ϫ 64. The rs-fcMRI run started with two dummy volumes before the onset of the fixation to allow for equilibration of T1 saturation effects, followed by 144 image volumes.
Resting-state functional connectivity MRI procedure. For the rs-fcMRI, fixation crosshairs were presented in the center of the screen using E-Prime (Psychology Software Tools) for the duration of the run, which lasted 360 s (6 min). Participants were instructed to maintain gaze on the fixation crosshairs and to minimize movement, with the following in-struction on the screen before the onset of the scan: "This session takes approximately 6 min. You do not need to perform any tasks during this session. Just fix your eyes at the crosshairs. It is very important that you do not move your head or feet during this time or at all throughout the scan."

Data analysis
Skin conductance response preprocessing. The SCR waveform was downsampled by averaging the data points in each 2.5 s bin to match the TR (2.5 s) of the EPI scan of imaging acquisition. Then the SCR waveform (without thresholding) was detrended to remove the typical slow linear decrease of SCL as a function of relaxation time, and bandpass filtered with the same frequency range (0.01-0.12) as in typical rs-fcMRI analysis. Another consideration related to the frequency range selection is to overlap the known HRV, another important index of ANS activity (Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology, 1996). The power spectral density of HRV in the low-frequency (LF) range (0.04 -0.15) and in the high-frequency range (0.15-0.40) corresponds to both sympathetic and parasympathetic activity. The power in the very LF range (Ͻ0.04) warrants further examination. Individual skin conductance time series during the resting-state scan were also examined and compared to the standardized SCR curve (Boucsein, 1992;Patterson et al., 2002;Bach et al., 2010) to identify and count the number of SCRs. Coherence analysis. Coherence analysis between the posterior cingulate cortex (PCC) seed (see seed region definition in functional connectivity analysis) and SCR was also conducted to depict the nature of the correlation between this region of interest (ROI) and SCR. Both PCC time course and SCR were detrended before computing coherence using the mscohere function in Matlab (MathWorks) to estimate the magnitude squared coherence function. FFT length, which determines the frequencies at which the coherence is estimated, was 129; sampling frequency was 1/2.5 Hz; the number of samples to use for each section was 18; and the number of samples by which the sections overlapped was 9 (50%). Following Fisher's r-to-z transformation, the data points of magnitude squared coherence against frequency from each subject were averaged and were then transformed back to r.
Regression analysis. General linear modeling (GLM) of the imaging data was conducted using statistical parametric mapping (SPM8; Wellcome Trust Centre for Neuroimaging, London, UK). The EPI scans were realigned to the first volume, timing corrected, coregistered to the T2 image, normalized to a standard template [Montreal Neurological Institute (MNI)], and spatially smoothed with an 8 ϫ 8 ϫ 8 mm full-width at half-maximum (FWHM) Gaussian kernel. Timing correction of the scan was to the first slice which occurred at 0 s. However, the SCR was an average of the samples collected from time point 0 to 2.5 s. GLM (Friston et al., 1994) was then conducted with the SCR time series as a predictor of the observed BOLD signals. The regressor was not generated as in standard fMRI data analysis by convolving default SPM basis function with delta functions, because there was no sequence of individual events. Lowfrequency drifts in signal were removed using a standard high-pass filter with a 128 s cutoff. Serial correlation was estimated using an autoregressive AR(1) model. Mean voxel value was used for global calculation and grand mean scaling was applied with global normalization to remove nonspecific noise (Van Dijk et al., 2010). Ventricle and white-matter signals were extracted using corresponding mask and entered as covariates. In addition, the six parameters generated during motion correction were also entered as covariates.
Subsequently, the relationship between SCR and BOLD was tested by employing a mixed effect model (Friston et al., 2005) implemented in SPM12 ␣ (Wellcome Trust Centre for Neuroimaging). The mixed effect model was used because unlike event-related fMRI, the number of the nonspecific SCRs could not be controlled; nonspecific SCRs varied considerably across subjects (mean ϭ 5.8, SD ϭ 5.7). The fixed-effect model was specified and estimated first using the data from all subjects. Then the mixed effect model was specified and estimated. The resultant voxelwise statistical maps were thresholded for significance using a cluster-size algorithm that protects against an inflation of the false-positive rate of multiple comparisons. An uncorrected p value of 0.05 for the height (intensity) threshold of each activated voxel and an extent threshold of k ϭ 120 were used. A Monte Carlo simulation of the brain volume of the current study was conducted to establish an appropriate voxel contiguity threshold (Slotnick and Schacter, 2004). Assuming an individual voxel type I error of p Ͻ 0.05, a cluster extent of 120 contiguous resampled voxels (2 ϫ 2 ϫ 2 mm 3 ) was indicated as necessary to correct for multiple voxel comparisons at p Ͻ 0.05. The same cluster-level threshold was applied to all contrasts. Statistical results were mapped onto the surface of the cerebral cortex.
Functional connectivity analysis. Functional connectivity, operationally defined as simple correlations between activation of brain areas (Friston et al., 1997), was computed using the traditional correlation analysis of the time courses among brain regions (Koshino et al., 2005). For each subject, time series volumes of rs-fcMRI scan images were preprocessed using the Data Processing Assistant for Resting-State fMRI (DPARSF) toolbox (Chao-Gan and Yu-Feng, 2010). The analysis included slice timing correction, realignment (to the first volume), coregistration, normalization (to the MNI space with unified segmentation on T2 images), spatial smoothing (using a 8 mm FWHM Gaussian kernel), and detrending (to remove the systematic drift or trend) and temporal filtering [bandpass, 0.01ϳ0.12, based on coherence analysis, which is wider than 0.01ϳ0.08 (Biswal et al., 1995), to reduce the effect of low frequency drift and high frequency physiological signal or noise]. To examine the connectivity of the default mode network, the time course of the PCC (left and right combined) was extracted using the automated anatomical labeling template (Tzourio-Mazoyer et al., 2002). Then the resting-state fMRI data analysis toolkit (REST) (Song et al., 2011) was used to calculate the voxelwise linear correlation between the mean time course of the PCC and the time course of each voxel in the whole brain, with six head motion parameters, global mean signal, white-matter signal, and CSF signal as covariates. To explore the effect of SCR, the second voxelwise connectivity analysis was performed by adding the SCR as a covariate to regress out the SCR effect. The correlation coefficients were transformed using Fisher's r-to-z transformations. A paired t test was conducted be- . Connectivity between PCC and other brain regions without regressed out SCR-related brain activity. The green color indicates voxels with positive connectivity to PCC, while the red color indicates voxels with anticorrelated connectivity. p Ͻ 0.001 and k Ͼ 120 resampled voxel size thresholds were used to plot this figure. tween the functional connectivity of PCC before and after regressing out the SCR effect. Psychophysiological interaction analysis. Psychophysiological interaction (PPI) analysis was conducted, setting PCC as the seed region, to examine whether fluctuations in arousal (the psychological factor, reflected by SCR) modulate connectivity between PCC and other key regions [such as the ventromedial prefrontal cortex (vmPFC)] in the RSNs. The raw spontaneous fluctuation (SF) wave of skin conductance was Figure 4. The connectivity changes in default mode and other brain networks before and after regressing SCR out. The red color indicates voxels with decreased anticorrelated connectivity, while the green color indicates voxels with decreased positive connectivity with PCC after regressing out SCR.  (Bach et al., 2009(Bach et al., , 2010(Bach et al., , 2011 was used to deconvolve the SCR signal, generating the onset and amplitude vectors of SCR. Because a fixed number of SCRs is assumed and only their amplitudes are estimated, the amplitude of "unnecessary" SFs is estimated to be close to zero. Events where the amplitude was below a threshold (in unit of S, varied across subjects) were excluded. The onset and amplitude vectors were used in a second GLM using the amplitude vectors as the modulator of the onset stick of SCR, convolved with the hemodynamic response function, but not SCR function, prepared for the PPI analysis. All the covariates used in the first GLM were also used here for the second GLM and the following PPI analysis.
From the first-level analysis, the BOLD signal of the PCC ROI was extracted by using the PCC mask (as in previous functional connectivity analysis), and the PPI variables were created using the F contrast of the SCR regressor. SF amplitude was not included because it does not depend on tonic arousal, while the number of SFs does. GLM was conducted with PPI regressors of (1) interaction of PCC by SCR, (2) main effect of PCC, and (3) main effect of SCR, corresponding to PPI.ppi, PPI.Y, and PPI.P in the design matrix. Here PPI.P (SCR) was treated as the psychological factor. The second-level group data analysis was then conducted using the mixed-effect model. The threshold was the same as in the previous GLM analysis (individual voxel p Ͻ 0.05 and extent k Ͼ 120 voxels). The significantly activated regions may indicate that (1) the contribution of PCC to those regions is altered by the experimental (psychological) context (here it is the SCR that is related to arousal), or (2) the response of those regions to SCR due to the contribution of PCC.

Coherence between the posterior cingulate cortex and SCR
Because PCC is the key node of the DMN (Greicius et al., 2003), we first examined the time course of the PCC ROI and SCR ( Fig.  1 A, B, D, E, from a single subject). The averaged power spectral density showed that SCR and PCC time courses share the same frequency range (Fig. 1C,F, averaged across all subjects). Coherence analysis, to identify variations of the two signals with similar spectral properties, showed that the activity of the PCC ROI and SCR covary (Fig. 1G, averaged across all subjects). This indicates a potentially important relationship between spontaneous BOLD fluctuations of PCC and SCR activity. The pattern of coherence between BOLD signal of PCC and SCR may reflect the underlying control of the ANS by the DMN and/or TPN, or brain response to ANS activity.

Brain activity as a function of nonspecific SCR
GLM showed that deactivation of the PCC (and precuneus) of DMN and activation of the ACC and AI of the TPN, were significantly correlated with nonspecific SCR (Fig. 2; Table 1). The activity of other regions of the DMN, such as subgenual ACC, and of the TPN, such as supplementary motor area, middle frontal gyrus, and inferior parietal lobule, also showed significant negative and positive correlations with SCR. More specifically, the findings indicate that increases in autonomic arousal relate to decreased DMN and increased TPN activity.

Functional connectivity change after regressing out SCR
Using PCC as the seed, the functional connectivity between the PCC and other brain regions showed a significant positive correlation with DMN and anticorrelation with TPN ( Fig. 3; Table 2). The magnitude of the anticorrelation between PCC (of the DMN) and frontoparietal regions (of the TPN) was stronger before in contrast to after SCR was regressed out (Fig. 4; Table 3). There was also greater PCC connectivity with other regions in the DMN (especially the medial prefrontal cortex) when SCR was not regressed out. However, under SCR modeling, there was no significant change of connectivity between PCC and AI, indicating that the intrinsic anticorrelation was not significantly affected by SCR, although the amplitudes of both regions were negatively and positively correlated with SCR, as shown in GLM analyses.

Psychophysiological interaction
Treating arousal (indexed by SCR) as the psychological context and using the BOLD signal of PCC as the physiological signal, the PPI analysis showed an enhanced positive connectivity between PCC and other brain regions in the DMN (e.g., vmPFC and PCC itself; Fig. 5, regions in green; Table 4), as well as an enhanced anticorrelation between PCC and regions in the TPN (e.g., ACC, precentral gyrus, areas near/along intraparietal sulcus). These findings are consistent with our functional connectivity analysis using the classic rs-fcMRI analytic method (Fox et al., 2005). The intrinsic connectivity between PCC and AI was not significantly modulated. Interestingly, the anticorrelation between PCC and visual cortex was significantly enhanced by SCR.

Discussion
Although complex approaches have been used to examine the functional fractionation of the DMN (Andrews- Hanna et al., 2010), and RSNs more generally (Deco et al., 2011), task-based comparisons to spontaneous resting activity (Andrews-Hanna et al., 2010) and computational simulations (Deco et al., 2011) cannot provide an objective index of mental state during rest. This limitation raises questions about possible psychological processes that may relate and/or contribute to RSN activity and connectivity. Some more popular hypotheses about the functional basis of RSNs have included the generation of spontaneous thoughts (Mason et al., 2007) and self-relevant mental simulations (Buckner et al., 2008), as well as predictions about and preparation for environmental demands (Deco et al., 2011). Previous associations of SCR with the TPN suggest that SCR has some putative contributions to consciousness and bodily awareness/interoception (Craig, 2009). Thus SCR, in relation to spontaneous brain activity, is an excellent index of homeostatic monitoring, some of Regions in green indicate that the contribution of PCC to those regions is increased by SCR, or the response of those regions to SCR is increased due to the contribution of PCC. Regions in red indicate an enhanced anticorrelation between PCC and these regions by SCR.
which may be conscious. In addition to being a good index of autonomic activity, the nonspecific SCRs during "rest" (in contrast to the specific SCR related to a stimulus) occur about 1-3 times per minute (Boucsein, 1992), an oscillation occurring within the typically used filter band of 0.01-0.08 Hz for rs-fcMRI analysis (Biswal et al., 1995). Our findings suggest that autonomic arousal, and the bodily and psychophysiological states that it reflects, have a critical relationship to the intrinsic properties of RSNs. Combinations of physiological recording and neuroimaging have not only permitted scientists to index emotion objectively, but also to identify the relationships between bodily arousal, brain activity, and possibly mental states (Harrison et al., 2010). With a few key exceptions (e.g., affective neuroscience), psychophysiological measures are most commonly used as a means to remove noise from the brain's signal in neuroimaging methodology (Critchley et al., 2011). In contrast to the usual effort to exclude psychophysiological measures from the brain activity as noise, we present evidence that general ANS activity is significantly related to spontaneous BOLD activity. The fact that activity of the key nodes of the DMN and the anticorrelated TPN is associated with SCR suggests that a significant portion of RSN activity may be linked to monitoring internal bodily and psychophysiological states (Thompson and Varela, 2001;Craig, 2002Craig, , 2009. Previous associations of SCR with the AI and ACC have led some to conclude that activity of DMN and TPN may reflect a dynamic relationship between externally and internally (especially in relation to interoception) focused attention (Critchley et al., 2011). The AI integrates high-order cognitive, sensory, and interoceptive signals, contributing to the conditions of core af- fective feeling and/or subjective awareness (Craig, 2009). The ACC, in direct coupling with the AI, functions to implement neurobehavioral responses consistent with subjective experience (Medford and Critchley, 2010). The AI and ACC are also implicated in the monitoring and control, respectively, of autonomic activity (Critchley et al., 2011). Thus, the observed positive correlation between SCR and the TPN, as well as the decrease in PCC connectivity after SCR-related activity is removed (i.e., identified as noise), suggest increased attentional coupling with interoception during rest. When external demands are limited, it seems that spontaneous activity of the brain, in part, may function to increase processing (both unconscious and conscious) of the state of the body. In contrast to approaches that consider bodily processes as noise in relation to activity of the brain, our findings provide support for an embodied mind (Thompson and Varela, 2001), a view that necessitates recognition of a dynamic relationship among mental states, bodily functions, and brain processes (Thompson and Varela, 2001;Harrison et al., 2010;Critchley et al., 2011;Deco et al., 2011). The idea that limited external demands may increase focus on internal bodily processes and psychophysiological state is somewhat consistent with a previous examination of SCR and taskrelated brain activity (Patterson et al., 2002). However, SCR may not reflect a single physiological or psychological process (Boucsein, 1992;Critchley, 2002). Patterson et al. (2002) have previously shown that changes in SCR are correlated with activation of DMN regions (e.g., vmPFC) that influence the sympathetic nervous system, independent of the tasks being performed. Consistent with their overall theory, our findings further showed that decreased activation and enhanced connectivity of vmPFC and PCC in DMN and activation of ACC/AI in TPN are related to SCR. The broad involvement of brain regions/networks suggests that SCR may reflect a multifaceted psychophysiological response and interaction. Given the present findings, along with others that suggest that autonomic activity can impact RSNs (Shmueli et al., 2007), as well as the potential relationship between autonomic activity and consciousness, important contributions from the body and mind to RSNs are potentially being overlooked in current analyses of spontaneous fluctuations of the brain.
Although we speculate that bodily arousal shapes neural responses of the RSN, it is also possible that activity of this network produces the arousal changes themselves, or a third factor, such as cognitive processes, drives both. Given that our methods are correlational, this question could not be answered in the current study and calls for further investigation. However, the present results may have implications for existing resting-state studies that have not taken arousal into account, especially those that compare DMN across health and disease states. For example, the abnormal connectivity patterns of DMN in patients with autism spectrum disorders might be related to deficits of ANS, creating the classic third variable problem.
There is also the possibility that motion is a contributing variable (i.e., motion may be a factor in some presently observed patterns). Our findings on this issue were mixed. In GLM, after the motion-related effects were regressed out (using the six motion correction parameters in the modeling), there was no activation associated with SCR in the precentral gyrus. Actually, SCR was negatively correlated with the activation of precentral gyrus and positively correlated with activation of the caudate, indicating the possibility of motor inhibition. However, our PPI analysis showed an enhanced anticorrelation between PCC and motor areas, which may suggest a contribution of motion, though it may also reflect top-down control.
At the very least, our findings, along with the findings of others (Shmueli et al., 2007;Birn et al., 2008), suggest considerable benefit of measuring and incorporating (rather than excluding) psychophysiological signals (especially ANS activity) in rs-fcMRI and fMRI. There are likely numerous sources driving the spontaneous fluctuations of the brain, some of which undoubtedly relate to neural activity (Mantini et al., 2007) and anatomical structure (Honey et al., 2009). However, methods which exclude signals from the body (Birn et al., 2008;Chang and Glover, 2009;Iacovella and Hasson, 2011), or rely mostly on reverse inference (Deco et al., 2011), are limiting the means by which we can understand the causal and functional basis of RSNs. SCR is not the only, nor necessarily the best, means to index psychophysiological states (Critchley et al., 2011). Although the present findings suggest an important role of the ANS in RSNs, the ANS has both sensory (afferent) and motor (efferent) subsystems, which are difficult to tease apart (Harrison et al., 2010). Further, SCR is but one of a variety of measures (e.g., electrodermal, cardiovascular, and blood pressure responses) to assess ANS activity. Given the ability of SCR to predict RSN activity, inclusion of multiple physiological measures may increase the amount of RSN variance that can be predicted. Using pattern classification analysis with multiple physiological indices may also help to identify the probable psychophysiological states of the subject during rest (Stephens et al., 2010). Such a program of research could provide important insight into the sources and functions of RSNs.
Finally, we would be remiss to ignore some additional questions raised by our findings. There are multiple factors that may underlie the pattern of correlations between SCR and rs-fcMRI, such as autonomic arousal, interoceptive awareness, volitional and cognitive intent, and perhaps others. A simplistic explanation for the present association between SCR fluctuations and the TPN may have to do with suppression of natural urges during rest. ACC and AI have previously been shown to help control and suppress spontaneous blinking (Lerner et al., 2009). Thus, our findings may simply reflect participants' efforts to remain still during the resting task. Even if true, this conclusion would still suggest a relation between physiological states during rest and the correlational patterns of RSNs. Alternatively, the present findings may relate to the emergence of subjective awareness, in part related to interoception, and subsequent implementation of consistent neurobehavioral responses (Medford and Critchley, 2010). If true, the latter explanation would relate to empirical evidence that the brain is one part of a dynamic homeostatic system, which includes the body and its environmental context (Thompson and Varela, 2001). This evidence advocates not only for more inclusive and objective means to measure the functional role of and contributions to RSNs, but also for consideration of the mind/brain/body relationship, both in the context of RSNs and in cognitive neuroscience more generally.