Top-Down versus Bottom-up Processing in the Human Brain: Distinct Directional Influences Revealed by Integrating SOBI and Granger Causality

Similar documents
DATA MANAGEMENT & TYPES OF ANALYSES OFTEN USED. Dennis L. Molfese University of Nebraska - Lincoln

From Single-trial EEG to Brain Area Dynamics

Contrasting single-trial ERPs between experimental manipulations: Improving differentiability by blind source separation

Beyond Blind Averaging: Analyzing Event-Related Brain Dynamics. Scott Makeig. sccn.ucsd.edu

Feasibility Study of the Time-variant Functional Connectivity Pattern during an Epileptic Seizure

Decisions Have Consequences

From single-trial EEG to brain area dynamics

CS/NEUR125 Brains, Minds, and Machines. Due: Friday, April 14

Est-ce que l'eeg a toujours sa place en 2019?

The neurolinguistic toolbox Jonathan R. Brennan. Introduction to Neurolinguistics, LSA2017 1

Sixth International Conference on Emerging trends in Engineering and Technology (ICETET'16) ISSN: , pp.

On the tracking of dynamic functional relations in monkey cerebral cortex

EEG changes accompanying learned regulation of 12-Hz EEG activity

Oscillations: From Neuron to MEG

Processed by HBI: Russia/Switzerland/USA

WAVELET ENERGY DISTRIBUTIONS OF P300 EVENT-RELATED POTENTIALS FOR WORKING MEMORY PERFORMANCE IN CHILDREN

Neural Correlates of Human Cognitive Function:

Biomedical Research 2013; 24 (3): ISSN X

Nature Neuroscience: doi: /nn Supplementary Figure 1. Large-scale calcium imaging in vivo.

The Sonification of Human EEG and other Biomedical Data. Part 3

Emotion Detection Using Physiological Signals. M.A.Sc. Thesis Proposal Haiyan Xu Supervisor: Prof. K.N. Plataniotis

A Brain Computer Interface System For Auto Piloting Wheelchair

Toward a noninvasive automatic seizure control system with transcranial focal stimulations via tripolar concentric ring electrodes

Simultaneous Real-Time Detection of Motor Imagery and Error-Related Potentials for Improved BCI Accuracy

Effects of Light Stimulus Frequency on Phase Characteristics of Brain Waves

Matrix Energetics Research Brainwaves and Heart waves Research on Matrix Energetics in Action

Human Brain Institute Russia-Switzerland-USA

USING AUDITORY SALIENCY TO UNDERSTAND COMPLEX AUDITORY SCENES

Event Related Potentials: Significant Lobe Areas and Wave Forms for Picture Visual Stimulus

Sum of Neurally Distinct Stimulus- and Task-Related Components.

Brain Computer Interface. Mina Mikhail

EEG in the ICU: Part I

ABSTRACT 1. INTRODUCTION 2. ARTIFACT REJECTION ON RAW DATA

DETECTION AND CORRECTION OF EYE BLINK ARTIFACT IN SINGLE CHANNEL ELECTROENCEPHALOGRAM (EEG) SIGNAL USING A SIMPLE k-means CLUSTERING ALGORITHM

Attention Response Functions: Characterizing Brain Areas Using fmri Activation during Parametric Variations of Attentional Load

Investigations in Resting State Connectivity. Overview

Amy Kruse, Ph.D. Strategic Analysis, Inc. LCDR Dylan Schmorrow USN Defense Advanced Research Projects Agency

EEG/ERP Control Process Separation

Competing Streams at the Cocktail Party

Multichannel Classification of Single EEG Trials with Independent Component Analysis

This presentation is the intellectual property of the author. Contact them for permission to reprint and/or distribute.

Neuro Q no.2 = Neuro Quotient

PSD Analysis of Neural Spectrum During Transition from Awake Stage to Sleep Stage

Power-Based Connectivity. JL Sanguinetti

HST 583 fmri DATA ANALYSIS AND ACQUISITION

Outline of Talk. Introduction to EEG and Event Related Potentials. Key points. My path to EEG

Robust time-varying multivariate coherence estimation: Application to electroencephalogram recordings during general anesthesia

Seizure onset can be difficult to asses in scalp EEG. However, some tools can be used to increase the seizure onset activity over the EEG background:

Music-induced Emotions and Musical Regulation and Emotion Improvement Based on EEG Technology

Principles of Science

Fusing Simultaneous EEG-fMRI by Linking Multivariate Classifiers

Assessing Functional Neural Connectivity as an Indicator of Cognitive Performance *

Visual Context Dan O Shea Prof. Fei Fei Li, COS 598B

Novel single trial movement classification based on temporal dynamics of EEG

The impact of numeration on visual attention during a psychophysical task; An ERP study

EEG Analysis on Brain.fm (Focus)

PARAFAC: a powerful tool in EEG monitoring

Submitted report on Sufi recordings at AAPB 2013 in Portland. Not for general distribution. Thomas F. Collura, Ph.D. July, 2013

Functional Connectivity and the Neurophysics of EEG. Ramesh Srinivasan Department of Cognitive Sciences University of California, Irvine

Automatic Identification and Removal of Ocular Artifacts from EEG using Wavelet Transform

Development of 2-Channel Eeg Device And Analysis Of Brain Wave For Depressed Persons

Discrimination between ictal and seizure free EEG signals using empirical mode decomposition

EEG-Rhythm Dynamics during a 2-back Working Memory Task and Performance

THE PENNSYLVANIA STATE UNIVERSITY SCHREYER HONORS COLLEGE DEPARTMENT OF INTERCOLLEGE PROGRAM PERSONALITY AND THE DEFAULT MODE NETWORK: AN EEG STUDY

Neurophysiologically Driven Image Triage: A Pilot Study

PCA Enhanced Kalman Filter for ECG Denoising

Supplemental Digital Content 1: Supplemental Results

Supplementary Figure 1

To link to this article: PLEASE SCROLL DOWN FOR ARTICLE

Rhythm and Rate: Perception and Physiology HST November Jennifer Melcher

Introduction to Computational Neuroscience

Word Length Processing via Region-to-Region Connectivity

1 Decoding attentional orientation from EEG spectra. Ramesh Srinivasan, Samuel Thorpe, Siyi Deng, Tom Lappas & Michael D'Zmura

Electroencephalographic Study of Essential Oils for Stress Relief

Title of Thesis. Study on Audiovisual Integration in Young and Elderly Adults by Event-Related Potential

Supplementary materials for: Executive control processes underlying multi- item working memory

Selection of Feature for Epilepsy Seizer Detection Using EEG

EDGE DETECTION. Edge Detectors. ICS 280: Visual Perception

CHAPTER 6 INTERFERENCE CANCELLATION IN EEG SIGNAL

1. Introduction

AUTOCORRELATION AND CROSS-CORRELARION ANALYSES OF ALPHA WAVES IN RELATION TO SUBJECTIVE PREFERENCE OF A FLICKERING LIGHT

Supplementary Figure 1. Example of an amygdala neuron whose activity reflects value during the visual stimulus interval. This cell responded more

John E. Richards* Department of Psychology, University of South Carolina, Columbia, SC 29208, United States

Brain and Cognition. Cognitive Neuroscience. If the brain were simple enough to understand, we would be too stupid to understand it

Frequency Tracking: LMS and RLS Applied to Speech Formant Estimation

Supplementary material

Classification of EEG signals in an Object Recognition task

Neurophysiologically Driven Image Triage: A Pilot Study

ISSN: (Online) Volume 3, Issue 7, July 2015 International Journal of Advance Research in Computer Science and Management Studies

the body and the front interior of a sedan. It features a projected LCD instrument cluster controlled

Lateral Geniculate Nucleus (LGN)

An Introduction to Translational Neuroscience Approaches to Investigating Autism

The EEG Analysis of Auditory Emotional Stimuli Perception in TBI Patients with Different SCG Score

Neurophysiology & EEG

Diagnosing Complicated Epilepsy: Mapping of the Epileptic Circuitry. Michael R. Sperling, M.D. Thomas Jefferson University Philadelphia, PA

Quiroga, R. Q., Reddy, L., Kreiman, G., Koch, C., Fried, I. (2005). Invariant visual representation by single neurons in the human brain, Nature,

Proceedings of Meetings on Acoustics

EEG reveals divergent paths for speech envelopes during selective attention

Support Vector Machine Classification and Psychophysiological Evaluation of Mental Workload and Engagement of Intuition- and Analysis-Inducing Tasks

Transcription:

Top-Down versus Bottom-up Processing in the Human Brain: Distinct Directional Influences Revealed by Integrating SOBI and Granger Causality Akaysha C. Tang 1, Matthew T. Sutherland 1, Peng Sun 2, Yan Zhang 3, Masato Nakazawa 1, Amy Korzekwa 1, Zhen Yang 1, Mingzhou Ding 3 1 University of New Mexico, Department of Psychology, Albuquerque, NM, USA 2 University of New Mexico, Department of Electrical Engineering, Albuquerque, NM, USA 3 University of Florida, J. Crayton Pruitt Family Department of Biomedical Engineering, Gainesville, FL, USA akaysha@unm.edu, msuther@unm.edu pengsun@unm.edu, zhyang @unm.edu, maliszt@unm.edu, akorzekw@unm.edu, mding@bme.flu.edu Abstract. Top-down and bottom-up processing are two distinct yet highly interactive modes of neuronal activity underlying normal and abnormal human cognition. Here we characterize the dynamic processes that contribute to these two modes of cognitive operation. We used a blind source separation algorithm called second-order blind identification (SOBI [1]) to extract from high-density scalp EEG (128 channels) two components that index neuronal activity in two distinct local networks: one in the occipital lobe and one in the frontal lobe. We then applied Granger causality analysis to the SOBI-recovered neuronal signals from these two local networks to characterize feed-forward and feedback influences between them. With three repeated observations made at least one week apart, we show that feed-forward influence is dominated by alpha while feedback influence is dominated by theta band activity and that this directionselective dominance pattern is jointly modulated by situational familiarity and demand for visual processing. Keywords: electroencephalogram, second-order blind identification (SOBI), coherence, Granger causality, top-down, bottom-up, feed-forward, feedback. 1 Introduction Second-order blind identification (SOBI) [1] is an emerging signal processing technique that can be used to facilitate source analysis from high-density EEG. Similar to other ICA algorithms that have been applied to EEG data [2], [3], SOBI can be used to isolate and remove ocular artifact [4]. In our laboratory, we have conducted extensive investigations to demonstrate the utility of SOBI in aiding source analysis from high-density EEG. Specifically, we have shown that: (1) SOBI can correctly recover known noise sources (noisy sensors and artificially injected noise at known electrodes) and known neuronal sources (SI activation by median nerve

stimulation) [5]; (2) SOBI can increase signal to noise ratios leading to improved performance in single-trial ERP classification [6]; (3) SOBI can recover neuronal sources whose activations are correlated [7]; (4) SOBI can recover neuronal sources using EEG collected when the brain is in its default mode (i.e., the resting state) [8]; (5) SOBI can recover neuronal sources during free viewing of continuous streams of visual information [9]; and (6) SOBI can recover weak neuronal signals that temporally overlap with much stronger signals (e.g. signals associated with ipsilateral activation of primary somatosensory cortex) [10]. In this paper, we set out to achieve three goals. First, we seek to provide further validation for SOBI recovered neuronal sources by investigating whether the same neuronal sources can be recovered from repeated EEG measures that are obtained days and weeks apart. Second, we combine SOBI with Granger causality analysis to show distinct patterns of theta(θ)/alpha(α) contributions in the feed-forward and feedback influences between the frontal and occipital cortices. Third, we investigate how such asymmetrical influence between the frontal and occipital cortices is modulated by sensory processing and by situational familiarity. 2 Methods Eight right-handed subjects volunteered to participate in the present study. All subjects were free of any history of neurological or psychological disorders. The experimental procedures were conducted in accordance with the Human Research Review Committee at the University of New Mexico. Each subject was tested in three sessions at Week 0, Week 1, and Week 4 or later. Up to 7 min of continuous 128- channel EEG data were collected at 1000 Hz during: (1) eyes-closed resting ; (2) eyes-open resting ; (3) video-viewing (a silently played nature video); (4) listening to only the audio track of the video; and (5) forming mental images of scenes from the video. This paper limits the discussion to conditions 1-3. SOBI was applied to the continuous EEG data x(t), across all conditions to extract the continuous time course of activation from two types of neuronal components--- an anterior (A) and a posterior (P) component. For details on SOBI application, see [5]. Briefly, SOBI recovers the underlying sources, s(t), by minimizing the sum squared cross-correlations between s i (t) and s j (t + τ), across all pairs of sources and across multiple time delays, τs. A subset of SOBI-recovered components can be verified as neuronal sources via source localization using a forward model (e.g. BESA 5.0) [3]. Here we focused our analysis on two such neuronal components that correspond to focal regions within the frontal and occipital lobes. Feed-forward (FF) and feedback (FB) influences were quantified by Granger causality between the two components, reflecting long-distance directional influences between the frontal and occipital cortices. Granger causality analysis was carried out on the continuous time courses, s i (t), from the selected A and P components according to methods detailed in [11], [12]. As Granger causality can be decomposed into its frequency content, we computed Granger causality spectrum and measured power within the θ (4-7 Hz) and α (8-14 Hz) bands using a moving window of 30-sec with 5-

sec increments. Power in the θ and α bands from the A and P components were also computed as indicators of synchronization within the local networks. 3 Results Reliable Extraction and Identification of Neuronal Components from Repeated Measures made Weeks Apart. In all 8 subjects, across all 3 sessions, we were able to recover SOBI components that corresponded to two distinct neuronal sources, one localized to a rather focal region within the frontal cortex, in or near anterior cingulate cortex (ACC) and the other to focal regions within the occipital lobe (occipital gyrus). Repeated-measure ANOVA revealed no statistically significant differences in the location of the corresponding ECD models across the 3 recording sessions. As no session-to-session difference was found, the averaged locations across the 3 sessions are shown in Fig. 1. ECDs for each of the 8 subjects are superimposed in the figure revealing a tight clustering of ECDs across subjects. This result demonstrates that SOBI can reliably recover components that correspond to anatomically well defined brain regions even when the recording sessions were made weeks apart. It is important to emphasize that the recovery of these two neuronal sources was achieved without imposing constraints of fixation or use of event-related stimulation paradigms. Instead, subjects were allowed to freely move or blink their eyes as needed during the recording conditions. No segment of the EEG data was excluded prior to SOBI application. These unique features of SOBI processing have non-trivial Fig. 1. Equivalent current dipole (ECD) locations for the SOBI recovered A and P components. implications for the study of mental disorders and the study of early development or aging where subjects are often unable to conform to typical experimental constraints.

Theoretically, this result implies that SOBI s ability to recover anatomically welldefined neuronal sources does not depend upon the use of an event-related stimulation paradigm. Thus, fast electrical brain activity in the default mode [13] can be investigated in terms of neuronal signals originating from specific, focal cortical areas. In comparison to default mode activity revealed by fmri, the default mode activity revealed with SOBI and EEG will offer millisecond temporal resolution, allowing for the characterization of default mode brain dynamics within a new temporal domain. Fig. 2. Median power spectra of two SOBI-neuronal components as a function of repeated exposures to the same experimental situation. Session 1: week 0; Session 2: week 1; Session 3: week 4+. Local Network Synchrony Shows Distinct Patterns of Change across 3 Repeated Exposures to the same Experimental Situation. For each of the 3 recording sessions, power spectra from the component time courses were computed for ~5-min segments during which the subjects had their eyes-closed (red), eyes-open (blue), or viewed a nature video (green), respectively (Fig. 2). The anterior component had peak power within the θ band while the posterior component had peak power within the α band, indicated by a significant main effect of Region in the θ-to-α ratio (F[1,7] = 52.12, p < 0.001, partial η 2 = 0.88). This is consistent with the well established fact that the posterior and anterior parts of the brain are major sources of α and θ generators, respectively.

Power spectra in these two components were differentially modulated by sessions and experimental conditions [interaction effect: Region x Session (contrast coefficients: 1, -1, 0) x Condition (1, 0, -1), F(1,7) = 3.52, p = 0.05, 1-tailed, partial η 2 =0.33)]. For the P component, the power spectra revealed a systematic effect of session and experimental condition. Across the 3 repeated exposures to the same experimental conditions, peak α power decreased as the testing situation became increasingly familiar. Across the 3 experimental conditions, the highest peak α power was associated with the eyes closed condition and the peak α power was successively reduced when the demand for visual processing increased from the eyes-closed to the eyes-open and video-viewing conditions. This latter observation is consistent with the known observation that visual processing suppresses α band activity. In contrast, for the anterior component, the power spectra showed a relative insensitivity to repeated exposures to the same experimental conditions and little modulation by the eyesclosed, eyes-open, and video-viewing conditions. Differential Modulation of θ/α Contribution to Feed-Forward and Feedback Influences by Situational Familiarity and Visual Processing. FF(posterior-to-anterior) and FB (anterior-to-posterior) influences were measured by Granger causality in the θ and α band activity separately. FF and FB Granger causality measures were plotted as a function of time (Fig. 3). For the FF influence, when the eyes were closed, α band Fig. 3. Theta dominance over alpha in the anterior-to-posterior feedback (lower) influence and its reversal in the posterior-to-anterior feed-forward influence (upper) from a single-subject. activity clearly dominated as indicated by the α waveforms (black) having greater area underneath the curve than the θ waveforms (grey). This α dominance was clearly reduced when the eyes were open and was further reduced to nearly non-existent when the subjects viewed a video. For the FB influence, the pattern of α dominance over θ was reversed showing uniform θ dominance over α across all 3 experimental conditions.

Using the area underneath the curve as a dependent measure, we summarize results from all 8 subjects across all 3 recording sessions in Fig. 4. To determine whether θ and α band activity contribute differentially to the FF and FB influences and how such differential contributions are modulated by situational familiarity and sensory processing, we performed an ANOVA on the θ/α ratio. The θ/α ratio differed significantly between the FF and FB influences with a greater ratio for FB influence than for the FF influence (main effect of Direction, F[1,7] = 34.64, p < 0.001, partial η 2 = 0.83), i.e. a θ dominance in FB influence. This can be seen by the higher measures for the θ band activity than the α band activity for the FB influences in most of the 9 conditions and clear reversal or reduction of this θ dominance in the FF influence (Fig. 4). Fig. 4. Cumulative Granger Causality (area underneath the curve in Fig. 3) in the θ and α band as a function of situational familiarity (repeated sessions) and a function of visual processing (eyes-closed, eyes-open, video-viewing). This reversal of θ dominance from FF and FB influences was significantly modulated by the familiarity of the situation [Direction x Session (contrast coefficients: 1, -1, 0), F(1,7) = 11.97, p=0.005, 1-tailed, partial η 2 =0.63]. The reversal is more prominent when the situation was novel (Week 0) than when it became more familiar (Week 1 and 4+). This is best seen in the case of eyes-closed condition. The magnitude of reversal is clearly reduced from Week 0 in comparison to Week 4+. For the eyes-open condition, the θ dominance was reversed in Week 0 and 1 and reduced in Week 4+. For the video-viewing condition, the reversal of θ dominance does not appear to be influenced by the increasing situational familiarity. These

patterns indicate that the FF/FB contrast is dependent upon the amount of visual information processing involved. When the subjects were engaged in visual perception during video-viewing, θ dominance in the FB influence and θ-α balance in the feed-forward influence are maintained across recording sessions. This visual processing-dependent effect is supported by a significant 3-way interaction [Direction x Session (1, -1, 0) x Condition (1, 0, -1), F(1,7) = 7.13, p = 0.02, 1-tailed partial η 2 =0.63]. Within Week0 when the recording situation was novel (which is comparable to most studies that do not deal with the issue of task familiarity), θ dominance in the FB influence was maintained despite varying demand for visual processing. In contrast, the α dominance in the FF influence in the case of eyes-closed condition was reduced by increasing demand for sensory processing. In fact, visual processing was accompanied not only by a reduction in α but an increase in θ band activity in the FF influence. We speculate that this increase in θ band activity serves to match the θ- dominance in the FB influence to mediate the dynamic two-way communication between the posterior and anterior parts of the brain. 4 Discussion We analyzed EEG data collected from 8 subjects in three sessions that were weeks apart, each including a period of resting with eyes-closed, resting with eyes-open, and visual perception while free viewing a nature video. We extracted neuronal signals from focal brain regions within the frontal and occipital lobes and showed that such extraction can be achieved under free viewing conditions and from recordings made weeks apart. As many intervening events must have taken place during the intersession intervals, the reliable extraction of the neuronal sources raises the possibility that such a wide range of variations may be overcome by the use of SOBI in longitudinal experimental designs necessary for developmental and aging studies. Applying Granger causality analysis to the time courses of the frontal and occipital SOBI components, we presented evidence indicating distinct patterns of θ/α band activity in the FF and FB influences between the two components, with a θ dominance characterizing the FB influence and an α dominance in the FF influence. By comparing the feed-forward and feedback influences under varying degrees of situational familiarity (sessions) and under conditions of varying degrees of visual processing (eyes-closed, eyes-open, and video viewing), we presented evidence that the balance in θ-α band activity between the FF and FB influences is modulated by two factors. First, situational familiarity can reduce the degrees of θ and α dominance in the FB and FF influences, respectively (as in the case of eyesclosed). Second, the amount of sensory processing increases the θ band contribution and decreases α band contribution to FF influence but has little effect on FB influence. Finally, situational familiarity and sensory processing jointly determine the θ-α balance. Increasing familiarity and increasing visual processing both increases θ band contribution to FF influence. In contrast, for FB influences, increasing familiarity decreases θ band contribution when there is little demand for visual

processing (eyes-closed) and has no effect on θ band contribution when there is high demand for visual processing (visual). Together, these findings demonstrate a novel non-invasive approach to the assessment of top-down and bottom-up influences in the human brain. These findings may particularly benefit those clinicians and researchers who are interested in how bottom-up and top-down influences interact in both diseased and normal brains. Future work will extend this analysis to networks involving more functionally distinct brain regions. Acknowledgments. Supported by a grant from Sandia National Labs to A.C.T. References 1. Belouchrani, A., Abed-Meraim, K., Cardoso, J.-F., Moulines, E.: A blind source separation technique using second-order statistics. IEEE Trans. Signal Process. 5 (1997) 434 444 2. Bell, A.J., and Sejnowski, T.J.: An information-maximization approach to blind separation and blind deconvolution. Neural Computation 1 (1995) 1129-1159 3. Hyvarinen, A., and Oja E.: A fast fixed-point algorithm for independent component analysis. Neural Computation 9 (1997) 1483-1492 4. Joyce, C.A., Gorodnitsky, I.F., Kutas, M.: Automatic removal of eye movement and blink artifacts from EEG data using blind component separation. Psychophysiology 41 (2004) 313 325 5. Tang, A.C., Sutherland, M.T., McKinney, C.J.: Validation of SOBI components from high density EEG. NeuroImage 25 (2005) 539 553 6. Tang, A.C., Sutherland, M.T., and Wang, Y.: Contrasting single-trial ERPs between experimental manipulations: Improving differentiability by blind source separation. NeuroImage, 29 (2006) 335-346 7. Tang, A.C., Liu, J.Y., and Sutherland, M.T.: Recovery of correlated neuronal sources from EEG: The good and bad ways of using SOBI. NeuroImage 28 (2005) 507-519 8. Sutherland, M.T., and Tang, A.C.: Blind source separation can recover systematically distributed neuronal sources from resting EEG. In: Eurasip Proceedings of the Second International Symposium on Communications, Control, and Signal Processing (ISCCSP 2006). Marrakech, Morocco; March 13-15; (2006) http://www.eurasip.org/content/eusipco/isccsp06/defevent/papers/cr1307.pdf 9. Tang, A.C., Sutherland, M.T., McKinney, C.J., Liu, J.Y., Wang, Y., Parra, L.C., Gerson, A.D., and Sajda, P.: Classifying single-trial ERPs from visual and frontal cortex during free viewing. In: IEEE Proceedings of the International Joint Conference on Neural Networks (IJCNN 2006), Vancouver, BC, Canada; July 16-21; (2006) 1376-1383 10. Sutherland, M.T., and Tang, A.C.: Reliable detection of bilateral activation in human primary somatosensory cortex by unilateral median nerve stimulation. NeuroImage, 33 (2006) 1042-1054 11. Ding, M.: Short-window spectral analysis of cortical event-related potentials by adaptive multivariate autoregressive modeling: Data preprocessing, model validation, and variability assessment. Biological Cybernetics 83 (2000) 35 45 12. Ding, M., Chen, Y., and Bressler, S.L.: Granger causality: Basic theory and application to neuroscience. In: Winterhalder, M., Schelter, B., Timmer, J. (eds.): Handbook of Time Series Analysis. Wiley, (2006) 13. Gusnard, D.A., and Raichle, M.E.: Searching for a baseline: Functional imaging and the resting human brain. Nature Review Neuroscience, 102 (2001) 685-694