Short timescale abnormalities in the states of spontaneous synchrony in the functional neural networks in Alzheimer's disease

Comparison between EEG and MEG of static and dynamic resting-state networks

The characterisation of resting-state networks (RSNs) using neuroimaging techniques has significantly contributed to our understanding of the organisation of brain activity. Prior work has demonstrated the electrophysiological basis of RSNs and their dynamic nature, revealing transient activations of brain networks with millisecond timescales. While previous research has confirmed the comparability of RSNs identified by electroencephalography (EEG) to those identified by magnetoencephalography (MEG) and functional magnetic resonance imaging (fMRI), most studies have utilised static analysis techniques, ignoring the dynamic nature of brain activity. Often, these studies use high-density EEG systems, which limit their applicability in clinical settings. Addressing these gaps, our research studies RSNs using medium-density EEG systems (61 sensors), comparing both static and dynamic brain network features to those obtained from a high-density MEG system (306 sensors). We assess the qualitative and quantitative comparability of EEG-derived RSNs to those from MEG, including their ability to capture age-related effects, and explore the reproducibility of dynamic RSNs within and across the modalities. Our findings suggest that both MEG and EEG offer comparable static and dynamic network descriptions, albeit with MEG offering some increased sensitivity and reproducibility. Such RSNs and their comparability across the two modalities remained consistent qualitatively but not quantitatively when the data were reconstructed without subject-specific structural MRI images.

Dynamic functional connectivity and gene expression correlates in temporal lobe epilepsy: insights from hidden markov models

BACKGROUD

Temporal lobe epilepsy (TLE) is associated with abnormal dynamic functional connectivity patterns, but the dynamic changes in brain activity at each time point remain unclear, as does the potential molecular mechanisms associated with the dynamic temporal characteristics of TLE.

METHODS

Resting-state functional magnetic resonance imaging (rs-fMRI) was acquired for 84 TLE patients and 35 healthy controls (HCs). The data was then used to conduct HMM analysis on rs-fMRI data from TLE patients and an HC group in order to explore the intricate temporal dynamics of brain activity in TLE patients with cognitive impairment (TLE-CI). Additionally, we aim to examine the gene expression profiles associated with the dynamic modular characteristics in TLE patients using the Allen Human Brain Atlas (AHBA) database.

RESULTS

Five HMM states were identified in this study. Compared with HCs, TLE and TLE-CI patients exhibited distinct changes in dynamics, including fractional occupancy, lifetimes, mean dwell time and switch rate. Furthermore, transition probability across HMM states were significantly different between TLE and TLE-CI patients (p < 0.05). The temporal reconfiguration of states in TLE and TLE-CI patients was associated with several brain networks (including the high-order default mode network (DMN), subcortical network (SCN), and cerebellum network (CN). Furthermore, a total of 1580 genes were revealed to be significantly associated with dynamic brain states of TLE, mainly enriched in neuronal signaling and synaptic function.

CONCLUSIONS

This study provides new insights into characterizing dynamic neural activity in TLE. The brain network dynamics defined by HMM analysis may deepen our understanding of the neurobiological underpinnings of TLE and TLE-CI, indicating a linkage between neural configuration and gene expression in TLE.

Data-driven MEG analysis to extract fMRI resting-state networks

The electrophysiological basis of resting-state networks (RSN) is still under debate. In particular, no principled mechanism has been determined that is capable of explaining all RSN equally well. While magnetoencephalography (MEG) and electroencephalography are the methods of choice to determine the electrophysiological basis of RSN, no standard analysis pipeline of RSN yet exists. In this article, we compare the two main existing data-driven analysis strategies for extracting RSNs from MEG data and introduce a third approach. The first approach uses phase-amplitude coupling to determine the RSN. The second approach extracts RSN through an independent component analysis of the Hilbert envelope in different frequency bands, while the third new approach uses a singular value decomposition instead. To evaluate these approaches, we compare the MEG-RSN to the functional magnetic resonance imaging (fMRI)-RSN from the same subjects. Overall, it was possible to extract RSN with MEG using all three techniques, which matched the group-specific fMRI-RSN. Interestingly the new approach based on SVD yielded significantly higher correspondence to five out of seven fMRI-RSN than the two existing approaches. Importantly, with this approach, all networks-except for the visual network-had the highest correspondence to the fMRI networks within one frequency band. Thereby we provide further insights into the electrophysiological underpinnings of the fMRI-RSNs. This knowledge will be important for the analysis of the electrophysiological connectome.

osl-dynamics, a toolbox for modeling fast dynamic brain activity

Neural activity contains rich spatiotemporal structure that corresponds to cognition. This includes oscillatory bursting and dynamic activity that span across networks of brain regions, all of which can occur on timescales of tens of milliseconds. While these processes can be accessed through brain recordings and imaging, modeling them presents methodological challenges due to their fast and transient nature. Furthermore, the exact timing and duration of interesting cognitive events are often a priori unknown. Here, we present the OHBA Software Library Dynamics Toolbox (osl-dynamics), a Python-based package that can identify and describe recurrent dynamics in functional neuroimaging data on timescales as fast as tens of milliseconds. At its core are machine learning generative models that are able to adapt to the data and learn the timing, as well as the spatial and spectral characteristics, of brain activity with few assumptions. osl-dynamics incorporates state-of-the-art approaches that can be, and have been, used to elucidate brain dynamics in a wide range of data types, including magneto/electroencephalography, functional magnetic resonance imaging, invasive local field potential recordings, and electrocorticography. It also provides novel summary measures of brain dynamics that can be used to inform our understanding of cognition, behavior, and disease. We hope osl-dynamics will further our understanding of brain function, through its ability to enhance the modeling of fast dynamic processes.

Post-stroke upper limb recovery is correlated with dynamic resting-state network connectivity

Motor recovery is still limited for people with stroke especially those with greater functional impairments. In order to improve outcome, we need to understand more about the mechanisms underpinning recovery. Task-unbiased, blood flow-independent post-stroke neural activity can be acquired from resting brain electrophysiological recordings and offers substantial promise to investigate physiological mechanisms, but behaviourally relevant features of resting-state sensorimotor network dynamics have not yet been identified. Thirty-seven people with subcortical ischaemic stroke and unilateral hand paresis of any degree were longitudinally evaluated at 3 weeks (early subacute) and 12 weeks (late subacute) after stroke. Resting-state magnetoencephalography and clinical scores of motor function were recorded and compared with matched controls. Magnetoencephalography data were decomposed using a data-driven hidden Markov model into 10 time-varying resting-state networks. People with stroke showed statistically significantly improved Action Research Arm Test and Fugl-Meyer upper extremity scores between 3 weeks and 12 weeks after stroke (both P < 0.001). Hidden Markov model analysis revealed a primarily alpha-band ipsilesional resting-state sensorimotor network which had a significantly increased life-time (the average time elapsed between entering and exiting the network) and fractional occupancy (the occupied percentage among all networks) at 3 weeks after stroke when compared with controls. The life-time of the ipsilesional resting-state sensorimotor network positively correlated with concurrent motor scores in people with stroke who had not fully recovered. Specifically, this relationship was observed only in ipsilesional rather in contralesional sensorimotor network, default mode network or visual network. The ipsilesional sensorimotor network metrics were not significantly different from controls at 12 weeks after stroke. The increased recruitment of alpha-band ipsilesional resting-state sensorimotor network at subacute stroke served as functionally correlated biomarkers exclusively in people with stroke with not fully recovered hand paresis, plausibly reflecting functional motor recovery processes.

Impaired long-range excitatory time scale predicts abnormal neural oscillations and cognitive deficits in Alzheimer's disease

Alzheimer's disease (AD) is the most common form of dementia, progressively impairing memory and cognition. While neuroimaging studies have revealed functional abnormalities in AD, how these relate to aberrant neuronal circuit mechanisms remains unclear. Using magnetoencephalography imaging we documented abnormal local neural synchrony patterns in patients with AD. To identify abnormal biophysical mechanisms underlying these abnormal electrophysiological patterns, we estimated the parameters of a spectral graph-theory model (SGM). SGM is an analytic model that describes how long-range fiber projections in the brain mediate the excitatory and inhibitory activity of local neuronal subpopulations. The long-range excitatory time scale was associated with greater deficits in global cognition and was able to distinguish AD patients from controls with high accuracy. These results demonstrate that long-range excitatory time scale of neuronal activity, despite being a global measure, is a key determinant in the spatiospectral signatures and cognition in AD.

Altered brain dynamic in major depressive disorder: state and trait features

Temporal neural synchrony disruption can be linked to a variety of symptoms of major depressive disorder (MDD), including mood rigidity and the inability to break the cycle of negative emotion or attention biases. This might imply that altered dynamic neural synchrony may play a role in the persistence and exacerbation of MDD symptoms. Our study aimed to investigate the changes in whole-brain dynamic patterns of the brain functional connectivity and activity related to depression using the hidden Markov model (HMM) on resting-state functional magnetic resonance imaging (rs-fMRI) data. We compared the patterns of brain functional dynamics in a large sample of 314 patients with MDD (65.9% female; age (mean ± standard deviation): 35.9 ± 13.4) and 498 healthy controls (59.4% female; age: 34.0 ± 12.8). The HMM model was used to explain variations in rs-fMRI functional connectivity and averaged functional activity across the whole-brain by using a set of six unique recurring states. This study compared the proportion of time spent in each state and the average duration of visits to each state to assess stability between different groups. Compared to healthy controls, patients with MDD showed significantly higher proportional time spent and temporal stability in a state characterized by weak functional connectivity within and between all brain networks and relatively strong averaged functional activity of regions located in the somatosensory motor (SMN), salience (SN), and dorsal attention (DAN) networks. Both proportional time spent and temporal stability of this brain state was significantly associated with depression severity. Healthy controls, in contrast to the MDD group, showed proportional time spent and temporal stability in a state with relatively strong functional connectivity within and between all brain networks but weak averaged functional activity across the whole brain. These findings suggest that disrupted brain functional synchrony across time is present in MDD and associated with current depression severity.

New Therapeutics in Alzheimer's Disease Longitudinal Cohort study (NTAD): study protocol

INTRODUCTION

With the pressing need to develop treatments that slow or stop the progression of Alzheimer's disease, new tools are needed to reduce clinical trial duration and validate new targets for human therapeutics. Such tools could be derived from neurophysiological measurements of disease.

METHODS AND ANALYSIS

The New Therapeutics in Alzheimer's Disease study (NTAD) aims to identify a biomarker set from magneto/electroencephalography that is sensitive to disease and progression over 1 year. The study will recruit 100 people with amyloid-positive mild cognitive impairment or early-stage Alzheimer's disease and 30 healthy controls aged between 50 and 85 years. Measurements of the clinical, cognitive and imaging data (magnetoencephalography, electroencephalography and MRI) of all participants will be taken at baseline. These measurements will be repeated after approximately 1 year on participants with Alzheimer's disease or mild cognitive impairment, and clinical and cognitive assessment of these participants will be repeated again after approximately 2 years. To assess reliability of magneto/electroencephalographic changes, a subset of 30 participants with mild cognitive impairment or early-stage Alzheimer's disease will also undergo repeat magneto/electroencephalography 2 weeks after baseline. Baseline and longitudinal changes in neurophysiology are the primary analyses of interest. Additional outputs will include atrophy and cognitive change and estimated numbers needed to treat each arm of simulated clinical trials of a future disease-modifying therapy.

ETHICS AND DATA STATEMENT

The study has received a favourable opinion from the East of England Cambridge Central Research Ethics Committee (REC reference 18/EE/0042). Results will be disseminated through internal reports, peer-reviewed scientific journals, conference presentations, website publication, submission to regulatory authorities and other publications. Data will be made available via the Dementias Platform UK Data Portal on completion of initial analyses by the NTAD study group.

A review of MEG dynamic brain network research

The dynamic description of neural networks has attracted the attention of researchers for dynamic networks may carry more information compared with resting-state networks. As a non-invasive electrophysiological data with high temporal and spatial resolution, magnetoencephalogram (MEG) can provide rich information for the analysis of dynamic functional brain networks. In this review, the development of MEG brain network was summarized. Several analysis methods such as sliding window, Hidden Markov model, and time-frequency based methods used in MEG dynamic brain network studies were discussed. Finally, the current research about multi-modal brain network analysis and their applications with MEG neurophysiology, which are prospected to be one of the research directions in the future, were concluded.

Resting state fast brain dynamics predict interindividual variability in motor performance

Motor learning features rapid enhancement during practice then offline post-practice gains with the reorganization of related brain networks. We hypothesised that fast transient, sub-second variations in magnetoencephalographic (MEG) network activity during the resting-state (RS) reflect early learning-related plasticity mechanisms and/or interindividual motor variability in performance. MEG RS activity was recorded before and 20 min after motor learning. Hidden Markov modelling (HMM) of MEG power envelope signals highlighted 8 recurrent topographical states. For two states, motor performance levels were associated with HMM temporal parameters both in pre- and post-learning resting-state sessions. However, no association emerged with offline changes in performance. These results suggest a trait-like relationship between spontaneous transient neural dynamics at rest and interindividual variations in motor abilities. On the other hand, transient RS dynamics seem not to be state-dependent, i.e., modulated by learning experience and reflect neural plasticity, at least on the short timescale.

Microstates and power envelope hidden Markov modeling probe bursting brain activity at different timescales

State modeling of whole-brain electroencephalography (EEG) or magnetoencephalography (MEG) allows to investigate transient, recurring neurodynamical events. Two widely-used techniques are the microstate analysis of EEG signals and hidden Markov modeling (HMM) of MEG power envelopes. Both reportedly lead to similar state lifetimes on the 100 ms timescale, suggesting a common neural basis. To investigate whether microstates and power envelope HMM states describe the same neural dynamics, we used simultaneous MEG/EEG recordings at rest and compared the spatial signature and temporal activation dynamics of microstates and power envelope HMM states obtained separately from EEG and MEG. Results showed that microstates and power envelope HMM states differ both spatially and temporally. Microstates reflect sharp events of neural synchronization, whereas power envelope HMM states disclose network-level activity with 100–200 ms lifetimes. Further, MEG microstates do not correspond to the canonical EEG microstates but are better interpreted as split HMM states. On the other hand, both MEG and EEG HMM states involve the (de)activation of similar functional networks. Microstate analysis and power envelope HMM thus appear sensitive to neural events occurring over different spatial and temporal scales. As such, they represent complementary approaches to explore the fast, sub-second scale bursting electrophysiological dynamics in spontaneous human brain activity.

Age of onset modulates resting-state brain network dynamics in Friedreich Ataxia

This magnetoencephalography (MEG) study addresses (i) how Friedreich ataxia (FRDA) affects the sub‐second dynamics of resting‐state brain networks, (ii) the main determinants of their dynamic alterations, and (iii) how these alterations are linked with FRDA‐related changes in resting‐state functional brain connectivity (rsFC) over long timescales. For that purpose, 5 min of resting‐state MEG activity were recorded in 16 FRDA patients (mean age: 27 years, range: 12–51 years; 10 females) and matched healthy subjects. Transient brain network dynamics was assessed using hidden Markov modeling (HMM). Post hoc median‐split, nonparametric permutations and Spearman rank correlations were used for statistics. In FRDA patients, a positive correlation was found between the age of symptoms onset (ASO) and the temporal dynamics of two HMM states involving the posterior default mode network (DMN) and the temporo‐parietal junctions (TPJ). FRDA patients with an ASO <11 years presented altered temporal dynamics of those two HMM states compared with FRDA patients with an ASO > 11 years or healthy subjects. The temporal dynamics of the DMN state also correlated with minute‐long DMN rsFC. This study demonstrates that ASO is the main determinant of alterations in the sub‐second dynamics of posterior associative neocortices in FRDA patients and substantiates a direct link between sub‐second network activity and functional brain integration over long timescales.

When the time is right: Temporal dynamics of brain activity in healthy aging and dementia

Brain activity and communications are complex phenomena that dynamically unfold over time. However, in contrast with the large number of studies reporting neuroanatomical differences in activation relative to young adults, changes of temporal dynamics of neural activity during normal and pathological aging have been grossly understudied and are still poorly known. Here, we synthesize the current state of knowledge from MEG and EEG studies that aimed at specifying the effects of healthy and pathological aging on local and network dynamics, and discuss the clinical and theoretical implications of these findings. We argue that considering the temporal dynamics of brain activations and networks could provide a better understanding of changes associated with healthy aging, and the progression of neurodegenerative disease. Recent research has also begun to shed light on the association of these dynamics with other imaging modalities and with individual differences in cognitive performance. These insights hold great potential for driving new theoretical frameworks and development of biomarkers to aid in identifying and treating age-related cognitive changes.

Abnormal meta-state activation of dynamic brain networks across the Alzheimer spectrum

The characterization of the distinct dynamic functional connectivity (dFC) patterns that activate in the brain during rest can help to understand the underlying time-varying network organization. The presence and behavior of these patterns (known as meta-states) have been widely studied by means of functional magnetic resonance imaging (fMRI). However, modalities with high-temporal resolution, such as electroencephalography (EEG), enable the characterization of fast temporally evolving meta-state sequences. Mild cognitive impairment (MCI) and dementia due to Alzheimer's disease (AD) have been shown to disrupt spatially localized activation and dFC between different brain regions, but not much is known about how they affect meta-state network topologies and their network dynamics. The main hypothesis of the study was that MCI and dementia due to AD alter normal meta-state sequences by inducing a loss of structure in their patterns and a reduction of their dynamics. Moreover, we expected that patients with MCI would display more flexible behavior compared to patients with dementia due to AD. Thus, the aim of the current study was twofold: (i) to find repeating, distinctly organized network patterns (meta-states) in neural activity; and (ii) to extract information about meta-state fluctuations and how they are influenced by MCI and dementia due to AD. To accomplish these goals, we present a novel methodology to characterize dynamic meta-states and their temporal fluctuations by capturing aspects based on both their discrete activation and the continuous evolution of their individual strength. These properties were extracted from 60-s resting-state EEG recordings from 67 patients with MCI due to AD, 50 patients with dementia due to AD, and 43 cognitively healthy controls. First, the instantaneous amplitude correlation (IAC) was used to estimate instantaneous functional connectivity with a high temporal resolution. We then extracted meta-states by means of graph community detection based on recurrence plots (RPs), both at the individual- and group-level. Subsequently, a diverse set of properties of the continuous and discrete fluctuation patterns of the meta-states was extracted and analyzed. The main novelty of the methodology lies in the usage of Louvain GJA community detection to extract meta-states from IAC-derived RPs and the extended analysis of their discrete and continuous activation. Our findings showed that distinct dynamic functional connectivity meta-states can be found on the EEG time-scale, and that these were not affected by the oscillatory slowing induced by MCI or dementia due to AD. However, both conditions displayed a loss of meta-state modularity, coupled with shorter dwell times and higher complexity of the meta-state sequences. Furthermore, we found evidence that meta-state sequencing is not entirely random; it shows an underlying structure that is partially lost in MCI and dementia due to AD. These results show evidence that AD progression is associated with alterations in meta-state switching, and a degradation of dynamic brain flexibility.

Alterations in resting-state network dynamics along the Alzheimer's disease continuum

Human brain activity is intrinsically organized into resting-state networks (RSNs) that transiently activate or deactivate at the sub-second timescale. Few neuroimaging studies have addressed how Alzheimer's disease (AD) affects these fast temporal brain dynamics, and how they relate to the cognitive, structural and metabolic abnormalities characterizing AD. We aimed at closing this gap by investigating both brain structure and function using magnetoencephalography (MEG) and hybrid positron emission tomography-magnetic resonance (PET/MR) in 10 healthy elders, 10 patients with subjective cognitive decline (SCD), 10 patients with amnestic mild cognitive impairment (aMCI) and 10 patients with typical Alzheimer's disease with dementia (AD). The fast activation/deactivation state dynamics of RSNs were assessed using hidden Markov modeling (HMM) of power envelope fluctuations at rest measured with MEG. Correlations were sought between temporal properties of HMM states and participants' cognitive test scores, whole hippocampal grey matter volume and regional brain glucose metabolism. The posterior default-mode network (DMN) was less often activated and for shorter durations in AD patients than matched healthy elders. No significant difference was found in patients with SCD or aMCI. The time spent by participants in the activated posterior DMN state did not correlate significantly with cognitive scores, nor with the whole hippocampal volume. However, it correlated positively with the regional glucose consumption in the right dorsolateral prefrontal cortex (DLPFC). AD patients present alterations of posterior DMN power activation dynamics at rest that identify an additional electrophysiological correlate of AD-related synaptic and neural dysfunction. The right DLPFC may play a causal role in the activation of the posterior DMN, possibly linked to the occurrence of mind wandering episodes. As such, these data might suggest a neural correlate of the decrease in mind wandering episodes reported in pathological aging.

Changes in electrophysiological static and dynamic human brain functional architecture from childhood to late adulthood

This magnetoencephalography study aimed at characterizing age-related changes in resting-state functional brain organization from mid-childhood to late adulthood. We investigated neuromagnetic brain activity at rest in 105 participants divided into three age groups: children (6-9 years), young adults (18-34 years) and healthy elders (53-78 years). The effects of age on static resting-state functional brain integration were assessed using band-limited power envelope correlation, whereas those on transient functional brain dynamics were disclosed using hidden Markov modeling of power envelope activity. Brain development from childhood to adulthood came with (1) a strengthening of functional integration within and between resting-state networks and (2) an increased temporal stability of transient (100-300 ms lifetime) and recurrent states of network activation or deactivation mainly encompassing lateral or medial associative neocortical areas. Healthy aging was characterized by decreased static resting-state functional integration and dynamic stability within the primary visual network. These results based on electrophysiological measurements free of neurovascular biases suggest that functional brain integration mainly evolves during brain development, with limited changes in healthy aging. These novel electrophysiological insights into human brain functional architecture across the lifespan pave the way for future clinical studies investigating how brain disorders affect brain development or healthy aging.

Tau pathology in early Alzheimer's disease is linked to selective disruptions in neurophysiological network dynamics

Understanding the role of Tau protein aggregation in the pathogenesis of Alzheimer's disease is critical for the development of new Tau-based therapeutic strategies to slow or prevent dementia. We tested the hypothesis that Tau pathology is associated with functional organization of widespread neurophysiological networks. We used electro-magnetoencephalography with [¹⁸F]AV-1451 PET scanning to quantify Tau-dependent network changes. Using a graph theoretical approach to brain connectivity, we quantified nodal measures of functional segregation, centrality, and the efficiency of information transfer and tested them against levels of [¹⁸F]AV-1451. Higher Tau burden in early Alzheimer's disease was associated with a shift away from the optimal small-world organization and a more fragmented network in the beta and gamma bands, whereby parieto-occipital areas were disconnected from the anterior parts of the network. Similarly, higher Tau burden was associated with decreases in both local and global efficiency, especially in the gamma band. The results support the translational development of neurophysiological "signatures" of Alzheimer's disease, to understand disease mechanisms in humans and facilitate experimental medicine studies.

Cortical electrophysiological evidence for individual-specific temporal organization of brain functional networks

The human brain has been demonstrated to rapidly and continuously form and dissolve networks on a subsecond timescale, offering effective and essential substrates for cognitive processes. Understanding how the dynamic organization of brain functional networks on a subsecond level varies across individuals is, therefore, of great interest for personalized neuroscience. However, it remains unclear whether features of such rapid network organization are reliably unique and stable in single subjects and, therefore, can be used in characterizing individual networks. Here, we used two sets of 5‐min magnetoencephalography (MEG) resting data from 39 healthy subjects over two consecutive days and modeled the spontaneous brain activity as recurring networks fast shifting between each other in a coordinated manner. MEG cortical maps were obtained through source reconstruction using the beamformer method and subjects' temporal structure of recurring networks was obtained via the Hidden Markov Model. Individual organization of dynamic brain activity was quantified with the features of the network‐switching pattern (i.e., transition probability and mean interval time) and the time‐allocation mode (i.e., fractional occupancy and mean lifetime). Using these features, we were able to identify subjects from the group with significant accuracies (~40% on average in 0.5–48 Hz). Notably, the default mode network displayed a distinguishable pattern, being the least frequently visited network with the longest duration for each visit. Together, we provide initial evidence suggesting that the rapid dynamic temporal organization of brain networks achieved in electrophysiology is intrinsic and subject specific.

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.