Recent studies have provided promising evidence that neuroimaging data can predict treatment outcomes for patients with major depressive disorder (MDD). As most of these studies had small sample sizes, a meta-analysis is warranted to identify the most robust findings and imaging modalities, and to compare predictive outcomes obtained in magnetic resonance imaging (MRI) and studies using clinical and demographic features. We conducted a literature search from database inception to July 22, 2023, to identify studies using pretreatment clinical or brain MRI features to predict treatment outcomes in patients with MDD. Two meta-analyses were conducted on clinical and MRI studies, respectively. The meta-regression was employed to explore the effects of covariates and compare the predictive performance between clinical and MRI groups, as well as across MRI modalities and intervention subgroups. Meta-analysis of 13 clinical studies yielded an area under the curve (AUC) of 0.73, while in 44 MRI studies, the AUC was 0.89. MRI studies showed a higher sensitivity than clinical studies (0.78 vs. 0.62, Z = 3.42, P = 0.001). In MRI studies, resting-state functional MRI (rsfMRI) exhibited a higher specificity than task-based fMRI (tbfMRI) (0.79 vs. 0.69, Z = -2.86, P = 0.004). No significant differences in predictive performance were found between structural and functional MRI, nor between different interventions. Of note, predictive MRI features for treatment outcomes in studies using antidepressants were predominantly located in the limbic and default mode networks, while studies of electroconvulsive therapy (ECT) were restricted mainly to the limbic network. Our findings suggest a promise for pretreatment brain MRI features to predict MDD treatment outcomes, outperforming clinical features. While tasks in tbfMRI studies differed, those studies overall had less predictive utility than rsfMRI data. Overlapping but distinct network-level measures predicted antidepressants and ECT outcomes. Future studies are needed to predict outcomes using multiple MRI features, and to clarify whether imaging features predict outcomes generally or differ depending on treatments.