Mapping human brain networks with cortico-cortical evoked potentials

CARLA: Adjusted common average referencing for cortico-cortical evoked potential data

Human brain connectivity can be mapped by single pulse electrical stimulation during intracranial EEG measurements. The raw cortico-cortical evoked potentials (CCEP) are often contaminated by noise. Common average referencing (CAR) removes common noise and preserves response shapes but can introduce bias from responsive channels. We address this issue with an adjusted, adaptive CAR algorithm termed "CAR by Least Anticorrelation (CARLA)". CARLA was tested on simulated CCEP data and real CCEP data collected from four human participants. In CARLA, the channels are ordered by increasing mean cross-trial covariance, and iteratively added to the common average until anticorrelation between any single channel and all re-referenced channels reaches a minimum, as a measure of shared noise. We simulated CCEP data with true responses in 0-45 of 50 total channels. We quantified CARLA's error and found that it erroneously included 0 (median) truly responsive channels in the common average with ≤42 responsive channels, and erroneously excluded ≤2.5 (median) unresponsive channels at all responsiveness levels. On real CCEP data, signal quality was quantified with the mean R2 between all pairs of channels, which represents inter-channel dependency and is low for well-referenced data. CARLA re-referencing produced significantly lower mean R2 than standard CAR, CAR using a fixed bottom quartile of channels by covariance, and no re-referencing. CARLA minimizes bias in re-referenced CCEP data by adaptively selecting the optimal subset of non-responsive channels. It showed high specificity and sensitivity on simulated CCEP data and lowered inter-channel dependency compared to CAR on real CCEP data.

Causal Cortical and Thalamic Connections in the Human Brain

The brain's functional architecture is intricately shaped by causal connections between its cortical and subcortical structures. Here, we studied 27 participants with 4864 electrodes implanted across the anterior, mediodorsal, and pulvinar thalamic regions, and the cortex. Using data from electrical stimulation procedures and a data-driven approach informed by neurophysiological standards, we dissociated three unique spectral patterns generated by the perturbation of a given brain area. Among these, a novel waveform emerged, marked by delayed-onset slow oscillations in both ipsilateral and contralateral cortices following thalamic stimulations, suggesting a mechanism by which a thalamic site can influence bilateral cortical activity. Moreover, cortical stimulations evoked earlier signals in the thalamus than in other connected cortical areas suggesting that the thalamus receives a copy of signals before they are exchanged across the cortex. Our causal connectivity data can be used to inform biologically-inspired computational models of the functional architecture of the brain.

Intracortical mechanisms of single pulse electrical stimulation (SPES) evoked excitations and inhibitions in humans

Cortico-cortical evoked potentials (CCEPs) elicited by single-pulse electric stimulation (SPES) are widely used to assess effective connectivity between cortical areas and are also implemented in the presurgical evaluation of epileptic patients. Nevertheless, the cortical generators underlying the various components of CCEPs in humans have not yet been elucidated. Our aim was to describe the laminar pattern arising under SPES evoked CCEP components (P1, N1, P2, N2, P3) and to evaluate the similarities between N2 and the downstate of sleep slow waves. We used intra-cortical laminar microelectrodes (LMEs) to record CCEPs evoked by 10 mA bipolar 0.5 Hz electric pulses in seven patients with medically intractable epilepsy implanted with subdural grids. Based on the laminar profile of CCEPs, the latency of components is not layer-dependent, however their rate of appearance varies across cortical depth and stimulation distance, while the seizure onset zone does not seem to affect the emergence of components. Early neural excitation primarily engages middle and deep layers, propagating to the superficial layers, followed by mainly superficial inhibition, concluding in a sleep slow wave-like inhibition and excitation sequence.

Causal Cortical and Thalamic Connections in the Human Brain

The brain's functional architecture is intricately shaped by causal connections between its cortical and subcortical structures. Here, we studied 27 participants with 4864 electrodes implanted across the anterior, mediodorsal, and pulvinar thalamic regions, and the cortex. Using data from electrical stimulation procedures and a data-driven approach informed by neurophysiological standards, we dissociated three unique spectral patterns generated by the perturbation of a given brain area. Among these, a novel waveform emerged, marked by delayed-onset slow oscillations in both ipsilateral and contralateral cortices following thalamic stimulations, suggesting a mechanism by which a thalamic site can influence bilateral cortical activity. Moreover, cortical stimulations evoked earlier signals in the thalamus than in other connected cortical areas suggesting that the thalamus receives a copy of signals before they are exchanged across the cortex. Our causal connectivity data can be used to inform biologically-inspired computational models of the functional architecture of the brain.

Effects of transcranial magnetic stimulation on the human brain recorded with intracranial electrocorticography

Transcranial magnetic stimulation (TMS) is increasingly used as a noninvasive technique for neuromodulation in research and clinical applications, yet its mechanisms are not well understood. Here, we present the neurophysiological effects of TMS using intracranial electrocorticography (iEEG) in neurosurgical patients. We first evaluated safety in a gel-based phantom. We then performed TMS-iEEG in 22 neurosurgical participants with no adverse events. We next evaluated intracranial responses to single pulses of TMS to the dorsolateral prefrontal cortex (dlPFC) (N = 10, 1414 electrodes). We demonstrate that TMS is capable of inducing evoked potentials both locally within the dlPFC and in downstream regions functionally connected to the dlPFC, including the anterior cingulate and insular cortex. These downstream effects were not observed when stimulating other distant brain regions. Intracranial dlPFC electrical stimulation had similar timing and downstream effects as TMS. These findings support the safety and promise of TMS-iEEG in humans to examine local and network-level effects of TMS with higher spatiotemporal resolution than currently available methods.

Cortico-cortical evoked potentials of language tracts in minimally invasive glioma surgery guided by Penfield stimulation

OBJECTIVE

We investigated the feasibility of recording cortico-cortical evoked potentials (CCEPs) in patients with low- and high-grade glioma. We compared CCEPs during awake and asleep surgery, as well as those stimulated from the functional Broca area and recorded from the functional Wernicke area (BtW), and vice versa (WtB). We also analyzed CCEP properties according to tumor location, histopathology, and aphasia.

METHODS

We included 20 patients who underwent minimally invasive surgery in an asleep-awake-asleep setting. Strip electrode placement was guided by classical Penfield stimulation of positive language sites and fiber tracking of the arcuate fascicle. CCEPs were elicited with alternating monophasic single pulses of 1.1 Hz frequency and recorded as averaged signals. Intraoperatively, there was no post-processing of the signal.

RESULTS

Ninety-seven CCEPs from 19 patients were analyzed. There was no significant difference in CCEP properties when comparing awake versus asleep, nor BtW versus WtB. CCEP amplitude and latency were affected by tumor location and histopathology. CCEP features after tumor resection correlated with short- and long-term postoperative aphasia.

CONCLUSION

CCEP recordings are feasible during minimally invasive surgery. CCEPs might be surrogate markers for altered connectivity of the language tracts.

SIGNIFICANCE

This study may guide the incorporation of CCEPs into intraoperative neurophysiological monitoring.

Thalamic feedback shapes brain responses evoked by cortical stimulation in mice and humans

Cortical stimulation with single pulses is a common technique in clinical practice and research. However, we still do not understand the extent to which it engages subcortical circuits which contribute to the associated evoked potentials (EPs). Here we find that cortical stimulation generates remarkably similar EPs in humans and mice, with a late component similarly modulated by the subject's behavioral state. We optogenetically dissect the underlying circuit in mice, demonstrating that the late component of these EPs is caused by a thalamic hyperpolarization and rebound. The magnitude of this late component correlates with the bursting frequency and synchronicity of thalamic neurons, modulated by the subject's behavioral state. A simulation of the thalamo-cortical circuit highlights that both intrinsic thalamic currents as well as cortical and thalamic GABAergic neurons contribute to this response profile. We conclude that the cortical stimulation engages cortico-thalamo-cortical circuits highly preserved across different species and stimulation modalities.

Galvanic vs. pulsatile effects on decision-making networks: reshaping the neural activation landscape

Objective. Primarily due to safety concerns, biphasic pulsatile stimulation (PS) is the present standard for electrical excitation of neural tissue with a diverse set of applications. While pulses have been shown to be effective to achieve functional outcomes, they have well-known deficits. Due to recent technical advances, galvanic stimulation (GS), delivery of current for extended periods of time (>1 s), has re-emerged as an alternative to PS.Approach. In this paper, we use a winner-take-all decision-making cortical network model to investigate differences between pulsatile and GS in the context of a perceptual decision-making task.Main results. Based on previous work, we hypothesized that GS would produce more spatiotemporally distributed, network-sensitive neural responses, while PS would produce highly synchronized activation of a limited group of neurons. Our results in-silico support these hypotheses for low-amplitude GS but deviate when galvanic amplitudes are large enough to directly activate or block nearby neurons.Significance. We conclude that with careful parametrization, GS could overcome some limitations of PS to deliver more naturalistic firing patterns in the group of targeted neurons.

Subcallosal cingulate deep brain stimulation evokes two distinct cortical responses via differential white matter activation

Subcallosal cingulate (SCC) deep brain stimulation (DBS) is an emerging therapy for refractory depression. Good clinical outcomes are associated with the activation of white matter adjacent to the SCC. This activation produces a signature cortical evoked potential (EP), but it is unclear which of the many pathways in the vicinity of SCC is responsible for driving this response. Individualized biophysical models were built to achieve selective engagement of two target bundles: either the forceps minor (FM) or cingulum bundle (CB). Unilateral 2 Hz stimulation was performed in seven patients with treatment-resistant depression who responded to SCC DBS, and EPs were recorded using 256-sensor scalp electroencephalography. Two distinct EPs were observed: a 120 ms symmetric response spanning both hemispheres and a 60 ms asymmetrical EP. Activation of FM correlated with the symmetrical EPs, while activation of CB was correlated with the asymmetrical EPs. These results support prior model predictions that these two pathways are predominantly activated by clinical SCC DBS and provide first evidence of a link between cortical EPs and selective fiber bundle activation.

Temporally organized representations of reward and risk in the human brain

The value and uncertainty associated with choice alternatives constitute critical features relevant for decisions. However, the manner in which reward and risk representations are temporally organized in the brain remains elusive. Here we leverage the spatiotemporal precision of intracranial electroencephalography, along with a simple card game designed to elicit the unfolding computation of a set of reward and risk variables, to uncover this temporal organization. Reward outcome representations across wide-spread regions follow a sequential order along the anteroposterior axis of the brain. In contrast, expected value can be decoded from multiple regions at the same time, and error signals in both reward and risk domains reflect a mixture of sequential and parallel encoding. We further highlight the role of the anterior insula in generalizing between reward prediction error and risk prediction error codes. Together our results emphasize the importance of neural dynamics for understanding value-based decisions under uncertainty.

Thalamic stimulation induced changes in effective connectivity

Deep brain stimulation (DBS) is a viable treatment for a variety of neurological conditions, however, the mechanisms through which DBS modulates large-scale brain networks are unresolved. Clinical effects of DBS are observed over multiple timescales. In some conditions, such as Parkinson's disease and essential tremor, clinical improvement is observed within seconds. In many other conditions, such as epilepsy, central pain, dystonia, neuropsychiatric conditions or Tourette syndrome, the DBS related effects are believed to require neuroplasticity or reorganization and often take hours to months to observe. To optimize DBS parameters, it is therefore essential to develop electrophysiological biomarkers that characterize whether DBS settings are successfully engaging and modulating the network involved in the disease of interest. In this study, 10 individuals with drug resistant epilepsy undergoing intracranial stereotactic EEG including a thalamus electrode underwent a trial of repetitive thalamic stimulation. We evaluated thalamocortical effective connectivity using single pulse electrical stimulation, both at baseline and following a 145 Hz stimulation treatment trial. We found that when high frequency stimulation was delivered for >1.5 hours, the evoked potentials measured from remote regions were significantly reduced in amplitude and the degree of modulation was proportional to the strength of baseline connectivity. When stimulation was delivered for shorter time periods, results were more variable. These findings suggest that changes in effective connectivity in the network targeted with DBS accumulate over hours of DBS. Stimulation evoked potentials provide an electrophysiological biomarker that allows for efficient data-driven characterization of neuromodulation effects, which could enable new objective approaches for individualized DBS optimization.

Effective connectivity relates seizure outcome to electrode placement in responsive neurostimulation

Responsive neurostimulation is a closed-loop neuromodulation therapy for drug resistant focal epilepsy. Responsive neurostimulation electrodes are placed near ictal onset zones so as to enable detection of epileptiform activity and deliver electrical stimulation. There is no standard approach for determining the optimal placement of responsive neurostimulation electrodes. Clinicians make this determination based on presurgical tests, such as MRI, EEG, magnetoencephalography, ictal single-photon emission computed tomography and intracranial EEG. Currently functional connectivity measures are not being used in determining the placement of responsive neurostimulation electrodes. Cortico-cortical evoked potentials are a measure of effective functional connectivity. Cortico-cortical evoked potentials are generated by direct single-pulse electrical stimulation and can be used to investigate cortico-cortical connections in vivo. We hypothesized that the presence of high amplitude cortico-cortical evoked potentials, recorded during intracranial EEG monitoring, near the eventual responsive neurostimulation contact sites is predictive of better outcomes from its therapy. We retrospectively reviewed 12 patients in whom cortico-cortical evoked potentials were obtained during stereoelectroencephalography evaluation and subsequently underwent responsive neurostimulation therapy. We studied the relationship between cortico-cortical evoked potentials, the eventual responsive neurostimulation electrode locations and seizure reduction. Directional connectivity indicated by cortico-cortical evoked potentials can categorize stereoelectroencephalography electrodes as either receiver nodes/in-degree (an area of greater inward connectivity) or projection nodes/out-degree (greater outward connectivity). The follow-up period for seizure reduction ranged from 1.3-4.8 years (median 2.7) after responsive neurostimulation therapy started. Stereoelectroencephalography electrodes closest to the eventual responsive neurostimulation contact site tended to show larger in-degree cortico-cortical evoked potentials, especially for the early latency cortico-cortical evoked potentials period (10-60 ms period) in six out of 12 patients. Stereoelectroencephalography electrodes closest to the responsive neurostimulation contacts (≤5 mm) also had greater significant out-degree in the early cortico-cortical evoked potentials latency period than those further away (≥10 mm) (P < 0.05). Additionally, significant correlation was noted between in-degree cortico-cortical evoked potentials and greater seizure reduction with responsive neurostimulation therapy at its most effective period (P < 0.05). These findings suggest that functional connectivity determined by cortico-cortical evoked potentials may provide additional information that could help guide the optimal placement of responsive neurostimulation electrodes.

Localizing seizure onset zone by a cortico-cortical evoked potentials-based machine learning approach in focal epilepsy

OBJECTIVE

We aimed to develop a new approach for identifying the localization of the seizure onset zone (SOZ) based on corticocortical evoked potentials (CCEPs) and to compare the connectivity patterns in patients with different clinical phenotypes.

METHODS

Fifty patients who underwent stereoelectroencephalography and CCEP procedures were included. Logistic regression was used in the model, and six CCEP metrics were input as features: root mean square of the first peak (N1RMS) and second peak (N2RMS), peak latency, onset latency, width duration, and area.

RESULTS

The area under the curve (AUC) for localizing the SOZ ranged from 0.88 to 0.93. The N1RMS values in the hippocampus sclerosis (HS) group were greater than that of the focal cortical dysplasia (FCD) IIa group (p < 0.001), independent of the distance between the recorded and stimulated sites. The sensitivity of localization was higher in the seizure-free group than in the non-seizure-free group (p = 0.036).

CONCLUSIONS

This new method can be used to predict the SOZ localization in various focal epilepsy phenotypes.

SIGNIFICANCE

This study proposed a machine-learning approach for localizing the SOZ. Moreover, we examined how clinical phenotypes impact large-scale abnormality of the epileptogenic networks.

Single-pulse electrical stimulation artifact removal using the novel matching pursuit-based artifact reconstruction and removal method (MPARRM)

Objective.Single-pulse electrical stimulation (SPES) has been widely used to probe effective connectivity. However, analysis of the neural response is often confounded by stimulation artifacts. We developed a novel matching pursuit-based artifact reconstruction and removal method (MPARRM) capable of removing artifacts from stimulation-artifact-affected electrophysiological signals.Approach.To validate MPARRM across a wide range of potential stimulation artifact types, we performed a bench-top experiment in which we suspended electrodes in a saline solution to generate 110 types of real-world stimulation artifacts. We then added the generated stimulation artifacts to ground truth signals (stereoelectroencephalography signals from nine human subjects recorded during a receptive speech task), applied MPARRM to the combined signal, and compared the resultant denoised signal with the ground truth signal. We further applied MPARRM to artifact-affected neural signals recorded from the hippocampus while performing SPES on the ipsilateral basolateral amygdala in nine human subjects.Main results.MPARRM could remove stimulation artifacts without introducing spectral leakage or temporal spread. It accommodated variable stimulation parameters and recovered the early response to SPES within a wide range of frequency bands. Specifically, in the early response period (5-10 ms following stimulation onset), we found that the broadband gamma power (70-170 Hz) of the denoised signal was highly correlated with the ground truth signal (R=0.98±0.02, Pearson), and the broadband gamma activity of the denoised signal faithfully revealed the responses to the auditory stimuli within the ground truth signal with94%±1.47%sensitivity and99%±1.01%specificity. We further found that MPARRM could reveal the expected temporal progression of broadband gamma activity along the anterior-posterior axis of the hippocampus in response to the ipsilateral amygdala stimulation.Significance.MPARRM could faithfully remove SPES artifacts without confounding the electrophysiological signal components, especially during the early-response period. This method can facilitate the understanding of the neural response mechanisms of SPES.

Physiological and pathological neuronal connectivity in the living human brain based on intracranial EEG signals: the current state of research

Over the past decades, studies of human brain networks have received growing attention as the assessment and modelling of connectivity in the brain is a topic of high impact with potential application in the understanding of human brain organization under both physiological as well as various pathological conditions. Under specific diagnostic settings, human neuronal signal can be obtained from intracranial EEG (iEEG) recording in epilepsy patients that allows gaining insight into the functional organisation of living human brain. There are two approaches to assess brain connectivity in the iEEG-based signal: evaluation of spontaneous neuronal oscillations during ongoing physiological and pathological brain activity, and analysis of the electrophysiological cortico-cortical neuronal responses, evoked by single pulse electrical stimulation (SPES). Both methods have their own advantages and limitations. The paper outlines available methodological approaches and provides an overview of current findings in studies of physiological and pathological human brain networks, based on intracranial EEG recordings.

Direct cortical stimulation induces short-term plasticity of neural oscillations in humans

Patterned brain stimulation is commonly employed as a tool for eliciting plasticity in brain circuits and treating neuropsychiatric disorders. Although widely used in clinical settings, there remains a limited understanding of how stimulation-induced plasticity influences neural oscillations and their interplay with the underlying baseline functional architecture. To address this question, we applied 15 minutes of 10Hz focal electrical simulation, a pattern identical to 'excitatory' repetitive transcranial magnetic stimulation (rTMS), to 14 medically-intractable epilepsy patients undergoing intracranial electroencephalographic (iEEG). We quantified the spectral features of the cortico-cortical evoked potential (CCEPs) in these patients before and after stimulation. We hypothesized that for a given region the temporal and spectral components of the CCEP predicted the location and degree of stimulation-induced plasticity. Across patients, low frequency power (alpha and beta) showed the broadest change, while the magnitude of change was stronger in high frequencies (beta and gamma). Next we demonstrated that regions with stronger baseline evoked spectral responses were more likely to undergo plasticity after stimulation. These findings were specific to a given frequency in a specific temporal window. Post-stimulation power changes were driven by the interaction between direction of change in baseline power and temporal window of change. Finally, regions exhibiting early increases and late decreases in evoked baseline power exhibited power changes after stimulation and were independent of stimulation location. Together, these findings that time-frequency baseline features predict post-stimulation plasticity effects demonstrate properties akin to Hebbian learning in humans and extend this theory to the temporal and spectral window of interest. These findings can help improve our understanding of human brain plasticity and lead to more effective brain stimulation techniques.

Significance Statement

Brain stimulation is increasingly used to treat neuropsychiatric disorders by inducing changes in neural activity at specific brain regions. Despite their effectiveness, how these changes occur, specifically in the spectral domain, is unknown. To better understand how brain oscillations change after patterned stimulation, we performed focused stimulation in epilepsy patients and measured intracranial brain recordings. We found strong and predictable changes in brain oscillations (plasticity) after patterned stimulation. Specifically, low frequencies showing widespread effects and high frequencies exhibiting a greater magnitude of change. These changes were directly related to the temporal and spectral structure of brain responses prior to stimulation. Our study reveals that baseline brain activity patterns can predict how stimulation will induce plasticity in the spectral domain. These findings can help improve our understanding of human brain plasticity and lead to more effective brain stimulation techniques.

Highlights

We applied 15 minutes of repetitive 10Hz focal electrical stimulation and assessed the evoked brain-wide spectral changes with intracranial EEG.10Hz stimulation induced short-term plasticity in low frequency alpha evoked power broadly across regions and time windows and high frequency (beta, gamma) power specifically in early evoked time windows (10-50ms).Across patients, frequency bands, and time windows, brain regions with stronger baseline evoked power were more likely to undergo greater spectral changes after 10Hz stimulation.Post-stimulation spectral changes were specific; that is, for a given frequency band in a specific time window, baseline evoked power predicted post-stimulation change in the same frequency band and time window.Post-stimulation spectral change was driven by an interaction between direction of change and temporal window of baseline power; that is, regions exhibiting baseline evoked early (10-100ms) increases and late (100-200ms) decreases in power correlated with observed post-stimulation spectral changes.These results were independent of stimulation location.

Signatures of Electrical Stimulation Driven Network Interactions in the Human Limbic System

Stimulation-evoked signals are starting to be used as biomarkers to indicate the state and health of brain networks. The human limbic network, often targeted for brain stimulation therapy, is involved in emotion and memory processing. Previous anatomic, neurophysiological, and functional studies suggest distinct subsystems within the limbic network (Rolls, 2015). Studies using intracranial electrical stimulation, however, have emphasized the similarities of the evoked waveforms across the limbic network. We test whether these subsystems have distinct stimulation-driven signatures. In eight patients (four male, four female) with drug-resistant epilepsy, we stimulated the limbic system with single-pulse electrical stimulation. Reliable corticocortical evoked potentials (CCEPs) were measured between hippocampus and the posterior cingulate cortex (PCC) and between the amygdala and the anterior cingulate cortex (ACC). However, the CCEP waveform in the PCC after hippocampal stimulation showed a unique and reliable morphology, which we term the "limbic Hippocampus-Anterior nucleus of the thalamus-Posterior cingulate, HAP-wave." This limbic HAP-wave was visually distinct and separately decoded from the CCEP waveform in ACC after amygdala stimulation. Diffusion MRI data show that the measured end points in the PCC overlap with the end points of the parolfactory cingulum bundle rather than the parahippocampal cingulum, suggesting that the limbic HAP-wave may travel through fornix, mammillary bodies, and the anterior nucleus of the thalamus (ANT). This was further confirmed by stimulating the ANT, which evoked the same limbic HAP-wave but with an earlier latency. Limbic subsystems have unique stimulation-evoked signatures that may be used in the future to help network pathology diagnosis.SIGNIFICANCE STATEMENT The limbic system is often compromised in diverse clinical conditions, such as epilepsy or Alzheimer's disease, and characterizing its typical circuit responses may provide diagnostic insight. Stimulation-evoked waveforms have been used in the motor system to diagnose circuit pathology. We translate this framework to limbic subsystems using human intracranial stereo EEG (sEEG) recordings that measure deeper brain areas. Our sEEG recordings describe a stimulation-evoked waveform characteristic to the memory and spatial subsystem of the limbic network that we term the "limbic HAP-wave." The limbic HAP-wave follows anatomic white matter pathways from hippocampus to thalamus to the posterior cingulum and shows promise as a distinct biomarker of signaling in the human brain memory and spatial limbic network.

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.

An orexigenic subnetwork within the human hippocampus

Only recently have more specific circuit-probing techniques become available to inform previous reports implicating the rodent hippocampus in orexigenic appetitive processing1-4. This function has been reported to be mediated at least in part by lateral hypothalamic inputs, including those involving orexigenic lateral hypothalamic neuropeptides, such as melanin-concentrating hormone5,6. This circuit, however, remains elusive in humans. Here we combine tractography, intracranial electrophysiology, cortico-subcortical evoked potentials, and brain-clearing 3D histology to identify an orexigenic circuit involving the lateral hypothalamus and converging in a hippocampal subregion. We found that low-frequency power is modulated by sweet-fat food cues, and this modulation was specific to the dorsolateral hippocampus. Structural and functional analyses of this circuit in a human cohort exhibiting dysregulated eating behaviour revealed connectivity that was inversely related to body mass index. Collectively, this multimodal approach describes an orexigenic subnetwork within the human hippocampus implicated in obesity and related eating disorders.

Reliability and Validity of Transcranial Magnetic Stimulation-Electroencephalography Biomarkers

Noninvasive brain stimulation and neuroimaging have revolutionized human neuroscience with a multitude of applications, including diagnostic subtyping, treatment optimization, and relapse prediction. It is therefore particularly relevant to identify robust and clinically valuable brain biomarkers linking symptoms to their underlying neural mechanisms. Brain biomarkers must be reproducible (i.e., have internal reliability) across similar experiments within a laboratory and be generalizable (i.e., have external reliability) across experimental setups, laboratories, brain regions, and disease states. However, reliability (internal and external) is not alone sufficient; biomarkers also must have validity. Validity describes closeness to a true measure of the underlying neural signal or disease state. We propose that these metrics, reliability and validity, should be evaluated and optimized before any biomarker is used to inform treatment decisions. Here, we discuss these metrics with respect to causal brain connectivity biomarkers from coupling transcranial magnetic stimulation (TMS) with electroencephalography (EEG). We discuss controversies around TMS-EEG stemming from the multiple large off-target components (noise) and relatively weak genuine brain responses (signal), as is unfortunately often the case in noninvasive human neuroscience. We review the current state of TMS-EEG recordings, which consist of a mix of reliable noise and unreliable signal. We describe methods for evaluating TMS-EEG biomarkers, including how to assess internal and external reliability across facilities, cognitive states, brain networks, and disorders and how to validate these biomarkers using invasive neural recordings or treatment response. We provide recommendations to increase reliability and validity, discuss lessons learned, and suggest future directions for the field.

Quantitative approaches to guide epilepsy surgery from intracranial EEG

Over the past 10 years, the drive to improve outcomes from epilepsy surgery has stimulated widespread interest in methods to quantitatively guide epilepsy surgery from intracranial EEG (iEEG). Many patients fail to achieve seizure freedom, in part due to the challenges in subjective iEEG interpretation. To address this clinical need, quantitative iEEG analytics have been developed using a variety of approaches, spanning studies of seizures, interictal periods, and their transitions, and encompass a range of techniques including electrographic signal analysis, dynamical systems modeling, machine learning and graph theory. Unfortunately, many methods fail to generalize to new data and are sensitive to differences in pathology and electrode placement. Here, we critically review selected literature on computational methods of identifying the epileptogenic zone from iEEG. We highlight shared methodological challenges common to many studies in this field and propose ways that they can be addressed. One fundamental common pitfall is a lack of open-source, high-quality data, which we specifically address by sharing a centralized high-quality, well-annotated, multicentre dataset consisting of >100 patients to support larger and more rigorous studies. Ultimately, we provide a road map to help these tools reach clinical trials and hope to improve the lives of future patients.

Quantifying trial-by-trial variability during cortico-cortical evoked potential mapping of epileptogenic tissue

OBJECTIVE

Measuring cortico-cortical evoked potentials (CCEPs) is a promising tool for mapping epileptic networks, but it is not known how variability in brain state and stimulation technique might impact the use of CCEPs for epilepsy localization. We test the hypotheses that (1) CCEPs demonstrate systematic variability across trials and (2) CCEP amplitudes depend on the timing of stimulation with respect to endogenous, low-frequency oscillations.

METHODS

We studied 11 patients who underwent CCEP mapping after stereo-electroencephalography electrode implantation for surgical evaluation of drug-resistant epilepsy. Evoked potentials were measured from all electrodes after each pulse of a 30 s, 1 Hz bipolar stimulation train. We quantified monotonic trends, phase dependence, and standard deviation (SD) of N1 (15-50 ms post-stimulation) and N2 (50-300 ms post-stimulation) amplitudes across the 30 stimulation trials for each patient. We used linear regression to quantify the relationship between measures of CCEP variability and the clinical seizure-onset zone (SOZ) or interictal spike rates.

RESULTS

We found that N1 and N2 waveforms exhibited both positive and negative monotonic trends in amplitude across trials. SOZ electrodes and electrodes with high interictal spike rates had lower N1 and N2 amplitudes with higher SD across trials. Monotonic trends of N1 and N2 amplitude were more positive when stimulating from an area with higher interictal spike rate. We also found intermittent synchronization of trial-level N1 amplitude with low-frequency phase in the hippocampus, which did not localize the SOZ.

SIGNIFICANCE

These findings suggest that standard approaches for CCEP mapping, which involve computing a trial-averaged response over a .2-1 Hz stimulation train, may be masking inter-trial variability that localizes to epileptogenic tissue. We also found that CCEP N1 amplitudes synchronize with ongoing low-frequency oscillations in the hippocampus. Further targeted experiments are needed to determine whether phase-locked stimulation could have a role in localizing epileptogenic tissue.

Subsets of cortico-cortical evoked potentials propagate as traveling waves

Emerging evidence suggests that the temporal dynamics of cortico-cortical evoked potentials (CCEPs) may be used to characterize the patterns of information flow between and within brain networks. At present, however, the spatiotemporal dynamics of CCEP propagation cortically and subcortically are incompletely understood. We hypothesized that CCEPs propagate as an evoked traveling wave emanating from the site of stimulation. To elicit CCEPs, we applied single-pulse stimulation to stereoelectroencephalography (SEEG) electrodes implanted in 21 adult patients with intractable epilepsy. For each robust CCEP, we measured the timing of the maximal descent in evoked local field potentials and broadband high-gamma power (70-150 Hz) envelopes relative to the distance between the recording and stimulation contacts using three different metrics (i.e., Euclidean distance, path length, geodesic distance), representing direct, subcortical, and transcortical propagation, respectively. Many evoked responses to single-pulse electrical stimulation appear to propagate as traveling waves (~17-30%), even in the sparsely sampled, three-dimensional SEEG space. These results provide new insights into the spatiotemporal dynamics of CCEP propagation.

Mapping effective connectivity of human amygdala subdivisions with intracranial stimulation

The primate amygdala is a complex consisting of over a dozen nuclei that have been implicated in a host of cognitive functions, individual differences, and psychiatric illnesses. These functions are implemented through distinct connectivity profiles, which have been documented in animals but remain largely unknown in humans. Here we present results from 25 neurosurgical patients who had concurrent electrical stimulation of the amygdala with intracranial electroencephalography (electrical stimulation tract-tracing; es-TT), or fMRI (electrical stimulation fMRI; es-fMRI), methods providing strong inferences about effective connectivity of amygdala subdivisions with the rest of the brain. We quantified functional connectivity with medial and lateral amygdala, the temporal order of these connections on the timescale of milliseconds, and also detail second-order effective connectivity among the key nodes. These findings provide a uniquely detailed characterization of human amygdala functional connectivity that will inform functional neuroimaging studies in healthy and clinical populations.

Time-frequency signatures evoked by single-pulse deep brain stimulation to the subcallosal cingulate

Precision targeting of specific white matter bundles that traverse the subcallosal cingulate (SCC) has been linked to efficacy of deep brain stimulation (DBS) for treatment resistant depression (TRD). Methods to confirm optimal target engagement in this heterogenous region are now critical to establish an objective treatment protocol. As yet unexamined are the time-frequency features of the SCC evoked potential (SCC-EP), including spectral power and phase-clustering. We examined these spectral features—evoked power and phase clustering—in a sample of TRD patients (n = 8) with implanted SCC stimulators. Electroencephalogram (EEG) was recorded during wakeful rest. Location of electrical stimulation in the SCC target region was the experimental manipulation. EEG was analyzed at the surface level with an average reference for a cluster of frontal sensors and at a time window identified by prior study (50–150 ms). Morlet wavelets generated indices of evoked power and inter-trial phase clustering. Enhanced phase clustering at theta frequency (4–7 Hz) was observed in every subject and was significantly correlated with SCC-EP magnitude, but only during left SCC stimulation. Stimulation to dorsal SCC evinced stronger phase clustering than ventral SCC. There was a weak correlation between phase clustering and white matter density. An increase in evoked delta power (2–4 Hz) was also coincident with SCC-EP, but was less consistent across participants. DBS evoked time-frequency features index mm-scale changes to the location of stimulation in the SCC target region and correlate with structural characteristics implicated in treatment optimization. Results also imply a shared generative mechanism (inter-trial phase clustering) between evoked potentials evinced by electrical stimulation and evoked potentials evinced by auditory/visual stimuli and behavioral tasks. Understanding how current injection impacts downstream cortical activity is essential to building new technologies that adapt treatment parameters to individual differences in neurophysiology.

Anesthetic modulations dissociate neuroelectric characteristics between sensory-evoked and spontaneous activities across bilateral rat somatosensory cortical laminae

Spontaneous neural activity has been widely adopted to construct functional connectivity (FC) amongst distant brain regions. Although informative, the functional role and signaling mechanism of the resting state FC are not intuitive as those in stimulus/task-evoked activity. In order to bridge the gap, we investigated anesthetic modulation of both resting-state and sensory-evoked activities. We used two well-studied GABAergic anesthetics of varying dose (isoflurane: 0.5–2.0% and α-chloralose: 30 and 60 mg/kg∙h) and recorded changes in electrophysiology using a pair of laminar electrode arrays that encompass the entire depth of the bilateral somatosensory cortices (S1fl) in rats. Specifically, the study focused to describe how varying anesthesia conditions affect the resting state activities and resultant FC between bilateral hemispheres in comparison to those obtained by evoked responses. As results, isoflurane decreased the amplitude of evoked responses in a dose-dependent manner mostly due to the habituation of repetitive responses. However, α-chloralose rather intensified the amplitude without exhibiting habituation. No such diverging trend was observed for the spontaneous activity, in which both anesthetics increased the signal power. For α-chloralose, overall FC was similar to that obtained with the lowest dose of isoflurane at 0.5% while higher doses of isoflurane displayed increased FC. Interestingly, only α-chloralose elicited relatively much greater increases in the ipsi-stimulus evoked response (i.e., in S1fl ipsilateral to the stimulated forelimb) than those associated with the contra-stimulus response, suggesting enhanced neuronal excitability. Taken together, the findings demonstrate modulation of the FC profiles by anesthesia is highly non-linear, possibly with a distinct underlying mechanism that affects either resting state or evoked activities differently. Further, the current study warrants thorough investigation of the basal neuronal states prior to the interpretation of resting state FC and evoked activities for accurate understanding of neural signal processing and circuitry.

Movement-dependent electrical stimulation for volitional strengthening of cortical connections in behaving monkeys

Correlated activity of neurons can lead to long-term strengthening or weakening of the connections between them. In addition, the behavioral context, imparted by execution of physical movements or the presence of a reward, can modulate the plasticity induced by Hebbian mechanisms. In the present study, we have combined behavior and induced neuronal correlations to strengthen connections in the motor cortex of adult behaving monkeys. Correlated activity was induced using an electrical-conditioning protocol in which stimuli gated by voluntary movements were used to produce coactivation of neurons at motor-cortical sites involved in those movements. Delivery of movement-dependent stimulation resulted in small increases in the strength of associated cortical connections immediately after conditioning. Remarkably, when paired with further repetition of the movements that gated the conditioning stimuli, there were substantially larger gains in the strength of cortical connections, which occurred in a use-dependent manner, without delivery of additional conditioning stimulation. In the absence of such movements, little change was observed in the strength of motor-cortical connections. Performance of the motor behavior in the absence of conditioning also did not produce any changes in connectivity. Our results show that combining movement-gated stimulation with further natural use of the "conditioned" pathways after stimulation ends can produce use-dependent strengthening of connections in adult primates, highlighting an important role for behavior in cortical plasticity. Our data also provide strong support for combining movement-gated stimulation with use-dependent physical rehabilitation for strengthening connections weakened by a stroke or spinal cord injury.

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.

Difference in Network Effects of Pulsatile and Galvanic Stimulation

Biphasic pulsatile stimulation is the present standard for neural prosthetic use, and it is used to understand connectivity and functionality of the brain in brain mapping studies. While pulses have been shown to drive behavioral changes, such as biasing decision making, they have deficits. For example, cochlear implants restore hearing but lack the ability to restore pitch perception. Recent work shows that pulses produce artificial synchrony in networks of neurons and non-linear changes in firing rate with pulse amplitude. Studies also show galvanic stimulation, delivery of current for extended periods of time, produces more naturalistic behavioral responses than pulses. In this paper, we use a winner-take-all decision-making network model to investigate differences between pulsatile and galvanic stimulation at the single neuron and network level while accurately modeling the effects of pulses on neurons for the first time. Results show pulses bias spike timing and make neurons more resistive to natural network inputs than galvanic stimulation at an equivalent current amplitude. Clinical Relevance- This establishes that pulsatile stimulation may disrupt natural spike timing and network-level interactions while certain parameterizations of galvanic stimulation avoid these effects and can drive network firing more naturally.

Phase-based causality analysis with partial mutual information from mixed embedding

Instantaneous phases extracted from multivariate time series can retain information about the relationships between the underlying mechanisms that generate the series. Although phases have been widely used in the study of nondirectional coupling and connectivity, they have not found similar appeal in the study of causality. Herein, we present a new method for phase-based causality analysis, which combines ideas from the mixed embedding technique and the information-theoretic approach to causality in coupled oscillatory systems. We then use the introduced method to investigate causality in simulated datasets of bivariate, unidirectionally paired systems from combinations of Rössler, Lorenz, van der Pol, and Mackey-Glass equations. We observe that causality analysis using the phases can capture the true causal relation for coupling strength smaller than the analysis based on the amplitudes can capture. On the other hand, the causality estimation based on the phases tends to have larger variability, which is attributed more to the phase extraction process than the actual phase-based causality method. In addition, an application on real electroencephalographic data from an experiment on elicited human emotional states reinforces the usefulness of phases in causality identification.

Improved Speech Hearing in Noise with Invasive Electrical Brain Stimulation

Speech perception in noise is a challenging everyday task with which many listeners have difficulty. Here, we report a case in which electrical brain stimulation of implanted intracranial electrodes in the left planum temporale (PT) of a neurosurgical patient significantly and reliably improved subjective quality (up to 50%) and objective intelligibility (up to 97%) of speech in noise perception. Stimulation resulted in a selective enhancement of speech sounds compared with the background noises. The receptive fields of the PT sites whose stimulation improved speech perception were tuned to spectrally broad and rapidly changing sounds. Corticocortical evoked potential analysis revealed that the PT sites were located between the sites in Heschl's gyrus and the superior temporal gyrus. Moreover, the discriminability of speech from nonspeech sounds increased in population neural responses from Heschl's gyrus to the PT to the superior temporal gyrus sites. These findings causally implicate the PT in background noise suppression and may point to a novel potential neuroprosthetic solution to assist in the challenging task of speech perception in noise.SIGNIFICANCE STATEMENT Speech perception in noise remains a challenging task for many individuals. Here, we present a case in which the electrical brain stimulation of intracranially implanted electrodes in the planum temporale of a neurosurgical patient significantly improved both the subjective quality (up to 50%) and objective intelligibility (up to 97%) of speech perception in noise. Stimulation resulted in a selective enhancement of speech sounds compared with the background noises. Our local and network-level functional analyses placed the planum temporale sites in between the sites in the primary auditory areas in Heschl's gyrus and nonprimary auditory areas in the superior temporal gyrus. These findings causally implicate planum temporale in acoustic scene analysis and suggest potential neuroprosthetic applications to assist hearing in noise.

Simultaneous stereo-EEG and high-density scalp EEG recordings to study the effects of intracerebral stimulation parameters

BACKGROUND

Cortico-cortical evoked potentials (CCEPs) recorded by stereo-electroencephalography (SEEG) are a valuable tool to investigate brain reactivity and effective connectivity. However, invasive recordings are spatially sparse since they depend on clinical needs. This sparsity hampers systematic comparisons across-subjects, the detection of the whole-brain effects of intracortical stimulation, as well as their relationships to the EEG responses evoked by non-invasive stimuli.

OBJECTIVE

To demonstrate that CCEPs recorded by high-density electroencephalography (hd-EEG) provide additional information with respect SEEG alone and to provide an open, curated dataset to allow for further exploration of their potential.

METHODS

The dataset encompasses SEEG and hd-EEG recordings simultaneously acquired during Single Pulse Electrical Stimulation (SPES) in drug-resistant epileptic patients (N = 36) in whom stimulations were delivered with different physical, geometrical, and topological parameters. Differences in CCEPs were assessed by amplitude, latency, and spectral measures.

RESULTS

While invasively and non-invasively recorded CCEPs were generally correlated, differences in pulse duration, angle and stimulated cortical area were better captured by hd-EEG. Further, intracranial stimulation evoked site-specific hd-EEG responses that reproduced the spectral features of EEG responses to transcranial magnetic stimulation (TMS). Notably, SPES, albeit unperceived by subjects, elicited scalp responses that were up to one order of magnitude larger than the responses typically evoked by sensory stimulation in awake humans.

CONCLUSIONS

CCEPs can be simultaneously recorded with SEEG and hd-EEG and the latter provides a reliable descriptor of the effects of SPES as well as a common reference to compare the whole-brain effects of intracortical stimulation to those of non-invasive transcranial or sensory stimulations in humans.

Intracranial electroencephalographic biomarker predicts effective responsive neurostimulation for epilepsy prior to treatment

OBJECTIVE

Despite the overall success of responsive neurostimulation (RNS) therapy for drug-resistant focal epilepsy, clinical outcomes in individuals vary significantly and are hard to predict. Biomarkers that indicate the clinical efficacy of RNS-ideally before device implantation-are critically needed, but challenges include the intrinsic heterogeneity of the RNS patient population and variability in clinical management across epilepsy centers. The aim of this study is to use a multicenter dataset to evaluate a candidate biomarker from intracranial electroencephalographic (iEEG) recordings that predicts clinical outcome with subsequent RNS therapy.

METHODS

We assembled a federated dataset of iEEG recordings, collected prior to RNS implantation, from a retrospective cohort of 30 patients across three major epilepsy centers. Using ictal iEEG recordings, each center independently calculated network synchronizability, a candidate biomarker indicating the susceptibility of epileptic brain networks to RNS therapy.

RESULTS

Ictal measures of synchronizability in the high-γ band (95-105 Hz) significantly distinguish between good and poor RNS responders after at least 3 years of therapy under the current RNS therapy guidelines (area under the curve = .83). Additionally, ictal high-γ synchronizability is inversely associated with the degree of therapeutic response.

SIGNIFICANCE

This study provides a proof-of-concept roadmap for collaborative biomarker evaluation in federated data, where practical considerations impede full data sharing across centers. Our results suggest that network synchronizability can help predict therapeutic response to RNS therapy. With further validation, this biomarker could facilitate patient selection and help avert a costly, invasive intervention in patients who are unlikely to benefit.

Local and distant cortical responses to single pulse intracranial stimulation in the human brain are differentially modulated by specific stimulation parameters

BACKGROUND

Electrical neuromodulation via direct electrical stimulation (DES) is an increasingly common therapy for a wide variety of neuropsychiatric diseases. Unfortunately, therapeutic efficacy is inconsistent, likely due to our limited understanding of the relationship between the massive stimulation parameter space and brain tissue responses.

OBJECTIVE

To better understand how different parameters induce varied neural responses, we systematically examined single pulse-induced cortico-cortico evoked potentials (CCEP) as a function of stimulation amplitude, duration, brain region, and whether grey or white matter was stimulated.

METHODS

We measured voltage peak amplitudes and area under the curve (AUC) of intracranially recorded stimulation responses as a function of distance from the stimulation site, pulse width, current injected, location relative to grey and white matter, and brain region stimulated (N = 52, n = 719 stimulation sites).

RESULTS

Increasing stimulation pulse width increased responses near the stimulation location. Increasing stimulation amplitude (current) increased both evoked amplitudes and AUC nonlinearly. Locally (<15 mm), stimulation at the boundary between grey and white matter induced larger responses. In contrast, for distant sites (>15 mm), white matter stimulation consistently produced larger responses than stimulation in or near grey matter. The stimulation location-response curves followed different trends for cingulate, lateral frontal, and lateral temporal cortical stimulation.

CONCLUSION

These results demonstrate that a stronger local response may require stimulation in the grey-white boundary while stimulation in the white matter could be needed for network activation. Thus, stimulation parameters tailored for a specific anatomical-functional outcome may be key to advancing neuromodulatory therapy.

Neuromonitoring of the language pathways using cortico-cortical evoked potentials: a systematic review and meta-analysis

Cortico-cortical evoked potentials (CCEPs) are a surge in activity of one cortical zone caused by stimulation of another cortical zone. Recording of CCEP may be a useful method of intraoperative monitoring of the brain pathways, particularly of the language-related tracts. We aimed to conduct a systematic review and meta-analysis, dedicated to the clinical question: Does the CCEP recording effectively predict the postoperative speech deficits in neurosurgical patients? We conducted language-restricted PubMed, Google Scholar, Scopus, and Cochrane database search for eligible studies of CCEP published until March 2021. There were 4 articles (3 case series and 1 case report), which met our inclusion/exclusion criteria. A total of 32 patients (30 cases of tumors and 2 cavernomas) included in the analysis were divided into two cohorts - quantitative and qualitative, in accordance with the method of evaluating changes in the amplitude of CCEP after the lesion resection and postoperative alterations in speech function. Quantitative variables were studied using the Spearman rank correlation coefficient. Categorical variables were compared in groups by Fisher's exact test. We found a strong positive correlation between the decrease in the N1 wave amplitude and the severity of postoperative speech deficits (quantitative cohort: r = 0.57, p = 0.01; qualitative cohort: p = 0.02). Thus, the CCEP method using the N1 wave amplitude as a marker enables to effectively predict postoperative speech outcomes. Nevertheless, the low level of evidence for the included works indicated the necessity for additional research on this issue.

Network connectivity predicts effectiveness of responsive neurostimulation in focal epilepsy

Responsive neurostimulation is a promising treatment for drug-resistant focal epilepsy; however, clinical outcomes are highly variable across individuals. The therapeutic mechanism of responsive neurostimulation likely involves modulatory effects on brain networks; however, with no known biomarkers that predict clinical response, patient selection remains empiric. This study aimed to determine whether functional brain connectivity measured non-invasively prior to device implantation predicts clinical response to responsive neurostimulation therapy. Resting-state magnetoencephalography was obtained in 31 participants with subsequent responsive neurostimulation device implantation between 15 August 2014 and 1 October 2020. Functional connectivity was computed across multiple spatial scales (global, hemispheric, and lobar) using pre-implantation magnetoencephalography and normalized to maps of healthy controls. Normalized functional connectivity was investigated as a predictor of clinical response, defined as percent change in self-reported seizure frequency in the most recent year of clinic visits relative to pre-responsive neurostimulation baseline. Area under the receiver operating characteristic curve quantified the performance of functional connectivity in predicting responders (≥50% reduction in seizure frequency) and non-responders (<50%). Leave-one-out cross-validation was furthermore performed to characterize model performance. The relationship between seizure frequency reduction and frequency-specific functional connectivity was further assessed as a continuous measure. Across participants, stimulation was enabled for a median duration of 52.2 (interquartile range, 27.0-62.3) months. Demographics, seizure characteristics, and responsive neurostimulation lead configurations were matched across 22 responders and 9 non-responders. Global functional connectivity in the alpha and beta bands were lower in non-responders as compared with responders (alpha, p_fdr < 0.001; beta, p_fdr < 0.001). The classification of responsive neurostimulation outcome was improved by combining feature inputs; the best model incorporated four features (i.e. mean and dispersion of alpha and beta bands) and yielded an area under the receiver operating characteristic curve of 0.970 (0.919-1.00). The leave-one-out cross-validation analysis of this four-feature model yielded a sensitivity of 86.3%, specificity of 77.8%, positive predictive value of 90.5%, and negative predictive value of 70%. Global functional connectivity in alpha band correlated with seizure frequency reduction (alpha, P = 0.010). Global functional connectivity predicted responder status more strongly, as compared with hemispheric predictors. Lobar functional connectivity was not a predictor. These findings suggest that non-invasive functional connectivity may be a candidate personalized biomarker that has the potential to predict responsive neurostimulation effectiveness and to identify patients most likely to benefit from responsive neurostimulation therapy. Follow-up large-cohort, prospective studies are required to validate this biomarker. These findings furthermore support an emerging view that the therapeutic mechanism of responsive neurostimulation involves network-level effects in the brain.

Physiologically informed neuromodulation

The rapid evolution of neuromodulation techniques includes an increasing amount of research into stimulation paradigms that are guided by patients' neurophysiology, to increase efficacy and responder rates. Treatment personalisation and target engagement have shown to be effective in fields such as Parkinson's disease, and closed-loop paradigms have been successfully implemented in cardiac defibrillators. Promising avenues are being explored for physiologically informed neuromodulation in psychiatry. Matching the stimulation frequency to individual brain rhythms has shown some promise in transcranial magnetic stimulation (TMS). Matching the phase of those rhythms may further enhance neuroplasticity, for instance when combining TMS with electroencephalographic (EEG) recordings. Resting-state EEG and event-related potentials may be useful to demonstrate connectivity between stimulation sites and connected areas. These techniques are available today to the psychiatrist to diagnose underlying sleep disorders, epilepsy, or lesions as contributing factors to the cause of depression. These technologies may also be useful in assessing the patient's brain network status prior to deciding on treatment options. Ongoing research using invasive recordings may allow for future identification of mood biomarkers and network structure. A core limitation is that biomarker research may currently be limited by the internal heterogeneity of psychiatric disorders according to the current DSM-based classifications. New approaches are being developed and may soon be validated. Finally, care must be taken when incorporating closed-loop capabilities into neuromodulation systems, by ensuring the safe operation of the system and understanding the physiological dynamics. Neurophysiological tools are rapidly evolving and will likely define the next generation of neuromodulation therapies.

Six-dimensional dynamic tractography atlas of language connectivity in the developing brain

During a verbal conversation, our brain moves through a series of complex linguistic processing stages: sound decoding, semantic comprehension, retrieval of semantically coherent words, and overt production of speech outputs. Each process is thought to be supported by a network consisting of local and long-range connections bridging between major cortical areas. Both temporal and extratemporal lobe regions have functional compartments responsible for distinct language domains, including the perception and production of phonological and semantic components. This study provides quantitative evidence of how directly connected inter-lobar neocortical networks support distinct stages of linguistic processing across brain development. Novel six-dimensional tractography was used to intuitively visualize the strength and temporal dynamics of direct inter-lobar effective connectivity between cortical areas activated during each linguistic processing stage. We analysed 3401 non-epileptic intracranial electrode sites from 37 children with focal epilepsy (aged 5-20 years) who underwent extra-operative electrocorticography recording. Principal component analysis of auditory naming-related high-gamma modulations determined the relative involvement of each cortical area during each linguistic processing stage. To quantify direct effective connectivity, we delivered single-pulse electrical stimulation to 488 temporal and 1581 extratemporal lobe sites and measured the early cortico-cortical spectral responses at distant electrodes. Mixed model analyses determined the effects of naming-related high-gamma co-augmentation between connecting regions, age, and cerebral hemisphere on the strength of effective connectivity independent of epilepsy-related factors. Direct effective connectivity was strongest between extratemporal and temporal lobe site pairs, which were simultaneously activated between sentence offset and verbal response onset (i.e. response preparation period); this connectivity was approximately twice more robust than that with temporal lobe sites activated during stimulus listening or overt response. Conversely, extratemporal lobe sites activated during overt response were equally connected with temporal lobe language sites. Older age was associated with increased strength of inter-lobar effective connectivity especially between those activated during response preparation. The arcuate fasciculus supported approximately two-thirds of the direct effective connectivity pathways from temporal to extratemporal auditory language-related areas but only up to half of those in the opposite direction. The uncinate fasciculus consisted of <2% of those in the temporal-to-extratemporal direction and up to 6% of those in the opposite direction. We, for the first time, provided an atlas which quantifies and animates the strength, dynamics, and direction specificity of inter-lobar neural communications between language areas via the white matter pathways. Language-related effective connectivity may be strengthened in an age-dependent manner even after the age of 5.

Neurobehavioural comorbidities of epilepsy: towards a network-based precision taxonomy

Cognitive and behavioural comorbidities are prevalent in childhood and adult epilepsies and impose a substantial human and economic burden. Over the past century, the classic approach to understanding the aetiology and course of these comorbidities has been through the prism of the medical taxonomy of epilepsy, including its causes, course, characteristics and syndromes. Although this 'lesion model' has long served as the organizing paradigm for the field, substantial challenges to this model have accumulated from diverse sources, including neuroimaging, neuropathology, neuropsychology and network science. Advances in patient stratification and phenotyping point towards a new taxonomy for the cognitive and behavioural comorbidities of epilepsy, which reflects the heterogeneity of their clinical presentation and raises the possibility of a precision medicine approach. As we discuss in this Review, these advances are informing the development of a revised aetiological paradigm that incorporates sophisticated neurobiological measures, genomics, comorbid disease, diversity and adversity, and resilience factors. We describe modifiable risk factors that could guide early identification, treatment and, ultimately, prevention of cognitive and broader neurobehavioural comorbidities in epilepsy and propose a road map to guide future research.

Temporal order of signal propagation within and across intrinsic brain networks

We studied the temporal dynamics of activity within and across functional MRI (fMRI)-derived nodes of intrinsic resting-state networks of the human brain using intracranial electroencephalography (iEEG) and repeated single-pulse electrical stimulation (SPES) in neurosurgical subjects implanted with intracranial electrodes. We stimulated and recorded from 2,133 and 2,372 sites, respectively, in 29 subjects. We found that N1 and N2 segments of the evoked responses are associated with intra- and internetwork communications, respectively. In a separate cognitive experiment, evoked electrophysiological responses to visual target stimuli occurred with less temporal separation across pairs of electrodes that were located within the same fMRI-defined resting-state networks compared with those located across different resting-state networks. Our results suggest intranetwork prior to internetwork information processing at the subsecond timescale.

A brain atlas of axonal and synaptic delays based on modelling of cortico-cortical evoked potentials

Epilepsy presurgical investigation may include focal intracortical single-pulse electrical stimulations with depth electrodes, which induce cortico-cortical evoked potentials at distant sites because of white matter connectivity. Cortico-cortical evoked potentials provide a unique window on functional brain networks because they contain sufficient information to infer dynamical properties of large-scale brain connectivity, such as preferred directionality and propagation latencies.

 Here, we developed a biologically informed modelling approach to estimate the neural physiological parameters of brain functional networks from the cortico-cortical evoked potentials recorded in a large multicentric database. Specifically, we considered each cortico-cortical evoked potential as the output of a transient stimulus entering the stimulated region, which directly propagated to the recording region. Both regions were modelled as coupled neural mass models, the parameters of which were estimated from the first cortico-cortical evoked potential component, occurring before 80 ms, using dynamic causal modelling and Bayesian model inversion. This methodology was applied to the data of 780 patients with epilepsy from the F-TRACT database, providing a total of 34 354 bipolar stimulations and 774 445 cortico-cortical evoked potentials. The cortical mapping of the local excitatory and inhibitory synaptic time constants and of the axonal conduction delays between cortical regions was obtained at the population level using anatomy-based averaging procedures, based on the Lausanne2008 and the HCP-MMP1 parcellation schemes, containing 130 and 360 parcels, respectively. To rule out brain maturation effects, a separate analysis was performed for older (>15 years) and younger patients (<15 years). In the group of older subjects, we found that the cortico-cortical axonal conduction delays between parcels were globally short (median = 10.2 ms) and only 16% were larger than 20 ms. This was associated to a median velocity of 3.9 m/s. Although a general lengthening of these delays with the distance between the stimulating and recording contacts was observed across the cortex, some regions were less affected by this rule, such as the insula for which almost all efferent and afferent connections were faster than 10 ms. Synaptic time constants were found to be shorter in the sensorimotor, medial occipital and latero-temporal regions, than in other cortical areas. Finally, we found that axonal conduction delays were significantly larger in the group of subjects younger than 15 years, which corroborates that brain maturation increases the speed of brain dynamics.

To our knowledge, this study is the first to provide a local estimation of axonal conduction delays and synaptic time constants across the whole human cortex in vivo, based on intracerebral electrophysiological recordings.

Electrocorticography reveals thalamic control of cortical dynamics following traumatic brain injury

The return of consciousness after traumatic brain injury (TBI) is associated with restoring complex cortical dynamics; however, it is unclear what interactions govern these complex dynamics. Here, we set out to uncover the mechanism underlying the return of consciousness by measuring local field potentials (LFP) using invasive electrophysiological recordings in patients recovering from TBI. We found that injury to the thalamus, and its efferent projections, on MRI were associated with repetitive and low complexity LFP signals from a highly structured phase space, resembling a low-dimensional ring attractor. But why do thalamic injuries in TBI patients result in a cortical attractor? We built a simplified thalamocortical model, which connotes that thalamic input facilitates the formation of cortical ensembles required for the return of cognitive function and the content of consciousness. These observations collectively support the view that thalamic input to the cortex enables rich cortical dynamics associated with consciousness.

Closed-loop neuromodulation in an individual with treatment-resistant depression

Deep brain stimulation is a promising treatment for neuropsychiatric conditions such as major depression. It could be optimized by identifying neural biomarkers that trigger therapy selectively when symptom severity is elevated. We developed an approach that first used multi-day intracranial electrophysiology and focal electrical stimulation to identify a personalized symptom-specific biomarker and a treatment location where stimulation improved symptoms. We then implanted a chronic deep brain sensing and stimulation device and implemented a biomarker-driven closed-loop therapy in an individual with depression. Closed-loop therapy resulted in a rapid and sustained improvement in depression. Future work is required to determine if the results and approach of this n-of-1 study generalize to a broader population.

Effective connectivity differs between focal cortical dysplasia types I and II

OBJECTIVE

To determine whether brain connectivity differs between focal cortical dysplasia (FCD) types I and II.

METHODS

We compared cortico-cortical evoked potentials (CCEPs) as measures of effective brain connectivity in 25 FCD patients with drug-resistant focal epilepsy who underwent intracranial evaluation with stereo-electroencephalography (SEEG). We analyzed the amplitude and latency of CCEP responses following ictal-onset single-pulse electrical stimulation (iSPES).

RESULTS

In comparison to FCD type II, patients with type I demonstrated significantly larger responses in the electrodes near the ictal-onset zone (<50 mm). These findings persisted when controlling for the location of the epileptogenic zone, as noted in patients with temporal lobe epilepsies, as well as controlling for seizure type, as noted in patients with focal to bilateral tonic-clonic seizures (FBTCS). In type II, the root mean square (RMS) of CCEP responses dropped substantially from the early segment (10-60 ms) to the middle and late segments (60-600 ms). The middle and late CCEP latency segments showed the largest differences between FCD types I and II.

SIGNIFICANCE

Focal cortical dysplasia type I may have a greater degree of cortical hyperexcitability as compared with FCD type II. In addition, FCD type II displays a more restrictive area of hyperexcitability in both temporal and spatial domains. In patients with FBTCS and type I FCD, the increased amplitudes of RMS in the middle and late CCEP periods appear consistent with the cortico-thalamo-cortical network involvement of FBTCS. The notable differences in degree and extent of hyperexcitability may contribute to the different postsurgical seizure outcomes noted between these two pathological substrates.

Basis profile curve identification to understand electrical stimulation effects in human brain networks

Brain networks can be explored by delivering brief pulses of electrical current in one area while measuring voltage responses in other areas. We propose a convergent paradigm to study brain dynamics, focusing on a single brain site to observe the average effect of stimulating each of many other brain sites. Viewed in this manner, visually-apparent motifs in the temporal response shape emerge from adjacent stimulation sites. This work constructs and illustrates a data-driven approach to determine characteristic spatiotemporal structure in these response shapes, summarized by a set of unique "basis profile curves" (BPCs). Each BPC may be mapped back to underlying anatomy in a natural way, quantifying projection strength from each stimulation site using simple metrics. Our technique is demonstrated for an array of implanted brain surface electrodes in a human patient. This framework enables straightforward interpretation of single-pulse brain stimulation data, and can be applied generically to explore the diverse milieu of interactions that comprise the connectome.

Impact of gyral geometry on cortical responses to surface electrical stimulation: insights from experimental and modeling studies

Objective.Invasive simultaneous stimulation and recording from intracranial electrodes and microwire arrays were used to investigate direct cortical responses to single pulses of electrical stimulation in humans.Approach.Microwire contacts measured surface potentials in cortical microdomains at a distance of 2-6 mm from the intracranial electrode. Direct cortical responses to stimulation (<20 ms) consisted of a larger surface negative potentials.Main results. The latencies of these responses were directly or inversely correlated with distances between the intracranial electrode and microwire contacts. We hypothesize that surface negative potentials reflected local synchronous depolarization of apical dendrites of pyramidal neurons in cortical microdomains in the superficial cortical layer and resulted from the activation of gray matter axons that delivered excitatory inputs to apical dendrites after cortical stimulation. We further hypothesized that the positive or inverse distance-latency correlations of the recorded negative responses were measured depending on whether activation of neurons originated at one (crown) or multiple (crown, lip, bank) sites throughout the gyrus simultaneously. The inverse distance-latency correlations then reflected the spatiotemporal superposition of different nearby sources of neuronal recruitment in the gyrus. To prove this hypothesis, we built an anatomically informed and biophysically realistic cortical network model and simulated early responses of cortical neurons to electrical stimulation in this cortical network model. The model simulations yielded negative potentials in simulated microdomains in the cortical model consistent with those recorded from humans. The model predicted sensitivity of cortical responses to the alignment of the stimulating electrode and microwire array with respect to the cortical gyrus and confirmed that gyral geometry has a major impact on direct neuronal recruitment, the timing, and the time course of neuronal activation in cortical microdomains.Significance.In this work, we demonstrated how the high-resolution forward network models can be used for better understanding and detailed prediction of cortical stimulation effects. Accurate predictive modeling tools are needed for the progress of brain stimulation therapies.

Intracerebral electrical stimulations of the temporal lobe: A stereoelectroencephalography study

The functional anatomy of the anteromesial portion of the temporal lobe and its involvement in epilepsy can be explored by means of intracerebral electrical stimulations. Here, we aimed to expand the knowledge of its physiological and pathophysiological symptoms by conducting the first large-sample systematic analysis of 1529 electrical stimulations of this anatomical region. We retrospectively analysed all clinical manifestations induced by intracerebral electrical stimulations in 173 patients with drug-resistant focal epilepsy with at least one electrode implanted in this area. We found that high-frequency stimulations were more likely to evoke electroclinical manifestations (p < .0001) and also provoked 'false positive' seizures. Multimodal symptoms were associated with EEG electrical modification (after discharge) (p < .0001). Visual symptoms were not associated with after discharge (p = .0002) and were mainly evoked by stimulation of the hippocampus (p = .009) and of the parahippocampal gyrus (p = .0212). 'False positive seizures' can be evoked by stimulation of the hippocampus, parahippocampal gyrus and amygdala, likely due to their intrinsic low epileptogenic threshold. Visual symptoms evoked in the hippocampus and parahippocampal gyrus, without EEG changes, are physiological symptoms and suggest involvement of these areas in the visual ventral stream. Our findings provide meaningful guidance in the interpretation of intracranial EEG studies of the temporal lobe.