Physical model of coherent potentials measured with different electrode recording site sizes

Comparison of EEG Signal Characteristics of Subdural and Depth Electrodes

OBJECTIVES

Our study aimed to compare signal characteristics of subdural electrodes (SDE) and depth stereo EEG placed within a 5-mm vicinity in patients with drug-resistant epilepsy. We report how electrode design and placement collectively affect signal content from a shared source between these electrode types.

METHODS

In subjects undergoing invasive intracranial EEG evaluation at a surgical epilepsy center from 2012 to 2018, stereo EEG and SDE electrode contacts placed within a 5-mm vicinity were identified. Of these, 24 contacts (12 pairs) met our criteria for signal-to-noise ratio and data availability for final analysis. We used Welch method to analyze the correlation of power spectral densities of EEG segments, root mean square of 1-second windows, and fast-Fourier transform to calculate coherence across conventional frequency bands.

RESULTS

We observed a median distance of 3.7 mm between the electrode contact pairs. Time-aware analysis highlighted the coherence's strength primarily in the high-gamma band, where the median (r) was 0.889. In addition, the median power ratios between the SDE and stereo EEG signal was 1.99. This ratio decreased from high-gamma to infra-low frequencies, with medians of 2.07 and 0.97, respectively. The power spectral densities for the stereo EEG and SDE electrodes demonstrated a strong correlation, with a median correlation coefficient (r) of 0.99 and an interquartile range from 0.915 to 0.996.

CONCLUSIONS

Signals captured by standard subdural and depth (intracranial EEG) electrodes within a 5-mm radius exhibit band-specific coherence and are not identical. The association was most pronounced in the high-gamma band, with coherence decreasing with lower frequencies. Our findings underscore the combined effects of electrode size, design, placement, preferred bandwidth, and the nature of the activity source on signal recording. Particularly, SDE employed herein may offer advantages for high-frequency signals, but the impact of electrode size on recordings necessitates careful consideration in context-specific situations.

SIGNIFICANCE

The findings relate to surgical epilepsy care and may inform the design of brain-computer interface.

Electrode Surface Area Impacts Measurement of High Frequency Oscillations in Human Intracranial EEG

OBJECTIVE

High-frequency oscillations (HFOs) are a promising prognostic biomarker of surgical outcome in patients with epilepsy. Their rates of occurrence and morphology have been studied extensively using recordings from electrodes of various geometries. While electrode size is a potential confounding factor in HFO studies, it has largely been disregarded due to a lack of consistent evidence. Therefore, we designed an experiment to directly test the impact of electrode size on HFO measurement.

METHODS

We first simulated HFO measurement using a lumped model of the electrode-tissue interaction. Then eight human subjects were each implanted with a high-density 8x8 grid of subdural electrodes. After implantation, the electrode sizes were altered using a technique recently developed by our group, enabling intracranial EEG recordings for three different electrode surface areas from a static brain location. HFOs were automatically detected in the data and their characteristics were calculated.

RESULTS

The human subject measurements were consistent with the model. Specifically, HFO rate measured per area of tissue decreased significantly as electrode surface area increased. The smallest electrodes recorded more fast ripples than ripples. Amplitude of detected HFOs also decreased as electrode surface area increased, while duration and peak frequency were unaffected.

CONCLUSION

These results suggest that HFO rates measured using electrodes of different surface areas cannot be compared directly.

SIGNIFICANCE

This has significant implications for HFOs as a tool for surgical planning, particularly for individual patients implanted with electrodes of multiple sizes and comparisons of HFO rate made across patients and studies.

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

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

Flexible, high-resolution cortical arrays with large coverage capture microscale high-frequency oscillations in patients with epilepsy

OBJECTIVE

Effective surgical treatment of drug-resistant epilepsy depends on accurate localization of the epileptogenic zone (EZ). High-frequency oscillations (HFOs) are potential biomarkers of the EZ. Previous research has shown that HFOs often occur within submillimeter areas of brain tissue and that the coarse spatial sampling of clinical intracranial electrode arrays may limit the accurate capture of HFO activity. In this study, we sought to characterize microscale HFO activity captured on thin, flexible microelectrocorticographic (μECoG) arrays, which provide high spatial resolution over large cortical surface areas.

METHODS

We used novel liquid crystal polymer thin-film μECoG arrays (.76-1.72-mm intercontact spacing) to capture HFOs in eight intraoperative recordings from seven patients with epilepsy. We identified ripple (80-250 Hz) and fast ripple (250-600 Hz) HFOs using a common energy thresholding detection algorithm along with two stages of artifact rejection. We visualized microscale subregions of HFO activity using spatial maps of HFO rate, signal-to-noise ratio, and mean peak frequency. We quantified the spatial extent of HFO events by measuring covariance between detected HFOs and surrounding activity. We also compared HFO detection rates on microcontacts to simulated macrocontacts by spatially averaging data.

RESULTS

We found visually delineable subregions of elevated HFO activity within each μECoG recording. Forty-seven percent of HFOs occurred on single 200-μm-diameter recording contacts, with minimal high-frequency activity on surrounding contacts. Other HFO events occurred across multiple contacts simultaneously, with covarying activity most often limited to a .95-mm radius. Through spatial averaging, we estimated that macrocontacts with 2-3-mm diameter would only capture 44% of the HFOs detected in our μECoG recordings.

SIGNIFICANCE

These results demonstrate that thin-film microcontact surface arrays with both highresolution and large coverage accurately capture microscale HFO activity and may improve the utility of HFOs to localize the EZ for treatment of drug-resistant epilepsy.

A novel method for dynamically altering the surface area of intracranial EEG electrodes

Objective.Intracranial electroencephalogram (iEEG) plays a critical role in the treatment of neurological diseases, such as epilepsy and Parkinson's disease, as well as the development of neural prostheses and brain computer interfaces. While electrode geometries vary widely across these applications, the impact of electrode size on iEEG features and morphology is not well understood. Some insight has been gained from computer simulations, as well as experiments in which signals are recorded using electrodes of different sizes concurrently in different brain regions. Here, we introduce a novel method to record from electrodes of different sizes in the exact same location by changing the size of iEEG electrodes after implantation in the brain.Approach.We first present a theoretical model and anin vitrovalidation of the method. We then report the results of anin vivoimplementation in three human subjects with refractory epilepsy. We recorded iEEG data from three different electrode sizes and compared the amplitudes, power spectra, inter-channel correlations, and signal-to-noise ratio (SNR) of interictal epileptiform discharges, i.e. epileptic spikes.Main Results.We found that iEEG amplitude and power decreased as electrode size increased, while inter-channel correlation did not change significantly with electrode size. The SNR of epileptic spikes was generally highest in the smallest electrodes, but 39% of spikes had maximal SNR in larger electrodes. This likely depends on the precise location and spatial spread of each spike.Significance.Overall, this new method enables multi-scale measurements of electrical activity in the human brain that can facilitate our understanding of neurophysiology, treatment of neurological disease, and development of novel technologies.

Brain rhythms define distinct interaction networks with differential dependence on anatomy

Cognitive functions are subserved by rhythmic neuronal synchronization across widely distributed brain areas. In 105 area pairs, we investigated functional connectivity (FC) through coherence, power correlation, and Granger causality (GC) in the theta, beta, high-beta, and gamma rhythms. Between rhythms, spatial FC patterns were largely independent. Thus, the rhythms defined distinct interaction networks. Importantly, networks of coherence and GC were not explained by the spatial distributions of the strengths of the rhythms. Those networks, particularly the GC networks, contained clear modules, with typically one dominant rhythm per module. To understand how this distinctiveness and modularity arises on a common anatomical backbone, we correlated, across 91 area pairs, the metrics of functional interaction with those of anatomical projection strength. Anatomy was primarily related to coherence and GC, with the largest effect sizes for GC. The correlation differed markedly between rhythms, being less pronounced for the beta and strongest for the gamma rhythm.

Cellular Classes in the Human Brain Revealed In Vivo by Heartbeat-Related Modulation of the Extracellular Action Potential Waveform

Determining cell types is critical for understanding neural circuits but remains elusive in the living human brain. Current approaches discriminate units into putative cell classes using features of the extracellular action potential (EAP); in absence of ground truth data, this remains a problematic procedure. We find that EAPs in deep structures of the brain exhibit robust and systematic variability during the cardiac cycle. These cardiac-related features refine neural classification. We use these features to link bio-realistic models generated from in vitro human whole-cell recordings of morphologically classified neurons to in vivo recordings. We differentiate aspiny inhibitory and spiny excitatory human hippocampal neurons and, in a second stage, demonstrate that cardiac-motion features reveal two types of spiny neurons with distinct intrinsic electrophysiological properties and phase-locking characteristics to endogenous oscillations. This multi-modal approach markedly improves cell classification in humans, offers interpretable cell classes, and is applicable to other brain areas and species.

Adolescents Exhibit Late Maturation of Long-Range Beta Coherences in Affective Processing

We investigated whether intra-/interhemispheric long-range beta coherences mirror developmental changes in affective functional integration during adolescence. Electroencephalogram data were gathered from 15 young adolescents, 16 old adolescents, and 16 young adults during viewing affective pictures. The results indicated that both positive and negative pictures induced greater intra- and interhemispheric long-distance beta coherences than neutral pictures. However, opposite results were observed between young and old adolescents in terms of negative phase differences. Old adolescents exhibited greater beta coherences for positive and negative pictures than both young adolescents and young adults, but there was no difference between the groups for neutral pictures. These observations suggest that long-range beta coherence might reflect the late maturation of affective functional integration in adolescents.

Differences in MEG and EEG power-law scaling explained by a coupling between spatial coherence and frequency: a simulation study

Electrophysiological signals (electroencephalography, EEG, and magnetoencephalography, MEG), as many natural processes, exhibit scale-invariance properties resulting in a power-law (1/f) spectrum. Interestingly, EEG and MEG differ in their slopes, which could be explained by several mechanisms, including non-resistive properties of tissues. Our goal in the present study is to estimate the impact of space/frequency structure of source signals as a putative mechanism to explain spectral scaling properties of neuroimaging signals. We performed simulations based on the summed contribution of cortical patches with different sizes (ranging from 0.4 to 104.2 cm²). Small patches were attributed signals of high frequencies, whereas large patches were associated with signals of low frequencies, on a logarithmic scale. The tested parameters included i) the space/frequency structure (range of patch sizes and frequencies) and ii) the amplitude factor c parametrizing the spatial scale ratios. We found that the space/frequency structure may cause differences between EEG and MEG scale-free spectra that are compatible with real data findings reported in previous studies. We also found that below a certain spatial scale, there were no more differences between EEG and MEG, suggesting a limit for the resolution of both methods.Our work provides an explanation of experimental findings. This does not rule out other mechanisms for differences between EEG and MEG, but suggests an important role of spatio-temporal structure of neural dynamics. This can help the analysis and interpretation of power-law measures in EEG and MEG, and we believe our results can also impact computational modeling of brain dynamics, where different local connectivity structures could be used at different frequencies.

Fast Oscillatory Commands from the Motor Cortex Can Be Decoded by the Spinal Cord for Force Control

Oscillations in the beta and gamma bands (13-30 Hz; 35-70 Hz) have often been observed in motor cortical outputs that reach the spinal cord, acting on motoneurons and interneurons. However, the frequencies of these oscillations are above the muscle force frequency range. A current view is that the transformation of the motoneuron pool inputs into force is linear. For this reason possible roles for these oscillations are unclear, since if this transformation is linear, the high frequencies in the motoneuron inputs (e.g., 20 Hz from pyramidal tract neurons) would be filtered out by the muscle and have no effect on force control. A biologically inspired mathematical model of the neuromuscular system was used to investigate the impact of high-frequency cortical oscillatory activity on force control. The model simulation results evidenced that a typical motoneuron pool has a nonlinear behavior that enables the decoding of a high-frequency oscillatory input. An input at a single frequency (e.g., beta band) leads to an increase in the steady-state force generated by the muscle. When the input oscillation was amplitude modulated at a given low frequency, the force oscillated at this frequency. In both cases, the mechanism relies on the recruitment and derecruitment of motor units in response to the oscillatory descending drive. Therefore, the results from this study suggest a potential role in force control for cortical oscillations at frequencies at or above the beta band, despite the low-pass behavior of the muscles.

SIGNIFICANCE STATEMENT

The role of cortical oscillations in motor control has been a long-standing question, one view being that they are an epiphenomenon. Fast oscillations are known to reach the spinal cord, and hence they have been thought to affect muscle behavior. However, experimental limitations have hampered further advances to explain how they could influence muscle force. An approach for such a challenge was adopted in the present research: to study the problem through computer simulations of an advanced biologically compatible mathematical model. Using such a model, we found that the well-known mechanism of recruitment and derecruitment of the spinal cord motoneurons can allow the muscle to respond to cortical oscillations, suggesting that these oscillations are not epiphenomena in motor control.

A multi-channel, flex-rigid ECoG microelectrode array for visual cortical interfacing

High-density electrocortical (ECoG) microelectrode arrays are promising signal-acquisition platforms for brain-computer interfaces envisioned, e.g., as high-performance communication solutions for paralyzed persons. We propose a multi-channel microelectrode array capable of recording ECoG field potentials with high spatial resolution. The proposed array is of a 150 mm2 total recording area; it has 124 circular electrodes (100, 300 and 500 µm in diameter) situated on the edges of concentric hexagons (min. 0.8 mm interdistance) and a skull-facing reference electrode (2.5 mm2 surface area). The array is processed as a free-standing device to enable monolithic integration of a rigid interposer, designed for soldering of fine-pitch SMD-connectors on a minimal assembly area. Electrochemical characterization revealed distinct impedance spectral bands for the 100, 300 and 500 µm-type electrodes, and for the array's own reference. Epidural recordings from the primary visual cortex (V1) of an awake Rhesus macaque showed natural electrophysiological signals and clear responses to standard visual stimulation. The ECoG electrodes of larger surface area recorded signals with greater spectral power in the gamma band, while the skull-facing reference electrode provided higher average gamma power spectral density (γPSD) than the common average referencing technique.

Consequences of EEG electrode position error on ultimate beamformer source reconstruction performance

Inaccuracy of EEG electrode coordinates forms an error term in forward model generation and ultimate source reconstruction performance. This error arises from the combination of both intrinsic measurement noise of the digitization apparatus and manual coregistration error when selecting corresponding points on anatomical MRI volumes. A common assumption is that such an error would lead only to displacement of localized sources. Here, we measured electrode positions on a 3D-printed full-scale replica head, using three different techniques: a fringe projection 3D scanner, a novel "Flying Triangulation" 3D sensor, and a traditional electromagnetic digitizer. Using highly accurate fringe projection data as ground truth, the Flying Triangulation sensor had a mean error of 1.5 mm while the electromagnetic digitizer had a mean error of 6.8 mm. Then, again using the fringe projection as ground truth, individual EEG simulations were generated, with source locations across the brain space and a range of sensor noise levels. The simulated datasets were then processed using a beamformer in conjunction with the electrode coordinates registered with the Flying Triangulation and electromagnetic digitizer methods. The beamformer's output SNR was severely degraded with the digitizer-based positions but less severely with the Flying Triangulation coordinates. Therefore, the seemingly innocuous error in electrode registration may result in substantial degradation of beamformer performance, with output SNR penalties up to several decibels. In the case of low-SNR signals such as deeper brain structures or gamma band sources, this implies that sensor coregistration accuracy could make the difference between successful detection of such activity or complete failure to resolve the source.

Microscale inhomogeneity of brain tissue distorts electrical signal propagation

Interpretations of local field potentials (LFPs) are typically shaped on an assumption that the brain is a homogenous conductive milieu. However, microscale inhomogeneities including cell bodies, dendritic structures, axonal fiber bundles and blood vessels are unequivocally present and have different conductivities and permittivities than brain extracellular fluid. To determine the extent to which these obstructions affect electrical signal propagation on a microscale, we delivered electrical stimuli intracellularly to individual cells while simultaneously recording the extracellular potentials at different locations in a rat brain slice. As compared with relatively unobstructed paths, signals were attenuated across frequencies when fiber bundles were in between the stimulated cell and the extracellular electrode. Across group of cell bodies, signals were attenuated at low frequencies, but facilitated at high frequencies. These results show that LFPs do not reflect a democratic representation of neuronal contributions, as certain neurons may contribute to the LFP more than others based on the local extracellular environment surrounding them.