# Hyperlipidemia May Synergize with Hypomethylation in Establishing Trained Immunity and Promoting Inflammation in NASH and NAFLD

Authors: Charles I. V. Drummer; Fatma Saaoud; Yu Sun; Diana Atar; Keman Xu; Yifan Lu; Ying Shao; Candice Johnson; Lu Liu; Huimin Shen; Nirag C. Jhala; Xiaohua Jiang; Hong Wang; Xiaofeng Yang
Journal: Journal of Immunology Research
 (2021)

---

## Abstract

We performed a panoramic analysis on both human nonalcoholic steatohepatitis (NASH) microarray data and microarray/RNA-seq data from various mouse models of nonalcoholic fatty liver disease NASH/NAFLD with total 4249 genes examined and made the following findings: (i) human NASH and NAFLD mouse models upregulate both cytokines and chemokines; (ii) pathway analysis indicated that human NASH can be classified into metabolic and immune NASH; methionine- and choline-deficient (MCD)+high-fat diet (HFD), glycine N-methyltransferase deficient (GNMT-KO), methionine adenosyltransferase 1A deficient (MAT1A-KO), and HFCD (high-fat-cholesterol diet) can be classified into inflammatory, SAM accumulation, cholesterol/mevalonate, and LXR/RXR-fatty acid β-oxidation NAFLD, respectively; (iii) canonical and noncanonical inflammasomes play differential roles in the pathogenesis of NASH/NAFLD; (iv) trained immunity (TI) enzymes are significantly upregulated in NASH/NAFLD; HFCD upregulates TI enzymes more than cytokines, chemokines, and inflammasome regulators; (v) the MCD+HFD is a model with the upregulation of proinflammatory cytokines and canonical and noncanonical inflammasomes; however, the HFCD is a model with upregulation of TI enzymes and lipid peroxidation enzymes; and (vi) caspase-11 and caspase-1 act as upstream master regulators, which partially upregulate the expressions of cytokines, chemokines, canonical and noncanonical inflammasome pathway regulators, TI enzymes, and lipid peroxidation enzymes. Our findings provide novel insights on the synergies between hyperlipidemia and hypomethylation in establishing TI and promoting inflammation in NASH and NAFLD progression and novel targets for future therapeutic interventions for NASH and NAFLD, metabolic diseases, transplantation, and cancers.


---

## Body

## 1. Introduction

Several metabolic diseases significantly drive the development of cardiovascular disease, nonalcoholic fatty liver disease (NAFLD) [1], and nonalcoholic steatohepatitis (NASH) [2]. NASH is characterized by macrovascular steatosis, hepatocellular ballooning, lobular inflammation, and pericellular fibrosis. NAFLD constitutes a major health concern, and NAFLD prevalence is estimated between 25% and 30% in the general population [3]. However, the exact underlying inflammatory mechanisms that trigger the transition from fatty liver to NASH remain unclear.Innate immune system inflammation is considered a landmark of NASH pathogenesis. Mouse models of NASH have demonstrated that increased danger/pathogen-associated molecular pattern (DAMP/PAMP) receptor and Toll-like receptor (TLR) signaling [4] are associated with disease progression. Moreover, an intracellular DAMP/PAMP receptor and the inflammasome serve as a metabolic stress sensor and bridge metabolic distress to the initiation of pyroptosis, a proinflammatory programmed cell death [4–7]. In addition, inflammasomes serve as sensors of metabolic stresses for vascular, liver, and other inflammations [5, 6, 8–15]. Inflammasome pathways can be classified into the canonical caspase-1 or the noncanonical caspase-11-dependent inflammasome pathways [16]. The canonical inflammasome pathway has been shown to play a significant role in NASH and NAFLD pathogenesis [17, 18]. In addition to the molecular mechanisms of innate immunity, cellular mechanisms also play a significant role in NASH development.The secretome can be defined as the portion of total protein secreted by cells to the extracellular space to maintain tissue homeostasis [19–22]. Treg may use innate immunity-related secretomes such as caspase-1-gasdermin D- (GSDMD-) dependent secretome [23] and caspase-4 (humans)/caspase-11 (mice)-GSDMD-dependent secretome [24] to carry out their functions. However, it remains unknown whether all the cytokines and chemokines in canonically and noncanonically defined secretomes play roles in the progression of NAFLD [25].Innate immune cells can develop an exacerbated immunologic responses and long-term inflammatory phenotypes following brief exposure to endogenous or exogenous insults (the first challenge), which leads to an enhanced inflammatory response after a second challenge, which is known as trained immunity (TI) [4, 26–28]. TI is not only important for host defense and vaccine response but also for chronic inflammatory processes such as cardiovascular inflammation, EC activation [26, 27, 29], and liver injuries [30]. However, the role TI plays in the development of NASH/NAFLD has yet to be determined.Oxidation of cholesterol and phospholipids can lead to lipid peroxidation and cell death [31–35]. Furthermore, lipid peroxidation drives canonical and noncanonical inflammasome-GSDMD-mediated pyroptosis via a glutathione peroxidase 4- (GPX4-) suppressed mechanism [33]. GPX4 is highly expressed in the liver [36], and hepatic ferroptosis [37] initiates inflammation in NASH via GPX4-suppressed manner [38]. However, the role of differential lipid peroxidation enzyme expression in NASH/NAFLD progression remains unclear.DNA methylation plays a significant role in the progression of NAFLD. Previous studies have shown that there are differences in DNA methylation levels of multiple genes between liver samples from NAFLD patients with those from healthy individuals [39, 40], and further studies reveled that hypomethylation of different genes was more than the hypermethylation in advanced NAFLD versus mild NAFLD indicating that genes related to steatohepatitis, fibrosis, and tumorigenesis may be demethylated as NAFLD processes and the methylation levels of these genes may reflect the severity of NAFLD [41].To fill in important knowledge gaps, we performed a panoramic database mining analysis on both human NASH microarray data and microarray data from various NAFLD mouse models and examined a total of 4249 genes by using the strategy we pioneered [9, 42, 43]. We made significant findings, which provide novel insights on the roles of proinflammatory cytokines and chemokines, canonical and noncanonical inflammasome pathways, TI enzymes, and lipid peroxidation in promoting NASH/NAFLD progression.


## 2. Materials and Methods

### 2.1. Expression profile of innatomic and secretomic genes, canonical and noncanonical inflammasome pathway regulators, trained immune genes, and lipid peroxidation genes in patients’ nonalcoholic steatohepatitis (NASH) microarray data and microarray/RNA-seq data from various mouse models of nonalcoholic fatty liver disease

The 10 microarray datasets were collected from the National Institutes of Health- (NIH-) National Center for Biotechnology Information- (NCBI-) Gene Expression Omnibus (GEO) database and analyzed with an online software GEO2R (Figures1(a) and 1(b)). One dataset was from the high-fat diet (HFD) and methionine- and choline-deficient (MCD) diet model (GSE35961). One dataset was from the high-fat-cholesterol diet- (HFCD-) induced diet model. One dataset was from the glycine N-methyltransferase- (GNMT-) KO genetic model, (4) the liver-specific methionine adenosyltransferase 1A- (MAT1A-) KO (GSE63027) genetic model. There were two human NASH microarrays (GSE63067 and GSE17470) [49]. For the mechanism studies, we collected one dataset from caspase-11-KO mice (GSE115094) (Figures 1(a) and 1(b)). One dataset was from caspase-1-KO mice (GSE32515). One human dataset was for trained immunity pathway (GSE24187). For datasets not formatted for GEO2R, analysis was performed using DESeq2 in R Studio as described by Mistry et al. (2015) (Figure 1(c)). In brief, the expression data obtained from NCBI GEO database was converted to an expression set R script element. Differential gene expression was analyzed using the DESeq2 library. The numbers and detailed information of the 7 GEO datasets are listed in Supplementary Tables 2A and 2B.Figure 1
The flowchart of our data mining analyses. (a) Workflow included three parts: (1) identify genes and signaling pathways shared by microarrays from several mouse models of nonalcoholic fatty liver disease (NAFLD), (2) identify genes and signaling pathways common in mouse models of NAFLD and patients with nonalcoholic steatohepatitis (NASH), and (3) identify significant inflammatory mechanisms underlying the pathogenesis and progression of NASH/NAFLD. (b) Process automations using Microsoft Excel Macros, which have significantly facilitated the database mining processes comparing to that of our previous database mining papers. (c) GEO datasets without the GEO2R function were analyzed using DESeq2 library in R Studio. R script code used in C: Meeta Mistry, C. Titus Brown, jessicalumian, & tug65470 (2021, July 14), tug65470/msu_ngs2015: DESeq2 analysis of microarray and RNA-seq datasets (version v1.0.0), Zenodo, http://doi.org/10.5281/zenodo.5102949.
(a)(b)(c)


### 2.2. Ingenuity Pathway Analysis (IPA)

IPA was used to characterize the clinical relevance and molecular and cellular functions related to the identified genes. Differentially expressed genes were identified and uploaded into IPA for analysis. Core and canonical pathway analysis was used to identify molecular and cellular pathways [14, 15].


### 2.3. Statistical Analysis

We applied a statistical method similar to a meta-analysis [14, 27]. GEO dataset integrity used 7 housekeeping genes (Supplementary Tables 1 and 2) [44]. The target genes with expression changes more than 1.5-fold were defined as upregulated genes, while genes with their expression decreased more than 1.5-fold were defined as downregulated genes. Fold change calculations were performed in Excel as previously described by our group. Briefly, differential gene expression data were imported into Excel and then analyzed by four Excel Macros: (1) organize the DEG data from the DESeq2 and GEO2R, (2) filter data for significant p values, (3) filter data for foldchange±1.5, and (4) retrieve expression data for various genes of interest (Figure 1(c)).


## 2.1. Expression profile of innatomic and secretomic genes, canonical and noncanonical inflammasome pathway regulators, trained immune genes, and lipid peroxidation genes in patients’ nonalcoholic steatohepatitis (NASH) microarray data and microarray/RNA-seq data from various mouse models of nonalcoholic fatty liver disease

The 10 microarray datasets were collected from the National Institutes of Health- (NIH-) National Center for Biotechnology Information- (NCBI-) Gene Expression Omnibus (GEO) database and analyzed with an online software GEO2R (Figures1(a) and 1(b)). One dataset was from the high-fat diet (HFD) and methionine- and choline-deficient (MCD) diet model (GSE35961). One dataset was from the high-fat-cholesterol diet- (HFCD-) induced diet model. One dataset was from the glycine N-methyltransferase- (GNMT-) KO genetic model, (4) the liver-specific methionine adenosyltransferase 1A- (MAT1A-) KO (GSE63027) genetic model. There were two human NASH microarrays (GSE63067 and GSE17470) [49]. For the mechanism studies, we collected one dataset from caspase-11-KO mice (GSE115094) (Figures 1(a) and 1(b)). One dataset was from caspase-1-KO mice (GSE32515). One human dataset was for trained immunity pathway (GSE24187). For datasets not formatted for GEO2R, analysis was performed using DESeq2 in R Studio as described by Mistry et al. (2015) (Figure 1(c)). In brief, the expression data obtained from NCBI GEO database was converted to an expression set R script element. Differential gene expression was analyzed using the DESeq2 library. The numbers and detailed information of the 7 GEO datasets are listed in Supplementary Tables 2A and 2B.Figure 1
The flowchart of our data mining analyses. (a) Workflow included three parts: (1) identify genes and signaling pathways shared by microarrays from several mouse models of nonalcoholic fatty liver disease (NAFLD), (2) identify genes and signaling pathways common in mouse models of NAFLD and patients with nonalcoholic steatohepatitis (NASH), and (3) identify significant inflammatory mechanisms underlying the pathogenesis and progression of NASH/NAFLD. (b) Process automations using Microsoft Excel Macros, which have significantly facilitated the database mining processes comparing to that of our previous database mining papers. (c) GEO datasets without the GEO2R function were analyzed using DESeq2 library in R Studio. R script code used in C: Meeta Mistry, C. Titus Brown, jessicalumian, & tug65470 (2021, July 14), tug65470/msu_ngs2015: DESeq2 analysis of microarray and RNA-seq datasets (version v1.0.0), Zenodo, http://doi.org/10.5281/zenodo.5102949.
(a)(b)(c)


## 2.2. Ingenuity Pathway Analysis (IPA)

IPA was used to characterize the clinical relevance and molecular and cellular functions related to the identified genes. Differentially expressed genes were identified and uploaded into IPA for analysis. Core and canonical pathway analysis was used to identify molecular and cellular pathways [14, 15].


## 2.3. Statistical Analysis

We applied a statistical method similar to a meta-analysis [14, 27]. GEO dataset integrity used 7 housekeeping genes (Supplementary Tables 1 and 2) [44]. The target genes with expression changes more than 1.5-fold were defined as upregulated genes, while genes with their expression decreased more than 1.5-fold were defined as downregulated genes. Fold change calculations were performed in Excel as previously described by our group. Briefly, differential gene expression data were imported into Excel and then analyzed by four Excel Macros: (1) organize the DEG data from the DESeq2 and GEO2R, (2) filter data for significant p values, (3) filter data for foldchange±1.5, and (4) retrieve expression data for various genes of interest (Figure 1(c)).


## 3. Results

### 3.1. Human NASH and NAFLD Mouse Models Upregulate Both Cytokines and Chemokines Classified in the Innate Immune Database and Canonical Secretome

At least 11 NAFLD mouse models have been established [45, 46] including diet-induced and genetic-induced models [47, 48], but no model accurately reflects the disease’s progression in humans [3]. We include four mouse models that best represent the metabolic and inflammatory characteristics of human NAFLD [3, 46–48].We hypothesized that diet- and genetically induced NAFLD mouse models can share common pathways critical to the development of the disease. To test this, we collected microarray datasets from NAFLD mouse models and human NASH (Supplementary Table2a). The datasets we analyzed include (1) the high-fat diet (HFD) and methionine- and choline-deficient (MCD) diet model (GSE35961) [49], (2) the high-fat-cholesterol diet- (HFCD-) induced diet model [50], (3) the glycine N-methyltransferase- (GNMT-) KO genetic model, (4) the liver-specific methionine adenosyltransferase 1A- (MAT1A-) KO (GSE63027) [51] genetic model, and (5) two human NASH microarrays (GSE63067 [51] and GSE17470 [52]). The qualities of the datasets were assessed using a few housekeeping genes (Supplementary Tables 1a and 1b).Cytokines and chemokines [53–55] play significant roles in promoting the onset and progression of NAFLD [55]; however, their expression changes in human NASH and NAFLD mouse models remained poorly characterized. Cytokines and chemokines in the liver are soluble mediators of local and systemic inflammation (Figure 2(a)) [8, 43, 56, 57]. To analyze liver inflammation, we used IPA and classified 53 cytokines and chemokines from 1376 innate immune genes [58, 59] from the Innate Immunity Database [58, 60]. Human NASH upregulated six (11.3%) and four (7.5%) out of 53 innate immune cytokines and chemokines, respectively (Figure 2(b)). Of note, upregulations of C-X-C motif chemokine ligand 10 (CXCL10) and CKLF-like MARVEL transmembrane domain-containing 2 (CMTM2) were shared in these two human NASH datasets. MCD+HFD and MAT1A-KO upregulated 13 (24.5%) and 6 (11.3%) out of 53 innate immune cytokines and chemokines, respectively (Figure 2(c)). MCD+HFD, MAT1A-KO, GNMT-KO, and HFCD downregulated 9 (17%), 3 (5.7%), 1 (1.9%), and 12 (22.6%) out of 53 innate immune cytokines and chemokines, respectively.Figure 2
The cytokine and chemokine expressions are increased in both human nonalcoholic steatohepatitis (NASH) and mouse models of nonalcoholic fatty liver disease (NAFLD). Ingenuity Pathway Analysis results showed that 53 cytokines and chemokines were sorted out from 1376 innate immune genes (innatome, from the Innate Immune Database (https://www.innatedb.com/); also, see our paper (PMID: 33628851)) and 123 cytokines and chemokines were sorted out from 2641 canonical secretome (with signal peptide) genes (from the Human Protein Atlas database (https://www.proteinatlas.org/); also, see our paper (PMID: 32179051)). (a) A schematic figure showed the pathogenic effects of cytokines released from NASH on different organs. Two human NASH datasets (GSE63067 and GSE17470) and four mouse models of NAFLD datasets (GSE35961, GSE63027, GSE63027, and GSE53381) were analyzed. (b, c) The 53 innatome cytokines and chemokines were analyzed in human NASH and mouse models of NAFLD. (d) Venn diagram showed the significant overlapped regulated innatome cytokines and chemokines in human NASH and mouse models of NAFLD. The cytokine and chemokine expressions are increased in both human nonalcoholic steatohepatitis (NASH) and mouse models of nonalcoholic fatty liver disease (NAFLD). (e, f) The 123 canonical secretomic cytokines and chemokines were analyzed in human NASH and mouse models of NAFLD. (g) Venn diagram showed the overlapped significant regulated canonical secretomic cytokines and chemokines in human NASH and mouse models of NAFLD. The cytokine and chemokine expressions are increased in both human nonalcoholic steatohepatitis (NASH) and mouse models of nonalcoholic fatty liver disease (NAFLD). (h) The 16 cytokines and chemokines (innatome and secretome) were upregulated (upregulated in GSE63067 or GSE17470) in human NASH. PMID: 34084175. (i) The 16 cytokines and chemokines (innatome and secretome) were upregulated (GSE35961 or GSE63027) in mouse models of NAFLD. Secretome gene list detailed in our previous publication PMID: 34084175.
(a)(b)(c)(d)(e)(f)(g)(h)(i)To determine whether NAFLD mouse models share innate immune cytokines and chemokines with human NASH, we performed Venn diagram analysis. Six innate immune cytokines and chemokines including CXCL10, tumor necrosis factor (ligand) superfamily, member 13b (TNFSF13B), lymphotoxin beta (LTB, TNFSF3), TNFSF10, interleukin 1 receptor antagonist (IL-1RN), and nicotinamide phosphoribosyltransferase (NAMPT) were shared between human NASH and NAFLD mouse models (Figure2(d)). Two NAFLD mouse models upregulated 22 (78.6%) specific innate immune cytokines and chemokines; and human NASH upregulated two (25%) specific cytokines and chemokines C-C motif chemokine ligands 20 (CCL20) and CMTM2.To ensure a comprehensive analysis of liver inflammation in human NASH and NAFLD mouse models, we collect 2641 canonical secretome genes [61]. The 123 cytokines and chemokines were classified with IPA out of 2641 human canonical secretome genes. Human NASH upregulated 4 (3.3%) and 11 (8.9%) out of 123 canonical secretomic cytokines and chemokines (Figure 2(e)). Of note, upregulations of CXCL10 and thymic stromal lymphopoietin (TSLP) were shared in the two human NASH datasets. MCD+HFD, MAT1A-KO, GNMT-KO, and HFCD models upregulated 21 (17.1%), 7 (5.7%), 0, and 3 (2.4%) (Figure 2(f)) and downregulated 10 (8.1%), 4 (3.3%), 1 (0.8%), and 10 (8.1%) out of 123 canonical secretomic cytokines and chemokines, respectively.To determine whether NAFLD mouse models share canonical secretomic cytokines and chemokines with human NASH, we performed Venn diagram analysis. Seven canonical secretomic cytokines and chemokines including TSLP, CXCL10, CCL28, CCL25, CXCL14, IL21, and IL1RN were shared between NAFLD mouse models and human NASH (Figure2(g)). By comparison, three NAFLD mouse models upregulated 32 (82.1%) specific canonical secretomic cytokines and chemokines; and human NASH upregulated 8 (53.3%) specific cytokines and chemokines including erythropoietin (EPO), IL17F, CCL8, IL27, CCL21, CXCL12, CCL15, and CCL20.Caspase-1-deficient mice are protected from high fat-induced hepatic steatosis, inflammation, and early fibrogenesis [62], suggesting that pyroptosis-related caspase-1/inflammasomes [6] promote NASH progression. However, an important question remains whether noncanonical inflammasome pathway-related caspase-4/11 [63] plays any role in progression of NAFLD.To determine whether NASH/NAFLD upregulated cytokines and chemokines classified into noncanonical secretomes, we collected four types of secretomes with total 11,385 proteins including canonical secretome (2641 proteins) [61], caspase-1-dependent noncanonical secretome (961 proteins) [23], caspase-4-dependent noncanonical secretome (1223 proteins) [24], and exosome secretome (6560 proteins) [25]. We combined human NASH and NAFLD mouse upregulated cytokines and chemokines and performed Venn diagram analysis. Among 15 human NASH-upregulated cytokines and chemokines, two cytokines CMTM2 and LTB were nonclassified, 12 out of 15 (80%) were classified in the canonical secretome, TNFSF10 was classified in canonical and exosome secretomes, and NAMPT was classified in canonical and caspase-4 secretomes (Figure 2(h)). Of note, since 123 cytokines and chemokines identified with the IPA are derived from the canonical secretome, it was expected that the high percentage (80%) NASH-upregulated cytokines and chemokines were from the canonical secretome. Among 26 NAFLD mouse-upregulated cytokines and chemokines, one was the nonclassified, 15 were classified in the canonical secretome, seven were classified in canonical and exosome secretomes, NAMPT was classified in canonical and caspase-4 secretomes, and two cytokines were classified in canonical, exosomes, and caspase-4 secretomes (Figure 2(i)).Taken together, our results have demonstrated innate immune secretomic cytokines and chemokines suggesting that human NASH and NAFLD mouse models share innate immune mechanisms and there is a significant role of exosomes and caspase-4 secretomic in driving liver and systemic inflammations.


### 3.2. Human NASH Can Be Classified into Metabolic and Immune NASH; MCD+HFD, GNMT-KO, MAT1A-KO, and HFCD Can Be Classified into Inflammatory, SAM Accumulation, Cholesterol/Mevalonate, and LXR/RXR-Fatty Acidβ-Oxidation NAFLD, Respectively

To determine the signaling pathways mediating the transcriptomic changes in human NASH and NAFLD mouse models, we performed IPA analysis for the top 20 pathways for significantly modulated genes in each dataset. Human NASH dataset GSE17470 had top upregulated pathways including phospholipases, liver X receptor (LXR)/retinoid X receptor (RXR) activation, fatty acidβ-oxidation I, and leukocyte extravasation, among which leukocyte extravasation was shared by the second human NASH and MCD+HFD, LXR/RXR activation was shared by HFCD, and fatty acid β-oxidation was shared by GNMT-KO and HFCD (Table 1, Supplementary Figure 1).Table 1
Ingenuity Pathway Analysis (IPA) for top 20 pathways of significantly modulated genes in each dataset of human NASH studies and mouse models of NAFLD. Leukocyte extravasation was shared by both human NASH studies and in the MCD+HFD mouse model; LXR/RXR activation was shared by HFCD model of NAFLD; and fatty acidβ-oxidation was shared by the GNMT-KO and HFCD models of NAFLD. IL-8 signaling was shared by the second human NASH study and the MCD+HFD model. Fatty acid β-oxidation I and stearate biosynthesis I were shared by the first human NASH study and GNMT-KO and HFCD models. Nicotine degradation II was shared by the first human NASH study and the MAT1A-KO model.
(a)HumanMouseGEO IDGSE17470GSE63067GSE35961GSE63027GSE63027GSE53381ComparisonNASH vs. healthyNASH vs. healthyMCD+HFDGNMT-KOMAT1A-KOHFCDPathwaysSerotonin degradation↑↑Ethanol degradation II↑Synaptic long-term depression↑Noradrenaline and adrenaline degradation↑Sperm motility↑Phospholipases↑LXR/RXR activation↑↑Fatty acidβ-oxidation I↑↑↑Nicotine degradation II↑↑↑eNOS signaling↑Leukocyte extravasation signaling↑↑↑p70S6K signaling↑Retinol biosynthesis↑Stearate biosynthesis I (animals)↑↑Histamine degradation↑Aldosterone signaling in epithelial cells↑HIPPO signaling↑cAMP-mediated signaling↑Synaptic long-term potentiation↑Fcγ receptor-mediated phagocytosis in macrophages and monocytes↑↑IL-8 signaling↑↑Rac signaling↑NF-κB activation by viruses↑Salvage pathways of pyrimidine ribonucleotides↑Gαq signaling↑Pyridoxal 5′-phosphate salvage pathway↑GP6 signaling pathway↑Colorectal cancer metastasis signaling↑ILK signaling↑Integrin signaling↑mTOR signaling↑Phospholipase C signaling↑Cardiac hypertrophy signaling (enhanced)↑↑ERK5 signaling↑Sphingosine-1-phosphate signaling↑Small cell lung cancer signaling↑Gα12/13 signaling↑Signaling by rho family GTPases↑Superpathway of cholesterol biosynthesis↑PPAR signalingNicotine degradation III↑Melatonin degradation ISuperpathway of melatonin degradationXenobiotic metabolism CAR signaling pathwayXenobiotic metabolism PXR signaling pathwayCholesterol biosynthesis I↑Cholesterol biosynthesis II (via 24,25-dihydrolanosterol)Cholesterol biosynthesis III (via desmosterol)↑(b)HumanMouseGEO IDGSE17470GSE63067GSE35961GSE63027GSE63027GSE53381ComparisonNASH vs. healthyNASH vs. healthyMCD+HFD vs. NCDGNMT-KO vs. WTMAT1A-KO vs. WTHFCD vs. NCDPathwaysSuperpathway of geranylgeranyl diphosphate biosynthesis I (via mevalonate)Estrogen biosynthesis↑Acetone degradation I (to methylglyoxal)↑Mevalonate pathway I↑Tryptophan degradation III (eukaryotic)Antioxidant action of vitamin C↑Oxidative phosphorylation↑Glutaryl-CoA degradation↑Valine degradation I↑Isoleucine degradation I↑RhoA signaling↑Mitotic roles of polo-like kinase↑Triacylglycerol biosynthesis↑TCA cycle II (eukaryotic)↑Role of CHK proteins in cell cycle checkpoint control↑Sumoylation pathway↑Apelin adipocyte signaling pathway↑↑Glutathione-mediated detoxification↑Tetrahydrofolate salvage from 5,10-methenyltetrahydrofolate↑CMP-N-acetylneuraminate biosynthesis I (eukaryotes)↑EIF2 signaling↑tRNA charging↑Folate transformations I↑Endocannabinoid cancer inhibition pathway↑Pyrimidine ribonucleotide interconversion↑Pyrimidine ribonucleotide de novo biosynthesis↑NAD salvage pathway II↑PPARα/RXRα activation↑NRF2-mediated oxidative stress response↑↑Glutathione redox reactions I↑Pancreatic adenocarcinoma signaling↑Estrogen-mediated S-phase entry↑Type II diabetes mellitus signaling↑PCP pathway↑PDGF signaling↑NF-κB signaling↑Neuroinflammation signaling pathway↑Tec kinase signaling↑B cell receptor signaling↑Osteoarthritis pathway↑Th17 activation pathway↑Role of pattern recognition receptors in recognition of bacteria and viruses↑Apelin endothelial signaling pathway↑Cholecystokinin/gastrin-mediated signaling↑T cell exhaustion signaling pathway↑Endothelin-1 signaling↑B cell activating factor signaling↑RANK signaling in osteoclasts↑Toll-like receptor signaling↑Renin-angiotensin signaling↑The GNMT-KO (GSE63027) had top upregulated pathways including fatty acidβ-oxidation I, stearate biosynthesis I, antioxidant action of vitamin C, oxidative phosphorylation, glutaryl-CoA degradation kinase, triacylglycerol biosynthesis, and tricarboxylic acid cycle (TCA) cycle II, among which fatty acid β-oxidation I and stearate biosynthesis I were shared with the first human NASH (GSE17470). The MAT1A-KO (GSE63027) had top upregulated pathways including nicotine degradation II, superpathway of cholesterol biosynthesis, cholesterol biosynthesis I, cholesterol biosynthesis III, mevalonate pathway I, and NAD salvage pathway II, among which nicotine degradation II was shared with the first human NASH.The HFCD model (GSE53381) has only four upregulated pathways including LXR/RXR activation and fatty acidβ-oxidation I, among which LXR/RXR activation was shared with the first human NASH, fatty acid β-oxidation II was shared with the first human NASH and GNMT-KO, and nicotine degradation II was shared with MAT1A-KO.The second human NASH (GSE63067) had top upregulated pathways including NF-κB signaling and neuroinflammation (Table 1). The MCD+HFD model (GSE35961) had top upregulated pathways including leukocyte extravasation, antibody Fc fragment gamma (Fcγ) receptor-mediated phagocytosis, IL-8 signaling, Rac signaling, and glucose 6-phosphate (G6P) signaling, among which IL-8 signaling was shared with the second human NASH.Taken together, our results have demonstrated that the two human NASH datasets have diversified signaling pathways, which allow us to classify the first one into LRX/RXR activation-related metabolic NASH and the second NASH into NF-κB signaling [64]/T helper cell 17 (Th17) [56, 65] inflammatory NASH. IPA of upregulated genes in each mouse models of NASH and human NASH identified top pathways for significantly modulated genes. The major differences between metabolic NASH and immune NASH include the following: first, numerous inflammation and immune pathways are identified on the top pathway list, and second, metabolic reprogramming in any given mouse models affects four metabolic pathways including (1) increased glycolysis, (ii) glutaminolysis, (iii) increased accumulation of tricarboxylic acid cycle metabolites and acetyl-coenzyme A production, and (iv) increased mevalonate synthesis. Furthermore, this classification provides a criterion for mouse model stratification and selection, providing novel insight of which mouse models best represent various human cases of NASH/NAFLD.


### 3.3. Canonical and Noncanonical Inflammasome Pathways Play Differential Roles in the Pathogenesis of NASH/NAFLD

To examine the differential roles of caspase-1- and caspase-11-dependent pyroptosis, we collected 90 canonical and 14 noncanonical inflammasome pathway regulator genes and determined the expression changes of these inflammasome regulators in four NAFLD mouse models.MCD+HFD upregulated 18 (20%) including caspase-1 and downregulated 11 (12.2%) out of 90 canonical inflammasome pathway regulators. MAT1A-KO upregulated 5 (5.6%) and downregulated 4 (4.4%) out of 90 canonical inflammasome pathway regulators. GNMT-KO upregulated 1 (1.1%) and downregulated 2 (2.2%) out of 90 canonical inflammasome pathway regulators. HFCD upregulated 7 (7.8%) and downregulated 26 (28.9%) out of 90 canonical inflammasome pathway regulators (Figure3(b)). In addition, MCD+HFD downregulated 6 out of 14 (42.8%) noncanonical inflammasome pathway regulators. MAT1A-KO upregulated 3 (21.4%) and downregulated 3 (21.4%) out of 14 noncanonical inflammasome pathway regulators. GNMT-KO downregulated 1 out of 14 (7.1%) noncanonical inflammasome pathway regulators. HFCD upregulated 1 (7.1%) and downregulated 9 (64.3%) out of 14 noncanonical inflammasome pathway regulators (Figure 3(c)).Figure 3
Mouse pyroptosis pathway regulators. (a) Schematic representation of canonical vs. noncanonical pyroptosis. (b) Mouse canonical inflammasome pathway regulators (90 genes, KEGG pathway hsa04621) were differentially expressed in NASH/NAFLD models. The MCD+HFD model upregulated 18 out of 90 (20%) canonical inflammasome pathway regulators including caspase-1 and downregulated 11 out of 90 (12.2%) canonical inflammasome pathway regulators. The MAT1A-KO model upregulated 5 out of 90 (5.6%) and downregulated 4 out of 90 (4.4%) canonical inflammasome pathway regulators. The GNMT-KO model upregulated 1 out of 90 (1.1%) and downregulated 2 out of 90 (2.2%) canonical inflammasome pathway regulators. The HFCD model upregulated 7 out of 90 (7.8%) and downregulated 26 out of 90 (28.9%) canonical inflammasome pathway regulators. (c) Mouse noncanonical inflammasome pathway regulators (14 genes, KEGG pathway hsa04621) were differentially expressed in NASH/NAFLD. The MCD+HFD model upregulated zero out of 14 (0%) and downregulated 6 out of 14 (42.8%) noncanonical inflammasome pathway regulators. The MAT1A-KO model upregulated 3 out of 14 (21.4%) and downregulated 3 out of 14 (21.4%) noncanonical inflammasome pathway regulators. The GNMT-KO model upregulated 0 out of 14 (0%) and downregulated 1 out of 14 (7.1%) noncanonical inflammasome pathway regulators. The HFCD model upregulated 1 out of 14 (7.1%) and downregulated 9 out of 14 (64.3%) noncanonical inflammasome pathway regulators.
(a)(b)(c)Taken together, these results have demonstrated that the MCD+HFD upregulates more canonical inflammasome regulators than the other three NAFLD models. Conversely, the MAT1A-KO upregulates more noncanonical inflammasome regulators including caspase-4 than other models, suggesting that canonical and noncanonical inflammasome pathways are regulated differentially in various NAFLD models. The HFCD downregulates more than upregulates canonical and noncanonical inflammasome regulators, suggesting that HFCD model uses more focused set of regulators in modulating the expression of these regulators.


### 3.4. Trained Immunity (TI) Enzymes Are Significantly Upregulated in Human NASH and NAFLD Mouse Models; the HFD Upregulates TI Enzymes More Than Cytokines, Chemokines, and Canonical and Noncanonical Inflammasome Regulators; and Statins Promote rather than Suppress TI Enzyme Expressions

We hypothesized that the expressions of TI enzymes are modulated in NASH and NAFLD to establish metabolic reprogramming and achieve innate immune memory [66]. Three metabolic pathways play significant roles in establishing TI (Figure 4). To examine this, we collected 71 glycolysis enzyme genes, 24 acetyl-CoA generation enzyme genes, and 10 mevalonate genes [27, 29, 30, 67, 68].Figure 4
Schematic representation of fundamental mechanisms in trained immunity. Danger-associated molecular patterns (DAMPs) can bind to their corresponding pattern recognition receptors (PRRs) and alter metabolic pathways including increased glycolysis, increased acetyl-CoA generation, and mevalonate synthesis in cholesterol pathway, leading to histone modifications that enable chromatin regions to be more open for transcription. Increased gene expression of proinflammatory cytokines/chemokines and enhanced proinflammatory innate immune response against pathogens during secondary exposure. Figure created with Smart Servier Medical Art and Omnigraffle.The first human NASH upregulated 25 (35.2%) and downregulated 5 (7%) out of 71 glycolysis enzymes genes (Table2(a)). The second human NASH upregulated 5 out of 71 (7%) genes. The MCD+HFD upregulated 20 (28.2%) and downregulated 24 (33.8%) out of 71 genes. The MAT1A-KO upregulated 6 (8.5%) and downregulated 13 (18.3%) out of 71 genes. The GNMT-KO model upregulated 7 (9.9%) and downregulated 4 (5.6%) out of 71 genes. The HFCD upregulated 17 (23.9%) and downregulated 8 (11.3%) out of 71 genes.Table 2
(a) The 61 out of 71 mRNAs (85.9%) of glycolysis pathway enzyme genes are regulated in human and mouse NASH/NAFLD. Glycolysis pathway enzyme gene expression (PMID: 31153039) for human NASH studies (GSE17470 and GSE63067), MCD+HFD (GSE35961), HFC (GSE53381), MAT1A-KO (GSE63027), and GNMT-KO (GSE63027). The first human NASH study (GSE17470) upregulated 25 out of 71 (35.2%) and downregulated 5 out of 71 (7%) glycolysis enzymes. The second human NASH study (GSE63067) upregulated 5 out of 71 (7%) and downregulated 0 out of 71 (0%) glycolysis enzymes. The MCD+HFD model upregulated 20 out of 71 (28.2%) and downregulated 24 out of 71 (33.8%) glycolysis enzymes. The MAT1A-KO model upregulated 6 out of 71 (8.5%) and downregulated 13 out of 71 (18.3%) glycolysis enzymes. The GNMT-KO model upregulated 7 out of 71 (9.9%) and downregulated 4 out of 71 (5.6%) glycolysis enzymes. The HFCD model upregulated 17 out of 71 (23.9%) and downregulated 8 out of 71 (11.3%) glycolysis enzymes. The 71 glycolysis pathway genes are shown in Supplementary Table3a. (b) The 13 out of 24 mRNA (54.2%) of acetyl-CoA generating enzyme are regulated in human and mouse NASH/NAFLD. Acetyl-CoA generating pathway enzyme gene expression (PMID: 31153039) for human NASH studies (GSE17470 and GSE63067), MCD+HFD (GSE35961), HFC (GSE53381), MAT1A-KO (GSE63027), and GNMT-KO (GSE63027). The first human NASH study upregulated 4 out of 24 (16.6%) and downregulated 0 out of 24 (0%) acetyl-CoA generation enzymes. The second human NASH study upregulated 1 out of 24 (4.2%) and downregulated 0 out of 24 (0%) acetyl-CoA generation enzymes. The MCD+HFD model upregulated 3 out of 24 (12.5%) and downregulated 7 out of 24 (29.2%) acetyl-CoA generation enzymes. The MAT1A-KO model upregulated 1 out of 24 (4.2%) and downregulated 4 out of 24 (16.7%) acetyl-CoA generation enzymes. The GNMT-KO model upregulated 2 out of 24 (8.3%) and downregulated 0 out of 24 (0%) acetyl-CoA generation enzymes. The HFCD model upregulated 7 out of 24 (29.2%) and downregulated 2 out of 24 (8.3%) acetyl-CoA generation enzymes. The 24 acetyl-CoA generating enzyme genes are shown in Supplementary Table 3b. (c) The 10 out of 10 mRNAs (100%) of mevalonate pathway enzyme genes are regulated in human and mouse NASH/NAFLD. Mevalonate pathway enzyme gene expression (PMID 31153039 and PMID 29328908) for human NASH studies (GSE17470 and GSE63067), MCD+HFD (GSE35961), HFC (GSE53381), MAT1A-KO (GSE63027), and GNMT-KO (GSE63027). The first human NASH study upregulated 2 out of 10 (20%) and downregulated 3 out of 10 (30%) mevalonate pathway enzymes. The second human NASH study upregulated 1 out of 10 (10%) and downregulated 0 out of 10 (0%) mevalonate pathway enzymes. The MCD+HFD model upregulated 7 out of 10 (70%) and downregulated 1 out of 10 (10%) mevalonate pathway enzymes. The MAT1A-KO model upregulated 4 out of 10 (40%) and downregulated 1 out of 10 (10%) mevalonate pathway enzymes. The GNMT-KO model upregulated 3 out of 10 (30%) and downregulated 0 out of 24 (0%) mevalonate pathway enzymes. The HFCD model upregulated 7 out of 10 (70%) and downregulated 0 out of 10 (0%) mevalonate pathway enzymes. The 10 mevalonate pathway enzyme genes are shown in Supplementary Table 3C.
(a)PatientsMouseGEO IDGSE17470GSE63067GSE35961GSE63027GSE63027GSE53381ComparisonNASH vs. healthyNASH vs. healthyMCD+HFD vs. NCDMAT1A-KO vs. WTGNMT-KO vs. WTHFCD vs. NCDGene symbolACSS14.47-2.67ACSS2-2.55-1.623.45ADH72.761.641.481.39AHD1A7.97ADH1B2.82ADH1C4.09ALDH1A36.34ALDH1B14.381.622.24-1.63ALDH26.33-1.70-1.501.59ALDH3A22.032.05-2.875.86ALDH3B12.49ALDH7A13.81-4.071.65ALDOA3.69ALDOC2.211.641.323.96DLAT-2.861.40-2.341.401.36ENO33.934.11FBP23.13HK12.73HK22.10HKDC1-1.855.361.60LDHA-1.79LDHB4.99PCK27.054.38-1.33PDHA11.91PDHB-1.39-1.72PFKFB13.641.631.62PFKFB42.691.62PFKP7.30-2.27PGAM11.44PGK1-6.501.75PGM11.682.701.71-1.64PKM5.27-2.76ADH441.77-2.88-1.51-2.01ADH54.73-1.331.45ADH65.42-4.71ADPGK-2.07-1.82ALDH3B2-1.67ALDH9A1-1.351.92ALDOB2.051.50-1.661.30BPGM-1.32-1.371.32DLD2.58-2.021.68ENO1-1.551.36ENO2-2.13FBP12.24-2.231.49G6PC-5.06-1.59GALM2.16-1.721.37GAPDH-1.61GCK17.03-33.70GPI1.95HK32.71-3.37LDHC-2.04PANK1-1.391.77PCK1-1.88-1.45PFKFB21.74PFKFB31.982-1.65-2.48PFKL1.738-1.37PFKM-1.98-1.49-1.31PGAM222.09-2.02PGM21.52PKLR1.738-1.80SLC2A29.160-1.761.36% up25/71 (35.2%)5/71 (7%)20/71 (28.2%)6/71 (8.5%)7/71 (9.9%)17/71 (23.9%)% down5/71 (7%)0/71 (0%)24/71 (33.8%)13/71 (18.3%)4/71 (5.6%)8/71 (11.3%)(b)HumanMouseGEO IDGSE17470GSE63067GSE35961GSE63027GSE63027GSE53381ComparisonNASH vs. healthyNASH vs. healthyMCD+HFD vs. NCDMAT1A-KO vs. WTGNMT-KO vs. WTHFCD vs. NCDGene symbolBDH12.331.89ACSS14.47-2.67ALDH26.33-1.70-1.501.59ACSS2-2.55-1.623.45ACAA22.83-1.511.381.49HADH1.501.37ADH1B2.82GOT11.79ACO1-1.411.45ACLY-2.40-1.361.75GLS3.27-2.40IDH12.83-1.40-1.451.60GLUD1-1.55% up4/24 (16.6%)1/24 (4.16%)3/24 (12.5%)1/24 (4.2%)2/24 (8.3%)7/24 (29.2%)% down0/24 (0%)0/24 (0%)7/24 (29.2%)4/24 (16.7%)0/24 (0%)2/24 (8.3%)(c)HumanMouseGEO IDGSE17470GSE63067GSE35961GSE63027GSE63027GSE53381ComparisonNASH vs. healthyNASH vs. healthyMCD+HFD vs. NCDMAT1A-KO vs. WTGNMT-KO vs. WTHFCD vs. NCDGene symbolIDI12.7781.5710.974IDI2-1.655HMGCS1-2.4401.8952.647.177MVD-1.7954.0031.3815.819MVK4.0965.4511.8583.895HMGCR2.0911.7642.663ACAT12.1351.3741.528PMVK1.453.366ACAT21.63HMGCS27.434-1.908-1.7381.681% up2/10 (20%)1/10 (10%)7/10 (70%)4/10 (40%)3/10 (30%)7/10 (70%)% down3/10 (30%)0/10 (10%)1/10 (10%)1/10 (10%)0/10 (0%)0/10 (0%)Furthermore, the first human NASH upregulated 4 out of 24 (16.6%) acetyl-CoA generation enzyme genes. The second human NASH upregulated 1 out of 24 (4.2%) genes. The MCD+HFD upregulated 3 (12.5%) and downregulated 7 (29.2%) out of 24 genes. The MAT1A-KO upregulated 1 (4.2%) and downregulated 4 (16.7%) out of 24 genes. The HFCD upregulated 7 (29.2%) and downregulated 2 (8.3%) out of 24 genes (Table2(b)).In addition, the first human NASH upregulated 2 (20%) and downregulated 3 (30%) out of 10 mevalonate pathway enzyme genes. The second human NASH upregulated 1 out of 10 (10%) gene. The MCD+HFD upregulated 7 (70%) and downregulated 1 (10%) out of 10 genes. The MAT1A-KO upregulated 4 (40%) and downregulated 1 (10%) out of 10 genes. The GNMT-KO upregulated 3 out of 10 (30%) genes. The HFCD upregulated 7 out of 10 (70%) genes (Table2(c)).Statins (β-hydroxy β-methylglutaryl-CoA (HMG-CoA) reductase inhibitors) block cholesterol synthesis and mevalonate generation, preventing TI induction [69]; and statin treatments do not revert the TI phenotype in patients with familiar hypercholesterolemia [70]. We hypothesized that statins may modulate TI enzyme expression in NAFLD models. For this, we collected two microarray datasets of statin-treated human primary hepatocytes [71].Among the 2004 differentially expressed genes, 16 genes were overlapped with 99 TI enzyme genes (16.2%) [27], of which 15 (15.2%) genes were upregulated and one (1%) gene was downregulated by atorvastatin (Figures 5(a)–5(c)). In addition, among 2448 differentially expressed genes in rosuvastatin-treated human primary hepatocytes, 19 genes were overlapped with 99 TI enzyme genes (19.2%), of which 15 (15.2%) genes were upregulated and 4 (4%) genes were downregulated by rosuvastatin. Then, we used Venn diagram to determine the causative effects of statins in TI enzyme genes upregulated in NAFLD models. The results showed that the atorvastatin treatment overlapped with 13 TI enzyme genes upregulated in human and mouse NASH/NAFLD models. Similarly, rosuvastatin treatment overlapped with 15 TI enzyme genes upregulated in human and mouse NAFLD models. These results have demonstrated that TI inhibitors, statins, promote the expression of TI enzyme genes not only in mevalonate pathways but also in glycolysis pathways and acetyl-CoA generation pathway.Figure 5
Trained immunity enzyme genes were differentially expressed in statin treatment (HMG-CoA reductase inhibitor) in the mevalonate synthesis pathway and modulated in human and mouse NASH/NAFLD. (a) Venn diagram showed the overlapped trained immunity enzyme genes (PMID: 31153039) with differentially expressed genes (DEG) in statin (atorvastatin and rosuvastatin) treatment (GSE24187). (b) List of trained immunity enzyme genes differentially expressed in atorvastatin and rosuvastatin and upregulated in human and mouse NASH/NAFLD models. (c) Venn diagram showed that the upregulated trained immunity enzyme genes in human and mouse NASH/NAFLD overlapped with the differentially expressed genes in atorvastatin and rosuvastatin treatment.
(a)(b)(c)


### 3.5. Expression Analysis of Lipid Peroxidation Enzymes Indicates That There Are Two Representative NASH/NAFLD Models: The MCD+HFD Model Is Proinflammatory Cytokine- and Canonical and Noncanonical Inflammasome-Upregulated Model; in Contrast, the HFCD Model Is Lipid Peroxidation Enzyme- and TI Enzyme-Upregulated Model

Arachidonic acid (AA) undergoes oxidative metabolism under enzymatic or nonenzymatic mechanisms. Cyclooxygenase (COX), lipoxygenase (LOX), and cytochrome p450s (CYPs) mediated enzymatic peroxidation of AA to the respective metabolites. The 12-lipoxygenase (12-LO) and 15-lipoxygenase (15-LO) mediated conversion of AA to 12- and 15-hydroperoxyeicosatetraenoic acid (HpETE) and the downstream glutathione peroxidase- (GPx-) mediated conversion to the respective hydroxyeicosatetraenoic acids (12- and 15-HETE). ROS-induced lipid peroxidation (LPO) of 12- and 15-HpETE results in the generation of 4-hydroxynonenal (4-HNE) and 4-hydroxydodecadeinal (4-HDDE). These reactive aldehydes interact with and inactivate GPx, leading to an increased rate of 12-HpETE and 15-HpETE peroxidation [72]. The nonenzymatic peroxidation mediates conversion of AA to 4-HNE and MDA metabolites [73] (Figure 6).Figure 6
Schematic representation of arachidonic acid metabolism and lipid peroxidation. Cyclooxygenase (COX), lipoxygenase (LOX), and cytochrome p450s (CYPs) mediated enzymatic peroxidation of arachidonic acid (AA) to the respective metabolites. The 12-lipoxygenase (12-LO) and 15-lipoxygenase (15-LO) mediated conversion of AA to 12- and 15-hydroperoxyeicosatetraenoic acid (HpETE) and the downstream glutathione peroxidase- (GPx-) mediated conversion to the respective hydroxyeicosatetraenoic acids (12- and 15-HETE). Reactive oxygen species- (ROS-) induced lipid peroxidation (LPO) of 12- and 15-HpETE results in the generation of 4-hydroxynonenal (4-HNE) and 4-hydroxydodecadeinal (4-HDDE). These reactive aldehydes interact with and inactivate GPx, leading to an increased rate of 12-HpETE and 15-HpETE peroxidation (PMID: 29610056). The nonenzymatic peroxidation mediates conversion of AA to 4-HNE and MDA metabolites (PMID: 24999379).We hypothesized that downregulation of lipid peroxidation antioxidant enzymes CYPs and GPXs is inversely associated with upregulation of cytokines and chemokines in NAFLD and that upregulation of proinflammatory lipid peroxidation enzymes COXs and LOXs [74] is associated with upregulations of TI enzymes in NAFLD. To examine these hypotheses, we collected 26 COX lipid peroxidation enzymes, 44 cytochrome P450 CYP lipid peroxidation antioxidant (anti-inflammatory) enzymes [75, 76], 26 AA metabolism enzymes, and 8 hepatic glutathione peroxidase antioxidant (anti-inflammatory) enzymes [77].The MCD+HFD upregulated 2 (7.7%) and downregulated 9 (34.6%) COXs (Table3(a)). The MAT1A-KO upregulated 2 (7.7%) and downregulated 1 (3.8%) out of 26 COXs. The GNMT-KO upregulated 3 (11.5%) COXs. The HFCD upregulated 9 (34.6%) and downregulated 1 (3.8%) COXs.Table 3
(a) Cyclooxygenase (COX) lipid peroxidation enzymes were differentially expressed in mouse NASH/NAFLD models. The MCD+HFD model upregulated 2 out of 26 (7.7%) and downregulated 9 out of 26 (34.6%) COXs. The MAT1A-KO model upregulated 2 out of 26 (7.7%) and downregulated 1 out of 26 (3.8%) COXs. The GNMT-KO model upregulated 3 out of 26 (11.5%) and downregulated 0 out of 26 (0%) COXs. The HFCD model upregulated 9 out of 26 (34.6%) and downregulated 1 out of 26 (3.8%) COXs. (b) Hepatocyte glutathione peroxidase lipid antioxidant enzymes were differentially expressed in mouse NASH/NAFLD models. The MCD+HFD model upregulated 4 out of 8 (50%) and downregulated 0 out of 8 (0%) GPXs. The MAT1A-KO model upregulated 2 out of 8 (25%) and downregulated 0 out of 8 (0%) GPXs. The GNMT-KO model upregulated 1 out of 8 (12.5%) and downregulated 0 out of 8 (0%) GPXs. The HFCD model upregulated 0 out of 8 (0%) and downregulated 1 out of 8 (12.5%) GPXs. (c) Cytochrome P450 CYP lipid peroxidation antioxidant (anti-inflammatory) enzymes were differentially expressed in mouse NASH/NAFLD models. The MCD+HFD model upregulated 6 out of 44 (13.6%) and downregulated 13 out of 44 (29.5%) CYP lipid peroxidation antioxidant enzymes (CYP-LPAEs). The MAT1A-KO model upregulated 4 out of 44 (9.1%) and downregulated 8 out of 44 (18.2%) CYP-LPAEs. The GNMT-KO model upregulated 5 out of 44 (11.4%) and downregulated 8 out of 44 (18.2%) CYP-LPAEs. The HFCD model upregulated 22 out of 44 (50%) and downregulated 1 out of 44 (2.3%) CYP-LPAEs. (d) Arachidonic acid metabolism enzymes were differentially expressed in mouse NASH/NAFLD models. The MCD+HFD model upregulated 4 out of 26 (15.4%) and downregulated 11 out of 26 (42.3%) arachidonic acid metabolism enzymes. The MAT1A-KO model upregulated 4 out of 26 (15.4%) and downregulated 6 out of 26 (23.1%) arachidonic acid metabolism enzymes. The GNMT-KO model upregulated 4 out of 26 (15.4%) and downregulated 7 out of 26 (26.9%) arachidonic acid metabolism enzymes. The HFCD model upregulated 24 out of 26 (92.3%) and downregulated 0 out of 26 (0%) arachidonic acid metabolism enzymes.
(a)
COX lipid peroxidation enzymes
GEO IDGSE35961GSE63027GSE63027GSE53381ComparisonMCD+HFD vs. NCDMAT1A-KO vs. WTGNMT-KO vs. WTHFCD vs. NCDGene symbolCOX101.391COX111.388COX141.303COX151.372COX171.3021.364COX18-1.3681.477COX191.438COX4I1-1.314COX5A-1.584COX5B-1.5241.344COX6A1-1.5981.387COX6A21.378-2.695COX6B1-1.437COX6B22.1351.9231.409COX6C-1.334COX7A1-1.7301.439COX7A2-1.313COX7A2LCOX7C1.356COX8C-1.761% up2/26 (7.7%)2/26 (7.7%)3/26 (11.5%)9/26 (34.6%)% down9/26 (34.6%)1/26 (3.8%)0/26 (0%)1/26 (3.8%)(b)
Hepatocyte glutathione peroxidase lipid antioxidant enzyme
GEO IDGSE35961GSE63027GSE63027GSE53381ComparisonMCD+HFD vs. NCDMAT1A-KO vs. WTGNMT-KO vs. WTHFCD vs. NCDGene symbolGPX38.8572.387GPX41.8321.590GPX73.7632.196GPX82.266-1.693% up4/8 (50%)2/8 (25%)1/8 (12.5%)0/8 (0%)% down0/8 (0%)0/8 (0%)0/8 (0%)1/8 (12.5%)(c)
Cytochrome P450 CYP lipid peroxidation antioxidant (anti-inflammatory) enzymes
GEO IDGSE35961GSE63027GSE63027GSE53381ComparisonMCD+HFD vs. NCDMAT1A-KO vs. WTGNMT-KO vs. WTHFCD vs. NCDGene symbolCYP11A1-2.023CYP11B2-3.851-1.459CYP17A1-2.7817.210CYP1B13.0801.496CYP20A12.293CYP21A2CYP2B10-3.16012.086CYP2B139.4553.6487.655CYP2B948.3036.8114.811CYP2C29-4.350-2.577-1.7361.663CYP2C37-2.233-3.070-2.3951.834CYP2C38-1.531-2.7842.6942.499CYP2C392.7711.751CYP2C401.830CYP2C54-113.732-1.679-3.6811.907CYP2C551.681-2.0854.526CYP2C651.413CYP2F1CYP2J111.3471.621CYP2J2CYP2J5-4.063-1.5691.442CYP2R1-2.226-2.0321.441CYP2S1-1.756CYP2U1-6.007-1.9901.473CYP2W1CYP3A7CYP4A102.235CYP4A11CYP4A12A-7.4363.6641.6702.931CYP4A12B2.789CYP4A143.575-13.39216.11114.158CYP4A31-2.77313.473CYP4B1-1.598-1.619CYP4F13-2.8111.525CYP4F18% up6/44 (13.6%)4/44 (9.1%)5/44 (11.4%)22/44 (50%)% down13/44 (29.5%)8/44 (18.2%)8/44 (18.2%)1/44 (2.3%)(d)
Arachidonic acid metabolism enzymes
GEO IDGSE35961GSE63027GSE63027GSE53381MCD+HFD vs. NCDMAT1A-KO vs. WTGNMT-KO vs. WTHFCD vs. NCDGene symbolCYP2C37-2.233-3.070-2.3951.834CYP2C29-4.350-2.577-1.7361.663CYP2C54-113.732-1.679-3.6811.907CYP2J5-4.063-1.5691.442CYP4F13-2.8111.525CYP2U1-6.007-1.9901.473PLA2G4C-3.499LTA4H-2.953CYP 20-1.8181.432CYP4A12A-7.4363.6641.6702.931CYP2C38-1.531-2.7842.6942.499CYP2C651.413PLA2G61.429ALOX12B1.439PTGES21.508CYP2J111.3471.621CYP2C401.830CYP4A102.235CYP4A12B2.789CYP2C551.681-2.0854.526CYP2B10-3.16012.086CYP4A31-2.77313.473CYP4F18CYP2C392.7711.751CYP4A143.575-13.39216.11114.158CYP2B139.4553.6487.655CYP2B948.3036.8114.811% up4/26 (15.4%)4/26 (15.4%)4/26 (15.4%)24/26 (92.3%)% down11/26 (42.3%)6/26 (23.1%)7/26 (26.9%)0/26 (0%)The MCD+HFD upregulated 6 (13.6%) and downregulated 13 (29.5%) CYP lipid peroxidation antioxidant enzymes (CYP-LPAEs) (Table3(c)). The MAT1A-KO upregulated 4 (9.1%) and downregulated 8 (18.2%) CYP-LPAEs. The GNMT-KO upregulated 5 (11.4%) and downregulated 8 (18.2%) CYP-LPAEs. The HFCD upregulated 22 (50%) and downregulated 1 (2.3%) CYP-LPAEs (Figure 6(c)).The MCD+HFD upregulated 4 (15.4%) and downregulated 11 (42.3%) AA metabolism enzyme genes (Table3(d)). The MAT1A-KO upregulated 4 (15.4%) and downregulated 6 (23.1%) genes. The GNMT-KO upregulated 4 (15.4%) and downregulated 7 (26.9%) genes. The HFCD upregulated 24 (92.3%) genes (Figure 6(d)).The MCD+HFD upregulated 4 out of 8 (50%) GPXs (Table3(b)). The MAT1A-KO upregulated 2 (25%) GPXs. The GNMT-KO upregulated 1 out of 8 (12.5%) GPXs. The HFCD downregulated 1 out of 8 (12.5%) GPXs (Figure 6(e)).Taken together, our results have demonstrated that,first, proinflammatory MCD+HFD upregulates 7.7% COXs, 13.6% CYP-LPAEs, 15.4% AA metabolism enzymes, and 50% GPXs and downregulates 34.6% COXs, 29.5% CYP-LPAEs, and 42.3% AA metabolism enzymes, which are in big contrast with that of less inflammatory HFCD with upregulations of 34.6% COXs, 50% CYP-LPAEs, and 92.3% AA metabolism enzymes and downregulation of 3.8% COXs, 2.3% CYP-LPAEs, and 12.5% GPXs; second, significant upregulation of lipid peroxidation enzymes in HFCD is positively correlated with significant upregulation of TI enzymes and is, however, inversely associated with low-level upregulation of canonical inflammasome regulators (6.7%), noncanonical inflammasome regulators (7.1%), and 2.4% secretomic cytokines; third, no upregulation of GPXs in HFCD (Table 3(b)) is well correlated with other’s findings that GPX4 decreases lipid peroxidation and inflammation, and caspase-11-dependent pyroptosis mediates atherosclerosis [78], endotoxic shock [79], and septic death in GPX4-KO mice [33]; and fourth, the MCD+HFD and HFCD are the two representative models: the MCD+HFD is a model of NASH/NAFLD with upregulation of proinflammatory cytokines and canonical and noncanonical inflammasome regulators. In contrast, the HFCD is a model with upregulation of lipid peroxidation enzymes and TI enzymes.


### 3.6. As Upstream Master Regulators, Caspase-11 and Caspase-1 Partially Upregulate the Expressions of Cytokines, Chemokines, Canonical and Noncanonical Inflammasome Pathway Regulators, TI Enzymes, and Lipid Peroxidation Enzymes

It has been reported that caspase-11 mediates hepatocyte pyroptosis and promotes the progression of MCD diet-induced NASH/NAFLD [80]; and inflammasome-gasdermin D-pyroptosis [81] may regulate the progression of NASH/NAFLD [82]. In addition to promoting inflammation by enzymatic cleavage of pre-IL-1β, pre-IL-18, and N-terminal GSDMD, caspase-1 also cleaves as many as 114 substrates, suggesting that caspase-1 may regulate inflammation indirectly via all the protein substrate-mediated signaling [15], 38 interaction protein-mediated signaling [15], and GSDMD-secretome-mediated signaling [23]. However, an important question remains whether caspase-11 modulates the expression of cytokines, chemokines, and TI enzymes. We hypothesized that caspase-11 promotes inflammatory cytokines and chemokines upregulated in NASH/NAFLD and modulates the expression of canonical and noncanonical inflammasome pathway genes and TI enzymes. To examine this hypothesis, we collected all the cytokines and chemokines, canonical and noncanonical inflammasome pathway genes, and TI enzymes upregulated in NASH/NAFLD and crossed with caspase-11-KO and caspase-1-KO datasets (Supplementary Table 2b).Among 39 innatomic and secretomic cytokines and chemokines upregulated in NASH/NAFLD, caspase-11 deficiency upregulated 18 (46.2%) and downregulated 8 (20.5%) NASH/NAFLD-upregulated cytokines and chemokines. In comparison, caspase-1 deficiency upregulated one cytokine CXCL14 (2.6%) and downregulated one cytokine LTB (2.6%) NASH/NAFLD-upregulated cytokines and chemokines (Table4(a)). These results suggest that capsase-11 deficiency promotes the upregulation of most cytokines and chemokines induced by NASH/NAFLD. In addition, caspase-11 deficiency upregulated 5 (17.9%) and downregulated 5 (17.9%) canonical inflammasome regulators induced by NASH/NAFLD. However, caspase-1 deficiency did not modulate the expression of canonical inflammasome regulators induced by NASH/NAFLD (Table 4(b)). Furthermore, caspase-11 deficiency upregulated 6 (8.8%) and downregulated 31 (45.6%) TI enzymes upregulated in NASH/NAFLD. In comparison, caspase-1 deficiency upregulated 4 (5.9%) and downregulated 6 (8.8%) TI enzymes upregulated in NASH/NAFLD (Table 4(c)). Finally, caspase-11 deficiency upregulated 13 (28.3%) and downregulated 6 (13%) lipid peroxidation enzymes upregulated in NASH/NAFLD. In comparison, caspase-1 deficiency upregulated 3 (6.5%) and downregulated 6 (13%) lipid peroxidation enzymes upregulated in NASH/NAFLD (Table 4(d)).Table 4
(a) The cytokines and chemokines upregulated in NASH/NAFLD were differentially expressed in caspase-11 and caspase-1 knockout. Among 28 cytokines and chemokines upregulated in NASH/NAFLD, caspase-11 deficiency upregulated 18 out of 28 (64.3%) and downregulated 8 out of 28 (28.6%) NASH/NAFLD-upregulated cytokines and chemokines. Caspase-1 deficiency upregulated one cytokine CXCL14 out of 39 (2.6%) and downregulated one cytokine LTB out of 39 (2.6%) NASH/NAFLD-upregulated cytokines and chemokines. (b) The upregulated canonical inflammasome regulators in NASH/NAFLD were differentially expressed in caspase-11 knockout. Caspase-11 deficiency upregulated 5 out of 28 (17.9%) and downregulated five out of 28 (17.9%) canonical inflammasome regulators induced by NASH/NAFLD. Caspase-1 deficiency did not modulate the expression of canonical inflammasome regulators induced by NASH/NAFLD. (c) The upregulated trained immunity enzyme genes in NASH/NAFLD were differentially expressed in caspase-11 knockout and caspase-1 knockout mouse models. Caspase-11 deficiency upregulated 6 out of 68 (8.8%) and downregulated 31 out of 68 (45.6%) trained immunity enzymes upregulated in NASH/NAFLD. Caspase-1 deficiency upregulated 4 out of 68 (5.9%) and downregulated 6 out of 68 (8.8%) trained immunity enzymes upregulated in NASH/NAFLD. (d) The upregulated lipid peroxidation enzyme genes in NASH/NAFLD were differentially expressed in caspase-11 knockout and caspase-1 knockout mouse models. Caspase-11 deficiency upregulated 13 out of 46 (28.3%) and downregulated 6 out of 46 (13%) lipid peroxidation enzymes upregulated in NASH/NAFLD. Caspase-1 deficiency upregulated 3 out of 46 (6.5%) and downregulated 6 out of 46 (13%) lipid peroxidation enzymes upregulated in NASH/NAFLD.
(a)GEO IDGSE115094GSE32515ComparisonCaspase-11-KO vs. WTCaspase-1-KO vs. WTGene symbolCXCL101.98CCL205.07EDN13.45CXCL162.50CKLF2.83CCL54.81TIMP110.52CCL24.75Ccl75.99TNF2.72TSLP7.42CCL82.70IL274.15PF44.82CCL1713.15CCL33.58CCL47.09CXCL54.50CXCL143.219DKK 3-5.54CX3CL1-4.49CXCL12-3.50CSF1-3.04TNFSF10-2.23NAMPT-2.08WNT5A-1.79SPP1-5.87LTB-1.551% up18/39 (46.2%)1/39 (2.6%)% down8/39 (20.5)1/39 (2.6%)(b)Canonical inflammasome pathway regulators upregulated in NASHGEO IDGSE115094GSE32515ComparisonCaspase-11-KO vs. WTCaspase-1-KO vs. WTGene symbolTXN21.749PYCARD2.469Naip51.602Naip21.895CYBA3.599PRKCD-1.643ANTXR2-3.523CTSB-2.328NAMPT-2.083RIPK1-1.754% up5/28 (17.9%)0/28 (0%)% down5/28 (17.9%)0/28 (0%)(c)Trained immunityGEO IDGSE115094GSE32515PathwaysComparisonCaspase-11-KO vs. WTCaspase-1-KO vs. WTGlycolysisLDHB2.520PFKFB31.844ALDOA1.883FBP14.254ALDOA1.883ADH71.624PFKFB11.511PFKFB21.606HK1-2.535HK2-4.554HKDC1-1.550PDHA1-3.277PFKFB4-1.570PFKP-1.544PGAM1-2.484ENO1-1.689DLAT-1.770GALM-1.700PANK1-1.892PFKL-4.133PGM2-1.663PKLR-3.863-2.354ALDH1B1-2.881ALDH2-2.171ALDH3A2-1.840DLAT-1.770ACSS2-1.836ADH4-2.019GCK-3.334SLC2A2-1.789ACSS1-3.206ALDH1B1-2.881ALDH2-2.171ALDH3A2-1.840MevalonatePMVK1.699IDI1-3.248HMGCS1-1.818HMGCR-1.894ACAT1-1.917IDI1-3.248Acetyl-CoA generationGOT1-1.6841.701BDH1-1.964ACO1-1.601ACLY-2.203-2.33% up6/68 (8.8%)4/68 (5.9%)% down31/68 (45.6%)6/68 (8.8%)(d)Lipid peroxidationGEO IDGSE115094GSE32515EnzymesComparisonCaspase-11-KO vs. WTCaspase-1-KO vs. WTCOXsCOX145.806COX181.862COX5B6.309COX6A15.582COX6B212.677COX7A111.346COX7C11.363COX194.005COX15-3.650CYPsCYP2R14.459CYP2B107.4684.070CYP2C3915.441CYP4A311.831CYP2U1-2.833CYP1B1-2.058CYP20A1-2.466CYP2B9-3.075CYP2C37-6.050CYP2C38-3.200CYP2C54-2.609CYP2C55-6.187CYP4F18-1.575GPXGPX46.696GPX73.769GPX83.769GPX3-1.549AA metabolismPLA2G6-2.293% up13/46 (28.3%)3/46 (6.5%)% down6/46 (13%)6/46 (13%)Taken together, our results have demonstrated that,first, caspase-11 and caspase-1 partially modulate the expressions of cytokines, chemokines, canonical and noncanonical inflammasome pathway regulators, TI enzymes, and lipid peroxidation enzymes; second, in the context of NASH/NAFLD, caspase-11 has stronger capacities than caspase-1 in modulating inflammatory gene expressions; and third, experimental conditions in caspase-11 deficiency [79] and caspase-1 deficiency [83] were not exactly the same as NASH/NAFLD. Since we used inflammatory gene lists upregulated in the four NASH/NAFLD models for the analyses, our results are at least partially relevant to NASH/NAFLD pathology. Future work on transcriptomic analysis of caspase-11 deficiency and caspase-1 deficiency in NASH/NAFLD models, respectively, is needed to verify these findings presented here.


## 3.1. Human NASH and NAFLD Mouse Models Upregulate Both Cytokines and Chemokines Classified in the Innate Immune Database and Canonical Secretome

At least 11 NAFLD mouse models have been established [45, 46] including diet-induced and genetic-induced models [47, 48], but no model accurately reflects the disease’s progression in humans [3]. We include four mouse models that best represent the metabolic and inflammatory characteristics of human NAFLD [3, 46–48].We hypothesized that diet- and genetically induced NAFLD mouse models can share common pathways critical to the development of the disease. To test this, we collected microarray datasets from NAFLD mouse models and human NASH (Supplementary Table2a). The datasets we analyzed include (1) the high-fat diet (HFD) and methionine- and choline-deficient (MCD) diet model (GSE35961) [49], (2) the high-fat-cholesterol diet- (HFCD-) induced diet model [50], (3) the glycine N-methyltransferase- (GNMT-) KO genetic model, (4) the liver-specific methionine adenosyltransferase 1A- (MAT1A-) KO (GSE63027) [51] genetic model, and (5) two human NASH microarrays (GSE63067 [51] and GSE17470 [52]). The qualities of the datasets were assessed using a few housekeeping genes (Supplementary Tables 1a and 1b).Cytokines and chemokines [53–55] play significant roles in promoting the onset and progression of NAFLD [55]; however, their expression changes in human NASH and NAFLD mouse models remained poorly characterized. Cytokines and chemokines in the liver are soluble mediators of local and systemic inflammation (Figure 2(a)) [8, 43, 56, 57]. To analyze liver inflammation, we used IPA and classified 53 cytokines and chemokines from 1376 innate immune genes [58, 59] from the Innate Immunity Database [58, 60]. Human NASH upregulated six (11.3%) and four (7.5%) out of 53 innate immune cytokines and chemokines, respectively (Figure 2(b)). Of note, upregulations of C-X-C motif chemokine ligand 10 (CXCL10) and CKLF-like MARVEL transmembrane domain-containing 2 (CMTM2) were shared in these two human NASH datasets. MCD+HFD and MAT1A-KO upregulated 13 (24.5%) and 6 (11.3%) out of 53 innate immune cytokines and chemokines, respectively (Figure 2(c)). MCD+HFD, MAT1A-KO, GNMT-KO, and HFCD downregulated 9 (17%), 3 (5.7%), 1 (1.9%), and 12 (22.6%) out of 53 innate immune cytokines and chemokines, respectively.Figure 2
The cytokine and chemokine expressions are increased in both human nonalcoholic steatohepatitis (NASH) and mouse models of nonalcoholic fatty liver disease (NAFLD). Ingenuity Pathway Analysis results showed that 53 cytokines and chemokines were sorted out from 1376 innate immune genes (innatome, from the Innate Immune Database (https://www.innatedb.com/); also, see our paper (PMID: 33628851)) and 123 cytokines and chemokines were sorted out from 2641 canonical secretome (with signal peptide) genes (from the Human Protein Atlas database (https://www.proteinatlas.org/); also, see our paper (PMID: 32179051)). (a) A schematic figure showed the pathogenic effects of cytokines released from NASH on different organs. Two human NASH datasets (GSE63067 and GSE17470) and four mouse models of NAFLD datasets (GSE35961, GSE63027, GSE63027, and GSE53381) were analyzed. (b, c) The 53 innatome cytokines and chemokines were analyzed in human NASH and mouse models of NAFLD. (d) Venn diagram showed the significant overlapped regulated innatome cytokines and chemokines in human NASH and mouse models of NAFLD. The cytokine and chemokine expressions are increased in both human nonalcoholic steatohepatitis (NASH) and mouse models of nonalcoholic fatty liver disease (NAFLD). (e, f) The 123 canonical secretomic cytokines and chemokines were analyzed in human NASH and mouse models of NAFLD. (g) Venn diagram showed the overlapped significant regulated canonical secretomic cytokines and chemokines in human NASH and mouse models of NAFLD. The cytokine and chemokine expressions are increased in both human nonalcoholic steatohepatitis (NASH) and mouse models of nonalcoholic fatty liver disease (NAFLD). (h) The 16 cytokines and chemokines (innatome and secretome) were upregulated (upregulated in GSE63067 or GSE17470) in human NASH. PMID: 34084175. (i) The 16 cytokines and chemokines (innatome and secretome) were upregulated (GSE35961 or GSE63027) in mouse models of NAFLD. Secretome gene list detailed in our previous publication PMID: 34084175.
(a)(b)(c)(d)(e)(f)(g)(h)(i)To determine whether NAFLD mouse models share innate immune cytokines and chemokines with human NASH, we performed Venn diagram analysis. Six innate immune cytokines and chemokines including CXCL10, tumor necrosis factor (ligand) superfamily, member 13b (TNFSF13B), lymphotoxin beta (LTB, TNFSF3), TNFSF10, interleukin 1 receptor antagonist (IL-1RN), and nicotinamide phosphoribosyltransferase (NAMPT) were shared between human NASH and NAFLD mouse models (Figure2(d)). Two NAFLD mouse models upregulated 22 (78.6%) specific innate immune cytokines and chemokines; and human NASH upregulated two (25%) specific cytokines and chemokines C-C motif chemokine ligands 20 (CCL20) and CMTM2.To ensure a comprehensive analysis of liver inflammation in human NASH and NAFLD mouse models, we collect 2641 canonical secretome genes [61]. The 123 cytokines and chemokines were classified with IPA out of 2641 human canonical secretome genes. Human NASH upregulated 4 (3.3%) and 11 (8.9%) out of 123 canonical secretomic cytokines and chemokines (Figure 2(e)). Of note, upregulations of CXCL10 and thymic stromal lymphopoietin (TSLP) were shared in the two human NASH datasets. MCD+HFD, MAT1A-KO, GNMT-KO, and HFCD models upregulated 21 (17.1%), 7 (5.7%), 0, and 3 (2.4%) (Figure 2(f)) and downregulated 10 (8.1%), 4 (3.3%), 1 (0.8%), and 10 (8.1%) out of 123 canonical secretomic cytokines and chemokines, respectively.To determine whether NAFLD mouse models share canonical secretomic cytokines and chemokines with human NASH, we performed Venn diagram analysis. Seven canonical secretomic cytokines and chemokines including TSLP, CXCL10, CCL28, CCL25, CXCL14, IL21, and IL1RN were shared between NAFLD mouse models and human NASH (Figure2(g)). By comparison, three NAFLD mouse models upregulated 32 (82.1%) specific canonical secretomic cytokines and chemokines; and human NASH upregulated 8 (53.3%) specific cytokines and chemokines including erythropoietin (EPO), IL17F, CCL8, IL27, CCL21, CXCL12, CCL15, and CCL20.Caspase-1-deficient mice are protected from high fat-induced hepatic steatosis, inflammation, and early fibrogenesis [62], suggesting that pyroptosis-related caspase-1/inflammasomes [6] promote NASH progression. However, an important question remains whether noncanonical inflammasome pathway-related caspase-4/11 [63] plays any role in progression of NAFLD.To determine whether NASH/NAFLD upregulated cytokines and chemokines classified into noncanonical secretomes, we collected four types of secretomes with total 11,385 proteins including canonical secretome (2641 proteins) [61], caspase-1-dependent noncanonical secretome (961 proteins) [23], caspase-4-dependent noncanonical secretome (1223 proteins) [24], and exosome secretome (6560 proteins) [25]. We combined human NASH and NAFLD mouse upregulated cytokines and chemokines and performed Venn diagram analysis. Among 15 human NASH-upregulated cytokines and chemokines, two cytokines CMTM2 and LTB were nonclassified, 12 out of 15 (80%) were classified in the canonical secretome, TNFSF10 was classified in canonical and exosome secretomes, and NAMPT was classified in canonical and caspase-4 secretomes (Figure 2(h)). Of note, since 123 cytokines and chemokines identified with the IPA are derived from the canonical secretome, it was expected that the high percentage (80%) NASH-upregulated cytokines and chemokines were from the canonical secretome. Among 26 NAFLD mouse-upregulated cytokines and chemokines, one was the nonclassified, 15 were classified in the canonical secretome, seven were classified in canonical and exosome secretomes, NAMPT was classified in canonical and caspase-4 secretomes, and two cytokines were classified in canonical, exosomes, and caspase-4 secretomes (Figure 2(i)).Taken together, our results have demonstrated innate immune secretomic cytokines and chemokines suggesting that human NASH and NAFLD mouse models share innate immune mechanisms and there is a significant role of exosomes and caspase-4 secretomic in driving liver and systemic inflammations.


## 3.2. Human NASH Can Be Classified into Metabolic and Immune NASH; MCD+HFD, GNMT-KO, MAT1A-KO, and HFCD Can Be Classified into Inflammatory, SAM Accumulation, Cholesterol/Mevalonate, and LXR/RXR-Fatty Acidβ-Oxidation NAFLD, Respectively

To determine the signaling pathways mediating the transcriptomic changes in human NASH and NAFLD mouse models, we performed IPA analysis for the top 20 pathways for significantly modulated genes in each dataset. Human NASH dataset GSE17470 had top upregulated pathways including phospholipases, liver X receptor (LXR)/retinoid X receptor (RXR) activation, fatty acidβ-oxidation I, and leukocyte extravasation, among which leukocyte extravasation was shared by the second human NASH and MCD+HFD, LXR/RXR activation was shared by HFCD, and fatty acid β-oxidation was shared by GNMT-KO and HFCD (Table 1, Supplementary Figure 1).Table 1
Ingenuity Pathway Analysis (IPA) for top 20 pathways of significantly modulated genes in each dataset of human NASH studies and mouse models of NAFLD. Leukocyte extravasation was shared by both human NASH studies and in the MCD+HFD mouse model; LXR/RXR activation was shared by HFCD model of NAFLD; and fatty acidβ-oxidation was shared by the GNMT-KO and HFCD models of NAFLD. IL-8 signaling was shared by the second human NASH study and the MCD+HFD model. Fatty acid β-oxidation I and stearate biosynthesis I were shared by the first human NASH study and GNMT-KO and HFCD models. Nicotine degradation II was shared by the first human NASH study and the MAT1A-KO model.
(a)HumanMouseGEO IDGSE17470GSE63067GSE35961GSE63027GSE63027GSE53381ComparisonNASH vs. healthyNASH vs. healthyMCD+HFDGNMT-KOMAT1A-KOHFCDPathwaysSerotonin degradation↑↑Ethanol degradation II↑Synaptic long-term depression↑Noradrenaline and adrenaline degradation↑Sperm motility↑Phospholipases↑LXR/RXR activation↑↑Fatty acidβ-oxidation I↑↑↑Nicotine degradation II↑↑↑eNOS signaling↑Leukocyte extravasation signaling↑↑↑p70S6K signaling↑Retinol biosynthesis↑Stearate biosynthesis I (animals)↑↑Histamine degradation↑Aldosterone signaling in epithelial cells↑HIPPO signaling↑cAMP-mediated signaling↑Synaptic long-term potentiation↑Fcγ receptor-mediated phagocytosis in macrophages and monocytes↑↑IL-8 signaling↑↑Rac signaling↑NF-κB activation by viruses↑Salvage pathways of pyrimidine ribonucleotides↑Gαq signaling↑Pyridoxal 5′-phosphate salvage pathway↑GP6 signaling pathway↑Colorectal cancer metastasis signaling↑ILK signaling↑Integrin signaling↑mTOR signaling↑Phospholipase C signaling↑Cardiac hypertrophy signaling (enhanced)↑↑ERK5 signaling↑Sphingosine-1-phosphate signaling↑Small cell lung cancer signaling↑Gα12/13 signaling↑Signaling by rho family GTPases↑Superpathway of cholesterol biosynthesis↑PPAR signalingNicotine degradation III↑Melatonin degradation ISuperpathway of melatonin degradationXenobiotic metabolism CAR signaling pathwayXenobiotic metabolism PXR signaling pathwayCholesterol biosynthesis I↑Cholesterol biosynthesis II (via 24,25-dihydrolanosterol)Cholesterol biosynthesis III (via desmosterol)↑(b)HumanMouseGEO IDGSE17470GSE63067GSE35961GSE63027GSE63027GSE53381ComparisonNASH vs. healthyNASH vs. healthyMCD+HFD vs. NCDGNMT-KO vs. WTMAT1A-KO vs. WTHFCD vs. NCDPathwaysSuperpathway of geranylgeranyl diphosphate biosynthesis I (via mevalonate)Estrogen biosynthesis↑Acetone degradation I (to methylglyoxal)↑Mevalonate pathway I↑Tryptophan degradation III (eukaryotic)Antioxidant action of vitamin C↑Oxidative phosphorylation↑Glutaryl-CoA degradation↑Valine degradation I↑Isoleucine degradation I↑RhoA signaling↑Mitotic roles of polo-like kinase↑Triacylglycerol biosynthesis↑TCA cycle II (eukaryotic)↑Role of CHK proteins in cell cycle checkpoint control↑Sumoylation pathway↑Apelin adipocyte signaling pathway↑↑Glutathione-mediated detoxification↑Tetrahydrofolate salvage from 5,10-methenyltetrahydrofolate↑CMP-N-acetylneuraminate biosynthesis I (eukaryotes)↑EIF2 signaling↑tRNA charging↑Folate transformations I↑Endocannabinoid cancer inhibition pathway↑Pyrimidine ribonucleotide interconversion↑Pyrimidine ribonucleotide de novo biosynthesis↑NAD salvage pathway II↑PPARα/RXRα activation↑NRF2-mediated oxidative stress response↑↑Glutathione redox reactions I↑Pancreatic adenocarcinoma signaling↑Estrogen-mediated S-phase entry↑Type II diabetes mellitus signaling↑PCP pathway↑PDGF signaling↑NF-κB signaling↑Neuroinflammation signaling pathway↑Tec kinase signaling↑B cell receptor signaling↑Osteoarthritis pathway↑Th17 activation pathway↑Role of pattern recognition receptors in recognition of bacteria and viruses↑Apelin endothelial signaling pathway↑Cholecystokinin/gastrin-mediated signaling↑T cell exhaustion signaling pathway↑Endothelin-1 signaling↑B cell activating factor signaling↑RANK signaling in osteoclasts↑Toll-like receptor signaling↑Renin-angiotensin signaling↑The GNMT-KO (GSE63027) had top upregulated pathways including fatty acidβ-oxidation I, stearate biosynthesis I, antioxidant action of vitamin C, oxidative phosphorylation, glutaryl-CoA degradation kinase, triacylglycerol biosynthesis, and tricarboxylic acid cycle (TCA) cycle II, among which fatty acid β-oxidation I and stearate biosynthesis I were shared with the first human NASH (GSE17470). The MAT1A-KO (GSE63027) had top upregulated pathways including nicotine degradation II, superpathway of cholesterol biosynthesis, cholesterol biosynthesis I, cholesterol biosynthesis III, mevalonate pathway I, and NAD salvage pathway II, among which nicotine degradation II was shared with the first human NASH.The HFCD model (GSE53381) has only four upregulated pathways including LXR/RXR activation and fatty acidβ-oxidation I, among which LXR/RXR activation was shared with the first human NASH, fatty acid β-oxidation II was shared with the first human NASH and GNMT-KO, and nicotine degradation II was shared with MAT1A-KO.The second human NASH (GSE63067) had top upregulated pathways including NF-κB signaling and neuroinflammation (Table 1). The MCD+HFD model (GSE35961) had top upregulated pathways including leukocyte extravasation, antibody Fc fragment gamma (Fcγ) receptor-mediated phagocytosis, IL-8 signaling, Rac signaling, and glucose 6-phosphate (G6P) signaling, among which IL-8 signaling was shared with the second human NASH.Taken together, our results have demonstrated that the two human NASH datasets have diversified signaling pathways, which allow us to classify the first one into LRX/RXR activation-related metabolic NASH and the second NASH into NF-κB signaling [64]/T helper cell 17 (Th17) [56, 65] inflammatory NASH. IPA of upregulated genes in each mouse models of NASH and human NASH identified top pathways for significantly modulated genes. The major differences between metabolic NASH and immune NASH include the following: first, numerous inflammation and immune pathways are identified on the top pathway list, and second, metabolic reprogramming in any given mouse models affects four metabolic pathways including (1) increased glycolysis, (ii) glutaminolysis, (iii) increased accumulation of tricarboxylic acid cycle metabolites and acetyl-coenzyme A production, and (iv) increased mevalonate synthesis. Furthermore, this classification provides a criterion for mouse model stratification and selection, providing novel insight of which mouse models best represent various human cases of NASH/NAFLD.


## 3.3. Canonical and Noncanonical Inflammasome Pathways Play Differential Roles in the Pathogenesis of NASH/NAFLD

To examine the differential roles of caspase-1- and caspase-11-dependent pyroptosis, we collected 90 canonical and 14 noncanonical inflammasome pathway regulator genes and determined the expression changes of these inflammasome regulators in four NAFLD mouse models.MCD+HFD upregulated 18 (20%) including caspase-1 and downregulated 11 (12.2%) out of 90 canonical inflammasome pathway regulators. MAT1A-KO upregulated 5 (5.6%) and downregulated 4 (4.4%) out of 90 canonical inflammasome pathway regulators. GNMT-KO upregulated 1 (1.1%) and downregulated 2 (2.2%) out of 90 canonical inflammasome pathway regulators. HFCD upregulated 7 (7.8%) and downregulated 26 (28.9%) out of 90 canonical inflammasome pathway regulators (Figure3(b)). In addition, MCD+HFD downregulated 6 out of 14 (42.8%) noncanonical inflammasome pathway regulators. MAT1A-KO upregulated 3 (21.4%) and downregulated 3 (21.4%) out of 14 noncanonical inflammasome pathway regulators. GNMT-KO downregulated 1 out of 14 (7.1%) noncanonical inflammasome pathway regulators. HFCD upregulated 1 (7.1%) and downregulated 9 (64.3%) out of 14 noncanonical inflammasome pathway regulators (Figure 3(c)).Figure 3
Mouse pyroptosis pathway regulators. (a) Schematic representation of canonical vs. noncanonical pyroptosis. (b) Mouse canonical inflammasome pathway regulators (90 genes, KEGG pathway hsa04621) were differentially expressed in NASH/NAFLD models. The MCD+HFD model upregulated 18 out of 90 (20%) canonical inflammasome pathway regulators including caspase-1 and downregulated 11 out of 90 (12.2%) canonical inflammasome pathway regulators. The MAT1A-KO model upregulated 5 out of 90 (5.6%) and downregulated 4 out of 90 (4.4%) canonical inflammasome pathway regulators. The GNMT-KO model upregulated 1 out of 90 (1.1%) and downregulated 2 out of 90 (2.2%) canonical inflammasome pathway regulators. The HFCD model upregulated 7 out of 90 (7.8%) and downregulated 26 out of 90 (28.9%) canonical inflammasome pathway regulators. (c) Mouse noncanonical inflammasome pathway regulators (14 genes, KEGG pathway hsa04621) were differentially expressed in NASH/NAFLD. The MCD+HFD model upregulated zero out of 14 (0%) and downregulated 6 out of 14 (42.8%) noncanonical inflammasome pathway regulators. The MAT1A-KO model upregulated 3 out of 14 (21.4%) and downregulated 3 out of 14 (21.4%) noncanonical inflammasome pathway regulators. The GNMT-KO model upregulated 0 out of 14 (0%) and downregulated 1 out of 14 (7.1%) noncanonical inflammasome pathway regulators. The HFCD model upregulated 1 out of 14 (7.1%) and downregulated 9 out of 14 (64.3%) noncanonical inflammasome pathway regulators.
(a)(b)(c)Taken together, these results have demonstrated that the MCD+HFD upregulates more canonical inflammasome regulators than the other three NAFLD models. Conversely, the MAT1A-KO upregulates more noncanonical inflammasome regulators including caspase-4 than other models, suggesting that canonical and noncanonical inflammasome pathways are regulated differentially in various NAFLD models. The HFCD downregulates more than upregulates canonical and noncanonical inflammasome regulators, suggesting that HFCD model uses more focused set of regulators in modulating the expression of these regulators.


## 3.4. Trained Immunity (TI) Enzymes Are Significantly Upregulated in Human NASH and NAFLD Mouse Models; the HFD Upregulates TI Enzymes More Than Cytokines, Chemokines, and Canonical and Noncanonical Inflammasome Regulators; and Statins Promote rather than Suppress TI Enzyme Expressions

We hypothesized that the expressions of TI enzymes are modulated in NASH and NAFLD to establish metabolic reprogramming and achieve innate immune memory [66]. Three metabolic pathways play significant roles in establishing TI (Figure 4). To examine this, we collected 71 glycolysis enzyme genes, 24 acetyl-CoA generation enzyme genes, and 10 mevalonate genes [27, 29, 30, 67, 68].Figure 4
Schematic representation of fundamental mechanisms in trained immunity. Danger-associated molecular patterns (DAMPs) can bind to their corresponding pattern recognition receptors (PRRs) and alter metabolic pathways including increased glycolysis, increased acetyl-CoA generation, and mevalonate synthesis in cholesterol pathway, leading to histone modifications that enable chromatin regions to be more open for transcription. Increased gene expression of proinflammatory cytokines/chemokines and enhanced proinflammatory innate immune response against pathogens during secondary exposure. Figure created with Smart Servier Medical Art and Omnigraffle.The first human NASH upregulated 25 (35.2%) and downregulated 5 (7%) out of 71 glycolysis enzymes genes (Table2(a)). The second human NASH upregulated 5 out of 71 (7%) genes. The MCD+HFD upregulated 20 (28.2%) and downregulated 24 (33.8%) out of 71 genes. The MAT1A-KO upregulated 6 (8.5%) and downregulated 13 (18.3%) out of 71 genes. The GNMT-KO model upregulated 7 (9.9%) and downregulated 4 (5.6%) out of 71 genes. The HFCD upregulated 17 (23.9%) and downregulated 8 (11.3%) out of 71 genes.Table 2
(a) The 61 out of 71 mRNAs (85.9%) of glycolysis pathway enzyme genes are regulated in human and mouse NASH/NAFLD. Glycolysis pathway enzyme gene expression (PMID: 31153039) for human NASH studies (GSE17470 and GSE63067), MCD+HFD (GSE35961), HFC (GSE53381), MAT1A-KO (GSE63027), and GNMT-KO (GSE63027). The first human NASH study (GSE17470) upregulated 25 out of 71 (35.2%) and downregulated 5 out of 71 (7%) glycolysis enzymes. The second human NASH study (GSE63067) upregulated 5 out of 71 (7%) and downregulated 0 out of 71 (0%) glycolysis enzymes. The MCD+HFD model upregulated 20 out of 71 (28.2%) and downregulated 24 out of 71 (33.8%) glycolysis enzymes. The MAT1A-KO model upregulated 6 out of 71 (8.5%) and downregulated 13 out of 71 (18.3%) glycolysis enzymes. The GNMT-KO model upregulated 7 out of 71 (9.9%) and downregulated 4 out of 71 (5.6%) glycolysis enzymes. The HFCD model upregulated 17 out of 71 (23.9%) and downregulated 8 out of 71 (11.3%) glycolysis enzymes. The 71 glycolysis pathway genes are shown in Supplementary Table3a. (b) The 13 out of 24 mRNA (54.2%) of acetyl-CoA generating enzyme are regulated in human and mouse NASH/NAFLD. Acetyl-CoA generating pathway enzyme gene expression (PMID: 31153039) for human NASH studies (GSE17470 and GSE63067), MCD+HFD (GSE35961), HFC (GSE53381), MAT1A-KO (GSE63027), and GNMT-KO (GSE63027). The first human NASH study upregulated 4 out of 24 (16.6%) and downregulated 0 out of 24 (0%) acetyl-CoA generation enzymes. The second human NASH study upregulated 1 out of 24 (4.2%) and downregulated 0 out of 24 (0%) acetyl-CoA generation enzymes. The MCD+HFD model upregulated 3 out of 24 (12.5%) and downregulated 7 out of 24 (29.2%) acetyl-CoA generation enzymes. The MAT1A-KO model upregulated 1 out of 24 (4.2%) and downregulated 4 out of 24 (16.7%) acetyl-CoA generation enzymes. The GNMT-KO model upregulated 2 out of 24 (8.3%) and downregulated 0 out of 24 (0%) acetyl-CoA generation enzymes. The HFCD model upregulated 7 out of 24 (29.2%) and downregulated 2 out of 24 (8.3%) acetyl-CoA generation enzymes. The 24 acetyl-CoA generating enzyme genes are shown in Supplementary Table 3b. (c) The 10 out of 10 mRNAs (100%) of mevalonate pathway enzyme genes are regulated in human and mouse NASH/NAFLD. Mevalonate pathway enzyme gene expression (PMID 31153039 and PMID 29328908) for human NASH studies (GSE17470 and GSE63067), MCD+HFD (GSE35961), HFC (GSE53381), MAT1A-KO (GSE63027), and GNMT-KO (GSE63027). The first human NASH study upregulated 2 out of 10 (20%) and downregulated 3 out of 10 (30%) mevalonate pathway enzymes. The second human NASH study upregulated 1 out of 10 (10%) and downregulated 0 out of 10 (0%) mevalonate pathway enzymes. The MCD+HFD model upregulated 7 out of 10 (70%) and downregulated 1 out of 10 (10%) mevalonate pathway enzymes. The MAT1A-KO model upregulated 4 out of 10 (40%) and downregulated 1 out of 10 (10%) mevalonate pathway enzymes. The GNMT-KO model upregulated 3 out of 10 (30%) and downregulated 0 out of 24 (0%) mevalonate pathway enzymes. The HFCD model upregulated 7 out of 10 (70%) and downregulated 0 out of 10 (0%) mevalonate pathway enzymes. The 10 mevalonate pathway enzyme genes are shown in Supplementary Table 3C.
(a)PatientsMouseGEO IDGSE17470GSE63067GSE35961GSE63027GSE63027GSE53381ComparisonNASH vs. healthyNASH vs. healthyMCD+HFD vs. NCDMAT1A-KO vs. WTGNMT-KO vs. WTHFCD vs. NCDGene symbolACSS14.47-2.67ACSS2-2.55-1.623.45ADH72.761.641.481.39AHD1A7.97ADH1B2.82ADH1C4.09ALDH1A36.34ALDH1B14.381.622.24-1.63ALDH26.33-1.70-1.501.59ALDH3A22.032.05-2.875.86ALDH3B12.49ALDH7A13.81-4.071.65ALDOA3.69ALDOC2.211.641.323.96DLAT-2.861.40-2.341.401.36ENO33.934.11FBP23.13HK12.73HK22.10HKDC1-1.855.361.60LDHA-1.79LDHB4.99PCK27.054.38-1.33PDHA11.91PDHB-1.39-1.72PFKFB13.641.631.62PFKFB42.691.62PFKP7.30-2.27PGAM11.44PGK1-6.501.75PGM11.682.701.71-1.64PKM5.27-2.76ADH441.77-2.88-1.51-2.01ADH54.73-1.331.45ADH65.42-4.71ADPGK-2.07-1.82ALDH3B2-1.67ALDH9A1-1.351.92ALDOB2.051.50-1.661.30BPGM-1.32-1.371.32DLD2.58-2.021.68ENO1-1.551.36ENO2-2.13FBP12.24-2.231.49G6PC-5.06-1.59GALM2.16-1.721.37GAPDH-1.61GCK17.03-33.70GPI1.95HK32.71-3.37LDHC-2.04PANK1-1.391.77PCK1-1.88-1.45PFKFB21.74PFKFB31.982-1.65-2.48PFKL1.738-1.37PFKM-1.98-1.49-1.31PGAM222.09-2.02PGM21.52PKLR1.738-1.80SLC2A29.160-1.761.36% up25/71 (35.2%)5/71 (7%)20/71 (28.2%)6/71 (8.5%)7/71 (9.9%)17/71 (23.9%)% down5/71 (7%)0/71 (0%)24/71 (33.8%)13/71 (18.3%)4/71 (5.6%)8/71 (11.3%)(b)HumanMouseGEO IDGSE17470GSE63067GSE35961GSE63027GSE63027GSE53381ComparisonNASH vs. healthyNASH vs. healthyMCD+HFD vs. NCDMAT1A-KO vs. WTGNMT-KO vs. WTHFCD vs. NCDGene symbolBDH12.331.89ACSS14.47-2.67ALDH26.33-1.70-1.501.59ACSS2-2.55-1.623.45ACAA22.83-1.511.381.49HADH1.501.37ADH1B2.82GOT11.79ACO1-1.411.45ACLY-2.40-1.361.75GLS3.27-2.40IDH12.83-1.40-1.451.60GLUD1-1.55% up4/24 (16.6%)1/24 (4.16%)3/24 (12.5%)1/24 (4.2%)2/24 (8.3%)7/24 (29.2%)% down0/24 (0%)0/24 (0%)7/24 (29.2%)4/24 (16.7%)0/24 (0%)2/24 (8.3%)(c)HumanMouseGEO IDGSE17470GSE63067GSE35961GSE63027GSE63027GSE53381ComparisonNASH vs. healthyNASH vs. healthyMCD+HFD vs. NCDMAT1A-KO vs. WTGNMT-KO vs. WTHFCD vs. NCDGene symbolIDI12.7781.5710.974IDI2-1.655HMGCS1-2.4401.8952.647.177MVD-1.7954.0031.3815.819MVK4.0965.4511.8583.895HMGCR2.0911.7642.663ACAT12.1351.3741.528PMVK1.453.366ACAT21.63HMGCS27.434-1.908-1.7381.681% up2/10 (20%)1/10 (10%)7/10 (70%)4/10 (40%)3/10 (30%)7/10 (70%)% down3/10 (30%)0/10 (10%)1/10 (10%)1/10 (10%)0/10 (0%)0/10 (0%)Furthermore, the first human NASH upregulated 4 out of 24 (16.6%) acetyl-CoA generation enzyme genes. The second human NASH upregulated 1 out of 24 (4.2%) genes. The MCD+HFD upregulated 3 (12.5%) and downregulated 7 (29.2%) out of 24 genes. The MAT1A-KO upregulated 1 (4.2%) and downregulated 4 (16.7%) out of 24 genes. The HFCD upregulated 7 (29.2%) and downregulated 2 (8.3%) out of 24 genes (Table2(b)).In addition, the first human NASH upregulated 2 (20%) and downregulated 3 (30%) out of 10 mevalonate pathway enzyme genes. The second human NASH upregulated 1 out of 10 (10%) gene. The MCD+HFD upregulated 7 (70%) and downregulated 1 (10%) out of 10 genes. The MAT1A-KO upregulated 4 (40%) and downregulated 1 (10%) out of 10 genes. The GNMT-KO upregulated 3 out of 10 (30%) genes. The HFCD upregulated 7 out of 10 (70%) genes (Table2(c)).Statins (β-hydroxy β-methylglutaryl-CoA (HMG-CoA) reductase inhibitors) block cholesterol synthesis and mevalonate generation, preventing TI induction [69]; and statin treatments do not revert the TI phenotype in patients with familiar hypercholesterolemia [70]. We hypothesized that statins may modulate TI enzyme expression in NAFLD models. For this, we collected two microarray datasets of statin-treated human primary hepatocytes [71].Among the 2004 differentially expressed genes, 16 genes were overlapped with 99 TI enzyme genes (16.2%) [27], of which 15 (15.2%) genes were upregulated and one (1%) gene was downregulated by atorvastatin (Figures 5(a)–5(c)). In addition, among 2448 differentially expressed genes in rosuvastatin-treated human primary hepatocytes, 19 genes were overlapped with 99 TI enzyme genes (19.2%), of which 15 (15.2%) genes were upregulated and 4 (4%) genes were downregulated by rosuvastatin. Then, we used Venn diagram to determine the causative effects of statins in TI enzyme genes upregulated in NAFLD models. The results showed that the atorvastatin treatment overlapped with 13 TI enzyme genes upregulated in human and mouse NASH/NAFLD models. Similarly, rosuvastatin treatment overlapped with 15 TI enzyme genes upregulated in human and mouse NAFLD models. These results have demonstrated that TI inhibitors, statins, promote the expression of TI enzyme genes not only in mevalonate pathways but also in glycolysis pathways and acetyl-CoA generation pathway.Figure 5
Trained immunity enzyme genes were differentially expressed in statin treatment (HMG-CoA reductase inhibitor) in the mevalonate synthesis pathway and modulated in human and mouse NASH/NAFLD. (a) Venn diagram showed the overlapped trained immunity enzyme genes (PMID: 31153039) with differentially expressed genes (DEG) in statin (atorvastatin and rosuvastatin) treatment (GSE24187). (b) List of trained immunity enzyme genes differentially expressed in atorvastatin and rosuvastatin and upregulated in human and mouse NASH/NAFLD models. (c) Venn diagram showed that the upregulated trained immunity enzyme genes in human and mouse NASH/NAFLD overlapped with the differentially expressed genes in atorvastatin and rosuvastatin treatment.
(a)(b)(c)


## 3.5. Expression Analysis of Lipid Peroxidation Enzymes Indicates That There Are Two Representative NASH/NAFLD Models: The MCD+HFD Model Is Proinflammatory Cytokine- and Canonical and Noncanonical Inflammasome-Upregulated Model; in Contrast, the HFCD Model Is Lipid Peroxidation Enzyme- and TI Enzyme-Upregulated Model

Arachidonic acid (AA) undergoes oxidative metabolism under enzymatic or nonenzymatic mechanisms. Cyclooxygenase (COX), lipoxygenase (LOX), and cytochrome p450s (CYPs) mediated enzymatic peroxidation of AA to the respective metabolites. The 12-lipoxygenase (12-LO) and 15-lipoxygenase (15-LO) mediated conversion of AA to 12- and 15-hydroperoxyeicosatetraenoic acid (HpETE) and the downstream glutathione peroxidase- (GPx-) mediated conversion to the respective hydroxyeicosatetraenoic acids (12- and 15-HETE). ROS-induced lipid peroxidation (LPO) of 12- and 15-HpETE results in the generation of 4-hydroxynonenal (4-HNE) and 4-hydroxydodecadeinal (4-HDDE). These reactive aldehydes interact with and inactivate GPx, leading to an increased rate of 12-HpETE and 15-HpETE peroxidation [72]. The nonenzymatic peroxidation mediates conversion of AA to 4-HNE and MDA metabolites [73] (Figure 6).Figure 6
Schematic representation of arachidonic acid metabolism and lipid peroxidation. Cyclooxygenase (COX), lipoxygenase (LOX), and cytochrome p450s (CYPs) mediated enzymatic peroxidation of arachidonic acid (AA) to the respective metabolites. The 12-lipoxygenase (12-LO) and 15-lipoxygenase (15-LO) mediated conversion of AA to 12- and 15-hydroperoxyeicosatetraenoic acid (HpETE) and the downstream glutathione peroxidase- (GPx-) mediated conversion to the respective hydroxyeicosatetraenoic acids (12- and 15-HETE). Reactive oxygen species- (ROS-) induced lipid peroxidation (LPO) of 12- and 15-HpETE results in the generation of 4-hydroxynonenal (4-HNE) and 4-hydroxydodecadeinal (4-HDDE). These reactive aldehydes interact with and inactivate GPx, leading to an increased rate of 12-HpETE and 15-HpETE peroxidation (PMID: 29610056). The nonenzymatic peroxidation mediates conversion of AA to 4-HNE and MDA metabolites (PMID: 24999379).We hypothesized that downregulation of lipid peroxidation antioxidant enzymes CYPs and GPXs is inversely associated with upregulation of cytokines and chemokines in NAFLD and that upregulation of proinflammatory lipid peroxidation enzymes COXs and LOXs [74] is associated with upregulations of TI enzymes in NAFLD. To examine these hypotheses, we collected 26 COX lipid peroxidation enzymes, 44 cytochrome P450 CYP lipid peroxidation antioxidant (anti-inflammatory) enzymes [75, 76], 26 AA metabolism enzymes, and 8 hepatic glutathione peroxidase antioxidant (anti-inflammatory) enzymes [77].The MCD+HFD upregulated 2 (7.7%) and downregulated 9 (34.6%) COXs (Table3(a)). The MAT1A-KO upregulated 2 (7.7%) and downregulated 1 (3.8%) out of 26 COXs. The GNMT-KO upregulated 3 (11.5%) COXs. The HFCD upregulated 9 (34.6%) and downregulated 1 (3.8%) COXs.Table 3
(a) Cyclooxygenase (COX) lipid peroxidation enzymes were differentially expressed in mouse NASH/NAFLD models. The MCD+HFD model upregulated 2 out of 26 (7.7%) and downregulated 9 out of 26 (34.6%) COXs. The MAT1A-KO model upregulated 2 out of 26 (7.7%) and downregulated 1 out of 26 (3.8%) COXs. The GNMT-KO model upregulated 3 out of 26 (11.5%) and downregulated 0 out of 26 (0%) COXs. The HFCD model upregulated 9 out of 26 (34.6%) and downregulated 1 out of 26 (3.8%) COXs. (b) Hepatocyte glutathione peroxidase lipid antioxidant enzymes were differentially expressed in mouse NASH/NAFLD models. The MCD+HFD model upregulated 4 out of 8 (50%) and downregulated 0 out of 8 (0%) GPXs. The MAT1A-KO model upregulated 2 out of 8 (25%) and downregulated 0 out of 8 (0%) GPXs. The GNMT-KO model upregulated 1 out of 8 (12.5%) and downregulated 0 out of 8 (0%) GPXs. The HFCD model upregulated 0 out of 8 (0%) and downregulated 1 out of 8 (12.5%) GPXs. (c) Cytochrome P450 CYP lipid peroxidation antioxidant (anti-inflammatory) enzymes were differentially expressed in mouse NASH/NAFLD models. The MCD+HFD model upregulated 6 out of 44 (13.6%) and downregulated 13 out of 44 (29.5%) CYP lipid peroxidation antioxidant enzymes (CYP-LPAEs). The MAT1A-KO model upregulated 4 out of 44 (9.1%) and downregulated 8 out of 44 (18.2%) CYP-LPAEs. The GNMT-KO model upregulated 5 out of 44 (11.4%) and downregulated 8 out of 44 (18.2%) CYP-LPAEs. The HFCD model upregulated 22 out of 44 (50%) and downregulated 1 out of 44 (2.3%) CYP-LPAEs. (d) Arachidonic acid metabolism enzymes were differentially expressed in mouse NASH/NAFLD models. The MCD+HFD model upregulated 4 out of 26 (15.4%) and downregulated 11 out of 26 (42.3%) arachidonic acid metabolism enzymes. The MAT1A-KO model upregulated 4 out of 26 (15.4%) and downregulated 6 out of 26 (23.1%) arachidonic acid metabolism enzymes. The GNMT-KO model upregulated 4 out of 26 (15.4%) and downregulated 7 out of 26 (26.9%) arachidonic acid metabolism enzymes. The HFCD model upregulated 24 out of 26 (92.3%) and downregulated 0 out of 26 (0%) arachidonic acid metabolism enzymes.
(a)
COX lipid peroxidation enzymes
GEO IDGSE35961GSE63027GSE63027GSE53381ComparisonMCD+HFD vs. NCDMAT1A-KO vs. WTGNMT-KO vs. WTHFCD vs. NCDGene symbolCOX101.391COX111.388COX141.303COX151.372COX171.3021.364COX18-1.3681.477COX191.438COX4I1-1.314COX5A-1.584COX5B-1.5241.344COX6A1-1.5981.387COX6A21.378-2.695COX6B1-1.437COX6B22.1351.9231.409COX6C-1.334COX7A1-1.7301.439COX7A2-1.313COX7A2LCOX7C1.356COX8C-1.761% up2/26 (7.7%)2/26 (7.7%)3/26 (11.5%)9/26 (34.6%)% down9/26 (34.6%)1/26 (3.8%)0/26 (0%)1/26 (3.8%)(b)
Hepatocyte glutathione peroxidase lipid antioxidant enzyme
GEO IDGSE35961GSE63027GSE63027GSE53381ComparisonMCD+HFD vs. NCDMAT1A-KO vs. WTGNMT-KO vs. WTHFCD vs. NCDGene symbolGPX38.8572.387GPX41.8321.590GPX73.7632.196GPX82.266-1.693% up4/8 (50%)2/8 (25%)1/8 (12.5%)0/8 (0%)% down0/8 (0%)0/8 (0%)0/8 (0%)1/8 (12.5%)(c)
Cytochrome P450 CYP lipid peroxidation antioxidant (anti-inflammatory) enzymes
GEO IDGSE35961GSE63027GSE63027GSE53381ComparisonMCD+HFD vs. NCDMAT1A-KO vs. WTGNMT-KO vs. WTHFCD vs. NCDGene symbolCYP11A1-2.023CYP11B2-3.851-1.459CYP17A1-2.7817.210CYP1B13.0801.496CYP20A12.293CYP21A2CYP2B10-3.16012.086CYP2B139.4553.6487.655CYP2B948.3036.8114.811CYP2C29-4.350-2.577-1.7361.663CYP2C37-2.233-3.070-2.3951.834CYP2C38-1.531-2.7842.6942.499CYP2C392.7711.751CYP2C401.830CYP2C54-113.732-1.679-3.6811.907CYP2C551.681-2.0854.526CYP2C651.413CYP2F1CYP2J111.3471.621CYP2J2CYP2J5-4.063-1.5691.442CYP2R1-2.226-2.0321.441CYP2S1-1.756CYP2U1-6.007-1.9901.473CYP2W1CYP3A7CYP4A102.235CYP4A11CYP4A12A-7.4363.6641.6702.931CYP4A12B2.789CYP4A143.575-13.39216.11114.158CYP4A31-2.77313.473CYP4B1-1.598-1.619CYP4F13-2.8111.525CYP4F18% up6/44 (13.6%)4/44 (9.1%)5/44 (11.4%)22/44 (50%)% down13/44 (29.5%)8/44 (18.2%)8/44 (18.2%)1/44 (2.3%)(d)
Arachidonic acid metabolism enzymes
GEO IDGSE35961GSE63027GSE63027GSE53381MCD+HFD vs. NCDMAT1A-KO vs. WTGNMT-KO vs. WTHFCD vs. NCDGene symbolCYP2C37-2.233-3.070-2.3951.834CYP2C29-4.350-2.577-1.7361.663CYP2C54-113.732-1.679-3.6811.907CYP2J5-4.063-1.5691.442CYP4F13-2.8111.525CYP2U1-6.007-1.9901.473PLA2G4C-3.499LTA4H-2.953CYP 20-1.8181.432CYP4A12A-7.4363.6641.6702.931CYP2C38-1.531-2.7842.6942.499CYP2C651.413PLA2G61.429ALOX12B1.439PTGES21.508CYP2J111.3471.621CYP2C401.830CYP4A102.235CYP4A12B2.789CYP2C551.681-2.0854.526CYP2B10-3.16012.086CYP4A31-2.77313.473CYP4F18CYP2C392.7711.751CYP4A143.575-13.39216.11114.158CYP2B139.4553.6487.655CYP2B948.3036.8114.811% up4/26 (15.4%)4/26 (15.4%)4/26 (15.4%)24/26 (92.3%)% down11/26 (42.3%)6/26 (23.1%)7/26 (26.9%)0/26 (0%)The MCD+HFD upregulated 6 (13.6%) and downregulated 13 (29.5%) CYP lipid peroxidation antioxidant enzymes (CYP-LPAEs) (Table3(c)). The MAT1A-KO upregulated 4 (9.1%) and downregulated 8 (18.2%) CYP-LPAEs. The GNMT-KO upregulated 5 (11.4%) and downregulated 8 (18.2%) CYP-LPAEs. The HFCD upregulated 22 (50%) and downregulated 1 (2.3%) CYP-LPAEs (Figure 6(c)).The MCD+HFD upregulated 4 (15.4%) and downregulated 11 (42.3%) AA metabolism enzyme genes (Table3(d)). The MAT1A-KO upregulated 4 (15.4%) and downregulated 6 (23.1%) genes. The GNMT-KO upregulated 4 (15.4%) and downregulated 7 (26.9%) genes. The HFCD upregulated 24 (92.3%) genes (Figure 6(d)).The MCD+HFD upregulated 4 out of 8 (50%) GPXs (Table3(b)). The MAT1A-KO upregulated 2 (25%) GPXs. The GNMT-KO upregulated 1 out of 8 (12.5%) GPXs. The HFCD downregulated 1 out of 8 (12.5%) GPXs (Figure 6(e)).Taken together, our results have demonstrated that,first, proinflammatory MCD+HFD upregulates 7.7% COXs, 13.6% CYP-LPAEs, 15.4% AA metabolism enzymes, and 50% GPXs and downregulates 34.6% COXs, 29.5% CYP-LPAEs, and 42.3% AA metabolism enzymes, which are in big contrast with that of less inflammatory HFCD with upregulations of 34.6% COXs, 50% CYP-LPAEs, and 92.3% AA metabolism enzymes and downregulation of 3.8% COXs, 2.3% CYP-LPAEs, and 12.5% GPXs; second, significant upregulation of lipid peroxidation enzymes in HFCD is positively correlated with significant upregulation of TI enzymes and is, however, inversely associated with low-level upregulation of canonical inflammasome regulators (6.7%), noncanonical inflammasome regulators (7.1%), and 2.4% secretomic cytokines; third, no upregulation of GPXs in HFCD (Table 3(b)) is well correlated with other’s findings that GPX4 decreases lipid peroxidation and inflammation, and caspase-11-dependent pyroptosis mediates atherosclerosis [78], endotoxic shock [79], and septic death in GPX4-KO mice [33]; and fourth, the MCD+HFD and HFCD are the two representative models: the MCD+HFD is a model of NASH/NAFLD with upregulation of proinflammatory cytokines and canonical and noncanonical inflammasome regulators. In contrast, the HFCD is a model with upregulation of lipid peroxidation enzymes and TI enzymes.


## 3.6. As Upstream Master Regulators, Caspase-11 and Caspase-1 Partially Upregulate the Expressions of Cytokines, Chemokines, Canonical and Noncanonical Inflammasome Pathway Regulators, TI Enzymes, and Lipid Peroxidation Enzymes

It has been reported that caspase-11 mediates hepatocyte pyroptosis and promotes the progression of MCD diet-induced NASH/NAFLD [80]; and inflammasome-gasdermin D-pyroptosis [81] may regulate the progression of NASH/NAFLD [82]. In addition to promoting inflammation by enzymatic cleavage of pre-IL-1β, pre-IL-18, and N-terminal GSDMD, caspase-1 also cleaves as many as 114 substrates, suggesting that caspase-1 may regulate inflammation indirectly via all the protein substrate-mediated signaling [15], 38 interaction protein-mediated signaling [15], and GSDMD-secretome-mediated signaling [23]. However, an important question remains whether caspase-11 modulates the expression of cytokines, chemokines, and TI enzymes. We hypothesized that caspase-11 promotes inflammatory cytokines and chemokines upregulated in NASH/NAFLD and modulates the expression of canonical and noncanonical inflammasome pathway genes and TI enzymes. To examine this hypothesis, we collected all the cytokines and chemokines, canonical and noncanonical inflammasome pathway genes, and TI enzymes upregulated in NASH/NAFLD and crossed with caspase-11-KO and caspase-1-KO datasets (Supplementary Table 2b).Among 39 innatomic and secretomic cytokines and chemokines upregulated in NASH/NAFLD, caspase-11 deficiency upregulated 18 (46.2%) and downregulated 8 (20.5%) NASH/NAFLD-upregulated cytokines and chemokines. In comparison, caspase-1 deficiency upregulated one cytokine CXCL14 (2.6%) and downregulated one cytokine LTB (2.6%) NASH/NAFLD-upregulated cytokines and chemokines (Table4(a)). These results suggest that capsase-11 deficiency promotes the upregulation of most cytokines and chemokines induced by NASH/NAFLD. In addition, caspase-11 deficiency upregulated 5 (17.9%) and downregulated 5 (17.9%) canonical inflammasome regulators induced by NASH/NAFLD. However, caspase-1 deficiency did not modulate the expression of canonical inflammasome regulators induced by NASH/NAFLD (Table 4(b)). Furthermore, caspase-11 deficiency upregulated 6 (8.8%) and downregulated 31 (45.6%) TI enzymes upregulated in NASH/NAFLD. In comparison, caspase-1 deficiency upregulated 4 (5.9%) and downregulated 6 (8.8%) TI enzymes upregulated in NASH/NAFLD (Table 4(c)). Finally, caspase-11 deficiency upregulated 13 (28.3%) and downregulated 6 (13%) lipid peroxidation enzymes upregulated in NASH/NAFLD. In comparison, caspase-1 deficiency upregulated 3 (6.5%) and downregulated 6 (13%) lipid peroxidation enzymes upregulated in NASH/NAFLD (Table 4(d)).Table 4
(a) The cytokines and chemokines upregulated in NASH/NAFLD were differentially expressed in caspase-11 and caspase-1 knockout. Among 28 cytokines and chemokines upregulated in NASH/NAFLD, caspase-11 deficiency upregulated 18 out of 28 (64.3%) and downregulated 8 out of 28 (28.6%) NASH/NAFLD-upregulated cytokines and chemokines. Caspase-1 deficiency upregulated one cytokine CXCL14 out of 39 (2.6%) and downregulated one cytokine LTB out of 39 (2.6%) NASH/NAFLD-upregulated cytokines and chemokines. (b) The upregulated canonical inflammasome regulators in NASH/NAFLD were differentially expressed in caspase-11 knockout. Caspase-11 deficiency upregulated 5 out of 28 (17.9%) and downregulated five out of 28 (17.9%) canonical inflammasome regulators induced by NASH/NAFLD. Caspase-1 deficiency did not modulate the expression of canonical inflammasome regulators induced by NASH/NAFLD. (c) The upregulated trained immunity enzyme genes in NASH/NAFLD were differentially expressed in caspase-11 knockout and caspase-1 knockout mouse models. Caspase-11 deficiency upregulated 6 out of 68 (8.8%) and downregulated 31 out of 68 (45.6%) trained immunity enzymes upregulated in NASH/NAFLD. Caspase-1 deficiency upregulated 4 out of 68 (5.9%) and downregulated 6 out of 68 (8.8%) trained immunity enzymes upregulated in NASH/NAFLD. (d) The upregulated lipid peroxidation enzyme genes in NASH/NAFLD were differentially expressed in caspase-11 knockout and caspase-1 knockout mouse models. Caspase-11 deficiency upregulated 13 out of 46 (28.3%) and downregulated 6 out of 46 (13%) lipid peroxidation enzymes upregulated in NASH/NAFLD. Caspase-1 deficiency upregulated 3 out of 46 (6.5%) and downregulated 6 out of 46 (13%) lipid peroxidation enzymes upregulated in NASH/NAFLD.
(a)GEO IDGSE115094GSE32515ComparisonCaspase-11-KO vs. WTCaspase-1-KO vs. WTGene symbolCXCL101.98CCL205.07EDN13.45CXCL162.50CKLF2.83CCL54.81TIMP110.52CCL24.75Ccl75.99TNF2.72TSLP7.42CCL82.70IL274.15PF44.82CCL1713.15CCL33.58CCL47.09CXCL54.50CXCL143.219DKK 3-5.54CX3CL1-4.49CXCL12-3.50CSF1-3.04TNFSF10-2.23NAMPT-2.08WNT5A-1.79SPP1-5.87LTB-1.551% up18/39 (46.2%)1/39 (2.6%)% down8/39 (20.5)1/39 (2.6%)(b)Canonical inflammasome pathway regulators upregulated in NASHGEO IDGSE115094GSE32515ComparisonCaspase-11-KO vs. WTCaspase-1-KO vs. WTGene symbolTXN21.749PYCARD2.469Naip51.602Naip21.895CYBA3.599PRKCD-1.643ANTXR2-3.523CTSB-2.328NAMPT-2.083RIPK1-1.754% up5/28 (17.9%)0/28 (0%)% down5/28 (17.9%)0/28 (0%)(c)Trained immunityGEO IDGSE115094GSE32515PathwaysComparisonCaspase-11-KO vs. WTCaspase-1-KO vs. WTGlycolysisLDHB2.520PFKFB31.844ALDOA1.883FBP14.254ALDOA1.883ADH71.624PFKFB11.511PFKFB21.606HK1-2.535HK2-4.554HKDC1-1.550PDHA1-3.277PFKFB4-1.570PFKP-1.544PGAM1-2.484ENO1-1.689DLAT-1.770GALM-1.700PANK1-1.892PFKL-4.133PGM2-1.663PKLR-3.863-2.354ALDH1B1-2.881ALDH2-2.171ALDH3A2-1.840DLAT-1.770ACSS2-1.836ADH4-2.019GCK-3.334SLC2A2-1.789ACSS1-3.206ALDH1B1-2.881ALDH2-2.171ALDH3A2-1.840MevalonatePMVK1.699IDI1-3.248HMGCS1-1.818HMGCR-1.894ACAT1-1.917IDI1-3.248Acetyl-CoA generationGOT1-1.6841.701BDH1-1.964ACO1-1.601ACLY-2.203-2.33% up6/68 (8.8%)4/68 (5.9%)% down31/68 (45.6%)6/68 (8.8%)(d)Lipid peroxidationGEO IDGSE115094GSE32515EnzymesComparisonCaspase-11-KO vs. WTCaspase-1-KO vs. WTCOXsCOX145.806COX181.862COX5B6.309COX6A15.582COX6B212.677COX7A111.346COX7C11.363COX194.005COX15-3.650CYPsCYP2R14.459CYP2B107.4684.070CYP2C3915.441CYP4A311.831CYP2U1-2.833CYP1B1-2.058CYP20A1-2.466CYP2B9-3.075CYP2C37-6.050CYP2C38-3.200CYP2C54-2.609CYP2C55-6.187CYP4F18-1.575GPXGPX46.696GPX73.769GPX83.769GPX3-1.549AA metabolismPLA2G6-2.293% up13/46 (28.3%)3/46 (6.5%)% down6/46 (13%)6/46 (13%)Taken together, our results have demonstrated that,first, caspase-11 and caspase-1 partially modulate the expressions of cytokines, chemokines, canonical and noncanonical inflammasome pathway regulators, TI enzymes, and lipid peroxidation enzymes; second, in the context of NASH/NAFLD, caspase-11 has stronger capacities than caspase-1 in modulating inflammatory gene expressions; and third, experimental conditions in caspase-11 deficiency [79] and caspase-1 deficiency [83] were not exactly the same as NASH/NAFLD. Since we used inflammatory gene lists upregulated in the four NASH/NAFLD models for the analyses, our results are at least partially relevant to NASH/NAFLD pathology. Future work on transcriptomic analysis of caspase-11 deficiency and caspase-1 deficiency in NASH/NAFLD models, respectively, is needed to verify these findings presented here.


## 4. Discussion

In this study, we performed a panoramic database mining analysis on microarray data of both human NASH and NAFLD mouse models. We made the following significant findings: (i) human NASH and NAFLD mouse models upregulate both cytokines and chemokines and canonical secretome; (ii) pathway analysis indicated that human NASH can be classified into metabolic and immune NASH; MCD+HFD, GNMT-KO, MAT1A-KO, and HFCD can be classified into inflammatory NAFLD, SAM accumulation NAFLD, cholesterol/mevalonate NAFLD, and LXR/RXR-fatty acid β-oxidation NAFLD, respectively; (iii) canonical and noncanonical inflammasome pathways play differential roles in the pathogenesis of NASH/NAFLD; (iv) TI enzymes are significantly upregulated in human NASH and NAFLD mouse models; HFD upregulates TI enzymes more than cytokines, chemokines, and canonical and noncanonical inflammasome regulators; statins promote rather than suppress TI enzyme expression; (v) lipid peroxidation enzyme expression indicated that there are two representative NAFLD models: the MCD+HFD is a proinflammatory cytokine- and canonical and noncanonical inflammasome-upregulated model; in contrast, the HFCD is a lipid peroxidation enzyme- and TI enzyme-upregulated model; and (vi) caspase-11 and caspase-1 partially upregulate the expressions of cytokines, chemokines, canonical and noncanonical inflammasome pathway regulators, TI enzymes, and lipid peroxidation enzymes.We attempted to integrate all the findings presented here, our papers, and others’ reports and proposed a novel working model of multiple-hit TI model (Figure7) with enhanced inflammation for NASH/NAFLD development with a synergy between hyperlipidemia induced by HFD feeding and hypomethylation induced by nutrient and gene deficiencies related to methionine-homocysteine circle as we reported [84–91]. Hyperlipidemia and NAFLD are highly associated [1, 92–96]. HFCD model upregulated significantly TI enzymes and lipid peroxidation enzymes but did not significantly upregulate cytokines, chemokines, and canonical and noncanonical inflammasome regulators. Of note, most NAFLD mouse models have hyperlipidemia component (HFD) [45, 46, 97] (Figure 7(c)). Mechanistically, hyperlipidemia induced by HFD feeding acts as the first stimulation for TI [98, 99] (Figure 7(a)) via lysoPC stimulation [27], oxidized low-density lipoprotein (oxLDL) binding to CD36 [100], TLR4 [5]/TLR2 [101] and nucleotide-binding domain, leucine-rich-containing family, pyrin domain-containing 3 (NLRP3) inflammasomes [6, 26, 99, 102, 103] (Figure 7(b)).Figure 7
Our new integrated working model: multiple-hit trained immunity model with hyperlipidemia and hypomethylation (S-adenosylmethionine (SAM) decrease) for inflammation enhancement in various diet-induced and gene-deficient/genetic mouse models of nonalcoholic steatohepatitis (NASH)/nonalcoholic fatty liver disease (NAFLD). (a) Schematic presentation of trained immunity in NASH and NAFLD. Hyperlipidemia (high-fat diet) acts as the first stimuli to prime the innate immune cells and induce trained immunity. Hypomethylation that is caused by methionine-related nutrient deficiency and SAM synthase deficiency is regarded as second stimulation. This restimulation establishes trained immunity and metabolic reprogramming, leading to a large number of inflammatory cytokine secretion. (b) Epigenetic reprogramming pathway of trained immunity: hyperlipidemia and hyperglycemia inducers bind to their receptors in the cell membrane, increasing trained immunity-related glycolysis, acetyl-CoA production, and mevalonate pathways. Upregulation of three trained immunity pathways leads to increased trained immunity-related histone acetylation (H3K27ac and H3K14ac) and enhance proinflammatory innate immune response. (c) Overview figure shows the trained immunity model in detail. Abbreviations: Me: methylation; AC: acetylation; oxLDL: oxidized low-density lipoprotein; LPS: lipopolysaccharides; CD36: cluster of differentiation 36; TLR4: Toll-like receptor 4; NLRP3: Nod-like receptor family pyrin domain-containing 3; cROS: cytosolic reactive oxygen species; Akt: protein kinase B; mTOR: mechanistic target of rapamycin; ATP: adenosine triphosphate; GLUT1: glucose transporter 1; PTEN KO: phosphatase and tensin homology knockout; acetyl-CoA: acetyl coenzyme A; HAT: histone acetyltransferases; HMGCAR: 3-hydroxy-3-methylglutaryl CoA reductase; ACLY: ATP citrate lyase; H3K27: histone 3 lysine 27 acetylation; H314K: histone 3 lysine 14 acetylation; CBS: cystathionine-β-synthase; BHMT: betaine–homocysteine methyltransferase; MAT1A: methionine adenosyltransferase 1A; SAM: S-adenosylmethionine; SAH: S-adenosylhomocysteine; GNMT: glycine N-methyltransferase; PEMT: phosphatidylethanolamine N-methyltransferase; PC: phosphatidylcholine; CH3: methyl group.
(a)(b)(c)Others also reported that HFD and HFD plus HCD feeding promote the progression of NASH/NAFLD [45, 97], which are well correlated with our new working model that hyperlipidemic stimulations are functional as the first stimulation and second stimulation in TI to enhance innate immune responses. Very interestingly, we found that methionine deficiency diet, choline deficiency diet in MCD model, and MAT1A-KO all lead to S-adenosylmethionine (SAM) decrease or deficiency [104] in methionine-homocysteine cycle [91]. Furthermore, since glycine methylation is one of the reactions that contribute most to total transmethylation flux [105], GNMT1-KO lead to decreased glycine methylation [91] and global DNA hypomethylation [106]. Since SAM is cellular methyl donor we reported [85, 87–91, 107] and reviewed [84, 108], the SAM decrease/deficiency and weakened SAM function lead to decreased ratios of SAM/S-adenosylhomocysteine (SAH), cellular hypomethylation, decreased histone methylation, and consequently increased histone acetylation, presumably including TI-related histone 3 lysine 27 acetylation (H3K27ac) and H3K14ac, as we reviewed [26, 28] and reported [27, 29, 30, 109].Furthermore, increased histone acetylation by inhibiting histone deacetylase-2 (HDAC2) promotes macrophage infiltration and progression of NASH [110]. MAT1A and GNMT are relatively liver-specific enzymes; hypomethylation in the liver of MAT1A-KO and GNMT-KO can be striking. Hypomethylation may be an essential factor for liver injury but not essentially required for the establishment of TI since H3K4me3 also mediates TI [26] [111]. In support of our conclusion, it was reported that both hypermethylation and hypomethylation of lipid metabolism genes are found in obese patients with hypercholesterolemia [112]. Since hyperlipidemia can act alone to promote NASH/NAFLD progression and presumably the establishment of TI, both hyperlipidemia and hypomethylation act as second stimuli for TI (Figure 7(c)). As we indicated [27], the TI is a novel mechanism for qualifying any stimuli in the environments and underlying how chronic metabolic diseases [27, 28, 113] such as inflammation progression in NASH/NAFLD. Interestingly, we demonstrated that caspase-1 and caspase-11/4 not only serve as metabolic stress-derived DAMP sensors and inflammation initiators but also serve as upstream master regulators for at least partially upregulating inflammatory cytokines, chemokines, canonical and noncanonical inflammasome regulators, TI enzymes, and lipid peroxidation enzymes.One limitation of the current study is that, due to the low-throughput nature of verification techniques in all the research laboratories, we could not verify every result we identified [61, 114]. We acknowledge that carefully designed in vitro and in vivo experimental models will be needed to verify all the findings. Nevertheless, our findings provide novel insights on the roles of proinflammatory cytokines and chemokines [58, 59] and canonical secretome [61, 115, 116], canonical and noncanonical inflammasome pathways, TI enzymes, and lipid peroxidation enzymes in promoting NASH/NAFLD progression as well as novel targets for the future therapeutic interventions for NASH/NAFLD, metabolic diseases, transplantation, and cancers.



---

*Source: 3928323-2021-11-23.xml*