Migraine is a type of headache with multiple neurological symptoms. Prior neuroimaging studies in patients with migraine based on functional magnetic resonance imaging have found regional as well as network-level alterations in brain function. Here, we expand on prior studies by establishing whole-brain functional connectivity patterns in patients with migraine using dimensionality reduction techniques. We studied functional brain connectivity in 50 patients with episodic migraine and sex- and age-matched healthy controls. Using dimensionality reduction techniques that project high-dimensional functional connectivity onto low-dimensional representations (i.e., eigenvectors), we found significant between-group differences in the eigenvectors between patients with migraine and healthy controls, particularly in the sensory/motor and limbic cortices. Furthermore, we assessed between-group differences in subcortical connectivity with subcortical weighted manifolds defined by subcortico-cortical connectivity multiplied by cortical eigenvectors and revealed significant alterations in the amygdala. Finally, leveraging supervised machine learning, we moderately predicted headache frequency using cortical and subcortical functional connectivity features, again indicating that sensory and limbic regions play a particularly important role in predicting migraine frequency. Our study confirmed that migraine is a hierarchical disease of the brain that shows alterations along the sensory-limbic axis, and therefore, the functional connectivity in these areas could be a useful marker to investigate migraine symptomatology.