Next Article in Journal
Dengue Surveillance System in Brazil: A Qualitative Study in the Federal District
Next Article in Special Issue
Legionella Detection in Water Networks as per ISO 11731:2017: Can Different Filter Pore Sizes and Direct Placement on Culture Media Influence Laboratory Results?
Previous Article in Journal
Do the Walkability and Urban Leisure Amenities of Neighborhoods Affect the Body Mass Index of Individuals? Based on a Case Study in Seoul, South Korea
Previous Article in Special Issue
Two Drinking Water Outbreaks Caused by Wastewater Intrusion Including Sapovirus in Finland
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nine-Year Nationwide Environmental Surveillance of Hepatitis E Virus in Urban Wastewaters in Italy (2011–2019)

by
Marcello Iaconelli
1,
Giusy Bonanno Ferraro
1,
Pamela Mancini
1,
Elisabetta Suffredini
2,
Carolina Veneri
1,
Anna Rita Ciccaglione
3,
Roberto Bruni
3,
Simonetta Della Libera
1,
Francesco Bignami
1,
Massimo Brambilla
4,
Dario De Medici
2,
David Brandtner
5,
Pietro Schembri
6,
Stefania D’Amato
7 and
Giuseppina La Rosa
1,*
1
Department of Environment and Health - Istituto Superiore di Sanità, 00161 Rome, Italy
2
Department of Food Safety, Nutrition and Veterinary Public Health, Istituto Superiore di Sanità, 00161 Rome, Italy
3
Department Infectious Diseases, Istituto Superiore di Sanità, 00161 Rome, Italy
4
Consiglio per la ricerca in agricoltura e l’analisi dell’economia agraria (CREA), Research Centre for Engineering and Agri Food Processing, 24047 Treviglio, BG, Italy
5
Indipendent Researcher, 00156 Rome, Italy
6
Regional Department for Health Activities and Epidemiological Observatory of the Sicilian Region, 90146 Palermo, Italy
7
Ministry of Health, Directorate-General for Prevention, 00144 Rome, Italy
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2020, 17(6), 2059; https://doi.org/10.3390/ijerph17062059
Submission received: 27 February 2020 / Revised: 18 March 2020 / Accepted: 18 March 2020 / Published: 20 March 2020

Abstract

:
Hepatitis E virus (HEV) is an emerging causative agent of acute hepatitis worldwide. To provide insights into the epidemiology of HEV in Italy, a large-scale investigation was conducted into urban sewage over nine years (2011–2019), collecting 1374 sewage samples from 48 wastewater treatment plants located in all the 20 regions of Italy. Broadly reactive primers targeting the ORF1 and ORF2 regions were used for the detection and typing of HEV, followed by Sanger and next generation sequencing (NGS). Real-time RT-qPCR was also used to attempt quantification of positive samples. HEV RNA detection occurred in 74 urban sewage samples (5.4%), with a statistically significant higher frequency (7.1%) in central Italy. Fifty-six samples were characterized as G3 strains and 18 as G1. While the detection of G3 strains occurred in all the surveillance period, G1 strains were mainly detected in 2011–2012, and never in 2017–2019. Typing was achieved in 2 samples (3f subtype). Viral concentrations in quantifiable samples ranged from 1.2 × 103 g.c./L to 2.8 × 104 g.c./L. Our results suggest the considerable circulation of the virus in the Italian population, despite a relatively small number of notified cases, a higher occurrence in central Italy, and a noteworthy predominance of G3 strains.

1. Introduction

Hepatitis E virus (HEV) is one of the major causes of acute hepatitis in the world, commonly transmitted by the fecal–oral route through contaminated water and food. The World Health Organization (WHO) estimates that every year 20 million hepatitis E infections occur, resulting in over 3 million symptomatic cases (http://www.who.int/mediacentre/factsheets/fs280/en/).
HEV infection develops in most individuals as self-limiting acute hepatitis, with mortality rates of around 1%. However, some individuals may develop fulminant hepatic failure, a severe condition that is frequently fatal. This complication is particularly common when the infection occurs in pregnant women, where mortality rates rise to up to 25% [1]. Furthermore, immunocompromised patients may progress to chronic hepatitis following an HEV infection [2].
Large-scale outbreaks of acute hepatitis occur in regions where HEV is endemic, such as Africa, Asia, and the Middle East. In such countries, HEV is primarily a waterborne illness associated with large epidemics due to contaminated water and water supplies and poor sanitation conditions [1,3]. In contrast, hepatitis E occurs as sporadic cases in Europe, Japan, and the USA, where HEV is spread mainly by zoonotic transmission, and domestic pigs and wild boars are considered main reservoirs [4,5]. Foodborne HEV is most commonly transmitted to humans consuming raw or undercooked pork products [6]. However, virus transmission can also follow the consumption of the meat of animals such as wild boar, deer, rabbits, and camels [7,8] as well as professional exposure to pigs and pork [9]. In low endemicity countries, hepatitis E was considered initially as a travel-associated disease, but in the last few years, an increasing number of autochthonous cases have been reported [8]. In Europe, the incidence of confirmed cases has steadily increased over the last decade [6]. The recognition of a high anti-HEV seroprevalence in a significant proportion of healthy individuals from nonendemic countries suggests widespread infection, which most likely occurs at a subclinical or asymptomatic level [10,11].
HEV is a member of the family Hepeviridae, genus Orthohepevirus, which is classified into four species: Orthohepevirus A–D [12,13]. All the human HEV genotypes (G1 to G4) belong to Orthohepevirus A, along with other genotypes (G5 and G6) that were isolated from wild boar [12], and G7 HEV isolated from camelids and a patient with chronic hepatitis [14]. Genotypes 1 (G1) and 2 (G2) are mainly restricted to humans, while genotypes 3 (G3) and 4 (G4) are isolated in both humans and animals (pig, deer, wild boar, rabbit) in developing and industrialized countries. Genotypes 1 (detected in Asia and Africa), and 2 (detected in Mexico and Africa), include strains from areas where the virus is endemic. G3 strain isolation occurs mostly in North America, Europe, South America, and Japan [15]. Genotype 4 has mainly been identified in China, Taiwan, Japan, and Vietnam, but a few cases have also been reported in Europe, including Italy [16,17,18].
The purpose of this study was to assess the occurrence and diversity of HEV in the Italian population through molecular analysis of urban sewage samples collected at wastewater treatment plants (WTPs). Environmental surveillance through the monitoring of sewer systems has been demonstrated as a useful auxiliary tool in the study of the epidemiology of viral pathogens, and an essential component of early warning systems. Indeed, a number of studies on the occurrence of HEV in urban sewages have been conducted in European countries, like Norway, Portugal, France, Spain, Sweden, Switzerland, UK, Greece, and Italy [19,20,21,22,23,24,25,26,27,28]. In the present study, however, a significant effort to provide a robust estimate of HEV circulation at a country level was undertaken through the application of a large-scale sampling scheme (approximately 1400 samples), over an extended temporal (nine years, from 2011 to 2019), and spatial scale (48 wastewater treatment plants distributed over the 20 regions of the country).
This allowed us to: 1) monitor HEV trends over time, both in terms of virus occurrence, and of circulating genotypes, comparing these trends with the data reported for notified human cases; and 2) assess the differences in the geographical occurrence of HEV, and identify hotspots of virus excretion, correlating them with data from the clinical setting.

2. Materials and Methods

2.1. Sampling

In the period 2011 to 2019, a total of 1374 raw sewage samples were collected from 48 WTPs distributed over the entire country and covering all the 20 Italian administrative regions: 529 samples were taken from 8 Regions of Northern Italy (Emilia Romagna, Friuli Venezia Giulia, Liguria, Lombardy, Piedmont, Trentino-South Tyrol, Veneto, and Aosta Valley), 590 from 5 Regions of Central Italy (Abruzzo, Latium, Marche, Tuscany, and Umbria) and 255 from 7 Regions of Southern Italy (Apulia, Basilicata, Calabria, Campania, Molise, Sardinia, and Sicily). Figure 1 shows a geographic information system (GIS) map of the WTPs included in the study.
The environmental surveillance for HEV started within the framework of the 2011–2016 European Project PREDEMICS (Preparedness, prediction, and prevention of emerging zoonotic viruses with pandemic potential using multidisciplinary approaches - https://predemics.biomedtrain.eu/cms/). After that, the surveillance continued with the 2017–2019 national project “Hepatitis E, an emerging problem in food safety: ‘One Health’ approach for risk assessment” (CCM 2016 program) financed by the Italian Ministry of Health.
All the samples, were registered in the PostgreSQL database (available to recorded users at https://cosmos.bio.uniroma1.it) and received a univocal ID code. Before concentration, samples underwent spiking with a known quantity of murine norovirus (MNV-1), used as a control process.

2.2. Concentration and Treatment

Two different sewage treatment protocols were used during the 9-year surveillance; wastewater samples collected in 2011–2016 were analyzed, as previously described [29]. An aliquot of 20 mL of untreated wastewater (collected at the inlet collector, before treatment) were treated with 2 mL of 2.5 M glycine (pH 9.5), and incubated on ice for 30 min. The solution was then treated with 2.2 mL chloroform and centrifuged at 2300 g for 10 min, and the supernatant was recovered and measured. From 2017 onwards, with the start of the new project “Hepatitis E, an emerging problem in food safety: ‘One Health’ approach for risk assessment,” virus concentration was performed using the two-phase separation method (polyethylene glycol, (PEG)–dextran) according to the WHO guidelines for poliovirus environmental surveillance [30]. A sample of 250 mL of sewage underwent centrifugation at 4400 g for 30 min at 4 °C, and the resulting pellet was stored at 4 °C for the subsequent steps. The recovered supernatant, instead, was treated with 17 mL of 5 M sodium chloride, 142 mL of 29% polyethylene glycol 6000, and 20 mL of 22% dextran, and the mixture was transferred into a separation funnel and left overnight at 4 °C. After phase separation, about 10 mL of the lower phase, which included the interphase, was collected drop by drop into a sterile 50 mL centrifuge tube and was mixed with the previously saved pellet. Then, a 20% volume of chloroform was added to the recovered volume, thoroughly mixed by vortexing, and then centrifuged at 1000 g for 10 min at 4 °C. The supernatant was then recovered, and sample volumes recorded.
Viral nucleic acid extraction was carried out on 5 mL of chloroform-treated samples using the NucliSENS MiniMag (bioMerieux, Marcy l’Etoile, France) semi-automated extraction system, in accordance with the manufacturer’s instructions. Eluted RNA (100 μL) was stored in aliquots at −80 °C until molecular analysis.

2.3. Reverse Transcription Polymerase Chain Reaction Nested Amplification (Nested RT-PCR)

A nested RT-PCR assay, using universal primers targeting the methyltransferase (MTase) gene of the ORF1 region [31], was used. This assay, that amplifies a 172 bp region of all HEV genotypes, was proven to be the most sensitive in a previous study comparing five different protocols [32].
PCR amplification was performed in a T100 thermal cycler (BioRad, Hercules, CA, USA) using MyTaq one-step RT-PCR kit (Bioline, Memphis, TN, USA). PCR reactions were prepared in a 25 µL mixture containing 12.5 µL of one-step mix, 1 µL (10 pmol) of each primer, and 2 µL of the extracted RNA. PCR cycling conditions were: reverse transcription at 45 °C for 20 min, 95 °C for 1 min, followed by 40 cycles of 95 °C for 10 s, 51 °C for 10 s and 72 °C for 30 s, and a final step at 72 °C for 5 min. After the first round of PCR amplification, one µL of PCR product underwent the nested amplification (35 cycles), performed using MyTaq red mix kit (Bioline). The PCR cycle included 1 min of initial denaturation at 95 °C; followed by 35 cycles of 95 °C for 15 s, 48 °C for 15 s, 72 °C for 10 s with a final elongation step of 5 min at 72 °C. The PCR products were visualized on 1.5% agarose gel stained with GelRed® (Biotium, Fremont, CA, USA). Standard precautions were taken to prevent false-positive results.
All positive samples were reanalyzed using an assay targeting the ORF2 region [33,34,35], to allow a characterization of the strain, not only at the genotype level but also at the subtype level.

2.4. Sequencing and Molecular Characterization

PCR amplicons were purified using a Montage PCRm96 microwell filter plate (Millipore, Billerica, MA, USA) and were subjected to automated sequencing by BioFab Research (Rome, Italy), using the PCR amplification primers. As no ORF2 sequences could be obtained from the samples collected in 2017–2019 by Sanger sequencing, due to undetectable or unspecific PCR products, Next- Generation Sequencing (NGS) was also performed to sequence even low amounts of amplified specific fragments. All reactions of the ORF2 nested amplifications were pooled and sequenced by ultra-deep sequencing onto a MiSeq platform (Illumina, San Diego, CA, USA) as previously described [36]. The consensus sequence was constructed using Bowtie2 (Galaxy version 2.3.4.3+galaxy0, free software) as reads aligner and SSAKE (Galaxy version 0.010, MIT, Cambridge, USA) as contig assembler [37,38,39].
Phylogenetic relationships were studied by a neighbor-joining phylogenetic tree calculated on the maximum composite likelihood distances. Tree reliability was assessed through bootstrap analysis (1000 replicates) and considered significant linkage values >70%. The tree included the representative sequences belonging to genotypes 1 and 3 obtained from human and animal strains. All entries were reported with accession number, followed by host suffix (Sw, swine; Hu, human; WTP, wastewater treatment plant). Still unpublished sequences from clinical cases provided by the Italian National Reference Laboratory (NRL) for human Viral Hepatitis were reported as sample ID, followed by the year of isolation and administrative region origin.
Sequences were submitted to NCBI GenBank with the accession numbers LN624460 to LN624484, LN906799 to LN906809, and MN251630 to MN251657.

2.5. Real-Time RT-qPCR

Real-time RT-qPCR analysis was conducted as previously described [40] on samples that were positive according to nested RT-PCR. Amplifications were carried out with the RNA UltraSense one-step qRT-PCR system (LifeTechnologies, Carlsbad, CA, USA), using 1.25 μL of UltraSense enzyme mix and the following concentrations for primers and probe: 500 nM for forward primer JVHEVF (5′-GGTGGTTTCTGGGGTGAC-3′), 900 nM for reverse primer JVHEVR (5′-AGGGGTTGGTTGGATGAA-3′), and 250 nM for probe JVHEVP-MGB (5′-FAM-TGATTCTCAGCCCTTCGC-MGB-3′) [41,42]. The thermal profile was as follows: reverse transcription at 50 °C for 60 min, inactivation at 95 °C for 5 min, followed by 45 cycles of 15 s at 95 °C, 1 min at 60 °C, and 1 min at 65 °C. Analyses were performed in duplicate reactions (5 μL of RNA each) on a Quant Studio 12K flex (ThermoFisher Scientific, Waltham, MA, USA), and the average concentration of the two replicate reactions was used for quantification. Two negative PCR controls were included in each run. PCR inhibition was ruled out using an external amplification control (in vitro synthesized RNA) and amplifications were considered acceptable if inhibition was ≤50%. A linearized plasmid containing the target sequence was used to generate the standard curve (dynamic range 1 × 100 – 1 × 104 copies/μL); curves with a slope lying between −3.1 and −3.6, and a R2 ≥ 0.98 were used for quantification. Quantitative values in raw sewage (expressed in genome copies, g.c./L) were calculated based on the fraction of the sample subjected to RNA extraction compared to its total volume after sample concentration.

2.6. Statistical Analysis

Statistical analysis (chi-square test) was performed using the Minitab 17 statistical software version 2010 (Minitab Ltd., Coventry, UK).

3. Results

3.1. HEV Detection in Geographic Areas and Over Time

Results of HEV detection by nested PCR are reported in Table 1. Positive and negative controls, used to avoid the risk of false-positive or false-negative results, always gave the expected results. Overall, 74/1374 (5.4%, CI 95%; 4.3–6.7%) samples tested positive using the ORF1 assay: 23/529 (4.3%, CI 95%; 2.9–6.5%) in Northern Italy, 42/590 (7.1%, CI 95%; 5.3–9.5%) in Central Italy, and 9/255 (3.5%, CIm95%; 1.8–6.7%) in Southern Italy. A statistically significant difference (chi-square p value = 0.042) was present in HEV detection in the three geographic areas. No positive samples were detected in 7 regions (Basilicata, Campania, Liguria, Sardinia Tuscany, Umbria, and Valle d’Aosta), while in the remaining 13 regions HEV occurrence ranged from 3.0% (Veneto) to 14.2% (Molise).
Positive samples were found in all the years of the monitoring, ranging from 2.7% to 13.3% (Figure 2).

3.2. HEV Quantification in Positive Sample

Seventy-two samples underwent analysis by real-time RT-qPCR (Table 2). Thirty-two samples were below the detection limit, and most of the positive samples (31 out of 40) were below the analytical LOQ of 1.1 × 103 genome copies (g.c.)/L, with estimated concentrations ranging from ~20 to 9.9 × 102 g.c./L. Quantifiable samples ranged from 1.2 × 103 g.c./L to 2.8 × 104 g.c./L, with a median value of 3.2 × 103 g.c./L.

3.3. HEV Genotyping

Of the 74 positive samples, 56 were characterized as G3 and 18 as G1. While G3 sequences were detected in all the 9 years of monitoring (Figure 3), G1 sequences were mainly detected in 2011–2012 (13 of 18; 72%). Thereafter, only five G1 samples were found between 2013 and 2016, while in the triennium 2017–2019 only G3 sequences were detected.
A phylogenetic tree of the Mtase-ORF1 fragment obtained from the HEV strains detected in this study is shown in Figure 4A,B. The tree includes, in addition to the study sequences, clinical, environmental, and animal (swine) sequences previously detected in Italy. According to the sequencing results, the phylogenetic analyses showed that the strains were grouped into two main clusters, related to HEV G1 and G3. HEV G3 strains showed high sequence identity with human sequences detected in autochthonous cases (no travel history) and swine in Italy, and most of the G1 strains showed 100% nt identity with Italian strains detected in patients traveling in endemic areas.
An attempt to further characterize the positive samples was performed using broad-range primers targeting the ORF2 gene, but only 2 of these samples provided amplification in this region (ID sample 2222 and 2275, both collected from WTPs in Rome, in 2014 and 2015, respectively), and were characterized as subtype 3f by Sanger sequencing. NGS performed on pooled PCR amplicons collected in years 2017–2019 detected another G3f sequence.

4. Discussion

In Italy, the circulation of HEV has been demonstrated since the 1980s. Sporadic cases have been reported [18,32,35,43], related both to traveling to endemic countries and autochthonous. The latter cases were associated with consumption of game meat, dried pork liver sausage, and raw-dried wild boar sausages [44].
As for animal HEV, it has been demonstrated that HEV circulates largely within pig herds in Italy [45,46,47,48,49]. Genotype 3 is widespread among Italian pigs, even though genotype 4 was also reported in swine farms in Northern Italy [17]. HEV RNA has also been detected in pig slurry collected in Italian farms, showing a high degree of identity with strains detected in autochthonous HEV cases [50]. Moreover, different studies have demonstrated the presence of HEV RNA in liver samples from the wild boar population [40,51,52,53] and even in sheep serum and fecal samples [54].
The present study aimed to provide information on the occurrence and diversity of HEV in the Italian population through molecular analysis of urban sewage samples collected at WTPs. A number of studies have investigated HEV in urban sewages in nonendemic countries [19,23,26,28,55,56,57,58,59,60,61,62,63,64]. However, most of these studies are limited in time and conducted on a small number of samples. The present study, instead, describes a large-scale investigation, with a significant number of samples, a nationwide geographic coverage, and a nine-year observation period.
HEV RNA detection in sewage occurred throughout the years 2011–2019 (mean detection frequency 5.3%) with higher occurrence in 2014 (13.3%). In a study conducted following an HEV outbreak on a small French island, Miura et al. estimated that if 1–4% of inhabitants connected to a WTP were infected with HEV, the virus would be detectable in urban wastewater [57]. While it is not possible to directly infer these conclusions to WTPs of different capacities and size, our results on the occurrence of HEV in sewage suggests the considerable circulation of the virus in the Italian population, despite the relatively small number of notified cases. During the period from 2011 to the first semester of 2019 only 323 cases of acute viral hepatitis were reported in the Italian integrated epidemiological system for acute viral hepatitis (SEIEVA) managed by the Italian Institute of Health (ISS) (https://www.epicentro.iss.it/epatite/bollettino/Bollettino-4-marzo-2019.pdf and https://www.epicentro.iss.it/epatite/bollettino/Bollettino-n-5-novembre-2019.pdf), with a trend of autochthonous cases on the rise over the study period (7 in 2011 vs 45 in 2018). The discordance between the relatively constant HEV detection in Italian sewages and the low incidence of HEV infection in humans may be associated with nondetection of mild, subclinical, or asymptomatic cases, which, however, shed virus into local sewage systems and contributed to virus circulation. This phenomenon is recognized as the “surveillance pyramid” [65]; because of under and misdiagnoses, clinical surveillance is likely to only capture the tip of the iceberg of viral diseases, generally corresponding to the hospitalized or laboratory diagnosed cases, while the full extent of disease can only be captured through studies addressing incidence and etiology at a community level. In this sense, the monitoring of urban wastewaters provides a powerful tool for indirect assessment of the community health status.
Comparing to data available from other European countries, the percentages obtained in the present study are similar to those found in Norway and Portugal (8.0% and 6.7%, respectively) [19,20], but lower than those documented in France (22–25%) [21], Spain (13.5–43.5%) [22,23,24], Sweden (33.0%) [25], Switzerland (24.1–32.3%) [21], and UK (93.0%) [26]; on the other hand, no positive results were detected in samples from Greece [23,27]. All these studies, however, are difficult to compare because they were very different in terms of the period of sampling, sample size, sample volume and processing, and molecular methods used. In our study, during the 9-year surveillance period, a shift of the method for concentration of viruses took place, moving from the “glycine-chloroform” to the “two-phase separation” method recommended by WHO for poliovirus environmental surveillance, which is considered the gold standard for the concentration of viruses from sewage. Both methods were applied with the use of a control virus, that was added to the sample before virus concentration in order to monitor the efficiency of the entire procedure and ensure the quality of the analytical results. Indeed, no statistically significant difference was detected in the observed period based on the concentration method used (5.8% of positive samples with “glycine-chloroform” and 4.6% with the “two-phase separation”; p value = 0.33).
As for quantitative data, we found overall low HEV concentrations in raw sewage (range of concentrations ~20 to 2.8 × 104 g.c./L), with only nine samples exceeding 103 g.c./L and a median value of quantifiable samples of 3.2 × 103 g.c./L. Few studies have been performed on HEV quantitative levels in raw sewage, and comparison of available data is hindered by the diversity of analytical methods and by WTP capacity. In three studies carried out in Switzerland, Sweden, and Norway [19,21,25] most of the positive samples were below the corresponding LOQs and, overall, only two samples were quantified with concentrations of 7.8 × 104 c.g./L [21] and 4.5 × 104 c.g./L (25). Miura and colleagues [57] in their study performed during an HEV outbreak on an island, reported virus concentrations in inlet sewage up to 6.3 × 105 c.g./L after the outbreak, with values decreasing to 103 c.g./L within a month. Therefore, the quantitative data of the present study seem to confirm that, during the observation period, the detection of HEV occurred at concentrations lower than those presumably associated with relevant outbreaks in the populations.
Of all the positive samples detected, 75.6% were G3 and 24.3% G1. While G3 strains were consistently detected in the 9 years of surveillance, G1 were mainly detected in 2011 (6 samples) and 2012 (7 samples, outnumbering G3), and then sporadically, in 2013–2016. In the last three years, only G3 strains were detected. Genotype 1 is the most common cause of hepatitis E in regions with a lower socioeconomic status, where it is predominantly spread via contaminated potable water. In Italy, infections with HEV G1 have been documented in immigrants from countries at risk after returning from a stay in their country of origin, as well as in Italian people traveling from highly endemic countries [32,35,66]. The environmental data provided in this study reflect those from clinical surveillance. Indeed, of the total 285 cases documented in SEIEVA in 2011–2018 (https://www.epicentro.iss.it/epatite/bollettino/Bollettino-4-marzo-2019.pdf), 23.5% were imported (these are most probably caused by HEV G1), and 76.5% were autochthonous (caused by G3). The sporadic circulation of genotype 1, along with the predominance of genotype 3, is consistent with reports in the literature about other industrialized countries [22,23]. The predominance of G3 cases is likely to result in a higher concentration of G3 molecules in sewage samples compared to G1, leading to a predominance of G3 PCR amplification products and to the apparent absence of G1 in some years. Unfortunately, only two of the G3 positive samples could be typed by sequencing the ORF2 region. The possible explanations for the negative ORF2 detection in the vast majority of positive samples, besides the lower sensitivity of the ORF2 assay, include the possibility of RNA damage due to the long storage time before testing, despite the efforts for sample preservation. Moreover, being that the ORF2 assay is an amplification providing a longer template compared to ORF1 (>450 bps vs 172 bp), amplification from environmental samples, which often contain inactivated viruses and fragmented RNA, may be more difficult to achieve. The two positive samples were both characterized as G3f, which is also the predominant type reported in the same period in the clinical cases tested by the Italian NRL for human viral hepatitis (unpublished data). A G3f sequence was also detected by NGS in 2017–2019 pooled samples.
Geographically, HEV identification occurred in 13 Italian Regions out of 20. In the seven regions where the virus has not been identified, we cannot be sure whether such an outcome relates to the number of samples that could be collected. However, it should also be considered that the high variability of HEV occurrence in sewage, besides the kind of influent the WTP is fed with, strictly depends on the sampling volumes and the time of collection [67]. Indeed, due to the long surveillance period, few WTPs provided samples continuously over the nine years, while the vast majority provided samples in limited periods within this time, resulting in an unequal number of samples by WTP and by administrative region. Therefore, a detailed comparative analysis of the distribution of positive samples by WTP and administrative region was not feasible. Generally speaking, higher percentages of positive samples were detected in Central Italy (a difference found to be statistically significant) and, in particular, in Latium, Marche, and Abruzzo. These data were in line with clinical data; indeed, according to data from the Italian NRL for human viral hepatitis (ISS, Rome), 67% of the confirmed cases (IgM+ and/or HEV RNA +) of HEV in Italy are from Latium and Abruzzo (unpublished data).
Serological data fosters the marked heterogeneity of seroprevalence in the different Italian geographic areas. As a matter of fact, rates in the general population range from 0.12% to 17.8% [68,69,70,71,72,73]. The different sensitivity of the anti-HEV detection systems used in different studies (the older ones using kits with lower sensitivity) partially accounts for the observed heterogeneity. However, recent testing on 10,011 blood donors throughout Italy, using the same highly sensitive detection kit, showed that a true geographical prevalence heterogeneity exists, with very high rates of infection in some areas of Central Italy [74]. Those data confirmed a previous study whose results pointed at Abruzzo as a hyper-endemic region, with a seroprevalence rate of 49% among blood donors [75]. Abruzzo and Latium can be therefore considered hot spot regions for HEV circulation. It may be relevant that these administrative regions have traditional dishes based on raw/undercooked pork products, in particular liver sausages that, both raw and dried, may contain HEV [76]. Recently, an outbreak occurred in the bordering region Marche, linked to the consumption of a locally produced raw pork liver sausage [77].
The occurrence of HEV in urban wastewater raises the question of whether sewage-contaminated surface waters can be responsible for HEV waterborne transmission. A recent review by Fenaux et al. [67] highlighted the need to investigate the possible transmission of the hepatitis E virus by water in high-income countries. Indeed, HEV has been detected in different water matrices even in these countries; moreover, HEV survival in water is similar to other enteric viruses since it seems that the virus has similar resistance to inactivating factors like heat, UV, and chlorine [67]. In Italy, HEV has been recently detected in river waters receiving wastewater discharges [78]. Moreover, HEV strains have been isolated from shellfish, marine waters, and underwater sewage discharges [79,80].

5. Conclusions

This surveillance of human wastewaters allowed us to generate data on the circulation of HEV in the general population in Italy. HEV RNA, genotype 1 and 3, were detected at the WTP entrance, with the highest rates reported from the central regions of Italy. Viral sequences were genetically close to those reported in clinical specimens. A marked predominance of G3 strains was detected as expected, however, G1 strains were also identified, except in years 2016–2019, pointing to the circulation of this genotype also in nonendemic areas.

Author Contributions

Conceptualization, G.L.R.; study methodology, G.L.R. and M.I.; samples analysis, G.B.F., P.M., C.V., F.B., S.D.L., and E.S.; analysis, D.B.; statistical analysis, M.B.; comparison clinical/environmental molecular data, A.R.C. and R.B.; writing, G.L.R., E.S., and M.I.; project administration and management, D.D.M., P.S., and S.D. All authors have read and agreed to the published version of the manuscript.

Funding

The research leading to these results has received partial funding from the European Union’s Seventh Framework Programme for research, technological development, and demonstration under Grant Agreement no. 278433-PREDEMICS, and partially from the project CCM 2016, “Hepatitis E, an emerging problem in food safety: ‘One Health’ approach for risk assessment.”.

Acknowledgments

We thank Herbert W. Virgin, Washington University (St. Louis, Missouri, United States) for providing the MNV strain used as sample process control. For urban wastewater sample collection, we gratefully acknowledge the assistance of the following members of the Italian Wastewater Network: Meris Poli (Eco Center AG, Bolzano); Margherita De Ceglia, Marco Simonetti (SMAT, Castiglione T.se); Francesca Pizza, Marco Blazina (MM SpA, Milano); Patrizia Albertelli (ARPAE, Bologna); Alesandro Orsini, Giuseppe Micciché, Tania Giachini, Federico Barducci (ACEA SpA, Roma); Giuseppe Gala (ACA SpA, Pescara); Paolo Di Marzio, Piero Coletti (Consorzio Bonifica Centro, Chieti); Jlenia Cupaiolo (ARAP Abruzzo, Avezzano, Sulmona, Vasto, San Salvo), Maria Teresa Gallo (ARPAC, Napoli); Giuseppe Siragusa, Antonella Cicala (AMAP, Palermo).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Perez-Gracia, M.T.; Garcia, M.; Suay, B.; Mateos-Lindemann, M.L. Current Knowledge on Hepatitis E. J. Clin. Transl. Hepatol. 2015, 3, 117–126. [Google Scholar]
  2. Mohamed, O.E.; Jones, J.; Osman, H.; Huissoon, A.P. Unexplained abnormal liver function in patients with primary antibody deficiency: Could it be chronic hepatitis E infection? J. Clin. Pathol. 2018, 71, 81–83. [Google Scholar] [CrossRef]
  3. Purcell, R.H.; Emerson, S.U. Hepatitis E: An emerging awareness of an old disease. J. Hepatol. 2008, 48, 494–503. [Google Scholar] [CrossRef]
  4. Caballero-Gómez, J.; Jiménez-Ruiz, S.; Lopez-Lopez, P.; Vicente, J.; Risalde, M.A.; Cano-Terriza, D.; Frias, M.; Barasona, J.A.; Rivero, A.; García-Bocanegra, I.; et al. Emergent subtype of hepatitis E virus genotype 3 in wild boar in Spain. Transbound. Emerg. Dis. 2019, 66, 1803–1808. [Google Scholar] [CrossRef]
  5. Montone, A.M.; De Sabato, L.; Suffredini, E.; Alise, M.; Zaccherini, A.; Volzone, P.; Di Maro, O.; Neola, B.; Capuano, F.; Di Bartolo, I. Occurrence of HEV-RNA in Italian Regional Pork and Wild Boar Food Products. Food Environ. Virol. 2019, 11, 420–426. [Google Scholar] [CrossRef]
  6. EFSA. BIOHAZ Panel (EFSA Panel on Biological Hazards). Scientific Opinion on the public health risks associated with hepatitis E virus (HEV) as a food-borne pathogen. 2017, 15, 4886–4889. [Google Scholar]
  7. Sridhar, S.; Teng, J.L.L.; Chiu, T.H.; Lau, S.K.P.; Woo, P.C.Y. Hepatitis E Virus Genotypes and Evolution: Emergence of Camel Hepatitis E Variants. Int. J. Mol. Sci. 2017, 18, 869. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Ruggeri, F.M.; Di Bartolo, I.; Ponterio, E.; Angeloni, G.; Trevisani, M.; Ostanello, F. Zoonotic transmission of hepatitis E virus in industrialized countries. New Microbiol. 2013, 36, 331–344. [Google Scholar] [PubMed]
  9. Hoan, N.X.; Huy, P.X.; Sy, B.T.; Meyer, C.G.; Son, T.V.; Binh, M.T.; Giang, D.P.; Tu Anh, D.; Bock, C.T.; Wang, B.; et al. High Hepatitis E virus (HEV) Positivity Among Domestic Pigs and Risk of HEV Infection of Individuals Occupationally Exposed to Pigs and Pork Meat in Hanoi, Vietnam. Open Forum Infect. Dis. 2019, 6, ofz306. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Aggarwal, R. Hepatitis E: Historical, contemporary and future perspectives. J. Gastroenterol. Hepatol. 2011, 26 (Suppl. 1), 72–82. [Google Scholar] [CrossRef]
  11. Emerson, S.U.; Purcell, R.H. Hepatitis E virus. Rev. Med. Virol. 2003, 13, 145–154. [Google Scholar] [CrossRef] [PubMed]
  12. Smith, D.B.; Simmonds, P.; Izopet, J.; Oliveira-Filho, E.F.; Ulrich, R.G.; Johne, R.; Koenig, M.; Jameel, S.; Harrison, T.J.; Meng, X.J.; et al. Proposed reference sequences for hepatitis E virus subtypes. J. Gen. Virol. 2016, 97, 537–542. [Google Scholar] [CrossRef] [PubMed]
  13. Purdy, M.A.; Harrison, T.J.; Jameel, S.; Meng, X.J.; Okamoto, H.; Van der Poel, W.H.M.; Smith, D.B.; Ictv Report Consortium. ICTV Virus Taxonomy Profile: Hepeviridae. J. Gen. Virol. 2017, 98, 2645–2646. [Google Scholar] [CrossRef] [PubMed]
  14. Lee, G.H.; Tan, B.H.; Teo, E.C.; Lim, S.G.; Dan, Y.Y.; Wee, A.; Aw, P.P.; Zhu, Y.; Hibberd, M.L.; Tan, C.K.; 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–357. [Google Scholar] [CrossRef] [Green Version]
  15. Blasco-Perrin, H.; Abravanel, F.; Blasco-Baque, V.; Peron, J.M. Hepatitis E, the neglected one. Liver Int. 2016, 36 (Suppl. 1), 130–134. [Google Scholar] [CrossRef]
  16. Zhang, W.; Shen, Q.; Mou, J.; Gong, G.; Yang, Z.; Cui, L.; Zhu, J.; Ju, G.; Hua, X. Hepatitis E virus infection among domestic animals in eastern China. Zoonoses Public Health 2008, 55, 291–298. [Google Scholar] [CrossRef]
  17. Monne, I.; Ceglie, L.; Di Martino, G.; Natale, A.; Zamprogna, S.; Morreale, A.; Rampazzo, E.; Cattoli, G.; Bonfanti, L. Hepatitis E virus genotype 4 in a pig farm, Italy, 2013. Epidemiol. Infect. 2015, 143, 529–533. [Google Scholar] [CrossRef]
  18. Garbuglia, A.R.; Scognamiglio, P.; Petrosillo, N.; Mastroianni, C.M.; Sordillo, P.; Gentile, D.; La Scala, P.; Girardi, E.; Capobianchi, M.R. Hepatitis E virus genotype 4 outbreak, Italy, 2011. Emerg. Infect. Dis. 2013, 19, 110–114. [Google Scholar] [CrossRef] [Green Version]
  19. Myrmel, M.; Lange, H.; Rimstad, E. A 1-Year Quantitative Survey of Noro-, Adeno-, Human Boca-, and Hepatitis E Viruses in Raw and Secondarily Treated Sewage from Two Plants in Norway. Food Environ. Virol. 2015, 7, 213–223. [Google Scholar] [CrossRef]
  20. Matos, A.; Mesquita, J.; Gonçalves, D.; Joana Abreu-Silva Luxo, C.; Nascimento, M. First detection and molecular characterization of hepatitis E virus in water from wastewater treatment plants in Portugal. Ann. Agric. Environ. Med. 2018, 25, 364–367. [Google Scholar] [CrossRef] [Green Version]
  21. Masclaux, F.G.; Hotz, P.; Friedli, D.; Savova-Bianchi, D.; Oppliger, A. High occurrence of hepatitis E virus in samples from wastewater treatment plants in Switzerland and comparison with other enteric viruses. Water Res. 2013, 47, 5101–5109. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Rodriguez-Manzano, J.; Miagostovich, M.; Hundesa, A.; Clemente-Casares, P.; Carratala, A.; Buti, M.; Jardi, R.; Girones, R. Analysis of the evolution in the circulation of HAV and HEV in eastern Spain by testing urban sewage samples. J. Water Health 2010, 8, 346–354. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Clemente-Casares, P.; Pina, S.; Buti, M.; Jardi, R.; Martin, M.; Bofill-Mas, S.; Girones, R. Hepatitis E virus epidemiology in industrialized countries. Emerg. Infect. Dis. 2003, 9, 448–454. [Google Scholar] [CrossRef] [PubMed]
  24. Rusinol, M.; Fernandez-Cassi, X.; Timoneda, N.; Carratala, A.; Abril, J.F.; Silvera, C.; Figueras, M.J.; Gelati, E.; Rodo, X.; Kay, D.; et al. Evidence of viral dissemination and seasonality in a Mediterranean river catchment: Implications for water pollution management. J. Environ. Manag. 2015, 159, 58–67. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Wang, H.; Sikora, P.; Rutgersson, C.; Lindh, M.; Brodin, T.; Bjorlenius, B.; Larsson, D.G.J.; Norder, H. Differential removal of human pathogenic viruses from sewage by conventional and ozone treatments. Int. J. Hyg. Environ. Health 2018, 221, 479–488. [Google Scholar] [CrossRef] [PubMed]
  26. Smith, D.B.; Paddy, J.O.; Simmonds, P. The use of human sewage screening for community surveillance of hepatitis E virus in the UK. J. Med. Virol. 2016, 88, 915–918. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Kokkinos, P.A.; Ziros, P.G.; Mpalasopoulou, A.; Galanis, A.; Vantarakis, A. Molecular detection of multiple viral targets in untreated urban sewage from Greece. Virol. J. 2011, 8, 195. [Google Scholar] [CrossRef] [Green Version]
  28. La Rosa, G.; Pourshaban, M.; Iaconelli, M.; Vennarucci, V.S.; Muscillo, M. Molecular detection of hepatitis E virus in sewage samples. Appl. Environ. Microbiol. 2010, 76, 5870–5873. [Google Scholar] [CrossRef] [Green Version]
  29. La Rosa, G.; Della Libera, S.; Iaconelli, M.; Ciccaglione, A.R.; Bruni, R.; Taffon, S.; Equestre, M.; Alfonsi, V.; Rizzo, C.; Tosti, M.E.; et al. Surveillance of hepatitis A virus in urban sewages and comparison with cases notified in the course of an outbreak, Italy 2013. BMC Infect. Dis. 2014, 14, 419. [Google Scholar] [CrossRef] [Green Version]
  30. World Health Organization. Guidelines for Environmental Surveillance of Poliovirus Circulation. World Health Organization. 2003. Available online: https://apps.who.int/iris/handle/10665/67854 (accessed on 27 February 2020).
  31. Fogeda, M.; Avellon, A.; Cilla, C.G.; Echevarria, J.M. Imported and autochthonous hepatitis E virus strains in Spain. J. Med. Virol. 2009, 81, 1743–1749. [Google Scholar] [CrossRef]
  32. La Rosa, G.; Fratini, M.; Muscillo, M.; Iaconelli, M.; Taffon, S.; Equestre, M.; Chionne, P.; Madonna, E.; Pisani, G.; Bruni, R.; et al. Molecular characterisation of human hepatitis E virus from Italy: Comparative analysis of five reverse transcription-PCR assays. Virol. J. 2014, 11, 72. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Okamoto, H.; Takahashi, M.; Nishizawa, T.; Fukai, K.; Muramatsu, U.; Yoshikawa, A. Analysis of the complete genome of indigenous swine hepatitis E virus isolated in Japan. Biochem. Biophys. Res. Commun. 2001, 289, 929–936. [Google Scholar] [CrossRef] [PubMed]
  34. Shrestha, S.M.; Shrestha, S.; Tsuda, F.; Nishizawa, T.; Gotanda, Y.; Takeda, N.; Okamoto, H. Molecular investigation of hepatitis E virus infection in patients with acute hepatitis in Kathmandu, Nepal. J. Med. Virol. 2003, 69, 207–214. [Google Scholar] [CrossRef] [PubMed]
  35. La Rosa, G.; Muscillo, M.; Vennarucci, V.S.; Garbuglia, A.R.; La, S.P.; Capobianchi, M.R. Hepatitis E Virus In Italy: Molecular Analysis Of Travel-Related And Autochthonous Cases. J. Gen. Virol. 2011, 92, 1617–1626. [Google Scholar] [CrossRef] [PubMed]
  36. Iaconelli, M.; Valdazo-Gonzalez, B.; Equestre, M.; Ciccaglione, A.R.; Marcantonio, C.; Della Libera, S.; La Rosa, G. Molecular characterization of human adenoviruses in urban wastewaters using next generation and Sanger sequencing. Water Res. 2017, 121, 240–247. [Google Scholar] [CrossRef]
  37. Cuccuru, G.; Orsini, M.; Pinna, A.; Sbardellati, A.; Soranzo, N.; Travaglione, A.; Uva, P.; Zanetti, G.; Fotia, G. Orione, a web-based framework for NGS analysis in microbiology. Bioinformatics 2014, 30, 1928–1929. [Google Scholar] [CrossRef] [Green Version]
  38. Langmead, B.; Salzberg, S.L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 2012, 9, 357–359. [Google Scholar] [CrossRef] [Green Version]
  39. Warren, R.L.; Sutton, G.G.; Jones, S.J.; Holt, R.A. Assembling millions of short DNA sequences using SSAKE. Bioinformatics 2007, 23, 500–501. [Google Scholar] [CrossRef] [Green Version]
  40. Di Pasquale, S.; De Santis, P.; La Rosa, G.; Di Domenico, K.; Iaconelli, M.; Micarelli, G.; Martini, E.; Bilei, S.; De Medici, D.; Suffredini, E. Quantification and genetic diversity of Hepatitis E virus in wild boar (Sus scrofa) hunted for domestic consumption in Central Italy. Food Microbiol. 2019, 82, 194–201. [Google Scholar] [CrossRef]
  41. Garson, J.A.; Ferns, R.B.; Grant, P.R.; Ijaz, S.; Nastouli, E.; Szypulska, R.; Tedder, R.S. Minor groove binder modification of widely used TaqMan probe for hepatitis E virus reduces risk of false negative real-time PCR results. J. Virol. Methods 2012, 186, 157–160. [Google Scholar] [CrossRef]
  42. 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]
  43. Candido, A.; Taffon, S.; Chionne, P.; Pisani, G.; Madonna, E.; Dettori, S.; Hamza, A.; Valdarchi, C.; Bruni, R.; Ciccaglione, A.R. Diagnosis of HEV infection by serological and real-time PCR assays: A study on acute non-A-C hepatitis collected from 2004 to 2010 in Italy. BMC Res. Notes 2012, 5, 297. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Marcantonio, C.; Pezzotti, P.; Bruni, R.; Taliani, G.; Chionne, P.; Madonna, E.; Villano, U.; Pisani, G.; Equestre, M.; Dell’Orso, L.; et al. Incidence of hepatitis E virus infection among blood donors in a high endemic area of Central Italy. J. Viral Hepat. 2019, 26, 506–512. [Google Scholar] [CrossRef] [PubMed]
  45. Monini, M.; Di Bartolo, I.; Ianiro, G.; Angeloni, G.; Magistrali, C.F.; Ostanello, F.; Ruggeri, F.M. Detection and molecular characterization of zoonotic viruses in swine fecal samples in Italian pig herds. Arch. Virol. 2015, 160, 2547–2556. [Google Scholar] [CrossRef] [PubMed]
  46. 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, 133. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Di Bartolo, I.; Ponterio, E.; Castellini, L.; Ostanello, F.; Ruggeri, F.M. Viral and antibody HEV prevalence in swine at slaughterhouse in Italy. Vet. Microbiol. 2011, 149, 330–338. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Martelli, F.; Toma, S.; Di Bartolo, I.; Caprioli, A.; Ruggeri, F.M.; Lelli, D.; Bonci, M.; Ostanello, F. Detection of Hepatitis E Virus (HEV) in Italian pigs displaying different pathological lesions. Res. Vet. Sci. 2010, 88, 492–496. [Google Scholar] [CrossRef]
  49. Di Bartolo, I.; Martelli, F.; Inglese, N.; Pourshaban, M.; Caprioli, A.; Ostanello, F.; Ruggeri, F.M. Widespread diffusion of genotype 3 hepatitis E virus among farming swine in Northern Italy. Vet. Microbiol. 2008, 132, 47–55. [Google Scholar] [CrossRef] [Green Version]
  50. La Rosa, G.; Della Libera, S.; Brambilla, M.; Bisaglia, C.; Pisani, G.; Ciccaglione, A.R.; Bruni, R.; Taffon, S.; Equestre, M.; Iaconelli, M. Hepatitis E Virus (Genotype 3) in Slurry Samples from Swine Farming Activities in Italy. Food Environ. Virol. 2016, 9, 219–229. [Google Scholar] [CrossRef]
  51. Martelli, F.; Caprioli, A.; Zengarini, M.; Marata, A.; Fiegna, C.; Di Bartolo, I.; Ruggeri, F.M.; Delogu, M.; Ostanello, F. Detection of hepatitis E virus (HEV) in a demographic managed wild boar (Sus scrofa scrofa) population in Italy. Vet. Microbiol. 2008, 126, 74–81. [Google Scholar] [CrossRef]
  52. Serracca, L.; Battistini, R.; Rossini, I.; Mignone, W.; Peletto, S.; Boin, C.; Pistone, G.; Ercolini, R.; Ercolini, C. Molecular Investigation on the Presence of Hepatitis E Virus (HEV) in Wild Game in North-Western Italy. Food Environ. Virol. 2015, 7, 206–212. [Google Scholar] [CrossRef] [PubMed]
  53. De Sabato, L.; Ostanello, F.; De Grossi, L.; Marcario, A.; Franzetti, B.; Monini, M.; Di Bartolo, I. Molecular survey of HEV infection in wild boar population in Italy. Transbound. Emerg. Dis. 2018, 65, 1749–1756. [Google Scholar] [CrossRef] [PubMed]
  54. Sarchese, V.; Di, P.F.; Melegari, I.; Palombieri, A.; Sanchez, S.B.; Arbuatti, A.; Ciuffetelli, M.; Marsilio, F.; Martella, V.; Di, M.B. Hepatitis E virus in sheep in Italy. Transbound. Emerg. Dis. 2019, 66, 1120–1125. [Google Scholar] [CrossRef]
  55. Clemente-Casares, P.; Rodriguez-Manzano, J.; Girones, R. Hepatitis E virus genotype 3 and sporadically also genotype 1 circulate in the population of Catalonia, Spain. J. Water Health 2009, 7, 664–673. [Google Scholar] [CrossRef] [PubMed]
  56. Kamel, A.H.; Ali, M.A.; El-Nady, H.G.; Deraz, A.; Aho, S.; Pothier, P.; Belliot, G. Presence of enteric hepatitis viruses in the sewage and population of Greater Cairo. Clin. Microbiol. Infect. 2011, 17, 1182–1185. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Miura, T.; Lhomme, S.; Le Saux, J.C.; Le, M.P.; Guillois, Y.; Couturier, E.; Izopet, J.; Abranavel, F.; Le Guyader, F.S. Detection of Hepatitis E Virus in Sewage After an Outbreak on a French Island. Food Environ. Virol. 2016, 8, 194–199. [Google Scholar] [CrossRef] [Green Version]
  58. Ahmad, T.; Anjum, S.; Sadaf Zaidi, N.U.; Ali, A.; Waqas, M.; Afzal, M.S.; Arshad, N. Frequency of hepatitis E and Hepatitis A virus in water sample collected from Faisalabad, Pakistan. Ann. Agric. Environ. Med. 2015, 22, 661–664. [Google Scholar] [CrossRef]
  59. Beji-Hamza, A.; Hassine-Zaafrane, M.; Khelifi-Gharbi, H.; Della Libera, S.; Iaconelli, M.; Muscillo, M.; Petricca, S.; Ciccaglione, A.R.; Bruni, R.; Taffon, S.; et al. Hepatitis E virus genotypes 1 and 3 in wastewater samples in Tunisia. Arch. Virol. 2015, 160, 183–189. [Google Scholar] [CrossRef]
  60. Martinez Wassaf, M.G.; Pisano, M.B.; Barril, P.A.; Elbarcha, O.C.; Pinto, M.A.; Mendes de, O.J.; DiGiusto, P.; Nates, S.V.; Re, V.E. First detection of hepatitis E virus in Central Argentina: Environmental and serological survey. J. Clin. Virol. 2014, 61, 334–339. [Google Scholar] [CrossRef]
  61. Idolo, A.; Serio, F.; Lugoli, F.; Grassi, T.; Bagordo, F.; Guido, M.; Privitera, G.; Lobreglio, G.; De Donno, A. Identification of HEV in symptom-free migrants and environmental samples in Italy. J. Viral Hepat. 2013, 20, 438–443. [Google Scholar] [CrossRef]
  62. Ishida, S.; Yoshizumi, S.; Ikeda, T.; Miyoshi, M.; Goto, A.; Matsubayashi, K.; Ikeda, H. Detection and molecular characterization of hepatitis E virus in clinical, environmental and putative animal sources. Arch. Virol. 2012, 157, 2363–2368. [Google Scholar] [CrossRef] [PubMed]
  63. Iaconelli, M.; Muscillo, M.; Della Libera, S.; Fratini, M.; Meucci, L.; De Ceglia, M.; Giacosa, D.; La Rosa, G. One-year Surveillance of Human Enteric Viruses in Raw and Treated Wastewaters, Downstream River Waters, and Drinking Waters. Food Environ. Virol. 2017, 9, 79–88. [Google Scholar] [CrossRef] [PubMed]
  64. Ram, D.; Manor, Y.; Gozlan, Y.; Schwartz, E.; Ben-Ari, Z.; Mendelson, E.; Mor, O. Hepatitis E Virus Genotype 3 in Sewage and Genotype 1 in Acute Hepatitis Cases, Israel. Am. J. Trop. Med. Hyg. 2016, 95, 216–220. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Nieuwenhuijse, D.F.; Koopmans, M.P. Metagenomic Sequencing for Surveillance of Food- and Waterborne Viral Diseases. Front. Microbiol. 2017, 8, 230. [Google Scholar] [CrossRef] [Green Version]
  66. Bradanini, L.; Youkee, D.; Fabris, P.; Romano, L.; Brunetti, E.; Giordani, M.T. Acute hepatitis E virus infection in a migrant population in North East Italy: A retrospective analysis. Travel Med. Infect. Dis. 2017, 20, 37–42. [Google Scholar] [CrossRef]
  67. Fenaux, H.; Chassaing, M.; Berger, S.; Gantzer, C.; Bertrand, I.; Schvoerer, E. Transmission of hepatitis E virus by water: An issue still pending in industrialized countries. Water Res. 2019, 151, 144–157. [Google Scholar] [CrossRef]
  68. Puttini, C.; Riccio, M.L.; Redi, D.; Tordini, G.; Cenerini, M.; Romanello, F.; De, L.A.; Carmellini, M.; Fossombroni, V.; Cusi, M.G.; et al. Seroprevalence of hepatitis E virus (HEV) infection in blood donors and renal transplant recipients: A retrospective study from central Italy. Infez. Med. 2015, 23, 253–256. [Google Scholar]
  69. Zuin, M.; Caserta, C.; Romano, L.; Mele, A.; Zanetti, A.; Cannatelli, R.; Giorgini, A.; Tagliacarne, C.; Amante, A.; Marcucci, F.; et al. Seroepidemiology of HEV and HAV in two populations with different socio-economic levels and hygienic/sanitary conditions. Eur. J. Clin. Microbiol. Infect. Dis 2017, 36, 479–485. [Google Scholar] [CrossRef]
  70. Zanetti, A.R.; Dawson, G.J. Hepatitis type E in Italy: A seroepidemiological survey. Study Group of Hepatitis E. J. Med. Virol. 1994, 42, 318–320. [Google Scholar] [CrossRef]
  71. Coppola, R.C.; Masia, G.; Romano, L.; Tanzi, E.; Zanetti, A.R. Epidemiology and prevention of enterically transmitted hepatitis in Italy. Res. Virol. 1998, 149, 271–276. [Google Scholar] [CrossRef]
  72. Caruso, C.; Peletto, S.; Rosamilia, A.; Modesto, P.; Chiavacci, L.; Sona, B.; Balsamelli, F.; Ghisetti, V.; Acutis, P.L.; Pezzoni, G.; et al. Hepatitis E Virus: A Cross-Sectional Serological and Virological Study in Pigs and Humans at Zoonotic Risk within a High-Density Pig Farming Area. Transbound. Emerg. Dis. 2017, 64, 1443–1453. [Google Scholar] [CrossRef] [PubMed]
  73. Scotto, G.; Martinelli, D.; Centra, M.; Querques, M.; Vittorio, F.; Delli, C.P.; Tartaglia, A.; Campanale, F.; Bulla, F.; Prato, R.; et al. Epidemiological and clinical features of HEV infection: A survey in the district of Foggia (Apulia, Southern Italy). Epidemiol. Infect. 2014, 142, 287–294. [Google Scholar] [CrossRef] [PubMed]
  74. Spada, E.; Pupella, S.; Pisani, G.; Bruni, R.; Chionne, P.; Madonna, E.; Villano, U.; Simeoni, M.; Fabi, S.; Marano, G.; et al. A nationwide retrospective study on prevalence of hepatitis E virus infection in Italian blood donors. Blood Transfus 2018, 16, 413–421. [Google Scholar] [CrossRef] [PubMed]
  75. Lucarelli, C.; Spada, E.; Taliani, G.; Chionne, P.; Madonna, E.; Marcantonio, C.; Pezzotti, P.; Bruni, R.; La, R.G.; Pisani, G.; et al. High prevalence of anti-hepatitis E virus antibodies among blood donors in central Italy, February to March 2014. Eurosurveillance 2016, 21, 30299. [Google Scholar] [CrossRef] [PubMed]
  76. Di Bartolo, I.; Angeloni, G.; Ponterio, E.; Ostanello, F.; Ruggeri, F.M. Detection of hepatitis E virus in pork liver sausages. Int. J. Food Microbiol. 2015, 193, 29–33. [Google Scholar] [CrossRef] [PubMed]
  77. Alfonsi, V.; Romano, L.; Ciccaglione, A.R.; La Rosa, G.; Bruni, R.; Zanetti, A.; Della Libera, S.; Iaconelli, M.; Bagnarelli, P.; Capobianchi, M.R.; et al. Hepatitis E in Italy: 5 years of national epidemiological, virological and environmental surveillance, 2012 to 2016. Eurosurveillance 2018, 23, 1700517. [Google Scholar] [CrossRef]
  78. Iaconelli, M.; Purpari, G.; Della Libera, S.; Petricca, S.; Guercio, A.; Ciccaglione, A.R.; Bruni, R.; Taffon, S.; Equestre, M.; Fratini, M.; et al. Hepatitis A and E Viruses in Wastewaters, in River Waters, and in Bivalve Molluscs in Italy. Food Environ. Virol. 2015, 7, 316–324. [Google Scholar] [CrossRef]
  79. La Rosa, G.; Proroga, Y.T.R.; De, M.D.; Capuano, F.; Iaconelli, M.; Della, L.S.; Suffredini, E. First Detection of Hepatitis E Virus in Shellfish and in Seawater from Production Areas in Southern Italy. Food Environ. Virol. 2018, 10, 127–131. [Google Scholar] [CrossRef]
  80. Purpari, G.; Macaluso, G.; Di, B.S.; Gucciardi, F.; Mira, F.; Di, M.P.; Lastra, A.; Petersen, E.; La, R.G.; Guercio, A. Molecular characterization of human enteric viruses in food, water samples, and surface swabs in Sicily. Int. J. Infect. Dis. 2019, 80, 66–72. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Geographic information system (GIS) map of the WTPs included in the study.
Figure 1. Geographic information system (GIS) map of the WTPs included in the study.
Ijerph 17 02059 g001
Figure 2. The distribution of positive samples (blue columns) vs total samples (orange columns) over the sampling period (2011–2019). The percentage of positive samples for each year is indicated above the absolute number.
Figure 2. The distribution of positive samples (blue columns) vs total samples (orange columns) over the sampling period (2011–2019). The percentage of positive samples for each year is indicated above the absolute number.
Ijerph 17 02059 g002
Figure 3. The number of HEV genotype 1 (blue columns) and genotype 3 (orange columns) positive samples detected during the study period (2011–2019).
Figure 3. The number of HEV genotype 1 (blue columns) and genotype 3 (orange columns) positive samples detected during the study period (2011–2019).
Ijerph 17 02059 g003
Figure 4. (top) Neighbor-joining phylogenetic tree of HEV-G3 sequences. (bottom) Neighbor-joining phylogenetic tree of HEV-G1 sequences.
Figure 4. (top) Neighbor-joining phylogenetic tree of HEV-G3 sequences. (bottom) Neighbor-joining phylogenetic tree of HEV-G1 sequences.
Ijerph 17 02059 g004aIjerph 17 02059 g004b
Table 1. Nested RT-PCR results according to Region and geographic area.
Table 1. Nested RT-PCR results according to Region and geographic area.
AreaRegionNumber of Tested SamplesNumber of Positive Samples (%)
NorthEmilia Romagna352 (5.7)
Friuli Venezia Giulia111 (9.1)
Liguria250 (0)
Lombardy1628 (4.9)
Piedmont825 (6.1)
Trentino-South Tyrol1445 (3.5)
Veneto662 (3.0)
Aosta Valley40 (0)
subtotal 52923 (4.3)
CenterTuscany240 (0)
Umbria50 (0)
Latium27020 (7.4)
Marche141 (7.1)
Abruzzo a27721 (7.2)
subtotal 59042 (7.1)
SouthApulia402 (5.0)
Basilicata120 (0)
Calabria251 (4.0)
Campania890 (0)
Molise355 (14.2)
Sardinia240 (0)
Sicily301 (3.3)
subtotal 2559 (3.5)
Total 137474 (5.4)
a Abruzzo is officially a region of Southern Italy, because of its historic association with the Kingdom of the Two Sicilies, but is geographically located in Central Italy, bordering Lazio and Marche administrative regions.
Table 2. Genotyping and HEV concentration in positive samples.
Table 2. Genotyping and HEV concentration in positive samples.
Sample IDCitySampling Date (day/month/year)RegionAreaGenotype Concentration (g.c./L)
1480Rome18/05/2011LatiumCenterG37.9 × 102
1484Rome 18/05/2011LatiumCenterG3nt
1498Rome20/06/2011LatiumCenterG39.2 × 102
1536Bari01/08/2011ApuliaSouthG1<LOD
1542Palermo08/03/2011SicilySouthG1<LOD
1567Venice27/09/2011VenetoNorthG1nt
1577Milan10/09/2011LombardyNorthG3nt
1581Milan21/09/2011LombardyNorthG1<LOD
1588Milan16/11/2011LombardyNorthG37.3 × 102
1648Ancona18/11/2011MarcheCenterG3<LOD
1708Bari06/12/2011ApuliaSouthG35.9 × 102
1729Brescia01/06/2012LombardyNorthG1<LOD
1733Rome31/05/2012LatiumCenterG19.9 × 102
1734 Bologna22/11/2011Emilia RomagnaNorthG12.3 × 103
1739Venice 23/04/2012VenetoNorthG13.7 × 101
1745Trento 07/02/2012Trentino Alto Adige NorthG31.3 × 102
1746Trento 07/03/2012Trentino Alto Adige NorthG1<LOD
1747Trento 04/04/2012Trentino Alto Adige NorthG12.0 × 102
1750Trento 18/04/2012Trentino Alto Adige NorthG3<LOD
1762Turin24/05/2012PiedmontNorthG1<LOD
1763Turin20/12/2011PiedmontNorthG12.0 × 101
1785Bolzano01/05/2012Trentino Alto Adige NorthG3<LOD
1792Reggio Calabria08/04/2012CalabriaSouthG1nt
1927Milan14/03/2013LombardyNorthG11.2 × 103
1929Milan26/06/2013LombardyNorthG35.5 × 101
2020Collegno06/03/2013PiedmontNorthG33.1 × 101
2124Trieste02/10/2013Friuli Venezia GiuliaNorthG3nt
2151Pescara03/04/2012AbruzzoCenterG31.9 × 102
2208Milan28/08/2014LombardyNorthG32.9 × 101
2212Campobasso18/03/2014MoliseSouthG3<LOD
2214Campobasso10/02/2014MoliseSouthG3<LOD
2215Campobasso06/08/2014MoliseSouthG31.4 × 103
2219Rome16/07/2014LatiumCenterG17.8 × 101
2220Rome26/07/2014LatiumCenterG3<LOD
2222Rome04/08/2014LatiumCenterG1nt
2226Campobasso20/10/2014MoliseSouthG3nt
2267Campobasso30/03/2015MoliseSouthG1<LOD
2275Rome10/07/2015LatiumCenterG3<LOD
2320San Salvo27/01/2016AbruzzoCenterG19.0 × 101
2345Vasto03/02/2016AbruzzoCenterG37.5 × 102
2353Pescara15/03/2016AbruzzoCenterG33.4 × 102
2359Montesilvano15/03/2016AbruzzoCenterG31.0 × 102
2390Rome09/05/2016LatiumCenterG34.4 × 102
2393Rome30/05/2016LatiumCenterG3<LOD
2399Rome25/07/2016LatiumCenterG32.7 × 101
2415Pescara11/07/2016AbruzzoCenterG33.8 × 102
2420Pescara26/09/2016AbruzzoCenterG32.8 × 102
2428Chieti08/07/2016Abruzzo CenterG3<LOD
2431Chieti26/08/2016Abruzzo CenterG3<LOD
2433Chieti29/09/2016Abruzzo CenterG31.3 × 102
2459Pescara12/12/2016AbruzzoCenterG3<LOD
2487Sulmona26/06/2017AbruzzoCenterG3nt
2492Vasto19/04/2017AbruzzoCenterG32.5 × 102
2501Rome 18/07/2017LatiumCenterG38.0 × 101
2547Rome22/12/2017LatiumCenterG32.0 × 102
2625Bologna08/08/2017Emilia RomagnaNorthG36.3 × 101
2592Rome25/01/2018LatiumCenterG3<LOD
2596Rome22/02/2018LatiumCenterG3nt
2602Rome21/03/2018LatiumCenterG3<LOD
2646Vasto14/03/2018AbruzzoCenterG34.3 × 101
2696Rome27/06/2018LatiumCenterG32.0 × 102
2729Turin26/04/2018PiedmontNorthG36.4 × 102
2733Turin13/06/2018PiedmontNorthG3<LOD
2805Rome17/04/2018LatiumCenterG33.2 × 103
2812Rome01/12/2018LatiumCenterG3<LOD
2842Avezzano04/07/2018AbruzzoCenterG3<LOD
2846San Salvo09/07/2018AbruzzoCenterG35.3 × 103
2852Vasto12/07/2018AbruzzoCenterG3<LOD
2891Milan10/07/2018LombardyNorthG32.8 × 103
3051Vasto 17/04/2019AbruzzoCenterG36.9 × 103
3053Vasto 17/06/2019AbruzzoCenterG32.8 × 104
3093Pescara15/05/2019AbruzzoCenterG39.5 × 103
3154Rome01/05/2019LatiumCenterG35.2 × 102
3161San Salvo31/10/2019AbruzzoCenterG3<LOD
nt = not tested due to insufficient RNA. Quantitative values in bold were above the analytical LOQ (1.1 × 103 g.c./L); other values are provided as estimated concentrations.

Share and Cite

MDPI and ACS Style

Iaconelli, M.; Bonanno Ferraro, G.; Mancini, P.; Suffredini, E.; Veneri, C.; Ciccaglione, A.R.; Bruni, R.; Della Libera, S.; Bignami, F.; Brambilla, M.; et al. Nine-Year Nationwide Environmental Surveillance of Hepatitis E Virus in Urban Wastewaters in Italy (2011–2019). Int. J. Environ. Res. Public Health 2020, 17, 2059. https://doi.org/10.3390/ijerph17062059

AMA Style

Iaconelli M, Bonanno Ferraro G, Mancini P, Suffredini E, Veneri C, Ciccaglione AR, Bruni R, Della Libera S, Bignami F, Brambilla M, et al. Nine-Year Nationwide Environmental Surveillance of Hepatitis E Virus in Urban Wastewaters in Italy (2011–2019). International Journal of Environmental Research and Public Health. 2020; 17(6):2059. https://doi.org/10.3390/ijerph17062059

Chicago/Turabian Style

Iaconelli, Marcello, Giusy Bonanno Ferraro, Pamela Mancini, Elisabetta Suffredini, Carolina Veneri, Anna Rita Ciccaglione, Roberto Bruni, Simonetta Della Libera, Francesco Bignami, Massimo Brambilla, and et al. 2020. "Nine-Year Nationwide Environmental Surveillance of Hepatitis E Virus in Urban Wastewaters in Italy (2011–2019)" International Journal of Environmental Research and Public Health 17, no. 6: 2059. https://doi.org/10.3390/ijerph17062059

APA Style

Iaconelli, M., Bonanno Ferraro, G., Mancini, P., Suffredini, E., Veneri, C., Ciccaglione, A. R., Bruni, R., Della Libera, S., Bignami, F., Brambilla, M., De Medici, D., Brandtner, D., Schembri, P., D’Amato, S., & La Rosa, G. (2020). Nine-Year Nationwide Environmental Surveillance of Hepatitis E Virus in Urban Wastewaters in Italy (2011–2019). International Journal of Environmental Research and Public Health, 17(6), 2059. https://doi.org/10.3390/ijerph17062059

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop