Distinguishing between direct and indirect directional couplings in large oscillator networks: Partial or non-partial phase analyses?

Inferring connectivity of an oscillatory network via the phase dynamics reconstruction

We review an approach for reconstructing oscillatory networks' undirected and directed connectivity from data. The technique relies on inferring the phase dynamics model. The central assumption is that we observe the outputs of all network nodes. We distinguish between two cases. In the first one, the observed signals represent smooth oscillations, while in the second one, the data are pulse-like and can be viewed as point processes. For the first case, we discuss estimating the true phase from a scalar signal, exploiting the protophase-to-phase transformation. With the phases at hand, pairwise and triplet synchronization indices can characterize the undirected connectivity. Next, we demonstrate how to infer the general form of the coupling functions for two or three oscillators and how to use these functions to quantify the directional links. We proceed with a different treatment of networks with more than three nodes. We discuss the difference between the structural and effective phase connectivity that emerges due to high-order terms in the coupling functions. For the second case of point-process data, we use the instants of spikes to infer the phase dynamics model in the Winfree form directly. This way, we obtain the network's coupling matrix in the first approximation in the coupling strength.

Ordinal methods for a characterization of evolving functional brain networks

Ordinal time series analysis is based on the idea to map time series to ordinal patterns, i.e., order relations between the values of a time series and not the values themselves, as introduced in 2002 by C. Bandt and B. Pompe. Despite a resulting loss of information, this approach captures meaningful information about the temporal structure of the underlying system dynamics as well as about properties of interactions between coupled systems. This-together with its conceptual simplicity and robustness against measurement noise-makes ordinal time series analysis well suited to improve characterization of the still poorly understood spatiotemporal dynamics of the human brain. This minireview briefly summarizes the state-of-the-art of uni- and bivariate ordinal time-series-analysis techniques together with applications in the neurosciences. It will highlight current limitations to stimulate further developments, which would be necessary to advance characterization of evolving functional brain networks.

Network structure from a characterization of interactions in complex systems

Many natural and man-made complex dynamical systems can be represented by networks with vertices representing system units and edges the coupling between vertices. If edges of such a structural network are inaccessible, a widely used approach is to identify them with interactions between vertices, thereby setting up a functional network. However, it is an unsolved issue if and to what extent important properties of a functional network on the global and the local scale match those of the corresponding structural network. We address this issue by deriving functional networks from characterizing interactions in paradigmatic oscillator networks with widely-used time-series-analysis techniques for various factors that alter the collective network dynamics. Surprisingly, we find that particularly key constituents of functional networks—as identified with betweenness and eigenvector centrality—coincide with ground truth to a high degree, while global topological and spectral properties—clustering coefficient, average shortest path length, assortativity, and synchronizability—clearly deviate. We obtain similar concurrences for an empirical network. Our findings are of relevance for various scientific fields and call for conceptual and methodological refinements to further our understanding of the relationship between structure and function of complex dynamical systems.

Partial event coincidence analysis for distinguishing direct and indirect coupling in functional network construction

Correctly identifying interaction patterns from multivariate time series presents an important step in functional network construction. In this context, the widespread use of bivariate statistical association measures often results in a false identification of links because strong similarity between two time series can also emerge without the presence of a direct interaction due to intermediate mediators or common drivers. In order to properly distinguish such direct and indirect links for the special case of event-like data, we present here a new generalization of event coincidence analysis to a partial version thereof, which is aimed at excluding possible transitive effects of indirect couplings. Using coupled chaotic systems and stochastic processes on two generic coupling topologies (star and chain configuration), we demonstrate that the proposed methodology allows for the correct identification of indirect interactions. Subsequently, we apply our partial event coincidence analysis to multi-channel EEG recordings to investigate possible differences in coordinated alpha band activity among macroscopic brain regions in resting states with eyes open (EO) and closed (EC) conditions. Specifically, we find that direct connections typically correspond to close spatial neighbors while indirect ones often reflect longer-distance connections mediated via other brain regions. In the EC state, connections in the frontal parts of the brain are enhanced as compared to the EO state, while the opposite applies to the posterior regions. In general, our approach leads to a significant reduction in the number of indirect connections and thereby contributes to a better understanding of the alpha band desynchronization phenomenon in the EO state.

Impact of local network characteristics on network reconstruction

When a network is inferred from data, two types of errors can occur: false positive and false negative conclusions about the presence of links. We focus on the influence of local network characteristics on the probability α of false positive conclusions, and on the probability β of false negative conclusions, in the case of networks of coupled oscillators. We demonstrate that false conclusion probabilities are influenced by local connectivity measures such as the shortest path length and the detour degree, which can also be estimated from the inferred network when the true underlying network is not known a priori. These measures can then be used for quantification of the confidence level of link conclusions, and for improving the network reconstruction via advanced concepts of link weights thresholding.

The Human Organism as an Integrated Interaction Network: Recent Conceptual and Methodological Challenges

The field of Network Physiology aims to advance our understanding of how physiological systems and sub-systems interact to generate a variety of behaviors and distinct physiological states, to optimize the organism's functioning, and to maintain health. Within this framework, which considers the human organism as an integrated network, vertices are associated with organs while edges represent time-varying interactions between vertices. Likewise, vertices may represent networks on smaller spatial scales leading to a complex mixture of interacting homogeneous and inhomogeneous networks of networks. Lacking adequate analytic tools and a theoretical framework to probe interactions within and among diverse physiological systems, current approaches focus on inferring properties of time-varying interactions-namely strength, direction, and functional form-from time-locked recordings of physiological observables. To this end, a variety of bivariate or, in general, multivariate time-series-analysis techniques, which are derived from diverse mathematical and physical concepts, are employed and the resulting time-dependent networks can then be further characterized with methods from network theory. Despite the many promising new developments, there are still problems that evade from a satisfactory solution. Here we address several important challenges that could aid in finding new perspectives and inspire the development of theoretic and analytical concepts to deal with these challenges and in studying the complex interactions between physiological systems.

Evaluation of Granger Causality Measures for Constructing Networks from Multivariate Time Series

Granger causality and variants of this concept allow the study of complex dynamical systems as networks constructed from multivariate time series. In this work, a large number of Granger causality measures used to form causality networks from multivariate time series are assessed. These measures are in the time domain, such as model-based and information measures, the frequency domain, and the phase domain. The study aims also to compare bivariate and multivariate measures, linear and nonlinear measures, as well as the use of dimension reduction in linear model-based measures and information measures. The latter is particular relevant in the study of high-dimensional time series. For the performance of the multivariate causality measures, low and high dimensional coupled dynamical systems are considered in discrete and continuous time, as well as deterministic and stochastic. The measures are evaluated and ranked according to their ability to provide causality networks that match the original coupling structure. The simulation study concludes that the Granger causality measures using dimension reduction are superior and should be preferred particularly in studies involving many observed variables, such as multi-channel electroencephalograms and financial markets.

Impact of type of intracranial EEG sensors on link strengths of evolving functional brain networks

Objective and Approach: Investigating properties of evolving functional brain networks has become a valuable tool to characterize the complex dynamics of the epileptic brain. Such networks are usually derived from electroencephalograms (EEG) recorded with sensors implanted chronically into deeper structures of the brain and/or placed onto the cortex. It is still unclear, however, whether the use of different sensors for an identification of network nodes affects properties of functional brain networks. We address this question by investigating properties of links of such networks that we characterize by assessing interactions in multi-sensor, multi-day EEG data recorded from 49 epilepsy patients during presurgical evaluation. These data allow us to study the impact of different types of sensors together with the impact of various physiologic and pathophysiologic activities on the properties of links.

MAIN RESULTS

We observe that different types of sensors differently impact on spatial means and temporal fluctuations of link strengths. Moreover, the impact depends on the relative anatomical location of sensors with respect to location and extent of sources of the prevailing activities.

SIGNIFICANCE

Type and location of sensors should be considered when constructing networks.

Reconstructing networks of pulse-coupled oscillators from spike trains

We present an approach for reconstructing networks of pulse-coupled neuronlike oscillators from passive observation of pulse trains of all nodes. It is assumed that units are described by their phase response curves and that their phases are instantaneously reset by incoming pulses. Using an iterative procedure, we recover the properties of all nodes, namely their phase response curves and natural frequencies, as well as strengths of all directed connections.

Neural Cross-Frequency Coupling Functions

Although neural interactions are usually characterized only by their coupling strength and directionality, there is often a need to go beyond this by establishing the functional mechanisms of the interaction. We introduce the use of dynamical Bayesian inference for estimation of the coupling functions of neural oscillations in the presence of noise. By grouping the partial functional contributions, the coupling is decomposed into its functional components and its most important characteristics-strength and form-are quantified. The method is applied to characterize the δ-to-α phase-to-phase neural coupling functions from electroencephalographic (EEG) data of the human resting state, and the differences that arise when the eyes are either open (EO) or closed (EC) are evaluated. The δ-to-α phase-to-phase coupling functions were reconstructed, quantified, compared, and followed as they evolved in time. Using phase-shuffled surrogates to test for significance, we show how the strength of the direct coupling, and the similarity and variability of the coupling functions, characterize the EO and EC states for different regions of the brain. We confirm an earlier observation that the direct coupling is stronger during EC, and we show for the first time that the coupling function is significantly less variable. Given the current understanding of the effects of e.g., aging and dementia on δ-waves, as well as the effect of cognitive and emotional tasks on α-waves, one may expect that new insights into the neural mechanisms underlying certain diseases will be obtained from studies of coupling functions. In principle, any pair of coupled oscillations could be studied in the same way as those shown here.