Contamination Scenario Matters when Using Viral and Bacterial Human-Associated Genetic Markers as Indicators of a Health Risk in Untreated Sewage-Impacted Recreational Waters

Assessing human-source microbial contamination of stormwater in the U.S

Stormwater capture and use (SCU) projects have the potential to provide a significant portion of municipal water demand. However, uncertainty about the degree of microbial contamination in stormwater and the required treatment is a barrier for implementation of SCU projects. Stormwater runoff could become contaminated with human fecal matter in areas with deteriorating infrastructure where raw wastewater exfiltrates from sewer networks to stormwater collection networks, homeless encampments exist, or sanitary sewer overflows (SSOs) occur. Estimation of human fecal contamination can inform selection of stormwater treatment targets. This study investigates stormwater microbial contamination originating from human fecal matter using observed detections and concentrations of human microbial source tracking (MST) markers and potentially human-infectious pathogens (PHIPs). First, a systematic review complied measurements of human MST markers in wet and dry weather stormwater flows and influent wastewater. In addition, measurements of viral pathogens (e.g., adenoviruses, norovirus GI+GII, and enteroviruses) and protozoan pathogens (e.g., Giardia lamblia and Cryptosporidium parvum) in wet weather flows and influent wastewater were assessed. Human MST marker and PHIP data were statistically analyzed and applied to estimate a human fecal contamination analog (HFCA) which is an estimate of the amount of human fecal matter based on relative concentrations of microbial contaminants in stormwater compared to municipal wastewater. Human MST-based HFCAs in wet and dry weather flows ranged from <10-7.0 to 10-1.5 (median = 10-4.5) and 10-12 to 10-2.6 (median = 10-7.0), respectively. PHIP-based HFCAs in wet weather flows ranged from ∼10-8 to 10-0.14. Estimates of human MST-based HFCAs are more reliable than PHIP-based HFCAs because the current PHIP datasets are generally limited by the number of data points, percent detection, variability observed within the statistical distributions, and the geographical span of sampling locations. The use of human MST-based HFCAs is recommended to guide the selection of stormwater treatment process trains that are protective of public health based on the intended end use. Application of HFCA 10-1 (i.e., sewage dilution 10-1) remains a reasonable conservative estimate of human fecal contamination in stormwater to inform selection of pathogen log reduction targets based on the data presently available.

The effect of sewage source on HF183 risk-based threshold estimation for recreational water quality management

Host-associated fecal indicator measurements can be coupled with quantitative microbial risk assessment to develop risk-based thresholds for recreational use of potential sewage-contaminated waters. These assessments require information on the relative concentrations of indicators and pathogens in discharged sewage, typically based on data collected from wastewater treatment plant influent samples. However, most untreated sewage releases occur from within the collection system itself (i.e. compromised sewer laterals, compromised gravity and force mains, sanitary sewer overflows), where these relationships may differ. This study therefore analyzed the concentrations of a selected reference pathogen (norovirus) and fecal indicator (HF183) in sewage samples from upper and lower segments of gravity sewage collection systems, wastewater pumpstations, and the influent and effluent of treatment plants, to characterize variability in their relative concentrations. Norovirus detection rates were lower and more variable in upper collection system samples due to the smaller population represented; whereas, HF183 was routinely detected at all sites with higher concentrations in the collection system compared to treatment plant influent, resulting in variable comparative relationships across sample locations (types). Mean HF183:NoV ratios ranged from 1.0 × 105 for sewer lateral samples to 7 × 10° for force main samples. Results were used to develop risk-based thresholds for HF183 based on estimated recreational exposure to norovirus following a release from each potential sewage source, with higher thresholds for treatment facility influent compared to forced mains, or effluent. Consequently, this approach can allow for the rapid application of potential risk-based thresholds for recreational water quality applications based on different types of sewage discharge events.

Risk-based critical concentrations of enteric pathogens for recreational water criteria and recommended minimum sample volumes for routine water monitoring

Concerns are rising about the contamination of recreational waters from human and animal waste, along with associated risks to public health. However, existing guidelines for managing pathogens in these environments have not yet fully integrated risk-based pathogen-specific criteria, which, along with recent advancements in indicators and markers, are essential to improve the protection of public health. This study aimed to establish risk-based critical concentration benchmarks for significant enteric pathogens, i.e., norovirus, rotavirus, adenovirus, Cryptosporidium spp., Giardia lamblia, Campylobacter jejuni, Salmonella spp., and Escherichia coli O157:H7. Applying a 0.036 risk benchmark to both marine and freshwater environments, the study identified the lowest critical concentrations for children, who are the most susceptible group. Norovirus, C. jejuni, and Cryptosporidium presented lowest median critical concentrations for virus, bacteria, and protozoa, respectively: 0.74 GC, 1.73 CFU, and 0.39 viable oocysts per 100 mL in freshwater for children. These values were then used to determine minimum sample volumes corresponding to different recovery rates for culture method, digital polymerase chain reaction and quantitative PCR methods. The results indicate that for children, norovirus required the largest sample volumes of freshwater and marine water (52.08 to 178.57 L, based on the 5th percentile with a 10 % recovery rate), reflecting its low critical concentration and high potential for causing illness. In contrast, adenovirus and rotavirus required significantly smaller volumes (approximately 0.24 to 1.33 L). C. jejuni and Cryptosporidium, which required the highest sampling volumes for bacteria and protozoa, needed 1.72 to 11.09 L and 4.17 to 25.51 L, respectively. Additionally, the presented risk-based framework could provide a model for establishing pathogen thresholds, potentially guiding the creation of extensive risk-based criteria for various pathogens in recreational waters, thus aiding public health authorities in decision-making, strengthening pathogen monitoring, and improving water quality testing accuracy for enhanced health protection.

Two risk assessments: Evaluating the use of indicator HF183 Bacteroides versus pathogen measurements for modelling recreational illness risks in an urban watershed

The purpose of this study was to evaluate the performance of HF183 Bacteroides for estimating pathogen exposures during recreational water activities. We compared the use of Bacteroides-based exposure assessment to exposure assessment that relied on pathogen measurements. We considered two types of recreational water sites: those impacted by combined sewer overflows (CSOs) and those not impacted by CSOs. Samples from CSO-impacted and non-CSO-impacted urban creeks were analysed by quantitative polymerase chain reaction (qPCR) for HF183 Bacteroides and eight human gastrointestinal pathogens. Exposure assessment was conducted two ways for each type of site (CSO-impacted vs. non-CSO impacted): 1) by estimating pathogen concentrations from HF183 Bacteroides concentrations using published ratios of HF183 to pathogens in sewage and 2) by estimating pathogen concentrations from qPCR measurements. QMRA (quantitative microbial risk assessment) was then conducted for swimming, wading, and fishing exposures. Overall, mean risk estimates varied from 0.27 to 53 illnesses per 1,000 recreators depending on exposure assessment, site, activity, and norovirus dose-response model. HF183-based exposure assessment identified CSO-impacted sites as higher risk, and the recommended HF183 risk-based threshold of 525 genomic copies per 100 mL was generally protective of public health at the CSO-impacted sites but was not as protective at the non-CSO-impacted sites. In the context of our urban watershed, HF183-based exposure assessment over- and under-estimated risk relative to exposure assessment based on pathogen measurements, and the etiology of predicted pathogen-specific illnesses differed significantly. Across all sites, the HF183 model overestimated risk for norovirus, adenovirus, and Campylobacter jejuni, and it underestimated risk for E. coli and Cryptosporidium. To our knowledge, this study is the first to directly compare health risk estimates using HF183 and empirical pathogen measurements from the same waterways. Our work highlights the importance of site-specific hazard identification and exposure assessment to decide whether HF183 is applicable for monitoring risk.

A multi-marker assessment of sewage contamination in streams using human-associated indicator bacteria, human-specific viruses, and pharmaceuticals

Human sewage contaminates waterways, delivering excess nutrients, pathogens, chemicals, and other toxic contaminants. Contaminants and various sewage indicators are measured to monitor and assess water quality, but these analytes vary in their representation of sewage contamination and the inferences about water quality they support. We measured the occurrence and concentration of multiple microbiological (n = 21) and chemical (n = 106) markers at two urban stream locations in Milwaukee, Wisconsin, USA over two years. Five-day composite water samples (n = 98) were collected biweekly, and sewage influent samples (n = 25) were collected monthly at a Milwaukee, WI water reclamation facility. We found the vast majority of markers were not sensitive enough to detect sewage contamination. To compare analytes for monitoring applications, five consistently detected human sewage indicators were used to evaluate temporal patterns of sewage contamination, including microbiological (pepper mild mottle virus, human Bacteroides, human Lachnospiraceae) and chemical (acetaminophen, metformin) markers. The proportion of human sewage in each stream was estimated using the mean influent concentration from the water reclamation facility and the mean concentration of all stream samples for each sewage indicator marker. Estimates of instream sewage pollution varied by marker, differing by up to two orders of magnitude, but four of the five sewage markers characterized Underwood Creek (mean proportions of human sewage ranged 0.0025 % - 0.075 %) as less polluted than Menomonee River (proportions ranged 0.013 % - 0.14 %) by an order of magnitude more. Chemical markers correlated with each other and yielded higher estimates of sewage pollution than microbial markers, which exhibited greater temporal variability. Transport, attenuation, and degradation processes can influence chemical and microbial markers differently and cause variation in human sewage estimates. Given the range of potential human and ecological health effects of human sewage contamination, robust characterization of sewage contamination that uses multiple lines of evidence supports monitoring and research applications.

Site-specific risk-based threshold (RBT) concentrations for sewage-associated markers in estuarine swimming waters

This study establishes site-specific risk-based threshold (RBT) concentrations for sewage-associated markers, including Bacteroides HF183 (HF183), Lachnospiraceae Lachno3 (Lachno3), cross-assembly phage (CrAssphage), and pepper mild mottle virus (PMMoV), utilizing quantitative microbial risk assessment (QMRA) for recreational estuarine waters (EW). The QMRA model calculates a RBT concentration corresponding to a selected target illness risk for ingestion of EW contaminated with untreated sewage. RBT concentrations were estimated considering site-specific decay rates and concentrations of markers and reference pathogen (human norovirus; HNoV), aiding in the identification of high-risk days during the swimming season. Results indicated varying RBT concentrations for fresh (Day 0) and aged (Days 1 to 10) sewage contamination scenarios over 10 days. HF183 exhibited the highest RBT concentration (26,600 gene copis (GC)/100 mL) initially but decreased rapidly with aging (2570 to 3120 GC/100 mL on Day 10) depending on the decay rates, while Lachno3 and CrAssphage remained relatively stable. PMMoV, despite lower initial RBT (3920 GC/100 mL), exhibited increased RBT (4700 to 6440 GC/100 mL) with aging due to its slower decay rate compared to HNoV. Sensitivity analysis revealed HNoV concentrations as the most influential parameter. Comparison of marker concentrations in estuarine locations with RBT concentrations showed instances of marker exceedance, suggesting days of potential higher risks. The observed discrepancies between bacterial and viral marker concentrations in EW highlight the need for optimized sample concentration method and simultaneous measurement of multiple markers for enhanced risk predictions. Future research will explore the utility of multiple markers in risk management. Overall, this study contributes to better understanding human health risks in recreational waters, aiding regulators, and water quality managers in effective decision-making for risk prioritization and mitigation strategies.

Occurrence, transmission and risks assessment of pathogens in aquatic environments accessible to humans

Pathogens are ubiquitously detected in various natural and engineered water systems, posing potential threats to public health. However, it remains unclear which human-accessible waters are hotspots for pathogens, how pathogens transmit to these waters, and what level of health risk associated with pathogens in these environments. This review collaboratively focuses and summarizes the contamination levels of pathogens on the 5 water systems accessible to humans (natural water, drinking water, recreational water, wastewater, and reclaimed water). Then, we showcase the pathways, influencing factors and simulation models of pathogens transmission and survival. Further, we compare the health risk levels of various pathogens through Quantitative Microbial Risk Assessment (QMRA), and assess the limitations of water-associated QMRA application. Pathogen levels in wastewater are consistently higher than in other water systems, with no significant variation for Cryptosporidium spp. among five water systems. Hydraulic conditions primarily govern the transmission of pathogens into human-accessible waters, while environmental factors such as temperature impact pathogens survival. The median and mean values of computed public health risk levels posed by pathogens consistently surpass safety thresholds, particularly in the context of recreational waters. Despite the highest pathogens levels found in wastewater, the calculated health risk is significantly lower than in other water systems. Except pathogens concentration, variables like the exposure mode, extent, and frequency are also crucial factors influencing the public health risk in water systems. This review shares valuable insights to the more accurate assessment and comprehensive management of public health risk in human-accessible water environments.

Risk of Gastroenteritis from Swimming at a Wastewater-Impacted Tropical Beach Varies across Localized Scales

Population growth and changing climate are expected to increase human exposure to pathogens in tropical coastal waters. We examined microbiological water quality in three rivers within 2.3 km of each other that impact a Costa Rican beach and in the ocean outside their plumes during the rainy and dry seasons. We performed quantitative microbial risk assessment (QMRA) to predict the risk of gastroenteritis associated with swimming and the amount of pathogen reduction needed to achieve safe conditions. Recreational water quality criteria based on enterococci were exceeded in >90% of river samples but in only 13% of ocean samples. Multivariate analysis grouped microbial observations by subwatershed and season in river samples but only by subwatershed in the ocean. The modeled median risk from all pathogens in river samples was between 0.345 and 0.577, 10-fold above the U.S. Environmental Protection Agency (U.S. EPA) benchmark of 0.036 (36 illnesses/1,000 swimmers). Norovirus genogroup I (NoVGI) contributed most to risk, but adenoviruses raised risk above the threshold in the two most urban subwatersheds. The risk was greater in the dry compared to the rainy season, due largely to the greater frequency of NoVGI detection (100% versus 41%). Viral log10 reduction needed to ensure safe swimming conditions varied by subwatershed and season and was greatest in the dry season (3.8 to 4.1 dry; 2.7 to 3.2 rainy). QMRA that accounts for seasonal and local variability of water quality contributes to understanding the complex influences of hydrology, land use, and environment on human health risk in tropical coastal areas and can contribute to improved beach management. IMPORTANCE This holistic investigation of sanitary water quality at a Costa Rican beach assessed microbial source tracking (MST) marker genes, pathogens, and indicators of sewage. Such studies are still rare in tropical climates. Quantitative microbial risk assessment (QMRA) found that rivers impacting the beach consistently exceeded the U.S. EPA risk threshold for gastroenteritis of 36/1,000 swimmers. The study improves upon many QMRA studies by measuring specific pathogens, rather than relying on surrogates (indicator organisms or MST markers) or estimating pathogen concentrations from the literature. By analyzing microbial levels and estimating the risk of gastrointestinal illness in each river, we were able to discern differences in pathogen levels and human health risks even though all rivers were highly polluted by wastewater and were located less than 2.5 km from one another. This variability on a localized scale has not, to our knowledge, previously been demonstrated.

Distribution of human fecal marker genes and their association with pathogenic viruses in untreated wastewater determined using quantitative PCR

Quantitative microbial risk assessment (QMRA) of human health risks using human fecal marker genes (HFMGs) is an useful water quality management tool. To inform accurate QMRA analysis, generation of probability distribution functions for HFMGs, and reference pathogenic viruses can be improved by input of correlation and ratios based upon measurement of HFMGs and gene copies (GC) of pathogenic viruses in untreated wastewater. The concentrations of four HFMGs (Bacteroides HF183, Lachnospiraceae Lachno3, CrAssphage and pepper mild mottle virus (PMMoV)), and GC of three reference pathogenic viruses human adenovirus 40/41 (HAdV 40/41), human norovirus GI + GII HNoV GI + GII and enterovirus (EV) were measured in untreated wastewater samples collected over a period of 12 months from two wastewater treatment plants in Sydney, Australia using quantitative polymerase chain reaction (qPCR) and reverse transcription qPCR (RT-qPCR). Over the course of the study, the GC of potential pathogenic viruses were 3-5 orders of magnitude lower than HFMGs in untreated wastewater. The GC of pathogenic viruses were highly variable over the course of the study, which contrasted with the concentrations of HFMGs that were quite stable with little variation observed within and between WWTPs. Among the HFMGs, HF183, CrAssphage and PMMoV correlated well with pathogenic virus GC, whereas weak or negative correlations were observed between Lachno3 and pathogenic virus GC. While the two assessed WWTPs had dissimilar population service sizes, the ratios between log₁₀ transformed pathogenic virus GC and HFMGs demonstrated similar central tendency and variability for the same combinations between WWTP A and WWTP B with no difference between the WWTPs. This suggests the widespread presence of these HFMGs in both populations serviced by these two WWTPs. The observed correlation and ratios of HFMGs and GC of reference pathogenic viruses can contribute to improved QMRA of human health risks in environmental waters subject to fresh sewer overflows.

Highly variable removal of pathogens, antibiotic resistance genes, conventional fecal indicators and human-associated fecal source markers in a pilot-scale stormwater biofilter operated under realistic stormflow conditions

Green stormwater infrastructure systems, such as biofilters, provide many water quality and other environmental benefits, but their ability to remove human pathogens and antibiotic resistance genes (ARGs) from stormwater runoff is not well documented. In this study, a field scale biofilter in Southern California (USA) was simultaneously evaluated for the breakthrough of a conservative tracer (bromide), conventional fecal indicators, bacterial and viral human-associated fecal source markers (HF183, crAssphage, and PMMoV), ARGs, and bacterial and viral pathogens. When challenged with a 50:50 mixture of untreated sewage and stormwater (to mimic highly contaminated storm flow) the biofilter significantly removed (p < 0.05) 14 of 17 microbial markers and ARGsin descending order of concentration reduction: ermB (2.5 log(base 10) reduction) > Salmonella (2.3) > adenovirus (1.9) > coliphage (1.5) > crAssphage (1.2) > E. coli (1.0) ∼ 16S rRNA genes (1.0) ∼ fecal coliform (1.0) ∼ intl1 (1.0) > Enterococcus (0.9) ∼ MRSA (0.9) ∼ sul1 (0.9) > PMMoV (0.7) > Entero1A (0.5). No significant removal was observed for GenBac3, Campylobacter, and HF183. From the bromide data, we infer that 0.5 log-units of attenuation can be attributed to the dilution of incoming stormwater with water stored in the biofilter; removal above this threshold is presumably associated with non-conservative processes, such as physicochemical filtration, die-off, and predation. Our study documents high variability (>100-fold) in the removal of different microbial contaminants and ARGs by a field-scale stormwater biofilter operated under transient flow and raises further questions about the utility of human-associated fecal source markers as surrogates for pathogen removal.

Modeling Untreated Wastewater Evolution and Swimmer Illness for Four Wastewater Infrastructure Scenarios in the San Diego-Tijuana (US/MX) Border Region

The popular beaches of the San Diego-Tijuana (US/MX) border region are often impacted by untreated wastewater sourced from Mexico-via the Tijuana River Estuary (TJRE) and San Antonio de los Buenos outfall at the Pt. Bandera (SAB/PTB) shoreline, leading to impacted beaches and human illness. The US-Mexico-Canada trade agreement will fund border infrastructure projects reducing untreated wastewater discharges. However, estimating project benefits such as reduced human illness and beach impacts is challenging. We develop a coupled hydrodynamic, norovirus (NoV) pathogen, and swimmer illness risk model with the wastewater sources for the year 2017. The model is used to evaluate the reduction in human illness and beach impacts under baseline conditions and three infrastructure diversion scenarios which (Scenario A) reduce SAB/PTB discharges and moderately reduce TJRE inflows or (Scenarios B, C) strongly reduce TJRE in inflows only. The model estimates shoreline untreated wastewater and NoV concentrations, and the number of NoV ill swimmers at Imperial Beach CA. In the Baseline, the percentage of swimmers becoming ill is 3.8% over 2017, increasing to 4.5% during the tourist season (Memorial to Labor Day) due to south-swell driven SAB/PTB plumes. Overall, Scenario A provides the largest reduction in ill swimmers and beach impacts for the tourist season and full year. The 2017 tourist season TJRE inflows were not representative of those in 2020, yet, Scenario A likely still provides the greatest benefit in other years. This methodology can be applied to other coastal regions with wastewater inputs.