Sparse weightings for collapsing inverse solutions to cortical parcellations optimize M/EEG source reconstruction accuracy

Characterization of Second-Order Mixing Effects in Reconstructed Cross-Spectra of Random Neural Fields

Functional connectivity in electroencephalography (EEG) and magnetoencephalography (MEG) data is commonly assessed by using measures that are insensitive to instantaneously interacting sources and as such would not give rise to false positive interactions caused by instantaneous mixing of true source signals (first-order mixing). Recent studies, however, have drawn attention to the fact that such measures are still susceptible to instantaneous mixing from lagged sources (i.e. second-order mixing) and that this can lead to a large number of false positive interactions. In this study we relate first- and second-order mixing effects on the cross-spectra of reconstructed source activity to the properties of the resolution operators that are used for the reconstruction. We derive two identities that relate first- and second-order mixing effects to the transformation properties of measurement and source configurations and exploit them to establish several basic properties of signal mixing. First, we provide a characterization of the configurations that are maximally and minimally sensitive to second-order mixing. It turns out that second-order mixing effects are maximal when the measurement locations are far apart and the sources coincide with the measurement locations. Second, we provide a description of second-order mixing effects in the vicinity of the measurement locations in terms of the local geometry of the point-spread functions of the resolution operator. Third, we derive a version of Lagrange's identity for cross-talk functions that establishes the existence of a trade-off between the magnitude of first- and second-order mixing effects. It also shows that, whereas the magnitude of first-order mixing is determined by the inner product of cross-talk functions, the magnitude of second-order mixing is determined by a generalized cross-product of cross-talk functions (the wedge product) which leads to an intuitive geometric understanding of the trade-off. All results are derived within the general framework of random neural fields on cortical manifolds.

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.

Strengthening of alpha synchronization is a neural correlate of cognitive transfer

Cognitive training can lead to improvements in both task-specific strategies and general capacities, such as visuo-spatial working memory (VSWM). The latter emerge slowly and linearly throughout training, in contrast to strategy where changes typically occur within the first days of training. Changes in strategy and capacity have not been separated in prior neuroimaging studies. Here, we used a within-participants design with dense temporal sampling to capture the time dynamics of neural mechanisms associated with change in capacity. In four participants, neural activity was recorded with magnetoencephalography on seven occasions over two months of visuo-spatial working memory training. During scanning, the participants performed a trained visuo-spatial working memory task, a transfer task, and a control task. First, we extracted an individual visuo-spatial working memory-load-dependent synchronization network for each participant. Next, we identified linear changes over time in the network, congruent with the temporal dynamics of capacity change. Three out of four participants showed a gradual strengthening of alpha synchronization. Strengthening of the same connections was also found in the transfer task but not in the control task. This suggests that cognitive transfer occurs through slow, gradual strengthening of alpha synchronization between cortical regions that are vital for both the trained task and the transfer task.

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.

PSIICOS projection optimality for EEG and MEG based functional coupling detection

Functional connectivity is crucial for cognitive processes in the healthy brain and serves as a marker for a range of neuropathological conditions. Non-invasive exploration of functional coupling using temporally resolved techniques such as MEG allows for a unique opportunity of exploring this fundamental brain mechanism. The indirect nature of MEG measurements complicates the estimation of functional coupling due to the volume conduction and spatial leakage effects. In the previous work (Ossadtchi et al., 2018), we introduced PSIICOS, a method that for the first time allowed us to suppress the volume conduction effect and yet retain information about functional networks whose nodes are coupled with close to zero or zero mutual phase lag. In this paper, we demonstrate analytically that the PSIICOS projection is optimal in achieving a controllable trade-off between suppressing mutual spatial leakage and retaining information about zero- or close to zero-phase coupled networks. We also derive an alternative solution using the regularization-based inverse of the mutual spatial leakage matrix and show its equivalence to the original PSIICOS. We then discuss how PSIICOS solution to the functional connectivity estimation problem can be incorporated into the conventional source estimation framework. Instead of sources, the unknowns are the elementary dyadic networks and their activation time series are formalized by the corresponding source-space cross-spectral coefficients. This view on connectivity estimation as a regression problem opens up new opportunities for formulating a set of principled estimators based on the rich intuition accumulated in the neuroimaging community.

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.

Identifying good practices for detecting inter-regional linear functional connectivity from EEG

Aggregating voxel-level statistical dependencies between multivariate time series is an important intermediate step when characterising functional connectivity (FC) between larger brain regions. However, there are numerous ways in which voxel-level data can be aggregated into inter-regional FC, and the advantages of each of these approaches are currently unclear. In this study we generate ground-truth data and compare the performances of various pipelines that estimate directed and undirected linear phase-to-phase FC between regions. We test the ability of several existing and novel FC analysis pipelines to identify the true regions within which connectivity was simulated. We test various inverse modelling algorithms, strategies to aggregate time series within regions, and connectivity metrics. Furthermore, we investigate the influence of the number of interactions, the signal-to-noise ratio, the noise mix, the interaction time delay, and the number of active sources per region on the ability of detecting phase-to-phase FC. Throughout all simulated scenarios, lowest performance is obtained with pipelines involving the absolute value of coherency. Further, the combination of dynamic imaging of coherent sources (DICS) beamforming with directed FC metrics that aggregate information across multiple frequencies leads to unsatisfactory results. Pipelines that show promising results with our simulated pseudo-EEG data involve the following steps: (1) Source projection using the linearly-constrained minimum variance (LCMV) beamformer. (2) Principal component analysis (PCA) using the same fixed number of components within every region. (3) Calculation of the multivariate interaction measure (MIM) for every region pair to assess undirected phase-to-phase FC, or calculation of time-reversed Granger Causality (TRGC) to assess directed phase-to-phase FC. We formulate recommendations based on these results that may increase the validity of future experimental connectivity studies. We further introduce the free ROIconnect plugin for the EEGLAB toolbox that includes the recommended methods and pipelines that are presented here. We show an exemplary application of the best performing pipeline to the analysis of EEG data recorded during motor imagery.

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.

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) Val158Met and brain-derived neurotrophic factor (BDNF) Val66Met. 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.

Harmoni: a Method for Eliminating Spurious Interactions due to the Harmonic Components in Neuronal Data

Cross-frequency synchronization (CFS) has been proposed as a mechanism for integrating spatially and spectrally distributed information in the brain. However, investigating CFS in Magneto- and Electroencephalography (MEG/EEG) is hampered by the presence of spurious neuronal interactions due to the non-sinusoidal waveshape of brain oscillations. Such waveshape gives rise to the presence of oscillatory harmonics mimicking genuine neuronal oscillations. Until recently, however, there has been no methodology for removing these harmonics from neuronal data. In order to address this long-standing challenge, we introduce a novel method (called HARMOnic miNImization - Harmoni) that removes the signal components which can be harmonics of a non-sinusoidal signal. Harmoni's working principle is based on the presence of CFS between harmonic components and the fundamental component of a non-sinusoidal signal. We extensively tested Harmoni in realistic EEG simulations. The simulated couplings between the source signals represented genuine and spurious CFS and within-frequency phase synchronization. Using diverse evaluation criteria, including ROC analyses, we showed that the within- and cross-frequency spurious interactions are suppressed significantly, while the genuine activities are not affected. Additionally, we applied Harmoni to real resting-state EEG data revealing intricate remote connectivity patterns which are usually masked by the spurious connections. Given the ubiquity of non-sinusoidal neuronal oscillations in electrophysiological recordings, Harmoni is expected to facilitate novel insights into genuine neuronal interactions in various research fields, and can also serve as a steppingstone towards the development of further signal processing methods aiming at refining within- and cross-frequency synchronization in electrophysiological recordings.

Principles and open questions in functional brain network reconstruction

Graph theory is now becoming a standard tool in system-level neuroscience. However, endowing observed brain anatomy and dynamics with a complex network representation involves often covert theoretical assumptions and methodological choices which affect the way networks are reconstructed from experimental data, and ultimately the resulting network properties and their interpretation. Here, we review some fundamental conceptual underpinnings and technical issues associated with brain network reconstruction, and discuss how their mutual influence concurs in clarifying the organization of brain function.

Reduced evoked activity and cortical oscillations are correlated with anisometric amblyopia and impairment of visual acuity

Amblyopia is a developmental disorder associated with abnormal visual experience during early childhood commonly arising from strabismus and/or anisometropia and leading to dysfunctions in visual cortex and to various visual deficits. The different forms of neuronal activity that are attenuated in amblyopia have been only partially characterized. In electrophysiological recordings of healthy human brain, the presentation of visual stimuli is associated with event-related activity and oscillatory responses. It has remained poorly understood whether these forms of activity are reduced in amblyopia and whether possible dysfunctions would arise from lower- or higher-order visual areas. We recorded neuronal activity with magnetoencephalography (MEG) from anisometropic amblyopic patients and control participants during two visual tasks presented separately for each eye and estimated neuronal activity from source-reconstructed MEG data. We investigated whether event-related and oscillatory responses would be reduced for amblyopia and localized their cortical sources. Oscillation amplitudes and evoked responses were reduced for stimuli presented to the amblyopic eye in higher-order visual areas and in parietal and prefrontal cortices. Importantly, the reduction of oscillation amplitudes but not that of evoked responses was correlated with decreased visual acuity in amblyopia. These results show that attenuated oscillatory responses are correlated with visual deficits in anisometric amblyopia.

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.

Genuine cross-frequency coupling networks in human resting-state electrophysiological recordings

Phase synchronization of neuronal oscillations in specific frequency bands coordinates anatomically distributed neuronal processing and communication. Typically, oscillations and synchronization take place concurrently in many distinct frequencies, which serve separate computational roles in cognitive functions. While within-frequency phase synchronization has been studied extensively, less is known about the mechanisms that govern neuronal processing distributed across frequencies and brain regions. Such integration of processing between frequencies could be achieved via cross-frequency coupling (CFC), either by phase–amplitude coupling (PAC) or by n:m-cross–frequency phase synchrony (CFS). So far, studies have mostly focused on local CFC in individual brain regions, whereas the presence and functional organization of CFC between brain areas have remained largely unknown. We posit that interareal CFC may be essential for large-scale coordination of neuronal activity and investigate here whether genuine CFC networks are present in human resting-state (RS) brain activity. To assess the functional organization of CFC networks, we identified brain-wide CFC networks at mesoscale resolution from stereoelectroencephalography (SEEG) and at macroscale resolution from source-reconstructed magnetoencephalography (MEG) data. We developed a novel, to our knowledge, graph-theoretical method to distinguish genuine CFC from spurious CFC that may arise from nonsinusoidal signals ubiquitous in neuronal activity. We show that genuine interareal CFC is present in human RS activity in both SEEG and MEG data. Both CFS and PAC networks coupled theta and alpha oscillations with higher frequencies in large-scale networks connecting anterior and posterior brain regions. CFS and PAC networks had distinct spectral patterns and opposing distribution of low- and high-frequency network hubs, implying that they constitute distinct CFC mechanisms. The strength of CFS networks was also predictive of cognitive performance in a separate neuropsychological assessment. In conclusion, these results provide evidence for interareal CFS and PAC being 2 distinct mechanisms for coupling oscillations across frequencies in large-scale brain networks.

Genuine interareal cross-frequency coupling (CFC) can be identified from human resting state activity using magnetoencephalography, stereoelectroencephalography, and novel network approaches. CFC couples slow theta and alpha oscillations to faster oscillations across brain regions.

Neuronal correlates of full and partial visual conscious perception

Stimuli may induce only partial consciousness-an intermediate between null and full consciousness-where the presence but not identity of an object can be reported. The differences in the neuronal basis of full and partial consciousness are poorly understood. We investigated if evoked and oscillatory activity could dissociate full from partial conscious perception. We recorded human cortical activity with magnetoencephalography (MEG) during a visual perception task in which stimulus could be either partially or fully perceived. Partial consciousness was associated with an early increase in evoked activity and theta/low-alpha-band oscillations while full consciousness was also associated with late evoked activity and beta-band oscillations. Full from partial consciousness was dissociated by stronger evoked activity and late increase in theta oscillations that were localized to higher-order visual regions and posterior parietal and prefrontal cortices. Our results reveal both evoked activity and theta oscillations dissociate partial and full consciousness.

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 locking value revisited: teaching new tricks to an old dog

OBJECTIVE

Despite the increase in calculation power over the last few decades, the estimation of brain connectivity is still a tedious task. The high computational cost of the algorithms escalates with the square of the number of signals evaluated, usually in the range of thousands. In this work we propose a re-formulation of a widely used algorithm that allows the estimation of whole brain connectivity in much smaller times.

APPROACH

We start from the original implementation of phase locking value (PLV) and re-formulated it in a computationally very efficient way. What is more, this formulation stresses its strong similarity with coherence, which we used to introduce two new metrics insensitive to zero lag synchronization: the imaginary part of PLV (iPLV) and its corrected counterpart (ciPLV).

MAIN RESULTS

The new implementation of PLV avoids some highly CPU-expensive operations and achieves a 100-fold speedup over the original algorithm. The new derived metrics were highly robust in the presence of volume conduction. Moreover, ciPLV proved capable of ignoring zero-lag connectivity, while correctly estimating nonzero-lag connectivity.

SIGNIFICANCE

Our implementation of PLV makes it possible to calculate whole-brain connectivity in much shorter times. The results of the simulations using ciPLV suggest that this metric is ideal to measure synchronization in the presence of volume conduction or source leakage effects.

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.

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.

Whole-Brain Source-Reconstructed MEG-Data Reveal Reduced Long-Range Synchronization in Chronic Schizophrenia

Current theories of schizophrenia (ScZ) posit that the symptoms and cognitive dysfunctions arise from a dysconnection syndrome. However, studies that have examined this hypothesis with physiological data at realistic time scales are so far scarce. The current study employed a state-of-the-art approach using Magnetoencephalography (MEG) to test alterations in large-scale phase synchronization in a sample of n = 16 chronic ScZ patients, 10 males and n = 19 healthy participants, 10 males, during a perceptual closure task. We identified large-scale networks from source reconstructed MEG data using data-driven analyses of neuronal synchronization. Oscillation amplitudes and interareal phase-synchronization in the 3–120 Hz frequency range were estimated for 400 cortical parcels and correlated with clinical symptoms and neuropsychological scores. ScZ patients were characterized by a reduction in γ-band (30–120 Hz) oscillation amplitudes that was accompanied by a pronounced deficit in large-scale synchronization at γ-band frequencies. Synchronization was reduced within visual regions as well as between visual and frontal cortex and the reduction of synchronization correlated with elevated clinical disorganization. Accordingly, these data highlight that ScZ is associated with a profound disruption of transient synchronization, providing critical support for the notion that core aspect of the pathophysiology arises from an impairment in coordination of distributed neural activity.

Adaptive cortical parcellations for source reconstructed EEG/MEG connectomes

There is growing interest in the rich temporal and spectral properties of the functional connectome of the brain that are provided by Electro- and Magnetoencephalography (EEG/MEG). However, the problem of leakage between brain sources that arises when reconstructing brain activity from EEG/MEG recordings outside the head makes it difficult to distinguish true connections from spurious connections, even when connections are based on measures that ignore zero-lag dependencies. In particular, standard anatomical parcellations for potential cortical sources tend to over- or under-sample the real spatial resolution of EEG/MEG. By using information from cross-talk functions (CTFs) that objectively describe leakage for a given sensor configuration and distributed source reconstruction method, we introduce methods for optimising the number of parcels while simultaneously minimising the leakage between them. More specifically, we compare two image segmentation algorithms: 1) a split-and-merge (SaM) algorithm based on standard anatomical parcellations and 2) a region growing (RG) algorithm based on all the brain vertices with no prior parcellation. Interestingly, when applied to minimum-norm reconstructions for EEG/MEG configurations from real data, both algorithms yielded approximately 70 parcels despite their different starting points, suggesting that this reflects the resolution limit of this particular sensor configuration and reconstruction method. Importantly, when compared against standard anatomical parcellations, resolution matrices of adaptive parcellations showed notably higher sensitivity and distinguishability of parcels. Furthermore, extensive simulations of realistic networks revealed significant improvements in network reconstruction accuracies, particularly in reducing false leakage-induced connections. Adaptive parcellations therefore allow a more accurate reconstruction of functional EEG/MEG connectomes.

•

We introduce adaptive cortical parcellation algorithms for E/MEG source estimation.

•

Optimum number, size and locations of parcels are found based on cross-talk functions

•

Algorithms yielded ∼70 distinguishable parcels regardless of the starting point.

•

Parcel resolution matrices were notably improved compared to anatomical atlases.

•

Network reconstruction accuracies of simulated connectomes improved significantly.

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

The role of cortical beta oscillations in time estimation

Estimation of time is central to perception, action, and cognition. Human functional magnetic resonance imaging (fMRI) and positron emission topography (PET) have revealed a positive correlation between the estimation of multi-second temporal durations and neuronal activity in a circuit of sensory and motor areas, prefrontal and temporal cortices, basal ganglia, and cerebellum. The systems-level mechanisms coordinating the collective neuronal activity in these areas have remained poorly understood. Synchronized oscillations regulate communication in neuronal networks and could hence serve such coordination, but their role in the estimation and maintenance of multi-second time intervals has remained largely unknown. We used source-reconstructed magnetoencephalography (MEG) to address the functional significance of local neuronal synchronization, as indexed by the amplitudes of cortical oscillations, in time-estimation. MEG was acquired during a working memory (WM) task where the subjects first estimated and then memorized the durations, or in the contrast condition, the colors of dynamic visual stimuli. Time estimation was associated with stronger beta (β, 14 - 30 Hz) band oscillations than color estimation in sensory regions and attentional cortical structures that earlier have been associated with time processing. In addition, the encoding of duration information was associated with strengthened gamma- (γ, 30 - 120 Hz), and the retrieval and maintenance with alpha- (α, 8 - 14 Hz) band oscillations. These data suggest that β oscillations may provide a mechanism for estimating short temporal durations, while γ and α oscillations support their encoding, retrieval, and maintenance in memory. Hum Brain Mapp 37:3262-3281, 2016. © 2016 Wiley Periodicals, Inc.

Cortical localization of phase and amplitude dynamics predicting access to somatosensory awareness

Neural dynamics leading to conscious sensory perception have remained enigmatic in despite of large interest. Human functional magnetic resonance imaging (fMRI) studies have revealed that a co-activation of sensory and frontoparietal areas is crucial for conscious sensory perception in the several second time-scale of BOLD signal fluctuations. Electrophysiological recordings with magneto- and electroencephalography (MEG and EEG) and intracranial EEG (iEEG) have shown that event related responses (ERs), phase-locking of neuronal activity, and oscillation amplitude modulations in sub-second timescales are greater for consciously perceived than for unperceived stimuli. The cortical sources of ER and oscillation dynamics predicting the conscious perception have, however, remained unclear because these prior studies have utilized MEG/EEG sensor-level analyses or iEEG with limited neuroanatomical coverage. We used a somatosensory detection task, magnetoencephalography (MEG), and cortically constrained source reconstruction to identify the cortical areas where ERs, local poststimulus amplitudes and phase-locking of neuronal activity are predictive of the conscious access of somatosensory information. We show here that strengthened ERs, phase-locking to stimulus onset (SL), and induced oscillations amplitude modulations all predicted conscious somatosensory perception, but the most robust and widespread of these was SL that was sustained in low-alpha (6-10 Hz) band. The strength of SL and to a lesser extent that of ER predicted conscious perception in the somatosensory, lateral and medial frontal, posterior parietal, and in the cingulate cortex. These data suggest that a rapid phase-reorganization and concurrent oscillation amplitude modulations in these areas play an instrumental role in the emergence of a conscious percept.

Relationship of fast- and slow-timescale neuronal dynamics in human MEG and SEEG

A growing body of evidence suggests that the neuronal dynamics are poised at criticality. Neuronal avalanches and long-range temporal correlations (LRTCs) are hallmarks of such critical dynamics in neuronal activity and occur at fast (subsecond) and slow (seconds to hours) timescales, respectively. The critical dynamics at different timescales can be characterized by their power-law scaling exponents. However, insight into the avalanche dynamics and LRTCs in the human brain has been largely obtained with sensor-level MEG and EEG recordings, which yield only limited anatomical insight and results confounded by signal mixing. We investigated here the relationship between the human neuronal dynamics at fast and slow timescales using both source-reconstructed MEG and intracranial stereotactical electroencephalography (SEEG). Both MEG and SEEG revealed avalanche dynamics that were characterized parameter-dependently by power-law or truncated-power-law size distributions. Both methods also revealed robust LRTCs throughout the neocortex with distinct scaling exponents in different functional brain systems and frequency bands. The exponents of power-law regimen neuronal avalanches and LRTCs were strongly correlated across subjects. Qualitatively similar power-law correlations were also observed in surrogate data without spatial correlations but with scaling exponents distinct from those of original data. Furthermore, we found that LRTCs in the autonomous nervous system, as indexed by heart-rate variability, were correlated in a complex manner with cortical neuronal avalanches and LRTCs in MEG but not SEEG. These scalp and intracranial data hence show that power-law scaling behavior is a pervasive but neuroanatomically inhomogeneous property of neuronal dynamics in central and autonomous nervous systems.

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.