Occurrence of human-associated Bacteroidetes genetic source tracking markers in raw and treated wastewater of municipal and domestic origin and comparison to standard and alternative indicators of faecal pollution.

This was a detailed investigation of the seasonal occurrence, dynamics, removal and resistance of human-associated genetic Bacteroidetes faecal markers (GeBaM) compared with ISO-based standard faecal indicator bacteria (SFIB), human-specific viral faecal markers and one human-associated Bacteroidetes phage in raw and treated wastewater of municipal and domestic origin. Characteristics of the selected activated sludge wastewater treatment plants (WWTPs) from Austria and Germany were studied in detail (WWTPs, n = 13, connected populations from 3 to 49000 individuals), supported by volume-proportional automated 24-h sampling and chemical water quality analysis. GeBaM were consistently detected in high concentrations in raw (median log10 8.6 marker equivalents (ME) 100 ml(-1)) and biologically treated wastewater samples (median log10 6.2-6.5 ME 100 ml(-1)), irrespective of plant size, type and time of the season (n = 53-65). GeBaM, Escherichia coli, and enterococci concentrations revealed the same range of statistical variability for raw (multiplicative standard deviations s* = 2.3-3.0) and treated wastewater (s* = 3.7-4.5), with increased variability after treatment. Clostridium perfringens spores revealed the lowest variability for raw wastewater (s* = 1.5). In raw wastewater correlations among microbiological parameters were only detectable between GeBaM, C. perfringens and JC polyomaviruses. Statistical associations amongst microbial parameters increased during wastewater treatment. Two plants with advanced treatment were also investigated, revealing a minimum log10 5.0 (10th percentile) reduction of GeBaM in the activated sludge membrane bioreactor, but no reduction of the genetic markers during UV irradiation (254 nm). This study highlights the potential of human-associated GeBaM to complement wastewater impact monitoring based on the determination of SFIB. In addition, human-specific JC polyomaviruses and adenoviruses seem to be a valuable support if highly specific markers are needed.

Introduction

Contamination of aquatic systems by wastewater of human origin can pose a serious threat to public health because it frequently contains high numbers of intestinal pathogens (Stevens et al., 2009). Appropriate disposal systems combined with efficient faecal pollution monitoring techniques for municipal and domestic wastewater are thus essential for safeguarding our water resources. Wastewater treatment plants (WWTPs) based on primary (mechanical), secondary (biological), and tertiary (enhanced biological and chemical) treatment are designed to remove organic carbon (C), nitrogen (N) and phosphorus (P) out of wastewater to a great extent. Although providing a first essential barrier, conventional WWTPs are not built to sufficiently remove microbial faecal loads to support the safe use of effluent wastewater for human related activities, such as recreational purposes or irrigation. Disinfection of wastewater effluents has not yet become a common standard in most regions of the world, and such advanced treatment is often restricted to the discharge of wastewater into sensitive aquatic areas. Rainfall events may also lead to a bypass of WWTPs (i.e. combined sewer overflows) and the contamination of water resources with raw wastewater (Molina et al., 2014; Shibata et al., 2014; Tryland et al., 2014).

Routine monitoring of microbial faecal pollution in the aquatic environment is still based on the selective cultivation of standard faecal indicator bacteria (SFIB), including Escherichia coli and intestinal enterococci (ISO, 2005). Without doubt, water quality testing based on the application of SFIB has contributed to a fundamental improvement in water safety management since the end of the 19th century (Tallon et al., 2005). However, the application of SFIB has also recently been subjected to increasing criticism (Ishii and Sadowsky, 2008). Several studies suggested that SFIB in aquatic habitats also originate from non-enteric compartments, such as soil, sediment and algae (Byappanahalli et al., 2012; Desmarais et al., 2001; Whitman et al., 2003). In addition, typing of isolated SFIB strains, to recover information about their origin, requires the formation of unrealistically large catchment-specific strain libraries and does currently not provide a feasible option for monitoring applications (Domingo et al., 2007). These limitations obviously call for additional indicators and tools to complement the existing standard methods to obtain a more detailed and certain view on the existing faecal pollution patterns to support MST and risk assessment (Harwood et al., 2014).

Amongst the vast number of alternative parameters (Hagedorn et al., 2011; Wuertz et al., 2011), PCR-based assays for the analysis of human-associated genetic Bacteroidetes faecal markers (GeBaM) have gained increasing popularity in the field of faecal pollution analysis and microbial source tracking (MST) during recent years (Harwood et al., 2014). Quantitative PCR (qPCR)-based GeBaM assays for general-, human-, wastewater-, or animal-associated faecal sources have been developed (Kildare et al., 2007; Layton et al., 2006; Reischer et al., 2006; Shanks et al., 2009). Several evaluation studies including various aquatic environments successfully demonstrated the value of GeBaM diagnostics (Boehm et al., 2009; Reischer et al., 2011; Ridley et al., 2014; Riedel et al., 2014; Sauer et al., 2011; Tambalo et al., 2012). However, the application of qPCR-based GeBaM assays is not yet standardized. It requires careful study design and background information on the catchment to create unbiased results and to recognize methodical limits (Boehm et al., 2013; Reischer et al., 2011).

A useful parameter for the analysis of general- or host associated microbial faecal pollution in water has to fulfil several basic performance criteria, including source-sensitivity and source-specificity (Wuertz et al., 2011). Considerable effort has been dedicated to sensitivity and specificity testing of GeBaM qPCR assays during recent years, most frequently based on individual sampling strategies covering various sources of animal and human excreta or wastewater (Ahmed et al., 2013; Boehm et al., 2013; Keity et al., 2012; Reischer et al., 2013, 2011; Riedel et al., 2014; Shanks et al., 2009). Emphasis has also been put on sampling techniques, DNA extraction, and PCR quantification procedures (Cankar et al., 2006; Karlen et al., 2007; Shanks et al., 2012; Siefring et al., 2008; Stoeckel et al., 2009). However, information on the occurrence of GeBaM in wastewater regarding the characteristics of the disposal system (combined and separate sewer systems), its seasonal variability, and its relationship to standard and alternative faecal indicators is scarce (Srinivasan et al., 2011).

The aim of this study was to investigate the prevalence and abundance of human-associated GeBaM by qPCR determination in raw and treated wastewater of well-characterized municipal wastewater treatment plants over one year. Emphasis was put on the selection of municipal WWTPs with primary, secondary, and tertiary treatment, as such systems are representative for the situation in Austria and the Central European Region (CER). Small domestic WWTP (dWWTPs) were also included in our investigation, as they are frequently implemented in remote areas, where the connection to municipal sewer systems is not possible. Although advanced treatment was not the main focus of this study, the investigation of UV disinfection at one selected WWTP was included, as such treatment is becoming increasingly important. The Taqman HF183 qPCR assay (Haugland et al., 2010) and the BacHUM qPCR assay (Kildare et al., 2007) were selected for the determination of human-associated GeBaM concentrations, following recommendations of recent evaluation studies (Boehm et al., 2013; Layton et al., 2013; Reischer et al., 2013). To support methodical cross-comparisons, cultivation-based SFIB using ISO standard methods and viral faecal markers for human-specific faecal pollution were simultaneously determined. Among these, JC polyomavirus (JCPyV) as well as human adenoviruses (HAdV), which have been used as human faecal viral indicators and highly specific MST tools (Bofill-Mas et al., 2000; Pina et al., 1998), and bacteriophages infecting Bacteroides thetaiotaomicron, which have been proposed as a human-associated faecal indicator, were tested. Raw and treated wastewater was also analysed by chemical standard methods to support treatment plant characterisation and comparison of elimination characteristics between microbial and chemical parameters.

Materials and methods

Selection criteria and parameters to characterize the sewer disposal systems and WWTPs

The Danube Region and other parts of the CER are defined as sensitive areas with respect to water bodies. In terms of nutrients, strict discharge limits for WWTPs according to the EU urban wastewater directive have been established (EC, 1991). Wastewater disposal is caused by very small (i.e., a few inhabitants) up to large treatment systems (>100 000 persons connected), as a mix of rural areas and large cities characterize this region. Activated sludge treatment is the common process to treat the wastewater. In Austria, WWTP effluent concentrations are restricted by certain removal performance targets related to the influent load (%) and by a maximum effluent concentration (mg L−1). These limits depend on the plant size (AEV, 1996). Only municipal WWTPs providing data for a basic characterization over the investigation period were selected. Essential information on WWTP design, including design capacity (p.e., population equivalent), actual average loading inhabitants connected, type of treatment (mechanical (M), carbon removal (C), nitrification (N), denitrification (D), phosphorus removal (P)), advanced treatment available (UV irradiation, membrane filtration), removal efficiency of nutrients (C, N and P), and sludge age, were required. To characterize the raw wastewater quantity and quality over the investigation period, volumetric flow rate (Q), chemical oxygen demand (COD) or biological oxygen demand (BOD), total nitrogen (TN), total phosphorus (TP), temperature (°C), and pH were requested. In addition to these parameters, ammonium (NH4), nitrate (NO3), and suspended solids (SS) were included for the treated wastewater. All of these flow and chemical data were provided on a daily basis because they were available from representative sampling by automated and cooled sampling devices (24 h proportional-flow sampling). All chemical information provided was cross-checked by our own investigations of the investigated samples (see below).

In contrast to municipal WWTP, data availability for dWWTPs (<150 p.e.) was very low. The basic requirements were information on the type of treatment system, the total number of persons connected, and the effluent concentrations of COD, NH4 and pH. Only information based on grab sample analysis was available.

Sampling for chemical and microbiological analysis

60 influent and 60 effluent 24 h volume-proportional composite samples were recovered by fixed installed and cooled (4 °C) automatic sampling devices (cf. supplemental material) from raw and treated municipal wastewater. Further background information on automated sampling in WWTP can also be found in Mayer et al. (Mayer et al., 2015). 2 L of samples were collected during the automated sampling process in sterilized 2 L glass bottles and immediately transferred to the laboratory at 5 ± 3 °C for further analysis within 8 h. To cover seasonal variations, samples were taken in 4- to 6-week intervals over one annual cycle. Following a homogenization by manual shaking, 1 L of the sample volume of each sample was used for chemical and microbiological analysis.

In summer and fall (2013) 16 grab samples were taken at the effluent discharge unit of the selected dWWTPs. 2 L of samples were collected in sterilized 2 L glass bottles and immediately transferred to the laboratory at 5 ± 3 °C for further analysis within 8 h. Following a homogenization by manual shaking, a volume of 1 L was used for chemical and microbiological investigations each.

To evaluate the effect of UV irradiation on GeBaM and other microbiological indicators, 10 additional pairs of 1 L grab samples (before and after the UV system) were taken at the WWTP5 throughout the season (UV-2000, Trojan Technologies, Canada; 48 UV lamps).

Chemical analysis

For COD, BOD5, TP, TN, pH, SS, and conductivity analysis, the first preparation step included the homogenization of the sample by manual agitation. A pre-filtration step, applying a 0.45 µm membrane (sterilized cellulose-nitrate filter), was needed to analyse the dissolved parameters PO4–P, NH4–N, NO2–N, NOx-N. All selected parameters were performed according to standardized methods, as given in detail in the additional material section (Table S2 additional material).

Quantification of genetic Bacteroidetes markers (GeBaM) by qPCR

For DNA extraction, polycarbonate membrane filtration (0.2 µm Millipore, Isopore Membrane Filter - GTTP) followed by phenol/chloroform extraction of 20 ml influent, 50 ml effluent and 1500 ml membrane-filtrated effluent, was used as previously described (Griffiths et al., 2000; Reischer et al., 2006). Every filtration event was monitored by a minimum of two blank filter controls and every extraction event (equals 11 samples) was also followed by a blank extraction control. The controls were processed in exactly the same way as the samples. Cells were lysed with a FastPrepR-24 Instrument (MP Biomedicals Inc., Irvine, USA) with a speed setting of 6 for 30 s. The extracted DNA was dissolved in 10 mM TRIS HCl pH = 8 and stored at −80 °C not longer than 21 days before PCR testing. Respective 16S-rRNA-gene markers for AllBac (Layton et al., 2006), BacHUM (Kildare et al., 2007), and HF183 TaqMan (Haugland et al., 2010) were quantified by qPCR. The rotor-discs and 96-well plates were loaded with the mastermix and the sample DNA using a Qiagility Roboter (Qiagen, Hilden, Germany). The measurements were subsequently performed on a Rotorgene Q Cycler (Qiagen, HILDEN, Germany). For the AllBac qPCR assay we used 2.5 µl of the respective DNA sample dilution, 600 nmol L−1 primer AllBac296f, 600 nmol L−1 primer AllBac412r, 25 nmol L−1 TaqMan MGB probe AllBac375Bhqr (Layton et al., 2006), 0.4 g L−1 bovine serum albumin (Roche Diagnostics, Mannheim, Germany), 7.5 µl of iQ Supermix (Bio-Rad, Hercules, USA) in a total reaction volume of 15 µl; additionally, 5 mmol L−1 MgCl2 was added to obtain a total Mg2+ concentration of 8 mmol L−1 (Layton et al., 2006). For the BacHUM assay we used 2.5 µl of the respective DNA sample dilution, 400 nmol L−1 primer BacHUM-160f, 400 nmol L−1 primer BacHUM-241r, 80 nmol L−1 TaqMan MGB probeeBacHUM-193p (Kildare et al., 2007), 0.4 g L−1 bovine serum albumin, and 7.5 µl of iQ Supermix in a total reaction volume of 15 ml. For the HF183 TaqMan assay, we used 2.5 µl of the respective DNA sample-dilution, 100 nmol L−1 primer HF183, 100 nmol L−1 primer BFD REV, 80 nmol L−1 TaqMan MGB probe eBFD FAM (Haugland et al., 2010) 0.4 g l L−1 bovine serum albumin, and 7.5 ml of iQ Supermix in a total reaction volume of 15 ml. The PCR program for the AllBac assay was initial 95 °C for 3 min and for the cycles of 95 °C for 30 s and 60 °C for 45 s; for the BacHUM assay, initial 95 °C for 3 min and for the cycles of 95 °C for 15 s and 60 °C for 1 min; for the HF183 TaqMan assay, initial 95 °C for 3 min and for the cycles of 95 °C for 15 s and 60 °C for 30 s. The real-time data were collected during the primer annealing step at 60 °C. Quantification was based on plasmid standard dilutions and given as molecular equivalent targets per volume (ME vo−1) as previously described (Reischer et al., 2006). The respective plasmid stock (cf. supplemental material) for each assay was diluted in an unspecific 500 ng ml−1 poly(dI-dC) (Roche Diagnostics, Mannheim, Germany) background to avoid adsorption of plasmid DNA to reaction vials at low plasmid concentrations. A total of at least seven tenfold serial dilutions of plasmid standard (100–106 gene copies) were run in each qPCR experiment.

Each DNA sample was analysed in two dilution steps (10- and 100-fold dilution) with each dilution in duplicate to check for a possible PCR inhibition. In case the quotient of the concentration ratio (after correction with the dilution factor) was in between the range of 0.5–2, inhibition was assumed to be absent (i.e. correspondence of results from different dilutions). No signs of PCR inhibition could be detected for any of the applied qPCR assays. Inhibition was also analysed by dilution plot analysis (cf. supplemental material) in which the concentrations of the two dilutions are shown on the x vs. y axes. Dilution plot analysis resulted in a very good correlation (R2) and a 1:1 slope amongst the concentrations, as revealed for the different dilutions (cf. Fig. S1, supplementary material). Every qPCR run also included a no template control in duplicates (only mastermix, primers and probe) to monitor for possible contaminations. All qPCR runs in this study revealed a calculated PCR efficiency of between 90% and 105% and coefficient of determination (R2) higher than 0.99 for the calibration curve. The no-template controls were consistently negative. The range of quantification was linked to the calibration curve covering concentrations in between 102 - 107 targets per reaction. To introduce a process control for the filtration and extraction step, an E. coli cell standard carrying a diagnostic gene fragment for qPCR detection at the chromosome was applied. The so called defined genetic target number cell standard (DeTaCs) was directly spiked into 50% of the collected wastewater samples as previously described (Kaiblinger, 2008) at concentrations of log10 7.0 (details on the design, production and detection of the DeTaCs are given in detail in the supplementary material section). The low variability of the concentrations obtained from the DeTaCs spikes by qPCR further proved the reliability and reproducibility of the filtration and extraction process. The multiplicative standard deviation for the compiled data set resulted in a value of s* = 2.1.

Cultivation-based enumeration of standard faecal indicator bacteria (SFIB)

Cultivation-based enumeration of SFIB (i.e., E. coli, Enterococci and Clostridium perfringens spores) was performed in the frame of our ISO 17025 accreditation. Before analysis, the samples were homogenized in an ultrasonic bath for 5 min (Bandelin Sonorex, RK 100, 35 kHz; Berlin, Germany). For membrane filtration, appropriate dilutions were performed (Farnleitner et al., 2010; Vierheilig et al., 2013). Enumeration of presumptive E. coli was based on the ISO standard 16649-1 (ISO, 2001a), using the chromogenic TBX agar (Oxoid, Thermo Fisher Scientific Inc., United Kingdom) and incubation at 44 ± 0.5 °C for 24 ± 0.5 h. Enumeration of Enterococci was based on the ISO standard 7899-2 (ISO, 2000), using Slanetz–Bartley medium (Oxoid) and incubation at 36 ± 2 °C for 44 ± 4 h. For quantification of C. perfringens spores, 5 ml influent and 15 ml effluent were pasteurized at 60 ± 2 °C for 15 min. C. perfringens was analysed according to the established ISO method 14189 (ISO, 2013), based on selective cultivation using TSC agar (Scharlau, Spain) at 44 ± 0.5 °C for 21 ± 3 h and subsequent identification of colonies by acid phosphatase reaction (Ryzinska-Paier et al., 2011). For quality assurance (media control as well as performance control of the methods), the reference strains E. coli NCTC 9001, Enterococcus faecalis NCTC 775 and Clostridium perfringens NCTC 8237 were used.

Quantification of human-specific viral faecal indicators by qPCR

To concentrate the desired viral DNA, 50 ml of influent and 500 ml of effluent were used. For membrane bioreactor 5000–10,000 ml of effluent grab samples, were used for the skimmed milk flocculation process as established by Calgua (Calgua et al., 2013a, 2008). Viral concentrates were resuspended in 1 ml of phosphate buffer (1:2, v/v of Na2HPO4 0.2 M and NaH2PO4 0.2 M) at pH 7.5. A control spike (adenovirus type 35) was also added as a process control to aliquoted composite influent and aliquoted composite effluent samples. The HAdV35 spiking stock was prepared in A549 cells and quantified by plaque assays and quantitative PCR (cf. for more details on plaque assay see supplemental material). A GC number of 105 HAdV35 were used for spiking. The qPCR assay used to quantify the spiked HAdV in the samples showed a limit of detection of 11 GC per reaction. According to the protocol of concentration applied, numbers lower than 102–103 GC/L were under the sample limit of detection. All process controls tested in this study were above the given sample limit of detection and well within the range of recovery rates, as previously reported for this specific protocol of concentration (Calgua et al., 2013b). The process of treatment for all samples investigated was thus considered acceptable (for more information see supplemental material) Ultrapure water (ISO 3696 grade 1) produced from tap water was used as negative control of the process. The water treatment system consists of a pretreatment cartridge followed by a reverse osmosis unit, UV irradiation chamber, a deionizing polishing cartridge and a 0.2 µm membrane filter (ELGA Medica-R 7, Veolia, Water Systems, Austria). Viral DNA was extracted from all samples using the QIAamp Viral RNA kit (Qiagen, Inc.). Nucleic acid eluates (holding time < 48 h at 21 °C). were sent at room temperature to the laboratory in Barcelona for quantification by qPCR. Specific real-time qPCR assays were used to quantify HAdV and JCPyV as previously described in detail (Bofill-Mas et al., 2006; Hernroth et al., 2002; Pal et al., 2006). Amplifications were performed in a 25-µl reaction mixture containing 10 µl of DNA and 15 µl of TaqMan Environmental Master Mix (Life Technologies) and 0.9 µM of each primer (AdF and AdR) and 0,225 µM of fluorogenic probe (AdP1) for HAdV detection. After activation of the AmpliTaq Gold for 10 min at 95 °C, 40 cycles (15 s at 95 °C and 1 min at 60 °C) were performed in a Stratagene Mx3000P detection system. DNA suspensions were analysed in duplicates using un-diluted and ten-fold dilutions to analyse environmental samples, whereas each dilution of standard DNA suspensions (cf. supplemental material) from 100 to 107 gene copies was run in triplicates. In all of the qPCRs carried out, the amount of DNA was defined as the mean of the data obtained. A non-template control (NTC) and a non-amplification control (NAC, DNA polymerase was omitted from the reaction to monitor background signal and probe stability) were added to each run. Tests for inhibition were performed according to the procedure described in section 2.5. To verify negative results known amounts of targeted viral DNA was added to negative samples and analysis was repeated.

Enumeration of human-associated bacteriophages infecting B. thetaiotaomicron

Bacteriophages infecting B. thetaiotaomicron were enumerated according to the standard method ISO 10705-4 (ISO, 2001b) as described for Bacteroides fragilis RYC2056 or HSP40 infecting phages. The host strain applied for the phage analyses was B. thetaiotaomicron (GA17), kindly provided by Prof. Maite Muniesa, University of Barcelona, Spain. To reduce the high concentration of background flora, samples were filtered through a low protein-binding membrane (0.2 µm; Minisart 16534, Sartorius). For enumeration of the phages, 1 ml of the host strain (inoculum culture) was added to tubes containing 2.5 ml of semisolid agar and aliquots of the samples to be tested, gently mixed, and poured onto solid agar plates. The plates were incubated under anaerobic conditions (AnaeroGen AN0025A, Oxoid) at 36 ± 2 °C for 21 ± 3 h. The results are expressed as the number of plaque-forming units (pfu) per sample volume.

Data analysis and statistics

All microbial data are expressed as log10 (x+1), after having performed all the needed calculations on the untransformed data. Reductions were calculated as log10 (effluent) minus log10 (influent). Microbial loads were calculated as numbers per inhabitant and day. To achieve this, the respective microbiological concentrations were multiplied by the amount of discharge and divided by the number of connected people. Visual and statistical data were analysed with Visplore 2.0 (Piringer et al., 2010) (VRVis GmbH, Austria, Vienna) and Sigma Plot 11.0 (SPSS Inc., Chicago, USA). To account for multiple testing, statistical significance levels were corrected according to Bonferroni. All graphs were prepared using Sigma Plot 13.0, Visplore 2.0 and CorelDraw X5 (Corel, Canada). To support correct comparisons of the variability of log-normal distributed variables, the multiplicative standard deviation s* was calculated for the recovered results according to Limpert et al. (2001). The multiplicative standard deviations s* ranged from 1.5 to 6.4 and from 1.5 to 1.7 for microbiological and chemical parameters, respectively. A s* > 1 is considered indicative to apply multiplicative standard deviation statistics (further details see Limpert et al., 2001).

Results

Characteristics of selected municipal WWTPs and chemical wastewater quality

Five municipal activated sludge WWTPs (WWTP2-6) in the metropolitan area of Vienna, Austria, with design capacities ranging from 20000 to 140000 p.e. and average loadings of 6600–78400 p.e. were selected. The number of connected inhabitants ranged from approx. 2000 up to approx. 31000 (Table 1). The catchments could be described as a mix of rural and urbanized areas. Sewers were constructed as combined systems, with a pressure pipe as inflowing sewer in case of WWTP4. Industrial influence in the catchment was evaluated as low to moderate. A potential impact on raw wastewater quality due to seasonal events relating to crop harvest and processing (i.e., wine production in the catchments of WWTPs 2, 4, 5, 6) and tourism (i.e., summer tourism in the catchment of WWTP5) could not be excluded. The average flow at the influent varied from 1600 m3 per day at WWTP5 to 18700 m3 per day at WWTP4 (Table 1).

Table 1

WWTP: wastewater treatment plant; PE: population equivalent; M: primary treatment: mechanical treatment step; C secondary treatment: biological carbon (C) removal; N,D,P: tertiary treatment: nutrient removal including nitrification (N): denitrification (D) and phosphorus removal (P); COD: chemical oxygen demand; TP: total phosphorus; TN: total nitrogen; NO3–N: nitrate nitrogen; NH4–N: ammonium nitrogen; infl.: influent; effl.: effluent of WWTP; * Annual mean COD load (kg a−1) divided by a COD load per person of 110 g COD d−1;** during summer.

WWTP	Design capacity [PE]	Actual average loading (PE)*	Inhabitants connected	Type	Advanced treatment/specifics	COD [mg L−1] median		TP [mg L−1] median		TN [mg L−1] median		NO3–N [mg L−1] median		NH4–N [mg L−1] median
						infl	effl	infl	effl	infl	effl	effl		effl
2	40,000	49,000	23,500	M, C, N, P	overloaded	466	36	9.3	0.7	58	25	21		0.3
3	23,000	16,300	10,800	M, C, N, D, P	–	553	16	8.8	0.3	55	9	4		0.1
4	140,500	78,400	30,800	M, C, N, D, P	–	441	16	5	0.3	43	9	4		0.1
5	20,000	6,600	2,100	M, C, N, D, P	UV disinfection**	443	16	5.5	0.2	38	2	0.4		0.2
6	45,000	20,000	15,000	M, C, N, D, P		555	24	8.5	0.7	63	14	11		0.4
7	21,000	17,200	5,900	M, C, N, D, P	membrane system	905	19	9	0.5	47	5	3		0.1

All WWTPs were using an activated sludge process with mechanical treatment and carbon removal. WWTPs 3–6 also performed nitrification and denitrification, while WWTP2 was overloaded and not designed for denitrification, resulting in lower nitrogen removal rates compared to other WWTPs (Table 1). The average sludge age ranged from 8 to 57 days (Table S1, additional material). P removal was achieved at all WWTPs by chemical precipitation with iron and/or aluminium salts, which is required for sensitive areas in the European Union (EC, 1991). WWTP5 was also equipped with additional UV-disinfection (UV-2000, Trojan Technologies, Canada; 48 UV lamps, operational parameters: UV Transmission UVT 10 mm: 65%, max discharge: 135 m3 h−1) at the effluent, which was operated only during summer. One additional activated sludge plant (WWTP7) was selected in the rural area of Bavaria, Germany. WWTP7 is a membrane bioreactor with a system of ultra-filtration membranes submerged in the aeration tank. The three vacuum rotation membrane units, each having a membrane surface of 2264 m2 and a pore size of approximately 38 nm, are used to separate the activated sludge flocs from the treated wastewater by means of a pressure difference. (physical solid-liquid separation process, Table 1).

In raw wastewater of WWTPs 2–7, median COD, TN, and TP showed values ranging from 440 to 900 mg L−1, from 40 to 60 mg L−1, and from 5 to 9 mg L−1, respectively (Table 1). Observed elimination rates were 94–98%, 77–93% (except WWTP2 with ~50%), and 90–96% for COD, TN, and TP, respectively, (Table S1, additional material). The median water temperature at the wastewater effluents was 12.5 °C–15.5 °C (further details in Table 1 and Table S1).

Characteristics of selected domestic WWTPs and chemical wastewater quality

Eight small dWWTPs (numbers 8–15) with a design capacity of 6–130 p.e., were selected in the metropolitan area of Vienna, Austria (Table S1, additional material). Two types of dWWTPs were encountered, including Dr. Renner® technology (also known as Gallé wastewater technique®) and Putox® technology (also known as Purator®). Two dWWTPs were linked to little taverns in the mountainous area of Vienna. Up to 50 persons were contributing their excreta to these sewer systems. Faecal load were strongly fluctuating. One dWWTP was localized at a horse barn with approximately 25 persons as permanent faecal sources. The rest of the dWWTPs was connected to individual households, with up to 5 contributing persons. Determined COD in the treated wastewater of dWWTPs 8–14 ranged from 15 to 56 mg L−1; the settle-able solids were generally less than 0.1 mg L−1. Ammonium and pH yielded values of 0.2–1.7 mg L−1 and 6.9–7.6, respectively, in the treated wastewater (see Table S1 additional material for more details).

Does wastewater from different municipal WWTPs show differences in GeBaM and SFIB concentrations?

One of the aims was to evaluate whether human-associated genetic Bacteroidetes faecal markers (GeBaM) and SFIB concentrations in raw and treated wastewater show significant differences with regard to the investigated municipal disposal systems or background conditions. The GeBaM BacHUM, HF183 TaqMan and AllBac as well as the cultivation-based parameter E. coli, enterococci and C. perfringens spores were selected as test parameters (for details see Table S3, additional material, n = 6–10 per WWTP). Statistical comparisons were performed between all the individually investigated municipal WWTPs, covering all WWTP2 to WWTP6 combinations (type 1 comparisons). Comparisons were also performed for the concentrations of microorganisms in wastewater from cool vs. warm seasons for the lumped results from WWTP 2–6 (type 2 comparisons). The results of type 1 and type 2 comparisons revealed no significant differences (Mann-Whitney Rank Sum Test, p < 0.05, Bonferroni-corrected). Hence, the results from WWTP2 to WWTP6 were pooled for further analysis (section 3.4). The results from WWTP7 (membrane filtration) and WWTP5 (UV-disinfection) are shown separately in sections 3.4.4. and 3.6.

Occurrence of microbial indicators in raw and biological treated wastewater from municipal systems

Prevalence and abundance of GeBaM and comparison with SFIB and human viral faecal markers

All investigated GeBaM markers (n = 37–45 per parameter and influent or effluent samples) showed 100% occurrence in raw and biologically treated wastewater with primary, secondary and tertiary treatments. The human-associated faecal markers BacHUM and HF183 TaqMan revealed remarkably similar concentrations, with medians of log10 8.6, log10 8.6 (raw) and log10 6.5, log10 6.4 (treated) ME per 100 ml wastewater, respectively (Fig. 1). The AllBac maker showed concentrations one order of magnitude higher, with medians of log10 9.6 (raw) and log10 7.8 (treated) ME per 100 ml wastewater.

Fig. 1

Concentration of standard faecal indicators and genetic microbial source tracking markers in raw (in) and treated (out) wastewater in pooled data from WWTP 2-6. Data shown are a pooled set. AllBac: genetic faecal marker for the total Bacteroidetes populations; BacHum, HF183: genetic faecal marker for human-associated Bacteroidetes populations; C. perf: Clostridium perfringens spores, Ent: enterococci, BtioPh: bacteriophages infecting Bacteroides thetaiotaomicron, HAdV: human adenovirus, JCPyV: JC polyomavirus, nq: number of quantifiable samples, s*: estimated multiplicative standard deviation, PR (%): Prevalence of investigated markers. Boxes cover the 25th to 75th percentile, line within the boxes, median; whiskers the 10th to 90th percentile.

SFIB (n = 47–49 per parameter and influent or effluent samples) also proved 100% prevalent in the investigated raw and treated wastewater, but with concentrations 2 to 3 orders of magnitude lower compared to GeBaM. The median concentrations for E. coli, Enterococci, and C. perfringens spores were log10 6.6, log10 6.1, log10 4.8 (raw) and log10 4.2, log10 3.8, and log10 3.4 (treated) CFU per 100 ml of wastewater, respectively (Fig. 1). The analysed bacteriophages and human viruses could not be detected in all samples. The prevalence of the human faecal-associated bacteriophage BtioPh (n = 24, 35 influent or effluent samples) was 97% in raw and treated wastewater. Prevalence rates for the human viruses HAdV and JCPyV (n = 24–47 per parameter and influent or effluent samples) were 92%, 98% (raw) and 85%, 60% (treated), respectively. Median concentrations for BtioPh, HAdV, and JCPyV were log10 4.5, log10 3.9, log10 4.2 (raw) and log10 2.5, log10 2.5, and log10 2.1 (treated) PFU or GC per 100 ml of wastewater (Fig. 1).

Variability of GeBaM concentrations in wastewater and comparisons with SFIB and human viral markers

In raw wastewater, the s* were remarkably low for the genetic marker concentrations of the BacHUM and HF183 TaqMan assay and were comparable with the variability of E. coli and enterococci, ranging from s* = 2.3 to s* = 3.0 (Fig. 1). C. perfringens spores proved to be the most constantly occurring indicator in raw wastewater (s* = 1.5).

A general increase in the variability of indicator concentrations between influent and effluent samples, irrespective of the considered parameter and WWTP, was obvious (p < 0.05, n = 9, Kruskal-Wallis). The multiplicative standard deviation s* increased by an average factor of 1.5 (range 0.9–2.0) during wastewater treatment (Fig. 1). In treated wastewater the variability of GeBaM concentrations matched the variability of SFIB concentrations as well (s* = 3.6 to 5.5).

The concentration variability of viral markers was higher, ranging from s* = 3.8 to 6.4 in raw wastewater and s* = 4.4 to 6.4 in treated wastewater. The investigated chemical parameters showed a statistical variability in the range of C. perfringens spores in raw wastewater (s* = 1.5, 1.8, 1.6 and 1.7 for COD, BOD, TN and TP, respectively). In treated wastewater the chemical parameters revealed lower variability compared to the microbiological parameter (s* = 1.6, 2.4, 3.1, 2.4 for COD, BOD, TN, and TP, respectively).

Establishing GeBaM loads per connected person and day

Medians for the calculated faecal marker loads AllBac, BacHUM, and HF183 TaqMan resulted in log10 13.2, log10 12.2, and log10 12.2 (raw wastewater) and log10 11.5, log1010.1, log10 9.9 (treated wastewater) ME per connected persons and day, respectively (Fig. S2, supplementary material). Quantitative relationships and statistical variability between GeBaM, SFIB and human viral faecal marker loads were similar to the obtained relationships regarding the concentrations (details are shown in Fig. S2, supplementary material).

Achieved microbiological reductions by wastewater treatment

C. perfringens revealed significantly lower reductions compared with the other SFIB and GeBaM. HAdV achieved lower reductions compared with E. coli (Kruskall Wallis, p < 0.001) during wastewater treatment (including primary, secondary and tertiary treatment in WWTPs 2–6). Some basic trends were obvious. E. coli achieved the highest reduction and showed a 10th-percentile value of -log101.6. GeBaM, enterococci, and bacteriophage BtioPh revealed very similar 10th-percentile values that ranged from -log101.0 to -log101.2 (Fig. 3). The lowest 10th-percentile reduction values were achieved by the human-specific viral faecal indicator HAdV and the bacterial faecal indicator C. perfringens spores, at -log10 0.5 and -log10 0.9 reductions, respectively (Fig. 3).

Fig. 3

Concentration of standard faecal indicators and genetic microbial source tracking markers in pooled data of municipal (M) WWTPs 2-6 effluents versus pooled data from domestic (D) dWWTPs 8-14 effluents. AllBac: genetic faecal marker for the total Bacteroidetes populations; BacHum, HF183: genetic faecal markers for human-associated Bacteroidetes populations; C. perf: Clostridium perfringens spores, Ent: Enterococci, nq: number of quantifiable samples, s*: estimated multiplicative standard deviation, PR (%): Prevalence of investigated markers. Boxes cover the 25th to 75th percentile, line within the boxes, median; whiskers the 10th to 90th percentile.

The reduction of GeBaM at the activated sludge membrane bioreactor (WWTP7) revealed a median 2.8 to 3.6 orders of magnitude increase in treatment efficacy of GeBaM compared to the conventional activated sludge treatment plants WWTPs 2–6 (Fig. S3, supplementary material).

Elucidating the relationships amongst GeBaM and other microbial/chemical variables

Except for total suspended solids (SS), statistical analysis of the pooled data set from WWTPs 2–6, including information from raw and treated wastewater, resulted in significant relationships amongst all parameters (correlation coefficients rho = 0.73–0.95, p < 0.0045; Table 2). Due to the inhomogeneous distribution of the raw vs. the treated wastewater data, such correlation analysis with pooled data led to a statistically biased relationship. Separate analysis was thus performed for the raw and treated wastewater data. Viral data were not included in this detailed correlation analysis because the replicate number was considered too low. Results for viral data are thus only given exemplarily for the raw waster data set, where appropriate.

Table 2

Correlation table showing the Spearman rank coefficient for the pooled influent and effluent (in & out), the influent (in) and the effluent (out) data of WWTP 2-6. AllBac: genetic faecal marker for the total Bacteroidetes populations; BacHum, HF183: genetic faecal markers for human-associated Bacteroidetes populations; C. perf: Clostridium perfringens spores, Ent: Enterococci; COD: chemical oxygen demand; TN: total nitrogen; NH4–N: ammonium nitrogen; TP: total phosphorus; TSS: total suspended solids.

AllBac																* = significant	Bonfer roni P = 0,0045
0.85*		BacHum
0.25	0.62*
0.81*		0.95*		HF183												in & out
0.19	0.45	0.80*	0.84*													in	out
0.77*		0.83*		0.78*		E. coli
0.11	0.42*	0.06	0.59*	0.08	0.58*
0.76*		0.83*		0.78*		0.86*		Ent
0.24	0.17	0.14	0.49*	0.04	0.46	0.39	0.52*
0.73*		0.82*		0.77*		0.80*		0.82*		C.perf.
0.26	-0.06	0.41*	0.15	0.22	0.18	0.39	-0.01	0.39	0.22
0.75*		0.81*		0.76*		0.80*		0.82*		0.89*		COD
0.29	0.03	0.36	0.17	0.19	0.23	0.43*	-0.03	0.39	0.13	0.65*	0.47*
0.60*		0.71*		0.67*		0.63*		0.65*		0.77*		0.80*		TN
0.27	0.01	0.51*	0.14	0.47	0.01	0.38	0.07	0.46*	0.08	0.62*	0.36	0.62*	0.66*
0.74*		0.83*		0.79*		0.77*		0.82*		0.81*		0.85*		0.71*		NH4
0.29	-0.04	0.48*	0.22	0.32	0.21	0.01	0.15	0.43	0.20	0.60*	-0.04	0.49*	0.32	0.73*	0.37
0.73*		0.80*		0.77*		0.78*		0.78*		0.84*		0.90*		0.82*		0.87*		TP
0.10	0.04	0.44	0.03	0.37	0.04	0.28	0.04	0.24	0.07	0.49*	0.28	0.58*	0.59*	0.76*	0.83*	0.61*	0.42
0.06		0.23		0.24		-0.03		0.20		0.52*		0.34		0.22		-0.03		0.13		TSS
n.a	0.06	n.a	0.23	n.a	0.24	n.a	-0.03	n.a	0.20	n.a	0.52*	n.a	0.34	n.a	0.22	n.a	-0.03	n.a	0.13

Correlation analysis for the raw wastewater data from catchments WWTPs 2–6 indicated a tight relationship between the human-associated faecal markers BacHUM and the HF183 TaqMan (rho = 0.80, p < 0.0045). A significant relationship became also obvious between the concentrations of the faecal marker HF183 TaqMan and the human-specific JC polyomavirus (rho = 0.45, p < 0.0045). In sharp contrast, the AllBac marker did not show a discernible relationship with the human-associated GeBaM or SFIB (Table 2). Interesting but non-significant correlation coefficients were obtained amongst the SFIB (rho = 0.39, p > 0.0045). Remarkably, a relationship between the human-associated GeBaM and the SFIB with the biological oxygen demand, the nitrogen content, and the phosphorus content in raw wastewater became evident. This was indicated by a range of significant correlations, including BacHUM, E. coli and enterococci with one or several components of the chemical parameters (rho = 0.46–0.51, p < 0.0045). Amongst the microbiological parameters, C. perfringens spores had the most pronounced relationship with the chemical quality characteristics of raw wastewater (rho = 0.49–0.65, p < 0.0045). Except for SS, a general interrelationship between all investigated chemical variables was obvious (rho = 0.49–0.76, p < 0.0045).

A different situation could be found for the data set of the treated wastewater. Except for C. perfringens spores, an increased tendency in associations amongst the microbiological variables during the treatment process was observed (i.e. 7 out of the 10 pairwise comparisons became significant after treatment, Table 2). In this respect, a correlation between the genetic AllBac and the human-associated BacHUM marker was discernible after the treatment (rho = 0.62, p < 0.0045). In contrast, relationships between the human-associated GeBaM and SFIB with the chemical quality characteristics of the wastewater disappeared during the treatment process (i.e. no significant correlation out of 25 pair-wise comparisons, Table 2). Again, the exception were C. perfringens spores, showing significant relationships with the COD and SS (Table 2).

Occurrence of GeBaM and SFIB in treated wastewater of small domestic WWTPs and comparison to municipal WWTPs

GeBaM in treated wastewater from dWWTPs showed 100% prevalence. The human-associated faecal markers BacHUM and HF183 TaqMan revealed similar concentrations, with medians of log10 6.3 and log10 6.2 ME per 100 ml wastewater (Fig. 3). The AllBac maker showed concentrations two orders of magnitude higher, with a median of log10 8.0 ME per 100 ml wastewater. FIB markers in treated domestic wastewater also resulted in 100% prevalence. The medians for E. coli, enterococci, and C. perfringens spores were log10 3.9, log10 3.6 and log10 3.9 CFU per 100 ml treated wastewater, respectively (Fig. 3). JCPyV were detected in 3 of 6 dWWTPs (log10 2.0-log10 3.0 ME per 100 ml of wastewater) evaluated that treated wastewater from 6 to 130 PE.

GeBaM and SFIB in treated wastewater had very similar concentrations for both the small domestic and the municipal WWTPs (Fig. 3). No differences in concentration could be detected (Man Whitney, p < 0.008, n from WWTPs randomly adjusted to number of dWWTPs). Information on the variability of results is reflected by the boxplots and the multiplicative standard deviations statistics (s*) directly given at Fig. 3.

Observed reductions due to UV irradiation (254 nm) (WWTP5)

After UV irradiation, the observed reductions for enterococci, E. coli, and somatic coliphages, given as 5th-percentile values (i.e., only 5% of the values showed a lower reduction), were log10 3.4, log10 3.0, and log 2.7, respectively (Fig. S4). C. perfringens spores were only slightly inactivated, resulting in a 5th-percentile reduction of log10 0.69. In contrast, no statistically significant reduction of GeBaM was detectable (one-way ANOVA, p < 0.05).

Discussion

Human-associated genetic Bacteroidetes faecal markers (GeBaM) were consistently detected in high concentrations in the investigated samples from raw and biological treated wastewater. The size of the studied wastewater systems varied over 4 orders of magnitude, with populations ranging from as few as 3 individuals up to 49000 inhabitants connected. Statistical analysis also demonstrated that GeBaM concentrations did not reveal differences regarding the type of the wastewater system or the time of the season investigated. Our results thus provide strong empirical evidence of the ubiquitous and abundant occurrence of GeBaM in raw and biologically treated wastewater, regardless whether the wastewater is derived from single households, larger settlements, or towns. Information on the quantitative occurrence and dynamics of GeBaM in wastewater along the wastewater and sanitation pathway has been limiting so far. The few studies available have focussed on individual samples or single systems (Ervin et al., 2013; Silkie and Nelson, 2009; Srinivasan et al., 2011; Stapleton et al., 2009). To our knowledge, our results provide the first comprehensive information on the occurrence and dynamics of GeBaM in raw and biologically treated wastewater from several well-characterized wastewater systems and treatment plants. Characterisation was facilitated by standard physicochemical and chemical analysis of raw and treated wastewater.

The selected systems were predominantly influenced by wastewater from households. No signs of significant influence from agriculture or industry could be found, and chemical analysis did not show any deviations from quality characteristics as expected for raw wastewater of municipal or domestic origin (Gujer, 2002). Data on the discharge dynamics also indicated that large rain events did not happen during the seasonal sampling campaigns (Table S1). A relevant influence on the wastewater quality due to surface runoff in the catchment area, potentially leading to strong dilution effects or to the input of faecal material from non-human sources, was thus not expected.

This study further supports the fact that faecal pollution based on GeBaM qPCR quantification can be performed with at least equal precision compared with traditional ISO-based cultivation techniques (Stapleton et al., 2009). The determined concentrations of GeBaM and SFIB (E.coli and enterococci) indicated equal statistical variability in raw and treated wastewater. This finding is of special interest regarding the current evaluation of GeBaM as a potentially new means to complement routine water quality testing (Betancourt and Fujioka, 2006; McQuaig et al., 2012; Molina et al., 2014). The statistical variability of GeBaM and SFIB concentrations was lowest for raw wastewater and, interestingly, increased during biological wastewater treatment. It can be assumed, that the raw wastewater (apparently) underlies sufficient mixing in the sewer channel to balance potential differences of input concentrations from the connected households. Diurnal variations were accounted for by the volume-proportional 24 h automated sampling. Interestingly, C. perfringens showed a statistical variability in the range of the measured chemical parameters, which was far lower than that of the rest of the microbiological parameters. A very low variability of C. perfringens in water has been reported previously (Byamukama et al., 2005).

It has to be mentioned that the statistical comparison of variability was supported by the implementation of two methodical innovations. To obtain representative samples from the influent and effluent of WWTPs, an automated 24-h volume-proportional and cooled sample was taken (Mayer et al., 2015). In contrast to this study, most studies dealing with microbiological investigations rely on randomly collected grab samples. In addition, the multiplicative standard deviation s* was introduced to the field of pollution microbiology to obtain an appropriate measure of statistical variability for log-normal distributed parameters (Limpert et al., 2001).

Unlike chemical load calculations, load calculations for microbial source tracking markers are very rare (Wilkes et al., 2014, 2013). The presented results suggest the future use of human-associated GeBaM loads as a valuable metric to estimate the impact of municipal and domestic wastewater input into the environment. The established median loads for raw and biological treated wastewater of approximately 1012 and 1010 molecule equivalents of human-associated GeBaM per person and day, respectively, demonstrate the sensitivity of GeBaM as a general measure of faecal pollution from municipal wastewater. Assuming a defecation rate of 100 g–1000 g of faecal excrement per person and day (Cummings et al., 1992; Geldreich, 1978), the estimated median load for raw wastewater can be converted back to a concentration range of 109 to 1010 human-associated GeBaM per g of faeces. This estimated range of GeBaM concentration compares well with concentrations in human faeces measured by qPCR (Haugland et al., 2010; Kildare et al., 2007).

The recovered GeBaM concentrations were in good agreement with previously reported levels (Reischer et al., 2013; Silkie and Nelson, 2009; Stapleton et al., 2009). GeBaM had concentrations at least two orders of magnitude higher than SFIB (cf. Fig. 1). Given the reported occurrence of intestinal microbiota in intestinal systems and human faeces, the dominance of GeBaM over SFIB is well known and expected (Ley et al., 2008; Reischer et al., 2007). This quantitative dominance of GeBaM in raw and biological treated wastewater is of high practical importance, regarding the sensitivity of molecular faecal pollution detection in comparison with cultivation-based standards. SFIB enumeration in water requires only minimal processing efforts. Samples are either directly analysed (MPN procedures) or subjected to membrane filtration before cultivation is started (ISO, 2000, 2005). PCR analysis involves several additional manipulation steps, including nucleic acid extraction, purification, and partial analysis of the extracted volumes (Ervin et al., 2013). Molecular detection methods thus have to apply higher sampling volumes or have to focus on more abundant targets to achieve comparable sample limits of detection (SLOD). The highly abundant nature of GeBaM in wastewater supports equal to superior sensitivity in comparison to SFIB methods, without the need for largely increased sampling volumes. This fact is the basis for the generation of large comparative sampling sets to appropriately cover pollution dynamics in aquatic systems (Ervin et al., 2013; Reischer et al., 2008, 2011; Riedel et al., 2014).

A high statistical association between the concentrations of the two human-associated GeBaM assays, the HF183 Taqman and the BacHUM, was observed for raw and treated wastewater (cf. Table 2). A possible explanation for this tight relationship can be found in the nature of the targeted human-associated Bacteroidetes populations. The most widely used human-associated GeBaM assays (including the above-mentioned ones) still focus on the same or similar phylogenetic sequence targets originally described by Bernhard and Field (Bernhard and Field, 2000). Recent research indicates that the HF183 Taqman and the BacHUM assay target populations within the species of Bacteroides dorei (McLellan and Eren, 2014). Although these assays revealed quite different specificity and sensitivity characteristics in a recent multi-laboratory study (Layton et al., 2013), our data elucidate the redundant nature of the HF183 Taqman and BacHUM assays for the detection of human-associated faecal pollution along the pathway of wastewater disposal. However, the tight association between these two independently performed assays proved the analytical precision of the recovered results within our study.

Statistically significant correlations between concentrations of GeBaM and SFIB could not be detected, although significant associations between the concentrations of faecal indicators and chemical parameters became obvious. A remarkable exception for raw wastewater was the slight but significant relationship between the human-associated BacHUM and the human-specific JCPyV (rho = 0.45, p < 0.02), pointing to the human-associated faecal pollution indication capacity of these molecular targets. Additionally, C. perfringens showed a slight correlation with BacHUM (rho = 0.41, p < 0.02). This is in line with a recent study, where C. perfringens revealed to be associated with human and carnivorous fecal sources, rather than with faecal emissions form herbivorous animals (Vierheilig et al., 2013).

The process of biological wastewater treatment had an equalising effect on the relationship amongst the microbiological variables, whilst the correlation of microbiological to chemical parameters, with the exception of a few cases, totally disappeared. To our knowledge, such a shift of correlation between microbiological parameters during wastewater treatment has not been reported so far and needs further verification.

Although the prevalence and abundance of the human-specific viruses HAdV and JCPyV found in the Austrian WWTP are lower than the range previously reported (Bofill-Mas et al., 2006; Rusinol et al., 2014), the nonetheless high prevalence of HAdV (92%) and JCPyV (98%) in the raw wastewater of the investigated municipal disposal systems in all seasons, and the detection of JCPyV in small dWWTPs, supports the usefulness of these specific tools as markers to trace human faecal pollution from WWTPs. The low abundance and prevalence were probably due to differences in the applied protocol (for instance, 500 ml of effluent municipal and domestic wastewater were concentrated instead of the 10 L tested in other studies). JCPyV and HAdV have been described as highly stable in the environment and present in nearly 100% of raw wastewater samples with concentrations up to log10 7.0 ME 100 ml−1 (Bofill-Mas et al., 2013). JCPyV is a highly specific human marker excreted in urine, and the detection of HAdV has been recently described as particularly useful for the prediction of risk in bathing waters (Marion et al., 2014). These parameters seem to be suitable tools to complement GeBaM/SFIB-based surface water monitoring for selected sampling locations or situations when higher sampling volumes can be processed. HAdV and JCPyV have been successfully applied for the identification of the source of contamination in river catchments covering various geographical areas (Rusinol et al., 2014). These viruses can also support verification of MST results in situations when the specificity level from human-associated GeBaM is deemed insufficient (Reischer et al., 2013).

GeBaM, E. coli and enterococci revealed similar reductions rates in the representatively chosen municipal activated sludge WWTPs (Fig. 2). Only C. perfringens spores demonstrated a lower reduction, most likely due to its resistant nature (Vierheilig et al., 2013). These results clearly demonstrate that GeBaM emission from municipal WWTP can be expected at fairly constant concentrations, and treatment just eliminates 2 log orders of magnitude from raw wastewater. The activated sludge membrane bioreactor removed approximately 5.4 and 5.0 log orders (10th percentile) of magnitude of BacHUM and HF183 TaqMan from wastewater, which is in line with previously reported data on bacterial removal in a membrane bioreactor (van den Akker et al., 2014). Data on UV irradiation (254 nm) indicated no discernible effect on the PCR detectable concentrations of GeBaM in wastewater, which is in agreement with a recently published study (Chern et al., 2014). It should be mentioned that the effect of chlorination was not investigated. This type of disinfection is not applied in European WWTPs. This study did not focus on other important factors, such as mobility, persistence, or specificity of the investigated indicators and markers in the aquatic environment. which also have to be considered for the appropriate interpretation of monitoring results.

Fig. 2

Reduction of standard faecal indicators and genetic microbial source tracking markers in municipal WWTPs 2-6 during wastewater treatment (pooled data). AllBac: genetic faecal marker for the total Bacteroidetes populations; BacHum, HF183: genetic faecal marker for human-associated Bacteroidetes populations; Ent: enterococci; C. perf: Clostridium perfringens spores,; BtioPh: bacteriophages infecting Bacteroides thetaiotaomicron, HAdV: human adenovirus, JCPyV: JC polyomavirus. Boxes cover the 25th to 75th percentile, line within the boxes, median; whiskers the 10th to 90th percentile.

In our study the occurrence and the dynamics of GeBaM in point sources along the human wastewater pathway as expected in the Central European Region have been shown.

A comparison to bacterial, viral, and process parameters at treatment (including membrane bioreactor and UV irradiation) was performed and volume-proportional 24-h auto-sampling and multiplicative standard deviation statistics were introduced.

The results strongly support the application of human-associated GeBaM to complement faecal pollution monitoring programs in water resources based on E. coli and enterococci.

Human-associated GeBaM were consistently detected in high concentrations in raw and biologically treated wastewater samples, irrespective of plant size, type and time of the season.

Human-specific JC polyomaviruses and adenoviruses seem to be a valuable support if highly specific markers are needed for MST.

Appendix A. Supplementary data

Supplementary data related to this article can be found at http://dx.doi.org/10.1016/j.watres.2015.12.031.