Municipal Wastewater as a Microbial Surveillance ... - ACS Publications

Mar 29, 2018 - Department of Civil and Environmental Engineering, the University of Hawaii at Manoa, Honolulu, Hawaii 96822, United States. ‡ .... s...
0 downloads 6 Views 3MB Size
Article Cite This: Environ. Sci. Technol. XXXX, XXX, XXX−XXX

pubs.acs.org/est

Municipal Wastewater as a Microbial Surveillance Platform for Enteric Diseases: A Case Study for Salmonella and Salmonellosis T. Yan,*,† P. O’Brien,‡ J. M. Shelton,† A. C. Whelen,‡,§ and E. Pagaling†,# †

Department of Civil and Environmental Engineering, the University of Hawaii at Manoa, Honolulu, Hawaii 96822, United States State Laboratories Division, Hawaii Department of Health, Honolulu, Hawaii 96782, United States § Department of Microbiology, the University of Hawaii at Manoa, Honolulu, Hawaii 96822, United States ‡

S Supporting Information *

ABSTRACT: Municipal wastewater (MW) contains a conglomeration of human enteric microbiota from a community and, hence, represents a potential surveillance tool for gastrointestinal infectious disease burden at the community level. To evaluate this, the concentration of Salmonella in MW samples from Honolulu, Hawaii, was monitored over a 54-week period, which showed positive and significant linear and rank correlation with clinical salmonellosis case numbers over the same period. Salmonella isolates were obtained from the MW samples and then compared with clinical isolates obtained by the Hawaii Department of Health State Laboratories over the same period. The MW isolate collection contained 34 serotypes, and the clinical isolate collection contained 47 serotypes, 21 of which were shared between the two isolate collections, including nine of the 12 most commonly detected clinical serotypes. Most notably, nine Salmonella strains, including one outbreak-associated Paratyphi B strain and eight other clinically rare strains, were shared and concurrently detected between the MW and the clinical isolate collections, indicating the feasibility of using enteric pathogens in the MW as a timely indication of community enteric disease activity.



communities.16−18 We hypothesize that raw MW represents a unique source of enteric pathogen information for a community that is independent of the clinical data stream and potentially offers a source of time-sensitive disease information that is otherwise difficult to obtain or is completely inaccessible. In a previous study, we illustrated the use of MW to estimate the disease burden of enteropathogenic E. coli (EPEC) and the genetic diversity of eae gene in a community.19 This work generated otherwise inaccessible information on a nonreportable disease; however, the lack of clinical information on EPEC prevented cross validation of the community disease burden estimated by the MW-based approach. Salmonella enterica consists of a diverse group of serotypes that are ubiquitous in nature and is a common cause of human salmonellosis, which leads to considerable morbidity and mortality.20,21 For example, Salmonella enterica serotype Typhi (S. Typhi) and Paratyphi (S. Paratyphi) are the etiological agents of typhoid fever, which causes 17 to 21 million infections and more than 600 000 deaths worldwide each year.22 Although typhoid fever has been essentially eradicated in the U.S., a significant number of nontyphoidal salmonellosis cases (1.4 million per year) still occur23 and result in 15 000 hospitalizations and 400 deaths a year.21 Salmonellosis is a reportable

INTRODUCTION Microbial infectious diseases continue to be a major cause of human illness in the U.S. and throughout the world. The U.S. CDC estimated that, each year, enteric pathogens alone cause 47.8 million illnesses, resulting in 3037 deaths,1,2 and globally, over 13 million deaths were attributed to microbial infectious diseases.3 This problem is further exacerbated by emerging and resurgent microbial pathogens, including pathogens with heightened virulence and increasing antimicrobial resistance.4,5 To effectively combat this problem, it is critical to have comprehensive and timely surveillance capabilities for microbial infectious diseases. However, most of the existing disease surveillance approaches rely predominantly on clinical diagnostic data, which only cover a limited number of diseases and lack timeliness due to reporting biases and procedural delays.6 Other alternative or supplemental mechanisms, including syndromic surveillance7 and the emerging big data analytics,8 have yet to prove their ability to provide immediate and actionable public health value. Since municipal wastewater (MW) systems continuously collect fecal human waste from communities, MW is expected to contain a collection of human gastrointestinal microbiota. Previous studies have investigated the feasibility of using MW to monitor tobacco and illicit drug use9,10 and to detect contaminants of emerging concern in human communities.11 Other related studies have focused on detecting and/or inactivating pathogens12−15 and determining MW as the source of microbial pathogens in the environment and diseases in the © XXXX American Chemical Society

Received: January 10, 2018 Revised: March 29, 2018 Accepted: April 4, 2018

A

DOI: 10.1021/acs.est.8b00163 Environ. Sci. Technol. XXXX, XXX, XXX−XXX

Article

Environmental Science & Technology

random from the colonies on XLD plates and then examined on LIA and TSA slants for confirmation of Salmonella. Further confirmation was performed using a qPCR assay that targets the Salmonella-specific ttr locus.27 A total of 378 confirmed Salmonella isolates were obtained for the MW isolate collection and were stored in LB broth with 50% (v/v) glycerol at −80 °C for further analyses. The clinical Salmonella isolates were identified by community clinical laboratories and submitted to Hawaii State Laboratories where they were confirmed biochemically and serotyped for routine public health surveillance of salmonellosis.28 PFGE Analysis. The whole cell PFGE method and subsequent data analysis of the Salmonella isolates from the clinical and the MW isolate collections were conducted using the standard operating procedure from PulseNet PFGE.29 Briefly, the Salmonella isolates were grown on TSA with 5% defibrinated sheep blood (TSA-SB) and incubated overnight. Colonies were picked and resuspended in cell suspension buffer (100 mM Tris-HCl/100 mM EDTA, pH 8.0) to achieve a concentration of 0.4−0.45 on a Dade turbidity meter. 200 μL of this suspension was embedded in 0.5% Seakem Gold agarose (Lonza, Allendale, NJ). These PFGE plugs were then incubated in cell lysis buffer (50 mM Tris/50 mM EDTA, pH 8.0, and 1% Sarcosyl) with proteinase K (1 mg/mL) for at least 1.5 h. Restriction digestion of genomic DNA within the agarose plugs was conducted using Xbal (50 U/sample). Electrophoresis was conducted on a CHEF Mapper (Bio-Rad, Richmond, CA). Analysis was performed using BioNumerics version 5.10 (Applied Maths, Austin, TX), and serotype classification was based on clustering analysis with entries in the PulseNet database. Data Analysis. Two derivative parameters of CSalm were calculated, including total weekly Salmonella flux (JSalm = CSalm × Qweek) and the ratio between CSalm and CEnt in the same MW samples (RS/E). Log-transformed CSalm and CEnt were analyzed to identify outliers using the 1.5 interquartile range (IQR) method, which identified two log CSalm outliers (Weeks 4 and 38). After outlier removal, the log CSalm and log CEnt data were tested and confirmed to satisfy the normality assumption based on the Shapiro-Wilk test, and the Student’s t-test was conducted to determine whether statistically significant differences existed in the log CSalm and log CEnt values between the summer months and winter months. The clinical salmonellosis case number (N), which is typically directly used without log transformation in epidemiology studies,30 was also analyzed using the IQR method, which identified one outlier (Week 18). In preparation for computing Pearson’s product-moment correlation coefficients (r), log CSalm, log CEnt, log JSalm, log RS/E, and N underwent pairwise deletion of outliers and were subsequently also confirmed of normality based on the ShapiroWilk test. Pearson’s product-moment correlation coefficients (r) between N and log CSalm, log JSalm, log RS/E, and log CEnt were computed to detect underlying linear relationships, and the significance of the correlation coefficient was determined using the t distribution. The nonparametric Spearman’s rank correlation coefficient (ρ), which is commonly used in epidemiology studies for its robustness toward influential data points and outliers and does not require normality,31,32 was also computed using all data (including outliers). The relative abundance of each Salmonella serotype was calculated as the ratio between the number of isolates classified to the serotype and the total number of isolates in the respective isolate collection. The “Undertermined” serotype is the collection of

disease, and positive clinical specimens or isolates are sent to public health laboratories for characterization. Information on clinical cases of salmonellosis is available through an extensive genetic database of clinical Salmonella isolates called PulseNet,24 which permits validation of the proposed MWbased approach in comparison with the traditional clinical approach. In this study, a 54-week MW sampling campaign was conducted to collect composite raw MW samples each week at the Sand Island Wastewater Treatment Plant (SIWTP) in Honolulu, Hawaii. Salmonella concentrations in the MW samples were determined and subsequently compared with the reported clinical salmonellosis case numbers through correlation analysis. Salmonella isolates were also obtained from the MW samples, and their serotypes were determined by pulsed-field gel electrophoresis (PFGE) and compared with that of the clinical isolate collection. Finally, identical Salmonella PFGE patterns between the MW isolate collection and the clinical isolate collection were identified, and their temporal detection patterns were compared to reveal concurrent detection in both collections.



MATERIALS AND METHODS MW Sample and Flow Data Collection. MW sampling was conducted at the Sand Island Wastewater Treatment Plant (SIWTP) in Honolulu, Hawaii, over a 54-week period (April 27, 2010 to May 9, 2011). The SIWTP collects sanitary MW from the Honolulu urban center and surrounding areas where major employment centers are located and treats approximately 60% of the MW of the City and County of Honolulu. Each Sunday, raw untreated MW samples (1 L) were collected hourly using an autosampler placed upstream of the primary clarifiers, and 40 mL of homogenized hourly samples was mixed to make the daily composite samples. The samples were stored at 4 °C in the dark and transported back to the laboratory for processing within 24 h. The MW flow data were collected by the SIWTP using a Doppler flow meter and were used to calculate the total weekly MW flow (Qweek). Bacterial Enumeration in MW Samples. A five-tube most probable number (MPN) method modified from the EPA method 1682 for biosolids25 was used to enumerate Salmonella in the MW samples. Briefly, 10-fold serial dilutions (up to 10−4) of the raw MW samples were used to inoculate tryptic soy broth (TSB). After incubation at 35 °C for 24 h, 30 μL of each MPN tube was inoculated onto the modified semisolid Rappaport Vassiliadis (MSRV) media and incubated at room temperature for 1 h and then at 42 °C for 16−18 h. Potential Salmonella-positive tubes were indicated by a ring of motility around the inoculation site on the MSRV plates. Samples from the motility zone were then taken with a sterile loop and streaked onto xylose lysine deoxycholate (XLD) agar and incubated at 35 °C for 24 h. Colonies on the XLD agar that were red, pink, or colorless with black centers were then stab inoculated onto lysine iron agar (LIA) and tryptic soy agar (TSA) slants to confirm that they were Salmonella. The Salmonella concentrations in the MW samples (CSalm) were calculated using the five-tube MPN table and reported as MPN/100 mL. The concentration of a commonly used fecal indicator bacteria enterococci was also determined in the MW samples (CEnt) using the standard membrane filtration mEI agar method.26 Salmonella Isolation. From each of the 49 Salmonellapositive MW samples, five to ten isolates were picked at B

DOI: 10.1021/acs.est.8b00163 Environ. Sci. Technol. XXXX, XXX, XXX−XXX

Article

Environmental Science & Technology

Figure 1. Concentration dynamics of Salmonella (CSalm) and enterococci (CEnt) in daily composite MW samples in Honolulu over 54 weeks (from April 26, 2010 to May 8, 2011). The inset is a box plot showing distribution of the concentration data.

of salmonellosis disease burden, as typically more salmonellosis disease cases are reported in summer months than in winter months in Honolulu. Similar seasonality of salmonellosis is commonly observed in temperate regions of the U.S.,37 which is attributed to higher ambient temperatures38 and intensified precipitation39 in summer months. Correlation between MW Salmonella and Clinical Salmonellosis Cases. To further determine the relationship between Salmonella in MW and salmonellosis prevalence in the community, the CSalm in MW samples was compared with clinical salmonellosis case numbers (N) of the same weeks through correlation analyses (Table 1). Since the CSalm in the

all isolates that were not assigned to any particular serotype. All statistical analyses were conducted in R, and the default significance value cutoff is P ≤ 0.05 unless stated otherwise.



RESULTS AND DISCUSSION Salmonella Concentration in Raw MW Samples. The concentrations of Salmonella (CSalm) and enterococci (CEnt), a common fecal indicator bacteria used here as a reference, were determined in 54 raw MW samples from SIWTP collected over a 54-week sampling period (Figure 1). Salmonella was detected in 49 of the 54 raw MW samples, and the CSalm fluctuated significantly over the sampling period, as indicated by a wide concentration range (4−105.2MPN/100 mL) and a large coefficient of variation (5.54). The lower quartile, median, and upper quartiles of CSalm were 18, 46, and 256 MPN/100 mL, respectively, with an overall geometric mean of 89 MPN/ 100 mL. This is in contrast with the relatively stable CEnt in the same raw MW samples, which showed a much narrower concentration range (105.5−106.8CFU/100 mL) and a much smaller coefficient of variation (0.74). The large fluctuation of CSalm (compared to the relative stability of CEnt) suggests considerable variation in the number of Salmonella cells excreted into the MW system and potentially fluctuating salmonellosis disease prevalence over time. The CSalm values detected in the raw MW samples in Honolulu were in agreement with some previous estimations of Salmonella in MW (e.g., 102−104 cells/100 mL33) and previous reports from different countries, including 1.5−27 CFU/100 mL in French MW,34 93−104.0 MPN/100 mL in Finnish MW,35 and 4.5 × 105−2.4 × 106 MPN/100 mL in Mexican MW.36 The CSalm values in the summer months (May−Oct., Weeks 1−27; geometric mean: 102.4 MPN/100 mL) were significantly higher (t test, P < 0.001) than those in the winter months (Nov.−Apr., Weeks 28−54; geometric mean: 101.5 MPN/100 mL). In contrast, the CEnt values in the summer months (geometric mean: 106.0 CFU/100 mL) showed no significant difference (t test, P = 0.99) from those in the winter months (geometric mean: 106.0 CFU/100 mL). This seasonality of CSalm in the MW samples corresponds to the seasonality

Table 1. Correlation between Log-Transformed Salmonella Concentration (log CSalm), Salmonella Weekly Flux (log JSalm), Enterococci Concentration (log CEnt), Salmonella/ Enterococci Concentration Ratio (log RS/E), and Clinical Salmonellosis Case Number (N) over the Sampling Period clinical salmonellosis case number (N) Pearson’s r (n = 46) log log log log

CSalm JSalm RS/E CEnt

0.397 0.392 0.386 −0.109

(P (P (P (P

= = = =

0.006) 0.007) 0.008) 0.471)

Spearman’s ρ (n = 49) 0.329 0.325 0.326 −0.208

(P (P (P (P

= = = =

0.022) 0.023) 0.023) 0.150)

MW samples can also be affected by wastewater flow, which can be significantly impacted by fluctuation in water consumption and rainfall-derived infiltration and inflow (RDII),40 the weekly flux of Salmonella (JSalm) and the concentration ratio between Salmonella and enterococci (RS/E) were also calculated and compared with N (Figure S1). Pearson’s correlation coefficients showed that the log transformed CSalm, JSalm, and RS/E of the MW samples all exhibited positive and statistically significant correlation with N (r: 0.386−0.397, P: 0.006−0.008). In comparison, the CEnt of the MW samples did not exhibit either positive or significant Pearson’s correlation with N (r = −0.109, P = 0.471). The Spearman’s rank correlation coefficients, which were computed using all data including outliers, also showed C

DOI: 10.1021/acs.est.8b00163 Environ. Sci. Technol. XXXX, XXX, XXX−XXX

Article

Environmental Science & Technology both positive and significant correlation between log CSalm, log JSalm, and log RS/E of the MW samples and N. Again, log CEnt showed neither positive nor significant Spearman’s rank correlation coefficient with N. The presence of positive and significant correlation between log CSalm and the clinical salmonellosis disease case number N demonstrates the direct association between Salmonella abundance in the MW samples and salmonellosis disease prevalence in the community. The two derivative parameters of log CSalm (i.e., log JSalm and log RS/E) showed similar correlation with N as log CSalm, suggesting that water consumption changes, and RDII did not significantly impact the correlation. Indeed, an extreme wet peaking factor of 4 would only change log CSalm by a factor of 0.6, which is dwarfed by the large variation in the amount of Salmonella cells excreted by infected individuals (estimated to be 2.5 × 105 to 1 × 109 Salmonella cells per gram fecal material41). It should be noted that the removal of two outlier data points in log CSalm and one outlier in N was necessary in order to meet the normality assumption for Pearson’s correlation analysis. These outliers were very influential in computing Pearson’s r, as their inclusion in the analysis would have increased the P value of correlation coefficients and would no longer have made the correlation significant (data not shown). The outlier data in log CSalm and N may be attributed to measurement errors and several other factors. First, the concentration of Salmonella in patients’ feces can vary widely and cause significant fluctuation in the amount of cells discharged into the MW system by the same number of patients.41 Second, besides fecal waste from infected humans, other sources may also contribute Salmonella to the MW, including kitchen waste from contaminated foods1 and house pets.42,43 Third, the Salmonella discharged into the MW collection systems may undergo growth and decay both in the bulk MW and within biofilms on the sewer walls, further complicating Salmonella concentration dynamics in the MW. Finally, the salmonellosis case number N can be significantly impacted by patients’ reluctance to seek medical assistance and varying clinical practices of disease diagnosis; it was estimated that only about 1/30 of individuals with enteric diseases seek medical assistance while health practitioners can also vary in lab test ordering practices.6,7 MW Salmonella Serotypes. To understand the diversity of Salmonella in the MW samples, the 378 Salmonella isolates from the MW isolate collection were characterized by PFGE. They were grouped into 34 different serotypes, which exhibited significant variation in their relative abundance in the MW isolate collection (Figure 2). S. Derby dominated the MW collection, making up 20.9% of the relative abundance, which is considerably higher than the second most abundant serotype, S. Paratyphi B var. L(+) tartrate(+) (7.7%). Overall, the 15 major serotypes (defined as ≥2.0% relative abundance), which included many of the most frequently reported clinical serotypes such as S. Enteritidis, S. Typhimurium, S. Newport, S. Muenchen, and S. Infantis,30 collectively accounted for 85.4% of the MW isolate collection. The remaining 19 minor serotypes in the MW collection accounted for the remaining 14.6% and were generally not frequently detected in health clinics. A similar distribution of Salmonella serotypes in MW samples was also reported previously in Zaragoza, Spain.44 The serotype distribution of the MW isolate collection was compared with that of the clinical isolate collection of the same period. The clinical isolate collection included 338 Salmonella isolates from clinical cases reported in Hawaii, which were

Figure 2. Relative abundance of Salmonella isolates grouped into different serotypes in the MW and clinical isolate collections. The serotypes are ranked on the basis of their relative abundances in the MW isolate collection.

grouped into 47 serotypes. The MW isolate collection shared 21 serotypes in common with the clinical isolation collection, which included seven of the 11 clinical major serotypes (defined as ≥2.0% relative abundance). Although the MW collection provided good coverage of the major serotypes contained in the clinical isolate collection, there remained 26 serotypes in the clinical isolate collection that were absent in the MW collection. The majority of these serotypes (76.9%) were represented by only one (57.7%) or two (19.2%) clinical isolates. Therefore, the absence of these serotypes in the MW isolate collection could have resulted from the small numbers of Salmonella cells discharged into the MW by a few infected individuals. The limited sampling depth (≤10 isolates per sample) used in constructing the MW isolate collection may also explain this observation. The MW isolate collection also contained many serotypes that were not detected in the clinical isolate collection, including four major serotypes S. Mbandaka (6.3% relative abundance), S. Senftenberg (5.8%), S. Anatum (4.5%), and S. London (2.9%). Among them, only S. Anatum appeared in the top 20 most frequently reported clinical serotypes according to the U.S. CDC Salmonella reports.30,37 The presence of these serotypes in MW but absence in the clinical data could be explained by variable virulence, subclinical or asymptomatic shedding,45 and the presence of other environmental sources (for example, S. Senftenberg is common in D

DOI: 10.1021/acs.est.8b00163 Environ. Sci. Technol. XXXX, XXX, XXX−XXX

Article

Environmental Science & Technology marine sediment and water46). Alternatively, it is also possible that infections caused by these serotypes were not detected by the traditional clinical approach for disease surveillance, which would further indicate the potential benefits of including MW monitoring as part of an integrated surveillance strategy. Concurrent Detection of an Outbreak-Causing Salmonella Strain. To determine whether timely disease information is contained in the MW, PFGE patterns of the Salmonella isolates were compared. Twenty-four different PFGE patterns, which henceforth are referred to as 24 Salmonella strains grouped into 13 different serotypes (Table 2), were shared between the MW and the clinical isolate Table 2. Total Numbers of Unique Salmonella PFGE Patterns (Strains) of Different Serotypes That Were Shared between the Clinical and MW Collections and the Number of Those Showing Concurrent Detection Salmonella serotypes

shared PFGE patterns

concurrent detection

Typhimurium Enteritidis Muenchen Infantis Newport Paratyphi B Derby Agona Amager Uganda Javiana Chester Takoradi total

3 2 4 3 1 1 3 1 1 2 1 1 1 24

0 0 0 0 1 1 1 0 1 2 1 1 1 9

Figure 3. PFGE patterns of clinical (labeled N) and MW (labeled Y) Salmonella Paratyphi B isolates (A) and their temporal detection patterns in the MW and clinical isolate collections over the 54-week sampling period (B).

outbreaks47 and the detection of Cryptosporidium in MW during a waterborne outbreak in the UK.48 A previous study in Sweden also reported that Salmonella strains in sewage sludge were identical to clinical strains recovered from an outbreak that had occurred months to years beforehand.49 To the best of the authors’ knowledge, the present study is the first report where a bacterial pathogen strain was detected in the MW during its outbreak in the community. Later in the study period, there were three sporadic clinical cases caused by the same S. Paratyphi B strain, though no corresponding matches were detected in the MW collection, which again was likely to be caused by insufficient sampling depth used when establishing the MW collection. From week 47 to week 54, a total of 12 isolates belonging to S. Paratyphi B were detected in the MW collection, and in week 53, S. Paratyphi B isolates were dominant in the MW sample. However, no clinical isolates of S. Paratyphi B were detected during that period, suggesting possible low sensitivity (falsenegatives) detection of salmonellosis disease burden in the community by the traditional clinical approach. It is also possible that the lack of correspondence between the MW and the clinical isolate collections was due to the S. Paratyphi B isolates being genetically different strains, since PFGE may be incapable of differentiating individual strains within a certain highly genetically homogeneous serotype.50,51 Other factors that were discussed previously regarding the relationship between CSalm in the MW and community salmonellosis case numbers N (such as potential survival of Salmonella in sewer biofilm) might also have played a role. Concurrent Detection of Clinically Rare Salmonella Strains. Also concurrently detected in the MW and the clinical Salmonella collections were two S. Uganda strains (named S. Uganda 1 and S. Uganda 2) that differed by three PFGE bands (Figure 4A). In Week 21, a single clinical S. Uganda 1

collections. Nine of the 24 shared strains were detected in both isolate collections within 1 week from each other (i.e., concurrent detection), including one S. Paratyphi B strain that caused a major outbreak with 17 confirmed cases in Hawaii from March to May 2010, which coincided with the beginning of our MW sampling campaign by 4 weeks (Figure 3). The PFGE patterns of the Salmonella isolates from the MW and the clinical collections were indistinguishable (Figure 3A), which together with their concurrent detection (Figure 3B) indicate that they are the same strain and have a common origin. The temporal detection pattern showed that the S. Paratyphi B strain was the most frequently isolated strain in the MW samples during the first 4 weeks of the sampling campaign when the outbreak was ongoing. After clinical cases completely subsided in Week 5, detection of the S. Paratyphi B isolates in MW eventually disappeared after persisting at a low level for two additional weeks. Although the outbreak caused by S. Paratyphi B (Figure 3) was the only clinical outbreak during the sampling period and hence did not permit robust statistical inference, the concurrent detection of this Salmonella strain in both MW and clinical samples supports the potential use of MW for time-sensitive outbreak surveillance. The outbreak-causing S. Paratyphi B led to its numerical dominance within the Salmonella population in the MW samples (e.g., six isolates in Week 2). Although high levels of enteric pathogens are intuitively expected in MW during disease outbreaks, only a few previous studies have reported time-sensitive detection, which include the detection of the Hepatitis A virus and Norovirus in Sweden MW during E

DOI: 10.1021/acs.est.8b00163 Environ. Sci. Technol. XXXX, XXX, XXX−XXX

Article

Environmental Science & Technology

detection of clinical cases and also persisted in MW for three additional weeks. S. Uganda is rarely isolated from human clinical cases in the U.S., with only a median of 48 human isolates per year from 1993 to 2001,52 so the detection of four Uganda clinical isolates in the clinical collection of 338 isolates (1.2%) can be considered rare events. The concurrent detection of two rare Uganda strains in the MW and clinical collections supports the notion that MW-based microbial surveillance has the potential to achieve high sensitivity. The continuous isolation of the same Uganda strains in MW samples for up to 6 weeks after the reported clinical cases suggests that the strains continued to circulate within the community while evading clinical detection. Although it is also possible that the strains were simply surviving in the sewer system, this possibility is considered unlikely because of the observation of a large number of Uganda 1 isolates in the MW samples in Weeks 25 and 26, which was highly unlikely to be primarily caused by surviving in the sewer system. Concurrent detection of the clinical and MW Salmonella isolates also included six additional strains that belong to six different serotypes (Figure 5). The MW and clinical Salmonella isolates were indistinguishable by their PFGE patterns (Figure S2). The S. Amager, S. Chester, S. Takoradi, and S. Newport strains were each detected in one cluster over the entire 54week sampling period in both the MW collection and the clinical collection (Figure 5A−D). S. Javiana (Figure 5E) and S. Derby (Figure 5F) appeared in two separate clusters in the MW collection, though only one cluster for each serotype showed temporal correlation with clinical isolates. S. Newport and S. Javiana ranked third and fourth most abundant among laboratory-confirmed Salmonella infections, while the other serotypes ranked out of the top 20.37 Although the six serotypes

Figure 4. PFGE patterns of two Uganda strains in the clinical (labeled N) and MW (labeled Y) isolate libraries (A) and their detection patterns in the MW and clinical isolate collections over the 54-week sampling period (B).

isolate was detected, while corresponding S. Uganda 1 isolates were detected in the MW sample collected after the clinical detection and persisted in the MW for 7 weeks (Figure 4B). The detection of the S. Uganda 2 strain in three clinical cases occurred in Week 35, 14 weeks after the S. Uganda 1 cluster. Similar to the S. Uganda 1 detection, the S. Uganda 2 isolates were also detected in the MW sample collected after the

Figure 5. Concurrent detection of the Salmonella serotype Amager strain (A), Chester strain (B), Takoradi strain (C), Newport strain (D), Javiana strain (E), and Derby strain (F) in the clinical and MW isolate collections over the 54-week sampling period. F

DOI: 10.1021/acs.est.8b00163 Environ. Sci. Technol. XXXX, XXX, XXX−XXX

Article

Environmental Science & Technology

intensive and time-consuming would be needed in order to make the proposed MW-based microbial surveillance approach practically feasible for routine monitoring.

typically exhibit different clinical detection frequencies in health clinics nationwide, all six strains were only represented by a single clinical case over the sampling period, indicating the rarity of their occurrence in the local clinical isolate collection. The current detection of these clinically rare Salmonella strains in the MW samples provides strong evidence of their association and being from common sources. The concurrent detection of the nine Salmonella strains in the MW and clinical collections (Figures 3, 4, and 5) clearly demonstrates that there is a correlation between the occurrence of salmonellosis in the community and the presence of Salmonella in MW. The grouping of Salmonella isolates into different serotypes and strains was based on indistinguishable PFGE patterns, which is the prevailing method for bacterial strain subtyping for epidemiological investigations.28,29 It should be noted that the PFGE-based strain subtyping was reported to be inadequate for S. Enteritidis because of its high genetic homogeneity,50,51 which requires methods with higher discriminatory power (such as whole genome sequencing) for strain differentiation.53 However, similar limitations have not been reported for any of the serotypes to which the nine Salmonella strains belong; for example, the serotype S. Paratyphi B is known to be highly genetically diverse.54 Overall, this case study of Salmonella and salmonellosis provides multiple lines of evidence supporting the feasibility of using MW as a microbial surveillance platform for enteric diseases in a community and the possibility of using MW to monitor salmonellosis activity in a time-sensitive fashion. First of all, Salmonella concentration in Honolulu MW exhibited positive and statistically significant linear and rank correlation with clinical salmonellosis case numbers over the same time period. Second, the MW Salmonella isolate collection shared 21 serotypes with the clinical Salmonella isolate collection established over the same time period, including nine of the 12 most common clinical serotypes. Finally, and most significantly, nine Salmonella strains, including an outbreakassociated S. Paratyphi B strain and numerous clinically rare strains, were concurrently detected in health clinics and in the MW samples, supporting the feasibility of monitoring enteric pathogens in the MW as a timely indication of community enteric disease activity. Although enteric disease information on the community is undoubtedly embedded within the MW, to achieve the full potential of the proposed MW-based microbial surveillance approach, one would need to overcome several major technological challenges. Although the presence of a linear correlation between CSalm in MW and salmonellosis case number N in the community (Table 1) suggests that directly quantifying CSalm in the MW for the estimation of disease prevalence is possible, the large fluctuation of CSalm and the presence of data outliers can present practical challenges in data interpretation, especially for predictive applications. Further testing of multiple outbreaks caused by Salmonella and other enteric pathogens is required in order to establish statistical confidence for MW-based microbial disease outbreak detection. Many of the clinical Salmonella strains were either not detected at all or not concurrently detected in the MW samples, which was attributed primarily to low MW sampling frequency, low sampling depth, and potential biases introduced by the cultivation-based approach, which was necessary for obtaining individual MW strains for direct comparison with the clinical isolates. Future development of high throughput and high automation molecular methods that are significantly less labor



ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.est.8b00163. Additional information on log JSalm and log RS/E and clinical salmonellosis case number N over the 54-week sampling (Figure S1); the PFGE patterns of the Salmonella isolates belonging to the serotypes S. Amager, S. Chester, S. Takoraki, S. Newport, S. Javiana, and S. Derby (Figure S2) (PDF)



AUTHOR INFORMATION

Corresponding Author

*Phone: 808-956-6024; fax: 808-956-5014; e-mail: taoyan@ hawaii.edu. ORCID

T. Yan: 0000-0001-7500-7993 Present Address #

E.P.: Ecological Sciences, The James Hutton Institute, Craigiebuckler, Aberdeen, Scotland, UK, AB15 8QH. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS The authors thank Mr. Ken Tenno of the Water Quality Laboratory, Department of Environmental Services of the City and County of Honolulu, for collecting sanitary sewage samples and providing wastewater flow data, Dr. Qian Zhang and Dr. Kun Yang for providing laboratory assistance, and Dr. Alexandria Boehm of Stanford Univeristy for sharing the laboratory protocol for Salmonella enumeration. We also thank student investigators Ana-Melissa Kea from the University of Hawaii Science, Technology, Engineering and Mathematics (STEM) program and Sasha Canovali from the Windward Community College Institutional Development Award (IDeA) Network for Biomedical Research Excellence (INBRE) program for performing much of the PFGE testing. This work was financially supported by grants from the National Science Foundation (CBET-0964260 and CBET-1507979 to T.Y.).



REFERENCES

(1) Scallan, E.; Hoekstra, R. M.; Angulo, F. J.; Tauxe, R. V.; Widdowson, M. A.; Roy, S. L.; Jones, J. L.; Griffin, P. M. Foodborne illness acquired in the United States–major pathogens. Emerging Infect. Dis. 2011, 17 (1), 7−15. (2) Scallan, E.; Griffin, P. M.; Angulo, F. J.; Tauxe, R. V.; Hoekstra, R. M. Foodborne illness acquired in the United States–unspecified agents. Emerging Infect. Dis. 2011, 17 (1), 16−22. (3) Cohen, M. L. Changing patterns of infectious disease. Nature 2000, 406 (6797), 762−7. (4) Khabbaz, R. F.; Moseley, R. R.; Steiner, R. J.; Levitt, A. M.; Bell, B. P. Challenges of infectious diseases in the USA. Lancet 2014, 384 (9937), 53−63. (5) World Health Organization. The Global Epidemiology of Infectious Diseases; WHO: Geneva, 2004. (6) Bradley, C. A.; Rolka, H.; Walker, D.; Loonsk, J. BioSense: implementation of a National Early Event Detection and Situational Awareness System. MMWR supplements 2005, 54, 11−19. G

DOI: 10.1021/acs.est.8b00163 Environ. Sci. Technol. XXXX, XXX, XXX−XXX

Article

Environmental Science & Technology (7) Berger, M.; Shiau, R.; Weintraub, J. M. Review of syndromic surveillance: implications for waterborne disease detection. J. Epidemiol. Community health 2006, 60 (6), 543−550. (8) Joe Alper, R. Big Data and Analytics for Infectious Disease Research, Operations, and Policy: Proceedings of a Workshop; The National Academides of Sciences, Engineering and Medicine: Washington, DC, 2016. (9) Castiglioni, S.; Senta, I.; Borsotti, A.; Davoli, E.; Zuccato, E. A novel approach for monitoring tobacco use in local communities by wastewater analysis. Tobaco Control 2015, 24 (1), 38−42. (10) Kim, K. Y.; Lai, F. Y.; Kim, H.-Y.; Thai, P. K.; Mueller, J. F.; Oh, J.-E. The first application of wastewater-based drug epidemiology in five South Korean cities. Sci. Total Environ. 2015, 524-525, 440−446. (11) Venkatesan, A. K.; Halden, R. U. Wastewater treatment plants as chemical observatories to forecast ecological and human health risks of manmade chemicals. Sci. Rep. 2015, 4, 3731. (12) Broek, J. C.; Mom, C. P. Contamination of polder water with Salmonella paratyphi as a consequence of the discharge of sewage. Antonie van Leeuwenhoek 1953, 19 (4), 363−72. (13) Leclerc, H.; Catsaras, M.; Savage, C.; Eymard, C. The isolation of Salmonella from heavily polluted environments. II. Sewage. Annales de l’Institut Pasteur de Lille 1970, 21, 277−293. (14) Langeland, G. Salmonella spp. in the working environment of sewage treatment plants in Oslo, Norway. Appl. Environ. Microbiol. 1982, 43 (5), 1111−1115. (15) Carrington, E. G.; Harman, S. A.; Pike, E. B. Inactivation of Salmonella during anaerobic digestion of sewage sludge. J. Appl. Bacteriol. 1982, 53 (3), 331−4. (16) Shinohara, N.; Tanaka, H.; Saito, T.; Deguchi, J.; Soda, K.; Sugiyama, T.; Ishimaru, Y.; Sonoda, S. Surveillance for typhoid fever in Matsuyama city during 1974−1981 and detection of Salmonella typhi in sewage and river waters. Jpn. J. Med. Sci. Biol. 1983, 36 (3), 191−7. (17) Johnston, W. S.; Hopkins, G. F.; Maclachlan, G. K.; Sharp, J. C. Salmonella in sewage effluent and the relationship to animal and human disease in the north of Scotland. Vet. Rec. 1986, 119 (9), 201− 3. (18) Reilly, W. J.; Oboegbulem, S. I.; Munro, D. S.; Forbes, G. I. The epidemiological relationship between Salmonella isolated from poultry meat and sewage effluents at a long-stay hospital. Epidemiol. Infect. 1991, 106 (1), 1−10. (19) Yang, K.; Pagaling, E.; Yan, T. Estimating the Prevalence of Potential Enteropathogenic Escherichia coli and Intimin Gene Diversity in a Human Community by Monitoring Sanitary Sewage. Appl. Environ. Microbiol. 2014, 80 (1), 119−127. (20) Mead, P. S.; Slutsker, L.; Dietz, V.; McCaig, L. F.; Bresee, J. S.; Shapiro, C.; Griffin, P. M.; Tauxe, R. V. Food-related illness and death in the United States. Emerging Infect. Dis. 1999, 5 (5), 607−25. (21) Voetsch, A. C.; Van Gilder, T. J.; Angulo, F. J.; Farley, M. M.; Shallow, S.; Marcus, R.; Cieslak, P. R.; Deneen, V. C.; Tauxe, R. V. FoodNet estimate of the burden of illness caused by nontyphoidal Salmonella infections in the United States. Clin. Infect. Dis. 2004, 38 (Suppl 3), S127−S134. (22) Valdez, Y.; Ferreira, R. B.; Finlay, B. B. Molecular mechanisms of Salmonella virulence and host resistance. Curr. Top. Microbiol. Immunol. 2009, 337, 93−127. (23) Mead, P. S.; Slutsker, L.; Dietz, V.; McCaig, L. F.; Bresee, J. S.; Shapiro, C.; Griffin, P. M.; Tauxe, R. V. Food-related illness and death in the United States. Emerging Infect. Dis. 1999, 5 (5), 607−625. (24) Swaminathan, B.; Barrett, T. J.; Hunter, S. B.; Tauxe, R. V. PulseNet: the molecular subtyping network for foodborne bacterial disease surveillance, United States. Emerging Infect. Dis. 2001, 7 (3), 382−9. (25) Walters, S. P.; Gonzalez-Escalona, N.; Son, I.; Melka, D. C.; Sassoubre, L. M.; Boehm, A. B. Salmonella enterica diversity in central Californian coastal waterways. Appl. Environ. Microbiol. 2013, 79 (14), 4199−209. (26) USEPA. USEPA Method 1600: Enterococci in Water by Membrane Filtration Using membrane-Enterococcus Indoxyl-β-D-Glucoside Agar (mEI); USEPA Office of Water: Washington DC, 2002.

(27) Malorny, B.; Paccassoni, E.; Fach, P.; Bunge, C.; Martin, A.; Helmuth, R. Diagnostic real-time PCR for detection of Salmonella in food. Appl. Environ. Microbiol. 2004, 70, 7046−7052. (28) CDC. National Enteric Disease Surveillance: Salmonella Surveillance Overview; CDC Enteric Diseases Epidemiology Branch, Atlanta, GA, 2011. (29) CDC. Standard Operating Procedure for PulseNet PFGE of Escherichia coli O157:H7, Escherichia coli non-O157 (STEC), Salmonella serotypes, Shigella sonnei and Shigella flexneri; CDC: Atlanta, GA, 2013. (30) CDC. National Enteric Disease Surveillance: Salmonella Annual Report, 2015; CDC National Center for Emerging and Zoonotic Infectious Diseases: Atlanta, GA, 2017. (31) Varga, C.; Pearl, D. L.; McEwen, S. A.; Sargeant, J. M.; Pollari, F.; Guerin, M. T. Incidence, distribution, seasonality, and demographic risk factors of Salmonella Enteritidis human infections in Ontario, Canada, 2007−2009. BMC Infect. Dis. 2013, 13, 212. (32) de Klerk, N. Commentary: Spearman’s ’The proof and measurement of association between two things. Int. J. Epidemiol. 2010, 39 (5), 1159−61. (33) Bitton, G. Wastewater Microbiology, 3rd ed.; John Wiley & Sons: Hoboken, NJ, 2005. (34) Wery, N.; Lhoutellier, C.; Ducray, F.; Delgenes, J. P.; Godon, J. J. Behaviour of pathogenic and indicator bacteria during urban wastewater treatment and sludge composting, as revealed by quantitative PCR. Water Res. 2008, 42 (1−2), 53−62. (35) Koivunen, J.; Siitonen, A.; Heinonen-Tanski, H. Elimination of enteric bacteria in biological-chemical wastewater treatment and tertiary filtration units. Water Res. 2003, 37 (3), 690−8. (36) Jimenez-Cisneros, B. E.; Maya-Rendon, C.; Salgado-Velazquez, G. The elimination of helminth ova, faecal coliforms, Salmonella and protozoan cysts by various physicochemical processes in wastewater and sludge. Water Sci. Technol. 2001, 43 (12), 179−182. (37) CDC. National Enteric Disease Surveillance: Salmonella Annual Report, 2011; CDC National Center for Emerging and Zoonotic Infectious Diseases: Atlanta, GA, 2013. (38) Naumova, E. N.; Jagai, J. S.; Matyas, B.; DeMaria, A., Jr.; MacNeill, I. B.; Griffiths, J. K. Seasonality in six enterically transmitted diseases and ambient temperature. Epidemiol. Infect. 2007, 135 (2), 281−92. (39) Stephen, D. M.; Barnett, A. G. Effect of temperature and precipitation on salmonellosis cases in South-East Queensland, Australia: an observational study. BMJ. open 2016, 6 (2), e010204. (40) Shelton, J. M.; Kim, L.; Fang, J.; Ray, C.; Yan, T. Assessing the severity of rainfall-derived infiltration and inflow and sewer deterioration based on the flux stability of sewage markers. Environ. Sci. Technol. 2011, 45 (20), 8683−90. (41) Thomson, S. The numbers of pathogenic bacilli in faeces in intestinal diseases. J. Hyg. 1955, 53 (2), 217−24. (42) Imanishi, M.; Rotstein, D. S.; Reimschuessel, R.; Schwensohn, C. A.; Woody, D. H., Jr.; Davis, S. W.; Hunt, A. D.; Arends, K. D.; Achen, M.; Cui, J.; Zhang, Y.; Denny, L. F.; Phan, Q. N.; Joseph, L. A.; Tuite, C. C.; Tataryn, J. R.; Behravesh, C. B. Outbreak of Salmonella enterica serotype Infantis infection in humans linked to dry dog food in the United States and Canada, 2012. J. Am. Vet. Med. Assoc. 2014, 244 (5), 545−53. (43) Kolker, S.; Itsekzon, T.; Yinnon, A. M.; Lachish, T. Osteomyelitis due to Salmonella enterica subsp. arizonae: the price of exotic pets. Clin. Microbiol. Infect. 2012, 18 (2), 167−70. (44) Lafarga, M. A.; Castillo, J.; Navarro, M.; Gomez-Lus, R. Serotypes of Salmonella enterica in sewage effluent from Zaragoza. Comparison with clinical isolates. 1982−1989. Microbiologia 1991, 7 (1), 23−34. (45) Gunn, J. S.; Marshall, J. M.; Baker, S.; Dongol, S.; Charles, R. C.; Ryan, E. T. Salmonella chronic carriage: epidemiology, diagnosis, and gallbladder persistence. Trends Microbiol. 2014, 22 (11), 648−55. (46) Setti, I.; Rodriguez-Castro, A.; Pata, M. P.; Cadarso-Suarez, C.; Yacoubi, B.; Bensmael, L.; Moukrim, A.; Martinez-Urtaza, J. Characteristics and dynamics of Salmonella contamination along the H

DOI: 10.1021/acs.est.8b00163 Environ. Sci. Technol. XXXX, XXX, XXX−XXX

Article

Environmental Science & Technology coast of Agadir, Morocco. Appl. Environ. Microbiol. 2009, 75 (24), 7700−9. (47) Hellmer, M.; Paxeus, N.; Magnius, L.; Enache, L.; Arnholm, B.; Johansson, A.; Bergstrom, T.; Norder, H. Detection of pathogenic viruses in sewage provided early warnings of hepatitis A virus and norovirus outbreaks. Appl. Environ. Microbiol. 2014, 80 (21), 6771−81. (48) Chalmers, R. M.; Robinson, G.; Elwin, K.; Hadfield, S. J.; Thomas, E.; Watkins, J.; Casemore, D.; Kay, D. Detection of Cryptosporidium species and sources of contamination with Cryptosporidium hominis during a waterborne outbreak in north west Wales. J. Water Health 2010, 8 (2), 311−25. (49) Sahlstrom, L.; de Jong, B.; Aspan, A. Salmonella isolated in sewage sludge traced back to human cases of salmonellosis. Lett. Appl. Microbiol. 2006, 43 (1), 46−52. (50) Zheng, J.; Keys, C. E.; Zhao, S.; Meng, J.; Brown, E. W. Enhanced subtyping scheme for Salmonella enteritidis. Emerging Infect. Dis. 2007, 13 (12), 1932−5. (51) Olson, A. B.; Andrysiak, A. K.; Tracz, D. M.; Guard-Bouldin, J.; Demczuk, W.; Ng, L. K.; Maki, A.; Jamieson, F.; Gilmour, M. W. Limited genetic diversity in Salmonella enterica serovar Enteritidis PT13. BMC Microbiol. 2007, 7, 87. (52) CDC. Salmonella surveillance: annual summary, 2001; Centers for Disease Control and Prevention: Atlanta, GA, 2002. (53) Deng, X.; Shariat, N.; Driebe, E. M.; Roe, C. C.; Tolar, B.; Trees, E.; Keim, P.; Zhang, W.; Dudley, E. G.; Fields, P. I.; Engelthaler, D. M. Comparative analysis of subtyping methods against a wholegenome-sequencing standard for Salmonella enterica serotype Enteritidis. J. Clin. Microbiol. 2015, 53 (1), 212−8. (54) Prager, R.; Rabsch, W.; Streckel, W.; Voigt, W.; Tietze, E.; Tschape, H. Molecular properties of Salmonella enterica serotype paratyphi B distinguish between its systemic and its enteric pathovars. J. Clin. Microbiol. 2003, 41 (9), 4270−8.

I

DOI: 10.1021/acs.est.8b00163 Environ. Sci. Technol. XXXX, XXX, XXX−XXX