Skip to main content

Alterations of functional connectivity of the lateral habenula in subclinical depression and major depressive disorder

Abstract

Background

Major depressive disorder (MDD) is a common cause of disability and morbidity, affecting about 10% of the population worldwide. Subclinical depression (SD) can be understood as a precursor of MDD, and therefore provides an MDD risk indicator. The pathogenesis of MDD and SD in humans is still unclear, and the current diagnosis lacks accurate biomarkers and gold standards.

Methods

A total of 40 MDD, 34 SD, and 40 healthy control (HC) participants matched by age, gender, and education were included in this study. Resting-state functional magnetic resonance images (rs-fMRI) were used to analyze the functional connectivity (FC) of the posterior parietal thalamus (PPtha), which includes the lateral habenula, as the region of interest. Analysis of variance with the post hoc t-test test was performed to find significant differences in FC and clarify the variations in FC among the HC, SD, and MDD groups.

Results

Increased FC was observed between PPtha and the left inferior temporal gyrus (ITG) for MDD versus SD, and between PPtha and the right ITG for SD versus HC. Conversely, decreased FC was observed between PPtha and the right middle temporal gyrus (MTG) for MDD versus SD and MDD versus HC. The FC between PPtha and the middle frontal gyrus (MFG) in SD was higher than that in MDD and HC. Compared with the HC group, the FC of PPtha-ITG (left and right) increased in both the SD and MDD groups, PPtha-MTG (right) decreased in both the SD and MDD groups and PPtha-MFG (right) increased in the SD group and decreased in the MDD group.

Conclusion

Through analysis of FC measured by rs-fMRI, the altered FC between PPtha and several brain regions (right and left ITG, right MTG, and right MFG) has been identified in participants with SD and MDD. Different alterations in FC between PPtha and these regions were identified for patients with depression. These findings might provide insights into the potential pathophysiological mechanisms of SD and MDD, especially related to PPtha and the lateral habenula.

Peer Review reports

Introduction

Major depressive disorder (MDD) is characterized by marked and persistent depression and is chronic, recurrent, and difficult to cure. Additionally, MDD is associated with a high suicide rate. These characteristics produce a severe adverse effect on the quality of life of patients with MDD [1,2,3,4,5].

Subclinical depression (SD) is regarded as an early stage indicator or precursor of MDD [6,7,8] because the depressive symptoms of patients with SD are not severe or persistent enough to be diagnosed as MDD [8,9,10]. SD is reported to affect 32% of Chinese college students [8, 11] and 23–39% of European college students [8, 12]. The public’s understanding of the mechanisms of MDD and SD is inadequate, and the diagnosis of MDD and SD still relies on physician interviews and scale tests, with a notable lack of accurate biomarkers.

Resting-state functional magnetic resonance imaging (rs-fMRI) examines the temporal correlation of spontaneous fluctuations in blood oxygen level dependent (BOLD) signals across brain regions [6, 13]. Functional connectivity (FC) measures the correlations between brain regions by calculating the Pearson correlation coefficient. By analyzing the FC between a specific brain region and other brain regions, abnormal connections can be identified. Exploring the FC changes between MDD and SD patients is of great significance [6, 13, 14]. Meta-analysis of rs-fMRI in MDD shows altered functional activity in a wide range of brain regions in the frontal-parietal network, frontal-limbic network, and default mode network [15, 16]. FC analyses of rs-fMRI data have characterized the altered brain network architecture based on intrinsic neural activity in the absence of any specified cognitive or affective load in MDD [17, 18].

Studies have shown that the lateral habenula (LHb) has an important relationship with the pathogenesis of depression [7, 19]. The lateral habenula (LHb) is a small epithelial structure located between the medial thalamus and the third ventricle [20,21,22], connecting the limbic system of the forebrain and the monoamine nucleus of the midbrain. According to the Montreal Neurological Institute (MNI) coordinates, the LHb is in the posterior parietal thalamus (PPtha) [7]. It is the anti-reward center of the brain and is considered to mediate most of the people's negative emotions, such as fear, tension, and anxiety. In depression, the LHb is involved in negative motivational value and decision-making [23,24,25,26,27,28]. A special discharge mode called cluster discharge of the LHb is a sufficient condition for the occurrence of depression [19, 29,30,31].

The LHb has become an anatomical target for testing antidepressants such as ketamine in animals and humans [19, 29, 32,33,34,35]. Concrete evidence demonstrates that bursts within the LHb drive depression in rats [19, 36,37,38]. Deep brain stimulation(DBS) of the LHb has been successfully used to treat patients with refractory MDD [39]. Additionally, clinical studies of DBS targeting the habenula in humans are underway [39,40,41]. Although this line of research is promising, it remains at an early stage and not all observations can be adequately explained by current theories on how LHb dysfunction contributes to depressive states [5]. Moreover, the connections of FC with the LHb are unknown in humans, especially when the FC is calculated from rs-fMRI data in both MDD and SD patients. Despite evidence that developmental trajectories in psychiatric disease are related to habenula volume [42], the FC alterations related to LHb in SD, MDD and health control (HC) groups have not been studied.

In this study, we compare the FC connected with LHb among MDD, SD and HC groups, and aim to clarify the FC alterations in HC, SD, and MDD. This study includes the collection of rs-fMRI data of three groups, extraction of the PPtha as the region of interest (ROI), calculation of the FC between PPtha and all other regions in the whole brain, and statistical analysis of altered FC.

Materials and methods

Participants

All participants in this study were college students from Guangzhou Medical University. A total of 40 MDD (11 males, 29 females), 34 SD (11 males, 23 females) and 40 HC (21 males, 19 females) subjects were screened from about 1000 college students over the period 2018–2019.

The MDD subjects were screened using the Hamilton Depression (HAMD) scale [43]. SD and HC subjects were screened using the Beck Depression Inventory II (BDI-II) scale [44]. MDD subjects had a HAMD score of > 17, SD subjects had a BDI-II score > 13 (total score, 14–28), and HC subjects had a BDI-II score of < 5 [44, 45]. The three groups were matched by gender, age, and education. All diagnostic criteria were following the Diagnostic and Statistical Manual of Mental Disorders-IV (DSM-IV) [46, 47]. Moreover, all participants were required to meet the following criteria: right handedness, no visible lesions on MRI, no neurological diseases, and no alcohol or drug dependence.

MRI acquisition

MRI images were obtained using a 3-Tesla MRI scanner (Siemens, Erlangen, Germany). To control head movement and reduce MRI noise, foam pads and earplugs were used. The rs-fMRI images were obtained using an echo planar imaging sequence: echo time (TE) = 21 ms, repetition time (TR) = 2500 ms, flip angle (FA) = 90°, field of view (FOV) = 200 mm × 200 mm, matrix = 64 × 64, voxel size = 3.5 mm × 3.1 mm × 3.1 mm, and 42 slices with no gap. There were 200 time-points (or volumes) of rs-fMRI collected for each participant. The parameters of the T1-weighted images were as follows: TR = 2530 ms, TE = 2.34 ms, FA = 7°, FOV = 256 mm × 224 mm, and a slice thickness of 1.0 mm with no gap.

The participants were asked to stay relaxed and keep their eyes closed without sleeping. None of the participants showed significant structural damage based on routine MRI scans.

Study design and main procedures

As shown in Fig. 1A, this study included four steps: (1) image preprocessing; (2) calculation of FC; (3) statistical analysis; (4) identification of alterations in FC of the LHb. In brief, the rs-fMRI image preprocessing was conducted and then the FC was calculated. One-way analysis of variance (ANOVA) with the post hoc t-test [48,49,50] was carried out during the statistical analysis. Then alterations in FC of the LHb were identified.

Fig. 1
figure 1

The study design and main procedures. A The overview of the study procedures; B Rs-fMRI image preprocessing; C The ROI extraction and FC calculation between PPtha and other brain regions; D The statistical analysis procedure of FC. (ROI, region of interest; FC, functional connectivity; ANOVA, analysis of variance; rs-fMRI, resting-state functional magnetic resonance imaging; T1-WI, T1-weighted imaging; PPtha, posterior parietal thalamus; R, right; L, left; MDD, major depressive disorder; SD, subclinical depression; HC, healthy control; vs., versus)

Rs-fMRI image preprocessing

The rs-fMRI data were preprocessed using the Data Processing Assistant for Resting-State fMRI (DPARSF, http://rfmri.org/DPARSF) [51] based on Statistical Parametric Mapping (SPM12, http://www.fil.ion.ucl.ac.uk/spm). The substeps are given as follows (Fig. 1B):

  1. (1)

    To eliminate the unstable factors of MRI equipment and subjects, the first ten time points of each subject were removed.

  2. (2)

    The slice timing was corrected by using the middle slice as the reference.

  3. (3)

    Realignment was used to correct the head motion resulting from breathing, heartbeats, and uncontrolled factors during the scan. To control the head motion further, scrubbing was conducted for all rs-fMRI in image preprocessing. The framewise displacement (FD) type is “FD (Power)” [52], the FD threshold is 0.2 mm and each bad time point (FD > 0.2 mm) is removed as a separate regression value [52,53,54]. Finally, the patients with head motions > 2 mm and rotations > 2° were removed.

  4. (4)

    The T1 image was registered to the MNI space using the Diffeomorphic Anatomical Registration Through Exponentiated Lie Algebra (DARTEL) segmentation algorithm [55].

  5. (5)

    To control the head motion as little as possible, the nuisance covariates (i.e., Friston 24-parameter model, cerebrospinal fluid signal, and white matter signal) were regressed out.

  6. (6)

    The normalized T1 image parameters were used to match the functional image to the MNI space. The function images were resampled with a pixel size of 3 mm × 3 mm × 3 mm for normalization.

  7. (7)

    The images were filtered with the band-pass filter of 0.01–0.08 Hz.

  8. (8)

    Finally, a Gaussian kernel with a full-width at half maximum (FWHM) of 6 mm × 6 mm × 6 mm was used to smoothen the images.

Calculation of FC

A seed-to-whole-brain method was used to calculate the FC (Fig. 1C). Specifically, the Brainnetome Atlas [56] with 246 subregions was used to locate the ROI. The left PPtha (239th subregion) and right PPtha (240th subregion) were combined to get the whole PPtha mask. The PPtha mask was selected as the seed and the FC between PPtha and other voxels in the whole brain was calculated.

The seed-based method in the FC module of the DPABI software ([57], http://www.rfmri.org/dpabi) was used to calculate the FC. First, the average time series of PPtha were calculated. Second, the time series of each voxel in the whole brain except the PPtha were extracted. The Pearson correlation coefficient between the PPtha and each voxel was obtained. The r values were converted to z-scores by DPABI to get the z-maps. Finally, the Anatomical Automatic Labeling (AAL) atlas was further utilized to identify where the clusters with significant differences were located.

Statistical analysis

Demographic data from the MDD, SD, and HC groups were compared using the Statistical Package for the Social Sciences software, version 24.0 (SPSS, https://www.ibm.com/products/spss-statistics). One-way ANOVA analysis was performed to assess differences in age and educational level, and a chi-squared test was performed to assess gender differences.

The presence of significant age differences, root mean squared (RMS) framewise displacement (FD) (one-way ANOVA tests, p < 0.05), gender (chi-squared test, p < 0.05), and educational level (one-way ANOVA tests, p < 0.05) among the MDD, SD, and HC groups was tested. The difference in BDI-II score between the SD and HC groups was tested using a two-sample t-test; p < 0.001 indicates a significant statistical difference.

Figure 1D illustrates the statistical analysis procedure. The statistical analysis was conducted by DPABI software. The mean framewise displacement (i.e., Mean FD_Jenkinson) was taken as the covariate to control the impact of unnecessary head motion in the statistical analysis. One-way ANOVA (p < 0.05, Bonferroni corrected) in the statistical analysis module of DPABI software was performed to compare variables among the three groups (HC, SD, and MDD). A mask was built according to the results of ANOVA. Based on the mask, the inter-group differences were obtained by using the post hoc t-test. Two-tailed Gaussian random field (GRF) correction [58] (voxel threshold of p < 0.005 and cluster threshold of p < 0.05, and voxel threshold of p < 0.001 and cluster threshold of p < 0.05) was performed during the two-sample t-test (the post hoc test). Finally, the brain regions with significant FC connected with PPtha were identified.

Results

Demographic and clinical characteristics

A total of 40 MDD subjects, 34 SD subjects, and 40 HC subjects participated in the rs-fMRI process. Using the head motion standard of 2 mm translation and 2° rotation, 2 MDD subjects, 8 SD subjects, and 7 HC subjects were excluded. Finally, the data from 38 MDD subjects, 26 SD subjects, and 33 HC subjects were included in the statistical analysis. Demographic and clinical data are presented in Table 1. No significant differences in terms of age, gender, and educational level (p > 0.05) were found among the MDD, SD, and HC groups. All participants with MDD and SD showed higher symptom scores than those in the HC group (p < 0.001).

Table 1 Demographic and clinical characteristic of three groups of participants

The RMS FD for each participant were calculated and compared across groups using one-way ANOVA tests. The results presented in Table 1, do not show any significant differences in RMS FD among the three groups.

ANOVA results for MDD, SD, and HC groups

Table 2 and Fig. 2 show the one-way ANOVA analysis of the FC among the three groups. Significant group differences in FC were observed in Temporal_Inf_L (AAL), Superior Temporal Gyrus, Temporal_Pole_Mid_R (AAL), Middle Frontal Gyrus, Frontal_Sup_R (AAL), Medial Frontal Gyrus, Cerebellum Posterior Lobe, Cerebellum Anterior Lobe, Middle Temporal Gyrus, Precentral Gyrus, Supp_Motor_Area_R (AAL), Precuneus_R (AAL), Postcentral Gyrus, and Parietal_Sup_L (AAL).

Table 2 Brain regions with significantly different FC connected with the posterior parietal thalamus among the SD, MDD patients, and healthy controls
Fig. 2
figure 2

Brain regions with significant differences in functional connectivity connected with the PPtha among three groups (p < 0.05, z > 3.1) after ANOVA analysis

Altered FC connected with PPtha for MDD versus SD

The significantly different FC between the PPtha and other brain regions in the MDD and SD groups is shown in Table 3 and Fig. 3(A). The MDD group presents the significantly increased FC between the PPtha and left inferior temporal gyrus (ITG) (Temporal_Inf_L [AAL]) in comparison with SD. Moreover, the FC between the PPtha and right superior temporal gyrus (STG), middle temporal gyrus (MTG), middle frontal gyrus (MFG), dorsolateral superior frontal gyrus (Frontal_Sup_R [AAL]), and medial frontal gyrus is decreased in the MDD group.

Table 3 Brain regions with significantly different FC connected with the posterior parietal thalamus
Fig. 3
figure 3

Brain regions with significant differences in functional connectivity connected with the posterior parietal thalamus among three groups (cluster p < 0.05, z > 2.8). Color bars represent the t-value from group analysis. A MDD vs. SD. B MDD vs. HC. C SD vs. HC

Altered FC connected with PPtha for MDD versus HC

In the MDD group, the FC between the PPtha and right cerebellum posterior lobe and cerebellum anterior lobe is higher than in the HC group. However, in comparison with HC, the FC between the PPtha and right MTG, precentral gyrus, Supp_Motor_Area_R [AAL], and precuneus is lower in the MDD group. Detailed results are shown in Table 3 and Fig. 3B.

Altered FC connected with PPtha for SD versus HC

The FC between the PPtha and right ITG (Temporal_Inf_R [AAL]), MFG, and postcentral gyrus is higher in the SD group than that in the HC group. Conversely, the FC between the PPtha and left superior parietal gyrus (SPG, Parietal_Sup_L [AAL]) is lower, as shown in Table 3 and Fig. 3C.

Important FC alteration for HC, SD and MDD

By comparing the FC alterations in MDD vs. SD, MDD vs. HC and SD vs. HC, the potential tendency for depression have been identified. Specifically, increased FC was found in ITG in the comparison of MDD vs. SD (left side) and SD vs. HC (right side). Conversely, decreased FC of the MDD group was observed between the PPtha and right MTG in the comparison of MDD vs. SD and MDD vs. HC. Additionally, the FC between the PPtha and right MFG was higher in the SD group than in both the MDD and HC groups. The locations of the four identified regions, (left ITG, right ITG, right MFG, and right MTG) are shown in Fig. 4. The FC value for MFG.R, MTG.R, ITG.L and ITG.R for each group is shown in Fig. 5.

Fig. 4
figure 4

Four main brain regions characterizing the alterations of functional connectivity connected with the posterior parietal thalamus in HC, SD, and MDD groups (ITG.L, the left inferior temporal gyrus; ITG.R, the right inferior temporal gyrus; MTG.R, the middle frontal gyrus; MFG.R, the right middle temporal gyrus)

Fig. 5
figure 5

The FC value for the significant ROIs for each group (ITG.L, the left inferior temporal gyrus; ITG.R, the right inferior temporal gyrus; MTG.R, the middle frontal gyrus; MFG.R, the right middle temporal gyrus)

It should be noted that GRF correct method was used with two criteria: (1) voxel p < 0.005 and cluster p < 0.05; (2) voxel p < 0.001 and cluster p < 0.05. As shown in Table 3, the brain regions with significant differences were the same in the two cases. The voxel size of clusters with significant differences was small (e.g., 1 voxel in some cases) while voxel p < 0.001 and cluster p < 0.05, potentially suggesting that the correct standard was too strict.

Discussion

To the best of our knowledge, this is the first study comparing the FC of PPtha among patients with MDD, SD, and HC. Altered FC between PPtha and several regions was identified in all three pairwise group comparisons. The regions included the ITG, MTG, and MFG. The FC alteration tendency from HC to SD and MDD was identified by comparing the three groups. These findings provide a new perspective and thinking for exploring the pathophysiological mechanism of depression. The significance of these findings is discussed in the following subsections.

Importance of LHb and PPtha in depression

SD has been considered the precursor to MDD. The cluster discharge of the LHb in rats is sufficient for the occurrence of depression [19, 29, 31]. Undoubtedly, the FC between the LHb and other related regions is affected. The LHb is an important subregion of the PPtha. In this study, PPtha was chosen as the ROI for FC analysis to research which region the LHb influenced. We have used two steps to ensure the habenula is located within the PPtha ROI for all participants. First, according to the reported locations of the habenula and PPtha in MNI coordinates, we can know that the habenula is located in the PPtha region. In our previous study [7], the habenula was located in the PPtha region in the same way as used in the current study, and abnormal brain connections with PPtha had been identified. Second, after determining the PPtha ROI, one experienced neurologist confirmed that the habenula was located within the PPtha ROI for each participant via visual inspection of T1-W images and prior anatomical knowledge.

The PPtha in the Brainnetome Atlas was used as the ROI, instead of LHb, because the precise location of the LHb is difficult to determine. First, the small volume of the LHb was a limitation in some studies. The habenula volume in vivo is estimated to be closer to 18.5 mm3 per hemisphere based on structural MRI and postmortem measurements [59, 60]. This may be an underestimation due to the difficulty of resolving the lateral and anterior boundaries of the habenula. The LHb volume is smaller than 18.5 mm3 per hemisphere, which can be smaller than the voxel size of standard fMRI [20, 59, 61]. Second, 3 T fMRI with whole-brain coverage and appropriate repetition times could only be reliably achieved with a resolution ≥ 27mm3. In this case, the conventional fMRI voxel size is larger than the LHb. Finally, the smoothing kernels with 5–12 mm FWHM are larger than the LHb. The habenula signal is likely to be contaminated by adjacent structures, such as the medial dorsal thalamus or the epithalamic paraventricular nucleus [20].

Increased FC between PPtha and ITG in both MDD and SD

All three pairwise group comparisons showed abnormal FC between the PPtha and the temporal lobe. Specifically, the FC between the PPtha and ITG in the MDD group was higher than that in the SD group. FC in the SD group was higher than that in the HC group. It is speculated that the FC between the PPtha and ITG is gradually enhanced during the progression from HC to MDD.

Decreased FC between PPtha and MTG in MDD

The FC between the PPtha and right MTG was lower in the MDD group than that in the SD and HC groups. It is speculated that the FC between the PPtha and right MTG gradually decreases as depression progresses. One study reported a decreased connectivity between the prefrontal thalamus and the MTG in MDD patients [6], which may provide a basis for estimating the development stage of depression. Comparing the SD and HC groups, no significant differences in FC were found between the PPtha and MTG. This finding might suggest that the PPtha-MTG connection abnormality is common and consistent in both the SD and MDD groups.

The ITG and the MTG are crucial regions in the temporal cortex. Some research has reported that thalamo-temporal connectivity changes significantly in MDD patients [62]. Moreover, higher thalamo-temporal connectivity is associated with more severe depressive symptoms, suggesting an association with the core psychopathology of MDD [62, 63]. Numerous PET and fMRI studies have shown that some lateral parts of the temporal cortex play significant roles in aspects of affective processing and social cognition [64]. Other studies have considered that the overall function of the anterior temporal lobe may be the semantic processing, which is personally, socially, and emotionally relevant [65]. Additionally, modest evidence suggests the presence of small clusters of decreased connectivity between the frontal cortex and medial thalamus [62]. Our findings are consistent with those of other neuroimaging studies, and possibly reflect the emotional and somatosensory dysregulation in SD and MDD patients.

Abnormal FC between PPtha and MFG in SD

In the MDD and SD groups, and the SD and HC groups, the FC between the PPtha and frontal lobes showed an abnormal value. The FC between the PPtha and right MFG in the MDD group was lower than that in the SD group. Moreover, the FC in the SD group was higher than that in the HC group. The FC between the hippocampus and MFG was impaired in MDD patients [66, 67]. Moreover, altered FC between the hippocampus and MFG has been reported in older adults with SD [68]. The left MFG indicates increased spontaneous neural activity in SD patients [8]. MFG is involved in regulating the intensity of response to emotional stimuli, and changes in MFG functional activities may lead to inappropriate responses to emotional events [69]. The depressive symptoms occur as a result.

Notably, the MFG is a crucial part of the dorsolateral prefrontal cortex (DLPFC), which has been tightly linked to depression [8]. In depression, DLPFC is involved in mood regulation. The DLPFC activity is inhibited at rest but increases when symptoms are relieved [70, 71]. Most importantly, regions of the DLPFC have been the targets for repetitive transcranial magnetic stimulation (rTMS) in depression treatment [7, 72, 73].

The mechanism underlying the higher PPtha-MFG connection in the SD group compared with both the MDD and HC groups is unknown. Here, we speculate that the connection between the PPtha and MFG may be severely changed in MDD patients. Another possibility is that the MFG or the PPtha may have been impaired in the post-depression stage.

Limitations and future works

There are several limitations in this study. First, although this study included three groups of rs-fMRI data, the sample size is largely small. This factor limited the statistical significance. Second, we didn’t extract the accurate LHb as the ROI but used the PPtha region from the Brainnetome Atlas. This change may affect our findings. This study investigated the FC between the PPtha and ITG, MTG, and MFG for MDD and SD. The validity of our findings requires verification through further studies.

In future work, the sample size will be expanded. Subjects from the WU-Minn Human Connectome Project (HCP) Consortium’s 500 Subjects Release can be included [60, 71]. Furthermore, 246 × 246 matrices of functional brain networks can be constructed using the human Brainnetome Atlas [56]. The information of each node and its connections can be further analyzed from the viewpoint of brain networks. Finally, after calculating the FC connections among different brain regions, machine learning models could be constructed to classify the MDD and SD participants [74].

Conclusion

Through analyzing the FC measured by using rs-fMRI data, this study explored the altered FC between PPtha and other brain regions in SD and MDD subjects. The FC between the PPtha and ITG, right MTG, and right MFG showed abnormality in both MDD and SD subjects. Different alterations for FC between PPtha and these brain regions were present for HC, SD, and MDD: increased PPtha-ITG (left and right) FC; decreased PPtha-MTG (right) FC; initially increased and subsequently decreased PPtha-MFG (right). These findings might provide insight into potential pathophysiologic mechanisms of SD and MDD, especially related to the lateral habenula.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

Abbreviations

AAL:

Anatomical automatic labeling

ANOVA:

Analysis of variance

BDI-II:

Beck Depression Inventory-II

BOLD:

Blood oxygen level dependent

DARTEL:

Diffeomorphic Anatomical Registration Through Exponentiated Lie Algebra

DBS:

Deep brain stimulation

DLPFC:

Dorsolateral prefrontal cortex

DPARSF:

Data processing assistant for resting-state fMRI

DSM-IV:

Diagnostic and Statistical Manual of Mental Disorders-IV

EPI:

Echo planar imaging

FA:

Flip angle

FC:

Functional connectivity

FOV:

Field of view

FWHM:

Full-width at half maximum

GRF:

Gaussian random field

HAMD:

Hamilton Depression

HC:

Healthy control

ITG:

Inferior temporal gyrus

LHb:

Lateral habenula

MDD:

Major depressive disorder

MFG:

Middle frontal gyrus

MNI:

Montreal neurological institute

MTG:

Middle temporal gyrus

PPtha:

Posterior parietal thalamus

ROI:

Region of interest

rs-fMRI:

Resting-state functional magnetic resonance imaging

rTMS:

Repetitive transcranial magnetic stimulation

SD:

Subclinical depression

SPM:

Statistical Parametric Mapping

SPSS:

Statistical package for the social sciences

TE:

Echo time

TR:

Repetition time

References

  1. Fava M, Kendler KS. Major depressive disorder. Neuron. 2000;28(2):335–41. https://doi.org/10.1016/S0896-6273(00)00112-4.

    Article  CAS  PubMed  Google Scholar 

  2. Tsuno N, Besset A, Ritchie K. Sleep and depression. J Clin Psychiatry. 2005;66(10):1254–69. https://doi.org/10.4088/JCP.v66n1008.

    Article  PubMed  Google Scholar 

  3. American Psychiatric Association, AP. Diagnostic and statistical manual of mental disorders (DSM-5®). Washington: American Psychiatric Pub; 2013.

    Book  Google Scholar 

  4. Li M, Metzger CD, Li W, Safron A, van Tol MJ, Lord A, et al. Dissociation of glutamate and cortical thickness is restricted to regions subserving trait but not state markers in major depressive disorder. J Affect Disord. 2014;169:91–100. https://doi.org/10.1016/j.jad.2014.08.001.

    Article  CAS  PubMed  Google Scholar 

  5. Liu CH, Zhang GZ, Li B, Li M, Woelfer M, Walter M, Wang LH. Role of inflammation in depression relapse. J Neuroinflamm. 2019;16:90. https://doi.org/10.1186/s12974-019-1475-7.

    Article  Google Scholar 

  6. Kong QM, Qiao H, Liu CZ, Zhang P, Li K, et al. Aberrant intrinsic functional connectivity in thalamo-cortical networks in major depressive disorder. CNS Neurosci Ther. 2018;24:1063–72. https://doi.org/10.1111/cns.12831.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Zhu YK, Qi SL, Zhang B, He DN, Teng YY, Hu JN, et al. Connectome-based biomarkers predict sub-clinical depression and identify abnormal brain connections with the lateral habenula and thalamus. Front Psychiatry. 2019;10:371. https://doi.org/10.3389/fpsyt.2019.00371.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Zhang B, Qi SL, Liu S, Liu XY, Wei XH, Ming D. Altered spontaneous neural activity in the precuneus, middle and superior frontal gyri, and hippocampus in college students with subclinical depression. BMC Psychiatry. 2021;21(1):1–10. https://doi.org/10.1186/s12888-021-03292-1.

    Article  Google Scholar 

  9. Fergusson DM, Horwood LJ, Ridder EM, Beautrais AL. Subthreshold depression in adolescence and mental health outcomes in adulthood. Arch Gen Psychiatry. 2005;62(1):66–72. https://doi.org/10.1001/archpsyc.62.1.66.

    Article  PubMed  Google Scholar 

  10. Segal DL. Diagnostic and statistical manual of mental disorders (DSM-IV-TR). Corsini Encyclopedia Psychol. 2010;1–3. https://doi.org/10.1002/9780470479216.corpsy0271.

  11. Jiang L, Wang Y, Zhang Y, Li R, Wu H, Li C, Wu YL, Qian T. The reliability and validity of the center for epidemiologic studies depression scale (CES-D) for Chinese university students. Front Psychiatry. 2019;10:315. https://doi.org/10.3389/fpsyt.2019.00315.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Mikolajczyk RT, Maxwell AE, El Ansari W, Naydenova V, Stock C, Ilieva S, et al. Prevalence of depressive symptoms in university students from Germany, Denmark, Poland and Bulgaria. Soc Psychiatry Psychiatr Epidemiol. 2008;43(2):105–12. https://doi.org/10.1007/s00127-007-0282-0.

    Article  PubMed  Google Scholar 

  13. Sundermann B, Beverborg M, Pfleiderer B. Meta-analysis of resting-state fMRI in depression: generating spatial hypotheses for potential clinical applications. Peerj. 2017;2. https://doi.org/10.7287/peerj.preprints.412v1.

  14. Woodward ND, Heckers S. Mapping thalamocortical functional connectivity in chronic and early stages of psychotic disorders. Biol Psychiatry. 2016;79:1016–25. https://doi.org/10.1016/j.biopsych.2015.06.026.

    Article  PubMed  Google Scholar 

  15. Graham J, Salimi-Khorshidi G, Hagan C, Walsh N, Goodyer I, Lennox B. Meta-analytic evidence for neuroimaging models of depression: state or trait? J Affective Dis. 2013;151(2):423–31. https://doi.org/10.1016/j.jad.2013.07.002.

    Article  Google Scholar 

  16. Kaiser RH, Andrews-Hanna JR, Wager TD, Pizzagalli DA. Large-scale network dysfunction in major depressive disordera meta-analysis of resting-state functional connectivity. JAMA Psychiat. 2015;72(6):1–10. https://doi.org/10.1001/jamapsychiatry.2015.0071.

    Article  Google Scholar 

  17. Raichle ME. The restless brain: how intrinsic activity organizes brain function. Philos Trans R Soc Lond B. 2015;370:20140172. https://doi.org/10.1098/rstb.2014.0172.

    Article  Google Scholar 

  18. Drysdale AT, Grosenick L, Downar J, Dunlop K, Mansouri F, Meng Y, et al. Resting-state connectivity biomarkers define neurophysiological subtypes of depression. Nat Med. 2017;23:28–38. https://doi.org/10.1038/nm.4246.

    Article  CAS  PubMed  Google Scholar 

  19. Yang Y, Cui Y, Sang K, Dong Y, Ni Z, Ma S, et al. Ketamine blocks bursting in the lateral habenula to rapidly relieve depression. Nature. 2018;554(7692):317–22. https://doi.org/10.1038/nature25509.

    Article  CAS  PubMed  Google Scholar 

  20. Lawson RP, Drevets WC, Roiser JP. Defining the habenula in human neuroimaging studies. Neuroimage. 2013;64:722–7. https://doi.org/10.1016/j.neuroimage.2012.08.076.

    Article  PubMed  Google Scholar 

  21. Boulos LJ, Darcq E, Kieffer BL. Translating the habenula-from rodents to humans. Biol Psychiatry. 2016;81(4):296–305. https://doi.org/10.1016/j.biopsych.2016.06.003.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Pauli W, Nili A, Tyszka J. A high-resolution probabilistic in vivo atlas of human subcortical brain nuclei. Scientific Data. 2018;5:180063. https://doi.org/10.1038/sdata.2018.63.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Shumake J, Edwards E, Gonzalez-Lima F. Opposite metabolic changes in the habenula and ventral tegmental area of a genetic model of helpless behavior. Brain Res. 2003;963(1):274–81. https://doi.org/10.1016/S0006-8993(02)04048-9.

    Article  CAS  PubMed  Google Scholar 

  24. Matsumoto M, Hikosaka O. Lateral habenula as a source of negative reward signals in dopamine neurons. Nature. 2007;447:1111–5. https://doi.org/10.1038/nature05860.

    Article  CAS  PubMed  Google Scholar 

  25. Hikosaka O. The habenula: from stress evasion to value-based decision-making. Nat Rev Neurosci. 2010;11(7):503–13. https://doi.org/10.1038/nrn2866.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Shumake J, Gonzalez-Lima F. Functional opposition between habenula metabolism and the brain reward system. Front Hum Neurosci. 2013;7:662. https://doi.org/10.3389/fnhum.2013.00662.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Lawson RP, Seymour B, Loh E, Lutti A, Dolan RJ, Dayan P, Weiskopf N, Roiser JP. The habenula encodes negative motivational value associated with primary punishment in humans. Proc Natl Acad Sci U S A. 2014;111(32):11858–63. https://doi.org/10.1073/pnas.1323586111.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Stopper CM, Floresco SB. What’s better for me? Fundamental role for lateral habenula in promoting subjective decision biases. Nat Neurosci. 2014;17(1):33–5. https://doi.org/10.1038/nn.3587.

    Article  CAS  PubMed  Google Scholar 

  29. Li K, Zhou T, Liao L, Yang ZF, Wong C, Henn F, Malinow R, Hu HL. βCaMKII in lateral habenula mediates core symptoms of depression. Science. 2013;341(6149):1016–20. https://doi.org/10.1126/science.1240729.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Lawson RP, Nord CL, Seymour B, Thomas DL, Dayan P, Pilling S, Roiser JP. Disrupted habenula function in major depression. Mol Psychiatry. 2017;22(2):202–8. https://doi.org/10.1038/mp.2016.81.

    Article  CAS  PubMed  Google Scholar 

  31. Cui YH, Yang Y, Ni ZY, Dong YY, Cai GH, Alexandre Foncelle A, Ma SS, Sang KN, Tang SY, Li YZ, Shen Y, Berry H, Wu XH, Hu HL. Astroglial Kir41 in the lateral habenula drives neuronal bursts in depression. Nature. 2018;554:323–7. https://doi.org/10.1038/nature25752.

    Article  CAS  PubMed  Google Scholar 

  32. Shabel SJ, Wang C, Monk B, Aronson S, Malinow R. Stress transforms lateral habenula reward responses into punishment signals. Proc Natl Acad Sci USA. 2019;116:12488–93. https://doi.org/10.1073/pnas.1903334116.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Klein ME, Chandra J, Sheriff S, Malinow R. Opioid system is necessary but not sufficient for antidepressive actions of ketamine in rodents. Proc Natl Acad Sci USA. 2020;117:2656–62. https://doi.org/10.1073/pnas.1916570117.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Hu H, Cui Y, Yang Y. Circuits and functions of the lateral habenula in health and in disease. Nat Rev Neurosci. 2020;21(5):277–95. https://doi.org/10.1038/s41583-020-0292-4.

    Article  CAS  PubMed  Google Scholar 

  35. Lipina TV, Khrapova MV, Serykh A, Dubrovina NI, Petrova ES, Mikhnevich N, Starostina MV, Amstyslavskaja TG. The increased density of the habenular neurons, high impulsivity, aggression and resistant fear memory in Disc1-Q31 L genetic mouse model of depression. Behav Brain Res. 2020;392:112693. https://doi.org/10.1016/j.bbr.2020.112693.

    Article  CAS  PubMed  Google Scholar 

  36. Caldecott-Hazard S, Mazziotta J, Phelps M. Cerebral correlates of depressed behavior in rats, visualized using 14C–2-deoxyglucose autoradiography. J Neurosci. 1988;8(6):1951–61. https://doi.org/10.1523/JNEUROSCI.08-06-01951.1988.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Li B, Piriz J, Mirrione M, Chung C, Proulx CD, Schulz D, et al. Synaptic potentiation onto habenula neurons in the learned helplessness model of depression. Nature. 2011;470:535–9. https://doi.org/10.1038/nature09742.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Stamatakis AM, Stuber GD. Activation of lateral habenula inputs to the ventral midbrain promotes ehavioral avoidance. Nat Neurosci. 2012;15:1105–7. https://doi.org/10.1038/nn.3145.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Sartorius A, Kiening KL, Kirsch P, von Gall CC, Haberkorn U, Unterberg AW, Henn FA, Meyer-Lindenberg A. Remission of major depression under deep brain stimulation of the lateral habenula in a therapy-refractory patient. Biol Psychiatry. 2010;67:e9–11. https://doi.org/10.1016/j.biopsych.2009.08.027.

    Article  PubMed  Google Scholar 

  40. Dandekar MP, Fenoy AJ, Carvalho AF, Soares JC, Quevedo J. Deep brain stimulation for treatment-resistant depression: An integrative review of preclinical and clinical findings and translational implications. Mol Psychiatry. 2018;23:1094–112. https://doi.org/10.1038/mp.2018.2.

    Article  CAS  PubMed  Google Scholar 

  41. Nuno-Perez A, Tchenio A, Mameli M, Lecca S. Lateral Habenula Gone Awry in Depression: Bridging Cellular Adaptations with Therapeutics. Front Neurosci. 2018;12:485. https://doi.org/10.3389/fnins.2018.00485.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Germann J, Gouveia FV, Martinez RCR, et al. Fully Automated habenula segmentation provides robust and reliable volume estimation across large magnetic resonance imaging datasets, suggesting intriguing developmental trajectories in psychiatric disease. Biol Psychiatry Cogn Neurosci Neuroimaging. 2020;5(9):923–9. https://doi.org/10.1016/j.bpsc.2020.01.004.

    Article  PubMed  Google Scholar 

  43. Sharp R. The Hamilton Rating Scale for Depression. Occup Med. 2015;65(4):340. https://doi.org/10.1093/occmed/kqv043.

    Article  Google Scholar 

  44. Dozois DJA, Dobson KS, Ahnberg JL. A psychometric evaluation of the Beck Depression Inventory–II. Psychol Assess. 1998;10(2):83. https://doi.org/10.1037/1040-3590.10.2.83.

    Article  Google Scholar 

  45. Titov N, Dear BF, Mcmillan D, et al. Psychometric comparison of the PHQ-9 and BDI-II for measuring response during treatment of depression. Cogn Behav Ther. 2011;40(2):126–36. https://doi.org/10.1080/16506073.2010.550059.

    Article  PubMed  Google Scholar 

  46. American Psychiatric Association. Diagnostic and statistical manual of mental disorders (DSM-IV-TR). Washington: American Psychiatric Pub; 2000. p. 339–50.

    Google Scholar 

  47. Volkmar FR. DSM-IV. 2013. https://doi.org/10.1007/978-1-4419-1698-3_1444

  48. Chen T, Xu M, Tu J, et al. Relationship between Omnibus and Post-Hoc Tests: An Investigation of Performance of the F Test in ANOVA. Shanghai Arch Psychiatry. 2018;30(1):60–4. https://doi.org/10.11919/j.issn.1002-0829.218014.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Gurvich V, Naumova M. Logical contradictions in the One-way ANOVA and Tukey-Kramer multiple comparisons tests with more than two groups of observations. arXiv preprint arXiv. 2021;2104.07552. https://doi.org/10.3390/sym13081387.

  50. Yap YT, Gouwanda D, Gopalai AA, Zheng CY. The effect of asymmetrical gait induced by unilateral knee brace on the knee flexor and extensor muscles. Med Biol Eng Compu. 2021;59(3):711–20. https://doi.org/10.1007/s11517-021-02337-7.

    Article  Google Scholar 

  51. Yan CG, Zang YF. DPARSF: a MATLAB toolbox for“pipeline”data analysis of resting-state fMRI. Front Syst Neurosci. 2010;4:13. https://doi.org/10.3389/fnsys.2010.00013.

    Article  Google Scholar 

  52. Power JD, Schlaggar BL, Petersen SE. Recent progress and outstanding issues in motion correction in resting-state fMRI. Neuroimage. 2015;105:536–51. https://doi.org/10.1016/j.neuroimage.2014.10.044.

    Article  PubMed  Google Scholar 

  53. Yan CG, Cheung B, Kelly C, Colcombe S, Craddock RC, Di M, et al. A comprehensive assessment of regional variation in the impact of head micromovements on functional connectomics. Neuroimage. 2013;76:183–201. https://doi.org/10.1016/j.neuroimage.2013.03.004.

    Article  PubMed  Google Scholar 

  54. Satterthwaite TD, Elliott MA, Gerraty RT, et al. An improved framework for confound regression and filtering for control of motion artifact in the preprocessing of resting-state functional connectivity data. Neuroimage. 2013;64:240–56. https://doi.org/10.1016/j.neuroimage.2012.08.052.

    Article  PubMed  Google Scholar 

  55. Bergouignan L, Chupin M, Czechowska Y, Kinkingnehun S, Lemogne C, Bastard GL, et al. Can voxel based morphometry, manual segmentation and automated segmentation equally detect hippocampal volume differences in acute depression? Neuroimage. 2009;45(1):29–37. https://doi.org/10.1016/j.neuroimage.2008.11.006.

    Article  PubMed  Google Scholar 

  56. Fan L, Li H, Zhuo J, Zhang Y, Wang J, Chen L, et al. The human brainnetome atlas: a new brain atlas based on connectional architecture. Cereb Cortex. 2016;26:3508. https://doi.org/10.1093/cercor/bhw157.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Yan CG, Wang XD, Zuo XN, Zang YF. DPABI: Data Processing & Analysis for (Resting-State) Brain Imaging. Neuroinformatics. 2016;14:339–51. https://doi.org/10.1007/s12021-016-9299-4.

    Article  PubMed  Google Scholar 

  58. Chen X, Lu B, Yan CG. Reproducibility of r-fmri metrics on the impact of different strategies for multiple comparison correction and sample sizes. Hum Brain Mapp. 2018;39(1):300–18. https://doi.org/10.1002/hbm.23843.

    Article  PubMed  Google Scholar 

  59. Ranft K, Dobrowolny H, Krell D, Bielau H, Bogerts B, Bernstein HG. Evidence for structural abnormalities of the human habenular complex in affective disorders but not in schizophrenia. Psychol Med. 2010;40(04):557–67. https://doi.org/10.1017/S0033291709990821.

    Article  CAS  PubMed  Google Scholar 

  60. Ely BA, Xu J, Goodman WK, Lapidus KA, Gabbay V, Stern ER. Resting-state functional connectivity of the human habenula in healthy individuals: associations with subclinical depression. Hum Brain Mapp. 2016;37(7):2369–84. https://doi.org/10.1002/hbm.23179.

    Article  PubMed  PubMed Central  Google Scholar 

  61. Savitz JB, Nugent AC, Bogers W, Roiser JP, Bain EE, Neumeister A, et al. Habenula volume in bipolar disorder and major depressive disorder: a high-resolution magnetic resonance imaging study. Biol Psychiatry. 2011;69(4):336–43. https://doi.org/10.1016/j.biopsych.2010.09.027.

    Article  PubMed  Google Scholar 

  62. Brown EC, Clark DL, Hassel S, et al. Thalamocortical connectivity in major depressive disorder. J Affect Disord. 2017;217:125–31. https://doi.org/10.1016/j.jad.2017.04.004.

    Article  PubMed  Google Scholar 

  63. Jin CY, Qi SL, Teng YY, Li C, Yao YD, Ruan XH, Wei XH. Altered Degree Centrality of Brain Networks in Parkinson’s Disease With Freezing of Gait: A Resting-State Functional MRI Study. Front Neurol. 2021;12:743135. https://doi.org/10.3389/fneur.2021.743135.

    Article  PubMed  PubMed Central  Google Scholar 

  64. Blair R, Morris JS, Frith CD, et al. Dissociable neural responses to facial expressions of sadness and anger. Brain. 1999;122(Pt 5):883–93. https://doi.org/10.1093/brain/122.5.883.

    Article  PubMed  Google Scholar 

  65. Wong C, Gallate J. The function of the anterior temporal lobe: A review of the empirical evidence. Brain Res. 2012;1449(none):94–116. https://doi.org/10.1016/j.brainres.2012.02.017.

    Article  CAS  PubMed  Google Scholar 

  66. Ye T, Peng J, Nie B, Gao J, Liu JT, Li Y, et al. Altered functional connectivity of the dorsolateral prefrontal cortex in first-episode patients with major depressive disorder. Eur J Radiol. 2012;81(12):4035–40. https://doi.org/10.1016/j.ejrad.2011.04.058.

    Article  PubMed  Google Scholar 

  67. Cao X, Liu Z, Cheng X, et al. Disrupted resting-state functional connectivity of the hippocampus in medication-nave patients with major depressive disorder. J Affect Disord. 2012;141(2–3):194–203. https://doi.org/10.1016/j.jad.2012.03.002.

    Article  PubMed  Google Scholar 

  68. Zhu XY, Li R, Wang PY, Li J. Aberrant functional connectivity of the hippocampus in older adults with subthreshold depression. Psych J. 2014;3(4):245–53. https://doi.org/10.1002/pchj.60.

    Article  PubMed  Google Scholar 

  69. Niendam TA, Laird AR, Ray KL, Monica Dean Y, Glahn DC, Carter CS. Meta-analytic evidence for a superordinate cognitive control network subserving diverse executive functions. Cogn Affect Behav Neurosci. 2012;12(2):241–68. https://doi.org/10.3758/s13415-011-0083-5.

    Article  PubMed  PubMed Central  Google Scholar 

  70. Koenigs M, Grafmanb J. The functional neuroanatomy of depression: distinct roles for ventromedial and dorsolateral prefrontal cortex. Behav Brain Res. 2009;201:239–43. https://doi.org/10.1016/j.bbr.2009.03.004.

    Article  PubMed  PubMed Central  Google Scholar 

  71. Van Essen DC, Smith SM, Barch DM, Behrens TE, Yacoub E, Ugurbil K, et al. The WU-Minn Human Connectome Project: an overview. Neuroimage. 2013;80:62–79. https://doi.org/10.1016/j.neuroimage.2013.05.041.

    Article  PubMed  Google Scholar 

  72. Downar J, Daskalakis ZJ. New targets for rTMS in depression: a review of convergent evidence. Brain Stimul. 2013;6(3):231–40. https://doi.org/10.1016/j.brs.2012.08.006.

    Article  PubMed  Google Scholar 

  73. Meng H, Wang Y, Huang M, Lin W, Wang S, Zhang B. Chronic deep brain stimulation of the lateral habenula nucleus in a rat model of depression. Brain Res. 2011;1422:32–8. https://doi.org/10.1016/j.brainres.2011.08.041.

    Article  CAS  PubMed  Google Scholar 

  74. Gao S, Calhoun VD, Sui J. Machine learning in major depression: From classification to treatment outcome prediction. CNS Neurosci Ther. 2018;24(11):1037–52. https://doi.org/10.1111/cns.13048.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Not Applicable (NA)

Funding

This work was partly supported by the National Natural Science Foundation of China under Grant (Nos. 82072008, 81871846), the Fundamental Research Funds for the Central Universities (N2119010), and the Natural Science Foundation of Guangdong Province, P.R. China (2021A1515011288), and the Science and Technology Project of Guangzhou, P.R. China (202102010020).

Author information

Authors and Affiliations

Authors

Contributions

LY performed experiments and analyzed the data along with CJ and SQ. SQ, YY, and XW conceived the study, presented the results, and wrote the manuscript along with LY. XR collected and analyzed the data. YT and CL supervised the algorithm development and analyzed the data. All authors read and approved the final manuscript.

Author’s information

Not Applicable (NA).

Corresponding author

Correspondence to Shouliang Qi.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the medical ethics committee of Guangzhou First People's Hospital Affiliated to Guangzhou Medical University. All methods were carried out in accordance with the Declaration of Helsinki and relevant guidelines and regulations. All participants had signed the written informed consent.

Consent for publication

Not applicable (NA).

Competing interests

The authors have declared no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Yang, L., Jin, C., Qi, S. et al. Alterations of functional connectivity of the lateral habenula in subclinical depression and major depressive disorder. BMC Psychiatry 22, 588 (2022). https://doi.org/10.1186/s12888-022-04221-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12888-022-04221-6

Keywords