Kurtosis and skewness of high-frequency brain signals are altered in paediatric epilepsy

Neural decoding and feature selection methods for closed-loop control of avoidance behavior

Objective.Many psychiatric disorders involve excessive avoidant or defensive behavior, such as avoidance in anxiety and trauma disorders or defensive rituals in obsessive-compulsive disorders. Developing algorithms to predict these behaviors from local field potentials (LFPs) could serve as the foundational technology for closed-loop control of such disorders. A significant challenge is identifying the LFP features that encode these defensive behaviors.Approach.We analyzed LFP signals from the infralimbic cortex and basolateral amygdala of rats undergoing tone-shock conditioning and extinction, standard for investigating defensive behaviors. We utilized a comprehensive set of neuro-markers across spectral, temporal, and connectivity domains, employing SHapley Additive exPlanations for feature importance evaluation within Light Gradient-Boosting Machine models. Our goal was to decode three commonly studied avoidance/defensive behaviors: freezing, bar-press suppression, and motion (accelerometry), examining the impact of different features on decoding performance.Main results.Band power and band power ratio between channels emerged as optimal features across sessions. High-gamma (80-150 Hz) power, power ratios, and inter-regional correlations were more informative than other bands that are more classically linked to defensive behaviors. Focusing on highly informative features enhanced performance. Across 4 recording sessions with 16 subjects, we achieved an average coefficient of determination of 0.5357 and 0.3476, and Pearson correlation coefficients of 0.7579 and 0.6092 for accelerometry jerk and bar press rate, respectively. Utilizing only the most informative features revealed differential encoding between accelerometry and bar press rate, with the former primarily through local spectral power and the latter via inter-regional connectivity. Our methodology demonstrated remarkably low training/inference time and memory usage, requiring<310 ms for training,<0.051 ms for inference, and 16.6 kB of memory, using a single core of AMD Ryzen Threadripper PRO 5995WX CPU.Significance.Our results demonstrate the feasibility of accurately decoding defensive behaviors with minimal latency, using LFP features from neural circuits strongly linked to these behaviors. This methodology holds promise for real-time decoding to identify physiological targets in closed-loop psychiatric neuromodulation.

Neural Decoding and Feature Selection Techniques for Closed-Loop Control of Defensive Behavior

Objective

Many psychiatric disorders involve excessive avoidant or defensive behavior, such as avoidance in anxiety and trauma disorders or defensive rituals in obsessive-compulsive disorders. Developing algorithms to predict these behaviors from local field potentials (LFPs) could serve as foundational technology for closed-loop control of such disorders. A significant challenge is identifying the LFP features that encode these defensive behaviors.

Approach

We analyzed LFP signals from the infralimbic cortex and basolateral amygdala of rats undergoing tone-shock conditioning and extinction, standard for investigating defensive behaviors. We utilized a comprehensive set of neuro-markers across spectral, temporal, and connectivity domains, employing SHapley Additive exPlanations for feature importance evaluation within Light Gradient-Boosting Machine models. Our goal was to decode three commonly studied avoidance/defensive behaviors: freezing, bar-press suppression, and motion (accelerometry), examining the impact of different features on decoding performance.

Main results

Band power and band power ratio between channels emerged as optimal features across sessions. High-gamma (80-150 Hz) power, power ratios, and inter-regional correlations were more informative than other bands that are more classically linked to defensive behaviors. Focusing on highly informative features enhanced performance. Across 4 recording sessions with 16 subjects, we achieved an average coefficient of determination of 0.5357 and 0.3476, and Pearson correlation coefficients of 0.7579 and 0.6092 for accelerometry jerk and bar press rate, respectively. Utilizing only the most informative features revealed differential encoding between accelerometry and bar press rate, with the former primarily through local spectral power and the latter via inter-regional connectivity. Our methodology demonstrated remarkably low time complexity, requiring <110 ms for training and <1 ms for inference.

Significance

Our results demonstrate the feasibility of accurately decoding defensive behaviors with minimal latency, using LFP features from neural circuits strongly linked to these behaviors. This methodology holds promise for real-time decoding to identify physiological targets in closed-loop psychiatric neuromodulation.

Predicting Alzheimer's disease CSF core biomarkers: a multimodal Machine Learning approach

Introduction

Alzheimer's disease (AD) is a progressive neurodegenerative disorder. Current core cerebrospinal fluid (CSF) AD biomarkers, widely employed for diagnosis, require a lumbar puncture to be performed, making them impractical as screening tools. Considering the role of sleep disturbances in AD, recent research suggests quantitative sleep electroencephalography features as potential non-invasive biomarkers of AD pathology. However, quantitative analysis of comprehensive polysomnography (PSG) signals remains relatively understudied. PSG is a non-invasive test enabling qualitative and quantitative analysis of a wide range of parameters, offering additional insights alongside other biomarkers. Machine Learning (ML) gained interest for its ability to discern intricate patterns within complex datasets, offering promise in AD neuropathology detection. Therefore, this study aims to evaluate the effectiveness of a multimodal ML approach in predicting core AD CSF biomarkers.

Methods

Mild-moderate AD patients were prospectively recruited for PSG, followed by testing of CSF and blood samples for biomarkers. PSG signals underwent preprocessing to extract non-linear, time domain and frequency domain statistics quantitative features. Multiple ML algorithms were trained using four subsets of input features: clinical variables (CLINVAR), conventional PSG parameters (SLEEPVAR), quantitative PSG signal features (PSGVAR) and a combination of all subsets (ALL). Cross-validation techniques were employed to evaluate model performance and ensure generalizability. Regression models were developed to determine the most effective variable combinations for explaining variance in the biomarkers.

Results

On 49 subjects, Gradient Boosting Regressors achieved the best results in estimating biomarkers levels, using different loss functions for each biomarker: least absolute deviation (LAD) for the Aβ42, least squares (LS) for p-tau and Huber for t-tau. The ALL subset demonstrated the lowest training errors for all three biomarkers, albeit with varying test performance. Specifically, the SLEEPVAR subset yielded the best test performance in predicting Aβ42, while the ALL subset most accurately predicted p-tau and t-tau due to the lowest test errors.

Conclusions

Multimodal ML can help predict the outcome of CSF biomarkers in early AD by utilizing non-invasive and economically feasible variables. The integration of computational models into medical practice offers a promising tool for the screening of patients at risk of AD, potentially guiding clinical decisions.

Epileptic Tissue Localization through Skewness-Based Functional Connectivity in the High-Frequency Band of Intracranial EEG

Functional connectivity analysis of intracranial electroencephalography (iEEG) plays an important role in understanding the mechanism of epilepsy and seizure dynamics. However, existing connectivity analysis is only suitable for low-frequency bands below 80 Hz. High-frequency oscillations (HFOs) and high-frequency activity (HFA) in the high-frequency band (80-500 Hz) are thought to be specific biomarkers in epileptic tissue localization. However, the transience in duration and variability of occurrence time and amplitudes of these events pose a challenge for conducting effective connectivity analysis. To deal with this problem, we proposed skewness-based functional connectivity (SFC) in the high-frequency band and explored its utility in epileptic tissue localization and surgical outcome evaluation. SFC comprises three main steps. The first step is the quantitative measurement of amplitude distribution asymmetry between HFOs/HFA and baseline activity. The second step is functional network construction on the basis of rank correlation of asymmetry across time. The third step is connectivity strength extraction from the functional network. Experiments were conducted in two separate datasets which consist of iEEG recordings from 59 patients with drug-resistant epilepsy. Significant difference (p<0.001) in connectivity strength was found between epileptic and non-epileptic tissue. Results were quantified via the receiver operating characteristic curve and the area under the curve (AUC). Compared with low-frequency bands, SFC demonstrated superior performance. With respect to pooled and individual epileptic tissue localization for seizure-free patients, AUCs were 0.66 (95% confidence interval (CI): 0.63-0.69) and (0.63 95% CI 0.56-0.71), respectively. For surgical outcome classification, the AUC was 0.75 (95% CI 0.59-0.85). Therefore, SFC can act as a promising assessment tool in characterizing the epileptic network and potentially provide better treatment options for patients with drug-resistant epilepsy.

Interictal localization of the epileptogenic zone: Utilizing the observed resonance behavior in the spectral band of surrounding inhibition

Introduction

The gold standard for identification of the epileptogenic zone (EZ) continues to be the visual inspection of electrographic changes around seizures' onset by experienced electroencephalography (EEG) readers. Development of an epileptogenic focus localization tool that can delineate the EZ from analysis of interictal (seizure-free) periods is still an open question of great significance for improved diagnosis (e.g., presurgical evaluation) and treatment of epilepsy (e.g., surgical outcome).

Methods

We developed an EZ interictal localization algorithm (EZILA) based on novel analysis of intracranial EEG (iEEG) using a univariate periodogram-type power measure, a straight-forward ranking approach, a robust dimensional reduction method and a clustering technique. Ten patients with temporal and extra temporal lobe epilepsies, and matching the inclusion criteria of having iEEG recordings at the epilepsy monitoring unit (EMU) and being Engel Class I ≥12 months post-surgery, were recruited in this study.

Results

In a nested k-fold cross validation statistical framework, EZILA assigned the highest score to iEEG channels within the EZ in all patients (10/10) during the first hour of the iEEG recordings and up to their first typical clinical seizure in the EMU (i.e., early interictal period). To further validate EZILA's performance, data from two new (Engel Class I) patients were analyzed in a double-blinded fashion; the EZILA successfully localized iEEG channels within the EZ from interictal iEEG in both patients.

Discussion

Out of the sampled brain regions, iEEG channels in the EZ were most frequently and maximally active in seizure-free (interictal) periods across patients in specific narrow gamma frequency band (∼60-80 Hz), which we have termed focal frequency band (FFB). These findings are consistent with the hypothesis that the EZ may interictally be regulated (controlled) by surrounding inhibitory neurons with resonance characteristics within this narrow gamma band.

Using search trends to analyze web-based users' behavior profiles connected with COVID-19 in mainland China: infodemiology study based on hot words and Baidu Index

Background

Mainland China, the world's most populous region, experienced a large-scale coronavirus disease 2019 (COVID-19) outbreak in 2020 and 2021, respectively. Existing infodemiology studies have primarily concentrated on the prospective surveillance of confirmed cases or symptoms which met the criterion for investigators; nevertheless, the actual impact regarding COVID-19 on the public and subsequent attitudes of different groups towards the COVID-19 epidemic were neglected.

Methods

This study aimed to examine the public web-based search trends and behavior patterns related to COVID-19 outbreaks in mainland China by using hot words and Baidu Index (BI). The initial hot words (the high-frequency words on the Internet) and the epidemic data (2019/12/01-2021/11/30) were mined from infodemiology platforms. The final hot words table was established by two-rounds of hot words screening and double-level hot words classification. Temporal distribution and demographic portraits of COVID-19 were queried by search trends service supplied from BI to perform the correlation analysis. Further, we used the parameter estimation to quantitatively forecast the geographical distribution of COVID-19 in the future.

Results

The final English-Chinese bilingual table was established including six domains and 32 subordinate hot words. According to the temporal distribution of domains and subordinate hot words in 2020 and 2021, the peaks of searching subordinate hot words and COVID-19 outbreak periods had significant temporal correlation and the subordinate hot words in COVID-19 Related and Territory domains were reliable for COVID-19 surveillance. Gender distribution results showed that Territory domain (the male proportion: 67.69%; standard deviation (SD): 5.88%) and Symptoms/Symptom and Public Health (the female proportion: 57.95%, 56.61%; SD: 0, 9.06%) domains were searched more by male and female groups respectively. The results of age distribution of hot words showed that people aged 20-50 (middle-aged people) had a higher online search intensity, and the group of 20-29, 30-39 years old focused more on Media and Symptoms/Symptom (proportion: 45.43%, 51.66%; SD: 15.37%, 16.59%) domains respectively. Finally, based on frequency rankings of searching hot words and confirmed cases in Mainland China, the epidemic situation of provinces and Chinese administrative divisions were divided into 5 levels of early-warning regions. Central, East and South China regions would be impacted again by the COVID-19 in the future.

Automated Detection of Left Bundle Branch Block from ECG Signal Utilizing the Maximal Overlap Discrete Wavelet Transform with ANFIS

Left bundle branch block (LBBB) is a common disorder in the heart’s electrical conduction system that leads to the ventricles’ uncoordinated contraction. The complete LBBB is usually associated with underlying heart failure and other cardiac diseases. Therefore, early automated detection is vital. This work aimed to detect the LBBB through the QRS electrocardiogram (ECG) complex segments taken from the MIT-BIH arrhythmia database. The used data contain 2655 LBBB (abnormal) and 1470 normal signals (i.e., 4125 total signals). The proposed method was employed in the following steps: (i) QRS segmentation and filtration, (ii) application of the Maximal Overlapped Discrete Wavelet Transform (MODWT) on the ECG R wave, (iii) selection of the detailed coefficients of the MODWT (D2, D3, D4), kurtosis, and skewness as extracted features to be fed into the Adaptive Neuro-Fuzzy Inference System (ANFIS) classifier. The obtained results proved that the proposed method performed well based on the achieved sensitivity, specificity, and classification accuracies of 99.81%, 100%, and 99.88%, respectively (F-Score is equal to 0.9990). Our results showed that the proposed method was robust and effective and could be used in real clinical situations.

Implementation of a Convolutional Neural Network for Eye Blink Artifacts Removal From the Electroencephalography Signal

Electroencephalography (EEG) signals are disrupted by technical and physiological artifacts. One of the most common artifacts is the natural activity that results from the movement of the eyes and the blinking of the subject. Eye blink artifacts (EB) spread across the entire head surface and make EEG signal analysis difficult. Methods for the elimination of electrooculography (EOG) artifacts, such as independent component analysis (ICA) and regression, are known. The aim of this article was to implement the convolutional neural network (CNN) to eliminate eye blink artifacts. To train the CNN, a method for augmenting EEG signals was proposed. The results obtained from the CNN were compared with the results of the ICA and regression methods for the generated and real EEG signals. The results obtained indicate a much better performance of the CNN in the task of removing eye-blink artifacts, in particular for the electrodes located in the central part of the head.

Interictal high frequency background activity as a biomarker of epileptogenic tissue

High frequency oscillations (HFOs) are very brief events that are a well-established biomarker of the epileptogenic zone (EZ) but are rare and comprise only a tiny fraction of the total recorded EEG. We hypothesize that the interictal high frequency ‘background’ data, which has received little attention but represents the majority of the EEG record, also may contain additional, novel information for identifying the EZ. We analysed intracranial EEG (30–500 Hz frequency range) acquired from 24 patients who underwent resective surgery. We computed 38 quantitative features based on all usable, interictal data (63–307 h per subject), excluding all detected HFOs. We assessed association between each feature and the seizure onset zone (SOZ) and resected volume (RV) using logistic regression. A pathology score per channel was also created via principle component analysis and logistic regression, using hold-out-one-patient cross-validation to avoid in-sample training. Association of the pathology score with the SOZ and RV was quantified using an asymmetry measure. Many features were associated with the SOZ: 23/38 features had odds ratios >1.3 or <0.7 and 17/38 had odds ratios different than zero with high significance (P < 0.001/39, logistic regression with Bonferroni Correction). The pathology score, the rate of HFOs, and their channel-wise product were each strongly associated with the SOZ [median asymmetry ≥0.44, good surgery outcome patients; median asymmetry ≥0.40, patients with other outcomes; 95% confidence interval (CI) > 0.27 in both cases]. The pathology score and the channel-wise product also had higher asymmetry with respect to the SOZ than the HFO rate alone (median difference in asymmetry ≥0.18, 95% CI >0.05). These results support that the high frequency background data contains useful information for determining the EZ, distinct and complementary to information from detected HFOs. The concordance between the high frequency activity pathology score and the rate of HFOs appears to be a better biomarker of epileptic tissue than either measure alone.

Diagnostic added value of interictal magnetic source imaging in presurgical evaluation of persons with epilepsy: A prospective blinded study

BACKGROUND AND PURPOSE

In presurgical evaluation for epilepsy surgery, information is sourced from various imaging modalities to accurately localize the epileptogenic zone. Magnetoencephalography (MEG) is a newer noninvasive technique for localization. However, there is limited literature to evaluate if MEG provides additional advantage over the conventional imaging modalities in clinical decision making. The objective of this study was to assess the diagnostic added value of MEG in decision making before epilepsy surgery.

METHOD

This was a prospective observational study. Patients underwent 3 h of recording in a MEG scanner, and the resulting localizations were compared with other complimentary investigations. Added value of MEG (considered separately from high-density electroencephalography) was defined as the frequency of cases in which (i) the information provided by magnetic source imaging (MSI) avoided implantation of intracranial electrodes and the patient was directly cleared for surgery, and (ii) MSI indicated additional substrates for implantation of intracranial electrodes. Postoperative seizure freedom was used as the diagnostic reference by which to measure the localizing accuracy of MSI.

RESULTS

A total of 102 patients underwent epilepsy surgery. MEG provided nonredundant information, which contributed to deciding the course of surgery in 33% of the patients, and prevented intracranial recordings in 19%. A total of 76% of the patients underwent surgical resection in sublobes concordant with MSI localization, and the diagnostic odds ratio for good (Engel I) outcome in these patients was 2.3 (95% confidence interval 0.68, 7.86; p = 0.183) after long-term follow-up of 36 months.

CONCLUSION

Magnetic source imaging yields additional useful information which can significantly alter as well as improve the surgical strategy for persons with epilepsy.

Neuromagnetic high frequency spikes are a new and noninvasive biomarker for localization of epileptogenic zones

OBJECTIVE

One barrier hindering high frequency brain signals (HFBS, >80 Hz) from wide clinical applications is that the brain generates both pathological and physiological HFBS. This study was to find specific biomarkers for localizing epileptogenic zones (EZs).

METHODS

Twenty three children with drug-resistant epilepsy and age/sex matched healthy controls were studied with magnetoencephalography (MEG). High frequency oscillations (HFOs, > 4 oscillatory waveforms) and high frequency spikes (HFSs, > 1 spiky or sharp waveforms) in 80-250 Hz and 250-600 Hz bands were blindly detected with an artificial intelligence method and validated with visual inspection. The magnitude of HFOs and HFSs were quantified with spectral analyses. Sources of HFSs and HFOs were localized and compared with clinical EZs determined by invasive recordings and surgical outcomes.

RESULTS

HFOs in 80-250 Hz and 250-600 Hz were identified in both epilepsy patients (18/23, 12/23, respectively) and healthy controls (6/23, 4/23, respectively). HFSs in 80-250 Hz and 250-600 Hz were detected in patients (16/23, 11/23, respectively) but not in healthy controls. A combination of HFOs and HFSs localized EZs for 22 (22/23, 96%) patients.

CONCLUSIONS

The results indicate, for the first time, that HFSs are a newer and more specific biomarker than HFOs for localizing EZs because HFOs appeared in both epilepsy patients and healthy controls while HFSs appeared only in epilepsy patients.

Delineation of epileptogenic zones with high frequency magnetic source imaging based on kurtosis and skewness

BACKGROUND

Neuromagnetic high frequency brain signals (HFBS, > 80 Hz) are a new biomarker for localization of epileptogenic zones (EZs) for pediatric epilepsy.

METHODS

Twenty three children with drug-resistant epilepsy and age/sex matched healthy controls were studied with magnetoencephalography (MEG). Epileptic HFBS in 80-250 Hz and 250-600 Hz were quantitatively determined by comparing with normative controls in terms of kurtosis and skewness. Magnetic sources of epileptic HFBS were localized and then compared to clinical EZs determined by invasive recordings and surgical outcomes.

RESULTS

Kurtosis and skewness of HFBS were significantly elevated in epilepsy patients compared to healthy controls (p < 0,001 and p < 0.0001, respectively). Sources of elevated MEG signals in comparison to normative data were co-localized to EZs for 22 (22/23, 96 %) patients.

CONCLUSIONS

The results indicate, for the first time, that epileptic HFBS can be noninvasively quantified by measuring kurtosis and skewness in MEG data. Magnetic source imaging based on kurtosis and skewness can accurately localize EZs.

SIGNIFICANCE

Source imaging of kurtosis and skewness of MEG HFBS provides a novel way for preoperative localization of EZs for epilepsy surgery.

Distinguishing false and true positive detections of high frequency oscillations

Objective.

High frequency oscillations (HFOs) are a promising biomarker of tissue that instigates seizures. However, ambiguous data and random background fluctuations can cause any HFO detector (human or automated) to falsely label non-HFO data as an HFO (a false positive detection). The objective of this paper was to identify quantitative features of HFOs that distinguish between true and false positive detections.

Approach.

Feature selection was performed using background data in multi-day, interictal intracranial recordings from ten patients. We selected the feature most similar between randomly selected segments of background data and HFOs detected in surrogate background data (false positive detections by construction). We then compared these results with fuzzy clustering of detected HFOs in clinical data to verify the feature’s applicability. We validated the feature is sensitive to false versus true positive HFO detections by using an independent data set (six subjects) scored for HFOs by three human reviewers. Lastly, we compared the effect of redacting putative false positive HFO detections on the distribution of HFOs across channels and their association with seizure onset zone (SOZ) and resected volume (RV).

Main results.

Of the 15 analyzed features, the analysis selected only skewness of the curvature (skewCurve). The feature was validated in human scored data to be associated with distinguishing true and false positive HFO detections. Automated HFO detections with higher skewCurve were more focal based on entropy measures and had increased localization to both the SOZ and RV.

Significance.

We identified a quantitative feature of HFOs which helps distinguish between true and false positive detections. Redacting putative false positive HFO detections improves the specificity of HFOs as a biomarker of epileptic tissue.