Performance of different synchronization measures in real data: a case study on electroencephalographic signals

Effect of Eight-Week Transcranial Direct-Current Stimulation Combined with Lat Pull-Down Resistance Training on Improving Pull-Up Performance for Male College Students

The aim of this study was to investigate the effects and potential mechanisms of 8-week transcranial direct-current stimulation (tDCS) combined with resistance training (RT) on pull-up performance in male college students. Twenty-five male college students were randomly assigned to either RT combined with anodal tDCS stimulation (RT + tDCS) or RT alone (RT). Participants of both groups engaged in lat pull-down training programs for 8 weeks, with the RT + tDCS group receiving 20 min tDCS before each RT session. Pre- and post-intervention tests included pull-up endurance (number of repetitions), flexed arm circumference (FAC), and lat pull-down maximal voluntary isometric contraction (MVIC) peak force. During the pull-up endurance test, surface electromyography (sEMG) was recorded for the bicep brachii (BB), tricep brachii (TB), brachioradialis (BR), anterior deltoid (AD), middle deltoid (MD), posterior deltoid (PD), pectoralis major (PM), and latissimus dorsi (LD) muscles. Both groups demonstrated significant improvements in pull-up endurance and lat pull-down MVIC peak force after training, but no significant difference between the two groups was observed in the post-training test. Additionally, muscle activation of BR, PD, and PM decreased significantly in both groups, while the RT + tDCS group also demonstrated a significant reduction in TB coactivation after training. These findings suggest that eight weeks of tDCS combined with lat pull-down training and lat pull-down training alone can both significantly improve pull-up performance, which may be attributed to enhanced muscle contraction capacity. Although no significant training gains were found between the two training groups, the RT + tDCS group showed a significant decrease in TB coactivation and the enhancement of elbow flexion muscle contraction efficiency after training.

Reliable detection of directional couplings using cross-vector measures

Detecting directional couplings from time series is crucial in understanding complex dynamical systems. Various approaches based on reconstructed state-spaces have been developed for this purpose, including a cross-distance vector measure, which we introduced in our recent work. Here, we devise two new cross-vector measures that utilize ranks and time series estimates instead of distances. We analyze various deterministic and stochastic dynamics to compare our cross-vector approach against some established state-space-based approaches. We demonstrate that all three cross-vector measures can identify the correct coupling direction for a broader range of couplings for all considered dynamics. Among the three cross-vector measures, the rank-based variant performs the best. Comparing this novel measure to an established rank-based measure confirms that it is more noise-robust and less affected by linear cross-correlation. To extend this comparison to real-world signals, we combine both measures with the method of surrogates and analyze a database of electroencephalographic (EEG) recordings from epilepsy patients. This database contains signals from brain areas where the patients' seizures were detected first and signals from brain areas that were not involved in the seizure onset. A better discrimination between these signal classes is obtained by the cross-rank vector measure. Additionally, this measure proves to be robust to non-stationarity, as its results remain nearly unchanged when the analysis is repeated for the subset of EEG signals that were identified as stationary in previous work. These findings suggest that the cross-vector approach can serve as a valuable tool for researchers analyzing complex time series and for clinical applications.

The part and the whole: how single nodes contribute to large-scale phase-locking in functional EEG networks

OBJECTIVE

The application of signal analysis techniques to electroencephalographic (EEG) recordings from epilepsy patients shows that epilepsy involves not only altered neuronal synchronization but also the reorganization of functional EEG networks. This study aims to assess the large-scale phase-locking of such functional networks and how individual network nodes contribute to this collective dynamics.

METHODS

We analyze the EEG recorded before, during and after seizures from sixteen patients with pharmacoresistant focal-onset epilepsy. The data is filtered to low (4-30 Hz) and high (80-150 Hz) frequencies. We define the multivariate phase-locking measure and the univariate phase-locking contribution measure. Surrogate signals are used to estimate baseline results expected under the null hypothesis that the EEG is a correlated linear stochastic process.

RESULTS

On average, nodes from inside and outside the seizure onset zone (SOZ) increase and decrease, respectively, the large-scale phase-locking. This difference becomes most evident in a joint analysis of low and high frequencies.

CONCLUSIONS

Nodes inside and outside the SOZ play opposite roles for the large-scale phase-locking in functional EEG network in epilepsy patients.

SIGNIFICANCE

The application of the phase-locking contribution measure to EEG recordings from epilepsy patients can potentially help in localizing the SOZ.

Characterization of cardiorespiratory coupling via a variability-based multi-method approach: Application to postural orthostatic tachycardia syndrome

There are several mechanisms responsible for the dynamical link between heart period (HP) and respiration (R), usually referred to as cardiorespiratory coupling (CRC). Historically, diverse signal processing techniques have been employed to study CRC from the spontaneous fluctuations of HP and respiration (R). The proposed tools differ in terms of rationale and implementation, capturing diverse aspects of CRC. In this review, we classify the existing methods and stress differences with the aim of proposing a variability-based multi-method approach to CRC evaluation. Ten methodologies for CRC estimation, namely, power spectral decomposition, traditional and causal squared coherence,\;information transfer, cross-conditional entropy, mixed prediction, Shannon entropy of the latency between heartbeat and inspiratory/expiratory onset, conditional entropy of the phase dynamics, synchrogram-based analysis, pulse-respiration quotient, and joint symbolic dynamics, are considered. The ability of these techniques was exemplified over recordings acquired from patients suffering from postural orthostatic tachycardia syndrome (POTS) and healthy controls. Analyses were performed at rest in the supine position (REST) and during head-up tilt (HUT). Although most of the methods indicated that at REST, the CRC was lower in POTS patients and decreased more evidently during HUT in POTS, peculiar differences stressed the complementary value of the approaches. The multiple perspectives provided by the variability-based multi-method approach to CRC evaluation help the characterization of a pathological state and/or the quantification of the effect of a postural challenge. The present work stresses the need for the application of multiple methods to derive a more complete evaluation of the CRC in humans.

Generalized synchronization in the presence of dynamical noise and its detection via recurrent neural networks

Given two unidirectionally coupled nonlinear systems, we speak of generalized synchronization when the responder "follows" the driver. Mathematically, this situation is implemented by a map from the driver state space to the responder state space termed the synchronization map. In nonlinear times series analysis, the framework of the present work, the existence of the synchronization map amounts to the invertibility of the so-called cross map, which is a continuous map that exists in the reconstructed state spaces for typical time-delay embeddings. The cross map plays a central role in some techniques to detect functional dependencies between time series. In this paper, we study the changes in the "noiseless scenario" just described when noise is present in the driver, a more realistic situation that we call the "noisy scenario." Noise will be modeled using a family of driving dynamics indexed by a finite number of parameters, which is sufficiently general for practical purposes. In this approach, it turns out that the cross and synchronization maps can be extended to the noisy scenario as families of maps that depend on the noise parameters, and only for "generic" driver states in the case of the cross map. To reveal generalized synchronization in both the noiseless and noisy scenarios, we check the existence of synchronization maps of higher periods (introduced in this paper) using recurrent neural networks and predictability. The results obtained with synthetic and real-world data demonstrate the capability of our method.

Fractal Conditional Correlation Dimension Infers Complex Causal Networks

Determining causal inference has become popular in physical and engineering applications. While the problem has immense challenges, it provides a way to model the complex networks by observing the time series. In this paper, we present the optimal conditional correlation dimensional geometric information flow principle (oGeoC) that can reveal direct and indirect causal relations in a network through geometric interpretations. We introduce two algorithms that utilize the oGeoC principle to discover the direct links and then remove indirect links. The algorithms are evaluated using coupled logistic networks. The results indicate that when the number of observations is sufficient, the proposed algorithms are highly accurate in identifying direct causal links and have a low false positive rate.

Effect of fatigue on intermuscular EMG-EMG coupling during bench press exercise at 60% 1RM workload in males

Objective

To explore the neuromuscular control mechanism and quantifying the fatigue response during bench press exercise is important references to prescribe an appropriate exercise program. However, current literature struggles to provide a concrete conclusion on the changes of intermuscular EMG-EMG coupling between synergistic and antagonist muscles during the exercise. Thus, the current study was designed to reveal fatigue-related changes of intermuscular EMG-EMG coupling during bench press exercise.

Methods

Thirty-one healthy male participants performed a bench press exercise on the Smith machine at 60% One Repetition Maximum (1RM) workload to exhaustion, while surface electromyographic signals (sEMG) were collected from triceps brachii (TB), biceps brachii (BB), anterior deltoid (AD), posterior deltoid (PD), and pectoralis major (PM). Surface EMG signals were divided into the first half and second half of the bench press exercise. Phase synchronization index (PSI) was calculated between sEMG of synergistic muscle pairs AD-TB, AD-PM and antagonist muscle pairs BB-TB, AD-PD.

Results

EMG power of TB, AD, PD, PM muscles in alpha (8-12 Hz) frequency band and EMG power of each muscle in beta (15-35 Hz), and gamma (35-60 Hz) frequency bands were all increased during the second half of contraction compared with the first half of contraction. PSI of gamma frequency band was significantly decreased in BB-TB muscle pair while EMG-EMG coupling of AD-TB in gamma frequency band was significantly increased during the second half of contraction compared to the first half of contraction.

Conclusions

The results indicated a decrease of interconnection between synchronized cortical neurons and the motoneuron pool of BB and TB, and an increase of interconnection between AD-TB muscles during fatiguing bench press exercise at 60% 1RM workload. The changes of intermuscular coupling may be related to the supraspinal modulations to compensate for the decrease of muscle force as well as a result of unbalanced changes of agonist and antagonist muscle contractility.

Assessment of Fractal Synchronization during an Epileptic Seizure

In this paper, we define fractal synchronization (FS) based on the idea of stochastic synchronization and propose a mathematical apparatus for estimating FS. One major advantage of our proposed approach is that fractal synchronization makes it possible to estimate the aggregate strength of the connection on multiple time scales between two projections of the attractor, which are time series with a fractal structure. We believe that one of the promising uses of FS is the assessment of the interdependence of encephalograms. To demonstrate this approach in evaluating the cross-dependence between channels in a network of electroencephalograms, we evaluated the FS of encephalograms during an epileptic seizure. Fractal synchronization demonstrates the presence of desynchronization during an epileptic seizure.

A rigorous and versatile statistical test for correlations between stationary time series

In disciplines from biology to climate science, a routine task is to compute a correlation between a pair of time series and determine whether the correlation is statistically significant (i.e., unlikely under the null hypothesis that the time series are independent). This problem is challenging because time series typically exhibit autocorrelation and thus cannot be properly analyzed with the standard iid-oriented statistical tests. Although there are well-known parametric tests for time series, these are designed for linear correlation statistics and thus not suitable for the increasingly popular nonlinear correlation statistics. There are also nonparametric tests that can be used with any correlation statistic, but for these, the conditions that guarantee correct false positive rates are either restrictive or unclear. Here, we describe the truncated time-shift (TTS) test, a nonparametric procedure to test for dependence between 2 time series. We prove that this test correctly controls the false positive rate as long as one of the time series is stationary, a minimally restrictive requirement among current tests. The TTS test is versatile because it can be used with any correlation statistic. Using synthetic data, we demonstrate that this test performs correctly even while other tests suffer high false positive rates. In simulation examples, simple guidelines for parameter choices allow high statistical power to be achieved with sufficient data. We apply the test to datasets from climatology, animal behavior, and microbiome science, verifying previously discovered dependence relationships and detecting additional relationships.

An EEG-based marker of functional connectivity: detection of major depressive disorder

Major depressive disorder (MDD) is a prevalent psychiatric disorder globally. There are many assays for MDD, but rapid and reliable detection remains a pressing challenge. In this study, we present a fusion feature called P-MSWC, as a novel marker to construct brain functional connectivity matrices and utilize the convolutional neural network (CNN) to identify MDD based on electroencephalogram (EEG) signal. Firstly, we combine synchrosqueezed wavelet transform and coherence theory to get synchrosqueezed wavelet coherence. Then, we obtain the fusion feature by incorporating synchrosqueezed wavelet coherence value and phase-locking value, which outperforms conventional functional connectivity markers by comprehensively capturing the original EEG signal's information and demonstrating notable noise-resistance capabilities. Finally, we propose a lightweight CNN model that effectively utilizes the high-dimensional connectivity matrix of the brain, constructed using our novel marker, to enable more accurate and efficient detection of MDD. The proposed method achieves 99.92% accuracy on a single dataset and 97.86% accuracy on a combined dataset. Moreover, comparison experiments have shown that the performance of the proposed method is superior to traditional machine learning methods. Furthermore, visualization experiments reveal differences in the distribution of brain connectivity between MDD patients and healthy subjects, including decreased connectivity in the T7, O1, F8, and C3 channels of the gamma band. The results of the experiments indicate that the fusion feature can be utilized as a new marker for constructing functional brain connectivity, and the combination of deep learning and functional connectivity matrices can provide more help for the detection of MDD.

Response inhibition in premotor cortex corresponds to a complex reshuffle of the mesoscopic information network

Recent studies have explored functional and effective neural networks in animal models; however, the dynamics of information propagation among functional modules under cognitive control remain largely unknown. Here, we addressed the issue using transfer entropy and graph theory methods on mesoscopic neural activities recorded in the dorsal premotor cortex of rhesus monkeys. We focused our study on the decision time of a Stop-signal task, looking for patterns in the network configuration that could influence motor plan maturation when the Stop signal is provided. When comparing trials with successful inhibition to those with generated movement, the nodes of the network resulted organized into four clusters, hierarchically arranged, and distinctly involved in information transfer. Interestingly, the hierarchies and the strength of information transmission between clusters varied throughout the task, distinguishing between generated movements and canceled ones and corresponding to measurable levels of network complexity. Our results suggest a putative mechanism for motor inhibition in premotor cortex: a topological reshuffle of the information exchanged among ensembles of neurons.

Statistical Interdependence between Daily Precipitation and Extreme Daily Temperature in Regions of Mexico and Colombia

We study the statistical interdependence between daily precipitation and daily extreme temperature for regions of Mexico (14 climatic stations, period 1960-2020) and Colombia (7 climatic stations, period 1973-2020) using linear (cross-correlation and coherence) and nonlinear (global phase synchronization index, mutual information, and cross-sample entropy) synchronization metrics. The information shared between these variables is relevant and exhibits changes when comparing regions with different climatic conditions. We show that precipitation and temperature records from La Mojana are characterized by high persistence, while data from Mexico City exhibit lower persistence (less memory). We find that the information exchange and the level of coupling between the precipitation and temperature are higher for the case of the La Mojana region (Colombia) compared to Mexico City (Mexico), revealing that regions where seasonal changes are almost null and with low temperature gradients (less local variability) tend to display higher synchrony compared to regions where seasonal changes are very pronounced. The interdependence characterization between precipitation and temperature represents a robust option to characterize and analyze the collective dynamics of the system, applicable in climate change studies, as well as in changes not easily identifiable in future scenarios.

The causality measure of partial mutual information from mixed embedding (PMIME) revisited

The measure of partial mutual information from mixed embedding (PMIME) is an information theory-based measure to accurately identify the direct and directional coupling, termed Granger causality or simply causality, between the observed variables or subsystems of a high-dimensional dynamical and complex system, without any a priori assumptions about the nature of the coupling relationship. In its core, it is a forward selection procedure that aims to iteratively identify the lag-dependence structure of a given observed variable (response) to all the other observed variables (candidate drivers). This model-free approach is capable of detecting nonlinear interactions, abundantly present in real-world complex systems, and it was shown to perform well on multivariate time series of moderately high dimension. However, the PMIME presents some inefficiencies in its performance mainly when applied on strongly stochastic (linear or nonlinear) systems as it may falsely detect non-existent relationships. Moreover, and by construction, the measure cannot extract purely synergetic relationships present in a system. In the current work, the issue of false detections is addressed by introducing an improved resampling significance test and a procedure of rechecking the identified drivers (backward revision). Regarding the inability to detect synergetic relationships, the PMIME is further enhanced by checking pairs as candidate drivers for the response variable after having considered all drivers individually. The effects of these modifications are investigated in a systematic simulation study on properly designed systems involving strong stochasticity, regressor terms with synergetic effects, and a system dimension ranging from 3 to 30. The overall results of the simulations indicate that these modifications indeed improve the performance of PMIME and alleviate to a significant degree the issues of the original algorithm. Guidelines for balancing between accuracy and computational efficiency are also given, particularly relevant for real-world applications. Finally, the measure performance is investigated in the study of futures of various government bonds and stock market indices in the period around COVID-19 pandemic.

Spike-and-wave discharges of absence seizures in a sleep waves-constrained corticothalamic model

AIMS

Recurrent network activity in corticothalamic circuits generates physiological and pathological EEG waves. Many computer models have simulated spike-and-wave discharges (SWDs), the EEG hallmark of absence seizures (ASs). However, these models either provided detailed simulated activity only in a selected territory (i.e., cortical or thalamic) or did not test whether their corticothalamic networks could reproduce the physiological activities that are generated by these circuits.

METHODS

Using a biophysical large-scale corticothalamic model that reproduces the full extent of EEG sleep waves, including sleep spindles, delta, and slow (<1 Hz) waves, here we investigated how single abnormalities in voltage- or transmitter-gated channels in the neocortex or thalamus led to SWDs.

RESULTS

We found that a selective increase in the tonic γ-aminobutyric acid type A receptor (GABA-A) inhibition of first-order thalamocortical (TC) neurons or a selective decrease in cortical phasic GABA-A inhibition is sufficient to generate ~4 Hz SWDs (as in humans) that invariably start in neocortical territories. Decreasing the leak conductance of higher-order TC neurons leads to ~7 Hz SWDs (as in rodent models) while maintaining sleep spindles at 7-14 Hz.

CONCLUSION

By challenging key features of current mechanistic views, this simulated ictal corticothalamic activity provides novel understanding of ASs and makes key testable predictions.

A machine learning approach for real-time cortical state estimation

Objective.Cortical function is under constant modulation by internally-driven, latent variables that regulate excitability, collectively known as 'cortical state'. Despite a vast literature in this area, the estimation of cortical state remains relatively ad hoc, and not amenable to real-time implementation. Here, we implement robust, data-driven, and fast algorithms that address several technical challenges for online cortical state estimation.Approach. We use unsupervised Gaussian mixture models to identify discrete, emergent clusters in spontaneous local field potential signals in cortex. We then extend our approach to a temporally-informed hidden semi-Markov model (HSMM) with Gaussian observations to better model and infer cortical state transitions. Finally, we implement our HSMM cortical state inference algorithms in a real-time system, evaluating their performance in emulation experiments.Main results. Unsupervised clustering approaches reveal emergent state-like structure in spontaneous electrophysiological data that recapitulate arousal-related cortical states as indexed by behavioral indicators. HSMMs enable cortical state inferences in a real-time context by modeling the temporal dynamics of cortical state switching. Using HSMMs provides robustness to state estimates arising from noisy, sequential electrophysiological data.Significance. To our knowledge, this work represents the first implementation of a real-time software tool for continuously decoding cortical states with high temporal resolution (40 ms). The software tools that we provide can facilitate our understanding of how cortical states dynamically modulate cortical function on a moment-by-moment basis and provide a basis for state-aware brain machine interfaces across health and disease.

Temporal segmentation of EEG based on functional connectivity network structure

In the study of brain functional connectivity networks, it is assumed that a network is built from a data window in which activity is stationary. However, brain activity is non-stationary over sufficiently large time periods. Addressing the analysis electroencephalograph (EEG) data, we propose a data segmentation method based on functional connectivity network structure. The goal of segmentation is to ensure that within a window of analysis, there is similar network structure. We designed an intuitive and flexible graph distance measure to quantify the difference in network structure between two analysis windows. This measure is modular: a variety of node importance indices can be plugged into it. We use a reference window versus sliding window comparison approach to detect changes, as indicated by outliers in the distribution of graph distance values. Performance of our segmentation method was tested in simulated EEG data and real EEG data from a drone piloting experiment (using correlation or phase-locking value as the functional connectivity strength metric). We compared our method under various node importance measures and against matrix-based dissimilarity metrics that use singular value decomposition on the connectivity matrix. The results show the graph distance approach worked better than matrix-based approaches; graph distance based on partial node centrality was most sensitive to network structural changes, especially when connectivity matrix values change little. The proposed method provides EEG data segmentation tailored for detecting changes in terms of functional connectivity networks. Our study provides a new perspective on EEG segmentation, one that is based on functional connectivity network structure differences.

Generalized synchronization mediated by a flat coupling between structurally nonequivalent chaotic systems

Synchronization of chaotic systems is usually investigated for structurally equivalent systems typically coupled through linear diffusive functions. Here, we focus on a particular type of coupling borrowed from a nonlinear control theory and based on the optimal placement of a sensor-a device measuring the chosen variable-and an actuator-a device applying the actuating (control) signal to a variable's derivative-in the response system, leading to the so-called flat control law. We aim to investigate the dynamics produced by a response system that is flat coupled to a drive system and to determine the degree of generalized synchronization between them using statistical and topological arguments. The general use of a flat control law for getting generalized synchronization is discussed.

Weighted ordinal connection based functional network classification for schizophrenia disease detection using EEG signal

A brain connectivity network (BCN) is an advanced approach to examining brain functionality in various conditions. However, the predictability of the BCN is affected by the connectivity measure used for the network construction. Various connectivity measures available in the literature differ according to the domain of their working data. The application of random connectivity measures might result in an inefficient BCN that ultimately hampers its predictability. Therefore, selecting an appropriate functional connectivity metric is crucial in clinical as well as cognitive neuroscience. In parallel to this, an effective network identifier plays a vital role in distinguishing different brain states. Hence, the objective of this paper is two-fold, which includes identifying suitable connectivity measures and proposing an efficient network identifier. For this, the weighted BCN (WBCN) is constructed using multiple connectivity measures like correlation coefficient (r), coherence (COH), phase-locking value (PLV), and mutual information (MI) from electroencephalogram (EEG) signals. The most recent technique for feature extraction, i.e., weighted ordinal connections, has been applied to EEG-based BCN. EEG signals data has been taken from the schizophrenia disease database. Further, several classification algorithms such as k-nearest neighbours (KNN), support vector machine (SVM) with linear, radial basis function and polynomial kernels, random forest (RF), and 1D convolutional neural network (CNN1D) are used to classify the brain states based on extracted features. In classification, 90% accuracy is achieved by the CNN1D classifier with WBCN based on the coherence connectivity measure. The study also provides a structural analysis of the BCN.

Comparison of Bootstrap Methods for Estimating Causality in Linear Dynamic Systems: A Review

In this study, we present a thorough comparison of the performance of four different bootstrap methods for assessing the significance of causal analysis in time series data. For this purpose, multivariate simulated data are generated by a linear feedback system. The methods investigated are uncorrelated Phase Randomization Bootstrap (uPRB), which generates surrogate data with no cross-correlation between variables by randomizing the phase in the frequency domain; Time Shift Bootstrap (TSB), which generates surrogate data by randomizing the phase in the time domain; Stationary Bootstrap (SB), which calculates standard errors and constructs confidence regions for weakly dependent stationary observations; and AR-Sieve Bootstrap (ARSB), a resampling method based on AutoRegressive (AR) models that approximates the underlying data-generating process. The uPRB method accurately identifies variable interactions but fails to detect self-feedback in some variables. The TSB method, despite performing worse than uPRB, is unable to detect feedback between certain variables. The SB method gives consistent causality results, although its ability to detect self-feedback decreases, as the mean block width increases. The ARSB method shows superior performance, accurately detecting both self-feedback and causality across all variables. Regarding the analysis of the Impulse Response Function (IRF), only the ARSB method succeeds in detecting both self-feedback and causality in all variables, aligning well with the connectivity diagram. Other methods, however, show considerable variations in detection performance, with some detecting false positives and others only detecting self-feedback.

Usefulness and limitations of convergent cross sorting and continuity scaling methods for their application in simulated and real-world time series

Causality detection methods are valuable tools for detecting causal links in complex systems. The efficiency of continuity scaling (CS) and the convergent cross sorting (CSS) methods to detect causality was analysed. Usefulness and limitations of both methods in their application to simulated and real-world time series was explored under different scenarios. We find that CS is more robust and efficient than the CSS method for all simulated systems, even when increasing noise levels were considered. Both methods were not able to infer causality when time series with a marked difference in their main frequencies were analysed. Minimum time-series length required for the detection of a causal link depends on intrinsic system dynamics and on the method selected to detect it. Using simulated time series, only the CS method was capable to detect bidirectional causality. Causality detection, using the CS method, should at least include: (i) causality strength convergence analysis, (ii) statistical tests of significance, (iii) time-series standardization, and (iv) causality strength ratios as a strength indicator of relative causality between systems. Causality cannot be detected by either method in simulated time series that exhibit generalized synchronization.

EEG phase synchronization during absence seizures

Absence seizures-generalized rhythmic spike-and-wave discharges (SWDs) are the defining property of childhood (CAE) and juvenile (JAE) absence epilepsies. Such seizures are the most compelling examples of pathological neuronal hypersynchrony. All the absence detection algorithms proposed so far have been derived from the properties of individual SWDs. In this work, we investigate EEG phase synchronization in patients with CAE/JAE and healthy subjects to explore the possibility of using the wavelet phase synchronization index to detect seizures and quantify their disorganization (fragmentation). The overlap of the ictal and interictal probability density functions was high enough to preclude effective seizure detection based solely on changes in EEG synchronization. We used a machine learning classifier with the phase synchronization index (calculated for 1 s data segments with 0.5 s overlap) and the normalized amplitude as features to detect generalized SWDs. Using 19 channels (10-20 setup), we identified 99.2% of absences. However, the overlap of the segments classified as ictal with seizures was only 83%. The analysis showed that seizures were disorganized in approximately half of the 65 subjects. On average, generalized SWDs lasted about 80% of the duration of abnormal EEG activity. The disruption of the ictal rhythm can manifest itself as the disappearance of epileptic spikes (with high-amplitude delta waves persisting), transient cessation of epileptic discharges, or loss of global synchronization. The detector can analyze a real-time data stream. Its performance is good for a six-channel setup (Fp1, Fp2, F7, F8, O1, O2), which can be implemented as an unobtrusive EEG headband. False detections are rare for controls and young adults (0.03% and 0.02%, respectively). In patients, they are more frequent (0.5%), but in approximately 82% cases, classification errors are caused by short epileptiform discharges. Most importantly, the proposed detector can be applied to parts of EEG with abnormal EEG activity to quantitatively determine seizure fragmentation. This property is important because a previous study reported that the probability of disorganized discharges is eight times higher in JAE than in CAE. Future research must establish whether seizure properties (frequency, length, fragmentation, etc.) and clinical characteristics can help distinguish CAE and JAE.

Compensatory mechanisms of reduced interhemispheric EEG connectivity during sleep in patients with apnea

We performed a mathematical analysis of functional connectivity in electroencephalography (EEG) of patients with obstructive sleep apnea (OSA) (N = 10; age: 52.8 ± 13 years; median age: 49 years; male/female ratio: 7/3), compared with a group of apparently healthy participants (N = 15; age: 51.5 ± 29.5 years; median age: 42 years; male/female ratio: 8/7), based on the calculation of wavelet bicoherence from nighttime polysomnograms. Having observed the previously known phenomenon of interhemispheric synchronization deterioration, we demonstrated a compensatory increase in intrahemispheric connectivity, as well as a slight increase in the connectivity of the central and occipital areas for high-frequency EEG activity. Significant changes in functional connectivity were extremely stable in groups of apparently healthy participants and OSA patients, maintaining the overall pattern when comparing different recording nights and various sleep stages. The maximum variability of the connectivity was observed at fast oscillatory processes during REM sleep. The possibility of observing some changes in functional connectivity of brain activity in OSA patients in a state of passive wakefulness opens up prospects for further research. Developing the methods of hypnogram evaluation that are independent of functional connectivity may be useful for implementing a medical decision support system.

A Survey on Brain Effective Connectivity Network Learning

Human brain effective connectivity characterizes the causal effects of neural activities among different brain regions. Studies of brain effective connectivity networks (ECNs) for different populations contribute significantly to the understanding of the pathological mechanism associated with neuropsychiatric diseases and facilitate finding new brain network imaging markers for the early diagnosis and evaluation for the treatment of cerebral diseases. A deeper understanding of brain ECNs also greatly promotes brain-inspired artificial intelligence (AI) research in the context of brain-like neural networks and machine learning. Thus, how to picture and grasp deeper features of brain ECNs from functional magnetic resonance imaging (fMRI) data is currently an important and active research area of the human brain connectome. In this survey, we first show some typical applications and analyze existing challenging problems in learning brain ECNs from fMRI data. Second, we give a taxonomy of ECN learning methods from the perspective of computational science and describe some representative methods in each category. Third, we summarize commonly used evaluation metrics and conduct a performance comparison of several typical algorithms both on simulated and real datasets. Finally, we present the prospects and references for researchers engaged in learning ECNs.

Metastability as a candidate neuromechanistic biomarker of schizophrenia pathology

The disconnection hypothesis of schizophrenia proposes that symptoms of the disorder arise as a result of aberrant functional integration between segregated areas of the brain. The concept of metastability characterizes the coexistence of competing tendencies for functional integration and functional segregation in the brain, and is therefore well suited for the study of schizophrenia. In this study, we investigate metastability as a candidate neuromechanistic biomarker of schizophrenia pathology, including a demonstration of reliability and face validity. Group-level discrimination, individual-level classification, pathophysiological relevance, and explanatory power were assessed using two independent case-control studies of schizophrenia, the Human Connectome Project Early Psychosis (HCPEP) study (controls n = 53, non-affective psychosis n = 82) and the Cobre study (controls n = 71, cases n = 59). In this work we extend Leading Eigenvector Dynamic Analysis (LEiDA) to capture specific features of dynamic functional connectivity and then implement a novel approach to estimate metastability. We used non-parametric testing to evaluate group-level differences and a naïve Bayes classifier to discriminate cases from controls. Our results show that our new approach is capable of discriminating cases from controls with elevated effect sizes relative to published literature, reflected in an up to 76% area under the curve (AUC) in out-of-sample classification analyses. Additionally, our new metric showed explanatory power of between 81-92% for measures of integration and segregation. Furthermore, our analyses demonstrated that patients with early psychosis exhibit intermittent disconnectivity of subcortical regions with frontal cortex and cerebellar regions, introducing new insights about the mechanistic bases of these conditions. Overall, these findings demonstrate reliability and face validity of metastability as a candidate neuromechanistic biomarker of schizophrenia pathology.

Comparison of discretization strategies for the model-free information-theoretic assessment of short-term physiological interactions

This work presents a comparison between different approaches for the model-free estimation of information-theoretic measures of the dynamic coupling between short realizations of random processes. The measures considered are the mutual information rate (MIR) between two random processes X and Y and the terms of its decomposition evidencing either the individual entropy rates of X and Y and their joint entropy rate, or the transfer entropies from X to Y and from Y to X and the instantaneous information shared by X and Y. All measures are estimated through discretization of the random variables forming the processes, performed either via uniform quantization (binning approach) or rank ordering (permutation approach). The binning and permutation approaches are compared on simulations of two coupled non-identical Hènon systems and on three datasets, including short realizations of cardiorespiratory (CR, heart period and respiration flow), cardiovascular (CV, heart period and systolic arterial pressure), and cerebrovascular (CB, mean arterial pressure and cerebral blood flow velocity) measured in different physiological conditions, i.e., spontaneous vs paced breathing or supine vs upright positions. Our results show that, with careful selection of the estimation parameters (i.e., the embedding dimension and the number of quantization levels for the binning approach), meaningful patterns of the MIR and of its components can be achieved in the analyzed systems. On physiological time series, we found that paced breathing at slow breathing rates induces less complex and more coupled CR dynamics, while postural stress leads to unbalancing of CV interactions with prevalent baroreflex coupling and to less complex pressure dynamics with preserved CB interactions. These results are better highlighted by the permutation approach, thanks to its more parsimonious representation of the discretized dynamic patterns, which allows one to explore interactions with longer memory while limiting the curse of dimensionality.

Detecting Nonlinear Interactions in Complex Systems: Application in Financial Markets

Emerging or diminishing nonlinear interactions in the evolution of a complex system may signal a possible structural change in its underlying mechanism. This type of structural break may exist in many applications, such as in climate and finance, and standard methods for change-point detection may not be sensitive to it. In this article, we present a novel scheme for detecting structural breaks through the occurrence or vanishing of nonlinear causal relationships in a complex system. A significance resampling test was developed for the null hypothesis (H₀) of no nonlinear causal relationships using (a) an appropriate Gaussian instantaneous transform and vector autoregressive (VAR) process to generate the resampled multivariate time series consistent with H₀; (b) the modelfree Granger causality measure of partial mutual information from mixed embedding (PMIME) to estimate all causal relationships; and (c) a characteristic of the network formed by PMIME as test statistic. The significance test was applied to sliding windows on the observed multivariate time series, and the change from rejection to no-rejection of H₀, or the opposite, signaled a non-trivial change of the underlying dynamics of the observed complex system. Different network indices that capture different characteristics of the PMIME networks were used as test statistics. The test was evaluated on multiple synthetic complex and chaotic systems, as well as on linear and nonlinear stochastic systems, demonstrating that the proposed methodology is capable of detecting nonlinear causality. Furthermore, the scheme was applied to different records of financial indices regarding the global financial crisis of 2008, the two commodity crises of 2014 and 2020, the Brexit referendum of 2016, and the outbreak of COVID-19, accurately identifying the structural breaks at the identified times.

Exploration of sleep function connection and classification strategies based on sub-period sleep stages

Background

As a medium for developing brain-computer interface systems, EEG signals are complex and difficult to identify due to their complexity, weakness, and differences between subjects. At present, most of the current research on sleep EEG signals are single-channel and dual-channel, ignoring the research on the relationship between different brain regions. Brain functional connectivity is considered to be closely related to brain activity and can be used to study the interaction relationship between brain areas.

Methods

Phase-locked value (PLV) is used to construct a functional connection network. The connection network is used to analyze the connection mechanism and brain interaction in different sleep stages. Firstly, the entire EEG signal is divided into multiple sub-periods. Secondly, Phase-locked value is used for feature extraction on the sub-periods. Thirdly, the PLV of multiple sub-periods is used for feature fusion. Fourthly, the classification performance optimization strategy is used to discuss the impact of different frequency bands on sleep stage classification performance and to find the optimal frequency band. Finally, the brain function network is constructed by using the average value of the fusion features to analyze the interaction of brain regions in different frequency bands during sleep stages.

Results

The experimental results have shown that when the number of sub-periods is 30, the α (8-13 Hz) frequency band has the best classification effect, The classification result after 10-fold cross-validation reaches 92.59%.

Conclusion

The proposed algorithm has good sleep staging performance, which can effectively promote the development and application of an EEG sleep staging system.

Mechanisms of hemoglobin cycling in anemia patients treated with erythropoiesis-stimulating agents

Patients with renal anemia are frequently treated with erythropoiesis-stimulating agents (ESAs), which are dynamically dosed in order to stabilize blood hemoglobin levels within a specified target range. During typical ESA treatments, a fraction of patients experience hemoglobin 'cycling' periods during which hemoglobin levels periodically over- and undershoot the target range. Here we report a specific mechanism of hemoglobin cycling, whereby cycles emerge from the patient's delayed physiological response to ESAs and concurrent ESA dose adjustments. We introduce a minimal theoretical model that can explain dynamic hallmarks of observed hemoglobin cycling events in clinical time series and elucidates how physiological factors (such as red blood cell lifespan and ESA responsiveness) and treatment-related factors (such as dosing schemes) affect cycling. These results show that in general, hemoglobin cycling cannot be attributed to patient physiology or ESA treatment alone but emerges through an interplay of both, with consequences for the design of ESA treatment strategies.

Adaptation of Partial Mutual Information from Mixed Embedding to Discrete-Valued Time Series

A causality analysis aims at estimating the interactions of the observed variables and subsequently the connectivity structure of the observed dynamical system or stochastic process. The partial mutual information from mixed embedding (PMIME) is found appropriate for the causality analysis of continuous-valued time series, even of high dimension, as it applies a dimension reduction by selecting the most relevant lag variables of all the observed variables to the response, using conditional mutual information (CMI). The presence of lag components of the driving variable in this vector implies a direct causal (driving-response) effect. In this study, the PMIME is appropriately adapted to discrete-valued multivariate time series, called the discrete PMIME (DPMIME). An appropriate estimation of the discrete probability distributions and CMI for discrete variables is implemented in the DPMIME. Further, the asymptotic distribution of the estimated CMI is derived, allowing for a parametric significance test for the CMI in the DPMIME, whereas for the PMIME, there is no parametric test for the CMI and the test is performed using resampling. Monte Carlo simulations are performed using different generating systems of discrete-valued time series. The simulation suggests that the parametric significance test for the CMI in the progressive algorithm of the DPMIME is compared favorably to the corresponding resampling significance test, and the accuracy of the DPMIME in the estimation of direct causality converges with the time-series length to the accuracy of the PMIME. Further, the DPMIME is used to investigate whether the global financial crisis has an effect on the causality network of the financial world market.

Metastability, fractal scaling, and synergistic information processing: What phase relationships reveal about intrinsic brain activity

Dynamic functional connectivity (dFC) in resting-state fMRI holds promise to deliver candidate biomarkers for clinical applications. However, the reliability and interpretability of dFC metrics remain contested. Despite a myriad of methodologies and resulting measures, few studies have combined metrics derived from different conceptualizations of brain functioning within the same analysis - perhaps missing an opportunity for improved interpretability. Using a complexity-science approach, we assessed the reliability and interrelationships of a battery of phase-based dFC metrics including tools originating from dynamical systems, stochastic processes, and information dynamics approaches. Our analysis revealed novel relationships between these metrics, which allowed us to build a predictive model for integrated information using metrics from dynamical systems and information theory. Furthermore, global metastability - a metric reflecting simultaneous tendencies for coupling and decoupling - was found to be the most representative and stable metric in brain parcellations that included cerebellar regions. Additionally, spatiotemporal patterns of phase-locking were found to change in a slow, non-random, continuous manner over time. Taken together, our findings show that the majority of characteristics of resting-state fMRI dynamics reflect an interrelated dynamical and informational complexity profile, which is unique to each acquisition. This finding challenges the interpretation of results from cross-sectional designs for brain neuromarker discovery, suggesting that individual life-trajectories may be more informative than sample means.

Repetitive nociceptive stimulation increases spontaneous neural activation similar to nociception-induced activity in mouse insular cortex

Recent noninvasive neuroimaging technology has revealed that spatiotemporal patterns of cortical spontaneous activity observed in chronic pain patients are different from those in healthy subjects, suggesting that the spontaneous cortical activity plays a key role in the induction and/or maintenance of chronic pain. However, the mechanisms of the spontaneously emerging activities supposed to be induced by nociceptive inputs remain to be established. In the present study, we investigated spontaneous cortical activities in sessions before and after electrical stimulation of the periodontal ligament (PDL) by applying wide-field and two-photon calcium imaging to anesthetized GCaMP6s transgenic mice. First, we identified the sequential cortical activation patterns from the primary somatosensory and secondary somatosensory cortices to the insular cortex (IC) by PDL stimulation. We, then found that spontaneous IC activities that exhibited a similar spatiotemporal cortical pattern to evoked activities by PDL stimulation increased in the session after repetitive PDL stimulation. At the single-cell level, repetitive PDL stimulation augmented the synchronous neuronal activity. These results suggest that cortical plasticity induced by the repetitive stimulation leads to the frequent PDL stimulation-evoked-like spontaneous IC activation. This nociception-induced spontaneous activity in IC may be a part of mechanisms that induces chronic pain.

A comparison of muscle activation and concomitant intermuscular coupling of antagonist muscles among bench presses with different instability degrees in untrained men

The aim of this study was to analyze and compare the muscle activation and concomitant intermuscular coupling of antagonist muscles among bench presses with different instability degrees. Twenty-nine untrained male college students performed bench press exercises at an intensity of 60% 1 RM on three conditions: small unstable bench press with Smith machine (SBP), medium unstable bench press of free weight (FWBP), and large unstable bench press with increased instability by suspending the load with elastic bands (IIBP). One-way repeated measures analysis of variance was used to compare integrated EMG activity values of the biceps brachii (BB), posterior deltoid (PD), long head of the triceps brachii (TB), anterior deltoid (AD), upper portion of the pectoralis major (PM) muscles, and phase synchronization index (PSI) of BB-TB and PD-AD antagonist muscle pairs. A higher integrated EMG of BB muscle was found during bench press with a more unstable condition. IIBP showed a higher integrated EMG of prime movers (TB, AD, and PM) and stabilizing of BB than SBP and FWBP. PSI between muscle pairs of BB-TB in the gamma frequency band was higher in SBP than the other bench presses with unstable conditions, which may be related to the optimal “internal model” for antagonist muscles during bench press exercise. Therefore, IIBP training may be an effective accessory exercise to maintain a higher level of muscle activation across primary and stabilizing muscles with a lighter load for untrained men, while SBP may be a suitable bench press exercise for untrained participants who have not developed the neuromuscular adaptations necessary for correct stabilization of the elbow joint.

Disruption in structural-functional network repertoire and time-resolved subcortical fronto-temporoparietal connectivity in disorders of consciousness

Understanding recovery of consciousness and elucidating its underlying mechanism is believed to be crucial in the field of basic neuroscience and medicine. Ideas such as the global neuronal workspace (GNW) and the mesocircuit theory hypothesize that failure of recovery in conscious states coincide with loss of connectivity between subcortical and frontoparietal areas, a loss of the repertoire of functional networks states and metastable brain activation. We adopted a time-resolved functional connectivity framework to explore these ideas and assessed the repertoire of functional network states as a potential marker of consciousness and its potential ability to tell apart patients in the unresponsive wakefulness syndrome (UWS) and minimally conscious state (MCS). In addition, the prediction of these functional network states by underlying hidden spatial patterns in the anatomical network, that is so-called eigenmodes, was supplemented as potential markers. By analysing time-resolved functional connectivity from functional MRI data, we demonstrated a reduction of metastability and functional network repertoire in UWS compared to MCS patients. This was expressed in terms of diminished dwell times and loss of nonstationarity in the default mode network and subcortical fronto-temporoparietal network in UWS compared to MCS patients. We further demonstrated that these findings co-occurred with a loss of dynamic interplay between structural eigenmodes and emerging time-resolved functional connectivity in UWS. These results are, amongst others, in support of the GNW theory and the mesocircuit hypothesis, underpinning the role of time-resolved thalamo-cortical connections and metastability in the recovery of consciousness.

Rijke tube: A nonlinear oscillator

Dynamical systems theory has emerged as an interdisciplinary area of research to characterize the complex dynamical transitions in real-world systems. Various nonlinear dynamical phenomena and bifurcations have been discovered over the decades using different reduced-order models of oscillators. Different measures and methodologies have been developed theoretically to detect, control, or suppress the nonlinear oscillations. However, obtaining such phenomena experimentally is often challenging, time-consuming, and risky mainly due to the limited control of certain parameters during experiments. With this review, we aim to introduce a paradigmatic and easily configurable Rijke tube oscillator to the dynamical systems community. The Rijke tube is commonly used by the combustion community as a prototype to investigate the detrimental phenomena of thermoacoustic instability. Recent investigations in such Rijke tubes have utilized various methodologies from dynamical systems theory to better understand the occurrence of thermoacoustic oscillations and their prediction and mitigation, both experimentally and theoretically. The existence of various dynamical behaviors has been reported in single and coupled Rijke tube oscillators. These behaviors include bifurcations, routes to chaos, noise-induced transitions, synchronization, and suppression of oscillations. Various early warning measures have been established to predict thermoacoustic instabilities. Therefore, this review article consolidates the usefulness of a Rijke tube oscillator in terms of experimentally discovering and modeling different nonlinear phenomena observed in physics, thus transcending the boundaries between the physics and the engineering communities.

Spatiotemporal patterns of firearm acquisition in the United States in different presidential terms

This study develops mathematical tools and approaches to investigate spatiotemporal patterns of firearm acquisition in the U.S. complemented by hypothesis testing and statistical analysis. First, state-level and nation-level instant background check (BC) data are employed as proxy of firearm acquisition corresponding to 1999-2021. The relative-phase time-series of BC in each U.S. state is recovered and utilized to calculate the time-series of the U.S. states' synchronization degree. We reveal that U.S. states present a high-level degree of synchronization except in 2010-2011 and after 2018. Comparing these results with respect to a sitting U.S. president provides additional information: specifically, any two presidential terms are characterized by statistically different synchronization degrees except G. W. Bush's first term and B. H. Obama's second term. Next, to detail variations of BC, short-time Fourier transform, dimensionality reduction techniques, and diffusion maps are implemented within a time-frequency representation. Firearm acquisition in the high frequency band is described by a low-dimensional embedding, in the form of a plane with two embedding coordinates. Data points on the embedding plane identify separate clusters that signify state transitions in the original BC data with respect to different time windows. Through this analysis, we reveal that the frequency content of the BC data has a time-dependent characteristic. By comparing the diffusion map at hand with respect to a presidential term, we find that at least one of the embedding coordinates presents statistically significant variations between any two presidential terms except B. H. Obama's first term and D. J. Trump's pre-COVID term. The results point at a possible interplay between firearm acquisition in the U.S. and a presidential term.

May the 4C's be with you: an overview of complexity-inspired frameworks for analysing resting-state neuroimaging data

Competing and complementary models of resting-state brain dynamics contribute to our phenomenological and mechanistic understanding of whole-brain coordination and communication, and provide potential evidence for differential brain functioning associated with normal and pathological behaviour. These neuroscientific theories stem from the perspectives of physics, engineering, mathematics and psychology and create a complicated landscape of domain-specific terminology and meaning, which, when used outside of that domain, may lead to incorrect assumptions and conclusions within the neuroscience community. Here, we review and clarify the key concepts of connectivity, computation, criticality and coherence-the 4C's-and outline a potential role for metastability as a common denominator across these propositions. We analyse and synthesize whole-brain neuroimaging research, examined through functional magnetic imaging, to demonstrate that complexity science offers a principled and integrated approach to describe, and potentially understand, macroscale spontaneous brain functioning.

From Coherence to Multivariate Causal Estimators of EEG Connectivity

The paper concerns the development of methods of EEG functional connectivity estimation including short overview of the currently applied measures describing their advantages and flaws. Linear and non-linear, bivariate and multivariate methods are confronted. The performance of different connectivity measures in respect of robustness to noise, common drive effect and volume conduction is considered providing a guidance towards future developments in the field, which involve evaluation not only functional, but also effective (causal) connectivity. The time-varying connectivity measure making possible estimation of dynamical information processing in brain is presented. The methods of post-processing of connectivity results are considered involving application of advanced graph analysis taking into account community structure of networks and providing hierarchy of networks rather than the single, binary networks currently used.

Phase irregularity: A conceptually simple and efficient approach to characterize electroencephalographic recordings from epilepsy patients

The severe neurological disorder epilepsy affects almost 1% of the world population. For patients who suffer from pharmacoresistant focal-onset epilepsy, electroencephalographic (EEG) recordings are essential for the localization of the brain area where seizures start. Apart from the visual inspection of the recordings, quantitative EEG signal analysis techniques proved to be useful for this purpose. Among other features, regularity versus irregularity and phase coherence versus phase independence allowed characterizing brain dynamics from the measured EEG signals. Can phase irregularities also characterize brain dynamics? To address this question, we use the univariate coefficient of phase velocity variation, defined as the ratio of phase velocity standard deviation and the mean phase velocity. Beyond that, as a bivariate measure we use the classical mean phase coherence to quantify the degree of phase locking. All phase-based measures are combined with surrogates to test null hypotheses about the dynamics underlying the signals. In the first part of our analysis, we use the Rössler model system to study our approach under controlled conditions. In the second part, we use the Bern-Barcelona EEG database which consists of focal and nonfocal signals extracted from seizure-free recordings. Focal signals are recorded from brain areas where the first seizure EEG signal changes can be detected, and nonfocal signals are recorded from areas that are not involved in the seizure at its onset. Our results show that focal signals have less phase variability and more phase coherence than nonfocal signals. Once combined with surrogates, the mean phase velocity proved to have the highest discriminative power between focal and nonfocal signals. In conclusion, conceptually simple and easy to compute phase-based measures can help to detect features induced by epilepsy from EEG signals. This holds not only for the classical mean phase coherence but even more so for univariate measures of phase irregularity.

Halo Sport Transcranial Direct Current Stimulation Improved Muscular Endurance Performance and Neuromuscular Efficiency During an Isometric Submaximal Fatiguing Elbow Flexion Task

The present study examined the effects of transcranial direct current stimulation (tDCS) using Halo Sport on the time to exhaustion (TTE) in relation with muscle activities and corticomuscular coupling of agonist and antagonist muscles during a sustained isometric fatiguing contraction performed with the elbow flexors. Twenty healthy male college students were randomly assigned to tDCS group and control group. The two group participants performed two experimental sessions which consisted of pre-fatigue isometric maximal voluntary contraction (MVC), sustained submaximal voluntary contractions (30% maximal torque) performed to exhaustion, and post-fatigue MVC with the right elbow flexor muscles. Sham stimulation (90 s) and tDCS (20 min) were applied for control and tDCS group participants 20 min prior to the second session test, respectively. MVC strength in pre- and post-fatigue test, TTE, electroencephalogram (EEG), and electromyography (EMG) of biceps brachii (BB) and triceps brachii (TB) were recorded during the tests. It was found that tDCS using the Halo Sport device significantly increased TTE and thus improved muscular endurance performance. The improvement may be partly related to the improvement of neuromuscular efficiency as reflected by decrease of antagonistic muscle coactivation activities, which may be related to cortical originated central processing mechanism of neuromuscular activities.

Monitoring the Evolution of Asynchrony between Mean Arterial Pressure and Mean Cerebral Blood Flow via Cross-Entropy Methods

Cerebrovascular control is carried out by multiple nonlinear mechanisms imposing a certain degree of coupling between mean arterial pressure (MAP) and mean cerebral blood flow (MCBF). We explored the ability of two nonlinear tools in the information domain, namely cross-approximate entropy (CApEn) and cross-sample entropy (CSampEn), to assess the degree of asynchrony between the spontaneous fluctuations of MAP and MCBF. CApEn and CSampEn were computed as a function of the translation time. The analysis was carried out in 23 subjects undergoing recordings at rest in supine position (REST) and during active standing (STAND), before and after surgical aortic valve replacement (SAVR). We found that at REST the degree of asynchrony raised, and the rate of increase in asynchrony with the translation time decreased after SAVR. These results are likely the consequence of the limited variability of MAP observed after surgery at REST, more than the consequence of a modified cerebrovascular control, given that the observed differences disappeared during STAND. CApEn and CSampEn can be utilized fruitfully in the context of the evaluation of cerebrovascular control via the noninvasive acquisition of the spontaneous MAP and MCBF variability.

Directional couplings between the respiration and parasympathetic control of the heart rate during sleep and wakefulness in healthy subjects at different ages

Cardiorespiratory interactions are important, both for understanding the fundamental processes of functioning of the human body and for development of methods for diagnostics of various pathologies. The properties of cardiorespiratory interaction are determined by the processes of autonomic control of blood circulation, which are modulated by the higher nervous activity. We study the directional couplings between the respiration and the process of parasympathetic control of the heart rate in the awake state and different stages of sleep in 96 healthy subjects from different age groups. The detection of directional couplings is carried out using the method of phase dynamics modeling applied to experimental RR-intervals and the signal of respiration. We reveal the presence of bidirectional couplings between the studied processes in all age groups. Our results show that the coupling from respiration to the process of parasympathetic control of the heart rate is stronger than the coupling in the opposite direction. The difference in the strength of bidirectional couplings between the considered processes is most pronounced in deep sleep.

Connectivity Analysis for Multivariate Time Series: Correlation vs. Causality

The study of the interdependence relationships of the variables of an examined system is of great importance and remains a challenging task. There are two distinct cases of interdependence. In the first case, the variables evolve in synchrony, connections are undirected and the connectivity is examined based on symmetric measures, such as correlation. In the second case, a variable drives another one and they are connected with a causal relationship. Therefore, directed connections entail the determination of the interrelationships based on causality measures. The main open question that arises is the following: can symmetric correlation measures or directional causality measures be applied to infer the connectivity network of an examined system? Using simulations, we demonstrate the performance of different connectivity measures in case of contemporaneous or/and temporal dependencies. Results suggest the sensitivity of correlation measures when temporal dependencies exist in the data. On the other hand, causality measures do not spuriously indicate causal effects when data present only contemporaneous dependencies. Finally, the necessity of introducing effective instantaneous causality measures is highlighted since they are able to handle both contemporaneous and causal effects at the same time. Results based on instantaneous causality measures are promising; however, further investigation is required in order to achieve an overall satisfactory performance.

Classifying subclinical depression using EEG spectral and connectivity measures

Detecting depression on its early stages helps preventing the onset of severe depressive episodes. In this study, we propose an automatic classification pipeline to detect subclinical depression (i.e., dysphoria) through the electroencephalography (EEG) signal. To this aim, we recorded the EEG signals in resting condition from 26 female participants with dysphoria and 38 female controls. The EEG signals were processed to extract several spectral and functional connectivity features to feed a nonlinear Support Vector Machine (SVM) classifier embedded with a Recursive Feature Elimination (RFE) algorithm. Our recognition pipeline obtained a maximum classification accuracy of 83.91% in recognizing dysphoria patients with a combination of connectivity and spectral measures. Moreover, an accuracy of 76.11% was achieved with only the 4 most informative functional connections, suggesting a central role of cortical connectivity in the theta band for early depression recognition. The present study can facilitate the diagnosis of subclinical conditions of depression and may provide reliable indicators of depression for the clinical community.

Comparison of undirected frequency-domain connectivity measures for cerebro-peripheral analysis

Analyses of cerebro-peripheral connectivity aim to quantify ongoing coupling between brain activity (measured by MEG/EEG) and peripheral signals such as muscle activity, continuous speech, or physiological rhythms (such as pupil dilation or respiration). Due to the distinct rhythmicity of these signals, undirected connectivity is typically assessed in the frequency domain. This leaves the investigator with two critical choices, namely a) the appropriate measure for spectral estimation (i.e., the transformation into the frequency domain) and b) the actual connectivity measure. As there is no consensus regarding best practice, a wide variety of methods has been applied. Here we systematically compare combinations of six standard spectral estimation methods (comprising fast Fourier and continuous wavelet transformation, bandpass filtering, and short-time Fourier transformation) and six connectivity measures (phase-locking value, Gaussian-Copula mutual information, Rayleigh test, weighted pairwise phase consistency, magnitude squared coherence, and entropy). We provide performance measures of each combination for simulated data (with precise control over true connectivity), a single-subject set of real MEG data, and a full group analysis of real MEG data. Our results show that, overall, WPPC and GCMI tend to outperform other connectivity measures, while entropy was the only measure sensitive to bimodal deviations from a uniform phase distribution. For group analysis, choosing the appropriate spectral estimation method appears to be more critical than the connectivity measure. We discuss practical implications (sampling rate, SNR, computation time, and data length) and aim to provide recommendations tailored to particular research questions.

Windowed multiscale synchrony: modeling time-varying and scale-localized interpersonal coordination dynamics

Social interactions are pervasive in human life with varying forms of interpersonal coordination emerging and spanning different modalities (e.g. behaviors, speech/language, and neurophysiology). However, during social interactions, as in any dynamical system, patterns of coordination form and dissipate at different scales. Historically, researchers have used aggregate measures to capture coordination over time. While those measures (e.g. mean relative phase, cross-correlation, coherence) have provided a wealth of information about coordination in social settings, some evidence suggests that multiscale coordination may change over the time course of a typical empirical observation. To address this gap, we demonstrate an underutilized method, windowed multiscale synchrony, that moves beyond quantifying aggregate measures of coordination by focusing on how the relative strength of coordination changes over time and the scales that comprise social interaction. This method involves using a wavelet transform to decompose time series into component frequencies (i.e. scales), preserving temporal information and then quantifying phase synchronization at each of these scales. We apply this method to both simulated and empirical interpersonal physiological and neuromechanical data. We anticipate that demonstrating this method will stimulate new insights on the mechanisms and functions of synchrony in interpersonal contexts using neurophysiological and behavioral measures.

Evaluating the sensitivity of functional connectivity measures to motion artifact in resting-state fMRI data

Functional connectivity (FC) networks are typically inferred from resting-state fMRI data using the Pearson correlation between BOLD time series from pairs of brain regions. However, alternative methods of estimating functional connectivity have not been systematically tested for their sensitivity or robustness to head motion artifact. Here, we evaluate the sensitivity of eight different functional connectivity measures to motion artifact using resting-state data from the Human Connectome Project. We report that FC estimated using full correlation has a relatively high residual distance-dependent relationship with motion compared to partial correlation, coherence, and information theory-based measures, even after implementing rigorous methods for motion artifact mitigation. This disadvantage of full correlation, however, may be offset by higher test-retest reliability, fingerprinting accuracy, and system identifiability. FC estimated by partial correlation offers the best of both worlds, with low sensitivity to motion artifact and intermediate system identifiability, with the caveat of low test-retest reliability and fingerprinting accuracy. We highlight spatial differences in the sub-networks affected by motion with different FC metrics. Further, we report that intra-network edges in the default mode and retrosplenial temporal sub-networks are highly correlated with motion in all FC methods. Our findings indicate that the method of estimating functional connectivity is an important consideration in resting-state fMRI studies and must be chosen carefully based on the parameters of the study.

Super-Selective Reconstruction of Causal and Direct Connectivity With Application to in vitro iPSC Neuronal Networks

Despite advancements in the development of cell-based in-vitro neuronal network models, the lack of appropriate computational tools limits their analyses. Methods aimed at deciphering the effective connections between neurons from extracellular spike recordings would increase utility of in vitro local neural circuits, especially for studies of human neural development and disease based on induced pluripotent stem cells (hiPSC). Current techniques allow statistical inference of functional couplings in the network but are fundamentally unable to correctly identify indirect and apparent connections between neurons, generating redundant maps with limited ability to model the causal dynamics of the network. In this paper, we describe a novel mathematically rigorous, model-free method to map effective-direct and causal-connectivity of neuronal networks from multi-electrode array data. The inference algorithm uses a combination of statistical and deterministic indicators which, first, enables identification of all existing functional links in the network and then reconstructs the directed and causal connection diagram via a super-selective rule enabling highly accurate classification of direct, indirect, and apparent links. Our method can be generally applied to the functional characterization of any in vitro neuronal networks. Here, we show that, given its accuracy, it can offer important insights into the functional development of in vitro hiPSC-derived neuronal cultures.

Low-frequency component of photoplethysmogram reflects the autonomic control of blood pressure

The question of how much information the photoplethysmogram (PPG) signal contains on the autonomic regulation of blood pressure (BP) remains unsolved. This study aims to compare the low-frequency (LF) and high-frequency components of PPG and BP and assess their correlation with oscillations in interbeat (RR) intervals at similar frequencies. The PPG signal from the distal phalanx of the right index finger recorded using a reflective PPG sensor at green light, the BP signal from the left hand recorded using a Finometer, and RR intervals were analyzed. These signals were simultaneously recorded within 15 min in a supine resting condition in 17 healthy subjects (12 males and 5 females) aged 33 ± 9 years (mean ± SD). The study revealed the high coherence of LF components of PPG and BP with the LF component of RR intervals. The high-frequency components of these signals had low coherence. The analysis of the signal instantaneous phases revealed the presence of high-phase coherence between the LF components of PPG and BP. It is shown that the LF component of PPG is determined not only by local myogenic activity but also reflects the processes of autonomic control of BP.