Next Article in Journal
Characterization of Novel Transcripts in Pseudorabies Virus
Previous Article in Journal
The Use of Convalescent Sera in Immune-Electron Microscopy to Detect Non-Suspected/New Viral Agents
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 7
Issue 5
10.3390/v7052704
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hepatitis E Virus Genotype 3 Diversity: Phylogenetic Analysis and Presence of Subtype 3b in Wild Boar in Europe

by
Ariel Vina-Rodriguez
1,
Josephine Schlosser
1,
Dietmar Becher
2,
Volker Kaden
3,
Martin H. Groschup
1 and
Martin Eiden
1,*
1
Institute for Novel and Emerging Infectious Diseases, Friedrich-Loeffler-Institut, Südufer 10, 17493 Greifswald-Insel Riems, Germany
2
Micromun GmbH, Greifswald, Walther-Rathenau-Straße 49A, 17489 Greifswald, Germany
3
Institute of Infectology, Friedrich-Loeffler-Institut, Südufer 10, 17493 Greifswald-Insel Riems, Germany
*
Author to whom correspondence should be addressed.
Viruses 2015, 7(5), 2704-2726; https://doi.org/10.3390/v7052704
Submission received: 14 April 2015 / Revised: 13 May 2015 / Accepted: 20 May 2015 / Published: 22 May 2015
(This article belongs to the Section Animal Viruses)
Browse Figures
Versions Notes

Abstract

:
An increasing number of indigenous cases of hepatitis E caused by genotype 3 viruses (HEV-3) have been diagnosed all around the word, particularly in industrialized countries. Hepatitis E is a zoonotic disease and accumulating evidence indicates that domestic pigs and wild boars are the main reservoirs of HEV-3. A detailed analysis of HEV-3 subtypes could help to determine the interplay of human activity, the role of animals as reservoirs and cross species transmission. Although complete genome sequences are most appropriate for HEV subtype determination, in most cases only partial genomic sequences are available. We therefore carried out a subtype classification analysis, which uses regions from all three open reading frames of the genome. Using this approach, more than 1000 published HEV-3 isolates were subtyped. Newly recovered HEV partial sequences from hunted German wild boars were also included in this study. These sequences were assigned to genotype 3 and clustered within subtype 3a, 3i and, unexpectedly, one of them within the subtype 3b, a first non-human report of this subtype in Europe.

1. Introduction

The Hepatitis E virus (HEV) is a causative agent of acute hepatitis in developing countries in Asia, Africa and Latin America where it is transmitted primarily via contaminated drinking water. Sporadic cases of HEV are reported in developed countries, partially imported by travelers from endemic areas, but there are also an increasing number of reports of autochthonous HIV infections. The transmission route of most of the autochthonous infections in industrialized countries still remains unclear. Reports of transfusion and transplant related infections exist [1,2,3], but accumulating evidence suggests that hepatitis E is a zoonotic disease with domestic pigs and wild boars being the main reservoirs. Moreover, the consumption of undercooked meat products poses a risk for HEV infection [4,5,6,7].
Studies on HEV RNA detection in animals revealed, that HEV is ubiquitous in domestic pigs and wild boars throughout Europe [8]. This includes the United Kingdom [9], France [10], Germany [11,12,13,14,15], Hungary [16], Italy [17,18], The Netherlands [19,20], Belgium [21], Spain [22], Slovenia [23], Czech Republic [24] and Sweden [25]. HEV infection in farmed pigs affects up to 80%–100% of the animals worldwide and usually occurs at the age of 2–4 months [26].
The virion is approximately 27–34 nm in diameter and most likely icosahedral. HEV has a positive sense single-stranded RNA genome of approximately 7.2 kb, which contains a short 5' untranslated region (UTR), a short 3' UTR and three open reading frames (ORF1, ORF2 and ORF3) [27]. The ORF1 encodes for viral non-structural proteins carrying domains with methyl transferase, helicase and replicase activities [28]. The ORF2 codes for the viral capsid protein of about 660 amino acids. The ORF3 is almost completely overlapped by the ORF2 and codes for a small phosphoprotein of about 114 amino acids, which is putatively responsible for the virion egress from infected cells [29].
A new proposed consensus for the HEV classification [30] divides the Hepeviridae family in two genera: Orthohepevirus and Piscihepevirus. The latter includes only isolates from cutthroat trout so far. The genus Orthohepevirus is further subdivided into four species: Orthohepevirus A with isolates from human, pig, wild boar, deer, mongoose, rabbit and camel, and Orthohepevirus B, C and D with avian and other mammal isolates. Orthohepevirus A is subsequently divided into at least six genotypes (HEV-1, HEV-2, etc.). HEV-1 and 2 include exclusively human HEV strains, whereas HEV-3 and 4 can also infect other animal species, particularly domestic pigs and wild boar. The separation in anthropotropic (HEV-1 and -2) and enzoonotic (HEV-3 and -4) forms may have occurred more than 500 year ago [31].
HEV-1 is found in Asia and Africa [32], whereas HEV-2 was first isolated in Mexico [33] and later in Africa [34]. HEV-4 includes strains from sporadic human HEV cases in Asia [35]. HEV-3 was isolated initially from human cases in the USA [36] and has been detected in all continents including Europe [37].
Since 2001, 2703 human HEV cases [38] have been reported in Germany, which include an increasing number of non-travel associated autochthonous cases. The origin of infection remains unclear for most of the autochthonous cases, however, often the suspected HEV sources are domestic pigs and wild boars [39]. Similar increases are also reported in other European countries [40].
The role of different HEV-3 genetic variants in the evolution of the disease [1,7,39,41,42], the possibility of tracking the routes of infection and the influences of human activity on it [43,44,45,46] are currently under study. The direct comparison of isolates is still hampered by the limited number of complete genome (CG) sequences available. Due to this limitation, the subtyping scheme proposed by Lu et al. [32] has been commonly used and have been supported by epidemiological and/or statistical analysis [4,43,45,47,48] but questions have arisen, partially due to the lack of commonly accepted reference sequences for some subtypes [1]. In response to this, and to the increasing number of partial sequences, the subtyping of genotype 3 strains was actualized in order to provide an update of the subtyping scheme of HEV-3 and of the set of reference sequences. We carried out the classification of newly recovered HEV isolates from German wild boar and detected HEV subtype 3b strains for the first time in animals in Europe and, possibly, the first from a wild animal outside Japan.

2. Materials and Methods

2.1. Samples and RNA Extraction

Blood samples were collected from wild boar hunted in Mecklenburg-Western Pomerania during the seasons 1996/1997 (955 samples) and 2005/2006 (58 samples). Liver samples were collected in 2009/2010 from 134 animals hunted in the region of Greifswald and from another five from Western Pomerania. All samples were stored at −80 °C prior to their use. RNA was extracted with the RNeasy Mini Kit (QIAGEN, Hilden, Germany) according to the manufacturer’s protocol. A synthetic RNA (IC) was used as internal extraction control [49].

2.2. Primers and Probe Design

For primer and probe design, an alignment of 351 HEV sequences was constructed using the Vector NTI Advanced v.10 (Invitrogen, Carlsbad, CA, USA), BioEdit v.7.0.5.3 [50] and MEGA v6 [51] software. This alignment was manually curated using both the nucleotide and the deduced amino acid sequences. Very similar sequences were not included (more than 99% identity). HEV-1, -2, -3 and -4 genotypes were included (with preference to genotype 3), covering all subtypes, and including 131 CG (48 of them cited by Lu et al. [32]) and 65 German HEV sequences (the accession numbers are included in Supplementary Table S01). For genotyping and subtyping, four sets of nested degenerated primers were selected from this alignment, which target different regions of the genome. Previously published primers [14] were used to amplify an RNA-dependent RNA polymerase (RdRp) region. A novel diagnostic quantitative real-time RT-PCR assay (qRT-PCR) that targets ORF3 was also designed, which we already used in a recently published work [52]. Primer and probes used are listed in Table 1 (nucleotide positions refer to FJ705359, strain wbGER27, a German wild boar isolate [14]) and were included in the GenBank sequences entries.

2.3. PCR

The diagnostic/screening RT-qPCR was performed using the QuantiTec Probe RT-PCR kit (QIAGEN) in 25 µL reaction volume. In all reactions, the final concentration of each primer was 0.8 µM, and of the probe 0.1 µM if present. A volume of 5 µL of the RNA eluate was added. The reverse transcription (RT) was carried out at 50 °C for 30 min, followed by denaturation/activation at 95 °C for 15 min. DNA was amplified immediately with 45 cycles at 95 °C (10 s), 55 °C (25 s) and 72 °C (25 s).
Table
Table 1. Primer and probes used in this study. All nucleotide positions refer to FJ705359 (strain wbGER27, a German wild boar isolate). Abbreviations: open reading frame (ORF), hyper variable region (HVR), RNA dependent RNA polymerase (RdRp).
Table 1. Primer and probes used in this study. All nucleotide positions refer to FJ705359 (strain wbGER27, a German wild boar isolate). Abbreviations: open reading frame (ORF), hyper variable region (HVR), RNA dependent RNA polymerase (RdRp).
Region Name and Internal Length (nt)Primer NamePositionSequenceStepProduct Length (bp)
ORF1 -5´ This studyHEV.ORF1_F133–58CCCAYCAGTTYATWAAGGCTCCTGGCRT-PCR493
HEV.ORF1_R1497–525TGCARDGARTANARRGCNAYNCCNGTCTC
HEV.ORF1_F298–126AAYTCYGCCYTGGCGAATGCTGTGGTGGTnested PCR302
HEV.ORF1_R2377–399CCVCGRGTNGGRGCRGWRTACCA
HVR (for genotype 3) This studyHEV.HVR_F12069–2091TTYTCYCCTGGGCAYMTYTGGGART-PCR401
HEV.HVR_R12441–2469TTAACCARCCARTCACARTCYGAYTCAAA
HEV.HVR_F2a2135–2157ACYTGGTCHACATCTGGYTTYTCnested PCR263 or 293
HEV.HVR_F2b2165–2184TTYTCCCCYCCTGAGGCGGC
HEV.HVR_R22405–2427TACACCTTRGCSCCRTCRGGRTA
RdRp 280 (4312–4591) (Johne et al. 2010)HEV-cs4181–4203TCGCGCATCACMTTYTTCCARAART-PCR469
HEV-cas4628–4650GCCATGTTCCAGACDGTRTTCCA
HEV-csn4287–4311GTGCTCTGTTTGGCCCNTGGTTYMGnested PCR330
HEV-casn4592–4617CCAGGCTCACCRGARTGYTTCTTCCA
ORF3 225 (5205–5429) This studyHEV.ORF3_F15126–5145MGGKTRGAATGAATAACATGRT-PCR326 or 362
HEV.ORF3_R25430–5451GGCGCTGGGAYTGGTCRCGCCA
HEV.ORF3_R15467–5487CAGYTGGGGYAGRTCGACGRC
HEV.ORF3_F25182–5204GGGCTGTTCTGTTKYTGYTCYTCnested PCR219 or 269
HEV.ORF3_R25430–5451GGCGCTGGGAYTGGTCRCGCCA
HEV.ORF3_R2a5382–5401CGAGGGCGAGCTCCAGCCCC
modified Diagnostic-qPCR - This study and (Schlosser et al. 2014)HEV.Fa5278–5294GTGCCGGCGGTGGTTTCRT-qPCR81
HEV.Fb5278–5296GTGCCGGCGGTGGTTTCTG
HEV.R5340–5359GCGAAGGGGTTGGTTGGATG
HEV.P5300–5320FAM-TGACMGGGTTGATTCTCAGCC-BHQ1
ORF2 187 (6277–6488) This studyHEV.ORF2_F16205–6223CDGCNACYCGBTTYATGAART-PCR393
HEV.ORF2_R1a6573–6598GTKAGRGARAGCCAWAGYACATCATT
HEV.ORF2_R1b6573–6598GTRAGNGADAGCCACARRACATCATT
HEV.ORF2_F26276–6301GCBYTHACNYTRTTYAAYCTTGCTGAnested PCR241
HEV.ORF2_R26489–6517TGYTCRTGYTGRTTRTCRTARTCYTGDAT
The determination of the HEV RNA concentration was carried out using a standard curve according to a synthetic external calibrator. This calibrator encompassed the 81 bp sequence of the diagnostic qRT-PCR amplicon and included the T7 promoter sequence at the 5'-end for in vitro transcription. The RNA synthesis and DNA degradation were carried out by the Riboprobe® Combination System—T3/T7 RNA Polymerase (Promega Corporation’s, Madison, WI, USA); and the QIAamp Viral RNA Mini Kit (QIAGEN) kit was used for RNA isolation (without carrier RNA). The RNA concentration was estimated with the Quant-It™ RNA Assay Kit, Broad Range (Invitrogen) and confirmed by endpoint dilution PCR.
For genotyping, the initial RT-qPCR was performed with the QuantiTec SYBR Green RT-PCR kit (QIAGEN) in 25 µL reaction volume using 5 µL of the sample RNA. The thermal profile applied was: 30 min at 50 °C for RT, 15 min 95 °C denaturation/activation followed by 45 cycles of 95 °C for 10 s, 55 °C for 25 s, 72 °C for 25 s and 80 °C for 5 s (with fluorescence reading). A final dissociation curve generation step was also included. Two microliters of the resulting solution was added to 23 µL of the Maxima™ SYBR Green/ROX qPCR Master Mix kit (Fermentas, Canada) containing the primers and the PCR was carried out under similar conditions (45 cycles of: 95 °C for 10 s, 55 °C for 25 s, 72 °C for 25 s and 80 °C for 5 s with a final dissociation curve generation step).

2.4. Sequencing, Phylogenetic Analysis and Classification

RT-PCR or nested-PCR products were directly sequenced with the corresponding forward and reverse PCR primers using the BigDye Terminator v1.1 Cycle Sequencing Kit on the DNA sequencer “3130 Genetic Analyzer” (Applied Biosystems, Waltham, MA, USA).
The newly generated sequences were manually inserted in the multi-alignment previously used for primer design. This multi-alignment was updated with new HEV sequences (NCBI, 2014-12-15), up to more than 1400 sequences (Supplementary Table S01), mainly genotype 3 (near 1300 sequences). This included all HEV-3 sequences longer 1000 nt, all HEV NCBI nucleotide entries with the keyword “Germany”, the 135 sequences cited by Lu et al. [32] and other sequences from around the word, but particularly from other European countries. The evolutionary history was inferred by using the Maximum Likelihood method based on the Kimura 2-parameter model [53]. The trees with the highest log are shown. The percentage of trees (boostrap values for 500 replicates for the first CG tree—and 100 replicates for all others trees) in which the associated taxa clustered together is shown next to the branches. Initial tree(s) for the heuristic search were obtained by applying the Neighbor-Joining (NJ) method to a matrix of pairwise distances estimated using the Maximum Composite Likelihood (MCL) approach. A discrete Gamma distribution was used to model evolutionary rate differences among sites (5 categories). The trees were drawn to scale, with branch lengths measured in the number of substitutions per site. Codon positions included were 1st + 2nd + 3rd. Evolutionary analyses were conducted in MEGA6 [51]. The alignment in FASTA format and the auxiliary worksheet for classification, selection and automatic labeling of sequences in MEGA are provided in Supplementary Files S01, S02 (updated versions are planned to be available from the authors). The sheet enables a quick selection of all the sequences spanning a given genomic region, which can be set in alignment coordinates (nt position) or by referencing to the sequences M73218 (Burma) or FJ705359 (wbGER27). Starting with the CG, and followed by the longest sequences, we built phylogenetic trees, and labeled each sequence as a “reference” to be used in subsequence classifications only if the tree reproduced the same topology as the tree for CG and if the clade was supported by bootstrap values of more than 70%. That is: HEV-3 sequences were labeled as “reference” only if the subtype was unambiguously determined. We assumed that all sequences from one strain represent the same genome sequence, and if one of them was labeled “reference”, all the others were also labeled.

3. Results

3.1. HEV RNA Detection

HEV RNA was detected in 32 out of 955 blood samples from 1996/97 and three out of 58 blood samples from 2005/2006, which suggests a virus prevalence of about 3.4% and 5.2%, respectively. In addition, HEV RNA was found in 14 out of 134 wild boar derived liver samples from the Greifswald region, giving a prevalence rate of about 10.4%. Finally, two wild boar liver samples (WS03-09 and WS05-09) from individual hunts were also positive. All HEV RNA positive samples were re-tested with the PCR for genotyping, and partial sequences from 12 animals could be recovered and subjected to phylogenetic analysis.

3.2. Phylogenetic Analyses

A reference phylogenetic tree was constructed based on 166 CG sequences, including eight German HEV isolates and 98 strains of genotype 3 (Figure 1a, which corresponds to Figure 4 in [32]). The hypervariable region (HVR) (2146–2358 nt) was excluded from this analysis. All nucleotide positions refer to sequence M73218. This tree confirmed a good separation of the HEV-3 from all other HEV genotypes. The sequences clustered into four monophyletic groups: “3jab”, “3chi”, “3feg” and “rabbit”. A detailed overview of the HEV-3 clade is shown in Figure 1b.
Viruses 07 02704 g001a 1024Viruses 07 02704 g001b 1024
Figure 1. Molecular Phylogenetic analysis of 166 complete HEV genomes by Maximum Likelihood method based on the Kimura 2-parameter model. The percentage of trees (from 500 bootstrap replicates) in which the associated taxa clustered together is shown next to the branches when over 70%. Initial tree(s) for the heuristic search were obtained by applying the Neighbor-Joining method to a matrix of pairwise distances estimated using the Maximum Composite Likelihood (MCL) approach. A discrete Gamma distribution was used to model evolutionary rate differences among sites (five categories (+G, parameter = 0.5255)). The tree is drawn to scale, with branch lengths measured in the number of substitutions per site. All positions with less than 95% site coverage were eliminated. There were a total of 6868 positions in the final dataset. Evolutionary analyses were conducted in MEGA6. (a) Global view of the unrooted tree, which corresponds to Figure 4 in Lu et al. [32]; (b) Detailed view of the HEV-3 clade. (*)—references sequences cited by Lu; HRC-HE104—strain used in the HEV RNA WHO standard; wb—wild boar. See Supplementary Table 1 for more information.
Figure 1. Molecular Phylogenetic analysis of 166 complete HEV genomes by Maximum Likelihood method based on the Kimura 2-parameter model. The percentage of trees (from 500 bootstrap replicates) in which the associated taxa clustered together is shown next to the branches when over 70%. Initial tree(s) for the heuristic search were obtained by applying the Neighbor-Joining method to a matrix of pairwise distances estimated using the Maximum Composite Likelihood (MCL) approach. A discrete Gamma distribution was used to model evolutionary rate differences among sites (five categories (+G, parameter = 0.5255)). The tree is drawn to scale, with branch lengths measured in the number of substitutions per site. All positions with less than 95% site coverage were eliminated. There were a total of 6868 positions in the final dataset. Evolutionary analyses were conducted in MEGA6. (a) Global view of the unrooted tree, which corresponds to Figure 4 in Lu et al. [32]; (b) Detailed view of the HEV-3 clade. (*)—references sequences cited by Lu; HRC-HE104—strain used in the HEV RNA WHO standard; wb—wild boar. See Supplementary Table 1 for more information.
Viruses 07 02704 g001aViruses 07 02704 g001b
The analysis was continued only with HEV-3, excluding the rabbit sequences, which form a well-separated clade.
Comparisons amongst complete HEV-3 genomes display different levels of diversity within this genotype (Figure 2). The analysis of the frequency of corrected distances between sequences shows a possible separation around 0.14 substitutions per site. Sequences with lower differences belong to the same subtype. The graphic of the frequencies depicts two additional intermediary peaks, which are the basis for the definition of groups (3jab, 3chi and 3feg) and of major clades 3-I (3jab and 3chi) and 3-II (3feg) [54] (Figure 2).
The newly obtained tree (Figure 3a) for complete HEV-3 genomes further segregated the sequences in subtypes clades 3j, 3a, 3b, 3c, 3h, 3e and 3f. This classification was supported by bootstrap values of 99%–100%. Subtypes 3i and 3g were represented by individual isolates. Ninety-six of these sequences were marked as “reference”. Two strains were marked 3ef. German strains were grouped in the 3a, 3c, 3i, 3e and 3f subtypes. This Figure represents our best approximation to the true topology and displays the reference structure for all subsequent trees.
Viruses 07 02704 g002 1024
Figure 2. Frequency diagram of corrected distances between sequences from Figure 3a. It shows a possible separation of subtypes and suggests a possible definition of groups (3chi, 3jab and 3feg) and of major clades 3-I and 3-II (Norder, 2009). Pairwise distances were estimated using the Maximum Composite Likelihood (MCL) approach and grouped in intervals of 0.003. A discrete Gamma distribution was used to model evolutionary rate differences among sites (five categories). All positions with less than 95% site coverage were eliminated.
Figure 2. Frequency diagram of corrected distances between sequences from Figure 3a. It shows a possible separation of subtypes and suggests a possible definition of groups (3chi, 3jab and 3feg) and of major clades 3-I and 3-II (Norder, 2009). Pairwise distances were estimated using the Maximum Composite Likelihood (MCL) approach and grouped in intervals of 0.003. A discrete Gamma distribution was used to model evolutionary rate differences among sites (five categories). All positions with less than 95% site coverage were eliminated.
Viruses 07 02704 g002
Viruses 07 02704 g003 1024
Figure 3. Molecular phylogenetic analysis of: (a) CG-[HVR]: 97 complete HEV-3 genomes, excluding the HVR (M73218:2146-2358 nt); (b) ORF1.1860nt: 96 partial HEV-3 genome sequences spanning the region 93–1952 nt, which correspond to a partial sequence JQ807482(●) obtained in this study from the animal WS05-09; (c) ORF1.318nt: 124 partial HEV-3 genome sequences spanning the region Burma.M73218: 61–378 nt approximately correspond to Figure 2 in [32]. Three of the sequences obtained in this study are included here (●): JQ807489, JQ807482 and JQ807479 (Bugewitz strain). Detailed view of the sequences used in b and c can be seen in supplementary Figures S01 and S02.
Figure 3. Molecular phylogenetic analysis of: (a) CG-[HVR]: 97 complete HEV-3 genomes, excluding the HVR (M73218:2146-2358 nt); (b) ORF1.1860nt: 96 partial HEV-3 genome sequences spanning the region 93–1952 nt, which correspond to a partial sequence JQ807482(●) obtained in this study from the animal WS05-09; (c) ORF1.318nt: 124 partial HEV-3 genome sequences spanning the region Burma.M73218: 61–378 nt approximately correspond to Figure 2 in [32]. Three of the sequences obtained in this study are included here (●): JQ807489, JQ807482 and JQ807479 (Bugewitz strain). Detailed view of the sequences used in b and c can be seen in supplementary Figures S01 and S02.
Viruses 07 02704 g003
Viruses 07 02704 g004 1024
Figure 4. Molecular phylogenetic analysis of: (a) HVR.319nt: 126 partial HEV-3 genome sequences spanning the region 2094–2412 nt. The branches for sequences AF455784-3g and JN564006-3a are very large and have been truncated, including eight of the sequences obtained in this study (●); (b) RdRp.280nt: 147 partial HEV-3 genome sequences spanning the region 4284–4563 nt, including three of the sequences obtained in this study (●); (c) ORF2.187nt: 124 partial HEV-3 genome sequences spanning the region 6277–6463 nt, including five of the sequences obtained in this study (●), and stool pools from The Netherlands: 3c (▲) NLSW36 and NLSW105 and 3c+3f: (▼)NLSW20 and NLSW99 (see discussion). Detailed view of the sequences can be seen in supplementary Figures S07–S09.
Figure 4. Molecular phylogenetic analysis of: (a) HVR.319nt: 126 partial HEV-3 genome sequences spanning the region 2094–2412 nt. The branches for sequences AF455784-3g and JN564006-3a are very large and have been truncated, including eight of the sequences obtained in this study (●); (b) RdRp.280nt: 147 partial HEV-3 genome sequences spanning the region 4284–4563 nt, including three of the sequences obtained in this study (●); (c) ORF2.187nt: 124 partial HEV-3 genome sequences spanning the region 6277–6463 nt, including five of the sequences obtained in this study (●), and stool pools from The Netherlands: 3c (▲) NLSW36 and NLSW105 and 3c+3f: (▼)NLSW20 and NLSW99 (see discussion). Detailed view of the sequences can be seen in supplementary Figures S07–S09.
Viruses 07 02704 g004
Viruses 07 02704 g005 1024
Figure 5. Scheme of the alignment of the deduced amino acid sequences including the HVR, which is flanked by conserved regions (some sequences only partially span this region and gaps in the conserved flanks are not part of the alignment). Obtained from the Alignment Explorer of MEGA6 (screenshot).
Figure 5. Scheme of the alignment of the deduced amino acid sequences including the HVR, which is flanked by conserved regions (some sequences only partially span this region and gaps in the conserved flanks are not part of the alignment). Obtained from the Alignment Explorer of MEGA6 (screenshot).
Viruses 07 02704 g005
Each sequence in the alignment was analyzed by constructing a tree based on its exact (or near equal) length and including all the other sequences that span the region. More than 50 trees were analyzed, and as a result 68 more sequences were selected as “reference”, 48 of them longer than 1000 nt. Trees constructed for sequences more than 1500 nt length retained the same topology and support, with bootstrap values of at least 80% and each of the subtypes with more than one sequence. A representative tree based on 96 partial sequences spanning a region of 1860 nucleotide from position 93 to 1925 in the ORF1 (ORF1.1860nt) is depicted in Figure 3b.
Reduced sequence length maintains the same basic tree topology, but leads to a reduced or no bootstrap value support. This goes along with an increasing number of available sequences, better reflecting the HEV-3 diversity (Supplementary Figure S03). The tree based on a 318-nucleotide region from the 5’ section of the genome (position 61–378, ORF1.318nt), which includes the alignment of 124 partial sequences, is shown in Figure 3b. Newly recovered sequences from German wild boar isolates could be classified as subtype 3b (isolate WS03-09) and as subtype 3i (isolate WS05-09). Isolate WS34-10 was tentatively assigned to subtype 3b.
An overview of all the newly classified isolates from Mecklenburg-Western Pomerania and the corresponding accession numbers are depicted in Supplementary Table S02. A phylogenetic tree based on a 242-nucleotide region (positions 125–366 nt, ORF2.242nt) including 294 isolates is deposited at Supplementary Figure S03a.
The HVR (nucleotide position 2146–2358) was found to be particularly variable and was therefore manually aligned considering the deduced amino acid sequences. A scheme of a protein alignment of this region is shown in Figure 5. The longest stable HVR was found in strains of the subtype 3f, which were isolated from humans in France and Spanish pigs. In comparison, all other isolates contained shorter sequences without changing the coding frame. A tree based on a 319-nucleotide region (position 2094–2412) with 126 partial HEV-3 nucleotide sequences is shown in Figure 4a. This grouping assigned the strain WS03/09 to subtype 3b, the two isolates WS34-10 and WS35-10 were associated only weakly with subtype 3a and five isolates (WS05-09, WS21-10, 5160, 5304, 4322) clustered within subtype 3i.
A tree for 147 partial sequences within the RdRp region (280 nucleotide length, position 4284–4563, RdRp.280nt, Figure 4b, which approximately corresponds to Figure 6 presented by Lu et al. [32]) reproduced a similar topology compared to the CG. This tree confirmed the assignment of the strain WS03/09 to subtype 3b and a weak association of two isolates WS34-10 and WS28-10 to subtype 3a.
Finally, we carried out a phylogenetic analysis with sequences of the ORF3 (position 5180–5404, ORF3.225nt, Supplementary Figure S10/S11) and ORF2 (position 6277–6463, ORF2.187nt in Figure 4c, which partially correspond to Figure 5 in [32]) including 114 and 124 sequences, respectively. In both cases, monophyletic groups were not well supported leading to bootstrap values below 70% and in the case of the ORF3 without clear separation into subtypes. Thus, the subtypes were proposed (if possible) assuming the subtype of the nearest sequences. Using ORF3 derived sequences, we assigned three isolates to subtype 3a (isolate 8603, isolate 4701 and 4973), isolate WS03-09 to 3b, three isolates (WS28-10, WS34-10, WS35-10) weakly to 3b and two isolates to 3i (WS05-09, WS21-10) In the case of partial ORF2 derived sequences, two isolates segregated to subtype 3a (8603, 4701) and three isolates into subtype 3b (WS34-10, WS35-10, WS03-09).
In summary, twelve new isolates from Germany were analyzed using partial sequences from different regions of the genome (Supplementary Table S02). More than 1200 publicly available sequences, representing more than 1100 isolates, were subtyped accordingly and are listed in Supplementary Table S01. The classified strains are summarized by geographic region, subtype and host in Table 2.
Table
Table 2. Subtype distribution of the 1109 isolates of genotype 3, represented by the 1283 sequences analyzed. Summarized by continent, country, and source—shown in order: human‘pig‘wild animal‘other. For example 1’2’3’4 means: 1 human isolate; 2 isolates from pig; 3 isolates from wild animals; and 4 from others sources (environment samples, water, etc.). As some isolates have undetermined subtypes and are not shown, not all rows or columns sum the real total. Due to various types of biases in the selection of the biological samples by the original authors and in our selection of sequences, the only approximate analysis possible was the comparison of the proportion of subtypes by host in different countries or regions. An alternative, graphical view is available at Supplementary Figure S13, which can also be explored interactively.
Table 2. Subtype distribution of the 1109 isolates of genotype 3, represented by the 1283 sequences analyzed. Summarized by continent, country, and source—shown in order: human‘pig‘wild animal‘other. For example 1’2’3’4 means: 1 human isolate; 2 isolates from pig; 3 isolates from wild animals; and 4 from others sources (environment samples, water, etc.). As some isolates have undetermined subtypes and are not shown, not all rows or columns sum the real total. Due to various types of biases in the selection of the biological samples by the original authors and in our selection of sequences, the only approximate analysis possible was the comparison of the proportion of subtypes by host in different countries or regions. An alternative, graphical view is available at Supplementary Figure S13, which can also be explored interactively.
Origen33jab3chi 3feg Total
3a3b3d3j3h3L3i3c3e3ef3f3g3k
Africa 7 7
CMR 0’4 0’4
CAN 0’3 0’3
America 272 2 15 47
ARG 1 3’1 5’1
URY 10 10
BRA 0’1 0’1
CUB 7’6 7’6
CAN 2’11’1 3’2
MEX 0’1 0’1
USA 6’4’0’1 6’4’0’1
Asia 1867333 13151 114
CHN 1’10’3 1’4
JPN 10’3’230’12’23 1’2 4’3’4 46’20’29
KOR 0’2 0’2
KGZ 0’1 0’1
MNG 0’1 0’1 0’2
NZL 1’1 1 2’1
THA 4’1 4’1
Europe8592 10 2527732622089 940
ESP 0’0’0’1 0’0’0’1 0’2 7’27’0’13 7’29’0’15
FRA 1 3’1 3’14’1’0’1 24’9 35’12’0’1
GRC 1 1 2’3
ITA 1 0’1 0’2 1’3
AUT 0’4 1’1 1’5
NLD 1’10 16’29’12’2’0’1 5’33’0’1 24’74’1’2
DEU49’2’110’0’1 0’0’4 0’0’2140’6’1411’8’6215’1’11 82’17’61
HUN 0’8’3 3’12’2 2 5’20’5
CZE 2 1’41’4’2 4’8’2
SRB 0’4 0’1’2 0’5
SVN 0’30’3
GBR36 165268 57 507
SWE1 20’0’1 4’1 7’0’1’1
Total8104713520152527733932131031109

4. Discussion

The high sequence variability of HEV genomes is the central problem affecting the screening and diagnostic methods for the detection and quantification of the viral RNA. Therefore, the conserved ORF3 region offers a promising target for PCR assays (as reported by Jothikumar et al. [55]). Within this region, considering potential secondary structures, the high cg content, and trying to avoid the frequently observed non-specific signals with high Ct-values, we developed a new real-time RT-PCR assay (Table 1).
The HEV RNA detection rate of 3.4% in nearly 1000 sera collected from swine in 1996/97 and of 5.2% collected in 2005/2006 in the here presented study was similar to a previously reported 5.3% prevalence rate found in 189 samples collected in 1995/96 from the same region [12]. This indicates a constant circulation of HEV in this region. In addition, the prevalence of HEV in liver samples was 10.4%, which corresponds to the 14.9% [14] and 18.1% [13] found in wild boar derived liver samples from other regions in Germany.
The use of different and short genome regions for genotyping can lead to incongruences and provides insufficient evidence for establishing or refuting phylogenetic hypotheses [56]. Considering the given restriction, Lu et al. [32] proposed a comprehensive subtype scheme for the phylogenetic analysis of Hepatitis E virus, which has been commonly used. Nevertheless, Lu pointed to some incongruence in this scheme due to the use of different regions and to the small number of sequences within some subtypes that were available.
A major source of inconsistency during subtyping is the combined use of short sequences and the pooling of samples with a subsequent in silico concatenation of sequences. For example, Lu pointed out that the Arkell strain isolated from a pool of pig feces in Canada [57] is probably an artificial mixture of sequences, which could explain the inconsistent classification of this strain using different regions. Similarly, ORF1.242nt and ORF2.301nt sequences derived from pools of 20 to 60 pig faces [19] were used by Lu to define the 3c subtype (in the major clade I), but if only the ORF.148nt is analyzed, two of them distantly cluster 3f in the major clade II. Both 3c and 3f subtypes are common in Dutch pigs. Most of these inconsistences can be avoided using the original set (or a new consensus set) of reference sequences. Nevertheless, in general, this effect could appear without pooling of samples, due to co-infections and to true recombination between distant strains that are presumably rare events.
Another source of inconsistency is the lack, or insufficient number, of strains in some subtypes. The introduction of a new subtype based only on one single or a few short sequences can be error prone, due to laboratory artifacts, insufficient phylogenetic information, recombination, etc. For example, a new Hepatitis C virus (HCV) subtype (among other requisites) is created only when one complete genome (CG) and two other sufficiently informative sequences are available [58]. Basically, the CG will serve as a reference along the whole genome and the other two will determine the cluster, or prove the existence of a relatively recent common ancestor. This is an obvious problem within the HEV group 3chi: only 12 sequences longer than 1500 nt are available, from which only seven comprise CG (Supplementary Figure S12). Within the group 3chi, the best-described subtype appears to be subtype 3c (in the set of sequences we had already analyzed), with three CG and a large number of partial sequences. However, the 3i subtype is represented by only one CG, thus making it difficult to compare the sequences from different genomic regions that could be assigned to this subtype. Three CG were assigned to the subtype 3h, but they are highly divergent. It is important to note that the current poor structure of the group 3chi is not due to a rare detection of 3chi sequences, but rather to a relatively limited effort to obtain CG or nearly complete sequences. Most of the long HEV-3 sequences have been obtained in Japan (54 CG of 97 and 58 of 117 sequences longer 1500 nt), where 3chi apparently does not circulate.
Based on Lu’s classification scheme, we generated an updated phylogenetic tree with all newly available CG of genotype 3 and used the corresponding structure (Figure 3a) as template for subsequent classification of other strains based on partial sequences only. Our experimentally recovered partial sequences from different genomic regions, covering the 5’ ORF1, HVR, ORF3 and ORF2 (target regions selected in this study) and the RdRp regions [59] were originally selected to match that of the majority of the European sequences. The results show that partial sequences from our ORF1, HVR, RdRp and ORF2 regions generated trees with similar structures compared to the reference tree and can be used to subtype most sequences. In contrast, sequences from ORF3.225 are only partially suitable for classification up to the subtype level. In this context, the tree topology of Figure 3b and Figure 5 correspond to previously published trees (Figure 2 and Figure 6, Lu et al. [32])
The HVR is not a typical hypervariable region, but rather a genotypically diverse sequence [31]. The variability of this region has two components: (1) a higher mutation rate, and (2) insertions and/or deletions of one or two triplets or of much longer sequences (but maintaining the same reading frame) [60,61]. Taking this into account, it is almost trivial to find the right alignment manually and to decide whether the region should be included or not in the phylogenetic analyses. This alignment alone (Figure 5) allows an approximate reconstruction of the evolutionary history of the HEV genotypes and subtypes.
Based on the analysis of the 1652 nt part of the ORF2 not overlapping the ORF3 region, Purdy et al. [31] (including near 55 HEV-3 sequences) calculated the Time of the Most Recent Common Ancestor (TMRCA) of the four genotypes to be around the year 1475, and for HEV-3 and -4 around 1595. The TMRCA of major clades 3-I and 3-II (HEV-3 excluding the rabbit sequences) was determined near the year 1790, and for the clade 3-I (corresponding to 3chi and 3jab together) in 1865. Each group (3jab, 3chi and 3feg) has evolved roughly from 1900. While subtypes 3a and 3e have a TMRCA in 1945, the 3b and 3f+3ef are approximately 15 year older (without the 3ef, it is reduced for 3f to around 1960). Interestingly, the whole HEV-1 is only 100 years old. In a more recent analysis [43] using only the ORF2.301nt (including 208 HEV-3 sequences) the TMRCA for both major clades I and II together was dated to 1810 and for the clade I alone (3chi + 3jab groups) to 1895. Each group was correspondingly estimated: 3jab—1920; 3chi—1919; and 3feg—1889. For subtypes: 3a—1959; 3b—1944; 3f—1935; and 3e—1917 (include sequences not included by Purdy). The subtype 3d was the last separation, in 2002. Another study [45] shows compatible results and possible sources of minor discordance are discussed [43].
Not all possible methods of evolutionary tree reconstruction were thoroughly evaluated, but we noted that modeling evolutionary rate differences among sites have a major impact on the consistency of the results and that the tree generated with a ML method have longer internal branches and shorter terminal branches than with the NJ method, which is considered a good characteristic [62].
Partial sequences from twelve field isolates could be recovered and they all clustered within genotype 3. Three strains (4701, 4973, 8603) from the retrospective samples segregated to subtype 3a, which has been already detected in German autochthonous human infections [63], in wild boars around the city of Potsdam, in Brandenburg [14], as well as in human and pig samples from Bavaria [15]. Subtype 3a appears to be worldwide represented in samples from humans, pigs and wild animals (especially boar). In American samples, 3a could be the predominant and potential indigenous subtype (Supplementary Table S02), and import of USA pigs has been pointed to as a source of infection in South Korea and Japan [45].
Unexpectedly, the strain WS03-09 collected from an animal hunted in Western Pomerania clustered within subtype 3b for the four regions analyzed. In Europe, subtype 3b has not been detected in wild boars or domestic pig populations so far. This subtype probably originated from Japan [45] and has been mainly identified in humans, wild boars, domestic pigs and deer from that country. It has also been isolated from one Canadian pig and (reportedly) from humans and swine in Brazil [64]. In Europe, only one human isolate from France was grouped into subtype 3b [65]. We report here the first non-Japanese 3b isolate obtained from a wild animal.
WS 34/10, WS 35/10, and WS 28/10 could be classified 3jab, but could not be unambiguously subtyped and further investigations using longer sequences are needed to define if they cluster into existing subtypes (no -a or -b) or whether they define a new subtype within the 3jab group.
No sequences within subtype 3c were discovered in this study, although 3c appears to be specific for Central Europe and is the major subtype in Germany, Netherlands and recently United Kingdom, detected in humans, pigs and wild boars.
Sequences from animals WS 05/09, WS 21/10, 5160, 5304 and 4322 clustered within subtype 3i, which is closely related to the 3c, and could have similar distribution. Curiously, none of the 101 analyzed strains from The Netherlands was classified 3i, but 21 (all from wild board) out of 162 Germans strains were classified 3i. Until now, this subtype has been detected in Germany in only wild boars, but in Austria and Argentina has been also detected in humans.
Other subtypes were not detected, although especially subtype 3e and 3f are widely distributed in Europe (Table 2). The subtype 3e appears to be more widely distributed than 3c, including clusters of sequences from Japan and West Europe but it is more represented in Central Europe. In contrast, 3f sequences are more frequently found in Spain and France, and also found in other European countries. 3f has been also detected outside Europe in Thailand (two Japan patients were infected with this subtype after a trip to this country). Interesting, only one 3f strain was isolated from wild boar (from 94 total wild boar analyzed), but 204 were isolated from humans and domestic pigs (out of 998 from all subtypes). Finally, we recommend the use of partial sequences only when the obtained tree reproduces the same structure compared to the CG tree. Ideally, sequences with more than 1000 nt should be used for classification. In contrast, sequences below 200 nt should be avoided for subtyping. In particular, the commonly used ORF2.148nt, and the OFR2.171nt generate poorly structured trees. ORF3 sequences are sufficient for genotype, but not for subtype determination. HVR sequences should only be used for intra genotype comparisons, and alignments have to be checked manually, especially in the case of sequences with long insertions, which are impossible to be compared with the reference sequences. Do not define or modify subtypes based only on a single CG or only on short sequences (less than 1500 nt).

5. Conclusions

We designed RT-PCR assays for screening, quantification and genotyping of HEV-3 strains, and detected viral RNA in wild boar samples from Mecklenburg-Western Pomerania, Germany. Twelve strains clustered into subtypes 3a, 3i and, unexpectedly, also 3b, which is a common subtype in Japan, but has not been reported in animals in Europe. The phylogenetic trees based on our partial sequences of ORF1, RdRp, HVR and ORF2 regions reproduced similar topology as obtained from complete genome analysis and were useful for subtyping.
More than 30 different PCR fragments and the corresponding genomic regions have been used for genotyping and subtyping so far, which is a source of ambiguous subtyping schemes and inadequate classification. The presented study offers an updated set of reference sequences for the relatively simple and neutral subtype scheme proposed by Lu et al. [32], which could eliminates most of the existing incongruences and creates the basis for new hypotheses regarding the Hepatitis E epidemiology. A comprehensive subtyping of HEV-3 according to this classification scheme could enable a detailed view of the spread of HEV-3 strains among pigs, wild life and humans, and could allow determining the consequences of infections with different subtypes on humans and finally help limit the potential spread of the disease.

Acknowledgments

We thank Birke Boettcher and Gina Lucht for excellent technical assistance.

Author Contributions

Conceived and designed the experiments: AVR, MHG, ME. VK and DB performed sample collection and characterization. JS was involved in RNA isolation and sequencing. AVR, JS, ME performed data analysis. AVR and ME wrote the paper. AVR, JS, MHG, ME edited the manuscript.

Supplementary Information

Tables

Supplementary Table S01. Global classification of Hepatitis E virus isolates using complete genomes (CG) or partial sequences from open reading frame 1 (ORF1), hyper variable region (HVR), RNA dependent RNA polymerase (RdRp), open reading frame 3 (ORF3) and open reading frame 2 (ORF2). All nucleotide positions refer to FJ705359 (wbGer27). Abbreviations: (*) reference sequences as cited by Lu et al. [32]; (**) strain used in the HEV RNA WHO standard (WHO, 2013); 3b.JIO—virulent strain from Japan; the sequences selected as “reference” are marked with an +; wb—wild boar; 3-letter country code after ISO 3166-1 alpha-3. When the information was available, this country code reflects the country from where the infection was imported. German isolates derived from Federal Institute for Risk Assessment (FIRA), Robert-Koch-Institut (RKI), Friedrich-Loeffler-Institut (FLI), Germany, University of Regensburg, Institute of Medical Microbiology and Hygiene (UReg), Bernhard Nocht Institute for Tropical Medicine (BNI) and University of Veterinary Medicine Hannover (TiHo). Burma (Myanmar, MMR) was conserved due to the widespread use in HEV related literature. mac-hum—isolated from macaques original inoculated witch human feces, or from both.
Supplementary Table S02. Subtype designation and accession numbers of sequences from isolates recovered in this study. *tentative assignments.

Figures

Supplementary Figure S01. Detailed view of the sequences used in Figure 3b.
Supplementary Figure S02. Detailed view of the sequences used in Figure 3c.
Supplementary Figure S03. Molecular phylogenetic analysis of: (a) ORF1.242nt (left): 294 partial HEV-3 genome sequences spanning the region 125–366 nt, including all the sequences from Figure 2 in [32], a large number of European sequences, partial sequences from animals WS34-10: (●) JQ807489, and (●) JQ807482, and stool pools from Netherlands: 3c (▲) NLSW36 and NLSW105 and 3c + 3f: (▼)NLSW20 and NLSW99 (see discussion); (b) ORF2.280nt (center): 665 partial HEV-3 genome sequences spanning the region 6017–6296 nt including near 500 sequences from England and Wales (Ijaz et al. 2014), and from the same stool pools from The Netherlands: (▲) and (▼); (c) ORF2.241nt (right): 156 partial HEV-3 genome sequences spanning the region 5596–5836 nt. A detailed view of all the sequences used is given in Supplementary Figures S04, S05 and S06, respectively.
Supplementary Figure S04. Detailed view of the sequences used in Supplementary Figure S03a.
Supplementary Figure S05. Detailed view of the sequences used in Supplementary Figure S03b.
Supplementary Figure S06. Detailed view of the sequences used in Supplementary Figure S03c.
Supplementary Figure S07. Detailed view of the sequences used in Figure 4a.
Supplementary Figure S08. Detailed view of the sequences used in Figure 4b.
Supplementary Figure S09. Detailed view of the sequences used in Figure 4c.
Supplementary Figure S10. Molecular phylogenetic analysis of: (a) ORF3.225nt: 114 partial HEV-3 genome sequences spanning the region 5180–5404 nt, including sequences obtained in this study (●).
Supplementary Figure S11. Detailed view of the sequences used in in Supplementary Figure S10.
Supplementary Figure S12. Detailed view of the 3chi group from three partial genomic regions.
Supplementary Figure S13. Geographical distribution of subtypes by host (human, domestic pig, or wild animals). This information can be also interactively explored at: https://public.tableau.com/views/HEV-g3/Story3?:embed=y&:showTabs=y&:display_count=yes.

Supplementary Files

Supplementary File S01: HEV-A.fasta.zip. Multialignment in FASTA format, compressed.
Supplementary File S02: HEVsubtypingMEGAut.xlsx. Auxiliary worksheet for classification, selection and automatic labeling of sequences in MEGA.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Moal, V.; Gerolami, R.; Ferretti, A.; Purgus, R.; Devichi, P.; Burtey, S.; Colson, P. Hepatitis e virus of subtype 3i in chronically infected kidney transplant recipients in southeastern france. J. Clin. Microbiol. 2014, 52, 3967–3972. [Google Scholar] [CrossRef] [PubMed]
  2. Moal, V.; Ferretti, A.; Devichi, P.; Colson, P. Genome sequence of a hepatitis e virus of genotype 3e from a chronically infected kidney transplant recipient. Genome Announc. 2014, 2, e01156–e0115613. [Google Scholar] [CrossRef] [PubMed]
  3. Marano, G.; Vaglio, S.; Pupella, S.; Facco, G.; Bianchi, M.; Calizzani, G.; Candura, F.; Catalano, L.; Farina, B.; Lanzoni, M.; et al. Hepatitis E: An old infection with new implications. Blood Transfus. = Trasfus. Sangue 2015, 13, 6–20. [Google Scholar]
  4. Colson, P.; Borentain, P.; Queyriaux, B.; Kaba, M.; Moal, V.; Gallian, P.; Heyries, L.; Raoult, D.; Gerolami, R. Pig liver sausage as a source of hepatitis e virus transmission to humans. J. Infect. Dis. 2010, 202, 825–834. [Google Scholar] [CrossRef] [PubMed]
  5. Li, T.C.; Chijiwa, K.; Sera, N.; Ishibashi, T.; Etoh, Y.; Shinohara, Y.; Kurata, Y.; Ishida, M.; Sakamoto, S.; Takeda, N.; et al. Hepatitis e virus transmission from wild boar meat. Emerg. Infectious Dis. 2005, 11, 11–13. [Google Scholar] [CrossRef]
  6. Tei, S.; Kitajima, N.; Takahashi, K.; Mishiro, S. Zoonotic transmission of hepatitis e virus from deer to human beings. Lancet 2003, 362, 371–373. [Google Scholar] [CrossRef] [PubMed]
  7. Renou, C.; Pariente, A.; Cadranel, J.F.; Nicand, E.; Pavio, N. Clinically silent forms may partly explain the rarity of acute cases of autochthonous genotype 3c hepatitis e infection in france. J. Clin. Virol. 2011, 51, 139–141. [Google Scholar] [CrossRef] [PubMed]
  8. Lewis, H.C.; Wichmann, O.; Duizer, E. Transmission routes and risk factors for autochthonous hepatitis e virus infection in europe: A systematic review. Epidemiol. Infect. 2010, 138, 145–166. [Google Scholar] [CrossRef] [PubMed]
  9. Meader, E.; Thomas, D.; Salmon, R.; Sillis, M. Seroprevalence of hepatitis e virus in the UK farming population. Zoonoses Public Health 2010, 57, 504–509. [Google Scholar] [CrossRef] [PubMed]
  10. Kaba, M.; Davoust, B.; Marie, J.L.; Colson, P. Detection of hepatitis e virus in wild boar (sus scrofa) livers. Vet. J. 2010, 186, 259–261. [Google Scholar] [CrossRef] [PubMed]
  11. Preiss, J.C.; Plentz, A.; Engelmann, E.; Schneider, T.; Jilg, W.; Zeitz, M.; Duchmann, R. Autochthonous hepatitis e virus infection in germany with sequence similarities to other european isolates. Infection 2006, 34, 173–175. [Google Scholar] [CrossRef] [PubMed]
  12. Kaci, S.; Nöckler, K.; Johne, R. Detection of hepatitis e virus in archived german wild boar serum samples. Vet. Microbiol. 2008, 128, 380–385. [Google Scholar] [CrossRef] [PubMed]
  13. Adlhoch, C.; Wolf, A.; Meisel, H.; Kaiser, M.; Ellerbrok, H.; Pauli, G. High hev presence in four different wild boar populations in east and west germany. Vet. Microbiol. 2009, 139, 270–278. [Google Scholar] [CrossRef] [PubMed]
  14. Schielke, A.; Sachs, K.; Lierz, M.; Appel, B.; Jansen, A.; Johne, R. Detection of hepatitis e virus in wild boars of rural and urban regions in germany and whole genome characterization of an endemic strain. Virol. J. 2009, 6, e58. [Google Scholar] [CrossRef]
  15. Wenzel, J.J.; Preiß, J.; Schemmerer, M.; Huber, B.; Plentz, A.; Jilg, W. Detection of hepatitis e virus (hev) from porcine livers in southeastern germany and high sequence homology to human hev isolates. J. Clin. Virol. 2011, 52, 50–54. [Google Scholar] [CrossRef] [PubMed]
  16. Reuter, G.; Fodor, D.; Forgách, P.; Kátai, A.; Szucs, G. Characterization and zoonotic potential of endemic hepatitis e virus (hev) strains in humans and animals in hungary. J. Clin. Virol. 2009, 44, 277–281. [Google Scholar] [CrossRef] [PubMed]
  17. Ponterio, E.; di Bartolo, I.; Orru, G.; Liciardi, M.; Ostanello, F.; Ruggeri, F.M. Detection of serum antibodies to hepatitis e virus in domestic pigs in italy using a recombinant swine hev capsid protein. BMC Vet. Res. 2014, 10, e133. [Google Scholar] [CrossRef]
  18. Caruso, C.; Modesto, P.; Bertolini, S.; Peletto, S.; Acutis, P.L.; Dondo, A.; Robetto, S.; Mignone, W.; Orusa, R.; Ru, G.; et al. Serological and virological survey of hepatitis e virus in wild boar populations in northwestern italy: Detection of hev subtypes 3e and 3f. Arch. Virol. 2015, 160, 153–160. [Google Scholar] [CrossRef] [PubMed]
  19. Van der Poel, W.H.; Verschoor, F.; van der Heide, R.; Herrera, M.I.; Vivo, A.; Kooreman, M.; de Roda Husman, A.M. Hepatitis e virus sequences in swine related to sequences in humans, the netherlands. Emerg. Infect. Dis. 2001, 7, 970–976. [Google Scholar] [CrossRef] [PubMed]
  20. Rutjes, S.A.; Lodder, W.J.; Lodder-Verschoor, F.; van den Berg, H.H.; Vennema, H.; Duizer, E.; Koopmans, M.; de Roda Husman, A.M. Sources of hepatitis e virus genotype 3 in the netherlands. Emerg. Infect. Dis. 2009, 15, 381–387. [Google Scholar] [CrossRef] [PubMed]
  21. Hakze-van der Honing, R.W.; van Coillie, E.; Antonis, A.F.G.; van der Poel, W.H.M. First isolation of hepatitis e virus genotype 4 in europe through swine surveillance in the netherlands and belgium. PLoS ONE 2011, 6, e22673. [Google Scholar] [CrossRef] [PubMed]
  22. De Deus, N.; Seminati, C.; Pina, S.; Mateu, E.; Martín, M.; Segalés, J. Detection of hepatitis e virus in liver, mesenteric lymph node, serum, bile and faeces of naturally infected pigs affected by different pathological conditions. Vet. Microbiol. 2007, 119, 105–114. [Google Scholar] [CrossRef] [PubMed]
  23. Steyer, A.; Naglič, T.; Močilnik, T.; Poljšak-Prijatelj, M.; Poljak, M. Hepatitis e virus in domestic pigs and surface waters in slovenia: Prevalence and molecular characterization of a novel genotype 3 lineage. Infect. Genet. Evol. 2011, 11, 1732–1737. [Google Scholar] [CrossRef] [PubMed]
  24. Vasickova, P.; Psikal, I.; Widen, F.; Smitalova, R.; Bendova, J.; Pavlik, I.; Kralik, P. Detection and genetic characterisation of hepatitis e virus in czech pig production herds. Res. Vet. Sci. 2009, 87, 143–148. [Google Scholar] [CrossRef] [PubMed]
  25. Widén, F.; Sundqvist, L.; Matyi-Toth, A.; Metreveli, G.; Belák, S.; Hallgren, G.; Norder, H. Molecular epidemiology of hepatitis e virus in humans, pigs and wild boars in sweden. Epidemiol. Infect. 2011, 139, 361–371. [Google Scholar] [CrossRef] [PubMed]
  26. Meng, X.J. Recent advances in hepatitis e virus. J. Viral Hepat. 2010, 17, 153–161. [Google Scholar] [CrossRef] [PubMed]
  27. Tam, A.W.; Smith, M.M.; Guerra, M.E.; Huang, C.-C.C.; Bradley, D.W.; Fry, K.E.; Reyes, G.R. Hepatitis e virus (hev): Molecular cloning and sequencing of the full-length viral genome. Virology 1991, 185, 120–131. [Google Scholar] [CrossRef] [PubMed]
  28. Koonin, E.V.; Gorbalenya, A.E.; Purdy, M.A.; Rozanov, M.N.; Reyes, G.R.; Bradley, D.W. Computer-assisted assignment of functional domains in the nonstructural polyprotein of hepatitis e virus: Delineation of an additional group of positive-strand rna plant and animal viruses. Proc. Natl. Acad. Sci. USA 1992, 89, 8259–8263. [Google Scholar] [CrossRef] [PubMed]
  29. Okamoto, H. Hepatitis e virus cell culture models. Virus Res. 2011, 161, 65–77. [Google Scholar] [CrossRef] [PubMed]
  30. Smith, D.B.; Simmonds, P.; members of the International Committee on the Taxonomy of Viruses Hepeviridae Study Group; Jameel, S.; Emerson, S.U.; Harrison, T.J.; Meng, X.J.; Okamoto, H.; van der Poel, W.H.; Purdy, M.A. Consensus proposals for classification of the family hepeviridae. J. Gen. Virol. 2014, 95, 2223–2232. [Google Scholar] [CrossRef] [PubMed]
  31. Purdy, M.A.; Khudyakov, Y.E. Evolutionary history and population dynamics of hepatitis e virus. PLoS ONE 2010, 5, e14376. [Google Scholar] [CrossRef] [PubMed]
  32. Lu, L.; Li, C.; Hagedorn, C.H. Phylogenetic analysis of global hepatitis e virus sequences: Genetic diversity, subtypes and zoonosis. Rev. Med. Virol. 2006, 16, 5–36. [Google Scholar] [CrossRef] [PubMed]
  33. Huang, C.C.; Nguyen, D.; Fernandez, J.; Yun, K.Y.; Fry, K.E.; Bradley, D.W.; Tam, A.W.; Reyes, G.R. Molecular cloning and sequencing of the mexico isolate of hepatitis e virus (hev). Virology 1992, 191, 550–558. [Google Scholar] [CrossRef] [PubMed]
  34. Van Cuyck-Gandré, H.; Zhang, H.Y.; Tsarev, S.A.; Clements, N.J.; Cohen, S.J.; Caudill, J.D.; Buisson, Y.; Coursaget, P.; Warren, R.L.; Longer, C.F. Characterization of hepatitis e virus (hev) from algeria and chad by partial genome sequence. J. Med. Virol. 1997, 53, 340–347. [Google Scholar] [CrossRef] [PubMed]
  35. Wang, Y.; Ling, R.; Erker, J.C.; Zhang, H.; Li, H.; Desai, S.; Mushahwar, I.K.; Harrison, T.J. A divergent genotype of hepatitis e virus in chinese patients with acute hepatitis. J. Gen. Virol. 1999, 80, 169–177. [Google Scholar] [PubMed]
  36. Kwo, P.Y.; Schlauder, G.G.; Carpenter, H.A.; Murphy, P.J.; Rosenblatt, J.E.; Dawson, G.J.; Mast, E.E.; Krawczynski, K.; Balan, V. Acute hepatitis e by a new isolate acquired in the united states. Mayo Clin. Proc. Mayo Clin. 1997, 72, 1133–1136. [Google Scholar] [CrossRef]
  37. Aggarwal, R.; Naik, S. Epidemiology of hepatitis e: Current status. J. Gastroenterol. Hepatol. 2009, 24, 1484–1493. [Google Scholar] [CrossRef] [PubMed]
  38. Robert Koch-Institut. Available online: http://www.rki.de (accessed on 7 April 2015).
  39. Scharn, N.; Ganzenmueller, T.; Wenzel, J.J.; Dengler, R.; Heim, A.; Wegner, F. Guillain-barre syndrome associated with autochthonous infection by hepatitis e virus subgenotype 3c. Infection 2014, 42, 171–173. [Google Scholar] [CrossRef] [PubMed]
  40. Manka, P.; Bechmann, L.P.; Coombes, J.D.; Thodou, V.; Schlattjan, M.; Kahraman, A.; Syn, W.K.; Saner, F.; Gerken, G.; Baba, H.; et al. Hepatitis e virus infection as a possible cause of acute liver failure in europe. Clin. Gastroenterol. Hepatol. 2015. [Google Scholar] [CrossRef]
  41. Takahashi, K.; Okamoto, H.; Abe, N.; Kawakami, M.; Matsuda, H.; Mochida, S.; Sakugawa, H.; Suginoshita, Y.; Watanabe, S.; Yamamoto, K.; et al. Virulent strain of hepatitis e virus genotype 3, japan. Emerg. Infect. Dis. 2009, 15, 704–709. [Google Scholar] [CrossRef] [PubMed]
  42. Pfefferle, S.; Frickmann, H.; Gabriel, M.; Schmitz, N.; Gunther, S.; Schmidt-Chanasit, J. Fatal course of an autochthonous hepatitis e virus infection in a patient with leukemia in germany. Infection 2012, 40, 451–454. [Google Scholar] [CrossRef] [PubMed]
  43. Zehender, G.; Ebranati, E.; Lai, A.; Luzzago, C.; Paladini, S.; Tagliacarne, C.; Galli, C.; Galli, M.; Ciccozzi, M.; Zanetti, A.R.; et al. Phylogeography and phylodynamics of european genotype 3 hepatitis e virus. Infect. Genet. Evol. 2014, 25, 138–143. [Google Scholar] [CrossRef] [PubMed]
  44. Nakano, T.; Okano, H.; Kobayashi, M.; Ito, K.; Ohmori, S.; Nomura, T.; Kato, H.; Ayada, M.; Nakano, Y.; Akachi, S.; et al. Molecular epidemiology and genetic history of european-type genotype 3 hepatitis e virus indigenized in the central region of japan. Infect. Genet. Evol. 2012, 12, 1524–1534. [Google Scholar] [CrossRef] [PubMed]
  45. Nakano, T.; Takahashi, K.; Pybus, O.G.; Hashimoto, N.; Kato, H.; Okano, H.; Kobayashi, M.; Fujita, N.; Shiraki, K.; Takei, Y.; et al. New findings regarding the epidemic history and population dynamics of japan-indigenous genotype 3 hepatitis e virus inferred by molecular evolution. Liver Int. 2012, 32, 675–688. [Google Scholar] [CrossRef] [PubMed]
  46. Nakano, T.; Takahashi, K.; Arai, M.; Okano, H.; Kato, H.; Ayada, M.; Okamoto, H.; Mishiro, S. Identification of european-type hepatitis e virus subtype 3e isolates in Japanese wild boars: Molecular tracing of hev from swine to wild boars. Infect. Genet. Evol. 2013, 18, 287–298. [Google Scholar] [CrossRef] [PubMed]
  47. Colson, P.; Borentain, P.; Motte, A.; Lagrange, X.; Kaba, M.; Henry, M.; Tamalet, C.; Gerolami, R. First human cases of hepatitis e infection with genotype 3c strains. J. Clin. Virol. 2007, 40, 318–320. [Google Scholar] [CrossRef] [PubMed]
  48. Kamar, N.; Dalton, H.R.; Abravanel, F.; Izopet, J. Hepatitis e virus infection. Clin. Microbiol. Rev. 2014, 27, 116–138. [Google Scholar] [CrossRef] [PubMed]
  49. Hoffmann, B.; Depner, K.; Schirrmeier, H.; Beer, M. A universal heterologous internal control system for duplex real-time rt-pcr assays used in a detection system for pestiviruses. J. Virol. Methods 2006, 136, 200–209. [Google Scholar] [CrossRef] [PubMed]
  50. Hall, T.A. Bioedit: A user-friendly biological sequence alignment editor and analysis program for windows 95/98/nt. Nucl. Acids Symp. Ser. 1999, 41, 95–98. [Google Scholar]
  51. Tamura, K.; Stecher, G.; Peterson, D.; Filipski, A.; Kumar, S. Mega6: Molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evol. 2013, 30, 2725–2729. [Google Scholar] [CrossRef] [PubMed]
  52. Schlosser, J.; Eiden, M.; Vina-Rodriguez, A.; Fast, C.; Dremsek, P.; Lange, E.; Ulrich, R.G.; Groschup, M.H. Natural and experimental hepatitis e virus genotype 3-infection in european wild boar is transmissible to domestic pigs. Vet. Res. 2014, 45, e121. [Google Scholar] [CrossRef]
  53. Kimura, M. A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J. Mol. Evol. 1980, 16, 111–120. [Google Scholar] [CrossRef] [PubMed]
  54. Norder, H.; Sundqvist, L.; Magnusson, L.; Ostergaard Breum, S.; Löfdahl, M.; Larsen, L.E.; Hjulsager, C.K.; Magnius, L.; Böttiger, B.E.; Widen, F. Endemic hepatitis e in two nordic countries. Eurosurveillance 2009, 14, 1–9. [Google Scholar]
  55. Jothikumar, N.; Cromeans, T.L.; Robertson, B.H.; Meng, X.J.; Hill, V.R. A broadly reactive one-step real-time RT-PCR assay for rapid and sensitive detection of hepatitis e virus. J. Virol. Methods 2006, 131, 65–71. [Google Scholar] [CrossRef] [PubMed]
  56. Rokas, A.; Williams, B.L.; King, N.; Carroll, S.B. Genome-scale approaches to resolving incongruence in molecular phylogenies. Nature 2003, 425, 798–804. [Google Scholar] [CrossRef] [PubMed]
  57. Pei, Y.; Yoo, D. Genetic characterization and sequence heterogeneity of a canadian isolate of swine hepatitis e virus. J. Clin. Microbiol. 2002, 40, 4021–4029. [Google Scholar] [CrossRef] [PubMed]
  58. Simmonds, P.; Bukh, J.; Combet, C.; Deleage, G.; Enomoto, N.; Feinstone, S.; Halfon, P.; Inchauspe, G.; Kuiken, C.; Maertens, G.; et al. Consensus proposals for a unified system of nomenclature of hepatitis c virus genotypes. Hepatology 2005, 42, 962–973. [Google Scholar] [CrossRef] [PubMed]
  59. Johne, R.; Plenge-Bönig, A.; Hess, M.; Ulrich, R.G.; Reetz, J.; Schielke, A. Detection of a novel hepatitis e-like virus in faeces of wild rats using a nested broad-spectrum RT-PCR. J. Gen. Virol. 2010, 91, 750–758. [Google Scholar] [CrossRef] [PubMed]
  60. Purdy, M.A. Evolution of the hepatitis e virus polyproline region: Order from disorder. J. Virol. 2012, 86, 10186–10193. [Google Scholar] [CrossRef] [PubMed]
  61. Lhomme, S.; Abravanel, F.; Dubois, M.; Sandres-Saune, K.; Mansuy, J.M.; Rostaing, L.; Kamar, N.; Izopet, J. Characterization of the polyproline region of the hepatitis e virus in immunocompromised patients. J. Virol. 2014, 88, 12017–12025. [Google Scholar] [CrossRef] [PubMed]
  62. Yang, Z. Computational Molecular Evolution; Oxford University Press: Oxford, UK, 2006; pp. xvi, p. 357. [Google Scholar]
  63. Wichmann, O.; Schimanski, S.; Koch, J.; Kohler, M.; Rothe, C.; Plentz, A.; Jilg, W.; Stark, K. Phylogenetic and case-control study on hepatitis e virus infection in germany. J. Infect. Dis. 2008, 198, 1732–1741. [Google Scholar] [CrossRef] [PubMed]
  64. Lopes Dos Santos, D.R.; Lewis-Ximenez, L.L.; da Silva, M.F.; de Sousa, P.S.; Gaspar, A.M.; Pinto, M.A. First report of a human autochthonous hepatitis e virus infection in brazil. J. Clin. Virol. 2010, 47, 276–279. [Google Scholar] [CrossRef] [PubMed]
  65. Legrand-Abravanel, F.; Mansuy, J.M.; Dubois, M.; Kamar, N.; Peron, J.M.; Rostaing, L.; Izopet, J. Hepatitis e virus genotype 3 diversity, france. Emerg. Infect. Dis. 2009, 15, 110–114. [Google Scholar] [CrossRef] [PubMed]

Share and Cite

MDPI and ACS Style

Vina-Rodriguez, A.; Schlosser, J.; Becher, D.; Kaden, V.; Groschup, M.H.; Eiden, M. Hepatitis E Virus Genotype 3 Diversity: Phylogenetic Analysis and Presence of Subtype 3b in Wild Boar in Europe. Viruses 2015, 7, 2704-2726. https://doi.org/10.3390/v7052704

AMA Style

Vina-Rodriguez A, Schlosser J, Becher D, Kaden V, Groschup MH, Eiden M. Hepatitis E Virus Genotype 3 Diversity: Phylogenetic Analysis and Presence of Subtype 3b in Wild Boar in Europe. Viruses. 2015; 7(5):2704-2726. https://doi.org/10.3390/v7052704

Chicago/Turabian Style

Vina-Rodriguez, Ariel, Josephine Schlosser, Dietmar Becher, Volker Kaden, Martin H. Groschup, and Martin Eiden. 2015. "Hepatitis E Virus Genotype 3 Diversity: Phylogenetic Analysis and Presence of Subtype 3b in Wild Boar in Europe" Viruses 7, no. 5: 2704-2726. https://doi.org/10.3390/v7052704

APA Style

Vina-Rodriguez, A., Schlosser, J., Becher, D., Kaden, V., Groschup, M. H., & Eiden, M. (2015). Hepatitis E Virus Genotype 3 Diversity: Phylogenetic Analysis and Presence of Subtype 3b in Wild Boar in Europe. Viruses, 7(5), 2704-2726. https://doi.org/10.3390/v7052704

Article Metrics

Back to TopTop