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.

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.

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.

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.

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.

The Past, Present, and Future of the Brain Imaging Data Structure (BIDS).

The Brain Imaging Data Structure (BIDS) is a community-driven standard for the organization of data and metadata from a growing range of neuroscience modalities. This paper is meant as a history of how the standard has developed and grown over time. We outline the principles behind the project, the mechanisms by which it has been extended, and some of the challenges being addressed as it evolves. We also discuss the lessons learned through the project, with the aim of enabling researchers in other domains to learn from the success of BIDS.

Dynamic Visualization of Gyral and Sulcal Stereoelectroencephalographic Contacts in Humans.

Stereoelectroencephalography (SEEG) is a neurosurgical method to survey electrophysiological activity within the brain to treat disorders such as Epilepsy. In this stereotactic approach, leads are implanted through straight trajectories to survey both cortical and sub-cortical activity.Visualizing the recorded locations covering sulcal and gyral activity while staying true to the cortical architecture is challenging due to the folded, three-dimensional nature of the human cortex.To overcome this challenge, we developed a novel visualization concept, allowing investigators to dynamically morph between the subjects' cortical reconstruction and an inflated cortex representation. This inflated view, in which gyri and sulci are viewed on a smooth surface, allows better visualization of electrodes buried within the sulcus while staying true to the underlying cortical architecture.Clinical relevance- These visualization techniques might also help guide clinical decision-making when defining seizure onset zones or resections for patients undergoing SEEG monitoring for intractable epilepsy.

A comparison of neuroelectrophysiology databases.

As data sharing has become more prevalent, three pillars - archives, standards, and analysis tools - have emerged as critical components in facilitating effective data sharing and collaboration. This paper compares four freely available intracranial neuroelectrophysiology data repositories: Data Archive for the BRAIN Initiative (DABI), Distributed Archives for Neurophysiology Data Integration (DANDI), OpenNeuro, and Brain-CODE. The aim of this review is to describe archives that provide researchers with tools to store, share, and reanalyze both human and non-human neurophysiology data based on criteria that are of interest to the neuroscientific community. The Brain Imaging Data Structure (BIDS) and Neurodata Without Borders (NWB) are utilized by these archives to make data more accessible to researchers by implementing a common standard. As the necessity for integrating large-scale analysis into data repository platforms continues to grow within the neuroscientific community, this article will highlight the various analytical and customizable tools developed within the chosen archives that may advance the field of neuroinformatics.

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.

A Comparison of Neuroelectrophysiology Databases.

As data sharing has become more prevalent, three pillars - archives, standards, and analysis tools - have emerged as critical components in facilitating effective data sharing and collaboration. This paper compares four freely available intracranial neuroelectrophysiology data repositories: Data Archive for the BRAIN Initiative (DABI), Distributed Archives for Neurophysiology Data Integration (DANDI), OpenNeuro, and Brain-CODE. The aim of this review is to describe archives that provide researchers with tools to store, share, and reanalyze both human and non-human neurophysiology data based on criteria that are of interest to the neuroscientific community. The Brain Imaging Data Structure (BIDS) and Neurodata Without Borders (NWB) are utilized by these archives to make data more accessible to researchers by implementing a common standard. As the necessity for integrating large-scale analysis into data repository platforms continues to grow within the neuroscientific community, this article will highlight the various analytical and customizable tools developed within the chosen archives that may advance the field of neuroinformatics.

Electrical Stimulation of Temporal and Limbic Circuitry Produces Distinct Responses in Human Ventral Temporal Cortex.

The human ventral temporal cortex (VTC) is highly connected to integrate visual perceptual inputs with feedback from cognitive and emotional networks. In this study, we used electrical brain stimulation to understand how different inputs from multiple brain regions drive unique electrophysiological responses in the VTC. We recorded intracranial EEG data in 5 patients (3 female) implanted with intracranial electrodes for epilepsy surgery evaluation. Pairs of electrodes were stimulated with single-pulse electrical stimulation, and corticocortical evoked potential responses were measured at electrodes in the collateral sulcus and lateral occipitotemporal sulcus of the VTC. Using a novel unsupervised machine learning method, we uncovered 2-4 distinct response shapes, termed basis profile curves (BPCs), at each measurement electrode in the 11-500 ms after stimulation interval. Corticocortical evoked potentials of unique shape and high amplitude were elicited following stimulation of several regions and classified into a set of four consensus BPCs across subjects. One of the consensus BPCs was primarily elicited by stimulation of the hippocampus; another by stimulation of the amygdala; a third by stimulation of lateral cortical sites, such as the middle temporal gyrus; and the final one by stimulation of multiple distributed sites. Stimulation also produced sustained high-frequency power decreases and low-frequency power increases that spanned multiple BPC categories. Characterizing distinct shapes in stimulation responses provides a novel description of connectivity to the VTC and reveals significant differences in input from cortical and limbic structures.SIGNIFICANCE STATEMENT Disentangling the numerous input influences on highly connected areas in the brain is a critical step toward understanding how brain networks work together to coordinate human behavior. Single-pulse electrical stimulation is an effective tool to accomplish this goal because the shapes and amplitudes of signals recorded from electrodes are informative of the synaptic physiology of the stimulation-driven inputs. We focused on targets in the ventral temporal cortex, an area strongly implicated in visual object perception. By using a data-driven clustering algorithm, we identified anatomic regions with distinct input connectivity profiles to the ventral temporal cortex. Examining high-frequency power changes revealed possible modulation of excitability at the recording site induced by electrical stimulation of connected regions.

Canonical Response Parameterization: Quantifying the structure of responses to single-pulse intracranial electrical brain stimulation.

Single-pulse electrical stimulation in the nervous system, often called cortico-cortical evoked potential (CCEP) measurement, is an important technique to understand how brain regions interact with one another. Voltages are measured from implanted electrodes in one brain area while stimulating another with brief current impulses separated by several seconds. Historically, researchers have tried to understand the significance of evoked voltage polyphasic deflections by visual inspection, but no general-purpose tool has emerged to understand their shapes or describe them mathematically. We describe and illustrate a new technique to parameterize brain stimulation data, where voltage response traces are projected into one another using a semi-normalized dot product. The length of timepoints from stimulation included in the dot product is varied to obtain a temporal profile of structural significance, and the peak of the profile uniquely identifies the duration of the response. Using linear kernel PCA, a canonical response shape is obtained over this duration, and then single-trial traces are parameterized as a projection of this canonical shape with a residual term. Such parameterization allows for dissimilar trace shapes from different brain areas to be directly compared by quantifying cross-projection magnitudes, response duration, canonical shape projection amplitudes, signal-to-noise ratios, explained variance, and statistical significance. Artifactual trials are automatically identified by outliers in sub-distributions of cross-projection magnitude, and rejected. This technique, which we call "Canonical Response Parameterization" (CRP) dramatically simplifies the study of CCEP shapes, and may also be applied in a wide range of other settings involving event-triggered data.

A motor association area in the depths of the central sulcus.

Cells in the precentral gyrus directly send signals to the periphery to generate movement and are principally organized as a topological map of the body. We find that movement-induced electrophysiological responses from depth electrodes extend this map three-dimensionally throughout the gyrus. Unexpectedly, this organization is interrupted by a previously undescribed motor association area in the depths of the midlateral aspect of the central sulcus. This 'Rolandic motor association' (RMA) area is active during movements of different body parts from both sides of the body and may be important for coordinating complex behaviors.

Developmental trajectory of transmission speed in the human brain.

The structure of the human connectome develops from childhood throughout adolescence to middle age, but how these structural changes affect the speed of neuronal signaling is not well described. In 74 subjects, we measured the latency of cortico-cortical evoked responses across association and U-fibers and calculated their corresponding transmission speeds. Decreases in conduction delays until at least 30 years show that the speed of neuronal communication develops well into adulthood.

Dynamic Visualization of Gyral and Sulcal Stereoelectroencephalographic contacts in Humans.

Stereoelectroencephalography (SEEG) is a neurosurgical method to survey electrophysiological activity within the brain to treat disorders such as Epilepsy. In this stereotactic approach, leads are implanted through straight trajectories to survey both cortical and sub-cortical activity. Visualizing the recorded locations covering sulcal and gyral activity while staying true to the cortical architecture is challenging due to the folded, three-dimensional nature of the human cortex. To overcome this challenge, we developed a novel visualization concept, allowing investigators to dynamically morph between the subjects' cortical reconstruction and an inflated cortex representation. This inflated view, in which gyri and sulci are viewed on a smooth surface, allows better visualization of electrodes buried within the sulcus while staying true to the underlying cortical architecture. Clinical relevance­: These visualization techniques might also help guide clinical decision-making when defining seizure onset zones or resections for patients undergoing SEEG monitoring for intractable epilepsy.

Measuring brain beats: Cardiac-aligned fast functional magnetic resonance imaging signals.

Blood and cerebrospinal fluid (CSF) pulse and flow throughout the brain, driven by the cardiac cycle. These fluid dynamics, which are essential to healthy brain function, are characterized by several noninvasive magnetic resonance imaging (MRI) methods. Recent developments in fast MRI, specifically simultaneous multislice acquisition methods, provide a new opportunity to rapidly and broadly assess cardiac-driven flow, including CSF spaces, surface vessels and parenchymal vessels. We use these techniques to assess blood and CSF flow dynamics in brief (3.5 min) scans on a conventional 3 T MRI scanner in five subjects. Cardiac pulses are measured with a photoplethysmography (PPG) on the index finger, along with functional MRI (fMRI) signals in the brain. We, retrospectively, align the fMRI signals to the heartbeat. Highly reliable cardiac-gated fMRI temporal signals are observed in CSF and blood on the timescale of one heartbeat (test-retest reliability within subjects R2  > 50%). In blood vessels, a local minimum is observed following systole. In CSF spaces, the ventricles and subarachnoid spaces have a local maximum following systole instead. Slower resting-state scans with slice timing, retrospectively, aligned to the cardiac pulse, reveal similar cardiac-gated responses. The cardiac-gated measurements estimate the amplitude and phase of fMRI pulsations in the CSF relative to those in the arteries, an estimate of the local intracranial impedance. Cardiac aligned fMRI signals can provide new insights about fluid dynamics or diagnostics for diseases where these dynamics are important.

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.

A Practical Workflow for Organizing Clinical Intraoperative and Long-term iEEG Data in BIDS.

The neuroscience community increasingly uses the Brain Imaging Data Structure (BIDS) to organize data, extending from MRI to electrophysiology data. While automated tools and workflows are developed that help organize MRI data from the scanner to BIDS, these workflows are lacking for clinical intracranial EEG (iEEG data). We present a practical workflow on how to organize full clinical iEEG epilepsy data into BIDS. We present electrophysiological datasets recorded from twelve subjects who underwent intracranial monitoring followed by resective epilepsy surgery at the University Medical Center Utrecht, the Netherlands, and became seizure-free after surgery. These data include intraoperative electrocorticography recordings from six patients, long-term electrocorticography recordings from three patients and stereo-encephalography recordings from three patients. We describe the 6 steps in the pipeline that are essential to structure the data from these clinical iEEG recordings into BIDS and the challenges during this process. These proposed workflow enable centers performing clinical iEEG recordings to structure their data to improve accessibility, reusability and interoperability of clinical data.