Gene-metabolite network analysis in different nonalcoholic fatty liver disease phenotypes

We sought to identify common key regulators and build a gene-metabolite network in different nonalcoholic fatty liver disease (NAFLD) phenotypes. We used a high-fat diet (HFD), a methionine-choline-deficient diet (MCDD) and streptozocin (STZ) to establish nonalcoholic fatty liver (NAFL), nonalcoholic steatohepatitis (NASH) and NAFL+type 2 diabetes mellitus (T2DM) in rat models, respectively. Transcriptomics and metabolomics analyses were performed in rat livers and serum. A functional network-based regulation model was constructed using Cytoscape with information derived from transcriptomics and metabolomics. The results revealed that 96 genes, 17 liver metabolites and 4 serum metabolites consistently changed in different NAFLD phenotypes (>2-fold, P<0.05). Gene-metabolite network analysis identified ccl2 and jun as hubs with the largest connections to other genes, which were mainly involved in tumor necrosis factor, P53, nuclear factor-kappa B, chemokine, peroxisome proliferator activated receptor and Toll-like receptor signaling pathways. The specifically regulated genes and metabolites in different NAFLD phenotypes constructed their own networks, which were mainly involved in the lipid and fatty acid metabolism in HFD models, the inflammatory and immune response in MCDD models, and the AMPK signaling pathway and response to insulin in HFD+STZ models. Our study identified networks showing the general and specific characteristics in different NAFLD phenotypes, complementing the genetic and metabolic features in NAFLD with hepatic and extra-hepatic manifestations.