Background: Metabolic-associated fatty liver disease (MAFLD) is a major cause of liver disease worldwide and comprises non-alcoholic fatty liver (NAFL) and non-alcoholic steatohepatitis (NASH). Due to the high prevalence and poor prognosis of NASH, it is critical to identify and treat patients at risk. However, the aetiology and mechanisms remain largely unknown, warranting further analysis.
Methods: We first identified differential genes in NASH by single-cell analysis of the GSE129516 dataset and conducted expression profiling data analysis of the GSE184019 dataset from the Gene Expression Omnibus (GEO) database. Then single-cell trajectory reconstruction and analysis, immune gene score, cellular communication, key gene screening, functional enrichment analysis, and immune microenvironment analysis were carried out. Finally, cell experiments were performed to verify the role of key genes in NASH.
Results: We conducted transcriptome profiling of 30,038 single cells, including hepatocytes and non-hepatocytes from normal and steatosis adult mouse livers. Comparative analysis of hepatocytes and non-hepatocytes revealed pronounced heterogeneity as non-hepatocytes acted as major cell-communication hubs. The results showed that Hspa1b, Tfrc, Hmox1 and Map4k4 could effectively distinguish NASH tissues from normal samples. The results of scRNA-seq and qPCR indicated that the expression levels of hub genes in NASH were significantly higher than in normal cells or tissues. Further immune infiltration analysis showed significant differences in M2 macrophage distribution between healthy and metabolic-associated fatty liver samples.
Conclusions: Our results suggest that Hspa1b, Tfrc, Hmox1 and Map4k4 have huge prospects as diagnostic and prognostic biomarkers for NASH and may be potential therapeutic targets for NASH.
Keywords: GEO database; Metabolic-associated fatty liver disease; cell experiments; single-cell sequencing analysis.