Post-task responses following working memory and movement are driven by transient spectral bursts with similar characteristics.

The post-movement beta rebound has been studied extensively using magnetoencephalography (MEG) and is reliably modulated by various task parameters as well as illness. Our recent study showed that rebounds, which we generalise as "post-task responses" (PTRs), are a ubiquitous phenomenon in the brain, occurring across the cortex in theta, alpha, and beta bands. Currently, it is unknown whether PTRs following working memory are driven by transient bursts, which are moments of short-lived high amplitude activity, similar to those that drive the post-movement beta rebound. Here, we use three-state univariate hidden Markov models (HMMs), which can identify bursts without a priori knowledge of frequency content or response timings, to compare bursts that drive PTRs in working memory and visuomotor MEG datasets. Our results show that PTRs across working memory and visuomotor tasks are driven by pan-spectral transient bursts. These bursts have very similar spectral content variation over the cortex, correlating strongly between the two tasks in the alpha (R2 = .89) and beta (R2 = .53) bands. Bursts also have similar variation in duration over the cortex (e.g., long duration bursts occur in the motor cortex for both tasks), strongly correlating over cortical regions between tasks (R2 = .56), with a mean over all regions of around 300 ms in both datasets. Finally, we demonstrate the ability of HMMs to isolate signals of interest in MEG data, such that the HMM probability timecourse correlates more strongly with reaction times than frequency filtered power envelopes from the same brain regions. Overall, we show that induced PTRs across different tasks are driven by bursts with similar characteristics, which can be identified using HMMs. Given the similarity between bursts across tasks, we suggest that PTRs across the cortex may be driven by a common underlying neural phenomenon.

Dissecting unsupervised learning through hidden Markov modeling in electrophysiological data.

Unsupervised, data-driven methods are commonly used in neuroscience to automatically decompose data into interpretable patterns. These patterns differ from one another depending on the assumptions of the models. How these assumptions affect specific data decompositions in practice, however, is often unclear, which hinders model applicability and interpretability. For instance, the hidden Markov model (HMM) automatically detects characteristic, recurring activity patterns (so-called states) from time series data. States are defined by a certain probability distribution, whose state-specific parameters are estimated from the data. But what specific features, from all of those that the data contain, do the states capture? That depends on the choice of probability distribution and on other model hyperparameters. Using both synthetic and real data, we aim to better characterize the behavior of two HMM types that can be applied to electrophysiological data. Specifically, we study which differences in data features (such as frequency, amplitude, or signal-to-noise ratio) are more salient to the models and therefore more likely to drive the state decomposition. Overall, we aim at providing guidance for the appropriate use of this type of analysis on one- or two-channel neural electrophysiological data and an informed interpretation of its results given the characteristics of the data and the purpose of the analysis.NEW & NOTEWORTHY Compared with classical supervised methods, unsupervised methods of analysis have the advantage to be freer of subjective biases. However, it is not always clear what aspects of the data these methods are most sensitive to, which complicates interpretation. Focusing on the hidden Markov model, commonly used to describe electrophysiological data, we explore in detail the nature of its estimates through simulations and real data examples, providing important insights about what to expect from these models.

Mapping Interictal activity in epilepsy using a hidden Markov model: A magnetoencephalography study.

Epilepsy is a highly heterogeneous neurological disorder with variable etiology, manifestation, and response to treatment. It is imperative that new models of epileptiform brain activity account for this variability, to identify individual needs and allow clinicians to curate personalized care. Here, we use a hidden Markov model (HMM) to create a unique statistical model of interictal brain activity for 10 pediatric patients. We use magnetoencephalography (MEG) data acquired as part of standard clinical care for patients at the Children's Hospital of Philadelphia. These data are routinely analyzed using excess kurtosis mapping (EKM); however, as cases become more complex (extreme multifocal and/or polymorphic activity), they become harder to interpret with EKM. We assessed the performance of the HMM against EKM for three patient groups, with increasingly complicated presentation. The difference in localization of epileptogenic foci for the two methods was 7 ± 2 mm (mean ± SD over all 10 patients); and 94% ± 13% of EKM temporal markers were matched by an HMM state visit. The HMM localizes epileptogenic areas (in agreement with EKM) and provides additional information about the relationship between those areas. A key advantage over current methods is that the HMM is a data-driven model, so the output is tuned to each individual. Finally, the model output is intuitive, allowing a user (clinician) to review the result and manually select the HMM epileptiform state, offering multiple advantages over previous methods and allowing for broader implementation of MEG epileptiform analysis in surgical decision-making for patients with intractable epilepsy.

The relationship between frequency content and representational dynamics in the decoding of neurophysiological data.

Decoding of high temporal resolution, stimulus-evoked neurophysiological data is increasingly used to test theories about how the brain processes information. However, a fundamental relationship between the frequency spectra of the neural signal and the subsequent decoding accuracy timecourse is not widely recognised. We show that, in commonly used instantaneous signal decoding paradigms, each sinusoidal component of the evoked response is translated to double its original frequency in the subsequent decoding accuracy timecourses. We therefore recommend, where researchers use instantaneous signal decoding paradigms, that more aggressive low pass filtering is applied with a cut-off at one quarter of the sampling rate, to eliminate representational alias artefacts. However, this does not negate the accompanying interpretational challenges. We show that these can be resolved by decoding paradigms that utilise both a signal's instantaneous magnitude and its local gradient information as features for decoding. On a publicly available MEG dataset, this results in decoding accuracy metrics that are higher, more stable over time, and free of the technical and interpretational challenges previously characterised. We anticipate that a broader awareness of these fundamental relationships will enable stronger interpretations of decoding results by linking them more clearly to the underlying signal characteristics that drive them.

Effective psychological therapy for PTSD changes the dynamics of specific large-scale brain networks.

In posttraumatic stress disorder (PTSD), re-experiencing of the trauma is a hallmark symptom proposed to emerge from a de-contextualised trauma memory. Cognitive therapy for PTSD (CT-PTSD) addresses this de-contextualisation through different strategies. At the brain level, recent research suggests that the dynamics of specific large-scale brain networks play an essential role in both the healthy response to a threatening situation and the development of PTSD. However, very little is known about how these dynamics are altered in the disorder and rebalanced after treatment and successful recovery. Using a data-driven approach and fMRI, we detected recurring large-scale brain functional states with high temporal precision in a population of healthy trauma-exposed and PTSD participants before and after successful CT-PTSD. We estimated the total amount of time that each participant spent on each of the states while being exposed to trauma-related and neutral pictures. We found that PTSD participants spent less time on two default mode subnetworks involved in different forms of self-referential processing in contrast to PTSD participants after CT-PTSD (mtDMN+ and dmDMN+ ) and healthy trauma-exposed controls (only mtDMN+ ). Furthermore, re-experiencing severity was related to decreased time spent on the default mode subnetwork involved in contextualised retrieval of autobiographical memories, and increased time spent on the salience and visual networks. Overall, our results support the hypothesis that PTSD involves an imbalance in the dynamics of specific large-scale brain network states involved in self-referential processes and threat detection, and suggest that successful CT-PTSD might rebalance this dynamic aspect of brain function.

Dissociable Components of Information Encoding in Human Perception.

Brain decoding can predict visual perception from non-invasive electrophysiological data by combining information across multiple channels. However, decoding methods typically conflate the composite and distributed neural processes underlying perception that are together present in the signal, making it unclear what specific aspects of the neural computations involved in perception are reflected in this type of macroscale data. Using MEG data recorded while participants viewed a large number of naturalistic images, we analytically decomposed the brain signal into its oscillatory and non-oscillatory components, and used this decomposition to show that there are at least three dissociable stimulus-specific aspects to the brain data: a slow, non-oscillatory component, reflecting the temporally stable aspect of the stimulus representation; a global phase shift of the oscillation, reflecting the overall speed of processing of specific stimuli; and differential patterns of phase across channels, likely reflecting stimulus-specific computations. Further, we show that common cognitive interpretations of decoding analysis, in particular about how representations generalize across time, can benefit from acknowledging the multicomponent nature of the signal in the study of perception.

Hierarchical modelling of functional brain networks in population and individuals from big fMRI data.

A major goal of large-scale brain imaging datasets is to provide resources for investigating heterogeneous populations. Characterisation of functional brain networks for individual subjects from these datasets will have an enormous potential for prediction of cognitive or clinical traits. We propose for the first time a technique, Stochastic Probabilistic Functional Modes (sPROFUMO), that is scalable to UK Biobank (UKB) with expected 100,000 participants, and hierarchically estimates functional brain networks in individuals and the population, while allowing for bidirectional flow of information between the two. Using simulations, we show the model's utility, especially in scenarios that involve significant cross-subject variability, or require delineation of fine-grained differences between the networks. Subsequently, by applying the model to resting-state fMRI from 4999 UKB subjects, we mapped resting state networks (RSNs) in single subjects with greater detail than has been possible previously in UKB (>100 RSNs), and demonstrate that these RSNs can predict a range of sensorimotor and higher-level cognitive functions. Furthermore, we demonstrate several advantages of the model over independent component analysis combined with dual-regression (ICA-DR), particularly with respect to the estimation of the spatial configuration of the RSNs and the predictive power for cognitive traits. The proposed model and results can open a new door for future investigations into individualised profiles of brain function from big data.

Automatic decomposition of electrophysiological data into distinct nonsinusoidal oscillatory modes.

Neurophysiological signals are often noisy, nonsinusoidal, and consist of transient bursts. Extraction and analysis of oscillatory features (such as waveform shape and cross-frequency coupling) in such data sets remains difficult. This limits our understanding of brain dynamics and its functional importance. Here, we develop iterated masking empirical mode decomposition (itEMD), a method designed to decompose noisy and transient single-channel data into relevant oscillatory modes in a flexible, fully data-driven way without the need for manual tuning. Based on empirical mode decomposition (EMD), this technique can extract single-cycle waveform dynamics through phase-aligned instantaneous frequency. We test our method by extensive simulations across different noise, sparsity, and nonsinusoidality conditions. We find itEMD significantly improves the separation of data into distinct nonsinusoidal oscillatory components and robustly reproduces waveform shape across a wide range of relevant parameters. We further validate the technique on multimodal, multispecies electrophysiological data. Our itEMD extracts known rat hippocampal θ waveform asymmetry and identifies subject-specific human occipital α without any prior assumptions about the frequencies contained in the signal. Notably, it does so with significantly less mode mixing compared with existing EMD-based methods. By reducing mode mixing and simplifying interpretation of EMD results, itEMD will enable new analyses into functional roles of neural signals in behavior and disease.NEW & NOTEWORTHY We introduce a novel, data-driven method to identify oscillations in neural recordings. This approach is based on empirical mode decomposition and reduces mixing of components, one of its main problems. The technique is validated and compared with existing methods using simulations and real data. We show our method better extracts oscillations and their properties in highly noisy and nonsinusoidal datasets.

Within-cycle instantaneous frequency profiles report oscillatory waveform dynamics.

The nonsinusoidal waveform is emerging as an important feature of neuronal oscillations. However, the role of single-cycle shape dynamics in rapidly unfolding brain activity remains unclear. Here, we develop an analytical framework that isolates oscillatory signals from time series using masked empirical mode decomposition to quantify dynamical changes in the shape of individual cycles (along with amplitude, frequency, and phase) with instantaneous frequency. We show how phase-alignment, a process of projecting cycles into a regularly sampled phase grid space, makes it possible to compare cycles of different durations and shapes. "Normalized shapes" can then be constructed with high temporal detail while accounting for differences in both duration and amplitude. We find that the instantaneous frequency tracks nonsinusoidal shapes in both simulated and real data. Notably, in local field potential recordings of mouse hippocampal CA1, we find that theta oscillations have a stereotyped slow-descending slope in the cycle-wise average yet exhibit high variability on a cycle-by-cycle basis. We show how principal component analysis allows identification of motifs of theta cycle waveform that have distinct associations to cycle amplitude, cycle duration, and animal movement speed. By allowing investigation into oscillation shape at high temporal resolution, this analytical framework will open new lines of inquiry into how neuronal oscillations support moment-by-moment information processing and integration in brain networks.NEW & NOTEWORTHY We propose a novel analysis approach quantifying nonsinusoidal waveform shape. The approach isolates oscillations with empirical mode decomposition before waveform shape is quantified using phase-aligned instantaneous frequency. This characterizes the full shape profile of individual cycles while accounting for between-cycle differences in duration, amplitude, and timing. We validated in simulations before applying to identify a range of data-driven nonsinusoidal shape motifs in hippocampal theta oscillations.

Optimising network modelling methods for fMRI.

A major goal of neuroimaging studies is to develop predictive models to analyze the relationship between whole brain functional connectivity patterns and behavioural traits. However, there is no single widely-accepted standard pipeline for analyzing functional connectivity. The common procedure for designing functional connectivity based predictive models entails three main steps: parcellating the brain, estimating the interaction between defined parcels, and lastly, using these integrated associations between brain parcels as features fed to a classifier for predicting non-imaging variables e.g., behavioural traits, demographics, emotional measures, etc. There are also additional considerations when using correlation-based measures of functional connectivity, resulting in three supplementary steps: utilising Riemannian geometry tangent space parameterization to preserve the geometry of functional connectivity; penalizing the connectivity estimates with shrinkage approaches to handle challenges related to short time-series (and noisy) data; and removing confounding variables from brain-behaviour data. These six steps are contingent on each-other, and to optimise a general framework one should ideally examine these various methods simultaneously. In this paper, we investigated strengths and short-comings, both independently and jointly, of the following measures: parcellation techniques of four kinds (categorized further depending upon number of parcels), five measures of functional connectivity, the decision of staying in the ambient space of connectivity matrices or in tangent space, the choice of applying shrinkage estimators, six alternative techniques for handling confounds and finally four novel classifiers/predictors. For performance evaluation, we have selected two of the largest datasets, UK Biobank and the Human Connectome Project resting state fMRI data, and have run more than 9000 different pipeline variants on a total of ∼14000 individuals to determine the optimum pipeline. For independent performance validation, we have run some best-performing pipeline variants on ABIDE and ACPI datasets (∼1000 subjects) to evaluate the generalisability of proposed network modelling methods.

Post-stimulus beta responses are modulated by task duration.

Modulation of beta-band neural oscillations during and following movement is a robust marker of brain function. In particular, the post-movement beta rebound (PMBR), which occurs on movement cessation, has been related to inhibition and connectivity in the healthy brain, and is perturbed in disease. However, to realise the potential of the PMBR as a biomarker, its modulation by task parameters must be characterised and its functional role determined. Here, we used MEG to image brain electrophysiology during and after a grip-force task, with the aim to characterise how task duration, in the form of an isometric contraction, modulates beta responses. Fourteen participants exerted a 30% maximum voluntary grip-force for 2, 5 and 10 s. Our results showed that the amplitude of the PMBR is modulated by task duration, with increasing duration significantly reducing PMBR amplitude and increasing its time-to-peak. No variation in the amplitude of the movement related beta decrease (MRBD) with task duration was observed. To gain insight into what may underlie these trial-averaged results, we used a Hidden Markov Model to identify the individual trial dynamics of a brain network encompassing bilateral sensorimotor areas. The rapidly evolving dynamics of this network demonstrated similar variation with task parameters to the 'classical' rebound, and we show that the modulation of the PMBR can be well-described in terms of increased frequency of beta events on a millisecond timescale rather than modulation of beta amplitude during this time period. Our results add to the emerging picture that, in the case of a carefully controlled paradigm, beta modulation can be systematically controlled by task parameters and such control can reveal new information as to the processes that generate the average beta timecourse. These findings will support design of clinically relevant paradigms and analysis pipelines in future use of the PMBR as a marker of neuropathology.

Modelling subject variability in the spatial and temporal characteristics of functional modes.

Recent work has highlighted the scale and ubiquity of subject variability in observations from functional MRI data (fMRI). Furthermore, it is highly likely that errors in the estimation of either the spatial presentation of, or the coupling between, functional regions can confound cross-subject analyses, making accurate and unbiased representations of functional data essential for interpreting any downstream analyses. Here, we extend the framework of probabilistic functional modes (PFMs) (Harrison et al., 2015) to capture cross-subject variability not only in the mode spatial maps, but also in the functional coupling between modes and in mode amplitudes. A new implementation of the inference now also allows for the analysis of modern, large-scale data sets, and the combined inference and analysis package, PROFUMO, is available from git.fmrib.ox.ac.uk/samh/profumo. A new implementation of the inference now also allows for the analysis of modern, large-scale data sets. Using simulated data, resting-state data from 1000 subjects collected as part of the Human Connectome Project (Van Essen et al., 2013), and an analysis of 14 subjects in a variety of continuous task-states (Kieliba et al., 2019), we demonstrate how PFMs are able to capture, within a single model, a rich description of how the spatio-temporal structure of resting-state fMRI activity varies across subjects. We also compare the new PFM model to the well established independent component analysis with dual regression (ICA-DR) pipeline. This reveals that, under PFM assumptions, much more of the (behaviorally relevant) cross-subject variability in fMRI activity should be attributed to the variability in spatial maps, and that, after accounting for this, functional coupling between modes primarily reflects current cognitive state. This has fundamental implications for the interpretation of cross-sectional studies of functional connectivity that do not capture cross-subject variability to the same extent as PFMs.

The role of transient spectral 'bursts' in functional connectivity: A magnetoencephalography study.

Neural oscillations dominate electrophysiological measures of macroscopic brain activity and fluctuations in these rhythms offer an insightful window on cortical excitation, inhibition, and connectivity. However, in recent years the 'classical' picture of smoothly varying oscillations has been challenged by the idea that many 'oscillations' may actually be formed from the recurrence of punctate high-amplitude bursts in activity, whose spectral composition intersects the traditionally defined frequency ranges (e.g. alpha/beta band). This finding offers a new interpretation of measurable brain activity, however neither the methodological means to detect bursts, nor their link to other findings (e.g. connectivity) have been settled. Here, we use a new approach to detect bursts in magnetoencephalography (MEG) data. We show that a time-delay embedded Hidden Markov Model (HMM) can be used to delineate single-region bursts which are in agreement with existing techniques. However, unlike existing techniques, the HMM looks for specific spectral patterns in timecourse data. We characterise the distribution of burst duration, frequency of occurrence and amplitude across the cortex in resting state MEG data. During a motor task we show how the movement related beta decrease and post movement beta rebound are driven by changes in burst occurrence. Finally, we show that the beta band functional connectome can be derived using a simple measure of burst overlap, and that coincident bursts in separate regions correspond to a period of heightened coherence. In summary, this paper offers a new methodology for burst identification and connectivity analysis which will be important for future investigations of neural oscillations.

Temporally Unconstrained Decoding Reveals Consistent but Time-Varying Stages of Stimulus Processing.

In this article, we propose a method to track trial-specific neural dynamics of stimulus processing and decision making with high temporal precision. By applying this novel method to a perceptual template-matching task, we tracked representational brain states associated with the cascade of neural processing, from early sensory areas to higher order areas that are involved in integration and decision making. We address a major limitation of the traditional decoding approach: that it relies on consistent timing of these processes over trials. Using a TUDA approach, we found that the timing of the cognitive processes involved in perceptual judgments can vary considerably over trials. This revealed that the sequence of processing states was consistent for all subjects and trials, even when the timing of these states varied. Furthermore, we found that the specific timing of states on each trial was related to the quality of performance over trials. Altogether, this work not only highlights the serious pitfalls and misleading interpretations that result from assuming stimulus processing to be synchronous across trials but can also open important avenues to investigate learning and quantify plasticity.

Brain network dynamics are hierarchically organized in time.

The brain recruits neuronal populations in a temporally coordinated manner in task and at rest. However, the extent to which large-scale networks exhibit their own organized temporal dynamics is unclear. We use an approach designed to find repeating network patterns in whole-brain resting fMRI data, where networks are defined as graphs of interacting brain areas. We find that the transitions between networks are nonrandom, with certain networks more likely to occur after others. Further, this nonrandom sequencing is itself hierarchically organized, revealing two distinct sets of networks, or metastates, that the brain has a tendency to cycle within. One metastate is associated with sensory and motor regions, and the other involves areas related to higher order cognition. Moreover, we find that the proportion of time that a subject spends in each brain network and metastate is a consistent subject-specific measure, is heritable, and shows a significant relationship with cognitive traits.

Do the posterior midline cortices belong to the electrophysiological default-mode network?

The default-mode network (DMN) and its principal core hubs in the posterior midline cortices (PMC), i.e., the precuneus and the posterior cingulate cortex, play a critical role in the human brain structural and functional architecture. Because of their centrality, they are affected by a wide spectrum of brain disorders, e.g., Alzheimer's disease. Non-invasive electrophysiological techniques such as magnetoencephalography (MEG) are crucial to the investigation of the neurophysiology of the DMN and its alteration by brain disorders. However, MEG studies relying on band-limited power envelope correlation diverge in their ability to identify the PMC as a part of the DMN in healthy subjects at rest. Since these works were based on different MEG recording systems and different source reconstruction pipelines, we compared DMN functional connectivity estimated with two distinct MEG systems (Elekta, now MEGIN, and CTF) and two widely used reconstruction algorithms (Minimum Norm Estimation and linearly constrained minimum variance Beamformer). Our results identified the reconstruction method as the critical factor influencing PMC functional connectivity, which was significantly dampened by the Beamformer. On this basis, we recommend that future electrophysiological studies on the DMN should rely on Minimum Norm Estimation (or close variants) rather than on the classical Beamformer. Crucially, based on analytic knowledge about these two reconstruction algorithms, we demonstrated with simulations that this empirical observation could be explained by the existence of a spontaneous linear, approximately zero-lag synchronization structure between areas of the DMN or among multiple sources within the PMC. This finding highlights a novel property of the neural dynamics and functional architecture of a core human brain network at rest.

Tracking dynamic brain networks using high temporal resolution MEG measures of functional connectivity.

Fluctuations in functional interactions between brain regions typically occur at the millisecond time scale. Conventional connectivity metrics are not adequately time-resolved to detect such fast fluctuations in functional connectivity. At the same time, attempts to use conventional metrics in a time-resolved manner usually come with the selection of sliding windows of fixed arbitrary length. In the current work, we evaluated the use of high temporal resolution metrics of functional connectivity in conjunction with non-negative tensor factorisation to detect fast fluctuations in connectivity and temporally evolving subnetworks. To this end, we used the phase difference derivative, wavelet coherence, and we also introduced a new metric, the instantaneous amplitude correlation. In order to deal with the inherently noisy nature of magnetoencephalography data and large datasets, we make use of recurrence plots and we used pair-wise orthogonalisation to avoid spurious estimates of functional connectivity due to signal leakage. Firstly, metrics were evaluated in the context of dynamically coupled neural mass models in the presence and absence of delays and also compared to conventional static metrics with fixed sliding windows. Simulations showed that these high temporal resolution metrics outperformed conventional static connectivity metrics. Secondly, the sensitivity of the metrics to fluctuations in connectivity was analysed in post-movement beta rebound magnetoencephalography data, which showed time locked sensorimotor subnetworks that modulated with the post-movement beta rebound. Finally, sensitivity of the metrics was evaluated in resting-state magnetoencephalography, showing similar spatial patterns across metrics, thereby indicating the robustness of the current analysis. The current methods can be applied in cognitive experiments that involve fast modulations in connectivity in relation to cognition. In addition, these methods could also be used as input to temporal graph analysis to further characterise the rapid fluctuation in brain network topology.

Stable between-subject statistical inference from unstable within-subject functional connectivity estimates.

Spatial or temporal aspects of neural organization are known to be important indices of how cognition is organized. However, measurements and estimations are often noisy and many of the algorithms used are probabilistic, which in combination have been argued to limit studies exploring the neural basis of specific aspects of cognition. Focusing on static and dynamic functional connectivity estimations, we propose to leverage this variability to improve statistical efficiency in relating these estimations to behavior. To achieve this goal, we use a procedure based on permutation testing that provides a way of combining the results from many individual tests that refer to the same hypothesis. This is needed when testing a measure whose value is obtained from a noisy process, which can be repeated multiple times, referred to as replications. Focusing on functional connectivity, this noisy process can be: (a) computational, for example, when using an approximate inference algorithm for which different runs can produce different results or (b) observational, if we have the capacity to acquire data multiple times, and the different acquired data sets can be considered noisy examples of some underlying truth. In both cases, we are not interested in the individual replications but on the unobserved process generating each replication. In this note, we show how results can be combined instead of choosing just one of the estimated models. Using both simulations and real data, we show the benefits of this approach in practice.

A biophysical model of dynamic balancing of excitation and inhibition in fast oscillatory large-scale networks.

Over long timescales, neuronal dynamics can be robust to quite large perturbations, such as changes in white matter connectivity and grey matter structure through processes including learning, aging, development and certain disease processes. One possible explanation is that robust dynamics are facilitated by homeostatic mechanisms that can dynamically rebalance brain networks. In this study, we simulate a cortical brain network using the Wilson-Cowan neural mass model with conduction delays and noise, and use inhibitory synaptic plasticity (ISP) to dynamically achieve a spatially local balance between excitation and inhibition. Using MEG data from 55 subjects we find that ISP enables us to simultaneously achieve high correlation with multiple measures of functional connectivity, including amplitude envelope correlation and phase locking. Further, we find that ISP successfully achieves local E/I balance, and can consistently predict the functional connectivity computed from real MEG data, for a much wider range of model parameters than is possible with a model without ISP.

Unpacking Transient Event Dynamics in Electrophysiological Power Spectra.

Electrophysiological recordings of neuronal activity show spontaneous and task-dependent changes in their frequency-domain power spectra. These changes are conventionally interpreted as modulations in the amplitude of underlying oscillations. However, this overlooks the possibility of underlying transient spectral 'bursts' or events whose dynamics can map to changes in trial-average spectral power in numerous ways. Under this emerging perspective, a key challenge is to perform burst detection, i.e. to characterise single-trial transient spectral events, in a principled manner. Here, we describe how transient spectral events can be operationalised and estimated using Hidden Markov Models (HMMs). The HMM overcomes a number of the limitations of the standard amplitude-thresholding approach to burst detection; in that it is able to concurrently detect different types of bursts, each with distinct spectral content, without the need to predefine frequency bands of interest, and does so with less dependence on a priori threshold specification. We describe how the HMM can be used for burst detection and illustrate its benefits on simulated data. Finally, we apply this method to empirical data to detect multiple burst types in a task-MEG dataset, and illustrate how we can compute burst metrics, such as the task-evoked timecourse of burst duration.