Co-quantification of crAssphage increases confidence in wastewater-based epidemiology for SARS-CoV-2 in low prevalence areas.

Wastewater surveillance of SARS-CoV-2 RNA is increasingly being incorporated into public health efforts to respond to the COVID-19 pandemic. In order to obtain the maximum benefit from these efforts, approaches to wastewater monitoring need to be rapid, sensitive, and relatable to relevant epidemiological parameters. In this study, we present an ultracentrifugation-based method for the concentration of SARS-CoV-2 wastewater RNA and use crAssphage, a bacteriophage specific to the human gut, to help account for RNA loss during transit in the wastewater system and sample processing. With these methods, we were able to detect, and sometimes quantify, SARS-CoV-2 RNA from 20 mL wastewater samples within as little as 4.5 hours. Using known concentrations of bovine coronavirus RNA and deactivated SARS-CoV-2, we estimate recovery rates of approximately 7-12% of viral RNA using our method. Results from 24 sewersheds across Upstate New York during the spring and summer of 2020 suggested that stronger signals of SARS-CoV-2 RNA from wastewater may be indicative of greater COVID-19 incidence in the represented service area approximately one week in advance. SARS-CoV-2 wastewater RNA was quantifiable in some service areas with daily positives tests of less than 1 per 10,000 people or when weekly positive test rates within a sewershed were as low as 1.7%. crAssphage DNA concentrations were significantly lower during periods of high flow in almost all areas studied. After accounting for flow rate and population served, crAssphage levels per capita were estimated to be about 1.35 × 1011 and 2.42 × 108 genome copies per day for DNA and RNA, respectively. A negative relationship between per capita crAssphage RNA and service area size was also observed likely reflecting degradation of RNA over long transit times. Our results reinforce the potential for wastewater surveillance to be used as a tool to supplement understanding of infectious disease transmission obtained by traditional testing and highlight the potential for crAssphage co-detection to improve interpretations of wastewater surveillance data.

High-throughput wastewater analysis for substance use assessment in central New York during the COVID-19 pandemic.

Wastewater entering sewer networks represents a unique source of pooled epidemiological information. In this study, we coupled online solid-phase extraction with liquid chromatography-high resolution mass spectrometry to achieve high-throughput analysis of health and lifestyle-related substances in untreated municipal wastewater during the coronavirus disease 2019 (COVID-19) pandemic. Twenty-six substances were identified and quantified in influent samples collected from six wastewater treatment plants during the COVID-19 pandemic in central New York. Over a 12 week sampling period, the mean summed consumption rate of six major substance groups (i.e., antidepressants, antiepileptics, antihistamines, antihypertensives, synthetic opioids, and central nervous system stimulants) correlated with disparities in household income, marital status, and age of the contributing populations as well as the detection frequency of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA in wastewater and the COVID-19 test positivity in the studied sewersheds. Nontarget screening revealed the covariation of piperine, a nontarget substance, with SARS-CoV-2 RNA in wastewater collected from one of the sewersheds. Overall, this proof-of-the-concept study demonstrated the utility of high-throughput wastewater analysis for assessing the population-level substance use patterns during a public health crisis such as COVID-19.

Development and validation of rapid environmental DNA (eDNA) detection methods for bog turtle (Glyptemys muhlenbergii).

Bog turtles (Glyptemys muhlenbergii) are listed as Species of Greatest Conservation Need (SGCN) for wildlife action plans in every state it occurs and multi-state efforts are underway to better characterize extant populations and prioritize restoration efforts. However, traditional sampling methods can be ineffective due to the turtle's wetland habitat, small size, and burrowing nature. Molecular methods, such as qPCR, provide the ability to overcome this challenge by effectively quantifying minute amounts of turtle DNA left behind in its environment (eDNA). Developing such methods for bog turtles has proved difficult partly because of the high sequence similarity between bog turtles and closely-related, cohabitating species, most often wood turtles (Glyptemys insculpta). Additionally, substrates containing bog turtle eDNA are often rich in organics or other substances that frequently inhibit both DNA extraction and qPCR amplification. Here, we describe the development and validation of a qPCR assay, BT3, targeting the mitochondrial cytochrome oxidase I gene that correctly identifies bog turtles with 100% specificity and sensitivity when tested on 201 blood samples collected from six species over a wide geographic range. We also developed a full-process internal control employing a genetically modified strain of Caenorhabditis elegans to improve DNA extraction methods, limit false negative results due to qPCR inhibition, and measure total DNA recovery from each sample. Using the internal control, we found that DNA recovery varied by over an order of magnitude between samples and likely explains the lack of bog turtle detection in some cases. Methods presented herein are highly-specific and may offer a more cost effective, non-invasive tool to supplement bog turtle population assessments in the Eastern United States. Poor or differential DNA recovery, which remains unmeasured in the vast majority of eDNA studies, significantly reduced the ability to detect bog turtle in their natural environment.

Identification of Specialists and Abundance-Occupancy Relationships among Intestinal Bacteria of Aves, Mammalia, and Actinopterygii.

The coalescence of next-generation DNA sequencing methods, ecological perspectives, and bioinformatics analysis tools is rapidly advancing our understanding of the evolution and function of vertebrate-associated bacterial communities. Delineation of host-microbe associations has applied benefits ranging from clinical treatments to protecting our natural waters. Microbial communities follow some broad-scale patterns observed for macroorganisms, but it remains unclear how the specialization of intestinal vertebrate-associated communities to a particular host environment influences broad-scale patterns in microbial abundance and distribution. We analyzed the V6 region of 16S rRNA genes amplified from 106 fecal samples spanning Aves, Mammalia, and Actinopterygii (ray-finned fish). We investigated the interspecific abundance-occupancy relationship, where widespread taxa tend to be more abundant than narrowly distributed taxa, among operational taxonomic units (OTUs) within and among host species. In a separate analysis, we identified specialist OTUs that were highly abundant in a single host and rare in all other hosts by using a multinomial model without excluding undersampled OTUs a priori. We show that intestinal microbes in humans and other vertebrates display abundance-occupancy relationships, but because intestinal host-associated communities have undergone intense specialization, this trend is violated by a disproportionately large number of specialist taxa. Although it is difficult to distinguish the effects of dispersal limitations, host selection, historical contingency, and stochastic processes on community assembly, results suggest that intestinal bacteria can be shared among diverse hosts in ways that resemble the distribution of "free-living" bacteria in the extraintestinal environment.

Comparison of Sewage and Animal Fecal Microbiomes by Using Oligotyping Reveals Potential Human Fecal Indicators in Multiple Taxonomic Groups.

Most DNA-based microbial source tracking (MST) approaches target host-associated organisms within the order Bacteroidales, but the gut microbiota of humans and other animals contain organisms from an array of other taxonomic groups that might provide indicators of fecal pollution sources. To discern between human and nonhuman fecal sources, we compared the V6 regions of the 16S rRNA genes detected in fecal samples from six animal hosts to those found in sewage (as a proxy for humans). We focused on 10 abundant genera and used oligotyping, which can detect subtle differences between rRNA gene sequences from ecologically distinct organisms. Our analysis showed clear patterns of differential oligotype distributions between sewage and animal samples. Over 100 oligotypes of human origin occurred preferentially in sewage samples, and 99 human oligotypes were sewage specific. Sequences represented by the sewage-specific oligotypes can be used individually for development of PCR-based assays or together with the oligotypes preferentially associated with sewage to implement a signature-based approach. Analysis of sewage from Spain and Brazil showed that the sewage-specific oligotypes identified in U.S. sewage have the potential to be used as global alternative indicators of human fecal pollution. Environmental samples with evidence of prior human fecal contamination had consistent ratios of sewage signature oligotypes that corresponded to the trends observed for sewage. Our methodology represents a promising approach to identifying new bacterial taxa for MST applications and further highlights the potential of the family Lachnospiraceae to provide human-specific markers. In addition to source tracking applications, the patterns of the fine-scale population structure within fecal taxa suggest a fundamental relationship between bacteria and their hosts.

Development of rapid canine fecal source identification PCR-based assays.

The extent to which dogs contribute to aquatic fecal contamination is unknown despite the potential for zoonotic transfer of harmful human pathogens. We used genome fragment enrichment (GFE) to identify novel nonribosomal microbial genetic markers potentially useful for detecting dog fecal contamination with PCR-based methods in environmental samples. Of the 679 sequences obtained from GFE, we used 84 for the development of PCR assays targeting putative canine-associated genetic markers. Twelve genetic markers were shown to be prevalent among dog fecal samples and were rarely found in other animals. Three assays, DG3, DG37, and DG72, performed best in terms of specificity and sensitivity and were used for the development of SYBR Green and TaqMan quantitative PCR (qPCR) assays. qPCR analysis of 244 fecal samples collected from a wide geographic range indicated that marker concentrations were below limits of detection in noncanine hosts. As a proof-of-concept, these markers were detected in urban stormwater samples, suggesting a future application of newly developed methods for water quality monitoring.

Improved HF183 quantitative real-time PCR assay for characterization of human fecal pollution in ambient surface water samples.

Quantitative real-time PCR (qPCR) assays that target the human-associated HF183 bacterial cluster within members of the genus Bacteroides are among the most widely used methods for the characterization of human fecal pollution in ambient surface waters. In this study, we show that a current TaqMan HF183 qPCR assay (HF183/BFDrev) routinely forms nonspecific amplification products and introduce a modified TaqMan assay (HF183/BacR287) that alleviates this problem. The performance of each qPCR assay was compared in head-to-head experiments investigating limits of detection, analytical precision, predicted hybridization to 16S rRNA gene sequences from a reference database, and relative marker concentrations in fecal and sewage samples. The performance of the modified HF183/BacR287 assay is equal to or improves upon that of the original HF183/BFDrev assay. In addition, a qPCR chemistry designed to combat amplification inhibition and a multiplexed internal amplification control are included. In light of the expanding use of PCR-based methods that rely on the detection of extremely low concentrations of DNA template, such as qPCR and digital PCR, the new TaqMan HF183/BacR287 assay should provide more accurate estimations of human-derived fecal contaminants in ambient surface waters.

Enterococcus and Escherichia coli fecal source apportionment with microbial source tracking genetic markers--is it feasible?

Fecal pollution is measured in surface waters using culture-based measurements of enterococci and Escherichia coli bacteria. Source apportionment of these two fecal indicator bacteria is an urgent need for prioritizing remediation efforts and quantifying health risks associated with source-specific pathogens. There are a number of quantitative real-time PCR (QPCR) assays that estimate concentrations of source-associated genetic markers; however, their concentrations are not necessarily amenable to source apportionment because the markers may differ in prevalence across sources. Here we mathematically derive and test, under ideal conditions, a method that utilizes the ratios of fecal source-associated genetic markers and culture and molecular measurements of general fecal indicators to apportion enterococci and E. coli. The source contribution is approximately equal to the ratio of the source-associated and the general fecal indicator concentrations in a water sample divided by their ratio in the source material, so long as cross-reactivity is negligible. We illustrate the utility of the ratio method using samples consisting of mixtures of various fecal pollution sources. The results from the ratio method correlated well with the actual source apportionment in artificial samples. However, aging of contamination can confound source allocation predictions. In particular, culturable enterococci and E. coli, the organisms presently regulated in the United States and much of the world, decay at different rates compared to source-associated markers and as a result cannot be apportioned using this method. However, limited data suggest a similar decay rate between source-associated and QPCR-measured Enterococcus and E. coli genetic markers, indicating that apportionment may be possible for these organisms; however further work is needed to confirm.

Sensitive detection of sample interference in environmental qPCR.

Sample interference in environmental applications of quantitative PCR (qPCR) can prevent accurate estimations of molecular markers in the environment. We developed a spike-and-recovery approach using a mutant strain of Escherichia coli that contains a chromosomal insertion of a mutant GFP gene. The method was tested in water samples by separately reducing extraction efficiency or adding humic acids and ethanol, compounds that often contaminate environmental DNA extracts, and analyzing qPCR amplification of the spiked E. coli control and human fecal Bacteroides markers (HF183 and HF134). This approach, coupled with previously developed kinetic outlier detection (KOD) methods, allowed sensitive detection of PCR inhibition at much lower inhibitor concentrations than alternative approaches using Cq values or amplification efficiencies. Although HF183 was more sensitive to the effects of qPCR inhibitors than the E. coli control assay, KOD methods correctly identified inhibition of both control and HF183 assays in samples containing as little as 0.1 ng humic acids per reaction or 5% ethanol. Because sigmoidal modeling methods allow distinction of qPCR inhibition from poor DNA recovery, we were able to simultaneously identify qPCR-inhibited reactions and estimate recovery of nucleic acids in environmental samples using a single control assay. Since qPCR is currently used to estimate important water quality parameters that have serious economic and human health outcomes, these results are timely. While we demonstrate the methods in the context of water quality regulation, they will be useful in all areas of environmental research that use qPCR.

Genetic markers for rapid PCR-based identification of gull, Canada goose, duck, and chicken fecal contamination in water.

Avian feces contaminate waterways but contribute fewer human pathogens than human sources. Rapid identification and quantification of avian contamination would therefore be useful to prevent overestimation of human health risk. We used subtractive hybridization of PCR-amplified gull fecal 16S RNA genes to identify avian-specific fecal rRNA gene sequences. The subtracters were rRNA genes amplified from human, dog, cat, cow, and pig feces. Recovered sequences were related to Enterobacteriaceae (47%), Helicobacter (26%), Catellicoccus (11%), Fusobacterium (11%), and Campylobacter (5%). Three PCR assays, designated GFB, GFC, and GFD, were based on recovered sequence fragments. Quantitative PCR assays for GFC and GFD were developed using SYBR green. GFC detected down to 0.1 mg gull feces/100 ml (corresponding to 2 gull enterococci most probable number [MPN]/100 ml). GFD detected down to 0.1 mg chicken feces/100 ml (corresponding to 13 Escherichia coli MPN/100 ml). GFB and GFC were 97% and 94% specific to gulls, respectively. GFC cross-reacted with 35% of sheep samples but occurred at about 100,000 times lower concentrations in sheep. GFD was 100% avian specific and occurred in gulls, geese, chickens, and ducks. In the United States, Canada, and New Zealand, the three markers differed in their geographic distributions but were found across the range tested. These assays detected four important bird groups contributing to fecal contamination of waterways: gulls, geese, ducks, and chickens. Marker distributions across North America and in New Zealand suggest that they will have broad applicability in other parts of the world as well.

Differential decay of human faecal Bacteroides in marine and freshwater.

Genetic markers from Bacteroides and other faecal bacteria are being tested for inclusion in regulations to quantify aquatic faecal contamination and estimate public health risk. For the method to be used quantitatively across environments, persistence and decay of markers must be understood. We measured concentrations of contaminant molecular markers targeting Enterococcus and Bacteroides spp. in marine and freshwater microcosms spiked with human sewage and exposed to either sunlight or dark treatments. We used Bayesian statistics with a delayed Chick-Watson model to estimate kinetic parameters for target decay. DNA- and RNA-based targets decayed at approximately the same rate. Molecular markers persisted (could be detected) longer in marine water. Sunlight increased the decay rates of cultured indicators more than those of molecular markers; sunlight also limited persistence of molecular markers. Within each treatment, Bacteroides markers had similar decay profiles, but some Bacteroides markers significantly differed in decay rates. The role of extracellular DNA in persistence appeared unimportant in the microcosms. Because conditions were controlled, microcosms allowed the effects of specific environmental variables on marker persistence and decay to be measured. While marker decay profiles in more complex environments would be expected to vary from those observed here, the differences we measured suggest that water matrix is an important factor affecting quantitative source tracking and microbial risk assessment applications.