In ∼20 million people with drug-resistant epilepsy, focal seizures originating in dysfunctional brain networks will often evolve and spread to surrounding tissue, disrupting function in otherwise normal brain regions. To identify network control mechanisms that regulate seizure spread, we developed a novel tool for pinpointing brain regions that facilitate synchronization in the epileptic network. Our method measures the impact of virtually resecting putative control regions on synchronization in a validated model of the human epileptic network. By applying our technique to time-varying functional networks, we identified brain regions whose topological role is to synchronize or desynchronize the epileptic network. Our results suggest that greater antagonistic push-pull interaction between synchronizing and desynchronizing brain regions better constrains seizure spread. These methods, while applied here to epilepsy, are generalizable to other brain networks and have wide applicability in isolating and mapping functional drivers of brain dynamics in health and disease.

The epileptic network is characterized by pathologic, seizure-generating foci embedded in a web of structural and functional connections. Clinically, seizure foci are considered optimal targets for surgery. However, poor surgical outcome suggests a complex relationship between foci and the surrounding network that drives seizure dynamics. We developed a novel technique to objectively track seizure states from dynamic functional networks constructed from intracranial recordings. Each dynamical state captures unique patterns of network connections that indicate synchronized and desynchronized hubs of neural populations. Our approach suggests that seizures are generated when synchronous relationships near foci work in tandem with rapidly changing desynchronous relationships from the surrounding epileptic network. As seizures progress, topographical and geometrical changes in network connectivity strengthen and tighten synchronous connectivity near foci-a mechanism that may aid seizure termination. Collectively, our observations implicate distributed cortical structures in seizure generation, propagation and termination, and may have practical significance in determining which circuits to modulate with implantable devices.

Neural sensing and stimulation have been the backbone of neuroscience research, brain-machine interfaces and clinical neuromodulation therapies for decades. To-date, most of the neural stimulation systems have relied on sharp metal microelectrodes with poor electrochemical properties that induce extensive damage to the tissue and significantly degrade the long-term stability of implantable systems. Here, we demonstrate a flexible cortical microelectrode array based on porous graphene, which is capable of efficient electrophysiological sensing and stimulation from the brain surface, without penetrating into the tissue. Porous graphene electrodes show superior impedance and charge injection characteristics making them ideal for high efficiency cortical sensing and stimulation. They exhibit no physical delamination or degradation even after 1 million biphasic stimulation cycles, confirming high endurance. In in vivo experiments with rodents, same array is used to sense brain activity patterns with high spatio-temporal resolution and to control leg muscles with high-precision electrical stimulation from the cortical surface. Flexible porous graphene array offers a minimally invasive but high efficiency neuromodulation scheme with potential applications in cortical mapping, brain-computer interfaces, treatment of neurological disorders, where high resolution and simultaneous recording and stimulation of neural activity are crucial.

Human epilepsy patients suffer from spontaneous seizures, which originate in brain regions that also subserve normal function. Prior studies demonstrate focal, neocortical epilepsy is associated with dysfunction, several hours before seizures. How does the epileptic network perpetuate dysfunction during baseline periods? To address this question, we developed an unsupervised machine learning technique to disentangle patterns of functional interactions between brain regions, or subgraphs, from dynamic functional networks constructed from approximately 100 h of intracranial recordings in each of 22 neocortical epilepsy patients. Using this approach, we found: (1) subgraphs from ictal (seizure) and interictal (baseline) epochs are topologically similar, (2) interictal subgraph topology and dynamics can predict brain regions that generate seizures, and (3) subgraphs undergo slower and more coordinated fluctuations during ictal epochs compared to interictal epochs. Our observations suggest that seizures mark a critical shift away from interictal states that is driven by changes in the dynamical expression of strongly interacting components of the epileptic network.

Objective

Current brain stimulation paradigms are largely empirical rather than theoretical. An opportunity exists to improve upon their modest effectiveness in closed-loop control strategies with the development of theoretically grounded, model-based designs.

Approach

Inspired by this need, here we couple experimental data and mathematical modeling with a control-theoretic strategy for seizure termination. We begin by exercising a dynamical systems approach to model seizures (n = 94) recorded using intracranial EEG (iEEG) from 21 patients with medication-resistant, localization-related epilepsy.

Main results

Although each patient's seizures displayed unique spatial and temporal patterns, their evolution can be parsimoniously characterized by the same model form. Idiosyncracies of the model can inform individualized intervention strategies, specifically in iEEG samples with well-localized seizure onset zones. Temporal fluctuations in the spatial profiles of the oscillatory modes show that seizure onset marks a transition into a regime in which the underlying system supports prolonged rhythmic and focal activity. Based on these observations, we propose a control-theoretic strategy that aims to stabilize ictal activity using static output feedback for linear time-invariant switching systems. Finally, we demonstrate in silico that our proposed strategy allows us to dampen the emerging focal oscillatory sources using only a small set of electrodes.

Significance

Our integrative study informs the development of modulation and control algorithms for neurostimulation that could improve the effectiveness of implantable, closed-loop anti-epileptic devices.

Objective

Despite the overall success of responsive neurostimulation (RNS) therapy for drug-resistant focal epilepsy, clinical outcomes in individuals vary significantly and are hard to predict. Biomarkers that indicate the clinical efficacy of RNS-ideally before device implantation-are critically needed, but challenges include the intrinsic heterogeneity of the RNS patient population and variability in clinical management across epilepsy centers. The aim of this study is to use a multicenter dataset to evaluate a candidate biomarker from intracranial electroencephalographic (iEEG) recordings that predicts clinical outcome with subsequent RNS therapy.

Methods

We assembled a federated dataset of iEEG recordings, collected prior to RNS implantation, from a retrospective cohort of 30 patients across three major epilepsy centers. Using ictal iEEG recordings, each center independently calculated network synchronizability, a candidate biomarker indicating the susceptibility of epileptic brain networks to RNS therapy.

Results

Ictal measures of synchronizability in the high-γ band (95-105 Hz) significantly distinguish between good and poor RNS responders after at least 3 years of therapy under the current RNS therapy guidelines (area under the curve = .83). Additionally, ictal high-γ synchronizability is inversely associated with the degree of therapeutic response.

Significance

This study provides a proof-of-concept roadmap for collaborative biomarker evaluation in federated data, where practical considerations impede full data sharing across centers. Our results suggest that network synchronizability can help predict therapeutic response to RNS therapy. With further validation, this biomarker could facilitate patient selection and help avert a costly, invasive intervention in patients who are unlikely to benefit.

Optimizing direct electrical stimulation for the treatment of neurological disease remains difficult due to an incomplete understanding of its physical propagation through brain tissue. Here, we use network control theory to predict how stimulation spreads through white matter to influence spatially distributed dynamics. We test the theory's predictions using a unique dataset comprising diffusion weighted imaging and electrocorticography in epilepsy patients undergoing grid stimulation. We find statistically significant shared variance between the predicted activity state transitions and the observed activity state transitions. We then use an optimal control framework to posit testable hypotheses regarding which brain states and structural properties will efficiently improve memory encoding when stimulated. Our work quantifies the role that white matter architecture plays in guiding the dynamics of direct electrical stimulation and offers empirical support for the utility of network control theory in explaining the brain's response to stimulation.

Chronically implantable neurostimulation devices are becoming a clinically viable option for treating patients with neurological disease and psychiatric disorders. Neurostimulation offers the ability to probe and manipulate distributed networks of interacting brain areas in dysfunctional circuits. Here, we use tools from network control theory to examine the dynamic reconfiguration of functionally interacting neuronal ensembles during targeted neurostimulation of cortical and subcortical brain structures. By integrating multimodal intracranial recordings and diffusion-weighted imaging from patients with drug-resistant epilepsy, we test hypothesized structural and functional rules that predict altered patterns of synchronized local field potentials. We demonstrate the ability to predictably reconfigure functional interactions depending on stimulation strength and location. Stimulation of areas with structurally weak connections largely modulates the functional hubness of downstream areas and concurrently propels the brain towards more difficult-to-reach dynamical states. By using focal perturbations to bridge large-scale structure, function, and markers of behavior, our findings suggest that stimulation may be tuned to influence different scales of network interactions driving cognition.

Patients with drug-resistant epilepsy often require surgery to become seizure-free. While laser ablation and implantable stimulation devices have lowered the morbidity of these procedures, seizure-free rates have not dramatically improved, particularly for patients without focal lesions. This is in part because it is often unclear where to intervene in these cases. To address this clinical need, several research groups have published methods to map epileptic networks but applying them to improve patient care remains a challenge. In this study we advance clinical translation of these methods by: (i) presenting and sharing a robust pipeline to rigorously quantify the boundaries of the resection zone and determining which intracranial EEG electrodes lie within it; (ii) validating a brain network model on a retrospective cohort of 28 patients with drug-resistant epilepsy implanted with intracranial electrodes prior to surgical resection; and (iii) sharing all neuroimaging, annotated electrophysiology, and clinical metadata to facilitate future collaboration. Our network methods accurately forecast whether patients are likely to benefit from surgical intervention based on synchronizability of intracranial EEG (area under the receiver operating characteristic curve of 0.89) and provide novel information that traditional electrographic features do not. We further report that removing synchronizing brain regions is associated with improved clinical outcome, and postulate that sparing desynchronizing regions may further be beneficial. Our findings suggest that data-driven network-based methods can identify patients likely to benefit from resective or ablative therapy, and perhaps prevent invasive interventions in those unlikely to do so.

Bioresorbable silicon electronics technology offers unprecedented opportunities to deploy advanced implantable monitoring systems that eliminate risks, cost and discomfort associated with surgical extraction. Applications include postoperative monitoring and transient physiologic recording after percutaneous or minimally invasive placement of vascular, cardiac, orthopaedic, neural or other devices. We present an embodiment of these materials in both passive and actively addressed arrays of bioresorbable silicon electrodes with multiplexing capabilities, which record in vivo electrophysiological signals from the cortical surface and the subgaleal space. The devices detect normal physiologic and epileptiform activity, both in acute and chronic recordings. Comparative studies show sensor performance comparable to standard clinical systems and reduced tissue reactivity relative to conventional clinical electrocorticography (ECoG) electrodes. This technology offers general applicability in neural interfaces, with additional potential utility in treatment of disorders where transient monitoring and modulation of physiologic function, implant integrity and tissue recovery or regeneration are required.