Neuronal synchrony and critical bistability: Mechanistic biomarkers for localizing the epileptogenic network

OBJECTIVE

Postsurgical seizure freedom in drug-resistant epilepsy (DRE) patients varies from 30% to 80%, implying that in many cases the current approaches fail to fully map the epileptogenic zone (EZ). We aimed to advance a novel approach to better characterize epileptogenicity and investigate whether the EZ encompasses a broader epileptogenic network (EpiNet) beyond the seizure zone (SZ) that exhibits seizure activity.

METHODS

We first used computational modeling to test putative complex systems-driven and systems neuroscience-driven mechanistic biomarkers for epileptogenicity. We then used these biomarkers to extract features from resting-state stereoelectroencephalograms recorded from DRE patients and trained supervised classifiers to localize the SZ against gold standard clinical localization. To further explore the prevalence of pathological features in an extended brain network outside of the clinically identified SZ, we also used unsupervised classification.

RESULTS

Supervised SZ classification trained on individual features achieved accuracies of .6-.7 area under the receiver operating characteristic curve (AUC). Combining all criticality and synchrony features further improved the AUC to .85. Unsupervised classification discovered an EpiNet-like cluster of brain regions, in which 51% of brain regions were outside of the SZ. Brain regions in the EpiNet-like cluster engaged in interareal hypersynchrony and locally exhibited high-amplitude bistability and excessive inhibition, which was strikingly similar to the high seizure risk regime revealed by our computational modeling.

SIGNIFICANCE

The finding that combining biomarkers improves SZ localization accuracy indicates that the novel mechanistic biomarkers for epileptogenicity employed here yield synergistic information. On the other hand, the discovery of SZ-like brain dynamics outside of the clinically defined SZ provides empirical evidence of an extended pathophysiological EpiNet.

Rhythmicity of neuronal oscillations delineates their cortical and spectral architecture

Neuronal oscillations are commonly analyzed with power spectral methods that quantify signal amplitude, but not rhythmicity or 'oscillatoriness' per se. Here we introduce a new approach, the phase-autocorrelation function (pACF), for the direct quantification of rhythmicity. We applied pACF to human intracerebral stereoelectroencephalography (SEEG) and magnetoencephalography (MEG) data and uncovered a spectrally and anatomically fine-grained cortical architecture in the rhythmicity of single- and multi-frequency neuronal oscillations. Evidencing the functional significance of rhythmicity, we found it to be a prerequisite for long-range synchronization in resting-state networks and to be dynamically modulated during event-related processing. We also extended the pACF approach to measure 'burstiness' of oscillatory processes and characterized regions with stable and bursty oscillations. These findings show that rhythmicity is double-dissociable from amplitude and constitutes a functionally relevant and dynamic characteristic of neuronal oscillations.

Measuring digital intervention user experience with a novel ecological momentary assessment (EMA) method, CORTO

Digital interventions often suffer from low usage, which may reflect insufficient attention to user experience. Moreover, the existing evaluation methods have limited applicability in the remote study of user experience of complex interventions that have expansive content and that are used over an extensive period of time. To alleviate these challenges, we describe here a novel qualitative Ecological Momentary Assessment (EMA) method: the CORTO method (Contextual, One-item, Repeated, Timely, Open-ended). We used it to gather digital intervention user experience data from Finnish adults (n = 184) who lived with interview-confirmed major depressive disorder (MDD) and took part in a randomized controlled trial (RCT) that studied the efficacy of a novel 12-week game-based digital intervention for depression. A second dataset on user experience was gathered with retrospective interviews (n = 22). We inductively coded the CORTO method and retrospective interview data, which led to four user experience categories: (1) contextual use, (2) interaction-elicited emotional experience, (3) usability, and (4) technical issues. Then, we used the created user experience categories and Template Analysis to analyze both datasets together, and reported the results qualitatively. Finally, we compared the two datasets with each other. We found that the data generated with the CORTO method offered more insights into usability and technical categories than the interview data that particularly illustrated the contextual use. The emotional valence of the interview data was more positive compared with the CORTO data. Both the CORTO and interview data detected 55 % of the micro-level categories; 20 % of micro-level categories were only detected by the CORTO data and 25 % only by the interview data. We found that the during-intervention user experience measurement with the CORTO method can provide intervention-specific insights, and thereby further the iterative user-centered intervention development. Overall, these findings highlight the impact of evaluation methods on the categories and qualities of insights acquired in intervention research.

Infra-slow fluctuations in cortical potentials and respiration drive fast cortical EEG rhythms in sleeping and waking states

OBJECTIVE

Infra-slow fluctuations (ISF, 0.008-0.1 Hz) characterize hemodynamic and electric potential signals of human brain. ISFs correlate with the amplitude dynamics of fast (>1 Hz) neuronal oscillations, and may arise from permeability fluctuations of the blood-brain barrier (BBB). It is unclear if physiological rhythms like respiration drive or track fast cortical oscillations, and the role of sleep in this coupling is unknown.

METHODS

We used high-density full-band electroencephalography (EEG) in healthy human volunteers (N = 21) to measure concurrently the ISFs, respiratory pulsations, and fast neuronal oscillations during periods of wakefulness and sleep, and to assess the strength and direction of their phase-amplitude coupling.

RESULTS

The phases of ISFs and respiration were both coupled with the amplitude of fast neuronal oscillations, with stronger ISF coupling being evident during sleep. Phases of ISF and respiration drove the amplitude dynamics of fast oscillations in sleeping and waking states, with different contributions.

CONCLUSIONS

ISFs in slow cortical potentials and respiration together significantly determine the dynamics of fast cortical oscillations.

SIGNIFICANCE

We propose that these slow physiological phases play a significant role in coordinating cortical excitability, which is a fundamental aspect of brain function.

Critical-like Brain Dynamics in a Continuum from Second- to First-Order Phase Transition

The classic brain criticality hypothesis postulates that the brain benefits from operating near a continuous second-order phase transition. Slow feedback regulation of neuronal activity could, however, lead to a discontinuous first-order transition and thereby bistable activity. Observations of bistability in awake brain activity have nonetheless remained scarce and its functional significance unclear. Moreover, there is no empirical evidence to support the hypothesis that the human brain could flexibly operate near either a first- or second-order phase transition despite such a continuum being common in models. Here, using computational modeling, we found bistable synchronization dynamics to emerge through elevated positive feedback and occur exclusively in a regimen of critical-like dynamics. We then assessed bistability in vivo with resting-state MEG in healthy adults (7 females, 11 males) and stereo-electroencephalography in epilepsy patients (28 females, 36 males). This analysis revealed that a large fraction of the neocortices exhibited varying degrees of bistability in neuronal oscillations from 3 to 200 Hz. In line with our modeling results, the neuronal bistability was positively correlated with classic assessment of brain criticality across narrow-band frequencies. Excessive bistability was predictive of epileptic pathophysiology in the patients, whereas moderate bistability was positively correlated with task performance in the healthy subjects. These empirical findings thus reveal the human brain as a one-of-a-kind complex system that exhibits critical-like dynamics in a continuum between continuous and discontinuous phase transitions.SIGNIFICANCE STATEMENT In the model, while synchrony per se was controlled by connectivity, increasing positive local feedback led to gradually emerging bistable synchrony with scale-free dynamics, suggesting a continuum between second- and first-order phase transitions in synchrony dynamics inside a critical-like regimen. In resting-state MEG and SEEG, bistability of ongoing neuronal oscillations was pervasive across brain areas and frequency bands and was observed only with concurring critical-like dynamics as the modeling predicted. As evidence for functional relevance, moderate bistability was positively correlated with executive functioning in the healthy subjects, and excessive bistability was associated with epileptic pathophysiology. These findings show that critical-like neuronal dynamics in vivo involves both continuous and discontinuous phase transitions in a frequency-, neuroanatomy-, and state-dependent manner.

The influence of inter-regional delays in generating large-scale brain networks of phase synchronization

Large-scale networks of phase synchronization are considered to regulate the communication between brain regions fundamental to cognitive function, but the mapping to their structural substrates, i.e., the structure-function relationship, remains poorly understood. Biophysical Network Models (BNMs) have demonstrated the influences of local oscillatory activity and inter-regional anatomical connections in generating alpha-band (8-12 Hz) networks of phase synchronization observed with Electroencephalography (EEG) and Magnetoencephalography (MEG). Yet, the influence of inter-regional conduction delays remains unknown. In this study, we compared a BNM with standard "distance-dependent delays", which assumes constant conduction velocity, to BNMs with delays specified by two alternative methods accounting for spatially varying conduction velocities, "isochronous delays" and "mixed delays". We followed the Approximate Bayesian Computation (ABC) workflow, i) specifying neurophysiologically informed prior distributions of BNM parameters, ii) verifying the suitability of the prior distributions with Prior Predictive Checks, iii) fitting each of the three BNMs to alpha-band MEG resting-state data (N = 75) with Bayesian optimization for Likelihood-Free Inference (BOLFI), and iv) choosing between the fitted BNMs with ABC model comparison on a separate MEG dataset (N = 30). Prior Predictive Checks revealed the range of dynamics generated by each of the BNMs to encompass those seen in the MEG data, suggesting the suitability of the prior distributions. Fitting the models to MEG data yielded reliable posterior distributions of the parameters of each of the BNMs. Finally, model comparison revealed the BNM with "distance-dependent delays", as the most probable to describe the generation of alpha-band networks of phase synchronization seen in MEG. These findings suggest that distance-dependent delays might contribute to the neocortical architecture of human alpha-band networks of phase synchronization. Hence, our study illuminates the role of inter-regional delays in generating the large-scale networks of phase synchronization that might subserve the communication between regions vital to cognition.

The Development of Game-Based Digital Mental Health Interventions: Bridging the Paradigms of Health Care and Entertainment

Game elements are increasingly used to improve user engagement in digital mental health interventions, and specific game mechanics may yield therapeutic effects per se and thereby contribute to digital mental health intervention efficacy. However, only a few commercial game-based interventions are available. We suggest that the key challenge in their development reflects the tension between the 2 underlying paradigms, health care and entertainment, which have disparate goals and processes in digital development. We describe 3 approaches currently used to negotiate the 2 paradigms: the gamification of health care software, designing serious games, and purpose shifting existing entertainment games. We advanced an integrative framework to focus attention on 4 key themes in intervention development: target audience, engagement, mechanisms of action, and health-related effectiveness. On each theme, we show how the 2 paradigms contrast and can complement each other. Finally, we consider the 4 interdependent themes through the new product development phases from concept to production. Our viewpoint provides an integrative synthesis that facilitates the research, design, and development of game-based digital mental health interventions.

Brain criticality predicts individual levels of inter-areal synchronization in human electrophysiological data

Neuronal oscillations and their synchronization between brain areas are fundamental for healthy brain function. Yet, synchronization levels exhibit large inter-individual variability that is associated with behavioral variability. We test whether individual synchronization levels are predicted by individual brain states along an extended regime of critical-like dynamics - the Griffiths phase (GP). We use computational modelling to assess how synchronization is dependent on brain criticality indexed by long-range temporal correlations (LRTCs). We analyze LRTCs and synchronization of oscillations from resting-state magnetoencephalography and stereo-electroencephalography data. Synchronization and LRTCs are both positively linearly and quadratically correlated among healthy subjects, while in epileptogenic areas they are negatively linearly correlated. These results show that variability in synchronization levels is explained by the individual position along the GP with healthy brain areas operating in its subcritical and epileptogenic areas in its supercritical side. We suggest that the GP is fundamental for brain function allowing individual variability while retaining functional advantages of criticality.

Factors Affecting Digital Tool Use in Client Interaction According to Mental Health Professionals: Interview Study

BACKGROUND

Digital tools and interventions are being increasingly developed in response to the growing mental health crisis, and mental health professionals (MHPs) considerably influence their adoption in client practice. However, how MHPs use digital tools in client interaction is yet to be sufficiently understood, which poses challenges to their design, development, and implementation.

OBJECTIVE

This study aimed to create a contextual understanding of how MHPs use different digital tools in clinical client practice and what characterizes the use across tools.

METHODS

A total of 19 Finnish MHPs participated in semistructured interviews, and the data were transcribed, coded, and inductively analyzed.

RESULTS

We found that MHP digital tool use was characterized by 3 distinct functions: communication, diagnosis and evaluation, and facilitating therapeutic change. The functions were addressed using analog tools, digitized tools that mimic their analog counterparts, and digital tools that use the possibilities native to digital. The MHP-client communication included various media alongside face-to-face meetings, the MHPs increasingly used digitized tools in client evaluation, and the MHPs actively used digitized materials to facilitate therapeutic change. MHP tool use was generally characterized by adaptability-it was negotiated in client interactions. However, there was considerable variance in the breadth of MHPs' digital toolbox. The existing clinical practices emphasized MHP-client interaction and invited incremental rather than radical developments, which challenged the achievement of the scalability benefits expected from digital tools.

CONCLUSIONS

MHPs use digitized and digital tools in client practice. Our results contribute to the user-centered research, development, and implementation of new digital solutions in mental health care by classifying them according to their function and medium and describing how MHPs use and do not use them.

Genetic polymorphisms in COMT and BDNF influence synchronization dynamics of human neuronal oscillations

Neuronal oscillations, their inter-areal synchronization, and scale-free dynamics constitute fundamental mechanisms for cognition by regulating communication in neuronal networks. These oscillatory dynamics have large inter-individual variability that is partly heritable. We hypothesized that this variability could be partially explained by genetic polymorphisms in neuromodulatory genes. We recorded resting-state magnetoencephalography (MEG) from 82 healthy participants and investigated whether oscillation dynamics were influenced by genetic polymorphisms in catechol-O-methyltransferase (COMT) Val¹⁵⁸Met and brain-derived neurotrophic factor (BDNF) Val⁶⁶Met. Both COMT and BDNF polymorphisms influenced local oscillation amplitudes and their long-range temporal correlations (LRTCs), while only BDNF polymorphism affected the strength of large-scale synchronization. Our findings demonstrate that COMT and BDNF genetic polymorphisms contribute to inter-individual variability in neuronal oscillation dynamics. Comparison of these results to computational modeling of near-critical synchronization dynamics further suggested that COMT and BDNF polymorphisms influenced local oscillations by modulating the excitation-inhibition balance according to the brain criticality framework.

•

Human local oscillation dynamics is influenced by polymorphisms in COMT and BNDF

•

COMT and BDNF influence oscillation amplitudes and long-range temporal correlations

•

BDNF polymorphism affected the strength of large-scale synchronization

•

Framework of brain criticality links COMT and BDNF with local E/I-balance

Phase-Based Cortical Synchrony Is Affected by Prematurity

Inter-areal synchronization by phase–phase correlations (PPCs) of cortical oscillations mediates many higher neurocognitive functions, which are often affected by prematurity, a globally prominent neurodevelopmental risk factor. Here, we used electroencephalography to examine brain-wide cortical PPC networks at term-equivalent age, comparing human infants after early prematurity to a cohort of healthy controls. We found that prematurity affected these networks in a sleep state-specific manner, and the differences between groups were also frequency-selective, involving brain-wide connections. The strength of synchronization in these networks was predictive of clinical outcomes in the preterm infants. These findings show that prematurity affects PPC networks in a clinically significant manner, suggesting early functional biomarkers of later neurodevelopmental compromise that may be used in clinical or translational studies after early neonatal adversity.

Long range temporal correlations (LRTCs) in MEG-data during emerging psychosis: Relationship to symptoms, medication-status and clinical trajectory

Long-Range Temporal Correlations (LRTCs) index the capacity of the brain to optimally process information. Previous research has shown that patients with chronic schizophrenia present altered LRTCs at alpha and beta oscillations. However, it is currently unclear at which stage of schizophrenia aberrant LRTCs emerge. To address this question, we investigated LRTCs in resting-state magnetoencephalographic (MEG) recordings obtained from patients with affective disorders and substance abuse (clinically at low-risk of psychosis, CHR-N), patients at clinical high-risk of psychosis (CHR-P) (n = 115), as well as patients with a first episode (FEP) (n = 25). Matched healthy controls (n = 47) served as comparison group. LRTCs were obtained for frequencies from 4 to 40 Hz and correlated with clinical and neuropsychological data. In addition, we examined the relationship between LRTCs and transition to psychosis in CHR-P participants, and the relationship between LRTC and antipsychotic medication in FEP participants. Our results show that participants from the clinical groups have similar LRTCs to controls. In addition, LRTCs did not correlate with clinical and neurocognitive variables across participants nor did LRTCs predict transition to psychosis. Therefore, impaired LRTCs do not reflect a feature in the clinical trajectory of psychosis. Nevertheless, reduced LRTCs in the beta-band over posterior sensors of medicated FEP participants indicate that altered LRTCs may appear at the onset of the illness. Future studies are needed to elucidate the role of anti-psychotic medication in altered LRTCs.

Detrended fluctuation analysis in the presurgical evaluation of parietal lobe epilepsy patients

OBJECTIVE

To examine the usability of long-range temporal correlations (LRTCs) in non-invasive localization of the epileptogenic zone (EZ) in refractory parietal lobe epilepsy (RPLE) patients.

METHODS

We analyzed 10 RPLE patients who had presurgical MEG and underwent epilepsy surgery. We quantified LRTCs with detrended fluctuation analysis (DFA) at four frequency bands for 200 cortical regions estimated using individual source models. We correlated individually the DFA maps to the distance from the resection area and from cortical locations of interictal epileptiform discharges (IEDs). Additionally, three clinical experts inspected the DFA maps to visually assess the most likely EZ locations.

RESULTS

The DFA maps correlated with the distance to resection area in patients with type II focal cortical dysplasia (FCD) (p<0.05), but not in other etiologies. Similarly, the DFA maps correlated with the IED locations only in the FCD II patients. Visual analysis of the DFA maps showed high interobserver agreement and accuracy in FCD patients in assigning the affected hemisphere and lobe.

CONCLUSIONS

Aberrant LRTCs correlate with the resection areas and IED locations.

SIGNIFICANCE

This methodological pilot study demonstrates the feasibility of approximating cortical LRTCs from MEG that may aid in the EZ localization and provide new non-invasive insight into the presurgical evaluation of epilepsy.

Measurement of excitation-inhibition ratio in autism spectrum disorder using critical brain dynamics

Balance between excitation (E) and inhibition (I) is a key principle for neuronal network organization and information processing. Consistent with this notion, excitation-inhibition imbalances are considered a pathophysiological mechanism in many brain disorders including autism spectrum disorder (ASD). However, methods to measure E/I ratios in human brain networks are lacking. Here, we present a method to quantify a functional E/I ratio (fE/I) from neuronal oscillations, and validate it in healthy subjects and children with ASD. We define structural E/I ratio in an in silico neuronal network, investigate how it relates to power and long-range temporal correlations (LRTC) of the network's activity, and use these relationships to design the fE/I algorithm. Application of this algorithm to the EEGs of healthy adults showed that fE/I is balanced at the population level and is decreased through GABAergic enforcement. In children with ASD, we observed larger fE/I variability and stronger LRTC compared to typically developing children (TDC). Interestingly, visual grading for EEG abnormalities that are thought to reflect E/I imbalances revealed elevated fE/I and LRTC in ASD children with normal EEG compared to TDC or ASD with abnormal EEG. We speculate that our approach will help understand physiological heterogeneity also in other brain disorders.

Spectral and Anatomical Patterns of Large-Scale Synchronization Predict Human Attentional Capacity

The capacity of visual attention determines how many visual objects may be perceived at any moment. This capacity can be investigated with multiple object tracking (MOT) tasks, which have shown that it varies greatly between individuals. The neuronal mechanisms underlying capacity limits have remained poorly understood. Phase synchronization of cortical oscillations coordinates neuronal communication within the fronto-parietal attention network and between the visual regions during endogenous visual attention. We tested a hypothesis that attentional capacity is predicted by the strength of pretarget synchronization within attention-related cortical regions. We recorded cortical activity with magneto- and electroencephalography (M/EEG) while measuring attentional capacity with MOT tasks and identified large-scale synchronized networks from source-reconstructed M/EEG data. Individual attentional capacity was correlated with load-dependent strengthening of theta (3–8 Hz), alpha (8–10 Hz), and gamma-band (30–120 Hz) synchronization that connected the visual cortex with posterior parietal and prefrontal cortices. Individual memory capacity was also preceded by crossfrequency phase–phase and phase–amplitude coupling of alpha oscillation phase with beta and gamma oscillations. Our results show that good attentional capacity is preceded by efficient dynamic functional coupling and decoupling within brain regions and across frequencies, which may enable efficient communication and routing of information between sensory and attentional systems.

Comparison of Methods to Identify Modules in Noisy or Incomplete Brain Networks

Community structure, or "modularity," is a fundamentally important aspect in the organization of structural and functional brain networks, but their identification with community detection methods is confounded by noisy or missing connections. Although several methods have been used to account for missing data, the performance of these methods has not been compared quantitatively so far. In this study, we compared four different approaches to account for missing connections when identifying modules in binary and weighted networks using both Louvain and Infomap community detection algorithms. The four methods are "zeros," "row-column mean," "common neighbors," and "consensus clustering." Using Lancichinetti-Fortunato-Radicchi benchmark-simulated binary and weighted networks, we find that "zeros," "row-column mean," and "common neighbors" approaches perform well with both Louvain and Infomap, whereas "consensus clustering" performs well with Louvain but not Infomap. A similar pattern of results was observed with empirical networks from stereotactical electroencephalography data, except that "consensus clustering" outperforms other approaches on weighted networks with Louvain. Based on these results, we recommend any of the four methods when using Louvain on binary networks, whereas "consensus clustering" is superior with Louvain clustering of weighted networks. When using Infomap, "zeros" or "common neighbors" should be used for both binary and weighted networks. These findings provide a basis to accounting for noisy or missing connections when identifying modules in brain networks.

Dynamic large-scale network synchronization from perception to action

Sensory-guided actions entail the processing of sensory information, generation of perceptual decisions, and the generation of appropriate actions. Neuronal activity underlying these processes is distributed into sensory, fronto-parietal, and motor brain areas, respectively. How the neuronal processing is coordinated across these brain areas to support functions from perception to action remains unknown. We investigated whether phase synchronization in large-scale networks coordinate these processes. We recorded human cortical activity with magnetoencephalography (MEG) during a task in which weak somatosensory stimuli remained unperceived or were perceived. We then assessed dynamic evolution of phase synchronization in large-scale networks from source-reconstructed MEG data by using advanced analysis approaches combined with graph theory. Here we show that perceiving and reporting of weak somatosensory stimuli is correlated with sustained strengthening of large-scale synchrony concurrently in delta/theta (3-7 Hz) and gamma (40-60 Hz) frequency bands. In a data-driven network localization, we found this synchronization to dynamically connect the task-relevant, that is, the fronto-parietal, sensory, and motor systems. The strength and temporal pattern of interareal synchronization were also correlated with the response times. These data thus show that key brain areas underlying perception, decision-making, and actions are transiently connected by large-scale dynamic phase synchronization in the delta/theta and gamma bands.

Phase-lags in large scale brain synchronization: Methodological considerations and in-silico analysis

Architecture of phase relationships among neural oscillations is central for their functional significance but has remained theoretically poorly understood. We use phenomenological model of delay-coupled oscillators with increasing degree of topological complexity to identify underlying principles by which the spatio-temporal structure of the brain governs the phase lags between oscillatory activity at distant regions. Phase relations and their regions of stability are derived and numerically confirmed for two oscillators and for networks with randomly distributed or clustered bimodal delays, as a first approximation for the brain structural connectivity. Besides in-phase, clustered delays can induce anti-phase synchronization for certain frequencies, while the sign of the lags is determined by the natural frequencies and by the inhomogeneous network interactions. For in-phase synchronization faster oscillators always phase lead, while stronger connected nodes lag behind the weaker during frequency depression, which consistently arises for in-silico results. If nodes are in anti-phase regime, then a distance π is added to the in-phase trends. The statistics of the phases is calculated from the phase locking values (PLV), as in many empirical studies, and we scrutinize the method’s impact. The choice of surrogates do not affects the mean of the observed phase lags, but higher significance levels that are generated by some surrogates, cause decreased variance and might fail to detect the generally weaker coherence of the interhemispheric links. These links are also affected by the non-stationary and intermittent synchronization, which causes multimodal phase lags that can be misleading if averaged. Taken together, the results describe quantitatively the impact of the spatio-temporal connectivity of the brain to the synchronization patterns between brain regions, and to uncover mechanisms through which the spatio-temporal structure of the brain renders phases to be distributed around 0 and π.

Trial registration: South African Clinical Trials Register: http://www.sanctr.gov.za/SAClinicalbrnbspTrials/tabid/169/Default.aspx, then link to respiratory tract then link to tuberculosis, pulmonary; and TASK Applied Sciences Clinical Trials, AP-TB-201-16 (ALOPEXX): https://task.org.za/clinical-trials/.

Functional connectivity, and in particular, phase coupling between distant brain regions may be fundamental in regulating neuronal processing and communication. However, phase relationships between the nodes of the brain and how they are confined by its spatio-temporal structure, have been mostly overlooked. We use a model of oscillatory dynamics superimposed on the space-time structure defined by the connectome, and we analyze the possible regimes of synchronization. Limitations of data analysis are also considered and we show that the choice of the significance threshold for coherence does not essentially impact the statistics of the observed phase lags, although it is crucial for the right detection of statistically significant coherence. Analytical insights are obtained for networks with heterogeneous time-delays, based on the empirical data from the connectome, and these are confirmed by numerical simulations, which show in- or anti-phase synchronization depending on the frequency and the distribution of time-delays. Phase lags are shown to result from inhomogeneous network interactions, so that stronger connected nodes generally phase lag behind the weaker.

Hyperedge bundling: Data, source code, and precautions to modeling-accuracy bias to synchrony estimates

It has not been well documented that MEG/EEG functional connectivity graphs estimated with zero-lag-free interaction metrics are severely confounded by a multitude of spurious interactions (SI), i.e., the false-positive "ghosts" of true interactions [1], [2]. These SI are caused by the multivariate linear mixing between sources, and thus they pose a severe challenge to the validity of connectivity analysis. Due to the complex nature of signal mixing and the SI problem, there is a need to intuitively demonstrate how the SI are discovered and how they can be attenuated using a novel approach that we termed hyperedge bundling. Here we provide a dataset with software with which the readers can perform simulations in order to better understand the theory and the solution to SI. We include the supplementary material of [1] that is not directly relevant to the hyperedge bundling per se but reflects important properties of the MEG source model and the functional connectivity graphs. For example, the gyri of dorsal-lateral cortices are the most accurately modeled areas; the sulci of inferior temporal, frontal and the insula have the least modeling accuracy. Importantly, we found the interaction estimates are heavily biased by the modeling accuracy between regions, which means the estimates cannot be straightforwardly interpreted as the coupling between brain regions. This raise a red flag that the conventional method of thresholding graphs by estimate values is rather suboptimal: because the measured topology of the graph reflects the geometric property of source-model instead of the cortical interactions under investigation.

Ghost interactions in MEG/EEG source space: A note of caution on inter-areal coupling measures

When combined with source modeling, magneto- (MEG) and electroencephalography (EEG) can be used to study long-range interactions among cortical processes non-invasively. Estimation of such inter-areal connectivity is nevertheless hindered by instantaneous field spread and volume conduction, which artificially introduce linear correlations and impair source separability in cortical current estimates. To overcome the inflating effects of linear source mixing inherent to standard interaction measures, alternative phase- and amplitude-correlation based connectivity measures, such as imaginary coherence and orthogonalized amplitude correlation have been proposed. Being by definition insensitive to zero-lag correlations, these techniques have become increasingly popular in the identification of correlations that cannot be attributed to field spread or volume conduction. We show here, however, that while these measures are immune to the direct effects of linear mixing, they may still reveal large numbers of spurious false positive connections through field spread in the vicinity of true interactions. This fundamental problem affects both region-of-interest-based analyses and all-to-all connectome mappings. Most importantly, beyond defining and illustrating the problem of spurious, or "ghost" interactions, we provide a rigorous quantification of this effect through extensive simulations. Additionally, we further show that signal mixing also significantly limits the separability of neuronal phase and amplitude correlations. We conclude that spurious correlations must be carefully considered in connectivity analyses in MEG/EEG source space even when using measures that are immune to zero-lag correlations.

Fluctuations of the EEG-fMRI correlation reflect intrinsic strength of functional connectivity in default mode network

Both functional magnetic resonance imaging (fMRI) and electrophysiological recordings have revealed that resting-state functional connectivity is temporally variable in human brain. Combined full-band electroencephalography-fMRI (fbEEG-fMRI) studies have shown that infraslow (<.1 Hz) fluctuations in EEG scalp potential are correlated with the blood-oxygen-level-dependent (BOLD) fMRI signals and that also this correlation appears variable over time. Here, we used simultaneous fbEEG-fMRI to test the hypothesis that correlation dynamics between BOLD and fbEEG signals could be explained by fluctuations in the activation properties of resting-state networks (RSNs) such as the extent or strength of their activation. We used ultrafast magnetic resonance encephalography (MREG) fMRI to enable temporally accurate and statistically robust short-time-window comparisons of infra-slow fbEEG and BOLD signals. We found that the temporal fluctuations in the fbEEG-BOLD correlation were dependent on RSN connectivity strength, but not on the mean signal level or magnitude of RSN activation or motion during scanning. Moreover, the EEG-fMRI correlations were strongest when the intrinsic RSN connectivity was strong and close to the pial surface. Conversely, weak fbEEG-BOLD correlations were attributable to periods of less coherent or spatially more scattered intrinsic RSN connectivity, or RSN activation in deeper cerebral structures. The results thus show that the on-average low correlations between infra-slow EEG and BOLD signals are, in fact, governed by the momentary coherence and depth of the underlying RSN activation, and may reach systematically high values with appropriate source activities. These findings further consolidate the notion of slow scalp potentials being directly coupled to hemodynamic fluctuations.

Hyperedge bundling: A practical solution to spurious interactions in MEG/EEG source connectivity analyses

Inter-areal functional connectivity (FC), neuronal synchronization in particular, is thought to constitute a key systems-level mechanism for coordination of neuronal processing and communication between brain regions. Evidence to support this hypothesis has been gained largely using invasive electrophysiological approaches. In humans, neuronal activity can be non-invasively recorded only with magneto- and electroencephalography (MEG/EEG), which have been used to assess FC networks with high temporal resolution and whole-scalp coverage. However, even in source-reconstructed MEG/EEG data, signal mixing, or "source leakage", is a significant confounder for FC analyses and network localization. Signal mixing leads to two distinct kinds of false-positive observations: artificial interactions (AI) caused directly by mixing and spurious interactions (SI) arising indirectly from the spread of signals from true interacting sources to nearby false loci. To date, several interaction metrics have been developed to solve the AI problem, but the SI problem has remained largely intractable in MEG/EEG all-to-all source connectivity studies. Here, we advance a novel approach for correcting SIs in FC analyses using source-reconstructed MEG/EEG data. Our approach is to bundle observed FC connections into hyperedges by their adjacency in signal mixing. Using realistic simulations, we show here that bundling yields hyperedges with good separability of true positives and little loss in the true positive rate. Hyperedge bundling thus significantly decreases graph noise by minimizing the false-positive to true-positive ratio. Finally, we demonstrate the advantage of edge bundling in the visualization of large-scale cortical networks with real MEG data. We propose that hypergraphs yielded by bundling represent well the set of true cortical interactions that are detectable and dissociable in MEG/EEG connectivity analysis.

Functional integration across oscillation frequencies by cross-frequency phase synchronization

Neuronal oscillations and their inter-areal synchronization may be instrumental in regulating neuronal communication in distributed networks. Several lines of research have, however, shown that cognitive tasks engage neuronal oscillations simultaneously in multiple frequency bands that have distinct functional roles in cognitive processing. Gamma oscillations (30-120 Hz) are associated with bottom-up processing, while slower oscillations in delta (1-4 Hz), theta (4-7 Hz), alpha (8-14 Hz) and beta (14-30 Hz) frequency bands may have roles in executive or top-down controlling functions, although also other distinctions have been made. Identification of the mechanisms that integrate such spectrally distributed processing and govern neuronal communication among these networks is crucial for understanding how cognitive functions are achieved in neuronal circuits. Cross-frequency interactions among oscillations have been recognized as a likely candidate mechanism for such integration. We advance here the hypothesis that phase-phase synchronization of neuronal oscillations in two different frequency bands, cross-frequency phase synchrony (CFS), could serve to integrate, coordinate and regulate neuronal processing distributed into neuronal assemblies concurrently in multiple frequency bands. A trail of studies over the past decade has revealed the presence of CFS among cortical oscillations and linked CFS with roles in cognitive integration. We propose that CFS could connect fast and slow oscillatory networks and thereby integrate distributed cognitive functions such as representation of sensory information with attentional and executive functions.

Critical dynamics of endogenous fluctuations predict cognitive flexibility in the Go/NoGo task

Fluctuations with power-law scaling and long-range temporal correlations (LRTCs) are characteristic to human psychophysical performance. Systems operating in a critical state exhibit such LRTCs, but phenomenologically similar fluctuations and LRTCs may also be caused by slow decay of the system’s memory without the system being critical. Theoretically, criticality endows the system with the greatest representational capacity and flexibility in state transitions. Without criticality, however, slowly decaying system memory would predict inflexibility. We addressed these contrasting predictions of the ‘criticality’ and ‘long-memory’ candidate mechanisms of human behavioral LRTCs by using a Go/NoGo task wherein the commission errors constitute a measure of cognitive flexibility. Response time (RT) fluctuations in this task exhibited power-law frequency scaling, autocorrelations, and LRTCs. We show here that the LRTC scaling exponents, quantifying the strength of long-range correlations, were negatively correlated with the commission error rates. Strong LRTCs hence parallel optimal cognitive flexibility and, in line with the criticality hypothesis, indicate a functionally advantageous state. This conclusion was corroborated by a positive correlation between the LRTC scaling exponents and executive functions measured with the Rey-Osterrieth Complex Figure test. Our results hence support the notion that LRTCs arise from critical dynamics that is functionally significant for human cognitive performance.

Modular co-organization of functional connectivity and scale-free dynamics in the human brain

Scale-free neuronal dynamics and interareal correlations are emergent characteristics of spontaneous brain activity. How such dynamics and the anatomical patterns of neuronal connectivity are mutually related in brain networks has, however, remained unclear. We addressed this relationship by quantifying the network colocalization of scale-free neuronal activity-both neuronal avalanches and long-range temporal correlations (LRTCs)-and functional connectivity (FC) by means of intracranial and noninvasive human resting-state electrophysiological recordings. We found frequency-specific colocalization of scale-free dynamics and FC so that the interareal couplings of LRTCs and the propagation of neuronal avalanches were most pronounced in the predominant pathways of FC. Several control analyses and the frequency specificity of network colocalization showed that the results were not trivial by-products of either brain dynamics or our analysis approach. Crucially, scale-free neuronal dynamics and connectivity also had colocalized modular structures at multiple levels of network organization, suggesting that modules of FC would be endowed with partially independent dynamic states. These findings thus suggest that FC and scale-free dynamics-hence, putatively, neuronal criticality as well-coemerge in a hierarchically modular structure in which the modules are characterized by dense connectivity, avalanche propagation, and shared dynamic states.

Cross-frequency synchronization connects networks of fast and slow oscillations during visual working memory maintenance

Neuronal activity in sensory and fronto-parietal (FP) areas underlies the representation and attentional control, respectively, of sensory information maintained in visual working memory (VWM). Within these regions, beta/gamma phase-synchronization supports the integration of sensory functions, while synchronization in theta/alpha bands supports the regulation of attentional functions. A key challenge is to understand which mechanisms integrate neuronal processing across these distinct frequencies and thereby the sensory and attentional functions. We investigated whether such integration could be achieved by cross-frequency phase synchrony (CFS). Using concurrent magneto- and electroencephalography, we found that CFS was load-dependently enhanced between theta and alpha–gamma and between alpha and beta-gamma oscillations during VWM maintenance among visual, FP, and dorsal attention (DA) systems. CFS also connected the hubs of within-frequency-synchronized networks and its strength predicted individual VWM capacity. We propose that CFS integrates processing among synchronized neuronal networks from theta to gamma frequencies to link sensory and attentional functions.

DOI:http://dx.doi.org/10.7554/eLife.13451.001

Experience Drives Synchronization: The phase and Amplitude Dynamics of Neural Oscillations to Musical Chords Are Differentially Modulated by Musical Expertise

Musical expertise is associated with structural and functional changes in the brain that underlie facilitated auditory perception. We investigated whether the phase locking (PL) and amplitude modulations (AM) of neuronal oscillations in response to musical chords are correlated with musical expertise and whether they reflect the prototypicality of chords in Western tonal music. To this aim, we recorded magnetoencephalography (MEG) while musicians and non-musicians were presented with common prototypical major and minor chords, and with uncommon, non-prototypical dissonant and mistuned chords, while watching a silenced movie. We then analyzed the PL and AM of ongoing oscillations in the theta (4–8 Hz) alpha (8–14 Hz), beta- (14–30 Hz) and gamma- (30–80 Hz) bands to these chords. We found that musical expertise was associated with strengthened PL of ongoing oscillations to chords over a wide frequency range during the first 300 ms from stimulus onset, as opposed to increased alpha-band AM to chords over temporal MEG channels. In musicians, the gamma-band PL was strongest to non-prototypical compared to other chords, while in non-musicians PL was strongest to minor chords. In both musicians and non-musicians the long-latency (> 200 ms) gamma-band PL was also sensitive to chord identity, and particularly to the amplitude modulations (beats) of the dissonant chord. These findings suggest that musical expertise modulates oscillation PL to musical chords and that the strength of these modulations is dependent on chord prototypicality.

Phase and amplitude correlations in resting-state activity in human stereotactical EEG recordings

Inter-areal interactions of neuronal oscillations may be a key mechanism in the coordination of anatomically distributed neuronal processing. In humans, invasive stereo-electroencephalography (SEEG) is emerging as a reference method for electrophysiological recordings because of its excellent spatial and temporal resolution. It could thus be also considered an optimal method for mapping neuronal inter-areal interactions. However, the common bipolar (BP) referencing of SEEG data may both confuse signals from distinct sources and suppress true neuronal interactions whereas the alternative monopolar (MP) reference yields data contaminated by volume conduction. We advance here a novel referencing scheme for SEEG data where electrodes in grey matter are referenced to closest white-matter (CW) electrodes. Using a 22 subject cohort and these three referencing schemes, we observed that both inter-areal phase and amplitude correlations decayed as function of distance and frequency but remained significant and stable across distances up to 10cm. Furthermore, we found that deep and superficial cortical laminae exhibit distinct spectral profiles of oscillation power as well as distinct patterns of inter-areal phase and amplitude interactions. These effects were qualitatively similar in MP and CW but distorted with BP referencing. Importantly CW was not influenced by the apparent large-scale volume conduction inherent to MP. We thus demonstrate here that with CW referencing, the superior anatomical accuracy of SEEG can be leveraged to yield accurate quantification and qualitatively novel insight into phase and amplitude interactions in human brain activity.

Gamma Oscillations Underlie the Maintenance of Feature-Specific Information and the Contents of Visual Working Memory

Visual working memory (VWM) sustains information online as integrated object representations. Neuronal mechanisms supporting the maintenance of feature-specific information have remained unidentified. Synchronized oscillations in the gamma band (30-120 Hz) characterize VWM retention and predict task performance, but whether these oscillations are specific to memorized features and VWM contents or underlie general executive VWM functions is not known. In the present study, we investigated whether gamma oscillations reflect the maintenance of feature-specific information in VWM. Concurrent magneto- and electroencephalography was recorded while subjects memorized different object features or feature conjunctions in identical VWM experiments. Using a data-driven source analysis approach, we show that the strength, load-dependence, and source topographies of gamma oscillations in the visual cortex differentiate these memorized features. Load-dependence of gamma oscillations in feature-specific visual and prefrontal areas also predicts VWM accuracy. Furthermore, corroborating the hypothesis that gamma oscillations support the perceptual binding of feature-specific neuronal assemblies, we also show that VWM for color-location conjunctions is associated with stronger gamma oscillations than that for these features separately. Gamma oscillations hence support the maintenance of feature-specific information and reflect VWM contents. The results also suggest that gamma oscillations contribute to feature binding in the formation of memory representations.

Discovering oscillatory interaction networks with M/EEG: challenges and breakthroughs

The systems-level neuronal mechanisms that coordinate temporally, anatomically and functionally distributed neuronal activity into coherent cognitive operations in the human brain have remained poorly understood. Synchronization of neuronal oscillations may regulate network communication and could thus serve as such a mechanism. Evidence for this hypothesis, however, was until recently sparse, as methodological challenges limit the investigation of interareal interactions with non-invasive magneto- and electroencephalography (M/EEG) recordings. Nevertheless, recent advances in M/EEG source reconstruction and clustering methods support complete phase-interaction mappings that are essential for uncovering the large-scale neuronal assemblies and their functional roles. These data show that synchronization is a robust and behaviorally significant phenomenon in task-relevant cortical networks and could hence bind distributed neuronal processing to coherent cognitive states.

Infra-slow fluctuations in electrophysiological recordings, blood-oxygenation-level-dependent signals, and psychophysical time series

Converging electrophysiological and neuroimaging data show that mammalian brain dynamics are governed by spontaneous modulations of neuronal activity levels in cortical and subcortical structures. The time scales of these fluctuations form a continuum from seconds to tens and hundreds of seconds corresponding to slow (0.1-1Hz), infra-slow (0.01-0.1Hz), and "ultradian" (<0.01Hz) frequency bands, respectively. We focus here on the spontaneous neuronal dynamics in the infra-slow frequency band, infra-slow fluctuations (ISFs), and explore their electrophysiological substrates and behavioral correlates. Although electrophysiological ISFs and the associated infra-slow modulations of fast (here, >1Hz) neuronal activities have been recognized on numerous occasions since late 50's, a resurgence in interest towards this frequency band has been driven by a discovery that ISFs in blood-oxygenation-level dependent (BOLD) signals are correlated among specific constellations of brain regions, which constitute intrinsic connectivity networks and define the dynamic architecture of spontaneous brain activity at large. Importantly, electrophysiological and BOLD signal ISFs are directly correlated both with ISFs in amplitudes of fast neuronal activities and with ISFs in behavioral performance. Moreover, both electrophysiological and neuroimaging data suggest that the apparently scale-free ISFs may arise from more local quasi-periodic infra-slow oscillations with a contribution of time-scale-specific cellular-level mechanisms. We conclude that ISFs in electrophysiological recordings, BOLD signals, neuronal activity levels, and behavioral time series are likely to reflect the same underlying phenomenon; a superstructure of interacting and transiently oscillatory ISFs that regulate both the integration within and decoupling between concurrently active neuronal communities.

Functional roles of alpha-band phase synchronization in local and large-scale cortical networks

Alpha-frequency band (8-14 Hz) oscillations are among the most salient phenomena in human electroencephalography (EEG) recordings and yet their functional roles have remained unclear. Much of research on alpha oscillations in human EEG has focused on peri-stimulus amplitude dynamics, which phenomenologically support an idea of alpha oscillations being negatively correlated with local cortical excitability and having a role in the suppression of task-irrelevant neuronal processing. This kind of an inhibitory role for alpha oscillations is also supported by several functional magnetic resonance imaging and trans-cranial magnetic stimulation studies. Nevertheless, investigations of local and inter-areal alpha phase dynamics suggest that the alpha-frequency band rhythmicity may play a role also in active task-relevant neuronal processing. These data imply that inter-areal alpha phase synchronization could support attentional, executive, and contextual functions. In this review, we outline evidence supporting different views on the roles of alpha oscillations in cortical networks and unresolved issues that should be addressed to resolve or reconcile these apparently contrasting hypotheses.

Functional roles of alpha-band phase synchronization in local and large-scale cortical networks

Alpha-frequency band (8–14 Hz) oscillations are among the most salient phenomena in human electroencephalography (EEG) recordings and yet their functional roles have remained unclear. Much of research on alpha oscillations in human EEG has focused on peri-stimulus amplitude dynamics, which phenomenologically support an idea of alpha oscillations being negatively correlated with local cortical excitability and having a role in the suppression of task-irrelevant neuronal processing. This kind of an inhibitory role for alpha oscillations is also supported by several functional magnetic resonance imaging and trans-cranial magnetic stimulation studies. Nevertheless, investigations of local and inter-areal alpha phase dynamics suggest that the alpha-frequency band rhythmicity may play a role also in active task-relevant neuronal processing. These data imply that inter-areal alpha phase synchronization could support attentional, executive, and contextual functions. In this review, we outline evidence supporting different views on the roles of alpha oscillations in cortical networks and unresolved issues that should be addressed to resolve or reconcile these apparently contrasting hypotheses.

Roles of multiscale brain activity fluctuations in shaping the variability and dynamics of psychophysical performance

Spontaneous brain activity across many time scales influences sensory perception and human cognitive performance. Empirical insight into the underlying systems-level mechanisms has, however, remained fragmented. We review here recent studies on how wideband scale-free and scale-specific neuronal activity fluctuations together bias sensory processing and perceptual performance. We posit that these fluctuations constitute the neurophysiological foundation for both the trial-to-trial behavioral variability and the scaling laws governing psychophysical performance.

Optimization of an NLEO-based algorithm for automated detection of spontaneous activity transients in early preterm EEG

We propose here a simple algorithm for automated detection of spontaneous activity transients (SATs) in early preterm electroencephalography (EEG). The parameters of the algorithm were optimized by supervised learning using a gold standard created from visual classification data obtained from three human raters. The generalization performance of the algorithm was estimated by leave-one-out cross-validation. The mean sensitivity of the optimized algorithm was 97% (range 91-100%) and specificity 95% (76-100%). The optimized algorithm makes it possible to systematically study brain state fluctuations of preterm infants.

Graph properties of synchronized cortical networks during visual working memory maintenance

Oscillatory synchronization facilitates communication in neuronal networks and is intimately associated with human cognition. Neuronal activity in the human brain can be non-invasively imaged with magneto- (MEG) and electroencephalography (EEG), but the large-scale structure of synchronized cortical networks supporting cognitive processing has remained uncharacterized. We combined simultaneous MEG and EEG (MEEG) recordings with minimum-norm-estimate-based inverse modeling to investigate the structure of oscillatory phase synchronized networks that were active during visual working memory (VWM) maintenance. Inter-areal phase-synchrony was quantified as a function of time and frequency by single-trial phase-difference estimates of cortical patches covering the entire cortical surfaces. The resulting networks were characterized with a number of network metrics that were then compared between delta/theta- (3-6 Hz), alpha- (7-13 Hz), beta- (16-25 Hz), and gamma- (30-80 Hz) frequency bands. We found several salient differences between frequency bands. Alpha- and beta-band networks were more clustered and small-world like but had smaller global efficiency than the networks in the delta/theta and gamma bands. Alpha- and beta-band networks also had truncated-power-law degree distributions and high k-core numbers. The data converge on showing that during the VWM-retention period, human cortical alpha- and beta-band networks have a memory-load dependent, scale-free small-world structure with densely connected core-like structures. These data further show that synchronized dynamic networks underlying a specific cognitive state can exhibit distinct frequency-dependent network structures that could support distinct functional roles.

New vistas for alpha-frequency band oscillations

The amplitude of alpha-frequency band (8-14 Hz) activity in the human electroencephalogram is suppressed by eye opening, visual stimuli and visual scanning, whereas it is enhanced during internal tasks, such as mental calculation and working memory. alpha-Frequency band oscillations have hence been thought to reflect idling or inhibition of task-irrelevant cortical areas. However, recent data on alpha-amplitude and, in particular, alpha-phase dynamics posit a direct and active role for alpha-frequency band rhythmicity in the mechanisms of attention and consciousness. We propose that simultaneous alpha-, beta- (14-30 Hz) and gamma- (30-70 Hz) frequency band oscillations are required for unified cognitive operations, and hypothesize that cross-frequency phase synchrony between alpha, beta and gamma oscillations coordinates the selection and maintenance of neuronal object representations during working memory, perception and consciousness.

Development of the spontaneous activity transients and ongoing cortical activity in human preterm babies

Recent experimental studies have shown that developing cortex in several animals species, including humans, exhibits spontaneous intermittent activity that is believed to be crucial for the proper wiring of early brain networks. The present study examined the developmental changes in these spontaneous activity transients (SAT) and in other ongoing cortical activities in human preterm babies. Full-band electroencephalography (FbEEG) recordings were obtained from 16 babies at conceptional ages between 32.8 and 40 wk. We examined the SATs and the intervening ongoing cortical activities (inter-SAT; iSAT) with average waveforms, their variance and power, as well as with wavelet-based time-frequency analyses. Our results show, that the low frequency power and the variance of the average waveform of SAT decrease during development. There was a simultaneous increase in the activity at higher frequencies, with most pronounced increase at theta-alpha range (4-9 Hz). In addition to the overall increase, the activity at higher frequencies showed an increased grouping into bursts that are nested in the low frequency (0.5-1 Hz) waves. Analysis of the iSAT epochs showed a developmental increase in power at lower frequencies in quiet sleep. There was an increase in a wide range of higher frequencies (4-16 Hz), whereas the ratio of beta (16-30 Hz) and theta-alpha (4-9 Hz) range activity declined, indicating a preferential increase at theta-alpha range activity. Notably, SAT and iSAT activities remained distinct throughout the development in all measures used in our study. The present results are consistent with the idea that SAT and the other ongoing cortical activities are distinct functional entities. Recognition of these two basic mechanisms in the cortical activity in preterm human babies opens new rational approaches for an evaluation and monitoring of early human brain function.

Epileptogenic neocortical networks are revealed by abnormal temporal dynamics in seizure-free subdural EEG

Long-term video electroencephalographic (EEG) recording is currently a routine procedure in the presurgical evaluation of localization-related epilepsies. Cortical epileptogenic zone is usually localized from ictal recordings with intracranial electrodes, causing a significant burden to patients and health care. Growing literature suggests that epileptogenic networks exhibit aberrant dynamics also during seizure-free periods. We examined if neocortical epileptogenic regions can be circumscribed by quantifying local long-range temporal (auto-)correlations (LRTC) with detrended fluctuation analysis of seizure-free ongoing subdural EEG activity in 4 frequency bands in 5 patients. We show here with subdural EEG recordings that the LRTC are abnormally strong near the seizure onset area. This effect was most salient in neocortical oscillations in the beta frequency band (14-30 Hz). Moreover, lorazepam, a widely used antiepileptic drug, exerted contrasting effects on LRTC (n = 2): lorazepam attenuated beta-band LRTC near the epileptic focus, whereas it strengthened LRTC in other cortical areas. Our findings demonstrate that interictal neuronal network activity near the focus of seizure onset has pathologically strong intrinsic temporal correlations. The observed effect by lorazepam on beta-band activity suggests that the antiepileptic mechanism of benzodiazepines may be related to the normalization of LRTC within the epileptic focus. We propose that this method may become a promising candidate for routine invasive and noninvasive presurgical localization of epileptic foci.

Early neural correlates of conscious somatosensory perception

The cortical processing of consciously perceived and unperceived somatosensory stimuli is thought to be identical during the first 100-120 ms after stimulus onset. Thereafter, the electrophysiological correlates of conscious perception have been shown to be reflected in the N1 component of the evoked response as well as in later (>200 ms) nonstimulus-locked gamma-band (28-50 Hz) oscillatory activity. To evaluate more specifically the time course and correlation of neuronal oscillations with conscious perception, we recorded neuromagnetic responses to threshold-intensity somatosensory stimuli. We show here that cortical broadband activities phase locked to the subsequently perceived stimuli in somatosensory, frontal, and parietal regions as early as 30-70 ms from stimulus onset, whereas the phase locking to the unperceived stimuli was weak and primarily restricted to somatosensory regions. Such stimulus locking also preceded the perceived stimuli, indicating that the phase of ongoing cortical activities biases subsequent perception. Furthermore, the data show that the stimulus locking was present in the theta- (4-8 Hz), alpha- (8-14 Hz), beta- (14-28 Hz), and gamma- (28-40 Hz) frequency bands, of which the widespread alpha-band component was dominant for the consciously perceived stimuli but virtually unobservable for the unperceived stimuli. Our results show that the neural correlates of conscious perception are already found during the earliest stages of cortical processing from 30 to 150 ms after stimulus onset and suggest that alpha-frequency-band oscillations have a role in the neural mechanisms of sensory awareness.

Phase synchrony among neuronal oscillations in the human cortex

Synchronization of neuronal activity, often associated with network oscillations, is thought to provide a means for integrating anatomically distributed processing in the brain. Neuronal processing, however, involves simultaneous oscillations in various frequency bands. The mechanisms involved in the integration of such spectrally distributed processing have remained enigmatic. We demonstrate, using magnetoencephalography, that robust cross-frequency phase synchrony is present in the human cortex among oscillations with frequencies from 3 to 80 Hz. Continuous mental arithmetic tasks demanding the retention and summation of items in the working memory enhanced the cross-frequency phase synchrony among alpha (approximately 10 Hz), beta (approximately 20 Hz), and gamma (approximately 30-40 Hz) oscillations. These tasks also enhanced the "classical" within-frequency synchrony in these frequency bands, but the spatial patterns of alpha, beta, and gamma synchronies were distinct and, furthermore, separate from the patterns of cross-frequency phase synchrony. Interestingly, an increase in task load resulted in an enhancement of phase synchrony that was most prominent between gamma- and alpha-band oscillations. These data indicate that cross-frequency phase synchrony is a salient characteristic of ongoing activity in the human cortex and that it is modulated by cognitive task demands. The enhancement of cross-frequency phase synchrony among functionally and spatially distinct networks during mental arithmetic tasks posits it as a candidate mechanism for the integration of spectrally distributed processing.

Slow endogenous activity transients and developmental expression of K+-Cl- cotransporter 2 in the immature human cortex

Spontaneous transients of correlated activity are a characteristic feature of immature brain structures, where they are thought to be crucial for the establishment of precise neuronal connectivity. Studies on experimental animals have shown that this kind of early activity in cortical structures is composed of long-lasting, intermittent network events, which undergo a developmental decline that is closely paralleled by the maturation of GABAergic inhibition. In order to examine whether similar events occur in the immature human cortex, we performed direct current-coupled electroencephalography (EEG) recordings from sleeping preterm babies. We show now that much of the preterm EEG activity is confined to spontaneous, slow activity transients. These transients are characterized by a large voltage deflection that nests prominent oscillatory activity in several frequency bands covering the whole frequency spectrum of the preterm EEG (<0.1-30 Hz). The slow voltage deflections had an amplitude of up to 800 microV. Most of these 'giant' events originated in the temporo-occipital areas, with a maximum rate of about 8/min, and their occurrence as well as amplitude showed a decline by the time of normal birth. In age-matched fetal brain tissue, this decrease in the spontaneous activity transients was associated with a developmental up-regulation of the neuronal chloride extruder K+-Cl- cotransporter 2, a crucial molecule for the generation of inhibitory GABAergic Cl- currents. Our work indicates that slow endogenous activity transients in the immature human neocortex are mostly confined to the prenatal stage and appear to be terminated in parallel with the maturation of functional GABAergic inhibition.