A computational study of stimulus driven epileptic seizure abatement

This is the final version of the article. Available from Public Library of Science via the DOI in this record.Active brain stimulation to abate epileptic seizures has shown mixed success. In spike-wave (SW) seizures, where the seizure and background state were proposed to coexist, single-pulse stimulations have been suggested to be able to terminate the seizure prematurely. However, several factors can impact success in such a bistable setting. The factors contributing to this have not been fully investigated on a theoretical and mechanistic basis. Our aim is to elucidate mechanisms that influence the success of single-pulse stimulation in noise-induced SW seizures. In this work, we study a neural population model of SW seizures that allows the reconstruction of the basin of attraction of the background activity as a four dimensional geometric object. For the deterministic (noise-free) case, we show how the success of response to stimuli depends on the amplitude and phase of the SW cycle, in addition to the direction of the stimulus in state space. In the case of spontaneous noise-induced seizures, the basin becomes probabilistic introducing some degree of uncertainty to the stimulation outcome while maintaining qualitative features of the noise-free case. Additionally, due to the different time scales involved in SW generation, there is substantial variation between SW cycles, implying that there may not be a fixed set of optimal stimulation parameters for SW seizures. In contrast, the model suggests an adaptive approach to find optimal stimulation parameters patient-specifically, based on real-time estimation of the position in state space. We discuss how the modelling work can be exploited to rationally design a successful stimulation protocol for the abatement of SW seizures using real-time SW detection.This work was supported by the EPSRC (EP/K026992/1), the BBSRC, the DTC for Systems Biology (University of Manchester), and the Nanyang Technological University Singapore. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript

Phase Estimation and Phase Ambiguity Resolution by Message Passing

Several code-aided algorithms for phase estimation have recently been proposed. While some of them are ad-hoc, others are derived in a systematic way. The latter can be divided into two different classes: phase estimators derived from the expectation-maximization (EM) principle and estimators that are approximations of the sum-product message passing algorithm. In this paper, the main differences and similarities between these two classes of phase estimation algorithms are outlined and their performance and complexity is compared. Furthermore, an alternative criterion for phase ambiguity resolution is presented and compared to an EM based approach proposed earlier

Information capacity in the weak-signal approximation

We derive an approximate expression for mutual information in a broad class of discrete-time stationary channels with continuous input, under the constraint of vanishing input amplitude or power. The approximation describes the input by its covariance matrix, while the channel properties are described by the Fisher information matrix. This separation of input and channel properties allows us to analyze the optimality conditions in a convenient way. We show that input correlations in memoryless channels do not affect channel capacity since their effect decreases fast with vanishing input amplitude or power. On the other hand, for channels with memory, properly matching the input covariances to the dependence structure of the noise may lead to almost noiseless information transfer, even for intermediate values of the noise correlations. Since many model systems described in mathematical neuroscience and biophysics operate in the high noise regime and weak-signal conditions, we believe, that the described results are of potential interest also to researchers in these areas.Comment: 11 pages, 4 figures; accepted for publication in Physical Review

A hybrid feature selection approach for the early diagnosis of Alzheimer's disease

Objective. Recently, significant advances have been made in the early diagnosis of Alzheimer’s disease from EEG. However, choosing suitable measures is a challenging task. Among other measures, frequency Relative Power and loss of complexity have been used with promising results. In the present study we investigate the early diagnosis of AD using synchrony measures and frequency Relative Power on EEG signals, examining the changes found in different frequency ranges. Approach. We first explore the use of a single feature for computing the classification rate, looking for the best frequency range. Then, we present a multiple feature classification system that outperforms all previous results using a feature selection strategy. These two approaches are tested in two different databases, one containing MCI and healthy subjects (patients age: 71.9 ± 10.2, healthy subjects age: 71.7 ± 8.3), and the other containing Mild AD and healthy subjects (patients age: 77.6 ± 10.0; healthy subjects age: 69.4± 11.5). Main Results. Using a single feature to compute classification rates we achieve a performance of 78.33% for the MCI data set and of 97.56 % for Mild AD. Results are clearly improved using the multiple feature classification, where a classification rate of 95% is found for the MCI data set using 11 features, and 100% for the Mild AD data set using 4 features. Significance. The new features selection method described in this work may be a reliable tool that could help to design a realistic system that does not require prior knowledge of a patient's status. With that aim, we explore the standardization of features for MCI and Mild AD data sets with promising results

Disentangling cortical functional connectivity strength and topography reveals divergent roles of genes and environment

The human brain varies across individuals in its morphology, function, and cognitive capacities. Variability is particularly high in phylogenetically modern regions associated with higher order cognitive abilities, but its relationship to the layout and strength of functional networks is poorly understood. In this study we disentangled the variability of two key aspects of functional connectivity: strength and topography. We then compared the genetic and environmental influences on these two features. Genetic contribution is heterogeneously distributed across the cortex and differs for strength and topography. In heteromodal areas genes predominantly affect the topography of networks, while their connectivity strength is shaped primarily by random environmental influence such as learning. We identified peak areas of genetic control of topography overlapping with parts of the processing stream from primary areas to network hubs in the default mode network, suggesting the coordination of spatial configurations across those processing pathways. These findings provide a detailed map of the diverse contribution of heritability and individual experience to the strength and topography of functional brain architecture.Nanyang Technological UniversityPublished versionThis work was supported by the Medical University of Vienna, the Austrian Research Fund (FWF) [grants P 35189, P 34198, and I 3925-B27] in collaboration with the French National Research Agency (ANR), the Vienna Science and Technology Fund (WWTF) [LS20-065], the European Research Council Grant [866533-CORTIGRAD], the National Natural Science Foundation of China [Grant No. 81790652, No.81790650] and the NAM Advanced Biomedical Imaging Program [FY2016] between Nanyang Technological University, Singapore and Medical University of Vienna, Austria

Investigating model performance in language identification: beyond simple error statistics

Language development experts need tools that can automatically identify languages from fluent, conversational speech, and provide reliable estimates of usage rates at the level of an individual recording. However, language identification systems are typically evaluated on metrics such as equal error rate and balanced accuracy, applied at the level of an entire speech corpus. These overview metrics do not provide information about model performance at the level of individual speakers, recordings, or units of speech with different linguistic characteristics. Overview statistics may therefore mask systematic errors in model performance for some subsets of the data, and consequently, have worse performance on data derived from some subsets of human speakers, creating a kind of algorithmic bias. In the current paper, we investigate how well a number of language identification systems perform on individual recordings and speech units with different linguistic properties in the MERLIon CCS Challenge. The Challenge dataset features accented English-Mandarin code-switched child-directed speech.Comment: Accepted to Interspeech 2023, 5 pages, 5 figure

MERLIon CCS Challenge: A English-Mandarin code-switching child-directed speech corpus for language identification and diarization

To enhance the reliability and robustness of language identification (LID) and language diarization (LD) systems for heterogeneous populations and scenarios, there is a need for speech processing models to be trained on datasets that feature diverse language registers and speech patterns. We present the MERLIon CCS challenge, featuring a first-of-its-kind Zoom video call dataset of parent-child shared book reading, of over 30 hours with over 300 recordings, annotated by multilingual transcribers using a high-fidelity linguistic transcription protocol. The audio corpus features spontaneous and in-the-wild English-Mandarin code-switching, child-directed speech in non-standard accents with diverse language-mixing patterns recorded in a variety of home environments. This report describes the corpus, as well as LID and LD results for our baseline and several systems submitted to the MERLIon CCS challenge using the corpus.Comment: Accepted by Interspeech 2023, 5 pages, 2 figures, 3 table

Monte Carlo drift correction – quantifying the drift uncertainty of global climate models

Global climate models are susceptible to drift, causing spurious trends in output variables. Drift is often corrected using data from a control simulation. However, internal climate variability within the control simulation introduces uncertainty to the drift correction process. To quantify this drift uncertainty, we develop a probabilistic technique: Monte Carlo drift correction (MCDC). MCDC samples the standard error associated with drift in the control time series. We apply MCDC to an ensemble of global climate models from the Coupled Model Intercomparison Project Phase 6 (CMIP6). We find that drift correction partially addresses a problem related to drift: energy leakage. Nevertheless, the energy balance of several models remains suspect. We quantify the drift uncertainty of global quantities associated with the Earth's energy balance and thermal expansion of the ocean. When correcting drift in a cumulatively integrated energy flux, we find that it is preferable to integrate the flux before correcting the drift: an alternative method would be to correct the bias before integrating the flux, but this alternative method amplifies the drift uncertainty. Assuming that drift is linear likely leads to an underestimation of drift uncertainty. Time series with weak trends may be especially susceptible to drift uncertainty: for historical thermosteric sea level rise since the 1850s, the drift uncertainty can range from 3 to 24 mm, which is of comparable magnitude to the impact of omitting volcanic forcing in control simulations. Derived coefficients – such as the ocean's expansion efficiency of heat – can also be susceptible to drift uncertainty. When evaluating and analysing global climate model data that are susceptible to drift, researchers should consider drift uncertainty.</p

Recurrent Activity in Higher Order, Modality Non-Specific Brain Regions: A Granger Causality Analysis of Autobiographic Memory Retrieval

It has been proposed that the workings of the brain are mainly intrinsically generated recurrent neuronal activity, with sensory inputs as modifiers of such activity in both sensory and higher order modality non-specific regions. This is supported by the demonstration of recurrent neuronal activity in the visual system as a response to visual stimulation. In contrast recurrent activity has never been demonstrated before in higher order modality non-specific regions. Using magneto-encephalography and Granger causality analysis, we tested in a paralimbic network the hypothesis that stimulation may enhance causal recurrent interaction between higher-order, modality non-specific regions. The network includes anterior cingulate/medial prefrontal and posterior cingulate/medial parietal cortices together with pulvinar thalami, a network known to be effective in autobiographic memory retrieval and self-awareness. Autobiographic memory retrieval of previous personal judgments of visually presented words was used as stimuli. It is demonstrated that the prestimulus condition is characterized by causal, recurrent oscillations which are maximal in the lower gamma range. When retrieving previous judgments of visually presented adjectives, this activity is dramatically increased during the stimulus task as ascertained by Granger causality analysis. Our results confirm the hypothesis that stimulation may enhance causal interaction between higher order, modality non-specific brain regions, exemplified in a network of autobiographical memory retrieval