Granger causality analysis of fMRI BOLD signals is invariant to hemodynamic convolution but not downsampling

EEG-based functional and effective connectivity patterns during emotional episodes using graph theoretical analysis

Understanding the neural mechanisms underlying emotional processing is critical for advancing neuroscience and mental health interventions. This study examined these mechanisms by analyzing EEG connectivity patterns across different brain regions while participants evoked various emotions. After applying independent component analysis (ICA) to eliminate non-cortical activity, we assessed frequency-specific connectivity patterns using coherence, Granger causality, and graph theoretical measures to evaluate both functional and effective connectivity. Graph theoretical analysis revealed significant differences in connectivity between emotions across multiple frequency bands, particularly in the delta and beta bands. These results indicated modulations in key brain regions, such as the precentral, superior frontal, and temporal areas, suggesting that these regions play a crucial role in emotional processing. Coherence analysis demonstrated predominant alpha band activity across all emotions, with specific emotional states, such as fear, grief, and jealousy, exhibiting enhanced beta band activity. In terms of coherence strength, we observed that the gamma band was largely inactive, except for the emotion of sadness, which displayed significant activity in the right lobe, particularly in regions such as the supplementary motor area and the lingual gyrus. Additionally, Granger causality analysis highlighted that the beta and gamma bands were dominant across all emotional states, with minimal modulation observed in the theta band. Clustering coefficients from the graph analysis further revealed distinct patterns of connectivity in the delta and beta bands, with significant variations across different emotions, particularly in the temporal and frontal regions. These findings enhance our understanding of emotional processing and have practical applications in mental health, biomarker identification, and human-computer interaction.

Sex-dependent nonlinear Granger connectivity patterns of brain aging in healthy population

BACKGROUND

Brain aging is a complex process that involves functional alterations in multiple subnetworks and brain regions. However, most previous studies investigating aging-related functional connectivity (FC) changes using resting-state functional magnetic resonance images (rs-fMRIs) have primarily focused on the linear correlation between brain subnetworks, ignoring the nonlinear casual properties of fMRI signals.

METHODS

We introduced the neural Granger causality technique to investigate the sex-dependent nonlinear Granger connectivity (NGC) during aging on a publicly available dataset of 227 healthy participants acquired cross-sectionally in Leipzig, Germany.

RESULTS

Our findings indicate that brain aging may cause widespread declines in NGC at both regional and subnetwork scales. These findings exhibit high reproducibility across different network sparsities, demonstrating the efficacy of static and dynamic analysis strategies. Females exhibit greater heterogeneity and reduced stability in NGC compared to males during aging, especially the NGC between the visual network and other subnetworks. Besides, NGC strengths can well reflect the individual cognitive function, which may therefore work as a sensitive metric in cognition-related experiments for individual-scale or group-scale mechanism understanding.

CONCLUSION

These findings indicate that NGC analysis is a potent tool for identifying sex-dependent brain aging patterns. Our results offer valuable perspectives that could substantially enhance the understanding of sex differences in neurological diseases in the future, especially in degenerative disorders.

Fractal Conditional Correlation Dimension Infers Complex Causal Networks

Determining causal inference has become popular in physical and engineering applications. While the problem has immense challenges, it provides a way to model the complex networks by observing the time series. In this paper, we present the optimal conditional correlation dimensional geometric information flow principle (oGeoC) that can reveal direct and indirect causal relations in a network through geometric interpretations. We introduce two algorithms that utilize the oGeoC principle to discover the direct links and then remove indirect links. The algorithms are evaluated using coupled logistic networks. The results indicate that when the number of observations is sufficient, the proposed algorithms are highly accurate in identifying direct causal links and have a low false positive rate.

The subcortical brain regions influence the cortical areas during resting-state: an fMRI study

Introduction

Numerous modes or patterns of neural activity can be seen in the brain of individuals during the resting state. However, those functions do not persist long, and they are continuously altering in the brain. We have hypothesized that the brain activations during the resting state should themselves be responsible for this alteration of the activities.

Methods

Using the resting-state fMRI data of 63 healthy young individuals, we estimated the causality effects of each resting-state activation map on all other networks. The resting-state networks were identified, their causality effects on the other components were extracted, the networks with the top 20% of the causality were chosen, and the networks which were under the influence of those causal networks were also identified.

Results

Our results showed that the influence of each activation component over other components is different. The brain areas which showed the highest causality coefficients were subcortical regions, such as the brain stem, thalamus, and amygdala. On the other hand, nearly all the areas which were mostly under the causal effects were cortical regions.

Discussion

In summary, our results suggest that subcortical brain areas exert a higher influence on cortical regions during the resting state, which could help in a better understanding the dynamic nature of brain functions.

Causalized Convergent Cross Mapping and Its Implementation in Causality Analysis

Rooted in dynamic systems theory, convergent cross mapping (CCM) has attracted increased attention recently due to its capability in detecting linear and nonlinear causal coupling in both random and deterministic settings. One limitation with CCM is that it uses both past and future values to predict the current value, which is inconsistent with the widely accepted definition of causality, where it is assumed that the future values of one process cannot influence the past of another. To overcome this obstacle, in our previous research, we introduced the concept of causalized convergent cross mapping (cCCM), where future values are no longer used to predict the current value. In this paper, we focus on the implementation of cCCM in causality analysis. More specifically, we demonstrate the effectiveness of cCCM in identifying both linear and nonlinear causal coupling in various settings through a large number of examples, including Gaussian random variables with additive noise, sinusoidal waveforms, autoregressive models, stochastic processes with a dominant spectral component embedded in noise, deterministic chaotic maps, and systems with memory, as well as experimental fMRI data. In particular, we analyze the impact of shadow manifold construction on the performance of cCCM and provide detailed guidelines on how to configure the key parameters of cCCM in different applications. Overall, our analysis indicates that cCCM is a promising and easy-to-implement tool for causality analysis in a wide spectrum of applications.

Spectral dynamic causal modeling: A didactic introduction and its relationship with functional connectivity

We present a didactic introduction to spectral dynamic causal modeling (DCM), a Bayesian state-space modeling approach used to infer effective connectivity from noninvasive neuroimaging data. Spectral DCM is currently the most widely applied DCM variant for resting-state functional MRI analysis. Our aim is to explain its technical foundations to an audience with limited expertise in state-space modeling and spectral data analysis. Particular attention will be paid to cross-spectral density, which is the most distinctive feature of spectral DCM and is closely related to functional connectivity, as measured by (zero-lag) Pearson correlations. In fact, the model parameters estimated by spectral DCM are those that best reproduce the cross-correlations between all measurements-at all time lags-including the zero-lag correlations that are usually interpreted as functional connectivity. We derive the functional connectivity matrix from the model equations and show how changing a single effective connectivity parameter can affect all pairwise correlations. To complicate matters, the pairs of brain regions showing the largest changes in functional connectivity do not necessarily coincide with those presenting the largest changes in effective connectivity. We discuss the implications and conclude with a comprehensive summary of the assumptions and limitations of spectral DCM.

Exploring brain effective connectivity of early MCI with GRU_GC model on resting-state fMRI

BACKGROUND

Investigating the functional interactions between different brain regions and revealing the transmission of information by computing brain connectivity have great potential and significance in the diagnosis of early Mild Cognitive Impairment (EMCI).

METHODS

The Granger causality with Gate Recurrent Unit (GRU_GC) model is a suitable method that allows the detection of a nonlinear causal relationship and solves the limitation of fixed time lag, which cannot be detected by the classical Granger method. The model can transmit time series signals with any transmission delay length, and the time series can be screened and learned through the gate model.

RESULTS

The classification experiment of 89 EMCI and 73 neurologically healthy controls (HC) shows that the accuracy reached 87.88%. Compared with multivariate variables GC (MVGC) and Long Short-Term Memory-based GC (LSTM_GC), the GRU_GC significantly improved the estimation of brain connectivity communication. Constructing a difference network to explore the brain effective connectivity between EMCI and HC.

CONCLUSIONS

The GRU_GC can discover the abnormal brain regions, including the parahippocampal gyrus, the posterior cingulate gyrus. The method can be used in clinical applications as an effective brain connectivity analysis tool and provides auxiliary means for the medical diagnosis of EMCI.

Perineuronal nets in the rat medial prefrontal cortex alter hippocampal-prefrontal oscillations and reshape cocaine self-administration memories

The medial prefrontal cortex (mPFC) is a major contributor to relapse to cocaine in humans and to reinstatement behavior in rodent models of cocaine use disorder. Output from the mPFC is modulated by parvalbumin (PV)-containing fast-spiking interneurons, the majority of which are surrounded by perineuronal nets (PNNs). Here we tested whether chondroitinase ABC (ABC)- mediated removal of PNNs prevented the acquisition or reconsolidation of a cocaine self-administration memory. ABC injections into the dorsal mPFC prior to training attenuated the acquisition of cocaine self-administration. Also, ABC given 3 days prior to but not 1 hr after memory reactivation blocked cue-induced reinstatement. However, reduced reinstatement was present only in rats given a novel reactivation contingency, suggesting that PNNs are required for the updating of a familiar memory. In naive rats, ABC injections into mPFC did not alter excitatory or inhibitory puncta on PV cells but reduced PV intensity. Whole-cell recordings revealed a greater inter-spike interval 1 hr after ABC, but not 3 days later. In vivo recordings from the mPFC and dorsal hippocampus (dHIP) during novel memory reactivation revealed that ABC in the mPFC prevented reward-associated increases in beta and gamma activity as well as phase-amplitude coupling between the dHIP and mPFC. Together, our findings show that PNN removal attenuates the acquisition of cocaine self-administration memories and disrupts reconsolidation of the original memory when combined with a novel reactivation session. Further, reduced dHIP/mPFC coupling after PNN removal may serve as a key biomarker for how to disrupt reconsolidation of cocaine memories and reduce relapse.

Causalized convergent cross-mapping and its approximate equivalence with directed information in causality analysis

Convergent cross-mapping (CCM) has attracted increased attention recently due to its capability to detect causality in nonseparable systems under deterministic settings, which may not be covered by the traditional Granger causality. From an information-theoretic perspective, causality is often characterized as the directed information (DI) flowing from one side to the other. As information is essentially nondeterministic, a natural question is: does CCM measure DI flow? Here, we first causalize CCM so that it aligns with the presumption in causality analysis-the future values of one process cannot influence the past of the other, and then establish and validate the approximate equivalence of causalized CCM (cCCM) and DI under Gaussian variables through both theoretical derivations and fMRI-based brain network causality analysis. Our simulation result indicates that, in general, cCCM tends to be more robust than DI in causality detection. The underlying argument is that DI relies heavily on probability estimation, which is sensitive to data size as well as digitization procedures; cCCM, on the other hand, gets around this problem through geometric cross-mapping between the manifolds involved. Overall, our analysis demonstrates that cross-mapping provides an alternative way to evaluate DI and is potentially an effective technique for identifying both linear and nonlinear causal coupling in brain neural networks and other settings, either random or deterministic, or both.

Static and dynamic fMRI-derived functional connectomes represent largely similar information

Functional connectivity (FC) of blood oxygen level-dependent (BOLD) fMRI time series can be estimated using methods that differ in sensitivity to the temporal order of time points (static vs. dynamic) and the number of regions considered in estimating a single edge (bivariate vs. multivariate). Previous research suggests that dynamic FC explains variability in FC fluctuations and behavior beyond static FC. Our aim was to systematically compare methods on both dimensions. We compared five FC methods: Pearson's/full correlation (static, bivariate), lagged correlation (dynamic, bivariate), partial correlation (static, multivariate), and multivariate AR model with and without self-connections (dynamic, multivariate). We compared these methods by (i) assessing similarities between FC matrices, (ii) by comparing node centrality measures, and (iii) by comparing the patterns of brain-behavior associations. Although FC estimates did not differ as a function of sensitivity to temporal order, we observed differences between the multivariate and bivariate FC methods. The dynamic FC estimates were highly correlated with the static FC estimates, especially when comparing group-level FC matrices. Similarly, there were high correlations between the patterns of brain-behavior associations obtained using the dynamic and static FC methods. We conclude that the dynamic FC estimates represent information largely similar to that of the static FC.

Causal influence of brainstem response to transcutaneous vagus nerve stimulation on cardiovagal outflow

BACKGROUND

The autonomic response to transcutaneous auricular vagus nerve stimulation (taVNS) has been linked to the engagement of brainstem circuitry modulating autonomic outflow. However, the physiological mechanisms supporting such efferent vagal responses are not well understood, particularly in humans.

HYPOTHESIS

We present a paradigm for estimating directional brain-heart interactions in response to taVNS. We propose that our approach is able to identify causal links between the activity of brainstem nuclei involved in autonomic control and cardiovagal outflow.

METHODS

We adopt an approach based on a recent reformulation of Granger causality that includes permutation-based, nonparametric statistics. The method is applied to ultrahigh field (7T) functional magnetic resonance imaging (fMRI) data collected on healthy subjects during taVNS.

RESULTS

Our framework identified taVNS-evoked functional brainstem responses with superior sensitivity compared to prior conventional approaches, confirming causal links between taVNS stimulation and fMRI response in the nucleus tractus solitarii (NTS). Furthermore, our causal approach elucidated potential mechanisms by which information is relayed between brainstem nuclei and cardiovagal, i.e., high-frequency heart rate variability, in response to taVNS. Our findings revealed that key brainstem nuclei, known from animal models to be involved in cardiovascular control, exert a causal influence on taVNS-induced cardiovagal outflow in humans.

CONCLUSION

Our causal approach allowed us to noninvasively evaluate directional interactions between fMRI BOLD signals from brainstem nuclei and cardiovagal outflow.

Quantification of inter-brain coupling: A review of current methods used in haemodynamic and electrophysiological hyperscanning studies

Hyperscanning is a form of neuroimaging experiment where the brains of two or more participants are imaged simultaneously whilst they interact. Within the domain of social neuroscience, hyperscanning is increasingly used to measure inter-brain coupling (IBC) and explore how brain responses change in tandem during social interaction. In addition to cognitive research, some have suggested that quantification of the interplay between interacting participants can be used as a biomarker for a variety of cognitive mechanisms aswell as to investigate mental health and developmental conditions including schizophrenia, social anxiety and autism. However, many different methods have been used to quantify brain coupling and this can lead to questions about comparability across studies and reduce research reproducibility. Here, we review methods for quantifying IBC, and suggest some ways moving forward. Following the PRISMA guidelines, we reviewed 215 hyperscanning studies, across four different brain imaging modalities: functional near-infrared spectroscopy (fNIRS), functional magnetic resonance (fMRI), electroencephalography (EEG) and magnetoencephalography (MEG). Overall, the review identified a total of 27 different methods used to compute IBC. The most common hyperscanning modality is fNIRS, used by 119 studies, 89 of which adopted wavelet coherence. Based on the results of this literature survey, we first report summary statistics of the hyperscanning field, followed by a brief overview of each signal that is obtained from each neuroimaging modality used in hyperscanning. We then discuss the rationale, assumptions and suitability of each method to different modalities which can be used to investigate IBC. Finally, we discuss issues surrounding the interpretation of each method.

The arrow of time of brain signals in cognition: Potential intriguing role of parts of the default mode network

A promising idea in human cognitive neuroscience is that the default mode network (DMN) is responsible for coordinating the recruitment and scheduling of networks for computing and solving task-specific cognitive problems. This is supported by evidence showing that the physical and functional distance of DMN regions is maximally removed from sensorimotor regions containing environment-driven neural activity directly linked to perception and action, which would allow the DMN to orchestrate complex cognition from the top of the hierarchy. However, discovering the functional hierarchy of brain dynamics requires finding the best way to measure interactions between brain regions. In contrast to previous methods measuring the hierarchical flow of information using, for example, transfer entropy, here we used a thermodynamics-inspired, deep learning based Temporal Evolution NETwork (TENET) framework to assess the asymmetry in the flow of events, 'arrow of time', in human brain signals. This provides an alternative way of quantifying hierarchy, given that the arrow of time measures the directionality of information flow that leads to a breaking of the balance of the underlying hierarchy. In turn, the arrow of time is a measure of nonreversibility and thus nonequilibrium in brain dynamics. When applied to large-scale Human Connectome Project (HCP) neuroimaging data from close to a thousand participants, the TENET framework suggests that the DMN plays a significant role in orchestrating the hierarchy, that is, levels of nonreversibility, which changes between the resting state and when performing seven different cognitive tasks. Furthermore, this quantification of the hierarchy of the resting state is significantly different in health compared to neuropsychiatric disorders. Overall, the present thermodynamics-based machine-learning framework provides vital new insights into the fundamental tenets of brain dynamics for orchestrating the interactions between cognition and brain in complex environments.

Causal inference in cognitive neuroscience

Causal inference is a key step in many research endeavors in cognitive science and neuroscience, and particularly cognitive neuroscience. Statistical knowledge is sufficient for prediction and diagnosis, but causal knowledge is required for action and intervention. Most statistics courses and textbooks emphasize the difficulty of causal inference, focusing on the maxim that "correlation does not mean causation": there can be multiple causal possibilities, often many of them, consistent with given observed statistics. This paper focuses instead on the conceptual issues and assumptions that confront causal and other kinds of inference, primarily focusing on cognitive neuroscience. We connect inference methods with goals and challenges, and provide concrete guidance about how to select appropriate tools for the scientific task. This article is categorized under: Psychology > Theory and Methods Philosophy > Foundations of Cognitive Science.

Dynamical independence: Discovering emergent macroscopic processes in complex dynamical systems

We introduce a notion of emergence for macroscopic variables associated with highly multivariate microscopic dynamical processes. Dynamical independence instantiates the intuition of an emergent macroscopic process as one possessing the characteristics of a dynamical system "in its own right," with its own dynamical laws distinct from those of the underlying microscopic dynamics. We quantify (departure from) dynamical independence by a transformation-invariant Shannon information-based measure of dynamical dependence. We emphasize the data-driven discovery of dynamically independent macroscopic variables, and introduce the idea of a multiscale "emergence portrait" for complex systems. We show how dynamical dependence may be computed explicitly for linear systems in both time and frequency domains, facilitating discovery of emergent phenomena across spatiotemporal scales, and outline application of the linear operationalization to inference of emergence portraits for neural systems from neurophysiological time-series data. We discuss dynamical independence for discrete- and continuous-time deterministic dynamics, with potential application to Hamiltonian mechanics and classical complex systems such as flocking and cellular automata.

Static and dynamic functional connectomes represent largely similar information

Functional connectivity (FC) of blood-oxygen-level-dependent (BOLD) fMRI time series can be estimated using methods that differ in sensitivity to the temporal order of time points (static vs. dynamic) and the number of regions considered in estimating a single edge (bivariate vs. multivariate). Previous research suggests that dynamic FC explains variability in FC fluctuations and behavior beyond static FC. Our aim was to systematically compare methods on both dimensions. We compared five FC methods: Pearson's/full correlation (static, bivariate), lagged correlation (dynamic, bivariate), partial correlation (static, multivariate) and multivariate AR model with and without self-connections (dynamic, multivariate). We compared these methods by (i) assessing similarities between FC matrices, (ii) by comparing node centrality measures, and (iii) by comparing the patterns of brain-behavior associations. Although FC estimates did not differ as a function of sensitivity to temporal order, we observed differences between the multivariate and bivariate FC methods. The dynamic FC estimates were highly correlated with the static FC estimates, especially when comparing group-level FC matrices. Similarly, there were high correlations between the patterns of brain-behavior associations obtained using the dynamic and static FC methods. We conclude that the dynamic FC estimates represent information largely similar to that of the static FC.

Brain Connectivity Analysis in Distinct Footwear Conditions during Infinity Walk Using fNIRS

Gait and balance are an intricate interplay between the brain, nervous system, sensory organs, and musculoskeletal system. They are greatly influenced by the type of footwear, walking patterns, and surface. This exploratory study examines the effects of the Infinity Walk, pronation, and footwear conditions on brain effective connectivity patterns. A continuous-wave functional near-infrared spectroscopy device collected data from five healthy participants. A highly computationally efficient connectivity model based on the Grange causal relationship between the channels was applied to data to find the effective relationship between inter- and intra-hemispheric brain connectivity. Brain regions of interest (ROI) were less connected during the barefoot condition than during other complex walks. Conversely, the highest interconnectedness between ROI was observed while wearing flat insoles and medially wedged sandals, which is a relatively difficult type of footwear to walk in. No statistically significant (p-value <0.05) effect on connectivity patterns was observed during the corrected pronated posture. The regions designated as motoric, sensorimotor, and temporal became increasingly connected with difficult walking patterns and footwear conditions. The Infinity Walk causes effective bidirectional connections between ROI across all conditions and both hemispheres. Due to its repetitive pattern, the Infinity Walk is a good test method, particularly for neuro-rehabilitation and motoric learning experiments.

Heritability of human "directed" functional connectome

INTRODUCTION

The functional connectivity patterns in the brain are highly heritable; however, it is unclear how genetic factors influence the directionality of such "information flows." Studying the "directionality" of the brain functional connectivity and assessing how heritability modulates it can improve our understanding of the human connectome.

METHODS

Here, we investigated the heritability of "directed" functional connections using a state-space formulation of Granger causality (GC), in conjunction with blind deconvolution methods accounting for local variability in the hemodynamic response function. Such GC implementation is ideal to explore the directionality of functional interactions across a large number of networks. Resting-state functional magnetic resonance imaging data were drawn from the Human Connectome Project (total n = 898 participants). To add robustness to our findings, the dataset was randomly split into a "discovery" and a "replication" sample (each with n = 449 participants). The two cohorts were carefully matched in terms of demographic variables and other confounding factors (e.g., education). The effect of shared environment was also modeled.

RESULTS

The parieto- and prefronto-cerebellar, parieto-prefrontal, and posterior-cingulate to hippocampus connections showed the highest and most replicable heritability effects with little influence by shared environment. In contrast, shared environmental factors significantly affected the visuo-parietal and sensory-motor directed connectivity.

CONCLUSION

We suggest a robust role of heritability in influencing the directed connectivity of some cortico-subcortical circuits implicated in cognition. Further studies, for example using task-based fMRI and GC, are warranted to confirm the asymmetric effects of genetic factors on the functional connectivity within cognitive networks and their role in supporting executive functions and learning.

Brain effective connectivity and functional connectivity as markers of lifespan vascular exposures in middle-aged adults: The Bogalusa Heart Study

Introduction

Effective connectivity (EC), the causal influence that functional activity in a source brain location exerts over functional activity in a target brain location, has the potential to provide different information about brain network dynamics than functional connectivity (FC), which quantifies activity synchrony between locations. However, head-to-head comparisons between EC and FC from either task-based or resting-state functional MRI (fMRI) data are rare, especially in terms of how they associate with salient aspects of brain health.

Methods

In this study, 100 cognitively-healthy participants in the Bogalusa Heart Study aged 54.2 ± 4.3years completed Stroop task-based fMRI, resting-state fMRI. EC and FC among 24 regions of interest (ROIs) previously identified as involved in Stroop task execution (EC-task and FC-task) and among 33 default mode network ROIs (EC-rest and FC-rest) were calculated from task-based and resting-state fMRI using deep stacking networks and Pearson correlation. The EC and FC measures were thresholded to generate directed and undirected graphs, from which standard graph metrics were calculated. Linear regression models related graph metrics to demographic, cardiometabolic risk factors, and cognitive function measures.

Results

Women and whites (compared to men and African Americans) had better EC-task metrics, and better EC-task metrics associated with lower blood pressure, white matter hyperintensity volume, and higher vocabulary score (maximum value of p = 0.043). Women had better FC-task metrics, and better FC-task metrics associated with APOE-ε4 3–3 genotype and better hemoglobin-A1c, white matter hyperintensity volume and digit span backwards score (maximum value of p = 0.047). Better EC rest metrics associated with lower age, non-drinker status, and better BMI, white matter hyperintensity volume, logical memory II total score, and word reading score (maximum value of p = 0.044). Women and non-drinkers had better FC-rest metrics (value of p = 0.004).

Discussion

In a diverse, cognitively healthy, middle-aged community sample, EC and FC based graph metrics from task-based fMRI data, and EC based graph metrics from resting-state fMRI data, were associated with recognized indicators of brain health in differing ways. Future studies of brain health should consider taking both task-based and resting-state fMRI scans and measuring both EC and FC analyses to get a more complete picture of functional networks relevant to brain health.

Toward naturalistic neuroscience: Mechanisms underlying the flattening of brain hierarchy in movie-watching compared to rest and task

Identifying the functional specialization of the brain has moved from using cognitive tasks and resting state to using ecological relevant, naturalistic movies. We leveraged a large-scale neuroimaging dataset to directly investigate the hierarchical reorganization of functional brain activity when watching naturalistic films compared to performing seven cognitive tasks and resting. A thermodynamics-inspired whole-brain model paradigm revealed the generative underlying mechanisms for changing the balance in causal interactions between brain regions in different conditions. Paradoxically, the hierarchy is flatter for movie-watching, and the level of nonreversibility is significantly smaller in comparison to both rest and tasks, where the latter in turn have the highest levels of hierarchy and nonreversibility. The underlying mechanisms were revealed by the model-based generative effective connectivity (GEC). Naturalistic films could therefore provide a fast and convenient way to measure important changes in GEC (integrating functional and anatomical connectivity) found in, for example, neuropsychiatric disorders. Overall, this study demonstrates the benefits of moving toward a more naturalistic neuroscience.

The thalamus is the causal hub of intervention in patients with major depressive disorder: Evidence from the Granger causality analysis

Major depressive disorder (MDD) is the leading mental disorder and afflicts more than 350 million people worldwide. The underlying neural mechanisms of MDD remain unclear, hindering the accurate treatment. Recent brain imaging studies have observed functional abnormalities in multiple brain regions in patients with MDD, identifying core brain regions is the key to locating potential therapeutic targets for MDD. The Granger causality analysis (GCA) measures directional effects between brain regions and, therefore, can track causal hubs as potential intervention targets for MDD. We reviewed literature employing GCA to investigate abnormal brain connections in patients with MDD. The total degree of effective connections in the thalamus (THA) is more than twice that in traditional targets such as the superior frontal gyrus and anterior cingulate cortex. Altered causal connections in patients with MDD mainly included enhanced bottom-up connections from the thalamus to various cortical and subcortical regions and reduced top-down connections from these regions to the THA, indicating excessive uplink sensory information and insufficient downlink suppression information for negative emotions. We suggest that the thalamus is the most crucial causal hub for MDD, which may serve as the downstream target for non-invasive brain stimulation and medication approaches in MDD treatment.

Altered nonlinear Granger causality interactions in the large-scale brain networks of patients with schizophrenia

Objective.It has been demonstrated that schizophrenia (SZ) is characterized by functional dysconnectivity involving extensive brain networks. However, the majority of previous studies utilizing resting-state functional magnetic resonance imaging (fMRI) to infer abnormal functional connectivity (FC) in patients with SZ have focused on the linear correlation that one brain region may influence another, ignoring the inherently nonlinear properties of fMRI signals.Approach. In this paper, we present a neural Granger causality (NGC) technique for examining the changes in SZ's nonlinear causal couplings. We develop static and dynamic NGC-based analyses of large-scale brain networks at several network levels, estimating complicated temporal and causal relationships in SZ patients.Main results. We find that the NGC-based FC matrices can detect large and significant differences between the SZ and healthy control groups at both the regional and subnetwork scales. These differences are persistent and significantly overlapped at various network sparsities regardless of whether the brain networks were built using static or dynamic techniques. In addition, compared to controls, patients with SZ exhibited extensive NGC confusion patterns throughout the entire brain.Significance. These findings imply that the NGC-based FCs may be a useful method for quantifying the abnormalities in the causal influences of patients with SZ, hence shedding fresh light on the pathophysiology of this disorder.

Dynamic neural reconfiguration for distinct strategies during competitive social interactions

Information exchange between brain regions is key to understanding information processing for social decision-making, but most analyses ignore its dynamic nature. New insights on this dynamic might help us to uncover the neural correlates of social cognition in the healthy population and also to understand the malfunctioning neural computations underlying dysfunctional social behavior in patients with mental disorders. In this work, we used a multi-round bargaining game to detect switches between distinct bargaining strategies in a cohort of 76 healthy participants. These switches were uncovered by dynamic behavioral modeling using the hidden Markov model. Proposing a novel model of dynamic effective connectivity to estimate the information flow between key brain regions, we found a stronger interaction between the right temporoparietal junction (rTPJ) and the right dorsolateral prefrontal cortex (rDLPFC) for the strategic deception compared with the social heuristic strategies. The level of deception was associated with the information flow from the Brodmann area 10 to the rTPJ, and this association was modulated by the rTPJ-to-rDLPFC information flow. These findings suggest that dynamic bargaining strategy is supported by dynamic reconfiguration of the rDLPFC-and-rTPJ interaction during competitive social interactions.

Forecasting brain activity based on models of spatiotemporal brain dynamics: A comparison of graph neural network architectures

Comprehending the interplay between spatial and temporal characteristics of neural dynamics can contribute to our understanding of information processing in the human brain. Graph neural networks (GNNs) provide a new possibility to interpret graph-structured signals like those observed in complex brain networks. In our study we compare different spatiotemporal GNN architectures and study their ability to model neural activity distributions obtained in functional MRI (fMRI) studies. We evaluate the performance of the GNN models on a variety of scenarios in MRI studies and also compare it to a VAR model, which is currently often used for directed functional connectivity analysis. We show that by learning localized functional interactions on the anatomical substrate, GNN-based approaches are able to robustly scale to large network studies, even when available data are scarce. By including anatomical connectivity as the physical substrate for information propagation, such GNNs also provide a multimodal perspective on directed connectivity analysis, offering a novel possibility to investigate the spatiotemporal dynamics in brain networks.

Altered Effective Connectivity of Resting-State Networks by Tai Chi Chuan in Chronic Fatigue Syndrome Patients: A Multivariate Granger Causality Study

Numerous evidence has shown that patients with chronic fatigue syndrome (CFS) have changes in resting brain functional connectivity, but there is no study on the brain network effect of Tai Chi Chuan intervention in CFS. To explore the influence of Tai Chi Chuan exercise on the causal relationship between brain functional networks in patients with CFS, 21 patients with CFS and 19 healthy controls were recruited for resting-state functional magnetic resonance imaging (rs-fMRI) scanning and 36-item Short-Form Health Survey (SF-36) scale assessment before and after 1month-long training in Tai Chi Chuan. We extracted the resting brain networks using the independent component analysis (ICA) method, analyzed the changes of FC in these networks, conducted Granger causality analysis (GCA) on it, and analyzed the correlation between the difference causality value and the SF-36 scale. Compared to the healthy control group, the SF-36 scale scores of patients with CFS were lower at baseline. Meanwhile, the causal relationship between sensorimotor network (SMN) and default mode network (DMN) was weakened. The above abnormalities could be improved by Tai Chi Chuan training for 1 month. In addition, the correlation analyses showed that the causal relationship between SMN and DMN was positively correlated with the scores of Role Physical (RP) and Bodily Pain (BP) in CFS patients, and the change of causal relationship between SMN and DMN before and after training was positively correlated with the change of BP score. The findings suggest that Tai Chi Chuan is helpful to improve the quality of life for patients with CFS. The change of Granger causality between SMN and DMN may be a readout parameter of CFS. Tai Chi Chuan may promote the functional plasticity of brain networks in patients with CFS by regulating the information transmission between them.

Intrusive experiences in posttraumatic stress disorder: Treatment response induces changes in the directed functional connectivity of the anterior insula

•

Many causal paths were less influenced by the AI after effective therapy for PTSD.

•

Insular influences over the rest of the brain were found to be positively correlated with re-experiencing.

•

Re-experiencing was linked with changes in intrinsic networks’ spatial stability after treatment.

Background

One of the core features of posttraumatic stress disorder (PTSD) is re-experiencing trauma. The anterior insula (AI) has been proposed to play a crucial role in these intrusive experiences. However, the dynamic function of the AI in re-experiencing trauma and its putative modulation by effective therapy need to be specified.

Methods

Thirty PTSD patients were enrolled and exposed to traumatic memory reactivation therapy. Resting-state functional magnetic resonance imaging (fMRI) scans were acquired before and after treatment. To explore AI-directed influences over the rest of the brain, we referred to a mixed model using pre-/posttreatment Granger causality analysis seeded on the AI as a within-subject factor and treatment response as a between-subject factor. To further identify correlates of re-experiencing trauma, we investigated how intrusive severity affected (i) causality maps and (ii) the spatial stability of other intrinsic brain networks.

Results

We observed changes in AI-directed functional connectivity patterns in PTSD patients. Many within- and between-network causal paths were found to be less influenced by the AI after effective therapy. Insular influences were found to be positively correlated with re-experiencing symptoms, while they were linked with a stronger default mode network (DMN) and more unstable central executive network (CEN) connectivity.

Conclusion

We showed that directed changes in AI signaling to the DMN and CEN at rest may underlie the degree of re-experiencing symptoms in PTSD. A positive response to treatment further induced changes in network-to-network anticorrelated patterns. Such findings may guide targeted neuromodulation strategies in PTSD patients not suitably improved by conventional treatment.

Causal Analysis of Activity in Social Brain Areas During Human-Agent Conversation

This article investigates the differences in cognitive and neural mechanisms between human-human and human-virtual agent interaction using a dataset recorded in an ecologically realistic environment. We use Convergent Cross Mapping (CCM) to investigate functional connectivity between pairs of regions involved in the framework of social cognitive neuroscience, namely the fusiform gyrus, superior temporal sulcus (STS), temporoparietal junction (TPJ), and the dorsolateral prefrontal cortex (DLPFC)-taken as prefrontal asymmetry. Our approach is a compromise between investigating local activation in specific regions and investigating connectivity networks that may form part of larger networks. In addition to concording with previous studies, our results suggest that the right TPJ is one of the most reliable areas for assessing processes occurring during human-virtual agent interactions, both in a static and dynamic sense.

Statistical Perspective on Functional and Causal Neural Connectomics: A Comparative Study

Representation of brain network interactions is fundamental to the translation of neural structure to brain function. As such, methodologies for mapping neural interactions into structural models, i.e., inference of functional connectome from neural recordings, are key for the study of brain networks. While multiple approaches have been proposed for functional connectomics based on statistical associations between neural activity, association does not necessarily incorporate causation. Additional approaches have been proposed to incorporate aspects of causality to turn functional connectomes into causal functional connectomes, however, these methodologies typically focus on specific aspects of causality. This warrants a systematic statistical framework for causal functional connectomics that defines the foundations of common aspects of causality. Such a framework can assist in contrasting existing approaches and to guide development of further causal methodologies. In this work, we develop such a statistical guide. In particular, we consolidate the notions of associations and representations of neural interaction, i.e., types of neural connectomics, and then describe causal modeling in the statistics literature. We particularly focus on the introduction of directed Markov graphical models as a framework through which we define the Directed Markov Property—an essential criterion for examining the causality of proposed functional connectomes. We demonstrate how based on these notions, a comparative study of several existing approaches for finding causal functional connectivity from neural activity can be conducted. We proceed by providing an outlook ahead regarding the additional properties that future approaches could include to thoroughly address causality.

Altered effective connectivity from the pregenual anterior cingulate cortex to the laterobasal amygdala mediates the relationship between internet gaming disorder and loneliness

BACKGROUND

Individual with internet gaming disorder (IGD) often experience a high level of loneliness, and neuroimaging studies have demonstrated that amygdala function is associated with both IGD and loneliness. However, the neurobiological basis underlying these relationships remains unclear.

METHODS

In the current study, Granger causal analysis was performed to investigate amygdalar subdivision-based resting-state effective connectivity differences between 111 IGD subjects and 120 matched participants with recreational game use (RGUs). We further correlated neuroimaging findings with clinical measures. Mediation analysis was conducted to explore whether amygdalar subdivision-based effective connectivity mediated the relationship between IGD severity and loneliness.

RESULTS

Compared with RGUs, IGD subjects showed inhibitory effective connections from the left pregenual anterior cingulate cortex (pACC) to the left laterobasal amygdala (LBA) and from the right medial prefrontal cortex (mPFC) to the left LBA, as well as an excitatory effective connection from the left middle prefrontal gyrus (MFG) to the right superficial amygdala. Further analyses demonstrated that the left pACC-left LBA effective connection was negatively correlated with both Internet Addiction Test and UCLA Loneliness scores, and it mediated the relationship between the two.

CONCLUSION

IGD subjects and RGUs showed different connectivity patterns involving amygdalar subdivisions. These findings support a neurobiological mechanism for the relationship between IGD and loneliness, and suggest targets for therapeutic approaches that could be used to treat IGD.

Stimulus evoked causality estimation in stereo-EEG

Objective.Stereo-electroencephalography (SEEG) has recently gained importance in analyzing brain functions. Its high temporal resolution and spatial specificity make it a powerful tool to investigate the strength, direction, and spectral content of brain networks interactions, especially when these connections are stimulus-evoked. However, choosing the best approach to evaluate the flow of information is not trivial, due to the lack of validated methods explicitly designed for SEEG.Approach.We propose a novel non-parametric statistical test for event-related causality (ERC) assessment on SEEG recordings. Here, we refer to the ERC as the causality evoked by a particular part of the stimulus (a response window (RW)). We also present a data surrogation method to evaluate the performance of a causality estimation algorithm. We finally validated our pipeline using surrogate SEEG data derived from an experimentally collected dataset, and compared the most used and successful measures to estimate effective connectivity, belonging to the Geweke-Granger causality framework.Main results.Here we show that our workflow correctly identified all the directed connections in the RW of the surrogate data and prove the robustness of the procedure against synthetic noise with amplitude exceeding physiological-plausible values. Among the causality measures tested, partial directed coherence performed best.Significance.This is the first non-parametric statistical test for ERC estimation explicitly designed for SEEG datasets. The pipeline, in principle, can also be applied to the analysis of any type of time-varying estimator, if there exists a clearly defined RW.

Multivariate semi-blind deconvolution of fMRI time series

Whole brain estimation of the haemodynamic response function (HRF) in functional magnetic resonance imaging (fMRI) is critical to get insight on the global status of the neurovascular coupling of an individual in healthy or pathological condition. Most of existing approaches in the literature works on task-fMRI data and relies on the experimental paradigm as a surrogate of neural activity, hence remaining inoperative on resting-stage fMRI (rs-fMRI) data. To cope with this issue, recent works have performed either a two-step analysis to detect large neural events and then characterize the HRF shape or a joint estimation of both the neural and haemodynamic components in an univariate fashion. In this work, we express the neural activity signals as a combination of piece-wise constant temporal atoms associated with sparse spatial maps and introduce an haemodynamic parcellation of the brain featuring a temporally dilated version of a given HRF model in each parcel with unknown dilation parameters. We formulate the joint estimation of the HRF shapes and spatio-temporal neural representations as a multivariate semi-blind deconvolution problem in a paradigm-free setting and introduce constraints inspired from the dictionary learning literature to ease its identifiability. A fast alternating minimization algorithm, along with its efficient implementation, is proposed and validated on both synthetic and real rs-fMRI data at the subject level. To demonstrate its significance at the population level, we apply this new framework to the UK Biobank data set, first for the discrimination of haemodynamic territories between balanced groups (n=24 individuals in each) patients with an history of stroke and healthy controls and second, for the analysis of normal aging on the neurovascular coupling. Overall, we statistically demonstrate that a pathology like stroke or a condition like normal brain aging induce longer haemodynamic delays in certain brain areas (e.g. Willis polygon, occipital, temporal and frontal cortices) and that this haemodynamic feature may be predictive with an accuracy of 74 % of the individual's age in a supervised classification task performed on n=459 subjects.

A graph neural network framework for causal inference in brain networks

A central question in neuroscience is how self-organizing dynamic interactions in the brain emerge on their relatively static structural backbone. Due to the complexity of spatial and temporal dependencies between different brain areas, fully comprehending the interplay between structure and function is still challenging and an area of intense research. In this paper we present a graph neural network (GNN) framework, to describe functional interactions based on the structural anatomical layout. A GNN allows us to process graph-structured spatio-temporal signals, providing a possibility to combine structural information derived from diffusion tensor imaging (DTI) with temporal neural activity profiles, like that observed in functional magnetic resonance imaging (fMRI). Moreover, dynamic interactions between different brain regions discovered by this data-driven approach can provide a multi-modal measure of causal connectivity strength. We assess the proposed model's accuracy by evaluating its capabilities to replicate empirically observed neural activation profiles, and compare the performance to those of a vector auto regression (VAR), like that typically used in Granger causality. We show that GNNs are able to capture long-term dependencies in data and also computationally scale up to the analysis of large-scale networks. Finally we confirm that features learned by a GNN can generalize across MRI scanner types and acquisition protocols, by demonstrating that the performance on small datasets can be improved by pre-training the GNN on data from an earlier study. We conclude that the proposed multi-modal GNN framework can provide a novel perspective on the structure-function relationship in the brain. Accordingly this approach appears to be promising for the characterization of the information flow in brain networks.

Revisiting the global workspace orchestrating the hierarchical organization of the human brain

A central challenge in neuroscience is how the brain organizes the information necessary to orchestrate behaviour. Arguably, this whole-brain orchestration is carried out by a core subset of integrative brain regions, a 'global workspace', but its constitutive regions remain unclear. We quantified the global workspace as the common regions across seven tasks as well as rest, in a common 'functional rich club'. To identify this functional rich club, we determined the information flow between brain regions by means of a normalized directed transfer entropy framework applied to multimodal neuroimaging data from 1,003 healthy participants and validated in participants with retest data. This revealed a set of regions orchestrating information from perceptual, long-term memory, evaluative and attentional systems. We confirmed the causal significance and robustness of our results by systematically lesioning a generative whole-brain model. Overall, this framework describes a complex choreography of the functional hierarchical organization of the human brain.

Mindfulness-based stress reduction training modulates striatal and cerebellar connectivity

Mindfulness is a meditation practice frequently associated with changes in subjective evaluation of cognitive and sensorial experience, as well as with modifications of brain activity and morphometry. Aside from the anatomical localization of functional changes induced by mindfulness practice, little is known about changes in functional and effective functional magnetic resonance imaging (fMRI) connectivity. Here we performed a connectivity fMRI analysis in a group of healthy individuals participating in an 8-week mindfulness-based stress reduction (MBSR) training program. Data from both a "mind-wandering" and a "meditation" state were acquired before and after the MBSR course. Results highlighted decreased local connectivity after training in the right anterior putamen and insula during spontaneous mind-wandering and the right cerebellum during the meditative state. A further effective connectivity analysis revealed (a) decreased modulation by the anterior cingulate cortex over the anterior portion of the putamen, and (b) a change in left and right posterior putamen excitatory input and inhibitory output with the cerebellum, respectively. Results suggest a rearrangement of dorsal striatum functional and effective connectivity in response to mindfulness practice, with changes in cortico-subcortical-cerebellar modulatory dynamics. Findings might be relevant for the understanding of widely documented mindfulness behavioral effects, especially those related to pain perception.

Interrelations between dopamine and serotonin producing sites and regions of the default mode network

Recent functional magnetic resonance imaging (fMRI) studies showed that blood oxygenation level-dependent (BOLD) signal fluctuations in the default mode network (DMN) are functionally tightly connected to those in monoaminergic nuclei, producing dopamine (DA), and serotonin (5-HT) transmitters, in the midbrain/brainstem. We combined accelerated fMRI acquisition with spectral Granger causality and coherence analysis to investigate causal relationships between these areas. Both methods independently lead to similar results and confirm the existence of a top-down information flow in the resting-state condition, where activity in core DMN areas influences activity in the neuromodulatory centers producing DA/5-HT. We found that latencies range from milliseconds to seconds with high inter-subject variability, likely attributable to the resting condition. Our novel findings provide new insights into the functional organization of the human brain.

Association between emotional intelligence and effective brain connectome: A large-scale spectral DCM study

INTRODUCTION

Emotional Intelligence (EI) is a well-documented aspect of social and interpersonal functioning, but the underlying neural mechanisms for this capacity remain poorly understood. Here we used advanced brain connectivity techniques to explore the associations between EI and effective connectivity (EC) within four functional brain networks.

METHODS

The Mayer-Salovey-Caruso Emotional Intelligence Test (MSCEIT) was used to collect EI data from 55 healthy individuals (mean age = 30.56±8.3 years, 26 males). The MSCEIT comprises two area cores - experiential EI (T1) and strategic EI (T2). The T1 core included two sub-scales - perception of emotions (S1) and using emotions to facilitate thinking (S2), and the T2 core included two sub-scales - understanding of emotions (S3) and management of emotions (S4). All participants underwent structural and resting-state functional magnetic resonance imaging (rsfMRI) scans. The spectral dynamic causal modeling approach was implemented to estimate EC within four networks of interest - the default-mode network (DMN), dorsal attention network (DAN), control-execution network (CEN) and salience network (SN). The strength of EC within each network was correlated with the measures of EI, with correlations at p_FDR < 0.05 considered as significant.

RESULTS

There was no significant association between any of the measures of EI and EC strength within the DMN and DAN. For CEN, however, we found that there were significant negative associations between EC strength from the right anterior prefrontal cortex (RAPFC) to the left anterior prefrontal cortex (LAPFC) and both S2 and T1, and significant positive associations between EC strength from LAPFC to RAPFC and S2. EC strength from the right superior parietal cortex (SPC) to RAPFC also showed significant negative association with S4 and T2. For the SN, S3 showed significant negative association with EC strength from the right insula to RAPFC and significant positive association with EC strength from the left insula to dorsal anterior cingulate cortex (DACC).

CONCLUSIONS

We provide evidence that the negative ECs within the right hemisphere, and from the right to left hemisphere, and positive ECs within the left hemisphere and from the left to right hemisphere of CEN (involving bilateral frontal and right parietal region) and SN (involving right frontal, anterior cingulate and bilateral insula) play a significant role in regulating and processing emotions. These findings also suggest that measures of EC can be utilized as important biomarkers to better understand the underlying neural mechanisms of EI.

Fronto-cerebellar connectivity mediating cognitive processing speed

Processing speed is an important construct in understanding cognition. This study was aimed to control task specificity for understanding the neural mechanisms underlying cognitive processing speed. Forty young adult subjects performed attention tasks of two modalities (auditory and visual) and two levels of task rules (compatible and incompatible). Block-design fMRI captured BOLD signals during the tasks. Thirteen regions of interest were defined with reference to publicly available activation maps for processing speed tasks. Cognitive speed was derived from task reaction times, which yielded six sets of connectivity measures. Mixed-effect LASSO regression revealed six significant paths suggestive of a cerebello-frontal network predicting the cognitive speed. Among them, three are long range (two fronto-cerebellar, one cerebello-frontal), and three are short range (fronto-frontal, cerebello-cerebellar, and cerebello-thalamic). The long-range connections are likely to relate to cognitive control, and the short-range connections relate to rule-based stimulus-response processes. The revealed neural network suggests that automaticity, acting on the task rules and interplaying with effortful top-down attentional control, accounts for cognitive speed.

A Bayesian method for inference of effective connectivity in brain networks for detecting the Mozart effect

Several studies claim that listening to Mozart music affects cognition and can be used to treat neurological conditions like epilepsy. Research into this Mozart effect has not addressed how dynamic interactions between brain networks, i.e. effective connectivity, are affected. The Granger-causality analysis is often used to infer effective connectivity. First, we investigate if a new method, Bayesian topology identification, can be used as an alternative. Both methods are evaluated on simulation data, where the Bayesian method outperforms the Granger-causality analysis in the inference of connectivity graphs of dynamic networks, especially for short data lengths. In the second part, the Bayesian method is extended to enable the inference of changes in effective connectivity between groups of subjects. Next, we apply both methods to fMRI scans of 16 healthy subjects, who were scanned before and after the exposure to Mozart's sonata K448 at least 2 hours a day for 7 days. Here, we investigate if the effective connectivity of the subjects significantly changed after listening to Mozart music. The Bayesian method detected changes in effective connectivity between networks related to cognitive processing and control in the connection from the central executive to the superior sensori-motor network, in the connection from the posterior default mode to the fronto-parietal right network, and in the connection from the anterior default mode to the dorsal attention network. This last connection was only detected in a subgroup of subjects with a longer listening duration. Only in this last connection, an effect was found by the Granger-causality analysis.

Whole-brain estimates of directed connectivity for human connectomics

Connectomics is essential for understanding large-scale brain networks but requires that individual connection estimates are neurobiologically interpretable. In particular, a principle of brain organization is that reciprocal connections between cortical areas are functionally asymmetric. This is a challenge for fMRI-based connectomics in humans where only undirected functional connectivity estimates are routinely available. By contrast, whole-brain estimates of effective (directed) connectivity are computationally challenging, and emerging methods require empirical validation. Here, using a motor task at 7T, we demonstrate that a novel generative model can infer known connectivity features in a whole-brain network (>200 regions, >40,000 connections) highly efficiently. Furthermore, graph-theoretical analyses of directed connectivity estimates identify functional roles of motor areas more accurately than undirected functional connectivity estimates. These results, which can be achieved in an entirely unsupervised manner, demonstrate the feasibility of inferring directed connections in whole-brain networks and open new avenues for human connectomics.

Partial Similarity Reveals Dynamics in Brainstem-Midbrain Networks during Trigeminal Nociception

Imaging studies help us understand the important role of brainstem and midbrain regions in human trigeminal pain processing without solving the question of how these regions actually interact. In the current study, we describe this connectivity and its dynamics during nociception with a novel analytical approach called Partial Similarity (PS). We developed PS specifically to estimate the communication between individual hubs of the network in contrast to the overall communication within that network. Partial Similarity works on trial-to-trial variance of neuronal activity acquired with functional magnetic resonance imaging. It discovers direct communication between two hubs considering the remainder of the network as confounds. A similar method to PS is Representational Similarity, which works with ordinary correlations and does not consider any external influence on the communication between two hubs. Particularly the combination of Representational Similarity and Partial Similarity analysis unravels brainstem dynamics involved in trigeminal pain using the spinal trigeminal nucleus (STN)—the first relay station of peripheral trigeminal input—as a seed region. The combination of both methods can be valuable tools in discovering the network dynamics in fMRI and an important instrument for future insight into the nature of various neurological diseases like primary headaches.

Modeling time-varying brain networks with a self-tuning optimized Kalman filter

Brain networks are complex dynamical systems in which directed interactions between different areas evolve at the sub-second scale of sensory, cognitive and motor processes. Due to the highly non-stationary nature of neural signals and their unknown noise components, however, modeling dynamic brain networks has remained one of the major challenges in contemporary neuroscience. Here, we present a new algorithm based on an innovative formulation of the Kalman filter that is optimized for tracking rapidly evolving patterns of directed functional connectivity under unknown noise conditions. The Self-Tuning Optimized Kalman filter (STOK) is a novel adaptive filter that embeds a self-tuning memory decay and a recursive regularization to guarantee high network tracking accuracy, temporal precision and robustness to noise. To validate the proposed algorithm, we performed an extensive comparison against the classical Kalman filter, in both realistic surrogate networks and real electroencephalography (EEG) data. In both simulations and real data, we show that the STOK filter estimates time-frequency patterns of directed connectivity with significantly superior performance. The advantages of the STOK filter were even clearer in real EEG data, where the algorithm recovered latent structures of dynamic connectivity from epicranial EEG recordings in rats and human visual evoked potentials, in excellent agreement with known physiology. These results establish the STOK filter as a powerful tool for modeling dynamic network structures in biological systems, with the potential to yield new insights into the rapid evolution of network states from which brain functions emerge.

During normal behavior, brains transition between functional network states several times per second. This allows humans to quickly read a sentence, and a frog to catch a fly. Understanding these fast network dynamics is fundamental to understanding how brains work, but up to now it has proven very difficult to model fast brain dynamics for various methodological reasons. To overcome these difficulties, we designed a new Kalman filter (STOK) by innovating on previous solutions from control theory and state-space modelling. We show that STOK accurately models fast network changes in simulations and real neural data, making it an essential new tool for modelling fast brain networks in the time and frequency domain.

Effective connectivity of the right anterior insula in schizophrenia: The salience network and task-negative to task-positive transition

Triple network dysfunction theory of schizophrenia postulates that the interaction between the default-mode and the fronto-parietal executive network is disrupted by aberrant salience signals from the right anterior insula (rAI). To date, it is not clear how the proposed resting-state disruption translates to task-processing inefficiency in subjects with schizophrenia. Using a contiguous resting and 2-back task performance fMRI paradigm, we quantified the change in effective connectivity that accompanies rest-to-task state transition in 29 clinically stable patients with schizophrenia and 31 matched healthy controls. We found an aberrant task-evoked increase in the influence of the rAI to both executive (Cohen's d = 1.35, p = 2.8 × 10-6) and default-mode (Cohen's d = 1.22, p = 1.5 × 10-5) network regions occur in patients when compared to controls. In addition, the effective connectivity from middle occipital gyrus (dorsal visual cortex) to insula is also increased in patients as compared with healthy controls. Aberrant insula to executive network influence is pronounced in patients with more severe negative symptom burden. These findings suggest that control signals from rAI are abnormally elevated and directed towards both task-positive and task-negative brain regions, when task-related demands arise in schizophrenia. This aberrant, undiscriminating surge in salience signalling may disrupt contextually appropriate allocation of resources in the neuronal workspace in patients with schizophrenia.

The boundaries of state-space Granger causality analysis applied to BOLD simulated data: A comparative modelling and simulation approach

BACKGROUND

The analysis of connectivity has become a fundamental tool in human neuroscience. Granger Causality Mapping is a data-driven method that uses Granger Causality (GC) to assess the existence and direction of influence between signals, based on temporal precedence of information. More recently, a theory of Granger causality has been developed for state-space (SS-GC) processes, but little is known about its statistical validation and application on functional magnetic resonance imaging (fMRI) data.

NEW METHOD

We explored different multivariate computational frameworks to define the optimal combination for GC estimation. We hypothesized a new heuristic, combining SS-GC with a distinct statistical validation technique, Time Reversed Testing, validating it on synthetic data. We test its performance with a number of experimental parameters, including block structure, sampling frequency, noise and system mean pairwise correlation, using a statistical framework of binary classification.

RESULTS

We found that SS-GC with time reversed testing outperforms other frameworks. The results validate the application of SS-GC to generative models. When estimating reliable causal relations, SS-GC returns promising results, especially when considering synthetic data with a high impact of noise and sampling rate.

CONCLUSIONS

In this study, we empirically explored the boundaries of SS-GC with time reversed testing, a data-driven causality analysis framework with potential applicability to fMRI data.

Construction of functional brain connectivity networks from fMRI data with driving and modulatory inputs: an extended conditional Granger causality approach

We propose a numerical-based approach extending the conditional MVAR Granger causality (MVGC) analysis for the construction of directed connectivity networks in the presence of both exogenous/stimuli and modulatory inputs. The performance of the proposed scheme is validated using both synthetic stochastic data considering also the influence of haemodynamics latencies and a benchmark fMRI dataset related to the role of attention in the perception of visual motion. The particular fMRI dataset has been used in many studies to evaluate alternative model hypotheses using the Dynamic Causal Modelling (DCM) approach. Based on the use of the Bayes factor, we show that the obtained GC connectivity network compares well to a reference model that has been selected through DCM analysis among other candidate models. Thus, our findings suggest that the proposed scheme can be successfully used as a stand-alone or complementary to DCM approach to find directed causal connectivity patterns in task-related fMRI studies.

Kernel Granger Causality Based on Back Propagation Neural Network Fuzzy Inference System on fMRI Data

Granger causality (GC) is one of the most popular measures to investigate causality influence among brain regions and has been achieved significant results for exploring brain networks based on functional magnetic resonance imaging (fMRI). However, the predictors and order selection of conventional GC are based on linear models which result in such restrictions as poorly detection of nonlinearity and so on, in the application. This paper proposes a novel GC model called back propagation (BP) based kernel function Granger causality (BP_KFGC), in which symplectic geometry is used for embedding dimension and fuzzy inference system for predicting time series. The proposed method doesn't depend on the prediction of the vector auto-regression model, so that time series don't need to be wide-sense stationary as linear GC and kernel GC. In addition, it is a multivariate approach which is applicable to both linear and nonlinear systems and eliminates the effects of latent variables. The performance of the new method is evaluated and compared with linear GC, partial GC, neural network GC and kernel GC by simulated data with multiple adjustments to the nonlinearity. The results show that BP_KFGC outperforms the other four methods in detecting both linear and nonlinear causalities. Furthermore, we applied BP_KFGC to construct directed weight network (DWN) of Alzheimer's disease (AD) patients and health controls (HCs), and then nine graph-based features of DWN were used for classification by the classifier of support vector machine with radial basis kernel function. The accuracy of 95.89%, sensitivity of 93.31%, and specificity of 94.97% were achieved which may provide an auxiliary mean for the clinical diagnosis of AD.

Decoding Task-Specific Cognitive States with Slow, Directed Functional Networks in the Human Brain

Flexible functional interactions among brain regions mediate critical cognitive functions. Such interactions can be measured using functional magnetic resonance imaging (fMRI) data either with instantaneous (zero-lag) or lag-based (time-lagged) functional connectivity. Because the fMRI hemodynamic response is slow, and is sampled at a timescale (seconds) several orders of magnitude slower than the underlying neural dynamics (milliseconds), simulation studies have shown that lag-based fMRI functional connectivity, measured with approaches like Granger–Geweke causality (GC), provides spurious and unreliable estimates of underlying neural interactions. Experimental verification of this claim is challenging because neural ground truth connectivity is often unavailable concurrently with fMRI recordings. Here we demonstrate that, despite these widely held caveats, GC networks estimated from fMRI recordings contain useful information for classifying task-specific cognitive states. We estimated instantaneous and lag-based GC functional connectivity networks using fMRI data from 1000 participants (Human Connectome Project database). A linear classifier, trained on either instantaneous or lag-based GC, reliably discriminated among seven different task and resting brain states, with >80% cross-validation accuracy. With network simulations, we demonstrate that instantaneous and lag-based GC exploited interactions at fast and slow timescales, respectively, to achieve robust classification. With human fMRI data, instantaneous and lag-based GC identified complementary, task–core networks. Finally, variations in GC connectivity explained inter-individual variations in a variety of cognitive scores. Our findings show that instantaneous and lag-based methods reveal complementary aspects of functional connectivity in the brain, and suggest that slow, directed functional interactions, estimated with fMRI, may provide useful markers of behaviorally relevant cognitive states.

Model-based whole-brain effective connectivity to study distributed cognition in health and disease

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.

Converging Resting State Networks Unravels Potential Remote Effects of Transcranial Magnetic Stimulation for Major Depression

Despite being a commonly used protocol to treat major depressive disorder (MDD), the underlying mechanism of repetitive transcranial magnetic stimulation (rTMS) on dorsolateral prefrontal cortex (DLPFC) remains unclear. In the current study, we investigated the resting-state fMRI data of 100 healthy subjects by exploring three overlapping functional networks associated with the psychopathologically MDD-related areas (the nucleus accumbens, amygdala, and ventromedial prefrontal cortex). Our results showed that these networks converged at the bilateral DLPFC, which suggested that rTMS over DLPFC might improve MDD by remotely modulating the MDD-related areas synergistically. Additionally, they functionally converged at the DMPFC and bilateral insula which are known to be associated with MDD. These two areas could also be potential targets for rTMS treatment. Dynamic causal modelling (DCM) and Granger causality analysis (GCA) revealed that all pairwise connections among bilateral DLPFC, DMPFC, bilateral insula, and three psychopathologically MDD-related areas contained significant causality. The DCM results also suggested that most of the functional interactions between MDD-related areas and bilateral DLPFC, DMPFC, and bilateral insula can predominantly be explained by the effective connectivity from the psychopathologically MDD-related areas to the rTMS stimulation sites. Finally, we found the conventional functional connectivity to be a more representative measure to obtain connectivity parameters compared to GCA and DCM analysis. Our research helped inspecting the convergence of the functional networks related to a psychiatry disorder. The results identified potential targets for brain stimulation treatment and contributed to the optimization of patient-specific brain stimulation protocols.

Decreased directed functional connectivity in the psychedelic state

Neuroimaging studies of the psychedelic state offer a unique window onto the neural basis of conscious perception and selfhood. Despite well understood pharmacological mechanisms of action, the large-scale changes in neural dynamics induced by psychedelic compounds remain poorly understood. Using source-localised, steady-state MEG recordings, we describe changes in functional connectivity following the controlled administration of LSD, psilocybin and low-dose ketamine, as well as, for comparison, the (non-psychedelic) anticonvulsant drug tiagabine. We compare both undirected and directed measures of functional connectivity between placebo and drug conditions. We observe a general decrease in directed functional connectivity for all three psychedelics, as measured by Granger causality, throughout the brain. These data support the view that the psychedelic state involves a breakdown in patterns of functional organisation or information flow in the brain. In the case of LSD, the decrease in directed functional connectivity is coupled with an increase in undirected functional connectivity, which we measure using correlation and coherence. This surprising opposite movement of directed and undirected measures is of more general interest for functional connectivity analyses, which we interpret using analytical modelling. Overall, our results uncover the neural dynamics of information flow in the psychedelic state, and highlight the importance of comparing multiple measures of functional connectivity when analysing time-resolved neuroimaging data.