- 1Sino-Dutch Biomedical and Information Engineering School, Northeastern University, Shenyang, China
- 2Key Laboratory of Medical Image Computing of Northeastern University (Ministry of Education), Shenyang, China
- 3Department of Radiology, Wayne State University, Detroit, United States
- 4Department of Radiology, Guangzhou First People’s Hospital, School of Medicine, South China University of Technology, Guangzhou, China
Subclinical depression (SD) has been considered as the precursor to major depressive disorder. Accurate prediction of SD and identification of its etiological origin are urgent. Bursts within the lateral habenula (LHb) drive depression in rats, but whether dysfunctional LHb is associated with SD in human is unknown. Here we develop connectome-based biomarkers which predict SD and identify dysfunctional brain regions and connections. T1 weighted images and resting-state functional MRI (fMRI) data were collected from 34 subjects with SD and 40 healthy controls (HCs). After the brain is parcellated into 48 brain regions (246 subregions) using the human Brainnetome Atlas, the functional network of each participant is constructed by calculating the correlation coefficient for the time series of fMRI signals of each pair of subregions. Initial candidates of abnormal connections are identified by the two-sample t-test and input into Support Vector Machine models as features. A total of 24 anatomical-region-based models, 231 sliding-window-based models, and 100 random-selection-based models are built. The performance of these models is estimated through leave-one-out cross-validation and evaluated by measures of accuracy, sensitivity, confusion matrix, receiver operating characteristic curve, and the area under the curve (AUC). After confirming the region with the highest accuracy, subregions within the thalamus and connections associated with subregions of LHb are compared. It is found that five prediction models using connections of the thalamus, posterior superior temporal sulcus, cingulate gyrus, superior parietal lobule, and superior frontal gyrus achieve an accuracy >0.9 and an AUC >0.93. Among 90 abnormal connections associated with the thalamus, the subregion of the right posterior parietal thalamus where LHb is located has the most connections (n = 18), the left subregion only has 3 connections. In SD group, 10 subregions in the thalamus have significantly different node degrees with those in the HC group, while 8 subregions have lower degrees ( p < 0.01), including the one with LHb. These results implicate abnormal brain connections associated with the thalamus and LHb to be associated with SD. Integration of these connections by machine learning can provide connectome-based biomarkers to accurately diagnose SD.
Introduction
On the depression severity continuum, subclinical depression (SD) is a mild condition considered to be the precursor to major depressive disorder (MDD) (1, 2). Subjects with SD are very vulnerable to depression and are apt to generate suicide ideation (3, 4). The increasingly high incidence of SD among both college students and the elderly (estimated as high as 15%) clearly demonstrates the need for intensive investigation (5–7). Unfortunately, knowledge of neural substrates of SD is incomplete, making it difficult to identify reliable diagnostic biomarkers and take preventative treatments (8).
Some dysfunctional brain regions and connections have been evaluated in order to identify new biomarkers for SD. Via resting state fMRI (rs-fMRI), we have previously demonstrated that the altered spontaneous neuronal activity by measurement of amplitude of low-frequency fluctuations (ALFF) and disrupted functional connectivity (FC) are implicated in SD (9–11). We also found that SD presents the increased interhemispheric FC and cortical degree centrality, as well as decreased subcortical degree centrality. These measures differentiate SD subjects from healthy controls (HCs) (10–12). SD is characterized by changed FCs between subregions of the anterior cingulate cortex (ACC), increased FC of Hb within default model network regions, and decreased FC within salience network regions (8, 13). Kaiser et al. (14) demonstrated that there exists a high correlation between the neural activity of dorsal anterior cingulate cortex (dACC) and posterior cingulate cortex (PCC) in SD subjects, indicating that SD subjects are confronted with greater difficulty of shifting out of internally directed and ruminative thinking. Dedovic et al. (15) and Petrican et al. (16) reported the weaker functional dominance in dorsal attention network (DAN) [low connectivity between the superior parietal lobule (SPL) and the frontoparietal control network].
Network neuroscience explores interactions of different neurobiological element from an integrative perspective and is capable of providing with better predictive biomarkers for brain disorders by machine learning (17, 18). Machine learning is suitable for individual-level prediction from a prospective viewpoint, and it is a potentially powerful tool for precision psychiatry (19). For example, Support Vector Machine (SVM), as a typical method of machine learning, has been widely used to identify imaging biomarkers in diseases such as schizophrenia, major depression, bipolar disorder, etc. (20). For more information on machine learning and its application in psychiatry, one can refer to the comprehensive reviews (21–24). Recently, machine learning has proved useful to build connectome-based biomarkers for autism spectrum disorder, bipolar disorder, subtypes of depression, and schizophrenia (25–28). However, not many connectome-based biomarkers have been developed for SD.
Compared with SD, MDD has received more attention and significant breakthroughs have been achieved. For example, concrete evidence has demonstrated that bursts within the lateral habenula (LHb) drive depression in rats (29). As an evolutionary conserved epithalamic structure, LHb is involved in negative motivational value and decision-making (30–33). LHb is also considered to be the pathophysiological basement of MDD (34, 35). For more details on LHb, one can refer to these recent reviews (36–38). The deep brain stimulation of LHb has been successfully used to treat patients with refractory MDD (39). These findings on MDD may provide useful clues regarding SD.
LHb has been investigated by multimodal MRI in depressive and healthy subjects, but not in subjects with SD. LHb volume measured by high-resolution T1-weighted images decreases in depression, but not in posttraumatic stress disorder or schizophrenia (40, 41). Using task-based functional MRI (fMRI), Salas et al. (42) have shown that LHb is activated in response to negative reward prediction. It is worth noting that the fMRI study on LHb has several limitations. First, the habenula volume approximately ranges from 29 to 36 mm3 in each hemisphere based on structural MRI and postmortem measurement, which can be smaller than the voxel size of standard fMRI (40, 41, 43). The smoothing kernels [5–12 mm full width at half maximum (FWHM)] are larger in size than LHb. Second, the habenular signal is likely contaminated by adjacent structures, such as the medial dorsal thalamus or the epithalamic paraventricular nucleus (43).
Herein, connectome-based biomarkers are developed to predict subclinical depression through a machine learning algorithm and identify dysfunctional brain regions and connections. The method of predictive modeling used in our study is different with the traditional method of brain mapping (13). Predictive modeling integrates all brain data or features into a single prediction of outcome, making multiple comparisons unnecessary and increasing statistical power (18). Specifically, we parcellate the whole brain into 48 regions and 246 subregions using the latest human Brainnetome Atlas (44) and build large-scale resting-state functional brain networks using fMRI data. A two-sample t-test is used to identify initial candidate connections, and the resultant regional connections are input into SVM models as features. The performance of the predictive models is estimated by leave-one-out cross-validation. The node degree of subregions within the thalamus is compared between SD and HC groups. Connections linked with subregions of LHb are further investigated.
Materials and Methods
Participants
All the participants were enlisted from volunteers who had undergone health screening at Guangzhou Medical University from 2012 to 2014. The Beck Depression Inventory II (BDI-II) scale is administered to evaluate the depression symptom severity. Thirty-four subjects (11 males, 23 females) with BDI-II score >13 are placed into the SD group (BDI score mean ± SD: 22.58 ± 6.92) and 40 healthy controls (21 males, 19 females) are selected to match the SD group by age, sex, and education. According to the two-sample t-test, there is no significant difference for the age (years) between SD and HC groups (mean ± SD: 19.91 ± 1.64 vs. 19.70 ± 0.85, p = 0.50), neither for the education (years) (mean ± SD: 13.18 ± 0.58 vs. 13.08 ± 0.62, p = 0.47). By the chi-square test, there is no significant difference for the gender (p = 0.07). None of the participants fulfilled the criteria for MDD based on Diagnostic and Statistical Manual of Mental Disorders IV (DSM-IV). Other inclusion criteria for all participants include age ranging from 19 to 25 years, right-handedness, no visualized lesion on any MRI scans, no neurological illness, and no alcohol or drug dependence. The study is approved by the Medical Ethics Committee of Guangzhou First People’s Hospital of Guangzhou Medical University and is in accordance with the 1964 Helsinki declaration and its later amendments or comparable ethical standards. All participants signed a written informed consent in accordance with the Declaration of Helsinki (2000).
MRI Imaging Data Acquisition
All MRI images were acquired using one 3-Tesla MRI scanner (Siemens, Erlangen, Germany) with an eight-channel phase-array brain coil. Foam pads and headphone were utilized to minimize the head motion and reduce noise, respectively. As in our previous studies (10–12), high-resolution T1-weighted images were obtained with a standard magnetization prepared rapid gradient echo (MP-RAGE) sequence [repetition time (TR)/echo time (TE) 2,530/2.34 ms, flip angle (FA) 7°, field of view (FOV) 256 × 224 mm, slice thickness 1.0 mm]. The resting-state fMRI data were collected by one echo-planar imaging (EPI) sequence (TR/TE 2,500/21 ms, FA 90°, FOV 200 × 200 mm, matrix 64 × 64, 42 slices without gap, voxel size 3.5 × 3.1 × 3.1 mm). The images of 200 time points were collected, and the total amount of fMRI acquisition time is 500 s. During the resting-state fMRI scan, the participants were asked to relax, to close their eyes, not to think of anything in particular, and not to fall asleep. Wakefulness of participants has been confirmed immediately after the fMRI scanning session.
Study Design and Main Procedures
The study design and procedures are schematically shown in Figure 1. There are six steps for this study (Figure1A). After the first step of image preprocessing, functional brain networks for HCs and SDs are constructed. Then two-sample t-tests are used to identify potential dysfunctional connections. Three methods are proposed to further select connections from previously identified candidates. These selected connections are used to train and test the predictive models of SD. After excluding confounders such as the number of connections and p-value, dysfunctional brain regions and connections are determined by examining the models with high predictive accuracy. Finally the emphasis is placed on the dysfunctional thalamus and LHb. Abnormal connections associated with the thalamus and its subregions, including LHb, as well as the node degree of these subregions are characterized. These six steps are described in details below.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/441521/fpsyt-10-00371-HTML/image_m/fpsyt-10-00371-g001.jpg)
Figure 1 Study design and procedures. (A) Overview of the study procedures; (B) functional MRI (fMRI) image preprocessing; (C) construction of functional brain networks; (D) identification of dysfunctional connections; (E) connection selection and predictive models.
Functional MRI Image Preprocessing
As shown in Figure 1B, the T1-weighted and rs-fMRI data is preprocessed using the DPARSF toolbox (http://www.restfmri.net/forum/DPARSF) (45, 46) as follows. First, the initial 20 time points of raw fMRI data are removed in order to eliminate unstable factors. Second, the time layer correction and head movement correction are carried out. Third, the brain of each subject is registered to a normative template through spatial standardization. Fourth, a band-pass filtering of 0.01–0.1 Hz and Gaussian smoothing with 6 mm FWHM are implemented.
Construction of Functional Brain Networks
The procedure for constructing functional brain networks is shown in Figure 1C. First, the newly developed human Brainnetome Atlas is used to parcellate the whole brain into 48 brain regions (246 subregions). This atlas is four to five times as accurate as the traditional Brodmann map and has a more objective and accurate boundary (44). Each subregion represents a node in the constructed brain network. The time course of each subregion is calculated by averaging the time course of all voxels therein. The strength of functional connection or the connection weight (Wij), also identified as edge weight, is expressed as the Pearson correlation coefficient between the time courses of any two subregions (i, j). The correlation matrix is transformed into Z scores by applying Fisher’s r-to-Z transformation. For each individual, a weighted undirected network is obtained in the form of a 246 × 246 adjacency matrix (A). Given that it is controversial for interpreting negative correlation or functional connectivity (47, 48), the normalized absolute value of the matrix is used as done in previous studies (49, 50), such that 0 ≤ Wij ≤ 1 for all i and j.
Node degree and edge weight are used to determine whether a brain region (or subregion) is connected or dysfunctional in SD. The node degree (ki) refers to the number of connections that link this node to the rest of the network. For our weighted networks, the definition can be transformed as
where Wij is the strength of the connection between node i and node j, and N is the set of nodes in the network. The edge weight (Wij) is an important measure for evaluating the alteration in the strength of a connection in SD.
Identification of Dysfunctional Connections
As shown in Figure 1D, two-sample t-tests are performed to examine significant differences between edge weight in SD and HC groups (p < 0.05). For multiple comparisons, the false discovery rate (FDR) is controlled by the linear step-up procedure introduced by Benjamini and Hochberg (51). To avoid the information leakage, the two-sample t-test is carried out after leaving one out, not for all subjects. This step generated 74 different masks of abnormal connections. Based on each mask, the work in Connection Selection and Predictive Models is done. However, for the group study in Dysfunctional Thalamus and Lateral Habenula, the two-sample t-test is done for all subjects.
Connection Selection and Predictive Models
The study uses the library for support vector machines (LIBSVM) toolkit developed by Professor Lin of Taiwan University (https://www.csie.ntu.edu.tw/∼cjlin/libsvm/), which integrates many functions such as kernel selection, parameter adjustment, and prediction. For the training of SVM classification model, the radial basis kernel function (RBF) is used. This kernel function provides a good classification for samples with nonlinear relationship between labels and features as expressed below.
According to the recommendation from LIBSVM, the values of the optimal penalty coefficient C and the kernel function parameter γ are determined by the way of “grid-search” using cross-validation (52). After going through all the pairs of (C, γ) with C = 2−5, 2−3,…, 215 and γ = 2−15, 2−13,…, 23, the pair leading to the best cross-validation accuracy is found.
Three methods are proposed to select connections from identified candidates, and these selected connections are used to train and test the predictive SVM models of SD, as shown in Figure 1E. For the first method, the significantly altered connections associated with each brain region defined by the human Briannetome Atlas are used as input features. A total of 24 SVM models are built to predict SD, and they are named as anatomical-region-based models. The performance of these models is estimated through leave-one-out cross-validation using measures of accuracy, sensitivity, confusion matrix, receiver operating characteristic (ROC) curve, and the area under the curve (AUC). These models are ranked by accuracy. The brain regions leading to the models with an accuracy >0.90 are considered to be dysfunctional.
To determine whether models using connections associated with subregions not belonging to one specific anatomically well-defined brain region and with subregions that are anatomically nonadjacent can achieve comparable performance to the anatomical-region-based models, two more independent experiments are conducted. First, the method of sliding window with 16 subregions is employed to generate different input features and models. The reason why the number of subregions was set as 16 is that the thalamus leading to the predictive model of the highest ACC owns 16 subregions. As shown in the middle column of Figure 1E, to slide the window row by row throughout the adjacency matrix (246 × 246) will generate 231 windows (246 – 16 + 1 = 231). The models using the connections within each individual window as input features are named as sliding-window-based models. Second, a model is constructed using the functional connections within 16 randomly selected subregions as input features. A total of 100 similar models are generated and identified as random-selection-based models. The accuracy values of these three categories of models are compared.
Exclusion of Confounders
To estimate whether the performances of the anatomical-region-based models are dependent on the number of connections associated with the brain region and the p-value of these connections, their correlation coefficients are assessed. The distribution of connections in the model with the highest accuracy is investigated to explore whether these models with good performance are independent.
Dysfunctional Thalamus and Lateral Habenula
In order to further identify dysfunctional subregions and connections, the connections of brain regions achieving the highest predictive accuracy (thalamus) are examined. The number and p-value of connections associated with each of the 16 subregions are identified. Finally, the node degree of each subregions is compared between SD and HC groups.
Results
Anatomical-Region-Based Models
The 24 anatomical-region-based models ranked by the accuracy of predicting SD are presented in Figure 2A. The accuracy ranges from 0.65 to 0.92. The top five models used connections associated with the regions of thalamus, posterior superior temporal sulcus, cingulate gyrus, superior parietal lobule, and superior frontal gyrus. The accuracy of each of these five models is higher than 0.90. The anatomical locations are shown in Figure 2B. The ROC curves and the AUC values are shown in Figure 2C. The cingulate model achieves the highest AUC of 0.957. The thalamus model yields the second highest AUC of 0.943. The confusion matrices of the top five anatomical-region-based models are listed in Table 1. For the thalamus model, 31 out of 34 subjects with SD (91.2%, also defined as sensitivity) and 37 out of 40 HCs (92.5%, also defined as specificity) are predicted accurately. The posterior superior temporal sulcus model yields the highest specificity of 95.0%, and the posterior superior temporal sulcus model yields the highest sensitivity.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/441521/fpsyt-10-00371-HTML/image_m/fpsyt-10-00371-g002.jpg)
Figure 2 The performances of predictive models of subclinical depression (SD) and their comparison. (A) The prediction accuracy of 24 anatomical-region-based models; (B) the brain regions leading to the top five accuracy models; (C) receiver operating characteristic (ROC) curves and area under the curve (AUC) of the top five models; (D) comparison of the accuracy of models using connections with thalamus, 24 anatomical-region-based models, 231 sliding-window-based models, and 100 random-selection-based models.
Other Subregion Selection Strategies
The accuracy of the models using other subregion selection strategies is compared with that of the anatomical-region-based models, as shown in Figure 2D. No significant difference in accuracy is observed between the 231 sliding-window-based models and the 24 anatomical-region-based models (0.81 ± 0.06 vs. 0.80 ± 0.08). The accuracy of the 100 random-selection-based models is 0.45 ± 0.06, which is significantly lower than that of the anatomical-region-based models and the sliding-window-based models (p < 0.001). The top five anatomical-region-based models, and in particular the thalamus model, achieve extraordinarily higher accuracy, as compared with the other models. Two important features are worthy to be noted. First, the brain regions involved in the top five anatomical-region-based models are potentially dysfunctional due to SD. Second, the arbitrary anatomically adjacent subregions (obtained by the sliding window method) can generate comparable prediction accuracy with anatomically well-defined subregions (obtained by the anatomical-region-based method), but the randomly selected subregions cannot reliably predict SD.
Effect of the Number of Connections and the p-value
The number of connections associated with the 24 brain regions used for the predictive models of SD ranges from 52 to 240, as shown in Figure 3A. The top five models, which correspond to the regions of thalamus, posterior superior temporal sulcus, cingulate gyrus, superior parietal lobule, and superior frontal gyrus, have connections with the average number of 85, 120, 83, 222, and 131, respectively. Figure 3B shows the mean p-value of connections in individual brain regions, which ranges from 0.025 to 0.031.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/441521/fpsyt-10-00371-HTML/image_m/fpsyt-10-00371-g003.jpg)
Figure 3 Exclusion of confounders (the number of connections and the mean p-value). (A) The number of connections with significant difference between healthy controls (HCs) and SDs for each of 24 brain regions; (B) the p-value for significant difference of connection weight between HCs and SDs for each of 24 regions; (C) the relationship between the accuracy of prediction and the number of connections with significant difference (left part), between the accuracy of prediction and the mean of p-value for significant difference of connection weight (right part).
The dependence of the accuracy of the predictive models on the number of connections associated within individual brain regions and the mean p-value of connections are shown in Figure 3C. The accuracy of the predictive models is related to the number of connection associated with brain regions (r = 0.093), but there is no statistical significance (p = 0.665). When the number of connections associated with brain regions is randomly reduced to 90, as with the thalamus, the accuracy does not increase. For instance, after reducing the number of connections in middle frontal gyrus from 229 to 90, the accuracy decreases from 0.89 to 0.70 (the mean of 100 times random samples). The accuracy of the predictive models is negatively related to the mean p-value of connections associated with brain regions without statistical significance (r = 0.169, p = 0.429). The high value of accuracy of the top five models is due to neither the large number of connections, nor the small p-value of the connections.
Dysfunctional Connections With the Thalamus
Given that the thalamus is seen as one possible dysfunctional brain region of SD, connections associated with the thalamus are investigated, as shown in Figure 4. The number of connections between the thalamus and precuneus, insular gyrus, paracentral gyrus, and amygdala is higher than that of the other regions (11, 8, 8, and 8, respectively). However, the number of connections between the thalamus and itself, the posterior superior temporal sulcus, cingulate gyrus, superior parietal lobule, and superior frontal gyrus is only 6, 5, 3, 4, and 4, respectively. High accuracy values of the models using connections associated with the posterior superior temporal sulcus, cingulate gyrus, superior parietal lobule, and superior frontal gyrus are independent on the thalamus. These regions may also be impacted by SD.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/441521/fpsyt-10-00371-HTML/image_m/fpsyt-10-00371-g004.jpg)
Figure 4 The number of connections that present significant difference and connect the thalamus and each of 24 brain regions.
Subregions Within the Thalamus and Lateral Habenula
The distribution of the 90 significantly different connections associated with the thalamus among 16 subregions is shown in Figure 5A. There are 18 connections associated with the right posterior parietal thalamus (PPtha_r), much higher than those connected with the other regions. The significant asymmetry is observed, i.e., the right side has more connections than the left. Astonishingly, only two edges are connected to the left posterior parietal thalamus. The p-value of connections associated with PPtha_r is smaller than that of PPtha_l, as illustrated in the right part of Figure 5A. Based on Montreal Neurological Institute (MNI) coordinates, LHb is located in the posterior parietal thalamus (Figure 5B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/441521/fpsyt-10-00371-HTML/image_m/fpsyt-10-00371-g005.jpg)
Figure 5 Dysfunctional thalamus and lateral habenula. (A) The number of dysfunctional connections for16 subregions of thalamus (the left part) and the p-value of the dysfunctional connections with 16 subregions of thalamus (the right part); (B) the anatomical atlas of thalamus and LHb.
Node Degree of Subregions Within the Thalamus
The node degrees of 16 subregions within the thalamus are compared between SD and HC groups, as shown in Figure 6. Significant difference is found for 10 subregions. For eight subregions, the node degree of SD is significantly smaller than that of HC. The node degree of subregions on the right is higher than that of subregions on the left, for both SD and HC groups.
Discussions
Sophisticated connectome-based brain biomarkers permit the association of brain measures with both subjective experiences and objective behaviors, leading to a reconceptualization of diagnoses of mental illness. Herein, we have built several reliable brain biomarkers (>0.9 accuracy) that predict SD using abnormal functional connections as input features and SVM as the machine learning algorithm. We have found dysfunctional brain regions, especially the thalamus and LHb, which may be the etiological origin of SD. We have observed a reduction of the node degree for the right LHb in SD, but not for the left. The significance of these findings and the related advantages of this methodology are interpreted and discussed in the following subsections.
Reliable Biomarkers for Subclinical Depression Prediction
In this study, we have identified reliable brain biomarkers for SD prediction through the large-scale brain networks driven from resting state fMRI and a machine learning algorithm. Previously we had constructed biomarkers using the degree of centrality of different brain regions. The highest AUC was 0.82 for the right posterior parietal lobule (12). Here, the biomarkers are more reliable, and the highest AUC of 0.957 is achieved while using connections with the cingulate gyrus. Moreover, the arbitrary anatomically adjacent subregions (obtained by the sliding window method) and the anatomically well-defined subregions (obtained by the anatomical-region-based method) produce models with similar performances. These models present significantly higher accuracy than those driven by the randomly selected subregions. These results suggest that anatomical adjacency is important in the selection of feature (or connections) while building brain models or biomarkers. However, more sophisticated algorithms of feature selection, such as L1-regularized sparse canonical correlation analysis (L1-SCCA), may yield better biomarkers than anatomical adjacency (18, 25).
Dysfunctional Brain Regions in Subclinical Depression
Using the criterion of owing prediction accuracy greater than 0.90, we identified the top five regions associated with dysfunction in SD from 24 cortical and subcortical regions, defined by the human Brainnetome Atlas (44). These regions include the thalamus, posterior superior temporal sulcus, cingulate gyrus, superior parietal lobule, and superior frontal gyrus. Most of these regions had been reported in previous studies of SD. The related findings for each dysfunctional region are described below.
Not unexpectedly, given that the thalamus has multiple functions of relaying information between different subcortical regions and the cerebral cortex, the dysfunctional thalamus is identified in SD. It had been reported that two subtypes of depression had hyperconnectivity between the thalamic and frontostriatal network, resulting in symptoms related to reward processing, adaptive motor control, and action initiation (27, 53). Of particular importance, LHb, a small epithalamic structure, is believed to control reward and aversion processing. The importance of these observations will be discussed in detail below.
Distinct connectivity patterns are observed in subregions of the posterior cingulate cortex for SD (8). The anterior cingulate cortex (ACC) is an important component of reward circuitry, with abnormalities resulting in anhedonia (loss of interest/pleasure), a core symptom in MDD (54). Abnormal ACC is also linked to default model network (DMN, self-related thoughts), hyperconnectivity, and switching between the DMN and the central executive network (CEN, externally-focused cognition) (8, 14, 55, 56).
Previously we had reported that the superior parietal lobule (SPL, Brodmann area 7, BA 7) presented the decreased fractional ALFF (fALFF) (9). The SPL had been proposed to be the key component controlling the executive network and playing a critical role in working memory (57).
The superior frontal gyrus includes the dorsolateral prefrontal cortex (DLPFC) and the medial prefrontal cortex (MPFC). In depression, DLPFC is used for emotion adjustment, with the activity of DLPFC inhibited at rest but increased during symptom remission (58, 59). Our previous work had shown that the functional connectivity between SPL and DLPFC was reduced in SD (9). MPFC is an important component of DMN playing a crucial role in self-referential processing. A lack of DMN inhibition, i.e., self-focus, is a core issue of MDD (60). Most importantly, both regions of the superior frontal gyrus had been the targets for repetitive transcranial magnetic stimulation (rTMS) in depression treatment (61).
Dysfunctional Brain Regions Connected With the Thalamus
We have found that the dysfunctional thalamus in SD is mainly linked with the precuneus, insular gyrus, paracentral lobule, and amygdala. It is not surprising to observe the insula and amygdala because they are the neuroanatomical core of MDD pathology and closely related to anxiety (27). The precuneus is related to anhedonia, and the paracentral lobule (premotor) to anxiety. Positive connectivity between the LHb (an epithalamic structure) and the sensorimotor cortex had been reported by Ely et al. (13).
There is no overlap between the four regions with large number of abnormal connections associated with thalamus and the four regions (except thalamus) of the top five models ranked by accuracy. This is explained by the difference between machine learning and classical statistics which is discussed below.
Lateral Habenula—Beyond a Reasonable Doubt
Only one previous study had investigated the resting state functional connectivity of the LHb in SD (13). Herein, the SD group shows greater LHb connectivity with DMN and lower connectivity with the salience network, which is consistent with prior finding in MDD. Here we found a lateralized decrease of node degree in the subregion with LHb (at right) in SD. This finding is consistent with our previous finding of decreased subcortical degree centrality (12). We speculate that the decreased node degree of LHb corresponds to its hyperactivity and abnormal bursts. Given that LHb has an inhibitory effect on dopamine neurons, Hikosak (32) proposed that the hyperactivity of LHb results in hypoactivity of dopamine neurons, reducing motor activity in MDD. Hyperactivity of LHb could be the result of bursts. According to Yang et al. (29), LHb burst firing increases in depression and LHb bursts lead to depression in rats. Interestingly, Ely et al. (13) found that LHb connectivity increased in the left and decreased in the right. This may partly explain the decreased node degree of the right LHb.
Among many dysfunctional brain regions, which one is the most likely etiological origin of SD? Is it LHb, as in depression for rats (29)? Given the correlative nature of resting-state fMRI, it is difficult to establish a causal inference (58). Therefore, we do not know which brain region is the cause or consequence of SD. However, based on these evidences given in our study, we believe that LHb may be the origin of SD, beyond a reasonable doubt.
Machine Learning and Classical Statistics
In this study, we have used the classical statistical method, the two-sample t-test, to initially screen the candidate connections with significant difference between SD and HC groups. The selected connections are input into the SVM models as features. The two-sample t-test is actually used as a one feature selection algorithm. The method of using group tests has been proven to yield an inflated bias (62). Thus, we did not carry out strict multiple comparison corrections. More powerful feature selection or dimension reduction algorithms, such as L1-regularized sparse canonical correlation analysis, linear elastic-net, and minimum-redundancy maximum relevancy, can be assessed in the future (25, 63).
Moreover, we found that connections with small p-values do not always lead to high prediction accuracy in machine-learning-based models. This is consistent with previous studies, and originates from the essential difference between group difference and classification (62, 64, 65).
Limitations and Future Works
The sample size is relatively small even though a large population is screened. The prediction accuracy may decrease with the sample size for most disorders (62). The generalizability of these models needs to be validated in independent cohorts. The 50 subjects (25 with high SD scores and 25 with low scores) from the WU-Minn Human Connectome Project (HCP) Consortium’s 500 Subjects Release (66), previously used by Ely et al. (13), will be used in the future.
Only the criterion of BDI-II score >13 is relatively simple compared with the complicated clinical evaluation of SD. However, there is no consensus among researchers regarding how to combine different scales to define SD. Besides the widely applied BDI, the National Institutes of Health (NIH) Toolbox Negative Affect Survey Sadness Subscale, and the Achenbach Adult Self-Report (ASR) have been utilized (13, 16). Regression model of BDI-II score may tackle the issue induced by the fixed threshold.
We only distinguished subjects with SD from HC. However, it is probable that SD has several neurophysiological subtypes, as shown in depression (27). With an increased sample size, it may be possible to differentiate these subtypes. Moreover, generalizing our identified biomarkers to distinguish between SD and MDD cannot be done before the specific verification. The study done by Lawson et al. (67) demonstrated that habenula function is also disrupted in MDD. In our group, the data of MDD is being collected so that the verification will be carried out in the near future.
The voxel size of our fMRI images is 3.5 × 3.1 × 3.1 mm, which is larger than the volume of LHb, approximately 18.5 mm3 per hemisphere (41). As a consequence the precise location of LHb is difficult to determine. The fMRI data with the whole brain coverage and high resolution (2 mm isotropic) is available for 25 subjects with high SD from the HCP (66). Using those high-resolution data, the node degree for LHb and its lateralization may be investigated clearly.
Actually, we have identified the dysfunctional brain regions from the predictive models driven by dysfunctional connections and especially identified LHb as the possible etiological origin of SD. To find the predictive and dysfunctional functional connections is also important. However, there are too many predictive functional connections in the current work to explain one by one. For example, even only for the right posterior parietal thalamus, 18 abnormal functional connections are available. Compared with the connections, the findings on brain regions are more reliable. Once the number of predictive functional connections is reduced to be less than 20 as done by Yahata et al. (25) their explanations will be feasible and valuable.
Conclusion
Through integrating functional brain connections and SVM, this study has provided the connectome-based biomarkers for accurate prediction of SD. Abnormal brain connections with the thalamus and LHb are found to be implicated with SD. The right LHb in SD shows the decreased node degree comparing with HC group, but the left LHb does not. This evidence indicates that LHb may be the etiological origin of SD. The generated biomarkers can aid early diagnosis of SD. Furthermore, the identified dysfunctional brain connections and regions may help localize the etiological origin of SD and understand the pathogenesis of SD.
Ethics Statement
This study was carried out in accordance with the recommendations of the guidelines of the Declaration of Helsinki, the Medical Ethics Committee of Guangzhou First People’s Hospital of Guangzhou Medical University. All subjects gave written informed consent in accordance with the Declaration of Helsinki. The protocol was approved by the Medical Ethics Committee of Guangzhou First People’s Hospital of Guangzhou Medical University.
Author Contributions
SQ, JH, and XW designed and directed the study. YZ, SQ, BZ, DH, YT, and JH analyzed data. XW recruited participants and acquired data. YZ, SQ, and YT drafted the manuscript. All authors revised and approved the final version of the manuscript.
Funding
The Fundamental Research Funds for the Central Universities (N181904003, N172008008), the National Science Foundation of China (81871846), and the Science and Technology Planning Project of Guangzhou (201804010032).
Conflict of Interest Statement
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.
References
1. Lewinsohn PM, Klein DN, Durbin EC, Seeley JR, Rohde P. Family study of subthreshold depressive symptoms: risk factor for MDD? J Affect Disord (2003) 77(2):149–57. doi: 10.1016/S0165-0327(02)00106-4
2. Shankman SA, Lewinsohn PM, Klein DN, Small JW, Seeley JR, Altman SE. Subthreshold conditions as precursors for full syndrome disorders: a 15-year longitudinal study of multiple diagnostic classes. J Child Psychol Psychiatry (2009) 50(12):1485–94. doi: 10.1111/j.1469-7610.2009.02117.x
3. Fergusson DM, Horwood LJ, Ridder EM, Beautrais AL. Subthreshold depression in adolescence and mental health outcomes in adulthood. Arch Gen Psychiatry (2005) 62:66–72. doi: 10.1001/archpsyc.62.1.66
4. Cukrowicz KC, Schlegel EF, Smith PN, Jacobs MP, Van Orden KA, Paukert AP, et al. Suicide ideation among college students evidencing subclinical depression. J Am Coll Health (2011) 59(7):575–81. doi: 10.1080/07448481.2010.483710
5. Lavretsky H, Kumar A. Clinically significant non-major depression: old concepts, new insights. Am J Geriatr Psychiatry (2002) 10:239–55. doi: 10.1176/appi.ajgp.10.3.239
6. VanItallie TB. Subsyndromal depression in the elderly: underdiagnosed and untreated. Metab Clin Exp (2005) 54:39–44. doi: 10.1016/j.metabol.2005.01.012
7. Mikolajczyk RT, Maxwell AE, ElAnsari W, Naydenova V, Stock C, Ilieva S, et al. Prevalence of depressive symptoms in university students from Germany, Denmark, Poland and Bulgaria. Soc. Epidemiol Psychiatr Sci (2008) 43:105–12. doi: 10.1007/s00127-007-0282-0
8. Philippi CL, Motzkin JC, Pujara MS, Koenigs M. Subclinical depression severity is associated with distinct patterns of functional connectivity for subregions of anterior cingulate cortex. J Psychiatr Res (2015) 71:103–11. doi: 10.1016/j.jpsychires.2015.10.005
9. Wei X, Shen H, Ren J, Li X, Xu X, Yang R, et al. Altered resting-state connectivity in college students with nonclinical depressive symptoms. Plos One (2014) 9(12):e114603. doi: 10.1371/journal.pone.0114603
10. Wei X, Shen H, Ren J, Liu W, Yang R, Liu J, et al. Alteration of spontaneous neuronal activity in young adults with non-clinical depressive symptoms. Psychiat Res Neuroim (2015) 233(1):36–42. doi: 10.1016/j.pscychresns.2015.04.008
11. Wei X, Ren J, Liu W, Yang R, Xu X, Liu J, et al. Increased interhemispheric functional connectivity in college students with non-clinical depressive symptoms in resting state. Neurosci (2015) 589:67–72. doi: 10.1016/j.neulet.2015.01.034
12. Gao C, Liu W, Liu Y, Ruan X, Chen X, Liu L, et al. Decreased subcortical and increased cortical degree centrality in a nonclinical college student sample with subclinical depressive symptoms: a resting-state fMRI study. Front Hum Neurosci (2016) 10:1–9. doi: 10.3389/fnhum.2016.00617
13. 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. doi: 10.1002/hbm.23179
14. Kaiser RH, Andrews-Hanna JR, Spielberg JM, Warren SL, Sutton BP, Miller GA, et al. Distracted and down: neural mechanisms of affective interference in subclinical depression. Soc Cogn Affect Neurosci (2015) 10(5):654–63. doi: 10.1093/scan/nsu100
15. Dedovic K, Giebl S, Duchesne A, Lue SD, Andrews J, Efanov S, et al. Psychological, endocrine, and neural correlates of attentional bias in subclinical depression. Anxiety Stress Coping (2016) 29(5):479–96. doi: 10.1080/10615806.2015.1101457
16. Petrican R, Saverino C, Rosenbaum RS, Grady C. Inter-individual differences in the experience of negative emotion predict variations in functional brain architecture. Neuroimage (2015) 123:80–8. doi: 10.1016/j.neuroimage.2015.08.031
18. Woo CW, Chang LJ, Lindquist MA, Wager TD. Building better biomarkers: brain models in translational neuroimaging. Nat Neurosci (2017) 20(3):365–77. doi: 10.1038/nn.4478
19. Bzdok D, Meyer-Lindenberg A. Machine learning for precision psychiatry: opportunities and challenges. Biol Psychiatry Cogn Neurosci Neuroim (2018) 3(3):223. doi: 10.1016/j.bpsc.2017.11.007
20. 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(4):1140–52. doi: 10.1016/j.neubiorev.2012.01.004
21. Iniesta R, Stahl D, Mcguffin P. Machine learning, statistical learning and the future of biological research in psychiatry. Psychol Med (2016) 46:2455–65. doi: 10.1017/S0033291716001367
22. Vieira S, Pinaya WHL, Mechelli A. Using deep learning to investigate the neuroimaging correlates of psychiatric and neurological disorders: methods and applications. Neurosci Biobehav Rev (2017) 74:58–75. doi: 10.1016/j.neubiorev.2017.01.002
23. Gao S, Calhoun VD, Shui J. Machine learning in major depression: from classification to treatment outcome prediction. CNS Neurosci Ther (2018) 24:1037–52. doi: 10.1111/cns.13048
24. Dwyer DB, Falkai P, Koutsouleris N. Machine learning approaches for clinical psychology and psychiatry. Annu Rev Clin Psychol (2018) 14:91–118. doi: 10.1146/annurev-clinpsy-032816-045037
25. Yahata N, Morimoto J, Hashimoto R, Lisi G, Shibata K, Kawakubo Y, et al. A small number of abnormal brain connections predicts adult autism spectrum disorder. Nat Commun (2016) 7:11254. doi: 10.1038/ncomms11254
26. Librenzagarcia D, Kotzian BJ, Yang J, Mwangi B, Cao B, Lima LNP, et al. The impact of machine learning techniques in the study of bipolar disorder: a systematic review. Neurosci Biobehav Rev (2017) 80:538–54. doi: 10.1016/j.neubiorev.2017.07.004
27. 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(1):28–38. doi: 10.1038/nm.4246
28. Sui J, Qi S, van Erp TGM, Bustillo J, Jiang R, Lin D, et al. Multimodal neuromarkers in schizophrenia via cognition-guided MRI fusion. Nat Commun (2018) 9:3028. doi: 10.1038/s41467-018-05432-w
29. 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. doi: 10.1038/nature25509
30. Matsumoto M, Hikosaka O. Lateral habenula as a source of negative reward signals in dopamine neurons. Nature (2007) 447:1111–5. doi: 10.1038/nature05860
31. Lawson RP, Seymour B, Loh E, Lutti A, Dolan RJ, Dayan P, et al. The habenula encodes negative motivational value associated with primary punishment in humans. Proc Natl Acad Sci USA (2014) 111(32):11858–63. doi: 10.1073/pnas.1323586111
32. Hikosaka O. The habenula: from stress evasion to value-based decision-making. Nat Rev Neurosci (2010) 11(7):503–13. doi: 10.1038/nrn2866
33. 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. doi: 10.1038/nn.3587
34. 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. doi: 10.1038/nature09742
35. Stamatakis AM, Stuber GD. Activation of lateral habenula inputs to the ventral midbrain promotes behavioral avoidance. Nat Neurosci (2012) 15:1105–7. doi: 10.1038/nn.3145
36. Proulx CD, Hikosaka O, Malinow R. Reward processing by the lateral habenula in normal and depressive behaviors. Nat Neurosci (2014) 17:1146–52. doi: 10.1038/nn.3779
37. Benarroch EE. Habenula: recently recognized functions and potential clinical relevance. Neurology (2015) 85(11):992–1000. doi: 10.1212/WNL.0000000000001937
38. Boulos LJ, Darcq E, Kieffer BL. Translating the habenula—from rodents to humans. Biol Psychiatry (2016) 81(4):296–305. doi: 10.1016/j.biopsych.2016.06.003
39. Sartorius A, Kiening KL, Kirsch P, von Gall CC, Haberkorn U, Unterberg AW, et al. Remission of major depression under deep brain stimulation of the lateral habenula in a therapy-refractory patient. Biol Psychiatry (2010) 67:e9–e11. doi: 10.1016/j.biopsych.2009.08.027
40. 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. doi: 10.1016/j.biopsych.2010.09.027
41. 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. doi: 10.1017/S0033291709990821
42. Salas R, Baldwin P, De Biasi M, Montague R. BOLD responses to negative reward prediction errors in the human habenula. Front Hum Neurosci (2010) 36(4). doi: 10.3389/fnhum.2010.00036
43. Lawson RP, Drevets WC, Roiser JP. Defining the habenula in human neuroimaging studies. Neuroimage (2013) 64:722–7. doi: 10.1016/j.neuroimage.2012.08.076
44. 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. doi: 10.1093/cercor/bhw157
45. Yan C, Zang Y. DPARSF: a MATLAB toolbox for “pipeline” data analysis of resting-state fMRI. Front Syst Neurosci (2010) 4:13. doi: 10.3389/fnsys.2010.00013
46. Cohen JD, Daw N, Engelhardt B, Hasson U, Li K, Niv Y, et al. Computational approaches to fMRI analysis. Nat Neurosci (2017) 20:304–13. doi: 10.1038/nn.4499
47. Weissenbacher A, Kasess C, Gerstl F, Lanzenberger R, Moser E, Windischberger C. Correlations and anticorrelations in resting-state functional connectivity MRI: a quantitative comparison of preprocessing strategies. Neuroimage (2009) 47:1408–16. doi: 10.1016/j.neuroimage.2009.05.005
48. Murphy K, Birn RM, Handwerker DA, Jones TB, Bandettini PA. The impact of global signal regression on resting state correlations: are anti-correlated networks introduced? Neuroimage (2009) 44:893–905. doi: 10.1016/j.neuroimage.2008.09.036
49. Wang J, Zuo X, Gohel S, Milham MP, Biswal BB, He Y. Graph theoretical analysis of functional brain networks: test–retest evaluation on short- and long-term resting-state functional MRI data. PLoS One (2011) 6:e21976. doi: 10.1371/journal.pone.0021976
50. Wang J, Wang X, Xia M, Liao X, Evans A, He Y. GRETNA: a graph theoretical network analysis toolbox for imaging connectomics. Front Hum Neurosci (2015) 9:386. doi: 10.3389/fnhum.2015.00458
51. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Statist Soc B (1995) 57:289–300. doi: 10.2307/2346101
52. Chang CC, Lin CJ. LIBSVM: a library for support vector machines. ACM Trans Intell Syst Technol (2011) 2(27):1–27. doi: 10.1145/1961189.1961199
53. Ferenczi EA, Zalocusky KA, Liston C, Grosenick L, Warden MR, Amatya D, et al. Prefrontal cortical regulation of brainwide circuit dynamics and reward-related behavior. Science (2016) 351:aac9698. doi: 10.1126/science.aac9698
54. Pujara M, Koenigs M. Mechanisms of reward circuit dysfunction in psychiatric illness prefrontalestriatal interactions. Neuroscience (2014) 20:82e95. doi: 10.1177/1073858413499407
55. Whitfield-Gabrieli S, Ford JM. Default mode network activity and connectivity in psychopathology. Annu Rev Clin Psychol (2012) 8:49e76. doi: 10.1146/annurev-clinpsy-032511-143049
56. Hamilton JP, Furman DJ, Chang C, Thomason ME, Dennis E, Gotlib IH. Default-mode and task-positive network activity in major depressive disorder: implications for adaptive and maladaptive rumination. Biol Psychiatry (2011) 70:327e333. doi: 10.1016/j.biopsych.2011.02.003
57. Koenigs M, Barbey AK, Postle BR, Grafman J. Superior parietal cortex is critical for the manipulation of information in working memory. J Neurosci (2009) 29(47):14980–6. doi: 10.1523/JNEUROSCI.3706-09.2009
58. Koenigs M, Grafmanb J. The functional neuroanatomy of depression: distinct roles for ventromedial and dorsolateral prefrontal cortex. Behav Brain Res (2009) 201:239–43. doi: 10.1016/j.bbr.2009.03.004
59. Viviani R. Emotion regulation, attention to emotion, and the ventral attentional network. Front Hum Neurosci (2013) 7:746. doi: 10.3389/fnhum.2013.00746
60. Lemogne C, Delaveau P, Freton M, Guionnet S, Fossati P. Medial prefrontal cortex and the self in major depression. J Affect Disord (2012) 136(1–2):e1–e11. doi: 10.1016/j.jad.2010.11.034
61. Downar J, Daskalakis ZJ. New targets for rTMS in depression: a review of convergent evidence. Brain Stimul (2013) 6(3):231–40. doi: 10.1016/j.brs.2012.08.006
62. Arbabshirani MR, Plis S, Sui J, Calhoun VD. Single subject prediction of brain disorders in neuroimaging: promises and pitfalls. Neuroimage (2017) 145:137–65. doi: 10.1016/j.neuroimage.2016.02.079
63. Brown G, Pocock A, Zhao MJ, Lujan M. Conditional likelihood maximisation: a unifying framework for mutual information feature selection. J Mach Learn Res (2012) 13:27–66. doi: 10.1080/00207179.2012.669851
65. Bzdok D. Classical statistics and statistical learning in imaging neuroscience. Front Neurosci (2017) 11:543. doi: 10.3389/fnins.2017.00543
66. 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. doi: 10.1016/j.neuroimage.2013.05.041
Keywords: resting state functional MRI, brain network, subclinical depression, brain biomarker, functional connection, node degree
Citation: Zhu Y, Qi S, Zhang B, He D, Teng Y, Hu J and Wei X (2019) Connectome-Based Biomarkers Predict Subclinical Depression and Identify Abnormal Brain Connections With the Lateral Habenula and Thalamus. Front. Psychiatry 10:371. doi: 10.3389/fpsyt.2019.00371
Received: 07 December 2018; Accepted: 13 May 2019;
Published: 12 June 2019.
Edited by:
Jing Sui, Institute of Automation (CAS), ChinaReviewed by:
Rongtao Jiang, Institute of Automation (CAS), ChinaYegang Hu, Shanghai Mental Health Center (SMHC), China
Copyright © 2019 Zhu, Qi, Zhang, He, Teng, Hu and Wei. 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: Shouliang Qi, qisl@bmie.neu.edu.cn; Xinhua Wei, weixinhua@aliyun.com