Lead/Lag directionality is not generally equivalent to causality in nonlinear systems: Comparison of phase slope index and conditional mutual information.

Applications of causal techniques to neural time series have increased extensively over last decades, including a wide and diverse family of methods focusing on electroencephalogram (EEG) analysis. Besides connectivity inferred in defined frequency bands, there is a growing interest in the analysis of cross-frequency interactions, in particular phase and amplitude coupling and directionality. Some studies show contradicting results of coupling directionality from high frequency to low frequency signal components, in spite of generally considered modulation of a high-frequency amplitude by a low-frequency phase. We have compared two widely used methods to estimate the directionality in cross frequency coupling: conditional mutual information (CMI) and phase slope index (PSI). The latter, applied to infer cross-frequency phase-amplitude directionality from animal intracranial recordings, gives opposite results when comparing to CMI. Both metrics were tested in a numerically simulated example of unidirectionally coupled Rössler systems, which helped to find the explanation of the contradictory results: PSI correctly estimates the lead/lag relationship which, however, is not generally equivalent to causality in the sense of directionality of coupling in nonlinear systems, correctly inferred by using CMI with surrogate data testing.

Identification of Early Hippocampal Dynamics during Recognition Memory with Independent Component Analysis.

The hippocampus is generally considered to have relatively late involvement in recognition memory, its main electrophysiological signature being between 400 and 800 ms after stimulus onset. However, most electrophysiological studies have analyzed the hippocampus as a single responsive area, selecting only a single-site signal exhibiting the strongest effect in terms of amplitude. These classical approaches may not capture all the dynamics of this structure, hindering the contribution of other hippocampal sources that are not located in the vicinity of the selected site. We combined intracerebral electroencephalogram recordings from epileptic patients with independent component analysis during a recognition memory task involving the recognition of old and new images. We identified two sources with different responses emerging from the hippocampus: a fast one (maximal amplitude at ∼250 ms) that could not be directly identified from raw recordings and a latter one, peaking at ∼400 ms. The former component presented different amplitudes between old and new items in 6 out of 10 patients. The latter component had different delays for each condition, with a faster activation (∼290 ms after stimulus onset) for recognized items. We hypothesize that both sources represent two steps of hippocampal recognition memory, the faster reflecting the input from other structures and the latter the hippocampal internal processing. Recognized images evoking early activations would facilitate neural computation in the hippocampus, accelerating memory retrieval of complementary information. Overall, our results suggest that the hippocampal activity is composed of several sources with an early activation related to recognition memory.

Combining independent component analysis and source localization for improving spatial sampling of stereoelectroencephalography in epilepsy.

Stereoelectroencephalography is a powerful intracerebral EEG recording method for the presurgical evaluation of epilepsy. It consists in implanting depth electrodes in the patient's brain to record electrical activity and map the epileptogenic zone, which should be resected to render the patient seizure-free. Stereoelectroencephalography has high spatial accuracy and signal-to-noise ratio but remains limited in the coverage of the explored brain regions. Thus, the implantation might provide a suboptimal sampling of epileptogenic regions. We investigate the potential of improving a suboptimal stereoelectroencephalography recording by performing source localization on stereoelectroencephalography signals. We propose combining independent component analysis, connectivity measures to identify components of interest, and distributed source modelling. This approach was tested on two patients with two implantations each, the first failing to characterize the epileptogenic zone and the second giving a better diagnosis. We demonstrate that ictal and interictal source localization performed on the first stereoelectroencephalography recordings matches the findings of the second stereo-EEG exploration. Our findings suggest that independent component analysis followed by source localization on the topographies of interest is a promising method for retrieving the epileptogenic zone in case of suboptimal implantation.

Studying memory processes at different levels with simultaneous depth and surface EEG recordings.

Investigating cognitive brain functions using non-invasive electrophysiology can be challenging due to the particularities of the task-related EEG activity, the depth of the activated brain areas, and the extent of the networks involved. Stereoelectroencephalographic (SEEG) investigations in patients with drug-resistant epilepsy offer an extraordinary opportunity to validate information derived from non-invasive recordings at macro-scales. The SEEG approach can provide brain activity with high spatial specificity during tasks that target specific cognitive processes (e.g., memory). Full validation is possible only when performing simultaneous scalp SEEG recordings, which allows recording signals in the exact same brain state. This is the approach we have taken in 12 subjects performing a visual memory task that requires the recognition of previously viewed objects. The intracranial signals on 965 contact pairs have been compared to 391 simultaneously recorded scalp signals at a regional and whole-brain level, using multivariate pattern analysis. The results show that the task conditions are best captured by intracranial sensors, despite the limited spatial coverage of SEEG electrodes, compared to the whole-brain non-invasive recordings. Applying beamformer source reconstruction or independent component analysis does not result in an improvement of the multivariate task decoding performance using surface sensor data. By analyzing a joint scalp and SEEG dataset, we investigated whether the two types of signals carry complementary information that might improve the machine-learning classifier performance. This joint analysis revealed that the results are driven by the modality exhibiting best individual performance, namely SEEG.

Reconstruction and localization of auditory sources from intracerebral SEEG using independent component analysis.

Stereo-electroencephalography (SEEG) is the surgical implantation of electrodes in the brain to better localize the epileptic network in pharmaco-resistant epileptic patients. This technique has exquisite spatial and temporal resolution. Still, the number and the position of the electrodes in the brain is limited and determined by the semiology and/or preliminary non-invasive examinations, leading to a large number of unexplored brain structures in each patient. Here, we propose a new approach to reconstruct the activity of non-sampled structures in SEEG, based on independent component analysis (ICA) and dipole source localization. We have tested this approach with an auditory stimulation dataset in ten patients. The activity directly recorded from the auditory cortex served as ground truth and was compared to the ICA applied on all non-auditory electrodes. Our results show that the activity from the auditory cortex can be reconstructed at the single trial level from contacts as far as ∼40 mm from the source. Importantly, this reconstructed activity is localized via dipole fitting in the proximity of the original source. In addition, we show that the size of the confidence interval of the dipole fitting is a good indicator of the reliability of the result, which depends on the geometry of the SEEG implantation. Overall, our approach allows reconstructing the activity of structures far from the electrode locations, partially overcoming the spatial sampling limitation of intracerebral recordings.

Comparison of beamformer and ICA for dynamic connectivity analysis: A simultaneous MEG-SEEG study.

Magnetoencephalography (MEG) is a powerful tool for estimating brain connectivity with both good spatial and temporal resolution. It is particularly helpful in epilepsy to characterize non-invasively the epileptic networks. However, using MEG to map brain networks requires solving a difficult inverse problem that introduces uncertainty in the activity localization and connectivity measures. Our goal here was to compare independent component analysis (ICA) followed by dipole source localization and the linearly constrained minimum-variance beamformer (LCMV-BF) for characterizing regions with interictal epileptic activity and their dynamic connectivity. After a simulation study, we compared ICA and LCMV-BF results with intracerebral EEG (stereotaxic EEG, SEEG) recorded simultaneously in 8 epileptic patients, which provide a unique 'ground truth' to which non-invasive results can be confronted. We compared the signal time courses extracted applying ICA and LCMV-BF on MEG data to that of SEEG, both for the actual signals and the dynamic connectivity computed using cross-correlation (evolution of links in time). With our simulations, we illustrated the different effect of the temporal and spatial correlation among sources on the two methods. While ICA was more affected by the temporal correlation but robust against spatial configurations, LCMV-BF showed opposite behavior. Moreover, ICA seems more suited to retrieve the simulated networks. In case of real patient data, good MEG/SEEG correlation and good localization were obtained in 6 out of 8 patients. In 4 of them ICA had the best performance (higher correlation, lower localization distance). In terms of dynamic connectivity, the evolution in time of the cross-correlation links could be retrieved in 5 patients out of 6, however, with more variable results in terms of correlation and distance. In two patients LCMV-BF had better results than ICA. In one patient the two methods showed equally good outcomes, and in the remaining two patients ICA performed best. In conclusion, our results obtained by exploiting simultaneous MEG/SEEG recordings suggest that ICA and LCMV-BF have complementary qualities for retrieving the dynamics of interictal sources and their network interactions.

Site-dependent shaping of field potential waveforms.

The activity of neuron populations gives rise to field potentials (FPs) that extend beyond the sources. Their mixing in the volume dilutes the original temporal motifs in a site-dependent manner, a fact that has received little attention. And yet, it potentially rids of physiological significance the time-frequency parameters of individual waves (amplitude, phase, duration). This is most likely to happen when a single source or a local origin is erroneously assumed. Recent studies using spatial treatment of these signals and anatomically realistic modeling of neuron aggregates provide convincing evidence for the multisource origin and site-dependent blend of FPs. Thus, FPs generated in primary structures like the neocortex and hippocampus reach far and cross-contaminate each other but also, they add and even impose their temporal traits on distant regions. Furthermore, both structures house neurons that act as spatially distinct (but overlapped) FP sources whose activation is state, region, and time dependent, making the composition of so-called local FPs highly volatile and strongly site dependent. Since the spatial reach cannot be predicted without source geometry, it is important to assess whether waveforms and temporal motifs arise from a single source; otherwise, those from each of the co-active sources should be sought.

Magnetoencephalography can reveal deep brain network activities linked to memory processes.

Recording from deep neural structures such as hippocampus noninvasively and yet with high temporal resolution remains a major challenge for human neuroscience. Although it has been proposed that deep neuronal activity might be recordable during cognitive tasks using magnetoencephalography (MEG), this remains to be demonstrated as the contribution of deep structures to MEG recordings may be too small to be detected or might be eclipsed by the activity of large-scale neocortical networks. In the present study, we disentangled mesial activity and large-scale networks from the MEG signals thanks to blind source separation (BSS). We then validated the MEG BSS components using intracerebral EEG signals recorded simultaneously in patients during their presurgical evaluation of epilepsy. In the MEG signals obtained during a memory task involving the recognition of old and new images, we identified with BSS a putative mesial component, which was present in all patients and all control subjects. The time course of the component selectively correlated with stereo-electroencephalography signals recorded from hippocampus and rhinal cortex, thus confirming its mesial origin. This finding complements previous studies with epileptic activity and opens new possibilities for using MEG to study deep brain structures in cognition and in brain disorders.

Low-Power Lossless Data Compression for Wireless Brain Electrophysiology.

Wireless electrophysiology opens important possibilities for neuroscience, especially for recording brain activity in more natural contexts, where exploration and interaction are not restricted by the usual tethered devices. The limiting factor is transmission power and, by extension, battery life required for acquiring large amounts of neural electrophysiological data. We present a digital compression algorithm capable of reducing electrophysiological data to less than 65.5% of its original size without distorting the signals, which we tested in vivo in experimental animals. The algorithm is based on a combination of delta compression and Huffman codes with optimizations for neural signals, which allow it to run in small, low-power Field-Programmable Gate Arrays (FPGAs), requiring few hardware resources. With this algorithm, a hardware prototype was created for wireless data transmission using commercially available devices. The power required by the algorithm itself was less than 3 mW, negligible compared to the power saved by reducing the transmission bandwidth requirements. The compression algorithm and its implementation were designed to be device-agnostic. These developments can be used to create a variety of wired and wireless neural electrophysiology acquisition systems with low power and space requirements without the need for complex or expensive specialized hardware.

Brain Network Allostasis after Chronic Alcohol Drinking Is Characterized by Functional Dedifferentiation and Narrowing.

Alcohol use disorder (AUD) causes complex alterations in the brain that are poorly understood. The heterogeneity of drinking patterns and the high incidence of comorbid factors compromise mechanistic investigations in AUD patients. Here we used male Marchigian Sardinian alcohol-preferring (msP) rats, a well established animal model of chronic alcohol drinking, and a combination of longitudinal resting-state fMRI and manganese-enhanced MRI to provide objective measurements of brain connectivity and activity, respectively. We found that 1 month of chronic alcohol drinking changed the correlation between resting-state networks. The change was not homogeneous, resulting in the reorganization of pairwise interactions and a shift in the equilibrium of functional connections. We identified two fundamentally different forms of network reorganization. First is functional dedifferentiation, which is defined as a regional increase in neuronal activity and overall correlation, with a concomitant decrease in preferential connectivity between specific networks. Through this mechanism, occipital cortical areas lost their specific interaction with sensory-insular cortex, striatal, and sensorimotor networks. Second is functional narrowing, which is defined as an increase in neuronal activity and preferential connectivity between specific brain networks. Functional narrowing strengthened the interaction between striatal and prefrontocortical networks, involving the anterior insular, cingulate, orbitofrontal, prelimbic, and infralimbic cortices. Importantly, these two types of alterations persisted after alcohol discontinuation, suggesting that dedifferentiation and functional narrowing rendered persistent network states. Our results support the idea that chronic alcohol drinking, albeit at moderate intoxicating levels, induces an allostatic change in the brain functional connectivity that propagates into early abstinence.SIGNIFICANCE STATEMENT Excessive consumption of alcohol is positioned among the top five risk factors for disease and disability. Despite this priority, the transformations that the nervous system undergoes from an alcohol-naive state to a pathologic alcohol drinking are not well understood. In our study, we use an animal model with proven translational validity to study this transformation longitudinally. The results show that shortly after chronic alcohol consumption there is an increase in redundant activity shared by brain structures, and the specific communication shrinks to a set of pathways. This functional dedifferentiation and narrowing are not reversed immediately after alcohol withdrawal but persist during early abstinence. We causally link chronic alcohol drinking with an early and abstinence-persistent retuning of the functional equilibrium of the brain.

Extracting seizure onset from surface EEG with independent component analysis: Insights from simultaneous scalp and intracerebral EEG.

The success of stereoelectroencephalographic (SEEG) investigations depends crucially on the hypotheses on the putative location of the seizure onset zone. This information is derived from non-invasive data either based on visual analysis or advanced source localization algorithms. While source localization applied to interictal spikes recorded on scalp is the classical method, it does not provide unequivocal information regarding the seizure onset zone. Raw ictal activity contains a mixture of signals originating from several regions of the brain as well as EMG artifacts, hampering direct input to the source localization algorithms. We therefore introduce a methodology that disentangles the various sources contributing to the scalp ictal activity using independent component analysis and uses equivalent current dipole localization as putative locus of ictal sources. We validated the results of our analysis pipeline by performing long-term simultaneous scalp - intracerebral (SEEG) recordings in 14 patients and analyzing the wavelet coherence between the independent component encoding the ictal discharge and the SEEG signals in 8 patients passing the inclusion criteria. Our results show that invasively recorded ictal onset patterns, including low-voltage fast activity, can be captured by the independent component analysis of scalp EEG. The visibility of the ictal activity strongly depends on the depth of the sources. The equivalent current dipole localization can point to the seizure onset zone (SOZ) with an accuracy that can be as high as 10 mm for superficially located sources, that gradually decreases for deeper seizure generators, averaging at 47 mm in the 8 analyzed patients. Independent component analysis is therefore shown to have a promising SOZ localizing value, indicating whether the seizure onset zone is neocortical, and its approximate location, or located in mesial structures. That may contribute to a better crafting of the hypotheses used as basis of the stereo-EEG implantations.

Functional Interactions between Entorhinal Cortical Pathways Modulate Theta Activity in the Hippocampus.

Theta oscillations organize neuronal firing in the hippocampus during context exploration and memory formation. Recently, we have shown that multiple theta rhythms coexist in the hippocampus, reflecting the activity in their afferent regions in CA3 (Schaffer collaterals) and the entorhinal cortex layers II (EC-II, perforant pathway) and III (EC-III, temporoammonic pathway). Frequency and phase coupling between theta rhythms were modulated by the behavioral state, with synchronized theta rhythmicity preferentially occurring in tasks involving memory updating. However, information transmission between theta generators was not investigated. Here, we used source separation techniques to disentangle the current generators recorded in the hippocampus of rats exploring a known environment with or without a novel stimulus. We applied analytical tools based on Granger causality and transfer entropy to investigate linear and non-linear directed interactions, respectively, between the theta activities. Exploration in the novelty condition was associated with increased theta power in the generators with EC origin. We found a significant directed interaction from the Schaffer input over the EC-III input in CA1, and a bidirectional interaction between the inputs in the hippocampus originating in the EC, likely reflecting the connection between layers II and III. During novelty exploration, the influence of the EC-II over the EC-III generator increased, while the Schaffer influence decreased. These results associate the increase in hippocampal theta activity and synchrony during novelty exploration with an increase in the directed functional connectivity from EC-II to EC-III.

Changes in epileptogenicity biomarkers after stereotactic thermocoagulation.

OBJECTIVE: Stereo-electroencephalography (SEEG)-guided radiofrequency thermocoagulation (RF-TC) aims at modifying epileptogenic networks to reduce seizure frequency. High-frequency oscillations (HFOs), spikes, and cross-rate are quantifiable epileptogenic biomarkers. In this study, we sought to evaluate, using SEEG signals recorded before and after thermocoagulation, whether a variation in these markers is related to the therapeutic effect of this procedure and to the outcome of surgery. METHODS: Interictal segments of SEEG signals were analyzed in 38 patients during presurgical evaluation. We used an automatized method to quantify the rate of spikes, rate of HFOs, and cross-rate (a measure combining spikes and HFOs) before and after thermocoagulation. We analyzed the differences both at an individual level with a surrogate approach and at a group level with analysis of variance. We then evaluated the correlation between these variations and the clinical response to RF-TC and to subsequent resective surgery. RESULTS: After thermocoagulation, 19 patients showed a clinical improvement. At the individual level, clinically improved patients more frequently had a reduction in spikes and cross-rate in the epileptogenic zone than patients without clinical improvement (p = .002, p = .02). At a group level, there was a greater decrease of HFOs in epileptogenic and thermocoagulated zones in patients with clinical improvement (p < .05) compared to those with no clinical benefit. Eventually, a significant decrease of all the markers after RF-TC was found in patients with a favorable outcome of resective surgery (spikes, p = .026; HFOs, p = .03; cross-rate, p = .03). SIGNIFICANCE: Quantified changes in the rate of spikes, rate of HFOs, and cross-rate can be observed after thermocoagulation, and the reduction of these markers correlates with a favorable clinical outcome after RF-TC and with successful resective surgery. This may suggest that interictal biomarker modifications after RF-TC can be clinically used to predict the effectiveness of the thermocoagulation procedure and the outcome of resective surgery.

Different theta frameworks coexist in the rat hippocampus and are coordinated during memory-guided and novelty tasks.

Hippocampal firing is organized in theta sequences controlled by internal memory processes and by external sensory cues, but how these computations are coordinated is not fully understood. Although theta activity is commonly studied as a unique coherent oscillation, it is the result of complex interactions between different rhythm generators. Here, by separating hippocampal theta activity in three different current generators, we found epochs with variable theta frequency and phase coupling, suggesting flexible interactions between theta generators. We found that epochs of highly synchronized theta rhythmicity preferentially occurred during behavioral tasks requiring coordination between internal memory representations and incoming sensory information. In addition, we found that gamma oscillations were associated with specific theta generators and the strength of theta-gamma coupling predicted the synchronization between theta generators. We propose a mechanism for segregating or integrating hippocampal computations based on the flexible coordination of different theta frameworks to accommodate the cognitive needs.

In the brain, a vast number of neurons coordinate their activity to support complex cognitive processes. One of the best places to see this in action is the hippocampus, a brain structure with a key role in memory and navigation. The hippocampus shows waves of electrical activity, which represent the synchronized firing of large numbers of neurons. The hippocampus can generate multiple rhythms at once. The two main rhythms are theta and gamma. Theta waves are slow, with a frequency of about 8 Hertz. Gamma waves are faster with a frequency of up to 120 Hertz or even more. Theta waves are always present in the brains of freely moving animals, whereas gamma waves occur in brief bursts. These bursts usually correspond to a particular point on the theta wave. One burst may occur just before each peak of the theta wave, for example, whereas another burst may occur just after the peak. This separation enables individual bursts of gamma to carry different messages without them becoming mixed up. This is similar to how radio stations broadcast their signals at different carrier frequencies to avoid interference. By recording hippocampal activity in rats exploring a maze, Lopez-Madrona et al. now show that the hippocampus has not one, but three generators of theta waves. Having three sources of theta, each of which can be synchronized with gamma, provides a more versatile system for encoding and sending information. It also means that the three theta generators can vary the degree to which they coordinate their firing. This helps the brain combine or separate streams of information as required. By working together to create a single theta rhythm, for example, the three theta generators can help animals combine information stored in memory with incoming sensory input. How the coordination of theta rhythms in the hippocampus influences the activity of other brain regions involved in learning and memory remains unclear. However, uncoupling of theta and gamma waves seems to be an early sign of Alzheimer's disease and can also be seen in the brains of people with schizophrenia and other psychiatric disorders. Understanding how this process occurs could provide clues to the origin of these disorders.

Inferring correlations associated to causal interactions in brain signals using autoregressive models.

The specific connectivity of a neuronal network is reflected in the dynamics of the signals recorded on its nodes. The analysis of how the activity in one node predicts the behaviour of another gives the directionality in their relationship. However, each node is composed of many different elements which define the properties of the links. For instance, excitatory and inhibitory neuronal subtypes determine the functionality of the connection. Classic indexes such as the Granger causality (GC) quantifies these interactions, but they do not infer into the mechanism behind them. Here, we introduce an extension of the well-known GC that analyses the correlation associated to the specific influence that a transmitter node has over the receiver. This way, the G-causal link has a positive or negative effect if the predicted activity follows directly or inversely, respectively, the dynamics of the sender. The method is validated in a neuronal population model, testing the paradigm that excitatory and inhibitory neurons have a differential effect in the connectivity. Our approach correctly infers the positive or negative coupling produced by different types of neurons. Our results suggest that the proposed approach provides additional information on the characterization of G-causal connections, which is potentially relevant when it comes to understanding interactions in the brain circuits.

Slow-Wave Activity in the S1HL Cortex Is Contributed by Different Layer-Specific Field Potential Sources during Development.

Spontaneous correlated activity in cortical columns is critical for postnatal circuit refinement. We used spatial discrimination techniques to explore the late maturation of synaptic pathways through the laminar distribution of the field potential (FP) generators underlying spontaneous and evoked activities of the S1HL cortex in juvenile (P14-P16) and adult anesthetized rats. Juveniles exhibit an intermittent FP pattern resembling Up/Down states in adults, but with much reduced power and different laminar distribution. Whereas FPs in active periods are dominated by a layer VI generator in juveniles, in adults a developing multipart generator takes over, displaying current sinks in middle layers (III-V). The blockade of excitatory transmission in upper and middle layers of adults recovered the juvenile-like FP profiles. In addition to the layer VI generator, a gamma-specific generator in supragranular layers was the same in both age groups. While searching for dynamical coupling among generators in juveniles we found significant cross-correlation in ∼one-half of the tested pairs, whereas excessive coherence hindered their efficient separation in adults. Also, potentials evoked by tactile and electrical stimuli showed different short-latency dipoles between the two age groups, and the juveniles lacked the characteristic long latency UP state currents in middle layers. In addition, the mean firing rate of neurons was lower in juveniles. Thus, cortical FPs originate from different intra-columnar segments as they become active postnatally. We suggest that although some cortical segments are active early postnatally, a functional sensory-motor control relies on a delayed maturation and network integration of synaptic connections in middle layers.SIGNIFICANCE STATEMENT Early postnatal activity in the rodent cortex is mostly endogenous, whereas it becomes driven by peripheral input at later stages. The precise schedule for the maturation of synaptic pathways is largely unknown. We explored this in the somatosensory hindlimb cortex at an age when animals begin to use their limbs by uncovering the laminar distribution of the field potential generators underlying the dominant delta waves in juveniles and adults. Our results suggest that field potentials are mostly generated by a pathway in deep layers, whereas other pathways mature later in middle layers and take over in adults. We suggest that a functional sensory-motor control relies on a delayed maturation and network integration of synaptic connections in middle layers.

On the role of the entorhinal cortex in the effective connectivity of the hippocampal formation.

Inferring effective connectivity from neurophysiological data is a challenging task. In particular, only a finite (and usually small) number of sites are simultaneously recorded, while the response of one of these sites can be influenced by other sites that are not being recorded. In the hippocampal formation, for instance, the connections between areas CA1-CA3, the dentate gyrus (DG), and the entorhinal cortex (EC) are well established. However, little is known about the relations within the EC layers, which might strongly affect the resulting effective connectivity estimations. In this work, we build excitatory/inhibitory neuronal populations representing the four areas CA1, CA3, the DG, and the EC and fix their connectivities. We model the EC by three layers (LII, LIII, and LV) and assume any possible connection between them. Our results, based on Granger Causality (GC) and Partial Transfer Entropy (PTE) measurements, reveal that the estimation of effective connectivity in the hippocampus strongly depends on the connectivities between EC layers. Moreover, we find, for certain EC configurations, very different results when comparing GC and PTE measurements. We further demonstrate that causal links can be robustly inferred regardless of the excitatory or inhibitory nature of the connection, adding complexity to their interpretation. Overall, our work highlights the importance of a careful analysis of the connectivity methods to prevent unrealistic conclusions when only partial information about the experimental system is available, as usually happens in brain networks. Our results suggest that the combination of causality measures with neuronal modeling based on precise neuroanatomical tracing may provide a powerful framework to disambiguate causal interactions in the brain.