A machine-learning framework for robust and reliable prediction of short- and long-term treatment response in initially antipsychotic-naïve schizophrenia patients based on multimodal neuropsychiatric data.

The reproducibility of machine-learning analyses in computational psychiatry is a growing concern. In a multimodal neuropsychiatric dataset of antipsychotic-naïve, first-episode schizophrenia patients, we discuss a workflow aimed at reducing bias and overfitting by invoking simulated data in the design process and analysis in two independent machine-learning approaches, one based on a single algorithm and the other incorporating an ensemble of algorithms. We aimed to (1) classify patients from controls to establish the framework, (2) predict short- and long-term treatment response, and (3) validate the methodological framework. We included 138 antipsychotic-naïve, first-episode schizophrenia patients with data on psychopathology, cognition, electrophysiology, and structural magnetic resonance imaging (MRI). Perinatal data and long-term outcome measures were obtained from Danish registers. Short-term treatment response was defined as change in Positive And Negative Syndrome Score (PANSS) after the initial antipsychotic treatment period. Baseline diagnostic classification algorithms also included data from 151 matched controls. Both approaches significantly classified patients from healthy controls with a balanced accuracy of 63.8% and 64.2%, respectively. Post-hoc analyses showed that the classification primarily was driven by the cognitive data. Neither approach predicted short- nor long-term treatment response. Validation of the framework showed that choice of algorithm and parameter settings in the real data was successfully guided by results from the simulated data. In conclusion, this novel approach holds promise as an important step to minimize bias and obtain reliable results with modest sample sizes when independent replication samples are not available.

Pharmaco-electroencephalographic responses in the rat differ between active and inactive locomotor states.

Quantitative electroencephalography from freely moving rats is commonly used as a translational tool for predicting drug-effects in humans. We hypothesized that drug-effects may be expressed differently depending on whether the rat is in active locomotion or sitting still during recording sessions, and proposed automatic state-detection as a viable tool for estimating drug-effects free of hypo-/hyperlocomotion-induced effects. We aimed at developing a fully automatic and validated method for detecting two behavioural states: active and inactive, in one-second intervals and to use the method for evaluating ketamine, DOI, d-cycloserine, d-amphetamine, and diazepam effects specifically within each state. The developed state-detector attained high precision with more than 90% of the detected time correctly classified, and multiple differences between the two detected states were discovered. Ketamine-induced delta activity was found specifically related to locomotion. Ketamine and DOI suppressed theta and beta oscillations exclusively during inactivity. Characteristic gamma and high-frequency oscillations (HFO) enhancements of the NMDAR and 5HT2A modulators, speculated associated with locomotion, were profound and often largest during the inactive state. State-specific analyses, theoretically eliminating biases from altered occurrence of locomotion, revealed only few effects of d-amphetamine and diazepam. Overall, drug-effects were most abundant in the inactive state. In conclusion, this new validated and automatic locomotion state-detection method enables fast and reliable state-specific analysis facilitating discovery of state-dependent drug-effects and control for altered occurrence of locomotion. This may ultimately lead to better cross-species translation of electrophysiological effects of pharmacological modulations.

Is there a reasonable excuse for not providing post-operative analgesia when using animal models of peripheral neuropathic pain for research purposes?

INTRODUCTION: The induction of neuropathic pain-like behaviors in rodents often requires surgical intervention. This engages acute nociceptive signaling events that contribute to pain and stress post-operatively that from a welfare perspective demands peri-operative analgesic treatment. However, a large number of researchers avoid providing such care based largely on anecdotal opinions that it might interfere with model pathophysiology in the longer term. OBJECTIVES: To investigate effects of various peri-operative analgesic regimens encapsulating different mechanisms and duration of action, on the development of post-operative stress/welfare and pain-like behaviors in the Spared Nerve Injury (SNI)-model of neuropathic pain. METHODS: Starting on the day of surgery, male Sprague-Dawley rats were administered either vehicle (s.c.), carprofen (5.0mg/kg, s.c.), buprenorphine (0.1mg/kg s.c. or 1.0mg/kg p.o. in Nutella®), lidocaine/bupivacaine mixture (local irrigation) or a combination of all analgesics, with coverage from a single administration, and up to 72 hours. Post-operative stress and recovery were assessed using welfare parameters, bodyweight, food-consumption, and fecal corticosterone, and hindpaw mechanical allodynia was tested for assessing development of neuropathic pain for 28 days. RESULTS: None of the analgesic regimes compromised the development of mechanical allodynia. Unexpectedly, the combined treatment with 0.1mg/kg s.c. buprenorphine and carprofen for 72 hours and local irrigation with lidocaine/bupivacaine, caused severe adverse effects with peritonitis. This was not observed when the combination included a lower dose of buprenorphine (0.05mg/kg, s.c.), or when buprenorphine was administered alone (0.1mg/kg s.c. or 1.0mg/kg p.o.) for 72 hours. An elevated rate of wound dehiscence was observed especially in the combined treatment groups, underlining the need for balanced analgesia. Repeated buprenorphine injections had positive effects on body weight the first day after surgery, but depressive effects on food intake and body weight later during the first week. CONCLUSION: Post-operative analgesia does not appear to affect established neuropathic hypersensitivity outcome in the SNI model.

Sleep stability and transitions in patients with idiopathic REM sleep behavior disorder and patients with Parkinson's disease.

OBJECTIVE: Patients with idiopathic rapid eye movement (REM) sleep behavior disorder (iRBD) are at high risk of developing Parkinson's disease (PD). As wake/sleep-regulation is thought to involve neurons located in the brainstem and hypothalamic areas, we hypothesize that the neurodegeneration in iRBD/PD is likely to affect wake/sleep and REM/non-REM (NREM) sleep transitions. METHODS: We determined the frequency of wake/sleep and REM/NREM sleep transitions and the stability of wake (W), REM and NREM sleep as measured by polysomnography (PSG) in 27 patients with PD, 23 patients with iRBD, 25 patients with periodic leg movement disorder (PLMD) and 23 controls. Measures were computed based on manual scorings and data-driven labeled sleep staging. RESULTS: Patients with PD showed significantly lower REM stability than controls and patients with PLMD. Patients with iRBD had significantly lower REM stability compared with controls. Patients with PD and RBD showed significantly lower NREM stability and significantly more REM/NREM transitions than controls. CONCLUSIONS: We conclude that W, NREM and REM stability and transitions are progressively affected in iRBD and PD, probably reflecting the successive involvement of brain stem areas from early on in the disease. SIGNIFICANCE: Sleep stability and transitions determined by a data-driven approach could support the evaluation of iRBD and PD patients.

Data-driven modeling of sleep EEG and EOG reveals characteristics indicative of pre-Parkinson's and Parkinson's disease.

BACKGROUND: Manual scoring of sleep relies on identifying certain characteristics in polysomnograph (PSG) signals. However, these characteristics are disrupted in patients with neurodegenerative diseases. NEW METHOD: This study evaluates sleep using a topic modeling and unsupervised learning approach to identify sleep topics directly from electroencephalography (EEG) and electrooculography (EOG). PSG data from control subjects were used to develop an EOG and an EEG topic model. The models were applied to PSG data from 23 control subjects, 25 patients with periodic leg movements (PLMs), 31 patients with idiopathic REM sleep behavior disorder (iRBD) and 36 patients with Parkinson's disease (PD). The data were divided into training and validation datasets and features reflecting EEG and EOG characteristics based on topics were computed. The most discriminative feature subset for separating iRBD/PD and PLM/controls was estimated using a Lasso-regularized regression model. RESULTS: The features with highest discriminability were the number and stability of EEG topics linked to REM and N3, respectively. Validation of the model indicated a sensitivity of 91.4% and a specificity of 68.8% when classifying iRBD/PD patients. COMPARISON WITH EXISTING METHOD: The topics showed visual accordance with the manually scored sleep stages, and the features revealed sleep characteristics containing information indicative of neurodegeneration. CONCLUSIONS: This study suggests that the amount of N3 and the ability to maintain NREM and REM sleep have potential as early PD biomarkers. Data-driven analysis of sleep may contribute to the evaluation of neurodegenerative patients.

Automatic sleep classification using a data-driven topic model reveals latent sleep states.

BACKGROUND: The golden standard for sleep classification uses manual scoring of polysomnography despite points of criticism such as oversimplification, low inter-rater reliability and the standard being designed on young and healthy subjects. NEW METHOD: To meet the criticism and reveal the latent sleep states, this study developed a general and automatic sleep classifier using a data-driven approach. Spectral EEG and EOG measures and eye correlation in 1s windows were calculated and each sleep epoch was expressed as a mixture of probabilities of latent sleep states by using the topic model Latent Dirichlet Allocation. Model application was tested on control subjects and patients with periodic leg movements (PLM) representing a non-neurodegenerative group, and patients with idiopathic REM sleep behavior disorder (iRBD) and Parkinson's Disease (PD) representing a neurodegenerative group. The model was optimized using 50 subjects and validated on 76 subjects. RESULTS: The optimized sleep model used six topics, and the topic probabilities changed smoothly during transitions. According to the manual scorings, the model scored an overall subject-specific accuracy of 68.3 ± 7.44 (% µ ± σ) and group specific accuracies of 69.0 ± 4.62 (control), 70.1 ± 5.10 (PLM), 67.2 ± 8.30 (iRBD) and 67.7 ± 9.07 (PD). COMPARISON WITH EXISTING METHOD: Statistics of the latent sleep state content showed accordances to the sleep stages defined in the golden standard. However, this study indicates that sleep contains six diverse latent sleep states and that state transitions are continuous processes. CONCLUSIONS: The model is generally applicable and may contribute to the research in neurodegenerative diseases and sleep disorders.

Decreased sleep spindle density in patients with idiopathic REM sleep behavior disorder and patients with Parkinson's disease.

OBJECTIVE: To determine whether sleep spindles (SS) are potentially a biomarker for Parkinson's disease (PD). METHODS: Fifteen PD patients with REM sleep behavior disorder (PD+RBD), 15 PD patients without RBD (PD-RBD), 15 idiopathic RBD (iRBD) patients and 15 age-matched controls underwent polysomnography (PSG). SS were scored in an extract of data from control subjects. An automatic SS detector using a Matching Pursuit (MP) algorithm and a Support Vector Machine (SVM) was developed and applied to the PSG recordings. The SS densities in N1, N2, N3, all NREM combined and REM sleep were obtained and evaluated across the groups. RESULTS: The SS detector achieved a sensitivity of 84.7% and a specificity of 84.5%. At a significance level of α=1%, the iRBD and PD+RBD patients had a significantly lower SS density than the control group in N2, N3 and all NREM stages combined. At a significance level of α=5%, PD-RBD had a significantly lower SS density in N2 and all NREM stages combined. CONCLUSIONS: The lower SS density suggests involvement in pre-thalamic fibers involved in SS generation. SS density is a potential early PD biomarker. SIGNIFICANCE: It is likely that an automatic SS detector could be a supportive diagnostic tool in the evaluation of iRBD and PD patients.

Classification of iRBD and Parkinson's patients using a general data-driven sleep staging model built on EEG.

Sleep analysis is an important diagnostic tool for sleep disorders. However, the current manual sleep scoring is time-consuming as it is a crude discretization in time and stages. This study changes Esbroeck and Westover's [1] latent sleep staging model into a global model. The proposed data-driven method trained a topic mixture model on 10 control subjects and was applied on 10 other control subjects, 10 iRBD patients and 10 Parkinson's patients. In that way 30 topic mixture diagrams were obtained from which features reflecting distinct sleep architectures between control subjects and patients were extracted. Two features calculated on basis of two latent sleep states classified subjects as "control" or "patient" by a simple clustering algorithm. The mean sleep staging accuracy compared to classical AASM scoring was 72.4% for control subjects and a clustering of the derived features resulted in a sensitivity of 95% and a specificity of 80 %. This study demonstrates that frequency analysis of sleep EEG can be used for data-driven global sleep classification and that topic features separates iRBD and Parkinson's patients from control subjects.

Classification of iRBD and Parkinson's disease patients based on eye movements during sleep.

Patients suffering from the sleep disorder idiopathic rapid-eye-movement sleep behavior disorder (iRBD) have been observed to be in high risk of developing Parkinson's disease (PD). This makes it essential to analyze them in the search for PD biomarkers. This study aims at classifying patients suffering from iRBD or PD based on features reflecting eye movements (EMs) during sleep. A Latent Dirichlet Allocation (LDA) topic model was developed based on features extracted from two electrooculographic (EOG) signals measured as parts in full night polysomnographic (PSG) recordings from ten control subjects. The trained model was tested on ten other control subjects, ten iRBD patients and ten PD patients, obtaining a EM topic mixture diagram for each subject in the test dataset. Three features were extracted from the topic mixture diagrams, reflecting "certainty", "fragmentation" and "stability" in the timely distribution of the EM topics. Using a Naive Bayes (NB) classifier and the features "certainty" and "stability" yielded the best classification result and the subjects were classified with a sensitivity of 95 %, a specificity of 80% and an accuracy of 90 %. This study demonstrates in a data-driven approach, that iRBD and PD patients may exhibit abnorm form and/or timely distribution of EMs during sleep.

Gene expression profiles in relation to tension and dissociation in borderline personality disorder.

The biological underpinnings of borderline personality disorder (BPD) and its psychopathology including states of aversive tension and dissociation is poorly understood. Our goal was to examine transcriptional changes associated with states of tension or dissociation within individual patients in a pilot study. Dissociation is not only a critical symptom of BPD but has also been associated with higher risk for self-mutilation and depression. We conducted a whole blood gene expression profile analysis using quantitative PCR in 31 female inpatients with BPD. For each individual, two samples were drawn during a state of high tension and dissociation, while two samples were drawn at non-tension states. There was no association between gene expression and tension states. However, we could show that Interleukin-6 was positively correlated to dissociation scores, whereas Guanine nucleotide-binding protein G(s) subunit alpha isoforms, Mitogen-activated protein kinase 3 and 8, Guanine nucleotide-binding protein G(i) subunit alpha-2, Beta-arrestin-1 and 2, and Cyclic AMP-responsive element-binding protein were negatively correlated to dissociation. Our data point to a potential association of dissociation levels with the expression of genes involved in immune system regulation as well as cellular signalling/second-messenger systems. Major limitations of the study are the the possibly heterogeneous cell proportions in whole blood and the heterogeneous medication.

Separation of Parkinson's patients in early and mature stages from control subjects using one EOG channel.

In this study, polysomnographic left side EOG signals from ten control subjects, ten iRBD patients and ten Parkinson's patients were decomposed in time and frequency using wavelet transformation. A total of 28 features were computed as the means and standard deviations in energy measures from different reconstructed detail subbands across all sleep epochs during a whole night of sleep. A subset of features was chosen based on a cross validated Shrunken Centroids Regularized Discriminant Analysis, where the controls were treated as one group and the patients as another. Classification of the subjects was done by a leave-one-out validation approach using same method, and reached a sensitivity of 95%, a specificity of 70% and an accuracy of 86.7%. It was found that in the optimal subset of features, two hold lower frequencies reflecting the rapid eye movements and two hold higher frequencies reflecting EMG activity. This study demonstrates that both analysis of eye movements during sleep as well as EMG activity measured at the EOG channel hold potential of being biomarkers for Parkinson's disease.