Pair-wise Clustering of Large Scale Granger Causality Index Matrices for Revealing Communities

Automated diagnosis of HIV-associated neurocognitive disorders using large-scale Granger causality analysis of resting-state functional MRI

HIV-associated neurocognitive disorders (HAND) represent an important source of neurologic complications in individuals with HIV. The dynamic, often subclinical, course of HAND has rendered diagnosis, which currently depends on neuropsychometric (NP) evaluation, a challenge for clinicians. Here, we present evidence that functional brain connectivity, derived by large-scale Granger causality (lsGC) analysis of resting-state functional MRI (rs-fMRI) time-series, represents a potential biomarker to address this critical diagnostic need. Brain graph properties were used as features in machine learning tasks to 1) classify individuals as HIV⁺ or HIV^- and 2) to predict overall cognitive performance, as assessed by NP scores, in a 22-subject (13 HIV^-, 9 HIV⁺) cohort. Over nearly all seven brain parcellation templates considered, support vector machine (SVM) classifiers based on lsGC-derived brain graph features significantly outperformed those based on conventional Pearson correlation (PC)-derived features (p<0.05, Bonferroni-corrected). In a second task for which the objective was to predict the overall NP score of each subject, the lsGC-based SVM regressors consistently outperformed the PC-based regressors (p<0.05, Bonferroni-corrected) on nearly all templates. With the widely used Automated Anatomical Labeling (AAL90) template, it was determined that the brain regions that figured most strongly in the SVM classifiers included those of the default mode network (posterior cingulate cortex, angular gyrus) and basal ganglia (caudate nucleus), dysfunction in both of which have been observed in previous structural and functional analyses of HAND.

Multivariate Granger causality unveils directed parietal to prefrontal cortex connectivity during task-free MRI

While a large body of research has focused on the study of functional brain "connectivity", few investigators have focused on directionality of brain-brain interactions which, in spite of the mostly bidirectional anatomical substrates, cannot be assumed to be symmetrical. We employ a multivariate Granger Causality-based approach to estimating directed in-network interactions and quantify its advantages using extensive realistic synthetic BOLD data simulations to match Human Connectome Project (HCP) data specification. We then apply our framework to resting state functional MRI (rs-fMRI) data provided by the HCP to estimate the directed connectome of the human brain. We show that the functional interactions between parietal and prefrontal cortices commonly observed in rs-fMRI studies are not symmetrical, but consists of directional connectivity from parietal areas to prefrontal cortices rather than vice versa. These effects are localized within the same hemisphere and do not generalize to cross-hemispheric functional interactions. Our data are consistent with neurophysiological evidence that posterior parietal cortices involved in processing and integration of multi-sensory information modulate the function of more anterior prefrontal regions implicated in action control and goal-directed behaviour. The directionality of functional connectivity can provide an additional layer of information in interpreting rs-fMRI studies both in health and disease.

Exploring connectivity with large-scale Granger causality on resting-state functional MRI

BACKGROUND

Large-scale Granger causality (lsGC) is a recently developed, resting-state functional MRI (fMRI) connectivity analysis approach that estimates multivariate voxel-resolution connectivity. Unlike most commonly used multivariate approaches, which establish coarse-resolution connectivity by aggregating voxel time-series avoiding an underdetermined problem, lsGC estimates voxel-resolution, fine-grained connectivity by incorporating an embedded dimension reduction.

NEW METHOD

We investigate application of lsGC on realistic fMRI simulations, modeling smoothing of neuronal activity by the hemodynamic response function and repetition time (TR), and empirical resting-state fMRI data. Subsequently, functional subnetworks are extracted from lsGC connectivity measures for both datasets and validated quantitatively. We also provide guidelines to select lsGC free parameters.

RESULTS

Results indicate that lsGC reliably recovers underlying network structure with area under receiver operator characteristic curve (AUC) of 0.93 at TR=1.5s for a 10-min session of fMRI simulations. Furthermore, subnetworks of closely interacting modules are recovered from the aforementioned lsGC networks. Results on empirical resting-state fMRI data demonstrate recovery of visual and motor cortex in close agreement with spatial maps obtained from (i) visuo-motor fMRI stimulation task-sequence (Accuracy=0.76) and (ii) independent component analysis (ICA) of resting-state fMRI (Accuracy=0.86).

COMPARISON WITH EXISTING METHOD(S)

Compared with conventional Granger causality approach (AUC=0.75), lsGC produces better network recovery on fMRI simulations. Furthermore, it cannot recover functional subnetworks from empirical fMRI data, since quantifying voxel-resolution connectivity is not possible as consequence of encountering an underdetermined problem.

CONCLUSIONS

Functional network recovery from fMRI data suggests that lsGC gives useful insight into connectivity patterns from resting-state fMRI at a multivariate voxel-resolution.

Mutual Connectivity Analysis (MCA) Using Generalized Radial Basis Function Neural Networks for Nonlinear Functional Connectivity Network Recovery in Resting-State Functional MRI

We investigate the applicability of a computational framework, called mutual connectivity analysis (MCA), for directed functional connectivity analysis in both synthetic and resting-state functional MRI data. This framework comprises of first evaluating non-linear cross-predictability between every pair of time series prior to recovering the underlying network structure using community detection algorithms. We obtain the non-linear cross-prediction score between time series using Generalized Radial Basis Functions (GRBF) neural networks. These cross-prediction scores characterize the underlying functionally connected networks within the resting brain, which can be extracted using non-metric clustering approaches, such as the Louvain method. We first test our approach on synthetic models with known directional influence and network structure. Our method is able to capture the directional relationships between time series (with an area under the ROC curve = 0.92 ± 0.037) as well as the underlying network structure (Rand index = 0.87 ± 0.063) with high accuracy. Furthermore, we test this method for network recovery on resting-state fMRI data, where results are compared to the motor cortex network recovered from a motor stimulation sequence, resulting in a strong agreement between the two (Dice coefficient = 0.45). We conclude that our MCA approach is effective in analyzing non-linear directed functional connectivity and in revealing underlying functional network structure in complex systems.

Large-Scale Granger Causality Analysis on Resting-State Functional MRI

We demonstrate an approach to measure the information flow between each pair of time series in resting-state functional MRI (fMRI) data of the human brain and subsequently recover its underlying network structure. By integrating dimensionality reduction into predictive time series modeling, large-scale Granger Causality (lsGC) analysis method can reveal directed information flow suggestive of causal influence at an individual voxel level, unlike other multivariate approaches. This method quantifies the influence each voxel time series has on every other voxel time series in a multivariate sense and hence contains information about the underlying dynamics of the whole system, which can be used to reveal functionally connected networks within the brain. To identify such networks, we perform non-metric network clustering, such as accomplished by the Louvain method. We demonstrate the effectiveness of our approach to recover the motor and visual cortex from resting state human brain fMRI data and compare it with the network recovered from a visuomotor stimulation experiment, where the similarity is measured by the Dice Coefficient (DC). The best DC obtained was 0.59 implying a strong agreement between the two networks. In addition, we thoroughly study the effect of dimensionality reduction in lsGC analysis on network recovery. We conclude that our approach is capable of detecting causal influence between time series in a multivariate sense, which can be used to segment functionally connected networks in the resting-state fMRI.