Structural, geometric and genetic factors predict interregional brain connectivity patterns probed by electrocorticography

Identifying epileptogenic abnormality by decomposing intracranial EEG and MEG power spectra

BACKGROUND

Accurate identification of abnormal electroencephalographic (EEG) activity is pivotal for diagnosing and treating epilepsy. Recent studies indicate that decomposing brain activity into periodic (oscillatory) and aperiodic (trend across all frequencies) components can illuminate the drivers of spectral activity changes.

NEW METHODS

We analysed intracranial EEG (iEEG) data from 234 subjects, creating a normative map. This map was compared to a cohort of 63 patients with refractory focal epilepsy under consideration for neurosurgery. The normative map was computed using three approaches: (i) relative complete band power, (ii) relative band power with the aperiodic component removed, and (iii) the aperiodic exponent. Abnormalities were calculated for each approach in the patient cohort. We evaluated the spatial profiles, assessed their ability to localize abnormalities, and replicated the findings using magnetoencephalography (MEG).

RESULTS

Normative maps of relative complete band power and relative periodic band power exhibited similar spatial profiles, while the aperiodic normative map revealed higher exponent values in the temporal lobe. Abnormalities estimated through complete band power effectively distinguished between good and bad outcome patients. Combining periodic and aperiodic abnormalities enhanced performance, like the complete band power approach.

COMPARISON WITH EXISTING METHODS AND CONCLUSIONS

Sparing cerebral tissue with abnormalities in both periodic and aperiodic activity may result in poor surgical outcomes. Both periodic and aperiodic components do not carry sufficient information in isolation. The relative complete band power solution proved to be the most reliable method for this purpose. Future studies could investigate how cerebral location or pathology influences periodic or aperiodic abnormalities.

Utility of intracranial EEG networks depends on re-referencing and connectivity choice

Studies of intracranial EEG networks have been used to reveal seizure generators in patients with drug-resistant epilepsy. Intracranial EEG is implanted to capture the epileptic network, the collection of brain tissue that forms a substrate for seizures to start and spread. Interictal intracranial EEG measures brain activity at baseline, and networks computed during this state can reveal aberrant brain tissue without requiring seizure recordings. Intracranial EEG network analyses require choosing a reference and applying statistical measures of functional connectivity. Approaches to these technical choices vary widely across studies, and the impact of these technical choices on downstream analyses is poorly understood. Our objective was to examine the effects of different re-referencing and connectivity approaches on connectivity results and on the ability to lateralize the seizure onset zone in patients with drug-resistant epilepsy. We applied 48 pre-processing pipelines to a cohort of 125 patients with drug-resistant epilepsy recorded with interictal intracranial EEG across two epilepsy centres to generate intracranial EEG functional connectivity networks. Twenty-four functional connectivity measures across time and frequency domains were applied in combination with common average re-referencing or bipolar re-referencing. We applied an unsupervised clustering algorithm to identify groups of pre-processing pipelines. We subjected each pre-processing approach to three quality tests: (i) the introduction of spurious correlations; (ii) robustness to incomplete spatial sampling; and (iii) the ability to lateralize the clinician-defined seizure onset zone. Three groups of similar pre-processing pipelines emerged: common average re-referencing pipelines, bipolar re-referencing pipelines and relative entropy-based connectivity pipelines. Relative entropy and common average re-referencing networks were more robust to incomplete electrode sampling than bipolar re-referencing and other connectivity methods (Friedman test, Dunn-Šidák test P < 0.0001). Bipolar re-referencing reduced spurious correlations at non-adjacent channels better than common average re-referencing (Δ mean from machine ref = -0.36 versus -0.22) and worse in adjacent channels (Δ mean from machine ref = -0.14 versus -0.40). Relative entropy-based network measures lateralized the seizure onset hemisphere better than other measures in patients with temporal lobe epilepsy (Benjamini-Hochberg-corrected P < 0.05, Cohen's d: 0.60-0.76). Finally, we present an interface where users can rapidly evaluate intracranial EEG pre-processing choices to select the optimal pre-processing methods tailored to specific research questions. The choice of pre-processing methods affects downstream network analyses. Choosing a single method among highly correlated approaches can reduce redundancy in processing. Relative entropy outperforms other connectivity methods in multiple quality tests. We present a method and interface for researchers to optimize their pre-processing methods for deriving intracranial EEG brain networks.

Injectable Fluorescent Neural Interfaces for Cell-Specific Stimulating and Imaging

Building on current explorations in chronic optical neural interfaces, it is essential to address the risk of photothermal damage in traditional optogenetics. By focusing on calcium fluorescence for imaging rather than stimulation, injectable fluorescent neural interfaces significantly minimize photothermal damage and improve the accuracy of neuronal imaging. Key advancements including the use of injectable microelectronics for targeted electrical stimulation and their integration with cell-specific genetically encoded calcium indicators have been discussed. These injectable electronics that allow for post-treatment retrieval offer a minimally invasive solution, enhancing both usability and reliability. Furthermore, the integration of genetically encoded fluorescent calcium indicators with injectable bioelectronics enables precise neuronal recording and imaging of individual neurons. This shift not only minimizes risks such as photothermal conversion but also boosts safety, specificity, and effectiveness of neural imaging. Embracing these advancements represents a significant leap forward in biomedical engineering and neuroscience, paving the way for advanced brain-machine interfaces.

Group-common and individual-specific effects of structure-function coupling in human brain networks with graph neural networks

The human cerebral cortex is organized into functionally segregated but synchronized regions bridged by the structural connectivity of white matter pathways. While structure-function coupling has been implicated in cognitive development and neuropsychiatric disorders, studies yield inconsistent findings. The extent to which the structure-function coupling reflects reliable individual differences or primarily group-common characteristics remains unclear, at both the global and regional brain levels. By leveraging two independent, high-quality datasets, we found that the graph neural network accurately predicted unseen individuals' functional connectivity from structural connectivity, reflecting a strong structure-function coupling. This coupling was primarily driven by network topology and was substantially stronger than that of the linear models. Moreover, we observed that structure-function coupling was dominated by group-common effects, with subtle yet significant individual-specific effects. The regional group and individual effects of coupling were hierarchically organized across the cortex along a sensorimotor-association axis, with lower group and higher individual effects in association cortices. These findings emphasize the importance of considering both group and individual effects in understanding cortical structure-function coupling, suggesting insights into interpreting individual differences of the coupling and informing connectivity-guided therapeutics.

Discrete tactile feature comparison subprocess in human brain during a decision-making process

From sensory input to motor action, encoded sensory features flow sequentially along cortical networks for decision-making. Despite numerous studies probing the decision-making process, the subprocess that compares encoded sensory features before making a decision has not been fully elucidated in humans. In this study, we investigated sensory feature comparison by presenting two different tasks (a discrimination task, in which participants made decisions by comparing two sequential tactile stimuli; and a detection task, in which participants responded to the second tactile stimulus in two sequential stimuli) to epilepsy patients while recording electrocorticography (ECoG). By comparing tactile-specific gamma band (30-200 Hz) power between the two tasks, the decision-making process was divided into three subprocesses-categorization, comparison, and decision-consistent with a previous study (Heekeren et al., 2004). These subprocesses occurred sequentially in the dorsolateral prefrontal cortex, premotor cortex, secondary somatosensory cortex, and parietal lobe. Gamma power showed two different patterns of correlation with response time. In the inferior parietal lobule (IPL), there was a negative correlation. This means that as gamma power increased, response time decreased. In the secondary somatosensory cortex (S2), there was a positive correlation. Here, as gamma power increased, response time also increased. These results indicate that the IPL and S2 encode tactile feature comparison differently. Our connectivity analysis showed that the S2 transmitted tactile information to the IPL. Our findings suggest that multiple areas in the parietal lobe encode sensory feature comparison differently before making a decision.

Modular subgraphs in large-scale connectomes underpin spontaneous co-fluctuation events in mouse and human brains

Previous studies have adopted an edge-centric framework to study fine-scale network dynamics in human fMRI. To date, however, no studies have applied this framework to data collected from model organisms. Here, we analyze structural and functional imaging data from lightly anesthetized mice through an edge-centric lens. We find evidence of "bursty" dynamics and events - brief periods of high-amplitude network connectivity. Further, we show that on a per-frame basis events best explain static FC and can be divided into a series of hierarchically-related clusters. The co-fluctuation patterns associated with each cluster centroid link distinct anatomical areas and largely adhere to the boundaries of algorithmically detected functional brain systems. We then investigate the anatomical connectivity undergirding high-amplitude co-fluctuation patterns. We find that events induce modular bipartitions of the anatomical network of inter-areal axonal projections. Finally, we replicate these same findings in a human imaging dataset. In summary, this report recapitulates in a model organism many of the same phenomena observed in previously edge-centric analyses of human imaging data. However, unlike human subjects, the murine nervous system is amenable to invasive experimental perturbations. Thus, this study sets the stage for future investigation into the causal origins of fine-scale brain dynamics and high-amplitude co-fluctuations. Moreover, the cross-species consistency of the reported findings enhances the likelihood of future translation.

The involvement of the human prefrontal cortex in the emergence of visual awareness

Exploring the neural mechanisms of awareness is a fundamental task of cognitive neuroscience. There is an ongoing dispute regarding the role of the prefrontal cortex (PFC) in the emergence of awareness, which is partially raised by the confound between report- and awareness-related activity. To address this problem, we designed a visual awareness task that can minimize report-related motor confounding. Our results show that saccadic latency is significantly shorter in the aware trials than in the unaware trials. Local field potential (LFP) data from six patients consistently show early (200-300ms) awareness-related activity in the PFC, including event-related potential and high-gamma activity. Moreover, the awareness state can be reliably decoded by the neural activity in the PFC since the early stage, and the neural pattern is dynamically changed rather than being stable during the representation of awareness. Furthermore, the enhancement of dynamic functional connectivity, through the phase modulation at low frequency, between the PFC and other brain regions in the early stage of the awareness trials may explain the mechanism of conscious access. These results indicate that the PFC is critically involved in the emergence of awareness.

Macroscopic resting-state brain dynamics are best described by linear models

It is typically assumed that large networks of neurons exhibit a large repertoire of nonlinear behaviours. Here we challenge this assumption by leveraging mathematical models derived from measurements of local field potentials via intracranial electroencephalography and of whole-brain blood-oxygen-level-dependent brain activity via functional magnetic resonance imaging. We used state-of-the-art linear and nonlinear families of models to describe spontaneous resting-state activity of 700 participants in the Human Connectome Project and 122 participants in the Restoring Active Memory project. We found that linear autoregressive models provide the best fit across both data types and three performance metrics: predictive power, computational complexity and the extent of the residual dynamics unexplained by the model. To explain this observation, we show that microscopic nonlinear dynamics can be counteracted or masked by four factors associated with macroscopic dynamics: averaging over space and over time, which are inherent to aggregated macroscopic brain activity, and observation noise and limited data samples, which stem from technological limitations. We therefore argue that easier-to-interpret linear models can faithfully describe macroscopic brain dynamics during resting-state conditions.

Eigenmode-based approach reveals a decline in brain structure-function liberality across the human lifespan

While brain function is supported and constrained by the underlying structure, the connectome-based link estimated by current approaches is either relatively moderate or accompanied by high model complexity, with the essential principles underlying structure-function coupling remaining elusive. Here, by proposing a mapping method based on network eigendecomposition, we present a concise and strong correspondence between structure and function. We show that the explanation of functional connectivity can be significantly improved by incorporating interactions between different structural eigenmodes. We also demonstrate the pronounced advantage of the present mapping in capturing individual-specific information with simple implementation. Applying our methodology to the human lifespan, we find that functional diversity decreases with age, with functional interactions increasingly dominated by the leading functional mode. We also find that structure-function liberality weakens with age, which is driven by the decreases in functional components that are less constrained by anatomy, while the magnitude of structure-aligned components is preserved. Overall, our work enhances the understanding of structure-function coupling from a collective, connectome-oriented perspective and promotes a more refined identification of functional portions relevant to human aging, holding great potential for mechanistic insights into individual differences associated with cognition, development, and neurological disorders.

Controversies and progress on standardization of large-scale brain network nomenclature

Progress in scientific disciplines is accompanied by standardization of terminology. Network neuroscience, at the level of macroscale organization of the brain, is beginning to confront the challenges associated with developing a taxonomy of its fundamental explanatory constructs. The Workgroup for HArmonized Taxonomy of NETworks (WHATNET) was formed in 2020 as an Organization for Human Brain Mapping (OHBM)-endorsed best practices committee to provide recommendations on points of consensus, identify open questions, and highlight areas of ongoing debate in the service of moving the field toward standardized reporting of network neuroscience results. The committee conducted a survey to catalog current practices in large-scale brain network nomenclature. A few well-known network names (e.g., default mode network) dominated responses to the survey, and a number of illuminating points of disagreement emerged. We summarize survey results and provide initial considerations and recommendations from the workgroup. This perspective piece includes a selective review of challenges to this enterprise, including (1) network scale, resolution, and hierarchies; (2) interindividual variability of networks; (3) dynamics and nonstationarity of networks; (4) consideration of network affiliations of subcortical structures; and (5) consideration of multimodal information. We close with minimal reporting guidelines for the cognitive and network neuroscience communities to adopt.

Soft and stretchable organic bioelectronics for continuous intraoperative neurophysiological monitoring during microsurgery

In microneurosurgery, it is crucial to maintain the structural and functional integrity of the nerve through continuous intraoperative identification of neural anatomy. To this end, here we report the development of a translatable system leveraging soft and stretchable organic-electronic materials for continuous intraoperative neurophysiological monitoring. The system uses conducting polymer electrodes with low impedance and low modulus to record near-field action potentials continuously during microsurgeries, offers higher signal-to-noise ratios and reduced invasiveness when compared with handheld clinical probes for intraoperative neurophysiological monitoring and can be multiplexed, allowing for the precise localization of the target nerve in the absence of anatomical landmarks. Compared with commercial metal electrodes, the neurophysiological monitoring system allowed for enhanced post-operative prognoses after tumour-resection surgeries in rats. Continuous recording of near-field action potentials during microsurgeries may allow for the precise identification of neural anatomy through the entire procedure.

Intracranial EEG Structure-Function Coupling and Seizure Outcomes After Epilepsy Surgery

BACKGROUND AND OBJECTIVES

Surgery is an effective treatment for drug-resistant epilepsy, which modifies the brain's structure and networks to regulate seizure activity. Our objective was to examine the relationship between brain structure and function to determine the extent to which this relationship affects the success of the surgery in controlling seizures. We hypothesized that a stronger association between brain structure and function would lead to improved seizure control after surgery.

METHODS

We constructed functional and structural brain networks in patients with drug-resistant focal epilepsy by using presurgery functional data from intracranial EEG (iEEG) recordings, presurgery and postsurgery structural data from T1-weighted MRI, and presurgery diffusion-weighted MRI. We quantified the relationship (coupling) between structural and functional connectivity by using the Spearman rank correlation and analyzed this structure-function coupling at 2 spatial scales: (1) global iEEG network level and (2) individual iEEG electrode contacts using virtual surgeries. We retrospectively predicted postoperative seizure freedom by incorporating the structure-function connectivity coupling metrics and routine clinical variables into a cross-validated predictive model.

RESULTS

We conducted a retrospective analysis on data from 39 patients who met our inclusion criteria. Brain areas implanted with iEEG electrodes had stronger structure-function coupling in seizure-free patients compared with those with seizure recurrence (p = 0.002, d = 0.76, area under the receiver operating characteristic curve [AUC] = 0.78 [95% CI 0.62-0.93]). Virtual surgeries on brain areas that resulted in stronger structure-function coupling of the remaining network were associated with seizure-free outcomes (p = 0.007, d = 0.96, AUC = 0.73 [95% CI 0.58-0.89]). The combination of global and local structure-function coupling measures accurately predicted seizure outcomes with a cross-validated AUC of 0.81 (95% CI 0.67-0.94). These measures were complementary to other clinical variables and, when included for prediction, resulted in a cross-validated AUC of 0.91 (95% CI 0.82-1.0), accuracy of 92%, sensitivity of 93%, and specificity of 91%.

DISCUSSION

Our study showed that the strength of structure-function connectivity coupling may play a crucial role in determining the success of epilepsy surgery. By quantitatively incorporating structure-function coupling measures and standard-of-care clinical variables into presurgical evaluations, we may be able to better localize epileptogenic tissue and select patients for epilepsy surgery.

CLASSIFICATION OF EVIDENCE

This is a Class IV retrospective case series showing that structure-function mapping may help determine the outcome from surgical resection for treatment-resistant focal epilepsy.

Brain network communication: concepts, models and applications

Understanding communication and information processing in nervous systems is a central goal of neuroscience. Over the past two decades, advances in connectomics and network neuroscience have opened new avenues for investigating polysynaptic communication in complex brain networks. Recent work has brought into question the mainstay assumption that connectome signalling occurs exclusively via shortest paths, resulting in a sprawling constellation of alternative network communication models. This Review surveys the latest developments in models of brain network communication. We begin by drawing a conceptual link between the mathematics of graph theory and biological aspects of neural signalling such as transmission delays and metabolic cost. We organize key network communication models and measures into a taxonomy, aimed at helping researchers navigate the growing number of concepts and methods in the literature. The taxonomy highlights the pros, cons and interpretations of different conceptualizations of connectome signalling. We showcase the utility of network communication models as a flexible, interpretable and tractable framework to study brain function by reviewing prominent applications in basic, cognitive and clinical neurosciences. Finally, we provide recommendations to guide the future development, application and validation of network communication models.

Spatio-temporal evolution of human neural activity during visually cued hand movements

Making hand movements in response to visual cues is common in daily life. It has been well known that this process activates multiple areas in the brain, but how these neural activations progress across space and time remains largely unknown. Taking advantage of intracranial electroencephalographic (iEEG) recordings using depth and subdural electrodes from 36 human subjects using the same task, we applied single-trial and cross-trial analyses to high-frequency iEEG activity. The results show that the neural activation was widely distributed across the human brain both within and on the surface of the brain, and focused specifically on certain areas in the parietal, frontal, and occipital lobes, where parietal lobes present significant left lateralization on the activation. We also demonstrate temporal differences across these brain regions. Finally, we evaluated the degree to which the timing of activity within these regions was related to sensory or motor function. The findings of this study promote the understanding of task-related neural processing of the human brain, and may provide important insights for translational applications.

Co-evolving dynamics and topology in a coupled oscillator model of resting brain function

Dynamic models of ongoing BOLD fMRI brain dynamics and models of communication strategies have been two important approaches to understanding how brain network structure constrains function. However, dynamic models have yet to widely incorporate one of the most important insights from communication models: the brain may not use all of its connections in the same way or at the same time. Here we present a variation of a phase delayed Kuramoto coupled oscillator model that dynamically limits communication between nodes on each time step. An active subgraph of the empirically derived anatomical brain network is chosen in accordance with the local dynamic state on every time step, thus coupling dynamics and network structure in a novel way. We analyze this model with respect to its fit to empirical time-averaged functional connectivity, finding that, with the addition of only one parameter, it significantly outperforms standard Kuramoto models with phase delays. We also perform analyses on the novel time series of active edges it produces, demonstrating a slowly evolving topology moving through intermittent episodes of integration and segregation. We hope to demonstrate that the exploration of novel modeling mechanisms and the investigation of dynamics of networks in addition to dynamics on networks may advance our understanding of the relationship between brain structure and function.

Temporal stability of intracranial electroencephalographic abnormality maps for localizing epileptogenic tissue

OBJECTIVE

Identifying abnormalities on interictal intracranial electroencephalogram (iEEG), by comparing patient data to a normative map, has shown promise for the localization of epileptogenic tissue and prediction of outcome. The approach typically uses short interictal segments of approximately 1 min. However, the temporal stability of findings has not been established.

METHODS

Here, we generated a normative map of iEEG in nonpathological brain tissue from 249 patients. We computed regional band power abnormalities in a separate cohort of 39 patients for the duration of their monitoring period (.92-8.62 days of iEEG data, mean = 4.58 days per patient, >4800 hours recording). To assess the localizing value of band power abnormality, we computedD RS -a measure of how different the surgically resected and spared tissue was in terms of band power abnormalities-over time.

RESULTS

In each patient, theD RS value was relatively consistent over time. The medianD RS of the entire recording period separated seizure-free (International League Against Epilepsy [ILAE] = 1) and not-seizure-free (ILAE> 1) patients well (area under the curve [AUC] = .69). This effect was similar interictally (AUC = .69) and peri-ictally (AUC = .71).

SIGNIFICANCE

Our results suggest that band power abnormality D_RS, as a predictor of outcomes from epilepsy surgery, is a relatively robust metric over time. These findings add further support for abnormality mapping of neurophysiology data during presurgical evaluation.

The Interictal Suppression Hypothesis in focal epilepsy: network-level supporting evidence

Why are people with focal epilepsy not continuously having seizures? Previous neuronal signalling work has implicated gamma-aminobutyric acid balance as integral to seizure generation and termination, but is a high-level distributed brain network involved in suppressing seizures? Recent intracranial electrographic evidence has suggested that seizure-onset zones have increased inward connectivity that could be associated with interictal suppression of seizure activity. Accordingly, we hypothesize that seizure-onset zones are actively suppressed by the rest of the brain network during interictal states. Full testing of this hypothesis would require collaboration across multiple domains of neuroscience. We focused on partially testing this hypothesis at the electrographic network level within 81 individuals with drug-resistant focal epilepsy undergoing presurgical evaluation. We used intracranial electrographic resting-state and neurostimulation recordings to evaluate the network connectivity of seizure onset, early propagation and non-involved zones. We then used diffusion imaging to acquire estimates of white-matter connectivity to evaluate structure-function coupling effects on connectivity findings. Finally, we generated a resting-state classification model to assist clinicians in detecting seizure-onset and propagation zones without the need for multiple ictal recordings. Our findings indicate that seizure onset and early propagation zones demonstrate markedly increased inwards connectivity and decreased outwards connectivity using both resting-state (one-way ANOVA, P-value = 3.13 × 10-13) and neurostimulation analyses to evaluate evoked responses (one-way ANOVA, P-value = 2.5 × 10-3). When controlling for the distance between regions, the difference between inwards and outwards connectivity remained stable up to 80 mm between brain connections (two-way repeated measures ANOVA, group effect P-value of 2.6 × 10-12). Structure-function coupling analyses revealed that seizure-onset zones exhibit abnormally enhanced coupling (hypercoupling) of surrounding regions compared to presumably healthy tissue (two-way repeated measures ANOVA, interaction effect P-value of 9.76 × 10-21). Using these observations, our support vector classification models achieved a maximum held-out testing set accuracy of 92.0 ± 2.2% to classify early propagation and seizure-onset zones. These results suggest that seizure-onset zones are actively segregated and suppressed by a widespread brain network. Furthermore, this electrographically observed functional suppression is disproportionate to any observed structural connectivity alterations of the seizure-onset zones. These findings have implications for the identification of seizure-onset zones using only brief electrographic recordings to reduce patient morbidity and augment the presurgical evaluation of drug-resistant epilepsy. Further testing of the interictal suppression hypothesis can provide insight into potential new resective, ablative and neuromodulation approaches to improve surgical success rates in those suffering from drug-resistant focal epilepsy.

Intracranial electrophysiological and structural basis of BOLD functional connectivity in human brain white matter

While functional MRI (fMRI) studies have mainly focused on gray matter, recent studies have consistently found that blood-oxygenation-level-dependent (BOLD) signals can be reliably detected in white matter, and functional connectivity (FC) has been organized into distributed networks in white matter. Nevertheless, it remains unclear whether this white matter FC reflects underlying electrophysiological synchronization. To address this question, we employ intracranial stereotactic-electroencephalography (SEEG) and resting-state fMRI data from a group of 16 patients with drug-resistant epilepsy. We find that BOLD FC is correlated with SEEG FC in white matter, and this result is consistent across a wide range of frequency bands for each participant. By including diffusion spectrum imaging data, we also find that white matter FC from both SEEG and fMRI are correlated with white matter structural connectivity, suggesting that anatomical fiber tracts underlie the functional synchronization in white matter. These results provide evidence for the electrophysiological and structural basis of white matter BOLD FC, which could be a potential biomarker for psychiatric and neurological disorders.

Connectional Hierarchy in Human Brain Revealed by Individual Variability of Functional Network Edges

The human cerebral cortex is connected by intricate inter-areal wiring at the macroscale. The cortical hierarchy from primary sensorimotor to higher-order association areas is a unifying organizational principle across various neurobiological properties; however, previous studies have not clarified whether the connections between cortical regions exhibit a similar hierarchical pattern. Here, we identify a connectional hierarchy indexed by inter-individual variability of functional connectivity edges, which continuously progresses along a hierarchical gradient from within-network connections to between-network edges connecting sensorimotor and association networks. We found that this connectional hierarchy of variability aligns with both hemodynamic and electromagnetic connectivity strength and is constrained by structural connectivity strength. Moreover, the patterning of connectional hierarchy is related to inter-regional similarity in transcriptional and neurotransmitter receptor profiles. Using the Neurosynth cognitive atlas and cortical vulnerability maps in 13 brain disorders, we found that the connectional hierarchy of variability is associated with similarity networks of cognitive relevance and that of disorder vulnerability. Finally, we found that the prominence of this hierarchical gradient of connectivity variability declines during youth. Together, our results reveal a novel hierarchal organizational principle at the connectional level that links multimodal and multiscale human connectomes to individual variability in functional connectivity.

Communication dynamics in the human connectome shape the cortex-wide propagation of direct electrical stimulation

Communication between gray matter regions underpins all facets of brain function. We study inter-areal communication in the human brain using intracranial EEG recordings, acquired following 29,055 single-pulse direct electrical stimulations in a total of 550 individuals across 20 medical centers (average of 87 ± 37 electrode contacts per subject). We found that network communication models-computed on structural connectivity inferred from diffusion MRI-can explain the causal propagation of focal stimuli, measured at millisecond timescales. Building on this finding, we show that a parsimonious statistical model comprising structural, functional, and spatial factors can accurately and robustly predict cortex-wide effects of brain stimulation (R²=46% in data from held-out medical centers). Our work contributes toward the biological validation of concepts in network neuroscience and provides insight into how connectome topology shapes polysynaptic inter-areal signaling. We anticipate that our findings will have implications for research on neural communication and the design of brain stimulation paradigms.

High gamma activity distinguishes frontal cognitive control regions from adjacent cortical networks

Though the lateral frontal cortex is broadly implicated in cognitive control, functional MRI (fMRI) studies suggest fine-grained distinctions within this region. To examine this question electrophysiologically, we placed electrodes on the lateral frontal cortex in patients undergoing awake craniotomy for tumor resection. Patients performed verbal tasks with a manipulation of attentional switching, a canonical control demand. Power in the high gamma range (70-250 Hz) distinguished electrodes based on their location within a high-resolution fMRI network parcellation of the frontal lobe. Electrodes within the canonical fronto-parietal control network showed increased power in the switching condition, a result absent in electrodes within default mode, language and somato-motor networks. High gamma results contrasted with spatially distributed power decreases in the beta range (12-30 Hz). These results confirm the importance of fine-scale functional distinctions within the human frontal lobe, and pave the way for increased precision of functional mapping in tumor surgeries.

Longitudinal alterations of cortical structural-functional coupling in temporal lobe epilepsy

BACKGROUND AND PURPOSE

To investigate the longitudinal alterations of cortical structural-functional coupling (SF coupling) in patients with temporal lobe epilepsy (TLE) over a 2-year follow-up, thereby exploring the neuropathophysiological mechanisms of TLE.

METHODS

Twenty-eight TLE patients and 42 age- and gender-matched healthy controls (HCs) were recruited. We used resting-state functional MRI and diffusion-weighted imaging to estimate and compare SF coupling at the multiscale network level (whole-brain, modular, and regional levels). Then, we analyzed the relationships between the spatial patterns of SF coupling, the principal functional connectivity (FC) gradient, and the functional participation coefficient (PC). Finally, we related regional SF coupling changes between baseline and follow-up to the expression of regional TLE-specific genes.

RESULTS

Compared with HCs, TLE patients showed higher baseline SF couplings within the whole-brain, limbic, and default-mode modules. SF couplings within visual and dorsal attention modules were increased at follow-up compared to baseline. In all three groups, the spatial patterns of SF coupling aligned with the principal FC gradient and the functional PC. The longitudinal change in regional SF coupling in TLE patients was significantly positively correlated with the expression of the CUX2 gene.

CONCLUSIONS

Aberrant SF coupling was revealed in TLE and related to macroscale cortical hierarchies, functional segregation, and TLE-specific gene expression; these data help increase our understanding of the neuropathophysiological mechanisms underlying TLE.

Neural underpinnings of default mode network on empathy revealed by intracranial stereoelectroencephalography

BACKGROUND

Empathy is the ability to understand and share the feelings of others. It is fundamental to emotional intelligence and social iterations. Neuroimaging studies have demonstrated that empathy activates brain regions associated with the social cognition network.

AIM

To explore the neural underpinnings of empathy revealed by stereoelectroencephalography utilizing recurrence quantification analysis (RQA).

METHODS

This retrospective cohort included 38 epilepsy patients with stereoelectroencephalography implantation. RQA metrics were applied to parameterize the network organization of default mode network (DMN) brain regions. The relationships between DMN, seizure burden activity, and empathy, as measured using the Interpersonal Reactivity Index, were examined using partial least-square regression and mediation analysis.

RESULTS

RQA metrics with DMN (R²  = 0.75, P_Bonferroni  < 0.001) and its subsystems (medial temporal subsystem: R²  = 0.53, P_Bonferroni  < 0.001; core subsystem: R²  = 0.70, P_Bonferroni  < 0.001; dorsal medial subsystem: R²  = 0.48, P_Bonferroni  < 0.001) were positively correlated with empathy scores. Of 13 RQA metrics, the mean diagonal line length, entropy of the diagonal line lengths, trapping time, maximal vertical line length, and recurrence time of second type were found to be statistically higher in patient cohorts with reportedly high empathy. Furthermore, DMN characteristics (b path: F = 3.69, P = 0.04), rather than seizure burdens (direct effect: t = 0.33, P = 0.74, c' = - 0.007), mediated empathy status.

CONCLUSION

The present study used various RQA metrics to parameterize the network organization of DMN and determine the neural underpinning of DMN for empathy modulation.

Multimodal analysis demonstrating the shaping of functional gradients in the marmoset brain

The discovery of functional gradients introduce a new perspective in understanding the cortical spectrum of intrinsic dynamics, as it captures major axes of functional connectivity in low-dimensional space. However, how functional gradients arise and dynamically vary remains poorly understood. In this study, we investigated the biological basis of functional gradients using awake resting-state fMRI, retrograde tracing and gene expression datasets in marmosets. We found functional gradients in marmosets showed a sensorimotor-to-visual principal gradient followed by a unimodal-to-multimodal gradient, resembling functional gradients in human children. Although strongly constrained by structural wirings, functional gradients were dynamically modulated by arousal levels. Utilizing a reduced model, we uncovered opposing effects on gradient dynamics by structural connectivity (inverted U-shape) and neuromodulatory input (U-shape) with arousal fluctuations, and dissected the contribution of individual neuromodulatory receptors. This study provides insights into biological basis of functional gradients by revealing the interaction between structural connectivity and ascending neuromodulatory system.

Addressing spatial bias in intracranial EEG functional connectivity analyses for epilepsy surgical planning

Objective.To determine the effect of epilepsy on intracranial electroencephalography (EEG) functional connectivity, and the ability of functional connectivity to localize the seizure onset zone (SOZ), controlling for spatial biases.Approach.We analyzed intracranial EEG data from patients with drug-resistant epilepsy admitted for pre-surgical planning. We calculated intracranial EEG functional networks and determined whether changes in functional connectivity lateralized the SOZ using a spatial subsampling method to control for spatial bias. We developed a 'spatial null model' to localize the SOZ electrode using only spatial sampling information, ignoring EEG data. We compared the performance of this spatial null model against models incorporating EEG functional connectivity and interictal spike rates.Main results.About 110 patients were included in the study, although the number of patients differed across analyses. Controlling for spatial sampling, the average connectivity was lower in the SOZ region relative to the same anatomic region in the contralateral hemisphere. A model using intra-hemispheric connectivity accurately lateralized the SOZ (average accuracy 75.5%). A spatial null model incorporating spatial sampling information alone achieved moderate accuracy in classifying SOZ electrodes (mean AUC = 0.70, 95% CI 0.63-0.77). A model incorporating intracranial EEG functional connectivity and spike rate data further outperformed this spatial null model (AUC 0.78,p= 0.002 compared to spatial null model). However, a model incorporating functional connectivity without spike rate data did not significantly outperform the null model (AUC 0.72,p= 0.38).Significance.Intracranial EEG functional connectivity is reduced in the SOZ region, and interictal data predict SOZ electrode localization and laterality, however a predictive model incorporating functional connectivity without interictal spike rates did not significantly outperform a spatial null model. We propose constructing a spatial null model to provide an estimate of the pre-implant hypothesis of the SOZ, and to serve as a benchmark for further machine learning algorithms in order to avoid overestimating model performance because of electrode sampling alone.

Human electromagnetic and haemodynamic networks systematically converge in unimodal cortex and diverge in transmodal cortex

Whole-brain neural communication is typically estimated from statistical associations among electromagnetic or haemodynamic time-series. The relationship between functional network architectures recovered from these 2 types of neural activity remains unknown. Here, we map electromagnetic networks (measured using magnetoencephalography (MEG)) to haemodynamic networks (measured using functional magnetic resonance imaging (fMRI)). We find that the relationship between the 2 modalities is regionally heterogeneous and systematically follows the cortical hierarchy, with close correspondence in unimodal cortex and poor correspondence in transmodal cortex. Comparison with the BigBrain histological atlas reveals that electromagnetic–haemodynamic coupling is driven by laminar differentiation and neuron density, suggesting that the mapping between the 2 modalities can be explained by cytoarchitectural variation. Importantly, haemodynamic connectivity cannot be explained by electromagnetic activity in a single frequency band, but rather arises from the mixing of multiple neurophysiological rhythms. Correspondence between the two is largely driven by MEG functional connectivity at the beta (15 to 29 Hz) frequency band. Collectively, these findings demonstrate highly organized but only partly overlapping patterns of connectivity in MEG and fMRI functional networks, opening fundamentally new avenues for studying the relationship between cortical microarchitecture and multimodal connectivity patterns.

What is the relationship between functional network architectures inferred from electromagnetic and haemodynamic data? This study shows that superposition of electromagnetic networks at canonical frequency bands manifests as highly structured patterns of haemodynamic functional connectivity in the human brain, reflecting systematic variation in cytoarchitecture.

Toward an integrative neurovascular framework for studying brain networks

Brain functional connectivity based on the measure of blood oxygen level-dependent (BOLD) functional magnetic resonance imaging (fMRI) signals has become one of the most widely used measurements in human neuroimaging. However, the nature of the functional networks revealed by BOLD fMRI can be ambiguous, as highlighted by a recent series of experiments that have suggested that typical resting-state networks can be replicated from purely vascular or physiologically driven BOLD signals. After going through a brief review of the key concepts of brain network analysis, we explore how the vascular and neuronal systems interact to give rise to the brain functional networks measured with BOLD fMRI. This leads us to emphasize a view of the vascular network not only as a confounding element in fMRI but also as a functionally relevant system that is entangled with the neuronal network. To study the vascular and neuronal underpinnings of BOLD functional connectivity, we consider a combination of methodological avenues based on multiscale and multimodal optical imaging in mice, used in combination with computational models that allow the integration of vascular information to explain functional connectivity.

Advances in human intracranial electroencephalography research, guidelines and good practices

Since the second-half of the twentieth century, intracranial electroencephalography (iEEG), including both electrocorticography (ECoG) and stereo-electroencephalography (sEEG), has provided an intimate view into the human brain. At the interface between fundamental research and the clinic, iEEG provides both high temporal resolution and high spatial specificity but comes with constraints, such as the individual's tailored sparsity of electrode sampling. Over the years, researchers in neuroscience developed their practices to make the most of the iEEG approach. Here we offer a critical review of iEEG research practices in a didactic framework for newcomers, as well addressing issues encountered by proficient researchers. The scope is threefold: (i) review common practices in iEEG research, (ii) suggest potential guidelines for working with iEEG data and answer frequently asked questions based on the most widespread practices, and (iii) based on current neurophysiological knowledge and methodologies, pave the way to good practice standards in iEEG research. The organization of this paper follows the steps of iEEG data processing. The first section contextualizes iEEG data collection. The second section focuses on localization of intracranial electrodes. The third section highlights the main pre-processing steps. The fourth section presents iEEG signal analysis methods. The fifth section discusses statistical approaches. The sixth section draws some unique perspectives on iEEG research. Finally, to ensure a consistent nomenclature throughout the manuscript and to align with other guidelines, e.g., Brain Imaging Data Structure (BIDS) and the OHBM Committee on Best Practices in Data Analysis and Sharing (COBIDAS), we provide a glossary to disambiguate terms related to iEEG research.

Normative intracranial EEG maps epileptogenic tissues in focal epilepsy

Planning surgery for patients with medically refractory epilepsy often requires recording seizures using intracranial EEG. Quantitative measures derived from interictal intracranial EEG yield potentially appealing biomarkers to guide these surgical procedures; however, their utility is limited by the sparsity of electrode implantation as well as the normal confounds of spatiotemporally varying neural activity and connectivity. We propose that comparing intracranial EEG recordings to a normative atlas of intracranial EEG activity and connectivity can reliably map abnormal regions, identify targets for invasive treatment and increase our understanding of human epilepsy. Merging data from the Penn Epilepsy Center and a public database from the Montreal Neurological Institute, we aggregated interictal intracranial EEG retrospectively across 166 subjects comprising >5000 channels. For each channel, we calculated the normalized spectral power and coherence in each canonical frequency band. We constructed an intracranial EEG atlas by mapping the distribution of each feature across the brain and tested the atlas against data from novel patients by generating a z-score for each channel. We demonstrate that for seizure onset zones within the mesial temporal lobe, measures of connectivity abnormality provide greater distinguishing value than univariate measures of abnormal neural activity. We also find that patients with a longer diagnosis of epilepsy have greater abnormalities in connectivity. By integrating measures of both single-channel activity and inter-regional functional connectivity, we find a better accuracy in predicting the seizure onset zones versus normal brain (area under the curve = 0.77) compared with either group of features alone. We propose that aggregating normative intracranial EEG data across epilepsy centres into a normative atlas provides a rigorous, quantitative method to map epileptic networks and guide invasive therapy. We publicly share our data, infrastructure and methods, and propose an international framework for leveraging big data in surgical planning for refractory epilepsy.

Structural insight into the individual variability architecture of the functional brain connectome

Human cognition and behaviors depend upon the brain's functional connectomes, which vary remarkably across individuals. However, whether and how the functional connectome individual variability architecture is structurally constrained remains largely unknown. Using tractography- and morphometry-based network models, we observed the spatial convergence of structural and functional connectome individual variability, with higher variability in heteromodal association regions and lower variability in primary regions. We demonstrated that functional variability is significantly predicted by a unifying structural variability pattern and that this prediction follows a primary-to-heteromodal hierarchical axis, with higher accuracy in primary regions and lower accuracy in heteromodal regions. We further decomposed group-level connectome variability patterns into individual unique contributions and uncovered the structural-functional correspondence that is associated with individual cognitive traits. These results advance our understanding of the structural basis of individual functional variability and suggest the importance of integrating multimodal connectome signatures for individual differences in cognition and behaviors.

Network communication models narrow the gap between the modular organization of structural and functional brain networks

Structural and functional brain networks are modular. Canonical functional systems, such as the default mode network, are well-known modules of the human brain and have been implicated in a large number of cognitive, behavioral and clinical processes. However, modules delineated in structural brain networks inferred from tractography generally do not recapitulate canonical functional systems. Neuroimaging evidence suggests that functional connectivity between regions in the same systems is not always underpinned by anatomical connections. As such, direct structural connectivity alone would be insufficient to characterize the functional modular organization of the brain. Here, we demonstrate that augmenting structural brain networks with models of indirect (polysynaptic) communication unveils a modular network architecture that more closely resembles the brain's established functional systems. We find that diffusion models of polysynaptic connectivity, particularly communicability, narrow the gap between the modular organization of structural and functional brain networks by 20-60%, whereas routing models based on single efficient paths do not improve mesoscopic structure-function correspondence. This suggests that functional modules emerge from the constraints imposed by local network structure that facilitates diffusive neural communication. Our work establishes the importance of modeling polysynaptic communication to understand the structural basis of functional systems.

Social cognitive network neuroscience

Over the past three decades, research from the field of social neuroscience has identified a constellation of brain regions that relate to social cognition. Although these studies have provided important insights into the specific neural regions underlying social behavior, they may overlook the broader neural context in which those regions and the interactions between them are embedded. Network neuroscience is an emerging discipline that focuses on modeling and analyzing brain networks-collections of interacting neural elements. Because human cognition requires integrating information across multiple brain regions and systems, we argue that a novel social cognitive network neuroscience approach-which leverages methods from the field of network neuroscience and graph theory-can advance our understanding of how brain systems give rise to social behavior. This review provides an overview of the field of network neuroscience, discusses studies that have leveraged this approach to advance social neuroscience research, highlights the potential contributions of social cognitive network neuroscience to understanding social behavior and provides suggested tools and resources for conducting network neuroscience research.

Normative brain mapping of interictal intracranial EEG to localize epileptogenic tissue

The identification of abnormal electrographic activity is important in a wide range of neurological disorders, including epilepsy for localizing epileptogenic tissue. However, this identification may be challenging during non-seizure (interictal) periods, especially if abnormalities are subtle compared to the repertoire of possible healthy brain dynamics. Here, we investigate if such interictal abnormalities become more salient by quantitatively accounting for the range of healthy brain dynamics in a location-specific manner. To this end, we constructed a normative map of brain dynamics, in terms of relative band power, from interictal intracranial recordings from 234 participants (21 598 electrode contacts). We then compared interictal recordings from 62 patients with epilepsy to the normative map to identify abnormal regions. We proposed that if the most abnormal regions were spared by surgery, then patients would be more likely to experience continued seizures postoperatively. We first confirmed that the spatial variations of band power in the normative map across brain regions were consistent with healthy variations reported in the literature. Second, when accounting for the normative variations, regions that were spared by surgery were more abnormal than those resected only in patients with persistent postoperative seizures (t = -3.6, P = 0.0003), confirming our hypothesis. Third, we found that this effect discriminated patient outcomes (area under curve 0.75 P = 0.0003). Normative mapping is a well-established practice in neuroscientific research. Our study suggests that this approach is feasible to detect interictal abnormalities in intracranial EEG, and of potential clinical value to identify pathological tissue in epilepsy. Finally, we make our normative intracranial map publicly available to facilitate future investigations in epilepsy and beyond.

Local structure-function relationships in human brain networks across the lifespan

A growing number of studies have used stylized network models of communication to predict brain function from structure. Most have focused on a small set of models applied globally. Here, we compare a large number of models at both global and regional levels. We find that globally most predictors perform poorly. At the regional level, performance improves but heterogeneously, both in terms of variance explained and the optimal model. Next, we expose synergies among predictors by using pairs to jointly predict FC. Finally, we assess age-related differences in global and regional coupling across the human lifespan. We find global decreases in the magnitude of structure-function coupling with age. We find that these decreases are driven by reduced coupling in sensorimotor regions, while higher-order cognitive systems preserve local coupling with age. Our results describe patterns of structure-function coupling across the cortex and how this may change with age.

Electrophysiological foundations of the human default-mode network revealed by intracranial-EEG recordings during resting-state and cognition

Investigations using noninvasive functional magnetic resonance imaging (fMRI) have provided significant insights into the unique functional organization and profound importance of the human default mode network (DMN), yet these methods are limited in their ability to resolve network dynamics across multiple timescales. Electrophysiological techniques are critical to address these challenges, yet few studies have explored the neurophysiological underpinnings of the DMN. Here we investigate the electrophysiological organization of the DMN in a common large-scale network framework consistent with prior fMRI studies. We used intracranial EEG (iEEG) recordings, and evaluated intra- and cross-network interactions during resting-state and its modulation during a cognitive task involving episodic memory formation. Our analysis revealed significantly greater intra-DMN phase iEEG synchronization in the slow-wave (< 4 Hz), while DMN interactions with other brain networks was higher in the beta (12-30 Hz) and gamma (30-80 Hz) bands. Crucially, slow-wave intra-DMN synchronization was observed in the task-free resting-state and during both verbal memory encoding and recall. Compared to resting-state, slow-wave intra-DMN phase synchronization was significantly higher during both memory encoding and recall. Slow-wave intra-DMN phase synchronization increased during successful memory retrieval, highlighting its behavioral relevance. Finally, analysis of nonlinear dynamic causal interactions revealed that the DMN is a causal outflow network during both memory encoding and recall. Our findings identify frequency specific neurophysiological signatures of the DMN which allow it to maintain stability and flexibility, intrinsically and during task-based cognition, provide novel insights into the electrophysiological foundations of the human DMN, and elucidate network mechanisms by which it supports cognition.

A framework For brain atlases: Lessons from seizure dynamics

Brain maps, or atlases, are essential tools for studying brain function and organization. The abundance of available atlases used across the neuroscience literature, however, creates an implicit challenge that may alter the hypotheses and predictions we make about neurological function and pathophysiology. Here, we demonstrate how parcellation scale, shape, anatomical coverage, and other atlas features may impact our prediction of the brain’s function from its underlying structure. We show how network topology, structure–function correlation (SFC), and the power to test specific hypotheses about epilepsy pathophysiology may change as a result of atlas choice and atlas features. Through the lens of our disease system, we propose a general framework and algorithm for atlas selection. This framework aims to maximize the descriptive, explanatory, and predictive validity of an atlas. Broadly, our framework strives to provide empirical guidance to neuroscience research utilizing the various atlases published over the last century.

Neurophysiological Evidence for Cognitive Map Formation during Sequence Learning

Humans deftly parse statistics from sequences. Some theories posit that humans learn these statistics by forming cognitive maps, or underlying representations of the latent space which links items in the sequence. Here, an item in the sequence is a node, and the probability of transitioning between two items is an edge. Sequences can then be generated from walks through the latent space, with different spaces giving rise to different sequence statistics. Individual or group differences in sequence learning can be modeled by changing the time scale over which estimates of transition probabilities are built, or in other words, by changing the amount of temporal discounting. Latent space models with temporal discounting bear a resemblance to models of navigation through Euclidean spaces. However, few explicit links have been made between predictions from Euclidean spatial navigation and neural activity during human sequence learning. Here, we use a combination of behavioral modeling and intracranial encephalography (iEEG) recordings to investigate how neural activity might support the formation of space-like cognitive maps through temporal discounting during sequence learning. Specifically, we acquire human reaction times from a sequential reaction time task, to which we fit a model that formulates the amount of temporal discounting as a single free parameter. From the parameter, we calculate each individual's estimate of the latent space. We find that neural activity reflects these estimates mostly in the temporal lobe, including areas involved in spatial navigation. Similar to spatial navigation, we find that low-dimensional representations of neural activity allow for easy separation of important features, such as modules, in the latent space. Lastly, we take advantage of the high temporal resolution of iEEG data to determine the time scale on which latent spaces are learned. We find that learning typically happens within the first 500 trials, and is modulated by the underlying latent space and the amount of temporal discounting characteristic of each participant. Ultimately, this work provides important links between behavioral models of sequence learning and neural activity during the same behavior, and contextualizes these results within a broader framework of domain general cognitive maps.

Modularity maximization as a flexible and generic framework for brain network exploratory analysis

The modular structure of brain networks supports specialized information processing, complex dynamics, and cost-efficient spatial embedding. Inter-individual variation in modular structure has been linked to differences in performance, disease, and development. There exist many data-driven methods for detecting and comparing modular structure, the most popular of which is modularity maximization. Although modularity maximization is a general framework that can be modified and reparamaterized to address domain-specific research questions, its application to neuroscientific datasets has, thus far, been narrow. Here, we highlight several strategies in which the "out-of-the-box" version of modularity maximization can be extended to address questions specific to neuroscience. First, we present approaches for detecting "space-independent" modules and for applying modularity maximization to signed matrices. Next, we show that the modularity maximization frame is well-suited for detecting task- and condition-specific modules. Finally, we highlight the role of multi-layer models in detecting and tracking modules across time, tasks, subjects, and modalities. In summary, modularity maximization is a flexible and general framework that can be adapted to detect modular structure resulting from a wide range of hypotheses. This article highlights multiple frontiers for future research and applications.

Modular origins of high-amplitude cofluctuations in fine-scale functional connectivity dynamics

The topology of structural brain networks shapes brain dynamics, including the correlation structure of brain activity (functional connectivity) as estimated from functional neuroimaging data. Empirical studies have shown that functional connectivity fluctuates over time, exhibiting patterns that vary in the spatial arrangement of correlations among segregated functional systems. Recently, an exact decomposition of functional connectivity into frame-wise contributions has revealed fine-scale dynamics that are punctuated by brief and intermittent episodes (events) of high-amplitude cofluctuations involving large sets of brain regions. Their origin is currently unclear. Here, we demonstrate that similar episodes readily appear in silico using computational simulations of whole-brain dynamics. As in empirical data, simulated events contribute disproportionately to long-time functional connectivity, involve recurrence of patterned cofluctuations, and can be clustered into distinct families. Importantly, comparison of event-related patterns of cofluctuations to underlying patterns of structural connectivity reveals that modular organization present in the coupling matrix shapes patterns of event-related cofluctuations. Our work suggests that brief, intermittent events in functional dynamics are partly shaped by modular organization of structural connectivity.

The diversity and multiplexity of edge communities within and between brain systems

The human brain is composed of functionally specialized systems that support cognition. Recently, we proposed an edge-centric model for detecting overlapping communities. It remains unclear how these communities and brain systems are related. Here, we address this question using data from the Midnight Scan Club and show that all brain systems are linked via at least two edge communities. We then examine the diversity of edge communities within each system, finding that heteromodal systems are more diverse than sensory systems. Next, we cluster the entire cortex to reveal it according to the regions' edge-community profiles. We find that regions in heteromodal systems are more likely to form their own clusters. Finally, we show that edge communities are personalized. Our work reveals the pervasive overlap of edge communities across the cortex and their relationship with brain systems. Our work provides pathways for future research using edge-centric brain networks.

Distributed Subnetworks of Depression Defined by Direct Intracranial Neurophysiology

Major depressive disorder is a common and disabling disorder with high rates of treatment resistance. Evidence suggests it is characterized by distributed network dysfunction that may be variable across patients, challenging the identification of quantitative biological substrates. We carried out this study to determine whether application of a novel computational approach to a large sample of high spatiotemporal resolution direct neural recordings in humans could unlock the functional organization and coordinated activity patterns of depression networks. This group level analysis of depression networks from heterogenous intracranial recordings was possible due to application of a correlational model-based method for inferring whole-brain neural activity. We then applied a network framework to discover brain dynamics across this model that could classify depression. We found a highly distributed pattern of neural activity and connectivity across cortical and subcortical structures that was present in the majority of depressed subjects. Furthermore, we found that this depression signature consisted of two subnetworks across individuals. The first was characterized by left temporal lobe hypoconnectivity and pathological beta activity. The second was characterized by a hypoactive, but hyperconnected left frontal cortex. These findings have applications toward personalization of therapy.

Mapping gene transcription and neurocognition across human neocortex

Regulation of gene expression drives protein interactions that govern synaptic wiring and neuronal activity. The resulting coordinated activity among neuronal populations supports complex psychological processes, yet how gene expression shapes cognition and emotion remains unknown. Here, we directly bridge the microscale and macroscale by mapping gene expression patterns to functional activation patterns across the cortical sheet. Applying unsupervised learning to the Allen Human Brain Atlas and Neurosynth databases, we identify a ventromedial-dorsolateral gradient of gene assemblies that separate affective and perceptual domains. This topographic molecular-psychological signature reflects the hierarchical organization of the neocortex, including systematic variations in cell type, myeloarchitecture, laminar differentiation and intrinsic network affiliation. In addition, this molecular-psychological signature strengthens over neurodevelopment and can be replicated in two independent repositories. Collectively, our results reveal spatially covarying transcriptomic and cognitive architectures, highlighting the influence that molecular mechanisms exert on psychological processes.

Mesoscopic physiological interactions in the human brain reveal small-world properties

Cognition depends on rapid and robust communication between neural circuits spanning different brain areas. We investigated the mesoscopic network of cortico-cortical interactions in the human brain in an extensive dataset consisting of 6,024 h of intracranial field potential recordings from 4,142 electrodes in 48 subjects. We evaluated communication between brain areas at the network level across different frequency bands. The interaction networks were validated against known anatomical measurements and neurophysiological interactions in humans and monkeys. The resulting human brain interactome is characterized by a broad and spatially specific, dynamic, and extensive network. The physiological interactome reveals small-world properties, which we conjecture might facilitate efficient and reliable information transmission. The interaction dynamics correlate with the brain sleep/awake state. These results constitute initial steps toward understanding how the interactome orchestrates cortical communication and provide a reference for future efforts assessing how dysfunctional interactions may lead to mental disorders.

Cognition relies on rapid and robust communication between brain areas. Wang et al. leverage multi-day intracranial field potential recordings to characterize the human mesoscopic functional interactome. They validated the methods using monkey anatomical and physiological data. The human interactome reveals small-world properties and is modulated by sleep versus awake state.

Evaluating the sensitivity of functional connectivity measures to motion artifact in resting-state fMRI data

Functional connectivity (FC) networks are typically inferred from resting-state fMRI data using the Pearson correlation between BOLD time series from pairs of brain regions. However, alternative methods of estimating functional connectivity have not been systematically tested for their sensitivity or robustness to head motion artifact. Here, we evaluate the sensitivity of eight different functional connectivity measures to motion artifact using resting-state data from the Human Connectome Project. We report that FC estimated using full correlation has a relatively high residual distance-dependent relationship with motion compared to partial correlation, coherence, and information theory-based measures, even after implementing rigorous methods for motion artifact mitigation. This disadvantage of full correlation, however, may be offset by higher test-retest reliability, fingerprinting accuracy, and system identifiability. FC estimated by partial correlation offers the best of both worlds, with low sensitivity to motion artifact and intermediate system identifiability, with the caveat of low test-retest reliability and fingerprinting accuracy. We highlight spatial differences in the sub-networks affected by motion with different FC metrics. Further, we report that intra-network edges in the default mode and retrosplenial temporal sub-networks are highly correlated with motion in all FC methods. Our findings indicate that the method of estimating functional connectivity is an important consideration in resting-state fMRI studies and must be chosen carefully based on the parameters of the study.

Electrocorticography and stereo EEG provide distinct measures of brain connectivity: implications for network models

Brain network models derived from graph theory have the potential to guide functional neurosurgery, and to improve rates of post-operative seizure freedom for patients with epilepsy. A barrier to applying these models clinically is that intracranial EEG electrode implantation strategies vary by centre, region and country, from cortical grid & strip electrodes (Electrocorticography), to purely stereotactic depth electrodes (Stereo EEG), to a mixture of both. To determine whether models derived from one type of study are broadly applicable to others, we investigate the differences in brain networks mapped by electrocorticography and stereo EEG in a cohort of patients who underwent surgery for temporal lobe epilepsy and achieved a favourable outcome. We show that networks derived from electrocorticography and stereo EEG define distinct relationships between resected and spared tissue, which may be driven by sampling bias of temporal depth electrodes in patients with predominantly cortical grids. We propose a method of correcting for the effect of internodal distance that is specific to electrode type and explore how additional methods for spatially correcting for sampling bias affect network models. Ultimately, we find that smaller surgical targets tend to have lower connectivity with respect to the surrounding network, challenging notions that abnormal connectivity in the epileptogenic zone is typically high. Our findings suggest that effectively applying computational models to localize epileptic networks requires accounting for the effects of spatial sampling, particularly when analysing both electrocorticography and stereo EEG recordings in the same cohort, and that future network studies of epilepsy surgery should also account for differences in focality between resection and ablation. We propose that these findings are broadly relevant to intracranial EEG network modelling in epilepsy and an important step in translating them clinically into patient care.

Cortico-muscular functional network: an exploration of cortico-muscular coupling in hand movements

Objective. The main objective of this research was to study cortico-muscular, intra-cortical, and inter-muscular coupling. Herein, we established a cortico-muscular functional network (CMFN) to assess the network differences associated with making a fist, opening the hand, and wrist flexion.Approach. We used transfer entropy (TE) to calculate the causality between electroencephalographic and electromyographic data and established the TE connection matrix. We then applied graph theory to analyze the clustering coefficient, global efficiency, and small-world attributes of the CMFN. We also used Relief-F to extract the features of the TE connection matrix of the beta2 band for the different hand movements and observed high accuracy when this feature was used for action recognition.Main results. We found that the CMFN of the three actions in the beta band had small-world attributes, among which the beta2 band's small-world was stronger. Moreover, we found that the extracted features were mainly concentrated in the left frontal area, left motor area, occipital lobe, and related muscles, suggesting that the CMFN could be used to assess the coupling differences between the cortex and the muscles that are associated with different hand movements. Overall, our results showed that the beta2 (21-35 Hz) wave is the main information carrier between the cortex and the muscles, and the CMFN can be used in the beta2 band to assess cortico-muscular coupling.Significance. Our study preliminarily explored the CMFN associated with hand movements, providing additional insights regarding the transmission of information between the cortex and the muscles, thereby laying a foundation for future rehabilitation therapy targeting pathological cortical areas in stroke patients.

Overcoming false-positive gene-category enrichment in the analysis of spatially resolved transcriptomic brain atlas data

Transcriptomic atlases have improved our understanding of the correlations between gene-expression patterns and spatially varying properties of brain structure and function. Gene-category enrichment analysis (GCEA) is a common method to identify functional gene categories that drive these associations, using gene-to-category annotation systems like the Gene Ontology (GO). Here, we show that applying standard GCEA methodology to spatial transcriptomic data is affected by substantial false-positive bias, with GO categories displaying an over 500-fold average inflation of false-positive associations with random neural phenotypes in mouse and human. The estimated false-positive rate of a GO category is associated with its rate of being reported as significantly enriched in the literature, suggesting that published reports are affected by this false-positive bias. We show that within-category gene-gene coexpression and spatial autocorrelation are key drivers of the false-positive bias and introduce flexible ensemble-based null models that can account for these effects, made available as a software toolbox.

Pairwise maximum entropy model explains the role of white matter structure in shaping emergent co-activation states

A major challenge in neuroscience is determining a quantitative relationship between the brain's white matter structural connectivity and emergent activity. We seek to uncover the intrinsic relationship among brain regions fundamental to their functional activity by constructing a pairwise maximum entropy model (MEM) of the inter-ictal activation patterns of five patients with medically refractory epilepsy over an average of ~14 hours of band-passed intracranial EEG (iEEG) recordings per patient. We find that the pairwise MEM accurately predicts iEEG electrodes' activation patterns' probability and their pairwise correlations. We demonstrate that the estimated pairwise MEM's interaction weights predict structural connectivity and its strength over several frequencies significantly beyond what is expected based solely on sampled regions' distance in most patients. Together, the pairwise MEM offers a framework for explaining iEEG functional connectivity and provides insight into how the brain's structural connectome gives rise to large-scale activation patterns by promoting co-activation between connected structures.