BackgroundThere are emerging studies suggesting that non-alcoholic fatty liver disease (NAFLD) is a heterogeneous disease with multiple etiologies and molecular phenotypes. Fibrosis is the key process in NAFLD progression. In this study, we aimed to explore molecular phenotypes of NAFLD with a particular focus on the fibrosis phenotype and also aimed to explore the changes of macrophage subsets in the fibrosis subset of NAFLD.MethodsTo assess the transcriptomic alterations of key factors in NAFLD and fibrosis progression, we included 14 different transcriptomic datasets of liver tissues. In addition, two single-cell RNA sequencing (scRNA-seq) datasets were included to construct transcriptomic signatures that could represent specific cells. To explore the molecular subsets of fibrosis in NAFLD based on the transcriptomic features, we used a high-quality RNA-sequencing (RNA-seq) dataset of liver tissues from patients with NAFLD. Non-negative matrix factorization (NMF) was used to analyze the molecular subsets of NAFLD based on the gene set variation analysis (GSVA) enrichment scores of key molecule features in liver tissues.ResultsThe key transcriptomic signatures on NAFLD including non-alcoholic steatohepatitis (NASH) signature, fibrosis signature, non-alcoholic fatty liver (NAFL) signature, liver aging signature and TGF-β signature were constructed by liver transcriptome datasets. We analyzed two liver scRNA-seq datasets and constructed cell type-specific transcriptomic signatures based on the genes that were highly expressed in each cell subset. We analyzed the molecular subsets of NAFLD by NMF and categorized four main subsets of NAFLD. Cluster 4 subset is mainly characterized by liver fibrosis. Patients with Cluster 4 subset have more advanced liver fibrosis than patients with other subsets, or may have a high risk of liver fibrosis progression. Furthermore, we identified two key monocyte-macrophage subsets which were both significantly correlated with the progression of liver fibrosis in NAFLD patients.ConclusionOur study revealed the molecular subtypes of NAFLD by integrating key information from transcriptomic expression profiling and liver microenvironment, and identified a novel and distinct fibrosis subset of NAFLD. The fibrosis subset is significantly correlated with the profibrotic macrophages and M2 macrophage subset. These two liver macrophage subsets may be important players in the progression of liver fibrosis of NAFLD patients.