Assessing the coupling between local neural activity and global connectivity fluctuations: Application to human intracranial electroencephalography during a cognitive task.

Cognitive-relevant information is processed by different brain areas that cooperate to eventually produce a response. The relationship between local activity and global brain states during such processes, however, remains for the most part unexplored. To address this question, we designed a simple face-recognition task performed in patients with drug-resistant epilepsy and monitored with intracranial electroencephalography (EEG). Based on our observations, we developed a novel analytical framework (named "local-global" framework) to statistically correlate the brain activity in every recorded gray-matter region with the widespread connectivity fluctuations as proxy to identify concurrent local activations and global brain phenomena that may plausibly reflect a common functional network during cognition. The application of the local-global framework to the data from three subjects showed that similar connectivity fluctuations found across patients were mainly coupled to the local activity of brain areas involved in face information processing. In particular, our findings provide preliminary evidence that the reported global measures might be a novel signature of functional brain activity reorganization when a stimulus is processed in a task context regardless of the specific recorded areas.

Feed-forward information and zero-lag synchronization in the sensory thalamocortical circuit are modulated during stimulus perception.

The direction of functional information flow in the sensory thalamocortical circuit may play a role in stimulus perception, but, surprisingly, this process is poorly understood. We addressed this problem by evaluating a directional information 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 the presence or absence of a tactile stimulus. During stimulus presence, feed-forward information (VPL → S1) increased as a function of the stimulus amplitude, while pure feed-back information (S1 → VPL) was unaffected. In parallel, zero-lag interaction emerged with increasing stimulus amplitude, reflecting externally driven thalamocortical synchronization during stimulus processing. Furthermore, VPL → S1 information decreased during error trials. Also, VPL → S1 and zero-lag interaction decreased when monkeys were not required to report the stimulus presence. These findings provide evidence that both the direction of information flow and the instant synchronization in the sensory thalamocortical circuit play a role in stimulus perception.

Degenerate time-dependent network dynamics anticipate seizures in human epileptic brain.

Epileptic seizures are known to follow specific changes in brain dynamics. While some algorithms can nowadays robustly detect these changes, a clear understanding of the mechanism by which these alterations occur and generate seizures is still lacking. Here, we provide crossvalidated evidence that such changes are initiated by an alteration of physiological network state dynamics. Specifically, our analysis of long intracranial electroencephalography (iEEG) recordings from a group of 10 patients identifies a critical phase of a few hours in which time-dependent network states become less variable ("degenerate"), and this phase is followed by a global functional connectivity reduction before seizure onset. This critical phase is characterized by an abnormal occurrence of highly correlated network instances and is shown to be particularly associated with the activity of the resected regions in patients with validated postsurgical outcome. Our approach characterizes preseizure network dynamics as a cascade of 2 sequential events providing new insights into seizure prediction and control.

Low entropy map of brain oscillatory activity identifies spatially localized events: A new method for automated epilepsy focus prediction.

The spatial mapping of localized events in brain activity critically depends on the correct identification of the pattern signatures associated with those events. For instance, in the context of epilepsy research, a number of different electrophysiological patterns have been associated with epileptogenic activity. Motivated by the need to define automated seizure focus detectors, we propose a novel data-driven algorithm for the spatial identification of localized events that is based on the following rationale: the distribution of emerging oscillations during confined events across all recording sites is highly non-uniform and can be mapped using a spatial entropy function. By applying this principle to EEG recording obtained from 67 distinct seizure epochs, our method successfully identified the seizure focus on a group of ten drug-resistant temporal lobe epilepsy patients (average sensitivity: 0.94, average specificity: 0.90) together with its characteristic electrophysiological pattern signature. Cross-validation of the method outputs with postresective information revealed the consistency of our findings in long follow-up seizure-free patients. Overall, our methodology provides a reliable computational procedure that might be used as in both experimental and clinical domains to identify the neural populations undergoing an emerging functional or pathological transition.

The dynamics of human cognition: Increasing global integration coupled with decreasing segregation found using iEEG.

Cognitive processing requires the ability to flexibly integrate and process information across large brain networks. How do brain networks dynamically reorganize to allow broad communication between many different brain regions in order to integrate information? We record neural activity from 12 epileptic patients using intracranial EEG while performing three cognitive tasks. We assess how the functional connectivity between different brain areas changes to facilitate communication across them. At the topological level, this facilitation is characterized by measures of integration and segregation. Across all patients, we found significant increases in integration and decreases in segregation during cognitive processing, especially in the gamma band (50-90 Hz). We also found higher levels of global synchronization and functional connectivity during task execution, again particularly in the gamma band. More importantly, functional connectivity modulations were not caused by changes in the level of the underlying oscillations. Instead, these modulations were caused by a rearrangement of the mutual synchronization between the different nodes as proposed by the "Communication Through Coherence" Theory.

Task-driven intra- and interarea communications in primate cerebral cortex.

Neural correlations during a cognitive task are central to study brain information processing and computation. However, they have been poorly analyzed due to the difficulty of recording simultaneous single neurons during task performance. In the present work, we quantified neural directional correlations using spike trains that were simultaneously recorded in sensory, premotor, and motor cortical areas of two monkeys during a somatosensory discrimination task. Upon modeling spike trains as binary time series, we used a nonparametric Bayesian method to estimate pairwise directional correlations between many pairs of neurons throughout different stages of the task, namely, perception, working memory, decision making, and motor report. We find that solving the task involves feedforward and feedback correlation paths linking sensory and motor areas during certain task intervals. Specifically, information is communicated by task-driven neural correlations that are significantly delayed across secondary somatosensory cortex, premotor, and motor areas when decision making takes place. Crucially, when sensory comparison is no longer requested for task performance, a major proportion of directional correlations consistently vanish across all cortical areas.

Thalamocortical interactions shape hierarchical neural variability during stimulus perception.

The brain is organized hierarchically to process sensory signals. But, how do functional connections within and across areas contribute to this hierarchical order? We addressed this problem in the thalamocortical network, while monkeys detected vibrotactile stimulus. During this task, we quantified neural variability and directed functional connectivity in simultaneously recorded neurons sharing the cutaneous receptive field within and across VPL and areas 3b and 1. Before stimulus onset, VPL and area 3b exhibited similar fast dynamics while area 1 showed slower timescales. During the stimulus presence, inter-trial neural variability increased along the network VPL-3b-1 while VPL established two main feedforward pathways with areas 3b and 1 to process the stimulus. This lower variability of VPL and area 3b was found to regulate feedforward thalamocortical pathways. Instead, intra-cortical interactions were only anticipated by higher intrinsic timescales in area 1. Overall, our results provide evidence of hierarchical functional roles along the thalamocortical network.

A mixture of mobility and meteorological data provides a high correlation with COVID-19 growth in an infection-naive population: a study for Spanish provinces.

Introduction: We use Spanish data from August 2020 to March 2021 as a natural experiment to analyze how a standardized measure of COVID-19 growth correlates with asymmetric meteorological and mobility situations in 48 Spanish provinces. The period of time is selected prior to vaccination so that the level of susceptibility was high, and during geographically asymmetric implementation of non-pharmacological interventions. Methods: We develop reliable aggregated mobility data from different public sources and also compute the average meteorological time series of temperature, dew point, and UV radiance in each Spanish province from satellite data. We perform a dimensionality reduction of the data using principal component analysis and investigate univariate and multivariate correlations of mobility and meteorological data with COVID-19 growth. Results: We find significant, but generally weak, univariate correlations for weekday aggregated mobility in some, but not all, provinces. On the other hand, principal component analysis shows that the different mobility time series can be properly reduced to three time series. A multivariate time-lagged canonical correlation analysis of the COVID-19 growth rate with these three time series reveals a highly significant correlation, with a median R-squared of 0.65. The univariate correlation between meteorological data and COVID-19 growth is generally not significant, but adding its two main principal components to the mobility multivariate analysis increases correlations significantly, reaching correlation coefficients between 0.6 and 0.98 in all provinces with a median R-squared of 0.85. This result is robust to different approaches in the reduction of dimensionality of the data series. Discussion: Our results suggest an important effect of mobility on COVID-19 cases growth rate. This effect is generally not observed for meteorological variables, although in some Spanish provinces it can become relevant. The correlation between mobility and growth rate is maximal at a time delay of 2-3 weeks, which agrees well with the expected 5?10 day delays between infection, development of symptoms, and the detection/report of the case.

A comparative study between a power and a connectivity sEEG biomarker for seizure-onset zone identification in temporal lobe epilepsy.

BACKGROUND: Ictal stereo-encephalography (sEEG) biomarkers for seizure onset zone (SOZ) localization can be classified depending on whether they target abnormalities in signal power or functional connectivity between signals, and they may depend on the frequency band and time window at which they are estimated. NEW METHOD: This work aimed to compare and optimize the performance of a power and a connectivity-based biomarker to identify SOZ contacts from ictal sEEG recordings. To do so, we used a previously introduced power-based measure, the normalized mean activation (nMA), which quantifies the ictal average power activation. Similarly, we defined the normalized mean strength (nMS), to quantify the ictal mean functional connectivity of every contact with the rest. The optimal frequency bands and time windows were selected based on optimizing AUC and F2-score. RESULTS: The analysis was performed on a dataset of 67 seizures from 10 patients with pharmacoresistant temporal lobe epilepsy. Our results suggest that the power-based biomarker generally performs better for the detection of SOZ than the connectivity-based one. However, an equivalent performance level can be achieved when both biomarkers are independently optimized. Optimal performance was achieved in the beta and lower-gamma range for the power biomarker and in the lower- and higher-gamma range for connectivity, both using a 20 or 30 s period after seizure onset. CONCLUSIONS: The results of this study highlight the importance of this optimization step over frequency and time windows when comparing different SOZ discrimination biomarkers. This information should be considered when training SOZ classifiers on retrospective patients' data for clinical applications.

Inferring neural information flow from spiking data.

The brain can be regarded as an information processing system in which neurons store and propagate information about external stimuli and internal processes. Therefore, estimating interactions between neural activity at the cellular scale has significant implications in understanding how neuronal circuits encode and communicate information across brain areas to generate behavior. While the number of simultaneously recorded neurons is growing exponentially, current methods relying only on pairwise statistical dependencies still suffer from a number of conceptual and technical challenges that preclude experimental breakthroughs describing neural information flows. In this review, we examine the evolution of the field over the years, starting from descriptive statistics to model-based and model-free approaches. Then, we discuss in detail the Granger Causality framework, which includes many popular state-of-the-art methods and we highlight some of its limitations from a conceptual and practical estimation perspective. Finally, we discuss directions for future research, including the development of theoretical information flow models and the use of dimensionality reduction techniques to extract relevant interactions from large-scale recording datasets.

Human hippocampal theta oscillations reflect sequential dependencies during spatial planning.

Movement-related theta oscillations in rodent hippocampus coordinate 'forward sweeps' of location-specific neural activity that could be used to evaluate spatial trajectories online. This raises the possibility that increases in human hippocampal theta power accompany the evaluation of upcoming spatial choices. To test this hypothesis, we measured neural oscillations during a spatial planning task that closely resembles a perceptual decision-making paradigm. In this task, participants searched visually for the shortest path between a start and goal location in novel mazes that contained multiple choice points, and were subsequently asked to make a spatial decision at one of those choice points. We observed ~4-8 Hz hippocampal/medial temporal lobe theta power increases specific to sequential planning that were negatively correlated with subsequent decision speed, where decision speed was inversely correlated with choice accuracy. These results implicate the hippocampal theta rhythm in decision tree search during planning in novel environments.

Serotonin 5-HT1A, 5-HT2A and dopamine D2 receptors strongly influence prefronto-hippocampal neural networks in alert mice: Contribution to the actions of risperidone.

Atypical antipsychotic drugs (APDs) used to treat positive and negative symptoms in schizophrenia block serotonin receptors 5-HT2AR and dopamine receptors D2R and stimulate 5-HT1AR directly or indirectly. However, the exact cellular mechanisms mediating their therapeutic actions remain unresolved. We recorded neural activity in the prefrontal cortex (PFC) and hippocampus (HPC) of freely-moving mice before and after acute administration of 5-HT1AR, 5-HT2AR and D2R selective agonists and antagonists and atypical APD risperidone. We then investigated the contribution of the three receptors to the actions of risperidone on brain activity via statistical modeling and pharmacological reversal (risperidone + 5-HT1AR antagonist WAY-100635, risperidone + 5-HT2A/2CR agonist DOI, risperidone + D2R agonist quinpirole). Risperidone, 5-HT1AR agonism with 8-OH-DPAT, 5-HT2AR antagonism with M100907, and D2R antagonism with haloperidol reduced locomotor activity of mice that correlated with a suppression of neural spiking, power of theta and gamma oscillations in PFC and HPC, and reduction of PFC-HPC theta phase synchronization. By contrast, activation of 5-HT2AR with DOI enhanced high-gamma oscillations in PFC and PFC-HPC high gamma functional connectivity, likely related to its hallucinogenic effects. Together, power changes, regression modeling and pharmacological reversals suggest an important role of 5-HT1AR agonism and 5-HT2AR antagonism in risperidone-induced alterations of delta, beta and gamma oscillations, while D2R antagonism may contribute to risperidone-mediated changes in delta oscillations. This study provides novel insight into the neural mechanisms for widely prescribed psychiatric medication targeting the serotonin and dopamine systems in two regions involved in the pathophysiology of schizophrenia.

Detection of recurrent activation patterns across focal seizures: Application to seizure onset zone identification.

OBJECTIVE: We introduce a method that quantifies the consistent involvement of intracranially monitored regions in recurrent focal seizures. METHODS: We evaluated the consistency of two ictal spectral activation patterns (mean power change and power change onset time) in intracranial recordings across focal seizures from seven patients with clinically marked seizure onset zone (SOZ). We examined SOZ discrimination using both patterns in different frequency bands and periods of interest. RESULTS: Activation patterns were proved to be consistent across more than 80% of recurrent ictal epochs. In all patients, whole-seizure mean activations were significantly higher for SOZ than non-SOZ regions (P<0.05) while activation onset times were significantly lower for SOZ than for non-SOZ regions (P<0.001) in six patients. Alpha-beta bands (8-20Hz) achieved the highest patient-average effect size on the whole-seizure period while gamma band (20-70Hz) achieved the highest discrimination values between SOZ and non-SOZ sites near seizure onset (0-5s). CONCLUSIONS: Consistent spectral activation patterns in focal epilepsies discriminate the SOZ with high effect sizes upon appropriate selection of frequency bands and activation periods. SIGNIFICANCE: The present method may be used to improve epileptogenic identification as well as pinpoint additional regions that are functionally altered during ictal events.