Sample collection
A total of 388 rectal swabs were collected from pigs (Sus scrofa domestica) of different production stages (different age categories) from 25 pig farms located in four districts of Slovakia. These farms housed from ≤ 100 to ≥ 1000 animals. The majority of farms focused on two or three production stages of pigs, namely weaners, growers and fatteners (see below). Only two conventional closed farms had over 1000 pigs and all production categories including sows and piglets. Nine farms had less than 100 animals.
Samples were collected using swab applicators (Sarstedt AG & Co, Germany) from five production stages: suckling piglets (≤ 4 weeks, n=53), weaners (5–10 weeks, n=81), growers (11–16 weeks, n=67), fatteners (≥17 weeks, n=135) and sows (1–3 years, n=52). The health status of all pigs was evaluated by qualified veterinarians on each farm. All pigs were asymptomatic with no clinical signs observed.
HEV RNA extraction and reverse transcription
Rectal swabs were eluted in 1 ml of 0.01 mol/l PBS (Merck Millipore Corp., USA) for 30 min, vortexed at 2000 rev. min-1 for 3 min and centrifuged at 14, 000 x g for 5 min. Total RNA was extracted from 140 µl supernatant using the QIAamp® viral RNA mini kit (QIAGEN GmbH, Hilden, Germany) by robotic station (QIAcube GmbH, Hilden, Germany) according to the manufacturer’s instruction. Supernatants and extracted RNA were stored at -80°C. The cDNA was synthesized in a 20 µl reaction mix comprising 5 µl of extracted RNA, 0.5 mM dNTPs (Thermo Fisher Scientific, Inc., USA), 200 U RevertAid Premium reverse transcriptase with 1xRT buffer (Thermo Fisher Scientific, Inc., USA), 5 µM of gene specific reverse outer primers [17, 27] (Microsynth Austria, GmbH, Austria), 20 U RNase inhibitor (Invitrogen, Inc., USA) and molecular biology grade water (Merck, GmbH, Germany). The mix was incubated at 65°C for 5 min then chilled on ice. Subsequently, the mix was incubated at 50°C for 30 min to synthesise cDNA and finally at 85°C for 5 min to terminate the reaction.
Nested RT-PCR and sequencing
The detection of HEV was based on the amplification of a 287 bp fragment of methyltranspherase (MTase) in ORF1 using outer and inner primers [27] and a 348 bp fragment of capsid protein in ORF2 using primers by Meng et al. [5]. The PCR reaction mix (50 µl) was composed of 1x ThermoPol reaction buffer (New England Biolabs, Inc., USA), 0.2 mM dNTPs (Thermo Fisher Scientific, Inc., USA), 300 nM of outer primers, 1 U Taq DNA polymerase (New England Biolabs, Inc., USA), 4 µl cDNA and molecular biology grade water (Merck, GmbH, Germany). The first PCR was carried out under the following thermal profile: 1 cycle at 95°C for 1 min, and 35 cycles with denaturation at 95°C for 30 s, annealing at 55°C for 1 min, extension at 68°C for 1 min and final extension at 68°C for 5 min using Thermocycler C1000 (Bio-Rad Laboratories, Inc., USA). For the second PCR with inner primers a similar thermal profile was used. The size of PCR products was checked by electrophoresis in 2% agarose gel after staining with GelRedTM (Biotium, Inc., USA) and visualized by Gel Doc EZ imager (Bio-Rad Laboratories, Inc., USA). PCR products with the expected size of 287 bp and 348 bp were purified and sequenced by the Sanger method with the PCR primers by a commercial company (Microsynth Austria, GmbH, Austria).
Phylogenetic analysis of HEV sequences
Partial ORF1 (242 bp) and ORF2 (304 bp) sequences (primers were omitted) were edited and aligned using the programmes SeqMan, EditSeq and MegAlign (Lasergene, DNASTAR, Inc. USA). Sequences were first checked against the NCBI GenBank database using nucleotide BLAST (http://blast.ncbi.nlm.nih.gov/Blast.cgi).
Two representative phylogenetic trees were constructed based on the partial ORF1 and ORF2 nucleotide sequences. The model tests and p-distance calculations were performed by MEGA6 [28]. Maximum Likelihood phylogenetic analysis of partial ORF1 gene using the GTR+G+I (General time reversible model with Gamma distribution plus evolutionarily Invariable sites) model was used. For the ORF2 gene the K2+G+I (Kimura-2 parameter model with Gamma distribution plus evolutionarily Invariable sites) model was employed. Models with the lowest BIC scores (Bayesian Information Criterion) were used for phylogenetic analysis. The bootstrap support values of branches were calculated from 1000 replicates. All 45 Slovak HEV nucleotide sequences obtained in this study were submitted to NCBI GenBank database under accession numbers: MT408248-408292.
Statistical analysis
The incidence of HEV was statistically analysed by production categories (suckling piglets, weaners, growers, fatteners, sows) using the chi-square (χ2) test with confidence limits of 95%, P < 0.05 (statistically significant) or 99%, P < 0.01 (highly statistically significant) using GraphPad Prism 5 for Windows (GraphPad Software, USA).