Identification of interacting neural populations: methods and statistical considerations

As improved recording technologies have created new opportunities for neurophysiological investigation, emphasis has shifted from individual neurons to multiple populations that form circuits, and it has become important to provide evidence of cross-population coordinated activity. We review various methods for doing so, placing them in six major categories while avoiding technical descriptions and instead focusing on high-level motivations and concerns. Our aim is to indicate what the methods can achieve and the circumstances under which they are likely to succeed. Toward this end, we include a discussion of four cross-cutting issues: the definition of neural populations, trial-to-trial variability and Poisson-like noise, time-varying dynamics, and causality.

Algebraic approach to the Kuramoto model

We study the Kuramoto model with attractive sine coupling. We introduce a complex-valued matrix formulation whose argument coincides with the original Kuramoto dynamics. We derive an exact solution for the complex-valued model, which permits analytical insight into individual realizations of the Kuramoto model. The existence of a complex-valued form of the Kuramoto model provides a key demonstration that, in some cases, reformulations of nonlinear dynamics in higher-order number fields may provide tractable analytical approaches.

Spatial correlations in a finite-range Kuramoto model

We study spatial correlations between oscillator phases in the steady state of a Kuramoto model, in which phase oscillators that are randomly distributed in space interact with constant strength but within a limited range. Such a model could be relevant, for example, in the synchronization of gene expression oscillations in cells, where only oscillations of neighboring cells are coupled through cell-cell contacts. We analytically infer spatial phase-phase correlation functions from the known steady-state distribution of oscillators for the case of homogenous frequencies and show that these can contain information about the range and strength of interactions, provided the noise in the system can be estimated. We suggest a method for the latter, and also explore when correlations appear to be ergodic in this model, which would enable an experimental measurement of correlation functions to utilize temporal averages. Simulations show that our techniques also give qualitative results for the model with heterogenous frequencies. We illustrate our results by comparison with experimental data on genetic oscillations in the segmentation clock of zebrafish embryos.

TORUS GRAPHS FOR MULTIVARIATE PHASE COUPLING ANALYSIS

Angular measurements are often modeled as circular random variables, where there are natural circular analogues of moments, including correlation. Because a product of circles is a torus, a d-dimensional vector of circular random variables lies on a d-dimensional torus. For such vectors we present here a class of graphical models, which we call torus graphs, based on the full exponential family with pairwise interactions. The topological distinction between a torus and Euclidean space has several important consequences. Our development was motivated by the problem of identifying phase coupling among oscillatory signals recorded from multiple electrodes in the brain: oscillatory phases across electrodes might tend to advance or recede together, indicating coordination across brain areas. The data analyzed here consisted of 24 phase angles measured repeatedly across 840 experimental trials (replications) during a memory task, where the electrodes were in 4 distinct brain regions, all known to be active while memories are being stored or retrieved. In realistic numerical simulations, we found that a standard pairwise assessment, known as phase locking value, is unable to describe multivariate phase interactions, but that torus graphs can accurately identify conditional associations. Torus graphs generalize several more restrictive approaches that have appeared in various scientific literatures, and produced intuitive results in the data we analyzed. Torus graphs thus unify multivariate analysis of circular data and present fertile territory for future research.

A practical method for estimating coupling functions in complex dynamical systems

A foremost challenge in modern network science is the inverse problem of reconstruction (inference) of coupling equations and network topology from the measurements of the network dynamics. Of particular interest are the methods that can operate on real (empirical) data without interfering with the system. One such earlier attempt (Tokuda et al. 2007 Phys. Rev. Lett. 99, 064101. (doi:10.1103/PhysRevLett.99.064101)) was a method suited for general limit-cycle oscillators, yielding both oscillators' natural frequencies and coupling functions between them (phase equations) from empirically measured time series. The present paper reviews the above method in a way comprehensive to domain-scientists other than physics. It also presents applications of the method to (i) detection of the network connectivity, (ii) inference of the phase sensitivity function, (iii) approximation of the interaction among phase-coherent chaotic oscillators, and (iv) experimental data from a forced Van der Pol electric circuit. This reaffirms the range of applicability of the method for reconstructing coupling functions and makes it accessible to a much wider scientific community. This article is part of the theme issue 'Coupling functions: dynamical interaction mechanisms in the physical, biological and social sciences'.

Synchronization of Independent Neural Ensembles in Human EEG during Choice Tasks

During behavioral experiments, humans placed in a situation of having to choose between a more valuable but risky reward and a less valuable but guaranteed reward make their decisions in accordance with external situational factors and individual characteristics, such as inclination to risk or caution. In such situations, humans can be divided into “risk-inclined” and “risk-averse” (or “cautious”) subjects. In this work, characteristics of EEG rhythms, such as phase–phase relationships and time lags between rhythms, were studied in pairs of alpha–beta and theta–beta rhythms. Phase difference can also be expressed as a time lag. It has been suggested that statistically significant time lags between rhythms are due to the combined neural activity of anatomically separate, independent (in activation/inhibition processes) ensembles. The extents of synchronicity between rhythms were compared as percentages between risk-inclined and risk-averse subjects. The results showed that synchronicity in response to stimuli was more often observed in pairs of alpha–beta rhythms of risk-averse subjects compared with risk-inclined subjects during the choice of a more valuable but less probable reward. In addition, significant differences in the percentage ratio of alpha and beta rhythms were revealed between (i) cases of synchronization without long time lags and (ii) cases with long time lags between rhythms (from 0.08 to 0.1 s).

Aberrant MEG multi-frequency phase temporal synchronization predicts conversion from mild cognitive impairment-to-Alzheimer's disease

Many neuroimaging studies focus on a frequency-specific or a multi-frequency network analysis showing that functional brain networks are disrupted in patients with Alzheimer's disease (AD). Although those studies enriched our knowledge of the impact of AD in brain's functionality, our goal is to test the effectiveness of combining neuroimaging with network neuroscience to predict with high accuracy subjects with mild cognitive impairment (MCI) that will convert to AD. In this study, eyes-closed resting-state magnetoencephalography (MEG) recordings from 27 stable MCI (sMCI) and 27 progressive MCI (pMCI) from two scan sessions (baseline and follow-up after approximately 3 years) were projected via beamforming onto an atlas-based set of regions of interest (ROIs). Dynamic functional connectivity networks were constructed independently for the five classical frequency bands while a multivariate phase-based coupling metric was adopted. Thus, computing the distance between the fluctuation of functional strength of every pair of ROIs between the two conditions with dynamic time wrapping (DTW), a large set of features was extracted. A machine learning algorithm revealed 30 DTW-based features in the five frequency bands that can distinguish the sMCI from pMCI with absolute accuracy (100%). Further analysis of the selected links revealed that most of the connected ROIs were part of the default mode network (DMN), the cingulo-opercular (CO), the fronto-parietal and the sensorimotor network. Overall, our dynamic network multi-frequency analysis approach provides an effective framework of constructing a sensitive MEG-based connectome biomarker for the prediction of conversion from MCI to Alzheimer's disease.

Bayesian Estimation of Phase Dynamics Based on Partially Sampled Spikes Generated by Realistic Model Neurons

A dynamic system showing stable rhythmic activity can be represented by the dynamics of phase oscillators. This would provide a useful mathematical framework through which one can understand the system's dynamic properties. A recent study proposed a Bayesian approach capable of extracting the underlying phase dynamics directly from time-series data of a system showing rhythmic activity. Here we extended this method to spike data that otherwise provide only limited phase information. To determine how this method performs with spike data, we applied it to simulated spike data generated by a realistic neuronal network model. We then compared the estimated dynamics obtained based on the spike data with the dynamics theoretically derived from the model. The method successfully extracted the modeled phase dynamics, particularly the interaction function, when the amount of available data was sufficiently large. Furthermore, the method was able to infer synaptic connections based on the estimated interaction function. Thus, the method was found to be applicable to spike data and practical for understanding the dynamic properties of rhythmic neural systems.

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.

Statistics of natural binaural sounds

Binaural sound localization is usually considered a discrimination task, where interaural phase (IPD) and level (ILD) disparities at narrowly tuned frequency channels are utilized to identify a position of a sound source. In natural conditions however, binaural circuits are exposed to a stimulation by sound waves originating from multiple, often moving and overlapping sources. Therefore statistics of binaural cues depend on acoustic properties and the spatial configuration of the environment. Distribution of cues encountered naturally and their dependence on physical properties of an auditory scene have not been studied before. In the present work we analyzed statistics of naturally encountered binaural sounds. We performed binaural recordings of three auditory scenes with varying spatial configuration and analyzed empirical cue distributions from each scene. We have found that certain properties such as the spread of IPD distributions as well as an overall shape of ILD distributions do not vary strongly between different auditory scenes. Moreover, we found that ILD distributions vary much weaker across frequency channels and IPDs often attain much higher values, than can be predicted from head filtering properties. In order to understand the complexity of the binaural hearing task in the natural environment, sound waveforms were analyzed by performing Independent Component Analysis (ICA). Properties of learned basis functions indicate that in natural conditions soundwaves in each ear are predominantly generated by independent sources. This implies that the real-world sound localization must rely on mechanisms more complex than a mere cue extraction.

Detecting event-related changes of multivariate phase coupling in dynamic brain networks

Oscillatory phase coupling within large-scale brain networks is a topic of increasing interest within systems, cognitive, and theoretical neuroscience. Evidence shows that brain rhythms play a role in controlling neuronal excitability and response modulation (Haider B, McCormick D. Neuron 62: 171-189, 2009) and regulate the efficacy of communication between cortical regions (Fries P. Trends Cogn Sci 9: 474-480, 2005) and distinct spatiotemporal scales (Canolty RT, Knight RT. Trends Cogn Sci 14: 506-515, 2010). In this view, anatomically connected brain areas form the scaffolding upon which neuronal oscillations rapidly create and dissolve transient functional networks (Lakatos P, Karmos G, Mehta A, Ulbert I, Schroeder C. Science 320: 110-113, 2008). Importantly, testing these hypotheses requires methods designed to accurately reflect dynamic changes in multivariate phase coupling within brain networks. Unfortunately, phase coupling between neurophysiological signals is commonly investigated using suboptimal techniques. Here we describe how a recently developed probabilistic model, phase coupling estimation (PCE; Cadieu C, Koepsell K Neural Comput 44: 3107-3126, 2010), can be used to investigate changes in multivariate phase coupling, and we detail the advantages of this model over the commonly employed phase-locking value (PLV; Lachaux JP, Rodriguez E, Martinerie J, Varela F. Human Brain Map 8: 194-208, 1999). We show that the N-dimensional PCE is a natural generalization of the inherently bivariate PLV. Using simulations, we show that PCE accurately captures both direct and indirect (network mediated) coupling between network elements in situations where PLV produces erroneous results. We present empirical results on recordings from humans and nonhuman primates and show that the PCE-estimated coupling values are different from those using the bivariate PLV. Critically on these empirical recordings, PCE output tends to be sparser than the PLVs, indicating fewer significant interactions and perhaps a more parsimonious description of the data. Finally, the physical interpretation of PCE parameters is straightforward: the PCE parameters correspond to interaction terms in a network of coupled oscillators. Forward modeling of a network of coupled oscillators with parameters estimated by PCE generates synthetic data with statistical characteristics identical to empirical signals. Given these advantages over the PLV, PCE is a useful tool for investigating multivariate phase coupling in distributed brain networks.

Learning intermediate-level representations of form and motion from natural movies

We present a model of intermediate-level visual representation that is based on learning invariances from movies of the natural environment. The model is composed of two stages of processing: an early feature representation layer and a second layer in which invariances are explicitly represented. Invariances are learned as the result of factoring apart the temporally stable and dynamic components embedded in the early feature representation. The structure contained in these components is made explicit in the activities of second-layer units that capture invariances in both form and motion. When trained on natural movies, the first layer produces a factorization, or separation, of image content into a temporally persistent part representing local edge structure and a dynamic part representing local motion structure, consistent with known response properties in early visual cortex (area V1). This factorization linearizes statistical dependencies among the first-layer units, making them learnable by the second layer. The second-layer units are split into two populations according to the factorization in the first layer. The form-selective units receive their input from the temporally persistent part (local edge structure) and after training result in a diverse set of higher-order shape features consisting of extended contours, multiscale edges, textures, and texture boundaries. The motion-selective units receive their input from the dynamic part (local motion structure) and after training result in a representation of image translation over different spatial scales and directions, in addition to more complex deformations. These representations provide a rich description of dynamic natural images and testable hypotheses regarding intermediate-level representation in visual cortex.

Multivariate phase-amplitude cross-frequency coupling in neurophysiological signals

Phase-amplitude cross-frequency coupling (CFC)-where the phase of a low-frequency signal modulates the amplitude or power of a high-frequency signal-is a topic of increasing interest in neuroscience. However, existing methods of assessing CFC are inherently bivariate and cannot estimate CFC between more than two signals at a time. Given the increase in multielectrode recordings, this is a strong limitation. Furthermore, the phase coupling between multiple low-frequency signals is likely to produce a high rate of false positives when CFC is evaluated using bivariate methods. Here, we present a novel method for estimating the statistical dependence between one high-frequency signal and N low-frequency signals, termed multivariate phase-coupling estimation (PCE). Compared to bivariate methods, the PCE produces sparser estimates of CFC and can distinguish between direct and indirect coupling between neurophysiological signals-critical for accurately estimating coupling within multiscale brain networks.

Reconstructing phase dynamics of oscillator networks

We generalize our recent approach to the reconstruction of phase dynamics of coupled oscillators from data [B. Kralemann et al., Phys. Rev. E 77, 066205 (2008)] to cover the case of small networks of coupled periodic units. Starting from a multivariate time series, we first reconstruct genuine phases and then obtain the coupling functions in terms of these phases. Partial norms of these coupling functions quantify directed coupling between oscillators. We illustrate the method by different network motifs for three coupled oscillators and for random networks of five and nine units. We also discuss nonlinear effects in coupling.

Oscillatory phase coupling coordinates anatomically dispersed functional cell assemblies

Hebb proposed that neuronal cell assemblies are critical for effective perception, cognition, and action. However, evidence for brain mechanisms that coordinate multiple coactive assemblies remains lacking. Neuronal oscillations have been suggested as one possible mechanism for cell assembly coordination. Prior studies have shown that spike timing depends upon local field potential (LFP) phase proximal to the cell body, but few studies have examined the dependence of spiking on distal LFP phases in other brain areas far from the neuron or the influence of LFP-LFP phase coupling between distal areas on spiking. We investigated these interactions by recording LFPs and single-unit activity using multiple microelectrode arrays in several brain areas and then used a unique probabilistic multivariate phase distribution to model the dependence of spike timing on the full pattern of proximal LFP phases, distal LFP phases, and LFP-LFP phase coupling between electrodes. Here we show that spiking activity in single neurons and neuronal ensembles depends on dynamic patterns of oscillatory phase coupling between multiple brain areas, in addition to the effects of proximal LFP phase. Neurons that prefer similar patterns of phase coupling exhibit similar changes in spike rates, whereas neurons with different preferences show divergent responses, providing a basic mechanism to bind different neurons together into coordinated cell assemblies. Surprisingly, phase-coupling-based rate correlations are independent of interneuron distance. Phase-coupling preferences correlate with behavior and neural function and remain stable over multiple days. These findings suggest that neuronal oscillations enable selective and dynamic control of distributed functional cell assemblies.