Employment of digital gene expression profiling to identify potential pathogenic and therapeutic targets of fulminant hepatic failure
© Chen et al.; licensee BioMed Central. 2015
Received: 24 October 2014
Accepted: 5 January 2015
Published: 27 January 2015
The dysregulated cytokine metabolism and activity are crucial to the development of fulminant hepatic failure (FHF), and many different cytokines have been identified. However, the precise gene expression profile and their interactions association with FHF are yet to be further elucidated.
In this study, we detected the digital gene expression profile (DGEP) by high-throughput sequencing in normal and FHF mouse liver, and the candidate genes and potential targets for FHF therapy were verified. And the FHF mouse model was induced by D-Galactosamine (GalN)/lipopolysaccharide (LPS).
Totally 12727 genes were detected, and 3551 differentially expressed genes (DEGs) were obtained from RNA-seq data in FHF mouse liver. In FHF mouse liver, many of those DEGs were identified as differentially expressed in metabolic process, biosynthetic process, response to stimulus and response to stress, etc. Similarly, pathway enrichment analysis in FHF mouse liver showed that many significantly DEGs were also enriched in metabolic pathways, apoptosis, chemokine signaling pathways, etc. Considering the important role of nuclear factor-kappa B (NF-κB) in metabolic regulation and delicate balance between cell survival and death, several DEGs involved in NF-κB pathway were selected for experimental validation. As compared to normal control, NF-κBp65 and its inhibitory protein IκBα were both significantly increased, and NF-κB targeted genes including tumor necrosis factor α(TNFα), inducible nitric oxide synthase (iNOS), interleukin-1β, chemokines CCL3 and CCL4 were also increased in hepatic tissues of FHF. In addition, after NF-κB was successfully pre-blocked, there were significant alteration of hepatic pathological damage and mortality of FHF mouse model.
This study provides the globe gene expression profile of FHF mouse liver, and demonstrates the possibility of NF-κB gene as a potential therapeutic target for FHF.
KeywordsFulminant hepatic failure Differentially expressed genes Nuclear factor-kappa B Therapeutic target
Fulminant hepatic failure (FHF) is a severe, life-threatening disorder, which is characterized by a dramatic clinical syndrome resulting from massive hepatocyte death [1,2]. There is a large number of evidence showing that the dysregulated cytokine metabolism and activity are crucial to the progression of liver damage and the regulation of undesired inflammation is an attractive potential strategy in the treatment of FHF [3-5]. Though many different cytokines have been identified, the precise gene expression profile and their interactions in FHF are yet to be further elucidated.
In past decade, expression microarrays have been applied for genome-wide expression analysis for identification of differentially expressed genes (DEGs) associated with liver injury . Although these studies provided new information regarding gene expression patterns in FHF, interactions between genes and affected pathways may have been missed because of technical limitations. Recently, by the rapid development of a “second generation” of sequencing technologies, the digital gene expression profile (DGEP) method was made [7-9], and millions of expression tags can be measured simultaneously at a fraction of the cost of capillary sequencing used in the original serial analysis of gene expression method. A major advantage of DGEP, which overcomes the drawbacks of using predesigned hybridization probes as in microarray-based methods, is that it allows detection of expressed genes in a strand-specific manner; and DGEP method provides information on the polarity of the expressed transcripts, allowing detection of sense and antisense transcripts from the same gene region.
At present, DGEP experiments often involve measuring the relative amount of mRNA expressed in two or more experimental conditions. This is because altered levels of a specific sequence of mRNA suggest a changed need for the protein coded for by the mRNA, perhaps indicating a homeostatic response or a pathological condition . For example, if injured hepatic cells express higher levels of mRNA associated with a particular cytokine than normal cells do, it might be that this cytokine plays a role in hepatic injuries. Thus, the application of DGEP in transcriptome profiling of severe injured liver tissue is important to further understand the underlying mechanisms of FHF.
D-Galactosamine (GalN)/lipopolysaccharide (LPS)-induced liver failure in mice is dependent on macrophage-derived pro-inflammatory cytokines, including tumor necrosis factor (TNF)-α, interleukin (IL)-1β, IL-6 and interferon-γ, and accurately represents human FHF . Thus, in this study, we detected global gene expression profile in normal and FHF mouse liver by DGEP method, and showed the main FHF-associated hepatic biological processes, pathways, and key molecules potentially serving as important therapeutic target for FHF.
Materials and methods
Male BALB/c mice, 18–22 g, were obtained from Sichuan University (Chengdu, China). They were housed and cared for in rooms maintained at a constant temperature and humidity. Food and water were allowed ad libitum. Food was withdrawn the evening before the experiment. In this study, animal experiments involved in screening of DEGs in FHF mouse liver (Part. 1), and verification of some potential candidate genes related to the progression of FHF (Part. 2). In part. 1, a mouse model of FHF was established by intraperitoneal injections of D-GalN (800 mg/kg body weight; Sigma, Saint Louis, USA) and LPS (20 μg/kg body weight; Sigma), as described previously. Control mice were intraperitoneally injected with saline. Plasma and liver tissue samples were obtained 6 h after injection and were stored at -80°C until analysis. In part. 2, a potent NF-κB chemical inhibitor pyrrolidine dithiocarbamate (PDTC) was given by intraperitoneal injection to pre-block of NF-κB activation and expression, 12 h before D-GalN/LPS administration; and the plasma and liver tissue samples were also obtained 6 h after the construction of FHF mouse liver, and were stored at -80°C until analysis.
Sample preparation and solexa library construction
Liver samples were ground in liquid nitrogen, and total RNA was extracted with Trizol reagent (Invitrogen). The integrity of total RNA was evaluated using an Agilent Bioanalyzer 2100. Solexa libraries were constructed following the manufacturer’s standard according to the schematic as described previously.
Solexa sequencing and data analysis
The image files obtained from Illumina 1G sequencer were processed to produce sequence data. Then the high-quality reads were screened from the original raw data, and the adaptors were removed from each sequence. Subsequently, high quality clean tags were compared with RefSeq database and the expression level of each gene was normalized to transcripts per million (TPM). The significance of DGEP was analyzed as described previously. The expression pattern of these DEGs was visualized using the heat-map function in the R base package. In this study, to avoid the possible noise signal from high-throughput sequencing, the genes with an average TPM less than 1 were excluded; and the remained genes were used to calculate the fold changes and false discovery rate (FDR). And the absolute fold change no less than 1(|log2Ratio| ≥ 1) and FDR less than 0.001 were used to define the DEGs.
Gene ontology functional and pathway enrichment analysis for DEGs
The clustered genes were assigned to molecular function, cellular component and biological process based on Gene Ontology using the web tool DAVID (http://david.abcc.ncifcrf.gov/home.jsp). Hypergeometric test was used to select the enriched molecular function, cellular component and biological process in Gene Ontology for each cluster. The pathways associated with these gene clusters were analyzed by KEGG (www.genome.jp/kegg/pathway.html).
Detection of NF-κBp65 and its inhibitory protein IκBα in liver tissue
The protein expression levels of NF-κBp65 and its inhibitory protein IκBα in liver tissue were determined by western blot analysis, according to the standard manufacturer’s protocol. And the rabbit polyclonal anti-NF-κBp65 (Abcam, England) and anti-IκBα antibodies (Abcam, England) were both diluted at 1: 200; while rabbit polyclonal anti-GAPDH antibody (ProteinTech, USA) was diluted at 1:1000 as a control. Immunodetection was performed with the ECL-Plus kit (Pierce Biotechnology, USA), and immunoblot signals were quantified using Quantity One Software (Bio-Rad).
Biochemistry and plasma protein analysis
In order to assess the extent of liver injury, serum alanine aminotransferase (ALT), aspertate aminotransferase (AST), total bilirubin (TBil) and indirect bilirubin (IBil) were measured with automatic analyzer (HITACHI, Japan). Serum levels of TNF-α were measured by enzyme-linked immunosorbent assay (GenStar BioSolutions, China), with lower limits of detection of 7 pg/mL.
Determination of mRNA levels for candidate genes in liver tissue
Primer sequences used in the quantitative RT-PCR analysis
Liver tissue sections (paraffin embedded) were stained with hematoxylin-eosin using a standard protocol, and analyzed by light microscopy under 400-fold magnification. Morphological criteria including vacuolization, swollen cytoplasm with disrupted cell and organelle membranes, and lytic nuclear changes served to determine necrosis. Hepatocyte apoptosis was assessed by the TUNEL assay kit (Roche, Switzerland). The protein expression of TNFα and iNOS in liver tissue was detected by immunohistochemistry; and the percentage of positive cells and the positive staining intensity were checked and scored using the Axiotis score standard by observation of 5 randomly chosen fields at 400-fold magnification . The primary anti-TNFα and anti-iNOS antibodies were purchased from Cell Signaling Technology (USA). In this study, the liver histology was assessed by a single pathologist, unaware of the laboratory data.
All animal experiments were done in triplicate, and the results were expressed as mean ± SD. Statistical comparisons were made by 1-way analysis of variance using SPSS 17.0 (SPSS Inc, Chicago, IL). The differences were considered statistically significant if P value less than 0.05.
High-throughput sequencing data and DEGs in FHF mouse liver
The main biological processes and pathways altered in FHF mouse liver
A partial list of pathways significantly altered in FHF mouse liver
Protein processing in endoplasmic reticulum
Drug metabolism - other enzymes
Chemokine signaling pathway
Adipocytokine signaling pathway
Toll-like receptor signaling pathway
Pantothenate and CoA biosynthesis
Fructose and mannose metabolism
Ascorbate and aldarate metabolism
The DEGs regulated by nuclear factor κB in FHF mouse liver
Pre-block of NF-κB significant alleviates hepatic pathological damage and mortality of FHF mouse
In immunohistochemistry assay, the number of positive cells and intensity of the reaction of TNFα and iNOS were obvious in liver tissue of FHF (Figure 6D); and after NF-κB was pre-blocked, the expression of TNFα and iNOS in liver tissue were both significantly decreased (Figure 6D,E and F). In summary, our data suggested that suppression of NF-κB activation or expression mediated by chemical inhibitor PDTC could significantly improve liver function and reduce mortality of FHF mouse model, and NF-κB would be a potential new therapeutic target for FHF.
At present, the exact mechanisms of the pathophysiological events in the liver leading to FHF are not fully understood, there is still no available particularly effective therapy except for liver transplantation; however, the latter is also limited by a severe shortage of donor livers . Thus, there is an urgent need to explore new effective treatment measures. In past decade, some host genes or their encoded moleculars (including cytokines and chemokines) had been discovered that were associated with hepatocyte unnatural death and participated in the occurrence of FHF . However, the interactions of vast majority of cytokines/chemokines with each other remained unclear in the process of FHF.
Currently, detecting changes in gene expression of host and their interaction networks have become an increasingly important topic in FHF, and the genes differentially expressed in FHF are of foremost interest to both basic and clinical researches . With the assistance of the new-generation high throughput sequencing technology, DGE has be employed to economically and quickly capture the whole genome of the gene expression of certain tissue of a species under specific conditions . In this study, this technology was firstly used to identify DEGs of mouse liver tissue under different external stimuli conditions, and it was also the first reported attempt to reveal the global network of DEGs in FHF. It is worth to mention that our results have provided a variety of potential candidate targets for pathogenic mechanism and new drug development of FHF, and all raw data could be obtained in Additional file 1: Table S1.
Though hepatotropic virus infection is the main cause of FHF in clinical practice, there is still no hepatotropic virus induced FHF animal model available at present, which greatly limits the investigations of related molecular mechanisms and new treatment strategies of FHF. As we know, the activation of hepatic macrophages and its-associated inflammation are extremely common in patients with FHF, and this situation is largely related to endotoxin exposure . So, most scholars had to choose endotoxin-induced animal model to conduct related research work on FHF in past years, though it has defects in reflecting the role of hepatotropic virus in the pathogenesis of massive hepatic necrosis. Based on this present situation, the D-GalN/LPS-induced acute liver damage mouse model has been used to yield the global gene expression profile of FHF in present study; and a total of 3551 DEGs were successfully detected in liver tissue of FHF, including 1366 up-regulated and 2185 down-regulated DEGs. And majority of those detected DGEs were strongly enriched metabolic process, biosynthetic process, response to stimulus and response to stress. As we know, liver is a multifunctional organ involved in important metabolic functions, synthesis of plasma proteins and storage of lipids and retinol, and the failure of liver function could result in serious metabolic disorder of sugar, fat and protein. Thus the highly enriched metabolic process observed in this study was consistent with the objective fact. Although we did not conduct in-depth verification of the relevant DEGs on liver cellular metabolism, the global interaction map of corresponding DEGs reported in this study would provide important clues to future studies.
Excessive apoptosis and abnormal necrosis are prominent features of FHF . According to analysis of biological process, high to 579 and 365 DEGs involved in response to stimulus and response to stress, respectively. To shed more lights into the function roles of genes responsible for the progression of FHF, we further investigated the hepatocyte damage associated apoptosis and chemokine signaling pathways by the enrichment analysis of DEGs among normal and FHF liver samples, and found there were 39 DEGs in apoptosis pathway and 76 DEGs in chemokine signaling pathways. From the large number of DEGs in above two-mentioned pathways, it’s easy to see that inflammatory processes play a fundamental role in liver tissue damage, and the complexity of inflammatory cytokine network involved in the occurrence and development of FHF. In recent years, many in vitro and in vivo studies have examined the contribution of components of the NF-κB signaling pathways to the regulation of liver cell damage and repair, in particular, of liver cell death and regeneration in liver damage development [19,20]. However, the exact role of NF-κB in FHF in vivo is still controversial . In this study, according to the clues provided by DGEP experiment, we performed a cluster analysis of DEGs in NF-κB pathway, and there were 32 DEGs up-regulated or down-regulated in FHF. Additionally, the expression levels of NF-κBp65 and IκBα proteins were both increased significantly in liver tissue. So, we inferred that the activation of NF-κB mediated signaling pathways played important roles in the development of FHF. As the gene transcription of IκBα was regulated by NF-κB, the high expression of IκBα proteins in FHF could be explained by the highly activated NF-κBp65. In addition, the high expression of IκBα proteins in this experiment may also execute a negative feedback regulation of NF-κB overexpression in the development of FHF .
In this study, after the activation approach of NF-κB were pre-blocked, we found that D-GalN/ LPS intraperitoneal injection could not induce significantly liver damage and all mouse survived. So we believed that the overexpression of NF-κB or highly activated NF-κB should participate in the progression of liver damage to FHF.
Though the higher expression levels of NF-κB have been reported to protect cells from apoptosis and necrosis in some vitro experiments , the overexpression of NF-κB also would activate and enhance the host non-specific and specific immune response, and therefore result in tissue injury and organ dysfunction . As we know, the biological function of certain molecules would be affected by the changes of surrounding environment of cells, so the vast difference between in vitro and in vivo environment is likely to influence the functions of NF-κB. So, we believe that the overexpression of NF-κB should mediate or aggravate the liver cell necrosis in the process of FHF. Considering that oxidative stress also has been shown to aggravate the hepatotoxicity induced by LPS/D-GalN and PDTC is also an antioxidant, BAY-11-7082 (another inhibitor of NF-κB) also was used to verify the function of NF-κB, and the result was in good agreement with PDTC experiment (data unshown).
Previous studies showed NF-κB was able to bind promoter and enhancer regions containing κB sites, and TNFα, iNOS, IL-2, IL-6, IL-8, GM-CSF, ICAM-1, MHC-1, TNFβ, Toll -like receptors, chemokines and chemokine receptor (CR) were all regulated by NF-κB [25-27]. So in this study, the high-level expression of cytokines (such as TNFα and iNOS) and chemokines were associated with the over-activated NF-κB complex, and their direct and indirect damaging effects on liver cells also partially contributed to the pathogenic role of NF-κB activation .
This study is the first genome-wide effort to investigate the transcriptional changes in FHF development. Though it is impossible to discuss the relevance and impact of all genes and clusters that were differentially expressed in this study, present study still demonstrates the possibility of NF-κB gene as a potential therapeutic target for FHF.
All animal experimental procedures were approved by the Institutional Animal Care and Use Committee of Sichuan University.
This study was supported by the National Natural Science Foundation of China (No.81300319 and 30972622). We thank Dr. Bing Yang (BGI-Shenzhen, China) for his assistance with pathway enrichment analysis of DEGs.
- Ichai P, Samuel D. Etiology and prognosis of fulminant hepatitis in adults. Liver Transpl. 2008;14 Suppl 2:S67–79.PubMedView ArticleGoogle Scholar
- Graziadei IW. The clinical challenges of acute on chronic liver failure. Liver Int. 2011;31 Suppl 3:24–6.PubMedView ArticleGoogle Scholar
- Liu Q. Role of cytokines in the pathophysiology of acute-on-chronic liver failure. Blood Purif. 2009;28:331–41.PubMedView ArticleGoogle Scholar
- Rutherford A, Chung RT. Acute liver failure: mechanisms of hepatocyte injury and regeneration. Semin Liver Dis. 2008;28:167–74.PubMedView ArticleGoogle Scholar
- Dominik A, Stange J, Pfensig C, Borufka L, Weiss-Reining H, Eggert M. Reduction of elevated cytokine levels in acute/acute-on-chronic liver failure using super-large pore albumin dialysis treatment: an In vitro study. Ther Apher Dial. 2014;18:347–52.PubMedView ArticleGoogle Scholar
- Cobb JP, Laramie JM, Stormo GD, Morrissey JJ, Shannon WD, Qiu Y, et al. Sepsis gene expression profiling: murine splenic compared with hepatic responses determined by using complementary DNA microarrays. Crit Care Med. 2002;30:2711–21.PubMedView ArticleGoogle Scholar
- El-Metwally S, Hamza T, Zakaria M, Helmy M. Next-generation sequence assembly: four stages of data processing and computational challenges. PLoS Comput Biol. 2013;9:e1003345.PubMed CentralPubMedView ArticleGoogle Scholar
- Shendure J, Ji H. Next-generation DNA sequencing. Nat Biotechnol. 2008;26:1135–45.PubMedView ArticleGoogle Scholar
- Wang W, Qin Z, Feng Z, Wang X, Zhang X. Identifying differentially spliced genes from two groups of RNA-seq samples. Gene. 2013;518:164–70.PubMedView ArticleGoogle Scholar
- Huang Q, Lin B, Liu H, Ma X, Mo F, Yu W, et al. RNA-Seq analyses generate comprehensive transcriptomic landscape and reveal complex transcript patterns in hepatocellular carcinoma. PLoS One. 2011;6:e26168.PubMed CentralPubMedView ArticleGoogle Scholar
- Chen EQ, Gong DY, Leng XH, Bai L, Liu C, Wang LC, et al. Inhibiting the expression of hepatocyte nuclear factor 4 alpha attenuates lipopolysaccharide/D-galactosamine-induced fulminant hepatic failure in mice. Hepatobiliary Pancreat Dis Int. 2012;11:624–9.PubMedView ArticleGoogle Scholar
- Axiotis CA, Monteagudo C, Merino MJ, LaPorte N, Neumann RD. Immunohistochemical detection of P-glycoprotein in endometrial adenocarcinoma. Am J Pathol. 1991;138:799–806.PubMed CentralPubMedGoogle Scholar
- Merion RM. Current status and future of liver transplantation. Semin Liver Dis. 2010;30:411–21.PubMedView ArticleGoogle Scholar
- Sun S, Guo Y, Zhao G, Zhou X, Li J, Hu J, et al. Complement and the alternative pathway play an important role in LPS/D-GalN-induced fulminant hepatic failure. PLoS One. 2011;6:e26838.PubMed CentralPubMedView ArticleGoogle Scholar
- Duan ZP, Chin Y. Prospect of severe hepatitis and liver failure. Zhonghua Gan Zang Bing Za Zhi. 2009;17:10–1.PubMedGoogle Scholar
- Su Z, Ning B, Fang H, Hong H, Perkins R, Tong W, et al. Next-generation sequencing and its applications in molecular diagnostics. Expert Rev Mol Diagn. 2011;11:333–43.PubMedGoogle Scholar
- Yang Q, Shi Y, He J, Chen Z. The evolving story of macrophages in acute liver failure. Immunol Lett. 2012;147:1–9.PubMedView ArticleGoogle Scholar
- McGill MR, Jaeschke H. Apoptosis or necrosis in acetaminophen-induced acute liver failure? New insights from mechanistic biomarkers*. Crit Care Med. 2013;41:2653–4.PubMedView ArticleGoogle Scholar
- Papa S, Bubici C, Zazzeroni F, Franzoso G. Mechanisms of liver disease: cross-talk between the NF-kappaB and JNK pathways. Biol Chem. 2009;390:965–76.PubMed CentralPubMedView ArticleGoogle Scholar
- Robinson SM, Mann DA. Role of nuclear factor kappaB in liver health and disease. Clin Sci (Lond). 2010;118:691–705.View ArticleGoogle Scholar
- Wu YL, Lian LH, Wan Y, Nan JX. Baicalein inhibits nuclear factor-kappaB and apoptosis via c-FLIP and MAPK in D-GalN/LPS induced acute liver failure in murine models. Chem Biol Interact. 2010;188:526–34.PubMedView ArticleGoogle Scholar
- Campbell KJ, Perkins ND. Regulation of NF-kappaB function. Biochem Soc Symp. 2006;73:165–80.PubMedGoogle Scholar
- Bharti AC, Aggarwal BB. Nuclear factor-kappa B and cancer: its role in prevention and therapy. Biochem Pharmacol. 2002;64:883–8.PubMedView ArticleGoogle Scholar
- Jayakumar AR, Bethea JR, Tong XY, Gomez J, Norenberg MD. NF-kappaB in the mechanism of brain edema in acute liver failure: studies in transgenic mice. Neurobiol Dis. 2011;41:498–507.PubMed CentralPubMedView ArticleGoogle Scholar
- Erener S, Petrilli V, Kassner I, Minotti R, Castillo R, Santoro R, et al. Inflammasome-activated caspase 7 cleaves PARP1 to enhance the expression of a subset of NF-kappaB target genes. Mol Cell. 2014;46:200–11.View ArticleGoogle Scholar
- Zheng C, Yin Q, Wu H. Structural studies of NF-kappaB signaling. Cell Res. 2011;21:183–95.PubMed CentralPubMedView ArticleGoogle Scholar
- Baltathakis I, Alcantara O, Boldt DH. Expression of different NF-kappaB pathway genes in dendritic cells (DCs) or macrophages assessed by gene expression profiling. J Cell Biochem. 2001;83:281–90.PubMedView ArticleGoogle Scholar
- Brasier AR. The NF-kappaB regulatory network. Cardiovasc Toxicol. 2006;6:111–30.PubMedView ArticleGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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.