brainlife.io: a decentralized and open-source cloud platform to support neuroscience research.

Neuroscience is advancing standardization and tool development to support rigor and transparency. Consequently, data pipeline complexity has increased, hindering FAIR (findable, accessible, interoperable and reusable) access. brainlife.io was developed to democratize neuroimaging research. The platform provides data standardization, management, visualization and processing and automatically tracks the provenance history of thousands of data objects. Here, brainlife.io is described and evaluated for validity, reliability, reproducibility, replicability and scientific utility using four data modalities and 3,200 participants.

Task matters: Individual MEG signatures from naturalistic and neurophysiological brain states.

The discovery that human brain connectivity data can be used as a "fingerprint" to identify a given individual from a population, has become a burgeoning research area in the neuroscience field. Recent studies have identified the possibility to extract these brain signatures from the temporal rich dynamics of resting-state magneto encephalography (MEG) recordings. Nevertheless, it is still uncertain to what extent MEG signatures can serve as an indicator of human identifiability during task-related conduct. Here, using MEG data from naturalistic and neurophysiological tasks, we show that identification improves in tasks relative to resting-state, providing compelling evidence for a task dependent axis of MEG signatures. Notably, improvements in identifiability were more prominent in strictly controlled tasks. Lastly, the brain regions contributing most towards individual identification were also modified when engaged in task activities. We hope that this investigation advances our understanding of the driving factors behind brain identification from MEG signals.

Dynamic causal modelling shows a prominent role of local inhibition in alpha power modulation in higher visual cortex.

During resting-state EEG recordings, alpha activity is more prominent over the posterior cortex in eyes-closed (EC) conditions compared to eyes-open (EO). In this study, we characterized the difference in spectra between EO and EC conditions using dynamic causal modelling. Specifically, we investigated the role of intrinsic and extrinsic connectivity-within the visual cortex-in generating EC-EO alpha power differences over posterior electrodes. The primary visual cortex (V1) and the bilateral middle temporal visual areas (V5) were equipped with bidirectional extrinsic connections using a canonical microcircuit. The states of four intrinsically coupled subpopulations-within each occipital source-were also modelled. Using Bayesian model selection, we tested whether modulations of the intrinsic connections in V1, V5 or extrinsic connections (or a combination thereof) provided the best evidence for the data. In addition, using parametric empirical Bayes (PEB), we estimated group averages under the winning model. Bayesian model selection showed that the winning model contained both extrinsic connectivity modulations, as well as intrinsic connectivity modulations in all sources. The PEB analysis revealed increased extrinsic connectivity during EC. Overall, we found a reduction in the inhibitory intrinsic connections during EC. The results suggest that the intrinsic modulations in V5 played the most important role in producing EC-EO alpha differences, suggesting an intrinsic disinhibition in higher order visual cortex, during EC resting state.

Auditory driven gamma synchrony is associated with cortical thickness in widespread cortical areas.

OBJECTIVE: Gamma synchrony is a fundamental functional property of the cerebral cortex, impaired in multiple neuropsychiatric conditions (i.e. schizophrenia, Alzheimer's disease, stroke etc.). Auditory stimulation in the gamma range allows to drive gamma synchrony of the entire cortical mantle and to estimate the efficiency of the mechanisms sustaining it. As gamma synchrony depends strongly on the interplay between parvalbumin-positive interneurons and pyramidal neurons, we hypothesize an association between cortical thickness and gamma synchrony. To test this hypothesis, we employed a combined magnetoencephalography (MEG) - Magnetic Resonance Imaging (MRI) study. METHODS: Cortical thickness was estimated from anatomical MRI scans. MEG measurements related to exposure of 40 Hz amplitude modulated tones were projected onto the cortical surface. Two measures of cortical synchrony were considered: (a) inter-trial phase consistency at 40 Hz, providing a vertex-wise estimation of gamma synchronization, and (b) phase-locking values between primary auditory cortices and whole cortical mantle, providing a measure of long-range cortical synchrony. A correlation between cortical thickness and synchronization measures was then calculated for 72 MRI-MEG scans. RESULTS: Both inter-trial phase consistency and phase locking values showed a significant positive correlation with cortical thickness. For inter-trial phase consistency, clusters of strong associations were found in the temporal and frontal lobes, especially in the bilateral auditory and pre-motor cortices. Higher phase-locking values corresponded to higher cortical thickness in the frontal, temporal, occipital and parietal lobes. DISCUSSION AND CONCLUSIONS: In healthy subjects, a thicker cortex corresponds to higher gamma synchrony and connectivity in the primary auditory cortex and beyond, likely reflecting underlying cell density involved in gamma circuitries. This result hints towards an involvement of gamma synchrony together with underlying brain structure in brain areas for higher order cognitive functions. This study contributes to the understanding of inherent cortical functional and structural brain properties, which might in turn constitute the basis for the definition of useful biomarkers in patients showing aberrant gamma synchronization.

Brain simulation as a cloud service: The Virtual Brain on EBRAINS.

The Virtual Brain (TVB) is now available as open-source services on the cloud research platform EBRAINS (ebrains.eu). It offers software for constructing, simulating and analysing brain network models including the TVB simulator; magnetic resonance imaging (MRI) processing pipelines to extract structural and functional brain networks; combined simulation of large-scale brain networks with small-scale spiking networks; automatic conversion of user-specified model equations into fast simulation code; simulation-ready brain models of patients and healthy volunteers; Bayesian parameter optimization in epilepsy patient models; data and software for mouse brain simulation; and extensive educational material. TVB cloud services facilitate reproducible online collaboration and discovery of data assets, models, and software embedded in scalable and secure workflows, a precondition for research on large cohort data sets, better generalizability, and clinical translation.

rsHRF: A toolbox for resting-state HRF estimation and deconvolution.

The hemodynamic response function (HRF) greatly influences the intra- and inter-subject variability of brain activation and connectivity, and might confound the estimation of temporal precedence in connectivity analyses, making its estimation necessary for a correct interpretation of neuroimaging studies. Additionally, the HRF shape itself is a useful local measure. However, most algorithms for HRF estimation are specific for task-related fMRI data, and only a few can be directly applied to resting-state protocols. Here we introduce rsHRF, a Matlab and Python toolbox that implements HRF estimation and deconvolution from the resting-state BOLD signal. We first provide an overview of the main algorithm, practical implementations, and then demonstrate the feasibility and usefulness of rsHRF by validation experiments with a publicly available resting-state fMRI dataset. We also provide tools for statistical analyses and visualization. We believe that this toolbox may significantly contribute to a better analysis and understanding of the components and variability of BOLD signals.

How hot is the hot zone? Computational modelling clarifies the role of parietal and frontoparietal connectivity during anaesthetic-induced loss of consciousness.

In recent years, specific cortical networks have been proposed to be crucial for sustaining consciousness, including the posterior hot zone and frontoparietal resting state networks (RSN). Here, we computationally evaluate the relative contributions of three RSNs - the default mode network (DMN), the salience network (SAL), and the central executive network (CEN) - to consciousness and its loss during propofol anaesthesia. Specifically, we use dynamic causal modelling (DCM) of 10 min of high-density EEG recordings (N = 10, 4 males) obtained during behavioural responsiveness, unconsciousness and post-anaesthetic recovery to characterise differences in effective connectivity within frontal areas, the posterior 'hot zone', frontoparietal connections, and between-RSN connections. We estimate - for the first time - a large DCM model (LAR) of resting EEG, combining the three RSNs into a rich club of interconnectivity. Consistent with the hot zone theory, our findings demonstrate reductions in inter-RSN connectivity in the parietal cortex. Within the DMN itself, the strongest reductions are in feed-forward frontoparietal and parietal connections at the precuneus node. Within the SAL and CEN, loss of consciousness generates small increases in bidirectional connectivity. Using novel DCM leave-one-out cross-validation, we show that the most consistent out-of-sample predictions of the state of consciousness come from a key set of frontoparietal connections. This finding also generalises to unseen data collected during post-anaesthetic recovery. Our findings provide new, computational evidence for the importance of the posterior hot zone in explaining the loss of consciousness, highlighting also the distinct role of frontoparietal connectivity in underpinning conscious responsiveness, and consequently, suggest a dissociation between the mechanisms most prominently associated with explaining the contrast between conscious awareness and unconsciousness, and those maintaining consciousness.

Deletion of Autism Risk Gene Shank3 Disrupts Prefrontal Connectivity.

Mutations in the synaptic scaffolding protein SHANK3 are a major cause of autism and are associated with prominent intellectual and language deficits. However, the neural mechanisms whereby SHANK3 deficiency affects higher-order socio-communicative functions remain unclear. Using high-resolution functional and structural MRI in adult male mice, here we show that loss of Shank3 (Shank3B-/-) results in disrupted local and long-range prefrontal and frontostriatal functional connectivity. We document that prefrontal hypoconnectivity is associated with reduced short-range cortical projections density, and reduced gray matter volume. Finally, we show that prefrontal disconnectivity is predictive of social communication deficits, as assessed with ultrasound vocalization recordings. Collectively, our results reveal a critical role of SHANK3 in the development of prefrontal anatomy and function, and suggest that SHANK3 deficiency may predispose to intellectual disability and socio-communicative impairments via dysregulation of higher-order cortical connectivity.SIGNIFICANCE STATEMENT Mutations in the synaptic scaffolding protein SHANK3 are commonly associated with autism, intellectual, and language deficits. Previous research has linked SHANK3 deficiency to basal ganglia dysfunction, motor stereotypies, and social deficits. However, the neural mechanism whereby Shank3 gene mutations affects cortical functional connectivity and higher-order socio-communicative functions remain unclear. Here we show that loss of SHANK3 in mice results in largely disrupted functional connectivity and abnormal gray matter anatomy in prefrontal areas. We also show that prefrontal connectivity disruption is tightly linked to socio-communicative deficits. Our findings suggest that SHANK3 is a critical orchestrator of frontocortical function, and that disrupted connectivity of prefrontal areas may underpin socio-communicative impairments observed in SHANK3 mutation carriers.

The effect of global signal regression on DCM estimates of noise and effective connectivity from resting state fMRI.

The influence of global BOLD fluctuations on resting state functional connectivity in fMRI data remains a topic of debate, with little consensus. In this study, we assessed the effects of global signal regression (GSR) on effective connectivity within and between resting state networks (RSNs) - as estimated with dynamic causal modelling (DCM) for resting state fMRI (rsfMRI). DCM incorporates a forward (generative) model that quantifies the contribution of different types of noise (including global measurement noise), effective connectivity, and (neuro)vascular processes to functional connectivity measurements. DCM analyses were applied to two different designs; namely, longitudinal and cross-sectional designs. In the modelling of longitudinal designs, we considered four extensive longitudinal resting state fMRI datasets with a total number of 20 subjects. In the analysis of cross-sectional designs, we used rsfMRI data from 361 subjects from the Human Connectome Project. We hypothesized that (1) GSR would have no discernible impact on effective connectivity estimated with DCM, and (2) GSR would be reflected in the parameters representing global measurement noise. Additionally, we performed comparative analyses of information gain with and without GSR. Our results showed negligible to small effects of GSR on effective connectivity within small (separately estimated) RSNs. However, although the effect sizes were small, there was substantial to conclusive evidence for an effect of GSR on connectivity parameters. For between-network connectivity, we found two important effects: the effect of GSR on between-network effective connectivity (averaged over all connections) was negligible to small, while the effect of GSR on individual connections was non-negligible. In the cross-sectional (but not in the longitudinal) data, some connections showed substantial to conclusive evidence for an effect of GSR. Contrary to our expectations, we found either no effect (in the longitudinal designs) or a non-specific (cross-sectional design) effect of GSR on parameters characterising (global) measurement noise. Data without GSR were found to be more informative than data with GSR; however, in small resting state networks the precision of posterior estimates was greater after GSR. In conclusion, GSR is a minor concern in DCM studies; however, quantitative interpretation of between-network connections (as opposed to average between-network connectivity) and noise parameters should be treated with some caution. The Kullback-Leibler divergence of the posterior from the prior (i.e., information gain) - together with the precision of posterior estimates - might offer a useful measure to assess the appropriateness of GSR in resting state fMRI.

Effective connectivity modulations related to win and loss outcomes.

Previous studies have characterized the brain regions involved in encoding monetary reward and punishment outcomes. The question of how this information is integrated across brain regions has received less attention. Here, we investigated changes in effective connectivity related to the processing of positive and negative monetary outcomes using functional magnetic resonance imaging data from the Human Connectome Project. Specifically, subjects engaged in a card guessing game which could yield win, loss, or neutral outcomes. A general linear model was used to define a network of regions involved in win and loss outcome processing, including anterior insula, anterior cingulate cortex, and ventral striatum. Dynamic causal modelling (DCM) was implemented to study between-region couplings and outcome-related modulations thereof within this network. In addition, we explored the relation between effective connectivity patterns and choice behavior in the gambling task. Parametric empirical Bayesian modelling was conducted for group-level inferences of both DCM and the choice behavior. Behaviorally, both win and loss outcomes increased the probability of choice switches in subsequent gambles. In terms of connectivity, win outcomes were associated with increased extrinsic connectivity across the network, while loss outcomes featured a balance between increased and decreased extrinsic connectivity. Moreover, self-inhibitory connections tended to decrease for both win and loss outcomes. Interestingly, a substantial discrepancy was observed for occipital cortex connectivity, which was characterized by intrinsic disinhibition in loss but not in win trials. The observed differences in effective connectivity during the processing of positive and negative outcomes, despite similarities in average regional activity and choice behavior, highlight the value of exploring network dynamics in the context of incentive manipulations.

Modeling brain dynamics after tumor resection using The Virtual Brain.

Brain tumor patients scheduled for tumor resection often face significant uncertainty, as the outcome of neurosurgery is difficult to predict at the individual patient level. Recently, simulation of the activity of neural populations connected according to the white matter fibers, producing personalized brain network models, has been introduced as a promising tool for this purpose. The Virtual Brain provides a robust open source framework to implement these models. However, brain network models first have to be validated, before they can be used to predict brain dynamics. In prior work, we optimized individual brain network model parameters to maximize the fit with empirical brain activity. In this study, we extend this line of research by examining the stability of fitted parameters before and after tumor resection, and compare it with baseline parameter variability using data from healthy control subjects. Based on these findings, we perform the first "virtual neurosurgery", mimicking patient's actual surgery by removing white matter fibers in the resection mask and simulating again neural activity on this new connectome. We find that brain network model parameters are relatively stable over time in brain tumor patients who underwent tumor resection, compared with baseline variability in healthy control subjects. Concerning the virtual neurosurgery analyses, use of the pre-surgery model implemented on the virtually resected structural connectome resulted in improved similarity with post-surgical empirical functional connectivity in some patients, but negligible improvement in others. These findings reveal interesting avenues for increasing interactions between computational neuroscience and neuro-oncology, as well as important limitations that warrant further investigation.

Disambiguating the role of blood flow and global signal with partial information decomposition.

Global signal (GS) is an ubiquitous construct in resting state functional magnetic resonance imaging (rs-fMRI), associated to nuisance, but containing by definition most of the neuronal signal. Global signal regression (GSR) effectively removes the impact of physiological noise and other artifacts, but at the same time it alters correlational patterns in unpredicted ways. Performing GSR taking into account the underlying physiology (mainly the blood arrival time) has been proven to be beneficial. From these observations we aimed to: 1) characterize the effect of GSR on network-level functional connectivity in a large dataset; 2) assess the complementary role of global signal and vessels; and 3) use the framework of partial information decomposition to further look into the joint dynamics of the global signal and vessels, and their respective influence on the dynamics of cortical areas. We observe that GSR affects intrinsic connectivity networks in the connectome in a non-uniform way. Furthermore, by estimating the predictive information of blood flow and the global signal using partial information decomposition, we observe that both signals are present in different amounts across intrinsic connectivity networks. Simulations showed that differences in blood arrival time can largely explain this phenomenon, while using hemodynamic and calcium mouse recordings we were able to confirm the presence of vascular effects, as calcium recordings lack hemodynamic information. With these results we confirm network-specific effects of GSR and the importance of taking blood flow into account for improving de-noising methods. Additionally, and beyond the mere issue of data denoising, we quantify the diverse and complementary effect of global and vessel BOLD signals on the dynamics of cortical areas.

Objects seen as scenes: Neural circuitry for attending whole or parts.

Depending on our goals, we pay attention to the global shape of an object or to the local shape of its parts, since it's difficult to do both at once. This typically effortless process can be impaired in disease. However, it is not clear which cortical regions carry the information needed to constrain shape processing to a chosen global/local level. Here, novel stimuli were used to dissociate functional MRI responses to global and local shapes. This allowed identification of cortical regions containing information about level (independent from shape). Crucially, these regions overlapped part of the cortical network implicated in scene processing. As expected, shape information (independent of level) was mainly located in category-selective areas specialized for object- and face-processing. Regions with the same informational profile were strongly linked (as measured by functional connectivity), but were weak when the profiles diverged. Specifically, in the ventral-temporal-cortex (VTC) regions favoring level and shape were consistently separated by the mid-fusiform sulcus (MFS). These regions also had limited crosstalk despite their spatial proximity, thus defining two functional pathways within VTC. We hypothesize that object hierarchical level is processed by neural circuitry that also analyses spatial layout in scenes, contributing to the control of the spatial-scale used for shape recognition. Use of level information tolerant to shape changes could guide whole/part attentional selection but facilitate illusory shape/level conjunctions under impoverished vision.

Synergistic Information Transfer in the Global System of Financial Markets.

Uncovering dynamic information flow between stock market indices has been the topic of several studies which exploited the notion of transfer entropy or Granger causality, its linear version. The output of the transfer entropy approach is a directed weighted graph measuring the information about the future state of each target provided by the knowledge of the state of each driving stock market index. In order to go beyond the pairwise description of the information flow, thus looking at higher order informational circuits, here we apply the partial information decomposition to triplets consisting of a pair of driving markets (belonging to America or Europe) and a target market in Asia. Our analysis, on daily data recorded during the years 2000 to 2019, allows the identification of the synergistic information that a pair of drivers carry about the target. By studying the influence of the closing returns of drivers on the subsequent overnight changes of target indexes, we find that (i) Korea, Tokyo, Hong Kong, and Singapore are, in order, the most influenced Asian markets; (ii) US indices SP500 and Russell are the strongest drivers with respect to the bivariate Granger causality; and (iii) concerning higher order effects, pairs of European and American stock market indices play a major role as the most synergetic three-variables circuits. Our results show that the Synergy, a proxy of higher order predictive information flow rooted in information theory, provides details that are complementary to those obtained from bivariate and global Granger causality, and can thus be used to get a better characterization of the global financial system.

Information decomposition of multichannel EMG to map functional interactions in the distributed motor system.

The central nervous system needs to coordinate multiple muscles during postural control. Functional coordination is established through the neural circuitry that interconnects different muscles. Here we used multivariate information decomposition of multichannel EMG acquired from 14 healthy participants during postural tasks to investigate the neural interactions between muscles. A set of information measures were estimated from an instantaneous linear regression model and a time-lagged VAR model fitted to the EMG envelopes of 36 muscles. We used network analysis to quantify the structure of functional interactions between muscles and compared them across experimental conditions. Conditional mutual information and transfer entropy revealed sparse networks dominated by local connections between muscles. We observed significant changes in muscle networks across postural tasks localized to the muscles involved in performing those tasks. Information decomposition revealed distinct patterns in task-related changes: unimanual and bimanual pointing were associated with reduced transfer to the pectoralis major muscles, but an increase in total information compared to no pointing, while postural instability resulted in increased information, information transfer and information storage in the abductor longus muscles compared to normal stability. These findings show robust patterns of directed interactions between muscles that are task-dependent and can be assessed from surface EMG recorded during static postural tasks. We discuss directed muscle networks in terms of the neural circuitry involved in generating muscle activity and suggest that task-related effects may reflect gain modulations of spinal reflex pathways.

Dynamic causal modelling of fluctuating connectivity in resting-state EEG.

Functional and effective connectivity are known to change systematically over time. These changes might be explained by several factors, including intrinsic fluctuations in activity-dependent neuronal coupling and contextual factors, like experimental condition and time. Furthermore, contextual effects may be subject-specific or conserved over subjects. To characterize fluctuations in effective connectivity, we used dynamic causal modelling (DCM) of cross spectral responses over 1- min of electroencephalogram (EEG) recordings during rest, divided into 1-sec windows. We focused on two intrinsic networks: the default mode and the saliency network. DCM was applied to estimate connectivity in each time-window for both networks. Fluctuations in DCM connectivity parameters were assessed using hierarchical parametric empirical Bayes (PEB). Within-subject, between-window effects were modelled with a second-level linear model with temporal basis functions as regressors. This procedure was conducted for every subject separately. Bayesian model reduction was then used to assess which (combination of) temporal basis functions best explain dynamic connectivity over windows. A third (between-subject) level model was used to infer which dynamic connectivity parameters are conserved over subjects. Our results indicate that connectivity fluctuations in the default mode network and to a lesser extent the saliency network comprised both subject-specific components and a common component. For both networks, connections to higher order regions appear to monotonically increase during the 1- min period. These results not only establish the predictive validity of dynamic connectivity estimates - in virtue of detecting systematic changes over subjects - they also suggest a network-specific dissociation in the relative contribution of fluctuations in connectivity that depend upon experimental context. We envisage these procedures could be useful for characterizing brain state transitions that may be explained by their cognitive or neuropathological underpinnings.

Modulation of the spontaneous hemodynamic response function across levels of consciousness.

Functional imaging research has already contributed with several results to the study of neural correlates of consciousness. Apart from task-related activation derived in fMRI, PET based glucose metabolism rate or cerebral blood flow account for a considerable proportion of the study of brain activity under different levels of consciousness. Resting state functional connectivity MRI is playing a crucial role to explore the consciousness related functional integration, successfully complementing PET, another widely used neuroimaging technique. Here, spontaneous hemodynamic response is introduced to characterize resting state brain activity giving information on the local metabolism (neurovascular coupling), and useful to improve the time-resolved activity and connectivity measures based on BOLD fMRI. This voxel-wise measure is then used to investigate the loss of consciousness under Propofol anesthesia and unresponsive wakefulness syndrome. Changes in the hemodynamic response in precuneus and posterior cingulate are found to be a common principle underlying loss of consciousness in both conditions. The thalamus appears to be less obviously modulated by Propofol, compared with frontoparietal regions. However, a significant increase in spontaneous thalamic hemodynamic response was found in patients in unresponsive wakefulness syndrome compared with healthy controls. Our results ultimately show that anesthesia- or pathology-induced neurovascular coupling could be tracked by modulated spontaneous hemodynamic response derived from resting state fMRI.

The interplay between emotion regulation, emotional well-being, and cognitive functioning in brain tumor patients and their caregivers: An exploratory study.

OBJECTIVE: Brain tumor patients may suffer from a range of health-impairing problems reducing their quality of life. To identify potential targets for interventions, we examined the influence of different emotion regulation strategies on affective and cognitive functioning as indices of quality of life in patients and their caregivers in the early phase of treatment. METHODS: To this end, we conducted an exploratory longitudinal study on a small cohort, measuring emotion regulation, emotional well-being, and cognitive functioning on the day before each patient's tumor resection (28 patients and 11 caregivers) and several months after neurosurgery (22 patients and 10 caregivers). RESULTS: Results showed that emotion regulation strategies are relatively stable from preoperative to postoperative assessment. Nevertheless, several associations between emotion regulation strategies and quality of life indices were evident after tumor resection. In particular, our results were largely in line with previous research findings in healthy and other patient populations, corroborating the adaptive character of cognitive reappraisal, whereas suppression and expression of emotions were related to reduced cognitive and affective functioning, respectively. CONCLUSIONS: Based on these results, we suggest that further intervention or qualitative studies explore whether therapeutic interventions directed toward mastery of cognitive reappraisal techniques and appropriate expression of emotions could lead to improved long-term adjustment among brain tumor patients and their caregivers.

Estimating Directed Connectivity from Cortical Recordings and Reconstructed Sources.

In cognitive neuroscience, electrical brain activity is most commonly recorded at the scalp. In order to infer the contributions and connectivity of underlying neuronal sources within the brain, it is necessary to reconstruct sensor data at the source level. Several approaches to this reconstruction have been developed, thereby solving the so-called implicit inverse problem Michel et al. (Clin Neurophysiol 115:2195-2222, 2004). However, a unifying premise against which to validate these source reconstructions is seldom available. The dataset provided in this work, in which brain activity is simultaneously recorded on the scalp (non-invasively) by electroencephalography (EEG) and on the cortex (invasively) by electrocorticography (ECoG), can be of a great help in this direction. These multimodal recordings were obtained from a macaque monkey under wakefulness and sedation. Our primary goal was to establish the connectivity architecture between two sources of interest (frontal and parietal), and to assess how their coupling changes over the conditions. We chose these sources because previous studies have shown that the connections between them are modified by anaesthesia Boly et al. (J Neurosci 32:7082-7090, 2012). Our secondary goal was to evaluate the consistency of the connectivity results when analyzing sources recorded from invasive data (128 implanted ECoG sources) and source activity reconstructed from scalp recordings (19 EEG sensors) at the same locations as the ECoG sources. We conclude that the directed connectivity in the frequency domain between cortical sources reconstructed from scalp EEG is qualitatively similar to the connectivity inferred directly from cortical recordings, using both data-driven (directed transfer function) and biologically grounded (dynamic causal modelling) methods. Furthermore, the connectivity changes identified were consistent with previous findings Boly et al. (J Neurosci 32:7082-7090, 2012). Our findings suggest that inferences about directed connectivity based upon non-invasive electrophysiological data have construct validity in relation to invasive recordings.