Neuroimaging techniques are now widely used to study human cognition. The functional associations between brain areas have become a standard proxy to describe how cognitive processes are distributed across the brain network. Among the many analysis tools available, dynamic models of brain activity have been developed to overcome the limitations of original connectivity measures such as functional connectivity. This goes in line with the many efforts devoted to the assessment of directional interactions between brain areas from the observed neuroimaging activity. This opinion article provides an overview of our model-based whole-brain effective connectivity to analyze fMRI data, while discussing the pros and cons of our approach with respect to other established approaches. Our framework relies on the multivariate Ornstein-Uhlenbeck (MOU) process and is thus referred to as MOU-EC. Once tuned, the model provides a directed connectivity estimate that reflects the dynamical state of BOLD activity, which can be used to explore cognition. We illustrate this approach using two applications on task-evoked fMRI data. First, as a connectivity measure, MOU-EC can be used to extract biomarkers for task-specific brain coordination, understood as the patterns of areas exchanging information. The multivariate nature of connectivity measures raises several challenges for whole-brain analysis, for which machine-learning tools present some advantages over statistical testing. Second, we show how to interpret changes in MOU-EC connections in a collective and model-based manner, bridging with network analysis. Our framework provides a comprehensive set of tools that open exciting perspectives to study distributed cognition, as well as neuropathologies.

Significance: How do multiple neurons communicate to solve a cognitive task? To answer this question, we investigate spike-train directional correlations across five primate cortical areas simultaneously recorded during a somatosensory discrimination task. Correlations are inferred using a nonparametric procedure that models spike trains as Markovian binary series and dynamically estimates the directed information between every neuron pair at different delays. We find that information processing during the discrimination task can be described by intra- and interarea decision-driven delayed correlations, which are no longer found when a monkey receives both stimuli but does not perform the task.

Significance: The direction of information flow between brain circuits may be key in cognitive functions. We addressed this problem by evaluating a directional correlation measure between simultaneously recorded neurons from somatosensory thalamus (ventral posterolateral nucleus, VPL) and somatosensory cortex (S1) sharing the same cutaneous receptive field while monkeys judged presence or absence of tactile stimuli. During stimulus presence, feed-forward (VPL → S1) information increased as a function of stimulus amplitude, while feed-back (S1 → VPL) information was unaffected. Simultaneously, zero-lag interaction emerged with increasing stimulus amplitude, contributing to thalamocortical synchronization. Furthermore, VPL → S1 information decreased during error trials. Also, both VPL → S1 and zero-lag interactions decreased when monkeys were not required to report stimulus presence. Thus, directional and coordinated information in the thalamocortical circuit is associated with stimulus perception.

Two alternative exact characterizations of the minimum error probability of Bayesian M-ary hypothesis testing are derived. The first expression corresponds to the error probability of an induced binary hypothesis test and implies the tightness of the meta-converse bound by Polyanskiy et al.; the second expression is a function of an information-spectrum measure and implies the tightness of a generalized Verdú-Han lower bound. The formulas characterize the minimum error probability of several problems in information theory and help to identify the steps where existing converse bounds are loose.

Directed connectivity inference has become a cornerstone in neuroscience to analyze multivariate data from neuroimaging and electrophysiological techniques. Here we propose a nonparametric significance method to test the nonzero values of multivariate autoregressive model to infer interactions in recurrent networks. We use random permutations or circular shifts of the original time series to generate the null-hypothesis distributions. The underlying network model is the same as used in multivariate Granger causality, but our test relies on the autoregressive coefficients instead of error residuals. By means of numerical simulation over multiple network configurations, we show that this method achieves a good control of false positives (type 1 error) and detects existing pairwise connections more accurately than using the standard parametric test for the ratio of error residuals. In practice, our method aims to detect temporal interactions in real neuronal networks with nodes possibly exhibiting redundant activity. As a proof of concept, we apply our method to multiunit activity (MUA) recorded from Utah electrode arrays in a monkey and examine detected interactions between 25 channels. We show that during stimulus presentation our method detects a large number of interactions that cannot be solely explained by the increase in the MUA level.