Uncovering neglected subtypes and zoonotic transmission of Hepatitis E virus (HEV) in Brazil
Virology Journal volume 20, Article number: 83 (2023)
Hepatitis E virus (HEV) circulation in humans and swine has been extensively studied in South America over the last two decades. Nevertheless, only 2.1% of reported HEV strains are available as complete genome sequences. Therefore, many clinical, epidemiological, and evolutionary aspects of circulating HEV in the continent still need to be clarified. Here, we conducted a retrospective evolutionary analysis of one human case and six swine HEV strains previously reported in northeastern, southern, and southeastern Brazil. We obtained two complete and four nearly complete genomic sequences. Evolutionary analysis comparing the whole genomic and capsid gene sequences revealed high genetic variability. This included the circulation of at least one unrecognized unique South American subtype. Our results corroborate that sequencing the whole capsid gene could be used as an alternative for HEV subtype assignment in the absence of complete genomic sequences. Moreover, our results substantiate the evidence for zoonotic transmission by comparing a larger genomic fragment recovered from the sample of the autochthonous human hepatitis E case. Further studies should continuously investigate HEV genetic diversity and zoonotic transmission of HEV in South America.
Hepatitis E is a zoonotic liver disease infecting 20 million people and causing 44,000 deaths yearly . Hepatitis E virus (HEV) has been detected in a wide range of species and is the sole hepatotropic virus with documented zoonotic transmission . HEV members belong to the family Hepeviridae which is subdivided into Orthohepevirinae and Parahepevirinae subfamilies. The Orthohepevirinae subfamily is divided into four genera (Avihepevirus, Chirohepevirus, Paslahepevirus and Rocahepevirus) of which infect mammals and birds, while the Parahepevirinae subfamily has the Piscihepevirus genus of which infects trout and salmon. The members of the Paslahepevirus balayani species are divided into eight HEV genotypes. Genotypes 1–4 are associated with self-limited acute hepatitis with exclusively human-human (HEV-1 and HEV-2), or zoonotic transmission (HEV-3/HEV-4) [3, 4]. Genotypes 5–8 infect wild mammals [5,6,7] with one single report of HEV-7 in human, so far .
HEV subtyping criteria have been debated mainly because of their subjectivity and lack of consistency [9,10,11]. Although subtyping is not within the remit of the International Committee on Taxonomy of Viruses (ICTV), the ICTV Hepeviridae Study Group has tried to develop an agreed upon set of subtype reference sequences to improve clarity for researchers, epidemiologists and clinicians . Currently, viruses in genotype 3 of the species Paslahepevirus balayani are classified into 13 recognized subtypes. HEV unassigned sequences can be recognized as new proposed subtypes if at least three complete genomic sequences from epidemiologically unrelated source form a distinct phylogenetic group . The distribution of genotypes and subtypes differs between geographical regions and hosts [13, 14] and may be associated with different clinical outcomes [15, 16]. Although distance-based criteria for assigning new sequences to subtypes can be used, their limits are not clearly established. Therefore, sequences are assigned to subtypes according to phylogenetic position compared to reference sequences [10, 13].
In South America, due to the low number of partial and complete sequences (~ 340 entries compared to ~ 11,000 from Europe) deposited in GenBank, there is a lack of knowledge about HEV genetic variability (Supplementary Table 1). All previously reported sequences obtained clustered into genotype 3 and were detected in swine herds, pork products, human cases, experimental studies in primates or environmental samples. The few partial sequences obtained for Brazilian sequences indicate that they cluster into subtypes 3b, 3c, 3d, 3f, 3h and 3i. However, the absence of complete genomes prevents accurate analysis of HEV genetic diversity in Brazil. Therefore, in this study, we sequenced available HEV positive samples from previous studies to obtain complete genomic sequences.
Materials and methods
Sera and feces samples were obtained from six positive swine and one human sample previously reported in studies from the northeastern, southeastern, and southern regions in Brazil [17,18,19,20,21]. The study with human sampling was approved by the Fiocruz institutional ethical committee (protocol PO 307/06) . Sample data, including access numbers, hosts, geographical area, genotype and reported subtype are detailed in Table 1. RNA was extracted and purified from stool or serum samples using TRIzolTM Reagent (Invitrogen) following the manufacturer’s recommendations. cDNA was synthesized with the ImProm-II™ Reverse Transcription System (Promega). RT-PCR was performed using the SuperScript™ III One-Step RT-PCR System with Platinum™ Taq DNA Polymerase kit (Invitrogen) following the manufacturer’s instructions. Primers were designed based on template complete genomes of sequences genetically close to the 156 Brazilian partial HEV-3 sequences available in the GenBank.
The complete genome of PRsw1 and PEsw1 were obtained by Sanger sequencing, and the other positive samples were also initially sequenced by Sanger method (primers available in Supplementary Table 1). The failed runs were further tentatively sequenced using DNA-based High-throughput sequencing (HTS) using MiSeq (Illumina) . Sequencing libraries were prepared using Nextera XT Library Kit, and sequencing was performed using a paired-end strategy with a V3–150 cycles flow cell. Low-quality sequencing reads were removed using the default parameters of Trimmomatic 0.36 . Assembly was performed using Velvet 1.0.10 , default parameters. Sequence data were compared to all HEV genotype 3 complete or nearly complete genomes (genome size > 6.5 kb; n = 747) available in GenBank as of 28.03.2023, including the references for each HEV-3 subtype [12, 13]. Aiming to investigate the amount of unresolved to fully resolved trees, the phylogenetic signal was verified through likelihood mapping analysis of 10,000 random quartets . Then, phylogenetic reconstruction based on Maximum Likelihood was performed using IQ-Tree server (http://iqtree.cibiv.univie.ac.at/), through default settings, using the best-fit substitution model TIM2 + F + R4 with 1000 ultrafast bootstrap replicates.
Results and discussion
We were able to sequence two complete genomic sequences (PRsw1 and PEsw1) (Table 1) and, four nearly complete porcine-derived genomic sequence (RJsw1, RJsw2, RJsw3 ad RJsw4), and a fragment of approximately 1000 nucleotides of the first recorded human-derived sequence in Brazil (RJh1). (Fig. 1a). The two complete genomic sequences exhibited the typical features of an HEV genome and were 7215 and 7233 nucleotides length, respectively. Both HEV-3 genomes contain three partially overlapping open reading frames, with ORF3 overlapping ORF1 and ORF2 (Fig. 1b).
Because we were not able to obtain complete genomic sequences of all positive samples and considering the epidemiological importance of the genetic information, we performed further analyses with multiple datasets to compare with HEV subtype reference sequences. Accordingly, we compared the two complete genomic sequences obtained, the nearly complete genomic sequences and the complete capsid sequence of all six sequences from pigs to a set of HEV-3 reference sequences (Fig. 2a and b). Our analysis suggests that full genome sequences from South (PRsw1) and Southeastern (RJsw1), previously classified as HEV-3 subtype i, likely represent a new subtype with nucleotide divergence of 15.5% from subtype ‘m’ together with two other sequences found in Uruguay (MW596896 and MZ969073) (Fig. 3; Supplementary Table 2). Moreover, the nucleotide divergence of PEsw1 from the northeast region, previously classified as subtype 3f in our previous study , differed from the HEV-3f reference sequence at 12.0%, but remained in the same subtype classification. Interestingly, both nearly-complete-genome- and capsid-based analysis suggests that strain RJsw2 (from Southern region) clustered between the subtype 3i and a group of unassigned sequences represented by MF959764 (UN-V). Strains RJsw3 and RJsw4, from the same region, are related to subtype 3b (Fig. 2b). On the one hand, the robust results on the capsid-based phylogeny corroborate that ORF2 sequences can be alternatively used in the absence of full genomic sequences (Fig. 2c), as previously suggested [11, 20, 26], including with subgenomic regions of ORF1 and ORF2 optimized for robust phylogenetic inference and subtyping . By contrast, this illustrates the high genetic variability of the HEV strains circulating in Brazil and upholds the need to increase the number of complete genomic sequences from different samples from Latin America. Overall, the consistent topology derived from multiple inferences supports the robustness of the new proposed subtype.
To investigate the zoonotic origin of the first autochthonous human case of acute hepatitis E reported in Brazil, we compared the partial fragments of the capsid and the ORF1. Phylogenetic analyses using the ORF2 fragment, as well as the previously published ORF1 fragment, indicated that the sequence from the human patient (RJh1) clusters among those from pigs from southeastern region (RJsw3 and RJsw4) with nucleotide identities of 99.5 and 99.6%, respectively (Fig. 2d; Supplementary Fig. 1 and Supplementary Table 2). In the previous study , only small fragments of ORF1 and ORF2 were available, and by providing a larger fragment of ORF2, our data support the robustness of the evolutionary inferences and the evidence of zoonotic transmission. Therefore, regardless of whether strains RJsw3 and RJsw4 represent a new subtype or not, our results show strong evidence that the first reported case of HEV in Brazil is likely of zoonotic origin.
Noteworthy, our study is limited by the lack of success in obtaining larger fragments of four sequences which were probably related to the volume availability and conservation status of the samples. However, except for the human sample, the study provided a good amount of retrospective information that allowed robust analysis. Although our findings and those previously published support that the capsid-based analyses are robust, as recommended by the ICTV Hepeviridae Study Group, it would be necessary to obtain full genomic sequence before proposing a new subtype .
Finally, the genetic variability of HEV in Brazil, along with the occurrence of zoonotic transmission, illustrates the challenges and lack of knowledge of HEV epidemiology in South America. Therefore, further studies, including more representative sequences, are needed to investigate the genetic diversity of HEV in humans, animal reservoirs, food, and environmental samples to obtain an accurate picture of the genetic variability of HEV in South America and its implication in zoonotic transmission and pathogenicity.
The raw data supporting the conclusions in this paper will be made available by the authors without reservation.
Hepatitis E. [https://www.who.int/news-room/fact-sheets/detail/hepatitis-e.
Nan Y, Wu C, Zhao Q, Zhou EM. Zoonotic Hepatitis E virus: an ignored risk for Public Health. Front Microbiol. 2017;8:2396.
Kenney SP. The Current Host Range of Hepatitis E Viruses. Viruses 2019, 11.
Purdy MA, Drexler JF, Meng XJ, Norder H, Okamoto H, Van der Poel WHM, Reuter G, de Souza WM, Ulrich RG, Smith DB. ICTV Virus Taxonomy Profile: Hepeviridae 2022. J Gen Virol 2022, 103.
Takahashi M, Nishizawa T, Nagashima S, Jirintai S, Kawakami M, Sonoda Y, Suzuki T, Yamamoto S, Shigemoto K, Ashida K, et al. Molecular characterization of a novel hepatitis E virus (HEV) strain obtained from a wild boar in Japan that is highly divergent from the previously recognized HEV strains. Virus Res. 2014;180:59–69.
Woo PC, Lau SK, Teng JL, Tsang AK, Joseph M, Wong EY, Tang Y, Sivakumar S, Xie J, Bai R, et al. New hepatitis E virus genotype in camels, the Middle East. Emerg Infect Dis. 2014;20:1044–8.
Woo PC, Lau SK, Teng JL, Cao KY, Wernery U, Schountz T, Chiu TH, Tsang AK, Wong PC, Wong EY, Yuen KY. New Hepatitis E virus genotype in bactrian camels, Xinjiang, China, 2013. Emerg Infect Dis. 2016;22:2219–21.
Lee GH, Tan BH, Teo EC, Lim SG, Dan YY, Wee A, Aw PP, Zhu Y, Hibberd ML, Tan CK, et al. Chronic infection with Camelid Hepatitis E Virus in a liver transplant recipient who regularly consumes Camel meat and milk. Gastroenterology. 2016;150:355–357e353.
Wang S, Wei W, Luo X, Cai X. Genome-wide comparisons of phylogenetic similarities between partial genomic regions and the full-length genome in Hepatitis E virus genotyping. PLoS ONE. 2014;9:e115785.
Nicot F, Dimeglio C, Migueres M, Jeanne N, Latour J, Abravanel F, Ranger N, Harter A, Dubois M, Lameiras S, et al. Classification of the zoonotic Hepatitis E virus genotype 3 into distinct subgenotypes. Front Microbiol. 2020;11:634430.
Oliveira-Filho EF, Konig M, Thiel HJ. Genetic variability of HEV isolates: inconsistencies of current classification. Vet Microbiol. 2013;165:148–54.
Smith DB, Simmonds P, Izopet J, Oliveira-Filho EF, Ulrich RG, Johne R, Koenig M, Jameel S, Harrison TJ, Meng XJ, et al. Proposed reference sequences for hepatitis E virus subtypes. J Gen Virol. 2016;97:537–42.
Smith DB, Izopet J, Nicot F, Simmonds P, Jameel S, Meng XJ, Norder H, Okamoto H, van der Poel WHM, Reuter G, Purdy MA. Update: proposed reference sequences for subtypes of hepatitis E virus (species Orthohepevirus A). J Gen Virol. 2020;101:692–8.
De Sabato L, Suffredini E, Pasquale SD, La Rosa G, De Santis P, Giammarioli M, Vaccari G, Bartolo ID. Novel subtypes and unexpected heterogeneity of hepatitis E viral strains in wild boar captured in a small area in Central Italy. Transbound Emerg Dis. 2022;69:e2541–50.
Peeters M, Schenk J, De Somer T, Roskams T, Locus T, Klamer S, Subissi L, Suin V, Delwaide J, Starkel P, et al. Viral clade is associated with severity of symptomatic genotype 3 hepatitis E virus infections in Belgium, 2010–2018. J Hepatol. 2023;78:67–77.
Abravanel F, Dimeglio C, Castanier M, Peron JM, Kamar N, Lhomme S, Izopet J. Does HEV-3 subtype play a role in the severity of acute hepatitis E? Liver Int. 2020;40:333–7.
dos Santos DR, Vitral CL, de Paula VS, Marchevsky RS, Lopes JF, Gaspar AM, Saddi TM, Júnior NC, Guimarães FeR, Júnior JG, et al. Serological and molecular evidence of hepatitis E virus in swine in Brazil. Vet J. 2009;182:474–80.
Lopes Dos Santos DR, Lewis-Ximenez LL, da Silva MF, de Sousa PS, Gaspar AM, Pinto MA. First report of a human autochthonous hepatitis E virus infection in Brazil. J Clin Virol. 2010;47:276–9.
dos Santos DR, de Paula VS, de Oliveira JM, Marchevsky RS, Pinto MA. Hepatitis E virus in swine and effluent samples from slaughterhouses in Brazil. Vet Microbiol. 2011;149:236–41.
Oliveira-Filho EF, Dos Santos DR, Duraes-Carvalho R, da Silva A, de Lima GB, Batista Filho AFB, Pena LJ, Gil LH. Evolutionary study of potentially zoonotic hepatitis E virus genotype 3 from swine in Northeast Brazil. Mem Inst Oswaldo Cruz. 2019;114:e180585.
Gardinali NR, Barry AF, da Silva PF, de Souza C, Alfieri AF, Alfieri AA. Molecular detection and characterization of hepatitis E virus in naturally infected pigs from brazilian herds. Res Vet Sci. 2012;93:1515–9.
Sanger F. Determination of nucleotide sequences in DNA. Science. 1981;214:1205–10.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20.
Zerbino DR. Using the Velvet de novo assembler for short-read sequencing technologies. Curr Protoc Bioinformatics 2010, Chap. 11:Unit 11 15.
Schmidt HA, Strimmer K, Vingron M, von Haeseler A. TREE-PUZZLE: maximum likelihood phylogenetic analysis using quartets and parallel computing. Bioinformatics. 2002;18:502–4.
Purdy MA, Sue A. The effect of phylogenetic signal reduction on genotyping of hepatitis E viruses of the species Orthohepevirus A. Arch Virol. 2017;162:645–56.
EF Oliveira-Filho was supported by FACEPE (grant 1380 − 5.05/15) and CNPq (grant DCR- 0003-5.05/16) until 2018. RD-C is supported by the Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP), Brazil, grants 2019/01255-9 and 2021/03684-4 (Young Investigator Program). The authors also thank the National Laboratory for Scientific Computing (LNCC/MCTI, Brazil) for providing HPC resources of the Santos Dumont supercomputer (ID #45691, project “virusevolution”).
Open Access funding enabled and organized by Projekt DEAL.
Ethics approval and consent to participate
The research study involving human sampling was reviewed and approved in accordance with Fiocruz’s ethics committee (CEP-Fiocruz) under authorization no. PO 307/06. Animal sampling was conducted as approved by the Institutional Committee animal ethics committee (CEUA-Fiocruz) under authorization no. PO 0132/01.
Consent for publication
All the authors provided critical feedback and consent to the publication of the manuscript.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
dos Santos, D.R.L., Durães-Carvalho, R., Gardinali, N.R. et al. Uncovering neglected subtypes and zoonotic transmission of Hepatitis E virus (HEV) in Brazil. Virol J 20, 83 (2023). https://doi.org/10.1186/s12985-023-02047-6