- Research
- Open access
- Published:
Identification of significant modules and hub genes involved in hepatic encephalopathy using WGCNA
European Journal of Medical Research volume 27, Article number: 264 (2022)
Abstract
Background
Hepatic encephalopathy (HE) is a reversible syndrome of brain dysfunction caused by advanced liver disease. Weighted gene co-expression network analysis (WGCNA) could establish a robust co-expression network to identify the hub genes and underlying biological functions. This study was aimed to explore the potential therapeutic targets in HE by WGCNA.
Results
The green and brown modules were found to be significantly associated with the development of HE. Functional enrichment analyses suggested the neuroinflammation, neuroimmune, extracellular matrix (ECM), and coagulation cascade were involved in HE. CYBB and FOXO1 were calculated as hub genes, which were upregulated in the HE patients. Tamibarotene and vitamin E were suggested as possible drug candidates to alleviate HE.
Conclusions
It is the first time to analyze transcriptomic data of HE by WGCNA. Our study not only promoted the current understanding of neuroinflammation in HE, but also provided the first evidence that CYBB and FOXO1 played pivotal roles in the pathogenesis of HE, which might be potential biomarkers and therapeutic targets. Tamibarotene might be a novel drug compound against HE.
Background
Hepatic encephalopathy (HE) is a common but serious complication of advanced liver disease due to liver insufficiency, which is accompanied by a spectrum of neuropsychiatric manifestations and leads to high morbidity and mortality [1,2,3]. Therefore, effective therapeutic strategies are essential to improve the prognosis of HE.
In-depth understanding of pathophysiology in HE is the basis for effective therapeutic strategies. Hyperammonemia remains the central pathophysiological process in the pathogenesis of HE. However, previous studies have highlighted the limited predictive value of blood ammonia in HE, which suggested the existence and importance of other factors [4, 5]. Besides hyperammonemia, overwhelming evidence has demonstrated the important role of systemic and central inflammation, malnutrition, gut–liver–brain axis, and neurotransmitters during the development of HE [1, 5, 6].
During the past decades, an emerging role of neuroinflammation in the development of HE has been demonstrated [7, 8]. As a key pathogenic factor in HE, neuroinflammation can be induced by inflammatory cytokines, chemokines, and oxidative stress, which involves various cell types in brain including glial cells and peripherally recruited immune cells [1, 8, 9]. Meanwhile, several anti-inflammatory strategies were proven to be effective in ameliorating HE. Rifaximin protected against HE in cirrhotic patients by modulating gut-derived inflammation [10]. Babao Dan exerted anti-inflammatory effects and improved clinical minimal HE [11]. However, those available therapies were proved to be mainly against peripheral inflammation. Novel strategies were still lacking to directly alleviate the central inflammation. Moreover, most of anti-inflammatory therapies nonspecifically suppressed the inflammatory response, while targeted therapies were lacking to precisely impede the development of HE. Therefore, seeking for novel therapeutic targets against central inflammation is an important issue.
In recent years, a powerful method for systematic analysis called weighted gene co-expression network analysis (WGCNA) has been widely applied in bioinformatic analysis of various diseases [12]. The robust co-expression network could cluster genes with similar expression patterns into modules to identify the underlying biological functions. Moreover, the hub gene identification would reduce the bias brought by confounding factors, offering insight into potential therapeutic targets of the interested diseases. Therefore, in the current study, we used WGCNA to identify the clinically significant modules and hub genes of HE in cirrhotic patients, which might provide further evidence for novel therapeutic targets of HE.
Results
Construction of the HE co-expression network and identification of the modules
After the removal of the outlier sample (GSM1027458) by clustering analysis, a total of 30 samples were involved in the WGCNA, including 12 healthy control samples, 7 cirrhosis samples, and 12 HE samples (Additional file 1: Figure S2). A total of 16,416 genes were annotated, when 25% of the genes with the greatest variance were selected to construct the co-expression network. Based on the 4104 selected genes, the optimal soft-thresholding power β was selected as 4 to ensure the scale-free topology (Fig. 1A). When the power was 4, the scale-free topology (R2) was 0.93 (Fig. 1B). After the optimal soft-thresholding power β was selected, the adjacency matrix of selected genes was constructed and subsequently transformed to TOM. Thus, dissTOM of 4104 selected genes was obtained and shown by the network heatmap plot (Additional file 1: Figure S1). Then, the modules with a minimum size of 30 were identified by the dynamic tree cut method (Fig. 1C). All 4104 genes were divided into certain clusters, while no gene was assigned to the grey module. Moreover, the average linkage hierarchical clustering was performed based on the average distance between the modules. Finally, seven modules of co-expressed genes were obtained after merging similar modules with a merging threshold function at 0.25 (Fig. 1C).
Identification of the HE modules
After merging similar modules, the MEs of the seven obtained modules were calculated (Fig. 2C, D, Additional file 1: Figure S3A-E). To identify the clinically significant modules associated with HE, the Pearson’s correlations between the MEs of the seven modules and the clinical traits were obtained (Fig. 1D). According the module–trait relationships, the brown module (r = 0.7, P = 2e−05) and the green module (r = 0.68, P = 4e−05) were strongly associated with HE (Fig. 1D). Meanwhile, according to the correlations between the MEs and the interested clinical trait HE, a hierarchical clustering dendrogram and a heatmap were constructed, indicating that the brown module and the green module were tightly correlated to HE (Fig. 2A). Moreover, GS for HE was also calculated in each module. The mean absolute values of GS in each module were calculated and visualized, suggesting that the brown and green modules were the clinically significant modules with the highest GS for HE (Fig. 2B). Furthermore, the MM values were obtained, followed by the correlation analysis between the GS for HE and the MM for genes in the seven modules, respectively (Fig. 2E, F, Additional file 1: Figure S4A-E). With a comprehensive consideration of the correlation coefficients and the p-values, the brown module (r = 0.71, P = 2.6e−69) and the green module (r = 0.7, P = 3.6 e−105) were selected as clinically significant modules, which were consistent with the above results (Fig. 2E, F). Taken together, with three different approaches, we could conclude that the brown module and the green module were HE modules exhibiting the strongest associations with HE. Otherwise, the functional enrichment analysis of the other five modules did not achieve satisfactory results. Thus, we would focus on these two HE modules in the following investigations to further identify the biological functions and hub genes.
Functional enrichment analysis and hub genes identification of the brown module
Based on the clustered genes in the brown module, functional enrichment annotations were identified by GO enrichment analysis and KEGG pathway enrichment analysis. The top 10 significant terms of BP in GO enrichment analysis were visualized, suggesting that the genes in the brown module were associated with neuroinflammation- and neuroimmune-related biological processes (Fig. 3A). For KEGG pathway enrichment analysis, the top 10 enriched pathways were mainly related to neuroinflammation and neuroimmune as well (Fig. 3B). To identify the hub gene in the brown module, top 20 genes with the highest IC were extracted to construct the key network, which were defined as central genes. The connections among the central genes were visualized (Fig. 3C, Additional file 1: Table S1). Finally, CYBB was selected as the hub gene of the brown module as it had the highest degree of connections among 20 genes. Further analysis focusing on CYBB was performed. It was demonstrated that the expression level of CYBB was significantly upregulated in HE samples when compared to health control samples and cirrhosis samples, respectively (Fig. 3D). In addition, ROC analysis highlighted that CYBB could predict the diagnosis of HE (AUC = 0.78) (Fig. 3E). Meanwhile, CYBB was also capable of distinguishing between HE and cirrhosis (AUC = 0.974) (Fig. 3F).
Functional enrichment analysis and hub genes identification of the green module
GO enrichment analysis and KEGG pathway enrichment analysis were utilized to identify the functional enrichment annotations in the green module. The GO enrichment analysis showed that the green module was mainly associated with neuroinflammation-related biological processes (Fig. 4A). For KEGG pathway enrichment analysis, neuroinflammation-, ECM-, and coagulation-related pathways were significantly enriched in the green module (Fig. 4B). To identify the hub gene, the top 20 genes with the highest IC were visualized, which were defined as central genes (Fig. 4C, Additional file 1: Table S2). After constructing the key network, TUBA1C was the gene with the highest degree of connections. However, there was no significant difference among the expression levels of TUBA1C in normal, cirrhosis, and HE groups. Alternatively, FOXO1 was selected as the hub gene of the green module, which had the second highest degree of connections. An elevated expression level of FOXO1 in HE samples was observed compared to healthy control samples (Fig. 4D). Meanwhile, ROC analysis showed that the high expression of FOXO1 could well predict the existence of HE (AUC = 0.902) and distinguish HE from cirrhosis (AUC = 0.805) (Fig. 4E, F).
Identification of candidate drugs
Based on the central genes in the brown and the green modules, candidate drugs were identified, respectively. For both HE modules, top 10 drug molecules with significant p-values and q-values were shown, which might be potential therapeutic agents against HE. Tamibarotene was the top candidate drug of the brown module with the lowest p-value and q-value (Table 1). Meanwhile, VITAMIN E was identified as the top candidate drug of the green module (Table 2).
Changes in glial molecular markers
Overwhelming evidence has suggested the important role of glial cell activation during the pathogenesis of HE [8, 13, 14]. Therefore, after we obtained a thorough picture of the gene expression patterns in HE, further analysis was conducted to investigate the specific expression patterns of astrocyte activation and microglial polarization. Firstly, we adopted specific markers of reactive astrocyte subpopulations, which could represent the A1 and A2 subtypes, respectively [15, 16]. The expression levels of the astrocyte-specific markers among three groups were shown by line charts (Fig. 5A, B). We could observe that the A1 subtype-related markers (FKBP5, GBP2, PSMB8, and SRGN) were remarkably elevated in HE samples compared to healthy control samples and cirrhosis samples (Fig. 5A). On the other hand, the A2 subtype-related markers exhibited mixed patterns of changing (Fig. 5B). Only TGM1 showed a downward trend, while PTGS2, PTX3, and SLC10A6 did not show notable changes. Interestingly, B3GNT5 and EMP1 were upregulated in HE samples, which was probably a compensatory upregulation. Meanwhile, the polarization status of microglia was also investigated. The microglia-specific markers for M1 and M2 subtypes were selected to represent the subpopulations of microglia [16]. According to the line chart, it was obvious that most of the markers in the M1 group (HLA-DQB1, IL-6, and TNF) were upregulated in HE samples. Only NOS2 exhibited a slightly downward trend (Fig. 5C). As for the markers in the M2 group, only MRC1 was elevated in HE samples, especially when compared to cirrhosis samples. ARG1 and TGFB1 did not show notable changes among three groups (Fig. 5D).
Hub genes verification and GSEA in the rat model of HE
To verify the hub genes in the rat model, HE in rats was induced by BDL. Based on transcriptomic data of the rat brains, we compared the expression levels of the central genes in two HE modules. Four central genes in the brown module and one central gene in the green module were not identified in RNA sequencing that might be due to low expression level. According to the heatmaps, both CYBB and FOXO1 were significantly upregulated in the brains of BDL rats (Fig. 6A, B). The upregulation of CYBB and FOXO1 was further validated by RNA sequencing of the brains of BDL rats (Additional file 1: Figure S5A, B). Moreover, GSEA was performed to figure out the potential functions. The gene sets of dopaminergic synapse and spinocerebellar ataxia were enriched and downregulated, while glutathione metabolism and tyrosine metabolism were enriched as upregulated gene sets in BDL rats (Fig. 6C–F).
Discussion
In the current study, we found that neuroinflammation played a pivotal role in the pathogenesis of HE. CYBB and FOXO1 might be major regulators during the development from cirrhosis to HE. Moreover, microglia and astrocytes were activated to a pro-inflammatory and neurotoxic status, which facilitated the progression of HE. To the best of our knowledge, it was the first time to use WGCNA to analyze brain transcriptomic datasets from human, which emphasized the importance of neuroinflammation in HE. Furthermore, our study provided the first evidence that CYBB and FOXO1 might be potential biomarkers and therapeutic targets of HE in cirrhotic patients.
In the functional enrichment analysis, our results suggested a remarkable overlap of enriched biological processes and signaling pathways between two HE modules. Interestingly, we also observed that the brown module was negatively correlated with cirrhosis, while the green module was negatively correlated with healthy control. Taken together, we could conclude that there were similar underlying mechanisms during the development of cirrhosis and the development from cirrhosis to HE, including neuroinflammation, neuroimmune, ECM, and coagulation cascade.
According to our results, the neuroinflammation in HE was mainly associated with inflammatory cytokine release, macrophage activation, and inflammatory signaling pathways. Overwhelming evidence has suggested the vital role of neuroinflammation during the pathogenesis of HE [17, 18]. Our previous study also highlighted the importance inflammation by functional enrichment analysis based on lncRNAs in the serum transcriptome of HE, which showed similar results as those in the current study [6]. Moreover, Balzano et al. has reported the protective effects of anti-TNF-α treatment against inflammatory cytokine release and microglial activation in the rat brains [19]. We also found that VEGF-C overexpression could alleviate neuroinflammation by enhancing meningeal lymphatic drainage, which eventually improved HE in cirrhotic rats [8]. Therefore, it was obvious that neuroinflammation was a major contributor to the progression of HE. Anti-inflammatory treatments could be effective therapeutic strategies.
In addition to the significant role of neuroinflammation in HE, the importance of neuroimmune, ECM, and coagulation cascade should not be neglected. It was reported that immune paralysis, defined by decreased HLA-DR expression on monocytes, was associated with overt HE (OHE) in liver cirrhosis and ACLF [20, 21]. Moreover, elevated serum levels of ECM-associated proteins were observed in HE, including osteopontin and 7S domain of type IV collagen (4COL7S) [22, 23]. Otherwise, as coagulopathy was a common comorbidity of liver cirrhosis, disturbed coagulation in HE might be attributed to liver insufficiency [24, 25]. Interestingly, the serum level of 4COL7S was also correlated to coagulopathy in ALF, suggesting a complicated regulatory network underlying the pathophysiology of HE [23].
Following the functional enrichment analysis, we performed further analysis in HE modules. In the brown module, CYBB was regarded as the hub gene. CYBB, also known as gp91-phox, is the catalytic subunit of the phagocyte-like NADPH oxidase (NOX2) [26]. Genetic deletion and overexpression of CYBB could modulate the expression of NOX2 [27,28,29]. In another word, the expression level and catalytic activity of NOX2 largely depended on the expression of CYBB. The upregulation of CYBB we observed in our study suggested an upregulation of NOX2 in the brain of HE patients. Therefore, our results showed that brain NOX2 played an essential role in the pathogenesis of HE. Previous in vitro studies have demonstrated that ammonia could induce the activity of NOX family in astrocytes, which could be reversed by NOX inhibitor apocynin [30, 31]. Moreover, reduction of ammonia by ornithine phenylacetate also led to a decrease of brain NOX1 in BDL rats [32]. However, to the best of our knowledge, the direct relationship between NOX2 and HE has never been reported. NOX2 could mediate neuroinflammation in various CNS disorders, when it most prominently affected astrocytes and microglia via ROS generation [33,34,35,36]. Thus, we could speculate that brain NOX2 was an important positive regulator of oxidative stress in HE, which exacerbated neuroinflammation. Interestingly, our GSEA in BDL rats also highlighted an upregulation in glutathione metabolism, which was an essential biological process involved in oxidative stress. Furthermore, various kinds of NOX2 inhibitors have been developed [37, 38]. The protective effects of NOX2 inhibitors against neuroinflammation have been revealed, while the microglial cell was a common target [39,40,41]. Thus, NOX2 could be a potential therapeutic target of HE, which warranted further translational study.
In the green module, the transcription factor FOXO1 was identified as the hub gene, which was highly expressed in HE samples and BDL rats. As an important transcriptional activator, FOXO1 played a vital role in metabolism modulation, immune response, cell fate decisions, etc. [42]. In CNS disorders, FOXO1 could balance energy metabolism, activate neuroinflammation, induce autophagy and exert anti-oxidative effects [43,44,45]. However, both upstream and downstream regulatory mechanisms of FOXO1 in HE remain unclear. Previous study highlighted the importance of FOXO1 phosphorylation to modulate the transcriptional activity, when PI3K/Akt and MAPK signaling pathways were significantly enriched in the green module [46, 47]. For the downstream regulation of biological processes, FOXO1 was a transcriptional activator of GSDMD gene in microglia, which enhanced microglia pyroptosis and neuroinflammation [48]. On the other hand, the elevated FOXO1 in HE might be an anti-oxidative signal, when FOXO1 was considered as a sensing element in anti-oxidative signaling [42]. However, contradictory results have been observed for the role of FOXO1 in oxidative stress [42, 49,50,51]. Therefore, loss-of-function experiments were required to figure out the exact role of FOXO1 in HE. CUT&Tag profiling or chromatin immunoprecipitation sequencing might help to identify the transcriptional targets. Moreover, FOXO1 could directly regulate tyrosine hydroxylase in dopaminergic neurons, when GSEA in BDL rats suggested the alterations of dopaminergic synapse and tyrosine metabolism in HE [52]. Furthermore, potential translational value of FOXO1 should not be neglected. Targeted interventions might affect various related biological processes regulated by FOXO1. Several therapies antagonizing FOXO1 were under research in cancers and metabolic disorders [53]. However, for the clinical application of FOXO1-targeted therapy in HE, there is still a long way to go.
Additionally, tamibarotene and vitamin E were identified as the peak drug candidates of two HE modules, respectively. Overwhelming evidence has suggested that vitamin E was an anti-oxidant, which could alleviate oxidative stress and astrocyte swelling induced by ammonia in vitro [54, 55]. Animal studies also highlighted the anti-oxidative effects of vitamin E to protect against TAA-induced acute HE in the rat model [56, 57]. Furthermore, the combination of vitamin E with other vitamins and zinc could significantly improve HE in cirrhotic patients [58]. Tamibarotene, also known as Am80, was a retinoic acid receptor (RAR) agonist to treat acute promyelocytic leukemia (APL) [59]. Tamibarotene could protect dopaminergic neurons from microglia-mediated inflammatory damages [60]. Previous clinical studies have also reported the anti-inflammatory effects of tamibarotene in CNS disorders including Alzheimer’s disease and intracerebral hemorrhage [61, 62]. However, the effects of tamibarotene on HE remained unknown, which warranted further clinical trials.
Furthermore, it was worth mentioning that the vital role of glial cells in HE has already been uncovered, especially for microglia and astrocytes [8, 14]. Activated microglial cells were capable of inducing neurotoxic reactive astrocytes in neuroinflammation [15]. Thus, we also investigated the changes in glial markers, which showed a remarkable upregulation of M1- and A1-specific markers. On the other hand, the changes in M2- and A2-specific markers exhibited mixed patterns. In conclusion, our results revealed that glial cells were polarized to pro-inflammatory and neurotoxic statuses during the pathogenesis of HE, which further emphasized the importance of glial cells in HE.
The current study does have some limitations. Firstly, due to the lack of clinical characteristics, we could not identify the hub genes related to those more specific clinical parameters including clinical grading and the psychometric hepatic encephalopathy score (PHES). Secondly, gain- and loss-of-function experiments in animal models or cultured cells were also needed to further validate the importance of the hub genes in HE. Thirdly, the gene expression profiles for WGCNA were measured in post-mortem brain samples from human. There were confounding factors that would influence the consistency of sequencing results among the samples.
Conclusions
In our study, we have identified two clinically significant modules related to HE based on WGCNA, and confirmed the important role of neuroinflammation and neuroimmune in the development of HE. In addition, CYBB and FOXO1 were identified as hub genes and highly expressed in patients with HE. Tamibarotene might be a potential therapeutic agent against HE. Finally, we hope the current study could provide further evidence and ideas for future investigations in targeted therapy of HE.
Materials and methods
Dataset download and pre-processing
The HE expression profiles and the corresponding clinical information of the GSE41919 and GSE57193 datasets were systemically extracted from the GEO (http://www.ncbi.nlm.nih.gov/geo/) database. The probe IDs were converted into gene symbols. Expression values of genes corresponding to multiple probes were calculated by averaging the probe values. For the missing data in the gene expression matrices, imputed data were estimated by k-nearest neighbor (KNN) approach with a k = 10 [63]. The ComBat package was utilized to remove the batch effect and merge them into a single dataset for further analysis [64]. The merged dataset containing 31 samples consisted of three clinically defined groups: healthy control (n = 12), cirrhosis (n = 7), and HE (n = 12). According to cluster analysis, one sample (GSM1027458) was removed as an outlier (Additional file 1: Figure S2). Finally, a total of 30 samples were involved in the construction of co-expression network and further analysis, which included 12 healthy control samples, 7 cirrhosis samples, and 11 HE samples.
Construction of co-expression network
WGCNA R package were used to construct the weight gene co-expression network [12]. Instead of using all 16,416 annotated genes to construct the co-expression network, we selected a quarter of the genes with the greatest variance. A total of 4104 genes were finally involved in the WGCNA. Firstly, after the removal of the outlier, a correlation matrix was constructed according to Pearson’s correlations between all pairwise genes. A proper weighting coefficient β was selected to satisfy the scale-free topology. Then the adjacency matrix was constructed based on the correlation matrix with the weighting coefficient β. The adjacency matrix was subsequently transformed to topological overlap matrix (TOM). The average linkage hierarchical clustering was utilized to identify the genes with similar expression patterns and obtain hierarchical clustering tree based on TOM-based dissimilarity (dissTOM). Then the modules were identified by dynamic tree-cutting method with a minimum size of 30. Furthermore, the module eigengenes (MEs) of each module were calculated, which were defined as the first principal components and representative of the overall expression levels of the corresponding modules. Cluster analysis was performed based on the average distance between modules, which was determined by the Pearson’s correlation analysis between MEs of divided modules. According to the cluster analysis, the modules with high similarity were merged with a merging threshold function at 0.25 [65].
Identification of clinically significant modules
After merging the modules with high similarity based on the average distance between the MEs of each module, the MEs of the modules were calculated again. The Pearson’s correlation analysis was conducted between each clinical traits and the MEs of each module. P < 0.05 was considered as a significant correlation between the module and the clinical trait. Moreover, gene significance (GS) was determined by Pearson’s correlation analysis between the clinical traits and the expression levels of each gene in the modules, respectively [12]. The mean absolute values of GS in each module were calculated to represent the correlation between the module and the clinical trait in another way.
Identification of hub genes
Hub genes were defined as the genes highly interconnected to other genes of the module. To identify the hub genes, module membership (MM) was calculated to evaluate the correlation between genes and modules. MM was determined by the correlation between the expression levels of genes and the MEs of each module [5]. And the intramodular connectivity (IC) was determined by the sum of the correlation coefficients with other genes in the given molecule. A gene with a higher IC tended to have higher MM values, which suggested a more decisive role in the given module [66]. The top 20 genes with the highest IC were selected as central genes for further analysis. The central genes and the connections among them were subsequently imported into Cytoscape 3.8.0 to visualized the key networks. Finally, the hub gene was defined as the gene of the highest degree within the key networks.
Identification of candidate drugs
According to the top 20 genes with the highest IC in two HE modules, drug compounds recognition was performed, respectively, to identify the drugs with potential therapeutic effects on HE. Candidate drugs were predicted with the Drug Signatures Database (DSigDB) on the Enrichr platform (https://maayanlab.cloud/Enrichr/) [67, 68]. The final results were generated with p-values indicating the correlation of the gene with the predicted drug. The drugs were deemed as having potential therapeutic effects according to the p-values.
Functional enrichment analysis of clinically significant modules
Functional enrichment analysis including gene ontology (GO) enrichment analysis and Kyoto encyclopedia of genes and genomes (KEGG) pathway enrichment analysis were performed to investigate the biological functions and related signaling pathways of the genes within the clinically significant modules. Based on clusterProfiler package, central genes with MM > 0.6 were selected for functional enrichment analysis [69]. The cut-off criteria for identification of significantly enriched GO annotations and KEGG pathways were set as P < 0.05.
Transcriptomic analysis of BDL rats
The transcriptomic dataset of the brains from BDL rats (GSE149227) was downloaded from the GEO database [8]. FastQC software was used for quality control checks on raw fastq files. The reference genome of rats (RGSC 6.0/rn6) was acquired from the database at the University of California, Santa Cruz (https://genome.ucsc.edu/cgi-bin/hgTables). HISAT2 was utilized for the construction of genome index and the alignment of reads to the reference genome. The featureCounts command was used to calculate the number of fragments mapped to corresponding genes. Subsequently, based on the processed RNA-sequencing data, DESeq2 package was utilized to compare the expression levels of the central genes in the clinically significant modules between control group and BDL group [70]. The heatmaps were drawn by pheatmap package. Meanwhile, gene set enrichment analysis (GSEA) was also conducted by clusterProfiler package [69]. The gene set database consisted of the chemical and genetic perturbations of curated gene sets, biological process, cellular component, and molecular function in the Molecular Signatures Database. The number of permutations was set as 1000, while other parameters were set as default values. The cut-off criteria for identification of the statistically significant gene sets were set as P < 0.05. Figures of GSEA were drawn by enrichplot package.
Statistical analysis
R language 3.4.0 and GraphPad Prism 8.0.1 (San Diego, CA) were utilized for the statistical analysis. All the continuous variables were presented as the mean ± SEM. For two-group comparisons, unpaired Student’s t test was used, while one-way ANOVA followed by post hoc Tukey's test were used for multiple comparisons among three or more groups. The correlations were analyzed by the Pearson’s correlation analysis. Receiver operating characteristic (ROC) analysis was implemented using pROC package [71]. P < 0.05 indicated statistically significant differences.
Availability of data and materials
All the datasets involved in this study are available in the GEO (http://www.ncbi.nlm.nih.gov/geo/) repository.
References
Rose CF, Amodio P, Bajaj JS, Dhiman RK, Montagnese S, Taylor-Robinson SD, et al. Hepatic encephalopathy: Novel insights into classification, pathophysiology and therapy. J Hepatol. 2020;73(6):1526–47.
Felipo V. Hepatic encephalopathy: effects of liver failure on brain function. Nat Rev Neurosci. 2013;14(12):851–8.
Zhang G, Li Y, Zhang X, Huang L, Cheng Y, Shen W. Identifying mild hepatic encephalopathy based on multi-layer modular algorithm and machine learning. Front Neurosci. 2020;14: 627062.
Ong JP, Aggarwal A, Krieger D, Easley KA, Karafa MT, Van Lente F, et al. Correlation between ammonia levels and the severity of hepatic encephalopathy. Am J Med. 2003;114(3):188–93.
Jaffe A, Lim JK, Jakab SS. Pathophysiology of Hepatic Encephalopathy. Clin Liver Dis. 2020;24(2):175–88.
Wang K, Lu Y, Zhao Z, Zhang C. Bioinformatics-based analysis of lncRNA-mRNA interaction network of mild hepatic encephalopathy in cirrhosis. Comput Math Methods Med. 2021;2021:7777699.
Seyan AS, Hughes RD, Shawcross DL. Changing face of hepatic encephalopathy: role of inflammation and oxidative stress. World J Gastroenterol. 2010;16(27):3347–57.
Hsu SJ, Zhang C, Jeong J, Lee SI, McConnell M, Utsumi T, et al. Enhanced meningeal lymphatic drainage ameliorates neuroinflammation and hepatic encephalopathy in cirrhotic rats. Gastroenterology. 2021;160(4):1315-29.e13.
Liu R, Kang JD, Sartor RB, Sikaroodi M, Fagan A, Gavis EA, et al. Neuroinflammation in murine cirrhosis is dependent on the gut microbiome and is attenuated by fecal transplant. Hepatology (Baltimore, MD). 2020;71(2):611–26.
Patel VC, Lee S, McPhail MJW, Da Silva K, Guilly S, Zamalloa A, et al. Rifaximin-α reduces gut-derived inflammation and mucin degradation in cirrhosis and encephalopathy: RIFSYS randomised controlled trial. J Hepatol. 2022;76(2):332–42.
Lu B, Wu C, Azami NLB, Xie D, Zhao C, Xu W, et al. Babao Dan improves neurocognitive function by inhibiting inflammation in clinical minimal hepatic encephalopathy. Biomed Pharmacother. 2021;135:111084.
Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinform. 2008;9:559.
Zemtsova I, Görg B, Keitel V, Bidmon HJ, Schrör K, Häussinger D. Microglia activation in hepatic encephalopathy in rats and humans. Hepatology. 2011;54(1):204–15.
Braissant O, Rackayová V, Pierzchala K, Grosse J, McLin VA, Cudalbu C. Longitudinal neurometabolic changes in the hippocampus of a rat model of chronic hepatic encephalopathy. J Hepatol. 2019;71(3):505–15.
Liddelow SA, Guttenplan KA, Clarke LE, Bennett FC, Bohlen CJ, Schirmer L, et al. Neurotoxic reactive astrocytes are induced by activated microglia. Nature. 2017;541(7638):481–7.
Pan YB, Sun Y, Li HJ, Zhou LY, Zhang J, Feng DF. Transcriptome analyses reveal systematic molecular pathology after optic nerve crush. Front Cell Neurosci. 2021;15: 800154.
Butterworth RF. Hepatic encephalopathy in cirrhosis: pathology and pathophysiology. Drugs. 2019;79(Suppl 1):17–21.
Butterworth RF. The liver-brain axis in liver failure: neuroinflammation and encephalopathy. Nat Rev Gastroenterol Hepatol. 2013;10(9):522–8.
Balzano T, Dadsetan S, Forteza J, Cabrera-Pastor A, Taoro-Gonzalez L, Malaguarnera M, et al. Chronic hyperammonemia induces peripheral inflammation that leads to cognitive impairment in rats: Reversed by anti-TNF-α treatment. J Hepatol. 2020;73(3):582–92.
Wasmuth HE, Kunz D, Yagmur E, Timmer-Stranghöner A, Vidacek D, Siewert E, et al. Patients with acute on chronic liver failure display “sepsis-like” immune paralysis. J Hepatol. 2005;42(2):195–201.
Lin CY, Tsai IF, Ho YP, Huang CT, Lin YC, Lin CJ, et al. Endotoxemia contributes to the immune paralysis in patients with cirrhosis. J Hepatol. 2007;46(5):816–26.
Matsui A, Mochida S, Ohno A, Nagoshi S, Hirose T, Fujiwara K. Plasma osteopontin levels in patients with fulminant hepatitis. Hepatol Res. 2004;29(4):202–6.
Ugamura A, Chu PS, Nakamoto N, Taniki N, Ojiro K, Hibi T, et al. Liver fibrosis markers improve prediction of outcome in non-acetaminophen-associated acute liver failure. Hepatology communications. 2018;2(11):1331–43.
Loffredo L, Pastori D, Farcomeni A, Violi F. Effects of anticoagulants in patients with cirrhosis and portal vein thrombosis: a systematic review and meta-analysis. Gastroenterology. 2017;153(2):480-7.e1.
Drolz A, Schramm C, Seiz O, Groth S, Vettorazzi E, Horvatits T, et al. Risk factors associated with bleeding after prophylactic endoscopic variceal ligation in cirrhosis. Endoscopy. 2021;53(3):226–34.
Frazão JB, Thain A, Zhu Z, Luengo M, Condino-Neto A, Newburger PE. Regulation of CYBB gene expression in human phagocytes by a distant upstream NF-κB binding site. J Cell Biochem. 2015;116(9):2008–17.
Weaver CJ, Terzi A, Roeder H, Gurol T, Deng Q, Leung YF, et al. nox2/cybb deficiency affects zebrafish retinotectal connectivity. J Neurosci. 2018;38(26):5854–71.
Keller CW, Kotur MB, Mundt S, Dokalis N, Ligeon LA, Shah AM, et al. CYBB/NOX2 in conventional DCs controls T cell encephalitogenicity during neuroinflammation. Autophagy. 2021;17(5):1244–58.
Mighiu AS, Recalde A, Ziberna K, Carnicer R, Tomek J, Bub G, et al. Inducibility, but not stability, of atrial fibrillation is increased by NOX2 overexpression in mice. Cardiovasc Res. 2021;117(11):2354–64.
Reinehr R, Görg B, Becker S, Qvartskhava N, Bidmon HJ, Selbach O, et al. Hypoosmotic swelling and ammonia increase oxidative stress by NADPH oxidase in cultured astrocytes and vital brain slices. Glia. 2007;55(7):758–71.
Bobermin LD, Souza DO, Gonçalves CA, Quincozes-Santos A. Resveratrol prevents ammonia-induced mitochondrial dysfunction and cellular redox imbalance in C6 astroglial cells. Nutr Neurosci. 2018;21(4):276–85.
Balasubramaniyan V, Wright G, Sharma V, Davies NA, Sharifi Y, Habtesion A, et al. Ammonia reduction with ornithine phenylacetate restores brain eNOS activity via the DDAH-ADMA pathway in bile duct-ligated cirrhotic rats. Am J Physiol Gastrointest Liver Physiol. 2012;302(1):G145–52.
Hou L, Zhou X, Zhang C, Wang K, Liu X, Che Y, et al. NADPH oxidase-derived H(2)O(2) mediates the regulatory effects of microglia on astrogliosis in experimental models of Parkinson’s disease. Redox Biol. 2017;12:162–70.
Ravelli KG, Santos GD, Dos Santos NB, Munhoz CD, Azzi-Nogueira D, Campos AC, et al. Nox2-dependent Neuroinflammation in An EAE Model of Multiple Sclerosis. Transl Neurosci. 2019;10:1–9.
Hu CF, Wu SP, Lin GJ, Shieh CC, Hsu CS, Chen JW, et al. Microglial Nox2 plays a key role in the pathogenesis of experimental autoimmune encephalomyelitis. Front Immunol. 2021;12: 638381.
Simpson DSA, Oliver PL. ROS generation in microglia: understanding oxidative stress and inflammation in neurodegenerative disease. Antioxidants. 2020;9:8.
Surace MJ, Block ML. Targeting microglia-mediated neurotoxicity: the potential of NOX2 inhibitors. Cell Mol Life Sci. 2012;69(14):2409–27.
Gage MC, Thippeswamy T. Inhibitors of Src family kinases, inducible nitric oxide synthase, and nadph oxidase as potential CNS drug targets for neurological diseases. CNS Drugs. 2021;35(1):1–20.
Wang Q, Chu CH, Oyarzabal E, Jiang L, Chen SH, Wilson B, et al. Subpicomolar diphenyleneiodonium inhibits microglial NADPH oxidase with high specificity and shows great potential as a therapeutic agent for neurodegenerative diseases. Glia. 2014;62(12):2034–43.
Massey N, Puttachary S, Bhat SM, Kanthasamy AG, Charavaryamath C. HMGB1-RAGE Signaling Plays a Role in Organic Dust-Induced Microglial Activation and Neuroinflammation. Toxicol Sci. 2019;169(2):579–92.
Geng L, Fan LM, Liu F, Smith C, Li J. Nox2 dependent redox-regulation of microglial response to amyloid-β stimulation and microgliosis in aging. Sci Rep. 2020;10(1):1582.
Xing YQ, Li A, Yang Y, Li XX, Zhang LN, Guo HC. The regulation of FOXO1 and its role in disease progression. Life Sci. 2018;193:124–31.
Kim KW, Donato J Jr, Berglund ED, Choi YH, Kohno D, Elias CF, et al. FOXO1 in the ventromedial hypothalamus regulates energy balance. J Clin Investig. 2012;122(7):2578–89.
Deng S, Jin P, Sherchan P, Liu S, Cui Y, Huang L, et al. Recombinant CCL17-dependent CCR4 activation alleviates neuroinflammation and neuronal apoptosis through the PI3K/AKT/Foxo1 signaling pathway after ICH in mice. J Neuroinflammation. 2021;18(1):62.
Sadagurski M, Cheng Z, Rozzo A, Palazzolo I, Kelley GR, Dong X, et al. IRS2 increases mitochondrial dysfunction and oxidative stress in a mouse model of Huntington disease. J Clin Investig. 2011;121(10):4070–81.
Pan CW, Jin X, Zhao Y, Pan Y, Yang J, Karnes RJ, et al. AKT-phosphorylated FOXO1 suppresses ERK activation and chemoresistance by disrupting IQGAP1-MAPK interaction. EMBO J. 2017;36(8):995–1010.
Koh PO. Melatonin prevents the injury-induced decline of Akt/forkhead transcription factors phosphorylation. J Pineal Res. 2008;45(2):199–203.
Xu S, Wang J, Zhong J, Shao M, Jiang J, Song J, et al. CD73 alleviates GSDMD-mediated microglia pyroptosis in spinal cord injury through PI3K/AKT/Foxo1 signaling. Clin Transl Med. 2021;11(1): e269.
Ji L, Wang Q, Huang F, An T, Guo F, Zhao Y, et al. FOXO1 overexpression attenuates tubulointerstitial fibrosis and apoptosis in diabetic kidneys by ameliorating oxidative injury via TXNIP-TRX. Oxid Med Cell Longev. 2019;2019:3286928.
Li C, Sheng M, Lin Y, Xu D, Tian Y, Zhan Y, et al. Functional crosstalk between myeloid Foxo1-β-catenin axis and Hedgehog/Gli1 signaling in oxidative stress response. Cell Death Differ. 2021;28(5):1705–19.
Kyriazis ID, Hoffman M, Gaignebet L, Lucchese AM, Markopoulou E, Palioura D, et al. KLF5 is induced by FOXO1 and causes oxidative stress and diabetic cardiomyopathy. Circ Res. 2021;128(3):335–57.
Doan KV, Kinyua AW, Yang DJ, Ko CM, Moh SH, Shong KE, et al. FoxO1 in dopaminergic neurons regulates energy homeostasis and targets tyrosine hydroxylase. Nat Commun. 2016;7:12733.
Peng S, Li W, Hou N, Huang N. A Review of FoxO1-Regulated Metabolic Diseases and Related Drug Discoveries. Cells. 2020;9:1.
Rama Rao KV, Jayakumar AR, Norenberg MD. Role of oxidative stress in the ammonia-induced mitochondrial permeability transition in cultured astrocytes. Neurochem Int. 2005;47(1–2):31–8.
Sinke AP, Jayakumar AR, Panickar KS, Moriyama M, Reddy PV, Norenberg MD. NFkappaB in the mechanism of ammonia-induced astrocyte swelling in culture. J Neurochem. 2008;106(6):2302–11.
Mustafa HN, El Awdan SA, Hegazy GA. Protective role of antioxidants on thioacetamide-induced acute hepatic encephalopathy: biochemical and ultrastructural study. Tissue Cell. 2013;45(5):350–62.
Mostafa RE, Salama AAA, Abdel-Rahman RF, Ogaly HA. Hepato- and neuro-protective influences of biopropolis on thioacetamide-induced acute hepatic encephalopathy in rats. Can J Physiol Pharmacol. 2017;95(5):539–47.
Mousa N, Abdel-Razik A, Zaher A, Hamed M, Shiha G, Effat N, et al. The role of antioxidants and zinc in minimal hepatic encephalopathy: a randomized trial. Ther Adv Gastroenterol. 2016;9(5):684–91.
Shinagawa K, Yanada M, Sakura T, Ueda Y, Sawa M, Miyatake J, et al. Tamibarotene as maintenance therapy for acute promyelocytic leukemia: results from a randomized controlled trial. J Clin Oncol. 2014;32(33):3729–35.
Katsuki H, Kurimoto E, Takemori S, Kurauchi Y, Hisatsune A, Isohama Y, et al. Retinoic acid receptor stimulation protects midbrain dopaminergic neurons from inflammatory degeneration via BDNF-mediated signaling. J Neurochem. 2009;110(2):707–18.
Matsushita H, Hijioka M, Hisatsune A, Isohama Y, Shudo K, Katsuki H. A retinoic acid receptor agonist Am 80 rescues neurons, attenuates inflammatory reactions, and improves behavioral recovery after intracerebral hemorrhage in mice. J Cerebr Blood Flow Metab. 2011;31(1):222–34.
Fukasawa H, Nakagomi M, Yamagata N, Katsuki H, Kawahara K, Kitaoka K, et al. Tamibarotene: a candidate retinoid drug for Alzheimer’s disease. Biol Pharm Bull. 2012;35(8):1206–12.
García-Laencina PJ, Abreu PH, Abreu MH, Afonoso N. Missing data imputation on the 5-year survival prediction of breast cancer patients with unknown discrete values. Comput Biol Med. 2015;59:125–33.
Leek JT, Johnson WE, Parker HS, Jaffe AE, Storey JD. The sva package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics. 2012;28(6):882–3.
Zhang B, Horvath S. A general framework for weighted gene co-expression network analysis. Stat Appl Genet Mol Biol. 2005;4:17.
Long J, Huang S, Bai Y, Mao J, Wang A, Lin Y, et al. Transcriptional landscape of cholangiocarcinoma revealed by weighted gene coexpression network analysis. Briefings Bioinform. 2021;22:4.
Yoo M, Shin J, Kim J, Ryall KA, Lee K, Lee S, et al. DSigDB: drug signatures database for gene set analysis. Bioinformatics (Oxford, England). 2015;31(18):3069–71.
Chen EY, Tan CM, Kou Y, Duan Q, Wang Z, Meirelles GV, et al. Enrichr: interactive and collaborative HTML5 gene list enrichment analysis tool. BMC Bioinform. 2013;14:128.
Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.
Robin X, Turck N, Hainard A, Tiberti N, Lisacek F, Sanchez JC, et al. pROC: an open-source package for R and S+ to analyze and compare ROC curves. BMC Bioinformatics. 2011;12:77.
Acknowledgements
The authors gratefully acknowledge the GEO database available to us for this study.
Funding
This study was supported by National Natural Science Foundation of China (Nos. 81900550, 81970526, 82200630, 81770599), the Clinical Research Program of Ninth People's Hospital, Shanghai Jiao Tong University School of Medicine (JYLJ021).
Author information
Authors and Affiliations
Contributions
CH Zhang and GQ Luo analyzed the data and wrote the article. JY Lin participated in language polish and data check. ZF Zhao was responsible for data visualization. M Luo and HJ Li designed the study. All authors read and approved the final manuscript.
Corresponding authors
Ethics declarations
Ethics approval and informed consent
This study was approved by the ethics committee of Shanghai Ninth People’s Hospital.
Competing interests
The authors declare no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1.
Additional figures.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Zhang, C., Luo, G., Lin, J. et al. Identification of significant modules and hub genes involved in hepatic encephalopathy using WGCNA. Eur J Med Res 27, 264 (2022). https://doi.org/10.1186/s40001-022-00898-3
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s40001-022-00898-3