Background Neuroimaging studies on major depressive disorder (MDD) have identified an extensive range of brain structural abnormalities, but the exact neural mechanisms associated with MDD remain elusive. Most previous studies were performed with voxel- or surface-based morphometry which were univariate methods without considering spatial information across voxels/vertices. Methods Brain morphology was investigated using voxel-based morphometry (VBM) and source-based morphometry (SBM) in 1082 MDD patients and 990 healthy controls (HCs) from the REST-meta-MDD Consortium. We first examined group differences in regional grey matter (GM) volumes and structural covariance networks between patients and HCs. We then compared first-episode, drug-na...