Altered Electroencephalographic Resting-State Large-Scale Brain Network Dynamics in Euthymic Bipolar Disorder Patients

Background: Neuroimaging studies provided evidence for disrupted resting-state functional brain network activity in bipolar disorder (BD). Electroencephalographic (EEG) studies found altered temporal characteristics of functional EEG microstates during depressive episode within different affective disorders. Here we investigated whether euthymic patients with BD show deviant resting-state large-scale brain network dynamics as reflected by altered temporal characteristics of EEG microstates.

Methods: We used high-density EEG to explore between-group differences in duration, coverage, and occurrence of the resting-state functional EEG microstates in 17 euthymic adults with BD in on-medication state and 17 age- and gender-matched healthy controls. Two types of anxiety, state and trait, were assessed separately with scores ranging from 20 to 80.

Results: Microstate analysis revealed five microstates (A–E) in global clustering across all subjects. In patients compared to controls, we found increased occurrence and coverage of microstate A that did not significantly correlate with anxiety scores.

Conclusion: Our results provide neurophysiological evidence for altered large-scale brain network dynamics in BD patients and suggest the increased presence of A microstate to be an electrophysiological trait characteristic of BD.

Introduction

Bipolar disorder (BD) is a common and severe psychiatric disorder, with an important personal and societal burden (1, 2). The worldwide prevalence of BD is considered to range between 1 and 3% (3, 4). BD patients are frequently misdiagnosed and often identified at late stages of disease progression, which can lead to inadequate treatment (5) and worse functional prognosis (6). A better understanding of the underlying pathophysiology is needed to identify objective biomarkers of BD that would improve diagnostic and/or treatment stratification of patients.

Possible candidates for neurobiological biomarkers in BD could arise from the abnormalities of functional brain networks. Evidence from brain imaging studies consistently points to abnormalities in circuits implicated in emotion regulation and reactivity. Particularly, attenuated frontal and enhanced limbic activations are reported in BD patients (7–9). Interestingly, regions implicated in the pathophysiology of the disease, such as the inferior frontal gyrus, the medial prefrontal cortex (mPFC), and the amygdala present altered activation patterns even in unaffected first-degree relatives of BD patients (10), pointing toward brain alterations that could underlie disease vulnerability. Moreover, evidence from functional magnetic resonance imaging (fMRI) studies showed aberrant resting-state functional connectivity between frontal and meso-limbic areas in BD when compared to healthy controls (11). A recently developed functional neuroanatomic model of BD suggests, more specifically, decreased connectivity between ventral prefrontal networks and limbic brain regions including the amygdala (12, 13). The functional connectivity abnormalities in BD in brain areas associated with emotion processing were shown to vary with mood state. A resting-state functional connectivity study of emotion regulation networks demonstrated that subgenual anterior cingulate cortex (sgACC)-amygdala coupling is critically affected during mood episodes, and that functional connectivity of sgACC plays a pivotal role in mood normalization through its interactions with the ventrolateral PFC and posterior cingulate cortex (14). Nevertheless, although different fMRI metrics allowed to report deviant patterns of large-scale networks and altered resting-state functional connectivity (14, 15) in BD, the precise temporal dynamics of the functional brain networks at rest remain to be determined.

Large-scale neural networks dynamically and rapidly re-organize themselves to enable efficient functioning (16, 17). Fast dynamics of the resting-state large-scale neural networks can be studied on sub-second temporal scales with EEG microstate analysis (18–20). EEG microstates are defined as short periods (60–120 ms) of quasi-stable electric potential scalp topography (21, 22). Therefore, microstate analysis can cluster the scalp topographies of the resting-state EEG activity into the set of a few microstate classes including the four canonical classes A–D (20) and more recent additional ones (23, 24). Since each microstate class topography reflects a coherent neuronal activity (25, 20), the temporal characteristics, such as duration, occurrence, and coverage, may be linked to the expression of spontaneous mental states and be representative of the contents of consciousness (26, 27). Numerous studies reported abnormalities in temporal properties of resting-state EEG microstates in neuropsychiatric disorders (for review see 25, 20). Evidence from microstate studies suggests that altered resting-state brain network dynamics may represent a marker of risk to develop neuropsychiatric disorders (28–30), predict clinical variables of an illness (31), or help to assess the efficacy of a treatment (32, 33). Only two studies investigated resting-state EEG in BD patients (34, 35). These studies examined patients during a depressive episode within different affective disorders. Adaptive segmentation of resting-state EEG showed abnormal microstate topographies and reduced overall average microstate duration in patients that met criteria for unipolar or bipolar mood disorders or for dysthymia (34). Using a k-means cluster analysis, an increased occurrence of microstate A with depression as an effect related to the symptom severity was observed during a period of depression in unipolar and bipolar patients (35).

Trait markers of BD based on neurobiological findings can be considered as biomarkers of illness (36, 37). These trait markers of BD can be studied during the periods of remission, or euthymia. No microstate study on spontaneous activity, however, has been performed on euthymic BD patients to the best of our knowledge. Thus, the main goal of the current study was to explore group differences between euthymic patients with BD and healthy controls in terms of resting-state EEG microstate dynamics. We hypothesized that BD patients during remission will show altered temporal characteristics of EEG microstates such as duration, coverage, and occurrence.

Materials and Methods

Subjects

Data were collected from 17 euthymic adult patients with BD and 17 healthy control (HC) subjects. The patients were recruited from the Mood Disorders Unit at the Geneva University Hospital. A snowball convenience sampling was used for the selection of the BD patients. Control subjects were recruited by general advertisement. All subjects were clinically evaluated using clinical structured interview [DIGS: Diagnostic for Genetic Studies (38)]. BD was confirmed in the experimental group by the usual assessment of the specialized program, an interview with a psychiatrist, and a semi-structured interview and relevant questionnaires with a psychologist. Exclusion criteria for all participants were a history of head injury, current alcohol, or drug abuse. Additionally, a history of psychiatric or neurological illness and of any neurological comorbidity were exclusion criteria for controls and bipolar patients, respectively. Symptoms of mania and depression were evaluated using the Young Mania Rating Scale (YMRS) (39) and the Montgomery-Åsberg Depression Rating Scale (MADRS) (40), respectively. Participants were considered euthymic if they scored < 6 on YMRS and < 12 on MADRS at the time of the experiment, and were stable for at least 4 weeks before. All patients were medicated, receiving pharmacological therapy including antipsychotics, antidepressants, and mood stabilizers, and had to be under stable medication for at least 4 weeks. The experimental group included both BD I (n = 10) and BD II (n = 7) types.

To check for possible demographic or clinical differences between groups, subject characteristics such as age, education, or level of depression were compared between groups using independent t-tests. Anxiety is highly associated with BD (41, 42) and is a potential confounding variable when investigating microstate dynamics at rest. For example, decreased duration of EEG microstates at rest in patients with panic disorder has been reported (43). To check for possible differences in anxiety symptoms, all subjects were assessed with the State-Trait Anxiety Inventory (STAI) (44). Anxiety as an emotional state (state-anxiety) and anxiety as a personal characteristic (trait-anxiety) were evaluated separately. Scores of both state- and trait-anxiety range from 20 to 80, higher values indicating greater anxiety. The scores were compared between patients and controls using independent t-tests.

This study was carried out in accordance with the recommendations of the Ethics Committee for Human Research of the Geneva University Hospital, with written informed consent from all subjects. All subjects gave written informed consent in accordance with the Declaration of Helsinki. The protocol was approved by the Ethics Committee for Human Research of the Geneva University Hospital, Switzerland.

Electroencephalographic Recording and Pre-Processing

The EEG was recorded with a high density 256-channel system (EGI System 200; Electrical Geodesic Inc., OR, USA), sampling rate of 1 kHz, and Cz as acquisition reference. Subjects were sitting in a comfortable upright position and were instructed to stay as calm as possible, to keep their eyes closed and to relax for 6 min. They were asked to stay awake.

To remove muscular artifacts originating in the neck and face the data were reduced to 204 channels. Two to four minutes of EEG data were selected based on visual assessment of the artifacts and band-pass filtered between 1 and 40 Hz. Subsequently, in order to remove ballistocardiogram and oculo-motor artifacts, Infomax-based independent component analysis (45) was applied on all but one or two channels rejected due to abundant artifacts. Only components related to physiological noise, such as ballistocardiogram, saccadic eye movements, and eye blinking, were removed based on the waveform, topography, and time course of the component. The cleaned EEG recordings were down-sampled to 125 Hz and the previously identified noisy channels were interpolated using a three-dimensional spherical spline (46), and re-referenced to the average reference. All the preprocessing steps were done using MATLAB and the freely available Cartool Software 3.70 (https://sites.google.com/site/cartoolcommunity/home), programmed by Denis Brunet.

Electroencephalographic Data Analysis

To estimate the optimal set of topographies explaining the EEG signal, a standard microstate analysis was performed using k-means clustering (see Supplementary Figure 1). The polarity of the maps was ignored in this procedure (18, 47, 48). To determine the optimal number of clusters, we applied a meta-criterion that is a combination of seven independent optimization criteria (for details see 24). In order to improve the signal-to-noise ratio, only the data at the time points of the local maximum of the Global Field Power (GFP) were clustered (18, 22, 28, 49). The GFP is a scalar measure of the strength of the scalp potential field and is calculated as the standard deviation of all electrodes at a given time point (47, 48, 50). The cluster analysis was first computed at the individual level yielding one set of representative maps for each subject. Clustering at the group level followed, in which the individual representative maps of patients and controls were clustered separately. Then all participants’ representative maps were clustered at global level yielding one set of maps that represented the data of all subjects. This one set of global representative maps entered the fitting process and the presence of each global map in every subject was determined. This enabled us to compare groups in terms of the presence of these global maps in the original data, i.e., to compare the microstate temporal characteristics between patients and controls.

In order to retrieve the temporal characteristics of the microstates, the fitting procedure consisted in calculating the spatial correlation between every representative map identified at the global level across all subjects and the individual subject’s scalp potential map in every instant of the pre-processed EEG recording. Each continuous time point of the subject’s EEG (not only the GFP peaks) was then assigned to the microstate class of the highest correlation (winner-takes-all), again ignoring polarity (20, 24, 47, 51). Temporal smoothing parameters [window half size = 3, strength (Besag factor) = 10] ensured that the noise during low GFP did not artificially interrupt the temporal segments of stable topography (18, 47). These segments of stable topography assigned to a given microstate class were then evaluated. For each subject, three temporal parameters were calculated for each microstate class: i) occurrence, ii) coverage, and iii) duration. Occurrence indicates how many temporal segments of a given microstate class occur in 1 s. The coverage in percent represents the summed amount of time spent in a given microstate class as a portion of the whole analyzed period. The duration in milliseconds indicates the most common amount of time that a given microstate class is continuously present. The global explained variance for a specific microstate class was calculated by summing the squared spatial correlations between the representative map and its corresponding assigned scalp potential maps at each time point weighted by the GFP (48). Global explained variances of all microstate classes were then summed yielding the portion of data explained with the set of representative maps.

Microstate analysis was performed using the freely available Cartool Software 3.70, (https://sites.google.com/site/cartoolcommunity/home), programmed by Denis Brunet. Mann-Whitney U test was used to investigate group differences for temporal parameters of each microstate. Multiple comparisons were corrected using the false discovery rate (FDR) method (52).

Spearman’s rank correlations were calculated between the MADRS, YMRS, STAI-state, and STAI-trait scores and significant microstate parameters to check for possible relationships between symptoms and microstate dynamics. Eyes closed spontaneous EEG dynamics are mostly dominated by alpha activity (53). To provide an estimate of the absolute alpha power of each group, individual EEG data were Fourier transformed (Hanning windowing function, common averge reference). Mann-Whitney U test was used to investigate group differences for mean alpha (8–14 Hz) power across 204 channels. To evaluate the role of the alpha rhythm in the appearance of the microstates, the Spearman’s rank correlation was calculated between Alpha Power and those microstate parameters, in which significant group differences were found. All statistical evaluations were performed by the routines included in the program package Statistica’13 (1984–2018, TIBCO, Software Inc, Version 13.4.0.14).

Results

Clinical and Demographic Variables

There were no significant differences in age and level of education between the patient and the control groups. In both groups, very low mean scores on depression and mania symptoms were observed, which did not significantly differ between the two groups. BD patients showed higher scores on state and trait scales of the STAI. For all subject characteristics, see Table 1.

TABLE 1

Table 1 Subject characteristics.

Microstate Results

The meta-criterion used to determine the most dominant topographies revealed five resting-state microstate maps across patients, HCs, and all subjects, explaining 82.1, 83.1, and 82.2% of the global variance, respectively (Figure 1). The topographies resembled those previously reported as A, B, C, and D maps (20, 22, 25, 49) and one of the three recently identified additional maps (23). We labeled these scalp maps from A to E in accordance with the previous literature on microstates. The scalp topographies showed left posterior-right anterior orientation (map A), a right posterior-left anterior orientation (map B), an anterior-posterior orientation (map C), a fronto-central maximum (map D), and a parieto-occipital maximum (map E).

FIGURE 1

Figure 1 The five microstate topographies (A–E) identified in the group clusterings across patients (BD) and healthy controls (HC) and in the global clustering across all subjects (BD+HC). Note that similar five dominant microstate topographies were identified in all three clusterings.

Since some microstate parameters showed a non-homogeneity of variances in the two groups (Levene’s tests for the microstate C coverage and microstates A and C duration; p < 0.01), we decided to calculate Mann-Whitney U test to investigate group differences for temporal parameters of each microstate.

We found significant between-group differences for microstate classes A and B. Both microstates showed increased presence in patients in terms of occurrence and coverage. The two groups did not differ in any temporal parameter of microstates C, D, or E. The results of the temporal characteristics of each microstate are summarized in Table 2 and Figure 2.

TABLE 2

Table 2 Mann-Whitney U test for group comparisons in the investigated microstate parameters.

FIGURE 2

Figure 2 Temporal dynamics of electroencephalographic microstates in patients with bipolar disorder (BD) and in healthy controls (HC). In each subplot, the raw data are plotted on top of the boxplot showing the mean (▪), 95% confidence interval (box plot area), 1 standard deviation (whiskers), and significant differences (*). In all plots, x-axes represent the subject group; y-axes represent the occurrence (upper plots) or coverage (lower plots). Note significantly increased occurrence and coverage of the microstate A and B in the BD compared to HC group (FDR corrected p < 0.05).

Clinical Correlations

Correlations with clinical parameters were calculated for those microstate parameters, in which significant group differences were found. The results of Spearman’s rank correlation revealed a significant positive association between the coverage of the microstate B and the STAI-state (r = 0.40; p < 0.05) and STAI-trait (r = 0.54; p < 0.05) scores. The results of Spearman’s rank correlation revealed a significant positive association between the occurrence of the microstate B and the STAI-trait (r = 0.47; p < 0.05) scores. The results of Spearman’s rank correlation revealed no significant associations between the STAI-state or STAI-trait scores and the occurrence or coverage of the microstate A (all absolute r-values < 0.35).

The results of Spearman’s rank correlation revealed no significant associations between the MADRS and YMRS scores and the occurrence or coverage of the microstate A and B (all absolute r-values < 0.30).

Alpha Rhythm

The Mann-Whitney U test showed significantly decreased alpha power (p < 0.03, Z-value 2.7) in the BD compared to HC group (see Figure 3). The results of Spearman’s rank correlation revealed no significant associations between the alpha power and occurrence or coverage of microstates A and B (all absolute r-values < 0.40).

FIGURE 3

Figure 3 The alpha power in patients with bipolar disorder (BD) and healthy controls (HC). Raw data are plotted on top of each boxplot showing the mean (▪), 95% confidence interval (box plot area), and non-outlier range (whiskers). The x-axis represents the subject group; the y-axis represents the average alpha (8–14 Hz) power across 204 channels. Note significantly decreased alpha power in the BD compared to HC group (p < 0.03, Z value 2.7).

Discussion

Our study presents the first evidence for altered resting-state EEG microstate dynamics in euthymic patients with BD. Patients were stable and did not significantly differ in their depressive or manic symptomatology from HCs at the time of experiment. Despite this fact, they showed abnormally increased presence of microstates A and B, the latter correlating with the anxiety level.

In an earlier combined fMRI-EEG study the microstate A was associated with the auditory network (49). Moreover, generators of the functional EEG microstates were estimated in recent studies, where sources of the microstate A showed left-lateralized activity in the temporal lobe, insula, mPFC, and occipital gyri (23, 24).

In the fMRI literature as well, resting-state functional connectivity alterations of the insula (55), the auditory network (56), and the mPFC (57) were reported in BD patients. Verbal episodic memory deficits and language-related symptoms in BD patients were suggested to be associated with a diminished functional connectivity within the auditory/temporal gyrus and to be compensated by increased fronto-temporal functional connectivity (56). The mPFC was also identified as a major locus of shared abnormality in BD and schizophrenia (58), showing reduced default mode network connectivity from the mPFC to the hippocampus and fusiform gyrus, as well as increased connectivity between the mPFC and primary visual cortex in BD. Hypoconnectivity of the default mode network from the left posterior cingulate cortex to the bilateral mPFC and bilateral precuneus, and reduced salience connectivity of the left sgACC to the right inferior temporal gyrus in BD patients (57) was observed in unmedicated BD patients. In euthymic BD subjects compared to HCs, resting-state functional connectivity of the insula (59) and amygdala (60) to other brain regions was reported to be increased and decreased, respectively. In summary, the evidence from fMRI studies shows both hypoconnectivity (57, 58) and hyperconnectivity (56, 58, 59) pointing to complex alterations of functional resting-state networks. Our findings of increased presence of the microstate A in euthymic BD patients might be related to the hyperconnectivity of the underlying networks that involve the temporal lobe, insula, mPFC, and occipital gyri.

Anxiety symptoms were previously associated with greater severity and impairment in BD (41) and euthymic bipolar patients tend to present high residual level of anxiety (61), as it was the case here. No significant correlation was found between the increased anxiety scores and the increased occurrence or coverage of the microstate A. Our results therefore indicate that this alteration of microstate dynamics might represent a characteristic feature of BD that is not affected by anxiety. The demonstrated alterations in microstate A dynamics during clinical remission might reflect i) an impaired resting-state large-scale brain network dynamics as a trait characteristic of the disorder and/or ii) a compensatory mechanism needed for clinical stabilization of the disorder.

Our study is the first to examine EEG microstate dynamics of spontaneous activity in BD patients during remission. The here demonstrated group difference of BD patients vs. controls, is not congruent with the previously reported reduced duration of the EEG microstates during a depressive episode (34). The experimental group in that study was not restricted to bipolar patients, however, and included also patients who met the criteria for unipolar depression or dysthymia. Moreover, authors examined the overall microstate duration and did not examine distinct microstates separately. These and other aspects, such as different clustering methods used, make it difficult to compare our findings with that early evidence of disrupted microstate dynamics in depression. Nevertheless, in our recent resting-state study we showed positive associations of depressive symptoms with the occurrence of microstate A in a heterogenous group of patients with affective disorders (35).

The microstate B was previously associated with the visual network (23, 24, 49, 62) and thoughts related to the conscious experience of an episodic autobiographic memory, i.e., mental visualization of the scene (24). In our group of BD patients, we found an abnormally increased presence of microstate B that was associated with a higher anxiety. In particular, the occurence together with coverage and only the coverage were positively correlated with scores of trait and state anxiety, respectively. The observed change in microstate B dynamics might be, therefore, more related to a relatively stable disposition than to the actual emotional state. Previous studies also suggest that anxiety may influence visual processing (63, 64) and that connections between amygdala and visual cortex might underlie enhanced visual processing of emotionally salient stimuli in patients with social phobia (65). Our finding of increased presence of microstate B positively associated with anxiety level in euthymic BD patients is consistent with these observations. Additionally, a more regular appearance for microstate B and increased overall temporal dependencies among microstates were recently reported in mood and anxiety disorders, suggesting a decreased dynamicity in switching between different brain states in these psychiatric conditions (66). Another microstate study on anxiety disorders reported a decreased overall resting-state microstate duration in panic disorder (43). That early study, however, did not assess temporal characteristics of different microstates separately and it is therefore difficult to compare those findings with our observations. Further evidence is needed to determine, whether the increased presence of microstate B in our experimental group is a characteristic feature of BD or anxiety, or whether it is related to both conditions.

We found an unchanged duration but a higher occurence and coverage of A and B microstates in BD patients. In other words, an unchanged sustainability in time and still increased presence of these microstates in patients compared to HCs were observed. Possible explanation for this finding appears to be a redundance in activation of the sensory and autobiographic memory networks during spontaneous mentation in euthymic BD patients. Changes in A and B microstate occurence, duration, and coverage have been reported in several psychiatric conditions such as schizophrenia, dementia, narcolepsy, multiple sclerosis, panic disorder, etc. (for review see 20). In schizophrenia patients, both increased (30) and decreased (67, 68) durations of microstate B as well as increased occurrence of microstate A (68) were reported. Increases in duration and occurrence of microstates A and B were also observed in patients with multiple sclerosis, moreover predicting depression scores and other clinical variables (31). It was suggested that multiple sclerosis affects the “sensory” (visual, auditory) rather than the higher-order (salience, central executive) functional networks (20). Our findings of impaired dynamics in microstates A and B suggest a similar interpretation for the BD. Although the pathophysiology in multiple sclerosis differes from that in depression, the increased presence of A and B microstates might reflect aberrant temporal functioning of sensory-related and memory networks in both diseases. Evidence from fMRI studies points to topographical dysbalances between the default mode and sensorimotor networks in BD patients with opposing patterns in depression and mania (69). Cyclothymic and depressive temperaments were associated with opposite changes in the sensorimotor network variability in the resting state signal measured by fractional standard deviation of blood-oxygen-level dependent signal (70). Our findings of altered microstates A and B dynamics is consistent with this fMRI evidence of impaired sensorimotor network in affective disorders, and moreover suggests that neural correlates of these deficits are prominent even during the euthymic state in BD patients.

It is known that EEG microstate dynamics during eyes closed and eyes open resting-state is different. Manipulations of visual input showed increased occurrence and coverage of microstate B and shorter duration of microstate A in the eyes open condition in healthy young adults (71). Further studies are needed to investigate, whether the here observed abnormalities of the A and B microstates in BD patients are still present in the eyes open condition.

BD patients were previously shown to display lower alpha power as compared to HCs (72), as it was the case here. We failed, however, to find any significant correlation between the altered microstate dynamics and decreased alpha power. Our findings, therefore, further support the previously reported independence of microstates from EEG frequency power fluctuations (49).

In summary, results of the current study seem to indicate that dysfunctional activity of resting-state brain networks underlying microstates A and B is a detectable impairment in BD during an euthymic state. The presence of microstate A and B represents measures that might be implicated in clinical practice, although using these parameters for early identification of BD at individual level could prove challenging. If future studies confirm the same pattern in prodromal or vulnerable subjects, it could help detection of at-risk subjects and therefore the possiblility for early intervention. The present study has, however, some limitations. Our low sample size made it impossible to examine any potential influence of medication on the microstate parameters by comparing patients receiving a specific drug with those not receiving it. Possible effects of medication on our results should be therefore taken into account. Due to the same reason, it was not possible to examine any potential influence of subtypes of BD on microstate results.

Conclusions

Our study described altered EEG resting-state microstate temporal parameters in euthymic bipolar patients. Our findings provide an insight into the resting-state global brain network dynamics in BD. Since the increased presence of microstate A is not unique to BD patients, having been reported also in other psychiatric disorders (see 20), it might be considered only as a non-specific electrophysiological marker of BD. Moreover, studies examining possible interactions between microstate dynamics and BD symptoms are needed to better understand the dysfunction of large-scale brain network resting-state dynamics in this affective disorder.

Data Availability Statement

The raw data supporting the conclusions of this manuscript will be made available by the authors, without undue reservation, to any qualified researcher.

Ethics Statement

The studies involving human participants were reviewed and approved by Ethics Committee for Human Research of the Geneva University Hospital. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

AD designed the study, performed the analysis, and wrote the initial draft. J-MA, AD, and CP were responsible for clinical assessment. CM served as an advisor. CB collected the HD-EEG data and was responsible for the overall oversight of the study. All authors revised the manuscript.

Funding

This study was supported by the European Union Horizon 2020 research and innovation program under the Marie Skłodowska-Curie grant agreement No. 739939, by the Swiss National Science Foundation (grant No. 320030_184677), and by the National Centre of Competence in Research (NCCR) “SYNAPSY—The Synaptic Basis of Mental Diseases” (NCCR Synapsy Grant # “51NF40-185897). The funding sources had no role in the design, collection, analysis, or interpretation of the study.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

Special thanks go to Anne-Lise Kung, psychologist, for her involvement in clinical data collection.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpsyt.2019.00826/full#supplementary-material

Supplementary Figure 1 | Microstate analysis: (A) resting-state EEG from subsample of 16 out of 204 electrodes; (B) global field power (GFP) curve with the GFP peaks (vertical lines) in the same EEG period as shown in (A); (C) potential maps at successive GFP peaks, indicated in (B), from the first 1 s period of the recording; (D) set of five cluster maps best explaining the data as revealed by k-means clustering of the maps at the GFP peaks; (E) the original EEG recording shown in (A) with superimposed color-coded microstate segments. Note that each time point of the EEG recording was labelled with the cluster map, shown in (D), with which the instant map correlated best. The duration of segments, occurrence, and coverage for all microstates were computed on thus labeled EEG recording.

References

1. Cloutier M, Greene M, Guerin A, Touya M, Wu E. The economic burden of bipolar I disorder in the United States in 2015. J Affect Disord (2018) 226:45–51. doi: 10.1016/j.jad.2017.09.011

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Eaton WW, Alexandre P, Bienvenu OJ, Clarke D, Martins SS, Nestadt G, et al. The burden of mental disorders. In: Eaton WW, editor. Public Mental Health. Oxford, UK: Oxford University Press (2012). doi: 10.1093/acprof:oso/9780195390445.003.0001.

CrossRef Full Text | Google Scholar

3. Merikangas KR, Jin R, He JP, Kessler RC, Lee S, Sampson NA, et al. Prevalence and correlates of bipolar spectrum disorder in the World Mental health survey initiative. Arch Gen Psychiatry (2011) 68:241–51. doi: 10.1001/archgenpsychiatry.2011.12

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Ferrari AJ, Stockings E, Khoo JP, Erskine HE, Degenhardt L, Vos T, et al. The prevalence and burden of bipolar disorder: findings from the Global Burden of disease study 2013. Bipolar Disord. (2016) 18:440–50. doi: 10.1111/bdi.12423

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Hirschfeld RMA. Screening for bipolar disorder. Am J Managed Care (2007) 13:S164–9.

Google Scholar

6. Vieta E, Salagre E, Grande I, Carvalho AF, Fernandes BS, Berk M, et al. Early intervention in Bipolar disorder. Am J Psychiatry (2018) 175:411–26. doi: 10.1176/appi.ajp.2017.17090972

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Chen CH, Suckling J, Lennox BR, Ooi C, Bullmore ET. A quantitative meta-analysis of fMRI studies in bipolar disorder. Bipolar Disord (2011) 13:1–15. doi: 10.1111/j.1399-5618.2011.00893.x

CrossRef Full Text | Google Scholar

8. Houenou J, Frommberger J, Carde S, Glasbrenner M, Diener C, Leboyer M, et al. Neuroimaging-based markers of bipolar disorder: evidence from two meta-analyses. J Affect Disord (2011) 132:344–55. doi: 10.1016/j.jad.2011.03.016

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Kupferschmidt DA, Zakzanis KK. Toward a functional neuroanatomical signature of bipolar disorder: quantitative evidence from the neuroimaging literature. Psychiatry Res Neuroimaging (2011) 193:71–9. doi: 10.1016/j.pscychresns.2011.02.011

CrossRef Full Text | Google Scholar

10. Piguet C, Fodoulian L, Aubry JM, Vuilleumier P, Houenou J. Bipolar disorder: functional neuroimaging markers in relatives. Neurosci Biobehav. Rev (2015) 57:284–96. doi: 10.1016/j.neubiorev.2015.08.015

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Vargas C, López-Jaramillo C, Vieta E. A systematic literature review of resting state network-functional MRI in bipolar disorder. J Affect Disord (2013) 150:727–35. doi: 10.1016/j.jad.2013.05.083

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Strakowski SM, Adler CM, Almeida J, Altshuler LL, Blumberg HP, Chang KD, et al. The functional neuroanatomy of bipolar disorder: A consensus model. Bipolar Disord (2012) 14:313–25. doi: 10.1111/j.1399-5618.2012.01022.x

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Chase HW, Phillips ML. Elucidating neural network functional connectivity abnormalities in bipolar disorder: toward a harmonized methodological approach. Biol Psychiatry Cogn. Neurosci Neuroimaging (2016) 1:288–98. doi: 10.1016/j.bpsc.2015.12.006

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Rey G, Piguet C, Benders A, Favre S, Eickhoff SB, Aubry JM, et al. Resting-state functional connectivity of emotion regulation networks in euthymic and non-euthymic bipolar disorder patients. Eur Psychiatry (2016) 34:56–63. doi: 10.1016/j.eurpsy.2015.12.005

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Wang Y, Zhong S, Jia Y, Sun Y, Wang B, Liu T, et al. Disrupted resting-state functional connectivity in nonmedicated bipolar disorder. Radiology (2016) 280:529–36. doi: 10.1148/radiol.2016151641

PubMed Abstract | CrossRef Full Text | Google Scholar

16. de Pasquale. F, Corbetta M, Betti. V, Penna Della. Cortical cores in network dynamics. Neuroimage (2018) 180:370–82:S. doi: 10.1016/j.neuroimage.2017.09.063

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Bressler SL, Menon V. Large-scale brain networks in cognition: emerging methods and principles. Trends Cogn. Sci (2010) 14:277–90. doi: 10.1016/j.tics.2010.04.004

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Pascual-Marqui RD, Michel CM, Lehmann D. Segmentation of brain electrical activity into microstates; model estimation and validation. IEEE Trans Biomed Eng. (1995) 42:658–65. doi: 10.1109/10.391164

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Van De Ville D, Britz J, Michel CM. EEG microstate sequences in healthy humans at rest reveal scale-free dynamics. Proc Natl Acad Sci (2010) 107:18179–84. doi: 10.1073/pnas.1007841107

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Michel CM, Koenig T. EEG microstates as a tool for studying the temporal dynamics of whole-brain neuronal networks: a review. Neuroimage (2018) 180:577–93. doi: 10.1016/j.neuroimage.2017.11.062

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Lehmann D, Ozaki H, Pal I. EEG alpha map series: brain micro-states by space-oriented adaptive segmentation. Electroencephalogr. Clin Neurophysiol. (1987) 67:271–88. doi: 10.1016/0013-4694(87)90025-3

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Koenig T, Prichep L, Lehmann D, Sosa PV, Braeker E, Kleinlogel H, et al. Millisecond by millisecond, year by year: normative EEG microstates and developmental stages. Neuroimage (2002) 16:41–8. doi: 10.1006/nimg.2002.1070

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Custo A, Van De Ville D, Wells WM, Tomescu MI, Brunet D, Michel CM. Electroencephalographic resting-state networks: source localization of microstates. Brain Connectivity (2017) 7:671–82. doi: 10.1089/brain.2016.0476

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Bréchet L, Brunet D, Birot G, Gruetter R, Michel CM, Jorge J. Capturing the spatiotemporal dynamics of self-generated, task-initiated thoughts with EEG and fMRI. Neuroimage (2019) 194:82–92. doi: 10.1016/j.neuroimage.2019.03.029

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Khanna A, Pascual-Leone A, Michel CM, Farzan F. Microstates in resting-state EEG: Current status and future directions. Neurosci Biobehav. Rev (2015) 49:105–13. doi: 10.1016/j.neubiorev.2014.12.010

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Changeux JP, Michel CM. Mechanism of neural integration at the brain-scale Level. In: Grillner S, Graybiel AM, editors. Microcircuits. Cambridge, Mass, USA: MIT Press (2004). p. 347–70.

Google Scholar

27. Lehmann D. Brain Electric microstates and cognition: the atoms of thought. In: John ER, Harmony T, Prichep LS, Valdés-Sosa M, Valdés-Sosa PA, editors. Machinery of the Mind. Birkhäuser, Boston, MA, USA (1990). p. 209–24. doi: 10.1007/978-1-4757-1083-0_10

CrossRef Full Text | Google Scholar

28. Tomescu MI, Rihs TA, Becker R, Britz J, Custo A, Grouiller F, et al. Deviant dynamics of EEG resting state pattern in 22q11.2 deletion syndrome adolescents: a vulnerability marker of schizophrenia? Schizophr. Res (2014) 157:175–81. doi: 10.1016/j.schres.2014.05.036

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Tomescu MI, Rihs TA, Roinishvili M, Karahanoglu FI, Schneider M, Menghetti S, et al. Schizophrenia patients and 22q11.2 deletion syndrome adolescents at risk express the same deviant patterns of resting state EEG microstates: a candidate endophenotype of schizophrenia. Schizophr. Res Cogn. (2015) 2:159–65. doi: 10.1016/j.scog.2015.04.005

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Andreou C, Faber PL, Leicht G, Schoettle D, Polomac N, Hanganu-Opatz IL, et al. Resting-state connectivity in the prodromal phase of schizophrenia: Insights from EEG microstates. Schizophr. Res (2014) 152:513–20. doi: 10.1016/j.schres.2013.12.008

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Gschwind M, Hardmeier M, Van De Ville D, Tomescu MI, Penner I, Naegelin Y, et al. Fluctuations of spontaneous EEG topographies predict disease state in relapsing-remitting multiple sclerosis. NeuroImage Clin (2016) 12:466–77. doi: 10.1016/j.nicl.2016.08.008

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Atluri S, Wong W, Moreno S, Blumberger DM, Daskalakis ZJ, Farzan F. Selective modulation of brain network dynamics by seizure therapy in treatment-resistant depression. NeuroImage Clin (2018) 20:1176–90. doi: 10.1016/j.nicl.2018.10.015

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Sverak T, Albrechtova L, Lamos M, Rektorova I, Ustohal L. Intensive repetitive transcranial magnetic stimulation changes EEG microstates in schizophrenia: A pilot study. Schizophr. Res (2018) 193:451–2. doi: 10.1016/j.schres.2017.06.044

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Strik WK, Dierks T, Becker T, Lehmann D. Larger topographical variance and decreased duration of brain electric microstates in depression. J Neural Transm (1995) 99:213–22. doi: 10.1007/BF01271480

CrossRef Full Text | Google Scholar

35. Damborská A, Tomescu MI, Honzírková E, Barteček R, Hořínková J, Fedorová S, et al. EEG resting-state large-scale brain network dynamics are related to depressive symptoms. Front Psychiatry (2019) 548:10. doi: 10.3389/fpsyt.2019.00548

CrossRef Full Text | Google Scholar

36. Piguet C, Dayer A, Aubry JM. Biomarkers and vulnerability to bipolar disorders. Schweiz Arch Neurol Psychiatr (2016) 167:57–67. doi: 10.4414/sanp.2016.00387

CrossRef Full Text | Google Scholar

37. Berchio C, Piguet C, Michel CM, Cordera P, Rihs TA, Dayer AG, et al. Dysfunctional gaze processing in bipolar disorder. NeuroImage Clin (2017) 16:545–56. doi: 10.1016/j.nicl.2017.09.006

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Nurnberger JI Jr., Blehar MC, Kaufmann CA, York-Cooler C, Simpson SG, Harkavy-Friedman J, et al. Diagnostic interview for genetic studies: Rationale, unique features, and training. Arch Gen Psychiatry (1994) 51:849–59. doi: 10.1001/archpsyc.1994.03950110009002

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Young RC, Biggs JT, Ziegler VE, Meyer DA. A rating scale for mania: reliability, validity and sensitivity. Br J Psychiatry (1978) 133:429–35. doi: 10.1192/bjp.133.5.429

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Williams JBW, Kobak KA. Development and reliability of a structured interview guide for the Montgomery Asberg depression rating scale (sigma). Br J Psychiatry (2008) 192:52–8. doi: 10.1192/bjp.bp.106.032532

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Simon NM, Otto MW, Wisniewski SR, Fossey M, Sagduyu K, Frank E, et al. Anxiety disorder comorbidity in bipolar disorder patients: data from the first 500 participants in the Systematic Treatment Enhancement Program for Bipolar Disorder (STEP-BD). Am J Psychiatry (2004) 161:2222–9. doi: 10.1176/appi.ajp.161.12.2222

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Simon NM, Zalta AK, Otto MW, Ostacher MJ, Fischmann D, Chow CW, et al. The association of comorbid anxiety disorders with suicide attempts and suicidal ideation in outpatients with bipolar disorder. J Psychiatr Res (2007) 41:255–64. doi: 10.1016/j.jpsychires.2006.08.004

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Wiedemann G, Stevens A, Pauli P, Dengler W. Decreased duration and altered topography of electroencephalographic microstates in patients with panic disorder. Psychiatry Res Neuroimaging (1998) 84:37–48. doi: 10.1016/S0925-4927(98)00044-4

CrossRef Full Text | Google Scholar

44. Spielberger CD, Gorsuch RL, Lushene RE, (1970). Manual for the State-Trait Anxiety Inventory. Palo Alto, CA, USA: Consulting Psychologists Press.

Google Scholar

45. Jung T, Makeig S, Westerfield M, Townsend J, Courchesne E.c, et al. Removal of eye activity artifacts from visual event-related potentials in normal and clinical subjects. Clin Neurophysiol. (2000) 111:1745–58. doi: 10.1016/S1388-2457(00)00386-2

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Perrin F, Pernier J, Bertrand O, Echallier JF. Spherical splines for scalp potential and current density mapping. Electroencephalogr. Clin Neurophysiol. (1989) 72:184–7. doi: 10.1016/0013-4694(89)90180-6

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Brunet D, Murray M, Michel CM. Spatiotemporal analysis of multichannel EEG: CARTOOL. Comput Intell. Neurosci (2011) 2011, 813870. doi: 10.1155/2011/813870.

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Murray MM, Brunet D, Michel CM. Topographic ERP analyses: a step-by-step tutorial review. Brain Topogr. (2008) 20:249–64. doi: 10.1007/s10548-008-0054-5

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Britz J, Van De Ville D, Michel CM. BOLD correlates of EEG topography reveal rapid resting-state network dynamics. NeuroImage (2010) 52:1162–70. doi: 10.1016/j.neuroimage.2010.02.052

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Michel CM, Brandeis D, Skrandies W, Pascual R, Strik WK, Dierks T, et al. Global field power: a 'time-honoured' index for EEG/EP map analysis. Int J Psychophysiol. (1993) 15:1–2. doi: 10.1016/0167-8760(93)90088-7

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Santarnecchi E, Khanna AR, Musaeus CS, Benwell CSY, Davila P, Farzan F, et al. EEG microstate correlates of fluid intelligence and response to cognitive training. Brain Topogr. (2017) 30:502–20. doi: 10.1007/s10548-017-0565-z

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Benjamini Y. Discovering the false discovery rate. J R Stat Soc Ser B Stat Methodol (2010) 72:405–16. doi: 10.1111/j.1467-9868.2010.00746.x

CrossRef Full Text | Google Scholar

53. Niedermeyer E. Niedermeyer's electroencephalography: basic principles, clinical applications, and related fields. Philadelphia , PA, USA: Lippincott Williams & Wilkins (2011).

Google Scholar

54. Oldfield RC. The assessment and analysis of handedness: The Edinburgh inventory. Neuropsychologia (1971) 9:97–113. doi: 10.1016/0028-3932(71)90067-4

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Yin Z, Chang M, Wei S, Jiang X, Zhou Y, Cui L, et al. Decreased functional connectivity in insular subregions in depressive episodes of bipolar disorder and major depressive disorder. Front Neurosci (2018) 12, 842. doi: 10.3389/fnins.2018.00842.

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Reinke B, de Ven V, Matura S, Linden DEJ, Oertel-Knöchel V. Altered intrinsic functional connectivity in language-related brain regions in association with verbal memory performance in euthymic bipolar patients. Brain Sci (2013) 3:1357–73. doi: 10.3390/brainsci3031357

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Gong J, Chen G, Jia Y, Zhong S, Zhao L, Luo X, et al. Disrupted functional connectivity within the default mode network and salience network in unmedicated bipolar II disorder. Prog Neuro-Psychopharmacol. Biol Psychiatry (2019) 88:11–8. doi: 10.1016/j.pnpbp.2018.06.012

CrossRef Full Text | Google Scholar

58. Öngür D, Lundy M, Greenhouse I, Shinn AK, Menon V, Cohen BM, et al. Default mode network abnormalities in bipolar disorder and schizophrenia. Psychiatry Res Neuroimaging (2010) 183:59–68. doi: 10.1016/j.pscychresns.2010.04.008

CrossRef Full Text | Google Scholar

59. Minuzzi L, Syan SK, Smith M, Hall A, Hall GBC, Frey BN. Structural and functional changes in the somatosensory cortex in euthymic females with bipolar disorder. Aust New Z J Psychiatry (2018) 52:1075–83. doi: 10.1177/0004867417746001

CrossRef Full Text | Google Scholar

60. Li G, Liu P, Andari E, Zhang A, Zhang K. The role of amygdala in patients with euthymic bipolar disorder during resting state. Front Psychiatry (2018) 9, 445. doi: 10.3389/fpsyt.2018.00445

PubMed Abstract | CrossRef Full Text | Google Scholar

61. Albert U, Rosso G, Maina G, Bogetto F. Impact of anxiety disorder comorbidity on quality of life in euthymic bipolar disorder patients: differences between bipolar I and II subtypes. J Affect Disord (2008) 105:297–303. doi: 10.1016/j.jad.2007.05.020

PubMed Abstract | CrossRef Full Text | Google Scholar

62. Custo A, Vulliemoz S, Grouiller F, Van De Ville D, Michel C. EEG source imaging of brain states using spatiotemporal regression. Neuroimage (2014) 96:106–16. doi: 10.1016/j.neuroimage.2014.04.002

PubMed Abstract | CrossRef Full Text | Google Scholar

63. Phelps EA, Ling S, Carrasco M. Emotion facilitates perception and potentiates the perceptual benefits of attention. Psychol Sci (2006) 17:292–9. doi: 10.1111/j.1467-9280.2006.01701.x

PubMed Abstract | CrossRef Full Text | Google Scholar

64. Laretzaki G, Plainis S, Argyropoulos S, Pallikaris IG, Bitsios P. Threat and anxiety affect visual contrast perception. J Psychopharmacol (2010) 24:667–75. doi: 10.1177/0269881108098823

PubMed Abstract | CrossRef Full Text | Google Scholar

65. Goldin PR, Manber T, Hakimi S, Canli T, Gross JJ. Neural bases of social anxiety disorder: Emotional reactivity and cognitive regulation during social and physical threat. Arch Gen Psychiatry (2009) 66:170–80. doi: 10.1001/archgenpsychiatry.2008.525

PubMed Abstract | CrossRef Full Text | Google Scholar

66. Al Zoubi O, Mayeli A, Tsuchiyagaito A, Misaki M, Zotev V, Refai H, et al. EEG microstates temporal dynamics differentiate individuals with mood and anxiety disorders from healthy subjects. Front Hum Neurosci (2019) 13, 56. doi: 10.3389/fnhum.2019.00056

PubMed Abstract | CrossRef Full Text | Google Scholar

67. Nishida K, Morishima Y, Yoshimura M, Isotani T, Irisawa S, Jann K, et al. EEG microstates associated with salience and frontoparietal networks in frontotemporal dementia, schizophrenia and Alzheimer's disease. Clin Neurophysiol. (2013) 124:1106–14. doi: 10.1016/j.clinph.2013.01.005

PubMed Abstract | CrossRef Full Text | Google Scholar

68. Lehmann D, Faber PL, Galderisi S, Herrmann WM, Kinoshita T, Koukkou M, et al. EEG microstate duration and syntax in acute, medication-naïve, first-episode schizophrenia: a multi-center study. Psychiatry Res Neuroimaging (2005) 138:141–56. doi: 10.1016/j.pscychresns.2004.05.007

CrossRef Full Text | Google Scholar

69. Martino M, Magioncalda P, Huang Z, Conio B, Piaggio N, Duncan NW, et al. Contrasting variability patterns in the default mode and sensorimotor networks balance in bipolar depression and mania. Proc Natl Acad Sci USA (2016) 113:4824–9. doi: 10.1073/pnas.1517558113

PubMed Abstract | CrossRef Full Text | Google Scholar

70. Conio B, Magioncalda P, Martino M, Tumati S, Capobianco L, Escelsior A, et al. Opposing patterns of neuronal variability in the sensorimotor network mediate cyclothymic and depressive temperaments. Hum Brain Mapp. (2019) 40:1344–52. doi: 10.1002/hbm.24453

PubMed Abstract | CrossRef Full Text | Google Scholar

71. Seitzman BA, Abell M, Bartley SC, Erickson MA, Bolbecker AR, Hetrick WP. Cognitive manipulation of brain electric microstates. Neuroimage (2017) 146:533–43. doi: 10.1016/j.neuroimage.2016.10.002

PubMed Abstract | CrossRef Full Text | Google Scholar

72. Basar E, Güntekin B, Atagün I, Turp Gölbas B, Tülay E, Özerdem A. Brain's alpha activity is highly reduced in euthymic bipolar disorder patients. Cogn. Neurodynamics (2012) 6:11–20. doi: 10.1007/s11571-011-9172-y

CrossRef Full Text | Google Scholar