Background: Non-alcoholic fatty liver disease (NAFLD) stands as a predominant chronic liver ailment globally, yet its pathogenesis remains elusive. This study aims to identify Hub mitophagyrelated genes (MRGs), and explore the underlying pathological mechanisms through which these hub genes regulate NAFLD. Methods: A total of 3 datasets were acquired from the GEO database and integrated to identify differentially expressed genes (DEGs) in NAFLD and perform Gene Set Enrichment Analysis (GSEA). By intersecting DEGs with MRGs, mitophagy-related differentially expressed genes (MRDEGs) were obtained. Then, hub MRGs with diagnostic biomarker capability for NAFLD were screened and a diagnostic prediction model was constructed and assessed using Nomogram, Decision Curve Analysis (DCA), and ROC curves. Functional enrichment analysis was conducted on the identified hub genes to explore their biological significance. Additionally, regulatory networks were constructed using databases. NAFLD was stratified into high and low -risk groups based on the Riskscore from the diagnostic prediction model. Furthermore, single -sample gene set enrichment analysis (ssGSEA) and CIBERSORT algorithms were employed to analyze immune cell infiltration patterns and the relationship between Hub MRGs and immune cells. Results: The integrated dataset comprised 122 NAFLD samples and 31 control samples. After screening, 18 MRDEGs were identified. Subsequently, six hub MRGs (NR4A1, PPP2R2A, P4HA1, TUBB6, DUSP1, NAMPT) with diagnostic potential were selected through WGCNA, logistic regression, SVM, RF, and LASSO models, all significantly downregulated in NAFLD samples compared to the control group. A diagnostic prediction model based on these six genes demonstrated robust predictive performance. Functional enrichment analysis of the six hub genes revealed involvement in processes such as protein phosphorylation or dephosphorylation. Correlation analysis demonstrated a significant association between hub MRGs and infiltrating immune cells. Conclusion: We identified six hub MRGs in NAFLD and constructed a diagnostic prediction model based on these six genes, applicable for early NAFLD diagnosis. These genes may participate in regulating NAFLD progression through the modulation of mitophagy and immune activation. Our findings may contribute to subsequent clinical and basic research on NAFLD.