Next Article in Journal
Different Infectivity and Transmissibility of H5N8 and H5N1 High Pathogenicity Avian Influenza Viruses Isolated from Chickens in Japan in the 2021/2022 Season
Next Article in Special Issue
Simple Wastewater Preparation Protocol Applied to Monitor the Emergence of the Omicron 21L/BA.2 Variant by Genome Sequencing
Previous Article in Journal
Comparison of a Genotype 1 and a Genotype 2 Macaque Foamy Virus env Gene Indicates Distinct Infectivity and Cell-Cell Fusion but Similar Tropism and Restriction of Cell Entry by Interferon-Induced Transmembrane Proteins
Previous Article in Special Issue
New RT-PCR Assay for the Detection of Current and Future SARS-CoV-2 Variants
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Expanding the Pathogen Panel in Wastewater Epidemiology to Influenza and Norovirus

by
Rudolf Markt
1,2,*,
Fabian Stillebacher
1,
Fabiana Nägele
1,
Anna Kammerer
1,
Nico Peer
1,
Maria Payr
1,
Christoph Scheffknecht
3,
Silvina Dria
3,
Simon Draxl-Weiskopf
1,
Markus Mayr
1,
Wolfgang Rauch
4,
Norbert Kreuzinger
5,
Lukas Rainer
6,
Florian Bachner
6,
Martin Zuba
6,
Herwig Ostermann
6,
Nina Lackner
2,
Heribert Insam
1 and
Andreas Otto Wagner
1
1
Department of Microbiology, Universität Innsbruck, 6020 Innsbruck, Austria
2
Department of Health Sciences and Social Work, Carinthia University of Applied Sciences, 9020 Klagenfurt, Austria
3
Institut für Umwelt und Lebensmittelsicherheit des Landes Vorarlberg, 6900 Bregenz, Austria
4
Department of Infrastructure, Universität Innsbruck, 6020 Innsbruck, Austria
5
Institute for Water Quality and Resource Management, Technische Universität Wien, 1040 Vienna, Austria
6
Austrian National Public Health Institute, 1010 Vienna, Austria
*
Author to whom correspondence should be addressed.
Viruses 2023, 15(2), 263; https://doi.org/10.3390/v15020263
Submission received: 27 November 2022 / Revised: 1 January 2023 / Accepted: 12 January 2023 / Published: 17 January 2023
(This article belongs to the Special Issue Wastewater-Based Epidemiology (WBE) in COVID-19 Pandemics)

Abstract

:
Since the start of the 2019 pandemic, wastewater-based epidemiology (WBE) has proven to be a valuable tool for monitoring the prevalence of SARS-CoV-2. With methods and infrastructure being settled, it is time to expand the potential of this tool to a wider range of pathogens. We used over 500 archived RNA extracts from a WBE program for SARS-CoV-2 surveillance to monitor wastewater from 11 treatment plants for the presence of influenza and norovirus twice a week during the winter season of 2021/2022. Extracts were analyzed via digital PCR for influenza A, influenza B, norovirus GI, and norovirus GII. Resulting viral loads were normalized on the basis of NH4-N. Our results show a good applicability of ammonia-normalization to compare different wastewater treatment plants. Extracts originally prepared for SARS-CoV-2 surveillance contained sufficient genomic material to monitor influenza A, norovirus GI, and GII. Viral loads of influenza A and norovirus GII in wastewater correlated with numbers from infected inpatients. Further, SARS-CoV-2 related non-pharmaceutical interventions affected subsequent changes in viral loads of both pathogens. In conclusion, the expansion of existing WBE surveillance programs to include additional pathogens besides SARS-CoV-2 offers a valuable and cost-efficient possibility to gain public health information.

1. Introduction

Wastewater-based epidemiology (WBE) is a powerful tool for monitoring the spread of viral pathogens within a population connected to a sewage system, as shown during the COVID-19 pandemic [1,2,3,4,5]. Viruses such as poliovirus [6,7], hepatitis E virus [8], norovirus (NoV) [9,10], influenza virus [11,12], and many more have already been detected [13] and monitored in wastewater. In contrast to traditional individual-level testing, WBE offers several advantages for the cost-efficient early detection of local or regional outbreaks and long-term monitoring of pathogens at community level [14]. Hence, comparability among different wastewater treatment plants (WWTPs) is essential and can be achieved by using population size markers to normalize WBE data [15,16]. In this context, NH4-N concentration or NH4-N load is a valuable population size marker [16,17,18]. This allows an efficient pandemic or pathogen management, including the implementation and validation of non-pharmaceutical interventions (NPIs) [18]. Recently, the COVID-19 pandemic demonstrated the importance of WBE for pathogen surveillance and accelerated the world-wide implementation of this tool [19].
Although WBE is a cost-efficient way to track the spread of pathogens at community-level if calculated per single community member, the per-sample costs can be high. Sample logistics significantly contribute to the high costs, especially in rural regions [20]. Currently, wastewater samples are mainly prepared for the subsequent quantification of SARS-CoV-2. During this process, the genomic material of many viruses extracted as by-catch along with SARS-CoV-2 RNA is usually disregarded, although it could contain valuable information.
Influenza virus causes a seasonal epidemic (the flu season). The enveloped virus is usually spread through respiratory droplets and aerosols. Influenza A virus (IAV), influenza B virus (IBV), and influenza C virus infect humans, with IAV and IBV being responsible for most cases. IAV is normally responsible for epidemics and pandemics. IBV undergoes genetic and antigenic changes to a smaller extent and shows increased activity every 2–4 years [21]. Worldwide, influenza viruses result in up to five million cases of severe illness every year [22]. In the Northern Hemisphere, the typical flu season spans from October to May. The virus is transmitted mainly through the respiratory path, but its RNA may also be detected in the stool of symptomatic patients, with a detection rate ranging from 7.2% to 47% [23,24,25,26].
Noroviruses (NoVs) are a common cause of gastroenteritis and are usually spread through the fecal–oral route [27]. Several genogroups of the non-enveloped NoVs can be distinguished. The genogroups GI and GII are the most relevant to humans and are responsible for local and regional outbreaks. NoVs are a highly contagious group of viruses that have resulted in nearly 700 million annual cases worldwide [28], with infection commonly occurring during the winter months. The RNA of NoVs can be detected in the stool of patients with gastroenteritis with a detection rate ranging from 17.9% to 38.3% [29,30]. This RNA can still be detected several weeks after infection [31]. Another study found NoV in stool samples of 361 out of 2205 (16.4%) asymptomatic individuals [32].
Therefore, the aims of this retrospective study were (I) to detect and quantify two different viruses from wastewater extract as a by-catch of a WBE-based SARS-CoV-2 surveillance program, applying NH4-N as a population marker; (II) to compare the WBE data with incidence from inpatients; and (III) to examine the effect of COVID-19-related NPIs on the spread of the other pathogens. For this purpose, we searched for the enveloped, respiratory IAV and IBV known for their droplet-based transmission, and the non-enveloped, gastro-enteric NoV of GI and GII known for their fecal–oral transmission in already present wastewater extracts. The study includes samples from the SARS-CoV-2 surveillance program in Austria from the winter season of 2021/2022 [33].

2. Materials and Methods

2.1. Sample Collection and Sample Preparation

Approximately 250 mL of 24 h volume-proportionally-mixed influent was collected twice a week from 11 WWTPs in two federal states of Austria from September 2021 to April 2022 (Figure 1). Five WWTPs in Vorarlberg (V001–V005) covering ~62% of the population and six WWTPs in Salzburg (Z001, Z002, Z003, Z006, Z008, Z012) covering ~80% of the population were monitored. Raw samples were chilled to 4 °C [34] and processed within two days. Concentration via polyethylene glycol precipitation (PEG precipitation) and purification via the Monarch™ Total RNA Miniprep Kit (New England Biolabs, Ipswich, MA, USA) were conducted according to Markt et al. (2021) [34]. Residual RNA extracts after SARS-CoV-2 surveillance were stored at −80 °C until one-step reverse-transcriptase digital PCR analysis (RT-dPCR). Daily wastewater flow rate and ammonium–nitrogen concentration (NH4-N measured using standard commercial cuvette tests) were reported by the WWTPs.

2.2. Probe Based RT-Digital PCR-Detection of Influenza Virus and NoV

For the detection of IAV, primers and probes targeted the matrix proteins 1 and 2 in segment 7, and for IBV the nuclear export protein (NEP) and nonstructural protein 1 (NS1) in segment 8 of the genome, according the SARS-CoV-2/IAV/IBV multiplex assay published by the US Centers for Disease Control and Prevention [35]. Primers and probes targeting the ORF1–ORF2 junction were used for specific detection of NoV GI and GII following a duplex approach [36]. All primers and probes were obtained from Microsynth AG (Balgach, Switzerland). Quenchers and reporters were modified. Details on Primers and Probes can be found in Table S1. Positive control sequences (Table S2) for IAV, IBV, NoV GI and NoV GII were obtained from Eurofins Genomics AT (Vienna, Austria) as plasmids cloned into a standard plasmid vector (pEX-A128) (Table S2). In vitro amplification of the plasmid was performed using the NEBNext Ultra™ II Q5 Master Mix (New England Biolabs, Ipswich, MA, USA) using a standard vector primer pair (standard PEX primer) from Eurofins Genomics AT (Vienna, Austria). The thermocycling conditions for amplification of the positive controls were set according to the manufacturer’s protocol for the Q5 High-Fidelity DNA Polymerase. Agarose-gel purification of the resulting PCR product was performed using the Monarch DNA Gel Extraction Kit (New England Biolabs, Ipswich, MA, USA). Finally, the concentration of the in vitro amplified linear DNA was measured using a Quantus Fluorometer E6150 (Promega, Madison, WI, USA) and diluted to <1000 gcµL−1. The PCR products were used as positive controls and were integrated into the assay as linear dsDNA. To detect and quantify absolute RNA copy numbers, two duplex RT-digital PCR (RT-dPCR) assays were performed on a QIAcuity One 5-plex digital PCR system (QIAGEN GmbH, Hilden, Germany): one for influenza viruses and one for NoV detection. RT-dPCR of the 20 µL reaction mixture consisted of: 10 µL of Luna Universal Probe One-Step Reaction Mix (2x); 1 µL of Luna WarmStart RT Enzyme Mix (20x) (both New England Biolabs, Ipswich, MA, USA); 3.2 μL of mixed primers (forward and reverse, final concentration 0.8 μM for each target); 0.8 μL of mixed probes (final concentration 0.2 μM for each target); 0.5 µL of reference dye, and 4.5 μL of template. The PCR program was identical for both assays: An initial reverse transcription step at 55 °C for 10 min was followed by initial denaturation at 95 °C for 1 min. Subsequently, a total of 40 cycles at 95 °C for 10 s and 60 °C for 40 s were performed within 96-well Nanoplates with 8500 partitions per well (QIAGEN GmbH, Hilden, Germany). Fluorescence was read at the end of the 40 cycles. Gene copy numbers per positive partition were calculated according to the Poisson distribution. The concentration in the reaction mix was calculated based on the volume of positive and negative valid partitions. These steps were done using the internal software Qiacuity Software Suite (QIAGEN GmbH, Hilden, Germany). Subsequently, the results were compensated for the volume of the template in the dPCR reaction, the RNA extract eluted after purification, and the initial wastewater sample volume. The limit of detection for InfA and InfB was 1 gc per reaction, and for NoV GI and GII 3 gc per reaction.

2.3. Inpatients Infected with NoV or Influenza Virus and Anti-COVID-19 NPIs

Inpatients with primary or secondary infection with influenza virus (J09, J090, J091, J098, J10, J100, J101, J101, and J108 according to ICD-10-GM-2022) and norovirus (A081 according to ICD-10-GM-2022) were aggregated according to admission week and federal state. Six events increasing or decreasing the impact of NPIs on the spread of the COVID-19 pandemic were set from September 2021 to March 2022: (I) September 9—school opens after the summer holidays, (II) November 15—lockdown for the unvaccinated, (III) November 22nd—complete lockdown, (IV) December 12—partial opening, (V) January 11—Omicron measures, and (VI) March 5—complete opening.

2.4. Normalization of Wastewater Signal and Statistical Analysis

Viral load was defined as the concentration of gene copies detected multiplied by the daily wastewater flow rate. For comparison between WWTPs, viral loads were normalized to NH4-N load assuming a NH4-N excretion per person and day of 8 g N [37,38] resulting in mega (M), giga (G), or tera (T) gene copies per person per day (gcPE−1/d−1). To test for significant difference of reported NH4-N concentration between WWTPs, the Kruskal–Wallis rank sum test and the Wilcoxon rank sum as post hoc test with Bonferroni correction were performed.
Wastewater data were down-sampled from biweekly datapoints to weekly datapoints (mean loads of both datapoints), as inpatient numbers were resolved at a weekly basis [39]. To compare wastewater signal and inpatient number, a cross-correlation analysis and ranked correlation analysis (Spearman rank correlation test) between aggregated viral loads (sum of viral loads) per calendar week across all WWTPs and inpatient numbers were performed. The significance cut-off was set at α = 0.05 and at α = 0.005 for highly significant for all analyses. Smoothing was conducted using a LOWESS-smoothing with a span of 0.5.
Calculations, statistical analysis, and figures were prepared in RStudio© (2022.07.1 + 554) with R version 4.1.1. Data organization, normalization, and statistical analysis were carried out using R-base [40]. Further data organization was carried out using the packages tidyquant [41] and readxl [42]. Figures were prepared using the packages ggplot2 [43] and cowplot [44].

3. Results

Approximately 70% of the total population of the two Federal States Salzburg and Vorarlberg, or 0.7 Mio. inhabitants, are connected to the sewage system of the 11 WWTPs monitored for this study. The number of connected inhabitants per WWTP ranged from 321,000 (Z001) to 14,000 (Z012) as shown in Table 1. Compared with registered numbers, the population size marker NH4-N slightly overestimated the number of inhabitants with 113% ± 25% (mean ± SD) across all WWTPs. A conspicuous discrepancy was found for the WWTPs V003 (146%) and Z002 (170%). Moreover, highly significant differences were detected in the NH4-N concentration between WWTPs. The post hoc test showed that WWTPs Z002 and Z008 differed significantly from other WWTPs. Wastewater flow rate, NH4-N concentration, and the calculated population for each of the WWTPs are depicted in Figure S1.
In total we analyzed 538 RNA extracts. IAV, NoV GI, and NoV GII were detected at least once in each of the WWTPs, while IBV was not detected at all between September 2021 to March 2022. The normalized IAV load (Figure 2A) was three orders of magnitude lower than the normalized NoV GII load (Figure 2B). IAV was under the limit of detection for many WWTPs during autumn, but started to rise in the middle of December. The IAV signal still rose at the end of the measurement period in March 2022. On the contrary, NoV load in wastewater was highest during autumn. Comparing the two NoV groups, NoV GII was ten to a hundred times more abundant than NoV GI. Moreover, NoV GI and NoV GII did not show the same temporal dynamics. NoV GI was highest in the middle of September (calendar week 38–39, 2021), while NoV GII peaked four to five weeks later (calendar week 43–44, 2021).
Inpatient numbers for both viral infections were low during the measurement period compared with previous years (Figure 3). Additionally, it was not possible to distinguish between virus variants or genogroups for influenza or NoV, respectively. Although resolution of the data allowed WWTP-specific inpatient numbers, the granularity was very high at this level. Furthermore, patients from rural regions may change catchment area when hospitalized. Therefore, inpatient numbers were aggregated at a level which includes both federal states, together covering approx. 0.95 Mio. inhabitants. The number of hospitalized persons infected with NoV peaked in autumn 2021 (Figure 4B). Distinct from that dynamic, the number of hospitalized persons infected with the influenza virus started to rise during the Christmas Holidays of 2021 and continued to rise at the end of the measurement period in March 2022 (Figure 4A).
Cross-correlation analysis showed a significant correlation between the aggregated virus load from all WWTPs and the aggregated inpatient numbers using IAV viral load (Figure S2) and NoV GII viral load (Figure S3) at a weekly basis. Neither IAV nor NoV GII showed a positive or negative lag in their temporal dynamics in comparison with the timeseries from respective inpatient numbers integrating the complete analysis period. Hence, the rank correlation test was also highly significant for IAV and significant for NoV GII. No significant correlation was obtained by including NoV GI load into analysis.
NPIs due to the COVID-19 pandemic affected the spread of influenza and NoV in Austria. The normally occurring seasonal peak during the winter period was missed in 2020. Additionally, the onset for the winter season 2021/2022 was earlier and the inpatient number was lower for influenza infections in comparison with seasons before the COVID-19 pandemic (Figure 3). Effects of NPIs on the spread of both of the investigated pathogens were complex. Nevertheless, it could be shown that a change in NPIs strategy coincided with a change of the virus signal in wastewater and the number of inpatients (Figure 4). Hence, in temporal proximity to the lockdowns (calendar week 46–47, 2021), virus loads for IAV and NoV GII decreased in November 2021. After Omicron-related measures (calendar week 02, 2022), the signal of NoV from wastewater dropped again in January 2021. The signal drop of IAV after Omicron-related measures was temporally shifted, highlighting the complexity of this topic. In contrast to restrictions, relaxation of regulations came along with an increase in virus loads as can be seen for the partial opening (calendar week 50, 2021) in December 2021 and the complete opening (calendar week 09, 2022) in March 2022.

4. Discussion

Archived RNA extracts from SARS-CoV-2 surveillance could easily be used to retrospectively track the seasonal occurrence of infections with influenza and norovirus. To normalize pathogen concentration, we assumed a daily NH4-N load per person of 8 gN [37] and were able to emulate the number of inhabitants attached to the sewer system by 113% ± 25%. The conspicuous result for WWTP Z002 (overestimation of population by 170%) may be explained by a significantly high NH4-N concentration, causing an unexpectedly high daily NH4-N load. Investigations for WWTP Z002 revealed that increased NH4-N concentrations came from internal returns of process water (anaerobic sludge desiccation). NH4-N concentration of the second conspicuous WWTP, V003 (overestimation of population by 146%), was not significantly different from other WWTPs assuming normal, urban wastewater. However, investigations for WWTP V003 revealed that industrial discharges (from an industrial bakery and industrial meat processors) as well as internal returns of process water (anaerobic sludge desiccation) increased the NH4-N load. Consequently, we recommend normalization of viral loads for inter-WWTP comparability, but the applicability of the population-size marker NH4-N must be verified for each WWTP.
Non-pharmaceutical interventions during the COVID-19 pandemic affected the spread of influenza viruses [45,46] and NoV [47,48]. In comparison with the influenza and NoV epidemics of 2017–2019, the annual sharp increase in infected inpatients was missing for both pathogens in 2020 and 2021. In our study, IAV and NoV showed different prevalence, with concentrations of NoV being approximately 1000 times higher. IBV activity was low during the winter season 2021/2022 in the northern hemisphere [49] and was not found in wastewater samples. Although, IAV activity was low during the measurement period due to anti-COVID measures, quantifiable traces of genomic material were found in wastewater via RT-dPCR. Aside from low prevalence, several other factors could explain the low IAV concentration in wastewater. These include the droplet-based transmission route via the respiratory tract, the relatively low fecal shedding [24,25,26], and the lower stability of enveloped compared to non-enveloped viruses in wastewater [13]. Additionally, recent studies reported that IAV is predominantly found in the solid fraction of wastewater [11,12]. This fraction was neglected during sample preparation for this study, as the main aim of RNA extraction was the detection of SARS-CoV-2 genomic material. According to our experience, the purification of genomic material of SARS-CoV-2 from wastewater is more easily achieved by processing the liquid fraction. Nevertheless, the signal from the liquid fraction would be sufficient for a surveillance program regarding the spread of IAV. In contrast to IAV, the non-enveloped NoV causes waterborne outbreaks, is shed via feces in large quantities, and is relatively stable in wastewater [50], and thus is suitable for WBE surveillance from already existing extracts.
Correlation between viral loads in wastewater and inpatients were significant for both, IAV and NoV GII, which is in accordance with previous studies [9,10,11,12]. Inclusion of the GI wastewater signal reduced the correlation coefficient which indicates a difference in hospitalization rates between GII and GI. Clinical information on the matter, however, is sparsely available. Desai et al. (2012) [51] report in their review including 71,724 illnesses and 501 hospitalizations a significantly higher hospitalization rate in GII.4 cases than in non-GII.4 cases, showing that different genotypes can generally lead to different case severity. Cross-correlation analysis indicated that the wastewater signal coincided with inpatient numbers for IAV und NoV GII. Huang et al. (2022) [10] reported similar results for NoV GII comparing wastewater data with cases and outbreaks of acute gastroenteritis. Kazama et al. (2017) [9] reported similar results for NoV in comparison with clinical cases.
In contrast to our findings, Mercier et al. (2022) [11] reported a wastewater lead time for influenza of more than two weeks over inpatient testing and testing of institutional outbreaks. Discrepancies between our study and other findings regarding lead time may have several causes. One important aspect is that in our study infected inpatients were assigned to the admission week and not to the date of onset of symptoms or diagnosis. An additional limitation is the temporal resolution at a weekly basis which may be too coarse for lead/lag-time calculations for periods smaller than seven days. Concerning NoV, we observed further an interesting temporal shift of four to five weeks between NoV GI and GII. NoV activity normally increases during the cold months, with lower activity in summer [52]. However, NoV GI concentration was highest in September in the present study. Nordgren et al. (2009) reported a high concentration of NoV GI in wastewater during the summer in contrast to NoV GII [53]. Thus, NoV surveillance should also consider the warmer seasons.
To prove a causal relationship between single NPIs and pathogen spread, additional variables including measures-fatigue, NPI-awareness due to medial impact, weather conditions, and mobility would be necessary. Nevertheless, we were able to assign the onset of restrictive measures to a subsequent decrease and the relaxation of measures to a subsequent increase in wastewater signal for both IAV and NoV GII.

5. Conclusions

It could be shown that an existing WBE program for SARS-CoV-2 surveillance can easily be expanded to monitor influenza and norovirus. Concentrated and purified RNA extracts from the liquid fraction of wastewater for SARS-CoV-2 detection contains sufficient genomic material for the surveillance of both influenza virus and norovirus. The resulting data contain valuable public health information on the spread of viral pathogens. Further, normalization for population size is independent from the number of pathogens tracked. Thus, the effort for normalization remains constant, while the benefit increases with every additional pathogen detected. The molecular biological detection of additional pathogens as by-catch in existing WBE infrastructure offers great potential for pathogen surveillance.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/v15020263/s1, Table S1 Primers and Probes for IAV, IBV, NoV GI and NoV GII detection via RT-dPCR; Table S2 Linearized DNA positive control sequences after in vitro amplification; Figure S1 Flow rate, NH4-N concentration and calculated population size according to NH4-N normalization; Figure S2 Cross-correlation analysis for IVA and inpatients infected with influenza virus; Figure S3 Cross-correlation analysis for NoV GII and inpatients infected with NoVs.

Author Contributions

Conceptualization, R.M., H.I. and A.O.W.; Data curation, M.M., L.R., F.B., M.Z. and H.O.; Formal analysis, R.M. and N.L.; Funding acquisition, H.I. and A.O.W.; Investigation, F.S., F.N., A.K., N.P., M.P. and S.D.; Methodology, R.M. and W.R.; Project administration, A.O.W.; Resources, C.S., S.D.-W., W.R., N.K., L.R., F.B., M.Z. and H.O.; Software, M.M.; Supervision, R.M.; Validation, R.M.; Visualization, R.M. and N.L.; Writing—original draft, R.M., F.S. and N.L.; Writing—review & editing, R.M., F.S., N.L., H.I. and A.O.W. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the Austrian Federal Ministry of Education, Science and Research (project: Schulstandort-Monitoring). Publication was made possible by Publikationsfonds der Univerisität Innsbruck.

Institutional Review Board Statement

Ethical review and approval were waived for this study as they were not required. We used an existing secondary dataset maintained as database to collect the information needed. It includes anonymous information and it is not possible to link data in order to generate identifiable information. Data are properly anonymised. The information is legally accessible and appropriately protected by law.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data can be shared on request.

Acknowledgments

We would like to thank the wastewater treatment plants in Vorarlberg and Salzburg for their supply of wastewater samples and for providing the information on flowrate and NH4-N concentration. Furthermore, we are grateful to Arno Sorger and Michael Hirscher from W.H.U GmbH for the pre-processing of the raw samples from Salzburg.

Conflicts of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Sherchan, S.P.; Shahin, S.; Ward, L.M.; Tandukar, S.; Aw, T.G.; Schmitz, B.; Ahmed, W.; Kitajima, M. First detection of SARS-CoV-2 RNA in wastewater in North America: A study in Louisiana, USA. Sci. Total Environ. 2020, 743, 140621. [Google Scholar] [CrossRef] [PubMed]
  2. Ahmed, W.; Angel, N.; Edson, J.; Bibby, K.; Bivins, A.; O’Brien, J.W.; Choi, P.M.; Kitajima, M.; Simpson, S.L.; Li, J.; et al. First confirmed detection of SARS-CoV-2 in untreated wastewater in Australia: A proof of concept for the wastewater surveillance of COVID-19 in the community. Sci. Total Environ. 2020, 728, 138764. [Google Scholar] [CrossRef]
  3. Medema, G.; Heijnen, L.; Elsinga, G.; Italiaander, R.; Brouwer, A. Presence of SARS-Coronavirus-2 RNA in Sewage and Correlation with Reported COVID-19 Prevalence in the Early Stage of the Epidemic in The Netherlands. Environ. Sci. Technol. Lett. 2020, 7, 511–516. [Google Scholar] [CrossRef]
  4. Markt, R.; Endler, L.; Amman, F.; Schedl, A.; Penz, T.; Büchel-Marxer, M.; Grünbacher, D.; Mayr, M.; Peer, E.; Pedrazzini, M.; et al. Detection and abundance of SARS-CoV-2 in wastewater in Liechtenstein, and the estimation of prevalence and impact of the B.1.1.7 variant. J. Water Health 2022, 20, 114–125. [Google Scholar] [CrossRef] [PubMed]
  5. Graham, K.E.; Loeb, S.K.; Wolfe, M.K.; Catoe, D.; Sinnott-Armstrong, N.; Kim, S.; Yamahara, K.M.; Sassoubre, L.M.; Mendoza Grijalva, L.M.; Roldan-Hernandez, L.; et al. SARS-CoV-2 RNA in Wastewater Settled Solids Is Associated with COVID-19 Cases in a Large Urban Sewershed. Environ. Sci. Technol. 2021, 55, 488–498. [Google Scholar] [CrossRef]
  6. Pöyry, T.; Stenvik, M.; Hovi, T. Viruses in sewage waters during and after a poliomyelitis outbreak and subsequent nationwide oral poliovirus vaccination campaign in Finland. Appl. Environ. Microbiol. 1988, 54, 371–374. [Google Scholar] [CrossRef] [Green Version]
  7. Berchenko, Y.; Manor, Y.; Freedman, L.S.; Kaliner, E.; Grotto, I.; Mendelson, E.; Huppert, A. Estimation of polio infection prevalence from environmental surveillance data. Sci. Transl. Med. 2017, 9, eaaf6786. [Google Scholar] [CrossRef]
  8. Beyer, S.; Szewzyk, R.; Gnirss, R.; Johne, R.; Selinka, H.-C. Detection and Characterization of Hepatitis E Virus Genotype 3 in Wastewater and Urban Surface Waters in Germany. Food Environ. Virol. 2020, 12, 137–147. [Google Scholar] [CrossRef] [Green Version]
  9. Kazama, S.; Miura, T.; Masago, Y.; Konta, Y.; Tohma, K.; Manaka, T.; Liu, X.; Nakayama, D.; Tanno, T.; Saito, M.; et al. Environmental Surveillance of Norovirus Genogroups I and II for Sensitive Detection of Epidemic Variants. Appl. Environ. Microbiol. 2017, 83, e03406-16. [Google Scholar] [CrossRef] [Green Version]
  10. Huang, Y.; Zhou, N.; Zhang, S.; Yi, Y.; Han, Y.; Liu, M.; Han, Y.; Shi, N.; Yang, L.; Wang, Q.; et al. Norovirus detection in wastewater and its correlation with human gastroenteritis: A systematic review and meta-analysis. Environ. Sci. Pollut. Res. Int. 2022, 29, 22829–22842. [Google Scholar] [CrossRef]
  11. Mercier, E.; D’Aoust, P.M.; Thakali, O.; Hegazy, N.; Jia, J.-J.; Zhang, Z.; Eid, W.; Plaza-Diaz, J.; Kabir, M.P.; Fang, W.; et al. Municipal and neighbourhood level wastewater surveillance and subtyping of an influenza virus outbreak. Sci. Rep. 2022, 12, 15777. [Google Scholar] [CrossRef] [PubMed]
  12. Wolfe, M.K.; Duong, D.; Bakker, K.M.; Ammerman, M.; Mortenson, L.; Hughes, B.; Arts, P.; Lauring, A.S.; Fitzsimmons, W.J.; Bendall, E.; et al. Wastewater-Based Detection of Two Influenza Outbreaks. Environ. Sci. Technol. Lett. 2022, 9, 687–692. [Google Scholar] [CrossRef]
  13. Corpuz, M.V.A.; Buonerba, A.; Vigliotta, G.; Zarra, T.; Ballesteros, F.; Campiglia, P.; Belgiorno, V.; Korshin, G.; Naddeo, V. Viruses in wastewater: Occurrence, abundance and detection methods. Sci. Total Environ. 2020, 745, 140910. [Google Scholar] [CrossRef] [PubMed]
  14. Mao, K.; Zhang, K.; Du, W.; Ali, W.; Feng, X.; Zhang, H. The potential of wastewater-based epidemiology as surveillance and early warning of infectious disease outbreaks. Curr. Opin. Environ. Sci. Health 2020, 17, 1–7. [Google Scholar] [CrossRef]
  15. Thai, P.K.; O’Brien, J.W.; Banks, A.P.W.; Jiang, G.; Gao, J.; Choi, P.M.; Yuan, Z.; Mueller, J.F. Evaluating the in-sewer stability of three potential population biomarkers for application in wastewater-based epidemiology. Sci. Total Environ. 2019, 671, 248–253. [Google Scholar] [CrossRef] [PubMed]
  16. Maal-Bared, R.; Qiu, Y.; Li, Q.; Gao, T.; Hrudey, S.E.; Bhavanam, S.; Ruecker, N.J.; Ellehoj, E.; Lee, B.E.; Pang, X. Does normalization of SARS-CoV-2 concentrations by Pepper Mild Mottle Virus improve correlations and lead time between wastewater surveillance and clinical data in Alberta (Canada) comparing twelve SARS-CoV-2 normalization approaches. Sci. Total Environ. 2022, 856, 158964. [Google Scholar] [CrossRef]
  17. Sakarovitch, C.; Schlosser, O.; Courtois, S.; Proust-Lima, C.; Couallier, J.; Pétrau, A.; Litrico, X.; Loret, J.-F. Monitoring of SARS-CoV-2 in wastewater: What normalisation for improved understanding of epidemic trends? J. Water Health 2022, 20, 712–726. [Google Scholar] [CrossRef]
  18. Stephens, N.; Béen, F.; Savic, D. An Analysis of SARS-CoV-2 in Wastewater to Evaluate the Effectiveness of Nonpharmaceutical Interventions against COVID-19 in The Netherlands. ACS EST Water 2022, 2, 2158–2166. [Google Scholar] [CrossRef]
  19. Naughton, C.C.; Roman, F.A.; Alvarado, A.G.F.; Tariqi, A.Q.; Deeming, M.A.; Bibby, K.; Bivins, A.; Rose, J.B.; Medema, G.; Ahmed, W.; et al. Show us the Data: Global COVID-19 Wastewater Monitoring Efforts, Equity, and Gaps. medRxiv 2021. preprint. [Google Scholar] [CrossRef]
  20. Gwenzi, W. Wastewater, waste, and water-based epidemiology (WWW-BE): A novel hypothesis and decision-support tool to unravel COVID-19 in low-income settings? Sci. Total Environ. 2022, 806, 150680. [Google Scholar] [CrossRef]
  21. Stephenson, I.; Zambon, M. The epidemiology of influenza. Occup. Med. 2002, 52, 241–247. [Google Scholar] [CrossRef] [Green Version]
  22. Krammer, F.; Smith, G.J.D.; Fouchier, R.A.M.; Peiris, M.; Kedzierska, K.; Doherty, P.C.; Palese, P.; Shaw, M.L.; Treanor, J.; Webster, R.G.; et al. Influenza. Nat. Rev. Dis. Prim. 2018, 4, 3. [Google Scholar] [CrossRef]
  23. Yoo, S.J.; Moon, S.J.; Kuak, E.-Y.; Yoo, H.M.; Kim, C.K.; Chey, M.-J.; Shin, B.-M. Frequent detection of pandemic (H1N1) 2009 virus in stools of hospitalized patients. J. Clin. Microbiol. 2010, 48, 2314–2315. [Google Scholar] [CrossRef] [Green Version]
  24. Chan, M.C.W.; Lee, N.; Chan, P.K.S.; To, K.F.; Wong, R.Y.K.; Ho, W.-S.; Ngai, K.L.K.; Sung, J.J.Y. Seasonal influenza A virus in feces of hospitalized adults. Emerg. Infect. Dis. 2011, 17, 2038–2042. [Google Scholar] [CrossRef]
  25. Arena, C.; Amoros, J.P.; Vaillant, V.; Balay, K.; Chikhi-Brachet, R.; Varesi, L.; Arrighi, J.; Blanchon, T.; Carrat, F.; Hanslik, T.; et al. Simultaneous investigation of influenza and enteric viruses in the stools of adult patients consulting in general practice for acute diarrhea. Virol. J. 2012, 9, 116. [Google Scholar] [CrossRef] [Green Version]
  26. Minodier, L.; Charrel, R.N.; Ceccaldi, P.-E.; van der Werf, S.; Blanchon, T.; Hanslik, T.; Falchi, A. Prevalence of gastrointestinal symptoms in patients with influenza, clinical significance, and pathophysiology of human influenza viruses in faecal samples: What do we know? Virol. J. 2015, 12, 215. [Google Scholar] [CrossRef] [Green Version]
  27. Robilotti, E.; Deresinski, S.; Pinsky, B.A. Norovirus. Clin. Microbiol. Rev. 2015, 28, 134–164. [Google Scholar] [CrossRef] [Green Version]
  28. Pires, S.M.; Fischer-Walker, C.L.; Lanata, C.F.; Devleesschauwer, B.; Hall, A.J.; Kirk, M.D.; Duarte, A.S.R.; Black, R.E.; Angulo, F.J. Aetiology-Specific Estimates of the Global and Regional Incidence and Mortality of Diarrhoeal Diseases Commonly Transmitted through Food. PLoS ONE 2015, 10, e0142927. [Google Scholar] [CrossRef]
  29. Fretz, R.; Herrmann, L.; Christen, A.; Svoboda, P.; Dubuis, O.; Viollier, E.H.; Tanner, M.; Baumgartner, A. Frequency of Norovirus in stool samples from patients with gastrointestinal symptoms in Switzerland. Eur. J. Clin. Microbiol. Infect. Dis. 2005, 24, 214–216. [Google Scholar] [CrossRef]
  30. Lee, N.; Chan, M.C.W.; Wong, B.; Choi, K.W.; Sin, W.; Lui, G.; Chan, P.K.S.; Lai, R.W.M.; Cockram, C.S.; Sung, J.J.Y.; et al. Fecal viral concentration and diarrhea in norovirus gastroenteritis. Emerg. Infect. Dis. 2007, 13, 1399–1401. [Google Scholar] [CrossRef]
  31. Aoki, Y.; Suto, A.; Mizuta, K.; Ahiko, T.; Osaka, K.; Matsuzaki, Y. Duration of norovirus excretion and the longitudinal course of viral load in norovirus-infected elderly patients. J. Hosp. Infect. 2010, 75, 42–46. [Google Scholar] [CrossRef] [PubMed]
  32. Phillips, G.; Tam, C.C.; Rodrigues, L.C.; Lopman, B. Prevalence and characteristics of asymptomatic norovirus infection in the community in England. Epidemiol. Infect. 2010, 138, 1454–1458. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Amman, F.; Markt, R.; Endler, L.; Hupfauf, S.; Agerer, B.; Schedl, A.; Richter, L.; Zechmeister, M.; Bicher, M.; Heiler, G.; et al. Viral variant-resolved wastewater surveillance of SARS-CoV-2 at national scale. Nat. Biotechnol. 2022, 40, 1814–1822. [Google Scholar] [CrossRef] [PubMed]
  34. Markt, R.; Mayr, M.; Peer, E.; Wagner, A.O.; Lackner, N.; Insam, H. Detection and Stability of SARS-CoV-2 Fragments in Wastewater: Impact of Storage Temperature. Pathogens 2021, 10, 1215. [Google Scholar] [CrossRef]
  35. Centers for Disease Control and Prevention (CDC). Influenza SARS-CoV-2 (Flu SC2) Multiplex Assay Primers and Probes; Centers for Disease Control and Prevention: Atlanta, GA, USA, 2021.
  36. Liu, D.; Zhang, Z.; Wu, Q.; Tian, P.; Geng, H.; Xu, T.; Wang, D. Redesigned Duplex RT-qPCR for the Detection of GI and GII Human Noroviruses. Engineering 2020, 6, 442–448. [Google Scholar] [CrossRef]
  37. Been, F.; Rossi, L.; Ort, C.; Rudaz, S.; Delémont, O.; Esseiva, P. Population normalization with ammonium in wastewater-based epidemiology: Application to illicit drug monitoring. Environ. Sci. Technol. 2014, 48, 8162–8169. [Google Scholar] [CrossRef]
  38. Rauch, W.; Schenk, H.; Insam, H.; Markt, R.; Kreuzinger, N. Data modelling recipes for SARS-CoV-2 wastewater-based epidemiology. Environ. Res. 2022, 214, 113809. [Google Scholar] [CrossRef]
  39. Arabzadeh, R.; Grünbacher, D.M.; Insam, H.; Kreuzinger, N.; Markt, R.; Rauch, W. Data filtering methods for SARS-CoV-2 wastewater surveillance. Water Sci. Technol. 2021, 84, 1324–1339. [Google Scholar] [CrossRef]
  40. R Core Team. R: A Language and Environment for Statistical Computing; R Core Team: Vienna, Austria, 2021; Available online: https://www.R-project.org/ (accessed on 3 August 2022).
  41. Dancho, M.; Vaughan, D. tidyquant: Tidy Quantitative Financial Analysis. 2022. Available online: https://CRAN.R-project.org/package=tidyquant (accessed on 3 August 2022).
  42. Wickham, H.; Bryan, J. readxl: Read Excel Files. 2022. Available online: https://CRAN.R-project.org/package=readxl (accessed on 3 August 2022).
  43. Wickham, H. ggplot2: Elegant Graphics for Data Analysis; Springer: New York, NY, USA, 2016; ISBN 978-3-319-24277-4. [Google Scholar]
  44. Wilke, C.O. cowplot: Streamlined Plot Theme and Plot Annotations for ‘ggplot2’. 2020. Available online: https://CRAN.R-project.org/package=cowplot (accessed on 3 August 2022).
  45. Davis, W.W.; Mott, J.A.; Olsen, S.J. The role of non-pharmaceutical interventions on influenza circulation during the COVID-19 pandemic in nine tropical Asian countries. Influenza Other Respir. Viruses 2022, 16, 568–576. [Google Scholar] [CrossRef]
  46. Baker, R.E.; Park, S.W.; Yang, W.; Vecchi, G.A.; Metcalf, C.J.E.; Grenfell, B.T. The impact of COVID-19 nonpharmaceutical interventions on the future dynamics of endemic infections. Proc. Natl. Acad. Sci. USA 2020, 117, 30547–30553. [Google Scholar] [CrossRef]
  47. Bruggink, L.D.; Garcia-Clapes, A.; Tran, T.; Druce, J.D.; Thorley, B.R. Decreased incidence of enterovirus and norovirus infections during the COVID-19 pandemic, Victoria, Austraalia, 2020. Commun. Dis. Intell. (2018) 2020, 44, 1–8. [Google Scholar] [CrossRef]
  48. Lennon, R.P.; Griffin, C.; Miller, E.L.; Dong, H.; Rabago, D.; Zgierska, A.E. Norovirus Infections Drop 49% in the United States with Strict COVID-19 Public Health Interventions. Acta Med. Acad. 2020, 49, 278–280. [Google Scholar] [CrossRef] [PubMed]
  49. Han, S.; Zhang, T.; Lyu, Y.; Lai, S.; Dai, P.; Zheng, J.; Yang, W.; Zhou, X.; Feng, L. The Incoming Influenza Season—China, the United Kingdom, and the United States, 2021–2022. China CDC Wkly. 2021, 3, 1039–1045. [Google Scholar] [CrossRef] [PubMed]
  50. Kauppinen, A.; Miettinen, I.T. Persistence of Norovirus GII Genome in Drinking Water and Wastewater at Different Temperatures. Pathogens 2017, 6, 48. [Google Scholar] [CrossRef]
  51. Desai, R.; Hembree, C.D.; Handel, A.; Matthews, J.E.; Dickey, B.W.; McDonald, S.; Hall, A.J.; Parashar, U.D.; Leon, J.S.; Lopman, B. Severe outcomes are associated with genogroup 2 genotype 4 norovirus outbreaks: A systematic literature review. Clin. Infect. Dis. 2012, 55, 189–193. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Vega, E.; Barclay, L.; Gregoricus, N.; Williams, K.; Lee, D.; Vinjé, J. Novel surveillance network for norovirus gastroenteritis outbreaks, United States. Emerg. Infect. Dis. 2011, 17, 1389–1395. [Google Scholar] [CrossRef]
  53. Nordgren, J.; Matussek, A.; Mattsson, A.; Svensson, L.; Lindgren, P.-E. Prevalence of norovirus and factors influencing virus concentrations during one year in a full-scale wastewater treatment plant. Water Res. 2009, 43, 1117–1125. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Map of Austria and its political districts. Colored districts correspond to those containing the 6 WWTPs monitored in Salzburg (blue) and the 5 WWTPs in Vorarlberg (red), covering in total 0.7 Mio inhabitants.
Figure 1. Map of Austria and its political districts. Colored districts correspond to those containing the 6 WWTPs monitored in Salzburg (blue) and the 5 WWTPs in Vorarlberg (red), covering in total 0.7 Mio inhabitants.
Viruses 15 00263 g001
Figure 2. Normalized virus loads for (A) influenza A, (B) norovirus GII, and (C) norovirus GI per WWTP. More intense color corresponds to higher loads in wastewater. Each plot is scaled separately due to large discrepancy of absolute loads between pathogens. WWTP V001–V005 correspond to the federal state Vorarlberg and WWTPs Z001–Z012 correspond to the federal state Salzburg.
Figure 2. Normalized virus loads for (A) influenza A, (B) norovirus GII, and (C) norovirus GI per WWTP. More intense color corresponds to higher loads in wastewater. Each plot is scaled separately due to large discrepancy of absolute loads between pathogens. WWTP V001–V005 correspond to the federal state Vorarlberg and WWTPs Z001–Z012 correspond to the federal state Salzburg.
Viruses 15 00263 g002
Figure 3. Seasonal occurrence of (red line) influenza virus and (blue line) norovirus infections in inpatients aggregated per month across Austria. The dashed black line marks the date of the first clinical COVID-19 case in Austria.
Figure 3. Seasonal occurrence of (red line) influenza virus and (blue line) norovirus infections in inpatients aggregated per month across Austria. The dashed black line marks the date of the first clinical COVID-19 case in Austria.
Viruses 15 00263 g003
Figure 4. Comparison between viral loads in wastewater and number of infected inpatients from September 2021–March 2022 including time points for non-pharmaceutical intervention for IAV (A) and NoV GII (B). (Blue dots) Viral loads and (red dots) number of inpatients were aggregated at a weekly basis. Smoothing lines with 95% confidence bands correspond to LOWESS-smoothing. Six events increasing or decreasing the impact of NPIs on the spread of the COVID-19 pandemic were set from September 2021 to March 2022: September 9–school opens, (I) November 15—lockdown for those unvaccinated, (II) November 22nd—complete lockdown, (III) December 12—partial opening, (IV) January 11—Omicron measures, (V) March 5—complete opening.
Figure 4. Comparison between viral loads in wastewater and number of infected inpatients from September 2021–March 2022 including time points for non-pharmaceutical intervention for IAV (A) and NoV GII (B). (Blue dots) Viral loads and (red dots) number of inpatients were aggregated at a weekly basis. Smoothing lines with 95% confidence bands correspond to LOWESS-smoothing. Six events increasing or decreasing the impact of NPIs on the spread of the COVID-19 pandemic were set from September 2021 to March 2022: September 9–school opens, (I) November 15—lockdown for those unvaccinated, (II) November 22nd—complete lockdown, (III) December 12—partial opening, (IV) January 11—Omicron measures, (V) March 5—complete opening.
Viruses 15 00263 g004
Table 1. WWTPs monitored in this study including relative coverage corresponding to total population within the two federal states Vorarlberg and Salzburg, and the number of inhabitants connected to each sewage system. Calculated populations across all sampling points per WWTPs are stated as theoretical inhabitants after NH4-N normalization. Additionally, the ratio between theoretical and registered inhabitants was calculated in %.
Table 1. WWTPs monitored in this study including relative coverage corresponding to total population within the two federal states Vorarlberg and Salzburg, and the number of inhabitants connected to each sewage system. Calculated populations across all sampling points per WWTPs are stated as theoretical inhabitants after NH4-N normalization. Additionally, the ratio between theoretical and registered inhabitants was calculated in %.
Federal State
(Population Covered by the WWTPs Monitored)
WWTP IDInhabitants Connected According Civil RegisterTheoretical Population Calculated by NH4-N-Normalization
(Mean ± SD)
Calculated Population/
Registered Population
Vorarlberg
(~62% coverage)
V00166,96169,363 ± 13,908104%
V00257,46562,984 ± 13,152110%
V00351,36675,197 ± 14,846146%
V00441,69644,686 ± 8655107%
V00532,47132,373 ± 11,170100%
Salzburg
(~80% coverage)
Z001320,682273,966 ± 41,85185%
Z00238,60665,808 ± 10,666170%
Z00331,63838,317 ± 10,562121%
Z00622,44722,522 ± 3721100%
Z00820,37016,439 ± 174981%
Z01214,20016,491 ± 4078116%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Markt, R.; Stillebacher, F.; Nägele, F.; Kammerer, A.; Peer, N.; Payr, M.; Scheffknecht, C.; Dria, S.; Draxl-Weiskopf, S.; Mayr, M.; et al. Expanding the Pathogen Panel in Wastewater Epidemiology to Influenza and Norovirus. Viruses 2023, 15, 263. https://doi.org/10.3390/v15020263

AMA Style

Markt R, Stillebacher F, Nägele F, Kammerer A, Peer N, Payr M, Scheffknecht C, Dria S, Draxl-Weiskopf S, Mayr M, et al. Expanding the Pathogen Panel in Wastewater Epidemiology to Influenza and Norovirus. Viruses. 2023; 15(2):263. https://doi.org/10.3390/v15020263

Chicago/Turabian Style

Markt, Rudolf, Fabian Stillebacher, Fabiana Nägele, Anna Kammerer, Nico Peer, Maria Payr, Christoph Scheffknecht, Silvina Dria, Simon Draxl-Weiskopf, Markus Mayr, and et al. 2023. "Expanding the Pathogen Panel in Wastewater Epidemiology to Influenza and Norovirus" Viruses 15, no. 2: 263. https://doi.org/10.3390/v15020263

APA Style

Markt, R., Stillebacher, F., Nägele, F., Kammerer, A., Peer, N., Payr, M., Scheffknecht, C., Dria, S., Draxl-Weiskopf, S., Mayr, M., Rauch, W., Kreuzinger, N., Rainer, L., Bachner, F., Zuba, M., Ostermann, H., Lackner, N., Insam, H., & Wagner, A. O. (2023). Expanding the Pathogen Panel in Wastewater Epidemiology to Influenza and Norovirus. Viruses, 15(2), 263. https://doi.org/10.3390/v15020263

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