Identification of Autophagy-Related Genes in the Progression from Non-Alcoholic Fatty Liver to Non-Alcoholic Steatohepatitis

BACKGROUND: Non-alcoholic fatty liver disease (NAFLD) is the most common cause of chronic liver disease worldwide. Autophagy plays a vital role in NAFLD development and progression. We aimed to establish a novel autophagy-related gene (ARG) signature as a therapeutic target in NAFLD patients based on high-throughput sequencing data. METHODS: ARGs obtained from the HAMdb and high-sequencing data obtained from the Gene Expression Omnibus (GEO) database were analyzed to identify differentially expressed ARGs (DEARGs) between normal and NASH tissues. Then, gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were performed to explore potential biological and pathological functions of DEARGs. The protein–protein interaction (PPI) network of the DEARGs was established through the STRING website, and visualized by Cytoscape. In addition, hub genes were validated by an independent dataset GSE89632. Finally, we performed Gene Set Variation Analysis (GSVA) pathway-related analysis to identify the pivotal signaling pathways and genes for the progression of non-alcoholic fatty liver (NAFL) to non-alcoholic steatohepatitis (NASH). RESULTS: A total of 76 DEARGs were identified in the GSE126848 dataset, of which 45 genes were upregulated and 31 genes were downregulated. GO analysis showed that the biological functions of DEARGs focused primarily on autophagy, cellular response to external stimulus, fibroblast proliferation, late endosome, and ubiquitin protein ligase binding. KEGG pathway analysis showed that these DEARGs were mainly involved in the apoptosis, PI3K-Akt signaling pathway, and estrogen signaling pathway. Among DEGs, 9 most closely related genes were identified from the PPI network. Furthermore, NOS3, IGF1, VAMP8, FOS, and HMOX1 were verified in the GSE89632 dataset. At last, the MAPK signal pathway was identified as important pathway, and JUN was identified as a key gene involved in the progression from NAFL to NASH. CONCLUSION: This study may provide credible molecular biomarkers in terms of screening and diagnosis for NAFLD. Meanwhile, it also serves as a basis for exploring the molecular mechanisms underlying the progression of NAFL to NASH.