Brain tissue can be stimulated using wireless implants that comprise flexible substrates and electronics thereon, and include stimulating electrodes and recorder circuitry for monitoring brain activity.
People with epilepsy face the challenge of unpredictable seizures, significantly affecting their daily life quality. Even those resistant to medication spend less than 1% of their time experiencing a seizure, while the remainder is spent under constant seizure threat. Currently, people with epilepsy rely on daily anti-seizure medications that often have short-lived effects.
This document describes systems, methods, devices, and other techniques used for seizure detection, intervention, and prevention.
One aspect is a method for predicting and managing seizures in an individual, the method comprising obtaining at least one measurement of an impedance of brain tissue of the individual to determine an impedance signature of the brain tissue of the individual, processing the impedance signature of the brain tissue to predict whether the individual is likely to experience a seizure, and in response to predicting based on the impedance signature that the individual is likely to experience a seizure, initiating a remedial action that comprises at least one of alerting the individual or another user that the individual is likely to experience the seizure, applying an electrical stimulus to the brain tissue of the individual to reduce a likelihood of the seizure occurring or to reduce a severity of the seizure, or delivering a drug to the brain tissue of the individual to reduce a likelihood of the seizure occurring or to reduce a severity of the seizure.
Another aspect is a system comprising an implanted device configured to capture at least one measurement of an impedance of brain tissue of an individual and a computing device configured to obtain, from the implanted device, the at least one measurement of the impedance of the brain tissue of the individual and in response to a prediction that the individual is likely to experience a seizure, initiate a remedial action that comprises at least one of alerting the individual or another user that the individual is likely to experience the seizure, applying an electrical stimulus to the brain tissue of the individual to reduce a likelihood of the seizure occurring or to reduce a severity of the seizure, or delivering a drug to the brain tissue of the individual to reduce a likelihood of the seizure occurring or to reduce a severity of the seizure, wherein the prediction that the individual is likely to experience a seizure is determined by processing an impedance signature determined from the at least one measurement of the impedance of the brain tissue of the individual.
Like reference symbols in the various drawings indicate like elements.
This specification generally describes systems, methods, devices, and other techniques used for seizure detection, intervention, and prevention. In some implementations, systems and methods for brain impedance sensing for seizure detection, intervention, and prevention are disclosed. In some examples, brain impedance is periodically sampled at one or more locations over a long-term period. The detected impedance dynamics information is integrated into a system for seizure detection and forecasting. In some implementations, deep-learning models are used. In some implementations, adaptive therapy techniques are implemented and applied.
In some examples, the mechanisms underlying changes in extracellular space, impedance dynamics, and seizure generation are analyzed. Impedance measurements are then applied as measurements are detected to make seizure predictions. For example, to make predictions in real-time or near-real time.
In some aspects, systems and methods using brain impedance for seizure forecasting are disclosed. As described below, impedance has a circadian cycle. In some of the seizure forecasting methods disclosed herein, detected deviations from an individual patient's impedance circadian rhythm is used independently or in conjunction with other data input. Examples of other data input include iEEG data including spike detection, high frequency oscillations, or evoked potentials (induced by electrical or sensory stimulus). In some implementations, impedance measurements are used for sleep staging, clearing extracellular protein aggregates, and seizure detection and forecasting. In some implementations, data from wearables, data from implanted EEG, subscalp EEG, or scalp EEG (e.g. evoked potentials, high frequency oscillations, local field potentials, spike detection) and impedance can be combined into a holistic, real-time monitoring of brain excitability and seizure risk.
In some examples, impedance of brain tissue is measured using 2- or 4-point measurements. In some implementations, the measurements are taken local to a location where the electrical stimulation is applied. In some implementations, the measurements are taken from tissue that is not local to where the electrical stimulation is applied. For example, the network response may be measured from stimulation at one node and measured potential at a different node.
In some implementations, macro (mm scale), meso (0.1 mm scale), and/or micro (0.001 mm scale) electrodes are used to stimulate and record from contacts. In some implementations, a wide dynamic range for stimulations and sensing is used. In some particular examples the wide dynamic range includes values from 0.01 to 10,000 Hz. However, other values may also be used. In some implementations, the sampling frequency ranges from 100 Hz to 40,000 Hz for sensing. Other ranges and values can also be used. A range of stimulation currents for impedance can be evoked. For example, a range of stimulation currents from 0.05 to 5 mA can be used in some examples (other ranges are also possible). Additionally or alternatively, different stimulation paradigms can be applied. For example, white noise stimulation, sinusoidal current, square pulses, etc. In some examples the square pulses can range from 10 μs to 1000 μs. Any combination of the above can be used in different implementations, including dynamic implementations which selects values that are most accurately tuned to the implantable device and/or the patient.
In some implementations, the impedance of brain tissue is monitored over time a patient specific threshold or model is used to determine if a seizure is likely to occur. For example, a current impedance (e.g., from one reading a series of readings) can be compared to a longer term mean. An unexpected difference (e.g., a change not caused by change in brain state) is indicative that a seizure may occur. In some implementations, the threshold is a change in impedance value. In some of these implementations the threshold or model is further based on circadian rhythm. In one example, method brain state and impedance are monitored and changes in impedance that are anomalous to the current brain stare are indicative that a seizure is likely to occur. In some implementations, the threshold or model is further based on characteristics of the pulse used to stimulate the brain tissue for measuring tissue. For example, the pulse may use a square wave, sinusoidal wave, white noise simulation, etc. In some implementations, the threshold or model further consider a rate in change of impedance. In some examples, the velocity of the rate of change of impedance increases prior to a seizure. In some implementations, different time scales of impedance values are considered and compared. For example, impedance measurements collected over seconds, minutes, hours, days, weeks, months, years, etc. may considered and compared to detect patterns and features that are indicative of a user being likely to experience a seizure.
In some implementations, a patient specific classifier is built. In one example, an implantable device is implanted on a patient's brain tissue, where it continuously records impedance values that are used as training data for the patient specific classifier. These readings set baselines e.g., for expected sleep and wake state impedances as well as features observed before, after, and during a seizure. In some implementations, a patient self-records when a seizure occurred. In other implementations, impedance readings are processed to determine when a seizure may have occurred. In other implementations, further data is detected and processed to determine whether a seizure occurred. For example, other device may be used to measure brain activity. Various combinations of the data described above can be used to train the patient specific classifier. In some implementations, continuous learning techniques are used to further train the patient specific classifier or other classifier and/or deep learning models described herein.
In some implementations, the systems and method disclosed herein are used for spatial and temporal forecasting of seizures. In some implementations, the classifier disclosed herein predict when a seizure is going to occur and where in the brain the seizure is going to occur. In some implementations, the systems and methods disclosed herein are used to notify a user that the user has gone from a high risk of seizure state back to a low risk of seizure state. In some implementations, a patient in a high risk state is instructed to take a medication or another type of treatment may be applied. In other examples, the patient is instructed to put themselves in a safe location and position for experiencing a seizure.
In some implementations, systems and methods for electrical stimulation and linear response for tracking and modulating extracellular space and brain excitability are used.
Electrical stimulation and associated measurement of linear response (impedance) of tissue (e.g., brain tissue) is used to track and modulate the tissue extracellular space and excitability in a human brain. Impedance is a fundamental electrical property of brain tissue and important for understanding local field potentials and electrical stimulation. More specifically, brain impedance characterizes the resistance to ionic current flow through brain extracellular, intracellular, and vascular compartments.
In some implementations, brain impedance is used to measure the volume of the extracellular space; and then using the impedance measurements to guide manipulation of the size of the extracellular space therapeutically by stimulating brain tissue through the activation the glymphatic system based upon the brain impedance measurements, wherein the amount of brain impedance can be indicative of when to apply electrical stimulation to active the glymphatic system, which in turn sends fluid through the extracellular space to increase the size and volume of the extracellular space to wash away negative brain impurities. In some implementations, brain impedance measurements are used to determine the various levels of sleep stages a patient has undergone by measuring the glymphatic stages of the brain to determine how much deep sleep was done. In some examples, tying brain impedance measurements to sleep stages (and the lack thereof of sleep) and activating the glymphatic system to increase extracellular space at specified times, seizures can be predicted and prevented.
System and methods for classifying and modulating brain behavioral states are discussed in U.S. Pat. No. 11,647,962, entitled “SYSTEM AND METHOD FOR CLASSIFYING AND MODULATING BRAIN BEHAVIORAL STATES”, the disclosure of which is hereby incorporated by reference in its entirety. Materials and methods for using electrical stimulation to treat a mammal having a proteinopathy (e.g., neurodegenerative diseases) or at risk of developing a proteinopathy are discussed in U.S. Pat. No. 11,426,577, entitled “NEUROMODULATION TO MODULATE GLYMPHATIC CLEARANCE”, the disclosure of which is hereby incorporated by reference in its entirety. Therapy systems with quantified biomarker targeting, including for epilepsy treatment, and associated systems and methods are discussed in U.S. Patent Publication Number 20210346809, entitled “THERAPY SYSTEMS WITH QUANTIFIED BIOMARKER TARGETING, INCLUDING FOR EPILEPSY TREATMENT, AND ASSOCIATED SYSTEMS AND METHODS”, the disclosure of which is hereby incorporated by reference in its entirety.
The implanted device 102 is configured to capture impedance data 110 of brain tissue of the user U. The implanted device 102 comprises a plurality of electrodes that are surgically implanted in various locations with the user's brain. The implanted device 102 is configured to sample a brain's impedance by injecting current into the brain tissue and sensing the corresponding voltage signals. Further examples of implanted devices 102 are disclosed herein. In some implementations, the implanted device 102 is configured to continuously measure brain tissue. In some implementations, the brain tissue is densely sampled using the implanted device 102.
In some examples, the implantable device 102 is configured to interface with a cloud computing environment. For example, via the computing device 104 or other edge computing devices.
The computing device 104 comprises a user interface 108, and is configured to receive the impedance data 110 from the implanted device 102 and to present an alert via the user interface 108 based on a prediction that the user U (for example, a user associated with the computing device) is likely to experience a seizure, wherein the prediction is based on an impedance signature derived from the impedance data 110.
The seizure alert application 106 is an application configured to alert (e.g., via broadcasting visual and/or audio alarms) the user U when a seizure is predicted as likely to occur. In some implementations, the seizure alert application is configured to relay information from the implanted device to the cloud. In some implementations, the seizure alert application 106 is an application configured to receive a user input for indicating that the user is experiencing a seizure or recently experienced a seizure. In some of these implementations, the user input is used to determine when the user experience a seizure and impedance data from around the time the seizure occurred is collected and processed to determine an impedance signature that corresponds to the user being likely to experience a seizure.
In some implementations, the seizure alert application 106 runs a deep learning model to forecast whether a seizure is likely to occur at the computing device 104. For example, when the computing device 104 has sufficient computing capacity to run the deep learning model. In other examples, the seizure alert application 106 interfaces with a cloud computing environment to run the deep learning model.
In some implementations, the computing device 104 is a mobile computing device such as a tablet, smart phone, smart watch, etc. In some implementations, the computing device 104 is a specialized computing device designed and configured for seizure monitoring. In some implementations, the computing device 104 and the implanted device 102 are configured to allow bi-directional wireless streaming.
The operation 202 samples impedance measurements of brain tissue of a user to determine an impedance signature of the brain tissue of the user. In some implementations, sampling impedance measurements of brain tissue of a user is performed by, injecting an electrical current into the brain tissue with an implanted device, and sensing a voltage in the brain tissue with the implanted device in response to the applied electrical stimulation, and calculating impedance based on the electrical current and sensed voltage. In some implementations, the impedance signature includes an absolute value of observed impedance and a rate of change in impedance.
The operation 204 processes the impedance signature of the brain tissue to predict whether the user is likely to experience a seizure. In some of these implementations, an absolute value of observed impedance and/or a rate of change in impedance are compared to threshold values to predict whether a user is likely to experience a seizure. In some implementations, the thresholds are calibrated for a specific patient. In some implementations, a seizure is predicted as being likely to occur based on a detected a change in impedance of the brain tissue.
In some implementations, a deep learning model is used to predict whether the user is likely to experience a seizure. In some of these implementations, the primary architecture of the deep learning model is a bi-directional long-term memory model. In some implementations, training the machine learning model comprises aggregating model inputs by segmenting impedance data and determining rate of impedance change and time-frequency power signals for each segment, processing the model inputs with a parallel 1-D convolution network (CNN) to extract features, constructing feature representations to learn valid information through forward and/or backward propagation of the extracted features, applying an attention mechanism to redistribute weights of the feature representations to emphasize critical features of preictal stage; and train a full connection layer to output a preictal probability based on the feature representations with redistributed weights. In some examples, the preictal probability is compared to a threshold to determine whether a seizure is likely to occur. In some implementations, the deep learning model is trained using a cloud computing system.
In some implementations, the prediction of whether the user is likely to experience a seizure is further based on a detected brain state. In some of these implementations, the brain states include an awake state and one or more sleep states. Examples of sleep states include REM state, NREM state, etc. In some implementations, the method 200 further comprises measuring glymphatic stages of the brain tissue based on the impedance signature and determining levels of sleep stages for the patient based on the measured glymphatic stages.
The operation 206 provides an alert when it is predicted that the user is likely to experience the seizure. In some implementations, the alert is a visual and/or audio alter that is presented on a user computing device of the patient.
In some implementations, the method 200 further comprises applying a second electrical stimulation to manipulate an extra cellular space in the brain tissue in response to predicting that the seizure is likely to occur.
In some implementations, the operation 206 results in other actions, in addition to or instead of providing an alert. In some implementations, a treatment is applied or altered based on the prediction of whether a seizure will occur. Some implementations include providing an alert to patient alert the patient of a high risk period to allow the patient administers a prophylactic medication. In some implementations, an alert is provided to an implanted neurostimulator (e.g., the implanted device 102 shown in
In some implementations, the user is altered that they are likely to have a seizure within some period of time (i.e. 1 minute, 30 minutes, 3 hours). In some implementations, a user is alerted if they are at higher risk or lower risk of seizure that day than what is typical for that patient. In some implementations, a user is alterted to a higher or lower risk period in their infradian seizure cycle. Some implementations, include combinations thereof.
People with epilepsy (PWE) face the challenge of unpredictable seizures, significantly affecting their daily life quality. Even those resistant to medication spend less than 1% of their time experiencing a seizure, while the remainder is spent under constant seizure threat. PWE rely on daily anti-seizure medications that often have side effects. Seizure forecasting can transform epilepsy therapy and the lives of PWE by allowing preemptive administration of medications and electrical brain stimulation. A closed-loop system that integrates brain sensing, distributed computing, and neuromodulation therapy for seizure detection, intervention, and prevention is disclosed herein. As an integral component of the system, brain impedance is periodically sampled at multiple locations over a long-term period. Impedance has a close relationship with the brain's extracellular space (ECS), which plays a critical role in influencing neuronal excitability and seizures. There is a correlation between periodic changes in ECS, impedance oscillations, and behavioral state in human limbic system. A significant increase in impedance approximately 1-3 hours before a seizure is a useful biomarker indicative of increased seizure likelihoods. By integrating the impedance dynamics information into the system an automatic system for high performance seizure forecasting using deep-learning model and adaptive therapy with is implemented.
In some implementations, a deep-learning model for seizure forecasting is developed and applied. The deep learning model leverages impedance tracking and other biomarkers for automatic seizure prediction.
In some implementations, a clinically viable closed-loop neuromodulation system for seizure prediction, detection, and therapeutic intervention is developed. In some implementations, the system leverages distributed computing facilities, an epilepsy patient assist device. The impedance dynamics discussed herein are integrated to enhance seizure prediction and adaptive therapy.
In some implementations, the mechanisms underlying changes in extracellular space, impedance dynamics, and seizure generation are analyzed. Impedance measurements are applied to make seizure predictions.
Seizures occupy a relatively small fraction of time for most people with epilepsy (PWE), their apparent unpredictability is an overwhelming challenge of living with epilepsy. Even those with drug resistant epilepsy spend less than 1% of the time experiencing a seizure, with reminder of the time spent living under the threat of a potential seizure. Currently, PWE take anti-seizure medications daily that are often short-lived and infrequent. The potential for seizure forecasting to transform epilepsy therapy and improve the lives of PWE by enabling preemptive administration of medications and electrical brain stimulation.
A comprehensive approach encompassing biomarker detection, seizure forecasting/detection, intervention, and therapeutic stimulation/medicine delivery is disclosed. In some implementations, a system that integrates tracking brain activity, seizures and behavior in a closed-loop design is disclosed, enabling efficient delivery of epilepsy healthcare to patients in their natural home environment over an extended periods.
In some alternative examples, forecasting relies on biomarkers derived from local field potential or temporal structures of seizure occurrence.
In some examples, impedance as a novel biomarker in epilepsy is used for seizure forecasting. In some implementations, information about brain impedance temporal dynamics are leveraged. For example, fast impedance changes on millisecond timescales due to the opening of ion channels in active neuron membranes, and slow impedance changes occurring over seconds to minutes attributed to cell swelling, have been observed during seizures. There is a close relationship between impedance and the brain's extracellular space (ECS), which plays a role in influencing neuronal excitability and seizures. There is a strong correlation between periodic changes in ECS, impedance oscillations, and behavioral state in the human limbic system. A significant increase in impedance level and rate of impedance change several hours before a seizure is a biomarker indicative of high likelihood of seizure occurrence. As a component of the system, impedance is periodically sampled with high temporal resolution (about minute scale) at multiple brain locations over an extended period. In some implementations, a single brain region is sampled every 5-15 minutes. By integrating the impedance dynamics information into the system, an automatic and practical system for high performance seizure forecasting and adaptive therapy is developed.
In some implementations, the systems disclosed are used for seizure forecasting and prevention and advantages include improving the quality, safety, and effectiveness of seizure forecasting and intervention. Additionally, by incorporating information about impedance dynamics into the system, performance in seizure forecasting and adaptive therapy is improved, thereby allowing for the prevention of seizures before they occur. These forecasting models could use a machine learning method including deep learning, support vector machine (SVM) or random forest (RF). Example implementation 1 results.
In the example implementation, the following subjects and data was acquired. The pilot data came from a range of animal and human investigation into the temporal dynamics of brain impedance spectroscopy. Five human subjects diagnosed with drug-resistant bilateral mesial temporal lobe epilepsy were implanted with investigational Medtronic Summit RC+S™ devices targeting the bilateral anterior nucleus of the thalamus and bilateral mesial temporal structures. Each subject was implanted with four Platinum-Iridium (Pt—Ir) alloy leads, each containing four electrodes (contacts), resulting in 16 channels per subject. The leads were implanted in the left and right thalamus (THL), amygdala-hippocampus (AMG-HPC), and posterior hippocampus (post-HPC).
Electrical two-point monopolar impedance measurements were nonuniformly sampled every 5-15 minutes and transmitted to a cloud database through a wireless network. The effective impedance (Z) was calculated by injecting a single square-wave current pulse (I, 80 μs width), sensing the corresponding voltage (V), and applying the Ohm's low Z=V/I. This is equivalent to an impedance measured by injecting a 1 kHz sinusoidal current with an amplitude of 500 nA. Both the current injection and voltage sensing were delivered using the same electrodes.
In the example implementation, a cross-correlation estimation between seizures and impedance was made. The temporal dynamics of impedance around seizures were characterized by estimating the cross-correlation between seizure onset times and the measured impedance. To eliminate the slow impedance drift over an extended period, the impedance values were normalized by subtracting the daily impedance average. For each seizure, impedance values within a ±5-minute and a ±12-hour window surrounding the seizure onset were collected. The impedance measured across anatomical locations and across the five subjects by aligning at the onsets (
The significance of the observed impedance variation was assessed by calculating the baseline signal of impedance change. Time-matched interictal (days w/o seizures) surrogate data timestamps was generated according to the distribution of seizure onsets and repeated the same estimation procedure used for days with seizures. The baseline interictal data is from days without any seizures within ±12 hours (
In the example implementation, slow and very-slow seizure-related impedance changes were analyzed. Both the level and rate of change in impedance surrounding seizure onsets (
The results of very-slow impedance changes also suggest an increase in impedance level around seizure onsets (about 1 hour before and 4 hours after) within a ±12-hour window in hippocampus area (
In some implementations, an integrated seizure forecasting system with closed-loop design is disclosed. In some implementations, a seizure forecasting system utilizes impedance spectroscopy dynamics.
The implanted electrodes (A) are surgically implanted in various locations within the brain's limbic system. In some implementations, the implanted electrodes (A) are Pt—Ir electrodes. The brain's impedance is periodically sampled by injecting current into the brain tissue and sensing the corresponding voltage signals.
The implanted device (b) device manages the measurement of impedance among other signals and can optionally initiate preventive brain stimulation. In some examples, the system is agnostic to any implantable sensing and stimulation device with bi-directional wireless streaming for off-the-body analytics.
The electronic personal assistance device (EPAD) (C) is a handheld device (e.g., running on smartphone) that is configured to perform multiple functions, including relaying data from the implanted device to the cloud, broadcasting visual/audio alarms upon receiving a seizure risk alert, and receiving patient inputs. In some implementations, the EPAD (C) has sufficient edge computing capacity, and can run the deep-learning (DL) forecasting model directly.
The cloud co-processor (d) trains the deep learning (DL) model. The probability of the predicted seizure risk can also be computed in the cloud when the EPAD (C) does not have adequate computing bandwidth. In some implementations, an adaptation of Bi-directional Long Short-Term Memory (BiLS™) model for seizure forecasting is used. The original impedance data is segmented, and additional inputs such as the rate of impedance change, and time-frequency power signals of each segment are added to the model input. The input is then fed into a parallel 1-D convolution network (CNN) for feature extraction. The subsequent full connection layer constructs feature representation, which goes to the next BiLS™ layer to learn valid information from the noisy data through forward/backward propagation. An attention mechanism is then employed to redistribute the weights and emphasize the critical features of preictal stage from a long sequence of signals. Finally, a full connection layer outputs the preictal probability. The specific threshold for classification will be determined by forecasting performance requirements.
The monitoring station (e) receives estimated seizure risk and other related information from the cloud co-processor (d) in real-time (or near-real time), providing critical patient information. In some implementations, the monitoring station (e) is located at a designated medical center, where the station receives estimated seizure risk and other related information from the cloud co-processor (d) in real-time, providing critical patient information to the medical/research team. In some implementations, the monitoring station (e) may automatically send seizure alarm signals to EPAD (C).
In some implementations, a Medical/Research Team (F) reviews the automatic seizure alert signals and makes the final decision regarding sending warning signals to patients.
Impedance is an electrical property of brain tissue and plays role in shaping the characteristics of local field potentials, the extent of ephaptic coupling, and the volume of tissue activated by externally applied electrical brain stimulation. Brain impedance, sleep-wake behavioral state, and epileptiform activity were tracked in five people with epilepsy living in their natural environment using an investigational device. The example implementation identified impedance oscillations that span hours to weeks in the amygdala, hippocampus, and anterior nucleus thalamus. The impedance in these limbic brain regions exhibit multiscale cycles with ultradian (˜1.5-1.7 h), circadian (˜21.6-26.4 h), and infradian (˜20-33 d) periods. The ultradian and circadian period cycles are driven by sleep-wake state transitions between wakefulness, nonrapid eye movement (NREM) sleep, and rapid eye movement (REM) sleep. Limbic brain tissue impedance reaches a minimum value in NREM sleep, intermediate values in REM sleep, and rises through the day during wakefulness, reaching a maximum in the early evening before sleep onset. Infradian (˜20-33 d) impedance cycles were not associated with a distinct behavioral correlate. Brain tissue impedance is known to strongly depend on the extracellular space (ECS) volume, and the findings reported here are consistent with sleep-wake-dependent ECS volume changes recently observed in the rodent cortex related to the brain glymphatic system. In this example implementation, human limbic brain ECS changes during sleep-wake state transitions underlie the observed multiscale impedance cycles. Impedance is an electrophysiological biomarker that this implementation found useful for tracking ECS dynamics in human health, disease, and therapy.
Electrical impedance is an extrinsic tissue property that characterizes the frequency-dependent resistance to endogenously generated ionic currents and externally applied electrical stimulation currents in the brain. The impedance influences the spatiotemporal dynamics of extracellular ionic currents that give rise to local field potentials (LFPs) and determines the range of ephaptic coupling and the volume of tissue activated by electrical brain stimulation (EBS). The voltage generated in response to an externally applied current can be used to probe and track brain impedance, and at the scale of clinical EBS electrodes (˜mm2), the impedance depends on the electrode-tissue interface and ionic current conductance through the intracellular, vascular, and extracellular space (ECS) brain compartments. The impedance of rodent, nonhuman primate, and human brain tissue is primarily resistive over the frequency range of physiologic LFP and is strongly dependent on the ECS volume fraction. The interstitial conductive fluid filled ECS creates a complex electrical network of ionic current flow that is increasingly recognized for its role in sleep-dependent brain health, such as clearing neurotoxic metabolites and pathologic proteins associated with neurodegenerative disorders.
In this example implementation, it was hypothesized that if sleep-wake transitions drive changes in human brain ECS volume, and manifest as impedance cycles with lower impedance in NREM sleep, given the increased ECS volume compared with wakefulness and REM sleep. In this example implementation, an investigational implantable neural sensing and stimulation device was used to continuously track impedance, sleep-wake behavioral state, and epileptiform activity in anterior nucleus of thalamus (ANT), amygdala (AMG), and hippocampus (HPC) over several months in five people with drug-resistant mesial temporal lobe epilepsy living in their natural home environment (
This example implementation people with drug-resistant mesial temporal lobe epilepsy were investigated. Seven subjects with drug resistant mesial temporal lobe epilepsy were identified for the study. Five subjects ultimately met the full inclusion criteria with adequate numbers of baseline seizures and were implanted with the investigational neural sensing and stimulation device. The demographics and clinical information of the five subjects implanted with the Medtronic RC+S are provided in the table shown in
The epilepsy patient assistant (EPA) is a custom software application running on a handheld computer enabling bidirectional communication between the implanted device, wearables, and cloud computing resources. The EPA features include automated algorithms for continuous LFP data acquisition, electrical stimulation, impedance testing, LFP analysis and an interface for collecting patient interactions.
The investigational Medtronic RC+S is an investigational implantable device with bidirectional wireless communication capability, programmable 16-channel electrical stimulation, two-point impedance measurements, and continuous four-channel (selected bipolar pairs) wireless LFP streaming to a handheld computer and cloud environment
For lead and electrode contact localization (
In some implementations, the LFP is sensed using the RC+S device. The LFPs from left and right ANT, AMG, and HPC using four bipolar electrode pairs, sampled at a frequency of 250-500 Hz, were selected for continuous data streaming. The contacts used to create lead-specific bipolar recordings were selected by visual review of LFP data during seizures, resting wakefulness, and sleep-wake transitions. In some implementations, the ANT electrode contacts were not used for LFP recordings or further analysis.
The sensed LFP is analyzed. The LFP analysis was performed retrospectively. Previously validated automated algorithms were applied to the long-term intracranial LFP recordings to identify seizures, IES and classify wake-sleep state (Awake, REM, and NREM). The algorithm pipeline identifies seizures, IES and sleep-wake behavioral state for consecutive 30 s data segments. The sleep-wake classifications and seizure events were then synchronized with the impedance measurements for further analysis. Toolboxes and data for automated LFP analysis for detecting IES, seizures, and performing sleep-wake classification were developed and used.
Seizures and IES are detected from continuous LFP recordings. Interictal epileptiform spikes and seizures are electrographic biomarkers of pathologic, epileptogenic brain tissue and are readily identified by human visual review of LFP recordings. Continuous hippocampal LFP recordings from the five subjects were visually reviewed and IES and seizures labeled for subsequent use in training, validating, and testing automated detectors.
A validated algorithm for detecting IES transients from LFP recordings was used. The adaptive IES algorithm enables detecting IES in long-term LFP data with changing background activity commonly encountered in prolonged recordings spanning weeks of time. The training data was used to set a hypersensitive threshold for all subjects (
The intracranial LFP associated with epileptic seizures exhibits characteristic temporal and spectral evolution over a wide frequency range (
The temporal distribution of seizures was determined by plotting the circular histogram of seizure onset times for all verified seizures across all subjects. To compare the impedance profiles during seizures with those observed outside seizure events, surrogate nonseizure data were generated. These surrogate nonseizure data were carefully aligned to match the time-of-day distribution observed for actual seizures. The generation of surrogate nonseizure data involved simulating 24 h periods without any seizure activity and ensuring equivalent distributions of times for impedance comparison. Importantly, the surrogate data account for the natural circadian impedance cycles. To create the estimated probability distribution of seizures, a rejection method was used based on the geometric interpretation of probability distribution. This method allowed for modeling the probability distribution of seizure occurrences and evaluate the impedance patterns specific to seizure events against the baseline nonseizure impedance data.
In some examples, sleep-wake classification is made from continuous LFP recordings. Reliable sleep-wake classification was determined by training a subject-specific classifier on simultaneous LFP and gold standard sleep annotations from polysomnography (PSG). For each subject over the course of three nights in the hospital epilepsy monitoring unit, we recorded simultaneous scalp PSG and continuous intracranial LFP data. The PSG signals were scored into standard sleep categories (Awake and REM) and two NREM stages (N2, N3) using American Academy of Sleep Medicine 2012 scoring rules. A subject-specific sleep classifier (Naïve Bayes) was trained using the first night of data. The second night was used for validation and the third night for pseudoprospective testing. The Naïve Bayes classifier uses features extracted from the LFP that were determined to be useful for sleep-wake behavioral state classification.
Electrical impedance is measured and analyzed. The two-point monopolar electrical impedance was measured from all electrodes in AMG, HPC, and ANT using a square-wave current pulse (0.4 mA, 80 μs pulse width) over multiple months in five people with drug-resistant epilepsy. In some implementations, the AMG-HPC electrode contacts in seizure onset tissue and ANT electrode contacts used for therapeutic EBS were excluded from impedance analysis.
The RC+S two-electrode impedance method uses the same electrode contacts for delivering current stimulation and sensing the voltage response. The voltage response to the applied current pulse using the two-point method includes the electrode-tissue interface polarization, created by the electrical stimulation, in addition to the bulk tissue impedance. To address the electrode tissue interface contribution caused by electrolyte polarization, the four impedance calculation can be used. This method involves using one pair of electrodes to deliver the current and a separate pair of electrodes to measure the voltage response (
A composite model was used to compare the impedance measurements obtained with the two-point square-wave pulse current probe the RC+S device uses with the results obtained from the four-electrode and two-electrode measurement using sinusoidal stimulation currents.
In some implementations, saline/microbead impedance is measured. Benchtop experiments were conducted using saline and saline/microbead composites (
The two-point RC+S impedance measurement uses a single current pulse (0.4 mA, 80 μs pulse-width) to probe the tissue and calculates the impedance from Ohm's law Zeff=Vmeasured/Istim with the voltage measured at 70 μs, near the end of the 80 μs current test pulse (
The measured human brain impedance was analyzed. The impedance time series (sampled every 5-15 min) in the five subjects were resampled with a 30-min-long moving average window with a step of 10 min and bandpass filtered between 30 min to 100 d with low- and highpass finite impulse response, zero phase shift filters of 1,001st order.
Polar impedance plots with the histograms of the minimum and maximum impedance value recorded over each 24 h (day/night) cycle were tested for nonuniform circular distributions using the Kuiper (1960) test (*p, 0.05, **p, 0.01, ***p, 0.001).
The multiscale impedance cycle was analyzed. The continuous wavelet transform (CWT) was used to investigate brain impedance cycles with ultradian (<24 h), circadian (˜24 h), and infradian (>24 h) periodicities. The CWTs were implemented in MATLAB software (Morlet wavelets and L2 normalization). The Thomson F test multitaper scheme was used to test for significance of oscillations with a 5-day moving window used for ultradian and circadian cycles and a 100-day moving window (step of 5 days) for infradian cycles. The presence of ultradian, circadian, and infradian frequency band oscillations were identified using preprocessed impedance signals. The periodicity test in ultradian band for each 5 d signal segment was investigated using the time-half-band product (TW) set at three (TW=3), and the number of tapers (K) was five (K=5), resulting in the frequency resolution of 2 W=1.2 (or W=0.6) cycles/day. In ultradian band, frequency of interest (FOI) was chosen from 23/24=0.96 h/cycle to 6 h/cycle with an incremental step of 15 min; in infradian band, the FOI was from 2.5 d/cycle to 30 d/cycle with an incremental step of 15 min. At each fixed frequency, the F test was set at the level of the upper 99% (p=0.01) percentage point of the F distribution with 2 and 2 K−2=8 degrees of freedom under the null hypothesis of no spectral line, leading to the critical value of 8.65.
In some examples, the sleep-wake state dependence of impedance was analyzed. The three behavioral categories on a group and subject level were compared using a statistical two-tailed Mann-Whitney test (p<0.05, Bonferroni correction for multiple observations). The continuous behavioral state classifications were also reviewed for periods of daytime sleep (>30 min duration), and in four of five subjects' daytime naps were identified to analyze and comparison of sleep impedance during daytime naps and overnight sleep.
The LFP and impedance data was summarized. The RC+S provides streaming of four LFP channels (250-500 Hz sampling) and 16 monopolar impedance values (0.001-0.003 Hz sampling) to the patient's tablet computer for storage. When there is either cellular or wireless connectivity available the data are transferred from patient's tablet to a cloud database. The sampling frequencies were determined as a compromise between data quality, wireless data transmission, and device battery limitations. The LFP and impedance data were collected during an ongoing investigational device clinical trial that is testing different ANT EBS paradigms for drug-resistant epilepsy that include (1) no EBS baseline, (2) low-frequency EBS (2 and 7 Hz continuous), and (3) high-frequency EBS (125 Hz continuous, duty cycle, 1 min on and 5 min off, and responsive).
In one example, because of the acute effects on tissue impedance, the following was excluded (1) data from the initial 14 d after surgical implant, (2) data within 12 h before and after seizures, (3) data from electrode contacts located in seizure onset zone, (4) data from electrode contacts used for EBS, and (5) data from electrode contacts with high impedance (0.5000 V). This left 45.00% (36/80) of the electrode contacts (5 AMG, 14 HPC, 17 ANT) for subsequent analysis of impedance and association with sleep-wake behavioral state.
The complete dataset available for analysis included at least 150 d for each patient (M1, 828 d; M2, 549 d; M3, 346; M4, 549; M5, 150 d). In some examples, days with poor telemetry and, 70% data were excluded because of the possibility of missing seizures, leaving 1456 continuous 24 h epochs from 5 subjects (291±198 24 h epochs). The analysis of sleep-wake impedance changes was focused on the 30 d period before starting EBS. During the initial period without EBS there were 299 data segments (30-90 min durations) from the 36 electrode contacts without seizure activity for at least 24 h annotated across the five patients and three behavioral categories—Awake, NREM, REM (M1, 30/21/22; M2, 20/19/18; M3, 15/14/13; M4, 22/27/29; M5, 12/19/18) for analysis.
The LFP recordings from the five subjects were used as input to automated algorithms for detecting seizures and classifying sleep-wake state (
Saline/microbead composites were used to compare two-point RC+S impedance measurements that use a square-wave current pulse with both two-point and four-point impedance measurements using sinusoidal currents (
Impedance (e.g., human brain impedance) was measured in five people with drug-resistant mesial temporal lobe epilepsy over multiple months using the RC+S monopolar two-point current pulse method with periodic sampling (every 5-15 min) in AMG, ANT, and HPC. The impedance increases after surgery but reaches a stable value after ˜14 d (
Seizures effect brain impedance and ECS volume, and to address this possible confound we used a hypersensitive seizure detector combined with expert visual review to label all seizures in the LFP data. The seizures primarily occur in either the morning or late afternoon in all subjects M1-5 (
Group-level analysis (
In some implementations, multiscale impedance oscillations are analyzed. In addition to ˜24 h cycles captured with the polar plots (
In some implementations, correlations between impedance and behavioral states (wakefulness, NREM, and REM sleep) were analyzed. In total, 299 data segments from the 36 right-hemisphere electrode contacts without seizure activity for at least 24 h were identified across the five patients and three behavioral states. The impedance was lower in NREM sleep compared with wakefulness and REM sleep (
Four of the five subjects also took afternoon naps (nap onset 2:06 P.M.±1.2 h). Impedance in NREM sleep during naps was analyzed (M1, 7; M2, 1; M3, 7; M4, 3) and compared with the wakefulness impedance in the hour before the nap. The nap NREM sleep impedance is decreased when compared with wakefulness in HPC, AMG, and ANT (
The average baseline total impedance in the AMG, HPC, and ANT during wakefulness was 753.58±80.23 Ω, 892.74±97.27Ω, and 1268.28±63.17Ω, and during NREM sleep was reduced to 730.11±77.00 Ω, 874.30±99.49Ω, and 1246.61±54.35Ω, respectively (
The total impedance consists of the electrode-tissue impedance in series with the brain tissue impedance (
The human AMG, ANT, and HPC show complex temporal dynamics with multiscale oscillations spanning hours to weeks (
Sleep-wake state transitions appear to generate the observed impedance cycles with characteristic oscillations anchored by a highly stable, dominant 24 h cycle. Compared with the circadian cycles, the ultradian and infradian cycles show more variability (
In some examples, ECS volume changes drive the ultradian and circadian impedance cycles in human brain given the strong dependence of electrical impedance on ECS volume and recent demonstration of ECS volume expansion in rodent cortex during NREM sleep. The impedance data are indirect evidence that sleep-wake state transitions (wakefulness, NREM, and REM sleep) differentially modulate human ECS and is consistent with the ECS dynamics proposed for a human brain glymphatic system. The vascular and intracellular compartment conductances will have limited contributions to composite impedance with the current probe used here because of the high impedance of the blood-brain barrier and neuronal and glial cell membranes. The electrode-tissue interface contribution to the impedance is not expected to depend on the sleep-wake state and should not differentially contribute to impedance changes.
The interstitial-fluid-filled ECS at the physiological volume fraction (˜0.2) and complex cellular geometry provides a highly conductive network that may be near the percolation threshold enabling significant impedance changes with relatively small ECS volume changes.
It was estimated that the ECS volume change, ΔVECS, that would produce the observed circadian impedance change ΔZ. The impedance change is related to the change in effective specific brain conductivity ΔZ=1/4παΔσσeff where α is the cathode/anode stimulation electrode separation. The effective specific conductivity of a composite media of cells embedded in the ECS volume is Δσeff=(3/2) σisf ΔVecs, where σisf is the conductivity of the interstitial fluid. Usings σisf≅1:79 S/m for the conductivity of the ECS interstitial fluid and adjusting the ECS percolation threshold volume fraction to 0.07, a brain-specific conductivity of ˜0.250 S/m yields reasonable agreement with reported ECS volume fraction of ˜0.16.
In the real-time iontophoretic tetramethylammonium diffusion studies in mice the ECS volume in wakefulness was 0.136 and expanded to 0.227 in the anesthetic state, a AVECS≈0:09 which is a 67% relative volume change. In the composite conductivity model above this corresponds to a change in specific conductivity of Δσeff=0.24 S/m. A measured ΔZ≈30Ω in human AMG and HPC corresponds to a specific conductivity change of Δσeff=0.29 S/m from wakefulness to NREM sleep and a corresponding ΔVECS=0.10: Based on this model, the human ECS volume change associated with transition to NREM sleep is slightly larger than what is observed in mice.
In some examples, ECS volume dynamics underlie the observed sleep-wake state dependence and multiscale impedance cycles. Alternatively, sleep state-dependent change in astrocyte membrane impedance combined with a gap-junction-mediated astrocytic syncytium network could provide a low impedance pathway. Additionally, changes in ionic conductivity because of ion concentration changes in the interstitial fluid could have an effect. Given the previous two-photon imaging and real-time iontophoretic tetramethylammonium diffusion findings in mouse cortex during NREM and anesthesia, the observed ECS volume changes would have a large effect on impedance.
In some examples, the data segments selected were without seizures and brain electrodes with low IES rates were and not involved in seizure onset there is electrophysiological evidence that the brain regions analyzed are not entirely normal. There are IESs throughout the limbic network in all subjects (
Fast impedance changes, on millisecond timescales due to opening of ion channels in active neuron membrane, and slow impedance change, occurring over seconds to minutes attributed to cell swelling, have been observed during seizures. Very slow impedance changes occurring on timescales of minutes to hours are observed in the example implementation 3. This knowledge provides insights into the mechanism of seizure generation and be useful for seizure forecasting. In this example, seizure-related impedance change over 24 hours in humans with TLE were studied.
Five peoples with drug resistant TLE were implanted with bilateral electrodes in the thalamus, amygdala, and hippocampus. One patient was excluded from this analysis due to a prior epilepsy surgery resulting in inconsistent seizure-related impedance change. The impedance was sampled nonuniformly (every 5-15 minutes) over multiple months. Seizures were identified from both hemispheres, but mainly from the more epileptogenic left side. The distributions of seizure onset times were quantified with circular histograms. Surrogate data segments from days without seizures were generated and aligned to the same distributions of seizure onset as the data segments with seizures. Cross-correlation between impedance and seizures was estimated and compared to cross-correlation with surrogate seizures.
The estimated cross-correlation showed a peak in correlation in the amygdala-hippocampus and posterior hippocampus. In left thalamus, the mean value of normalized impedance related to real seizures was lower than that related to surrogate seizures about 0-10 hours before seizure onset. In left amygdala-hippocampus, the mean value of normalized impedance related to real seizures was higher than that related to surrogate ones about 0-10 hours before seizure onset. Significant difference was found in the left side of posterior hippocampus. The normalized impedance related to real seizures was significantly higher (mean values were separated by >2 SEM) than that related to the surrogate seizures ˜0-4 hours before seizure onset.
The peaks of cross-correlated impedance may suggest that seizure onsets could be phase locked to the circadian cycle of impedance. The impedance of thalamus is likely lowered, while that of hippocampus elevated hours before the seizure onset. The significant increase of impedance 0-4 hours before seizure onset in the left posterior hippocampus is useful for forecasting seizure alert.
In some implementations, the impedance describes the frequency dependent linear response voltage (0.001-10,000 Hz) generated by a test probe electric current. The impedance can be used a biomarker for cellular, extracellular matrix, and interstitial extracellular electrolyte fluid electrical interactions. The impedance frequency spectrum can be used to track behavioral state (sleep-wake), the glymphatic system function, and probability of seizure occurrence.
In some implementations, the temporal dynamics of impedance frequency dispersion (0.001-10,000 Hz) in different anatomical brain regions and their correlations are biomarkers are used for tracking behavioral state (sleep-wake), the glymphatic system function, and probability of seizure occurrence.
The computing device 1700 includes a processor 1702, a memory 1704, a storage device 1706, a high-speed interface 1708 connecting to the memory 1704 and multiple high-speed expansion ports 1710, and a low-speed interface 1712 connecting to a low-speed expansion port 1714 and the storage device 1706. Each of the processor 1702, the memory 1704, the storage device 1706, the high-speed interface 1708, the high-speed expansion ports 1710, and the low-speed interface 1712, are interconnected using various busses, and may be mounted on a common motherboard or in other manners as appropriate. The processor 1702 can process instructions for execution within the computing device 1700, including instructions stored in the memory 1704 or on the storage device 1706 to display graphical information for a GUI on an external input/output device, such as a display 1716 coupled to the high-speed interface 1708. In other implementations, multiple processors and/or multiple buses may be used, as appropriate, along with multiple memories and types of memory. Also, multiple computing devices may be connected, with each device providing portions of the necessary operations (for example, as a server bank, a group of blade servers, or a multi-processor system).
The memory 1704 stores information within the computing device 1700. In some implementations, the memory 1704 is a volatile memory unit or units. In some implementations, the memory 1704 is a non-volatile memory unit or units. The memory 1704 may also be another form of computer-readable medium, such as a magnetic or optical disk.
The storage device 1706 is capable of providing mass storage for the computing device 1700. In some implementations, the storage device 1706 may be or contain a computer-readable medium, such as a floppy disk device, a hard disk device, an optical disk device, or a tape device, a flash memory or other similar solid state memory device, or an array of devices, including devices in a storage area network or other configurations. The computer program product may also contain instructions that, when executed, perform one or more methods, such as those described above. The computer program product can also be tangibly embodied in a computer- or machine-readable medium, such as the memory 1704, the storage device 1706, or memory on the processor 1702.
The high-speed interface 1708 manages bandwidth-intensive operations for the computing device 1700, while the low-speed interface 1712 manages lower bandwidth-intensive operations. Such allocation of functions is exemplary only. In some implementations, the high-speed interface 1708 is coupled to the memory 1704, the display 1716 (for example, through a graphics processor or accelerator), and to the high-speed expansion ports 1710, which may accept various expansion cards (not shown). In the implementation, the low-speed interface 1712 is coupled to the storage device 1706 and the low-speed expansion port 1714. The low-speed expansion port 1714, which may include various communication ports (for example, USB, Bluetooth, Ethernet, wireless Ethernet) may be coupled to one or more input/output devices, such as a keyboard, a pointing device, a scanner, or a networking device such as a switch or router, for example, through a network adapter.
The computing device 1700 may be implemented in a number of different forms, as shown in the figure. For example, it may be implemented as a standard server 1720, or multiple times in a group of such servers. In addition, it may be implemented in a personal computer such as a laptop computer 1722. It may also be implemented as part of a rack server system 1724. Alternatively, components from the computing device 1700 may be combined with other components in a mobile device (not shown), such as a mobile computing device 1750. Each of such devices may contain one or more of the computing device 1700 and the mobile computing device 1750, and an entire system may be made up of multiple computing devices communicating with each other.
The mobile computing device 1750 includes a processor 1752, a memory 1764, an input/output device such as a display 1754, a communication interface 1766, and a transceiver 1768, among other components. The mobile computing device 1750 may also be provided with a storage device, such as a micro-drive or other device, to provide additional storage. Each of the processor 1752, the memory 1764, the display 1754, the communication interface 1766, and the transceiver 1768, are interconnected using various buses, and several of the components may be mounted on a common motherboard or in other manners as appropriate.
The processor 1752 can execute instructions within the mobile computing device 1750, including instructions stored in the memory 1764. The processor 1752 may be implemented as a chipset of chips that include separate and multiple analog and digital processors. The processor 1752 may provide, for example, for coordination of the other components of the mobile computing device 1750, such as control of user interfaces, applications run by the mobile computing device 1750, and wireless communication by the mobile computing device 1750.
The processor 1752 may communicate with a user through a control interface 1758 and a display interface 1756 coupled to the display 1754. The display 1754 may be, for example, a TFT (Thin-Film-Transistor Liquid Crystal Display) display or an OLED (Organic Light Emitting Diode) display, or other appropriate display technology. The display interface 1756 may comprise appropriate circuitry for driving the display 1754 to present graphical and other information to a user. The control interface 1758 may receive commands from a user and convert them for submission to the processor 1752. In addition, an external interface 1762 may provide communication with the processor 1752, so as to enable near area communication of the mobile computing device 1750 with other devices. The external interface 1762 may provide, for example, for wired communication in some implementations, or for wireless communication in other implementations, and multiple interfaces may also be used.
The memory 1764 stores information within the mobile computing device 1750. The memory 1764 can be implemented as one or more of a computer-readable medium or media, a volatile memory unit or units, or a non-volatile memory unit or units. An expansion memory 1774 may also be provided and connected to the mobile computing device 1750 through an expansion interface 1772, which may include, for example, a SIMM (Single In Line Memory Module) card interface. The expansion memory 1774 may provide extra storage space for the mobile computing device 1750, or may also store applications or other information for the mobile computing device 1750. Specifically, the expansion memory 1774 may include instructions to carry out or supplement the processes described above, and may include secure information also. Thus, for example, the expansion memory 1774 may be provide as a security module for the mobile computing device 1750, and may be programmed with instructions that permit secure use of the mobile computing device 1750. In addition, secure applications may be provided via the SIMM cards, along with additional information, such as placing identifying information on the SIMM card in a non-hackable manner.
The memory may include, for example, flash memory and/or NVRAM memory (non-volatile random access memory), as discussed below. The computer program product contains instructions that, when executed, perform one or more methods, such as those described above. The computer program product can be a computer- or machine-readable medium, such as the memory 1764, the expansion memory 1774, or memory on the processor 1752. In some implementations, the computer program product can be received in a propagated signal, for example, over the transceiver 1768 or the external interface 1762.
The mobile computing device 1750 may communicate wirelessly through the communication interface 1766, which may include digital signal processing circuitry where necessary. The communication interface 1766 may provide for communications under various modes or protocols, such as GSM voice calls (Global System for Mobile communications), SMS (Short Message Service), EMS (Enhanced Messaging Service), or MMS messaging (Multimedia Messaging Service), CDMA (code division multiple access), TDMA (time division multiple access), PDC (Personal Digital Cellular), WCDMA (Wideband Code Division Multiple Access), CDMA2000, or GPRS (General Packet Radio Service), among others. Such communication may occur, for example, through the transceiver 1768 using a radio-frequency. In addition, short-range communication may occur, such as using a Bluetooth, WiFi, or other such transceiver (not shown). In addition, a GPS (Global Positioning System) receiver module 1770 may provide additional navigation- and location-related wireless data to the mobile computing device 1750, which may be used as appropriate by applications running on the mobile computing device 1750.
The mobile computing device 1750 may also communicate audibly using an audio codec 1760, which may receive spoken information from a user and convert it to usable digital information. The audio codec 1760 may likewise generate audible sound for a user, such as through a speaker, for example, in a handset of the mobile computing device 1750. Such sound may include sound from voice telephone calls, may include recorded sound (e.g., voice messages, music files, etc.) and may also include sound generated by applications operating on the mobile computing device 1750.
The mobile computing device 1750 may be implemented in a number of different forms, as shown in the figure. For example, it may be implemented as a cellular telephone 1780. It may also be implemented as part of a smart-phone 1782, personal digital assistant, or other similar mobile device.
Various implementations of the systems and techniques described here can be realized in digital electronic circuitry, integrated circuitry, specially designed ASICs (application specific integrated circuits), computer hardware, firmware, software, and/or combinations thereof. These various implementations can include implementation in one or more computer programs that are executable and/or interpretable on a programmable system including at least one programmable processor, which may be special or general purpose, coupled to receive data and instructions from, and to transmit data and instructions to, a storage system, at least one input device, and at least one output device.
These computer programs (also known as programs, software, software applications or code) include machine instructions for a programmable processor, and can be implemented in a high-level procedural and/or object-oriented programming language, and/or in assembly/machine language. As used herein, the terms machine-readable medium and computer-readable medium refer to any computer program product, apparatus and/or device (e.g., magnetic discs, optical disks, memory, Programmable Logic Devices (PLDs)) used to provide machine instructions and/or data to a programmable processor, including a machine-readable medium that receives machine instructions as a machine-readable signal. The term machine-readable signal refers to any signal used to provide machine instructions and/or data to a programmable processor.
To provide for interaction with a user, the systems and techniques described here can be implemented on a computer having a display device (e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor) for displaying information to the user and a keyboard and a pointing device (e.g., a mouse or a trackball) by which the user can provide input to the computer. Other kinds of devices can be used to provide for interaction with a user as well; for example, feedback provided to the user can be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user can be received in any form, including acoustic, speech, or tactile input.
The systems and techniques described here can be implemented in a computing system that includes a back end component (e.g., as a data server), or that includes a middleware component (e.g., an application server), or that includes a front end component (e.g., a client computer having a graphical user interface or a Web browser through which a user can interact with an implementation of the systems and techniques described here), or any combination of such back end, middleware, or front end components. The components of the system can be interconnected by any form or medium of digital data communication (e.g., a communication network). Examples of communication networks include a local area network (LAN), a wide area network (WAN), and the Internet.
The computing system can include clients and servers. A client and server are generally remote from each other and typically interact through a communication network. The relationship of client and server arises by virtue of computer programs running on the respective computers and having a client-server relationship to each other.
Some of the examples described herein include or are defined by the following implementations.
Implementation A1 is a method for predicting and managing seizures in an individual, the method comprising obtaining at least one measurement of an impedance of brain tissue of the individual to determine an impedance signature of the brain tissue of the individual, processing the impedance signature of the brain tissue to predict whether the individual is likely to experience a seizure, and in response to predicting based on the impedance signature that the individual is likely to experience a seizure, initiating a remedial action that comprises at least one of alerting the individual or another user that the individual is likely to experience the seizure, applying an electrical stimulus to the brain tissue of the individual to reduce a likelihood of the seizure occurring or to reduce a severity of the seizure, or delivering a drug to the brain tissue of the individual to reduce a likelihood of the seizure occurring or to reduce a severity of the seizure.
Implementation A2 is the method of implementation A1, wherein sampling the impedance of the brain tissue of an individual includes: electrically stimulating the brain tissue with an implanted device that injects an electrical current into the brain tissue, and sensing a voltage in the brain tissue with the implanted device in response electrically stimulating the brain tissue, and calculating the impedance based on the electrical current and the sensed voltage.
Implementation A3 is the method of any of implementations A1 and A2, wherein a machine learning model is used to predict whether the individual is likely to experience a seizure.
Implementation A4 is the method of implementation A3, wherein the machine learning model is an unsupervised model such as a deep learning model
Implementation A5 is the method of any of implementations A3 and A4, wherein the machine learning algorithm is a support vector machine (SVM) or random forrest (RF).
Implementation A6 is the method of any of implementations A3-A5, wherein the machine learning model outputs the probability of an upcoming seizure, wherein the remedial action is initiated when the probability of the upcoming seizure exceeds a threshold.
Implementation A7 is the method of any of implementations A3-A6, wherein the machine learning model is a deep learning model.
Implementation A8 is the method of implementation A7, wherein the deep learning model is a bi-directional long-term memory model.
Implementation A9 is the method of any of implementations A7 and A8, wherein training the deep learning model comprises aggregating model inputs by segmenting impedance data and determining rate of impedance change and time-frequency power signals for each segment, processing the model inputs with a parallel 1-D convolution network (CNN) to extract features, constructing feature representations to learn valid information through forward and/or backward propagation of the extracted features, applying an attention mechanism to redistribute weights of the feature representations to emphasize critical features of preictal stage, and training a full connection layer to output a preictal probability based on the feature representations with redistributed weights.
Implementation A10 is the method of implementation A9, wherein the preictal probability is compared to a threshold to determine whether a seizure is likely to occur.
Implementation A11 is the method of any of implementations A7-A10, wherein the deep learning model is trained using a cloud computing system.
Implementation A12 is the method of any of implementations A1-A11, wherein processing the impedance signature of the brain tissue to predict whether the individual is likely to experience a seizure is further based on a detected brain state.
Implementation A13 is the method of implementation A12, wherein brain states include an awake state and one or more sleep states.
Implementation A14 is the method of any implementations A1-A13, the method further comprising measuring glymphatic stages of the brain tissue based on the impedance signature, and determining levels of sleep stages for the individual based on the measured glymphatic stages.
Implementation A15 is the method of any of implementations A1-A14, wherein the impedance signature includes an absolute value of observed impedance and a rate of change in impedance.
Implementation A16 is the method of any of implementations A1-A15, wherein prior altering the individual, sending the impedance signature to a monitoring station for review by a provider individual.
Implementation A17 is the method of any of implementations A1-16, wherein the individual is a mammal.
Implementation A18 is the method of any of implementations A1-A17, wherein the individual is a human.
Implementation A19 is the method of any of implementations A1-A18, wherein the impedance signature is based on measurements of impedance at a single location in a brain of the individual.
Implementation A20 is the method of any of implementations A1-A19, wherein the impedance signature is based on measurements of impedance at multiple locations in a brain of the individual.
Implementation A21 is the method of any of implementations A1-A20, wherein the impedance signature is based on measurements of impedance spanning a time interval of at least 1, 2, 3, 4, or 5 hours.
Implementation A22 is the method of any of implementations A1-A21, wherein the impedance signature is based on measurements of impedance spanning a time interval of at least 0.5, 1, 2, 3, 4, or 5 days.
Implementation A23 is the method of any implementations A1-A22, wherein the impedance signature is based on measurements of impedance spanning a time interval of at least 20-33 days.
Implementation A24 is the method of any of implementations A1-A23, wherein processing the impedance signature to predict whether the individual is likely to experience a seizure comprises comparing the impedance signature to a model impedance signature.
Implementation A25 is the method of implementation A24, wherein the model impedance signature is personalized to the individual based on recordings from the individual.
Implementation A26 is the method of any implementations A24 and A25, wherein the model impedance signature comprises a threshold impedance level, and processing the impedance signature comprises determining that at least one measurement of the impedance of the individual exceeds the threshold impedance level.
Implementation A27 is the method of implementation A26, wherein the threshold impedance level is adjusted based on a detected brain state of the individual, wherein the detected brain state is selected from an awake state, a rapid eye movement (REM) sleep state, a non-rapid eye movement (NREM) sleep state and a NREM category including drowsiness (N1), sleep (N2), and deep sleep (N3) or microstates within these NREM sleep categories.
Implementation A28 is the method of any of implementations A26 and A27, wherein the threshold impedance level is adjusted based on an ultradian rhythm of the individual, a circadian rhythm of the individual, or an infradian rhythm of the individual.
Implementation A29 is the method of any of implementations A24-A28, wherein the model impedance signature comprises a threshold rate of change of impedance, and processing the impedance signature comprises determining that a rate of change of impedance of the brain tissue of the individual is greater than the threshold rate of change of impedance indicated by the model impedance signature.
Implementation A30 is the method of implementation A29, wherein the threshold rate of change of impedance is adjusted based on a detected brain state of the individual, wherein the detected brain state is selected from an awake state, a rapid eye movement (REM) sleep state, a non-rapid eye movement (NREM) sleep state and a NREM category including drowsiness (N1), sleep (N2), and deep sleep (N3) or microstates within these NREM sleep categories.
Implementation A31 is the method of any of implementations A29 and A30, wherein the threshold rate of change of impedance is adjusted based on an ultradian rhythm of the individual, a circadian rhythm of the individual, or an infradian rhythm of the individual.
Implementation A32 is the method of any of implementations A1-A31, wherein the impedance describes a frequency dependent linear response voltage generated by a test probe electric current.
Implementation A33 is the method of implementation A32, wherein a frequency of the frequency dependent linear response voltage is between 0.001 hertz and 10,000 hertz.
Implementation A34 is the method of any of implementations A1-A33, wherein the impedance is used as a biomarker for at least one of: (a) cellular electrical interactions; (b) extracellular matrix electrical interactions; (c) interstitial extracellular electrolyte fluid electrical interactions; or (d) any combination of (a), (b), and (c).
Implementation A35 is the method of any of implementations A1-A34, wherein an impedance frequency spectrum is constructed using the at least one measurement of the impedance of the brain tissue.
Implementation A36 is the method of implementation A35, wherein the impedance frequency spectrum is used to track: a behavioral state of the individual; (b) a glymphatic system function of the individual; (c) a probability of seizure occurrence; or (d) any combination of (a), (b), and (c).
Implementation A37 is the system method any of implementations A1-A36, wherein temporal dynamics and correlations of impedance frequency dispersion in different anatomical brain regions are biomarkers used to track: (a) a behavioral state of the individual; (b) a glymphatic system function of the individual; (c) a probability of seizure occurrence; or (d) any combination of (a), (b), and (c).
Implementation A38 is the method of implementation A37, wherein the temporal dynamics and correlations of frequency dispersion are detected between 0.001 hertz and 10,000 hertz.
Implementation B1 is a system comprising an implanted device configured to capture at least one measurement of an impedance of brain tissue of an individual and a computing device configured to obtain, from the implanted device, the at least one measurement of the impedance of the brain tissue of the individual and in response to a prediction that the individual is likely to experience a seizure, initiate a remedial action that comprises at least one of alerting the individual or another user that the individual is likely to experience the seizure, applying an electrical stimulus to the brain tissue of the individual to reduce a likelihood of the seizure occurring or to reduce a severity of the seizure, or delivering a drug to the brain tissue of the individual to reduce a likelihood of the seizure occurring or to reduce a severity of the seizure, wherein the prediction that the individual is likely to experience a seizure is determined by processing an impedance signature determined from the at least one measurement of the impedance of the brain tissue of the individual.
Implementation B2 is the system of implementation B1, wherein the computing device operates an application that is configured to receive an individual input for indicating that the individual is experiencing a seizure or recently experienced a seizure.
Implementation B3 is the system of implementation B2, wherein the individual input is used to determine when the individual experience a seizure and impedance data from around a time the seizure occurred is collected and processed to determine an impedance signature that corresponds to the individual being likely to experience a seizure.
Implementation B4 is the system of any of implementations B1-B3, wherein the prediction is based on a detected change in impedance as compared to a threshold.
Implementation B5 is the system of implementation B4, wherein the threshold is patient specific.
Implementation B6 is the system of any of implementations B1-B5, wherein the alert is presented to an individual ten minutes before the seizure is expected to occur.
Implementation B7 is the system of any of implementations B1-B6, wherein the implantable device is configured to receive a two-point impedance measurement.
Implementation B8 is the system of any of implementations B1-B7, wherein the implantable device is configured to receive a four-point impedance measurement.
Implementation B9 is the system of any of implementations B1-B8, wherein the impedance describes a frequency dependent linear response voltage generated by a test probe electric current.
Implementation B10 is the system of implementation B9, wherein a frequency of the frequency dependent linear response voltage is between 0.001 hertz and 10,000 hertz.
Implementation B11 is the system of any of implementations B1-B10, wherein the impedance is used as a biomarker for at least one of: (a) cellular electrical interactions; (b) extracellular matrix electrical interactions; (c) interstitial extracellular electrolyte fluid electrical interactions; or (d) any combination of (a), (b), and (c).
Implementation B12 is the system of any of implementations B1-B11, wherein an impedance frequency spectrum is constructed using the at least one measurement of the impedance of the brain tissue.
Implementation B13 is the system of implementation B12, wherein the impedance frequency spectrum is used to track: a behavioral state of the individual; (b) a glymphatic system function of the individual; (c) a probability of seizure occurrence; or (d) any combination of (a), (b), and (c).
Implementation B14 is the system of any of implementations B1-B13, wherein temporal dynamics and correlations of impedance frequency dispersion in different anatomical brain regions are biomarkers used to track: (a) a behavioral state of the individual; (b) a glymphatic system function of the individual; (c) a probability of seizure occurrence; or (d) any combination of (a), (b), and (c).
Implementation B15 is the system of implementation B14, wherein the temporal dynamics and correlations of frequency dispersion are detected between 0.001 hertz and 10,000 hertz.
In situations in which the systems, methods, devices, and other techniques here collect personal information (e.g., context data) about users, or may make use of personal information, the users may be provided with an opportunity to control whether programs or features collect user information (e.g., information about a user's social network, social actions or activities, profession, a user's preferences, or a user's current location), or to control whether and/or how to receive content from the content server that may be more relevant to the user. In addition, certain data may be treated in one or more ways before it is stored or used, so that personally identifiable information is removed. For example, a user's identity may be treated so that no personally identifiable information can be determined for the user, or a user's geographic location may be generalized where location information is obtained (such as to a city, ZIP code, or state level), so that a particular location of a user cannot be determined. Thus, the user may have control over how information is collected about the user and used by a content server.
Although various implementations have been described in detail above, other modifications are possible. In addition, the logic flows depicted in the figures do not require the particular order shown, or sequential order, to achieve desirable results. In addition, other steps may be provided, or steps may be eliminated, from the described flows, and other components may be added to, or removed from, the described systems. Accordingly, other implementations are within the scope of the following claims.
This application claims the benefit of priority under 35 U.S.C. § 119(e) to U.S. provisional application Ser. No. 63/595,739 filed Nov. 2, 2023, the entire contents of which are incorporated by reference herein.
This invention was made with government support under NS095495 awarded by the National Institutes of Health. The government has certain rights in the invention.
Number | Date | Country | |
---|---|---|---|
63595739 | Nov 2023 | US |