Functional connectivity signatures of major depressive disorder: machine learning analysis of two multicenter neuroimaging studies.

The promise of machine learning has fueled the hope for developing diagnostic tools for psychiatry. Initial studies showed high accuracy for the identification of major depressive disorder (MDD) with resting-state connectivity, but progress has been hampered by the absence of large datasets. Here we used regular machine learning and advanced deep learning algorithms to differentiate patients with MDD from healthy controls and identify neurophysiological signatures of depression in two of the largest resting-state datasets for MDD. We obtained resting-state functional magnetic resonance imaging data from the REST-meta-MDD (N = 2338) and PsyMRI (N = 1039) consortia. Classification of functional connectivity matrices was done using support vector machines (SVM) and graph convolutional neural networks (GCN), and performance was evaluated using 5-fold cross-validation. Features were visualized using GCN-Explainer, an ablation study and univariate t-testing. The results showed a mean classification accuracy of 61% for MDD versus controls. Mean accuracy for classifying (non-)medicated subgroups was 62%. Sex classification accuracy was substantially better across datasets (73-81%). Visualization of the results showed that classifications were driven by stronger thalamic connections in both datasets, while nearly all other connections were weaker with small univariate effect sizes. These results suggest that whole brain resting-state connectivity is a reliable though poor biomarker for MDD, presumably due to disease heterogeneity as further supported by the higher accuracy for sex classification using the same methods. Deep learning revealed thalamic hyperconnectivity as a prominent neurophysiological signature of depression in both multicenter studies, which may guide the development of biomarkers in future studies.

Multimodal Data Integration Advances Longitudinal Prediction of the Naturalistic Course of Depression and Reveals a Multimodal Signature of Remission During 2-Year Follow-up.

BACKGROUND: The ability to predict the disease course of individuals with major depressive disorder (MDD) is essential for optimal treatment planning. Here, we used a data-driven machine learning approach to assess the predictive value of different sets of biological data (whole-blood proteomics, lipid metabolomics, transcriptomics, genetics), both separately and added to clinical baseline variables, for the longitudinal prediction of 2-year remission status in MDD at the individual-subject level. METHODS: Prediction models were trained and cross-validated in a sample of 643 patients with current MDD (2-year remission n = 325) and subsequently tested for performance in 161 individuals with MDD (2-year remission n = 82). RESULTS: Proteomics data showed the best unimodal data predictions (area under the receiver operating characteristic curve = 0.68). Adding proteomic to clinical data at baseline significantly improved 2-year MDD remission predictions (area under the receiver operating characteristic curve = 0.63 vs. 0.78, p = .013), while the addition of other omics data to clinical data did not yield significantly improved model performance. Feature importance and enrichment analysis revealed that proteomic analytes were involved in inflammatory response and lipid metabolism, with fibrinogen levels showing the highest variable importance, followed by symptom severity. Machine learning models outperformed psychiatrists' ability to predict 2-year remission status (balanced accuracy = 71% vs. 55%). CONCLUSIONS: This study showed the added predictive value of combining proteomic data, but not other omics data, with clinical data for the prediction of 2-year remission status in MDD. Our results reveal a novel multimodal signature of 2-year MDD remission status that shows clinical potential for individual MDD disease course predictions from baseline measurements.

Patient specific intracranial neural signatures of obsessions and compulsions in the ventral striatum.

Objective. Deep brain stimulation is a treatment option for patients with refractory obsessive-compulsive disorder. A new generation of stimulators hold promise for closed loop stimulation, with adaptive stimulation in response to biologic signals. Here we aimed to discover a suitable biomarker in the ventral striatum in patients with obsessive compulsive disorder using local field potentials.Approach.We induced obsessions and compulsions in 11 patients undergoing deep brain stimulation treatment using a symptom provocation task. Then we trained machine learning models to predict symptoms using the recorded intracranial signal from the deep brain stimulation electrodes.Main results.Average areas under the receiver operating characteristics curve were 62.1% for obsessions and 78.2% for compulsions for patient specific models. For obsessions it reached over 85% in one patient, whereas performance was near chance level when the model was trained across patients. Optimal performances for obsessions and compulsions was obtained at different recording sites.Significance. The results from this study suggest that closed loop stimulation may be a viable option for obsessive-compulsive disorder, but that intracranial biomarkers are patient and not disorder specific.Clinical Trial:Netherlands trial registry NL7486.

Probabilistic tractography recovers a rostrocaudal trajectory of connectivity variability in the human insular cortex.

The insular cortex of macaques has a wide spectrum of anatomical connections whose distribution is related to its heterogeneous cytoarchitecture. Although there is evidence of a similar cytoarchitectural arrangement in humans, the anatomical connectivity of the insula in the human brain has not yet been investigated in vivo. In the present work, we used in vivo probabilistic white-matter tractography and Laplacian eigenmaps (LE) to study the variation of connectivity patterns across insular territories in humans. In each subject and hemisphere, we recovered a rostrocaudal trajectory of connectivity variation ranging from the anterior dorsal and ventral insula to the dorsal caudal part of the long insular gyri. LE suggested that regional transitions among tractography patterns in the insula occur more gradually than in other brain regions. In particular, the change in tractography patterns was more gradual in the insula than in the medial premotor region, where a sharp transition between different tractography patterns was found. The recovered trajectory of connectivity variation in the insula suggests a relation between connectivity and cytoarchitecture in humans resembling that previously found in macaques: tractography seeds from the anterior insula were mainly found in limbic and paralimbic regions and in anterior parts of the inferior frontal gyrus, while seeds from caudal insular territories mostly reached parietal and posterior temporal cortices. Regions in the putative dysgranular insula displayed more heterogeneous connectivity patterns, with regional differences related to the proximity with either putative granular or agranular regions.

[Clinical course of COVID-19 in the Netherlands: an overview of 2607 patients in hospital during the first wave]. / Klinisch beloop van covid-19 in Nederland.

OBJECTIVE: To systematically collect clinical data from patients with a proven COVID-19 infection in the Netherlands. DESIGN: Data from 2579 patients with COVID-19 admitted to 10 Dutch centers in the period February to July 2020 are described. The clinical data are based on the WHO COVID case record form (CRF) and supplemented with patient characteristics of which recently an association disease severity has been reported. METHODS: Survival analyses were performed as primary statistical analysis. These Kaplan-Meier curves for time to (early) death (3 weeks) have been determined for pre-morbid patient characteristics and clinical, radiological and laboratory data at hospital admission. RESULTS: Total in-hospital mortality after 3 weeks was 22.2% (95% CI: 20.7% - 23.9%), hospital mortality within 21 days was significantly higher for elderly patients (> 70 years; 35, 0% (95% CI: 32.4% - 37.8%) and patients who died during the 21 days and were admitted to the intensive care (36.5% (95% CI: 32.1% - 41.3%)). Apart from that, in this Dutch population we also see a risk of early death in patients with co-morbidities (such as chronic neurological, nephrological and cardiac disorders and hypertension), and in patients with more home medication and / or with increased urea and creatinine levels. CONCLUSION: Early death due to a COVID-19 infection in the Netherlands appears to be associated with demographic variables (e.g. age), comorbidity (e.g. cardiovascular disease) but also disease char-acteristics at admission.

Structural neuroimaging biomarkers for obsessive-compulsive disorder in the ENIGMA-OCD consortium: medication matters.

No diagnostic biomarkers are available for obsessive-compulsive disorder (OCD). Here, we aimed to identify magnetic resonance imaging (MRI) biomarkers for OCD, using 46 data sets with 2304 OCD patients and 2068 healthy controls from the ENIGMA consortium. We performed machine learning analysis of regional measures of cortical thickness, surface area and subcortical volume and tested classification performance using cross-validation. Classification performance for OCD vs. controls using the complete sample with different classifiers and cross-validation strategies was poor. When models were validated on data from other sites, model performance did not exceed chance-level. In contrast, fair classification performance was achieved when patients were grouped according to their medication status. These results indicate that medication use is associated with substantial differences in brain anatomy that are widely distributed, and indicate that clinical heterogeneity contributes to the poor performance of structural MRI as a disease marker.

Individual prediction of psychotherapy outcome in posttraumatic stress disorder using neuroimaging data.

Trauma-focused psychotherapy is the first-line treatment for posttraumatic stress disorder (PTSD) but 30-50% of patients do not benefit sufficiently. We investigated whether structural and resting-state functional magnetic resonance imaging (MRI/rs-fMRI) data could distinguish between treatment responders and non-responders on the group and individual level. Forty-four male veterans with PTSD underwent baseline scanning followed by trauma-focused psychotherapy. Voxel-wise gray matter volumes were extracted from the structural MRI data and resting-state networks (RSNs) were calculated from rs-fMRI data using independent component analysis. Data were used to detect differences between responders and non-responders on the group level using permutation testing, and the single-subject level using Gaussian process classification with cross-validation. A RSN centered on the bilateral superior frontal gyrus differed between responders and non-responder groups (PFWE < 0.05) while a RSN centered on the pre-supplementary motor area distinguished between responders and non-responders on an individual-level with 81.4% accuracy (P < 0.001, 84.8% sensitivity, 78% specificity and AUC of 0.93). No significant single-subject classification or group differences were observed for gray matter volume. This proof-of-concept study demonstrates the feasibility of using rs-fMRI to develop neuroimaging biomarkers for treatment response, which could enable personalized treatment of patients with PTSD.

Individual Prediction of Behavioral Variant Frontotemporal Dementia Development Using Multivariate Pattern Analysis of Magnetic Resonance Imaging Data.

BACKGROUND: Patients with behavioral variant of frontotemporal dementia (bvFTD) initially may only show behavioral and/or cognitive symptoms that overlap with other neurological and psychiatric disorders. The diagnostic accuracy is dependent on progressive symptoms worsening and frontotemporal abnormalities on neuroimaging findings. Predictive biomarkers could facilitate the early detection of bvFTD. OBJECTIVE: To determine the prognostic accuracy of clinical and structural MRI data using a support vector machine (SVM) classification to predict the 2-year clinical follow-up diagnosis in a group of patients presenting late-onset behavioral changes. METHODS: Data from 73 patients were included and divided into probable/definite bvFTD (n = 18), neurological (n = 28), and psychiatric (n = 27) groups based on 2-year follow-up diagnosis. Grey-matter volumes were extracted from baseline structural MRI scans. SVM classifiers were used to perform three binary classifications: bvFTD versus neurological and psychiatric, bvFTD versus neurological, and bvFTD versus psychiatric group(s), and one multi-class classification. Classification performance was determined for clinical and neuroimaging data separately and their combination using 5-fold cross-validation. RESULTS: Accuracy of the binary classification tasks ranged from 72-82% (p < 0.001) with adequate sensitivity (67-79%), specificity (77-88%), and area-under-the-receiver-operator-curve (0.80-0.9). Multi-class accuracy ranged between 55-59% (p < 0.001). The combination of clinical and voxel-wise whole brain data showed the best performance overall. CONCLUSION: These results show the potential for automated early confirmation of diagnosis for bvFTD using machine learning analysis of clinical and neuroimaging data in a diverse and clinically relevant sample of patients.

Disentangling subgroups of participants recruiting shared as well as different brain regions for the execution of the verb generation task: A data-driven fMRI study.

The spatial pattern of task-related brain activity in fMRI studies might be expected to change according to several variables such as handedness and age. However this spatial heterogeneity might also be due to other unmodeled sources of inter-subject variability. Since group-level results reflect patterns of task-evoked brain activity common to most of the subjects in the sample, they could conceal the presence of subgroups recruiting other brain regions beyond the common pattern. To deal with these issues, data-driven methods can be used to detect the presence of sources of inter-subject variability that might be hard to identify and therefore model a priori. Here we assess the potential of Independent Component Analysis (ICA) to detect the presence of unexpected subgroups of participants. To this end, we acquired task-evoked fMRI data on 45 healthy adults using the verb generation (VGEN) task, in which participants are visually presented with the noun of an object of everyday use, and asked to covertly generate a verb describing the corresponding action. As expected, the task elicited activity in a temporo-parieto-frontal network typically found in previous VGEN experiments. We then quantified the contribution of every subject to nine task-related spatio-temporal processes identified by ICA. A cluster analysis of this quantity yielded three subgroups of participants. Differences between the three identified subgroups were distributed in left and right prefrontal, posterior parietal and extrastriate occipital regions. These results could not be explained by differences in sex, age or handedness across the participants. Furthermore, some regions where a significant difference was found between subgroups were not present in the group-level pattern of task-related activity. We discuss the potential application of this approach for characterizing brain activity in different subgroups of patients with neuropsychiatric or neurological conditions.

Increased Functional Connectivity Between Subcortical and Cortical Resting-State Networks in Autism Spectrum Disorder.

IMPORTANCE: Individuals with autism spectrum disorder (ASD) exhibit severe difficulties in social interaction, motor coordination, behavioral flexibility, and atypical sensory processing, with considerable interindividual variability. This heterogeneous set of symptoms recently led to investigating the presence of abnormalities in the interaction across large-scale brain networks. To date, studies have focused either on constrained sets of brain regions or whole-brain analysis, rather than focusing on the interaction between brain networks. OBJECTIVES: To compare the intrinsic functional connectivity between brain networks in a large sample of individuals with ASD and typically developing control subjects and to estimate to what extent group differences would predict autistic traits and reflect different developmental trajectories. DESIGN, SETTING, AND PARTICIPANTS: We studied 166 male individuals (mean age, 17.6 years; age range, 7-50 years) diagnosed as having DSM-IV-TR autism or Asperger syndrome and 193 typical developing male individuals (mean age, 16.9 years; age range, 6.5-39.4 years) using resting-state functional magnetic resonance imaging (MRI). Participants were matched for age, IQ, head motion, and eye status (open or closed) in the MRI scanner. We analyzed data from the Autism Brain Imaging Data Exchange (ABIDE), an aggregated MRI data set from 17 centers, made public in August 2012. MAIN OUTCOMES AND MEASURES: We estimated correlations between time courses of brain networks extracted using a data-driven method (independent component analysis). Subsequently, we associated estimates of interaction strength between networks with age and autistic traits indexed by the Social Responsiveness Scale. RESULTS: Relative to typically developing control participants, individuals with ASD showed increased functional connectivity between primary sensory networks and subcortical networks (thalamus and basal ganglia) (all t ≥ 3.13, P < .001 corrected). The strength of such connections was associated with the severity of autistic traits in the ASD group (all r ≥ 0.21, P < .0067 corrected). In addition, subcortico-cortical interaction decreased with age in the entire sample (all r ≤ -0.09, P < .012 corrected), although this association was significant only in typically developing participants (all r ≤ -0.13, P < .009 corrected). CONCLUSIONS AND RELEVANCE: Our results showing ASD-related impairment in the interaction between primary sensory cortices and subcortical regions suggest that the sensory processes they subserve abnormally influence brain information processing in individuals with ASD. This might contribute to the occurrence of hyposensitivity or hypersensitivity and of difficulties in top-down regulation of behavior.