BACKGROUND: Despite efforts to improve targeting accuracy of the dorsolateral prefrontal cortex (DLPFC) as a repetitive transcranial magnetic stimulation (rTMS) target for Major Depressive Disorder (MDD), the heterogeneity in clinical response remains unexplained. OBJECTIVE: We sought to compare the patterns of functional connectivity from the DLPFC treatment site in patients with MDD who were TMS responders to those who were TMS non-responders. METHODS: Baseline anatomical T1 magnetic resonance imaging (MRI), resting-state functional MRI, and diffusion weighted imaging scans were obtained from 37 participants before they underwent a course of rTMS to left Brodmann area 46. A novel machine learning method was utilized to identify brain regions associated with each item of the Beck's Depression Inventory II (BDI-II), and for 26 participants who underwent rTMS treatment over the left Brodmann area 46, identify regions differentiating rTMS responders and non-responders. RESULTS: Nine parcels of the Human Connectome Project Multimodal Parcellation Atlas matched to at least three items of the Beck's Depression Inventory II (BDI-II) as predictors of response to rTMS, with many in the temporal, parietal and cingulate cortices. Additionally, pre-treatment mapping for 17 items of the BDI-II demonstrated significant variability in symptom to parcel mapping. When parcels associated with symptom presence and symptom resolution were compared, 15 parcels were uniquely associated with resolution (potential targets), and 12 parcels were associated with both symptom presence and resolution (blockers or biomarkers). CONCLUSIONS: Machine learning approaches show promise for the development of pathoanatomical diagnosis and treatment algorithms for MDD. Prospective studies are required to facilitate clinical translation.