Disentangling causal webs in the brain using functional Magnetic Resonance Imaging: A review of current approaches

In the past two decades, functional Magnetic Resonance Imaging has been used to relate neuronal network activity to cognitive processing and behaviour. Recently this approach has been augmented by algorithms that allow us to infer causal links between component populations of neuronal networks. Multiple inference procedures have been proposed to approach this research question but so far, each method has limitations when it comes to establishing whole-brain connectivity patterns. In this work, we discuss eight ways to infer causality in fMRI research: Bayesian Nets, Dynamical Causal Modelling, Granger Causality, Likelihood Ratios, LiNGAM, Patel's Tau, Structural Equation Modelling, and Transfer Entropy. We finish with formulating some recommendations for the future directions in this area

Test-retest reliability of regression dynamic causal modeling

Regression dynamic causal modeling (rDCM) is a novel and computationally highly efficient method for inferring effective connectivity at the whole-brain level. While face and construct validity of rDCM have already been demonstrated, here we assessed its test-retest reliability-a test-theoretical property of particular importance for clinical applications-together with group-level consistency of connection-specific estimates and consistency of whole-brain connectivity patterns over sessions. Using the Human Connectome Project dataset for eight different paradigms (tasks and rest) and two different parcellation schemes, we found that rDCM provided highly consistent connectivity estimates at the group level across sessions. Second, while test-retest reliability was limited when averaging over all connections (range of mean intraclass correlation coefficient 0.24-0.42 over tasks), reliability increased with connection strength, with stronger connections showing good to excellent test-retest reliability. Third, whole-brain connectivity patterns by rDCM allowed for identifying individual participants with high (and in some cases perfect) accuracy. Comparing the test-retest reliability of rDCM connectivity estimates with measures of functional connectivity, rDCM performed favorably-particularly when focusing on strong connections. Generally, for all methods and metrics, task-based connectivity estimates showed greater reliability than those from the resting state. Our results underscore the potential of rDCM for human connectomics and clinical applications. Keywords: Connectomics; Effective connectivity; Generative model; Regression dynamic causal modeling; Test-retest reliability; rDC

Test-retest reliability of the magnetic mismatch negativity response to sound duration and omission deviants

Mismatch negativity (MMN) is a neurophysiological measure of auditory novelty detection that could serve as a translational biomarker of psychiatric disorders, such as schizophrenia. However, the replicability of its magnetoencephalographic (MEG) counterpart (MMNm) has been insufficiently addressed. In the current study, test-retest reliability of the MMNm response to both duration and omission deviants was evaluated over two MEG sessions in 16 healthy adults. MMNm amplitudes and latencies were obtained at both sensor- and source-level using a cortically-constrained minimum-norm approach. Intraclass correlations (ICC) were derived to assess stability of MEG responses over time. In addition, signal-to-noise ratios (SNR) and within-subject statistics were obtained in order to determine MMNm detectability in individual participants. ICC revealed robust values at both sensor- and source-level for both duration and omission MMNm amplitudes (ICC = 0.81-0.90), in particular in the right hemisphere, while moderate to strong values were obtained for duration MMNm and omission MMNm peak latencies (ICC = 0.74-0.88). Duration MMNm was robustly identified in individual participants with high SNR, whereas omission MMNm responses were only observed in half of the participants. Our data indicate that MMNm to unexpected duration changes and omitted sounds are highly reproducible, providing support for the use of MEG-parameters in basic and clinical research

Model specification and the reliability of fMRI results: Implications for longitudinal neuroimaging studies in psychiatry

Functional Magnetic Resonance Imagine (fMRI) is an important assessment tool in longitudinal studies of mental illness and its treatment. Understanding the psychometric properties of fMRI-based metrics, and the factors that influence them, will be critical for properly interpreting the results of these efforts. The current study examined whether the choice among alternative model specifications affects estimates of test-retest reliability in key emotion processing regions across a 6-month interval. Subjects (N = 46) performed an emotional-faces paradigm during fMRI in which neutral faces dynamically morphed into one of four emotional faces. Median voxelwise intraclass correlation coefficients (mvICCs) were calculated to examine stability over time in regions showing task-related activity as well as in bilateral amygdala. Four modeling choices were evaluated: a default model that used the canonical hemodynamic response function (HRF), a flexible HRF model that included additional basis functions, a modified CompCor (mCompCor) model that added corrections for physiological noise in the global signal, and a final model that combined the flexible HRF and mCompCor models. Model residuals were examined to determine the degree to which each pipeline met modeling assumptions. Results indicated that the choice of modeling approaches impacts both the degree to which model assumptions are met and estimates of test-retest reliability. ICC estimates in the visual cortex increased from poor (mvICC = 0.31) in the default pipeline to fair (mvICC = 0.45) in the full alternative pipeline - an increase of 45%. In nearly all tests, the models with the fewest assumption violations generated the highest ICC estimates. Implications for longitudinal treatment studies that utilize fMRI are discussed. © 2014 Fournier et al

Learning and comparing functional connectomes across subjects

Functional connectomes capture brain interactions via synchronized fluctuations in the functional magnetic resonance imaging signal. If measured during rest, they map the intrinsic functional architecture of the brain. With task-driven experiments they represent integration mechanisms between specialized brain areas. Analyzing their variability across subjects and conditions can reveal markers of brain pathologies and mechanisms underlying cognition. Methods of estimating functional connectomes from the imaging signal have undergone rapid developments and the literature is full of diverse strategies for comparing them. This review aims to clarify links across functional-connectivity methods as well as to expose different steps to perform a group study of functional connectomes

Bayesian population receptive field modelling

We introduce a probabilistic (Bayesian) framework and associated software toolbox for mapping population receptive fields (pRFs) based on fMRI data. This generic approach is intended to work with stimuli of any dimension and is demonstrated and validated in the context of 2D retinotopic mapping. The framework enables the experimenter to specify generative (encoding) models of fMRI timeseries, in which experimental manipulations enter a pRF model of neural activity, which in turns drives a nonlinear model of neurovascular coupling and Blood Oxygenation Level Dependent (BOLD) response. The neuronal and haemodynamic parameters are estimated together on a voxel-by-voxel or region-of-interest basis using a Bayesian estimation algorithm (variational Laplace). This offers several novel contributions to receptive field modelling. The variance / covariance of parameters are estimated, enabling receptive fields to be plotted while properly representing uncertainty about pRF size and location. Variability in the haemodynamic response across the brain is accounted for. Furthermore, the framework introduces formal hypothesis testing to pRF analysis, enabling competing models to be evaluated based on their model evidence (approximated by the variational free energy), which represents the optimal tradeoff between accuracy and complexity. Using simulations and empirical data, we found that parameters typically used to represent pRF size and neuronal scaling are strongly correlated, which should be taken into account when making inferences. We used the framework to compare the evidence for six variants of pRF model using 7T functional MRI data and we found a circular Difference of Gaussians (DoG) model to be the best explanation for our data overall. We hope this framework will prove useful for mapping stimulus spaces with any number of dimensions onto the anatomy of the brain.Comment: 30 pages, 10 figures. Code available at https://github.com/pzeidman/BayespR

Using resting-state DMN effective connectivity to characterize the neurofunctional architecture of empathy

Neuroimaging studies in social neuroscience have largely relied on functional connectivity (FC) methods to characterize the functional integration between different brain regions. However, these methods have limited utility in social-cognitive studies that aim to understand the directed information flow among brain areas that underlies complex psychological processes. In this study we combined functional and effective connectivity approaches to characterize the functional integration within the Default Mode Network (DMN) and its role in self-perceived empathy. Forty-two participants underwent a resting state fMRI scan and completed a questionnaire of dyadic empathy. Independent Component Analysis (ICA) showed that higher empathy scores were associated with an increased contribution of the medial prefrontal cortex (mPFC) to the DMN spatial mode. Dynamic causal modelling (DCM) combined with Canonical Variance Analysis (CVA) revealed that this association was mediated indirectly by the posterior cingulate cortex (PCC) via the right inferior parietal lobule (IPL). More specifically, in participants with higher scores in empathy, the PCC had a greater effect on bilateral IPL and the right IPL had a greater influence on mPFC. These results highlight the importance of using analytic approaches that address directed and hierarchical connectivity within networks, when studying complex psychological phenomena, such as empathy.- This study was funded by BIAL Foundation (Grant number 87/12); by the Portuguese Foundation for Science and Technology and the Portuguese Ministry of Education and Science through national funds and co-financed by FEDER through COMPETE2020 under the PT2020 Partnership Agreement (POCI-01-0145-FEDER-007653); by the postdoctoral scholarship UMINHO/BPD/18/2017 and by the Portuguese Foundation for Science Doctoral scholarship (PD/BD/105963/2014). This work was conducted at Psychology Research Centre (UID/PSI/01662/2013), University of Minho

An information-theoretic analysis of resting-state versus task fMRI

Resting-state fMRI is an increasingly popular alternative to task-based fMRI. However, a formal quantification of the amount of information provided by resting-state fMRI as opposed to active task conditions about neural responses is lacking. We conducted a systematic comparison of the quality of inferences derived from a resting-state and a task fMRI paradigm by means of Bayesian Data Comparison. In this framework, data quality is formally quantified in information-theoretic terms as the precision and amount of information provided by the data on the parameters of interest. Parameters of effective connectivity, estimated from the cross-spectral densities of resting-state- and task time series by means of dynamic causal modelling (DCM), were subjected to the analysis. Data from 50 individuals undergoing resting-state and a Theory-of-Mind task were compared, both datasets provided by the Human Connectome Project. A threshold of very strong evidence was reached in favour of the Theory-of-Mind task (>10 bits or natural units) regarding information gain, which could be attributed to the active task condition eliciting stronger effective connectivity. Extending these analyses to other tasks and cognitive systems will reveal whether the superior informative value of task-based fMRI observed here is case specific or a more general trend

Reliability of Graph Measures Derived from Resting-State MEG Data Using Source Space Functional Connectivity Analysis

The reliability of global graph measures derived from neuroimaging data is an important criterion for their use as biomarkers for neurological disorders. This study examined the reliability of the global efficiency (GE), characteristic path length (CPL), transitivity, and synchronizability of functional whole-brain and intra-hemispheric networks based on resting-state magnetoencephalography. Brain sources were reconstructed using atlas-based beamforming, and functional connectivity in six frequency bands was estimated using the debiased weighted phase lag index. An optimal threshold of 100% was chosen based on test-retest reliability of the measures. At this threshold, test-retest reliability of the GE, CPL, and transitivity was mostly fair to excellent except for in the delta band. However, test-retest reliability of the synchronizability was mostly poor to fair. There was no significant effect of gender on any graph measure. Overall, these results indicate that the GE, CPL, and transitivity in most of the frequency bands may be useful biomarkers

Bayesian Dynamic DAG Learning: Application in Discovering Dynamic Effective Connectome of Brain

Understanding the complex mechanisms of the brain can be unraveled by extracting the Dynamic Effective Connectome (DEC). Recently, score-based Directed Acyclic Graph (DAG) discovery methods have shown significant improvements in extracting the causal structure and inferring effective connectivity. However, learning DEC through these methods still faces two main challenges: one with the fundamental impotence of high-dimensional dynamic DAG discovery methods and the other with the low quality of fMRI data. In this paper, we introduce Bayesian Dynamic DAG learning with M-matrices Acyclicity characterization \textbf{(BDyMA)} method to address the challenges in discovering DEC. The presented dynamic causal model enables us to discover bidirected edges as well. Leveraging an unconstrained framework in the BDyMA method leads to more accurate results in detecting high-dimensional networks, achieving sparser outcomes, making it particularly suitable for extracting DEC. Additionally, the score function of the BDyMA method allows the incorporation of prior knowledge into the process of dynamic causal discovery which further enhances the accuracy of results. Comprehensive simulations on synthetic data and experiments on Human Connectome Project (HCP) data demonstrate that our method can handle both of the two main challenges, yielding more accurate and reliable DEC compared to state-of-the-art and baseline methods. Additionally, we investigate the trustworthiness of DTI data as prior knowledge for DEC discovery and show the improvements in DEC discovery when the DTI data is incorporated into the process