Dynamic tractography-based localization of spike sources and animation of spike propagations

Cortical and white matter substrates supporting visuospatial working memory

OBJECTIVE

In tasks involving new visuospatial information, we rely on working memory, supported by a distributed brain network. We investigated the dynamic interplay between brain regions, including cortical and white matter structures, to understand how neural interactions change with different memory loads and trials, and their subsequent impact on working memory performance.

METHODS

Patients undertook a task of immediate spatial recall during intracranial EEG monitoring. We charted the dynamics of cortical high-gamma activity and associated functional connectivity modulations in white matter tracts.

RESULTS

Elevated memory loads were linked to enhanced functional connectivity via occipital longitudinal tracts, yet decreased through arcuate, uncinate, and superior-longitudinal fasciculi. As task familiarity grew, there was increased high-gamma activity in the posterior inferior-frontal gyrus (pIFG) and diminished functional connectivity across a network encompassing frontal, parietal, and temporal lobes. Early pIFG high-gamma activity was predictive of successful recall. Including this metric in a logistic regression model yielded an accuracy of 0.76.

CONCLUSIONS

Optimizing visuospatial working memory through practice is tied to early pIFG activation and decreased dependence on irrelevant neural pathways.

SIGNIFICANCE

This study expands our knowledge of human adaptation for visuospatial working memory, showing the spatiotemporal dynamics of cortical network modulations through white matter tracts.

Epilepsy Networks and Their Surgical Relevance

Surgical epilepsy is a rapidly evolved field. As the understanding and concepts of epilepsy shift towards a network disorder, surgical outcomes may shed light on numerous components of these systems. This review documents the evolution of the understanding of epilepsy networks and examines the data generated by resective, ablative, neuromodulation, and invasive monitoring surgeries in epilepsy patients. As these network tools are better integrated into epilepsy practice, they may eventually inform surgical decisions and improve clinical outcomes.

High frequency oscillations in relation to interictal spikes in predicting postsurgical seizure freedom

We evaluate whether interictal spikes, epileptiform HFOs and their co-occurrence (Spike + HFO) were included in the resection area with respect to seizure outcome. We also characterise the relationship between high frequency oscillations (HFOs) and propagating spikes. We analysed intracranial EEG of 20 patients that underwent resective epilepsy surgery. The co-occurrence of ripples and fast ripples was considered an HFO event; the co-occurrence of an interictal spike and HFO was considered a Spike + HFO event. HFO distribution and spike onset were compared in cases of spike propagation. Accuracy in predicting seizure outcome was 85% for HFO, 60% for Spikes, and 79% for Spike + HFO. Sensitivity was 57% for HFO, 71% for Spikes and 67% for Spikes + HFO. Specificity was 100% for HFO, 54% for Spikes and 85% for Spikes + HFO. In 2/2 patients with spike propagation, the spike onset included the HFO area. Combining interictal spikes with HFO had comparable accuracy to HFO. In patients with propagating spikes, HFO rate was maximal at the onset of spike propagation.

Identifying sources of human interictal discharges with travelling wave and white matter propagation

Interictal epileptiform discharges have been shown to propagate from focal epileptogenic sources as travelling waves or through more rapid white matter conduction. We hypothesize that both modes of propagation are necessary to explain interictal discharge timing delays. We propose a method that, for the first time, incorporates both propagation modes to identify unique potential sources of interictal activity. We retrospectively analysed 38 focal epilepsy patients who underwent intracranial EEG recordings and diffusion-weighted imaging for epilepsy surgery evaluation. Interictal discharges were detected and localized to the most likely source based on relative delays in time of arrival across electrodes, incorporating travelling waves and white matter propagation. We assessed the influence of white matter propagation on distance of spread, timing and clinical interpretation of interictal activity. To evaluate accuracy, we compared our source localization results to earliest spiking regions to predict seizure outcomes. White matter propagation helps to explain the timing delays observed in interictal discharge sequences, underlying rapid and distant propagation. Sources identified based on differences in time of receipt of interictal discharges are often distinct from the leading electrode location. Receipt of activity propagating rapidly via white matter can occur earlier than more local activity propagating via slower cortical travelling waves. In our cohort, our source localization approach was more accurate in predicting seizure outcomes than the leading electrode location. Inclusion of white matter in addition to travelling wave propagation in our model of discharge spread did not improve overall accuracy but allowed for identification of unique and at times distant potential sources of activity, particularly in patients with persistent postoperative seizures. Since distant white matter propagation can occur more rapidly than local travelling wave propagation, combined modes of propagation within an interictal discharge sequence can decouple the commonly assumed relationship between spike timing and distance from the source. Our findings thus highlight the clinical importance of recognizing the presence of dual modes of propagation during interictal discharges, as this may be a cause of clinical mislocalization.

Connectome-based approaches in pediatric epilepsy surgery: "State-of-the art" and future perspectives

Modern epilepsy science has overcome the traditional interpretation of a strict region-specific origin of epilepsy, highlighting the involvement of wider patterns of altered neuronal circuits. In selected cases, surgery may constitute a valuable option to achieve both seizure freedom and neurocognitive improvement. Although epilepsy is now considered as a brain network disease, the most relevant literature concerning the "connectome-based" epilepsy surgery mainly refers to adults, with a limited number of studies dedicated to the pediatric population. In this review, the Authors summarized the main current available knowledge on the relevance of WM surgical anatomy in epilepsy surgery, the post-surgical modifications of brain structural connectivity and the related clinical impact of such modifications within the pediatric context. In the last part, possible implications and future perspectives of this approach have been discussed, especially concerning the optimization of surgical strategies and the predictive value of the epilepsy network analysis for planning tailored approaches, with the final aim of improving case selection, presurgical planning, intraoperative management, and postoperative results.

Sevoflurane-induced high-frequency oscillations, effective connectivity and intraoperative classification of epileptic brain areas

OBJECTIVE

To determine how sevoflurane anesthesia modulates intraoperative epilepsy biomarkers on electrocorticography, including high-frequency oscillation (HFO) effective connectivity (EC), and to investigate their relation to epileptogenicity and anatomical white matter.

METHODS

We studied eight pediatric drug-resistant focal epilepsy patients who achieved seizure control after invasive monitoring and resective surgery. We visualized spatial distributions of the electrocorticography biomarkers at an oxygen baseline, three time-points while sevoflurane was increasing, and at a plateau of 2 minimum alveolar concentration (MAC) sevoflurane. HFO EC was combined with diffusion-weighted imaging, in dynamic tractography.

RESULTS

Intraoperative HFO EC diffusely increased as a function of sevoflurane concentration, although most in epileptogenic sites (defined as those included in the resection); their ability to classify epileptogenicity was optimized at sevoflurane 2 MAC. HFO EC could be visualized on major white matter tracts, as a function of sevoflurane level.

CONCLUSIONS

The results strengthened the hypothesis that sevoflurane-activated HFO biomarkers may help intraoperatively localize the epileptogenic zone.

SIGNIFICANCE

Our results help characterize how HFOs at non-epileptogenic and epileptogenic networks respond to sevoflurane. It may be warranted to establish a normative HFO atlas incorporating the modifying effects of sevoflurane and major white matter pathways, as critical reference in epilepsy presurgical evaluation.

Intra- and inter-hemispheric network dynamics supporting object recognition and speech production

We built normative brain atlases that animate millisecond-scale intra- and inter-hemispheric white matter-level connectivity dynamics supporting object recognition and speech production. We quantified electrocorticographic modulations during three naming tasks using event-related high-gamma activity from 1,114 nonepileptogenic intracranial electrodes (i.e., non-lesional areas unaffected by epileptiform discharges). Using this electrocorticography data, we visualized functional connectivity modulations defined as significant naming-related high-gamma modulations occurring simultaneously at two sites connected by direct white matter streamlines on diffusion-weighted imaging tractography. Immediately after stimulus onset, intra- and inter-hemispheric functional connectivity enhancements were confined mainly across modality-specific perceptual regions. During response preparation, left intra-hemispheric connectivity enhancements propagated in a posterior-to-anterior direction, involving the left precentral and prefrontal areas. After overt response onset, inter- and intra-hemispheric connectivity enhancements mainly encompassed precentral, postcentral, and superior-temporal (STG) gyri. We found task-specific connectivity enhancements during response preparation as follows. Picture naming enhanced activity along the left arcuate fasciculus between the inferior-temporal and precentral/posterior inferior-frontal (pIFG) gyri. Nonspeech environmental sound naming augmented functional connectivity via the left inferior longitudinal and fronto-occipital fasciculi between the medial-occipital and STG/pIFG. Auditory descriptive naming task enhanced usage of the left frontal U-fibers, involving the middle-frontal gyrus. Taken together, the commonly observed network enhancements include inter-hemispheric connectivity optimizing perceptual processing exerted in each hemisphere, left intra-hemispheric connectivity supporting semantic and lexical processing, and inter-hemispheric connectivity for symmetric oral movements during overt speech. Our atlases improve the currently available models of object recognition and speech production by adding neural dynamics via direct intra- and inter-hemispheric white matter tracts.

Propofol disrupts alpha dynamics in functionally distinct thalamocortical networks during loss of consciousness

During propofol-induced general anesthesia, alpha rhythms measured using electroencephalography undergo a striking shift from posterior to anterior, termed anteriorization, where the ubiquitous waking alpha is lost and a frontal alpha emerges. The functional significance of alpha anteriorization and the precise brain regions contributing to the phenomenon are a mystery. While posterior alpha is thought to be generated by thalamocortical circuits connecting nuclei of the sensory thalamus with their cortical partners, the thalamic origins of the propofol-induced alpha remain poorly understood. Here, we used human intracranial recordings to identify regions in sensory cortices where propofol attenuates a coherent alpha network, distinct from those in the frontal cortex where it amplifies coherent alpha and beta activities. We then performed diffusion tractography between these identified regions and individual thalamic nuclei to show that the opposing dynamics of anteriorization occur within two distinct thalamocortical networks. We found that propofol disrupted a posterior alpha network structurally connected with nuclei in the sensory and sensory associational regions of the thalamus. At the same time, propofol induced a coherent alpha oscillation within prefrontal cortical areas that were connected with thalamic nuclei involved in cognition, such as the mediodorsal nucleus. The cortical and thalamic anatomy involved, as well as their known functional roles, suggests multiple means by which propofol dismantles sensory and cognitive processes to achieve loss of consciousness.

Integration of white matter architecture to stereo-EEG better describes epileptic spike propagation

OBJECTIVE

Stereo-electroencephalography (SEEG)-derived epilepsy networks are used to better understand a patient's epilepsy; however, a unimodal approach provides an incomplete picture. We combine tractography and SEEG to determine the relationship between spike propagation and the white matter architecture and to improve our understanding of spike propagation mechanisms.

METHODS

Probablistic tractography from diffusion imaging (dMRI) of matched subjects from the Human Connectome Project (HCP) was combined with patient-specific SEEG-derived spike propagation networks. Two regions-of-interest (ROIs) with a significant spike propagation relationship constituted a Propagation Pair.

RESULTS

In 56 of 59 patients, Propagation Pairs were more often tract-connected as compared to all ROI pairs (p < 0.01; d = -1.91). The degree of spike propagation between tract-connected ROIs was greater (39 ± 21%) compared to tract-unconnected ROIs (31 ± 18%; p < 0.0001). Within the same network, ROIs receiving propagation earlier were more often tract-connected to the source (59.7%) as compared to late receivers (25.4%; p < 0.0001).

CONCLUSIONS

Brain regions involved in spike propagation are more likely to be connected by white matter tracts. Between nodes, presence of tracts suggests a direct course of propagation, whereas the absence of tracts suggests an indirect course of propagation.

SIGNIFICANCE

We demonstrate a logical and consistent relationship between spike propagation and the white matter architecture.

Multiple Sources of Fast Traveling Waves during Human Seizures: Resolving a Controversy

During human seizures, organized waves of voltage activity rapidly sweep across the cortex. Two contradictory theories describe the source of these fast traveling waves: either a slowly advancing narrow region of multiunit activity (an ictal wavefront) or a fixed cortical location. Limited observations and different analyses prevent resolution of these incompatible theories. Here we address this disagreement by combining the methods and microelectrode array recordings (N = 11 patients, 2 females, N = 31 seizures) from previous human studies to analyze the traveling wave source. We find, inconsistent with both existing theories, a transient relationship between the ictal wavefront and traveling waves, and multiple stable directions of traveling waves in many seizures. Using a computational model that combines elements of both existing theories, we show that interactions between an ictal wavefront and fixed source reproduce the traveling wave dynamics observed in vivo We conclude that combining both existing theories can generate the diversity of ictal traveling waves.SIGNIFICANCE STATEMENT The source of voltage discharges that propagate across cortex during human seizures remains unknown. Two candidate theories exist, each proposing a different discharge source. Support for each theory consists of observations from a small number of human subject recordings, analyzed with separately developed methods. How the different, limited data and different analysis methods impact the evidence for each theory is unclear. To resolve these differences, we combine the unique, human microelectrode array recordings collected separately for each theory and analyze these combined data with a unified approach. We show that neither existing theory adequately describes the data. We then propose a new theory that unifies existing proposals and successfully reproduces the voltage discharge dynamics observed in vivo.

Diffusion tractography predicts propagated high-frequency activity during epileptic spasms

OBJECTIVE

Determine the structural networks that constrain propagation of ictal oscillations during epileptic spasm events, and compare observed propagation patterns across patients with successful or unsuccessful surgical outcomes.

METHODS

Subdural electrode recordings of 18 young patients (age 1-11 years) were analyzed during epileptic spasm events to determine ictal networks and quantify the amplitude and onset time of ictal oscillations across the cortical surface. Corresponding structural networks were generated with diffusion MRI tractography by seeding the cortical region associated with the earliest average oscillation onset time, and white matter pathways connecting active electrode regions within the ictal network were isolated. Properties of this structural network were used to predict oscillation onset times and amplitudes, and this relationship was compared across patients who did and did not achieve seizure freedom following resective surgery.

RESULTS

Onset propagation patterns were relatively consistent across each patients' spasm events. An electrode's average ictal oscillation onset latency was most significantly associated with the length of direct corticocortical tracts connecting to the area with the earliest average oscillation onset (p < .001, model R²  = 0.54). Moreover, patients demonstrating a faster propagation of ictal oscillation signals within the corticocortical network were more likely to have seizure recurrence following resective surgery (p = .039). Ictal oscillation amplitude was also associated with connecting tractography length and weighted fractional anisotropy (FA) measures along these pathways (p = .002/.030, model R²  = 0.31/0.25). Characteristics of analogous corticothalamic pathways did not show significant associations with ictal oscillation onset latency or amplitude.

SIGNIFICANCE

Spatiotemporal propagation patterns of high-frequency activity in epileptic spasms align with length and FA measures from onset-originating corticocortical pathways. Considering data in this individualized framework may help inform surgical decision making and expectations of surgical outcomes.

Temporally and functionally distinct large-scale brain network dynamics supporting task switching

Objective:

Our daily activities require frequent switches among competing responses at the millisecond time scale. We determined the spatiotemporal characteristics and functional significance of rapid, large-scale brain network dynamics during task switching.

Methods:

This cross-sectional study investigated patients with drug-resistant focal epilepsy who played a Lumosity cognitive flexibility training game during intracranial electroencephalography (iEEG) recording. According to a given task rule, unpredictably switching across trials, participants had to swipe the screen in the direction the stimulus was pointing or moving. Using this data, we described the spatiotemporal characteristics of iEEG high-gamma augmentation occurring more intensely during switch than repeat trials, unattributable to the effect of task rule (pointing or moving), within-stimulus congruence (the direction of stimulus pointing and moving was same or different in a given trial), or accuracy of an immediately preceding response. Diffusion-weighted imaging (DWI) tractography determined whether distant cortical regions showing enhanced activation during task switch trials were directly connected by white matter tracts. Trial-by-trial iEEG analysis deduced whether the intensity of task switch-related high-gamma augmentation was altered through practice and whether high-gamma amplitude predicted the accuracy of an upcoming response among switch trials.

Results:

The average number of completed trials during five-minute gameplay was 221.4 per patient (range: 171–285). Task switch trials increased the response times, whereas later trials reduced them. Analysis of iEEG signals sampled from 860 brain sites effectively elucidated the distinct spatiotemporal characteristics of task switch, task rule, and post-error-specific high-gamma modulations. Post-cue, task switch-related high-gamma augmentation was initiated in the right calcarine cortex after 260 ms, right precuneus after 330 ms, right entorhinal after 420 ms, and bilateral anterior middle-frontal gyri after 450 ms. DWI tractography successfully showed the presence of direct white matter tracts connecting the right visual areas to the precuneus and anterior middle-frontal regions but not between the right precuneus and anterior middle-frontal regions. Task-related high-gamma amplitudes in later trials were reduced in the calcarine, entorhinal and anterior middle-frontal regions, but increased in the precuneus. Functionally, enhanced post-cue precuneus high-gamma augmentation improved the accuracy of subsequent responses among switch trials.

Conclusions:

Our multimodal analysis uncovered two temporally and functionally distinct network dynamics supporting task switching. High-gamma augmentation in the visual-precuneus pathway may reflect the neural process facilitating an attentional shift to a given updated task rule. High-gamma activity in the visual-dorsolateral prefrontal pathway, rapidly reduced through practice, may reflect the cost of executing appropriate stimulus-response translation.

OUP accepted manuscript

Presurgical identification of the epileptogenic zone is a critical determinant of seizure control following surgical resection in epilepsy. Excitatory glutamate α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid receptor is a major component of neurotransmission. Although elevated α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid receptor levels are observed in surgically resected brain areas of patients with epilepsy, it remains unclear whether increased α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid receptor-mediated currents initiate epileptic discharges. We have recently developed the first PET tracer for α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid receptor, [¹¹C]K-2, to visualize and quantify the density of α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid receptors in living human brains. Here, we detected elevated [¹¹C]K-2 uptake in the epileptogenic temporal lobe of patients with mesial temporal lobe epilepsy. Brain areas with high [¹¹C]K-2 uptake are closely colocalized with the location of equivalent current dipoles estimated by magnetoencephalography or with seizure onset zones detected by intracranial electroencephalogram. These results suggest that epileptic discharges initiate from brain areas with increased α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid receptors, providing a biological basis for epileptic discharges and an additional non-invasive option to identify the epileptogenic zone in patients with mesial temporal lobe epilepsy.