Parameter-free extraction of the functional network topology from intracranial EEG recordings: a time-resolved study of graph properties in focal onset seizures

Rationale: Focal onset epileptic seizures are due to abnormal interactions between distributed brain areas. By estimating the cross-correlation matrix of multi-site intra-cerebral EEG recordings (iEEG), one can quantify these interactions. To assess the topology of the underlying functional network, the binary connectivity matrix has to be derived from the cross-correlation matrix by use of a threshold. Classically, a unique threshold is used that constrains the topology [1]. Our method aims to set the threshold in a data-driven way by separating genuine from random cross-correlation. We compare our approach to the fixed threshold method and study the dynamics of the functional topology. Methods: We investigate the iEEG of patients suffering from focal onset seizures who underwent evaluation for the possibility of surgery. The equal-time cross-correlation matrices are evaluated using a sliding time window. We then compare 3 approaches assessing the corresponding binary networks. For each time window: * Our parameter-free method derives from the cross-correlation strength matrix (CCS)[2]. It aims at disentangling genuine from random correlations (due to finite length and varying frequency content of the signals). In practice, a threshold is evaluated for each pair of channels independently, in a data-driven way. * The fixed mean degree (FMD) uses a unique threshold on the whole connectivity matrix so as to ensure a user defined mean degree. * The varying mean degree (VMD) uses the mean degree of the CCS network to set a unique threshold for the entire connectivity matrix. * Finally, the connectivity (c), connectedness (given by k, the number of disconnected sub-networks), mean global and local efficiencies (Eg, El, resp.) are computed from FMD, CCS, VMD, and their corresponding random and lattice networks. Results: Compared to FMD and VMD, CCS networks present: *topologies that are different in terms of c, k, Eg and El. *from the pre-ictal to the ictal and then post-ictal period, topological features time courses that are more stable within a period, and more contrasted from one period to the next. For CCS, pre-ictal connectivity is low, increases to a high level during the seizure, then decreases at offset. k shows a ‘‘U-curve’’ underlining the synchronization of all electrodes during the seizure. Eg and El time courses fluctuate between the corresponding random and lattice networks values in a reproducible manner. Conclusions: The definition of a data-driven threshold provides new insights into the topology of the epileptic functional networks.

Subsecond changes of global brain state in illusory multistable motion perception

This study explored transient changes in EEG microstates and spatial Omega complexity associated with changes in multistable perception. 21-channel EEG was recorded from 13 healthy subjects viewing an alternating dot pattern that induced illusory motion with ambiguous direction. Baseline epochs with stable motion direction were compared to epochs immediately preceding stimuli that were perceived with changed motion direction ('reference stimuli'). About 750 ms before reference stimuli, Omega complexity decreased as compared to baseline, and two of four classes of EEG microstates changed their probability of occurrence. About 300 ms before reference stimuli, Omega complexity increased and the previous deviations of EEG microstates were reversed. Given earlier results on Omega complexity and microstates, these sub-second EEG changes might parallel longer-lasting fluctuations in vigilance. Assumedly, the discontinuities of illusory motion thus occur during sub-second dips in arousal, and the following reconstruction of the illusion coincides with a state of relative over-arousal.

Decreased functional connectivity of EEG theta-frequency activity in first-episode, neuroleptic-naïve patients with schizophrenia: preliminary results

We explored and refined the hypothesis that during a first episode of acute schizophrenia a disorganization of brain functioning is present. A novel EEG measure was introduced, Global Field Synchronization (GFS), that estimates functional connectivity of brain processes in different EEG frequency bands. The measure was applied to EEG's from 11 never-treated, first-episode, young patients with an acute, positive, schizophrenic symptomatology and from 19 controls, residing in Bern, Switzerland. In comparison to age- and sex- matched controls, patients had significantly decreased GFS in the theta EEG frequency band, indicating a loosened functional connectivity of processes in this frequency. The result was confirmed in an independent, comparable patient group from Osaka, Japan (9 patients and 9 controls), thus making a total of 20 analyzed patients. Previous EEG research in healthy, awake subjects indicated a positive correlation of theta activity with memory functions. Thus, our result suggests a loss of mutual interdependence of memory functions in patients with acute schizophrenia, which agrees well with previous reports of working memory dysfunction in schizophrenia.

Topographic time-frequency decomposition of the EEG

Frequency-transformed EEG resting data has been widely used to describe normal and abnormal brain functional states as function of the spectral power in different frequency bands. This has yielded a series of clinically relevant findings. However, by transforming the EEG into the frequency domain, the initially excellent time resolution of time-domain EEG is lost. The topographic time-frequency decomposition is a novel computerized EEG analysis method that combines previously available techniques from time-domain spatial EEG analysis and time-frequency decomposition of single-channel time series. It yields a new, physiologically and statistically plausible topographic time-frequency representation of human multichannel EEG. The original EEG is accounted by the coefficients of a large set of user defined EEG like time-series, which are optimized for maximal spatial smoothness and minimal norm. These coefficients are then reduced to a small number of model scalp field configurations, which vary in intensity as a function of time and frequency. The result is thus a small number of EEG field configurations, each with a corresponding time-frequency (Wigner) plot. The method has several advantages: It does not assume that the data is composed of orthogonal elements, it does not assume stationarity, it produces topographical maps and it allows to include user-defined, specific EEG elements, such as spike and wave patterns. After a formal introduction of the method, several examples are given, which include artificial data and multichannel EEG during different physiological and pathological conditions.

Multichannel EEG fields during and without visual input: frequency domain model source locations and dimensional complexities

27-Channel EEG potential map series were recorded from 12 normals with closed and open eyes. Intracerebral dipole model source locations in the frequency domain were computed. Eye opening (visual input) caused centralization (convergence and elevation) of the source locations of the seven frequency bands, indicative of generalized activity; especially, there was clear anteriorization of α-2 (10.5–12 Hz) and β-2 (18.5–21 Hz) sources (α-2 also to the left). Complexity of the map series' trajectories in state space (assessed by Global Dimensional Complexity and Global OMEGA Complexity) increased significantly with eye opening, indicative of more independent, parallel, active processes. Contrary to PET and fMRI, these results suggest that brain activity is more distributed and independent during visual input than after eye closing (when it is more localized and more posterior).

Prestimulus EEG microstates influence visual event-related potential microstates in field maps with 47 channels

The influence of the immediate prestimulus EEG microstate (sub-second epoch of stable topography/map landscape) on the map landscape of visually evoked 47-channel event-related potential (ERP) microstates was examined using the frequent, non-target stimuli of a cognitive paradigm (12 volunteers). For the two most frequent prestimulus microstate classes (oriented left anterior-right posterior and right anterior-left posterior), ERP map series were selectively averaged. The post-stimulus ERP grand average map series was segmented into microstates; 10 were found. The centroid locations of positive and negative map areas were extracted as landscape descriptors. Significant differences (MANOVAs and t-tests) between the two prestimulus classes were found in four of the ten ERP microstates. The relative orientation of the two ERP microstate classes was the same as prestimulus in some ERP microstates, but reversed in others. — Thus, brain electric microstates at stimulus arrival influence the landscapes of the post-stimulus ERP maps and therefore, information processing; prestimulus microstate effects differed for different post-stimulus ERP microstates.

Single-dose piracetam effects on global complexity measures of human spontaneous multichannel EEG

Global complexity of 47-channel resting electroencephalogram (EEG) of healthy young volunteers was studied after intake of a single dose of a nootropic drug (piracetam, Nootropil® UCB Pharma) in 12 healthy volunteers. Four treatment levels were used: 2.4, 4.8, 9.6 g piracetam and placebo. Brain electric activity was assessed through Global Dimensional Complexity and Global Omega-Complexity as quantitative measures of the complexity of the trajectory of multichannel EEG in state space. After oral ingestion (1–1.5 h), both measures showed significant decreases from placebo to 2.4 g piracetam. In addition, Global Dimensional Complexity showed a significant return to placebo values at 9.6 g piracetam. The results indicate that a single dose of piracetam dose-dependently affects the spontaneous EEG in normal volunteers, showing effects at the lowest treatment level. The decreased EEG complexity is interpreted as increased cooperativity of brain functional processes.

EEG microstate duration and syntax in acute, medication-naive, first-episode schizophrenia: a multi-center study

In young, first-episode, productive, medication-naive patients with schizophrenia, EEG microstates (building blocks of mentation) tend to be shortened. Koenig et al. [Koenig, T., Lehmann, D., Merlo, M., Kochi, K., Hell, D., Koukkou, M., 1999. A deviant EEG brain microstate in acute, neuroleptic-naïve schizophrenics at rest. European Archives of Psychiatry and Clinical Neuroscience 249, 205–211] suggested that shortening concerned specific microstate classes. Sequence rules (microstate concatenations, syntax) conceivably might also be affected. In 27 patients of the above type and 27 controls, from three centers, multichannel resting EEG was analyzed into microstates using k-means clustering of momentary potential topographies into four microstate classes (A–D). In patients, microstates were shortened in classes B and D (from 80 to 70 ms and from 94 to 82 ms, respectively), occurred more frequently in classes A and C, and covered more time in A and less in B. Topography differed only in class B where LORETA tomography predominantly showed stronger left and anterior activity in patients. Microstate concatenation (syntax) generally were disturbed in patients; specifically, the class sequence A→C→D→A predominated in controls, but was reversed in patients (A→D→C→A). In schizophrenia, information processing in certain classes of mental operations might deviate because of precocious termination. The intermittent occurrence might account for Bleuler's “double bookkeeping.” The disturbed microstate syntax opens a novel physiological comparison of mental operations between patients and controls.

Spatio-temporal dynamics of alpha brain electric fields, and cognitive modes

Brian electric activity is viewed as sequences of momentary maps of potential distribution. Frequency-domain source modeling, estimation of the complexity of the trajectory of the mapped brain field distributions in state space, and microstate parsing were used as analysis tools. Input-presentation as well as task-free (spontaneous thought) data collection paradigms were employed. We found: Alpha EEG field strength is more affected by visualizing mentation than by abstract mentation, both input-driven as well as self-generated. There are different neuronal populations and brain locations of the electric generators for different temporal frequencies of the brain field. Different alpha frequencies execute different brain functions as revealed by canonical correlations with mentation profiles. Different modes of mentation engage the same temporal frequencies at different brain locations. The basic structure of alpha electric fields implies inhomogeneity over time — alpha consists of concatenated global microstates in the sub-second range, characterized by quasi-stable field topographies, and rapid transitions between the microstates. In general, brain activity is strongly discontinuous, indicating that parsing into field landscape-defined microstates is appropriate. Different modes of spontaneous and induced mentation are associated with different brain electric microstates; these are proposed as candidates for psychophysiological ``atoms of thought''.

EEG microstates associated with salience and frontoparietal networks in frontotemporal dementia, schizophrenia and Alzheimer's disease

OBJECTIVE: There are relevant links between resting-state fMRI networks, EEG microstate classes and psychopathological alterations in mental disorders associated with frontal lobe dysfunction. We hypothesized that a certain microstate class, labeled C and correlated with the salience network, was impaired early in frontotemporal dementia (FTD), and that microstate class D, correlated with the frontoparietal network, was impaired in schizophrenia. METHODS: We measured resting EEG microstate parameters in patients with mild FTD (n = 18), schizophrenia (n = 20), mild Alzheimer's disease (AD; n = 19) and age-matched controls (old n = 19, young n = 18) to investigate neuronal dynamics at the whole-brain level. RESULTS: The duration of class C was significantly shorter in FTD than in controls and AD, and the duration of class D was significantly shorter in schizophrenia than in controls, FTD and AD. Transition analysis showed a reversed sequence of activation of classes C and D in FTD and schizophrenia patients compared with that in controls, with controls preferring transitions from C to D, and patients preferring D to C. CONCLUSION: The duration and sequence of EEG microstates reflect specific aberrations of frontal lobe functions in FTD and schizophrenia. SIGNIFICANCE: This study highlights the importance of subsecond brain dynamics for understanding of psychiatric disorders.

Low resolution brain electromagnetic tomography (LORETA) functional imaging in acute, neuroleptic-naive, first-episode, productive schizophrenia

Functional imaging of brain electrical activity was performed in nine acute, neuroleptic-naive, first-episode, productive patients with schizophrenia and 36 control subjects. Low-resolution electromagnetic tomography (LORETA, three-dimensional images of cortical current density) was computed from 19-channel of electroencephalographic (EEG) activity obtained under resting conditions, separately for the different EEG frequencies. Three patterns of activity were evident in the patients: (1) an anterior, near-bilateral excess of delta frequency activity; (2) an anterior-inferior deficit of theta frequency activity coupled with an anterior-inferior left-sided deficit of alpha-1 and alpha-2 frequency activity; and (3) a posterior-superior right-sided excess of beta-1, beta-2 and beta-3 frequency activity. Patients showed deviations from normal brain activity as evidenced by LORETA along an anterior-left-to-posterior-right spatial axis. The high temporal resolution of EEG makes it possible to specify the deviations not only as excess or deficit, but also as inhibitory, normal and excitatory. The patients showed a dis-coordinated brain functional state consisting of inhibited prefrontal/frontal areas and simultaneously overexcited right parietal areas, while left anterior, left temporal and left central areas lacked normal routine activity. Since all information processing is brain-state dependent, this dis-coordinated state must result in inadequate treatment of (externally or internally generated) information.

Faces and emotions: brain electric field sources during covert emotional processing

Covert brain activity related to task-free, spontaneous (i.e. unrequested), emotional evaluation of human face images was analysed in 27-channel averaged event-related potential (ERP) map series recorded from 18 healthy subjects while observing random sequences of face images without further instructions. After recording, subjects self-rated each face image on a scale from “liked” to “disliked”. These ratings were used to dichotomize the face images into the affective evaluation categories of “liked” and “disliked” for each subject and the subjects into the affective attitudes of “philanthropists” and “misanthropists” (depending on their mean rating across images). Event-related map series were averaged for “liked” and “disliked” face images and for “philanthropists” and “misanthropists”. The spatial configuration (landscape) of the electric field maps was assessed numerically by the electric gravity center, a conservative estimate of the mean location of all intracerebral, active, electric sources. Differences in electric gravity center location indicate activity of different neuronal populations. The electric gravity center locations of all event-related maps were averaged over the entire stimulus-on time (450 ms). The mean electric gravity center for disliked faces was located (significant across subjects) more to the right and somewhat more posterior than for liked faces. Similar differences were found between the mean electric gravity centers of misanthropists (more right and posterior) and philanthropists. Our neurophysiological findings are in line with neuropsychological findings, revealing visual emotional processing to depend on affective evaluation category and affective attitude, and extending the conclusions to a paradigm without directed task.

Affective attitudes to face images associated with intracerebral EEG source location before face viewing

We investigated whether different, personality-related affective attitudes are associated with different brain electric field (EEG) sources before any emotional challenge (stimulus exposure). A 27-channel EEG was recorded in 15 subjects during eyes-closed resting. After recording, subjects rated 32 images of human faces for affective appeal. The subjects in the first (i.e., most negative) and fourth (i.e., most positive) quartile of general affective attitude were further analyzed. The EEG data (mean=25±4.8 s/subject) were subjected to frequency-domain model dipole source analysis (FFT-Dipole-Approximation), resulting in 3-dimensional intracerebral source locations and strengths for the delta–theta, alpha, and beta EEG frequency band, and for the full range (1.5–30 Hz) band. Subjects with negative attitude (compared to those with positive attitude) showed the following source locations: more inferior for all frequency bands, more anterior for the delta–theta band, more posterior and more right for the alpha, beta and 1.5–30 Hz bands. One year later, the subjects were asked to rate the face images again. The rating scores for the same face images were highly correlated for all subjects, and original and retest affective mean attitude was highly correlated across subjects. The present results show that subjects with different affective attitudes to face images had different active, cerebral, neural populations in a task-free condition prior to viewing the images. We conclude that the brain functional state which implements affective attitude towards face images as a personality feature exists without elicitors, as a continuously present, dynamic feature of brain functioning.

Face-elicited ERPs and affective attitude: brain electric microstate and tomography analyses

Objectives: Although behavioral studies have demonstrated that normative affective traits modulate the processing of facial and emotionally charged stimuli, direct electrophysiological evidence for this modulation is still lacking. Methods: Event-related potential (ERP) data associated with personal, traitlike approach- or withdrawal-related attitude (assessed post-recording and 14 months later) were investigated in 18 subjects during task-free (i.e. unrequested, spontaneous) emotional evaluation of faces. Temporal and spatial aspects of 27 channel ERP were analyzed with microstate analysis and low resolution electromagnetic tomography (LORETA), a new method to compute 3 dimensional cortical current density implemented in the Talairach brain atlas. Results: Microstate analysis showed group differences 132-196 and 196-272 ms poststimulus, with right-shifted electric gravity centers for subjects with negative affective attitude. During these (over subjects reliably identifiable) personality-modulated, face-elicited microstates, LORETA revealed activation of bilateral occipito-temporal regions, reportedly associated with facial configuration extraction processes. Negative compared to positive affective attitude showed higher activity right temporal; positive compared to negative attitude showed higher activity left temporo-parieto-occipital. Conclusions: These temporal and spatial aspects suggest that the subject groups differed in brain activity at early, automatic, stimulus-related face processing steps when structural face encoding (configuration extraction) occurs. In sum, the brain functional microstates associated with affect-related personality features modulate brain mechanisms during face processing already at early information processing stages.

Chronic schizophrenics with positive symptomatology have shortened EEG microstate durations

OBJECTIVE: In young, first-episode, never-treated schizophrenics compared with controls, (a) generally shorter durations of EEG microstates were reported (Koukkou et al., Brain Topogr 6 (1994) 251; Kinoshita et al., Psychiatry Res Neuroimaging 83 (1998) 58), and (b) specifically, shorter duration of a particular class of microstates (Koenig et al., Eur Arch Psychiatry Clin Neurosci 249 (1999) 205). We now examined whether older, chronic schizophrenic patients with positive symptomatology also show these characteristics. METHODS: Multichannel resting EEG (62.2 s/subject) from two subject groups, 14 patients (36.1+/-10.2 years old) and 13 controls (35.1+/-8.2 years old), all males, was analyzed into microstates using a global approach for microstate analysis that clustered the microstates into 4 classes (Koenig et al., 1999). RESULTS: (a) Hypothesis testing of general microstate shortening supported a trend (P=0.064). (b) Two-way repeated measure ANOVA (two subject groupsx4 microstate classes) showed a significant group effect for microstate duration. Posthoc tests revealed that a microstate class with brain electric field orientation from left central to right central-posterior had significantly shorter microstates in patients than controls (68.5 vs. 76.1 ms, P=0.034). CONCLUSIONS: The results were in line with the results from young, never-treated, productive patients, thus suggesting that in schizophrenic information processing, one class of mental operations might intermittently cause deviant mental constructs because of premature termination of processing.