Characterization of highly expressed novel hub genes in hepatitis E virus chronicity in rabbits: a bioinformatics and experimental analysis
BMC Veterinary Research volume 18, Article number: 239 (2022)
Hepatitis E virus (HEV), which is the leading cause of acute viral hepatitis worldwide, usually causes self-limited infections in common individuals. However, it can lead to chronic infection in immunocompromised individuals and its mechanisms remain unclear. Rabbits are the natural host of HEV, and chronic HEV infections have been observed in rabbits. Therefore, we aimed to investigate potential key genes in HEV chronicity process in rabbits. In this study, both bioinformatics and experimental analysis were performed to deepen the understanding of hub genes in HEV chronic infection in rabbits.
Ninety-four candidate differentially expressed genes (DEGs) and the pathways they enriched were identified to be related with HEV chronicity. A total of 10 hub genes were found by protein–protein interaction (PPI) network construction. Rabbits of group P (n = 4) which showed symptoms of chronic HEV infection were selected to be compared with HEV negative rabbits (group N, n = 6). By detecting the identified hub genes in groups P and N by real-time PCR, we found that the expressions of MX1, OAS2 and IFI44 were significantly higher in group P (P < 0.05).
In this work, we presented that MX1, OAS2 and IFI44 were significantly upregulated in HEV chronic infected rabbits, indicating that they may be involved in the pathogenesis of HEV chronicity.
Hepatitis E virus (HEV), which is the leading cause of acute viral hepatitis worldwide, can lead to 20 million HEV infections annually . HEV is a positive-sense, single-stranded RNA virus with a 7.2 kb genome and belongs to the Hepeviridae family and the Orthohepevirus genus . HEV has only one serotype, and can be mainly divided into 8 genotypes (HEV1-8) [1, 3], of which HEV1 and 2 only cause human infections and HEV3, 4, 7 and 8 have been recognized as zoonotic viruses [1, 4].
HEV usually causes self-limited HEV infections in common individuals . However, HEV can establish chronic infections in immunocompromised individuals, including patients with solid organ transplant (SOT) , human immunodeficiency virus (HIV) infection  and hematological malignancies . HEV chronic infections are mostly caused by HEV3 , and few cases have been observed with HEV4  and 7 infection . The specific mechanisms of HEV chronic infection are not well understood now and viral factors and host immunity may be related to HEV chronicity [5, 12]. It was demonstrated that a sustained stimulation of Th2 immunity was observed in chronic HEV infection animal models , and interferon-stimulated genes may be involved in the persistence of an HEV infection . Still, further studies are warranted to clarify the molecular mechanisms underlying HEV chronic infection.
With the development of high-throughput technologies and bioinformatic analysis, the pathogenesis of many complex diseases can be investigated [15,16,17]. Hub genes are crucial genes that have interactions with other genes in gene networks and play an important role in the network, which are defined as the genes with connectivity greater than 10 in the interaction network [18, 19]. However, so far studies regarding gene expression profiling of chronicity of HEV are rare. As the natural host of HEV, rabbit models have been used to investigate the clinical symptoms and pathogenesis of HEV infection in several studies [20, 21]. Therefore, in this study, we aimed to identify differentially expressed genes (DEGs) and potential key genes in HEV chronicity process and further demonstrate these genes in animal models.
Sample information processing and screening and identification of candidate differentially expressed genes (DEGs)
The results of ‘R-limma’ analysis showed that there were 94 DEGs, including 80 upregulated DEGs and 14 downregulated DEGs in HEV chronic infected samples compared with normal samples (Fig. 1). The screening criteria for differentially expressed genes were as follows: adjust P-value < 0.05 and log twofold change > 1. Therefore, those DEGs were identified as candidate DEGs in this study.
The enrichment analysis in HEV chronic infected samples
In the enrichment analysis, DAVID and Fun Rich software were used of the samples’ genes. The GO terms of targets of DEGs were shown in Fig. 2A. The DEGs enriched in the BP terms were mostly associated with type I interferon signaling pathway, defense response to virus, negative regulation of viral genome replication, response to virus and interferon-gamma-mediated signaling pathway; those in the CC category mostly included the cytosol, cytoplasm and mitochondrion; and those belonged to MF category were mainly related to 2’-5’-oligoadenylate synthetase activity, nucleotidyltransferase activity and transferase activity. The KEGG pathway enrichment analysis revealed that the targets of DEGs were enriched in pathways in hepatitis C, influenza A, measles 7, herpes simplex infection and RIG-I-like receptor signaling pathway. The top 5 results of the KEGG pathway enrichment analysis are shown in Fig. 2B.
Then the DEGs were uploaded to Funrich for further enrichment analysis. By the analysis, the results indicated that DEGs in HEV chronic infected samples were mainly enriched in interferon signaling pathway and immune-related pathways (Fig. 2C).
Protein–protein interaction (PPI) network construction and hub genes identification
PPI network which contained 549 nodes was then constructed by Cytoscape software on the basis of the STRING database to further screen the hub genes among the DEGs (Fig. 3A). MCODE was used to carry out the network data analysis to identify gene clusters, genes in the first gene cluster with the highest score were selected for BP enrichment analysis, and the results showed that these genes in the first gene cluster were mostly involved in defense response to virus and immune response functions. Then, 10 hub genes of HEV chronic infected samples were selected (Fig. 3B).
Pathogenicity of HEV chronic infection in rabbits
The duration of fecal virus shedding of P1, P2, P3 and P4 were 14th, 17th, 15th and 20th weeks, respectively. The fecal viral load of group P was 104–108 (Fig. 4A). Serum HEV antigen of P1, P2, P3 and P4 was continuously positive before 9th, 10th, 9th and 13th weeks, respectively (Fig. 4B). The anti-HEV antibody of P1, P2, P3 and P4 showed seroconversion for anti-HEV antibody at 11th, 13th, 11th and 10th week, respectively (Fig. 4C). Group N, the negative control, showed no symptoms of HEV infection during the study.
Identification of hub genes in chronic HEV infected rabbits
The rabbits of group P were euthanized 2 weeks after their virus fecal shedding ceased and the rabbits of group N were euthanized at the end of study. Their livers were collected and the RNA was extracted as described above. The hub genes identified in this study were quantified and compared between two groups (Fig. 5). The results showed that MX1, OAS2 and IFI44 were significantly up-regulated in HEV chronic infected rabbits (P < 0.05). Also, the expressions of these genes in positive dropouts (n = 7) were assessed (Supplementary Fig. 1), and the expressions of MX1, ISG15, RSAD2, OAS2, IFI44, IFI44L and DDX60 had significant differences in positive dropouts (P < 0.05).
So far, the validation of large independent datasets for HEV chronicity still remains rare. In this study, we aimed to identify potential key genes by the way of a series of bioinformatics analyses and experimental validation.
The bioinformatics analysis results showed that inflammatory- or immune response-related genes and pathways were involved in HEV chronicity. The enrichment analysis showed that DEGs were significantly enriched in the categories of cellular responses to interferon-β. KEGG pathway analysis showed that DEGs were significantly enriched in RIG-I-like receptor signaling pathway. Both innate and adaptive immune responses seem to play a robust role in HEV chronicity (. Studies have shown that chronic hepatitis E (CHE) patients had lower levels of IL-1 receptor agonist and TNF-α concentration , and their total counts of CD2, CD3 and CD4 T cells were also significantly lower .
Specifically, 10 potential hub genes were selected in this study, including MX1, ISG15, RSAD2, OAS2, IFIT1, IFIT2, IFIT3, IFI44, IFI44L and DDX60. Some of these interferon-stimulated genes (ISGs) have been shown to be related with HEV infection. To verify the expressions of these genes in HEV chronically infected rabbits, we performed the animal experiments. Chronic HEV infection was observed in rabbits of group P, which was manifested by continuous fecal virus shedding for more than 12 weeks, the HEV antigen positivity and seroconversion for anti-HEV antibody. By real-time RT-PCR, we found that the levels of MX1, OAS2 and IFI44 were significantly higher in group P. However, the expressions of remaining genes didn’t show significant difference. Therefore, these results indicated that MX1, OAS2 and IFI44may play key roles in the pathogenesis of HEV chronic infection of rabbits. MX1 has been reported to inhibit the multiplication of many RNA and DNA viruses and were found to be upregulated in HEV infected patients [23, 24]. OAS2 is an interferon-stimulated antiviral enzyme and may be involved in HEV infection . FI44 is an interferon-alfa inducible protein and is related with infection of several viruses, including hepatitis C virus (HCV), human papillomavirus (HPV) and human immunodeficiency virus (HIV) [26,27,28]. Our results showed that these gene was upregulated during the HEV chronic infection, indicating their potential role in HEV chronicity of rabbits.
In conclusion, a total of 94 DEGs and 10 hub genes were identified by bioinformatics analysis. These DEGs in HEV chronic infected samples were mainly enriched in interferon signaling pathway and immune-related pathways. By experimental reserach, the present study demonstrated that hub genes, including MX1, OAS2 and IFI44, were involved in the pathogenesis of chronic HEV infection, indicating that these genes have the potential to be biomarkers in HEV chronicity. Further studies to confirm the mechanisms of these key genes and pathways are needed.
Acquisition of microarray data
The microarray data obtained from the gene expression profiling for chronic HEV infection were downloaded from the Gene Expression Omnibus database (https://www.ncbi.nlm.nih.gov/gds/). The dataset GSE36539, based on the GPL6480 platform, consists of the results from 10 whole blood of patients with chronic hepatitis E infection and 6 whole blood of control patients.
Screening of candidate DEGs
To screen differentially expressed DEGs between the whole blood of chronic HEV infection patients and control patients, a method named “R-limma” package of Bioconductor (http://www.bioconductor.org/) and software R (version 3.4.4, https://www.r-project.org/) was introduced. Probe sets with a P value < 0.05 and |log 2 FC|> 1 were between two groups were identified as significant DEG sets.
Functional annotation and pathway enrichment analysis
In this study, DAVID website (http://david.ncifcrf.gov/home.jsp) and Funrich software (http://www.funrich.org/) were used to evaluate the enrichment of the DEGs into the Cell Components (CC), Molecular Functions (MF), and Biological Processes (BP) of Gene Ontology (GO) terms. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis was performed to analyze relevant signaling pathways of DEGs . P-value < 0.05 was regarded as statistically significant.
Protein–protein interaction (PPI) network construction
The PPI network of DEGs was constructed by STRING database (http://string-db.org). Then, the hub genes in this network were identified by analyzing the degree of connectivity through Cytoscape software (version 3.6.1).
Experimental design of animal study
Fifty 3-month-old Japanese white rabbits (weighing about 2–3 kg) obtained from the Department of Laboratory Animal Science of National Institutes for Food and Drug Control were randomly selected. Fecal and serum specimens of these rabbits were collected weekly for two consecutive weeks. All rabbits with positive fecal HEV RNA (n = 11) were selected for follow-up monitoring. Serum and fecal samples of the rabbits were collected weekly during the follow-up monitoring for HEV RNA, anti-HEV antibody, HEV antigen and alanine transaminase (ALT) level detection. Of all the 11 positive rabbits, 4 had fecal virus shedding for more than 12 weeks, which can be recognized as HEV chronic infection and classified into group P. Also, six 3-month-old Japanese white rabbits (weighing about 2–3 kg) with negative for fecal and serum HEV RNA and anti-HEV antibody were selected and classified into group N as a negative control. The animal experiments were approved by the Committee of Laboratory Animal Welfare and Ethics of National Institutes for Food and Drug Control (2022B012) and First Clinical Medical College, Shanxi Medical University (SYXK2019-0007).
Feces were homogenized in sterile phosphate-buffered saline (PBS) to make 10% fecal suspensions (wt/vol). The suspensions were then clarified by centrifugation. The dissection and collection of tissue samples was conducted following standard protocols. Serum, fecal suspension and tissue samples were stored at − 80 °C for detection.
Detection of anti-HEV, HEV antigen and HEV RNA
All rabbit serum samples were detected by ELISA for total anti-HEV antibody and HEV antigen according to the manufacturer’s instructions (Wantai, Biopharmaceutical, Beijing, China). Total RNA in fecal suspension and tissues was extracted by using QIAamp Viral RNA Mini Kit (Qiagen, Hilden, Germany), according to the protocol of the manufacturer. Real-time fluorescence quantitative PCR (RT-qPCR) in this study was performed by using the Taqman probe detection method and QuantiTect® Probe RT–PCR Kit (Qiagen, Hilden, Germany) .
Detection of hub genes in rabbits
The sequences of rabbit MX1, ISG15, IFIT2, RSAD2, OAS2, IFIT3, IFI44, IFI44L and DDX60 were obtained from GenBank (Accession numbers XM_008252473, XM_017340429, XM_002718372, XM_017350632, XM_002719777, XM_008270166, XM_017346152, XM_008265011, XM_008249585, respectively), as was the sequence for glyceraldehyde-3-phosphate dehydrogenase (GAPDH) (Accession number NM_001082253) (used as a housekeeping gene for normalization). IFIT1 was not detected due to the lack of GenBank sequence. Using these sequences, primers were designed for each of the cytokines (Table 1) for hub gene quantitation by real-time RT-PCR. Real-time PCR assay was performed on an ABI Prism 7500 Sequence Detection System (Applied Biosystems, Foster City, CA, USA) using QuantiTect SYBR Green PCR Kit (Qiagen, Hilden, Germany).
Statistical analysis was performed using the SPSS PASW Statistics v18.0 statistical software package (SPSS, Inc., USA, http://www.ibm.com/cn/). Data were compared using Student’s t-test or chi-squared test. A P-value of < 0.05 was considered significant.
Availability of data and materials
The data shown in this paper are available within the article. The raw data of RNA- sequencing are available from the NCBI database under accession number GSE36539.
Hepatitis E virus
Solid organ transplant
Human immunodeficiency virus
Differentially expressed gene
Gene ontology (GO)
Kyoto Encyclopedia of Genes and Genomes
Chronic hepatitis E
Hepatitis C virus
Nimgaonkar I, Ding Q, Schwartz RE, Ploss A. Hepatitis E virus: advances and challenges. Nat Rev Gastroenterol Hepatol. 2018;15(2):96–110.
Hoofnagle JH, Nelson KE, Purcell RH. Hepatitis E. N Engl J Med. 2012;367(13):1237–44.
Smith DB, Simmonds P, Members Of The International Committee On The Taxonomy Of Viruses Hepeviridae Study G, Jameel S, Emerson SU, Harrison TJ, et al. Consensus proposals for classification of the family Hepeviridae. J Gen Virol. 2015;96(Pt 5):1191–2.
Meng XJ. Expanding Host range and cross-species infection of Hepatitis E virus. PLoS Pathog. 2016;12(8):e1005695.
Thakur V, Ratho RK, Kumar S, Saxena SK, Bora I, Thakur P. Viral Hepatitis E and chronicity: a growing public health concern. Front Microbiol. 2020;11:577339.
Kamar N, Lhomme S, Abravanel F, Marion O, Peron JM, Alric L, et al. Treatment of HEV Infection in patients with a solid-organ transplant and chronic hepatitis. Viruses. 2016;8(8):222.
Kenfak-Foguena A, Schöni-Affolter F, Bürgisser P, Witteck A, Darling KE, Kovari H, et al. Hepatitis E virus seroprevalence and chronic infections in patients with HIV. Switzerland Emerg Infect Dis. 2011;17(6):1074–8.
Tavitian S, Péron JM, Huynh A, Mansuy JM, Ysebaert L, Huguet F, et al. Hepatitis E virus excretion can be prolonged in patients with hematological malignancies. J Clin Virol. 2010;49(2):141–4.
Kamar N, Rostaing L, Legrand-Abravanel F, Izopet J. How should hepatitis E virus infection be defined in organ-transplant recipients? Am J Transplant. 2013;13(7):1935–6.
Geng Y, Zhao C, Huang W, Harrison TJ, Zhang H, Geng K, et al. Detection and assessment of infectivity of hepatitis E virus in urine. J Hepatol. 2016;64(1):37–43.
Lee GH, Tan BH, Teo EC, Lim SG, Dan YY, Wee A, et al. Chronic infection with camelid hepatitis E virus in a liver transplant recipient who regularly consumes camel meat and milk. Gastroenterology. 2016;150(2):355-7.e3.
Narayanan S, Abutaleb A, Sherman KE, Kottilil S. Clinical features and determinants of chronicity in hepatitis E virus infection. J Viral Hepat. 2019;26(4):414–21.
Cao D, Cao QM, Subramaniam S, Yugo DM, Heffron CL, Rogers AJ, et al. Pig model mimicking chronic hepatitis E virus infection in immunocompromised patients to assess immune correlates during chronicity. Proc Natl Acad Sci U S A. 2017;114(27):6914–23.
Moal V, Textoris J, Ben Amara A, Mehraj V, Berland Y, Colson P, et al. Chronic hepatitis E virus infection is specifically associated with an interferon-related transcriptional program. J Infect Dis. 2013;207(1):125–32.
Shi X, Ren S, Zhang B, Guo S, He W, Yuan C, et al. Analysis of the role of Purα in the pathogenesis of Alzheimer’s disease based on RNA-seq and ChIP-seq. Sci Rep. 2021;11(1):12178.
Dasgupta U, Mishra GP, Dikshit HK, Mishra DC, Bosamia T, Roy A, et al. Comparative RNA-Seq analysis unfolds a complex regulatory network imparting yellow mosaic disease resistance in mungbean [Vigna radiata (L.) R. Wilczek]. PLoS One. 2021;16(1):e0244593.
Wang H, Deng Y, Peng S, Yan L, Xu H, Wang Q, et al. RNA-Seq based transcriptome analysis in oral lichen planus. Hereditas. 2021;158(1):39.
Yu D, Lim J, Wang X, Liang F, Xiao G. Enhanced construction of gene regulatory networks using hub gene information. BMC Bioinformatics. 2017;18(1):186.
Seo CH, Kim JR, Kim MS, Cho KH. Hub genes with positive feedbacks function as master switches in developmental gene regulatory networks. Bioinformatics. 2009;25(15):1898–904.
Wang L, Liu L, Wang L. An overview: rabbit hepatitis E virus (HEV) and rabbit providing an animal model for HEV study. Rev Med Virol. 2018;28(1):e1961.
Han J, Lei Y, Liu L, Liu P, Xia J, Zhang Y, et al. SPF rabbits infected with rabbit hepatitis E virus isolate experimentally showing the chronicity of hepatitis. PLoS One. 2014;9(6):e99861.
Kamar N, Selves J, Mansuy JM, Ouezzani L, Péron JM, Guitard J, et al. Hepatitis E virus and chronic hepatitis in organ-transplant recipients. N Engl J Med. 2008;358(8):811–7.
Zhang F, Qi Y, Harrison TJ, Luo B, Zhou Y, Li X, et al. Hepatitis E genotype 4 virus from feces of monkeys infected experimentally can be cultured in PLC/PRF/5 cells and upregulate host interferon-inducible genes. J Med Virol. 2014;86(10):1736–44.
Sadler AJ, Williams BR. Interferon-inducible antiviral effectors. Nat Rev Immunol. 2008;8(7):559–68.
Patil RN, Karpe YA. Uncovering the roles of miR-214 in hepatitis E virus replication. J Mol Biol. 2020;432(19):5322–42.
Honda Y, Kondo J, Maeda T, Yoshiyama Y, Yamada E, Shimizu YK, et al. Isolation and purification of a non-A, non-B hepatitis-associated microtubular aggregates protein. J Gen Virol. 1990;71(Pt 9):1999–2004.
Kaczkowski B, Rossing M, Andersen DK, Dreher A, Morevati M, Visser MA, et al. Integrative analyses reveal novel strategies in HPV11,-16 and -45 early infection. Sci Rep. 2012;2:515.
Power D, Santoso N, Dieringer M, Yu J, Huang H, Simpson S, et al. IFI44 suppresses HIV-1 LTR promoter activity and facilitates its latency. Virology. 2015;481:142–50.
Kanehisa M, Furumichi M, Tanabe M, Sato Y, Morishima K. KEGG: new perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res. 2017;45(D1):D353–61.
Liu P, Bu QN, Wang L, Han J, Du RJ, Lei YX, et al. Transmission of hepatitis E virus from rabbits to cynomolgus macaques. Emerg Infect Dis. 2013;19(4):559–65.
This work was supported by a grant from National Institutes for Food and Drug Control (2022C4) to M. Li and a grant from the Natural Science Foundation of Shanxi Province (20210302124092) to Y. Wang.
Ethics approval and consent to participate
We confirmed that all experimental protocols were approved by a Committee of Laboratory Animal Welfare and Ethics of National Institutes for Food and Drug Control and First Clinical Medical College, Shanxi Medical University. All methods were carried out in accordance with relevant guidelines and regulations. All methods are reported in accordance with ARRIVE guidelines (https://arriveguidelines.org) for the reporting of animal experiments. Rabbits in this study were sacrificed by overdose of sodium pentobarbital (65 mg/kg).
Consent for publication
The authors declare no conflict of interest exists.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Figure 1.The comparison of hub genes expressions in positive dropouts and negative rabbits. The expressions of MX1, ISG15, IFIT2, RSAD2, OAS2, IFIT3, IFI44, IFI44L and DDX60 were measured by real-time PCR. (*, Pvalue < 0.05).
About this article
Cite this article
Li, M., Wang, Y., Li, K. et al. Characterization of highly expressed novel hub genes in hepatitis E virus chronicity in rabbits: a bioinformatics and experimental analysis. BMC Vet Res 18, 239 (2022). https://doi.org/10.1186/s12917-022-03337-x
- Hepatitis E virus
- Chronic infection