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

Identifying effective connectivity parameters in simulated fMRI: a direct comparison of switching linear dynamic system, stochastic dynamic causal, and multivariate autoregressive models

abstract: The number and variety of connectivity estimation methods is likely to continue to grow over the coming decade. Comparisons between methods are necessary to prune this growth to only the most accurate and robust methods. However, the nature of connectivity is elusive with different methods potentially attempting to identify different aspects of connectivity. Commonalities of connectivity definitions across methods upon which base direct comparisons can be difficult to derive. Here, we explicitly define “effective connectivity” using a common set of observation and state equations that are appropriate for three connectivity methods: dynamic causal modeling (DCM), multivariate autoregressive modeling (MAR), and switching linear dynamic systems for fMRI (sLDSf). In addition while deriving this set, we show how many other popular functional and effective connectivity methods are actually simplifications of these equations. We discuss implications of these connections for the practice of using one method to simulate data for another method. After mathematically connecting the three effective connectivity methods, simulated fMRI data with varying numbers of regions and task conditions is generated from the common equation. This simulated data explicitly contains the type of the connectivity that the three models were intended to identify. Each method is applied to the simulated data sets and the accuracy of parameter identification is analyzed. All methods perform above chance levels at identifying correct connectivity parameters. The sLDSf method was superior in parameter estimation accuracy to both DCM and MAR for all types of comparisons.View the article as published at http://journal.frontiersin.org/article/10.3389/fnins.2013.00070/ful

Neuroimaging of structural pathology and connectomics in traumatic brain injury: Toward personalized outcome prediction.

Recent contributions to the body of knowledge on traumatic brain injury (TBI) favor the view that multimodal neuroimaging using structural and functional magnetic resonance imaging (MRI and fMRI, respectively) as well as diffusion tensor imaging (DTI) has excellent potential to identify novel biomarkers and predictors of TBI outcome. This is particularly the case when such methods are appropriately combined with volumetric/morphometric analysis of brain structures and with the exploration of TBI-related changes in brain network properties at the level of the connectome. In this context, our present review summarizes recent developments on the roles of these two techniques in the search for novel structural neuroimaging biomarkers that have TBI outcome prognostication value. The themes being explored cover notable trends in this area of research, including (1) the role of advanced MRI processing methods in the analysis of structural pathology, (2) the use of brain connectomics and network analysis to identify outcome biomarkers, and (3) the application of multivariate statistics to predict outcome using neuroimaging metrics. The goal of the review is to draw the community's attention to these recent advances on TBI outcome prediction methods and to encourage the development of new methodologies whereby structural neuroimaging can be used to identify biomarkers of TBI outcome

Connectivity Analysis in EEG Data: A Tutorial Review of the State of the Art and Emerging Trends

Understanding how different areas of the human brain communicate with each other is a crucial issue in neuroscience. The concepts of structural, functional and effective connectivity have been widely exploited to describe the human connectome, consisting of brain networks, their structural connections and functional interactions. Despite high-spatial-resolution imaging techniques such as functional magnetic resonance imaging (fMRI) being widely used to map this complex network of multiple interactions, electroencephalographic (EEG) recordings claim high temporal resolution and are thus perfectly suitable to describe either spatially distributed and temporally dynamic patterns of neural activation and connectivity. In this work, we provide a technical account and a categorization of the most-used data-driven approaches to assess brain-functional connectivity, intended as the study of the statistical dependencies between the recorded EEG signals. Different pairwise and multivariate, as well as directed and non-directed connectivity metrics are discussed with a pros-cons approach, in the time, frequency, and information-theoretic domains. The establishment of conceptual and mathematical relationships between metrics from these three frameworks, and the discussion of novel methodological approaches, will allow the reader to go deep into the problem of inferring functional connectivity in complex networks. Furthermore, emerging trends for the description of extended forms of connectivity (e.g., high-order interactions) are also discussed, along with graph-theory tools exploring the topological properties of the network of connections provided by the proposed metrics. Applications to EEG data are reviewed. In addition, the importance of source localization, and the impacts of signal acquisition and pre-processing techniques (e.g., filtering, source localization, and artifact rejection) on the connectivity estimates are recognized and discussed. By going through this review, the reader could delve deeply into the entire process of EEG pre-processing and analysis for the study of brain functional connectivity and learning, thereby exploiting novel methodologies and approaches to the problem of inferring connectivity within complex networks

Regression DCM for fMRI

The development of large-scale network models that infer the effective (directed) connectivity among neuronal populations from neuroimaging data represents a key challenge for computational neuroscience. Dynamic causal models (DCMs) of neuroimaging and electrophysiological data are frequently used for inferring effective connectivity but are presently restricted to small graphs (typically up to 10 regions) in order to keep model inversion computationally feasible. Here, we present a novel variant of DCM for functional magnetic resonance imaging (fMRI) data that is suited to assess effective connectivity in large (whole-brain) networks. The approach rests on translating a linear DCM into the frequency domain and reformulating it as a special case of Bayesian linear regression. This paper derives regression DCM (rDCM) in detail and presents a variational Bayesian inversion method that enables extremely fast inference and accelerates model inversion by several orders of magnitude compared to classical DCM. Using both simulated and empirical data, we demonstrate the face validity of rDCM under different settings of signal-to-noise ratio (SNR) and repetition time (TR) of fMRI data. In particular, we assess the potential utility of rDCM as a tool for whole-brain connectomics by challenging it to infer effective connection strengths in a simulated whole-brain network comprising 66 regions and 300 free parameters. Our results indicate that rDCM represents a computationally highly efficient approach with promising potential for inferring whole-brain connectivity from individual fMRI data

Parasympathetic-sympathetic causal interactions assessed by time-varying multivariate autoregressive modeling of electrodermal activity and heart-rate-variability

Objective: Most of the bodily functions are regulated by multiple interactions between the parasympathetic (PNS) and sympathetic (SNS) nervous system. In this study, we propose a novel framework to quantify the causal flow of information between PNS and SNS through the analysis of heart rate variability (HRV) and electrodermal activity (EDA) signals. Methods: Our method is based on a time-varying (TV) multivariate autoregressive model of EDA and HRV time-series and incorporates physiologically inspired assumptions by estimating the Directed Coherence in a specific frequency range. The statistical significance of the observed interactions is assessed by a bootstrap procedure purposely developed to infer causalities in the presence of both TV model coefficients and TV model residuals (i.e., heteroskedasticity). We tested our method on two different experiments designed to trigger a sympathetic response, i.e., a hand-grip task (HG) and a mental-computation task (MC). Results: Our results show a parasympathetic driven interaction in the resting state, which is consistent across different studies. The onset of the stressful stimulation triggers a cascade of events characterized by the presence or absence of the PNS-SNS interaction and changes in the directionality. Despite similarities between the results related to the two tasks, we reveal differences in the dynamics of the PNS-SNS interaction, which might reflect different regulatory mechanisms associated with different stressors. Conclusion: We estimate causal coupling between PNS and SNS through MVAR modeling of EDA and HRV time-series. Significance: Our results suggest promising future applicability to investigate more complex contexts such as affective and pathological scenarios

Doctor of Philosophy

dissertationMagnetic Resonance (MR) is a relatively risk-free and flexible imaging modality that is widely used for studying the brain. Biophysical and chemical properties of brain tissue are captured by intensity measurements in T1W (T1-Weighted) and T2W (T2-Weighted) MR scans. Rapid maturational processes taking place in the infant brain manifest as changes in co{\tiny }ntrast between white matter and gray matter tissue classes in these scans. However, studies based on MR image appearance face severe limitations due to the uncalibrated nature of MR intensity and its variability with respect to changing conditions of scan. In this work, we develop a method for studying the intensity variations between brain white matter and gray matter that are observed during infant brain development. This method is referred to by the acronym WIVID (White-gray Intensity Variation in Infant Development). WIVID is computed by measuring the Hellinger Distance of separation between intensity distributions of WM (White Matter) and GM (Gray Matter) tissue classes. The WIVID measure is shown to be relatively stable to interscan variations compared with raw signal intensity and does not require intensity normalization. In addition to quantification of tissue appearance changes using the WIVID measure, we test and implement a statistical framework for modeling temporal changes in this measure. WIVID contrast values are extracted from MR scans belonging to large-scale, longitudinal, infant brain imaging studies and modeled using the NLME (Nonlinear Mixed Effects) method. This framework generates a normative model of WIVID contrast changes with time, which captures brain appearance changes during neurodevelopment. Parameters from the estimated trajectories of WIVID contrast change are analyzed across brain lobes and image modalities. Parameters associated with the normative model of WIVID contrast change reflect established patterns of region-specific and modality-specific maturational sequences. We also detect differences in WIVID contrast change trajectories between distinct population groups. These groups are categorized based on sex and risk/diagnosis for ASD (Autism Spectrum Disorder). As a result of this work, the usage of the proposed WIVID contrast measure as a novel neuroimaging biomarker for characterizing tissue appearance is validated, and the clinical potential of the developed framework is demonstrated

PATH AND DIRECTION DISCOVERY IN INDIVIDUAL DYNAMIC MODELS: A REGULARIZED HYBRID UNIFIED STRUCTURAL EQUATION MODELING WITH LATENT VARIABLE

There recently has been growing interest in the study of psychological and neurologicalprocesses at an individual level. One goal in such endeavors is to construct person-specific dynamicassessments using time series techniques such as Vector Autoregressive (VAR) models. Within thepsychometric field, researchers have developed psychometric modeling frameworks to estimatedifferent variants of VAR models. These modeling frameworks estimate the dynamic relations (e.g.,temporal and contemporaneous) unpacked in a multivariate time series data. However, two problemsexist with current VAR specifications: 1) VAR models are restricted in that contemporaneousrelations are typically modeled either as undirected relations among residuals or directed relationsamong observed variables, but not both; 2) current estimation frameworks are limited by thereliance on stepwise model building procedures. This study adopts a new modeling approach, i.e.,LASSO regularized hybrid unified Structural Equation Model (SEM), for a global search andestimation of a more flexible VAR representation. The present study to our knowledge is the firstapplication of the recently developed regularized SEM technique to the estimation of a type of timeseries SEM, which points to a promising future for statistical learning in psychometric models.Doctor of Philosoph