Identifying Neural Drivers with Functional MRI: An Electrophysiological Validation

Whether functional magnetic resonance imaging (fMRI) allows the identification of neural drivers remains an open question of particular importance to refine physiological and neuropsychological models of the brain, and/or to understand neurophysiopathology. Here, in a rat model of absence epilepsy showing spontaneous spike-and-wave discharges originating from the first somatosensory cortex (S1BF), we performed simultaneous electroencephalographic (EEG) and fMRI measurements, and subsequent intracerebral EEG (iEEG) recordings in regions strongly activated in fMRI (S1BF, thalamus, and striatum). fMRI connectivity was determined from fMRI time series directly and from hidden state variables using a measure of Granger causality and Dynamic Causal Modelling that relates synaptic activity to fMRI. fMRI connectivity was compared to directed functional coupling estimated from iEEG using asymmetry in generalised synchronisation metrics. The neural driver of spike-and-wave discharges was estimated in S1BF from iEEG, and from fMRI only when hemodynamic effects were explicitly removed. Functional connectivity analysis applied directly on fMRI signals failed because hemodynamics varied between regions, rendering temporal precedence irrelevant. This paper provides the first experimental substantiation of the theoretical possibility to improve interregional coupling estimation from hidden neural states of fMRI. As such, it has important implications for future studies on brain connectivity using functional neuroimaging.     

Granger causality analysis of steady-state electroencephalographic signals during propofol-induced anaesthesia

Changes in conscious level have been associated with changes in dynamical integration and segregation among distributed brain regions. Recent theoretical developments emphasize changes in directed functional (i.e., causal) connectivity as reflected in quantities such as 'integrated information' and 'causal density'. Here we develop and illustrate a rigorous methodology for assessing causal connectivity from electroencephalographic (EEG) signals using Granger causality (GC). Our method addresses the challenges of non-stationarity and bias by dividing data into short segments and applying permutation analysis. We apply the method to EEG data obtained from subjects undergoing propofol-induced anaesthesia, with signals source-localized to the anterior and posterior cingulate cortices. We found significant increases in bidirectional GC in most subjects during loss-of-consciousness, especially in the beta and gamma frequency ranges. Corroborating a previous analysis we also found increases in synchrony in these ranges; importantly, the Granger causality analysis showed higher inter-subject consistency than the synchrony analysis. Finally, we validate our method using simulated data generated from a model for which GC values can be analytically derived. In summary, our findings advance the methodology of Granger causality analysis of EEG data and carry implications for integrated information and causal density theories of consciousness.     

Uncovering interactions in the frequency domain

Oscillatory activity plays a critical role in regulating biological processes at levels ranging from subcellular, cellular, and network to the whole organism, and often involves a large number of interacting elements. We shed light on this issue by introducing a novel approach called partial Granger causality to reliably reveal interaction patterns in multivariate data with exogenous inputs and latent variables in the frequency domain. The method is extensively tested with toy models, and successfully applied to experimental datasets, including (1) gene microarray data of HeLa cell cycle; (2) in vivo multi-electrode array (MEA) local field potentials (LFPs) recorded from the inferotemporal cortex of a sheep; and (3) in vivo LFPs recorded from distributed sites in the right hemisphere of a macaque monkey.     

Changes in Effective Connectivity by Propofol Sedation

Mechanisms of propofol-induced loss of consciousness remain poorly understood. Recent fMRI studies have shown decreases in functional connectivity during unconsciousness induced by this anesthetic agent. Functional connectivity does not provide information of directional changes in the dynamics observed during unconsciousness. The aim of the present study was to investigate, in healthy humans during an auditory task, the changes in effective connectivity resulting from propofol induced loss of consciousness. We used Dynamic Causal Modeling for fMRI (fMRI-DCM) to assess how causal connectivity is influenced by the anesthetic agent in the auditory system. Our results suggest that the dynamic observed in the auditory system during unconsciousness induced by propofol, can result in a mixture of two effects: a local inhibitory connectivity increase and a decrease in the effective connectivity in sensory cortices.     

A Multivariate Granger Causality Concept towards Full Brain Functional Connectivity

Detecting changes of spatially high-resolution functional connectivity patterns in the brain is crucial for improving the fundamental understanding of brain function in both health and disease, yet still poses one of the biggest challenges in computational neuroscience. Currently, classical multivariate Granger Causality analyses of directed interactions between single process components in coupled systems are commonly restricted to spatially low- dimensional data, which requires a pre-selection or aggregation of time series as a preprocessing step. In this paper we propose a new fully multivariate Granger Causality approach with embedded dimension reduction that makes it possible to obtain a representation of functional connectivity for spatially high-dimensional data. The resulting functional connectivity networks may consist of several thousand vertices and thus contain more detailed information compared to connectivity networks obtained from approaches based on particular regions of interest. Our large scale Granger Causality approach is applied to synthetic and resting state fMRI data with a focus on how well network community structure, which represents a functional segmentation of the network, is preserved. It is demonstrated that a number of different community detection algorithms, which utilize a variety of algorithmic strategies and exploit topological features differently, reveal meaningful information on the underlying network module structure.     

Comparing Families of Dynamic Causal Models

Mathematical models of scientific data can be formally compared using Bayesian model evidence. Previous applications in the biological sciences have mainly focussed on model selection in which one first selects the model with the highest evidence and then makes inferences based on the parameters of that model. This “best model” approach is very useful but can become brittle if there are a large number of models to compare, and if different subjects use different models. To overcome this shortcoming we propose the combination of two further approaches: (i) family level inference and (ii) Bayesian model averaging within families. Family level inference removes uncertainty about aspects of model structure other than the characteristic of interest. For example: What are the inputs to the system? Is processing serial or parallel? Is it linear or nonlinear? Is it mediated by a single, crucial connection? We apply Bayesian model averaging within families to provide inferences about parameters that are independent of further assumptions about model structure. We illustrate the methods using Dynamic Causal Models of brain imaging data.     

Transfer entropy—a model-free measure of effective connectivity for the neurosciences

Understanding causal relationships, or effective connectivity, between parts of the brain is of utmost importance because a large part of the brain’s activity is thought to be internally generated and, hence, quantifying stimulus response relationships alone does not fully describe brain dynamics. Past efforts to determine effective connectivity mostly relied on model based approaches such as Granger causality or dynamic causal modeling. Transfer entropy (TE) is an alternative measure of effective connectivity based on information theory. TE does not require a model of the interaction and is inherently non-linear. We investigated the applicability of TE as a metric in a test for effective connectivity to electrophysiological data based on simulations and magnetoencephalography (MEG) recordings in a simple motor task. In particular, we demonstrate that TE improved the detectability of effective connectivity for non-linear interactions, and for sensor level MEG signals where linear methods are hampered by signal-cross-talk due to volume conduction.     

Changes of Functional and Effective Connectivity in Smoking Replenishment on Deprived Heavy Smokers: A Resting-State fMRI Study

Previous researches have explored the changes of functional connectivity caused by smoking with the aid of fMRI. This study considers not only functional connectivity but also effective connectivity regarding both brain networks and brain regions by using a novel analysis framework that combines independent component analysis (ICA) and Granger causality analysis (GCA). We conducted a resting-state fMRI experiment in which twenty-one heavy smokers were scanned in two sessions of different conditions: smoking abstinence followed by smoking satiety. In our framework, group ICA was firstly adopted to obtain the spatial patterns of the default-mode network (DMN), executive-control network (ECN), and salience network (SN). Their associated time courses were analyzed using GCA, showing that the effective connectivity from SN to DMN was reduced and that from ECN/DMN to SN was enhanced after smoking replenishment. A paired t-test on ICA spatial patterns revealed functional connectivity variation in regions such as the insula, parahippocampus, precuneus, anterior cingulate cortex, supplementary motor area, and ventromedial/dorsolateral prefrontal cortex. These regions were later selected as the regions of interest (ROIs), and their effective connectivity was investigated subsequently using GCA. In smoking abstinence, the insula showed the increased effective connectivity with the other ROIs; while in smoking satiety, the parahippocampus had the enhanced inter-area effective connectivity. These results demonstrate our hypothesis that for deprived heavy smokers, smoking replenishment takes effect on both functional and effective connectivity. Moreover, our analysis framework could be applied in a range of neuroscience studies.     

Filtering noise for synchronised activity in multi-trial electrophysiology data using Wiener and Kalman filters

Novel approaches to effectively reduce noise in data recorded from multi-trial physiology experiments have been investigated using two-dimensional filtering methods, adaptive Wiener filtering and reduced update Kalman filtering. Test data based on signal and noise model consisting of different conditions of signal components mixed with noise have been considered with filtering effects evaluated using analysis of frequency coherence and of time-dependent coherence. Various situations that may affect the filtering results have been explored and reveal that Wiener and Kalman filtering can considerably improve the coherence values between two channels of multi-trial data and suppress uncorrelated components. We have extended our approach to experimental data: multi-electrode array (MEA) local field potential (LFPs) recordings from the inferotemporal cortex of sheep and LFP vs. electromyogram (LFP-EMG) recording data during resting tremor in Parkinson’s disease patients. Finally general procedures for implementation of these filtering techniques are described.     

A Brain Network Processing the Age of Faces

Age is one of the most salient aspects in faces and of fundamental cognitive and social relevance. Although face processing has been studied extensively, brain regions responsive to age have yet to be localized. Using evocative face morphs and fMRI, we segregate two areas extending beyond the previously established face-sensitive core network, centered on the inferior temporal sulci and angular gyri bilaterally, both of which process changes of facial age. By means of probabilistic tractography, we compare their patterns of functional activation and structural connectivity. The ventral portion of Wernicke''s understudied perpendicular association fasciculus is shown to interconnect the two areas, and activation within these clusters is related to the probability of fiber connectivity between them. In addition, post-hoc age-rating competence is found to be associated with high response magnitudes in the left angular gyrus. Our results provide the first evidence that facial age has a distinct representation pattern in the posterior human brain. We propose that particular face-sensitive nodes interact with additional object-unselective quantification modules to obtain individual estimates of facial age. This brain network processing the age of faces differs from the cortical areas that have previously been linked to less developmental but instantly changeable face aspects. Our probabilistic method of associating activations with connectivity patterns reveals an exemplary link that can be used to further study, assess and quantify structure-function relationships.     

Are all data types and connectivity models created equal? Validating common connectivity approaches with dispersal data

Aim

There is enormous interest in applying connectivity modelling to resistance surfaces for identifying corridors for conservation action. However, the multiple analytical approaches used to estimate resistance surfaces and predict connectivity across resistance surfaces have not been rigorously compared, and it is unclear what methods provide the best inferences about population connectivity. Using a large empirical data set on puma (Puma concolor), we are the first to compare several of the most common approaches for estimating resistance and modelling connectivity and validate them with dispersal data.

Location

Southern California, USA.

Methods

We estimate resistance using presence‐only data, GPS telemetry data from puma home ranges and genetic data using a variety of analytical methods. We model connectivity with cost distance and circuit theory algorithms. We then measure the ability of each data type and connectivity algorithm to capture GPS telemetry points of dispersing pumas.

Results

We found that resource selection functions based on GPS telemetry points and paths outperformed species distribution models when applied using cost distance connectivity algorithms. Point and path selection functions were not statistically different in their performance, but point selection functions were more sensitive to the transformation used to convert relative probability of use to resistance. Point and path selection functions and landscape genetics outperformed other methods when applied with cost distance; no methods outperformed one another with circuit theory.

Main conclusions

We conclude that path or point selection functions, or landscape genetic models, should be used to estimate landscape resistance for wildlife. In cases where resource limitations prohibit the collection of GPS collar or genetic data, our results suggest that species distribution models, while weaker, may still be sufficient for resistance estimation. We recommend the use of cost distance‐based approaches, such as least‐cost corridors and resistant kernels, for estimating connectivity and identifying functional corridors for terrestrial wildlife.

     

Integration of sleep homeostasis and navigation in Drosophila

During sleep, the brain undergoes dynamic and structural changes. In Drosophila, such changes have been observed in the central complex, a brain area important for sleep control and navigation. The connectivity of the central complex raises the question about how navigation, and specifically the head direction system, can operate in the face of sleep related plasticity. To address this question, we develop a model that integrates sleep homeostasis and head direction. We show that by introducing plasticity, the head direction system can function in a stable way by balancing plasticity in connected circuits that encode sleep pressure. With increasing sleep pressure, the head direction system nevertheless becomes unstable and a sleep phase with a different plasticity mechanism is introduced to reset network connectivity. The proposed integration of sleep homeostasis and head direction circuits captures features of their neural dynamics observed in flies and mice.     

BOLD Granger Causality Reflects Vascular Anatomy

A number of studies have tried to exploit subtle phase differences in BOLD time series to resolve the order of sequential activation of brain regions, or more generally the ability of signal in one region to predict subsequent signal in another region. More recently, such lag-based measures have been applied to investigate directed functional connectivity, although this application has been controversial. We attempted to use large publicly available datasets (FCON 1000, ADHD 200, Human Connectome Project) to determine whether consistent spatial patterns of Granger Causality are observed in typical fMRI data. For BOLD datasets from 1,240 typically developing subjects ages 7–40, we measured Granger causality between time series for every pair of 7,266 spherical ROIs covering the gray matter and 264 seed ROIs at hubs of the brain’s functional network architecture. Granger causality estimates were strongly reproducible for connections in a test and replication sample (n=620 subjects for each group), as well as in data from a single subject scanned repeatedly, both during resting and passive video viewing. The same effect was even stronger in high temporal resolution fMRI data from the Human Connectome Project, and was observed independently in data collected during performance of 7 task paradigms. The spatial distribution of Granger causality reflected vascular anatomy with a progression from Granger causality sources, in Circle of Willis arterial inflow distributions, to sinks, near large venous vascular structures such as dural venous sinuses and at the periphery of the brain. Attempts to resolve BOLD phase differences with Granger causality should consider the possibility of reproducible vascular confounds, a problem that is independent of the known regional variability of the hemodynamic response.     

Successful Reconstruction of a Physiological Circuit with Known Connectivity from Spiking Activity Alone

Identifying the structure and dynamics of synaptic interactions between neurons is the first step to understanding neural network dynamics. The presence of synaptic connections is traditionally inferred through the use of targeted stimulation and paired recordings or by post-hoc histology. More recently, causal network inference algorithms have been proposed to deduce connectivity directly from electrophysiological signals, such as extracellularly recorded spiking activity. Usually, these algorithms have not been validated on a neurophysiological data set for which the actual circuitry is known. Recent work has shown that traditional network inference algorithms based on linear models typically fail to identify the correct coupling of a small central pattern generating circuit in the stomatogastric ganglion of the crab Cancer borealis. In this work, we show that point process models of observed spike trains can guide inference of relative connectivity estimates that match the known physiological connectivity of the central pattern generator up to a choice of threshold. We elucidate the necessary steps to derive faithful connectivity estimates from a model that incorporates the spike train nature of the data. We then apply the model to measure changes in the effective connectivity pattern in response to two pharmacological interventions, which affect both intrinsic neural dynamics and synaptic transmission. Our results provide the first successful application of a network inference algorithm to a circuit for which the actual physiological synapses between neurons are known. The point process methodology presented here generalizes well to larger networks and can describe the statistics of neural populations. In general we show that advanced statistical models allow for the characterization of effective network structure, deciphering underlying network dynamics and estimating information-processing capabilities.     

Listen to Genes: Dealing with Microarray Data in the Frequency Domain

Background

We present a novel and systematic approach to analyze temporal microarray data. The approach includes normalization, clustering and network analysis of genes.

Methodology

Genes are normalized using an error model based uniform normalization method aimed at identifying and estimating the sources of variations. The model minimizes the correlation among error terms across replicates. The normalized gene expressions are then clustered in terms of their power spectrum density. The method of complex Granger causality is introduced to reveal interactions between sets of genes. Complex Granger causality along with partial Granger causality is applied in both time and frequency domains to selected as well as all the genes to reveal the interesting networks of interactions. The approach is successfully applied to Arabidopsis leaf microarray data generated from 31,000 genes observed over 22 time points over 22 days. Three circuits: a circadian gene circuit, an ethylene circuit and a new global circuit showing a hierarchical structure to determine the initiators of leaf senescence are analyzed in detail.

Conclusions

We use a totally data-driven approach to form biological hypothesis. Clustering using the power-spectrum analysis helps us identify genes of potential interest. Their dynamics can be captured accurately in the time and frequency domain using the methods of complex and partial Granger causality. With the rise in availability of temporal microarray data, such methods can be useful tools in uncovering the hidden biological interactions. We show our method in a step by step manner with help of toy models as well as a real biological dataset. We also analyse three distinct gene circuits of potential interest to Arabidopsis researchers.     

Environment-Related Adaptive Changes of Gut Commensal Microbiota Do not Alter Colonic Toll-Like Receptors but Modulate the Local Expression of Sensory-Related Systems in Rats

Pathogenic and protective roles have been attributed to gut commensal microbiota (GCM) in gastrointestinal inflammatory and functional disorders. We have shown that the adaptation to a new environment implies specific changes in the composition of GCM. Here we assessed if environment-related adaptive changes of GCM modulate the expression of colonic Toll-like receptors (TLRs) and sensory-related systems in rats. Adult male SD rats were maintained under different environmental conditions: barrier-breed-and-maintained, barrier-breed adapted to conventional conditions or conventional-breed-and-maintained. Fluorescent in situ hybridization and real-time quantitative PCR (qPCR) were used to characterize luminal ceco-colonic microbiota. Colonic expression of TLR2, TLR4, TLR5, and TLR7, cannabinoid receptors (CB1/CB2), μ-opioid receptor (MOR), transient receptor potential vanilloid (TRPV1, TRPV3, and TRPV4), protease-activated receptor 2 (PAR-2), and calcitonin gene-related peptide were quantified by RT-qPCR. CB1, CB2 and MOR expression, was evaluated also by immunohistochemistry. In rats, housing-related environmental conditions induce specific changes of GCM, without impact on the expression of TLR-dependent bacterial recognition systems. Expression of sensory-related markers (MOR, TRPV3, PAR-2, and CB2) decreased with the adaptation to a conventional environment, correlating with changes in Bacteroides spp., Lactobacillus spp., and Bifidobacterium spp. counts. This suggests an interaction between GCM and visceral sensory mechanisms, which might be part of the mechanisms underlying the beneficial effects of some bacterial groups on functional and inflammatory gastrointestinal disorders.     

Estimation of effective connectivity using multi-layer perceptron artificial neural network

Studies on interactions between brain regions estimate effective connectivity, (usually) based on the causality inferences made on the basis of temporal precedence. In this study, the causal relationship is modeled by a multi-layer perceptron feed-forward artificial neural network, because of the ANN’s ability to generate appropriate input–output mapping and to learn from training examples without the need of detailed knowledge of the underlying system. At any time instant, the past samples of data are placed in the network input, and the subsequent values are predicted at its output. To estimate the strength of interactions, the measure of “Causality coefficient” is defined based on the network structure, the connecting weights and the parameters of hidden layer activation function. Simulation analysis demonstrates that the method, called “CREANN” (Causal Relationship Estimation by Artificial Neural Network), can estimate time-invariant and time-varying effective connectivity in terms of MVAR coefficients. The method shows robustness with respect to noise level of data. Furthermore, the estimations are not significantly influenced by the model order (considered time-lag), and the different initial conditions (initial random weights and parameters of the network). CREANN is also applied to EEG data collected during a memory recognition task. The results implicate that it can show changes in the information flow between brain regions, involving in the episodic memory retrieval process. These convincing results emphasize that CREANN can be used as an appropriate method to estimate the causal relationship among brain signals.     

Detecting Genetic Association of Common Human Facial Morphological Variation Using High Density 3D Image Registration

Human facial morphology is a combination of many complex traits. Little is known about the genetic basis of common facial morphological variation. Existing association studies have largely used simple landmark-distances as surrogates for the complex morphological phenotypes of the face. However, this can result in decreased statistical power and unclear inference of shape changes. In this study, we applied a new image registration approach that automatically identified the salient landmarks and aligned the sample faces using high density pixel points. Based on this high density registration, three different phenotype data schemes were used to test the association between the common facial morphological variation and 10 candidate SNPs, and their performances were compared. The first scheme used traditional landmark-distances; the second relied on the geometric analysis of 15 landmarks and the third used geometric analysis of a dense registration of ∼30,000 3D points. We found that the two geometric approaches were highly consistent in their detection of morphological changes. The geometric method using dense registration further demonstrated superiority in the fine inference of shape changes and 3D face modeling. Several candidate SNPs showed potential associations with different facial features. In particular, one SNP, a known risk factor of non-syndromic cleft lips/palates, rs642961 in the IRF6 gene, was validated to strongly predict normal lip shape variation in female Han Chinese. This study further demonstrated that dense face registration may substantially improve the detection and characterization of genetic association in common facial variation.