Long-range temporal correlations in resting state alpha oscillations in major depressive disorder and obsessive-compulsive disorder

Introduction

Mental disorders are a significant concern in contemporary society, with a pressing need to identify biological markers. Long-range temporal correlations (LRTC) of brain rhythms have been widespread in clinical cohort studies, especially in major depressive disorder (MDD). However, research on LRTC in obsessive-compulsive disorder (OCD) is severely limited. Given the high co-occurrence of OCD and MDD, we conducted a comparative LRTC investigation. We assumed that the LRTC patterns will allow us to compare measures of brain cortical balance of excitation and inhibition in OCD and MDD, which will be useful in the area of differential diagnosis.

Methods

In this study, we used the 64-channel resting state EEG of 29 MDD participants, 26 OCD participants, and a control group of 37 volunteers. Detrended fluctuation analyzes was used to assess LRTC.

Results

Our results indicate that all scaling exponents of the three subject groups exhibited persistent LRTC of EEG oscillations. There was a tendency for LRTC to be higher in disorders than in controls, but statistically significant differences were found between the OCD and control groups in the entire frontal and left parietal occipital areas, and between the MDD and OCD groups in the middle and right frontal areas.

Discussion

We believe that these results indicate abnormalities in the inhibitory and excitatory neurotransmitter systems, predominantly affecting areas related to executive functions.

Accelerating spiking neural network simulations with PymoNNto and PymoNNtorch

Spiking neural network simulations are a central tool in Computational Neuroscience, Artificial Intelligence, and Neuromorphic Engineering research. A broad range of simulators and software frameworks for such simulations exist with different target application areas. Among these, PymoNNto is a recent Python-based toolbox for spiking neural network simulations that emphasizes the embedding of custom code in a modular and flexible way. While PymoNNto already supports GPU implementations, its backend relies on NumPy operations. Here we introduce PymoNNtorch, which is natively implemented with PyTorch while retaining PymoNNto's modular design. Furthermore, we demonstrate how changes to the implementations of common network operations in combination with PymoNNtorch's native GPU support can offer speed-up over conventional simulators like NEST, ANNarchy, and Brian 2 in certain situations. Overall, we show how PymoNNto's modular and flexible design in combination with PymoNNtorch's GPU acceleration and optimized indexing operations facilitate research and development of spiking neural networks in the Python programming language.

Empirical comparison of deep learning models for fNIRS pain decoding

Introduction

Pain assessment is extremely important in patients unable to communicate and it is often done by clinical judgement. However, assessing pain using observable indicators can be challenging for clinicians due to the subjective perceptions, individual differences in pain expression, and potential confounding factors. Therefore, the need for an objective pain assessment method that can assist medical practitioners. Functional near-infrared spectroscopy (fNIRS) has shown promising results to assess the neural function in response of nociception and pain. Previous studies have explored the use of machine learning with hand-crafted features in the assessment of pain.

Methods

In this study, we aim to expand previous studies by exploring the use of deep learning models Convolutional Neural Network (CNN), Long Short-Term Memory (LSTM), and (CNN-LSTM) to automatically extract features from fNIRS data and by comparing these with classical machine learning models using hand-crafted features.

Results

The results showed that the deep learning models exhibited favourable results in the identification of different types of pain in our experiment using only fNIRS input data. The combination of CNN and LSTM in a hybrid model (CNN-LSTM) exhibited the highest performance (accuracy = 91.2%) in our problem setting. Statistical analysis using one-way ANOVA with Tukey's (post-hoc) test performed on accuracies showed that the deep learning models significantly improved accuracy performance as compared to the baseline models.

Discussion

Overall, deep learning models showed their potential to learn features automatically without relying on manually-extracted features and the CNN-LSTM model could be used as a possible method of assessment of pain in non-verbal patients. Future research is needed to evaluate the generalisation of this method of pain assessment on independent populations and in real-life scenarios.

Multiscale co-simulation design pattern for neuroscience applications

Integration of information across heterogeneous sources creates added scientific value. Interoperability of data, tools and models is, however, difficult to accomplish across spatial and temporal scales. Here we introduce the toolbox Parallel Co-Simulation, which enables the interoperation of simulators operating at different scales. We provide a software science co-design pattern and illustrate its functioning along a neuroscience example, in which individual regions of interest are simulated on the cellular level allowing us to study detailed mechanisms, while the remaining network is efficiently simulated on the population level. A workflow is illustrated for the use case of The Virtual Brain and NEST, in which the CA1 region of the cellular-level hippocampus of the mouse is embedded into a full brain network involving micro and macro electrode recordings. This new tool allows integrating knowledge across scales in the same simulation framework and validating them against multiscale experiments, thereby largely widening the explanatory power of computational models.

The Locare workflow: representing neuroscience data locations as geometric objects in 3D brain atlases

Neuroscientists employ a range of methods and generate increasing amounts of data describing brain structure and function. The anatomical locations from which observations or measurements originate represent a common context for data interpretation, and a starting point for identifying data of interest. However, the multimodality and abundance of brain data pose a challenge for efforts to organize, integrate, and analyze data based on anatomical locations. While structured metadata allow faceted data queries, different types of data are not easily represented in a standardized and machine-readable way that allow comparison, analysis, and queries related to anatomical relevance. To this end, three-dimensional (3D) digital brain atlases provide frameworks in which disparate multimodal and multilevel neuroscience data can be spatially represented. We propose to represent the locations of different neuroscience data as geometric objects in 3D brain atlases. Such geometric objects can be specified in a standardized file format and stored as location metadata for use with different computational tools. We here present the Locare workflow developed for defining the anatomical location of data elements from rodent brains as geometric objects. We demonstrate how the workflow can be used to define geometric objects representing multimodal and multilevel experimental neuroscience in rat or mouse brain atlases. We further propose a collection of JSON schemas (LocareJSON) for specifying geometric objects by atlas coordinates, suitable as a starting point for co-visualization of different data in an anatomical context and for enabling spatial data queries.

In silico analyzes of the involvement of GPR55, CB1R and TRPV1: response to THC, contribution to temporal lobe epilepsy, structural modeling and updated evolution

Introduction

The endocannabinoid (eCB) system is named after the discovery that endogenous cannabinoids bind to the same receptors as the phytochemical compounds found in Cannabis. While endogenous cannabinoids include anandamide (AEA) and 2-arachidonoylglycerol (2-AG), exogenous phytocannabinoids include Δ-9 tetrahydrocannabinol (THC) and cannabidiol (CBD). These compounds finely tune neurotransmission following synapse activation, via retrograde signaling that activates cannabinoid receptor 1 (CB1R) and/or transient receptor potential cation channel subfamily V member 1 (TRPV1). Recently, the eCB system has been linked to several neurological diseases, such as neuro-ocular abnormalities, pain insensitivity, migraine, epilepsy, addiction and neurodevelopmental disorders. In the current study, we aim to: (i) highlight a potential link between the eCB system and neurological disorders, (ii) assess if THC exposure alters the expression of eCB-related genes, and (iii) identify evolutionary-conserved residues in CB1R or TRPV1 in light of their function.

Methods

To address this, we used several bioinformatic approaches, such as transcriptomic (Gene Expression Omnibus), protein–protein (STRING), phylogenic (BLASTP, MEGA) and structural (Phyre2, AutoDock, Vina, PyMol) analyzes.

Results

Using RNA sequencing datasets, we did not observe any dysregulation of eCB-related transcripts in major depressive disorders, bipolar disorder or schizophrenia in the anterior cingulate cortex, nucleus accumbens or dorsolateral striatum. Following in vivo THC exposure in adolescent mice, GPR55 was significantly upregulated in neurons from the ventral tegmental area, while other transcripts involved in the eCB system were not affected by THC exposure. Our results also suggest that THC likely induces neuroinflammation following in vitro application on mice microglia. Significant downregulation of TPRV1 occurred in the hippocampi of mice in which a model of temporal lobe epilepsy was induced, confirming previous observations. In addition, several transcriptomic dysregulations were observed in neurons of both epileptic mice and humans, which included transcripts involved in neuronal death. When scanning known interactions for transcripts involved in the eCB system (n = 13), we observed branching between the eCB system and neurophysiology, including proteins involved in the dopaminergic system. Our protein phylogenic analyzes revealed that CB1R forms a clade with CB2R, which is distinct from related paralogues such as sphingosine-1-phosphate, receptors, lysophosphatidic acid receptors and melanocortin receptors. As expected, several conserved residues were identified, which are crucial for CB1R receptor function. The anandamide-binding pocket seems to have appeared later in evolution. Similar results were observed for TRPV1, with conserved residues involved in receptor activation.

Conclusion

The current study found that GPR55 is upregulated in neurons following THC exposure, while TRPV1 is downregulated in temporal lobe epilepsy. Caution is advised when interpreting the present results, as we have employed secondary analyzes. Common ancestors for CB1R and TRPV1 diverged from jawless vertebrates during the late Ordovician, 450 million years ago. Conserved residues are identified, which mediate crucial receptor functions.

Enabling uncertainty estimation in neural networks through weight perturbation for improved Alzheimer’s disease classification

Background

The willingness to trust predictions formulated by automatic algorithms is key in a wide range of domains. However, a vast number of deep architectures are only able to formulate predictions without associated uncertainty.

Purpose

In this study, we propose a method to convert a standard neural network into a Bayesian neural network and estimate the variability of predictions by sampling different networks similar to the original one at each forward pass.

Methods

We combine our method with a tunable rejection-based approach that employs only the fraction of the data, i.e., the share that the model can classify with an uncertainty below a user-set threshold. We test our model in a large cohort of brain images from patients with Alzheimer's disease and healthy controls, discriminating the former and latter classes based on morphometric images exclusively.

Results

We demonstrate how combining estimated uncertainty with a rejection-based approach increases classification accuracy from 0.86 to 0.95 while retaining 75% of the test set. In addition, the model can select the cases to be recommended for, e.g., expert human evaluation due to excessive uncertainty. Importantly, our framework circumvents additional workload during the training phase by using our network “turned into Bayesian” to implicitly investigate the loss landscape in the neighborhood of each test sample in order to determine the reliability of the predictions.

Conclusion

We believe that being able to estimate the uncertainty of a prediction, along with tools that can modulate the behavior of the network to a degree of confidence that the user is informed about (and comfortable with), can represent a crucial step in the direction of user compliance and easier integration of deep learning tools into everyday tasks currently performed by human operators.

Improving the detection of sleep slow oscillations in electroencephalographic data

Study objectives

We aimed to build a tool which facilitates manual labeling of sleep slow oscillations (SOs) and evaluate the performance of traditional sleep SO detection algorithms on such a manually labeled data set. We sought to develop improved methods for SO detection.

Method

SOs in polysomnographic recordings acquired during nap time from ten older adults were manually labeled using a custom built graphical user interface tool. Three automatic SO detection algorithms previously used in the literature were evaluated on this data set. Additional machine learning and deep learning algorithms were trained on the manually labeled data set.

Results

Our custom built tool significantly decreased the time needed for manual labeling, allowing us to manually inspect 96,277 potential SO events. The three automatic SO detection algorithms showed relatively low accuracy (max. 61.08%), but results were qualitatively similar, with SO density and amplitude increasing with sleep depth. The machine learning and deep learning algorithms showed higher accuracy (best: 99.20%) while maintaining a low prediction time.

Conclusions

Accurate detection of SO events is important for investigating their role in memory consolidation. In this context, our tool and proposed methods can provide significant help in identifying these events.

Domain adaptation for EEG-based, cross-subject epileptic seizure prediction

The ability to predict the occurrence of an epileptic seizure is a safeguard against patient injury and health complications. However, a major challenge in seizure prediction arises from the significant variability observed in patient data. Common patient-specific approaches, which apply to each patient independently, often perform poorly for other patients due to the data variability. The aim of this study is to propose deep learning models which can handle this variability and generalize across various patients. This study addresses this challenge by introducing a novel cross-subject and multi-subject prediction models. Multiple-subject modeling broadens the scope of patient-specific modeling to account for the data from a dedicated ensemble of patients, thereby providing some useful, though relatively modest, level of generalization. The basic neural network architecture of this model is then adapted to cross-subject prediction, thereby providing a broader, more realistic, context of application. For accrued performance, and generalization ability, cross-subject modeling is enhanced by domain adaptation. Experimental evaluation using the publicly available CHB-MIT and SIENA data datasets shows that our multiple-subject model achieved better performance compared to existing works. However, the cross-subject faces challenges when applied to different patients. Finally, through investigating three domain adaptation methods, the model accuracy has been notably improved by 10.30% and 7.4% for the CHB-MIT and SIENA datasets, respectively.