Cortical Network Theory
Functional connectivity and graph theory in EEG
From signals to networks
Classical electroencephalography reads each channel as a window onto local cortex. Network neuroscience reframes the same recording as a system of interacting nodes: each electrode (or, better, each reconstructed cortical source) is a node, and a quantitative measure of statistical dependence between two nodes' time series becomes an edge. The collection of nodes and edges is a graph, and the mathematics of graph theory - originally developed for transport, social, and computer networks - can then be brought to bear on the brain. This shift matters because many neurological diseases are increasingly understood not as focal lesions but as disorders of distributed coordination: the lesion may be local, but the dysfunction is networked, and the most disabling consequences often arise from how a local abnormality reverberates through the wider system.
Two conceptual layers must be kept separate. Anatomical (structural) connectivity is the physical wiring - the white-matter tracts and cortico-cortical fibers - which scalp EEG cannot measure directly (it is the province of diffusion MRI tractography). What EEG estimates is the dynamic coordination running over that wiring, and this comes in two flavors that are routinely confused: functional connectivity and effective connectivity. Conflating them, or treating either as a direct readout of anatomy, is among the most common errors in the literature, so we treat them deliberately and then turn to the equally important question of what these estimates can and cannot be trusted to mean.
Functional versus effective connectivity
Functional connectivity is undirected statistical dependence: a symmetric measure of how much two signals covary, oscillate together, or share information, with no claim about who drives whom. Common EEG estimators include coherence and the imaginary part of coherency (which discards the zero-lag component to suppress instantaneous volume-conduction artifact), the phase-locking value (PLV), and the phase-lag index (PLI) together with its weighted variant (wPLI), which quantify the consistency of non-zero phase differences while being robust to common-source leakage. Amplitude-based measures such as the amplitude envelope correlation (AEC) (ideally with orthogonalization to remove zero-lag leakage) and information-theoretic mutual information round out the toolbox. The defining property is symmetry: by construction, connectivity from A to B equals connectivity from B to A.
Effective connectivity is directed causal influence: an asymmetric measure of how activity in one node shapes activity in another. Estimators include Granger causality and its spectral generalizations - the directed transfer function (DTF) and partial directed coherence (PDC) built on multivariate autoregressive (MVAR) models - and transfer entropy, a model-free, information-theoretic measure of directed information flow that, unlike linear Granger methods, captures nonlinear interactions. Model-based dynamic causal modeling (DCM) goes further still, fitting a generative biophysical (neural-mass) model to the data and using Bayesian inference to estimate the coupling parameters that best explain the observations, including their direction and even their excitatory or inhibitory character. Effective connectivity is what we want when the clinical question is directional - for example, which node leads at seizure onset, or whether top-down feedback is intact.
A single source seen by two scalp electrodes creates spurious zero-lag correlation that masquerades as connectivity - and the common average or linked-ears reference adds a further shared signal to every channel, inflating apparent coupling network-wide. This is why volume-conduction-robust metrics (imaginary coherency, PLI, wPLI, orthogonalized AEC) are preferred for scalp data, why source-space analysis is increasingly standard, and why any raw coherence finding at or near zero phase lag should be regarded with suspicion until artifact and reference effects are excluded.
| Property | Functional connectivity | Effective connectivity |
|---|---|---|
| Directionality | Undirected (symmetric) | Directed (asymmetric) |
| Question answered | Are A and B coordinated? | Does A drive B? |
| Typical EEG estimators | Coherence, imaginary coherency, PLI/wPLI, PLV, AEC, mutual information | Granger causality, DTF, PDC, transfer entropy, DCM |
| Model dependence | Largely model-free | Often model-based (autoregressive or generative) |
| Main pitfall | Volume-conduction / reference-driven zero-lag spuriousness | Sensitivity to model order, sampling, stationarity, and unobserved nodes |
The graph-theoretic vocabulary
Once edges are estimated and (usually) thresholded into a graph - whether binary or weighted - a compact set of metrics characterizes its topology. Degree is the number of edges incident on a node, its connectedness; in weighted graphs the analogous quantity is node strength, the sum of its edge weights. The degree distribution across the network reveals whether connectivity is homogeneous or dominated by a few richly connected nodes. Clustering coefficient quantifies local integration: the probability that a node's neighbours are themselves interconnected (the fraction of closed triangles), indexing segregated, specialized processing. Characteristic path length captures global integration: the average number of edges on the shortest path between node pairs, with shorter paths permitting more efficient communication across the whole network. Because path length is undefined for disconnected graphs, the closely related global efficiency (the average of the inverse shortest path lengths) is often preferred, and a complementary local efficiency can be computed within each node's neighbourhood.
The interplay of these two scales yields the concept of small-worldness. A small-world network combines high clustering (like a regular lattice) with short path length (like a random graph), achieving an efficient balance between segregation and integration. It is conventionally indexed by a ratio (often written sigma) that normalizes the observed clustering and path length against those of degree-matched random graphs; values meaningfully greater than one indicate small-world organization. The healthy brain consistently exhibits this regime across imaging modalities, and deviations from it - drift toward either a rigid, over-regular lattice or a randomized, inefficient topology - recur as candidate signatures of disease. A caveat worth carrying: the sigma metric is sensitive to network size, density, and the choice of null model, so apparent small-world differences between groups can be artifacts of differing edge density rather than genuine topological reorganization; many authors now report the more robust small-world propensity or compare networks at matched density.
Finally, hubs are high-degree, high-centrality nodes that disproportionately broker communication. Betweenness centrality counts how often a node lies on shortest paths between other nodes, and hubs with high betweenness are the network's critical relays; eigenvector centrality and related measures instead reward connection to other well-connected nodes. A further organizing principle is the rich club: the tendency of high-degree hubs to be densely interconnected with one another, forming a central backbone that integrates information across otherwise segregated modules. Hubs and the rich club confer efficiency but also vulnerability: targeted disruption of a hub damages global integration far more than removing a peripheral node, and rich-club regions are disproportionately implicated across a wide range of brain disorders. This economy-versus-robustness trade-off - the metabolic and wiring cost of hubs against the efficiency they buy and the fragility they introduce - is central to interpreting why certain lesions or seizure-onset zones have outsized network consequences.
Two opposing demands shape every brain network: segregation (local, specialized processing, indexed by clustering and local efficiency) and integration (rapid global communication, indexed by short path length / high global efficiency). Small-world topology, with a rich-club backbone linking its hubs, is the brain's solution to satisfying both at once, and a great deal of network pathology can be read as a shift in this balance.
Modules, the connectome, and dynamic connectivity
Two further ideas complete the working vocabulary. The first is community structure (modularity): real brain networks are not uniformly wired but partition into modules - subsets of nodes more densely connected among themselves than with the rest of the graph - that correspond roughly to functional systems (visual, sensorimotor, default-mode, frontoparietal control, and so on). Modularity confers robustness and supports specialized processing, and the nodes that link modules to one another, the connector hubs, are functionally distinct from provincial hubs that are richly connected only within their own module. This distinction matters clinically because connector hubs and the rich club are disproportionately the regions whose disruption fragments global communication, and they recur across the imaging literature as sites of selective vulnerability in disease. The second idea is simply terminological hygiene: the comprehensive map of these elements - nodes, edges, modules, hubs - is the connectome, and EEG contributes a functional, temporally resolved connectome that complements the structural connectome from diffusion MRI and the slower functional connectome from resting-state fMRI.
A worked conceptual example fixes why directed metrics matter. Suppose three sources participate in an ictal discharge and a clinician must decide which to target. A functional (undirected) analysis might show all three strongly coupled - true but useless, because coupling alone does not say who leads. A directed analysis using partial directed coherence or Granger causality may instead reveal that one node has high out-degree and the other two are predominantly receivers, identifying the driver of the network. The logic generalizes: undirected metrics describe the shape of coordination, directed metrics describe its flow, and clinical questions about seizure onset, surgical targets, or the integrity of top-down control are nearly always questions about flow. This is also why directed analyses are most trustworthy on intracranial data, where the volume-conduction confound that plagues scalp directionality is greatly attenuated and short conduction delays can be resolved.
Finally, connectivity is increasingly recognized as time-varying rather than static. The brain does not hold a fixed graph; it reconfigures its coupling on the timescale of cognition and of the transition into seizure. Dynamic (time-resolved) connectivity methods - sliding-window estimates, adaptive (Kalman-filtered) autoregressive models, hidden-Markov models that infer recurring connectivity states, and EEG microstate analysis that segments the scalp topography into a small alphabet of quasi-stable states lasting tens of milliseconds - aim to capture this flux. They reveal, for instance, that integration and segregation wax and wane spontaneously, that anaesthesia narrows the repertoire of visited network states, and that the approach to a seizure can be heralded by stereotyped shifts in connectivity state minutes before clinical onset. The price is statistical: shorter windows mean noisier estimates and more opportunities for the artifacts catalogued in the next section to masquerade as dynamics, so time-resolved findings demand even more methodological discipline than static ones.
Failure modes and the limits of connectivity estimation
Before celebrating these metrics it is essential to confront how easily they mislead, because a sophisticated reader is defined as much by knowing the traps as by knowing the tools. The deepest problem is volume conduction and source leakage. A single neural source projects to many electrodes simultaneously and instantaneously, manufacturing zero-lag correlations wherever it is seen by more than one sensor; even after projecting to source space, the inverse solution smears each true source across neighbouring reconstructed locations (the point-spread of the inverse operator), so that spatially adjacent sources show spurious coupling - the source-leakage problem. Leakage-robust metrics (imaginary coherency, PLI/wPLI, orthogonalized AEC) suppress the strictly zero-lag component, but they do so at a real cost: they are blind to genuine near-zero-lag synchronization, which is physiologically common, and they remain only partially protected against leakage that is not exactly at zero phase.
The reference electrode is a second, under-appreciated confounder unique to EEG. Any non-quiet reference injects its own signal into every channel with opposite sign, inflating apparent synchrony across the whole montage; the common average reference distributes a mixture of all channels into each, and a linked-ears or vertex reference imposes its activity globally. Different references can therefore produce qualitatively different connectivity graphs from identical data, and reference-free or source-space approaches (or the surface Laplacian / current source density) are increasingly required for credible results. Beyond these field-spread issues lie a host of estimation pitfalls: non-stationarity (most connectivity estimators assume the signal statistics are stable over the analysis epoch, which EEG routinely violates), sensitivity of MVAR-based directed measures to model order and to down-sampling and filtering, distortion of phase estimates by filtering, the field of view problem (unobserved or unmeasured nodes can create the illusion of a direct link between two recorded nodes that are in fact both driven by a hidden third), and the strong dependence of every graph metric on edge density, thresholding scheme, node count, epoch length, and signal-to-noise ratio. The unglamorous but correct conclusion is that two studies of the same disease can disagree purely because of pipeline differences, and that connectivity findings must be defended with matched null models, density-matched comparisons, leakage-robust metrics, and, wherever possible, replication.
1) Which estimator, and is it leakage-robust? 2) Electrode space or source space? 3) What reference, and was a reference-independent transform used? 4) Were comparisons made at matched edge density against degree-matched null models? 5) Was stationarity respected (or time-varying methods used)? 6) Could an unobserved common driver or muscle/eye artifact explain the edges? If these are unanswered, treat the topology as hypothesis-generating, not established.
Network disruption in disease
In epilepsy, network analysis has matured from a research curiosity into a framework that reshapes how seizures are conceived. The interictal epileptic network often shows increased local clustering and a drift toward regularization - a move away from the small-world optimum toward a more lattice-like, hyper-synchronizable topology that may favour seizure recruitment, though the literature is not unanimous and direction of change varies with frequency band, state, and method. The epileptogenic zone frequently behaves as a pathological hub, with abnormally high out-degree in directed (effective-connectivity) analyses at seizure onset; DTF and PDC have been used to identify the leading node of an ictal discharge, and such directed metrics - increasingly applied to intracranial (stereo-EEG) recordings, where volume conduction is far less severe - inform presurgical hypotheses about where to resect or where to place responsive-neurostimulation leads. Crucially, the network view explains a long-standing surgical puzzle: removing the apparent focus sometimes fails because the driver of the seizure network lies elsewhere, because critical hubs are spared, or because the network reorganizes around the resection. A growing body of work also examines interictal markers of the epileptogenic network, including high-frequency oscillations and the network's response to single-pulse electrical stimulation (cortico-cortical evoked potentials), as complements to ictal mapping.
In diffuse encephalopathy - metabolic, toxic, septic, anoxic, or related to sedation - the picture is broadly one of disintegration of global topology. As consciousness declines, EEG network studies frequently report falling global efficiency, lengthening path length, and a breakdown of the integrated hub and rich-club structure that supports normal cognition, alongside the familiar visual hallmark of background slowing. This dovetails with influential theories of consciousness - notably integrated information theory and the global-workspace framework - that emphasize the brain's capacity for integrated yet differentiated activity: anaesthesia and coma tend to be accompanied by a collapse of long-range integration and a retreat toward fragmented, locally clustered, or randomized networks. Quantitative network metrics, and complexity measures such as the perturbational complexity index derived from TMS-EEG, are accordingly being explored as graded markers of encephalopathy severity, of the level of consciousness, and of recovery potential, with the aspiration of being more sensitive and more reproducible than visual grading of background alone. These remain active research tools rather than validated bedside diagnostics, and their performance is bounded by exactly the methodological fragilities catalogued above.
The network framing extends beyond these two paradigms, though the evidence thins as one moves from acute to chronic disease and should be held more loosely. In dementia, and most studied in Alzheimer disease, the recurring EEG description is of a network drifting toward randomization with loss of small-world organization and reduced long-range, especially alpha-band, synchronization - a picture often summarized as functional disconnection, and consistent with the broader hypothesis that Alzheimer disease preferentially attacks highly connected hub regions (the disrupted-connectome or hub-vulnerability view). In disorders of consciousness following severe brain injury, graph and complexity measures derived from EEG are being investigated to help distinguish the unresponsive-wakefulness (vegetative) state from minimally conscious states and to detect covert awareness - residual, command-following cognition in behaviourally unresponsive patients - a clinically and ethically weighty application now the subject of active multicentre work. Across all of these, the same caution recurs: the direction and magnitude of reported network change depend heavily on the estimator, the band, the state, and the analysis pipeline, so convergence across methods and independent replication matter more than any single striking result.
Network metrics are genuinely powerful but methodologically fragile: results depend on the connectivity estimator, the thresholding and density choices, the number of nodes, epoch length, montage, reference, and artifact handling. Treat reported network changes as hypothesis-generating until reproduced with leakage-robust, source-space or intracranial methods, matched null models, and density-matched comparisons. The strongest current clinical traction is in presurgical seizure-network mapping with intracranial directed connectivity, where volume conduction is least problematic.
The trajectory of the field is clear. As source-localization and high-density recording improve, EEG connectivity is moving from electrode space toward anatomically interpretable source networks, and from static snapshots toward dynamic (time-varying) connectivity - including hidden-Markov and microstate approaches - that tracks how topology reconfigures across cognitive states and across the transition into seizure. Increasingly these EEG-derived networks are fused with structural connectomes in personalized whole-brain models to predict seizure spread and stimulation effects. For the clinician, the durable take-home is conceptual rather than computational: many of the disorders we are accustomed to localizing to a spot are better understood as diseases of a network, and the same EEG that yields a focus can, with disciplined and skeptical analysis, also yield the topology in which that focus is embedded.
1. A study reports strong zero-phase-lag coherence between two adjacent scalp electrodes and concludes the underlying regions are functionally connected. What is the most important objection?
2. Which combination of graph metrics best characterizes a small-world network, the topology typical of the healthy brain?
3. Two recorded EEG nodes show a strong apparent directed influence from A to B, but a hidden region C actually drives both. What estimation problem does this illustrate?