BackgroundRecent studies have shown that major depressive disorder (MDD) is associated with altered intrinsic functional connectivity (FC) of the thalamus; however, investigations of these alterations at a finer time scale and the level of thalamic subregions are still lacking.MethodsWe collected resting-state functional MRI data from 100 treatment-naïve, first-episode MDD patients and 99 age-, gender- and education-matched healthy controls (HCs). Seed-based whole-brain sliding window-based dFC analyses were performed for 16 thalamic subregions. Between-group differences in the mean and variance of dFC were determined using threshold-free cluster enhancement algorithm. For significant alterations, there relationships with clinical and neuropsychological variables were further examined via bivariate and multivariate correlation analyses.ResultsOf all thalamic subregions, only the left sensory thalamus (Stha) showed altered variance of dFC in the patients characterized by increases with the left inferior parietal lobule, left superior frontal gyrus, left inferior temporal gyrus, and left precuneus, and decreases with multiple frontal, temporal, parietal, and subcortical regions. These alterations accounted for, to a great extent, clinical, and neuropsychological characteristics of the patients as revealed by the multivariate correlation analysis. In addition, the bivariate correlation analysis revealed a positive correlation between the variance of dFC between the left Stha and right inferior temporal gurus/fusiform and childhood trauma questionnaires scores (r= 0.562,P< 0.001).ConclusionThese findings suggest that the left Stha is the most vulnerable thalamic subregion to MDD, whose dFC alterations may serve as potential biomarkers for the diagnosis of the disease.