- 1Department of Psychiatry, National Clinical Research Center for Mental Disorders, The Second Xiangya Hospital of Central South University, Changsha, China
- 2Department of Radiology, Tianjin Medical University General Hospital, Tianjin, China
- 3Department of Radiology, The Second Xiangya Hospital of Central South University, Changsha, China
- 4Hunan Key Laboratory of Children's Psychological Development and Brain Cognitive Science, Changsha, China
- 5Department of Psychiatry, The First Affiliated Hospital of China Medical University, Shenyang, China
- 6Department of Psychiatry, The Third People's Hospital of Foshan, Foshan, China
Background: Melancholic depression has been viewed as one severe subtype of major depressive disorder (MDD). However, it is unclear whether melancholic depression has distinct changes in brain imaging. We aimed to explore specific or distinctive alterations in melancholic MDD and whether the alterations could be used to separate melancholic MDD from non-melancholic MDD or healthy controls.
Materials and Methods: Thirty-one outpatients with melancholic MDD and thirty-three outpatients with non-melancholic MDD and thirty-two age- and gender-matched healthy controls were recruited. All participants were scanned by resting-state functional magnetic resonance imaging (fMRI). Imaging data were analyzed with the regional homogeneity (ReHo) and support vector machine (SVM) methods.
Results: Melancholic MDD patients exhibited lower ReHo in the right superior occipital gyrus/middle occipital gyrus than non-melancholic MDD patients and healthy controls. Merely for non-melancholic MDD patients, decreased ReHo in the right middle frontal gyrus was negatively correlated with the total HRSD-17 scores. SVM analysis results showed that a combination of abnormal ReHo in the right fusiform gyrus/cerebellum Crus I and the right superior occipital gyrus/middle occipital gyrus exhibited the highest accuracy of 83.05% (49/59), with a sensitivity of 90.32% (28/31), and a specificity of 75.00% (21/28) for discriminating patients with melancholic MDD from patients with non-melancholic MDD. And a combination of abnormal ReHo in the right fusiform gyrus/cerebellum VI and left postcentral gyrus/precentral gyrus exhibited the highest accuracy of 98.41% (62/63), with a sensitivity of 96.77% (30/31), and a specificity of 100.00%(32/32) for separating patients with melancholic MDD from healthy controls.
Conclusion: Our findings showed the distinctive ReHo pattern in patients with melancholic MDD and found brain area that may be associated with the pathophysiology of non-melancholic MDD. Potential imaging markers for discriminating melancholic MDD from non-melancholic MDD or healthy controls were reported.
Introduction
Major depressive disorder (MDD), a common, recurrent and disabling psychiatric disorder, severely leads to impaired psychosocial function and reduced quality of life (1). Melancholic depression has been viewed as one severe subtype of MDD. MDD is characterized by a persistent depressed mood, anxiety and dysphoria, psychomotor changes, alterations of motivation and social behavior, and sleep abnormalities (2). Despite of the common features of MDD, melancholic depression is characterized by pervasive anhedonia, unremitting apprehension and morbid statements, blunted emotional response, non-reactive mood, retardation, spontaneous agitation, reduced concentration, and impaired working memory (3). Although there is no consensus on whether melancholic depression is a categorically separate psychiatric disorder or a dimensionally severe subtype of MDD (4, 5), poorer cognitive performance was observed in melancholic MDD than non-melancholic MDD, particularly in visual learning and executive functions, like attention/working memory, reasoning/problem solving and processing speed (6, 7). The periods of cognitive recovery need longer time in melancholic MDD patients than the non-melancholic individuals (8). Besides, patients with melancholic MDD showed higher risk of relapse than patients with non-melancholic MDD (9). Thus, the diagnosis of melancholic MDD has a superior predictive validity for its treatment and prognosis (3).
Patients with MDD have been reported to show structural and functional abnormalities in brain imaging, like abnormal network homogeneity (NH) (10–12), functional connectivity (FC) (13–15), voxel-mirrored homotopic connectivity (VMHC) (16–18), low-frequency fluctuations (ALFF) (19, 20), voxel-based morphometry (VBM) (21, 22) and functional network connectivity (FNC) (23, 24). Meanwhile, structural or functional alterations in melancholic patients have been revealed in many previous studies, such as gray/white matter alterations (25), asymmetrical enlargement of CSF space in the Sylvian fissure region (26), decreased effective connectivity of cortical systems involved in attention and interoception compared with non-melancholic MDD patients and healthy controls (27), lower mean fractional anisotropy (FA) in the right ventral tegmental area-lateral orbitofrontal cortex (VTA-LOFC) connection than patients with non-melancholic MDD (28). Increased FA and decreased radial diffusivity (RD) were also identified in the right anterior limb of the internal capsule in melancholic MDD compared with health controls (29). A previous functional connectivity study reported that melancholic medication-free remitted MDD (rMDD) patients exhibited decreased subgenual cingulate cortex (SCC) connectivity with the parahippocampal gyrus and amygdala compared with non-melancholic rMDD patients (30). However, it is unclear which brain region is abnormal when one brain region shows abnormal functional connections to other brain regions.
Regional homogeneity (ReHo) method reflects the temporal homogeneity of the regional blood oxygen level-dependent (BOLD) signal that enables it to reveal the temporal homogeneity of neural activity. A previous Reho study believed that voxels within a functional brain region were more temporally homogeneous involving in a certain condition (31). So abnormal ReHo may be due to unsynchronized regional neural activity, and the ReHo method can be applied to find abnormal neural activity in the whole brain regions. Kendall's coefficient of concordance (KCC) can measure the similarity of multiple time series. KCC was applied to measure ReHo of the time series of one given voxel with those of its nearest voxels in a voxel-wise way (32, 33). Many previous studies have applied the ReHo method to study psychiatric disorders such as depression (34, 35), schizophrenia (36) and ADHD (37). However, the questions of whether melancholic MDD has specific or distinctive alterations relative to non-melancholic MDD in regional homogeneity and whether abnormal regional homogeneity could be used to distinguish melancholic MDD and non-melancholic MDD as well as healthy controls still remain unknown.
To address these questions, we applied the ReHo method to examine patients with melancholic MDD and patients with non-melancholic MDD. We hypothesized that abnormal ReHo would be observed in certain regions in patients with melancholic MDD at rest and abnormal regions might be distinct areas that could be used to discriminating melancholic MDD from non-melancholic MDD or healthy controls. To test our hypotheses, we compared ReHo of whole brain regions across patients with melancholic MDD, patients with non-melancholic MDD and healthy controls.
Methods
Participants
We recruited a total of 31 outpatients with melancholic MDD and 33 outpatients with non-melancholic MDD of the Second Xiangya Hospital of Central South University. They were all aged from 18 to 45 years old. The data were collected from May 4, 2014 to December 30, 2016. The diagnosis was independently ascertained by two psychiatrists, according to the Diagnostic and Statistical Manual of Mental Disorders, Fourth Edition (DSM-4). All the patients met the following inclusion criteria: (1) first major depressive episode with Hamilton Rating Scale for Depression (HRSD-17) (38) total scores ≥17 and HRSD item 1 score ≥ 2; (2) illness duration <12 months; (3) no history of antipsychotics and electroconvulsive therapy. The criteria for the features of melancholic MDD in DSM-4 were required as follows: (1) at least one of the following symptoms occurs during the most severe period of the current episode: loss of pleasure in all or almost all activities (pervasive anhedonia); lack of mood reactivity to usually pleasurable stimuli (does not feel much better, even temporally, when something good happens) (non-reactive mood); (2) three (or more) of the following symptoms: distinct quality of depressed mood (i.e., depressive mood experienced is qualitatively different from the feeling experienced when the loved one dies. Showing extreme despondency, despair, and/or morose mood or alleged empty mood); depression is regularly worse in the morning; early morning awakening (i.e., at least 2 h earlier than usual awakening) (HRSD item 6 ≥ 1); marked psychomotor agitation or retardation (HRSD items 8 or 9 ≥ 2); significant anorexia or weight loss (HRSD items 12 or 16 = 2); and excessive or inappropriate guilt (HRSD item 2 ≥ 2). The anhedonic states was assessed by using the Chinese version of Snaith-Hamilton Pleasure Scale (SHAPS-C) (39), the higher the score, the more severe the anhedonia; The Chinese version of Temporal Experience of Pleasure Scale (TEPS) was applied to capture the level of anticipatory and consummatory facets of pleasure, the lower the score, the greater the anhedonia. Non-melancholic MDD group consisted of patients who did not meet these criteria.
A total of 32 age- and gender-matched healthy controls were recruited from the community. Healthy controls would be excluded if they were related to the patients or had a family history of mental illness, especially their first-degree relatives. Besides, they would be ruled out if they had any neurological disorders, substance abuse or psychosis symptoms.
Exclusion criteria for all participants as follows: (1) other psychiatric disorders meeting DSM-4 diagnostic criteria; (2) any history of neurological disorders, severe physical illnesses, and substance abuse; (3) pregnancy; (4) abnormal cerebral structure after initial MRI scan; (5) any contraindications for MRI scan.
All participants were right-handed and Han Chinese with at least 9 years of education. The HRSD-17 was applied to determine the severity of depression; The Beck anxiety inventory (BAI) was administered to evaluate anxiety state. All participants were evaluated by the HRSD-17, SHAPS-C and BAI, and all patients completed TEPS evaluation.
The study was approved by the Medical Research Ethics Committee of the Second Xiangya Hospital of Central South University, China. The study was performed in line with the Helsinki Declaration. Each participant has completed an informed consent prior to enrollment.
Image Acquisition
Resting-state MRI data were obtained with a 3.0 T Siemens scanner (Germany) at the Second Xiangya Hospital of Central South University. All subjects were asked to lay supine, remain still, close their eyes and stay awake. Soft earplugs and foam pads were used to reduce scanner noise and head movement. Using the echo planar imaging (EPI) sequence, the resting-state functional images were obtained through the following parameters: repetition time/echo time (TR/TE) 2,500/25 ms, 39 slices, 64*64 matrix, 90°flip angle, 24 cm field of view, 3.5 mm slice thickness, no gap, and 200 volumes lasting for 500 s.
Data Preprocessing
Data preprocessing was conducted in Matlab (Mathworks) by using Data Processing Assistant for Resting-State fMRI (DPARSF) (40). Considering the instability of the initial MRI signal and the influences of subjects' adaptation time, the first 10 images were discarded. None of the participants had more than 2 mm of maximum displacement in x, y, or z axis or more than 2 of angular rotation after correction for slice timing and head motion. Then, the corrected imaging data were spatially normalized to the MNI space and got resampled to 3 × 3 × 3 mm3. After that, the fMRI data were temporally band-pass filtered (0.01–0.08 Hz) and linearly detrended. Several spurious covariates were also removed, such as the signal from the white matter-centered region and the ventricular seed-based region of interest (ROI) (Supplementary Figure 1) as well as the 24-head motion parameters obtained by rigid body correction. According to a previous study (41), global signal was retained when the resting-state functional connectivity (FC) data was pre-processed.
ReHo Analysis
Regional homogeneity (ReHo) analysis was conducted by in-house software REST (http://www.resting-fmri.sourceforge.net). Individual ReHo maps were generated by calculating the KCC of the time series of a given voxel with those of its nearest voxels (26 voxels) in a voxel-wise analysis. The formula for calculating the KCC value has been elaborated in a previous study (31). In order to reduce the influence of individual variation in the KCC value, we normalized the ReHo maps by dividing the KCC of each voxel by the averaged KCC of the whole brain. The resulting imaging data were then spatially smoothed with a Gaussian kernel of 4 mm full-width at half-maximum.
Statistical Analyses
Analysis of variance (ANOVA) was conducted to analyze group differences in age, years of education, HRSD-17 scores, BAI scores and SHAPS-C scores across the three groups by using SPSS19.0 (LSD between two group comparison), and we applied two-sample t-test to analyze whether there were group differences in the illness duration and the TEPS scores between melancholic MDD group and non-melancholic MDD group. A Chi-square test was performed to describe gender distributions. The significance level was set at p < 0.05.
Analysis of covariance (ANCOVA) was performed on individual whole brain ReHo maps across the three groups in a voxel-by-voxel manner to identify the group differences, followed by post-hoc t-tests. Age, years of education and framewise displacement were used as covariates. The results were FDR (false discovery rate) corrected at p < 0.05.
Correlation Analyses
The voxel-based mean ReHo values were extracted from the brain regions with abnormal ReHo by in-house software REST (http://www.resting-fmri.sourceforge.net). The correlations between abnormal ReHo and HRSD-17, BAI, SHARS, TEPS scores were determined using Pearson's correlation analyses with a threshold of Benjamini-Hochberg corrected p < 0.05.
Classification Analyses
In the LIBSVM software package (http://www.csie.ntu.edu.tw/~cjlin/libsvm/) in MATLAB, support vector machines (SVM) was applied to test the ability to separate melancholic MDD from non-melancholic MDD, melancholic MDD from healthy controls, and non-melancholic MDD from healthy controls by using the identified ReHo values in abnormal brain regions. The “leave-one-out” method was applied in the study.
Results
Demographic Characteristics and Clinical Information
Five patients with non-melancholic MDD were excluded due to excessive head movement. Thus, 31 patients with melancholic MDD, 28 patients with non-melancholic MDD and 32 healthy controls were finally included in the analyses. Age and gender did not significantly differ across the three groups, and no illness duration difference was observed between the two patient groups. Significant differences were observed across the three groups in years of education, HRSD-17 scores, BAI scores and SHAPS-C scores. The melancholic MDD group (p = 0.001) and the healthy control group (p = 0.01) showed significantly higher education level than that of the non-melancholic MDD group. The HRSD-17, BAI and SHAPS-C scores of the melancholic MDD group (p < 0.001) group and non-melancholic MDD group (p < 0.001) were significantly higher than those of the healthy control group. The non-melancholic group exhibited significantly lower BAI (p = 0.024) and SHAPS-C scores (p < 0.001) than those of the melancholic group, while no significant difference in HRSD-17 scores was observed between the two patient groups (p = 0.743). Part of the TEPS scores showed significant differences between the two patient groups. The TEPS total scores (p = 0.002), TEPS Abstract Anticipatory scores (p = 0.001) and TEPS Contextual Anticipatory scores (p = 0.001) of the melancholic group were significantly lower than those of the non-melancholic group. The TEPS Abstract Consummatory scores (p = 0.117) and TEPS Contextual Consummatory (p = 0.074) did not significantly differ between the two patient groups. More details of demographic and clinical data are presented in Table 1.
ReHo Differences Between Groups
By using ANCOVA, the brain regions mainly in the right middle frontal gyrus, right fusiform gyrus/cerebellum Crus I, bilateral cerebellum Crus II, right fusiform gyrus/cerebellum VI, left middle occipital gyrus/inferior occipital gyrus, bilateral superior occipital gyrus/middle occipital gyrus, bilateral postcentral gyrus/precentral gyrus and right middle temporal gyrus showed significant differences of ReHo values across the three groups (Figure 1). We reanalyzed ANCONA adding gender as a covariate since previous studies (42–44) have reported the gender difference in depression, which is more prone to females. Similar results were obtained (Supplementary Table 1).
Figure 1. Brain regions showing group differences of ReHo values across three groups. Color bar indicates F values from ANCOVA (age, years of education and framewise displacement as covariates). ReHo, regional homogeneity; ANCOVA, analysis of covariance.
Melancholic vs. Non-melancholic
Compared with non-melancholic group, melancholic group showed increased ReHo in the right middle frontal gyrus and decreased ReHo in the right fusiform gyrus/cerebellum Crus I and the right superior occipital gyrus/middle occipital gyrus (Figure 2, Table 2).
Figure 2. Statistical map depicts higher and lower ReHo of melancholic MDD patients compared with non-melancholic MDD patients. The threshold was set at p < 0.05. Blue denotes lower ReHo and red denotes higher ReHo. Color bar indicates T values from post-hoc t-tests. L, left side; R, right side; ReHo, regional homogeneity.
Melancholic vs. Healthy Controls
Compared with healthy control group, melancholic group exhibited increased ReHo in the bilateral cerebellum Crus II and decreased ReHo in the right fusiform gyrus/cerebellum VI, left middle occipital gyrus/inferior occipital gyrus, bilateral superior occipital gyrus/middle occipital gyrus, and bilateral postcentral gyrus/precentral gyrus (Figure 3, Table 2).
Figure 3. Statistical map depicts higher and lower ReHo of melancholic MDD patients compared with healthy controls. The threshold was set at p < 0.05. Blue denotes lower ReHo and red denotes higher ReHo. Color bar indicates T values from post-hoc t-tests. L, left side; R, right side; ReHo, regional homogeneity.
Non-melancholic vs. Healthy Controls
Compared with healthy control group, increased ReHo in the right middle temporal gyrus was found in non-melancholic group (Figure 4, Table 2).
Figure 4. Statistical map depicts higher ReHo of non-melancholic MDD patients compared with healthy controls. The threshold was set at p < 0.05. Red denotes higher ReHo. Color bar indicates T values from post-hoc t-tests. L, left side; R, right side; ReHo, regional homogeneity.
Correlations Between ReHo and Clinical Characteristics
No significant correlation between abnormal ReHo and clinical features was found in melancholic MDD patients. For non-melancholic MDD patients, decreased ReHo in the right middle frontal gyrus was negatively correlated with the total HRSD-17 scores (r = −0.527, p = 0.004, Benjamini-Hochberg correction p = 0.016).
Discrimination Patients With Melancholic MDD From Other Two Groups
Discriminating Patients With Melancholic MDD From Non-melancholic MDD
The SVM results showed that a combination of abnormal ReHo in the right fusiform gyrus/cerebellum Crus I and right superior occipital gyrus/middle occipital gyrus exhibited the highest accuracy of 83.05% (49/59), with a sensitivity of 90.32% (28/31), and a specificity of 75.00%(21/28) for discriminating patients with melancholic MDD from patients with non-melancholic MDD (Figures 5, 6).
Figure 5. Visualization of classifications through support vector machine (SVM) using the combination of ReHo values in the right fusiform gyrus/cerebellum Crus I and the right superior occipital gyrus/middle occipital gyrus to discriminate melancholic MDD and non-melancholic MDD. Left: SVM parameters result of 3D view. Right: dimension 1 and dimension 2 represent the ReHo values in the right fusiform gyrus/cerebellum Crus I and the right superior occipital gyrus/middle occipital gyrus, respectively. Red crosses represent non-melancholic MDD patients, and green crosses represent the melancholic MDD patients. MDD, major depression disorder.
Figure 6. The accuracy of using abnormal ReHo in different brain regions to classify patients with melancholic and non-melancholic MDD. 1, Right Middle Frontal Gyrus; 2, Right Fusiform Gyrus/Cerebellum Crus I; 3, Right Superior Occipital Gyrus/Middle Occipital Gyrus; 12, Right Middle Frontal Gyrus and Right Fusiform Gyrus/Cerebellum Crus I; 13, Right Middle Frontal Gyrus and Right Superior Occipital Gyrus/Middle Occipital Gyrus; 23, Right Fusiform Gyrus/Cerebellum Crus I and Right Superior Occipital Gyrus/Middle Occipital Gyrus; 123, Right Middle Frontal Gyrus and Right Fusiform Gyrus/Cerebellum Crus I and Right Superior Occipital Gyrus/Middle Occipital Gyrus.
Discriminating Patients With Melancholic MDD From Healthy Controls
The SVM results showed that a combination of abnormal ReHo in the right fusiform gyrus/cerebellum VI and left postcentral gyrus/precentral gyrus showed the highest accuracy of 98.41% (62/63), with a sensitivity of 96.77% (30/31), and a specificity of 100.00% (32/32) for distinguishing patients with melancholic MDD from healthy controls (Figure 7).
Figure 7. Visualization of classifications through support vector machine (SVM) using the combination of ReHo values in the right fusiform gyrus/cerebellum VI and left postcentral gyrus/precentral gyrus to discriminate melancholic MDD and from controls. Left: SVM parameters result of 3D view. Right: dimension 1 and dimension 2 represent the ReHo values in the right fusiform gyrus/cerebellum VI and left postcentral gyrus/precentral gyrus, respectively. Red crosses represent healthy controls, and green crosses represent the melancholic MDD patients. MDD, major depression disorder.
Discriminating Patients With Non-melancholic MDD From Healthy Controls
The SVM results showed that abnormal ReHo in the right middle temporal gyrus had an accuracy of 65% (39/60), with a sensitivity of 50% (14/28), and a specificity of 78.13% (28/32) for separating the non-melancholic MDD patients from healthy controls (Figure 8).
Figure 8. Visualization of classifications through support vector machine (SVM) using the increased ReHo values in the right middle temporal gyrus to discriminate non-melancholic MDD patients from healthy controls. Left: SVM parameters result of 3D view. Right: Classified map of the ReHo values in the right middle temporal gyrus. MDD, major depressive disorder.
Discussion
In the present study, we investigated the whole brain regional homogeneity in patients with melancholic MDD and non-melancholic MDD at rest. We observed that melancholic MDD patients exhibited decreased ReHo in the right superior occipital gyrus/middle occipital gyrus compared with non-melancholic MDD patients and healthy controls. Abnormal ReHo in the right middle temporal gyrus may play an important role in the pathophysiology of non-melancholic MDD. Merely for non-melancholic MDD patients, decreased ReHo in the right middle frontal gyrus was negatively correlated with the total HRSD-17 scores. In addition, the SVM analysis results showed that a combination of abnormal ReHo in the right fusiform gyrus/cerebellum Crus I and right superior occipital gyrus/middle occipital gyrus might be a potential imaging marker to separate melancholic MDD patients from non-melancholic MDD patients, and a combination of abnormal ReHo in the right fusiform gyrus/cerebellum VI and left postcentral gyrus/precentral gyrus might be a potential imaging marker to separate melancholic MDD patients from healthy controls.
Occipital cortex involves in consolidating information into visual working memory (45). A previous study has observed that the occipital bending in MDD patients was more common than healthy controls and enlargement of ventricle may exacerbate the natural curvature of the occipital regions (46). Another study has revealed the asymmetrical enlargement of CSF space in the Sylvian fissure region in melancholic MDD patients (26). That issue might indicate that abnormal structural or functional changes may occur in the occipital regions in melancholic MDD. Consistent with these studies, we found that melancholic MDD patients exhibited decreased ReHo in the right superior occipital gyrus/middle occipital gyrus compared with non-melancholic MDD patients and healthy controls. And that point can interpret working memory impairment in melancholic MDD (45). Hence, we can speculate that decreased ReHo in the right superior occipital gyrus/middle occipital gyrus might be stable and distinctive neurobiological characteristic of melancholic MDD.
The temporal lobe has been expounded to get involved in emotional regulation, process of memory, and social cognition (47, 48). Many previous studies have observed abnormal brain activity or functional connectivity in the temporal areas in MDD patients. Decreased NH was found in melancholic MDD in the right middle temporal gyrus and temporal pole (MTG/TP) compared with healthy controls (12). MDD patients showed lower NH in the right inferior temporal gyrus (ITG) than healthy controls (10). Remitted geriatric depression exhibited decreased ReHo in the right superior and middle temporal gyrus (49). MDD showed significantly increased ReHo in the right inferior temporal gyrus (50). In the present study, we found increased ReHo of the right middle temporal gyrus in non-melancholic MDD patients compared with healthy controls. Although we could not directly compare these above mentioned studies because of different composition of MDD patients, these studies indicated abnormal brain activity of temporal gyrus in non-melancholic MDD. We did not find abnormal ReHo in the temporal gyrus in patients with melancholic MDD compared with both healthy controls and non-melancholic MDD patients. We suspected that the differences may be that the non-melancholic patients were characterized by other main characteristics (like anxiety, atypical increased sleep, somatic symptoms, etc.) rather than melancholic traits. A voxel-based morphometry (VBM) study reported that anxiety was positive correlated with VBM abnormalities in the middle temporal gyrus (51). A ReHo study of MDD patients with somatic symptoms (somatic depression, SD) reported that abnormal ReHo in the frontal and temporal regions may be involved in the neural basis of SD (52). Some previous studies revealed that the structural alterations in the temporal gyrus might be associated with suicide attempt (53, 54). Although we did not further classify the subtypes of non-melancholic patients because of small sample size, the present study reminded us to examine this issue in the future study. Thus, we speculated that abnormal ReHo in the right middle temporal gyrus may play an important role in the pathophysiology of non-melancholic MDD.
Frontal lobe is associated with cognitive control (55). Abnormal structural and functional alterations in the frontal gyrus in patients MDD were reported in numerous previous studies (56–58). Melancholic MDD was reported to show cognitive impairment (like poorer processing speed and verbal fluency) than HCs and atypical depression (59). A distinct and poorer cognitive performance was also reported in melancholic MDD compared with non-melancholic MDD (6, 7). In present study, non-melancholic group showed decreased ReHo in the right middle frontal gyrus compared with melancholic group and decreased ReHo in the right middle frontal gyrus was negatively correlated with the total HRSD-17 scores in non-melancholic MDD patients. The significant differences in the frontal gyrus between two patient groups might interpret poorer cognitive performance in melancholic MDD compared with non-melancholic MDD (6).
Support vector machine (SVM) has been widely used in biomedical applications for diagnoses of psychiatric disorders like major depression (60, 61), schizophrenia (62) and others (63). The sensitivity or specificity of qualified diagnostic indicators should not be <60%, and more than 70% is conducive to the establishment of diagnostic indicators (64, 65). The SVM analysis results showed that a combination of abnormal ReHo in the right fusiform gyrus/cerebellum Crus I and the right superior occipital gyrus/middle occipital gyrus exhibited the highest accuracy of 83.05% (49/59), with a sensitivity of 90.32% (28/31), and a specificity of 75.00% (21/28) for discriminating patients with melancholic MDD from patients with non-melancholic MDD. Thus, the combination of ReHo in the right fusiform gyrus/cerebellum Crus I and the right superior occipital gyrus/middle occipital gyrus might be employed as a potential imaging marker for discriminating patients with melancholic MDD from patients with non-melancholic MDD. The melancholic MDD is characterized by pervasive anhedonia and non-reactive mood, etc. (as we mentioned in the introduction part) (3), and shows poorer cognitive performance than non-melancholic MDD (6, 7). In the present study, non-melancholic patients were composited of other subtype of MDD instead of melancholic subtype, and were characterized by other characteristics (such as anxiety, atypical increased sleep, somatic symptoms, etc.) rather than melancholic traits. This may be the main reason for the apparent difference between the two subtypes of MDD. Besides, to increase the accuracy, we could apply multimodal brain imaging data to discriminate the two subtypes of MDD in future studies according to a previous study (66), such as combining two types of functional data, or combining both anatomical and functional data. In the present study, the SVM results showed that a combination of abnormal ReHo in the right fusiform gyrus/cerebellum VI and left postcentral gyrus/precentral gyrus exhibited the highest accuracy of 98.41% (62/63), with a sensitivity of 96.77% (30/31), and a specificity of 100.00% (32/32) for separating patients with melancholic MDD from healthy controls. Thus, we speculated that the combination of abnormal ReHo in the right fusiform gyrus/cerebellum VI and left postcentral gyrus/precentral gyrus might be a potential imaging marker for separating melancholic MDD patients from healthy controls. When we used abnormal ReHo in the right middle temporal gyrus to separate non-melancholic MDD from healthy controls, we obtained a moderate accuracy of 65% (39/60), with a sensitivity of 50% (14/28) and a specificity of 78.13% (28/32), which might be weak for the classification due to the lower sensitivity. As we mentioned earlier, the different composition of non-melancholic patients may have an impact on this issue, and we can further classify these patients into different subtypes to explore them in future studies.
Some limitations exist in our study. First, our sample size is small. Second, we recruited both melancholic and non-melancholic MDD patients but we did not further classify the non-melancholic MDD patients due to small sample size. Therefore, we were unable to further explore the neurobiological differences between different subtypes of non-melancholic MDD.
Conclusion
Our findings showed the distinctive ReHo pattern in patients with melancholic MDD. Decreased ReHo in the right superior occipital gyrus/middle occipital gyrus might be a stable and distinctive neurobiological characteristic of melancholic MDD. Abnormal ReHo in the right middle temporal gyrus may play an important role in the pathophysiology of non-melancholic MDD. In addition, a combination of the ReHo values in the right fusiform gyrus/cerebellum Crus I and the right superior occipital gyrus/middle occipital gyrus might be employed as a potential imaging marker for separating patients with melancholic MDD from patients with non-melancholic MDD. And a combination of abnormal ReHo values in the right fusiform gyrus/cerebellum VI and left postcentral gyrus/precentral gyrus might be a potential imaging marker for separating melancholic MDD patients from healthy controls.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics Statement
The studies involving human participants were reviewed and approved by the Medical Research Ethics Committee of the Second Xiangya Hospital of Central South University, China. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
All authors listed have made a substantial, direct and intellectual contribution to the work, and approved it for publication.
Funding
This study was supported by grants from the National Key R&D Program of China (Grant no. 2016YFC1307100), National Natural Science Foundation of China (Grant no. 81771447), Natural Science Foundation of Hunan (Grant no. 2020JJ4784), Natural Science Foundation of Tianjin (Grant no. 18JCQNJC10900), Hunan Key Laboratory of Children's Psychological Development and Brain Cognitive Science (Grant no. 2019TP1032) and Hunan Provincial Innovation Foundation for Postgraduate (Grant no. CX20190159).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We thank all participants who served as research participants.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpsyt.2021.618805/full#supplementary-material
Supplementary Figure 1. The ventricular seed-based region of interest (ROI) in data preprocessing.
Supplementary Table 1. Significant ReHo differences across three groups*.
References
2. American Psychiatric Association. DSM-IV: Diagnostic and Statistical Manual of Mental Disorders 4th ed. Washington, DC: Am Psychiatric Assoc Press (1994).
3. Parker G, Fink M, Shorter E, Taylor MA, Akiskal H, Berrios G, et al. Issues for DSM-5: whither melancholia? The case for its classification as a distinct mood disorder. Am J Psychiatry. (2010) 167:745–7. doi: 10.1176/appi.ajp.2010.09101525
4. Parker G, Tavella G, Hadzi-Pavlovic D. Identifying and differentiating melancholic depression in a non-clinical sample. J Affect Disord. (2019) 243:194–200. doi: 10.1016/j.jad.2018.09.024
5. Day CV, Williams LM. Finding a biosignature for melancholic depression. Expert Rev Neurother. (2012) 12:835–47. doi: 10.1586/ern.12.72
6. Zaninotto L, Solmi M, Veronese N, Guglielmo R, Ioime L, Camardese G, et al. A meta-analysis of cognitive performance in melancholic versus non-melancholic unipolar depression. J Affect Disord. (2016) 201:15–24. doi: 10.1016/j.jad.2016.04.039
7. Bosaipo NB, Foss MP, Young AH, Juruena MF. Neuropsychological changes in melancholic and atypical depression: a systematic review. Neurosci Biobehav Rev. (2017) 73:309–25. doi: 10.1016/j.neubiorev.2016.12.014
8. Withall A, Harris LM, Cumming SR. A longitudinal study of cognitive function in melancholic and non-melancholic subtypes of major depressive disorder. J Affect Disord. (2010) 123:150–7. doi: 10.1016/j.jad.2009.07.012
9. Gupta R, Mirza T, Majeed MH, Seemüller F, Moeller HJ. Survival of melancholia: a retrospective study of patients with depressive disorders. East Asian Arch Psychiatry. (2020) 30:39–43. doi: 10.12809/eaap1884
10. Guo W, Liu F, Zhang J, Zhang Z, Yu L, Liu J, et al. Abnormal default-mode network homogeneity in first-episode, drug-naive major depressive disorder. PLoS ONE. (2014) 9:e91102. doi: 10.1371/journal.pone.0091102
11. Guo W, Liu F, Yu M, Zhang J, Zhang Z, Liu J, et al. Decreased regional activity and network homogeneity of the fronto-limbic network at rest in drug-naive major depressive disorder. Aust N Z J Psychiatry. (2015) 49:550–6. doi: 10.1177/0004867415577978
12. Cui X, Guo W, Wang Y, Yang TX, Yang XH, Wang Y, et al. Aberrant default mode network homogeneity in patients with first-episode treatment-naive melancholic depression. Int J Psychophysiol. (2017) 112:46–51. doi: 10.1016/j.ijpsycho.2016.12.005
13. Chen Y, Wang C, Zhu X, Tan Y, Zhong Y. Aberrant connectivity within the default mode network in first-episode, treatment-naive major depressive disorder. J Affect Disord. (2015) 183:49–56. doi: 10.1016/j.jad.2015.04.052
14. Kaiser RH, Andrews-Hanna JR, Wager TD, Pizzagalli DA. Large-scale network dysfunction in major depressive disorder: a meta-analysis of resting-state functional connectivity. JAMA Psychiatry. (2015) 72:603–11. doi: 10.1001/jamapsychiatry.2015.0071
15. Pan PM, Sato JR, Salum GA, Rohde LA, Gadelha A, Zugman A, et al. Ventral striatum functional connectivity as a predictor of adolescent depressive disorder in a longitudinal community-based sample. Am J Psychiatry. (2017) 174:1112–9. doi: 10.1176/appi.ajp.2017.17040430
16. Guo W, Liu F, Dai Y, Jiang M, Zhang J, Yu L, et al. Decreased interhemispheric resting-state functional connectivity in first-episode, drug-naive major depressive disorder. Prog Neuropsychopharmacol Biol Psychiatry. (2013) 41:24–9. doi: 10.1016/j.pnpbp.2012.11.003
17. Hermesdorf M, Sundermann B, Feder S, Schwindt W, Minnerup J, Arolt V, et al. Major depressive disorder: findings of reduced homotopic connectivity and investigation of underlying structural mechanisms. Hum Brain Mapp. (2016) 37:1209–17. doi: 10.1002/hbm.23097
18. Lai CH, Wu YT. Decreased inter-hemispheric connectivity in anterior sub-network of default mode network and cerebellum: significant findings in major depressive disorder. Int J Neuropsychopharmacol. (2014) 17:1935–42. doi: 10.1017/S1461145714000947
19. Guo WB, Liu F, Xue ZM, Xu XJ, Wu RR, Ma CQ, et al. Alterations of the amplitude of low-frequency fluctuations in treatment-resistant and treatment-response depression: a resting-state fMRI study. Prog Neuropsychopharmacol Biol Psychiatry. (2012) 37:153–60. doi: 10.1016/j.pnpbp.2012.01.011
20. Li W, Chen Z, Wu M, Zhu H, Gu L, Zhao Y, et al. Characterization of brain blood flow and the amplitude of low-frequency fluctuations in major depressive disorder: a multimodal meta-analysis. J Affect Disord. (2017) 210:303–11. doi: 10.1016/j.jad.2016.12.032
21. Bora E, Fornito A, Pantelis C, Yücel M. Gray matter abnormalities in major depressive disorder: a meta-analysis of voxel based morphometry studies. J Affect Disord. (2012) 138:9–18. doi: 10.1016/j.jad.2011.03.049
22. Gray JP, Müller VI, Eickhoff SB, Fox PT. Multimodal abnormalities of brain structure and function in major depressive disorder: a meta-analysis of neuroimaging studies. Am J Psychiatry. (2020) 177:422–34. doi: 10.1176/appi.ajp.2019.19050560
23. Zhi D, Calhoun VD, Lv L, Ma X, Ke Q, Fu Z, et al. Aberrant dynamic functional network connectivity and graph properties in major depressive disorder. Front Psychiatry. (2018) 9:339. doi: 10.3389/fpsyt.2018.00339
24. He H, Yu Q, Du Y, Vergara V, Victor TA, Drevets WC, et al. Resting-state functional network connectivity in prefrontal regions differs between unmedicated patients with bipolar and major depressive disorders. J Affect Disord. (2016) 190:483–93. doi: 10.1016/j.jad.2015.10.042
25. Soriano-Mas C, Hernandez-Ribas R, Pujol J, Urretavizcaya M, Deus J, Harrison BJ, et al. Cross-sectional and longitudinal assessment of structural brain alterations in melancholic depression. Biol Psychiatry. (2011) 69:318–25. doi: 10.1016/j.biopsych.2010.07.029
26. Pujol J, Cardoner N, Benlloch L, Urretavizcaya M, Deus J, Losilla JM, et al. CSF spaces of the Sylvian fissure region in severe melancholic depression. Neuroimage. (2002) 15:103–6. doi: 10.1006/nimg.2001.0928
27. Hyett MP, Breakspear MJ, Friston KJ, Guo CC, Parker GB. Disrupted effective connectivity of cortical systems supporting attention and interoception in melancholia. JAMA Psychiatry. (2015) 72:350–8. doi: 10.1001/jamapsychiatry.2014.2490
28. Bracht T, Horn H, Strik W, Federspiel A, Schnell S, Hofle O, et al. White matter microstructure alterations of the medial forebrain bundle in melancholic depression. J Affect Disord. (2014) 155:186–93. doi: 10.1016/j.jad.2013.10.048
29. Hyett MP, Perry A, Breakspear M, Wen W, Parker GB. White matter alterations in the internal capsule and psychomotor impairment in melancholic depression. PLoS ONE. (2018) 13:e0195672. doi: 10.1371/journal.pone.0195672
30. Workman CI, Lythe KE, McKie S, Moll J, Gethin JA, Deakin JF, et al. Subgenual cingulate-amygdala functional disconnection and vulnerability to melancholic depression. Neuropsychopharmacology. (2016) 41:2082–90. doi: 10.1038/npp.2016.8
31. Zang Y, Jiang T, Lu Y, He Y, Tian L. Regional homogeneity approach to fMRI data analysis. Neuroimage. (2004) 22:394–400. doi: 10.1016/j.neuroimage.2003.12.030
33. Tononi G, McIntosh AR, Russell DP, Edelman GM. Functional clustering: identifying strongly interactive brain regions in neuroimaging data. Neuroimage. (1998) 7:133–49. doi: 10.1006/nimg.1997.0313
34. Spati J, Hanggi J, Doerig N, Ernst J, Sambataro F, Brakowski J, et al. Prefrontal thinning affects functional connectivity and regional homogeneity of the anterior cingulate cortex in depression. Neuropsychopharmacology. (2015) 40:1640–8. doi: 10.1038/npp.2015.8
35. Guo WB, Sun XL, Liu L, Xu Q, Wu RR, Liu ZN, et al. Disrupted regional homogeneity in treatment-resistant depression: a resting-state fMRI study. Prog Neuropsychopharmacol Biol Psychiatry. (2011) 35:1297–302. doi: 10.1016/j.pnpbp.2011.02.006
36. Wang S, Zhang Y, Lv L, Wu R, Fan X, Zhao J, et al. Abnormal regional homogeneity as a potential imaging biomarker for adolescent-onset schizophrenia: a resting-state fMRI study and support vector machine analysis. Schizophr Res. (2018) 192:179–84. doi: 10.1016/j.schres.2017.05.038
37. Shang CY, Lin HY, Tseng WY, Gau SS. A haplotype of the dopamine transporter gene modulates regional homogeneity, gray matter volume, and visual memory in children with attention-deficit/hyperactivity disorder. Psychol Med. (2018) 48:2530–40. doi: 10.1017/S0033291718000144
38. Hamilton M. Development of a rating scale for primary depressive illness. Br J Soc Clin Psychol. (1967) 6:278–96. doi: 10.1111/j.2044-8260.1967.tb00530.x
39. Liu WH, Wang LZ, Zhu YH, Li MH, Chan RC. Clinical utility of the Snaith-Hamilton-Pleasure scale in the Chinese settings. BMC Psychiatry. (2012) 12:184. doi: 10.1186/1471-244X-12-184
40. Chao-Gan Y, Yu-Feng Z. DPARSF: a MATLAB toolbox for “pipeline” data analysis of resting-state fMRI. Front Syst Neurosci. (2010) 4:13. doi: 10.3389/fnsys.2010.00013
41. Hahamy A, Calhoun V, Pearlson G, Harel M, Stern N, Attar F, et al. Save the global: global signal connectivity as a tool for studying clinical populations with functional magnetic resonance imaging. Brain Connect. (2014) 4:395–403. doi: 10.1089/brain.2014.0244
42. Briceño EM, Rapport LJ, Kassel MT, Bieliauskas LA, Zubieta JK, Weisenbach SL, et al. Age and gender modulate the neural circuitry supporting facial emotion processing in adults with major depressive disorder. Am J Geriatr Psychiatry. (2015) 23:304–13. doi: 10.1016/j.jagp.2014.05.007
43. Yang X, Peng Z, Ma X, Meng Y, Li M, Zhang J, et al. Sex differences in the clinical characteristics and brain gray matter volume alterations in unmedicated patients with major depressive disorder. Sci Rep. (2017) 7:2515. doi: 10.1038/s41598-017-02828-4
44. Young KD, Bodurka J, Drevets WC. Functional neuroimaging of sex differences in autobiographical memory recall in depression. Psychol Med. (2017) 47:2640–52. doi: 10.1017/S003329171700112X
45. Makovski T, Lavidor M. Stimulating occipital cortex enhances visual working memory consolidation. Behav Brain Res. (2014) 275:84–7. doi: 10.1016/j.bbr.2014.09.004
46. Maller JJ, Thomson RH, Rosenfeld JV, Anderson R, Daskalakis ZJ, Fitzgerald PB. Occipital bending in depression. Brain. (2014) 137(Pt 6):1830–7. doi: 10.1093/brain/awu072
47. Beauregard M, Paquette V, Levesque J. Dysfunction in the neural circuitry of emotional self-regulation in major depressive disorder. Neuroreport. (2006) 17:843–6. doi: 10.1097/01.wnr.0000220132.32091.9f
48. Gallagher HL, Frith. CD. Functional imaging of 'theory of mind'. Trends Cogn Sci. (2003) 7:77–83. doi: 10.1016/S1364-6613(02)00025-6
49. Yuan Y, Zhang Z, Bai F, Yu H, Shi Y, Qian Y, et al. Abnormal neural activity in the patients with remitted geriatric depression: a resting-state functional magnetic resonance imaging study. J Affect Disord. (2008) 111:145–52. doi: 10.1016/j.jad.2008.02.016
50. Guo WB, Liu F, Xue ZM, Yu Y, Ma CQ, Tan CL, et al. Abnormal neural activities in first-episode, treatment-naive, short-illness-duration, and treatment-response patients with major depressive disorder: a resting-state fMRI study. J Affect Disord. (2011) 135:326–31. doi: 10.1016/j.jad.2011.06.048
51. Besteher B, Gaser C, Langbein K, Dietzek M, Sauer H, Nenadić I. Effects of subclinical depression, anxiety and somatization on brain structure in healthy subjects. J Affect Disord. (2017) 215:111–7. doi: 10.1016/j.jad.2017.03.039
52. Geng J, Yan R, Shi J, Chen Y, Mo Z, Shao J, et al. Altered regional homogeneity in patients with somatic depression: a resting-state fMRI study. J Affect Disord. (2019) 246:498–505. doi: 10.1016/j.jad.2018.12.066
53. Pan LA, Ramos L, Segreti A, Brent DA, Phillips ML. Right superior temporal gyrus volume in adolescents with a history of suicide attempt. Br J Psychiatry. (2015) 206:339–40. doi: 10.1192/bjp.bp.114.151316
54. McLellan Q, Wilkes TC, Swansburg R, Jaworska N, Langevin LM, MacMaster FP. History of suicide attempt and right superior temporal gyrus volume in youth with treatment-resistant major depressive disorder. J Affect Disord. (2018) 239:291–4. doi: 10.1016/j.jad.2018.07.030
55. Badre D, Nee DE. Frontal Cortex and the hierarchical control of behavior. Trends Cogn Sci. (2018) 22:170–88. doi: 10.1016/j.tics.2017.11.005
56. Suh JS, Schneider MA, Minuzzi L, MacQueen GM, Strother SC, Kennedy SH, et al. Cortical thickness in major depressive disorder: a systematic review and meta-analysis. Prog Neuropsychopharmacol Biol Psychiatry. (2019) 88:287–302. doi: 10.1016/j.pnpbp.2018.08.008
57. Zaremba D, Dohm K, Redlich R, Grotegerd D, Strojny R, Meinert S, et al. Association of brain cortical changes with relapse in patients with major depressive disorder. JAMA Psychiatry. (2018) 75:484–92. doi: 10.1001/jamapsychiatry.2018.0123
58. Satterthwaite TD, Cook PA, Bruce SE, Conway C, Mikkelsen E, Satchell E, et al. Dimensional depression severity in women with major depression and post-traumatic stress disorder correlates with fronto-amygdalar hypoconnectivty. Mol Psychiatry. (2016) 21:894–902. doi: 10.1038/mp.2015.149
59. Lin K, Xu G, Lu W, Ouyang H, Dang Y, Lorenzo-Seva U, et al. Neuropsychological performance in melancholic, atypical and undifferentiated major depression during depressed and remitted states: a prospective longitudinal study. J Affect Disord. (2014) 168:184–91. doi: 10.1016/j.jad.2014.06.032
60. Yu JS, Xue AY, Redei EE, Bagheri N. A support vector machine model provides an accurate transcript-level-based diagnostic for major depressive disorder. Transl Psychiatry. (2016) 6:e931. doi: 10.1038/tp.2016.198
61. Shan X, Liao R, Ou Y, Ding Y, Liu F, Chen J, et al. Metacognitive training modulates default-mode network homogeneity during 8-week olanzapine treatment in patients with schizophrenia. Front Psychiatry. (2020) 11:234. doi: 10.3389/fpsyt.2020.00234
62. Cao B, Cho RY, Chen D, Xiu M, Wang L, Soares JC, et al. Treatment response prediction and individualized identification of first-episode drug-naïve schizophrenia using brain functional connectivity. Mol Psychiatry. (2018) 25:906–13. doi: 10.1038/s41380-018-0106-5
63. Orru G, Pettersson-Yeo W, Marquand AF, Sartori G, Mechelli A. Using support vector machine to identify imaging biomarkers of neurological and psychiatric disease: a critical review. Neurosci Biobehav Rev. (2012) 36:1140–52. doi: 10.1016/j.neubiorev.2012.01.004
64. Gong Q, Wu Q, Scarpazza C, Lui S, Jia Z, Marquand A, et al. Prognostic prediction of therapeutic response in depression using high-field MR imaging. Neuroimage. (2011) 55:1497–503. doi: 10.1016/j.neuroimage.2010.11.079
65. Swets JA. Measuring the accuracy of diagnostic systems. Science. (1988) 240:1285–93. doi: 10.1126/science.3287615
Keywords: melancholic depression, non-melancholic depression, regional homogeneity, magnetic resonance imaging, resting state
Citation: Yan M, He Y, Cui X, Liu F, Li H, Huang R, Tang Y, Chen J, Zhao J, Xie G and Guo W (2021) Disrupted Regional Homogeneity in Melancholic and Non-melancholic Major Depressive Disorder at Rest. Front. Psychiatry 12:618805. doi: 10.3389/fpsyt.2021.618805
Received: 18 October 2020; Accepted: 25 January 2021;
Published: 16 February 2021.
Edited by:
André Schmidt, University of Basel, SwitzerlandReviewed by:
Xiufeng Xu, The First Affiliated Hospital of Kunming Medical University, ChinaJing Sui, Chinese Academy of Sciences, China
Copyright © 2021 Yan, He, Cui, Liu, Li, Huang, Tang, Chen, Zhao, Xie and Guo. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Wenbin Guo, guowenbin76@csu.edu.cn; Guangrong Xie, xiegr2000@126.com
†These authors have contributed equally to this work