Surrogate-assisted analysis of weighted functional brain networks

Time-series-analysis-based detection of critical transitions in real-world non-autonomous systems

Real-world non-autonomous systems are open, out-of-equilibrium systems that evolve in and are driven by temporally varying environments. Such systems can show multiple timescale and transient dynamics together with transitions to very different and, at times, even disastrous dynamical regimes. Since such critical transitions disrupt the systems' intended or desired functionality, it is crucial to understand the underlying mechanisms, to identify precursors of such transitions, and to reliably detect them in time series of suitable system observables to enable forecasts. This review critically assesses the various steps of investigation involved in time-series-analysis-based detection of critical transitions in real-world non-autonomous systems: from the data recording to evaluating the reliability of offline and online detections. It will highlight pros and cons to stimulate further developments, which would be necessary to advance understanding and forecasting nonlinear behavior such as critical transitions in complex systems.

The time-evolving epileptic brain network: concepts, definitions, accomplishments, perspectives

Epilepsy is now considered a network disease that affects the brain across multiple levels of spatial and temporal scales. The paradigm shift from an epileptic focus-a discrete cortical area from which seizures originate-to a widespread epileptic network-spanning lobes and hemispheres-considerably advanced our understanding of epilepsy and continues to influence both research and clinical treatment of this multi-faceted high-impact neurological disorder. The epileptic network, however, is not static but evolves in time which requires novel approaches for an in-depth characterization. In this review, we discuss conceptual basics of network theory and critically examine state-of-the-art recording techniques and analysis tools used to assess and characterize a time-evolving human epileptic brain network. We give an account on current shortcomings and highlight potential developments towards an improved clinical management of epilepsy.

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.

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.

FitzHugh-Nagumo oscillators on complex networks mimic epileptic-seizure-related synchronization phenomena

We study patterns of partial synchronization in a network of FitzHugh-Nagumo oscillators with empirical structural connectivity measured in human subjects. We report the spontaneous occurrence of synchronization phenomena that closely resemble the ones seen during epileptic seizures in humans. In order to obtain deeper insights into the interplay between dynamics and network topology, we perform long-term simulations of oscillatory dynamics on different paradigmatic network structures: random networks, regular nonlocally coupled ring networks, ring networks with fractal connectivities, and small-world networks with various rewiring probability. Among these networks, a small-world network with intermediate rewiring probability best mimics the findings achieved with the simulations using the empirical structural connectivity. For the other network topologies, either no spontaneously occurring epileptic-seizure-related synchronization phenomena can be observed in the simulated dynamics, or the overall degree of synchronization remains high throughout the simulation. This indicates that a topology with some balance between regularity and randomness favors the self-initiation and self-termination of episodes of seizure-like strong synchronization.

Von Interaktionen zu Interaktionsnetzwerken: Zeitabhängige funktionelle Netzwerke am Beispiel der Epilepsie

Das menschliche Gehirn ist ein komplexes Netzwerk aus interagierenden nichtstationären Subsystemen (Netzwerk von Netzwerken), deren komplizierte räumlich-zeitliche Dynamiken bis heute nur unzureichend verstanden sind. Dabei versprechen aktuelle Entwicklungen im Bereich der Zeitreihenanalyse sowie der Theorie komplexer Netzwerke neue und verbesserte Einblicke in die Dynamiken von Hirnnetzwerken auf verschiedenen räumlich-zeitlichen Skalen. Wir geben einen Überblick über diese Entwicklungen und besprechen am Beispiel zeitabhängiger epileptischer Hirnnetzwerke Fortschritte im Verständnis von Hirndynamiken, die über multiple Skalen hinweg variieren.

Precursors of seizures due to specific spatial-temporal modifications of evolving large-scale epileptic brain networks

Knowing when, where, and how seizures are initiated in large-scale epileptic brain networks remains a widely unsolved problem. Seizure precursors – changes in brain dynamics predictive of an impending seizure – can now be identified well ahead of clinical manifestations, but either the seizure onset zone or remote brain areas are reported as network nodes from which seizure precursors emerge. We aimed to shed more light on the role of constituents of evolving epileptic networks that recurrently transit into and out of seizures. We constructed such networks from more than 3200 hours of continuous intracranial electroencephalograms recorded in 38 patients with medication refractory epilepsy. We succeeded in singling out predictive edges and predictive nodes. Their particular characteristics, namely edge weight respectively node centrality (a fundamental concept of network theory), from the pre-ictal periods of 78 out of 97 seizures differed significantly from the characteristics seen during inter-ictal periods. The vast majority of predictive nodes were connected by most of the predictive edges, but these nodes never played a central role in the evolving epileptic networks. Interestingly, predictive nodes were entirely associated with brain regions deemed unaffected by the focal epileptic process. We propose a network mechanism for a transition into the pre-seizure state, which puts into perspective the role of the seizure onset zone in this transition and highlights the necessity to reassess current concepts for seizure generation and seizure prevention.

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.

Quantification of Graph Complexity Based on the Edge Weight Distribution Balance: Application to Brain Networks

The aim of this study was to introduce a novel global measure of graph complexity: Shannon graph complexity (SGC). This measure was specifically developed for weighted graphs, but it can also be applied to binary graphs. The proposed complexity measure was designed to capture the interplay between two properties of a system: the ‘information’ (calculated by means of Shannon entropy) and the ‘order’ of the system (estimated by means of a disequilibrium measure). SGC is based on the concept that complex graphs should maintain an equilibrium between the aforementioned two properties, which can be measured by means of the edge weight distribution. In this study, SGC was assessed using four synthetic graph datasets and a real dataset, formed by electroencephalographic (EEG) recordings from controls and schizophrenia patients. SGC was compared with graph density (GD), a classical measure used to evaluate graph complexity. Our results showed that SGC is invariant with respect to GD and independent of node degree distribution. Furthermore, its variation with graph size [Formula: see text] is close to zero for [Formula: see text]. Results from the real dataset showed an increment in the weight distribution balance during the cognitive processing for both controls and schizophrenia patients, although these changes are more relevant for controls. Our findings revealed that SGC does not need a comparison with null-hypothesis networks constructed by a surrogate process. In addition, SGC results on the real dataset suggest that schizophrenia is associated with a deficit in the brain dynamic reorganization related to secondary pathways of the brain network.

Surrogate-assisted identification of influences of network construction on evolving weighted functional networks

We aim at identifying factors that may affect the characteristics of evolving weighted networks derived from empirical observations. To this end, we employ various chains of analysis that are often used in field studies for a data-driven derivation and characterization of such networks. As an example, we consider fully connected, weighted functional brain networks before, during, and after epileptic seizures that we derive from multichannel electroencephalographic data recorded from epilepsy patients. For these evolving networks, we estimate clustering coefficient and average shortest path length in a time-resolved manner. Lastly, we make use of surrogate concepts that we apply at various levels of the chain of analysis to assess to what extent network characteristics are dominated by properties of the electroencephalographic recordings and/or the evolving weighted networks, which may be accessible more easily. We observe that characteristics are differently affected by the unavoidable referencing of the electroencephalographic recording, by the time-series-analysis technique used to derive the properties of network links, and whether or not networks were normalized. Importantly, for the majority of analysis settings, we observe temporal evolutions of network characteristics to merely reflect the temporal evolutions of mean interaction strengths. Such a property of the data may be accessible more easily, which would render the weighted network approach-as used here-as an overly complicated description of simple aspects of the data.

Long-term variability of importance of brain regions in evolving epileptic brain networks

We investigate the temporal and spatial variability of the importance of brain regions in evolving epileptic brain networks. We construct these networks from multiday, multichannel electroencephalographic data recorded from 17 epilepsy patients and use centrality indices to assess the importance of brain regions. Time-resolved indications of highest importance fluctuate over time to a greater or lesser extent, however, with some periodic temporal structure that can mostly be attributed to phenomena unrelated to the disease. In contrast, relevant aspects of the epileptic process contribute only marginally. Indications of highest importance also exhibit pronounced alternations between various brain regions that are of relevance for studies aiming at an improved understanding of the epileptic process with graph-theoretical approaches. Nonetheless, these findings may guide new developments for individualized diagnosis, treatment, and control.

Which Brain Regions are Important for Seizure Dynamics in Epileptic Networks? Influence of Link Identification and EEG Recording Montage on Node Centralities

Nodes in large-scale epileptic networks that are crucial for seizure facilitation and termination can be regarded as potential targets for individualized focal therapies. Graph-theoretical approaches based on centrality concepts can help to identify such important nodes, however, they may be influenced by the way networks are derived from empirical data. Here we investigate evolving functional epileptic brain networks during 82 focal seizures with different anatomical onset locations that we derive from multichannel intracranial electroencephalographic recordings from 51 patients. We demonstrate how the various methodological steps (from the recording montage via node and link inference to the assessment of node centralities) affect importance estimation and discuss their impact on the interpretability of findings in the context of pathophysiological aspects of seizure dynamics.

Are Structural Changes Induced by Lithium in the HIV Brain Accompanied by Changes in Functional Connectivity?

Lithium therapy has been shown to affect imaging measures of brain function and microstructure in human immunodeficiency virus (HIV)-infected subjects with cognitive impairment. The aim of this proof-of-concept study was to explore whether changes in brain microstructure also entail changes in functional connectivity. Functional MRI data of seven cognitively impaired HIV infected individuals enrolled in an open-label lithium study were included in the connectivity analysis. Seven regions of interest (ROI) were defined based on previously observed lithium induced microstructural changes measured by Diffusion Tensor Imaging. Generalized partial directed coherence (gPDC), based on time-variant multivariate autoregressive models, was used to quantify the degree of connectivity between the selected ROIs. Statistical analyses using a linear mixed model showed significant differences in the average node strength between pre and post lithium therapy conditions. Specifically, we found that lithium treatment in this population induced changes suggestive of increased strength in functional connectivity. Therefore, by exploiting the information about the strength of functional interactions provided by gPDC we can quantify the connectivity changes observed in relation to a given intervention. Furthermore, in conditions where the intervention is associated with clinical changes, we suggest that this methodology could enable an interpretation of such changes in the context of disease or treatment induced modulations in functional networks.

Incidental and intentional learning of verbal episodic material differentially modifies functional brain networks

Learning- and memory-related processes are thought to result from dynamic interactions in large-scale brain networks that include lateral and mesial structures of the temporal lobes. We investigate the impact of incidental and intentional learning of verbal episodic material on functional brain networks that we derive from scalp-EEG recorded continuously from 33 subjects during a neuropsychological test schedule. Analyzing the networks' global statistical properties we observe that intentional but not incidental learning leads to a significantly increased clustering coefficient, and the average shortest path length remains unaffected. Moreover, network modifications correlate with subsequent recall performance: the more pronounced the modifications of the clustering coefficient, the higher the recall performance. Our findings provide novel insights into the relationship between topological aspects of functional brain networks and higher cognitive functions.

On the discovery of group-consistent graph substructure patterns from brain networks

Complex networks constitute a recurring issue in the analysis of neuroimaging data. Recently, network motifs have been identified as patterns of interconnections since they appear in a significantly higher number than in randomized networks, in a given ensemble of anatomical or functional connectivity graphs. The current approach for detecting and enumerating motifs in brain networks requires a predetermined motif repertoire and can operate only with motifs of small size (consisting of few nodes). There is a growing interest in methodologies for frequent graph-based pattern mining in large graph datasets that can facilitate adaptive design of motifs. The results presented in this paper are based on the graph-based Substructure pattern mining (gSpan) algorithm and introduce a manifold of ways to exploit it for data-driven motif extraction in connectomics research. Functional connectivity graphs from electroencephalographic (EEG) recordings during resting state and mental calculations are used to demonstrate our approach. Relying on either time-invariant or time-evolving graphs, characteristic motifs associated with various frequency bands were derived and compared. With a suitable manipulation, the gSpan discovers motifs which are specific to performing mental arithmetics. Finally, the subject-dependent temporal signatures of motifs' appearance revealed the transient nature of the evolving functional connectivity (math-related motifs "come and go").