Feasibility of Stereo EEG Based Brain Computer Interfacing in An Adult and Pediatric Cohort.

Introduction: Stereoelectroencephalography (sEEG) is a mesoscale intracranial monitoring method which records from the brain volumetrically with depth electrodes. Implementation of sEEG in BCI has not been well-described across a diverse patient cohort. Methods: Across eighteen subjects, channels with high frequency broadband (HFB, 65-115Hz) power increases during hand, tongue, or foot movements during a motor screening task were provided real-time feedback based on these HFB power changes to control a cursor on a screen. Results: Seventeen subjects established successful control of the overt motor BCI, but only nine were able to control imagery BCI with ≥ 80% accuracy. In successful imagery BCI, HFB power in the two target conditions separated into distinct subpopulations, which appear to engage unique subnetworks of the motor cortex compared to cued movement or imagery alone. Conclusion: sEEG-based motor BCI utilizing overt movement and kinesthetic imagery is robust across patient ages and cortical regions with substantial differences in learning proficiency between real or imagined movement.

Distributed representations of prediction error signals across the cortical hierarchy are synergistic.

A relevant question concerning inter-areal communication in the cortex is whether these interactions are synergistic. Synergy refers to the complementary effect of multiple brain signals conveying more information than the sum of each isolated signal. Redundancy, on the other hand, refers to the common information shared between brain signals. Here, we dissociated cortical interactions encoding complementary information (synergy) from those sharing common information (redundancy) during prediction error (PE) processing. We analyzed auditory and frontal electrocorticography (ECoG) signals in five common awake marmosets performing two distinct auditory oddball tasks and investigated to what extent event-related potentials (ERP) and broadband (BB) dynamics encoded synergistic and redundant information about PE processing. The information conveyed by ERPs and BB signals was synergistic even at lower stages of the hierarchy in the auditory cortex and between auditory and frontal regions. Using a brain-constrained neural network, we simulated the synergy and redundancy observed in the experimental results and demonstrated that the emergence of synergy between auditory and frontal regions requires the presence of strong, long-distance, feedback, and feedforward connections. These results indicate that distributed representations of PE signals across the cortical hierarchy can be highly synergistic.

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.

Method for cycle detection in sparse, irregularly sampled, long-term neuro-behavioral timeseries: Basis pursuit denoising with polynomial detrending of long-term, inter-ictal epileptiform activity.

Numerous physiological processes are cyclical, but sampling these processes densely enough to perform frequency decomposition and subsequent analyses can be challenging. Mathematical approaches for decomposition and reconstruction of sparsely and irregularly sampled signals are well established but have been under-utilized in physiological applications. We developed a basis pursuit denoising with polynomial detrending (BPWP) model that recovers oscillations and trends from sparse and irregularly sampled timeseries. We validated this model on a unique dataset of long-term inter-ictal epileptiform discharge (IED) rates from human hippocampus recorded with a novel investigational device with continuous local field potential sensing. IED rates have well established circadian and multiday cycles related to sleep, wakefulness, and seizure clusters. Given sparse and irregular samples of IED rates from multi-month intracranial EEG recordings from ambulatory humans, we used BPWP to compute narrowband spectral power and polynomial trend coefficients and identify IED rate cycles in three subjects. In select cases, we propose that random and irregular sampling may be leveraged for frequency decomposition of physiological signals. Trial Registration: NCT03946618.

Functional Mapping of Movement and Speech Using Task-Based Electrophysiological Changes in Stereoelectroencephalography.

Introduction: Stereoelectroencephalography (sEEG) has become the predominant method for intracranial seizure localization. When imaging, semiology, and scalp EEG are not in full agreement or definitively localizing, implanted sEEG recordings are used to test candidate seizure onset zones (SOZs). Discovered SOZs may then be targeted for resection, laser ablation, or neurostimulation. If a SOZ is eloquent, resection and ablation are both contraindicated, so identifying functional representation is crucial for therapeutic decision making. Objective: We present a novel functional brain mapping technique that utilizes task-based electrophysiological changes in sEEG during behavioral tasks and test this in pediatric and adult patients. Methods: sEEG was recorded in twenty patients with epilepsy, aged 6-39 (12 female, 18 of 20 patients < 21 years old), who underwent implanted monitoring to identify seizure onset. Each performed 1) visually cued simple repetitive movements of the hand, foot, or tongue while electromyography was recorded, and 2) simple picture naming or verb generation speech tasks while audio was recorded. Broadband changes in the power spectrum of the sEEG were compared between behavior and rest. Results: Electrophysiological functional mapping of movement and/or speech areas was completed in all 20 patients. Eloquent representation was identified in both cortex and white matter, and generally corresponded to classically described functional anatomic organization as well as other clinical mapping results. Robust maps of brain activity were identified in healthy brain, regions of developmental or acquired structural abnormality, and SOZs. Conclusion: Task based electrophysiological mapping using broadband changes in the sEEG signal reliably identifies movement and speech representation in pediatric and adult epilepsy patients.

Mapping of the central sulcus using non-invasive ultra-high-density brain recordings.

Brain mapping is vital in understanding the brain's functional organization. Electroencephalography (EEG) is one of the most widely used brain mapping approaches, primarily because it is non-invasive, inexpensive, straightforward, and effective. Increasing the electrode density in EEG systems provides more neural information and can thereby enable more detailed and nuanced mapping procedures. Here, we show that the central sulcus can be clearly delineated using a novel ultra-high-density EEG system (uHD EEG) and somatosensory evoked potentials (SSEPs). This uHD EEG records from 256 channels with an inter-electrode distance of 8.6 mm and an electrode diameter of 5.9 mm. Reconstructed head models were generated from T1-weighted MRI scans, and electrode positions were co-registered to these models to create topographical plots of brain activity. EEG data were first analyzed with peak detection methods and then classified using unsupervised spectral clustering. Our topography plots of the spatial distribution from the SSEPs clearly delineate a division between channels above the somatosensory and motor cortex, thereby localizing the central sulcus. Individual EEG channels could be correctly classified as anterior or posterior to the central sulcus with 95.2% accuracy, which is comparable to accuracies from invasive intracranial recordings. Our findings demonstrate that uHD EEG can resolve the electrophysiological signatures of functional representation in the brain at a level previously only seen from surgically implanted electrodes. This novel approach could benefit numerous applications, including research, neurosurgical mapping, clinical monitoring, detection of conscious function, brain-computer interfacing (BCI), rehabilitation, and mental health.

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.

Clinical evaluation of a stereotactic system for single-stage deep brain stimulation surgery under general anesthesia: technical note.

OBJECTIVE: Conventional frame-based stereotactic systems have circumferential base frames, often necessitating deep brain stimulation (DBS) surgery in two stages: intracranial electrode insertion followed by surgical re-preparation and pulse generator implantation. Some patients do not tolerate awake surgery, underscoring the need for a safe alternative for asleep DBS surgery. A frame-based stereotactic system with a skull-mounted "key" in lieu of a circumferential base frame received US FDA clearance. The authors describe the system's application for single-stage, asleep DBS surgery in 8 patients at their institution and review its workflow and technical considerations. METHODS: Eight patients underwent DBS lead insertion and IPG implantation in a single surgical preparation under general anesthesia using the system. Postoperative CT imaging confirmed lead placement. RESULTS: Eight patients underwent implantation of 15 total leads targeting the ventral intermediate nucleus (4 patients), globus pallidus internus (GPi; 3 patients), and subthalamic nucleus (STN; 1 patient). Intraoperative microelectrode recording was conducted for GPi and STN targets. Postoperative CT imaging revealed a mean ± SD radial error of 1.24 ± 0.45 mm (n = 15 leads), without surgical complications. CONCLUSIONS: The stereotactic system facilitated safe and effective asleep, single-stage DBS surgery, maintaining traditional lead accuracy standards.

Acute to long-term characteristics of impedance recordings during neurostimulation in humans.

Objective.This study aims to characterize the time course of impedance, a crucial electrophysiological property of brain tissue, in the human thalamus (THL), amygdala-hippocampus, and posterior hippocampus over an extended period.Approach.Impedance was periodically sampled every 5-15 min over several months in five subjects with drug-resistant epilepsy using an investigational neuromodulation device. Initially, we employed descriptive piecewise and continuous mathematical models to characterize the impedance response for approximately three weeks post-electrode implantation. We then explored the temporal dynamics of impedance during periods when electrical stimulation was temporarily halted, observing a monotonic increase (rebound) in impedance before it stabilized at a higher value. Lastly, we assessed the stability of amplitude and phase over the 24 h impedance cycle throughout the multi-month recording.Main results.Immediately post-implantation, the impedance decreased, reaching a minimum value in all brain regions within approximately two days, and then increased monotonically over about 14 d to a stable value. The models accounted for the variance in short-term impedance changes. Notably, the minimum impedance of the THL in the most epileptogenic hemisphere was significantly lower than in other regions. During the gaps in electrical stimulation, the impedance rebound decreased over time and stabilized around 200 days post-implant, likely indicative of the foreign body response and fibrous tissue encapsulation around the electrodes. The amplitude and phase of the 24 h impedance oscillation remained stable throughout the multi-month recording, with circadian variation in impedance dominating the long-term measures.Significance.Our findings illustrate the complex temporal dynamics of impedance in implanted electrodes and the impact of electrical stimulation. We discuss these dynamics in the context of the known biological foreign body response of the brain to implanted electrodes. The data suggest that the temporal dynamics of impedance are dependent on the anatomical location and tissue epileptogenicity. These insights may offer additional guidance for the delivery of therapeutic stimulation at various time points post-implantation for neuromodulation therapy.

A platform for brain network sensing and stimulation with quantitative behavioral tracking: Application to limbic circuit epilepsy.

Temporal lobe epilepsy is a common neurological disease characterized by recurrent seizures. These seizures often originate from limbic networks and people also experience chronic comorbidities related to memory, mood, and sleep (MMS). Deep brain stimulation targeting the anterior nucleus of the thalamus (ANT-DBS) is a proven therapy, but the optimal stimulation parameters remain unclear. We developed a neurotechnology platform for tracking seizures and MMS to enable data streaming between an investigational brain sensing-stimulation implant, mobile devices, and a cloud environment. Artificial Intelligence algorithms provided accurate catalogs of seizures, interictal epileptiform spikes, and wake-sleep brain states. Remotely administered memory and mood assessments were used to densely sample cognitive and behavioral response during ANT-DBS. We evaluated the efficacy of low-frequency versus high-frequency ANT-DBS. They both reduced seizures, but low-frequency ANT-DBS showed greater reductions and better sleep and memory. These results highlight the potential of synchronized brain sensing and behavioral tracking for optimizing neuromodulation therapy.

High and low frequency anterior nucleus of thalamus deep brain stimulation: Impact on memory and mood in five patients with treatment resistant temporal lobe epilepsy.

High frequency anterior nucleus of the thalamus deep brain stimulation (ANT DBS) is an established therapy for treatment resistant focal epilepsies. Although high frequency-ANT DBS is well tolerated, patients are rarely seizure free and the efficacy of other DBS parameters and their impact on comorbidities of epilepsy such as depression and memory dysfunction remain unclear. The purpose of this study was to assess the impact of low vs high frequency ANT DBS on verbal memory and self-reported anxiety and depression symptoms. Five patients with treatment resistant temporal lobe epilepsy were implanted with an investigational brain stimulation and sensing device capable of ANT DBS and ambulatory intracranial electroencephalographic (iEEG) monitoring, enabling long-term detection of electrographic seizures. While patients received therapeutic high frequency (100 and 145 Hz continuous and cycling) and low frequency (2 and 7 Hz continuous) stimulation, they completed weekly free recall verbal memory tasks and thrice weekly self-reports of anxiety and depression symptom severity. Mixed effects models were then used to evaluate associations between memory scores, anxiety and depression self-reports, seizure counts, and stimulation frequency. Memory score was significantly associated with stimulation frequency, with higher free recall verbal memory scores during low frequency ANT DBS. Self-reported anxiety and depression symptom severity was not significantly associated with stimulation frequency. These findings suggest the choice of ANT DBS stimulation parameter may impact patients' cognitive function, independently of its impact on seizure rates.

Flexible multi-step hypothesis testing of human ECoG data using cluster-based permutation tests with GLMEs.

BACKGROUND: Time series analysis is critical for understanding brain signals and their relationship to behavior and cognition. Cluster-based permutation tests (CBPT) are commonly used to analyze a variety of electrophysiological signals including EEG, MEG, ECoG, and sEEG data without a priori assumptions about specific temporal effects. However, two major limitations of CBPT include the inability to directly analyze experiments with multiple fixed effects and the inability to account for random effects (e.g. variability across subjects). Here, we propose a flexible multi-step hypothesis testing strategy using CBPT with Linear Mixed Effects Models (LMEs) and Generalized Linear Mixed Effects Models (GLMEs) that can be applied to a wide range of experimental designs and data types. METHODS: We first evaluate the statistical robustness of LMEs and GLMEs using simulated data distributions. Second, we apply a multi-step hypothesis testing strategy to analyze ERPs and broadband power signals extracted from human ECoG recordings collected during a simple image viewing experiment with image category and novelty as fixed effects. Third, we assess the statistical power differences between analyzing signals with CBPT using LMEs compared to CBPT using separate t-tests run on each fixed effect through simulations that emulate broadband power signals. Finally, we apply CBPT using GLMEs to high-gamma burst data to demonstrate the extension of the proposed method to the analysis of nonlinear data. RESULTS: First, we found that LMEs and GLMEs are robust statistical models. In simple simulations LMEs produced highly congruent results with other appropriately applied linear statistical models, but LMEs outperformed many linear statistical models in the analysis of "suboptimal" data and maintained power better than analyzing individual fixed effects with separate t-tests. GLMEs also performed similarly to other nonlinear statistical models. Second, in real world human ECoG data, LMEs performed at least as well as separate t-tests when applied to predefined time windows or when used in conjunction with CBPT. Additionally, fixed effects time courses extracted with CBPT using LMEs from group-level models of pseudo-populations replicated latency effects found in individual category-selective channels. Third, analysis of simulated broadband power signals demonstrated that CBPT using LMEs was superior to CBPT using separate t-tests in identifying time windows with significant fixed effects especially for small effect sizes. Lastly, the analysis of high-gamma burst data using CBPT with GLMEs produced results consistent with CBPT using LMEs applied to broadband power data. CONCLUSIONS: We propose a general approach for statistical analysis of electrophysiological data using CBPT in conjunction with LMEs and GLMEs. We demonstrate that this method is robust for experiments with multiple fixed effects and applicable to the analysis of linear and nonlinear data. Our methodology maximizes the statistical power available in a dataset across multiple experimental variables while accounting for hierarchical random effects and controlling FWER across fixed effects. This approach substantially improves power leading to better reproducibility. Additionally, CBPT using LMEs and GLMEs can be used to analyze individual channels or pseudo-population data for the comparison of functional or anatomical groups of data.

Acute to long-term characteristics of impedance recordings during neurostimulation in humans.

Objective: This study aims to characterize the time course of impedance, a crucial electrophysiological property of brain tissue, in the human thalamus (THL), amygdala-hippocampus (AMG-HPC), and posterior hippocampus (post-HPC) over an extended period. Approach: Impedance was periodically sampled every 5-15 minutes over several months in five subjects with drug-resistant epilepsy using an experimental neuromodulation device. Initially, we employed descriptive piecewise and continuous mathematical models to characterize the impedance response for approximately three weeks post-electrode implantation. We then explored the temporal dynamics of impedance during periods when electrical stimulation was temporarily halted, observing a monotonic increase (rebound) in impedance before it stabilized at a higher value. Lastly, we assessed the stability of amplitude and phase over the 24-hour impedance cycle throughout the multi-month recording. Main results: Immediately post-implantation, the impedance decreased, reaching a minimum value in all brain regions within approximately two days, and then increased monotonically over about 14 days to a stable value. The models accounted for the variance in short-term impedance changes. Notably, the minimum impedance of the THL in the most epileptogenic hemisphere was significantly lower than in other regions. During the gaps in electrical stimulation, the impedance rebound decreased over time and stabilized around 200 days post-implant, likely indicative of the foreign body response and fibrous tissue encapsulation around the electrodes. The amplitude and phase of the 24-hour impedance oscillation remained stable throughout the multi-month recording, with circadian variation in impedance dominating the long-term measures. Significance: Our findings illustrate the complex temporal dynamics of impedance in implanted electrodes and the impact of electrical stimulation. We discuss these dynamics in the context of the known biological foreign body response of the brain to implanted electrodes. The data suggest that the temporal dynamics of impedance are dependent on the anatomical location and tissue epileptogenicity. These insights may offer additional guidance for the delivery of therapeutic stimulation at various time points post-implantation for neuromodulation therapy.

ER-detect: a pipeline for robust detection of early evoked responses in BIDS-iEEG electrical stimulation data.

Human brain connectivity can be measured in different ways. Intracranial EEG (iEEG) measurements during single pulse electrical stimulation provide a unique way to assess the spread of electrical information with millisecond precision. To provide a robust workflow to process these cortico-cortical evoked potential (CCEP) data and detect early evoked responses in a fully automated and reproducible fashion, we developed Early Response (ER)-detect. ER-detect is an open-source Python package and Docker application to preprocess BIDS structured iEEG data and detect early evoked CCEP responses. ER-detect can use three response detection methods, which were validated against 14-manually annotated CCEP datasets from two different sites by four independent raters. Results showed that ER-detect's automated detection performed on par with the inter-rater reliability (Cohen's Kappa of ~0.6). Moreover, ER-detect was optimized for processing large CCEP datasets, to be used in conjunction with other connectomic investigations. ER-detect provides a highly efficient standardized workflow such that iEEG-BIDS data can be processed in a consistent manner and enhance the reproducibility of CCEP based connectivity results.

Implications of intracranial hemorrhage associated with stereo-EEG.

OBJECTIVE: Intracranial hemorrhage (ICH) is a known complication during stereo-electroencephalography (sEEG) however true rates remain unknown. We provide a comprehensive review of ICH during sEEG regardless of clinical symptoms. Secondly, we analyzed sEEG recordings to identify electrographic correlates of ICH. METHODS: This is a retrospective study of patients undergoing sEEG between January 2016 and April 2022 at the Mayo Clinic in Rochester. We reviewed medical records and imaging studies to identify ICH. We analyzed ICH by type, electrode trajectories, timing, sEEG findings and outcomes. RESULTS: There were a total of 201 sEEG implants, of which 23 (11%) cases or 0.9% electrodes implanted had evidence of ICH. The majority of affected patients (82%) were either asymptomatic or had mild clinical neurological manifestations. In 90% of patients who proceeded with surgical treatments, outcomes were favorable. The most common sEEG finding in contacts in proximity of ICH was either focal slowing with interictal discharges or focal electrographic seizures. CONCLUSIONS: ICH associated with sEEG is likely under-reported in literature. We present electroencephalographic correlates of ICH that may aid identification of ICH in the course of performing sEEG monitoring. SIGNIFICANCE: Our data provides clinically relevant information on potential risks and outcomes of ICH. Furthermore, our findings aid identification of ICH during sEEG.

Impedance Rhythms in Human Limbic System.

The impedance is a fundamental electrical property of brain tissue, playing a crucial role in shaping the characteristics of local field potentials, the extent of ephaptic coupling, and the volume of tissue activated by externally applied electrical brain stimulation. We tracked brain impedance, sleep-wake behavioral state, and epileptiform activity in five people with epilepsy living in their natural environment using an investigational device. The study identified impedance oscillations that span hours to weeks in the amygdala, hippocampus, and anterior nucleus thalamus. The impedance in these limbic brain regions exhibit multiscale cycles with ultradian (∼1.5-1.7 h), circadian (∼21.6-26.4 h), and infradian (∼20-33 d) periods. The ultradian and circadian period cycles are driven by sleep-wake state transitions between wakefulness, nonrapid eye movement (NREM) sleep, and rapid eye movement (REM) sleep. Limbic brain tissue impedance reaches a minimum value in NREM sleep, intermediate values in REM sleep, and rises through the day during wakefulness, reaching a maximum in the early evening before sleep onset. Infradian (∼20-33 d) impedance cycles were not associated with a distinct behavioral correlate. Brain tissue impedance is known to strongly depend on the extracellular space (ECS) volume, and the findings reported here are consistent with sleep-wake-dependent ECS volume changes recently observed in the rodent cortex related to the brain glymphatic system. We hypothesize that human limbic brain ECS changes during sleep-wake state transitions underlie the observed multiscale impedance cycles. Impedance is a simple electrophysiological biomarker that could prove useful for tracking ECS dynamics in human health, disease, and therapy.SIGNIFICANCE STATEMENT The electrical impedance in limbic brain structures (amygdala, hippocampus, anterior nucleus thalamus) is shown to exhibit oscillations over multiple timescales. We observe that impedance oscillations with ultradian and circadian periodicities are associated with transitions between wakefulness, NREM, and REM sleep states. There are also impedance oscillations spanning multiple weeks that do not have a clear behavioral correlate and whose origin remains unclear. These multiscale impedance oscillations will have an impact on extracellular ionic currents that give rise to local field potentials, ephaptic coupling, and the tissue activated by electrical brain stimulation. The approach for measuring tissue impedance using perturbational electrical currents is an established engineering technique that may be useful for tracking ECS volume.

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.

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.