A Unified Estimation Framework for State-Related Changes in Effective Brain Connectivity

The research progress on effective connectivity in adolescent depression based on resting-state fMRI

Introduction

The brain's spontaneous neural activity can be recorded during rest using resting state functional magnetic resonance imaging (rs-fMRI), and intricate brain functional networks and interaction patterns can be discovered through correlation analysis. As a crucial component of rs-fMRI analysis, effective connectivity analysis (EC) may provide a detailed description of the causal relationship and information flow between different brain areas. It has been very helpful in identifying anomalies in the brain activity of depressed teenagers.

Methods

This study explored connectivity abnormalities in brain networks and their impact on clinical symptoms in patients with depression through resting state functional magnetic resonance imaging (rs-fMRI) and effective connectivity (EC) analysis. We first introduce some common EC analysis methods, discuss their application background and specific characteristics.

Results

EC analysis reveals information flow problems between different brain regions, such as the default mode network, the central executive network, and the salience network, which are closely related to symptoms of depression, such as low mood and cognitive impairment. This review discusses the limitations of existing studies while summarizing the current applications of EC analysis methods. Most of the early studies focused on the static connection mode, ignoring the causal relationship between brain regions. However, effective connection can reflect the upper and lower relationship of brain region interaction, and provide help for us to explore the mechanism of neurological diseases. Existing studies focus on the analysis of a single brain network, but rarely explore the interaction between multiple key networks.

Discussion

To do so, we can address these issues by integrating multiple technologies. The discussion of these issues is reflected in the text. Through reviewing various methods and applications of EC analysis, this paper aims to explore the abnormal connectivity patterns of brain networks in patients with depression, and further analyze the relationship between these abnormalities and clinical symptoms, so as to provide more accurate theoretical support for early diagnosis and personalized treatment of depression.

A Deep Dynamic Causal Learning Model to Study Changes in Dynamic Effective Connectivity During Brain Development

OBJECTIVE

Brain dynamic effective connectivity (dEC), characterizes the information transmission patterns between brain regions that change over time, which provides insight into the biological mechanism underlying brain development. However, most existing methods predominantly capture fixed or temporally invariant EC, leaving dEC largely unexplored.

METHODS

Herein we propose a deep dynamic causal learning model specifically designed to capture dEC. It includes a dynamic causal learner to detect time-varying causal relationships from spatio-temporal data, and a dynamic causal discriminator to validate these findings by comparing original and reconstructed data.

RESULTS

Our model outperforms established baselines in the accuracy of identifying dynamic causalities when tested on the simulated data. When applied to the Philadelphia Neurodevelopmental Cohort, the model uncovers distinct patterns in dEC networks across different age groups. Specifically, the evolution process of brain dEC networks in young adults is more stable than in children, and significant differences in information transfer patterns exist between them.

CONCLUSION

This study highlights the brain's developmental trajectory, where networks transition from undifferentiated to specialized structures with age, in accordance with the improvement of an individual's cognitive and information processing capability.

SIGNIFICANCE

The proposed model consists of the identification and verification of dynamic causality, utilizing the spatio-temporal fusing information from fMRI. As a result, it can accurately detect dEC and characterize its evolution over age.

A Survey on Brain Effective Connectivity Network Learning

Human brain effective connectivity characterizes the causal effects of neural activities among different brain regions. Studies of brain effective connectivity networks (ECNs) for different populations contribute significantly to the understanding of the pathological mechanism associated with neuropsychiatric diseases and facilitate finding new brain network imaging markers for the early diagnosis and evaluation for the treatment of cerebral diseases. A deeper understanding of brain ECNs also greatly promotes brain-inspired artificial intelligence (AI) research in the context of brain-like neural networks and machine learning. Thus, how to picture and grasp deeper features of brain ECNs from functional magnetic resonance imaging (fMRI) data is currently an important and active research area of the human brain connectome. In this survey, we first show some typical applications and analyze existing challenging problems in learning brain ECNs from fMRI data. Second, we give a taxonomy of ECN learning methods from the perspective of computational science and describe some representative methods in each category. Third, we summarize commonly used evaluation metrics and conduct a performance comparison of several typical algorithms both on simulated and real datasets. Finally, we present the prospects and references for researchers engaged in learning ECNs.

Improving EEG-Based Driver Distraction Classification Using Brain Connectivity Estimators

This paper discusses a novel approach to an EEG (electroencephalogram)-based driver distraction classification by using brain connectivity estimators as features. Ten healthy volunteers with more than one year of driving experience and an average age of 24.3 participated in a virtual reality environment with two conditions, a simple math problem-solving task and a lane-keeping task to mimic the distracted driving task and a non-distracted driving task, respectively. Independent component analysis (ICA) was conducted on the selected epochs of six selected components relevant to the frontal, central, parietal, occipital, left motor, and right motor areas. Granger-Geweke causality (GGC), directed transfer function (DTF), partial directed coherence (PDC), and generalized partial directed coherence (GPDC) brain connectivity estimators were used to calculate the connectivity matrixes. These connectivity matrixes were used as features to train the support vector machine (SVM) with the radial basis function (RBF) and classify the distracted and non-distracted driving tasks. GGC, DTF, PDC, and GPDC connectivity estimators yielded the classification accuracies of 82.27%, 70.02%, 86.19%, and 80.95%, respectively. Further analysis of the PDC connectivity estimator was conducted to determine the best window to differentiate between the distracted and non-distracted driving tasks. This study suggests that the PDC connectivity estimator can yield better classification accuracy for driver distractions.

Effective connectivity in the default mode network after paediatric traumatic brain injury

Children who experience a traumatic brain injury (TBI) are at elevated risk for a range of negative cognitive and neuropsychological outcomes. Identifying which children are at greatest risk for negative outcomes can be difficult due to the heterogeneity of TBI. To address this barrier, the current study applied a novel method of characterizing brain connectivity networks, Bayesian multi-subject vector autoregressive modelling (BVAR-connect), which used white matter integrity as priors to evaluate effective connectivity-the time-dependent relationship in functional magnetic resonance imaging (fMRI) activity between two brain regions-within the default mode network (DMN). In a prospective longitudinal study, children ages 8-15 years with mild to severe TBI underwent diffusion tensor imaging and resting state fMRI 7 weeks after injury; post-concussion and anxiety symptoms were assessed 7 months after injury. The goals of this study were to (1) characterize differences in positive effective connectivity of resting-state DMN circuitry between healthy controls and children with TBI, (2) determine if severity of TBI was associated with differences in DMN connectivity and (3) evaluate whether patterns of DMN effective connectivity predicted persistent post-concussion symptoms and anxiety. Healthy controls had unique positive connectivity that mostly emerged from the inferior temporal lobes. In contrast, children with TBI had unique effective connectivity among orbitofrontal and parietal regions. These positive orbitofrontal-parietal DMN effective connectivity patterns also differed by TBI severity and were associated with persisting behavioural outcomes. Effective connectivity may be a sensitive neuroimaging marker of TBI severity as well as a predictor of chronic post-concussion symptoms and anxiety.

On the characterization of cognitive tasks using activity-specific short-lived synchronization between electroencephalography channels

Accurate characterization of brain activity during a cognitive task is challenging due to the dynamically changing and the complex nature of the brain. The majority of the proposed approaches assume stationarity in brain activity and disregard the systematic timing organization among brain regions during cognitive tasks. In this study, we propose a novel cognitive activity recognition method that captures the activity-specific timing parameters from training data that elicits maximal average short-lived pairwise synchronization between electroencephalography signals. We evaluated the characterization power of the activity-specific timing parameter triplets in a motor imagery activity recognition framework. The activity-specific timing parameter triplets consist of latency of the maximally synchronized signal segments from activity onset Δt, the time lag between maximally synchronized signal segments τ, and the duration of the maximally synchronized signal segments w. We used cosine-based similarity, wavelet bi-coherence, phase-locking value, phase coherence value, linearized mutual information, and cross-correntropy to calculate the channel synchronizations at the specific timing parameters. Recognition performances as well as statistical analyses on both BCI Competition-III dataset IVa and PhysioNet Motor Movement/Imagery dataset, indicate that the inter-channel short-lived synchronization calculated using activity-specific timing parameter triplets elicit significantly distinct synchronization profiles for different motor imagery tasks and can thus reliably be used for cognitive task recognition purposes.

Detecting Dynamic Community Structure in Functional Brain Networks Across Individuals: A Multilayer Approach

OBJECTIVE

We present a unified statistical framework for characterizing community structure of brain functional networks that captures variation across individuals and evolution over time. Existing methods for community detection focus only on single-subject analysis of dynamic networks; while recent extensions to multiple-subjects analysis are limited to static networks.

METHOD

To overcome these limitations, we propose a multi-subject, Markov-switching stochastic block model (MSS-SBM) to identify state-related changes in brain community organization over a group of individuals. We first formulate a multilayer extension of SBM to describe the time-dependent, multi-subject brain networks. We develop a novel procedure for fitting the multilayer SBM that builds on multislice modularity maximization which can uncover a common community partition of all layers (subjects) simultaneously. By augmenting with a dynamic Markov switching process, our proposed method is able to capture a set of distinct, recurring temporal states with respect to inter-community interactions over subjects and the change points between them.

RESULTS

Simulation shows accurate community recovery and tracking of dynamic community regimes over multilayer networks by the MSS-SBM. Application to task fMRI reveals meaningful non-assortative brain community motifs, e.g., core-periphery structure at the group level, that are associated with language comprehension and motor functions suggesting their putative role in complex information integration. Our approach detected dynamic reconfiguration of modular connectivity elicited by varying task demands and identified unique profiles of intra and inter-community connectivity across different task conditions.

CONCLUSION

The proposed multilayer network representation provides a principled way of detecting synchronous, dynamic modularity in brain networks across subjects.

Intralobular and Interlobular Parietal Functional Network Correlated to MI-BCI Performance

Brain-computer interface (BCI) brings hope to patients suffering from neuromuscular diseases, by allowing the control of external devices using neural signals from the central nervous system. However, a portion of individuals was unable to operate BCI with high efficacy. This research aimed to study the brain-wide functional connectivity differences that contributed to BCI performance, and investigate the relationship between task-related connectivity strength and BCI performance. Functional connectivity was estimated using pairwise Pearson's correlation from the EEG of 48 subjects performing left or right hand motor imagery (MI) tasks. The classification accuracy of linear support vector machine (SVM) to distinguish both tasks were used to represent MI-BCI performance. The significant differences in connectivity strengths were examined using Welch's T-test. The association between accuracy and connection strength was studied using correlation model. Three intralobular and fourteen interlobular connections from the parietal lobe showed a correlation of 0.31 and -0.34 respectively. Results indicate that alpha wave connectivity from 8 Hz to 13 Hz was more related to classification performance compared to high-frequency waves. Subject-independent trial-based analysis shows that MI trials executed with stronger intralobular and interlobular parietal connections performed significantly better than trials with weaker connections. Further investigation from an independent MI dataset reveals several similar connections that were correlated with MI-BCI performance. The functional connectivity of the parietal lobe could potentially allow prediction of MI-BCI performance and enable implementation of neurofeedback training for users to improve the usability of MI-BCI.

Statistical model for dynamically-changing correlation matrices with application to brain connectivity

BACKGROUND

Recent studies have indicated that functional connectivity is dynamic even during rest. A common approach to modeling the dynamic functional connectivity in whole-brain resting-state fMRI is to compute the correlation between anatomical regions via sliding time windows. However, the direct use of the sample correlation matrices is not reliable due to the image acquisition and processing noises in resting-sate fMRI.

NEW METHOD

To overcome these limitations, we propose a new statistical model that smooths out the noise by exploiting the geometric structure of correlation matrices. The dynamic correlation matrix is modeled as a linear combination of symmetric positive-definite matrices combined with cosine series representation. The resulting smoothed dynamic correlation matrices are clustered into disjoint brain connectivity states using the k-means clustering algorithm.

RESULTS

The proposed model preserves the geometric structure of underlying physiological dynamic correlation, eliminates unwanted noise in connectivity and obtains more accurate state spaces. The difference in the estimated dynamic connectivity states between males and females is identified.

COMPARISON WITH EXISTING METHODS

We demonstrate that the proposed statistical model has less rapid state changes caused by noise and improves the accuracy in identifying and discriminating different states.

CONCLUSIONS

We propose a new regression model on dynamically changing correlation matrices that provides better performance over existing windowed correlation and is more reliable for the modeling of dynamic connectivity.

A Markov-Switching Model Approach to Heart Sound Segmentation and Classification

OBJECTIVE

We consider challenges in accurate segmentation of heart sound signals recorded under noisy clinical environments for subsequent classification of pathological events. Existing state-of-the-art solutions to heart sound segmentation use probabilistic models such as hidden Markov models (HMMs), which, however, are limited by its observation independence assumption and rely on pre-extraction of noise-robust features.

METHODS

We propose a Markov-switching autoregressive (MSAR) process to model the raw heart sound signals directly, which allows efficient segmentation of the cyclical heart sound states according to the distinct dependence structure in each state. To enhance robustness, we extend the MSAR model to a switching linear dynamic system (SLDS) that jointly model both the switching AR dynamics of underlying heart sound signals and the noise effects. We introduce a novel algorithm via fusion of switching Kalman filter and the duration-dependent Viterbi algorithm, which incorporates the duration of heart sound states to improve state decoding.

RESULTS

Evaluated on Physionet/CinC Challenge 2016 dataset, the proposed MSAR-SLDS approach significantly outperforms the hidden semi-Markov model (HSMM) in heart sound segmentation based on raw signals and comparable to a feature-based HSMM. The segmented labels were then used to train Gaussian-mixture HMM classifier for identification of abnormal beats, achieving high average precision of 86.1% on the same dataset including very noisy recordings.

CONCLUSION

The proposed approach shows noticeable performance in heart sound segmentation and classification on a large noisy dataset.

SIGNIFICANCE

It is potentially useful in developing automated heart monitoring systems for pre-screening of heart pathologies.

Statistical model-based approaches for functional connectivity analysis of neuroimaging data

We present recent literature on model-based approaches to estimating functional connectivity from neuroimaging data. In contrast to the typical focus on a particular scientific question, we reframe a wider literature in terms of the underlying statistical model used. We distinguish between directed versus undirected and static versus time-varying connectivity. There are numerous advantages to a model-based approach, including easily specified inductive bias, handling limited data scenarios, and building complex models from simpler building blocks.

Improved state change estimation in dynamic functional connectivity using hidden semi-Markov models

The study of functional brain networks has grown rapidly over the past decade. While most functional connectivity (FC) analyses estimate one static network structure for the entire length of the functional magnetic resonance imaging (fMRI) time series, recently there has been increased interest in studying time-varying changes in FC. Hidden Markov models (HMMs) have proven to be a useful modeling approach for discovering repeating graphs of interacting brain regions (brain states). However, a limitation lies in HMMs assuming that the sojourn time, the number of consecutive time points in a state, is geometrically distributed. This may encourage inaccurate estimation of the time spent in a state before switching to another state. We propose a hidden semi-Markov model (HSMM) approach for inferring time-varying brain networks from fMRI data, which explicitly models the sojourn distribution. Specifically, we propose using HSMMs to find each subject's most probable series of network states and the graphs associated with each state, while properly estimating and modeling the sojourn distribution for each state. We perform a simulation study, as well as an analysis on both task-based fMRI data from an anxiety-inducing experiment and resting-state fMRI data from the Human Connectome Project. Our results demonstrate the importance of model choice when estimating sojourn times and reveal their potential for understanding healthy and diseased brain mechanisms.

Identification of Temporal Transition of Functional States Using Recurrent Neural Networks from Functional MRI

Dynamic functional connectivity analysis provides valuable information for understanding brain functional activity underlying different cognitive processes. Besides sliding window based approaches, a variety of methods have been developed to automatically split the entire functional MRI scan into segments by detecting change points of functional signals to facilitate better characterization of temporally dynamic functional connectivity patterns. However, these methods are based on certain assumptions for the functional signals, such as Gaussian distribution, which are not necessarily suitable for the fMRI data. In this study, we develop a deep learning based framework for adaptively detecting temporally dynamic functional state transitions in a data-driven way without any explicit modeling assumptions, by leveraging recent advances in recurrent neural networks (RNNs) for sequence modeling. Particularly, we solve this problem in an anomaly detection framework with an assumption that the functional profile of one single time point could be reliably predicted based on its preceding profiles within a stable functional state, while large prediction errors would occur around change points of functional states. We evaluate the proposed method using both task and resting-state fMRI data obtained from the human connectome project and experimental results have demonstrated that the proposed change point detection method could effectively identify change points between different task events and split the resting-state fMRI into segments with distinct functional connectivity patterns.

Large-scale dynamic modeling of task-fMRI signals via subspace system identification

OBJECTIVE

We analyze task-based fMRI time series to produce large-scale dynamical models that are capable of approximating the observed signal with good accuracy.

APPROACH

We extend subspace system identification methods for deterministic and stochastic state-space models with external inputs. The dynamic behavior of the generated models is characterized using control-theoretic analysis tools. To validate their effectiveness, we perform a probabilistic inversion of the identified input-output relationships via joint state-input maximum likelihood estimation. Our experimental setup explores a large dataset generated using state-of-the-art acquisition and pre-processing methods from the Human Connectome Project.

MAIN RESULTS

We analyze both anatomically parcellated and spatially dense time series, and propose an efficient algorithm to address the high-dimensional optimization problem resulting from the latter. Our results enable the quantification of input-output transfer functions between each task condition and each region of the cortex, as exemplified by a motor task. Further, the identified models produce impulse response functions between task conditions and cortical regions that are compatible with typical hemodynamic response functions. We then extend subspace methods to account for multi-subject experimental configurations, identifying models that capture common dynamical characteristics across subjects. Finally, we show that system inversion via maximum-likelihood allows the time-of-occurrence of the task stimuli to be estimated from the observed outputs.

SIGNIFICANCE

The ability to produce dynamical input-output models might have an impact in the expanding field of neurofeedback. In particular, the models we produce allow the partial quantification of the effect of external task-related inputs on the metabolic response of the brain, conditioned on its current state. Such a notion provides a basis for leveraging control-theoretic approaches to neuromodulation and self-regulation in therapeutic applications.

A Bayesian Approach for Estimating Dynamic Functional Network Connectivity in fMRI Data

Dynamic functional connectivity, i.e., the study of how interactions among brain regions change dynamically over the course of an fMRI experiment, has recently received wide interest in the neuroimaging literature. Current approaches for studying dynamic connectivity often rely on ad-hoc approaches for inference, with the fMRI time courses segmented by a sequence of sliding windows. We propose a principled Bayesian approach to dynamic functional connectivity, which is based on the estimation of time varying networks. Our method utilizes a hidden Markov model for classification of latent cognitive states, achieving estimation of the networks in an integrated framework that borrows strength over the entire time course of the experiment. Furthermore, we assume that the graph structures, which define the connectivity states at each time point, are related within a super-graph, to encourage the selection of the same edges among related graphs. We apply our method to simulated task-based fMRI data, where we show how our approach allows the decoupling of the task-related activations and the functional connectivity states. We also analyze data from an fMRI sensorimotor task experiment on an individual healthy subject and obtain results that support the role of particular anatomical regions in modulating interaction between executive control and attention networks.

Sparse Estimation of Resting-State Effective Connectivity From fMRI Cross-Spectra

In functional magnetic resonance imaging (fMRI), functional connectivity is conventionally characterized by correlations between fMRI time series, which are intrinsically undirected measures of connectivity. Yet, some information about the directionality of network connections can nevertheless be extracted from the matrix of pairwise temporal correlations between all considered time series, when expressed in the frequency-domain as a cross-spectral density matrix. Using a sparsity prior, it then becomes possible to determine a unique directed network topology that best explains the observed undirected correlations, without having to rely on temporal precedence relationships that may not be valid in fMRI. Applying this method on simulated data with 100 nodes yielded excellent retrieval of the underlying directed networks under a wide variety of conditions. Importantly, the method did not depend on temporal precedence to establish directionality, thus reducing susceptibility to hemodynamic variability. The computational efficiency of the algorithm was sufficient to enable whole-brain estimations, thus circumventing the problem of missing nodes that otherwise occurs in partial-brain analyses. Applying the method to real resting-state fMRI data acquired with a high temporal resolution, the inferred networks showed good consistency with structural connectivity obtained from diffusion tractography in the same subjects. Interestingly, this agreement could also be seen when considering high-frequency rather than low-frequency connectivity (average correlation: r = 0.26 for f < 0.3 Hz, r = 0.43 for 0.3 < f < 5 Hz). Moreover, this concordance was significantly better (p < 0.05) than for networks obtained with conventional functional connectivity based on correlations (average correlation r = 0.18). The presented methodology thus appears to be well-suited for fMRI, particularly given its lack of explicit dependence on temporal lag structure, and is readily applicable to whole-brain effective connectivity estimation.

Estimating Dynamic Connectivity States in fMRI Using Regime-Switching Factor Models

We consider the challenges in estimating the state-related changes in brain connectivity networks with a large number of nodes. Existing studies use the sliding-window analysis or time-varying coefficient models, which are unable to capture both smooth and abrupt changes simultaneously, and rely on ad-hoc approaches to the high-dimensional estimation. To overcome these limitations, we propose a Markov-switching dynamic factor model, which allows the dynamic connectivity states in functional magnetic resonance imaging (fMRI) data to be driven by lower-dimensional latent factors. We specify a regime-switching vector autoregressive (SVAR) factor process to quantity the time-varying directed connectivity. The model enables a reliable, data-adaptive estimation of change-points of connectivity regimes and the massive dependencies associated with each regime. We develop a three-step estimation procedure: 1) extracting the factors using principal component analysis, 2) identifying connectivity regimes in a low-dimensional subspace based on the factor-based SVAR model, and 3) constructing high-dimensional state connectivity metrics based on the subspace estimates. Simulation results show that our estimator outperforms -means clustering of time-windowed coefficients, providing more accurate estimate of time-evolving connectivity. It achieves percentage of reduction in mean squared error by 60% when the network dimension is comparable to the sample size. When applied to the resting-state fMRI data, our method successfully identifies modular organization in the resting-statenetworksin consistencywith other studies. It further reveals changes in brain states with variations across subjects and distinct large-scale directed connectivity patterns across states.

Bayesian switching factor analysis for estimating time-varying functional connectivity in fMRI

There is growing interest in understanding the dynamical properties of functional interactions between distributed brain regions. However, robust estimation of temporal dynamics from functional magnetic resonance imaging (fMRI) data remains challenging due to limitations in extant multivariate methods for modeling time-varying functional interactions between multiple brain areas. Here, we develop a Bayesian generative model for fMRI time-series within the framework of hidden Markov models (HMMs). The model is a dynamic variant of the static factor analysis model (Ghahramani and Beal, 2000). We refer to this model as Bayesian switching factor analysis (BSFA) as it integrates factor analysis into a generative HMM in a unified Bayesian framework. In BSFA, brain dynamic functional networks are represented by latent states which are learnt from the data. Crucially, BSFA is a generative model which estimates the temporal evolution of brain states and transition probabilities between states as a function of time. An attractive feature of BSFA is the automatic determination of the number of latent states via Bayesian model selection arising from penalization of excessively complex models. Key features of BSFA are validated using extensive simulations on carefully designed synthetic data. We further validate BSFA using fingerprint analysis of multisession resting-state fMRI data from the Human Connectome Project (HCP). Our results show that modeling temporal dependencies in the generative model of BSFA results in improved fingerprinting of individual participants. Finally, we apply BSFA to elucidate the dynamic functional organization of the salience, central-executive, and default mode networks-three core neurocognitive systems with central role in cognitive and affective information processing (Menon, 2011). Across two HCP sessions, we demonstrate a high level of dynamic interactions between these networks and determine that the salience network has the highest temporal flexibility among the three networks. Our proposed methods provide a novel and powerful generative model for investigating dynamic brain connectivity.

Bayesian vector autoregressive model for multi-subject effective connectivity inference using multi-modal neuroimaging data

In this article a multi-subject vector autoregressive (VAR) modeling approach was proposed for inference on effective connectivity based on resting-state functional MRI data. Their framework uses a Bayesian variable selection approach to allow for simultaneous inference on effective connectivity at both the subject- and group-level. Furthermore, it accounts for multi-modal data by integrating structural imaging information into the prior model, encouraging effective connectivity between structurally connected regions. They demonstrated through simulation studies that their approach resulted in improved inference on effective connectivity at both the subject- and group-level, compared with currently used methods. It was concluded by illustrating the method on temporal lobe epilepsy data, where resting-state functional MRI and structural MRI were used. Hum Brain Mapp 38:1311-1332, 2017. © 2016 Wiley Periodicals, Inc.