Skip to main content
  • Research article
  • Open access
  • Published:

Functional connectivity between the thalamus and the primary somatosensory cortex in major depressive disorder: a resting-state fMRI study



Studies have confirmed that the thalamus and the primary somatosensory cortex (SI) are associated with cognitive function. These two brain regions are closely related in structure and function. The interactions between SI and the thalamus are of crucial significance for the cognitive process. Patients with major depressive disorder (MDD) have significant cognitive impairment. Based on these observations, we used resting-state functional magnetic resonance imaging (rs-fMRI) to investigate whether there is an abnormality in the SI-thalamic functional connection in MDD. Furthermore, we explored the clinical symptoms related to this abnormality.


We included 31 patients with first-episode major depressive disorder and 28 age-, gender-, and education-matched healthy controls (HC). The SI-thalamic functional connectivity was compared between the MDD and HC groups. The correlation analyses were performed between areas with abnormal connectivity and clinical characteristics.


Compared with healthy subjects, the MDD patients had enhanced functional connectivity between the thalamus and SI (p < 0.05, corrected). Brain areas with significantly different levels of connectivity had a negative correlation with the Assessment of Neuropsychological Status total score (r = − 0.383, p = 0.033), delayed memory score (r = − 0.376, p = 0.037) and two-digit continuous operation test score (r = − 0.369, p = 0.041) in MDD patients.


These results demonstrate that SI-thalamic functional connectivity is abnormal and associated with the core clinical symptoms in MDD. The SI-thalamic functional connectivity functions as a neurobiological feature and potential biomarker for MDD.

Peer Review reports


Major depressive disorder (MDD) is a prevalent psychiatric disorder [1]. MDD is characterized by affective disorder, vegetative symptoms, and cognitive function deficits. Clinicians and researchers have found that MDD patients are characterized by significantly poorer cognitive functions, including attention, executive functioning, episodic memory, and processing speed [2,3,4]. The evidence suggests that cognitive dysfunction is observed beyond the acute phase, and some of these impairments even persist during illness remission, which may lead to depressive relapse and dysfunction [5,6,7]. Researchers have proposed that functional and mood symptom recovery must include cognitive function recovery [8].

The thalamus is the hub of cortical-subcortical connections located deep in the brain and has been traditionally viewed as a simple relay station for neurons for all body sensations (except olfaction) to then project to the cerebral cortex [9]. The thalamus works as a tightly connected “central core” of brain regions for multiple tasks and behaviours [10]. Some evidence suggests that the thalamus combined with the dorsal anterior cingulate cortex and anterior insular compose the salience network, which has been described as having a central role in cognitive control [11, 12]. Recent research has shown that the thalamus can amplify and sustain cortical representations, and a new framework acknowledging thalamus-frontal circuits in cognition has been defined [13]. In disease conditions, increased thalamic functional connectivity was related to decreased cognitive function [14,15,16,17,18]. These findings suggest that the importance of the thalamus in cognitive activity is gradually being recognized.

The primary somatosensory cortex (SI) is located in the postcentral gyrus and structurally and functionally connected to the thalamus [19, 20]. SI is viewed as an acquisition and transformation sensory signal structure; meanwhile, involvement in controlling and modulating associatively learned behaviors is becoming apparent [21,22,23,24]. SI participated in the non-painful stimulation encoding process, and furthermore, SI maintained task-related tactile information in the late maintenance stage and contributed to the memory trace of a pain sensation [25,26,27,28]. A researcher even suggested SI as an “embodied mind” that makes the unconscious self exist in concert with the embodied facet of the self [29].

All of this evidence provides promising opportunities for understanding the importance of the thalamus and SI in cognitive function. These two brain regions are closely related in structure and function. The interactions between SI and the thalamus are of crucial significance for the cognitive process. In MDD, patients have cognitive impairments. Therefore, we hypothesized that abnormal SI-thalamic functional connectivity might be a neurobiological feature of MDD that is closely related to these clinical symptoms. Resting-state functional magnetic resonance imaging (rs-fMRI) is a powerful neuroimaging technique that enables researchers to measure spontaneous fluctuations in activity among distinct brain regions [30]. The method of using rs-fMRI to explore the brain’s intrinsic functional networks has been called resting-state functional connectivity (rs-FC) and has been used in numerous studies [31, 32]. Therefore, we used rs-FC as our experimental measure. This study may generate a new understanding of the underlying neurobiology depressive disorders.



Thirty-one patients with first-episode major depressive disorder were recruited from the First Hospital of Shanxi Medical University. At least two consultant psychiatrists diagnosed all patients according to the criteria for MDD in the “Diagnostic and Statistical Manual of Mental Disorders Fourth Edition (DSM-IV)” and the Chinese version of the Modified Structured Clinical Interview for DSM-IV TR Axis I Disorders Patient Edition (SCID-I/P). The patients were excluded if they met any of the following criteria: 1) had comorbid mental or neurological illnesses or personality disorders; 2) had used psychiatric drugs in the previous 2 weeks; 3) scored less than 17 on the 17-item Hamilton depressive scale (HAMD-17) administered by well-trained research assistants with backgrounds in psychology or psychiatry; 4) were pregnant; 5) were under 18 or over 60 years of age; 6) had a history of substance abuse or drug addiction; or 7) were left-handed or mixed-handed. We collected the patients’ general information, including name, gender, age, education level, occupation and family history of psychotic diseases.

Twenty-eight age-, gender-, and education-matched healthy controls were recruited for comparison. Participants were excluded if they met any of the following criteria: 1) had mental or neurological illnesses or personality disorders; 2) had a history of substance abuse or drug addiction; 3) had a family history of mental disorders; 4) were pregnant; 5) were under 18 or over 60 years of age; or 6) were left-handed or mixed-handed.

The Ethical Committee for Medicine of the First Hospital of Shanxi Medical University approved this study. Written informed consent was received from all participants prior to inclusion.

Clinical assessment

An increasing number of studies have shown that depression cannot be observed from the single point of view of depression but must be observed from the multiple perspectives of affective and somatization symptoms and cognition [33]. We used the HAMD-17 and the Snaith-Hamilton Pleasure Scale (SHAPS) to assess affective experience. We used the fatigue severity scale (FSS) to assess somatization symptoms. Finally, we used the Repeatable Battery for the Assessment of Neuropsychological Status (RBANS) and Continuous Performance Test (CPT) to assess cognition features. Our senior psychologist evaluated all the scales in this study.

The HAMD-17 is the most commonly used scale for assessing the severity of MDD, and there is a good degree of confidence in the scale’s validity and reliability [34].

The SHAPS is a 14-item scale, each item ranged from 0 to 3 (“strongly agree” to “disagree”), and it was filled out by the subject to evaluate the clinical absence of pleasure [35].

The FSS is a 9-item scale of fatigue severity, each item ranged from the lightest (i.e., 1) to the worst (i.e., 7), and evaluated clinically elevated fatigue [36].

The RBANS is a neuropsychological screening scale designed by Randolph in 1998 for ease of use and for rating neuropsychological function in 20 ~ 89-year-old people and has excellent reliability and validity [37]. The scale has 12 test items covering the five factors of immediate memory, visuospatial/constructional, language, attention, and delayed memory.

The CPT mainly examines the continuous attention level of the participants, the constant concentration on the response, and the level of arousal. For this task, different numbers are presented on the display screen in random order, asking the participants to identify when the same number appeared repeatedly and to react within a specified time. Based on the number of digits of the flashing figures, it is divided into two-digit, three-digit, and four-digit categories [38].

MRI acquisition and preprocessing of rs-fMRI data

MRI acquisition

A MAGNETOM Trio Tim 3.0 T scanner (Siemens Medical Solutions, Germany) with a 12-channel birdcage head coil located at the Shanxi Provincial People’s Hospital was used to acquire rs-fMRI. rs-fMRI was performed using an echo planar imaging (EPI) sequence with the following parameters: TR = 2000 ms, TE = 30 ms, flip angle = 70°, FOV = 24 × 24 cm, matrix = 64 × 64, slice gap = 2 mm, slice thickness = 2 mm, 6 min acquisition. During the resting functional scan, participants were instructed to keep their eyes closed and let their minds wander and not to fall asleep; all participants reported that they did not fall asleep.

Preprocessing of rs-fMRI data

Data preprocessing was conducted using the DPARSFA toolbox version 3.2, which was based on statistical parametric mapping 8 (SPM8) and the Resting-State fMRI Data Analysis Toolkit (REST) [39, 40]. The first 10 volumes of functional time points were discarded to allow the participants to adapt to the scanning noise. The remaining 170 volumes were preprocessed, which included the following: 1) slice timing, 2) realigning to reduce head motion (all head movements exceeding 2 mm were excluded), 3) spatial normalizing to the Montreal Neurological Institute (MNI) coordinate space with 3 × 3 × 3 mm, 4) spatial smoothing with a 6 × 6 × 6 full-width at half maximum (FWHM) kernel, 5) linear detrending, 6) temporal bandpass filtering (0.01–0.08 Hz), and 7) white matter signal, cerebrospinal fluid signal, global mean signal and six head motion parameters were used as covariates.

Regions of interest (ROI) were obtained with the WFU Pick Atlas 3.0.5, which automatically generates segmented atlas ROI templates in MNI space [41, 42]. The thalamus was defined as both “Thalamus_L” and “Thalamus_R”. SI was defined as “Brodmann areas 3, 1, 2”. Then, the image was resliced to the dimension of our functional images (voxel dimension: 3x3x3).

Connectivity between the thalamus and SI was examined using DPARSFA. The mean time course for the thalamus ROI was calculated by averaging the time course for all voxels within the thalamus ROI [43]. Calculations of the functional connection of the thalamus to the whole brain were made. The resulting correlation coefficients were transformed into z scores by Fisher’s z transform to create subject-specific maps [44]. Then, we performed a two-sample t-test across the two groups of subjects for the primary somatosensory cortex ROI.

Data analysis

All analyses were performed using IBM SPSS Statistics Version 23.0 (SPSS 23.0). Data are expressed as the means ± standard deviations (SD). Independent-sample t tests and χ2 tests were used to analyse the demographic data of the two groups. A two-sample t-test was conducted to compare the clinical symptoms scores across groups. P-values of 0.05 were considered statistically significant (two-tailed).

The rs-FC of the HC and MDD groups was contrasted using the DPARSFA toolbox version 3.2. The two-sample t-test was used to identify abnormalities between the two groups. To ensure that effects were not accounted for by other factors, such as age, sex and years of education, these variables were included as regressors of no interest. Multiple comparison corrections were determined by Monte Carlo simulation (1000 iterations) using the REST AlphaSim program [45]. Voxels with p < 0.01 and cluster size ≥4 voxels (108 mm3) were regarded as brain areas with a significant difference which confined within the SI mask, corresponding to corrected p < 0.05 [46]. From the preprocessed resting-state fMRI data, we extracted the mean time course of each of the statistically significant clusters.

Furthermore, we used Pearson’s correlation coefficients to evaluate any relationships between abnormal clusters and clinical features of MDD patients by SPSS 23.0. The results were considered statistically significant at p < 0.05.


Demographic data

As shown in Table 1, the mean age (mean ± SD) of the patients and the control group were 29.96 ± 9.73 and 26.79 ± 6.91 years, respectively. No differences in age, sex or years of education were observed between the patients and controls (p < 0.05). Table 1 also shows the HAMD-17, SHAPS, FSS, RBANS, and CPT scores for each group. All of these differences were significant between the patients and controls.

Table 1 Demographic and Clinical characteristics of participants

Functional connectivity analysis results

Compared with the HC group, the MDD group showed significant decreases in rs-FC between the thalamus and SI in two clusters. As shown in Fig. 1 and Table 2, both clusters were in the right central posterior gyrus: the first cluster had the maximally intense voxel at MNI coordinates x = 27, y = − 33, z = 75 (8 voxels, about 216 mm3, T = 3.6162, p < 0.05 AlphaSim correction), and the second cluster had the maximally intense voxel at MNI coordinates x = 45, y = − 30, z = 63 (8 voxels, about 216 mm3, T = 3.2803, p < 0.05 AlphaSim correction).

Fig. 1
figure 1

The images display the regions in the right central posterior gyrus that show increased functional connectivity from the thalamus in major depressive disorder (MDD) compared to healthy controls (HC) at rest. The colour bar represents the range of T values

Table 2 Differences in SI-thalamic functional connectivity between MDD and HC (MDD > HC)

Correlational analysis

As shown in Fig. 2, significantly different brain regions and patient clinical symptoms in the SHAPS total score (r = − 0.383, p = 0.033), delayed memory score (r = − 0.376, p = 0.037), and two-digit continuous operation test score (r = − 0.369, p = 0.041) showed significant negative correlations. The significantly different brain regions had no significant correlation with the other scores.

Fig. 2
figure 2

In MDD patients, brain areas significantly different levels of connectivity were negatively correlated with cognitive function, including SHAPS total score (r = − 0.383, p = 0.033), delayed memory score (r = − 0.376, p = 0.037), and two-digit continuous operation test score (r = − 0.369, p = 0.041) as shown in (a, b, and c), respectively


Analysis of fMRI has been widely performed to study the dysfunction of brain systems in patients with psychiatric disorders. SI-thalamic functional connectivity in MDD patients has not been investigated. Our study found enhanced SI-thalamic functional connectivity, which was associated with affective experience and cognition function.

The rs-FC between the thalamus and the whole brain has been previously examined, and there were relatively consistent findings. MDD subjects had altered thalamocortical connectivity characterized by an intricate pattern in the region that was associated with clinical symptoms; indeed, one of them indicated the thalamus is hyper-connected with SI [47,48,49,50]. The importance of the SI has been gradually discovered in schizophrenia patients. For example, its abnormal internal activity predicted changes in processing speed, and external connections with the thalamus, as part of the thalamocortical-cerebellar dysconnectivity associated with schizophrenia, have become a reliable neurobiological marker [51, 52]. Our study does not emphasize the functional abnormalities of these two independent brain regions and their relationship to the development of MDD. We proposed that the differences in the SI-thalamic functional connectivity are the characteristic changes, as the central components of the experience of MDD.

Anhedonia is the most important clinical manifestation of MDD. Our study showed that abnormal connections between the thalamus and SI leads to the loss of pleasure. Numerous studies have reported that the role of SI has always been controversial with regards to the emergence of pleasure. Early studies showed that SI played an essential role in tactile pleasure perception, while studies that have used physical stimuli to modulate touch pleasantness have failed to control for differences in physical intensity and the perceived power of the overstimulation [53,54,55]. The research has shown that pleasant experience comes from outside of SI, but it is closely related to SI and brain activity, which may involve the thalamus [56, 57]. Our study confirmed that hyper-connection between the thalamus and SI is closely related to emotional experience, and to a certain degree, provides an explanation for the lack of pleasure experienced in MDD patients. However, further precision in identifying the role of these areas requires careful exploration of brain neuronal activity.

The starting point of our research was that cognitive function was an essential function, and our study emphasized that SI-thalamic functional connectivity is indeed closely related to attention and memory in MDD. In our research, we not only stressed the reliability of our previous theoretical basis but also provided a more detailed analysis of the cognitive items related to the SI-thalamic functional connectivity in MDD patients.

The research highlights the importance of SI for memory in different diseases or conditions [58, 59]. However, there are still some contradictory conclusions where researchers have proposed that SI did not affect the metacognitive accuracy of either temporal or spatial tactile working memory [60]. Unfortunately, the impact of thalamus and SI connections on cognition has been less explored. Research has demonstrated that in schizophrenia, hyper-connectivity with the thalamus correlated positively with working memory [61]. Our study also showed that similar changes occurred in MDD patients. The structural connection between the thalamus and SI may partially explain why the functional connections between the two are related to cognition. Thalamocortical responses provide input to SI for establishing context and storing sensory memories (e.g., slowly changing body memories stored in layer 4 and sensory memories rapidly stored in layer 2/3) [62]. However, the pathophysiological process for enhanced functional connectivity has not been elucidated.

The thalamus has a function of adjusting attention [63]. The thalamus functions as a master regulator of functional cortical connectivity. Therefore, the construction of a directed arousal state is useful for attention, impacting how a cognitive process such as attentional control unfolds over cortical space and time [64,65,66]. In MDD patients, we tentatively proposed that the hyper-connectivity between the thalamus and SI may disrupt the balance of this network, allowing the time for concentration to be prolonged. However, a new study shows that thalamus engagement in a delay is not to relay specific rule information but to sustain existing cortical representations [67]. These studies prompted us to carefully consider the role of thalamus inattention, and the mechanism for the inattention needs further investigation.

Our study has several limitations. First, our study design did not allow us to investigate causality or the role of change in functional connectivity between the thalamus and SI on the development of the MDD patients’ cognitive impairment. Second, the thalamus could be divided into different subregions, but we did not subdivide when selecting regions of interest. Third, the heterogeneity of our patient population and a limited number of participants may lead to biased results. Fourth, long-term follow-up observations can lead to more profound experimental results. Fifth, we did not consider the impact of the length of the initial duration of the disease on brain function connectivity.


In conclusion, we have implicated enhanced SI-thalamic functional connectivity as a core feature of the pathophysiology underlying MDD. Furthermore, this abnormality is related to anhedonia, memory, and attention. The SI-thalamic functional connectivity may be a useful neural target for affective experience and cognitive interventions. Future research designed to track the effects of current treatments may develop more targeted pathophysiological therapies.



Continuous performance test


Diagnostic and Statistical Manual of Mental Disorders - Fourth Edition


17-item Hamilton depression rating scale


Healthy controls


Major depressive disorder


Montreal Neurological Institute


Assessment of Neuropsychological Status


Regions of interest


Resting-state functional connectivity


Resting-state functional magnetic resonance imaging


The Chinese version of the Modified Structured Clinical Interview for DSM-IV TR Axis I Disorders Patient Edition


Standard deviation


Snaith-Hamilton Pleasure Scale


Primary somatosensory cortex

SPSS 23.0:

IBM SPSS Statistics Version 23.0


  1. Kessler RC, Berglund P, Demler O, et al. Lifetime prevalence and age-of-onset distributions of DSM-IV disorders in the National Comorbidity Survey Replication. Arch Gen Psychiatry. 2005;62(6):593–602.

    Article  Google Scholar 

  2. Bora E, Harrison BJ, Yucel M, Pantelis C. Cognitive impairment in euthymic major depressive disorder: a meta-analysis. Psychol Med. 2013;43(10):2017–26.

    Article  CAS  PubMed  Google Scholar 

  3. Snyder HR. Major depressive disorder is associated with broad impairments on neuropsychological measures of executive function: a meta-analysis and review. Psychol Bull. 2013;139(1):81–132.

    Article  Google Scholar 

  4. Mak ADP, Lau DTY, Chan AKW, So SH, Leung O, Wong SL, et al. Cognitive impairment in treatment-Naïve bipolar II and unipolar depression. Sci Rep. 2018;8(1):1905.

    Article  Google Scholar 

  5. Polosa M, Lemogne C, Jardri R, Fossati P. Cognition-the core of major depressive disorder. L’Encephale. 2016;42(1 Suppl 1):1S3–11.

    Article  Google Scholar 

  6. Ahern E, Semkovska M. Cognitive functioning in the first-episode of major depressive disorder: a systematic review and meta-analysis. Neuropsychology. 2017;31(1):52–72.

    Article  PubMed  Google Scholar 

  7. Nierenberg AA, Husain MM, Trivedi MH, Fava M, Warden D, Wisniewski SR, et al. Residual symptoms after remission of major depressive disorder with citalopram and risk of relapse: a STAR* D report. Psychol Med. 2010;40(1):41–50.

    Article  CAS  Google Scholar 

  8. Trivedi MH, Greer TL. Cognitive dysfunction in unipolar depression: implications for treatment. J Affect Disord. 2014;152-154:19–27.

    Article  PubMed  Google Scholar 

  9. Levin HS, Eisenberg HM, Benton AL. Frontal lobe function and dysfunction. Frontal lobe function and dysfunction. Oxford: Oxford University Press; 1991.

  10. Modha DS, Mishkin M. Network architecture of the long-distance pathways in the macaque brain. Proc Natl Acad Sci U S A. 2010;107(30):13485–90.

    Article  CAS  Google Scholar 

  11. Peters SK, Dunlop K, Downar J. Cortico-striatal-thalamic loop circuits of the salience network: a central pathway in psychiatric disease and treatment. Front Syst Neurosci. 2016;10:104.

    Article  Google Scholar 

  12. Chang YT, Lu CH, Wu MK, Hsu SW, Huang CW, Chang WN, et al. Salience network and depressive severities in Parkinson’s disease with mild cognitive impairment: a structural covariance network analysis. Front Aging Neurosci. 2018;9:417.

    Article  Google Scholar 

  13. Parnaudeau S, Bolkan SS, Kellendonk C. The Mediodorsal thalamus: an essential partner of the prefrontal cortex for cognition. Biol Psychiatry. 2018;83(8):648–56.

    Article  PubMed  Google Scholar 

  14. Schoonheim MM, Hulst HE, Brandt RB, Strik M, Wink AM, Uitdehaag BM, et al. Thalamus structure and function determine severity of cognitive impairment in multiple sclerosis. Neurology. 2015;84(8):776–83.

    Article  Google Scholar 

  15. Tewarie P, Schoonheim MM, Schouten DI, Polman CH, Balk LJ, Uitdehaag BM, et al. Functional brain networks: linking thalamic atrophy to clinical disability in multiple sclerosis, a multimodal fMRI and MEG study. Hum Brain Mapp. 2015;36(2):603–18.

    Article  Google Scholar 

  16. Bisecco A, Rocca MA, Pagani E, Mancini L, Enzinger C, Gallo A, et al. Connectivity-based parcellation of the thalamus in multiple sclerosis and its implications for cognitive impairment: a multicenter study. Hum Brain Mapp. 2015;36(7):2809–25.

    Article  Google Scholar 

  17. Houtchens MK, Benedict RHB, Killiany R, Sharma J, Jaisani Z, Singh B, et al. Thalamic atrophy and cognition in multiple sclerosis. Neurology. 2007;69(12):1213–23.

    Article  CAS  Google Scholar 

  18. Koenig KA, Rao SM, Lowe MJ, Lin J, Sakaie KE, Stone L, et al. The role of the thalamus and hippocampus in episodic memory performance in patients with multiple sclerosis. Mult Scler J. 2018;1:1–11.

  19. Zhang D, Snyder AZ, Shimony JS, Fox MD, Raichle ME, et al. Noninvasive functional and structural connectivity mapping of the human Thalamocortical system. Cereb Cortex. 2010;20(5):1187–94.

    Article  Google Scholar 

  20. Zhang D, Snyder AZ, Fox MD, Sansbury MW, Shimony JS, Raichle ME, et al. Intrinsic functional relations between human cerebral cortex and thalamus. J Neurophysiol. 2008;100(4):1740–8.

    Article  Google Scholar 

  21. Penfield W, Boldrey E. Somatic motor and sensory representation in the cerebral cortex of man as studied by electrical stimulation. Brain. 1937;60(4):389–443.

    Article  Google Scholar 

  22. Simons D. Neuronal integration in the somatosensory whisker/barrel cortex. Cereb Cortex. 1995;11:263–97.

    Article  Google Scholar 

  23. Galvez R, Weiss C, Weible AP, Disterhoft JF, et al. Vibrissa-signaled eyeblink conditioning induces somatosensory cortical plasticity. J Neurosci the Official Journal of the Society for Neuroscience. 2006;26(22):6062–8.

    Article  CAS  Google Scholar 

  24. Chau LS, Davis AS, Galvez R. Neocortical synaptic proliferation following forebrain-dependent trace associative learning. Behav Neurosci. 2013;127(2):285–92.

    Article  Google Scholar 

  25. Tseng MT, Kong Y, Eippert F, Tracey I, et al. Determining the neural substrate for encoding a memory of human pain and the influence of anxiety. J Neurosci. 2017;37(49):11806–17.

    Article  CAS  Google Scholar 

  26. Zhao D, Zhou YD, Bodner M, Ku Y, et al. The causal role of the prefrontal cortex and somatosensory cortex in tactile working memory. Cereb Cortex. 2017;28(10):1–10.

    Article  Google Scholar 

  27. Hegner YL, Lutzenberger W, Leiberg S, Braun C, et al. The involvement of ipsilateral temporoparietal cortex in tactile pattern working memory as reflected in beta event-related desynchronization. NeuroImage. 2007;37(4):1362–70.

    Article  Google Scholar 

  28. Khoshnejad M, Roy M, Martinu K, Chen JI, Cohen-Adad J, Grondin S, et al. Brain processing of the temporal dimension of acute pain in short-term memory. Pain. 2017;158(10):2001–11.

    Article  Google Scholar 

  29. Schaefer M, Northoff G. Who am I: the conscious and the unconscious self. Front Hum Neurosci. 2017;11:126.

    Article  Google Scholar 

  30. Biswal B, Zerrin Yetkin F, Haughton VM, Hyde JS. Functional connectivity in the motor cortex of resting human brain using echo-planarmri. Magn Reson Med. 1995;34(4):537–41.

    Article  CAS  Google Scholar 

  31. Zuo XN, Kelly C, Martino AD, Mennes M, Margulies DS, Bangaru S, et al. Growing together and growing apart: regional and sex differences in the lifespan developmental trajectories of functional homotopy. J Neurosci the Official Journal of the Society for Neuroscience. 2010;30(45):15034–43.

    Article  CAS  Google Scholar 

  32. Gorges M, Müller HP, Lulé D, Ludolph AC, Pinkhardt EH, Kassubek J, et al. Functional connectivity within the default mode network is associated with saccadic accuracy in Parkinson's disease: a resting-state FMRI and videooculographic study. Brain Connectivity. 2013;3(3):265–72.

    Article  PubMed  Google Scholar 

  33. Harshaw C. Interoceptive dysfunction: toward an integrated framework for understanding somatic and affective disturbance in depression. Psychol Bull. 2015;141(2):311.

    Article  Google Scholar 

  34. Hamilton M. A rating scale for depression. J Neurol Neurosurg Psychiatry. 1960;23(1):56.

    Article  CAS  Google Scholar 

  35. Snaith RP, Hamilton M, Morley S, Humayan A, Hargreaves D, Trigwell P. A scale for the assessment of hedonic tone the snaith-Hamilton pleasure scale. Br J Psychiatry. 1995;167(1):99–103.

    Article  CAS  Google Scholar 

  36. Krupp LB, LaRocca NG, Muir-Nash J, Steinberg AD, et al. The fatigue severity scale: application to patients with multiple sclerosis and systemic lupus erythematosus. Arch Neurol. 1989;46(10):1121–3.

    Article  CAS  Google Scholar 

  37. Randolph C, Tierney MC, Mohr E, Chase TN, et al. The repeatable battery for the assessment of neuropsychological status (RBANS): preliminary clinical validity. J Clin Exp Neuropsychol. 1998;20(3):310–9.

    Article  CAS  Google Scholar 

  38. Mur M, Portella MJ, Martínezarán A, Pifarre J, Vieta E, et al. Long-term stability of cognitive impairment in bipolar disorder: a 2-year follow-up study of lithium-treated euthymic bipolar patients. J Clin Psychiatry. 2008;69(5):712–9.

    Article  CAS  Google Scholar 

  39. Yan C, Zang Y. DPARSF: a MATLAB toolbox for “pipeline” data analysis of resting-state fMRI. Front Syst Neurosci. 2010;4:13.

    Article  Google Scholar 

  40. Data Processing Assistant for Resting-State fMRI (DPARSF). Accessed 20 May 2017.

  41. Maldjian JA, Laurienti PJ, Kraft A, Burdette JH. An automated method for neuroanatomic and cytoarchitectonic atlas-based interrogation of fMRI data sets. NeuroImage. 2003;19(3):1233–39.

    Article  Google Scholar 

  42. Software Packages that are currently available for download from the ANSIR Laboratory at Wake Forest University School of Medicine. Accessed 22 May 2017.

  43. Salvador R, Suckling J, Schwarzbauer C, Bullmore E, et al. Undirected graphs of frequency-dependent functional connectivity in whole brain networks. Philos Trans R Soc Lond Ser B Biol Sci. 2005;360(1457):937–46.

    Article  Google Scholar 

  44. Liu Y, Yu C, Liang M, Li J, Tian L, Zhou Y, et al. Whole brain functional connectivity in the early blind. Brain. 2007;130(8):2085–96.

    Article  Google Scholar 

  45. Ledberg A, Åkerman S, Roland PE. Estimation of the probabilities of 3D clusters in functional brain images. NeuroImage. 1998;8(2):113–28.

    Article  CAS  Google Scholar 

  46. Simultaneous Inference for FMRI Data. Accessed 20 May 2017.

  47. Kong QM, Qiao H, Liu CZ, Zhang P, Li K, Wang L, et al. Aberrant intrinsic functional connectivity in thalamo-cortical networks in major depressive disorder. CNS Neurosci Ther. 2018;28:1–10.

    Article  Google Scholar 

  48. Kim K, Kim SW, Myung W, Han CE, Fava M, Mischoulon D, et al. Reduced orbitofrontal-thalamic functional connectivity related to suicidal ideation in patients with major depressive disorder. Sci Rep. 2017;7(1):15772;1–11.

  49. Zhu X, Helpman L, Papini S, Schneier F, Markowitz JC, Meter PE, et al. Altered resting state functional connectivity of fear and reward circuitry in comorbid PTSD and major depression. Depression and anxiety. 2017;34(7):641–50.

    Article  Google Scholar 

  50. Brown EC, Clark DL, Hassel S, MacQueen G, Ramasubbu R, et al. Thalamocortical connectivity in major depressive disorder. J Affect Disord. 2017;217:125–31.

    Article  Google Scholar 

  51. Woodward ND, Karbasforoushan H, Heckers S. Thalamocortical dysconnectivity in schizophrenia. Am J Psychiatry. 2012;169(10):1092–9.

    Article  Google Scholar 

  52. Krukow P, Jonak K, Karakuła-Juchnowicz H, Podkowiński A, Jonak K, Borys M, et al. Disturbed functional connectivity within the left prefrontal cortex and sensorimotor areas predicts impaired cognitive speed in patients with first-episode schizophrenia. Psychiatry Res Neuroimaging. 2018;275:28–35.

    Article  PubMed  Google Scholar 

  53. Löken LS, Wessberg J, McGlone F, Morrison I, Olausson H, et al. Coding of pleasant touch by unmyelinated afferents in humans. Nat Neurosci. 2009;12(5):547.

    Article  Google Scholar 

  54. Sel A, Calvo-Merino B, Tuettenberg S, Forster B, et al. When you smile, the world smiles at you: ERP evidence for self-expression effects on face processing. Soc Cogn Affect Neurosci. 2015;10(10):1316–22.

    Article  Google Scholar 

  55. Gazzola V, Spezio ML, Etzel JA, Castelli F, Adolphs R, Keysers C, et al. Primary somatosensory cortex discriminates affective significance in social touch. Proc Natl Acad Sci. 2012;109(25):E1657–66.

    Article  CAS  Google Scholar 

  56. Case LK, Laubacher CM, Olausson H, Wang B, Spagnolo PA, Bushnell MC, et al. Encoding of touch intensity but not pleasantness in human primary somatosensory cortex. J Neurosci. 2016;36(21):5850–60.

    Article  CAS  Google Scholar 

  57. McCabe C, Rolls ET, Bilderbeck A, Mcglone F, et al. Cognitive influences on the affective representation of touch and the sight of touch in the human brain. Soc Cogn Affect Neurosci. 2008;3(2):97–108.

    Article  Google Scholar 

  58. Hara Y, Ago Y, Higuchi M, Hasebe S, Nakazawa T, Hashimoto H, et al. Oxytocin attenuates deficits in social interaction but not recognition memory in a prenatal valproic acid-induced mouse model of autism. Horm Behav. 2017;96:130–6.

    Article  CAS  Google Scholar 

  59. Boutin A, Pinsard B, Boré A, Carrier FSM, Doyon J, et al. Transient synchronization of hippocampo-striato-thalamo-cortical networks during sleep spindle oscillations induces motor memory consolidation. NeuroImage. 2018;169:419–30.

    Article  PubMed  Google Scholar 

  60. Gogulski J, Zetter R, Nyrhinen M, Pertovaara A, Carlson S, et al. Neural substrate for metacognitive accuracy of tactile working memory. Cereb Cortex. 2017;27(11):5343–52.

    Article  Google Scholar 

  61. Giraldo-Chica M, Rogers BP, Damon SM, Landman BA, Woodward ND, et al. Prefrontal-thalamic anatomical connectivity and executive cognitive function in schizophrenia. Biol Psychiatry. 2018;83(6):509–17.

    Article  Google Scholar 

  62. Brecht M. The body model theory of somatosensory cortex. Neuron. 2017;94(5):985–92.

    Article  CAS  Google Scholar 

  63. Petersen SE, Robinson DL, Keys W. Pulvinar nuclei of the behaving rhesus monkey: visual responses and their modulation. J Neurophysiol. 1985;54(4):867–86.

    Article  CAS  Google Scholar 

  64. Nakajima M, Halassa MM. Thalamic control of functional cortical connectivity. Curr Opin Neurobiol. 2017;44:127–31.

    Article  CAS  Google Scholar 

  65. Schmitt LI, Halassa MM. Interrogating the mouse thalamus to correct human neurodevelopmental disorders. Mol Psychiatry. 2017;22(2):183.

    Article  CAS  Google Scholar 

  66. Duncan J. The structure of cognition: attentional episodes in mind and brain. Neuron. 2013;80(1):35–50.

    Article  CAS  Google Scholar 

  67. Schmitt LI, Wimmer RD, Nakajima M, Happ M, Mofakham S, Halassa MM, et al. Thalamic amplification of cortical connectivity sustains attentional control. Nature. 2017;545(7653):219.

    Article  CAS  Google Scholar 

Download references


We would like to thank all the participants in this study.


Design and data collection of this work were supported by grants from the National Natural Science Foundation of China (81471379 and 81601193) and the National Clinical Research Center on Mental Disorders (2015BAI13B02). Data analysis and interpretation were supported by the National Natural Science Foundation of China (81701345) and the Shanxi Provincial Science Foundation for Youth (2015021204). Writing the manuscript was supported by the National Key R&D Program of China (2016YFC1307103). We sincerely thank the patients and the healthy volunteers for their participation and all the medical staff involved in the collection of cases.

Availability of data and materials

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

Author information

Authors and Affiliations



ZKR designed the study and wrote the protocol. SN, WYF and ZAX were responsible for all aspects of patient recruitment and clinical diagnostic assessment. YCX, KLJ, LZF contributed to the collection of the data. ZAX, LPH, LGZ and KLJ performed the statistical analysis. KLJ and ZAX prepared the first draft. ZKR contributed to the conceptualization of the study and reviewed several manuscript drafts. All authors contributed to and approved the final manuscript.

Corresponding author

Correspondence to Kerang Zhang.

Ethics declarations

Ethics approval and consent to participate

The Ethical Committee for Medicine of the First Hospital of Shanxi Medical University approved this study. Written informed consent was received from all participants prior to inclusion.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

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 distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Kang, L., Zhang, A., Sun, N. et al. Functional connectivity between the thalamus and the primary somatosensory cortex in major depressive disorder: a resting-state fMRI study. BMC Psychiatry 18, 339 (2018).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: