Individualized localization and cortical surface-based registration of intracranial electrodes

A comparison of FreeSurfer-generated data with and without manual intervention

This paper examined whether FreeSurfer-generated data differed between a fully-automated, unedited pipeline and an edited pipeline that included the application of control points to correct errors in white matter segmentation. In a sample of 30 individuals, we compared the summary statistics of surface area, white matter volumes, and cortical thickness derived from edited and unedited datasets for the 34 regions of interest (ROIs) that FreeSurfer (FS) generates. To determine whether applying control points would alter the detection of significant differences between patient and typical groups, effect sizes between edited and unedited conditions in individuals with the genetic disorder, 22q11.2 deletion syndrome (22q11DS) were compared to neurotypical controls. Analyses were conducted with data that were generated from both a 1.5 tesla and a 3 tesla scanner. For 1.5 tesla data, mean area, volume, and thickness measures did not differ significantly between edited and unedited regions, with the exception of rostral anterior cingulate thickness, lateral orbitofrontal white matter, superior parietal white matter, and precentral gyral thickness. Results were similar for surface area and white matter volumes generated from the 3 tesla scanner. For cortical thickness measures however, seven edited ROI measures, primarily in frontal and temporal regions, differed significantly from their unedited counterparts, and three additional ROI measures approached significance. Mean effect sizes for edited ROIs did not differ from most unedited ROIs for either 1.5 or 3 tesla data. Taken together, these results suggest that although the application of control points may increase the validity of intensity normalization and, ultimately, segmentation, it may not affect the final, extracted metrics that FS generates. Potential exceptions to and limitations of these conclusions are discussed.

Combining task-evoked and spontaneous activity to improve pre-operative brain mapping with fMRI

Noninvasive localization of brain function is used to understand and treat neurological disease, exemplified by pre-operative fMRI mapping prior to neurosurgical intervention. The principal approach for generating these maps relies on brain responses evoked by a task and, despite known limitations, has dominated clinical practice for over 20years. Recently, pre-operative fMRI mapping based on correlations in spontaneous brain activity has been demonstrated, however this approach has its own limitations and has not seen widespread clinical use. Here we show that spontaneous and task-based mapping can be performed together using the same pre-operative fMRI data, provide complimentary information relevant for functional localization, and can be combined to improve identification of eloquent motor cortex. Accuracy, sensitivity, and specificity of our approach are quantified through comparison with electrical cortical stimulation mapping in eight patients with intractable epilepsy. Broad applicability and reproducibility of our approach are demonstrated through prospective replication in an independent dataset of six patients from a different center. In both cohorts and every individual patient, we see a significant improvement in signal to noise and mapping accuracy independent of threshold, quantified using receiver operating characteristic curves. Collectively, our results suggest that modifying the processing of fMRI data to incorporate both task-based and spontaneous activity significantly improves functional localization in pre-operative patients. Because this method requires no additional scan time or modification to conventional pre-operative data acquisition protocols it could have widespread utility.

Distribution, Amplitude, Incidence, Co-Occurrence, and Propagation of Human K-Complexes in Focal Transcortical Recordings

K-complexes (KCs) are thought to play a key role in sleep homeostasis and memory consolidation; however, their generation and propagation remain unclear. The commonly held view from scalp EEG findings is that KCs are primarily generated in medial frontal cortex and propagate parietally, whereas an electrocorticography (ECOG) study suggested dorsolateral prefrontal generators and an absence of KCs in many areas. In order to resolve these differing views, we used unambiguously focal bipolar depth electrode recordings in patients with intractable epilepsy to investigate spatiotemporal relationships of human KCs. KCs were marked manually on each channel, and local generation was confirmed with decreased gamma power. In most cases (76%), KCs occurred in a single location, and rarely (1%) in all locations. However, if automatically detected KC-like phenomena were included, only 15% occurred in a single location, and 27% occurred in all recorded locations. Locally generated KCs were found in all sampled areas, including cingulate, ventral temporal, and occipital cortices. Surprisingly, KCs were smallest and occurred least frequently in anterior prefrontal channels. When KCs occur on two channels, their peak order is consistent in only 13% of cases, usually from prefrontal to lateral temporal. Overall, the anterior-posterior separation of electrode pairs explained only 2% of the variance in their latencies. KCs in stages 2 and 3 had similar characteristics. These results open a novel view where KCs overall are universal cortical phenomena, but each KC may variably involve small or large cortical regions and spread in variable directions, allowing flexible and heterogeneous contributions to sleep homeostasis and memory consolidation.

Slow Spatial Recruitment of Neocortex during Secondarily Generalized Seizures and Its Relation to Surgical Outcome

Understanding the spatiotemporal dynamics of brain activity is crucial for inferring the underlying synaptic and nonsynaptic mechanisms of brain dysfunction. Focal seizures with secondary generalization are traditionally considered to begin in a limited spatial region and spread to connected areas, which can include both pathological and normal brain tissue. The mechanisms underlying this spread are important to our understanding of seizures and to improve therapies for surgical intervention. Here we study the properties of seizure recruitment-how electrical brain activity transitions to large voltage fluctuations characteristic of spike-and-wave seizures. We do so using invasive subdural electrode arrays from a population of 16 patients with pharmacoresistant epilepsy. We find an average delay of ∼30 s for a broad area of cortex (8 × 8 cm) to be recruited into the seizure, at an estimated speed of ∼4 mm/s. The spatiotemporal characteristics of recruitment reveal two categories of patients: one in which seizure recruitment of neighboring cortical regions follows a spatially organized pattern consistent from seizure to seizure, and a second group without consistent spatial organization of activity during recruitment. The consistent, organized recruitment correlates with a more regular, compared with small-world, connectivity pattern in simulation and successful surgical treatment of epilepsy. We propose that an improved understanding of how the seizure recruits brain regions into large amplitude voltage fluctuations provides novel information to improve surgical treatment of epilepsy and highlights the slow spread of massive local activity across a vast extent of cortex during seizure.

Dynamic encoding of speech sequence probability in human temporal cortex

Sensory processing involves identification of stimulus features, but also integration with the surrounding sensory and cognitive context. Previous work in animals and humans has shown fine-scale sensitivity to context in the form of learned knowledge about the statistics of the sensory environment, including relative probabilities of discrete units in a stream of sequential auditory input. These statistics are a defining characteristic of one of the most important sequential signals humans encounter: speech. For speech, extensive exposure to a language tunes listeners to the statistics of sound sequences. To address how speech sequence statistics are neurally encoded, we used high-resolution direct cortical recordings from human lateral superior temporal cortex as subjects listened to words and nonwords with varying transition probabilities between sound segments. In addition to their sensitivity to acoustic features (including contextual features, such as coarticulation), we found that neural responses dynamically encoded the language-level probability of both preceding and upcoming speech sounds. Transition probability first negatively modulated neural responses, followed by positive modulation of neural responses, consistent with coordinated predictive and retrospective recognition processes, respectively. Furthermore, transition probability encoding was different for real English words compared with nonwords, providing evidence for online interactions with high-order linguistic knowledge. These results demonstrate that sensory processing of deeply learned stimuli involves integrating physical stimulus features with their contextual sequential structure. Despite not being consciously aware of phoneme sequence statistics, listeners use this information to process spoken input and to link low-level acoustic representations with linguistic information about word identity and meaning.

Development of grouped icEEG for the study of cognitive processing

Invasive intracranial EEG (icEEG) offers a unique opportunity to study human cognitive networks at an unmatched spatiotemporal resolution. To date, the contributions of icEEG have been limited to the individual-level analyses or cohorts whose data are not integrated in any way. Here we discuss how grouped approaches to icEEG overcome challenges related to sparse-sampling, correct for individual variations in response and provide statistically valid models of brain activity in a population. By the generation of whole-brain activity maps, grouped icEEG enables the study of intra and interregional dynamics between distributed cortical substrates exhibiting task-dependent activity. In this fashion, grouped icEEG analyses can provide significant advances in understanding the mechanisms by which cortical networks give rise to cognitive functions.

Improved localization of implanted subdural electrode contacts on magnetic resonance imaging with an elastic image fusion algorithm in an invasive electroencephalography recording

BACKGROUND

Accurate projection of implanted subdural electrode contacts in presurgical evaluation of pharmacoresistant epilepsy cases by invasive electroencephalography is highly relevant. Linear fusion of computed tomography and magnetic resonance images may display the contacts in the wrong position as a result of brain shift effects.

OBJECTIVE

A retrospective study in 5 patients with pharmacoresistant epilepsy was performed to evaluate whether an elastic image fusion algorithm can provide a more accurate projection of the electrode contacts on the preimplantation magnetic resonance images compared with linear fusion.

METHODS

An automated elastic image fusion algorithm (AEF), a guided elastic image fusion algorithm (GEF), and a standard linear fusion algorithm were used on preoperative magnetic resonance images and postimplantation computed tomography scans. Vertical correction of virtual contact positions, total virtual contact shift, corrections of midline shift, and brain shifts caused by pneumocephalus were measured.

RESULTS

Both AEF and GEF worked well with all 5 cases. An average midline shift of 1.7 mm (SD, 1.25 mm) was corrected to 0.4 mm (SD, 0.8 mm) after AEF and to 0.0 mm (SD, 0 mm) after GEF. Median virtual distances between contacts and cortical surface were corrected by a significant amount, from 2.3 mm after linear fusion algorithm to 0.0 mm after AEF and GEF (P < .001). Mean total relative corrections of 3.1 mm (SD, 1.85 mm) after AEF and 3.0 mm (SD, 1.77 mm) after GEF were achieved. The tested version of GEF did not achieve a satisfying virtual correction of pneumocephalus.

CONCLUSION

The technique provided a clear improvement in fusion of preimplantation and postimplantation scans, although the accuracy is difficult to evaluate.

Automatic segmentation of deep intracerebral electrodes in computed tomography scans

BACKGROUND

Invasive monitoring of brain activity by means of intracerebral electrodes is widely practiced to improve pre-surgical seizure onset zone localization in patients with medically refractory seizures. Stereo-Electroencephalography (SEEG) is mainly used to localize the epileptogenic zone and a precise knowledge of the location of the electrodes is expected to facilitate the recordings interpretation and the planning of resective surgery. However, the localization of intracerebral electrodes on post-implant acquisitions is usually time-consuming (i.e., manual segmentation), it requires advanced 3D visualization tools, and it needs the supervision of trained medical doctors in order to minimize the errors. In this paper we propose an automated segmentation algorithm specifically designed to segment SEEG contacts from a thresholded post-implant Cone-Beam CT volume (0.4 mm, 0.4 mm, 0.8 mm). The algorithm relies on the planned position of target and entry points for each electrode as a first estimation of electrode axis. We implemented the proposed algorithm into DEETO, an open source C++ prototype based on ITK library.

RESULTS

We tested our implementation on a cohort of 28 subjects in total. The experimental analysis, carried out over a subset of 12 subjects (35 multilead electrodes; 200 contacts) manually segmented by experts, show that the algorithm: (i) is faster than manual segmentation (i.e., less than 1s/subject versus a few hours) (ii) is reliable, with an error of 0.5 mm ± 0.06 mm, and (iii) it accurately maps SEEG implants to their anatomical regions improving the interpretability of electrophysiological traces for both clinical and research studies. Moreover, using the 28-subject cohort we show here that the algorithm is also robust (error < 0.005 mm) against deep-brain displacements (< 12 mm) of the implanted electrode shaft from those planned before surgery.

CONCLUSIONS

Our method represents, to the best of our knowledge, the first automatic algorithm for the segmentation of SEEG electrodes. The method can be used to accurately identify the neuroanatomical loci of SEEG electrode contacts by a non-expert in a fast and reliable manner.

An open-source automated platform for three-dimensional visualization of subdural electrodes using CT-MRI coregistration

OBJECTIVE

Visualizing implanted subdural electrodes in three-dimensional (3D) space can greatly aid in planning, executing, and validating resection in epilepsy surgery. Coregistration software is available, but cost, complexity, insufficient accuracy, or validation limit adoption. We present a fully automated open-source application, based on a novel method using postimplant computerized tomography (CT) and postimplant magnetic resonance (MR) images, for accurately visualizing intracranial electrodes in 3D space.

METHODS

CT-MR rigid brain coregistration, MR nonrigid registration, and prior-based segmentation were carried out on seven patients. Postimplant CT, postimplant MR, and an external labeled atlas were then aligned in the same space. The coregistration algorithm was validated by manually marking identical anatomic landmarks on the postimplant CT and postimplant MR images. Following coregistration, distances between the center of the landmark masks on the postimplant MR and the coregistered CT images were calculated for all subjects. Algorithms were implemented in open-source software and translated into a "drag and drop" desktop application for Apple Mac OS X.

RESULTS

Despite postoperative brain deformation, the method was able to automatically align intrasubject multimodal images and segment cortical subregions, so that all electrodes could be visualized on the parcellated brain. Manual marking of anatomic landmarks validated the coregistration algorithm with a mean misalignment distance of 2.87 mm (standard deviation 0.58 mm)between the landmarks. Software was easily used by operators without prior image processing experience.

SIGNIFICANCE

We demonstrate an easy to use, novel platform for accurately visualizing subdural electrodes in 3D space on a parcellated brain. We rigorously validated this method using quantitative measures. The method is unique because it involves no preprocessing, is fully automated, and freely available worldwide. A desktop application, as well as the source code, are both available for download on the International Epilepsy Electrophysiology Portal (https://www.ieeg.org) for use and interactive refinement.

Synchronization of isolated downstates (K-complexes) may be caused by cortically-induced disruption of thalamic spindling

Sleep spindles and K-complexes (KCs) define stage 2 NREM sleep (N2) in humans. We recently showed that KCs are isolated downstates characterized by widespread cortical silence. We demonstrate here that KCs can be quasi-synchronous across scalp EEG and across much of the cortex using electrocorticography (ECOG) and localized transcortical recordings (bipolar SEEG). We examine the mechanism of synchronous KC production by creating the first conductance based thalamocortical network model of N2 sleep to generate both spontaneous spindles and KCs. Spontaneous KCs are only observed when the model includes diffuse projections from restricted prefrontal areas to the thalamic reticular nucleus (RE), consistent with recent anatomical findings in rhesus monkeys. Modeled KCs begin with a spontaneous focal depolarization of the prefrontal neurons, followed by depolarization of the RE. Surprisingly, the RE depolarization leads to decreased firing due to disrupted spindling, which in turn is due to depolarization-induced inactivation of the low-threshold Ca²⁺ current (I_T). Further, although the RE inhibits thalamocortical (TC) neurons, decreased RE firing causes decreased TC cell firing, again because of disrupted spindling. The resulting abrupt removal of excitatory input to cortical pyramidal neurons then leads to the downstate. Empirically, KCs may also be evoked by sensory stimuli while maintaining sleep. We reproduce this phenomenon in the model by depolarization of either the RE or the widely-projecting prefrontal neurons. Again, disruption of thalamic spindling plays a key role. Higher levels of RE stimulation also cause downstates, but by directly inhibiting the TC neurons. SEEG recordings from the thalamus and cortex in a single patient demonstrated the model prediction that thalamic spindling significantly decreases before KC onset. In conclusion, we show empirically that KCs can be widespread quasi-synchronous cortical downstates, and demonstrate with the first model of stage 2 NREM sleep a possible mechanism whereby this widespread synchrony may arise.

EEG in the most common stage of human sleep is dominated by K-complexes (KCs) and sleep spindles (bursts of 10–14 Hz oscillations) occupying the thalamus and cortex. Recently, we discovered that KCs are brief moments when the cortex becomes almost completely silent. Here, using recordings directly from the cortex of epileptic patients, we show that KCs can be quasi-synchronous across widespread cortical areas, and ask what mechanism could produce such a phenomenon. We created the first network model of realistic cortical and thalamic neurons, which spontaneously generate KCs as well as sleep spindles. We showed that the membrane channels in the reticular nucleus of the thalamus can be inactivated by excitatory inputs from the cortex, and this disrupts the spindle-generating network, which can trigger widespread cortical silence. The model prediction that thalamic spindle disruption occurs prior to KC was then observed in simultaneous recordings from the human thalamus and cortex. Understanding the cellular and network mechanisms whereby KCs arise is crucial to understanding its roles in maintaining sleep and consolidating memories.

Corticocortical evoked potentials reveal projectors and integrators in human brain networks

The cerebral cortex is composed of subregions whose functional specialization is largely determined by their incoming and outgoing connections with each other. In the present study, we asked which cortical regions can exert the greatest influence over other regions and the cortical network as a whole. Previous research on this question has relied on coarse anatomy (mapping large fiber pathways) or functional connectivity (mapping inter-regional statistical dependencies in ongoing activity). Here we combined direct electrical stimulation with recordings from the cortical surface to provide a novel insight into directed, inter-regional influence within the cerebral cortex of awake humans. These networks of directed interaction were reproducible across strength thresholds and across subjects. Directed network properties included (1) a decrease in the reciprocity of connections with distance; (2) major projector nodes (sources of influence) were found in peri-Rolandic cortex and posterior, basal and polar regions of the temporal lobe; and (3) major receiver nodes (receivers of influence) were found in anterolateral frontal, superior parietal, and superior temporal regions. Connectivity maps derived from electrical stimulation and from resting electrocorticography (ECoG) correlations showed similar spatial distributions for the same source node. However, higher-level network topology analysis revealed differences between electrical stimulation and ECoG that were partially related to the reciprocity of connections. Together, these findings inform our understanding of large-scale corticocortical influence as well as the interpretation of functional connectivity networks.

Localizing ECoG electrodes on the cortical anatomy without post-implantation imaging

INTRODUCTION

Electrocorticographic (ECoG) grids are placed subdurally on the cortex in people undergoing cortical resection to delineate eloquent cortex. ECoG signals have high spatial and temporal resolution and thus can be valuable for neuroscientific research. The value of these data is highest when they can be related to the cortical anatomy. Existing methods that establish this relationship rely either on post-implantation imaging using computed tomography (CT), magnetic resonance imaging (MRI) or X-Rays, or on intra-operative photographs. For research purposes, it is desirable to localize ECoG electrodes on the brain anatomy even when post-operative imaging is not available or when intra-operative photographs do not readily identify anatomical landmarks.

METHODS

We developed a method to co-register ECoG electrodes to the underlying cortical anatomy using only a pre-operative MRI, a clinical neuronavigation device (such as BrainLab VectorVision), and fiducial markers. To validate our technique, we compared our results to data collected from six subjects who also had post-grid implantation imaging available. We compared the electrode coordinates obtained by our fiducial-based method to those obtained using existing methods, which are based on co-registering pre- and post-grid implantation images.

RESULTS

Our fiducial-based method agreed with the MRI-CT method to within an average of 8.24 mm (mean, median = 7.10 mm) across 6 subjects in 3 dimensions. It showed an average discrepancy of 2.7 mm when compared to the results of the intra-operative photograph method in a 2D coordinate system. As this method does not require post-operative imaging such as CTs, our technique should prove useful for research in intra-operative single-stage surgery scenarios. To demonstrate the use of our method, we applied our method during real-time mapping of eloquent cortex during a single-stage surgery. The results demonstrated that our method can be applied intra-operatively in the absence of post-operative imaging to acquire ECoG signals that can be valuable for neuroscientific investigations.

Surface-based mixed effects multilevel analysis of grouped human electrocorticography

Electrocorticography (ECoG) in humans yields data with unmatched spatio-temporal resolution that provides novel insights into cognitive operations. However, the broader application of ECoG has been confounded by difficulties in accurately depicting individual data and performing statistically valid population-level analyses. To overcome these limitations, we developed methods for accurately registering ECoG data to individual cortical topology. We integrated this technique with surface-based co-registration and a mixed-effects multilevel analysis (MEMA) to control for variable cortical surface anatomy and sparse coverage across patients, as well as intra- and inter-subject variability. We applied this surface-based MEMA (SB-MEMA) technique to a face-recognition task dataset (n=22). Compared against existing techniques, SB-MEMA yielded results much more consistent with individual data and with meta-analyses of face-specific activation studies. We anticipate that SB-MEMA will greatly expand the role of ECoG in studies of human cognition, and will enable the generation of population-level brain activity maps and accurate multimodal comparisons.

Evoked effective connectivity of the human neocortex

The role of cortical connectivity in brain function and pathology is increasingly being recognized. While in vivo magnetic resonance imaging studies have provided important insights into anatomical and functional connectivity, these methodologies are limited in their ability to detect electrophysiological activity and the causal relationships that underlie effective connectivity. Here, we describe results of cortico-cortical evoked potential (CCEP) mapping using single pulse electrical stimulation in 25 patients undergoing seizure monitoring with subdural electrode arrays. Mapping was performed by stimulating adjacent electrode pairs and recording CCEPs from the remainder of the electrode array. CCEPs reliably revealed functional networks and showed an inverse relationship to distance between sites. Coregistration to Brodmann areas (BA) permitted group analysis. Connections were frequently directional with 43% of early responses and 50% of late responses of connections reflecting relative dominance of incoming or outgoing connections. The most consistent connections were seen as outgoing from motor cortex, BA6-BA9, somatosensory (SS) cortex, anterior cingulate cortex, and Broca's area. Network topology revealed motor, SS, and premotor cortices along with BA9 and BA10 and language areas to serve as hubs for cortical connections. BA20 and BA39 demonstrated the most consistent dominance of outdegree connections, while BA5, BA7, auditory cortex, and anterior cingulum demonstrated relatively greater indegree. This multicenter, large-scale, directional study of local and long-range cortical connectivity using direct recordings from awake, humans will aid the interpretation of noninvasive functional connectome studies.

A transition in brain state during propofol-induced unconsciousness

Rhythmic oscillations shape cortical dynamics during active behavior, sleep, and general anesthesia. Cross-frequency phase-amplitude coupling is a prominent feature of cortical oscillations, but its role in organizing conscious and unconscious brain states is poorly understood. Using high-density EEG and intracranial electrocorticography during gradual induction of propofol general anesthesia in humans, we discovered a rapid drug-induced transition between distinct states with opposite phase-amplitude coupling and different cortical source distributions. One state occurs during unconsciousness and may be similar to sleep slow oscillations. A second state occurs at the loss or recovery of consciousness and resembles an enhanced slow cortical potential. These results provide objective electrophysiological landmarks of distinct unconscious brain states, and could be used to help improve EEG-based monitoring for general anesthesia.

Rapid and efficient localization of depth electrodes and cortical labeling using free and open source medical software in epilepsy surgery candidates

Depth intracranial electrodes (IEs) placement is one of the most used procedures to identify the epileptogenic zone (EZ) in surgical treatment of drug resistant epilepsy patients, about 20–30% of this population. IEs localization is therefore a critical issue defining the EZ and its relation with eloquent functional areas. That information is then used to target the resective surgery and has great potential to affect outcome. We designed a methodological procedure intended to avoid the need for highly specialized medical resources and reduce time to identify the anatomical location of IEs, during the first instances of intracranial EEG recordings. This workflow is based on established open source software; 3D Slicer and Freesurfer that uses MRI and Post-implant CT fusion for the localization of IEs and its relation with automatic labeled surrounding cortex. To test this hypothesis we assessed the time elapsed between the surgical implantation process and the final anatomical localization of IEs by means of our proposed method compared against traditional visual analysis of raw post-implant imaging in two groups of patients. All IEs were identified in the first 24 H (6–24 H) of implantation using our method in 4 patients of the first group. For the control group; all IEs were identified by experts with an overall time range of 36 h to 3 days using traditional visual analysis. It included (7 patients), 3 patients implanted with IEs and the same 4 patients from the first group. Time to localization was restrained in this group by the specialized personnel and the image quality available. To validate our method; we trained two inexperienced operators to assess the position of IEs contacts on four patients (5 IEs) using the proposed method. We quantified the discrepancies between operators and we also assessed the efficiency of our method to define the EZ comparing the findings against the results of traditional analysis.

Dominant frequencies of resting human brain activity as measured by the electrocorticogram

The brain's spontaneous, intrinsic activity is increasingly being shown to reveal brain function, delineate large scale brain networks, and diagnose brain disorders. One of the most studied and clinically utilized types of intrinsic brain activity are oscillations in the electrocorticogram (ECoG), a relatively localized measure of cortical synaptic activity. Here we objectively characterize the types of ECoG oscillations commonly observed over particular cortical areas when an individual is awake and immobile with eyes closed, using a surface-based cortical atlas and cluster analysis. Both methods show that [1] there is generally substantial variability in the dominant frequencies of cortical regions and substantial overlap in dominant frequencies across the areas sampled (primarily lateral central, temporal, and frontal areas), [2] theta (4-8 Hz) is the most dominant type of oscillation in the areas sampled with a mode around 7 Hz, [3] alpha (8-13 Hz) is largely limited to parietal and occipital regions, and [4] beta (13-30 Hz) is prominent peri-Rolandically, over the middle frontal gyrus, and the pars opercularis. In addition, the cluster analysis revealed seven types of ECoG spectral power densities (SPDs). Six of these have peaks at 3, 5, 7 (narrow), 7 (broad), 10, and 17 Hz, while the remaining cluster is broadly distributed with less pronounced peaks at 8, 19, and 42 Hz. These categories largely corroborate conventional sub-gamma frequency band distinctions (delta, theta, alpha, and beta) and suggest multiple sub-types of theta. Finally, we note that gamma/high gamma activity (30+ Hz) was at times prominently observed, but was too infrequent and variable across individuals to be reliably characterized. These results should help identify abnormal patterns of ECoG oscillations, inform the interpretation of EEG/MEG intrinsic activity, and provide insight into the functions of these different oscillations and the networks that produce them. Specifically, our results support theories of the importance of theta oscillations in general cortical function, suggest that alpha activity is primarily related to sensory processing/attention, and demonstrate that beta networks extend far beyond primary sensorimotor regions.

Local cortical dynamics of burst suppression in the anaesthetized brain

Burst suppression is an electroencephalogram pattern that consists of a quasi-periodic alternation between isoelectric ‘suppressions’ lasting seconds or minutes, and high-voltage ‘bursts’. It is characteristic of a profoundly inactivated brain, occurring in conditions including hypothermia, deep general anaesthesia, infant encephalopathy and coma. It is also used in neurology as an electrophysiological endpoint in pharmacologically induced coma for brain protection after traumatic injury and during status epilepticus. Classically, burst suppression has been regarded as a ‘global’ state with synchronous activity throughout cortex. This assumption has influenced the clinical use of burst suppression as a way to broadly reduce neural activity. However, the extent of spatial homogeneity has not been fully explored due to the challenges in recording from multiple cortical sites simultaneously. The neurophysiological dynamics of large-scale cortical circuits during burst suppression are therefore not well understood. To address this question, we recorded intracranial electrocorticograms from patients who entered burst suppression while receiving propofol general anaesthesia. The electrodes were broadly distributed across cortex, enabling us to examine both the dynamics of burst suppression within local cortical regions and larger-scale network interactions. We found that in contrast to previous characterizations, bursts could be substantially asynchronous across the cortex. Furthermore, the state of burst suppression itself could occur in a limited cortical region while other areas exhibited ongoing continuous activity. In addition, we found a complex temporal structure within bursts, which recapitulated the spectral dynamics of the state preceding burst suppression, and evolved throughout the course of a single burst. Our observations imply that local cortical dynamics are not homogeneous, even during significant brain inactivation. Instead, cortical and, implicitly, subcortical circuits express seemingly different sensitivities to high doses of anaesthetics that suggest a hierarchy governing how the brain enters burst suppression, and emphasize the role of local dynamics in what has previously been regarded as a global state. These findings suggest a conceptual shift in how neurologists could assess the brain function of patients undergoing burst suppression. First, analysing spatial variation in burst suppression could provide insight into the circuit dysfunction underlying a given pathology, and could improve monitoring of medically-induced coma. Second, analysing the temporal dynamics within a burst could help assess the underlying brain state. This approach could be explored as a prognostic tool for recovery from coma, and for guiding treatment of status epilepticus. Overall, these results suggest new research directions and methods that could improve patient monitoring in clinical practice.

Electrode localization for planning surgical resection of the epileptogenic zone in pediatric epilepsy

PURPOSE

In planning for a potentially curative resection of the epileptogenic zone in patients with pediatric epilepsy, invasive monitoring with intracranial EEG is often used to localize the seizure onset zone and eloquent cortex. A precise understanding of the location of subdural strip and grid electrodes on the brain surface, and of depth electrodes in the brain in relationship to eloquent areas is expected to facilitate pre-surgical planning.

METHODS

We developed a novel algorithm for the alignment of intracranial electrodes, extracted from post-operative CT, with pre-operative MRI. Our goal was to develop a method of achieving highly accurate localization of subdural and depth electrodes, in order to facilitate surgical planning. Specifically, we created a patient-specific 3D geometric model of the cortical surface from automatic segmentation of a pre-operative MRI, automatically segmented electrodes from post-operative CT, and projected each set of electrodes onto the brain surface after alignment of the CT to the MRI. Also, we produced critical visualization of anatomical landmarks, e.g., vasculature, gyri, sulci, lesions, or eloquent cortical areas, which enables the epilepsy surgery team to accurately estimate the distance between the electrodes and the anatomical landmarks, which might help for better assessment of risks and benefits of surgical resection.

RESULTS

Electrode localization accuracy was measured using knowledge of the position of placement from 2D intra-operative photographs in ten consecutive subjects who underwent intracranial EEG for pediatric epilepsy. Average spatial accuracy of localization was 1.31 ± 0.69 mm for all 385 visible electrodes in the photos.

CONCLUSIONS

In comparison with previously reported approaches, our algorithm is able to achieve more accurate alignment of strip and grid electrodes with minimal user input. Unlike manual alignment procedures, our algorithm achieves excellent alignment without time-consuming and difficult judgements from an operator.

Neurophysiological investigation of spontaneous correlated and anticorrelated fluctuations of the BOLD signal

Analyses of intrinsic fMRI BOLD signal fluctuations reliably reveal correlated and anticorrelated functional networks in the brain. Because the BOLD signal is an indirect measure of neuronal activity and anticorrelations can be introduced by preprocessing steps, such as global signal regression, the neurophysiological significance of correlated and anticorrelated BOLD fluctuations is a source of debate. Here, we address this question by examining the correspondence between the spatial organization of correlated BOLD fluctuations and correlated fluctuations in electrophysiological high γ power signals recorded directly from the cortical surface of 5 patients. We demonstrate that both positive and negative BOLD correlations have neurophysiological correlates reflected in fluctuations of spontaneous neuronal activity. Although applying global signal regression to BOLD signals results in some BOLD anticorrelations that are not apparent in the ECoG data, it enhances the neuronal-hemodynamic correspondence overall. Together, these findings provide support for the neurophysiological fidelity of BOLD correlations and anticorrelations.

Speech-specific tuning of neurons in human superior temporal gyrus

How the brain extracts words from auditory signals is an unanswered question. We recorded approximately 150 single and multi-units from the left anterior superior temporal gyrus of a patient during multiple auditory experiments. Against low background activity, 45% of units robustly fired to particular spoken words with little or no response to pure tones, noise-vocoded speech, or environmental sounds. Many units were tuned to complex but specific sets of phonemes, which were influenced by local context but invariant to speaker, and suppressed during self-produced speech. The firing of several units to specific visual letters was correlated with their response to the corresponding auditory phonemes, providing the first direct neural evidence for phonological recoding during reading. Maximal decoding of individual phonemes and words identities was attained using firing rates from approximately 5 neurons within 200 ms after word onset. Thus, neurons in human superior temporal gyrus use sparse spatially organized population encoding of complex acoustic-phonetic features to help recognize auditory and visual words.

Recursive grid partitioning on a cortical surface model: an optimized technique for the localization of implanted subdural electrodes

Object

Precise localization of subdural electrodes (SDEs) is essential for the interpretation of data from intracranial electrocorticography recordings. Blood and fluid accumulation underneath the craniotomy flap leads to a nonlinear deformation of the brain surface and of the SDE array on postoperative CT scans and adversely impacts the accurate localization of electrodes located underneath the craniotomy. Older methods that localize electrodes based on their identification on a postimplantation CT scan with coregistration to a preimplantation MR image can result in significant problems with accuracy of the electrode localization. The authors report 3 novel methods that rely on the creation of a set of 3D mesh models to depict the pial surface and a smoothed pial envelope. Two of these new methods are designed to localize electrodes, and they are compared with 6 methods currently in use to determine their relative accuracy and reliability.

Methods

The first method involves manually localizing each electrode using digital photographs obtained at surgery. This is highly accurate, but requires time intensive, operator-dependent input. The second uses 4 electrodes localized manually in conjunction with an automated, recursive partitioning technique to localize the entire electrode array. The authors evaluated the accuracy of previously published methods by applying the methods to their data and comparing them against the photograph-based localization. Finally, the authors further enhanced the usability of these methods by using automatic parcellation techniques to assign anatomical labels to individual electrodes as well as by generating an inflated cortical surface model while still preserving electrode locations relative to the cortical anatomy.

Results

The recursive grid partitioning had the least error compared with older methods (672 electrodes, 6.4-mm maximum electrode error, 2.0-mm mean error, p < 10−18). The maximum errors derived using prior methods of localization ranged from 8.2 to 11.7 mm for an individual electrode, with mean errors ranging between 2.9 and 4.1 mm depending on the method used. The authors also noted a larger error in all methods that used CT scans alone to localize electrodes compared with those that used both postoperative CT and postoperative MRI. The large mean errors reported with these methods are liable to affect intermodal data comparisons (for example, with functional mapping techniques) and may impact surgical decision making.

Conclusions

The authors have presented several aspects of using new techniques to visualize electrodes implanted for localizing epilepsy. The ability to use automated labeling schemas to denote which gyrus a particular electrode overlies is potentially of great utility in planning resections and in corroborating the results of extraoperative stimulation mapping. Dilation of the pial mesh model provides, for the first time, a sense of the cortical surface not sampled by the electrode, and the potential roles this “electrophysiologically hidden” cortex may play in both eloquent function and seizure onset.

Rapid fragmentation of neuronal networks at the onset of propofol-induced unconsciousness

The neurophysiological mechanisms by which anesthetic drugs cause loss of consciousness are poorly understood. Anesthetic actions at the molecular, cellular, and systems levels have been studied in detail at steady states of deep general anesthesia. However, little is known about how anesthetics alter neural activity during the transition into unconsciousness. We recorded simultaneous multiscale neural activity from human cortex, including ensembles of single neurons, local field potentials, and intracranial electrocorticograms, during induction of general anesthesia. We analyzed local and global neuronal network changes that occurred simultaneously with loss of consciousness. We show that propofol-induced unconsciousness occurs within seconds of the abrupt onset of a slow (<1 Hz) oscillation in the local field potential. This oscillation marks a state in which cortical neurons maintain local patterns of network activity, but this activity is fragmented across both time and space. Local (<4 mm) neuronal populations maintain the millisecond-scale connectivity patterns observed in the awake state, and spike rates fluctuate and can reach baseline levels. However, neuronal spiking occurs only within a limited slow oscillation-phase window and is silent otherwise, fragmenting the time course of neural activity. Unexpectedly, we found that these slow oscillations occur asynchronously across cortex, disrupting functional connectivity between cortical areas. We conclude that the onset of slow oscillations is a neural correlate of propofol-induced loss of consciousness, marking a shift to cortical dynamics in which local neuronal networks remain intact but become functionally isolated in time and space.

Localization of dense intracranial electrode arrays using magnetic resonance imaging

Intracranial electrode arrays are routinely used in the pre-surgical evaluation of patients with medically refractory epilepsy, and recordings from these electrodes have been increasingly employed in human cognitive neurophysiology due to their high spatial and temporal resolution. For both researchers and clinicians, it is critical to localize electrode positions relative to the subject-specific neuroanatomy. In many centers, a post-implantation MRI is utilized for electrode detection because of its higher sensitivity for surgical complications and the absence of radiation. However, magnetic susceptibility artifacts surrounding each electrode prohibit unambiguous detection of individual electrodes, especially those that are embedded within dense grid arrays. Here, we present an efficient method to accurately localize intracranial electrode arrays based on pre- and post-implantation MR images that incorporates array geometry and the individual's cortical surface. Electrodes are directly visualized relative to the underlying gyral anatomy of the reconstructed cortical surface of individual patients. Validation of this approach shows high spatial accuracy of the localized electrode positions (mean of 0.96 mm ± 0.81 mm for 271 electrodes across 8 patients). Minimal user input, short processing time, and utilization of radiation-free imaging are strong incentives to incorporate quantitatively accurate localization of intracranial electrode arrays with MRI for research and clinical purposes. Co-registration to a standard brain atlas further allows inter-subject comparisons and relation of intracranial EEG findings to the larger body of neuroimaging literature.

Independence of early speech processing from word meaning

We combined magnetoencephalography (MEG) with magnetic resonance imaging and electrocorticography to separate in anatomy and latency 2 fundamental stages underlying speech comprehension. The first acoustic-phonetic stage is selective for words relative to control stimuli individually matched on acoustic properties. It begins ∼60 ms after stimulus onset and is localized to middle superior temporal cortex. It was replicated in another experiment, but is strongly dissociated from the response to tones in the same subjects. Within the same task, semantic priming of the same words by a related picture modulates cortical processing in a broader network, but this does not begin until ∼217 ms. The earlier onset of acoustic-phonetic processing compared with lexico-semantic modulation was significant in each individual subject. The MEG source estimates were confirmed with intracranial local field potential and high gamma power responses acquired in 2 additional subjects performing the same task. These recordings further identified sites within superior temporal cortex that responded only to the acoustic-phonetic contrast at short latencies, or the lexico-semantic at long. The independence of the early acoustic-phonetic response from semantic context suggests a limited role for lexical feedback in early speech perception.