Non alcoholic fatty liver disease (NAFLD) is the most common metabolic disease in humans, affecting the majority of individuals. In the current investigation, we aim to identify potential key genes linked with NAFLD through bioinformatics analyses of next generation sequencing (NGS) dataset. NGS dataset of GSE135251 from the Gene Expression Omnibus (GEO) database were retrieved. Differentially expressed genes (DEGs) were obtained by DESeq2 package. g:Profiler database was further used to identify the potential gene ontology (GO) and REACTOME pathways. Protein-protein interaction (PPI) network was constructed using the Hippie interactome database. miRNet and NetworkAnalyst databases were used to establish a miRNA-hub gene regulatory network and TF-hub gene regulatory network for the hub genes. Hub genes were verified based on receiver operating characteristic curve (ROC) analysis. Totally, 951 DEGs were identified including 476 up regulated genes and 475 down regulated genes screened in NAFLD and normal control. GO showed that DEGs were significantly enhanced for signaling and regulation of biological quality. REACTOME pathway analysis revealed that DEGs were enriched in signaling by interleukins and extracellular matrix organization. ESR2, JUN, PTN, PTGER3, CEBPB, IKBKG, HSPA8, SFN, CDKN1A and E2F1 were indicated as hub genes from PPI network, miRNA-hub gene regulatory network and TF-hub gene regulatory network. Furthermore, ROC analysis revealed that ESR2, JUN, PTN, PTGER3, CEBPB, IKBKG, HSPA8, SFN, CDKN1A and E2F1 might serve as diagnostic biomarkers in NAFLD. The current investigation provided insights into the molecular mechanism of NAFLD that might be useful in further investigations.