State-dependent effects of responsive neurostimulation depend on seizure localization.

Brain-responsive neurostimulation is firmly ensconced among treatment options for drug-resistant focal epilepsy, but over a quarter of patients treated with the RNS System do not experience meaningful seizure reduction. Initial titration of RNS therapy is typically similar for all patients, raising the possibility that treatment response might be enhanced by consideration of patient-specific variables. Indeed, small, single-center studies have yielded preliminary evidence that RNS System effectiveness depends on the brain state during which stimulation is applied. The generalizability of these findings remains unclear, however, and it is unknown whether state-dependent effects of responsive neurostimulation are also stratified by location of the seizure onset zone where stimulation is delivered. We aimed to determine whether state-dependent effects of the RNS System are evident in the large, diverse, multi-center cohort of RNS System clinical trial participants and to test whether these effects differ between mesiotemporal and neocortical epilepsies. Eighty-one of 256 patients who were treated with the RNS System across 31 centers during clinical trials met criteria for inclusion in this retrospective study. Risk states were defined in relation to phases of daily and multi-day cycles of interictal epileptiform activity that are thought to determine seizure likelihood. We found that the probabilities of risk state transitions depended on the stimulation parameter being changed, the starting seizure risk state, and the stimulated brain region. Changes in two commonly adjusted stimulation parameters, charge density and stimulation frequency, produced opposite effects on risk state transitions depending on seizure localization. Greater variance in acute risk state transitions was explained by state-dependent responsive neurostimulation for bipolar stimulation for neocortical epilepsies and for monopolar stimulation for mesiotemporal epilepsies. Variability in effectiveness of RNS System therapy across individuals may relate, at least partly, to the fact that current treatment paradigms do not account fully for fluctuations in brain states or locations of simulation sites. State-dependence of electrical brain stimulation may inform development of next-generation closed-loop devices that can detect changes in brain state and deliver adaptive, localization-specific patterns of stimulation to maximize therapeutic effects.

Intracranial electroencephalographic biomarker predicts effective responsive neurostimulation for epilepsy prior to treatment.

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.

Early seizure spread and epilepsy surgery: A systematic review.

OBJECTIVE: A fundamental question in epilepsy surgery is how to delineate the margins of cortex that must be resected to result in seizure freedom. Whether and which areas showing seizure activity early in ictus must be removed to avoid postoperative recurrence of seizures is an area of ongoing research. Seizure spread dynamics in the initial seconds of ictus are often correlated with postoperative outcome; there is neither a consensus definition of early spread nor a concise summary of the existing literature linking seizure spread to postsurgical seizure outcomes. The present study is intended to summarize the literature that links seizure spread to postoperative seizure outcome and to provide a framework for quantitative assessment of early seizure spread. METHODS: A systematic review was carried out according to Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines. A Medline search identified clinical studies reporting data on seizure spread measured by intracranial electrodes, having at least 10 subjects and reporting at least 1-year postoperative outcome in the English literature from 1990 to 2019. Studies were evaluated regarding support for a primary hypothesis: Areas of early seizure spread represent cortex with seizure-generating potential. RESULTS: The search yielded 4562 studies: 15 studies met inclusion criteria and 7 studies supported the primary hypothesis. The methods and metrics used to describe seizure spread were heterogenous. The timeframe of seizure spread associated with seizure outcome ranged from 1-14 seconds, with large, well-designed, retrospective studies pointing to 3-10 seconds as most likely to provide meaningful correlates of postoperative seizure freedom. SIGNIFICANCE: The complex correlation between electrophysiologic seizure spread and the potential for seizure generation needs further elucidation. Prospective cohort studies or trials are needed to evaluate epilepsy surgery targeting cortex involved in the first 3-10 seconds of ictus.

Virtual resection predicts surgical outcome for drug-resistant epilepsy.

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.

Correction: Subgraphs of functional brain networks identify dynamical constraints of cognitive control.

[This corrects the article DOI: 10.1371/journal.pcbi.1006234.].

Subgraphs of functional brain networks identify dynamical constraints of cognitive control.

Brain anatomy and physiology support the human ability to navigate a complex space of perceptions and actions. To maneuver across an ever-changing landscape of mental states, the brain invokes cognitive control-a set of dynamic processes that engage and disengage different groups of brain regions to modulate attention, switch between tasks, and inhibit prepotent responses. Current theory posits that correlated and anticorrelated brain activity may signify cooperative and competitive interactions between brain areas that subserve adaptive behavior. In this study, we use a quantitative approach to identify distinct topological motifs of functional interactions and examine how their expression relates to cognitive control processes and behavior. In particular, we acquire fMRI BOLD signal in twenty-eight healthy subjects as they perform two cognitive control tasks-a Stroop interference task and a local-global perception switching task using Navon figures-each with low and high cognitive control demand conditions. Based on these data, we construct dynamic functional brain networks and use a parts-based, network decomposition technique called non-negative matrix factorization to identify putative cognitive control subgraphs whose temporal expression captures distributed network structures involved in different phases of cooperative and competitive control processes. Our results demonstrate that temporal expression of the subgraphs fluctuate alongside changes in cognitive demand and are associated with individual differences in task performance. These findings offer insight into how coordinated changes in the cooperative and competitive roles of cognitive systems map trajectories between cognitively demanding brain states.

Beyond modularity: Fine-scale mechanisms and rules for brain network reconfiguration.

The human brain is in constant flux, as distinct areas engage in transient communication to support basic behaviors as well as complex cognition. The collection of interactions between cortical and subcortical areas forms a functional brain network whose topology evolves with time. Despite the nontrivial dynamics that are germane to this networked system, experimental evidence demonstrates that functional interactions organize into putative brain systems that facilitate different facets of cognitive computation. We hypothesize that such dynamic functional networks are organized around a set of rules that constrain their spatial architecture - which brain regions may functionally interact - and their temporal architecture - how these interactions fluctuate over time. To objectively uncover these organizing principles, we apply an unsupervised machine learning approach called non-negative matrix factorization to time-evolving, resting state functional networks in 20 healthy subjects. This machine learning approach automatically groups temporally co-varying functional interactions into subgraphs that represent putative topological modes of dynamic functional architecture. We find that subgraphs are stratified based on both the underlying modular organization and the topographical distance of their strongest interactions: while many subgraphs are largely contained within modules, others span between modules and are expressed differently over time. The relationship between dynamic subgraphs and modular architecture is further highlighted by the ability of time-varying subgraph expression to explain inter-individual differences in module reorganization. Collectively, these results point to the critical role that subgraphs play in constraining the topography and topology of functional brain networks. More broadly, this machine learning approach opens a new door for understanding the architecture of dynamic functional networks during both task and rest states, and for probing alterations of that architecture in disease.

Modeling and interpreting mesoscale network dynamics.

Recent advances in brain imaging techniques, measurement approaches, and storage capacities have provided an unprecedented supply of high temporal resolution neural data. These data present a remarkable opportunity to gain a mechanistic understanding not just of circuit structure, but also of circuit dynamics, and its role in cognition and disease. Such understanding necessitates a description of the raw observations, and a delineation of computational models and mathematical theories that accurately capture fundamental principles behind the observations. Here we review recent advances in a range of modeling approaches that embrace the temporally-evolving interconnected structure of the brain and summarize that structure in a dynamic graph. We describe recent efforts to model dynamic patterns of connectivity, dynamic patterns of activity, and patterns of activity atop connectivity. In the context of these models, we review important considerations in statistical testing, including parametric and non-parametric approaches. Finally, we offer thoughts on careful and accurate interpretation of dynamic graph architecture, and outline important future directions for method development.

Emerging Frontiers of Neuroengineering: A Network Science of Brain Connectivity.

Neuroengineering is faced with unique challenges in repairing or replacing complex neural systems that are composed of many interacting parts. These interactions form intricate patterns over large spatiotemporal scales and produce emergent behaviors that are difficult to predict from individual elements. Network science provides a particularly appropriate framework in which to study and intervene in such systems by treating neural elements (cells, volumes) as nodes in a graph and neural interactions (synapses, white matter tracts) as edges in that graph. Here, we review the emerging discipline of network neuroscience, which uses and develops tools from graph theory to better understand and manipulate neural systems from micro- to macroscales. We present examples of how human brain imaging data are being modeled with network analysis and underscore potential pitfalls. We then highlight current computational and theoretical frontiers and emphasize their utility in informing diagnosis and monitoring, brain-machine interfaces, and brain stimulation. A flexible and rapidly evolving enterprise, network neuroscience provides a set of powerful approaches and fundamental insights that are critical for the neuroengineer's tool kit.

Crowdsourcing seizure detection: algorithm development and validation on human implanted device recordings.

There exist significant clinical and basic research needs for accurate, automated seizure detection algorithms. These algorithms have translational potential in responsive neurostimulation devices and in automatic parsing of continuous intracranial electroencephalography data. An important barrier to developing accurate, validated algorithms for seizure detection is limited access to high-quality, expertly annotated seizure data from prolonged recordings. To overcome this, we hosted a kaggle.com competition to crowdsource the development of seizure detection algorithms using intracranial electroencephalography from canines and humans with epilepsy. The top three performing algorithms from the contest were then validated on out-of-sample patient data including standard clinical data and continuous ambulatory human data obtained over several years using the implantable NeuroVista seizure advisory system. Two hundred teams of data scientists from all over the world participated in the kaggle.com competition. The top performing teams submitted highly accurate algorithms with consistent performance in the out-of-sample validation study. The performance of these seizure detection algorithms, achieved using freely available code and data, sets a new reproducible benchmark for personalized seizure detection. We have also shared a 'plug and play' pipeline to allow other researchers to easily use these algorithms on their own datasets. The success of this competition demonstrates how sharing code and high quality data results in the creation of powerful translational tools with significant potential to impact patient care.

Dynamic Network Drivers of Seizure Generation, Propagation and Termination in Human Neocortical Epilepsy.

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.

Hippocampal network activity forecasts epileptic seizures.

Seizures in people with epilepsy were long thought to occur at random, but recent methods for seizure forecasting enable estimation of the likelihood of seizure occurrence over short horizons. These methods rely on days-long cyclical patterns of brain electrical activity and other physiological variables that determine seizure likelihood and that require measurement through long-term, multimodal recordings. In this retrospective cohort study of 15 adults with bitemporal epilepsy who had a device that provides chronic intracranial recordings, functional connectivity of hippocampal networks fluctuated in multiday cycles with patterns that mirrored cycles of seizure likelihood. A functional connectivity biomarker of seizure likelihood derived from 90-s recordings of background hippocampal activity generalized across individuals and forecasted 24-h seizure likelihood as accurately as cycle-based models requiring months-long baseline recordings. Larger, prospective studies are needed to validate this approach, but our results have the potential to make reliable seizure forecasts accessible to more people with epilepsy.

Closed-loop neurostimulation for the treatment of psychiatric disorders.

Despite increasing prevalence and huge personal and societal burden, psychiatric diseases still lack treatments which can control symptoms for a large fraction of patients. Increasing insight into the neurobiology underlying these diseases has demonstrated wide-ranging aberrant activity and functioning in multiple brain circuits and networks. Together with varied presentation and symptoms, this makes one-size-fits-all treatment a challenge. There has been a resurgence of interest in the use of neurostimulation as a treatment for psychiatric diseases. Initial studies using continuous open-loop stimulation, in which clinicians adjusted stimulation parameters during patient visits, showed promise but also mixed results. Given the periodic nature and fluctuations of symptoms often observed in psychiatric illnesses, the use of device-driven closed-loop stimulation may provide more effective therapy. The use of a biomarker, which is correlated with specific symptoms, to deliver stimulation only during symptomatic periods allows for the personalized therapy needed for such heterogeneous disorders. Here, we provide the reader with background motivating the use of closed-loop neurostimulation for the treatment of psychiatric disorders. We review foundational studies of open- and closed-loop neurostimulation for neuropsychiatric indications, focusing on deep brain stimulation, and discuss key considerations when designing and implementing closed-loop neurostimulation.

Closed-Loop Neurostimulation for Biomarker-Driven, Personalized Treatment of Major Depressive Disorder.

Deep brain stimulation involves the administration of electrical stimulation to targeted brain regions for therapeutic benefit. In the context of major depressive disorder (MDD), most studies to date have administered continuous or open-loop stimulation with promising but mixed results. One factor contributing to these mixed results may stem from when the stimulation is applied. Stimulation administration specific to high-symptom states in a personalized and responsive manner may be more effective at reducing symptoms compared to continuous stimulation and may avoid diminished therapeutic effects related to habituation. Additionally, a lower total duration of stimulation per day is advantageous for reducing device energy consumption. This protocol describes an experimental workflow using a chronically implanted neurostimulation device to achieve closed-loop stimulation for individuals with treatment-refractory MDD. This paradigm hinges on determining a patient-specific neural biomarker that is related to states of high symptoms and programming the device detectors, such that stimulation is triggered by this read-out of symptom state. The described procedures include how to obtain neural recordings concurrent with patient symptom reports, how to use these data in a state-space model approach to differentiate low- and high-symptom states and corresponding neural features, and how to subsequently program and tune the device to deliver closed-loop stimulation therapy.

Intracranial electrical stimulation of corticolimbic sites modulates arousal in humans.

BACKGROUND: Humans routinely shift their sleepiness and wakefulness levels in response to emotional factors. The diversity of emotional factors that modulates sleep-wake levels suggests that the ascending arousal network may be intimately linked with networks that mediate mood. Indeed, while animal studies have identified select limbic structures that play a role in sleep-wake regulation, the breadth of corticolimbic structures that directly modulates arousal in humans remains unknown. OBJECTIVE: We investigated whether select regional activation of the corticolimbic network through direct electrical stimulation can modulate sleep-wake levels in humans, as measured by subjective experience and behavior. METHODS: We performed intensive inpatient stimulation mapping in two human participants with treatment resistant depression, who underwent intracranial implantation with multi-site, bilateral depth electrodes. Stimulation responses of sleep-wake levels were measured by subjective surveys (i.e. Stanford Sleepiness Scale and visual-analog scale of energy) and a behavioral arousal score. Biomarker analyses of sleep-wake levels were performed by assessing spectral power features of resting-state electrophysiology. RESULTS: Our findings demonstrated three regions whereby direct stimulation modulated arousal, including the orbitofrontal cortex (OFC), subgenual cingulate (SGC), and, most robustly, ventral capsule (VC). Modulation of sleep-wake levels was frequency-specific: 100Hz OFC, SGC, and VC stimulation promoted wakefulness, whereas 1Hz OFC stimulation increased sleepiness. Sleep-wake levels were correlated with gamma activity across broad brain regions. CONCLUSIONS: Our findings provide evidence for the overlapping circuitry between arousal and mood regulation in humans. Furthermore, our findings open the door to new treatment targets and the consideration of therapeutic neurostimulation for sleep-wake disorders.

Decoding naturalistic affective behaviour from spectro-spatial features in multiday human iEEG.

The neurological basis of affective behaviours in everyday life is not well understood. We obtained continuous intracranial electroencephalography recordings from the human mesolimbic network in 11 participants with epilepsy and hand-annotated spontaneous behaviours from 116 h of multiday video recordings. In individual participants, binary random forest models decoded affective behaviours from neutral behaviours with up to 93% accuracy. Both positive and negative affective behaviours were associated with increased high-frequency and decreased low-frequency activity across the mesolimbic network. The insula, amygdala, hippocampus and anterior cingulate cortex made stronger contributions to affective behaviours than the orbitofrontal cortex, but the insula and anterior cingulate cortex were most critical for differentiating behaviours with observable affect from those without. In a subset of participants (N = 3), multiclass decoders distinguished amongst the positive, negative and neutral behaviours. These results suggest that spectro-spatial features of brain activity in the mesolimbic network are associated with affective behaviours of everyday life.

Long-term brain network reorganization predicts responsive neurostimulation outcomes for focal epilepsy.

Responsive neurostimulation (RNS) devices, able to detect imminent seizures and to rapidly deliver electrical stimulation to the brain, are effective in reducing seizures in some patients with focal epilepsy. However, therapeutic response to RNS is often slow, is highly variable, and defies prognostication based on clinical factors. A prevailing view holds that RNS efficacy is primarily mediated by acute seizure termination; yet, stimulations greatly outnumber seizures and occur mostly in the interictal state, suggesting chronic modulation of brain networks that generate seizures. Here, using years-long intracranial neural recordings collected during RNS therapy, we found that patients with the greatest therapeutic benefit undergo progressive, frequency-dependent reorganization of interictal functional connectivity. The extent of this reorganization scales directly with seizure reduction and emerges within the first year of RNS treatment, enabling potential early prediction of therapeutic response. Our findings reveal a mechanism for RNS that involves network plasticity and may inform development of next-generation devices for epilepsy.

Effects of anterior thalamic nuclei stimulation on hippocampal activity: Chronic recording in a patient with drug-resistant focal epilepsy.

Implanted neurostimulation devices are gaining traction as palliative treatment options for certain forms of drug-resistant epilepsy, but clinical utility of these devices is hindered by incomplete mechanistic understanding of their therapeutic effects. Approved devices for anterior thalamic nuclei deep brain stimulation (ANT DBS) are thought to work at a network level, but limited sensing capability precludes characterization of neurophysiological effects outside the thalamus. Here, we describe a patient with drug-resistant temporal lobe epilepsy who was implanted with a responsive neurostimulation device (RNS System), involving hippocampal and ipsilateral temporal neocortical leads, and subsequently received ANT DBS. Over 1.5 years, RNS System electrocorticography enabled multiscale characterization of neurophysiological effects of thalamic stimulation. In brain regions sampled by the RNS System, ANT DBS produced acute, phasic, frequency-dependent responses, including suppression of hippocampal low frequency local field potentials. ANT DBS modulated functional connectivity between hippocampus and neocortex. Finally, ANT DBS progressively suppressed hippocampal epileptiform activity in relation to the extent of hippocampal theta suppression, which informs stimulation parameter selection for ANT DBS. Taken together, this unique clinical scenario, involving hippocampal recordings of unprecedented chronicity alongside ANT DBS, sheds light on the therapeutic mechanism of thalamic stimulation and highlights capabilities needed in next-generation devices.

Evidence of state-dependence in the effectiveness of responsive neurostimulation for seizure modulation.

BACKGROUND: An implanted device for brain-responsive neurostimulation (RNS® System) is approved as an effective treatment to reduce seizures in adults with medically-refractory focal epilepsy. Clinical trials of the RNS System demonstrate population-level reduction in average seizure frequency, but therapeutic response is highly variable. HYPOTHESIS: Recent evidence links seizures to cyclical fluctuations in underlying risk. We tested the hypothesis that effectiveness of responsive neurostimulation varies based on current state within cyclical risk fluctuations. METHODS: We analyzed retrospective data from 25 adults with medically-refractory focal epilepsy implanted with the RNS System. Chronic electrocorticography was used to record electrographic seizures, and hidden Markov models decoded seizures into fluctuations in underlying risk. State-dependent associations of RNS System stimulation parameters with changes in risk were estimated. RESULTS: Higher charge density was associated with improved outcomes, both for remaining in a low seizure risk state and for transitioning from a high to a low seizure risk state. The effect of stimulation frequency depended on initial seizure risk state: when starting in a low risk state, higher stimulation frequencies were associated with remaining in a low risk state, but when starting in a high risk state, lower stimulation frequencies were associated with transition to a low risk state. Findings were consistent across bipolar and monopolar stimulation configurations. CONCLUSION: The impact of RNS on seizure frequency exhibits state-dependence, such that stimulation parameters which are effective in one seizure risk state may not be effective in another. These findings represent conceptual advances in understanding the therapeutic mechanism of RNS, and directly inform current practices of RNS tuning and the development of next-generation neurostimulation systems.

Distributed Subnetworks of Depression Defined by Direct Intracranial Neurophysiology.

Major depressive disorder is a common and disabling disorder with high rates of treatment resistance. Evidence suggests it is characterized by distributed network dysfunction that may be variable across patients, challenging the identification of quantitative biological substrates. We carried out this study to determine whether application of a novel computational approach to a large sample of high spatiotemporal resolution direct neural recordings in humans could unlock the functional organization and coordinated activity patterns of depression networks. This group level analysis of depression networks from heterogenous intracranial recordings was possible due to application of a correlational model-based method for inferring whole-brain neural activity. We then applied a network framework to discover brain dynamics across this model that could classify depression. We found a highly distributed pattern of neural activity and connectivity across cortical and subcortical structures that was present in the majority of depressed subjects. Furthermore, we found that this depression signature consisted of two subnetworks across individuals. The first was characterized by left temporal lobe hypoconnectivity and pathological beta activity. The second was characterized by a hypoactive, but hyperconnected left frontal cortex. These findings have applications toward personalization of therapy.