Critical Evaluation of CrAssphage as a Molecular Marker for Human-Derived Wastewater Contamination in the Aquatic Environment.

The discharge of human-derived wastewater represents a major threat to water quality with the potential for waterborne disease outbreaks mainly associated with enteric viruses. To prevent illnesses, indicators associated with fecal contamination are monitored in polluted areas, however, their prevalence often does not correlate well with viral pathogens. In this study, we used crAssphage, a recently discovered human-specific gut-associated bacteriophage, for the surveillance of wastewater-derived viral contamination. Untreated and treated wastewater, surface water, sediment and mussel samples were collected monthly over 1 year from the Conwy River and estuary (UK) and were analyzed for crAssphage marker by quantitative PCR. This is the first long-term catchment-to-coast scale study of environmental crAssphage concentrations. CrAssphage was detected in all sample types and showed no distinct seasonal pattern. CrAssphage concentrations were 2 × 105-109 genome copies (gc)/L in all untreated wastewater influent and 107-108 gc/L in secondary treated effluent samples, 3 × 103 gc/L-3 × 107 gc/L in surface water samples (94% positive) and 2 × 102-104 gc/g sediment (68% positive) and mussel digestive tissue (79% positive). CrAssphage concentrations were 1-5 log10 higher than human enteric virus titers (norovirus, sapovirus, adenovirus, polyomavirus). Our results indicate that crAssphage is well suited to tracking human wastewater contamination and pollution risk assessment in aquatic environments.

Introduction

Enteric viruses are the most common etiologic agents of gastroenteritis globally. They are discharged with treated and untreated wastewater into the aquatic environment where they are able to persist for long periods of time (Kotwal and Cannon 2014). Hence, these viruses are often responsible for waterborne and foodborne illnesses due to the use of contaminated recreational water, and the consumption of polluted water and shellfish (Radin 2014; Rodríguez-Lázaro et al. 2012). Over 150 human pathogenic viruses, including noroviruses (NoV), sapoviruses (SaV), rotaviruses, hepatitis A/E viruses, adenoviruses (AdV), enteroviruses and polyomaviruses, have been identified in watercourses (Rodríguez-Lázaro et al. 2012; Tran et al. 2015). As the detection and surveillance of all pathogenic viral strains is not feasible, indicators are often used for tracking wastewater contamination in the environment (Symonds and Breitbart 2015).

The crAssphage sensu stricto genome (RefSeq accession number NC_024711), hereafter referred to as crAssphage, was first assembled from human fecal microbiomes within the last 5 years (Dutilh et al. 2014). Further analysis of human metagenomes showed that crAssphage is a highly abundant virus in the human gut, especially in individuals living in industrialized areas (Honap et al. 2018; Stachler and Bibby 2014). Data mining of virome datasets subsequently revealed that there is a diverse group of crAss-related phages present in the human gut, potentially representing a new family of viruses (Yutin et al. 2017). The first isolated representative of this family has just been discovered, phiCrAss001, infecting the human gut bacterium Bacteroides intestinalis (Shkoporov et al. 2018). Following the discovery of the prevalence of crAss-like phages in the human gut, their usefulness as a tool for source tracking of human fecal pollution was recognized (Stachler and Bibby 2014). Quantitative PCR (qPCR) assays have been developed for the quantification of the crAssphage in stool samples (Cinek et al. 2018; Liang et al. 2018; Stachler et al. 2017). Assay specificity showed great variation, suggesting that some regions of the crAssphage genome may be similar to other bacteriophage genome segments, whereas other regions are indeed human-specific (Stachler et al. 2018). CrAssphage prevalence has been investigated in the aquatic environment and found at high concentrations in untreated wastewater, wastewater-contaminated stream and in stormwater (Ahmed et al. 2018; García-Aljaro et al. 2017; Stachler et al. 2018). The concentration correlated with rainfall frequency in a polluted stream in a 30-day monitoring trial (Stachler et al. 2018). Nonetheless, crAssphage has not been used for long-term monitoring of wastewater-contaminated environments.

The aim of the research described here was to assess the usefulness of a crAssphage marker specifically as an indicator for wastewater-associated viral contamination in the aquatic environment. CrAssphage marker concentrations were monitored in treated and untreated wastewater, riverine and estuarine water, and sediment and shellfish samples collected monthly over 1 year along the Conwy River and estuary, UK. CrAssphage marker concentrations were compared with common enteric viral contaminants, namely NoV genogroups GI and GII (NoVGI and NoVGII), SaV genogroup GI, AdV and JC polyomavirus (JCV).

Materials and Methods

Sample Types and Process

We investigated the presence and concentration of crAssphage in nucleic acid extracted from concentrated wastewater, surface water, sediment and shellfish (blue mussel, Mytilus edulis) samples that were known to be contaminated with enteric viruses. The sampling schedule, locations and sample processing are detailed in Farkas et al (2018a) and summarized in Table 1. Thirteen samples were taken of each sample type once per month over a 1-year cycle (August 2016–August 2017).

Table 1

Percentile of virus-positive samples in wastewater (WW), surface water (SW), sediment (Sed) and mussel (SF) samples collected in the Conwy catchment and estuary between August 2016 and August 2017

Sample code	Site	Sample type	NoVGI (%)	NoVGII (%)	SaV (%)	AdV (%)	JCV (%)	CrAssphage (%)
GI	Ganol	WW influent	62	69	38	92	100	100
BI	Betws-y-Coed	WW influent	31	62	15	85	85	100a
BE	Betws-y-Coed	WW effluent (AC)	31	38	0	77	85	100a
LI	Llanrwst	WW influent	15	54	15	85	85	100a
LE	Llanrwst	WW effluent (BF)	31	62	15	92	100	100
TI	Tal-y-Bont	WW influent	31	62	38	100	92	100a
TE	Tal-y-Bont	WW effluent (AC)	45	46	15	92	85	100a
SW1	Betws-y-Coed–Llugwy River	River water	15	23	8	85	54	92
SW2	Betws-y-Coed–Conwy River	River water	0	8	8	85	54	92a
SW3	Llanrwst	Estuarine water	15	8	0	92	92	92
SW4	Conwy	Estuarine water	8	15	0	92	62	92
Sed1	Deganwy beach	Estuarine sediment	31	15	0	77	46	62
Sed2	Morfa beach	Estuarine sediment	8	0	0	85	46	75a
Sed4	Conwy	Estuarine sediment	15	8	0	85	23	62
SF1	Deganwy beach	Mussel	15	15	0	85	31	77
SF2	Morfa beach	Mussel	23	8	0	85	38	83a

AC activated sludge treatment, BF biofilter; n = 13, an = 12. NoV GI and GII Norovirus GI and GII, SaV sapovirus, AdV adenovirus and JCV JC polyomavirus data were adopted from Farkas et al. (2018a)

The water and wastewater samples, 10 L and 1 L, respectively, were concentrated using tangential flow ultrafiltration followed by beef extract elution and polyethylene glycol (PEG) precipitation (Farkas et al. 2018). Viruses were eluted from 10 g of sediment using beef extract solution and precipitated with PEG (Farkas et al. 2017). The mussel digestive tissue samples, 2 g each, were processed as described in the ISO-TS 15216 standard (International Organizatoin for Standardization 2013). Nucleic acids were extracted from the concentrates using the MiniMag NucliSens Nucleic Acid Extraction System (BioMérieux SA, France). Assay efficiency was assessed using mengovirus as a process control, which gave > 10% recovery values (Farkas et al. 2018).

CrAssphage Quantification

CrAssphage DNA was quantified using the CPQ_056 TaqMan primers and probe, which have shown great accuracy and little cross-reaction with animal litter (Ahmed et al. 2018; Stachler et al. 2017). The qPCR assays were carried out in a QuantStudio® Flex 6 Real-Time PCR System (Applied Biosystems, USA). The 20 µL reaction mix contained 1x KAPA Probe Force qPCR mix (KAPA Biosystems, USA) with 10 pmol of the forward, 10 pmol of the reverse primers, 5 pmol of the probe, 1 µg bovine serum albumin, 50 nM ROX reference dye. The sample volume was 2 µL for wastewater influent, 4 µL for wastewater effluent and surface water, and 8 µL for sediment and mussel samples in each reaction. For quantification, dilution series of a plasmid DNA carrying the target sequence were used. Non-template controls (molecular-grade water) were added to each reaction plate. Amplification was carried out using the following thermal cycling conditions: 98 °C for 5 min, then 40 cycles of 95 °C for 15 s, 60 °C for 1 min. The assay efficiency was 90–110%.

Data Analysis

The viral concentrations were expressed as gc/L of wastewater and surface water or gc/g of sediment or mussel digestive tissue. Linear regressions and Spearman’s rank correlation coefficients (r) were calculated between viral concentrations using two-tailed 95% confidence intervals in SigmaPlot 13.0 (Systat Software Inc., US).

Results and Discussion

CrAssphage qPCR

The main goal of this study was to evaluate the usefulness of a novel human-associated phage, the crAssphage, as a viral indicator for wastewater contamination. We, therefore, analyzed crAssphage concentrations in wastewater influent and effluent samples released into the river, and in river and estuarine water, sediment and mussel samples contaminated with wastewater, and correlated the concentrations of the crAssphage with enteric viral titres in the samples using the qPCR results. In our study, the non-template controls were negative in each assay suggesting no cross-contamination. The LOD was 2 gc/reaction and samples below that concentration were considered negative. The LOQ was 20 gc/reaction. Dilution of samples did not affect qPCR results indicating no inhibition in the wastewater effluent, sediment or mussel matrices. We utilized a well-established primer and probe set that has previously been used for the quantification of crAssphage in stool, wastewater, storm drain outfall and surface water (Ahmed et al. 2018; Stachler et al. 2017, 2018). The qPCR assay has previously shown good human specificity, with little cross-reaction with dog, gull and poultry litter (Ahmed et al. 2018; Stachler et al. 2017). Overall, the qPCR assay used for crAssphage quantification appeared to be accurate and sensitive and hence suitable for the analysis of environmental samples.

CrAssphage Concentrations in Wastewater and in the Aquatic Environment

All wastewater influent and effluent samples were positive for crAssphage with no seasonal patterns discernible (Table 1; Table S1). In the wastewater influent samples, crAssphage concentrations varied between 2.2 × 105 gc/L and 1.2 × 109 gc/L, with the lowest concentrations observed in the samples collected at Llanrwst (serving approx. 4000 inhabitants) and the highest concentrations in samples derived from the Ganol WWTP (serving approx. 82,000 inhabitants). These concentrations are lower than previously observed in untreated wastewater in Tampa, Florida (109–1010 gc/L) (Ahmed et al. 2018) and in Catalonia, Spain (108–109 gc/L) (García-Aljaro et al. 2017). The differences in crAssphage concentrations in untreated wastewater may be due to the differences in geographic viral distribution (Stachler and Bibby 2014), WWTP sizes and the level of urbanization, which has been shown to affect crAssphage abundance among inhabitants (Honap et al. 2018).

CrAssphage concentrations in wastewater effluent samples ranged between 1.7 × 105 and 2.0 × 108 gc/L. In general, we observed a very low (less than 1-log10) reduction of crAssphage during wastewater treatment at Llanrwst (using biofilter treatment), whereas the concentrations were reduced by up to 2-log10 at the Betws-y-Coed and Tal-y-Bont WWTPs (using activated sludge treatment). This suggests that there is a difference in the performance of the different types of wastewater treatment methods. To date, no quantitative studies have assessed the prevalence of crAssphage in treated wastewater.

CrAssphage was detected in 92% of all river and estuarine water samples (Table 1; Table S1). The concentration of crAssphage was approx. 1–3 log10 lower in the surface water relative to the treated wastewater and varied between 3.0 × 103 gc/L and 2.5 × 107 gc/L. The highest concentrations were observed in the SW3 samples (3.0 × 106 gc/L on average) derived approximately 200 m downstream of the Llanrwst wastewater discharge point. These concentrations are consistent with the previously reported crAssphage concentrations in wastewater-polluted surface water (104–106 gc/L) (Stachler et al. 2018). No seasonal patterns in the crAssphage concentrations in the surface water samples were observed.

The sediment and mussel samples had more negatives than the wastewater and surface water samples (Table 1 and S1) with concentrations up to 1.9 × 104 gc/g sediment or mussel digestive tissue. No seasonality or significant differences were observed among the sediment and mussel sample crAssphage concentrations. Detection of crAssphage in sediment or in mussel samples has not previously been reported.

Comparison of Enteric Virus and CrAssphage Concentrations in Different Sample Types

CrAssphage concentrations in wastewater, surface water, sediment and mussel samples were compared with the human-specific NoVGI, NoVGII, SaV, AdV and JCV concentrations of the same samples, which were determined in a previous study (Farkas et al. 2018). The enteric viruses were less prevalent and were present at lower concentrations than crAssphage. The virus-positive wastewater samples, both influent and effluent wastewater samples, had approx. 2 log10 higher crAssphage than enteric virus titers (Fig. 1). A weak correlation (r < 0.5) among viral titers was observed in wastewater influent except for the crAssphage/JCV and the AdV/SaV pairs, which showed moderate correlation (r = 0.50 and 0.55; Table 2 and S2). The reduction rate of all tested viruses in the WWTPs were very similar, resulting in approximately 1 log10 reduction during biofilter treatment and 1–2 log10 reduction during activated sludge treatment. Similar reduction trends have been reported in other studies on enteric virus and coliphage titers in activated sludge and biofilter-treated wastewater (Campos et al. 2016; Hewitt et al. 2011; Kitajima et al. 2014; Sidhu et al. 2017). In wastewater effluent, strong correlation was noted for crAssphage/JCV (r = 0.67; Table 2). The crAssphage and JCV showed a linear correlation in the untreated wastewater samples, however, such correlation was only observed in the wastewater effluents that underwent biofilter treatment (R2 = 0.81; Figure S1).

Fig. 1

Box plot showing the median concentration (min–max) of crAssphage (blue), human adenovirus (AdV; green), JC polyomavirus (JCV; gray) norovirus GI (NoVGI; purple) and GII (NoVGII, yellow) and sapovirus (SaV; red) determined in (A) wastewater influent (GI, BI, LI and TI) and effluent (BE, LE, TE) and in (B) surface water (SW1-4), sediment (Sed1, 2, 4) and mussel (SF1, 2). Sample concentration below the limit of quantification was considered 100 gc/l or gc/g. Circles: outliers, starts extreme outliers

Table 2

Spearman correlation between crAssphage and human enteric virus titers

Sample type	Viruses	R	p
Wastewater influent	CrAssphage–AdV	0.31*	0.03
	CrAssphage–JCV	0.50***	< 0.01
	CrAssphage–NoVGI	0.26	0.07
	CrAssphage–NoVGII	0.34*	0.02
	CrAssphage–SaV	0.265	0.07
Wastewater effluent	CrAssphage–AdV	0.13	0.45
	CrAssphage–JCV	0.67***	< 0.01
	CrAssphage–NoVGI	0.20	0.23
	CrAssphage–NoVGII	0.43**	< 0.01
	CrAssphage–SaV	0.21	0.21
Surface water	CrAssphage–AdV	− 0.02	0.88
	CrAssphage–JCV	0.49***	< 0.01
	CrAssphage–NoVGI	0.08	0.55
	CrAssphage–NoVGII	0.07	0.64
	CrAssphage–SaV	− 0.07	0.63
Sediment	CrAssphage–AdV	0.37*	0.02
	CrAssphage–JCV	0.09	0.61
	CrAssphage–NoVGI	− 0.08	0.63
	CrAssphage–NoVGII	0.24	0.14

R: correlation coefficient. Further data are shown in Supplementary Table 2. *p < 0.05; **p < 0.01; ***p < 0.001. No SaV was detected in the sediment

In the surface water samples, crAssphage was more prevalent than the enteric viruses. Only AdV was detected in a similar number of samples to crAssphage (85–92%), whereas the detection rates of the JCV (54–92%) and the NoVs and SaV (0–23%) were lower (Table 1). The concentration of crAssphage was 1–2 log10 higher than the concentration of the AdV and JCV, and 2–5 log10 higher than that observed for NoVs and SaV (Fig. 1). In the sediment and mussel samples, the detection rate of crAssphage (62–75%) was slightly lower than the detection rate of AdV (77–85%) but higher than observed for the JCV (23–46%) and NoVs (0–31%; Table 1). The concentration of crAssphage in sediment was similar to the concentration of AdV and JCV in sediment and slightly higher in mussel (Fig. 1). Weak or no correlation was observed among viruses in surface water, sediment and shellfish samples (Table S2). A linear correlation between crAssphage and JCV was only observed in the highly polluted surface water sample (SW3), where the highest viral titers were observed (Figure S1). The lack of correlation between viral titers is probably due to the low copy numbers and the low number of positive samples.

Environmental Implications and Future Research

In this study, the concentration of crAssphage in wastewater and in wastewater-contaminated environments was tracked over a period of 1 year for the first time. Similar to AdV and JCV, which are considered markers for human wastewater pollution (Rachmadi et al. 2016; Rames et al. 2016), crAssphage incidence displayed no seasonal patterns and was prevalent in all sample types. These findings suggest that the use of crAssphage as a fecal contamination indicator enables the assessment of pollutant transport. This was the first time crAssphage detection was addressed in sediment and in bivalve shellfish. CrAssphage was frequently detected in mussel and sediment, collected at the wastewater-contaminated areas, demonstrating that crAssphage may be a suitable fecal pollution indicator in those matrices.

Future research should explore the usefulness of qPCR assays targeting different regions of the crAssphage genome to enhance specificity, explore human specificity and exploit the recently discovered new crAss-like genomes to identify geographic differences in diversity and assess correlations between the concentrations of crAss-like phages and a wide range of enteric viruses and other viral and bacterial indicators. The prevalence of crAssphage in waters not affected by wastewater contamination also needs to be established. Furthermore, with the very recent discovery of a crAss-like phage isolate (Guerin et al. 2018) and the successful in vitro culturing of that strain (Shkoporov et al. 2018) future approaches should be focused on comparative description of crAssphage, F-RNA phage and AdV/JCV infectivity, survival and persistence in wastewater and in the aquatic environment, including surface and groundwater, and in bivalve shellfish and fresh produce. Once established, crAssphage may be a suitable marker for routine bivalve shellfish testing in commercial fisheries.

Below is the link to the electronic supplementary material.

Supplementary material 1 (DOCX 61 KB)