Skip to main content

Characterization of highly expressed novel hub genes in hepatitis E virus chronicity in rabbits: a bioinformatics and experimental analysis



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.

Peer Review reports


Hepatitis E virus (HEV), which is the leading cause of acute viral hepatitis worldwide, can lead to 20 million HEV infections annually [1]. 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 [2]. 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 [5]. However, HEV can establish chronic infections in immunocompromised individuals, including patients with solid organ transplant (SOT) [6], human immunodeficiency virus (HIV) infection [7] and hematological malignancies [8]. HEV chronic infections are mostly caused by HEV3 [9], and few cases have been observed with HEV4 [10] and 7 infection [11]. 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 [13], and interferon-stimulated genes may be involved in the persistence of an HEV infection [14]. 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.

Fig. 1
figure 1

Volcano plotting of candidate DEGs between the whole blood of chronic HEV infection patients and control patients by bioinformatics analysis (P value < 0.05 and |log 2 FC|> 1)

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.

Fig. 2
figure 2

The enrichment analysis in HEV chronic infected samples. A GO terms of targets of DEGs. B Significantly infuenced KEGG pathways by enrichment analysis. C Enrichment analysis by Funrich

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).

Fig. 3
figure 3

The PPI network construction and hub genes identification. A PPI network was constructed on the basis of the STRING database. B 10 hub genes of HEV chronic infected samples were identified

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.

Fig. 4
figure 4

The characteristics of HEV chronic infection in rabbits. A The fecal viral load levels of groups P and N. B Serum HEV antigen dynamic changes of groups P and N. C Anti-HEV antibody levels of groups P and N

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).

Fig. 5
figure 5

The comparison of hub genes expressions in rabbits. The expressions of MX1, ISG15, IFIT2, RSAD2, OAS2, IFIT3, IFI44, IFI44L and DDX60 were measured by real-time PCR. (*, P value < 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 ([12]. Studies have shown that chronic hepatitis E (CHE) patients had lower levels of IL-1 receptor agonist and TNF-α concentration [12], and their total counts of CD2, CD3 and CD4 T cells were also significantly lower [22].

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 [25]. 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 ( 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 ( and software R (version 3.4.4, 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 ( and Funrich software ( 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 [29]. 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 ( 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).

Sample collection

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) [30].

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).

Table 1 Primer information for hub genes detection

Statistical analysis

Statistical analysis was performed using the SPSS PASW Statistics v18.0 statistical software package (SPSS, Inc., USA, 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


Cell components


Molecular functions


Biological processes


Gene ontology (GO)


Kyoto Encyclopedia of Genes and Genomes


Phosphate-buffered saline


Chronic hepatitis E


Hepatitis C virus


Human papillomavirus


Alanine transaminase


  1. Nimgaonkar I, Ding Q, Schwartz RE, Ploss A. Hepatitis E virus: advances and challenges. Nat Rev Gastroenterol Hepatol. 2018;15(2):96–110.

    Article  Google Scholar 

  2. Hoofnagle JH, Nelson KE, Purcell RH. Hepatitis E. N Engl J Med. 2012;367(13):1237–44.

    Article  CAS  Google Scholar 

  3. 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.

    Article  CAS  Google Scholar 

  4. Meng XJ. Expanding Host range and cross-species infection of Hepatitis E virus. PLoS Pathog. 2016;12(8):e1005695.

    Article  Google Scholar 

  5. 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.

    Article  Google Scholar 

  6. 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.

    Article  Google Scholar 

  7. 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.

    Article  Google Scholar 

  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.

    Article  Google Scholar 

  9. 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.

    Article  CAS  Google Scholar 

  10. 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.

    Article  Google Scholar 

  11. 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.

    Article  Google Scholar 

  12. 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.

  13. 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.

    Article  CAS  Google Scholar 

  14. 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.

    Article  CAS  Google Scholar 

  15. 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.

    Article  CAS  Google Scholar 

  16. 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.

    Article  CAS  Google Scholar 

  17. 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.

    Article  CAS  Google Scholar 

  18. 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.

    Article  Google Scholar 

  19. 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.

    Article  CAS  Google Scholar 

  20. 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.

  21. 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.

    Article  Google Scholar 

  22. 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.

    Article  CAS  Google Scholar 

  23. 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.

    Article  CAS  Google Scholar 

  24. Sadler AJ, Williams BR. Interferon-inducible antiviral effectors. Nat Rev Immunol. 2008;8(7):559–68.

    Article  CAS  Google Scholar 

  25. Patil RN, Karpe YA. Uncovering the roles of miR-214 in hepatitis E virus replication. J Mol Biol. 2020;432(19):5322–42.

    Article  CAS  Google Scholar 

  26. 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.

    Article  CAS  Google Scholar 

  27. 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.

    Article  Google Scholar 

  28. 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.

    Article  CAS  Google Scholar 

  29. 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.

    Article  CAS  Google Scholar 

  30. 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.

    Article  Google Scholar 

Download references


Not applicable.


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.

Author information

Authors and Affiliations



Experiment and data processing: M. Li, Y. Wang, H. Lan, K. Li. Wrote the manuscript: M. Li and Y. Wang. Revised the manuscript: M. Li and C. Zhou. All the authors read and approved the final manuscript.

Corresponding authors

Correspondence to Manyu Li or Cheng Zhou.

Ethics declarations

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 ( for the reporting of animal experiments. Rabbits in this study were sacrificed by overdose of sodium pentobarbital (65 mg/kg).

Consent for publication

Not applicable.

Competing interests

The authors declare no conflict of interest exists.

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:

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).

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 The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

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).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: