Skip to main content

Rapidly measured indicators of recreational water quality and swimming-associated illness at marine beaches: a prospective cohort study

Abstract

Introduction

In the United States and elsewhere, recreational water quality is monitored for fecal indicator bacteria to help prevent swimming-associated illnesses. Standard methods to measure these bacteria take at least 24 hours to obtain results. Molecular approaches such as quantitative polymerase chain reaction (qPCR) can estimate these bacteria faster, in under 3 hours. Previously, we demonstrated that measurements of the fecal indicator bacteria Enterococcus using qPCR were associated with gastrointestinal (GI) illness among swimmers at freshwater beaches. In this paper, we report on results from three marine beach sites.

Methods

We interviewed beach-goers and collected water samples at marine beaches affected by treated sewage discharges in Mississippi in 2005, and Rhode Island and Alabama in 2007. Ten to twelve days later, we obtained information about gastrointestinal, respiratory, eye, ear and skin symptoms by telephone. We tested water samples for fecal indicator organisms using qPCR and other methods.

Results

We enrolled 6,350 beach-goers. The occurrence of GI illness among swimmers was associated with a log10-increase in exposure to qPCR-determined estimates of fecal indicator organisms in the genus Enterococcus (AOR = 2.6, 95% CI 1.3-5.1) and order Bacteroidales (AOR = 1.9, 95% CI 1.3-2.9). Estimates of organisms related to Clostridium perfringens and a subgroup of organisms in the genus Bacteroides were also determined by qPCR in 2007, as was F+ coliphage, but relationships between these indicators and illness were not statistically significant.

Conclusions

This study provides the first evidence of a relationship between gastrointestinal illness and estimates of fecal indicator organisms determined by qPCR at marine beaches.

Peer Review reports

Background

It is usually impractical to test recreational waters directly for the many and diverse pathogenic microorganisms associated with human derived sewage. As a result, recreational waters are often monitored for fecal indicator bacteria. Fecal indicator bacteria such as Enterococcus spp. or Escherichia coli are ordinarily harmless microbes that are commonly found in sewage and other sources of fecal contamination [1]. These fecal indicator bacteria have been statistically associated with gastrointestinal (GI) illness in recreational waters [2, 3]. Standard methods for measuring water quality involve growing fecal indicator bacteria in culture which requires at least 24 hours and may result in incorrect assessments of water quality by beach managers and regulators [4]. Molecular methods, such as quantitative polymerase chain reaction (qPCR), have the ability to detect fecal indicator organisms much faster by targeting and measuring specific genetic markers [5]. We previously demonstrated that the fecal indicator bacteria Enterococcus spp. estimated by qPCR was well-associated with gastrointestinal illness among swimmers at freshwater beaches [6, 7].

Previous studies found different associations between fecal indicator bacteria measured by culture and swimming-associated illness in marine and fresh waters [8, 9], possibly due to differences in the fate and transport of indicator and pathogenic microorganisms [8, 10, 11]. In this report, we extend our research to three marine beaches in the continental United States where we examine the relationships between swimming-associated illnesses and fecal indicator bacteria determined by alternative rapid methods.

Methods

Beach sites and health survey

We conducted studies at three marine beaches affected by treated sewage discharge from nearby Publicly Owned Treatment Works (POTW). In 2005, we studied Edgewater Beach in Biloxi, Mississippi, and in 2007 we studied Goddard Beach in Goddard Memorial State Park in West Warwick, Rhode Island and Fairhope Municipal Beach in Fairhope, Alabama (Figure 1). Each beach site was located within 7 miles or less of a treated sewage discharge outfall from facilities that served populations of at least 15,000. Based on historical records, each beach site showed variability in water quality, but were generally in compliance with local and federal water quality guidelines.

Figure 1
figure 1

Marine beach sites.

Data collection procedures for the health survey have been described previously [6, 7]. In brief, we conducted surveys on weekends and holidays between May and September. Upon arrival beachgoers were provided a pamphlet describing the study. After reviewing the pamphlet, interviewers approached the household group or individual and offered them the opportunity to enroll in the study. On most days, all beachgoers arriving between approximately 11 AM and 4 PM were offered enrollment. Respondents were ineligible if they had completed the study in the previous 30 days, or if there was no adult (18 years of age or older) household member present. Households provided verbal consent and completed an enrollment questionnaire consisting of demographic information, swimming exposures in the previous two weeks, and the presence of underlying health conditions, such as chronic diarrhea, asthma, or skin conditions. As they left the beach for the day, participants completed a questionnaire to ascertain the extent and duration of their contact with water and other activities during their visit to the beach such as contact with sand and food consumption. Ten to twelve days following the beach visit, we telephoned participants and asked about the occurrence of new gastrointestinal, skin, respiratory, eye, or ear symptoms.

The study procedures, questionnaires, protocols and consent process were reviewed and approved by the Institutional Review Board of the Centers for Disease Control and Prevention.

Illness definitions

We considered the following health endpoints consistent with those we previously reported [6, 7].

"Gastrointestinal illness" (GI illness) was defined as any of the following: (1) diarrhea (three or more loose stools in a 24-hour period); (2) vomiting; (3) nausea and stomachache; (4) nausea or stomachache, and interference with regular activities (missed regular activities as a result of the illness).

"Upper respiratory illness" (URI) was defined as any 2 of the following: sore throat, cough, runny nose, cold, or fever.

"Rash" was defined as a rash or itchy skin.

"Eye irritations" were defined as either eye infection or watery eye.

"Earache" was defined as earache, ear infection, or runny ears.

Diarrhea was also considered as a stand alone outcome because it is a commonly used definition of gastroenteritis in population-based surveillance [12, 13].

Participants ill within 3 days before their beach visit were excluded from analysis of the health outcome related to their baseline symptoms.

Water sample collection and analysis

Protocols used for water sample collection have been described [5]. Briefly, we collected two 1-liter water samples at 8:00 AM, 11:00 AM, and 3:00 PM along 3 transects perpendicular to the shoreline. At each transect, we collected one sample in waist-high water (1 m deep) and one in shin-high water (0.3 m deep). Transects were located at least 60 m apart and encompassed the swimming area. Following collection, samples were placed in coolers and maintained on ice at 1 to 4°C. At each water sampling time we recorded environmental conditions, including air and water temperature, cloud cover, rainfall, wind speed and direction, wave height, number of people (on the beach and in the water), boats, animals (number and type, on the beach and in the water), tide stage, and debris.

Water samples were tested for total Enterococcus spp. [5] and total Bacteroidales spp. [14], hereafter referred to as Enterococcus and Bacteroidales respectively, by qPCR using previously published protocols [5, 14]. In 2007, we added qPCR tests for subgroups of Bacteroides[15], and Clostridium spp. or "Clostridium perfringens group" [16], hereafter referred to as "fecal Bacteroides" and "Clostridium", respectively. In 2007, we also included a novel, faster test for F+ (male-specific) coliphage based on a culture and latex agglutination assay (CLAT assay), which also distinguishes F+ RNA coliphage and F+ DNA coliphage [17]. F+ coliphage was also evaluated using a 24-hour spot test according to EPA Method 1601 [18]. Samples were also tested for Enterococcus spp. using EPA Method 1600 [19], a culture-based method. Results for Enterococcus measured by EPA Method 1600 are reported in colony forming units (CFU) per 100 ml sample and for F+ coliphage in most probable number (MPN) per 100 ml. Samples for EPA Method 1600 and qPCR analysis were filtered within 6 hours of collection. The filters were held at -20°C and shipped overnight to EMSL Analytical (Westmount NJ) on dry ice where DNA extraction and qPCR analyses for Enterococcus, Bacteroidales and fecal Bacteroides were conducted. Frozen DNA extracts were sent from EMSL Analytical to the US EPA in Cincinnati where qPCR analysis was conducted for Clostridium.

Primer and probe sequences used for the Enterococcus[5], Bacteroidales[14] and fecal Bacteroides[15] qPCR assays were described previously. Primers used for the amplification of Clostridium were those of the "Clostridium perfringens group" assay targeting about 34 Clostridium species as reported by Rinttilä et. al. [16]. Additional details regarding sample processing, DNA extraction and reaction conditions for qPCR analyses, are provided in additional file 1. In brief, following filtration, DNA was extracted, and polymerase chain reaction (PCR) amplification was carried out using the TaqMan PCR product detection system. The reactions were performed in a thermal cycling instrument (Smart-Cycler System, Cepheid, Sunnyvale, CA) except for Clostridium which was performed on a Model 7900 DNA thermal cycler (Applied Biosystems, Foster City, CA). Both instruments automated the detection and quantitative measurement of the fluorescent signals produced by TaqMan probe degradation during each cycle of amplification.

PCR cycle threshold (CT) measurements of the test sample DNA extracts were compared with those of similarly prepared extracts from calibrator samples containing a known quantity of the target organism cells. Ratios of the target sequences in the test and calibrator samples were converted to estimates of calibrator cell equivalents (CCE) in the test samples [20]. PCR measurements of salmon testes DNA, that was added to the extraction buffer as a source of reference target sequences, were used to estimate the relative efficiency of total DNA recovery from the water sample filters compared to the calibrator samples and to identify potential PCR inhibition [21]. Five-fold dilutions of the water filter and calibration sample extracts were analyzed and water filter extracts giving salmon DNA assay CT values that were > 3 CT units higher than the mean values from the calibration extracts were reanalyzed after additional 5-fold dilutions. If both dilutions failed the salmon CT criterion, the sample was excluded and results were replaced with mean of valid samples collected at the same location, depth and time. Salmon DNA assays were performed in separate reaction tubes.

Quantitation of fecal indicator bacteria by qPCR

Two basic approaches were used to quantify CCE: "delta delta-CT" (CCEΔΔ ) [5–7, 14], and "delta-CT" (CCE) [14]. CCE values were determined using only test sample and batch-mean calibration sample target organism assay CT values (CCEΔ) and also after corrections using CT values from the salmon reference assays (CCEΔΔ method). See previously published manuscripts [5, 14, 22] and additional file 1 for a discussion and description of these calculations. The CCEΔΔ calculation provides quantitative adjustment for partial inhibition [5, 22], but there is some evidence that the salmon reference assay may over correct the CCE quantitation due to a higher sensitivity to matrix inhibitory effects whereas CCEΔ may lead to underestimations [14]. Therefore, both calculation methods were used to determine whether health effects associations were substantially affected by the calculation approach. Additional details on the CCE calculations are provided in the additional file 1.

The lower detection limit was defined as the upper 95% CT bound of the Y-intercept from the pooled standard curve data that was generated from repeated analyses of serially diluted genomic DNA extracts from the calibrator bacterial strains during the study period. Target sequence concentrations in these genomic DNA extracts were determined as previously described [22]. CT values were restricted at this upper bound for all CCE calculations. One-half the calculated CCE was used for non-detects where there was no detection after 45 cycles. Results are reported in qPCR CCE per 100 ml of original sample.

Swimming exposure

Our primary definition of swimming was "body immersion", defined as immersion to the waist or higher. Previously, we observed similar risks of illness for those who immersed their body and those who immersed their head [7]. Non-swimmers were those who reported no water contact.

We used the mean of the log10 fecal indicator organism estimates to represent exposure. We created separate exposure indices based on all samples (daily average), representing an estimate of the overall water quality, and for morning samples (8:00 AM, 11:00 AM) to evaluate whether morning water quality measures were associated with illness. For categorical presentations, indicator groups were established according to quintiles (for indicator bacteria), presence/absence (F+coliphage by CLAT) or at the median (F+ coliphage by SPOT). Non-swimmers were considered unexposed to waterborne fecal indicator organisms.

Statistical analysis

We used logistic regression models to quantify and describe the relationship between estimates of fecal indicator organisms and the risk of illness among swimmers. The predictor of interest was the estimate of fecal indicator organisms. Factor variables representing "beach" were included in all models to control for differences in baseline illness. Robust estimates of variance were used to account for non-independence of observations within households [23–26]. Covariates which could plausibly affect the relationship between water quality and illness, or those which were associated with health outcomes were considered for inclusion in regression models. These included age, sex, race, contact with animals, other swimming in the past 1-week, contact with other persons with diarrhea, distance traveled to the beach, frequency of visits to the beach under study, any other chronic illnesses (GI, skin, asthma), digging in sand, use of insect repellent and sunscreen, and consumption of raw or undercooked meat. We accounted for the following environmental measures: precipitation since 3:00 PM the previous day, bather density, dogs, birds and other animals on the beach, air and water temperature, wind speed, tide stage, wave height, cloud cover, wind direction and boat density, which were often highly correlated, by reducing them to summary measures using principal components analysis. The first two principal components explained 53% of the variability and were characterized by air and water temperature (first component), and wave height, windspeed and precipitation (second component). These two components were included as covariates in regression models. When one or more of the environmental measurements were missing, we used best-subset regression to impute the principal components [27].

For each analysis, the set of covariates was reduced through a change-in-estimate procedure [28], where the parameter of interest was the regression coefficient for the fecal indicator organism, with a criterion of a 5% change in the coefficient. To evaluate differences in the indicator-illness relationship across beaches and age groups, we used likelihood ratio tests to compare full models with multiplicative interaction terms between beach or age group and water quality (which allowed slopes to differ) to restricted models constrained to a single slope. When data were sufficient, we conducted separate analysis for those 10 years of age and under, consistent with our previous report [6]. Stratified analyses for older age groups were not conducted due to small sample size, and infrequent swimming and reporting of illness.

Adjusted Odds Ratios (AORs) estimated from logistic regression models were used to represent the degree of association between fecal indicator estimates and risk of illness. An AOR of 1 indicates no association or a completely flat slope. AORs with a 95% confidence bound including 1 were considered not statistically significant. For graphical presentations, adjusted probabilities of illness were predicted from logistic regression models holding covariates constant at their mean value.

We used Stata version 10.1 for data analysis [29].

Results

Subject recruitment and respondent characteristics

We conducted interviews on 70 study days at the three beaches (Table 1). At Edgewater Beach, data collection was stopped several days early due to the effects of Hurricane Katrina. A total of 9,069 beach-goers were offered enrollment. Of these, 1,715 (19%) refused to participate or were ineligible. Of those who agreed to participate, 6,350 (78%) completed the telephone interview and were eligible for analysis. Selected characteristics of the study population are shown in Table 1. Seventy-five percent of children age 5-10 immersed their body compared to only 26% of those over 55 (data not shown). Swimming was also associated with male gender, non-white race, less frequent visits to the beach, the absence of chronic illnesses and less frequent consumption of raw or undercooked meat (data not shown).

Table 1 Enrollment and selected respondent characteristics

Water quality

We collected and tested a total of 1,242 water samples on study days (Table 2). In terms of overall water quality, Fairhope Beach had the highest geometric mean Enterococcus CFU (21 CFU/100 ml) and Goddard Beach the lowest (4 CFU/100 ml). Individual samples of Enterococcus CFU ranged from below detection to 3,000 CFU/100 ml. Overall, 142 samples (11.5%) exceeded 104 CFU/100 ml, the EPA recommended single sample maximum for marine beaches [30]. Individual samples most frequently exceeded the single sample maximum at Fairhope Beach (99 samples, 23%), followed by Edgewater Beach (32, 8.5%) and Goddard Beach (11 samples, 3%). Enterococcus CCE estimated by qPCR were higher than Enterococcus CFU. For all indicators measured by qPCR, estimates of CCEΔΔ were higher than CCEΔ. Bacteroidales CCE were highest among the indicators.

Table 2 Fecal indicator bacteria estimates at marine beaches

Fifty-six percent (100/222) of samples at Fairhope Beach and 65% (203/425) of samples at Goddard Beach were positive for F+ coliphage by the 24 hour SPOT test. Fewer samples were positive for F+ coliphage by the CLAT assay. At Fairhope Beach, 4% (8/228) and 14% (14/224) of samples indicated the presence of F+ RNA and F+ DNA coliphage, respectively. At Goddard Beach, 8% (31/425) and 9% (37/423) of samples tested indicated the presence of F+ RNA and F+ DNA coliphage, respectively.

Illness and water quality

Bacterial indicators of water quality and illness

The number and crude (unadjusted) percentage of respondents reporting illness among non-swimmers and swimmers for quintiles of exposure based on daily averages are shown for all subjects in Table 3 for qPCR CCEΔΔ and for Enterococcus CFU. For Enterococcus CCE, Bacteroidales CCE and Enterococcus CFU, the crude cumulative incidence of both GI illness and diarrhea increased with increasing levels of exposure and peaked in the highest exposure categories. On days when the daily geometric mean of Bacteroidales CCEΔΔ was highest (over 3,530 CCE/100 ml), approximately 12% of swimmers reported GI illness compared to 6% among non-swimmers and 4% among swimmers on days when Bacteroidales CCEΔΔ was low (< 542 CCE/100 ml). Similar patterns were noted for Enterococcus CFU and CCE.

Table 3 Illness among non-swimmers and swimmers by quintiles of daily average indicator exposures

AORs of illness among swimmers with respect to indicator density are shown for all subjects (Table 4). The risks of both GI illness and diarrhea were significantly associated with exposure to Enterococcus and Bacteroidales CCEΔ and CCEΔΔ (Figures 2 and 3). A log10 increase in the daily average of Enterococcus and Bacteroidales CCE was associated with an approximate doubling or greater (AORs≈2 or more) in the risk of GI illness (Table 4). Morning estimates of Bacteroidales and Enterococcus CCE (8:00 and/or 11:00 AM) were also significantly associated with both GI illness and diarrhea. Associations between GI illness and fecal Bacteroidales were generally positive but not statistically significant. Clostridium CCEΔΔ showed a strong statistically significant association with GI illness (AOR = 1.94, 95% CI 1.16-3.23, Table 4), but the relationship with Clostridium CCEΔ was weaker (AOR = 1.81, 95% CI 0.96-3.39), Table 4).

Table 4 Adjusted Odds Ratios for illness risk among swimmers for a 1 log10 increase in indicator density1
Figure 2
figure 2

Estimated probability of GI illness among swimmers as a function of daily averages of Enterococcus and Bacteroidales qPCR CCE ΔΔ .

Figure 3
figure 3

Estimated probability of Diarrhea among swimmers as a function of daily averages of Enterococcus and Bacteroidales qPCR CCE ΔΔ .

Enterococcus CFU were also positively associated with GI illness, but the associations were not statistically significant. None of the associations with non-enteric illnesses were consistently positive (data for rash, earache and eye irritations see additional file 2, Table S1).

Children 10 years of age and under were examined separately (see additional file 2, Table S2) but there was no evidence of an increased susceptibility to illness with exposure to fecal indicator bacteria. Statistical models which allowed a separate slope for children showed no improvement over models with a single slope for all subjects.

Associations between GI illness and qPCR CCE appeared to strengthen with more intense or prolonged water exposure, although results were based on few subjects. For example, among those exposed to water greater than 90 minutes (N = 751), AORs were 6.4 (95% CI: 1.2-33) and 7.14 (95% CI: 1.4-37) for associations between Enterococcus CCEΔΔ and GI illness and diarrhea, respectively. Among subjects who swallowed water (N = 632) the AOR for GI illness with respect to the daily average Enterococcus CCEΔΔ was 8.9 (95% CI: 2.2-37).

Associations between GI illness and qPCR CCE were robust to many ways of handling and interpreting the data. Allowing qPCR CT values to extend to either 40 or 45 for CCE calculations and using different approaches to impute results for non-detected samples [31] such as a maximum likelihood, regression on order statistics, or Kaplan-Meier estimates did not affect the results or the interpretation (data not shown). Excluding those with swimming contact in the previous 1-week reduced the sample size, but also had little affect on the results (for example: AOR = 2.4, 95% CI 1.04-5.5 and AOR = 2.3, 95% CI 1.3-3.8 for associations between GI illness and Enterococcus and Bacteroidales CCEΔΔ, respectively). Restricting the definition of swimmers to those who immersed their head also had little effect on the results as did using averages of water quality approximately specific to the swimmer's estimated time of day and location of swimming exposure (data not shown).

There was no evidence of systematic differences in the relationships between GI illness and diarrhea and Enterococcus or Bacteroidales CCE exposure across the three beaches. Likelihood ratio tests comparing a full model which allowed the slope to vary by beach was no better than a model that allowed a common slope (for example: p = 0.51 and 0.90 for GI illness associations with Enterococcus and Bacteroidales CCEΔΔ, respectively).

F+ Coliphage and Illness

The odds of GI illness was higher among swimmers compared to non-swimmers on days when F+ DNA or F+ RNA coliphage was detected by the CLAT assay or F+ coliphage was detected by the SPOT assay (Table 5). Positive associations between GI illness and F+ RNA and F+ DNA coliphage measured by the CLAT assay among swimmers were not statistically significant (Table 5), although statistical power was limited due to relatively few positive results. F+ coliphage measured by the SPOT assay was also not associated with GI illness among swimmers (Table 5). Other illnesses did not show a relationship with the presence of coliphage (data not shown).

Table 5 GI illness by swimming status and presence of F+ coliphage by 24-hour SPOT Assay and 5-hour CLAT Assay

Illness below and above current indicator criteria

Eleven of the 70 days studied exceeded the Enterococcus CFU guideline geometric mean value of 35 CFU per 100. On these days, we observed elevated odds of diarrhea, respiratory illness and earache among swimmers compared to non-swimmers (Table 6). For all categories of illness, there were no differences in the odds of illness among swimmers on days when Enterococcus exceeded 35 CFU/100 ml compared to days when Enterococcus was below 35 CFU/100 ml (Table 6).

Table 6 Illness by exposure to Enterococcus CFU, above and below EPA guidelines

Discussion

This is the first study to demonstrate a relationship between rapid, molecular measures of fecal indicator organisms and swimming-associated illness at marine beaches. As we observed at freshwater beaches, of the various symptom categories those most consistently associated with fecal indicator bacteria were symptoms of GI illness. Although daily averaged samples generally showed the strongest associations with illness, the six morning samples collected at either 8:00 AM or 11:00 AM also were associated with illness occurrence.

The different CCE calculation approaches had little effect on the interpretation of the results. The finding that associations were unaffected by the type of CCE calculation and that the associations strengthened among those with increased exposure (swallowing water, and among those exposure longer than 90 minutes) reinforced the validity of the results.

Average water quality as measured by Enterococcus CFU at these three marine beaches was relatively good. Most other epidemiology studies at marine beaches have reported higher levels of Enterococcus CFU (although direct comparisons are difficult due to differences in the way average indicator densities were reported) [32–34]. A series of randomized studies in the United Kingdom, however, observed exposure-response associations at marine beaches impacted by sewage with low levels of Enterococcus CFU exposures [35]. Here, we report an association between Enteroccocus and Bacteroidales CCE and GI illness despite relatively low levels of culturable Enterococcus. This is consistent with our previous findings at freshwater beaches and provides further evidence that fecal indicators estimated by qPCR are a sensitive marker of poor water quality and subsequent health effects among swimmers at beach sites impacted by treated sewage discharge. Compared to the freshwater beach sites we previously studied [6], average Enterococcus CCE were slightly higher and average Enterococcus CFU were lower. The relationship between Enterococcus qPCR CCE and Enterococcus CFU can be influenced by many factors and may vary from beach to beach [36]. Environmental factors, such as sunlight, may affect the differential persistence of the qPCR signal and culturable organisms [37–39].

The estimate of the association between Enterococcus CCEΔΔ and GI illness among swimmers (AOR = 2.59, 95%CI 1.29-5.11) was stronger than we reported at freshwater beaches (AOR = 1.26, 95% CI 95% CI 1.05-1.51) [6]. Although the 95% confidence bounds of the AORs for marine and freshwater overlap, such an overlap does not necessarily indicate no statistical difference [40, 41]. Our sample size at the marine beaches was smaller and as a result, the 95% confidence bounds are wider and the AOR is less precise. Previous studies [8] also observed a higher illness rate and slope for the association between culturable Enterococcus and GI illness in marine water which the authors suggested may have been due to a greater die off of indicators relative to pathogens in marine waters. We anticipate more detailed comparisons of these results will be the focus of future reports which should also consider baseline (non-swimmer) illness. We observed positive relationships between gastrointestinal symptoms and other fecal indicators, Clostridium, fecal Bacteroides CCE, as well as F+ coliphage, but since associations were not consistently statistically significant among swimmers we cannot make conclusions regarding these indicators.

Clostridium perfringens has been suggested for use as an indicator of fecal contamination in tropical environments [42], but it has not been previously associated with illness in epidemiological studies [3]. Male-specific coliphage also showed some evidence of an association with GI illness at a marine site in Mission Bay, California [32] but this was based on few positive samples. Levels of F+ RNA coliphage were associated with gastrointestinal illnesses among canoeists in freshwater rivers [43].

Limitations

This study was conducted at beaches in a temperate climate with nearby treated point source sewage discharges and results may not be directly applicable to sites affected by fecal contamination from other types of sources or sites with different climates. Human viruses such as norovirus have been identified as likely contributors to excess GI illness among swimmers at beaches impacted by treated sewage discharges [44]. Some human pathogens, particularly human enteric viruses, are unlikely to be associated with non-human sources of fecal indicator bacteria. Results from some recent studies support this limitation. A recently conducted randomized study in Flordia at a sub-tropical beach without known sources of sewage showed no association between fecal indicator bacteria levels and GI illness (including Enterococcus, and three Bacteroidales markers measured by qPCR [34]). This study did observe an association between skin rash and culturable Enterococcus[34, 45], which they suggested may have been marking conditions favorable to pathogens which could cause skin infections [34]. A study at a California marine beach site also failed to find any robust associations between fecal indicator organisms and swimming-associated illness, including Enterococcus and Bacteroidales measured by qPCR [32]. This study, conducted at a beach that was impacted by widespread and diffuse sources of fecal contamination but no point source of contamination [46]. Source tracking studies identified birds as the predominant source of fecal pollution [47]. The results we report from these temperate marine beach sites also may not be applicable to a humid tropical environmental where fecal indicator bacteria have been reported to persist or accummulate in the soil [48].

We used an observational cohort study design which has been used by numerous others both historically and recently to evaluate the risks associated with recreational water exposures (for example [32, 33, 49–52]). This design allowed us to establish associations between average water quality and the subsequent risk of illness among those exposed. One disadvantage of this approach is that such measures of average water quality may not represent an individual swimmer's true exposure. However, in our study, swimmers were exposed for long durations (an average of 74 minutes, with 25% exposed more than 2 hours) and entered the water at multiple time periods and locations, making attempts to characterize individual exposure difficult and impractical. This prolonged and discontinuous exposure may explain why attempts to create time and location specific averages of water quality measures for each swimmer showed no advantage over broader characterizations of water quality. An advantage of the study design we used is the ability to field the study over a wide range of study days which resulted in varying average water quality conditions and enrollment of a large study population. Furthermore, our sampling design allowed us to demonstrate that estimates Enterococcus and Bacteroidales CCE in six morning samples were associated with subsequent GI illness among swimmers exposed on that day.

Some of the health endpoints were non-specific, and may have been affected by recall bias. Broad endpoints accounted for the diverse range of symptoms potentially associated with recreational water exposure but such broad symptoms may obscure more specific effects of water quality and swimming exposure. The associations Enterococcus and Bacteroidales CCE and gastrointestinal symptoms, however, were robust to different definitions (e.g., diarrhea). While swimmers may have been more likely to recall illness than non-swimmers, it is unlikely such a recall bias would occur among swimmers at varying levels of water quality.

Conclusion

Rapid, molecular measures of water quality were associated with illnesses among swimmers at marine beach sites located in a temperate climate with nearby treated sewage discharges. The results provide further evidence that such indicators are a marker of GI illness risk at such beach sites.

Abbreviations

QPCR:

Quantitative polymerase chain reaction

AOR:

Adjusted odds ratio

P:

P-value

CI:

Confidence interval

GI:

Gastrointestinal

URI:

Upper respiratory illness

POTW:

Publicly owned treatment works

CCE:

Calibrator cell equivalents

CCEΔΔ :

Calibrator cell equivalents, delta-delta cycle threshold calculation

CCEΔ :

Calibrator cell equivalents, delta cycle threshold calculation

CT:

Cycle Threshold

CFU:

Colony forming units

CLAT:

Culture and latex agglutination assay for F+ coliphage

SPOT:

24-hour spot test for F+ coliphage.

References

  1. Gerba CP: Indicator Microorganisms. Environmental Microbiology. Edited by: Maier RM, Pepper IL, Gerba CP. 2000, San Diego, CA: Academic Press, 491-503.

    Google Scholar 

  2. Pruss A: Review of epidemiological studies on health effects from exposure to recreational water. Int J Epidemiol. 1998, 27: 1-9. 10.1093/ije/27.1.1.

    Article  CAS  Google Scholar 

  3. Wade TJ, Pai N, Eisenberg JN, Colford JMJ: Do U.S. Environmental Protection Agency Water Quality Guidelines for Recreational Waters Prevent Gastrointestinal Illness? A Systematic Review and Meta-analysis. Environmental Health Perspectives. 2003, 111: 1102-1109.

    Article  Google Scholar 

  4. Kim JH, Grant SB: Public Mis-Notification of Coastal Water Quality: A Probabilistic Evaluation of Posting Errors at Huntington Beach, California. Environmental Science and Technology. 2004, 38: 2497-2504. 10.1021/es034382v.

    Article  CAS  Google Scholar 

  5. Haugland RA, Siefring SC, Wymer LJ, Brenner KP, Dufour AP: Comparison of Enterococcus measurements in freshwater at two recreational beaches by quantitative polymerase chain reaction and membrane filter culture analysis. Water Res. 2005, 39: 559-568. 10.1016/j.watres.2004.11.011.

    Article  CAS  Google Scholar 

  6. Wade TJ, Calderon RL, Brenner KP, Sams E, Beach M, Haugland R, Wymer L, Dufour AP: High sensitivity of children to swimming-associated gastrointestinal illness: results using a rapid assay of recreational water quality. Epidemiology. 2008, 19: 375-383. 10.1097/EDE.0b013e318169cc87.

    Article  Google Scholar 

  7. Wade TJ, Calderon RL, Sams E, Beach M, Brenner KP, Williams AH, Dufour AP: Rapidly measured indicators of recreational water quality are predictive of swimming-associated gastrointestinal illness. Environ Health Perspect. 2006, 114: 24-28. 10.1289/ehp.8273.

    Article  Google Scholar 

  8. Dufour A: Health effects criteria for fresh recreational waters. Tech Rep EPA-600-1-84-004, US Environmental Protection Agency. 1984

    Google Scholar 

  9. Cabelli V: Health effects criteria for marine recreational waters. Tech Rep USEPA Report Number: EPA-600/1-80-031, US Environmental Protection Agency. 1983

    Google Scholar 

  10. Vaccaro RF, Briggs MP, Carey CL, Ketchum BH: Viability of Escherichia coli in sea water. Am J Public Health. 1950, 40: 1257-1266. 10.2105/AJPH.40.10.1257.

    Article  CAS  Google Scholar 

  11. Nasser AM, Telser L, Nitzan Y: Effect of sunlight on the infectivity of Cryptosporidium parvum in seawater. Can J Microbiol. 2007, 53: 1101-1105. 10.1139/W07-043.

    Article  CAS  Google Scholar 

  12. Herikstad H, Yang S, Van Gilder TJ, Vugia D, Hadler J, Blake P, Deneen V, Shiferaw B, Angulo FJ: A population-based estimate of the burden of diarrhoeal illness in the United States: FoodNet, 1996-7. Epidemiol Infect. 2002, 129: 9-17.

    CAS  Google Scholar 

  13. Imhoff B, Morse D, Shiferaw B, Hawkins M, Vugia D, Lance-Parker S, Hadler J, Medus C, Kennedy M, Moore MR, Van Gilder T: Burden of self-reported acute diarrheal illness in FoodNet surveillance areas, 1998-1999. Clin Infect Dis. 2004, 38 (Suppl 3): S219-226. 10.1086/381590.

    Article  Google Scholar 

  14. Siefring S, Varma M, Atikovic E, Wymer L, Haugland RA: Improved real-time PCR assays for the detection of fecal indicator bacteria in surface waters with different instrument and reagent systems. Journal of Water Health. 2008, 6: 225-237. 10.2166/wh.2008.022.

    Article  CAS  Google Scholar 

  15. Converse RR, Blackwood AD, Kirs M, Griffith JF, Noble RT: Rapid QPCR-based assay for fecal Bacteroides spp. as a tool for assessing fecal contamination in recreational waters. Water Res. 2009, 43: 4828-4837. 10.1016/j.watres.2009.06.036.

    Article  CAS  Google Scholar 

  16. Rinttilä T, Kassinen A, Malinen E, Krogius L, Palva A: Development of an extensive set of 16 S rDNA-targeted primers for quantification of pathogenic and indigenous bacteria in faecal samples by real-time PCR. J Appl Microbiol. 2004, 97: 1166-1177. 10.1111/j.1365-2672.2004.02409.x.

    Article  Google Scholar 

  17. Love DC, Sobsey MD: Simple and rapid F+ coliphage culture, latex agglutination, and typing assay to detect and source track fecal contamination. Appl Environ Microbiol. 2007, 73: 4110-4118. 10.1128/AEM.02546-06.

    Article  CAS  Google Scholar 

  18. USEPA: Method 1601: Male-specific (F+) and Somatic Coliphage in Water by Two-step Enrichment Procedure. Tech Rep EPA 821-R-01-030. 2001

    Google Scholar 

  19. USEPA: Method 1600: Membrane Filter Test for Enterococci in Water. Tech Rep EPA-821-R-97-004, Office of Water. 1997

    Google Scholar 

  20. Applied Biosystems: User Bulletin 2: ABI PRISM 7700 Sequence Detection System. Tech rep, Applied Biosystems Corporation. 1997

    Google Scholar 

  21. Brinkman NE, Haugland RA, Wymer LJ, Byappanahalli M, Whitman RL, vesper SJ: Evaluation of a rapid, quantitative real-time PCR method for enumeration of pathogenic Candida cells in water. Appl Environ Microbiol. 2003, 69: 1775-1782. 10.1128/AEM.69.3.1775-1782.2003.

    Article  CAS  Google Scholar 

  22. Chern EC, Brenner KP, Wymer L, Haugland RA: Comparison of Fecal Indicator Bacteria Densities in Marine Recreational Waters by QPCR. Water Quality Exposure and Health. 2009, 1: 203-214. 10.1007/s12403-009-0019-2.

    Article  Google Scholar 

  23. Huber P: The behavior of maximum likelihood estimates under non-standard conditions. Proceedings of the Fifth Berkeley Symposium on Mathematical Statistics and Probability. 1967, Berkeley, CA: University of California Press, 1: 221-233.

    Google Scholar 

  24. White H: A heteroskedasticity-consistent covariance matrix estimator and a direct test for heteroskedasticity. Econometrica. 1980, 48: 817-830. 10.2307/1912934.

    Article  Google Scholar 

  25. Rogers W: Regression standard errors in clustered samples. Stata Technical Bulletin. 1993, 13: 19-23.

    Google Scholar 

  26. Williams RL: A note on robust variance estimation for cluster-correlated data. Biometrics. 2000, 56: 645-6. 10.1111/j.0006-341X.2000.00645.x.

    Article  CAS  Google Scholar 

  27. Little R, Rubin D: Statistical Analysis with Missing Data. 2002, New York: Wiley, 2

    Google Scholar 

  28. Greenland S, Rothman KJ: Modern epidemiology. 1998, Philadelphia, PA:: Lippincott-Raven

    Google Scholar 

  29. StataCorp: Stata Statistical Software: Release 10.1. College Station, Texas: StataCorp LP. 2008

    Google Scholar 

  30. USEPA: Ambient Water Quality for Bacteria-1986. Tech Rep EPA-440/5-84-002, Office of Water Regulations and Standards. 1986

    Google Scholar 

  31. Helsel DR: More than obvious: better methods for interpreting nondetect data. Environ Sci Technol. 2005, 39: 419A-423A. 10.1021/es053368a.

    Article  CAS  Google Scholar 

  32. Colford JM, Wade TJ, Schi KC, Wright CC, Griffith JF, Sandhu SK, Burns S, Sobsey M, Lovelace G, Weisberg SB: Water quality indicators and the risk of illness at beaches with nonpoint sources of fecal contamination. Epidemiology. 2007, 18: 27-35. 10.1097/01.ede.0000249425.32990.b9.

    Article  Google Scholar 

  33. Cabelli VJ, Dufour AP, McCabe LJ, Levin MA: Swimming-associated gastroenteritis and water quality. American Journal of Epidemiology. 1982, 115: 606-616.

    CAS  Google Scholar 

  34. Sinigalliano CD, Fleisher JM, Gidley ML, Solo-Gabriele HM, Shibata T, Plano LR, Elmir SM, Wanless D, Bartkowiak J, Boiteau R, Withum K, Abdelzaher AM, He G, Ortega C, Zhu X, Wright ME, Kish J, Hollenbeck J, Scott T, Backer LC, Fleming LE: Traditional and molecular analyses for fecal indicator bacteria in non-point source subtropical recreational marine waters. Water Research. 2010, 44: 3763-3772. 10.1016/j.watres.2010.04.026.

    Article  CAS  Google Scholar 

  35. Fleisher JM, Jones F, Kay D, Stanwell-Smith R, Wyer M, Morano R: Water and non-water-related risk factors for gastroenteritis among bathers exposed to sewage-contaminated marine waters. International Journal of Epidemiology. 1993, 22: 698-708. 10.1093/ije/22.4.698.

    Article  CAS  Google Scholar 

  36. Telech JW, Brenner KP, Haugland R, Sams E, Dufour AP, Wymer L, Wade TJ: Modeling Enterococcus densities measured by quantitative polymerase chain reaction and membrane filtration using environmental conditions at four Great Lakes beaches. Water Res. 2009, 43: 4947-4955. 10.1016/j.watres.2009.07.002.

    Article  CAS  Google Scholar 

  37. Walters SP, Yamahara KM, Boehm AB: Persistence of nucleic acid markers of health-relevant organisms in seawater microcosms: implications for their use in assessing risk in recreational waters. Water Res. 2009, 43: 4929-4939. 10.1016/j.watres.2009.05.047.

    Article  CAS  Google Scholar 

  38. Bae S, Wuertz S: Rapid decay of host-specific fecal Bacteroidales cells in seawater as measured by quantitative PCR with propidium monoazide. Water Res. 2009, 43: 4850-4859. 10.1016/j.watres.2009.06.053.

    Article  CAS  Google Scholar 

  39. Boehm AB, Griffith J, McGee C, Edge TA, Solo-Gabriele HM, Whitman R, Cao Y, Getrich M, Jay JA, Ferguson D, Goodwin KD, Lee CM, Madison M, Weisberg SB: Faecal indicator bacteria enumeration in beach sand: a comparison study of extraction methods in medium to coarse sands. J Appl Microbiol. 2009, 107: 1740-1750. 10.1111/j.1365-2672.2009.04440.x.

    Article  CAS  Google Scholar 

  40. Mulla ZD, Cole SR: Re: "epidemiology of salmonellosis in California, 1990-1999: morbidity, mortality, and hospitalization costs". Am J Epidemiol. 2004, 159: 104-10.1093/aje/kwh007.

    Article  Google Scholar 

  41. Wolfe R, Hanley J: If we're so different, why do we keep overlapping? When 1 plus 1 doesn't make 2. CMAJ. 2002, 166: 65-66.

    Google Scholar 

  42. Fujioka RS: Monitoring coastal marine waters for spore-forming bacteria of faecal and soil origin to determine point from non-point source pollution. Water Sci Technol. 2001, 44: 181-188.

    CAS  Google Scholar 

  43. Lee JV, Dawson SR, Ward S, Surman SB, Neal KR: Bacteriophages are a better indicator of illness rates than bacteria amongst users of a white water course fed by a lowland river. Water Science and Technology. 1997, 35: 165-170. 10.1016/S0273-1223(97)00253-9.

    Article  CAS  Google Scholar 

  44. Soller JA, Bartrand T, Ashbolt NJ, Ravenscroft J, Wade TJ: Estimating the primary etiologic agents in recreational freshwaters impacted by human sources of faecal contamination. Water Res. 2010, 44: 4736-4747. 10.1016/j.watres.2010.07.064.

    Article  CAS  Google Scholar 

  45. Fleisher JM, Fleming LE, Solo-Gabriele HM, Kish JK, Sinigalliano CD, Plano L, abd John D, Wang SME, Withum K, Shibata T, Gidley ML, Abdelzaher A, He G, Ortega C, Zhu X, Wright M, Hollenbeck J, Backer LC: The BEACHES Study: health effects and exposures from non-point source microbial contaminants in subtropical recreational marine waters. International Journal of Epidemiology. 2010, Advance Access published June 3, 2010

    Google Scholar 

  46. Schiff K, Kinney P: Tracking sources of bacterial contamination in stormwater discharges to Mission Bay, California. Water Environment Research. 2001, 73: 534-542. 10.2175/106143001X139605.

    Article  CAS  Google Scholar 

  47. Gruber SJ, Kay LM, Kolb R, Henry K: Mission Bay bacterial source identification study. Stormwater. 2005, 6: 40-51.

    Google Scholar 

  48. Fujioka R: Monitoring coastal marine waters for spore-forming bacteria of faecal and soil origin to determine point from non-point source pollution. Water Science & Technology. 2001, 44: 181-188.

    CAS  Google Scholar 

  49. McBride G, Salmondo C, Bandaranayake D, Turner S, Lewis G, Till D: Health effects of marine bathing in New Zealand. International Journal of Environmental Health Research. 1998, 8: 173-189. 10.1080/09603129873462.

    Article  CAS  Google Scholar 

  50. Marion JW, Lee J, Lemeshow S, Buckley TJ: Association of gastrointestinal illness and recreational water exposure at an inland U.S. beach. Water Res. 2010, 44: 4796-4804. 10.1016/j.watres.2010.07.065.

    Article  CAS  Google Scholar 

  51. Haile RW, Witte JS, Gold M, Cressey R, McGee C, Millikan RC, Glasser A, Harawa N, Ervin C, Harmon P, Harper J, Dermand J, Alamillo J, Barrett K, Nides M, Wang G: The health effects of swimming in ocean water contaminated by storm drain runoff. Epidemiology. 1999, 10: 355-363. 10.1097/00001648-199907000-00004.

    Article  CAS  Google Scholar 

  52. Corbett SJ, Rubin GL, Curry GK, Kleinbaum DG: The health effects of swimming at Sydney beaches. The Sydney Beach Users Study Advisory Group. American Journal of Public Health. 1993, 83: 1701-1706. 10.2105/AJPH.83.12.1701.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

The authors would like to gratefully acknowledge the following individuals and agencies for their support and assistance. Ann H. Williams, Lucas Neas, Pauline Mendola, Mike Ray, (US EPA, Chapel Hill, North Carolina), for their assistance in administering the study. Westat Inc., including Karen Della Torre, Kurt Patrizi, Robert Clickner, Sharon Jasim-Hanif and many others. Gunther Craun, Charles Poole, Abdel El-Shaawari, Marty Kanarek, John Standridge and John Archer for their helpful comments regarding conduct of the study, analysis of data and reporting of study results. US EPA field monitors Joel Hansel, Edward Hudgens, Holly Wirick, Greg Hellyer, Danelle Lobdell, Beth Leamond, James Scott, Scott Rhoney, Markey Johnson. Bob Mitchell (Fairhope, Alabama) for providing a site for weather gauge installation. The EMSL Laboratory in Westmont, NJ, the Severn-Trent Laboratory in Mobile, Alabama, and the BAL Laboratory in Cranston, RI with special thanks to Darlene Capuano. The authors dedicate this manuscript to the memory of Dr. Rebecca Calderon.

This document has been reviewed in accordance with U.S. Environmental Protection Agency policy and approved for publication. Mention of trade names or commercial products does not constitute endorsement or recommendation for use. The findings and conclusions in this report are those of the authors and do not represent the views of the Centers for Disease Control and Prevention or the EPA.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Timothy J Wade.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

TW led study planning and implementation, data analysis, result interpretation and manuscript preparation. ES led survey and environmental data collection. KB led water sampling design and microbiological analysis. CR conducted microbiological analyses and coliphage testing. EC and RH conducted qPCR analysis and interpreted qPCR results. DL provided expertise on the rapid coliphage CLAT assay. QL conducted and provided input on statistical analyses. RN provided expertise on the Fecal Bacteroides assay. MB conceptualized the study, designed questionnaires and developed study materials. AD conceptualized the study, selected beach sites, and interpreted the results. All authors contributed to writing the manuscript and all authors read and approved the final manuscript.

Electronic supplementary material

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Wade, T.J., Sams, E., Brenner, K.P. et al. Rapidly measured indicators of recreational water quality and swimming-associated illness at marine beaches: a prospective cohort study. Environ Health 9, 66 (2010). https://doi.org/10.1186/1476-069X-9-66

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1476-069X-9-66

Keywords