New Performance Metrics for Quantitative ... - ACS Publications

Aug 16, 2013 - New Performance Metrics for Quantitative Polymerase Chain. Reaction-Based Microbial Source Tracking Methods. Dan Wang,*. ,†. Hyatt C...
1 downloads 0 Views 1MB Size
Letter pubs.acs.org/journal/estlcu

New Performance Metrics for Quantitative Polymerase Chain Reaction-Based Microbial Source Tracking Methods Dan Wang,*,† Hyatt C. Green,‡ Orin C. Shanks,‡ and Alexandria B. Boehm† †

Environmental and Water Studies, Department of Civil and Environmental Engineering, Stanford University, Stanford, California 94305, United States ‡ Office of Research and Development, National Risk Management Research Laboratory, U.S. Environmental Protection Agency, Cincinnati, Ohio 45268, United States S Supporting Information *

ABSTRACT: The ability to select the best quantitative microbial source tracking method for a particular application is paramount. Binary sensitivity and specificity metrics are not adequate to describe the performance of the quantitative methods because the estimates depend on the amount of material tested and the limit of detection. We introduce a new framework for comparing the performance of quantitative fecal source identification methods.



INTRODUCTION Microbial pollution of surface waters is a primary public health concern in the United States and worldwide.1,2 Pollution can originate from a number of fecal and nonfecal sources,3 and different sources can pose different levels of human health risk.4 Microbial source tracking (MST) methods seek to differentiate fecal pollution sources to enhance water quality management efforts. MST method performance is most commonly based on estimates of sensitivity and specificity obtained through the analysis of a set of samples of known fecal origin,5−10 where sensitivity (true positive rate) measures the proportion of actual positives that are correctly identified and specificity (true negative rate) describes the proportion of negatives that are correctly identified as such. However, new methods, such as quantitative polymerase chain reaction (qPCR), provide information beyond binary presence versus absence designations, which should be embraced in the evaluation of method performance. Although sensitivity and specificity are widely used to evaluate the performance of qPCR MST methods, it is difficult to reconcile results across studies because of a lack of standardization in performance testing protocols, namely the quantity of reference fecal material added to each qPCR reaction (denoted as “unit of measure” by Ervin et al.11) and the definition of the limit of detection (LOD) that has been used to classify a sample as being either positive or negative. First, it is well-documented that a larger quantity of fecal reference material added to qPCR test reaction mixtures can increase sensitivity estimates7 and decrease specificity estimates.11 This problem is further complicated because the unit © 2013 American Chemical Society

of measure used to assess sensitivity and specificity often varies between and within studies. In addition, the unit of measure has been expressed using various measurement scales such as DNA mass,6−8,12−15 wet or dry fecal mass,16−20 Escherichia coli or enterococci,21 or total Bacteroidales.5 Second, different definitions of LOD lead to dramatic shifts in sensitivity and specificity estimates.22,23 The primary challenge to determine the LOD of a qPCR assay is the stochastic effect observed at low copy numbers of a target DNA marker.24,25 In this study, we evaluate the performance of three MST qPCR assays using a range (6 orders of magnitude) of fecal reference material to illustrate the shortcomings of the traditional sensitivity−specificity approach and to introduce alternative metrics that circumvent the limitations described above. The requirement to test the reference fecal samples across a range of testing concentrations could be prohibitively expensive. Hence, we illustrate how an in silico calculation11,22,23 using data from a single test concentration can achieve the same objective.



MATERIALS AND METHODS MST Assays and Fecal Reference Sample Collection. Three human-associated assays (HumM2,13 HF183Taqman,6 and BsteriF16) with varying sensitivity and specificity22 were Received: Revised: Accepted: Published: 20

July 30, 2013 August 15, 2013 August 16, 2013 August 16, 2013 dx.doi.org/10.1021/ez400022t | Environ. Sci. Technol. Lett. 2014, 1, 20−25

Environmental Science & Technology Letters

Letter

the Cq value at each dilution level (Shapiro-Wilk normality test; p > 0.05) using the Cq_LOD value corresponding to the LOD [for example, 10 copies per reaction is equivalent to a Cq_LOD of 36.02 for the HF183Taqman assay (Figure 1, middle panel)].

tested following previously published protocols (lab 2 in ref 10) for five fecal reference sources (human, canine, cattle, chicken, and swine) that have been reported to react with humanassociated methods.22 DNA from individual animal fecal samples collected from different states was used to create a DNA composite for each source group: Oregon humans (n = 6), Florida (n = 9) and Wyoming (n = 10) canines, Nebraska cattle (n = 9), Kentucky chicken (n = 9), and Georgia swine (n = 9). DNA was extracted from each individual fecal sample using the DNA-EZ kit (GeneRite, North Brunswick, NJ), and DNA concentrations were quantified using the Quant-iT PicoGreen dsDNA kit (Invitrogen, Grand Island, NY). For each animal composite, an equal amount of DNA mass from each individual was combined, serially diluted, and used as qPCR template at concentrations ranging from 10 to 10−4 ng of DNA per reaction for human fecal extract and from 100 to 10−4 ng of DNA per reaction for other animals. The Oregon human composite did not have enough volume for all three assays, so a different human composite [California (n = 12), described by Boehm et al.10] was tested for BsteriF1 assay. No inhibition was observed for these composites according to the results from their serial dilutions.26 All process and no template controls were negative. qPCR Standard Curve. Decimal dilutions of an artificial plasmid standard27 ranging from 105 to 101 copies per reaction were used to generate calibration curves for the three assays. Six replicated qPCRs were tested at each dilution level. Standard curves were fit by weighted linear least-squares regression (WLS)28 to address the heteroscedascticity;24,25 i.e., reactions with smaller amount of plasmids have greater measurement variance in the quantification cycle, Cq.29 The weights used were the reciprocals of the variance, calculated at concentrations of 10 and 102 copies per reaction as well as the group variance for concentrations higher than 102 copies per reaction (calculated from the centralized Cq values, namely Cq values from which the mean Cq value of each concentration level was subtracted), because the F-test confirmed that the variance at concentrations higher than 102 copies per reaction was consistent but the variance at 10 and 102 copies per reaction was significantly larger (p < 0.05). All statistical analyses in this study were conducted using the base and stats packages in R (http://www.r-project.org/). In Silico Calculation. An in silico approach was developed using data from a single high-concentration sample that shows no signs of inhibition26,30 to predict results for other dilutions. The mean marker concentration in 10 (human) or 100 (other animals) ng of DNA per reaction was calculated from the standard curve using the mean Cq value of replicated qPCRs. The mean marker concentrations for other dilution levels were calculated in silico by dividing the mean concentration measured for 10 or 100 ng of DNA per reaction by appropriate dilution factors. The mean marker concentration was then translated to the mean Cq value by substitution into the standard curve. The standard deviation of this mean Cq value was extracted from the linear interpolation of the standard deviation estimated from the standard curve (standard deviation of HF183Taqman illustrated in Figure 1 of the Supporting Information). The 95% confidence interval (CI) of mean marker concentration at each dilution level was determined by the upper and lower boundaries of the Cq values (mean Cq ± 2 × standard deviation). The in silico probability that a dilution level will yield a measurement above the LOD was calculated according to the normal distribution of

Figure 1. Performance of the HF183Taqman assay. The middle panel illustrates the relationship between the target marker concentration and DNA mass per reaction established from qPCR (×) and in silico calculation (the solid line is the mean concentration, and dashed lines are the 95% confidence intervals of the DNA marker). The bottom and top panels show the probability that the observed concentration will be higher than the limit of detection (LOD) of 10 or 100 marker copies per reaction. The equivalent threshold (ET) for a source is the quantity of fecal material per reaction corresponding to the mean number of target marker copies per reaction (solid line in the middle panel) at the LOD, which also corresponds to the probability of 0.5 where half of the observed concentration is above the LOD. The false positive index (FPI) is the difference between the log10(ET) of the target and a nontarget source. It should be noted that the largest amount of fecal material tested is 100 ng of DNA per reaction following recommendations of the manufacturer of the qPCR instrument (Applied Biosystems).

This probability is the chance that the measured marker concentration equals ≥10 copies per reaction or Cq ≤ Cq_LOD. Thus, this probability can be calculated using the z-score in the standard normal distribution as p (Z ≤ z-score), where z-score is Cq_LOD − mean Cq standard deviation of Cq 21

dx.doi.org/10.1021/ez400022t | Environ. Sci. Technol. Lett. 2014, 1, 20−25

Environmental Science & Technology Letters



Letter

RESULTS AND DISCUSSION For both the target and nontarget sources, results from the in silico calculation at each dilution level matched the actual qPCR measurement on the corresponding dilution level. The 95% CI of the in silico calculation encompassed the actual qPCR measurements, even for the heteroscedascticity region (Figure 1, middle panel). The fraction of replicated qPCRs that will have a concentration higher than the LOD fell on the probability curve calculated from the in silico approach (Figure 1, top and bottom panels). Figure 1 illustrates that the probability that a measurement is above the LOD increases as the amount of DNA added to the qPCR reaction increases. The figure also illustrates that, given a fixed amount of DNA in a qPCR reaction, the probability that the measurement is above the LOD decreases as the LOD increases (Figure 1, top panel vs bottom panel). This illustrates why it is difficult to directly compare samples that are not tested with the same amount of DNA (or other unit of measure based on E. coli, enterococci, or Bacteroidales, for example) per reaction or the same LOD definition. In addition, Figure 1 shows that different sources actually have the same serial dilution pattern because qPCR detection is a function of the target marker copy numbers in the reaction. However, nontarget sources (chicken and cattle) shift to the right, relative to the target source (human), on the axis of DNA mass added to each reaction. This is because the prevalence of the qPCR marker in nontarget sources is typically lower than that in the target source; therefore, it takes more DNA from the nontarget sources to provide the same amount of target molecular marker to be detected by qPCR. In the study, each dilution level was assigned a probability that it would be detected above the LOD. This approach considers the “Poisson effect”-induced stochastic behavior observed at low copy numbers of a target DNA marker,24,25 namely, a large variance in Cq values around the LOD (Figure 1, middle panel). In fact, it is the measurement variance in the stochastic range that generates the S-shaped probability curve.24 The smaller the variance, the narrower the S-shaped curve will be (in Figure 1, the top panel with an LOD of 100 copies per reaction has a narrower S-shaped curve than the bottom panel with an LOD of 10 copies per reaction). The measurement variance has not been considered by other researchers from a review of the literature largely because the target source has typically been tested at a concentration above the stochastic range. However, many nontarget sources may yield results in the stochastic range where a sample has a probability between 0 and 100% to be above LOD. In the literature,7−9,21 such samples have been treated as having either 0 or 100% probability of being above the LOD, which biases the estimated specificity. Because of the influence that the quantity of fecal reference material tested in a qPCR reaction has on traditional specificity and sensitivity estimates, we propose two new performance metrics: the equivalent threshold (ET) and false positive index (FPI). ET is defined as the quantity of fecal reference material per reaction that corresponds to the LOD of an assay. Figure 1 indicates that the ET of a particular fecal source is the quantity of fecal material per reaction, for example, DNA mass per reaction that provides the amount of target marker (mean concentration, solid line in the middle panel) equal to the LOD, which also corresponds to the probability of 0.5 at which half of the observed concentration is above the LOD (top and bottom panels for two different LODs). Fecal sources with

varying prevalence of an assay’s target marker will have different ET values because it takes different amounts of fecal material to provide the amount of target marker per reaction at the LOD (Figure 1). Because the DNA mass per reaction (or other measurement scales) and the mean concentration of the target marker of a particular sample follow a linear relationship (Figure 1, solid lines in the middle panel), ET values can be calculated for different LOD values by scaling a previously calculated ET. For example, if an LOD of 100 marker copies per reaction is used (Figure 1, top panel), a new ET can be calculated by scaling the ET value based on an LOD of 10 marker copies per reaction by 10 (Table 1 and Figure 1, bottom panels). Given this relationship, the ET values of a fecal source from any two LOD definitions at the same measurement scale (e.g., DNA mass) can be readily converted from one to another. The ET value for a target source (e.g., human in this study) is closely related to what has been previously defined as sensitivity because ET value indicates the marker prevalence in the target source feces. The ET value, based on DNA mass, is similar for HF183Taqman and BsteriF1 for human fecal composites [although there are two different human composites, they had almost identical ET values for HF183Taqman (data not shown)], but ∼1 order of magnitude higher for HumM2 (Table 1 and Figure 2 of the Supporting Information). This is an indication that HF183 and BsteriF1 are ∼1 order of magnitude more sensitive for human feces than HumM2. When the ET value is defined on other measurement scales (fecal mass, FIB, Bacteroidales, etc.),11 the difference in ET values between two assays for the same source material will not vary from the difference when DNA mass is used as the scale,10 providing that DNA extraction does not preferentially yield more of the target marker for one assay versus another. Therefore, the relative ET of two assays is typically independent of measurement scale. The false positive index (FPI) of an assay is defined as the difference between the log10 transformed ET values of the target source and a nontarget source. It is closely related to what previously has been defined as specificity or cross reactivity. FPI is independent of the LOD definition (Table 1 and Figure 1). The FPI of chicken feces for the HF183Taqman assay (1.22 order of magnitude difference) and the FPI of canine feces for the BsteriF1 assay (0.35 and 1.30 order of magnitude differences for the two canine composites from different states) are 3.62

specificity has been reported to vary between groups of hosts5,13 or even between different subpopulations of the same animal host.7,31 Therefore, the conventional practice of combining representative samples from all nontarget sources to calculate one specificity value is problematic. We recommend only combining representative samples from closely related animal populations to calculate a corresponding range of ET and FPI.32 This will therefore allow comparisons of ET and FPI across animal groups from different geographic areas or from different times, as well as with different diets, ages, and health status.12 In summary, we introduce new performance metrics for quantitative MST methods designed to address current limitations in conventional sensitivity and specificity approaches. The conventional framework can be biased by the amount of reference fecal material tested in qPCR amplifications and LOD definition. We recommend providing the following information when evaluating assay performance: (1) Provide the measurement variance or standard deviation for Cq values generated from the standard curve experiments. (2) Test samples from a reference fecal pollution source by qPCR at a single high DNA concentration. Report the qPCR results in marker copies per reaction and DNA mass (other scales such as fecal mass, FIB, and Bacteroidales are useful as well) added to the reaction. Provision of this information could allow qPCR results from different studies to be converted and compared using the ET and FPI metrics outlined in this report assuming that potential differences in protocols across studies such as DNA extraction method, source of qPCR reagents, thermal cycling instrumentation, and the method for determining a Cq threshold do not bias the estimation of the mean marker concentration in a sample. This approach may also be used to compare different quantitative technology platforms such as digital PCR, microarray, and next-generation sequencing.

Note that the equivalent threshold is based on DNA mass per reaction. HF183Taqman amplified the human, chicken, and cattle composites. BsteriF1 amplified the human, canine, and swine composites. HumM2 amplified the human composite only when the total amount of DNA was ≤100 ng per reaction. bThe weighted-least-square (WLS) standard curve for each assay. cThe equivalent threshold is defined as the DNA mass per reaction that corresponds to the LOD. The ET value of human composites for the three assays indicates their relative sensitivity. BsteriF1 is slightly more sensitive than HF183, and both are ∼1 order of magnitude more sensitive than HumM2. dThe false positive index is defined as the difference in log10 transformed ET values of the target source (human) and a nontarget source. The FPI indicates the specificity of an assay. The larger the FPI is, the higher the specificity. The FPI is independent of the LOD. eNA indicates that there is no need to evaluate specificity for the target human source.

NAe >4.74 >4.74 1.30 0.35 NAe >4.42 NAe

1.22

>4.42

0.35 4.03 × 10−3 NAe 1.80 × 10−3 >5.42 >1.00 × 102 5.10 >1.00 × 102 1.22 6.38 × 10−2 NAe 3.83 × 10−3

4.03 × 10 1.80 × 10 >1.00 × 10 4.80 × 10 6.38 × 10 3.83 × 10

ET [DNA (ng/reaction)] for an RDL of 10 copies per reactionc FPI relative to humand ET [DNA (ng/reaction)] for an RDL of 100 copies per reactionc FPI relative to humand

Letter



ASSOCIATED CONTENT

S Supporting Information *

Figures S1 and S2. This material is available free of charge via the Internet at http://pubs.acs.org.



AUTHOR INFORMATION

Corresponding Author

*Department of Civil and Environmental Engineering, 473 Via Ortega, Room M08, Stanford University, Stanford, CA 94305. Telephone: (650) 725-3025. E-mail: [email protected]. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS We acknowledge Yiping Cao for providing comments on the manuscript. Funding for this project has been provided in part through an agreement with the California State Water Resources Control Board. H.C.G. was supported by the National Research Council. The contents of this document do not necessarily reflect the views and policies of the California State Water Resources Control Board and USEPA, nor does mention of trade names or commercial products constitute endorsement or recommendation for use.

a

>4.62 >1.00 × 102 NAe 2.38 × 10−2 >5.74 >1.00 × 102 5.18 >1.00 × 102

2.38 × 10 >1.00 × 10 2.75 × 10 3.60 × 10

1.30 3.60 × 10−2

>1.00 × 102

others

−3

human swine canine 2

−3 −4

canine 1

−4

human others human

−4

chicken

−3

cattle

2

others

2

HumM2 (slope, −3.48; intercept, 42.06; R2, 0.9966)b BsteriF1 (slope, −3.58; intercept, 40.42; R2, 0.9759)b HF183 (slope, −3.47; intercept, 39.49; R2, 0.9916)b

Table 1. Comparison of the Equivalent Threshold (ET) and False Positive Index (FPI) for Limits of Detection (LODs) of 10 and 100 Marker Copies per Reaction for Different Source Compositesa

Environmental Science & Technology Letters

23

dx.doi.org/10.1021/ez400022t | Environ. Sci. Technol. Lett. 2014, 1, 20−25

Environmental Science & Technology Letters



Letter

large alpine catchment area. Lett. Appl. Microbiol. 2007, 44 (4), 351− 356. (17) Reischer, G. H.; Kasper, D. C.; Steinborn, R.; Mach, R. L.; Farnleitner, A. H. Quantitative PCR method for sensitive detection of ruminant fecal pollution in freshwater and evaluation of this method in Alpine karstic regions. Appl. Environ. Microbiol. 2006, 72 (8), 5610− 5614. (18) Aslan, A.; Rose, J. B. Evaluation of the host specificity of Bacteroides thetaiotaomicron α-1-6, mannanase gene as a sewage marker. Lett. Appl. Microbiol. 2013, 56 (1), 51−56. (19) Bernhard, A. E.; Field, K. G. A PCR assay to discriminate human and ruminant feces on the basis of host differences in BacteriodesPrevotella genes encoding 16S rRNA. Appl. Environ. Microbiol. 2000, 66 (10), 4571−4574. (20) Lee, C. S.; Lee, J. Evaluation of new gyrB-based real-time PCR system for the detection of B. fragilis as an indicator of human-specific fecal contamination. J. Microbiol. Methods 2010, 82 (3), 311−318. (21) Green, H. C.; Dick, L. K.; Gilpin, B.; Samadpour, M.; Field, K. G. Genetic markers for rapid PCR-based identification of gull, Canada goose, duck, and chicken fecal contamination in water. Appl. Environ. Microbiol. 2012, 78 (2), 503−510. (22) Layton, B. A.; Cao, Y.; Ebentier, D. L.; Hanley, K. T.; Ballesté, E.; Brandão, J.; Byappanahalli, M. N.; Converse, R.; Farnleitner, A. H.; Gentry-Shields, J.; Gidley, M. L.; Gourmelon, M.; Lee, C. S.; Lee, J.; Lozach, S.; Madi, T.; Meijer, W.; Noble, R. T.; Peed, L.; Reischer, G. H.; Rodrigues, R.; Rose, J.; Schriewer, A.; Sinigalliano, C.; Srinivasan, S.; Stewart, J. R.; Van De Werfhorst, L. C.; Wang, D.; Whitman, R. L.; Wuertz, S.; Jay, J. A.; Holden, P. A.; Boehm, A. B.; Shanks, O. C.; Griffith, J. Performance of human fecal anaerobe-associated PCRbased assays in a multi-laboratory method evaluation study. Water Res. 2013, 10.1016/j.watres.2013.05.060. (23) Raith, M.; Kelty, C. A.; Griffith, J. F.; Schriewer, A.; Wuertz, S.; Mieszkin, S.; Gourmelon, M.; Reischer, G.; Farnleitner, A. H.; Ervin, J.; Holden, P. A.; Ebentier, D. L.; Jay, J. A.; Wang, D.; Boehm, A. B.; Aw, T.; Rose, J. B.; Balleste, E.; Meijer, W. G.; Sivaganesan, M.; Shanks, O. C. Comparison of PCR and quantitative real-time PCR methods for the characterization of ruminant and cattle fecal pollution sources. Water Res. 2013, 10.1016/j.watres.2013.03.061. (24) Final report and executive summaries from the AOAC international presidential task force on best practices in microbiological methodology. Report 223-01-2464; AOAC International: Gaithersburg, MD, 2006. (25) Walsh, P. S.; Erlich, H. A.; Higuchi, R. Preferential PCR amplification of alleles: Mechanisms and solutions. PCR Methods Appl. 1992, 1 (4), 241−250. (26) Rajal, V. B.; McSwain, B. S.; Thompson, D. E.; Leutenegger, C. M.; Kildare, B. J.; Wuertz, S. Validation of hollow fiber ultrafiltration and real-time PCR using bacteriophage PP7 as surrogate for the quantification of viruses from water samples. Water Res. 2007, 41 (7), 1411−1422. (27) Wang, D.; Farnleitner, A. H.; Field, K. G.; Green, H. C.; Shanks, O. C.; Boehm, A. B. Enterococcus and Escherichia coli fecal source apportionment with microbial source tracking genetic markers: Is it feasible? Water Res. 2013, 10.1016/j.watres.2013.02.058. (28) Strutz, T. In Data Fitting and Uncertainty: A Practical Introduction to Weighted Least Squares and Beyond; Vieweg+Teubner, 2010; Chapter 3. (29) Bustin, S. A.; Benes, V.; Garson, J. A.; Hellemans, J.; Huggett, J.; Kubista, M.; Mueller, R.; Nolan, T.; Pfaffl, M. W.; Shipley, G. L. The MIQE guidelines: Minimum information for publication of quantitative real-time PCR experiments. Clin. Chem. 2009, 55 (4), 611−622. (30) Haugland, R. A.; Siefring, S.; Lavender, J.; Varma, M. Influences of sample interference and interference controls on quantification of enterococci fecal indicator bacteria in surface water samples by the qPCR method. Water Res. 2012, 46 (18), 5989−6001. (31) Shanks, O. C.; Kelty, C. A.; Archibeque, S.; Jenkins, M.; Newton, R. J.; McLellan, S. L.; Huse, S. M.; Sogin, M. L. Community

REFERENCES

(1) National summary of impaired waters and TMDL information; U.S. Environmenal Protection Agency: Washington, DC, 2011. (2) Guidelines for safe recreational water environments, Volume 1: Coastal and fresh waters; World Health Organization: Geneva, 2003. (3) Byappanahalli, M. N.; Nevers, M. B.; Korajkic, A.; Staley, Z. R.; Harwood, V. J. Enterococci in the Environment. Microbiol. Mol. Biol. Rev. 2012, 76 (4), 685−706. (4) Soller, J. A.; Schoen, M. E.; Bartrand, T.; Ravenscroft, J. E.; Ashbolt, N. J. Estimated human health risks from exposure to recreational waters impacted by human and non-human sources of faecal contamination. Water Res. 2010, 44 (16), 4674−4691. (5) Kildare, B. J.; Leutenegger, C. M.; McSwain, B. S.; Bambic, D. G.; Rajal, V. B.; Wuertz, S. 16S rRNA-based assays for quantitative detection of universal, human-, cow-, and dog-specific fecal Bacteroidales: A Bayesian approach. Water Res. 2007, 41 (16), 3701−3715. (6) Haugland, R. A.; Varma, M.; Sivaganesan, M.; Kelty, C.; Peed, L.; Shanks, O. C. Evaluation of genetic markers from the 16S rRNA gene V2 region for use in quantitative detection of selected Bacteroidales species and human fecal waste by qPCR. Syst. Appl. Microbiol. 2010, 33 (6), 348−357. (7) Shanks, O. C.; White, K.; Kelty, C. A.; Hayes, S.; Sivaganesan, M.; Jenkins, M.; Varma, M.; Haugland, R. A. Performance Assessment PCR-Based Assays Targeting Bacteroidales Genetic Markers of Bovine Fecal Pollution. Appl. Environ. Microbiol. 2010, 76 (5), 1359−1366. (8) Shanks, O. C.; White, K.; Kelty, C. A.; Sivaganesan, M.; Blannon, J.; Meckes, M.; Varma, M.; Haugland, R. A. Performance of PCRbased assays targeting Bacteroidales genetic markers of human fecal pollution in sewage and fecal samples. Environ. Sci. Technol. 2010, 44 (16), 6281−6288. (9) Van De Werfhorst, L. C.; Sercu, B.; Holden, P. A. Comparison of the Host Specificities of Two Bacteroidales Quantitative PCR Assays Used for Tracking Human Fecal Contamination. Appl. Environ. Microbiol. 2011, 77 (17), 6258−6260. (10) Boehm, A. B.; Van De Werfhorst, L. C.; Griffith, J. F.; Holden, P. A.; Jay, J. A.; Shanks, O. C.; Wang, D.; Weisberg, S. B. Performance of forty-one microbial source tracking methods: A twenty-seven lab evaluation study. Water Res. 2013, 10.1016/j.watres.2012.12.046. (11) Ervin, J. S.; Russell, T. L.; Layton, B. A.; Yamahara, K. M.; Wang, D.; Sassoubre, L. M.; Cao, Y.; Kelty, C. A.; Sivaganesan, M.; Boehm, A. B.; Holden, P.; Weisberg, S. B.; Shanks, O. C. Characterization of fecal concentrations in human and other animal sources by physical, cutlure-based, and quantitative real-time PCR methods. Water Res. 2013, 10.1016/j.watres.2013.02.060. (12) Shanks, O. C.; Atikovic, E.; Blackwood, A. D.; Lu, J. R.; Noble, R. T.; Domingo, J. S.; Seifring, S.; Sivaganesan, M.; Haugland, R. A. Quantitative PCR for detection and enumeration of genetic markers of bovine fecal pollution. Appl. Environ. Microbiol. 2008, 74 (3), 745− 752. (13) Shanks, O. C.; Kelty, C. A.; Sivaganesan, M.; Varma, M.; Haugland, R. A. Quantitative PCR for genetic markers of human fecal pollution. Appl. Environ. Microbiol. 2009, 75 (17), 5507−5513. (14) Yampara-Iquise, H.; Zheng, G.; Jones, J. E.; Carson, C. A. Use of a Bacteroides thetaiotaomicron-specific α-1-6, mannanase quantitative PCR to detect human faecal pollution in water. J. Appl. Microbiol. 2008, 105 (5), 1686−1693. (15) Reischer, G.; Ebdon, J. E.; Bauer, J.; Schuster, N.; Ahmed, W.; Åström, J.; Blanch, A.; Blöschl, G.; Byamukama, D.; Coakley, T.; Ferguson, C.; Goshu, G.; Ko, G.; De Roda Husman, A. M.; Mushi, D.; Poma, R.; Pradhan, B.; Rajal, V.; Schade, M.; Sommer, R.; Taylor, H. D.; Toth, E.; Vrajmasu, V.; Wuertz, S.; Mach, R. L.; Farnleitner, A. H. Performance characteristics of qPCR assays targeting human- and ruminant-associated Bacteroidetes for microbial source tracking across sixteen countries on six continents. Environ. Sci. Technol. 2013, 47, 8548−8556. (16) Reischer, G. H.; Kasper, D. C.; Steinborn, R.; Farnleitner, A. H.; Mach, R. L. A quantitative real-time PCR assay for the highly sensitive and specific detection of human faecal influence in spring water from a 24

dx.doi.org/10.1021/ez400022t | Environ. Sci. Technol. Lett. 2014, 1, 20−25

Environmental Science & Technology Letters

Letter

structures of fecal bacteria in cattle from different animal feeding operations. Appl. Environ. Microbiol. 2011, 77 (9), 2992−3001. (32) Wang, D.; Silkie, S. S.; Nelson, K. L.; Wuertz, S. Estimating true human and animal host source contribution in quantitative microbial source tracking using the Monte Carlo method. Water Res. 2010, 44 (16), 4760−4775.

25

dx.doi.org/10.1021/ez400022t | Environ. Sci. Technol. Lett. 2014, 1, 20−25