Systematic review, meta-analysis, and risk ... - ACS Publications

Can we swim yet? Systematic review, meta-analysis, and risk assessment of aging sewage in. 1 surface waters. 2. 3. Alexandria B. Boehm*, Katherine E...
0 downloads 0 Views 668KB Size
Subscriber access provided by Kaohsiung Medical University

Characterization of Natural and Affected Environments

Can we swim yet? Systematic review, meta-analysis, and risk assessment of aging sewage in surface waters Alexandria B Boehm, Katherine E Graham, and Wiley Jennings Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.8b01948 • Publication Date (Web): 06 Aug 2018 Downloaded from http://pubs.acs.org on August 6, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 39

1 2 3 4 5 6 7 8 9 10 11 12

Environmental Science & Technology

Can we swim yet? Systematic review, meta-analysis, and risk assessment of aging sewage in surface waters Alexandria B. Boehm*, Katherine E. Graham, Wiley C. Jennings Department of Civil and Environmental Engineering, Stanford University, Stanford, California, USA 94305 Submitted to Environmental Science & Technology

*

Corresponding author: 473 Via Ortega, Room 189, Dept. of Civil & Environmental Engineering , Stanford University, Stanford, CA 994305. [email protected], (650) 724-9128 (tel)

ACS Paragon Plus Environment

1

Environmental Science & Technology

13 14

Page 2 of 39

TOC Art

15 16

ACS Paragon Plus Environment

2

Page 3 of 39

Environmental Science & Technology

17

Abstract

18

This study investigated the risk of gastrointestinal illness associated with swimming in surface

19

waters with aged sewage contamination. First, a systematic review compiled 333 first order

20

decay rate constants (k) for human norovirus and its surrogates feline calicivirus and murine

21

norovirus, Salmonella, Campylobacter, Escherichia coli O157:H7, Giardia, and

22

Cryptosporidium, and human-associated indicators in surface water. A meta-analysis

23

investigated effects of sunlight, temperature, and water matrix on k. There was a relatively large

24

number of k for bacterial pathogens and some human-associated indicators (n>40), fewer for

25

protozoan and (n=14-22), and few for human norovirus and its Caliciviridae surrogates (n=2-4).

26

Average k ranked: Campylobacter>human-associated markers>Salmonella> E. coli O157:H7 >

27

norovirus and its surrogates >Giardia>Cryptosporidium. Compiled k values were used in a

28

quantitative microbial risk assessment (QMRA) to simulate gastrointestinal illness risk

29

associated with swimming in water with aged sewage contamination. The QMRA used human-

30

associated fecal indicator HF183 as an index for the amount of sewage present and thereby

31

provided insight into how risk relates to HF183 concentrations in surface water. Because

32

exposure to norovirus contributed the majority of risk, and HF183 k is greater than norovirus k,

33

the risk associated with exposure to a fixed HF183 concentration increases with the age of

34

contamination. Swimmer exposure to sewage after it has aged ~3 days results in median risks

35

less than 30/1000. A risk-based water quality threshold for HF183 in surface waters that takes

36

into account uncertainty in contamination age is derived.

37

ACS Paragon Plus Environment

3

Environmental Science & Technology

Page 4 of 39

38

Introduction. Surface waters, including fresh, estuarine, and marine waters, serve as drinking

39

water sources, sites for recreation, sources of food, and essential organism habitat. They are

40

susceptible to microbial pollution1–5 from runoff6, animal feces7, and sewage discharges8.

41

Microbial pollutants include pathogenic protozoa, viruses, and bacteria. Microbial pollution

42

around the world is routinely assessed using concentrations of fecal indicator bacteria (FIB)

43

including Escherichia coli and enterococci as proxies for pathogenic organisms9. FIB presence in

44

bathing waters is linked quantitatively to gastrointestinal illness risk in swimmers when the FIB

45

source is runoff, sewage, or wastewater effluent10–14. However, FIB can be found in a variety of

46

sources including animal feces and environmental reservoirs15–17 making them non-ideal

47

indicators. Therefore, new indicators that are fecal-source associated have been recently

48

developed. For example, human-associated fecal indicators like HF183 and HumM218–21 are

49

highly abundant in, and specific to, human feces.

50 51

Once introduced to surface waters, microbial pollutants are transported and dispersed, and

52

removed from the water column by a variety of processes including inactivation22. Inactivation

53

of microbial pollutants in surface waters can occur via light and dark pathways. Inactivation via

54

dark pathways is caused by lack of nutrients, stress, senescence, and exposure to extracellular

55

biocidal compounds as well as grazing by bacterivorous organisms22,23. Inactivation can also

56

occur via light pathways where photons from the sun induce die-off directly (for example UVB

57

damaging genomic DNA) or indirectly via reactive species (generated when photons interact

58

with sensitizers like humic acids)24,25.

59

ACS Paragon Plus Environment

4

Page 5 of 39

Environmental Science & Technology

60

Several recent studies have compiled inactivation data on Escherichia coli and enterococci, as

61

well as human-associated indicators and pathogens in surface waters26–29. However, a systematic

62

screening of the literature on pathogen and human-associated indicators suggested that available

63

data are considerably more extensive than reported in these studies. A quantitative synthesis of

64

available inactivation data of waterborne pathogens and human indicators would facilitate the

65

modeling of these organisms in surface waters by providing a range of relevant decay rate

66

constants30. Such data would also be useful in quantitative microbial risk assessment (QMRA)

67

scenarios where fecal contamination may be aged31.

68 69

QMRAs have been used extensively to explore how risk of enteric illness varies with bathing

70

water exposure to microbial pollution. For example, QMRA was used to show that swimming in

71

water containing enterococci from different fecal sources, including gull feces, cow feces, and

72

sewage, is associated with different risks32 owing to the fact that different pathogens, with

73

different pathogenicity, are present in the different fecal sources. Other bathing water QMRA

74

applications allowed calculation of risk-based water quality thresholds for human-33 and gull-

75

associated indicators34 for recreational waters; when concentrations of the human- and gull-

76

associated indicators are below the risk-based thresholds, then water is safe for swimming. In

77

these QMRAs, however, the fecal source was assumed to be “unaged”, meaning that the ratio of

78

indicators to pathogens in the fecal source material, an important model input, was assumed to be

79

the same as that in the surface water. In reality, once introduced to the environment, the indicator

80

and pathogens may decay at different rates, thus altering their ratios.

81

ACS Paragon Plus Environment

5

Environmental Science & Technology

Page 6 of 39

82

In this study, a novel systematic review and meta-analysis of decay rate constants for key

83

pathogens and human-associated indicators was conducted. Target-specific decay rate constants

84

were summarized as statistical distributions so they may be readily used in a variety of modeling

85

applications. We subsequently used the decay rate constants in a QMRA to investigate how

86

exposure to fecal contamination from untreated sewage of different ages during swimming

87

affects the resultant gastrointestinal illness risk. The QMRA provides a framework for

88

identifying a risk-based water quality threshold for the human-associated fecal marker HF183.

89 90

Materials and Methods.

91

Systematic review. The systematic review and meta-analysis followed PRISMA guidelines35.

92

The goal of the review was to compile from the peer reviewed literature quantitative information

93

on the decay of waterborne pathogens and human-associated indicators of fecal pollution in

94

surface waters under environmentally-relevant conditions. Pathogens included in the review

95

were human norovirus, Campylobacter, Salmonella, Giardia, Cryptosporidium, and E. coli

96

O157:H7. Human norovirus surrogates in the Caliciviridae family, feline calicivirus and murine

97

norovirus, were also included. Human-associated fecal indicators included were BacHum-UCD,

98

HumM2, and HF183. These organisms were included as they were central to our planned QMRA

99

application.

100 101

Web of Science core collection (search field = topic), Scopus (search field = article title, abstract,

102

keyword), and PubMed (search field = all fields) were searched in August 2017 (Tables 1 and

103

S1). The search terms were “(X) AND (water OR seawater OR stormwater) AND (die-off OR

104

persistence OR survival OR inactivat* OR decay)” where X is the target-specific text (Table 1).

ACS Paragon Plus Environment

6

Page 7 of 39

Environmental Science & Technology

105

Identified articles were assembled and duplicates were removed. Details of the review process,

106

which involved two independent full text reviews of papers, are provided in the supporting

107

information (SI). The inclusion criteria were the paper: (1) contains quantitative data on the

108

decay of the target of interest in raw (unaltered) surface water, (2) is in English, (3) is not a

109

review, presents primary data, and is peer reviewed, (4) does not contain data solely on

110

disinfection treatments such as addition of oxidants or SODIS, (5) included data from decay

111

experiments where the temperature is greater than 4°C and less than 30°C, and (6) describes

112

methods to enumerate the target that are logical and justifiable.

113 114

Decay rate constants were extracted from papers by a single reviewer. First-order decay rate

115

constants (k), in units per day (d-1), calculated from natural log (ln)-transformed concentration

116

data as used in Chick’s law36, were sought. If a study presented k values, then they were

117

extracted from the paper along with any reported errors and model fit values (R2 and/or root

118

mean square error (RMSE)), and unit conversions were applied where appropriate. If a study

119

reported decay model parameters from a model that was not a first-order model (for example a

120

shoulder log-linear model, or biphasic model), then we extracted those reported model

121

parameters and any associated errors and model fit values. If a study only reported T90 or T99

122

values (time to 90% of 99% reduction in concentration, respectively), then they were converted

123

to first order decay rate constants assuming Chick’s law applied. If no first order decay rate

124

constant was reported by the study authors, then Plot Digitizer

125

(http://plotdigitizer.sourceforge.net) was used to digitize the concentration times series appearing

126

in graphs within the publication. To be clear, this included data from studies that only reported

127

decay model parameters from other types (not first-order log-linear) of decay models. k was then

ACS Paragon Plus Environment

7

Environmental Science & Technology

Page 8 of 39

128

calculated as the regression slope of ln(C/Co) versus time (in days) using linear least-squares

129

regression in R. In this formulation C is the concentration at time t, and Co is the concentration at

130

the start of the experiment at t=0. k and its associated error as well as model fit parameters were

131

recorded. In carrying out the linear regression, values reported at or below the detection limit

132

were included if and only if they were not preceded by other consecutive values at or below the

133

detection limit; the value directly reported by the author was used in these cases.

134 135

Once all data were compiled, datasets and model parameters were examined to assess whether a

136

non-linear model was needed to describe decay. The goodness of log-linear model fit to the data

137

(R2 and RMSE), and the number of data points that appeared to “deviate” from the log-linear

138

model were considered. In general, if R2 values were greater than 0.7 and RMSE was relatively

139

small (~1 ln unit), only one data point visually deviated from a straight line fit between time and

140

ln(C/Co), or the non-log-linear model fit was no better than the log-linear fit, then a log-linear

141

curve fit was deemed acceptable.

142 143

In addition to extracting information on the decay of the target of interest, whether the

144

experiment was conducted in (1) freshwater, estuarine water, or seawater and (2) natural sunlight

145

or the dark was also noted. If an experiment was reportedly carried out in sunlight, but at a depth

146

in the water column greater than ~25 cm or in a container that was opaque to UVA and UVB, the

147

experiment was categorized as carried out in the dark. This is justified by previous work

148

suggesting that such experiments do not receive sufficient photons to be deemed affected by

149

sunlight37. The temperature at which the experiment was conducted was also recorded. If a range

150

of temperatures was provided, the mean of the reported range was used. Finally, the method used

ACS Paragon Plus Environment

8

Page 9 of 39

Environmental Science & Technology

151

for target enumeration was noted (culture, microscopy, quantitative PCR (QPCR), reverse-

152

transcription QPCR (RT-QPCR), or propidium monoazide QPCR (PMA-QPCR)).

153 154

Twenty percent of the papers from which data were extracted by a single reviewer were

155

randomly chosen for a second round of data extraction by a different reviewer. Data extracted by

156

the two reviewers were compared to ensure consistency. A single author conducted detailed

157

review of all datasets to identify missing data, data outliers, and data entry mistakes.

158 159

Meta-analysis. Statistical distributions were fit to target-specific k values. Goodness of fit was

160

assessed by visual inspection of residual and Q-Q plots. This yielded log-normal fits for all

161

organisms/targets with the number of k values n>14. For congruity, log-normal distributions

162

were also used when n