Identification of m6A-associated autophagy genes in non-alcoholic fatty liver

Background Studies had shown that autophagy was closely related to nonalcoholic fat liver disease (NAFLD), while N6-methyladenosine (m6A) was involved in the regulation of autophagy. However, the mechanism of m6A related autophagy in NAFLD was unclear. Methods The NAFLD related datasets were gained via the Gene Expression Omnibus (GEO) database, and we also extracted 232 autophagy-related genes (ARGs) and 37 m6A. First, differentially expressed ARGs (DE-ARGs) and differentially expressed m6A (DE-m6A) were screened out by differential expression analysis. DE-ARGs associated with m6A were sifted out by Pearson correlation analysis, and the m6A-ARGs relationship pairs were acquired. Then, autophagic genes in m6A-ARGs pairs were analyzed for machine learning algorithms to obtain feature genes. Further, we validated the relationship between feature genes and NAFLD through quantitative real-time polymerase chain reaction (qRT-PCR), Western blot (WB). Finally, the immuno-infiltration analysis was implement, and we also constructed the TF-mRNA and drug-gene networks. Results There were 19 DE-ARGs and four DE-m6A between NAFLD and normal samples. The three m6A genes and five AGRs formed the m6A-ARGs relationship pairs. Afterwards, genes obtained from machine learning algorithms were intersected to yield three feature genes (TBK1, RAB1A, and GOPC), which showed significant positive correlation with astrocytes, macrophages, smooth muscle, and showed significant negative correlation with epithelial cells, and endothelial cells. Besides, qRT-PCR and WB indicate that TBK1, RAB1A and GOPC significantly upregulated in NAFLD. Ultimately, we found that the TF-mRNA network included FOXP1-GOPC, ATF1-RAB1A and other relationship pairs, and eight therapeutic agents such as R-406 and adavosertib were predicted based on the TBK1. Conclusion The study investigated the potential molecular mechanisms of m6A related autophagy feature genes (TBK1, RAB1A, and GOPC) in NAFLD through bioinformatic analyses and animal model validation. However, it is critical to note that these findings, although consequential, demonstrate correlations rather than cause-and-effect relationships. As such, more research is required to fully elucidate the underlying mechanisms and validate the clinical relevance of these feature genes.