Dissecting Transient Burst Events
Authors: Catharina Zich, Andrew J. Quinn, Lydia C. Mardell, Nick S. Ward, Sven Bestmann
Journal: Trends in Cognitive Sciences (2020)

Increasing efforts are being made to understand the role of intermittent, transient, high-power burst events of neural activity. These events have a temporal, spectral, and spatial domain. Here, we argue that considering all three domains is crucial to fully reveal the functional relevance of these events in health and disease. Increasing efforts are being made to understand the role of intermittent, transient, high-power burst events of neural activity. These events have a temporal, spectral, and spatial domain. Here, we argue that considering all three domains is crucial to fully reveal the functional relevance of these events in health and disease. Neural activity recorded from the scalp [magneto-/electroencephalogram (M/EEG)], cortical surface [electrocorticogram (EcoG)], or inside the brain [local field potential (LFP)] is traditionally analysed by averaging tens to hundreds of trials. The trial-wise average of spectrograms commonly shows periods of sustained high or low spectral power. The underlying single trial activity may manifest as oscillations (i.e., sustained rhythmic fluctuations of synchronous spiking activity) with periods of high or low amplitude. However, as recently demonstrated, sustained high spectral power in the averaged spectrogram can also arise from the accumulation of burst-like events across trial [1.Shin H. et al.The rate of transient beta frequency events predicts behavior across tasks and species.Elife. 2017; 6: 1-31Crossref Scopus (112) Google Scholar, 2.Tinkhauser G. et al.The modulatory effect of adaptive deep brain stimulation on beta bursts in Parkinson's disease.Brain. 2017; 140: 1053-1067Crossref PubMed Scopus (221) Google Scholar, 3.Feingold J. et al.Bursts of beta oscillation differentiate postperformance activity in the striatum and motor cortex of monkeys performing movement tasks.Proc. Natl. Acad. Sci. U. S. A. 2015; 112: 13687-13692Crossref PubMed Scopus (184) Google Scholar, 4.Lundqvist M. et al.Gamma and beta bursts underlie working memory.Neuron. 2016; 90: 152-164Abstract Full Text Full Text PDF PubMed Scopus (356) Google Scholar, 5.Sherman M.A. et al.Neural mechanisms of transient neocortical beta rhythms: converging evidence from humans, computational modeling, monkeys, and mice.Proc. Natl. Acad. Sci. U. S. A. 2016; 113: E4885-E4894Crossref PubMed Scopus (200) Google Scholar, 6.Jones S.R. When brain rhythms aren't 'rhythmic': implication for their mechanisms and meaning.Curr. Opin. Neurobiol. 2016; 40: 72-80Crossref PubMed Scopus (123) Google Scholar]. Burst-like events are intermittent, transient periods of synchronous spiking activity, the generator of which may or may not be rhythmic [7.van Ede F. et al.Neural oscillations: sustained rhythms or transient burst-events?.Trends Neurosci. 2018; 41: 415-417Abstract Full Text Full Text PDF PubMed Scopus (74) Google Scholar]. Accordingly, the underlying mechanism of high or low spectral power in the averaged spectrogram can be due to differences in event amplitude, or other event characteristics. Describing events provides an untapped opportunity to expand our understanding of brain function in health and disease. Events have a temporal, spectral, and spatial domain. However, most work so far has focussed on temporal event characteristics, such as event amplitude (traditionally seen as a temporal event characteristic), event duration, or event interval time (Figure 1A ). Events are thereby compressed to a singular dimension, by reducing the spatial and spectral domains of their underlying signals. However, this domain reduction removes potentially relevant aspects of the data. The spatial domain is commonly reduced by extracting the time series from a single spatial location, or summarising the time series of several spatial locations through their average or a linear combination [e.g., first principle component (PC)]. Similarly, the spectral domain is often reduced by selecting a single peak frequency, averaging the signal or amplitude envelope within a specified frequency band, or by fitting a state-wise frequency profile from a time-delay embedded or autoregressive-based Hidden Markov Model (HMM) [8.Quinn A.J. et al.Unpacking transient event dynamics in electrophysiological power spectra.Brain Topogr. 2019; 32: 1020-1034Crossref PubMed Scopus (20) Google Scholar]. Rarely, both the temporal and spectral characteristics of events are analysed together (Figure 1B [1.Shin H. et al.The rate of transient beta frequency events predicts behavior across tasks and species.Elife. 2017; 6: 1-31Crossref Scopus (112) Google Scholar]), but the spatial characteristics (Figure 1C) have not yet been incorporated. Given that differences in spectral power can be caused by changes in a variety of event characteristics (e.g., event duration and/or increase in event spatial width; Figure 1), consideration of all three domains is likely to be necessary for disclosing the underlying mechanisms of differences in population activity measures, such as spectral power. Domain reduction can have several unwanted consequences. First, different domain reduction methods likely result in differences in event characteristics. Focussing on the temporal domain, Figure 2A,B illustrates how event duration, event interval time, and event onset differ across common domain reduction approaches. Importantly, these differences are not attributable to differences in how events are identified (e.g., threshold). Furthermore, domain reduction can obscure interactions between domains, which may lead to ambiguous inferences and conclusions drawn from these analyses. For example, if events change their location in the spatial and/or spectral domain over the course of the event, domain reduction can lead to misestimation of the event duration. The second event in Figure 2B illustrates this point. As evident from Figure 2Biii, the event has a temporospatial gradient. Therefore, reducing the spatial domain can result in an underestimation of the actual event duration or missing the event altogether. Thus, it becomes apparent that domain reduction can hide important information about mass neural signals. Therefore, even when only one domain is of interest, considering all domains is most likely necessary to accurately describe each individual domain. This increases sensitivity to detect differences between conditions, individuals, or groups, and to minimise spurious inference. Considering all domains simultaneously yields a more detailed description of the underlying signal and allows interactions across domains to be characterised. These interactions can have many forms, and here we illustrate four such interactions between domains (Figure 2C) with a more detailed discussion of the interactions depicted in Figure 2Cii,iii. As shown in Figure 2Cii,iii, the spatial location of the event shifts from medial to lateral over the course of the event, while at the same time the upper frequency boundary of the event shifts to a lower frequency. Such a temporal–spectral–spatial interaction appears plausible given previous work outside the burst literature. For example, temporal–spatial interactions (travelling waves) have been described at similar temporal and spatial scales as burst events [9.Muller L. et al.Cortical travelling waves: mechanisms and computational principles.Nat. Rev. Neurosci. 2018; 19: 255-268Crossref PubMed Scopus (194) Google Scholar], and sensorimotor beta activity travels in a consistent anatomical direction [10.Rubino D. et al.Propagating waves mediate information transfer in the motor cortex.Nat. Neurosci. 2006; 9: 1549-1557Crossref PubMed Scopus (308) Google Scholar]. As another example, β peak frequency varies within different locations of the sensorimotor network (i.e., spectral–spatial interaction [11.Salmelin R. Hari R. Spatiotemporal characteristics of rhythmic neuromagnetic activity related to thumb movement.Neuroscience. 1994; 60: 537-550Crossref PubMed Scopus (616) Google Scholar]). These examples support the idea of interactions across the three domains. They illustrate how multidomain analysis can reveal gradients that would otherwise remain obscure but that are relevant for revealing the mechanistic role of events. Analysis of all three domains (with their five dimensions: temporal–spectral–spatial [x–y–z]) benefits from high signal-to-noise recordings of neural activity, which have become more readily accessible through recent advances in acquisition techniques of mass neural activity, such as large-array, high-density ECoG, high-precision MEG, or optically pumped magnetometers. Parallel development in signal processing, particularly effective spatial leakage correction, precise beamformer, high-resolution time-frequency analysis, and nonstationary approaches for event identification constitute additional advances for analysing events precisely across domains. These advances notwithstanding, challenges, such as varying signal-to-noise across time points, frequencies, and spatial locations, tools for analysing multidimensional event characteristics effectively across events and individuals, and, for M/EEG data, accuracy of spatial localisation (forward and inverse modelling), remain. Furthermore, the resolution of the temporal and spectral domains constitutes a compromise and depends on the choice of 'epoch' length in the time–frequency analysis. In summary, intermittent, transient, events have a temporal, spectral, and spatial domain. Considering all three domains will yield a more accurate description of individual event characteristics and their interactions. This is likely to underpin identification of more precise fingerprints of mechanisms underlying differences in population activity measures, such as spectral power. In turn, this can guide the development of rational and mechanistically grounded treatments that target neurological and psychiatric conditions. While we use sensorimotor β activity here as a showcase for discussing the impact of domain reduction and opportunities to consider all three domains, these points likely extend to other brain areas, frequency ranges, and pathologies. C.Z. is supported by the Brain Research UK (201718-13, 201617-03). A.J.Q. is supported by the NIHR Oxford Health Biomedical Research Centre and a Wellcome Trust Strategic Award (098369/Z/12/Z). L.C.M. is supported by the Medical Research Council (MR/N013867/1). The Wellcome Centre for Human Neuroimaging and The Wellcome Centre for Integrative Neuroimaging are supported by core funding from the Wellcome Trust (203147/Z/16/Z and 203139/Z/16/Z, respectively).