Investigating the deregulation of metabolic tasks via Minimum Network Enrichment Analysis (MiNEA) as applied to nonalcoholic fatty liver disease using mouse and human omics data

Nonalcoholic fatty liver disease (NAFLD) is associated with metabolic syndromes spanning a wide spectrum of diseases, from simple steatosis to the more complex nonalcoholic steatohepatitis. To identify the deregulation that occurs in metabolic processes at the molecular level that give rise to these various NAFLD phenotypes, algorithms such as pathway enrichment analysis (PEA) can be used. These analyses require the use of predefined pathway maps, which are composed of reactions describing metabolic processes/subsystems. Unfortunately, the annotation of the metabolic subsystems can differ depending on the pathway database used, making these approaches subject to biases associated with different pathway annotations, and these methods cannot capture the balancing of cofactors and byproducts through the complex nature and interactions of genome-scale metabolic networks (GEMs). Here, we introduce a framework entitled Minimum Network Enrichment Analysis (MiNEA) that is applied to GEMs to generate all possible alternative minimal networks (MiNs), which are possible and feasible networks composed of all the reactions pertaining to various metabolic subsystems that can synthesize a target metabolite. We applied MiNEA to investigate deregulated MiNs and to identify key regulators in different NAFLD phenotypes, such as a fatty liver and liver inflammation, in both humans and mice by integrating condition-specific transcriptomics data from liver samples. We identified key deregulations in the synthesis of cholesteryl esters, cholesterol, and hexadecanoate in both humans and mice, and we found that key regulators of the hydrogen peroxide synthesis network were regulated differently in humans and mice. We further identified which MiNs demonstrate the general and specific characteristics of the different NAFLD phenotypes. MiNEA is applicable to any GEM and to any desired target metabolite, making MiNEA flexible enough to study condition-specific metabolism for any given disease or organism. Author Summary This work aims to introduce a network-based enrichment analysis using metabolic networks and transcriptomics data. Previous pathways/subsystems enrichment methods use predefined gene annotations of metabolic processes and gene annotations can differ based on different resources and can produce bias in pathways definitions. Thus, we introduce a framework, Minimum Network Enrichment Analysis (MiNEA), which first finds all possible minimal-size networks for a given metabolic process/task and then identifies deregulated minimal networks using deregulated genes between two conditions. MiNEA also identifies the deregulation in key reactions that are overlapped across all possible minimal-size networks. We applied MiNEA to identify deregulated metabolic tasks and their synthesis networks in the steatosis and nonalcoholic steatohepatitis (NASH) disease using a metabolic network and transcriptomics data of mouse and human liver samples. We identified key regulators of NASH form the synthesis networks of hydrogen peroxide and ceramide in both humans and mice. We also identified opposite deregulation in NASH for the phosphatidylserine synthesis network between humans and mice. MiNEA finds synthesis networks for a given target metabolite and due to this it is flexible to study deregulation in different phenotypes. MiNEA can be widely applicable for studying context-specific metabolism for any organism because the metabolic networks and context-specific gene expression data are now available for many organisms.