Major depressive disorder (MDD) is a heterogeneous clinical syndrome, involving affective, cognitive, and somatic symptoms. Neuronal underpinnings and mechanisms of MDD still remain poorly understood. Exploring the electrophysiological basis of MDD is important to delineate the heterogeneity of MDD. Neuronal oscillations are fundamental for brain function allowing dynamic routing of activity in densely interconnected brain networks. We investigated whether large-scale oscillatory neuronal connectivity could reflect the heterogeneous MDD symptoms. Resting-state magnetoencephalography (MEG) and magnetic resonance imaging (MRI) were collected from 171 MDD patients to identify aberrant neuronal oscillation dynamics in patients. MEG data were preprocessed, source-reconstructed in individual cortical anatomy, and expressed as cortical parcel time series collapsed to 200 parcels of the Schaefer atlas. Inter-areal phase synchronization (PS) was estimated with a weighted phase lag index. We then estimated the correlation of oscillatory PS networks with individual symptom data.Preliminary results indicate that MDD patients exhibited significant correlations with symptoms in a frequency- and symptom-specific manner. More specifically we found PS networks in delta, high alpha, and beta frequencies to exhibit symptom correlations. The PS networks in different frequencies were also localized to anatomically different circuits.These results indicate that oscillatory PS networks reflect the heterogeneity of MDD symptom profiles.