Functional brain networks of patients with epilepsy exhibit pronounced multiscale periodicities, which correlate with seizure onset.

Epileptic seizure detection and prediction by using noninvasive measurements such as scalp EEG signals or invasive, intracranial recordings, has been at the heart of epilepsy studies for at least three decades. To this end, the most common approach has been to consider short-length recordings (several seconds to a few minutes) around a seizure, aiming to identify significant changes that occur before or during seizures. An inherent assumption in this approach is the presence of a relatively constant EEG activity in the interictal period, which is interrupted by seizure occurrence. Here, we examine this assumption by using long-duration scalp EEG data (21-94 hr) in nine patients with epilepsy, based on which we construct functional brain networks. Our results reveal that these networks vary over time in a periodic fashion, exhibiting multiple peaks at periods ranging between 1 and 24 hr. The effects of seizure onset on the functional brain network properties were found to be considerably smaller in magnitude compared to the changes due to these inherent periodic cycles. Importantly, the properties of the identified network periodic components (instantaneous phase) were found to be strongly correlated to seizure onset, especially for the periodicities around 3 and 5 hr. These correlations were found to be largely absent between EEG signal periodicities and seizure onset, suggesting that higher specificity may be achieved by using network-based metrics. In turn, this implies that more robust seizure detection and prediction can be achieved if longer term underlying functional brain network periodic variations are taken into account.

Graph Theoretical Characteristics of EEG-Based Functional Brain Networks in Patients With Epilepsy: The Effect of Reference Choice and Volume Conduction.

It is well-established that both volume conduction and the choice of recording reference (montage) affect the correlation measures obtained from scalp EEG, both in the time and frequency domains. As a result, a number of correlation measures have been proposed aiming to reduce these effects. In our previous work, we have showed that scalp-EEG based functional brain networks in patients with epilepsy exhibit clear periodic patterns at different time scales and that these patterns are strongly correlated to seizure onset, particularly at shorter time scales (around 3 and 5 h), which has important clinical implications. In the present work, we use the same long-duration clinical scalp EEG data (multiple days) to investigate the extent to which the aforementioned results are affected by the choice of reference choice and correlation measure, by considering several widely used montages as well as correlation metrics that are differentially sensitive to the effects of volume conduction. Specifically, we compare two standard and commonly used linear correlation measures, cross-correlation in the time domain, and coherence in the frequency domain, with measures that account for zero-lag correlations: corrected cross-correlation, imaginary coherence, phase lag index, and weighted phase lag index. We show that the graphs constructed with corrected cross-correlation and WPLI are more stable across different choices of reference. Also, we demonstrate that all the examined correlation measures revealed similar periodic patterns in the obtained graph measures when the bipolar and common reference (Cz) montage were used. This includes circadian-related periodicities (e.g., a clear increase in connectivity during sleep periods as compared to awake periods), as well as periodicities at shorter time scales (around 3 and 5 h). On the other hand, these results were affected to a large degree when the average reference montage was used in combination with standard cross-correlation, coherence, imaginary coherence, and PLI, which is likely due to the low number of electrodes and inadequate electrode coverage of the scalp. Finally, we demonstrate that the correlation between seizure onset and the brain network periodicities is preserved when corrected cross-correlation and WPLI were used for all the examined montages. This suggests that, even in the standard clinical setting of EEG recording in epilepsy where only a limited number of scalp EEG measurements are available, graph-theoretic quantification of periodic patterns using appropriate montage, and correlation measures corrected for volume conduction provides useful insights into seizure onset.

Efficient algorithms for counting and reporting segregating sites in genomic sequences.

The number of segregating sites provides an indicator of the degree of DNA sequence variation that is present in a sample, and has been of great interest to the biological, pharmaceutical and medical professions. In this paper, we first provide linear- and expected-sublinear-time algorithms for finding all the segregating sites of a given set of DNA sequences. We also describe a data structure for tracking segregating sites in a set of sequences, such that every time the set is updated with the insertion of a new sequence or removal of an existing one, the segregating sites are updated accordingly without the need to re-scan the entire set of sequences.

Unsupervised detection and removal of muscle artifacts from scalp EEG recordings using canonical correlation analysis, wavelets and random forests.

OBJECTIVE: This paper proposes supervised and unsupervised algorithms for automatic muscle artifact detection and removal from long-term EEG recordings, which combine canonical correlation analysis (CCA) and wavelets with random forests (RF). METHODS: The proposed algorithms first perform CCA and continuous wavelet transform of the canonical components to generate a number of features which include component autocorrelation values and wavelet coefficient magnitude values. A subset of the most important features is subsequently selected using RF and labelled observations (supervised case) or synthetic data constructed from the original observations (unsupervised case). The proposed algorithms are evaluated using realistic simulation data as well as 30min epochs of non-invasive EEG recordings obtained from ten patients with epilepsy. RESULTS: We assessed the performance of the proposed algorithms using classification performance and goodness-of-fit values for noisy and noise-free signal windows. In the simulation study, where the ground truth was known, the proposed algorithms yielded almost perfect performance. In the case of experimental data, where expert marking was performed, the results suggest that both the supervised and unsupervised algorithm versions were able to remove artifacts without affecting noise-free channels considerably, outperforming standard CCA, independent component analysis (ICA) and Lagged Auto-Mutual Information Clustering (LAMIC). CONCLUSION: The proposed algorithms achieved excellent performance for both simulation and experimental data. Importantly, for the first time to our knowledge, we were able to perform entirely unsupervised artifact removal, i.e. without using already marked noisy data segments, achieving performance that is comparable to the supervised case. SIGNIFICANCE: Overall, the results suggest that the proposed algorithms yield significant future potential for improving EEG signal quality in research or clinical settings without the need for marking by expert neurophysiologists, EMG signal recording and user visual inspection.

Epileptic seizure onset correlates with long term EEG functional brain network properties.

We investigated the correlation of epileptic seizure onset times with long term EEG functional brain network properties. To do so, we constructed binary functional brain networks from long-term, multichannel electroencephalographic data recorded from nine patients with epilepsy. The corresponding network properties were quantified using the average network degree. It was found that the network degree (as well as other network properties such as the network efficiency and clustering coefficient) exhibited large fluctuations over time; however, it also exhibited specific periodic temporal structure over different time scales (1.5hr-24hr periods) that was consistent across subjects. We investigated the correlation of the phases of these network periodicities with the seizure onset by using circular statistics. The results showed that the instantaneous phases of the 3.5hr, 5.5hr, 12hr and 24hr network degree periodic components are not uniformly distributed, suggesting that functional network properties are related to seizure generation and occurrence.

Automatic detection and removal of muscle artifacts from scalp EEG recordings in patients with epilepsy.

Automatic detection and removal of muscle artifacts plays an important role in long-term scalp electroencephalography (EEG) monitoring, and muscle artifact detection algorithms have been intensively investigated. This paper proposes an algorithm for automatic muscle artifacts detection and removal using canonical correlation analysis (CCA) and wavelet transform (WT) in epochs from long-term EEG recordings. The proposed method first performs CCA analysis and then conducts wavelet decomposition on the canonical components within a specific frequency range and selects a subset of the wavelet coefficients for subsequent processing. A set of features, including the mean of wavelet coefficients and the canonical component autocorrelation values, are extracted from the above analysis and subsequently used as input in a random forest (RF) classifier. The RF classifier produces a similarity measure between observations and selects a subset of the most important features by comparing the original data with a set of synthetic data that is constructed based on the latter. The RF predictor output is finally used in combination with unsupervised clustering algorithms to discriminate between contaminated and non-contaminated EEG epochs. The proposed method is evaluated in epochs of 30 min from scalp EEG recordings obtained from three patients with epilepsy and yields a sensitivity of 71% and 80%, as well as a specificity of 81% and 85% for k-means and spectral clustering, respectively.

Periodicity in functional brain networks: application to scalp EEG from epilepsy patients.

Seizure detection and prediction studies using scalp- or intracranial-EEG measurements often focus on short-length recordings around the occurrence of the seizure, normally ranging between several seconds and up to a few minutes before and after the event. The underlying assumption in these studies is the presence of a relatively constant EEG activity in the interictal period, that is presumably interrupted by the occurrence of a seizure, at the time the seizure starts or slightly earlier. In this study, we put this assumption under test, by examining long-duration scalp EEG recordings, ranging between 22 and 72 hours, of five patients with epilepsy. For each patient, we construct functional brain networks, by calculating correlations between the scalp electrodes, and examine how these networks vary in time. The results suggest not only that the network varies over time, but it does so in a periodic fashion, with periods ranging between 11 and 25 hours.

Epigenetics of human T cells during the G0-->G1 transition.

We investigated functional epigenetic changes that occur in primary human T lymphocytes during entry into the cell cycle and mapped these at the single-nucleosome level by ChIP-chip on tiling arrays for chromosomes 1 and 6. We show that nucleosome loss and flanking active histone marks define active transcriptional start sites (TSSs). Moreover, these signatures are already set at many inducible genes in quiescent cells prior to cell stimulation. In contrast, there is a dearth of the inactive histone mark H3K9me3 at the TSS, and under-representation of H3K9me2 and H3K9me3 defines the body of active genes. At the DNA level, cytosine methylation (meC) is enriched for nucleosomes that remain at the TSS, whereas in general there is a dearth of meC at TSSs. Furthermore, a drop in meC also marks 3' transcription termination, and a peak of meC occurs at stop codons. This mimics the 3' nucleosomal distribution in yeast, which we show does not occur in human T cells.