Enhancing Multi-Center Generalization of Machine Learning-Based Depression Diagnosis From Resting-State fMRI

Diagnostic machine learning applications on clinical populations using functional near infrared spectroscopy: a review

Functional near-infrared spectroscopy (fNIRS) and its interaction with machine learning (ML) is a popular research topic for the diagnostic classification of clinical disorders due to the lack of robust and objective biomarkers. This review provides an overview of research on psychiatric diseases by using fNIRS and ML. Article search was carried out and 45 studies were evaluated by considering their sample sizes, used features, ML methodology, and reported accuracy. To our best knowledge, this is the first review that reports diagnostic ML applications using fNIRS. We found that there has been an increasing trend to perform ML applications on fNIRS-based biomarker research since 2010. The most studied populations are schizophrenia (n = 12), attention deficit and hyperactivity disorder (n = 7), and autism spectrum disorder (n = 6) are the most studied populations. There is a significant negative correlation between sample size (>21) and accuracy values. Support vector machine (SVM) and deep learning (DL) approaches were the most popular classifier approaches (SVM = 20) (DL = 10). Eight of these studies recruited a number of participants more than 100 for classification. Concentration changes in oxy-hemoglobin (ΔHbO) based features were used more than concentration changes in deoxy-hemoglobin (ΔHb) based ones and the most popular ΔHbO-based features were mean ΔHbO (n = 11) and ΔHbO-based functional connections (n = 11). Using ML on fNIRS data might be a promising approach to reveal specific biomarkers for diagnostic classification.

Spectral Graph Neural Network-Based Multi-Atlas Brain Network Fusion for Major Depressive Disorder Diagnosis

Major Depressive Disorder (MDD) imposes a substantial burden within the healthcare domain, impacting millions of individuals worldwide. Functional Magnetic Resonance Imaging (fMRI) has emerged as a promising tool for the objective diagnosis of MDD, enabling the investigation of functional connectivity patterns in the brain associated with this disorder. However, most existing methods focus on a single brain atlas, which limits their ability to capture the complex, multi-scale nature of functional brain networks. To address these limitations, we propose a novel multi-atlas fusion method that incorporates early and late fusion in a unified framework. Our method introduces the concept of the holistic Functional Connectivity Network (FCN), which captures both intra-atlas relationships within individual atlases and inter-regional relationships between atlases with different brain parcellation scales. This comprehensive representation enables the identification of potential disease-related patterns associated with MDD in the early stage of our framework. Moreover, by decoding the holistic FCN from various perspectives through multiple spectral Graph Convolutional Neural Networks and fusing their results with decision-level ensembles, we further improve the performance of MDD diagnosis. Our approach is easily implemented with minimal modifications to existing model structures and demonstrates a robust performance across different baseline models. Our method, evaluated on public resting-state fMRI datasets, surpasses the current multi-atlas fusion methods, enhancing the accuracy of MDD diagnosis. The proposed novel multi-atlas fusion framework provides a more reliable MDD diagnostic technique. Experimental results show our approach outperforms both single- and multi-atlas-based methods, demonstrating its effectiveness in advancing MDD diagnosis.

Analyzing 20 years of Resting-State fMRI Research: Trends and collaborative networks revealed

Resting-state functional magnetic resonance imaging (rs-fMRI), initially proposed by Biswal et al. in 1995, has emerged as a pivotal facet of neuroimaging research. Its ability to examine brain activity during the resting state without the need for explicit tasks or stimuli has made it an integral component of brain imaging studies. In recent years, rs-fMRI has witnessed substantial growth and found widespread application in the investigation of functional connectivity within the brain. To delineate the developmental trajectory of rs-fMRI over the past two decades, we conducted a comprehensive analysis using bibliometric tool Citespace. Our analysis encompassed publication trends, authorship networks, institutional affiliations, international collaborations, as well as emergent themes in references and keywords. Our study reveals a remarkable increase in the volume of rs-fMRI publications over the past two decades, underscoring the burgeoning interest and potential within this field. Harvard University stands out as the institution with the highest number of research papers published in the realm of RS-fMRI, while the United States holds the highest overall influence in this domain. The recent emergence of keywords such as "machine learning" and "default mode," coupled with citation surges in reference to rs-fMRI, have paved new avenues for research within this field. Our study underscores the critical importance of integrating machine learning techniques into rs-fMRI investigations, offering valuable insights into brain function and disease diagnosis. These findings hold profound significance for the field of neuroscience and may furnish insights for future research employing rs-fMRI as a diagnostic tool for a wide array of neurological disorders, thus emphasizing its pivotal role and potential as a tool for investigating brain functionality.

The diagnostic performance of machine learning based on resting-state functional magnetic resonance imaging data for major depressive disorders: a systematic review and meta-analysis

Objective

Machine learning (ML) has been widely used to detect and evaluate major depressive disorder (MDD) using neuroimaging data, i.e., resting-state functional magnetic resonance imaging (rs-fMRI). However, the diagnostic efficiency is unknown. The aim of the study is to conduct an updated meta-analysis to evaluate the diagnostic performance of ML based on rs-fMRI data for MDD.

Methods

English databases were searched for relevant studies. The Quality Assessment of Diagnostic Accuracy Studies (QUADAS-2) was used to assess the methodological quality of the included studies. A random-effects meta-analytic model was implemented to investigate the diagnostic efficiency, including sensitivity, specificity, diagnostic odds ratio (DOR), and area under the curve (AUC). Regression meta-analysis and subgroup analysis were performed to investigate the cause of heterogeneity.

Results

Thirty-one studies were included in this meta-analysis. The pooled sensitivity, specificity, DOR, and AUC with 95% confidence intervals were 0.80 (0.75, 0.83), 0.83 (0.74, 0.82), 14.00 (9, 22.00), and 0.86 (0.83, 0.89), respectively. Substantial heterogeneity was observed among the studies included. The meta-regression showed that the leave-one-out cross-validation (loocv) (sensitivity: p < 0.01, specificity: p < 0.001), graph theory (sensitivity: p < 0.05, specificity: p < 0.01), n > 100 (sensitivity: p < 0.001, specificity: p < 0.001), simens equipment (sensitivity: p < 0.01, specificity: p < 0.001), 3.0T field strength (Sensitivity: p < 0.001, specificity: p = 0.04), and Beck Depression Inventory (BDI) (sensitivity: p = 0.04, specificity: p = 0.06) might be the sources of heterogeneity. Furthermore, the subgroup analysis showed that the sample size (n > 100: sensitivity: 0.71, specificity: 0.72, n < 100: sensitivity: 0.81, specificity: 0.79), the different levels of disease evaluated by the Hamilton Depression Rating Scale (HDRS/HAMD) (mild vs. moderate vs. severe: sensitivity: 0.52 vs. 0.86 vs. 0.89, specificity: 0.62 vs. 0.78 vs. 0.82, respectively), the depression scales in patients with comparable levels of severity. (BDI vs. HDRS/HAMD: sensitivity: 0.86 vs. 0.87, specificity: 0.78 vs. 0.80, respectively), and the features (graph vs. functional connectivity: sensitivity: 0.84 vs. 0.86, specificity: 0.76 vs. 0.78, respectively) selected might be the causes of heterogeneity.

Conclusion

ML showed high accuracy for the automatic diagnosis of MDD. Future studies are warranted to promote the potential use of these classification algorithms in clinical settings.

Understanding and predicting future relapse in depression from resting state functional connectivity and self-referential processing

BACKGROUND

The recurrent nature of Major Depressive Disorder (MDD) asks for a better understanding of mechanisms underlying relapse. Previously, self-referential processing abnormalities have been linked to vulnerability for relapse. We investigated whether abnormalities in self-referential cognitions and functioning of associated brain-networks persist upon remission and predict relapse.

METHODS

Remitted recurrent MDD patients (n = 48) and never-depressed controls (n = 23) underwent resting-state fMRI scanning at baseline and were additionally assessed for their implicit depressed self-associations and ruminative behaviour. A template-based dual regression approach was used to investigate between-group differences in default mode, cingulo-opercular and frontoparietal network resting-state functional connectivity (RSFC). Additional prediction of relapse status at 18-month follow-up was investigated within patients using both regression analyses and machine learning classifiers.

RESULTS

Remitted patients showed higher rumination, but no implicit depressed self-associations or RSFC abnormalities were observed between patients and controls. Nevertheless, relapse was related to i) baseline RSFC between the ventral default mode network and the precuneus, dorsomedial frontal gyrus, and inferior occipital lobe, ii) implicit self-associations, and iii) uncontrollability of ruminative thinking, when controlled for depressive symptomatology. Moreover, preliminary machine learning classifiers demonstrated that RSFC within the investigated networks predicted relapse on an individual basis.

CONCLUSIONS

Remitted MDD patients seem to be commonly characterized by abnormal rumination, but not by implicit self-associations or abnormalities in relevant brain networks. Nevertheless, relapse was predicted by self-related cognitions and default mode RSFC during remission, suggesting that variations in self-relevant processing play a role in the complex dynamics associated with the vulnerability to developing recurrent depressive episodes.

CLINICAL TRIAL REGISTRATION

Netherlands Trial Register, August 18, 2015, trial number NL53205.042.15.

Sampling inequalities affect generalization of neuroimaging-based diagnostic classifiers in psychiatry

BACKGROUND

The development of machine learning models for aiding in the diagnosis of mental disorder is recognized as a significant breakthrough in the field of psychiatry. However, clinical practice of such models remains a challenge, with poor generalizability being a major limitation.

METHODS

Here, we conducted a pre-registered meta-research assessment on neuroimaging-based models in the psychiatric literature, quantitatively examining global and regional sampling issues over recent decades, from a view that has been relatively underexplored. A total of 476 studies (n = 118,137) were included in the current assessment. Based on these findings, we built a comprehensive 5-star rating system to quantitatively evaluate the quality of existing machine learning models for psychiatric diagnoses.

RESULTS

A global sampling inequality in these models was revealed quantitatively (sampling Gini coefficient (G) = 0.81, p < .01), varying across different countries (regions) (e.g., China, G = 0.47; the USA, G = 0.58; Germany, G = 0.78; the UK, G = 0.87). Furthermore, the severity of this sampling inequality was significantly predicted by national economic levels (β =  - 2.75, p < .001, R²_adj = 0.40; r =  - .84, 95% CI: - .41 to - .97), and was plausibly predictable for model performance, with higher sampling inequality for reporting higher classification accuracy. Further analyses showed that lack of independent testing (84.24% of models, 95% CI: 81.0-87.5%), improper cross-validation (51.68% of models, 95% CI: 47.2-56.2%), and poor technical transparency (87.8% of models, 95% CI: 84.9-90.8%)/availability (80.88% of models, 95% CI: 77.3-84.4%) are prevailing in current diagnostic classifiers despite improvements over time. Relating to these observations, model performances were found decreased in studies with independent cross-country sampling validations (all p < .001, BF₁₀ > 15). In light of this, we proposed a purpose-built quantitative assessment checklist, which demonstrated that the overall ratings of these models increased by publication year but were negatively associated with model performance.

CONCLUSIONS

Together, improving sampling economic equality and hence the quality of machine learning models may be a crucial facet to plausibly translating neuroimaging-based diagnostic classifiers into clinical practice.

Evaluation of Risk of Bias in Neuroimaging-Based Artificial Intelligence Models for Psychiatric Diagnosis: A Systematic Review

IMPORTANCE

Neuroimaging-based artificial intelligence (AI) diagnostic models have proliferated in psychiatry. However, their clinical applicability and reporting quality (ie, feasibility) for clinical practice have not been systematically evaluated.

OBJECTIVE

To systematically assess the risk of bias (ROB) and reporting quality of neuroimaging-based AI models for psychiatric diagnosis.

EVIDENCE REVIEW

PubMed was searched for peer-reviewed, full-length articles published between January 1, 1990, and March 16, 2022. Studies aimed at developing or validating neuroimaging-based AI models for clinical diagnosis of psychiatric disorders were included. Reference lists were further searched for suitable original studies. Data extraction followed the CHARMS (Checklist for Critical Appraisal and Data Extraction for Systematic Reviews of Prediction Modeling Studies) and PRISMA (Preferred Reporting Items for Systematic Reviews and Meta-analyses) guidelines. A closed-loop cross-sequential design was used for quality control. The PROBAST (Prediction Model Risk of Bias Assessment Tool) and modified CLEAR (Checklist for Evaluation of Image-Based Artificial Intelligence Reports) benchmarks were used to systematically evaluate ROB and reporting quality.

FINDINGS

A total of 517 studies presenting 555 AI models were included and evaluated. Of these models, 461 (83.1%; 95% CI, 80.0%-86.2%) were rated as having a high overall ROB based on the PROBAST. The ROB was particular high in the analysis domain, including inadequate sample size (398 of 555 models [71.7%; 95% CI, 68.0%-75.6%]), poor model performance examination (with 100% of models lacking calibration examination), and lack of handling data complexity (550 of 555 models [99.1%; 95% CI, 98.3%-99.9%]). None of the AI models was perceived to be applicable to clinical practices. Overall reporting completeness (ie, number of reported items/number of total items) for the AI models was 61.2% (95% CI, 60.6%-61.8%), and the completeness was poorest for the technical assessment domain with 39.9% (95% CI, 38.8%-41.1%).

CONCLUSIONS AND RELEVANCE

This systematic review found that the clinical applicability and feasibility of neuroimaging-based AI models for psychiatric diagnosis were challenged by a high ROB and poor reporting quality. Particularly in the analysis domain, ROB in AI diagnostic models should be addressed before clinical application.

Ensemble graph neural network model for classification of major depressive disorder using whole-brain functional connectivity

Major depressive disorder (MDD) is characterized by impairments in mood and cognitive functioning, and it is a prominent source of global disability and stress. A functional magnetic resonance imaging (fMRI) can aid clinicians in their assessments of individuals for the identification of MDD. Herein, we employ a deep learning approach to the issue of MDD classification. Resting-state fMRI data from 821 individuals with MDD and 765 healthy controls (HCs) is employed for investigation. An ensemble model based on graph neural network (GNN) has been created with the goal of identifying patients with MDD among HCs as well as differentiation between first-episode and recurrent MDDs. The graph convolutional network (GCN), graph attention network (GAT), and GraphSAGE models serve as a base models for the ensemble model that was developed with individual whole-brain functional networks. The ensemble's performance is evaluated using upsampling and downsampling, along with 10-fold cross-validation. The ensemble model achieved an upsampling accuracy of 71.18% and a downsampling accuracy of 70.24% for MDD and HC classification. While comparing first-episode patients with recurrent patients, the upsampling accuracy is 77.78% and the downsampling accuracy is 71.96%. According to the findings of this study, the proposed GNN-based ensemble model achieves a higher level of accuracy and suggests that our model produces can assist healthcare professionals in identifying MDD.

A systematic review on the potential use of machine learning to classify major depressive disorder from healthy controls using resting state fMRI measures

BACKGROUND

Major Depressive Disorder (MDD) is a psychiatric disorder characterized by functional brain deficits, as documented by resting-state functional magnetic resonance imaging (rs-fMRI) studies.

AIMS

In recent years, some studies used machine learning (ML) approaches, based on rs-fMRI features, for classifying MDD from healthy controls (HC). In this context, this review aims to provide a comprehensive overview of the results of these studies.

DESIGN

The studies research was performed on 3 online databases, examining English-written articles published before August 5, 2022, that performed a two-class ML classification using rs-fMRI features. The search resulted in 20 eligible studies.

RESULTS

The reviewed studies showed good performance metrics, with better performance achieved when the dataset was restricted to a more homogeneous group in terms of disease severity. Regions within the default mode network, salience network, and central executive network were reported as the most important features in the classification algorithms.

LIMITATIONS

The small sample size together with the methodological and clinical heterogeneity limited the generalizability of the findings.

CONCLUSIONS

In conclusion, ML applied to rs-fMRI features can be a valid approach to classify MDD and HC subjects and to discover features that can be used for additional investigation of the pathophysiology of the disease.

Multi-centre classification of functional neurological disorders based on resting-state functional connectivity

•

Using machine learning on multi-centre data, FND patients were successfully classified with an accuracy of 72%.

•

The angular- and supramarginal gyri, cingular- and insular cortex, and the hippocampus were the most discriminant regions.

•

To provide diagnostic utility, future studies must include patients with similar symptoms but different diagnoses.

Background

Patients suffering from functional neurological disorder (FND) experience disabling neurological symptoms not caused by an underlying classical neurological disease (such as stroke or multiple sclerosis). The diagnosis is made based on reliable positive clinical signs, but clinicians often require additional time- and cost consuming medical tests and examinations. Resting-state functional connectivity (RS FC) showed its potential as an imaging-based adjunctive biomarker to help distinguish patients from healthy controls and could represent a “rule-in” procedure to assist in the diagnostic process. However, the use of RS FC depends on its applicability in a multi-centre setting, which is particularly susceptible to inter-scanner variability. The aim of this study was to test the robustness of a classification approach based on RS FC in a multi-centre setting.

Methods

This study aimed to distinguish 86 FND patients from 86 healthy controls acquired in four different centres using a multivariate machine learning approach based on whole-brain resting-state functional connectivity. First, previously published results were replicated in each centre individually (intra-centre cross-validation) and its robustness across inter-scanner variability was assessed by pooling all the data (pooled cross-validation). Second, we evaluated the generalizability of the method by using data from each centre once as a test set, and the data from the remaining centres as a training set (inter-centre cross-validation).

Results

FND patients were successfully distinguished from healthy controls in the replication step (accuracy of 74%) as well as in each individual additional centre (accuracies of 73%, 71% and 70%). The pooled cross validation confirmed that the classifier was robust with an accuracy of 72%. The results survived post-hoc adjustment for anxiety, depression, psychotropic medication intake, and symptom severity. The most discriminant features involved the angular- and supramarginal gyri, sensorimotor cortex, cingular- and insular cortex, and hippocampal regions. The inter-centre validation step did not exceed chance level (accuracy below 50%).

Conclusions

The results demonstrate the applicability of RS FC to correctly distinguish FND patients from healthy controls in different centres and its robustness against inter-scanner variability. In order to generalize its use across different centres and aim for clinical application, future studies should work towards optimization of acquisition parameters and include neurological and psychiatric control groups presenting with similar symptoms.

An insight into diagnosis of depression using machine learning techniques: a systematic review

BACKGROUND

In this modern era, depression is one of the most prevalent mental disorders from which millions of individuals are affected today. The symptoms of depression are heterogeneous and often coincide with other disorders such as bipolar disorder, Parkinson's, schizophrenia, etc. It is a serious mental illness that may lead to other health problems if left untreated. Currently, identifying individuals with depression is totally based on the expertise of the clinician's experience. In order to assist clinicians in identifying the characteristics and classifying depressed people, different types of data modalities and machine learning techniques have been incorporated by researchers in this field. This study aims to find the answers to some important questions related to the trend of publications, data modality, machine learning models, dataset usage, pre-processing techniques and feature extraction and selection techniques that are prevalent and guide the direction of future research on depression diagnosis.

METHODS

This systematic review was conducted using a broad range of articles from two major databases: IEEE Xplore and PubMed. Studies ranging from the years 2011 to April 2021 were retrieved from the databases resulting in a total of 590 articles (53 articles from the IEEE Xplore database and 537 articles from the PubMed database). Out of those, the articles which satisfied the defined inclusion criteria were investigated for further analysis.

RESULTS

A total of 135 articles were identified and analysed for this review. High growth in the number of publications has been observed in recent years. Furthermore, significant diversity in the use of data modalities and machine learning classifiers has also been noted in this study. fMRI data with an SVM classifier was found to be the most popular choice among researchers. In most of the studies, data scarcity and small sample size, particularly for neuroimaging data are major concerns. The use of identical data pre-processing tools for similar data modalities can be seen. This study also provides statistical analysis of the current framework with respect to the modality, machine learning classifier, sample size and accuracy by applying one-way ANOVA and the Tukey - Kramer test.

CONCLUSION

The results indicate that an effective fusion of machine learning techniques with a potential data modality has a promising future for assisting clinicians in automatic depression diagnosis.

Resting-state brain activity can predict target-independent aptitude in fMRI-neurofeedback training

Neurofeedback (NF) aptitude, which refers to an individual's ability to change brain activity through NF training, has been reported to vary significantly from person to person. The prediction of individual NF aptitudes is critical in clinical applications to screen patients suitable for NF treatment. In the present study, we extracted the resting-state functional brain connectivity (FC) markers of NF aptitude, independent of NF-targeting brain regions. We combined the data from fMRI-NF studies targeting four different brain regions at two independent sites (obtained from 59 healthy adults and six patients with major depressive disorder) to collect resting-state fMRI data associated with aptitude scores in subsequent fMRI-NF training. We then trained the multiple regression models to predict the individual NF aptitude scores from the resting-state fMRI data using a discovery dataset from one site and identified six resting-state FCs that predicted NF aptitude. Subsequently, the reproducibility of the prediction model was validated using independent test data from another site. The identified FC model revealed that the posterior cingulate cortex was the functional hub among the brain regions and formed predictive resting-state FCs, suggesting that NF aptitude may be involved in the attentional mode-orientation modulation system's characteristics in task-free resting-state brain activity.

AUTOMATED RECOGNITION OF DEPRESSION FROM FEWER-SHOT LEANING IN RESTING-STATE fMRI WITH ReHo USING DEEP CONVOLUTIONAL NEURAL NETWORK

Up to now, there is still the absence of research about depression recognition using resting-state functional magnetic resonance imaging (rest_fMRI) and deep learning. Previous studies have shown that regional homogeneity (ReHo) of rest_fMRI (rest_ReHo_fMRI) is a characterization of the functional synchronization of adjacent voxels in brain regions, and the mental and behavioral abnormalities in depression are due to an imbalance of ReHo synchronization in some brain functional areas. Accordingly, this paper presents a method for depression recognition using rest_ReHo_fMRI. First, the rest_ReHo_fMRI is extracted from the preprocessed rest-fMRI by calculation. Then, deep convolutional networks (such as VGG16) pretrained on ImageNet are used to automatically complete extracting the classification features from rest_ReHo_fMRI. Finally, the Kernel Extreme Learning Machine (KELM) was used to classify the depression. The results of the test set show that the proposed method achieves 89.07% in sensitivity and 89.74% in specificity. This study suggests that features of rest_ReHo_fMRI can be used as biomarkers to distinguish depression from normal people.

Evidence of shared and distinct functional and structural brain signatures in schizophrenia and autism spectrum disorder

Schizophrenia (SZ) and autism spectrum disorder (ASD) share considerable clinical features and intertwined historical roots. It is greatly needed to explore their similarities and differences in pathophysiologic mechanisms. We assembled a large sample size of neuroimaging data (about 600 SZ patients, 1000 ASD patients, and 1700 healthy controls) to study the shared and unique brain abnormality of the two illnesses. We analyzed multi-scale brain functional connectivity among functional networks and brain regions, intra-network connectivity, and cerebral gray matter density and volume. Both SZ and ASD showed lower functional integration within default mode and sensorimotor domains, but increased interaction between cognitive control and default mode domains. The shared abnormalties in intra-network connectivity involved default mode, sensorimotor, and cognitive control networks. Reduced gray matter volume and density in the occipital gyrus and cerebellum were observed in both illnesses. Interestingly, ASD had overall weaker changes than SZ in the shared abnormalities. Interaction between visual and cognitive regions showed disorder-unique deficits. In summary, we provide strong neuroimaging evidence of the convergent and divergent changes in SZ and ASD that correlated with clinical features.

Modeling Heterogeneous Brain Dynamics of Depression and Melancholia Using Energy Landscape Analysis

Our current understanding of melancholic depression is shaped by its position in the depression spectrum. The lack of consensus on how it should be treated-whether as a subtype of depression, or as a distinct disorder altogethe-interferes with the recovery of suffering patients. In this study, we analyzed brain state energy landscape models of melancholic depression, in contrast to healthy and non-melancholic energy landscapes. Our analyses showed significant group differences on basin energy, basin frequency, and transition dynamics in several functional brain networks such as basal ganglia, dorsal default mode, and left executive control networks. Furthermore, we found evidences suggesting the connection between energy landscape characteristics (basin characteristics) and depressive symptom scores (BDI-II and SHAPS). These results indicate that melancholic depression is distinguishable from its non-melancholic counterpart, not only in terms of depression severity, but also in brain dynamics.

Common Brain Networks Between Major Depressive-Disorder Diagnosis and Symptoms of Depression That Are Validated for Independent Cohorts

Large-scale neuroimaging data acquired and shared by multiple institutions are essential to advance neuroscientific understanding of pathophysiological mechanisms in psychiatric disorders, such as major depressive disorder (MDD). About 75% of studies that have applied machine learning technique to neuroimaging have been based on diagnoses by clinicians. However, an increasing number of studies have highlighted the difficulty in finding a clear association between existing clinical diagnostic categories and neurobiological abnormalities. Here, using resting-state functional magnetic resonance imaging, we determined and validated resting-state functional connectivity related to depression symptoms that were thought to be directly related to neurobiological abnormalities. We then compared the resting-state functional connectivity related to depression symptoms with that related to depression diagnosis that we recently identified. In particular, for the discovery dataset with 477 participants from 4 imaging sites, we removed site differences using our recently developed harmonization method and developed a brain network prediction model of depression symptoms (Beck Depression Inventory-II [BDI] score). The prediction model significantly predicted BDI score for an independent validation dataset with 439 participants from 4 different imaging sites. Finally, we found 3 common functional connections between those related to depression symptoms and those related to MDD diagnosis. These findings contribute to a deeper understanding of the neural circuitry of depressive symptoms in MDD, a hetero-symptomatic population, revealing the neural basis of MDD.