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.

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.

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.

Evaluating functional brain organization in individuals and identifying contributions to network overlap.

Individual differences in the spatial organization of resting state networks have received increased attention in recent years. Measures of individual-specific spatial organization of brain networks and overlapping network organization have been linked to important behavioral and clinical traits and are therefore potential biomarker targets for personalized psychiatry approaches. To better understand individual-specific spatial brain organization, this paper addressed three key goals. First, we determined whether it is possible to reliably estimate weighted (non-binarized) resting state network maps using data from only a single individual, while also maintaining maximum spatial correspondence across individuals. Second, we determined the degree of spatial overlap between distinct networks, using test-retest and twin data. Third, we systematically tested multiple hypotheses (spatial mixing, temporal switching, and coupling) as candidate explanations for why networks overlap spatially. To estimate weighted network organization, we adopt the Probabilistic Functional Modes (PROFUMO) algorithm, which implements a Bayesian framework with hemodynamic and connectivity priors to supplement optimization for spatial sparsity/independence. Our findings showed that replicable individual-specific estimates of weighted resting state networks can be derived using high quality fMRI data within individual subjects. Network organization estimates using only data from each individual subject closely resembled group-informed network estimates (which was not explicitly modeled in our individual-specific analyses), suggesting that cross-subject correspondence was largely maintained. Furthermore, our results confirmed the presence of spatial overlap in network organization, which was replicable across sessions within individuals and in monozygotic twin pairs. Intriguingly, our findings provide evidence that network overlap is indicative of linear additive coupling. These results suggest that regions of network overlap concurrently process information from all contributing networks, potentially pointing to the role of overlapping network organization in the integration of information across multiple brain systems.

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.

Reduced coupling between offline neural replay events and default mode network activation in schizophrenia.

Schizophrenia is characterized by an abnormal resting state and default mode network brain activity. However, despite intense study, the mechanisms linking default mode network dynamics to neural computation remain elusive. During rest, sequential hippocampal reactivations, known as 'replay', are played out within default mode network activation windows, highlighting a potential role of replay-default mode network coupling in memory consolidation and model-based mental simulation. Here, we test a hypothesis of reduced replay-default mode network coupling in schizophrenia, using magnetoencephalography and a non-spatial sequence learning task designed to elicit off-task (i.e. resting state) neural replay. Participants with a diagnosis of schizophrenia (n = 28, mean age 28.2 years, range 20-40, 6 females, 13 not taking antipsychotic medication) and non-clinical control participants (n = 29, mean age 28.1 years, range 18-45, 6 females, matched at group level for age, intelligence quotient, gender, years in education and working memory) underwent a magnetoencephalography scan both during task completion and during a post-task resting state session. We used neural decoding to infer the time course of default mode network activation (time-delay embedding hidden Markov model) and spontaneous neural replay (temporally delayed linear modelling) in resting state magnetoencephalography data. Using multiple regression, we then quantified the extent to which default mode network activation was uniquely predicted by replay events that recapitulated the learned task sequences (i.e. 'task-relevant' replay-default mode network coupling). In control participants, replay-default mode network coupling was augmented following sequence learning, an augmentation that was specific for replay of task-relevant (i.e. learned) state transitions. This task-relevant replay-default mode network coupling effect was significantly reduced in schizophrenia (t(52) = 3.93, P = 0.018). Task-relevant replay-default mode network coupling predicted memory maintenance of learned sequences (ρ(52) = 0.31, P = 0.02). Importantly, reduced task-relevant replay-default mode network coupling in schizophrenia was not explained by differential replay or altered default mode network dynamics between groups nor by reference to antipsychotic exposure. Finally, task-relevant replay-default mode network coupling during rest correlated with stimulus-evoked default mode network modulation as measured in a separate task session. In the context of a proposed functional role of replay-default mode network coupling, our findings shed light on the functional significance of default mode network abnormalities in schizophrenia and provide for a consilience between task-based and resting state default mode network findings in this disorder.

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.

Understanding Harmonic Structures Through Instantaneous Frequency.

The analysis of harmonics and non-sinusoidal waveform shape in time-series data is growing in importance. However, a precise definition of what constitutes a harmonic is lacking. In this paper, we propose a rigorous definition of when to consider signals to be in a harmonic relationship based on an integer frequency ratio, constant phase, and a well-defined joint instantaneous frequency. We show this definition is linked to extrema counting and Empirical Mode Decomposition (EMD). We explore the mathematics of our definition and link it to results from analytic number theory. This naturally leads to us to define two classes of harmonic structures, termed strong and weak, with different extrema behaviour. We validate our framework using both simulations and real data. Specifically, we look at the harmonic structures in shallow water waves, the FitzHugh-Nagumo neuronal model, and the non-sinusoidal theta oscillation in rat hippocampus local field potential data. We further discuss how our definition helps to address mode splitting in nonlinear time-series decomposition methods. A clear understanding of when harmonics are present in signals will enable a deeper understanding of the functional roles of non-sinusoidal oscillations.

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.

Transient beta activity and cortico-muscular connectivity during sustained motor behaviour.

Neural oscillations are thought to play a central role in orchestrating activity states between distant neural populations. For example, during isometric contraction, 13-30 Hz beta activity becomes phase coupled between the motor cortex and the contralateral muscle. This and related observations have led to the proposal that beta activity and connectivity sustain stable cognitive and motor states - or the 'status quo' - in the brain. Recently, however, beta activity at the single-trial level has been shown to be short-lived - though so far this has been reported for regional beta activity in tasks without sustained motor demands. Here, we measured magnetoencephalography (MEG) and electromyography (EMG) in 18 human participants performing a sustained isometric contraction (gripping) task. If cortico-muscular beta connectivity is directly responsible for sustaining a stable motor state, then beta activity within single trials should be (or become) sustained in this context. In contrast, we found that motor beta activity and connectivity with the downstream muscle were transient. Moreover, we found that sustained motor requirements did not prolong beta-event duration in comparison to rest. These findings suggest that neural synchronisation between the brain and the muscle involves short 'bursts' of frequency-specific connectivity, even when task demands - and motor behaviour - are sustained.

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.

Multi-dynamic modelling reveals strongly time-varying resting fMRI correlations.

The activity of functional brain networks is responsible for the emergence of time-varying cognition and behaviour. Accordingly, time-varying correlations (Functional Connectivity) in resting fMRI have been shown to be predictive of behavioural traits, and psychiatric and neurological conditions. Typically, methods that measure time varying Functional Connectivity (FC), such as sliding windows approaches, do not separately model when changes occur in the mean activity levels from when changes occur in the FC, therefore conflating these two distinct types of modulation. We show that this can bias the estimation of time-varying FC to appear more stable over time than it actually is. Here, we propose an alternative approach that models changes in the mean brain activity and in the FC as being able to occur at different times to each other. We refer to this method as the Multi-dynamic Adversarial Generator Encoder (MAGE) model, which includes a model of the network dynamics that captures long-range time dependencies, and is estimated on fMRI data using principles of Generative Adversarial Networks. We evaluated the approach across several simulation studies and resting fMRI data from the Human Connectome Project (1003 subjects), as well as from UK Biobank (13301 subjects). Importantly, we find that separating fluctuations in the mean activity levels from those in the FC reveals much stronger changes in FC over time, and is a better predictor of individual behavioural variability.

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.

Mild traumatic brain injury impairs the coordination of intrinsic and motor-related neural dynamics.

Mild traumatic brain injury (mTBI) poses a considerable burden on healthcare systems. Whilst most patients recover quickly, a significant number suffer from sequelae that are not accompanied by measurable structural damage. Understanding the neural underpinnings of these debilitating effects and developing a means to detect injury, would address an important unmet clinical need. It could inform interventions and help predict prognosis. Magnetoencephalography (MEG) affords excellent sensitivity in probing neural function and presents significant promise for assessing mTBI, with abnormal neural oscillations being a potential specific biomarker. However, growing evidence suggests that neural dynamics are (at least in part) driven by transient, pan-spectral bursting and in this paper, we employ this model to investigate mTBI. We applied a Hidden Markov Model to MEG data recorded during resting state and a motor task and show that previous findings of diminished intrinsic beta amplitude in individuals with mTBI are largely due to the reduced beta band spectral content of bursts, and that diminished beta connectivity results from a loss in the temporal coincidence of burst states. In a motor task, mTBI results in diminished burst amplitude, altered modulation of burst probability during movement, and a loss in connectivity in motor networks. These results suggest that, mechanistically, mTBI disrupts the structural framework underlying neural synchrony, which impairs network function. Whilst the damage may be too subtle for structural imaging to see, the functional consequences are detectable and persist after injury. Our work shows that mTBI impairs the dynamic coordination of neural network activity and proposes a potent new method for understanding mTBI.

Revealing the Dynamic Nature of Amplitude Modulated Neural Entrainment With Holo-Hilbert Spectral Analysis.

Patterns in external sensory stimuli can rapidly entrain neuronally generated oscillations observed in electrophysiological data. Here, we manipulated the temporal dynamics of visual stimuli with cross-frequency coupling (CFC) characteristics to generate steady-state visual evoked potentials (SSVEPs). Although CFC plays a pivotal role in neural communication, some cases reporting CFC may be false positives due to non-sinusoidal oscillations that can generate artificially inflated coupling values. Additionally, temporal characteristics of dynamic and non-linear neural oscillations cannot be fully derived with conventional Fourier-based analyses mainly due to trade off of temporal resolution for frequency precision. In an attempt to resolve these limitations of linear analytical methods, Holo-Hilbert Spectral Analysis (HHSA) was investigated as a potential approach for examination of non-linear and non-stationary CFC dynamics in this study. Results from both simulation and SSVEPs demonstrated that temporal dynamic and non-linear CFC features can be revealed with HHSA. Specifically, the results of simulation showed that the HHSA is less affected by the non-sinusoidal oscillation and showed possible cross frequency interactions embedded in the simulation without any a priori assumptions. In the SSVEPs, we found that the time-varying cross-frequency interaction and the bidirectional coupling between delta and alpha/beta bands can be observed using HHSA, confirming dynamic physiological signatures of neural entrainment related to cross-frequency coupling. These findings not only validate the efficacy of the HHSA in revealing the natural characteristics of signals, but also shed new light on further applications in analysis of brain electrophysiological data with the aim of understanding the functional roles of neuronal oscillation in various cognitive functions.

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.

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.

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.

Temporally delayed linear modelling (TDLM) measures replay in both animals and humans.

There are rich structures in off-task neural activity which are hypothesized to reflect fundamental computations across a broad spectrum of cognitive functions. Here, we develop an analysis toolkit - temporal delayed linear modelling (TDLM) - for analysing such activity. TDLM is a domain-general method for finding neural sequences that respect a pre-specified transition graph. It combines nonlinear classification and linear temporal modelling to test for statistical regularities in sequences of task-related reactivations. TDLM is developed on the non-invasive neuroimaging data and is designed to take care of confounds and maximize sequence detection ability. Notably, as a linear framework, TDLM can be easily extended, without loss of generality, to capture rodent replay in electrophysiology, including in continuous spaces, as well as addressing second-order inference questions, for example, its temporal and spatial varying pattern. We hope TDLM will advance a deeper understanding of neural computation and promote a richer convergence between animal and human neuroscience.

Adapting non-invasive human recordings along multiple task-axes shows unfolding of spontaneous and over-trained choice.

Choices rely on a transformation of sensory inputs into motor responses. Using invasive single neuron recordings, the evolution of a choice process has been tracked by projecting population neural responses into state spaces. Here, we develop an approach that allows us to recover similar trajectories on a millisecond timescale in non-invasive human recordings. We selectively suppress activity related to three task-axes, relevant and irrelevant sensory inputs and response direction, in magnetoencephalography data acquired during context-dependent choices. Recordings from premotor cortex show a progression from processing sensory input to processing the response. In contrast to previous macaque recordings, information related to choice-irrelevant features is represented more weakly than choice-relevant sensory information. To test whether this mechanistic difference between species is caused by extensive over-training common in non-human primate studies, we trained humans on >20,000 trials of the task. Choice-irrelevant features were still weaker than relevant features in premotor cortex after over-training.