Search

We describe OHBA Software Library for the analysis of electrophysiology data (osl-ephys). This toolbox builds on top of the widely used MNE-Python package and provides unique analysis tools for magneto−/electro-encephalography (M/EEG) sensor and source space analysis, which can be used modularly. In particular, it facilitates processing large amounts of data using batch parallel processing, with high standards for reproducibility through a config API and log keeping, and efficient quality assurance by producing HTML processing reports. It also provides new functionality for doing coregistration, source reconstruction and parcellation in volumetric space, allowing for an alternative pipeline that avoids the need for surface-based processing, e.g., through the use of Fieldtrip. Here, we introduce osl-ephys by presenting examples applied to a publicly available M/EEG data (the multimodal faces dataset). osl-ephys is open-source software distributed on the Apache License and available as a Python package through PyPi and GitHub.

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.

Despite advances in the field of dynamic connectivity, fixed sliding window approaches for the detection of fluctuations in functional connectivity are still widely used. The use of conventional connectivity metrics in conjunction with a fixed sliding window comes with the arbitrariness of the chosen window lengths. In this paper we use multivariate autoregressive and neural mass models with a priori defined ground truths to systematically analyze the sensitivity of conventional metrics in combination with different window lengths to detect genuine fluctuations in connectivity for various underlying state durations. Metrics of interest are the coherence, imaginary coherence, phase lag index, phase locking value and the amplitude envelope correlation. We performed analysis for two nodes and at the network level. We demonstrate that these metrics show indeed higher variability for genuine temporal fluctuations in connectivity compared to a static connectivity state superimposed by noise. Overall, the error of the connectivity estimates themselves decreases for longer state durations (order of seconds), while correlations of the connectivity fluctuations with the ground truth was higher for longer state durations. In general, metrics, in combination with a sliding window, perform poorly for very short state durations. Increasing the SNR of the system only leads to a moderate improvement. In addition, at the network level, only longer window widths were sufficient to detect plausible resting state networks that matched the underlying ground truth, especially for the phase locking value, amplitude envelope correlation and coherence. The length of these longer window widths did not necessarily correspond to the underlying state durations. For short window widths resting state network connectivity patterns could not be retrieved. We conclude that fixed sliding window approaches for connectivity can detect modulations of connectivity, but mostly if the underlying dynamics operate on moderate to slow timescales. In practice, this can be a drawback, as state durations can vary significantly in empirical data.

Complex thought and behavior arise through dynamic recruitment of large-scale brain networks. The signatures of this process may be observable in electrophysiological data; yet robust modeling of rapidly changing functional network structure on rapid cognitive timescales remains a considerable challenge. Here, we present one potential solution using Hidden Markov Models (HMMs), which are able to identify brain states characterized by engaging distinct functional networks that reoccur over time. We show how the HMM can be inferred on continuous, parcellated source-space Magnetoencephalography (MEG) task data in an unsupervised manner, without any knowledge of the task timings. We apply this to a freely available MEG dataset in which participants completed a face perception task, and reveal task-dependent HMM states that represent whole-brain dynamic networks transiently bursting at millisecond time scales as cognition unfolds. The analysis pipeline demonstrates a general way in which the HMM can be used to do a statistically valid whole-brain, group-level task analysis on MEG task data, which could be readily adapted to a wide range of task-based studies.