EMD: Empirical Mode Decomposition and Hilbert-Huang Spectral Analyses in Python

Comparison of measurement of integrated relaxation pressure by esophageal manometry with analysis of swallowing sounds with artificial intelligence in patients with achalasia

BACKGROUND

Esophageal motility disorders are mainly evaluated with high-resolution manometry (HRM) which is a time-consuming and uncomfortable procedure with potential adverse events. Acoustic characterization of the swallowing has the potential to be an alternative noninvasive procedure.

METHODS

We compared the findings on HRM and swallowing sounds in 43 patients who were referred for evaluation of dysphagia. The sound analysis was done with empirical mode decomposition method and with artificial intelligence (AI) and the estimated integrated relaxation pressure (IRP) from a two-layer neural network method was compared to measured IRP on HRM. The model then was tested in five patients.

KEY RESULTS

IRP was estimated with high accuracy using the model developed with two-layer neural network method.

CONCLUSIONS & INFERENCES

The analysis of acoustic properties of swallowing has the potential to be used for evaluation of esophageal motility disorders, this needs to be further evaluated in larger studies.

Age-related differences in how the shape of alpha and beta oscillations change during reaction time tasks

While the shape of cortical oscillations is increasingly recognised to be physiologically and functionally informative, its relevance to the aging motor system has not been established. We therefore examined the shape of alpha and beta band oscillations recorded at rest, as well as during performance of simple and go/no-go reaction time tasks, in 33 young (23.3 ± 2.9 years, 27 females) and 27 older (60.0 ± 5.2 years, 23 females) adults. The shape of individual oscillatory cycles was characterised using a recently developed pipeline involving empirical mode decomposition, before being decomposed into waveform motifs using principal component analysis. This revealed four principal components that were uniquely influenced by task and/or age. These described specific dimensions of shape and tended to be modulated during the reaction phase of each task. Our results suggest that although oscillation shape is task-dependent, the nature of this effect is altered by advancing age, possibly reflecting alterations in cortical activity. These outcomes demonstrate the utility of this approach for understanding the neurophysiological effects of ageing.

Exploring the impact of tonal inventory on speech perception across languages: a study of MMN responses in tonal language speakers

Previous research on the perception of segmental features of languages has established a correlation between the phoneme inventory of a language and its speakers' perceptual abilities, as indexed by discrimination tasks and Mismatch Negativity (MMN). Building on this background, the current study elucidated the relationship between perceptual ability and tonal inventory by utilizing two tonal languages. Two groups of participants were included in the present experiment: Mandarin speakers and Hakka-Mandarin speakers. Onset latency analysis revealed a significant difference in the Mandarin syllable condition, with Hakka-Mandarin speakers demonstrating earlier MMN latency than Mandarin speakers. This suggests a more efficient auditory processing mechanism in Hakka-Mandarin speakers. Both groups, however, showed similar MMN latency in the Hakka syllable condition. The interaction between language background and syllable type indicates that other factors, such as syllable sonority, also influence MMN responses. These findings highlight the importance of considering multiple phonemic inventories and syllable characteristics in studies of tonal perception.

Cell-substrate distance fluctuations of confluent cells enable fast and coherent collective migration

Collective cell migration is an emergent phenomenon, with long-range cell-cell communication influenced by various factors, including transmission of forces, viscoelasticity of individual cells, substrate interactions, and mechanotransduction. We investigate how alterations in cell-substrate distance fluctuations, cell-substrate adhesion, and traction forces impact the average velocity and temporal-spatial correlation of confluent monolayers formed by either wild-type (WT) MDCKII cells or zonula occludens (ZO)-1/2-depleted MDCKII cells (double knockdown [dKD]) representing highly contractile cells. The data indicate that confluent dKD monolayers exhibit decreased average velocity compared to less contractile WT cells concomitant with increased substrate adhesion, reduced traction forces, a more compact shape, diminished cell-cell interactions, and reduced cell-substrate distance fluctuations. Depletion of basal actin and myosin further supports the notion that short-range cell-substrate interactions, particularly fluctuations driven by basal actomyosin, significantly influence the migration speed of the monolayer on a larger length scale.

Dopaminergic versus anticholinergic treatment effects on physiologic complexity of hand tremor in Parkinson's disease: A randomized crossover study

AIMS

Parkinsonian tremor (PT) is regulated by numerous neurophysiological components across multiple temporospatial scales. The dynamics of tremor fluctuation are thus highly complex. This study aimed to explore the effects of different medications on tremor complexity, and how the underlying factors contribute to such tremor complexity.

METHODS

In this study, 66 participants received a 2-mg dose of benzhexol or a pre-determined dose of levodopa at two study visits in a randomized order. Before and after taking the medications, tremor fluctuation was recorded using surface electromyography electrodes and accelerometers in resting, posture, and weighting conditions with and without a concurrent cognitive task. Tremor complexity was quantified using multiscale entropy.

RESULTS

Tremor complexity in resting (p = 0.002) and postural condition (p < 0.0001) was lower when participants were performing a cognitive task compared to a task-free condition. After taking levodopa and benzhexol, participants had increased (p = 0.02-0.03) and decreased (p = 0.03) tremor complexity compared to pre-medication state, respectively. Tremor complexity and its changes as induced by medications were significantly correlated with clinical ratings and their changes (β = -0.23 to -0.39; p = 0.002-0.04), respectively.

CONCLUSION

Tremor complexity may be a promising marker to capture the pathophysiology underlying the development of PT, aiding the characterization of the effects medications have on PT regulation.

Waveform-based classification of dentate spikes

Synchronous excitatory discharges from the entorhinal cortex (EC) to the dentate gyrus (DG) generate fast and prominent patterns in the hilar local field potential (LFP), called dentate spikes (DSs). As sharp-wave ripples in CA1, DSs are more likely to occur in quiet behavioral states, when memory consolidation is thought to take place. However, their functions in mnemonic processes are yet to be elucidated. The classification of DSs into types 1 or 2 is determined by their origin in the lateral or medial EC, as revealed by current source density (CSD) analysis, which requires recordings from linear probes with multiple electrodes spanning the DG layers. To allow the investigation of the functional role of each DS type in recordings obtained from single electrodes and tetrodes, which are abundant in the field, we developed an unsupervised method using Gaussian mixture models to classify such events based on their waveforms. Our classification approach achieved high accuracies (> 80%) when validated in 8 mice with DG laminar profiles. The average CSDs, waveforms, rates, and widths of the DS types obtained through our method closely resembled those derived from the CSD-based classification. As an example of application, we used the technique to analyze single-electrode LFPs from apolipoprotein (apo) E3 and apoE4 knock-in mice. We observed that the latter group, which is a model for Alzheimer's disease, exhibited wider DSs of both types from a young age, with a larger effect size for DS type 2, likely reflecting early pathophysiological alterations in the EC-DG network, such as hyperactivity. In addition to the applicability of the method in expanding the study of DS types, our results show that their waveforms carry information about their origins, suggesting different underlying network dynamics and roles in memory processing.

Coordinating brain-distributed network activities in memory resistant to extinction

Certain memories resist extinction to continue invigorating maladaptive actions. The robustness of these memories could depend on their widely distributed implementation across populations of neurons in multiple brain regions. However, how dispersed neuronal activities are collectively organized to underpin a persistent memory-guided behavior remains unknown. To investigate this, we simultaneously monitored the prefrontal cortex, nucleus accumbens, amygdala, hippocampus, and ventral tegmental area (VTA) of the mouse brain from initial recall to post-extinction renewal of a memory involving cocaine experience. We uncover a higher-order pattern of short-lived beta-frequency (15-25 Hz) activities that are transiently coordinated across these networks during memory retrieval. The output of a divergent pathway from upstream VTA glutamatergic neurons, paced by a slower (4-Hz) oscillation, actuates this multi-network beta-band coactivation; its closed-loop phase-informed suppression prevents renewal of cocaine-biased behavior. Binding brain-distributed neural activities in this temporally structured manner may constitute an organizational principle of robust memory expression.

Interdependence of "What" and "When" in the Brain

From a brain's-eye-view, when a stimulus occurs and what it is are interrelated aspects of interpreting the perceptual world. Yet in practice, the putative perceptual inferences about sensory content and timing are often dichotomized and not investigated as an integrated process. We here argue that neural temporal dynamics can influence what is perceived, and in turn, stimulus content can influence the time at which perception is achieved. This computational principle results from the highly interdependent relationship of what and when in the environment. Both brain processes and perceptual events display strong temporal variability that is not always modeled; we argue that understanding-and, minimally, modeling-this temporal variability is key for theories of how the brain generates unified and consistent neural representations and that we ignore temporal variability in our analysis practice at the peril of both data interpretation and theory-building. Here, we review what and when interactions in the brain, demonstrate via simulations how temporal variability can result in misguided interpretations and conclusions, and outline how to integrate and synthesize what and when in theories and models of brain computation.

Waveform-based classification of dentate spikes

Synchronous excitatory discharges from the entorhinal cortex (EC) to the dentate gyrus (DG) generate fast and prominent patterns in the hilar local field potential (LFP), called dentate spikes (DSs). As sharp-wave ripples in CA1, DSs are more likely to occur in quiet behavioral states, when memory consolidation is thought to take place. However, their functions in mnemonic processes are yet to be elucidated. The classification of DSs into types 1 or 2 is determined by their origin in the lateral or medial EC, as revealed by current source density (CSD) analysis, which requires recordings from linear probes with multiple electrodes spanning the DG layers. To allow the investigation of the functional role of each DS type in recordings obtained from single electrodes and tetrodes, which are abundant in the field, we developed an unsupervised method using Gaussian mixture models to classify such events based on their waveforms. Our classification approach achieved high accuracies (> 80%) when validated in 8 mice with DG laminar profiles. The average CSDs, waveforms, rates, and widths of the DS types obtained through our method closely resembled those derived from the CSD-based classification. As an example of application, we used the technique to analyze single-electrode LFPs from apolipoprotein (apo) E3 and apoE4 knock-in mice. We observed that the latter group, which is a model for Alzheimer's disease, exhibited wider DSs of both types from a young age, with a larger effect size for DS type 2, likely reflecting early pathophysiological alterations in the EC-DG network, such as hyperactivity. In addition to the applicability of the method in expanding the study of DS types, our results show that their waveforms carry information about their origins, suggesting different underlying network dynamics and roles in memory processing.

Microscopic defect dynamics during a brittle-to-ductile transition

Deformation of all materials necessitates the collective propagation of various microscopic defects. On Earth, fracturing gives way to crystal-plastic deformation with increasing depth resulting in a "brittle-to-ductile" transition (BDT) region that is key for estimating the integrated strength of tectonic plates, constraining the earthquake cycle, and utilizing deep geothermal resources. Here, we show that the crossing of a BDT in marble during deformation experiments in the laboratory is accompanied by systematic increase in the frequency of acoustic emissions suggesting a profound change in the mean size and propagation velocity of the active defects. We further identify dominant classes of emitted waveforms using unsupervised learning methods and show that their relative activity systematically changes as the rocks cross the brittle-ductile transition. As pressure increases, long-period signals are suppressed and short-period signals become dominant. At higher pressures, signals frequently come in avalanche-like patterns. We propose that these classes of waveforms correlate with individual dominant defect types. Complex mixed-mode events indicate that interactions between the defects are common over the whole pressure range, in agreement with postmortem microstructural observations. Our measurements provide unique, real-time data of microscale dynamics over a broad range of pressures (10 to 200 MPa) and can inform micromechanical models for semi-brittle deformation.

High resolution seafloor thermometry for internal wave and upwelling monitoring using Distributed Acoustic Sensing

Temperature is an essential oceanographic variable (EOV) that still today remains coarsely resolved below the surface and near the seafloor. Here, we gather evidence to confirm that Distributed Acoustic Sensing (DAS) technology can convert tens of kilometer-long seafloor fiber-optic telecommunication cables into dense arrays of temperature anomaly sensors having millikelvin (mK) sensitivity, thus allowing to monitor oceanic processes such as internal waves and upwelling with unprecedented detail. Notably, we report high-resolution observations of highly coherent near-inertial and super-inertial internal waves in the NW Mediterranean sea, offshore of Toulon, France, having spatial extents of a few kilometers and producing maximum thermal anomalies of more than 5 K at maximum absolute rates of more than 1 K/h. We validate our observations with in-situ oceanographic sensors and an alternative optical fiber sensing technology. Currently, DAS only provides temperature changes estimates, however practical solutions are outlined to obtain continuous absolute temperature measurements with DAS at the seafloor. Our observations grant key advantages to DAS over established temperature sensors, showing its transformative potential for the description of seafloor temperature fluctuations over an extended range of spatial and temporal scales, as well as for the understanding of the evolution of the ocean in a broad sense (e.g. physical and ecological). Diverse ocean-oriented fields could benefit from the potential applications of this fast-developing technology.

A holo-spectral EEG analysis provides an early detection of cognitive decline and predicts the progression to Alzheimer's disease

Aims

Our aim was to differentiate patients with mild cognitive impairment (MCI) and Alzheimer's disease (AD) from cognitively normal (CN) individuals and predict the progression from MCI to AD within a 3-year longitudinal follow-up. A newly developed Holo-Hilbert Spectral Analysis (HHSA) was applied to resting state EEG (rsEEG), and features were extracted and subjected to machine learning algorithms.

Methods

A total of 205 participants were recruited from three hospitals, with CN (n = 51, MMSE > 26), MCI (n = 42, CDR = 0.5, MMSE ≥ 25), AD1 (n = 61, CDR = 1, MMSE < 25), AD2 (n = 35, CDR = 2, MMSE < 16), and AD3 (n = 16, CDR = 3, MMSE < 16). rsEEG was also acquired from all subjects. Seventy-two MCI patients (CDR = 0.5) were longitudinally followed up with two rsEEG recordings within 3 years and further subdivided into an MCI-stable group (MCI-S, n = 36) and an MCI-converted group (MCI-C, n = 36). The HHSA was then applied to the rsEEG data, and features were extracted and subjected to machine-learning algorithms.

Results

(a) At the group level analysis, the HHSA contrast of MCI and different stages of AD showed augmented amplitude modulation (AM) power of lower-frequency oscillations (LFO; delta and theta bands) with attenuated AM power of higher-frequency oscillations (HFO; beta and gamma bands) compared with cognitively normal elderly controls. The alpha frequency oscillation showed augmented AM power across MCI to AD1 with a reverse trend at AD2. (b) At the individual level of cross-sectional analysis, implementation of machine learning algorithms discriminated between groups with good sensitivity (Sen) and specificity (Spec) as follows: CN elderly vs. MCI: 0.82 (Sen)/0.80 (Spec), CN vs. AD1: 0.94 (Sen)/0.80 (Spec), CN vs. AD2: 0.93 (Sen)/0.90 (Spec), and CN vs. AD3: 0.75 (Sen)/1.00 (Spec). (c) In the longitudinal MCI follow-up, the initial contrasted HHSA between MCI-S and MCI-C groups showed significantly attenuated AM power of alpha and beta band oscillations. (d) At the individual level analysis of longitudinal MCI groups, deploying machine learning algorithms with the best seven features resulted in a sensitivity of 0.9 by the support vector machine (SVM) classifier, with a specificity of 0.8 yielded by the decision tree classifier.

Conclusion

Integrating HHSA into EEG signals and machine learning algorithms can differentiate between CN and MCI as well as also predict AD progression at the MCI stage.

Improved HHT-microstate analysis of EEG in nicotine addicts

Background

Substance addiction is a chronic disease which causes great harm to modern society and individuals. At present, many studies have applied EEG analysis methods to the substance addiction detection and treatment. As a tool to describe the spatio-temporal dynamic characteristics of large-scale electrophysiological data, EEG microstate analysis has been widely used, which is an effective method to study the relationship between EEG electrodynamics and cognition or disease.

Methods

To study the difference of EEG microstate parameters of nicotine addicts at each frequency band, we combine an improved Hilbert Huang Transformation (HHT) decomposition with microstate analysis, which is applied to the EEG of nicotine addicts.

Results

After using improved HHT-Microstate method, we notice that there is significant difference in EEG microstates of nicotine addicts between viewing smoke pictures group (smoke) and viewing neutral pictures group (neutral). Firstly, there is a significant difference in EEG microstates at full-frequency band between smoke and neutral group. Compared with the FIR-Microstate method, the similarity index of microstate topographic maps at alpha and beta bands had significant differences between smoke and neutral group. Secondly, we find significant class × group interactions for microstate parameters at delta, alpha and beta bands. Finally, the microstate parameters at delta, alpha and beta bands obtained by the improved HHT-microstate analysis method are selected as features for classification and detection under the Gaussian kernel support vector machine. The highest accuracy is 92% sensitivity is 94% and specificity is 91%, which can more effectively detect and identify addiction diseases than FIR-Microstate and FIR-Riemann methods.

Conclusion

Thus, the improved HHT-Microstate analysis method can effectively identify substance addiction diseases and provide new ideas and insights for the brain research of nicotine addiction.

Reduction or enhancement? Repetition effects on early brain potentials during visual word recognition are frequency dependent

Most studies on word repetition have demonstrated that repeated stimuli yield reductions in brain activity. Despite the well-known repetition reduction effect, some literature reports repetition enhancements in electroencephalogram (EEG) activities. However, although studies of object and face recognition have consistently demonstrated both repetition reduction and enhancement effects, the results of repetition enhancement effects were not consistent in studies of visual word recognition. Therefore, the present study aimed to further investigate the repetition effect on the P200, an early event-related potential (ERP) component that indexes the coactivation of lexical candidates during visual word recognition. To achieve a high signal-to-noise ratio, EEG signals were decomposed into various modes by using the Hilbert-Huang transform. Results demonstrated a repetition enhancement effect on P200 activity in alpha-band oscillation and that lexicality and orthographic neighborhood size would influence the magnitude of the repetition enhancement effect on P200. These findings suggest that alpha activity during visual word recognition might reflect the coactivation of orthographically similar words in the early stages of lexical processing. Meantime, there were repetition reduction effects on ERP activities in theta-delta band oscillation, which might index that the lateral inhibition between lexical candidates would be omitted in repetition.

A Frequency Estimation Scheme Based on Gaussian Average Filtering Decomposition and Hilbert Transform: With Estimation of Respiratory Rate as an Example

Frequency estimation plays a critical role in vital sign monitoring. Methods based on Fourier transform and eigen-analysis are commonly adopted techniques for frequency estimation. Because of the nonstationary and time-varying characteristics of physiological processes, time-frequency analysis (TFA) is a feasible way to perform biomedical signal analysis. Among miscellaneous approaches, Hilbert-Huang transform (HHT) has been demonstrated to be a potential tool in biomedical applications. However, the problems of mode mixing, unnecessary redundant decomposition and boundary effect are the common deficits that occur during the procedure of empirical mode decomposition (EMD) or ensemble empirical mode decomposition (EEMD). The Gaussian average filtering decomposition (GAFD) technique has been shown to be appropriate in several biomedical scenarios and can be an alternative to EMD and EEMD. This research proposes the combination of GAFD and Hilbert transform that is termed the Hilbert-Gauss transform (HGT) to overcome the conventional drawbacks of HHT in TFA and frequency estimation. This new method is verified to be effective for the estimation of respiratory rate (RR) in finger photoplethysmography (PPG), wrist PPG and seismocardiogram (SCG). Compared with the ground truth values, the estimated RRs are evaluated to be of excellent reliability by intraclass correlation coefficient (ICC) and to be of high agreement by Bland-Altman analysis.

Low frequency ambient noise dynamics and trends in the Indian Ocean, Cape Leeuwin, Australia

Examination of 18 years of nearly continuous low frequency deep ocean ambient noise offshore Cape Leeuwin, Australia, finds evidence of a decreasing nonlinear trend suggestive of long-term cyclic dynamics. The nonlinear trend is found to be consistent with trends in oceanographic sea surface temperature, which are thought to drive changes in Antarctic sea ice extent. Assessment of oscillatory dynamics finds causal covariance between ambient noise levels and Indian Ocean sea surface temperature dipoles. Dynamics of annual ambient noise and Antarctic sea ice extent are examined suggesting a phase-locked relationship revealing nonlinear characteristics of the presumed dependence. Collectively, the hypotheses that deep water ambient noise dynamics in the Indian Ocean are influenced by Antarctic sea ice extent and melt dynamics and that linear models do not fully capture long-term ambient noise trends and dynamics are supported.

An Evaluation of Non-Contact Photoplethysmography-Based Methods for Remote Respiratory Rate Estimation

The respiration rate (RR) is one of the physiological signals deserving monitoring for assessing human health and emotional states. However, traditional devices, such as the respiration belt to be worn around the chest, are not always a feasible solution (e.g., telemedicine, device discomfort). Recently, novel approaches have been proposed aiming at estimating RR in a less invasive yet reliable way, requiring the acquisition and processing of contact or remote Photoplethysmography (contact reference and remote-PPG, respectively). The aim of this paper is to address the lack of systematic evaluation of proposed methods on publicly available datasets, which currently impedes a fair comparison among them. In particular, we evaluate two prominent families of PPG processing methods estimating Respiratory Induced Variations (RIVs): the first encompasses methods based on the direct extraction of morphological features concerning the RR; and the second group includes methods modeling respiratory artifacts adopting, in the most promising cases, single-channel blind source separation. Extensive experiments have been carried out on the public BP4D+ dataset, showing that the morphological estimation of RIVs is more reliable than those produced by a single-channel blind source separation method (both in contact and remote testing phases), as well as in comparison with a representative state-of-the-art Deep Learning-based approach for remote respiratory information estimation.

Opening Questions in Visual Working Memory

In this reflective piece on visual working memory, I depart from the laboriously honed skills of writing a review. Instead of integrating approaches, synthesizing evidence, and building a cohesive perspective, I scratch my head and share niggles and puzzlements. I expose where my scholarship and understanding are stumped by findings and standard views in the literature.

Characterizing the effect of demographics, cardiorespiratory factors, and inter-subject variation on maternal heart rate variability in pregnancy with statistical modeling: a retrospective observational analysis

Pregnancy complications are associated with insufficient adaptation of the maternal autonomic nervous system to the physiological demands of pregnancy. Consequently, assessing maternal heart rate variability (mHRV)-which reflects autonomic regulation-is a promising tool for detecting early deterioration in maternal health. However, before mHRV can be used to screen for complications, an understanding of the factors influencing mHRV during healthy pregnancy is needed. In this retrospective observational study, we develop regression models to unravel the effects of maternal demographics (age, body mass index (BMI), gestational age (GA), and parity), cardiorespiratory factors (heart rate and breathing rate), and inter-subject variation on mHRV. We develop these models using two datasets which are comprised of, respectively, single measurements in 290 healthy pregnant women and repeated measurements (median = 8) in 29 women with healthy pregnancies. Our most consequential finding is that between one-third and two-thirds of the variation in mHRV can be attributed to inter-subject variability. Additionally, median heart rate dominantly affects mHRV (p < 0.001), while BMI and parity have no effect. Moreover, we found that median breathing rate, age, and GA all impact mHRV (p < 0.05). These results suggest that personalized, long-term monitoring would be necessary for using mHRV for obstetric screening.

Adult-born dentate granule cells promote hippocampal population sparsity

The dentate gyrus (DG) gates neocortical information flow to the hippocampus. Intriguingly, the DG also produces adult-born dentate granule cells (abDGCs) throughout the lifespan, but their contribution to downstream firing dynamics remains unclear. Here, we show that abDGCs promote sparser hippocampal population spiking during mnemonic processing of novel stimuli. By combining triple-(DG-CA3-CA1) ensemble recordings and optogenetic interventions in behaving mice, we show that abDGCs constitute a subset of high-firing-rate neurons with enhanced activity responses to novelty and strong modulation by theta oscillations. Selectively activating abDGCs in their 4-7-week post-birth period increases sparsity of hippocampal population patterns, whereas suppressing abDGCs reduces this sparsity, increases principal cell firing rates and impairs novel object recognition with reduced dimensionality of the network firing structure, without affecting single-neuron spatial representations. We propose that adult-born granule cells transiently support sparser hippocampal population activity structure for higher-dimensional responses relevant to effective mnemonic information processing.

Robust estimation of 1/f activity improves oscillatory burst detection

Neural oscillations often occur as transient bursts with variable amplitude and frequency dynamics. Quantifying these effects is important for understanding brain-behaviour relationships, especially in continuous datasets. To robustly measure bursts, rhythmical periods of oscillatory activity must be separated from arrhythmical background 1/f activity, which is ubiquitous in electrophysiological recordings. The Better OSCillation (BOSC) framework achieves this by defining a power threshold above the estimated background 1/f activity, combined with a duration threshold. Here we introduce a modification to this approach called fBOSC, which uses a spectral parametrisation tool to accurately model background 1/f activity in neural data. fBOSC (which is openly available as a MATLAB toolbox) is robust to power spectra with oscillatory peaks and can also model non-linear spectra. Through a series of simulations, we show that fBOSC more accurately models the 1/f power spectrum compared with existing methods. fBOSC was especially beneficial where power spectra contained a 'knee' below ~.5-10 Hz, which is typical in neural data. We also found that, unlike other methods, fBOSC was unaffected by oscillatory peaks in the neural power spectrum. Moreover, by robustly modelling background 1/f activity, the sensitivity for detecting oscillatory bursts was standardised across frequencies (e.g., theta- and alpha-bands). Finally, using openly available resting state magnetoencephalography and intracranial electrophysiology datasets, we demonstrate the application of fBOSC for oscillatory burst detection in the theta-band. These simulations and empirical analyses highlight the value of fBOSC in detecting oscillatory bursts, including in datasets that are long and continuous with no distinct experimental trials.

Good scientific practice in EEG and MEG research: Progress and perspectives

Good scientific practice (GSP) refers to both explicit and implicit rules, recommendations, and guidelines that help scientists to produce work that is of the highest quality at any given time, and to efficiently share that work with the community for further scrutiny or utilization. For experimental research using magneto- and electroencephalography (MEEG), GSP includes specific standards and guidelines for technical competence, which are periodically updated and adapted to new findings. However, GSP also needs to be regularly revisited in a broader light. At the LiveMEEG 2020 conference, a reflection on GSP was fostered that included explicitly documented guidelines and technical advances, but also emphasized intangible GSP: a general awareness of personal, organizational, and societal realities and how they can influence MEEG research. This article provides an extensive report on most of the LiveMEEG contributions and new literature, with the additional aim to synthesize ongoing cultural changes in GSP. It first covers GSP with respect to cognitive biases and logical fallacies, pre-registration as a tool to avoid those and other early pitfalls, and a number of resources to enable collaborative and reproducible research as a general approach to minimize misconceptions. Second, it covers GSP with respect to data acquisition, analysis, reporting, and sharing, including new tools and frameworks to support collaborative work. Finally, GSP is considered in light of ethical implications of MEEG research and the resulting responsibility that scientists have to engage with societal challenges. Considering among other things the benefits of peer review and open access at all stages, the need to coordinate larger international projects, the complexity of MEEG subject matter, and today's prioritization of fairness, privacy, and the environment, we find that current GSP tends to favor collective and cooperative work, for both scientific and for societal reasons.

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.

Robust muscle force prediction using NMFSEMD denoising and FOS identification

In this paper, an aliasing noise restraint technique and a system identification-based surface electromyography (sEMG)-force prediction model are proposed to realize a type of robust sEMG and muscle force prediction. For signal denoising, a novel non-negative matrix factorization screening empirical mode decomposition (NMFSEMD) and a fast orthogonal search (FOS)-based muscle force prediction model are developed. First, the NMFSEMD model is used to screen the empirical mode decomposition (EMD) results into the noisy intrinsic mode functions (IMF). Then, the noise matrix is computed using IMF translation and superposition, and the matrix is used as the input of NMF to obtain the denoised IMF. Furthermore, the reconstruction outcome of the NMFSEMD method can be used to estimate the denoised sEMG. Finally, a new sEMG muscle force prediction model, which considers a kind of candidate function in derivative form, is constructed, and a data-training-based linear weighted model is obtained. Extensive experimental results validate the suggested method's correction: after the NMFSEMD denoising of raw sEMG signal, the signal-noise ratio (SNR) can be improved by about 15.0 dB, and the energy percentage (EP) can be greater than 90.0%. Comparing with the muscle force prediction models using the traditional pretreatment and LSSVM, and the NMFSEMD plus LSSVM-based method, the mean square error (MSE) of our approach can be reduced by at least 1.2%.

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.

Emotional Activity Is Negatively Associated With Cognitive Load in Multimedia Learning: A Case Study With EEG Signals

We aimed to investigate the relationship between emotional activity and cognitive load during multimedia learning from an emotion dynamics perspective using electroencephalography (EEG) signals. Using a between-subjects design, 42 university students were randomly assigned to two video lecture conditions (color-coded vs. grayscale). While the participants watched the assigned video, their EEG signals were recorded. After processing the EEG signals, we employed the correlation-based feature selector (CFS) method to identify emotion-related subject-independent features. We then put these features into the Isomap model to obtain a one-dimensional trajectory of emotional changes. Next, we used the zero-crossing rate (ZCR) as the quantitative characterization of emotional changes ZCR EC . Meanwhile, we extracted cognitive load-related features to analyze the degree of cognitive load (CLI). We employed a linear regression fitting method to study the relationship between ZCR EC and CLI. We conducted this study from two perspectives. One is the frequency domain method (wavelet feature), and the other is the non-linear dynamic method (entropy features). The results indicate that emotional activity is negatively associated with cognitive load. These findings have practical implications for designing video lectures for multimedia learning. Learning material should reduce learners' cognitive load to keep their emotional experience at optimal levels to enhance learning.

Investigating the effect of Mindfulness-Based Stress Reduction on stress level and brain activity of college students

Financial constraints usually hinder students, especially those in low-middle income countries (LMICs), from seeking mental health interventions. Hence, it is necessary to identify effective, affordable and sustainable counter-stress measures for college students in the LMICs context. This study examines the sustained effects of mindfulness practice on the psychological outcomes and brain activity of students, especially when they are exposed to stressful situations. Here, we combined psychological and electrophysiological methods (EEG) to investigate the sustained effects of an 8-week-long standardized Mindfulness-Based Stress Reduction (MBSR) intervention on the brain activity of college students. We found that the Test group showed a decrease in negative emotional states after the intervention, compared to the no statistically significant result of the Control group, as indicated by the Perceived Stress Scale (PSS) (33% reduction in the negative score) and Depression, Anxiety, Stress Scale (DASS-42) scores (nearly 40% reduction of three subscale scores). Spectral analysis of EEG data showed that this intervention is longitudinally associated with increased frontal and occipital lobe alpha band power. Additionally, the increase in alpha power is more prevalent when the Test group was being stress-induced by cognitive tasks, suggesting that practicing MBSR might enhance the practitioners’ tolerance of negative emotional states. In conclusion, MBSR intervention led to a sustained reduction of negative emotional states as measured by both psychological and electrophysiological metrics, which supports the adoption of MBSR as an effective and sustainable stress-countering approach for students in LMICs.

The oscillatory effects of rhythmic median nerve stimulation

Entrainment of brain oscillations can be achieved using rhythmic non-invasive brain stimulation, and stimulation of the motor cortex at a frequency associated with sensorimotor inhibition can impair motor responses. Despite the potential for therapeutic application, these techniques do not lend themselves to use outside of a clinical setting. Here, the aim was to investigate whether rhythmic median nerve stimulation (MNS) could be used to entrain oscillations related to sensorimotor inhibition. MEG data were recorded from 20 participants during 400 trials, where for each trial 10 pulses of MNS were delivered either rhythmically or arrhythmically at 12 or 20 Hz. Our results demonstrate a frequency specific increase in relative amplitude in the contralateral somatosensory cortex during rhythmic but not arrhythmic stimulation. This was coupled with an increase in inter-trial phase coherence at the same frequency, suggesting that the oscillations synchronised with the pulses of MNS. The results show that 12 and 20 Hz rhythmic peripheral nerve stimulation can produce entrainment. Rhythmic MNS resulted in synchronous firing of neuronal populations within the contralateral somatosensory cortex meaning these neurons were engaged in processing of the afferent input. Therefore, MNS could prove therapeutically useful in disorders associated with hyperexcitability within the sensorimotor cortices.

In sync with your child: The potential of parent-child electroencephalography in developmental research

Healthy interaction between parent and child is foundational for the child's socioemotional development. Recently, an innovative paradigm shift in electroencephalography (EEG) research has enabled the simultaneous measurement of neural activity in caregiver and child. This dual-EEG or hyperscanning approach, termed parent-child dual-EEG, combines the strength of both behavioral observations and EEG methods. In this review, we aim to inform on the potential of dual-EEG in parents and children (0-6 years) for developmental researchers. We first provide a general overview of the dual-EEG technique and continue by reviewing the first empirical work on the emerging field of parent-child dual-EEG, discussing the limited but fascinating findings on parent-child brain-to-behavior and brain-to-brain synchrony. We then continue by providing an overview of dual-EEG analysis techniques, including the technical challenges and solutions one may encounter. We finish by discussing the potential of parent-child dual-EEG for the future of developmental research. The analysis of multiple EEG data is technical and challenging, but when performed well, parent-child EEG may transform the way we understand how caregiver and child connect on a neurobiological level. Importantly, studying objective physiological measures of parent-child interactions could lead to the identification of novel brain-to-brain synchrony markers of interaction quality.

Representation Learning for EEG-Based Biometrics Using Hilbert–Huang Transform

A promising approach to overcome the various shortcomings of password systems is the use of biometric authentication, in particular the use of electroencephalogram (EEG) data. In this paper, we propose a subject-independent learning method for EEG-based biometrics using Hilbert spectrograms of the data. The proposed neural network architecture treats the spectrogram as a collection of one-dimensional series and applies one-dimensional dilated convolutions over them, and a multi-similarity loss was used as the loss function for subject-independent learning. The architecture was tested on the publicly available PhysioNet EEG Motor Movement/Imagery Dataset (PEEGMIMDB) with a 14.63% Equal Error Rate (EER) achieved. The proposed approach’s main advantages are subject independence and suitability for interpretation via created spectrograms and the integrated gradients method.

Cognitive Outcome Prediction in Infants With Neonatal Hypoxic-Ischemic Encephalopathy Based on Functional Connectivity and Complexity of the Electroencephalography Signal

Impaired neurodevelopmental outcome, in particular cognitive impairment, after neonatal hypoxic-ischemic encephalopathy is a major concern for parents, clinicians, and society. This study aims to investigate the potential benefits of using advanced quantitative electroencephalography analysis (qEEG) for early prediction of cognitive outcomes, assessed here at 2 years of age. EEG data were recorded within the first week after birth from a cohort of twenty infants with neonatal hypoxic-ischemic encephalopathy (HIE). A proposed regression framework was based on two different sets of features, namely graph-theoretical features derived from the weighted phase-lag index (WPLI) and entropies metrics represented by sample entropy (SampEn), permutation entropy (PEn), and spectral entropy (SpEn). Both sets of features were calculated within the noise-assisted multivariate empirical mode decomposition (NA-MEMD) domain. Correlation analysis showed a significant association in the delta band between the proposed features, graph attributes (radius, transitivity, global efficiency, and characteristic path length) and entropy features (Pen and SpEn) from the neonatal EEG data and the cognitive development at age two years. These features were used to train and test the tree ensemble (boosted and bagged) regression models. The highest prediction performance was reached to 14.27 root mean square error (RMSE), 12.07 mean absolute error (MAE), and 0.45 R-squared using the entropy features with a boosted tree regression model. Thus, the results demonstrate that the proposed qEEG features show the state of brain function at an early stage; hence, they could serve as predictive biomarkers of later cognitive impairment, which could facilitate identifying those who might benefit from early targeted intervention.

Impacts of Color Coding on Programming Learning in Multimedia Learning: Moving Toward a Multimodal Methodology

In the present study, we tested the effectiveness of color coding on the programming learning of students who were learning from video lectures. Effectiveness was measured using multimodal physiological measures, combining eye tracking and electroencephalography (EEG). Using a between-subjects design, 42 university students were randomly assigned to two video lecture conditions (color-coded vs. grayscale). The participants' eye tracking and EEG signals were recorded while watching the assigned video, and their learning performance was subsequently assessed. The results showed that the color-coded design was more beneficial than the grayscale design, as indicated by smaller pupil diameter, shorter fixation duration, higher EEG theta and alpha band power, lower EEG cognitive load, and better learning performance. The present findings have practical implications for designing slide-based programming learning video lectures; slides should highlight the format of the program code using color coding.

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.

Multitaper estimates of phase-amplitude coupling

Phase-amplitude coupling (PAC) is the association of the amplitude of a high-frequency oscillation with the phase of a low-frequency oscillation. In neuroscience, this relationship provides a mechanism by which neural activity might be coordinated between distant regions. The dangers and pitfalls of assessing PAC with commonly used statistical measures have been well-documented. The limitations of these measures include: (1) response to non-oscillatory, high-frequency, broad-band activity, (2) response to high-frequency components of the low-frequency oscillation, (3) adhoc selection of analysis frequency-intervals, and (4) reliance upon data shuffling to assess statistical significance.Objective.To address issues (1)-(4) by introducing a nonparametric multitaper estimator of PAC.Approach.In this work, a multitaper PAC estimator is proposed that addresses these issues. Specifically, issue (1) is addressed by replacing the analytic signal envelope estimator computed using the Hilbert transform with a multitaper estimator that down-weights non-sinusoidal activity using a classical, multitaper super-resolution technique. Issue (2) is addressed by replacing coherence between the low-frequency and high-frequency components in a standard PAC estimator with multitaper partial coherence, while issue (3) is addressed with a physical argument regarding meaningful neural oscillation. Finally, asymptotic statistical assessment of the multitaper estimator is introduced to address issue (4).Main results.Multitaper estimates of PAC are introduced. Their efficacy is demonstrated in simulation and on human intracranial recordings obtained from epileptic patients.Significance.This work facilitates a more informative statistical assessment of PAC, a phenomena exhibited by many neural systems, and provides a basis upon which further nonparametric multitaper-related methods can be developed.

Aperiodic sleep networks promote memory consolidation

Hierarchical synchronization of sleep oscillations establishes communication pathways to support memory reactivation, transfer, and consolidation. From an information-theoretical perspective, oscillations constitute highly structured network states that provide limited information-coding capacity. Recent findings indicate that sleep oscillations occur in transient bursts that are interleaved with aperiodic network states, which were previously considered to be random noise. We argue that aperiodic activity exhibits unique and variable spatiotemporal patterns, providing an ideal information-rich neurophysiological substrate for imprinting new mnemonic patterns onto existing circuits. We discuss novel avenues in conceptualizing and quantifying aperiodic network states during sleep to further understand their relevance and interplay with sleep oscillations in support of memory consolidation.