Alterations in eating behavior characterized eating disorders (ED). The genetic factors shared between ED diagnoses have been underexplored. The present study performed a genome-wide association study in individuals with disordered eating behaviors in the Mexican population, blood methylation quantitative trait loci (blood-meQTL), summary data-based Mendelian randomization (SMR) analysis, and in silico function prediction by different algorithms. The analysis included a total of 1803 individuals. We performed a genome-wide association study and blood-meQTL analysis by logistic and linear regression. In addition, we analyzed in silico functional variant prediction, phenome-wide, and multi-tissue expression quantitative trait loci. The genome-wide association study identified 44 single-nucleotide polymorphisms (SNP) associated at a nominal value and seven blood-meQTL at a genome-wide threshold. The SNPs show enrichment in genome-wide associations of the metabolic and immunologic domains. In the in silico analysis, the SNP rs10419198 (-value = 4.85 × 10) located on an enhancer mark could change the expression of in blood, adipocytes, and brain areas that regulate food intake. Additionally, we found an association of DNA methylation levels of (-value = 6.76 × 10) and (-value = 5.73 × 10) by SMR analysis. The present study supports the previous associations of genetic variation in the metabolic domain with ED.