Applicability of crAssphage as a performance indicator for viral reduction during activated sludge wastewater treatment
Waterborne infections continue to have far-reaching public health and socioeconomic consequences in both the developed and developing worlds. WHO estimates that unsafe water, sanitation, and hygiene cause ~2 million deaths annually, mainly related to infectious diarrhea (WHO ). While viral pathogens are found in water, most countries still use classic fecal indicator bacteria (FIB) which have well-known shortcomings, including insufficiently reflecting viral risk to human health. Many reasons contribute to this, including their increased susceptibility to wastewater/water treatment, sensitivity to disinfectants, low tolerance to environmental conditions, and co-occurrence in animal species (Boehm et al. ; Harwood et al. ; Payment and Locas ). Enteric viruses are the most prevalent causative agents of gastroenteritis worldwide. Over 150 human pathogenic viruses have been detected in watercourses (Fong and Lipp ; Rodríguez-Lázaro et al. ). Thus, it is not practical to test the water samples for all of the enteric viruses; thus, surrogate indicators are still needed. Viral fecal pollution indicators have previously been suggested, not yet been extensively utilized for regulatory uses. These previously discovered markers are divided into two categories: human pathogens and bacteriophages. Human pathogens formerly considered as viral water quality indicators include human adenovirus (HAdV), human polyomavirus (HPyV), and Aichi virus 1 (AiV-1) (Albinana-Gimenez et al. ; Hamza et al. ; Kitajima et al. ). These viral indicators have the benefit of being very specific to humans, but they are limited by low and unpredictable quantities in wastewater. Bacteriophages have also been proposed as indicators of water quality. These phage-based approaches fulfill the criteria for an ideal viral water quality indicator, such as higher concentrations in wastewater than many human pathogenic viruses and rapid and easy culturability than human viral pathogens (Grabow ). Limited specificity to human fecal waste and lower concentrations than other recently found viral targets are potential obstacles to the use of coliphage as an indicator (Grabow ; Jofre et al. ). CrAssphage was identified by metagenomic analysis and was claimed to be the most prevalent virus in the human gut (Dutilh et al. ) before being proven to be globally dispersed (Edwards et al. ). It was highly abundant in the USA and Europe compared to Africa and Asia (Stachler and Bibby ). Further metagenomic analysis revealed that crAssphage is highly specific to human fecal material and was proposed for human fecal source identification (Stachler and Bibby ). However, previous research has identified crAssphage in seagull, dog, chicken, cat, and cow feces at lower quantities than in human sewage (Ahmed et al. , ; Stachler and Bibby ). Recent studies have also effectively identified crAssphage in various water matrices impacted with human fecal pollution including river water (Ballesté et al. ; Farkas et al. ), lake (Ahmed et al. ), stormwater (Ahmed et al. ), and seawater (Sala-Comorera et al. ; Sangkaew et al. ), showing that crAssphage may be used to identify viral contamination by municipal wastewater. The presence of crAssphage in sewage-impacted waters has also been linked to a higher risk for human health (Crank et al. ). Despite the fact that crAssphage has been studied extensively as a human fecal marker, few studies have yet been performed to assess crAssphage as a process indicator in conventional activated sludge wastewater treatment facilities (Tandukar et al. ; Wu et al. ). Also, to our knowledge, no data are available on crAssphage in the Egyptian environment. Thus, the primary objectives of the present study were to assess crAssphage removal during activated sludge wastewater treatment and the suitability of crAssphage as a viral process indicator. Over 1-year study, the occurrence and abundance of crAssphage in influent and effluent samples of two WWTPs in Greater Cairo were determined. Moreover, its association with human enteric viruses including HAdV, HPyV, and bocaviruses that showed high dissemination in the Egyptian environment before was demonstrated. HAdV can cause a variety of diseases including gastrointestinal, respiratory, and urinary infections. HAdV is frequently identified in a variety of water matrices (Bofill-Mas et al. ; Hamza et al. , ; Hewitt et al. ; Pina et al. ). Thus, it has been considered as an indicator of human fecal contamination in water. HPyV usually does not produce symptoms in healthy people, but it may cause severe infections in immunocompromised people. HPyV is found in wastewater across the world, and several studies have proposed HPyV as a viral fecal contamination indicator (Albinana-Gimenez et al. ; Bofill-Mas et al. ). HBoV has been isolated from stool samples collected from patients with gastroenteritis and respiratory tract samples (Allander ; Rizk et al. ; Weissbrich et al. ). Also, different studies showed that HBoV was highly abundant in environmental water samples (Blinkova et al. ; Hamza et al. ). 

Study sites and sampling A total of 46 sewage samples were collected anonymously from two wastewater treatment facilities: WWTP-A and WWTP-B, located in Greater Cairo. Samples were taken monthly as grab samples over a one-year study course between 08/2018 and 07/2019. The designed capacities of these WWTPs are 330,000 m 3 /day for WWTP-A and 600,000 m 3 /day for WWTP-B. The populations served by the WWTPs are approximately 1,320,000 for WWTP-A and 2,200,000 for WWTP-B. Activated sludge is implemented in all WWTPs as a secondary treatment process. Five-liter samples were collected from both the influent and effluent. Samples were collected in sterile bottles and transported within 1 h to the laboratory for analysis. Virus concentration Virus concentration was performed employing the virus adsorption elution method reported earlier by USEPA ( ). In brief, samples were processed by adding a final concentration of 0.05 M MgCl 2 . The pH was then adjusted to 3.5 with 1 N HCl . . Then samples were filtrated by a negatively charged HA nitrocellulose membrane with 0.45 m pore size and 142 mm diameter. Prior to the viral recovery using 70 ml of organic elution buffer (3% beef extract, 0.05 M glycine, pH 9.4), the membrane was washed with 0.5 mM H 2 SO 4 , pH 4. The eluates were subjected to an organic flocculation technique for viral re-concentration. DNA extraction Viral DNA was extracted from 200 µl of concentrated suspension using QIAamp DNA Blood Mini Kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions. Sterile nuclease-free water was included in each set of extractions as a negative control to monitor cross-contamination. Since environmental samples may have PCR inhibitors which can lead to underestimation of viral concentrations, the frequency of positive samples, murine norovirus (MNV), was added to the samples during extraction as an exogenous control to identify the occurrence of PCR inhibition. A comparison of the Ct value of the MNV to that of the negative control showed no inhibitory effect (data not shown). Quantification of viral genome by qPCR In this current study, HBoVs, HAdV, and HPyV were included as human viruses and crAssphage was tested as an indicator virus. Table contains a list of all the primers utilized in the current study. The quantification methodology for HBoV-1 targets the NP1 gene, according to Hamza et al. ( ). The quantification of HBoV-2, -3, and -4 employed a single-sense primer, whereas qPCR for HBoV-2 and -4 used the same antisense primer, according to Kantola et al. ( ). DNA standards of HBoVs were prepared according to Hamza et al. ( ). HAdV qPCR assay was used by Heim et al. ( ), and HPyV qPCR was used according to Biel et al. ( ). The DNA standards of HAdV and HPyV were prepared according to Hamza et al. ( ). CrAssphage concentrations were determined using the CPQ_56 assay developed by Stachler et al. ( ). TaqMan probe assay was used for the quantification of all viruses except HBoV-2/4, and 3 SYBR green qPCR assay was conducted. TaqMan real-time qPCR reactions were performed in a total volume of 20 µl containing 1 × (10 µl) Quantitect probe PCR kit (Qiagen, Hilden, Germany), 0.5 µM for both forward and reverse primers, 0.2 µM Taqman probe, and 2 µl DNA template. The qPCR program was 95 °C for 15 min as the initial activation step for HotStart Taq DNA Polymerase and 45 cycles of 2-step cycling for 15 s at 94 °C and 1 min at 60 °C. HBoV-2/4 and 3 SYBR green assays were conducted using Maxima SYBR Green qPCR Master Mix Kit (Thermo Scientifc). The PCR conditions were 10 min initial denaturation step at 95 °C, 45 cycles of denaturation at 95 °C for 15 s, and annealing extension at 60 °C for 1 min. Amplification was followed by one cycle of melting curve analysis. Dissociation was carried out from 60 to 95 °C with a temperature ramp of 0.05 °C/s. Analysis indicated a melting peak 81.5 °C ± 0.3 °C for HBoV 2/4 and 80 °C ± 0.2 °C for HBoV-3. In order to exclude data of cross-contamination, negative controls (NTC) were included in each run as nuclease-free water. All NTCs were negative throughout the qPCRs. The amplification and data analysis were performed using Rotorgene 6000. Statistical analysis The viral concentrations were expressed as gc/l of wastewater. Kruskal–Wallis test was used for multiple comparison procedures to determine possible significant variations in the concentrations of crAssphage and human enteric viruses. Human viruses concentrations were normalized as the ratios over crAssphage concentrations to evaluate differential fate. Wilcox test was used to compare the ratio of enteric viruses over crAssphage from influents and effluents. Spearman’s rank correlation coefficients ( r ) were calculated between viral concentrations using two-tailed 95% confidence intervals. 

A total of 46 sewage samples were collected anonymously from two wastewater treatment facilities: WWTP-A and WWTP-B, located in Greater Cairo. Samples were taken monthly as grab samples over a one-year study course between 08/2018 and 07/2019. The designed capacities of these WWTPs are 330,000 m 3 /day for WWTP-A and 600,000 m 3 /day for WWTP-B. The populations served by the WWTPs are approximately 1,320,000 for WWTP-A and 2,200,000 for WWTP-B. Activated sludge is implemented in all WWTPs as a secondary treatment process. Five-liter samples were collected from both the influent and effluent. Samples were collected in sterile bottles and transported within 1 h to the laboratory for analysis. Virus concentration Virus concentration was performed employing the virus adsorption elution method reported earlier by USEPA ( ). In brief, samples were processed by adding a final concentration of 0.05 M MgCl 2 . The pH was then adjusted to 3.5 with 1 N HCl . . Then samples were filtrated by a negatively charged HA nitrocellulose membrane with 0.45 m pore size and 142 mm diameter. Prior to the viral recovery using 70 ml of organic elution buffer (3% beef extract, 0.05 M glycine, pH 9.4), the membrane was washed with 0.5 mM H 2 SO 4 , pH 4. The eluates were subjected to an organic flocculation technique for viral re-concentration. DNA extraction Viral DNA was extracted from 200 µl of concentrated suspension using QIAamp DNA Blood Mini Kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions. Sterile nuclease-free water was included in each set of extractions as a negative control to monitor cross-contamination. Since environmental samples may have PCR inhibitors which can lead to underestimation of viral concentrations, the frequency of positive samples, murine norovirus (MNV), was added to the samples during extraction as an exogenous control to identify the occurrence of PCR inhibition. A comparison of the Ct value of the MNV to that of the negative control showed no inhibitory effect (data not shown). Quantification of viral genome by qPCR In this current study, HBoVs, HAdV, and HPyV were included as human viruses and crAssphage was tested as an indicator virus. Table contains a list of all the primers utilized in the current study. The quantification methodology for HBoV-1 targets the NP1 gene, according to Hamza et al. ( ). The quantification of HBoV-2, -3, and -4 employed a single-sense primer, whereas qPCR for HBoV-2 and -4 used the same antisense primer, according to Kantola et al. ( ). DNA standards of HBoVs were prepared according to Hamza et al. ( ). HAdV qPCR assay was used by Heim et al. ( ), and HPyV qPCR was used according to Biel et al. ( ). The DNA standards of HAdV and HPyV were prepared according to Hamza et al. ( ). CrAssphage concentrations were determined using the CPQ_56 assay developed by Stachler et al. ( ). TaqMan probe assay was used for the quantification of all viruses except HBoV-2/4, and 3 SYBR green qPCR assay was conducted. TaqMan real-time qPCR reactions were performed in a total volume of 20 µl containing 1 × (10 µl) Quantitect probe PCR kit (Qiagen, Hilden, Germany), 0.5 µM for both forward and reverse primers, 0.2 µM Taqman probe, and 2 µl DNA template. The qPCR program was 95 °C for 15 min as the initial activation step for HotStart Taq DNA Polymerase and 45 cycles of 2-step cycling for 15 s at 94 °C and 1 min at 60 °C. HBoV-2/4 and 3 SYBR green assays were conducted using Maxima SYBR Green qPCR Master Mix Kit (Thermo Scientifc). The PCR conditions were 10 min initial denaturation step at 95 °C, 45 cycles of denaturation at 95 °C for 15 s, and annealing extension at 60 °C for 1 min. Amplification was followed by one cycle of melting curve analysis. Dissociation was carried out from 60 to 95 °C with a temperature ramp of 0.05 °C/s. Analysis indicated a melting peak 81.5 °C ± 0.3 °C for HBoV 2/4 and 80 °C ± 0.2 °C for HBoV-3. In order to exclude data of cross-contamination, negative controls (NTC) were included in each run as nuclease-free water. All NTCs were negative throughout the qPCRs. The amplification and data analysis were performed using Rotorgene 6000. 

Virus concentration was performed employing the virus adsorption elution method reported earlier by USEPA ( ). In brief, samples were processed by adding a final concentration of 0.05 M MgCl 2 . The pH was then adjusted to 3.5 with 1 N HCl . . Then samples were filtrated by a negatively charged HA nitrocellulose membrane with 0.45 m pore size and 142 mm diameter. Prior to the viral recovery using 70 ml of organic elution buffer (3% beef extract, 0.05 M glycine, pH 9.4), the membrane was washed with 0.5 mM H 2 SO 4 , pH 4. The eluates were subjected to an organic flocculation technique for viral re-concentration. 

Viral DNA was extracted from 200 µl of concentrated suspension using QIAamp DNA Blood Mini Kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions. Sterile nuclease-free water was included in each set of extractions as a negative control to monitor cross-contamination. Since environmental samples may have PCR inhibitors which can lead to underestimation of viral concentrations, the frequency of positive samples, murine norovirus (MNV), was added to the samples during extraction as an exogenous control to identify the occurrence of PCR inhibition. A comparison of the Ct value of the MNV to that of the negative control showed no inhibitory effect (data not shown). 

In this current study, HBoVs, HAdV, and HPyV were included as human viruses and crAssphage was tested as an indicator virus. Table contains a list of all the primers utilized in the current study. The quantification methodology for HBoV-1 targets the NP1 gene, according to Hamza et al. ( ). The quantification of HBoV-2, -3, and -4 employed a single-sense primer, whereas qPCR for HBoV-2 and -4 used the same antisense primer, according to Kantola et al. ( ). DNA standards of HBoVs were prepared according to Hamza et al. ( ). HAdV qPCR assay was used by Heim et al. ( ), and HPyV qPCR was used according to Biel et al. ( ). The DNA standards of HAdV and HPyV were prepared according to Hamza et al. ( ). CrAssphage concentrations were determined using the CPQ_56 assay developed by Stachler et al. ( ). TaqMan probe assay was used for the quantification of all viruses except HBoV-2/4, and 3 SYBR green qPCR assay was conducted. TaqMan real-time qPCR reactions were performed in a total volume of 20 µl containing 1 × (10 µl) Quantitect probe PCR kit (Qiagen, Hilden, Germany), 0.5 µM for both forward and reverse primers, 0.2 µM Taqman probe, and 2 µl DNA template. The qPCR program was 95 °C for 15 min as the initial activation step for HotStart Taq DNA Polymerase and 45 cycles of 2-step cycling for 15 s at 94 °C and 1 min at 60 °C. HBoV-2/4 and 3 SYBR green assays were conducted using Maxima SYBR Green qPCR Master Mix Kit (Thermo Scientifc). The PCR conditions were 10 min initial denaturation step at 95 °C, 45 cycles of denaturation at 95 °C for 15 s, and annealing extension at 60 °C for 1 min. Amplification was followed by one cycle of melting curve analysis. Dissociation was carried out from 60 to 95 °C with a temperature ramp of 0.05 °C/s. Analysis indicated a melting peak 81.5 °C ± 0.3 °C for HBoV 2/4 and 80 °C ± 0.2 °C for HBoV-3. In order to exclude data of cross-contamination, negative controls (NTC) were included in each run as nuclease-free water. All NTCs were negative throughout the qPCRs. The amplification and data analysis were performed using Rotorgene 6000. 

The viral concentrations were expressed as gc/l of wastewater. Kruskal–Wallis test was used for multiple comparison procedures to determine possible significant variations in the concentrations of crAssphage and human enteric viruses. Human viruses concentrations were normalized as the ratios over crAssphage concentrations to evaluate differential fate. Wilcox test was used to compare the ratio of enteric viruses over crAssphage from influents and effluents. Spearman’s rank correlation coefficients ( r ) were calculated between viral concentrations using two-tailed 95% confidence intervals. 

Detection rates of crAssphage and human viruses Over a one-year study, all viruses could be detected in the tested wastewater samples at different frequencies (Table ). HboV-2/4, HBoV-3, and crAssphage were the most frequently detected in influent samples of WWTPs. Influent samples were positive for at least 5 out of six viruses. No clear seasonal pattern was observed for neither the human pathogenic viruses nor the indicators. In effluent samples, there was a slight difference between the detection rates of human enteric viruses, except for HBoV-1; it was only detected in seven and four samples of WWTP-A and B, respectively. CrAssphage was identified in 100% ( n = 23) of effluent samples (Table ). Concentrations of crAssphage and human viruses The concentration of crAssphage in influent samples was significantly higher than those of HAdV, HPyV, and HBoVs (ANOVA), p < 0.0001. In wastewater influent samples, the concentration of crAssphage ranged from 1.45E + 04 to 1.02E + 08 gc/l in WWTP-A and from 3.51E + 05 to 2.39E + 08 gc/l in WWTP-B (Fig. ). Regarding human viruses, HBoVs were detected in effluents wastewater samples at concentration orders of magnitude lower than HAdV and HPyV (Fig. ). Similarly, in effluent samples of WWTPs, the concentration of crAssphage was significantly higher than HBoVs and HPyV, ranging from 1.25E + 04 to 6.29E + 06 gc/l in WWTP-A and 7.49E + 04 to 7.88E + 06 gc/l for WWTP-B (Kruskal–Wallis test, P = 0.0001) (Fig. ). However, no significant difference between crAssphage concentration and HAdV in effluent wastewater samples was identified. Additionally, in Fig. , the overall annual viral concentration is compared between influent and effluent samples of WWTPs. CrAssphage and human enteric viruses concentrations were relatively stable during the study course. Viral reduction during the treatment process The annual mean reduction of all tested viruses was relatively similar, varying between ~ 1 ± 0.64 log 10 for HBoVs, 0.84 ± 0.5 log 10 for HAdV, 1.1 ± 0.8log 10 for HPyV, and 1.32 ± 0.7log 10 for crAssphage. No significant difference between viral reduction was observed. Figure shows ratios of human pathogenic viruses concentrations at influent and effluent samples, normalized over crAssphage concentrations. The ratio showed a slight increase from influents to effluents. These ratios were used to assess the differences in the fate of crAssphage and other human viruses during the wastewater treatment process. Only samples with both targets within the quantifiable range were considered in pair comparison. CrAssphage correlation with human viruses Spearman’s rank correlation coefficients ( r ) were determined between human viruses and crAssphage concentrations in influent and effluent wastewater samples. As noted in Table , a strong positive correlation ( P = 0.001) was found between crAssphage and HPyV in influent samples. Also, a significant correlation ( P < 0.05) was detected between crAssphage and both HAdV and HPyV in the treated samples. 

Over a one-year study, all viruses could be detected in the tested wastewater samples at different frequencies (Table ). HboV-2/4, HBoV-3, and crAssphage were the most frequently detected in influent samples of WWTPs. Influent samples were positive for at least 5 out of six viruses. No clear seasonal pattern was observed for neither the human pathogenic viruses nor the indicators. In effluent samples, there was a slight difference between the detection rates of human enteric viruses, except for HBoV-1; it was only detected in seven and four samples of WWTP-A and B, respectively. CrAssphage was identified in 100% ( n = 23) of effluent samples (Table ). Concentrations of crAssphage and human viruses The concentration of crAssphage in influent samples was significantly higher than those of HAdV, HPyV, and HBoVs (ANOVA), p < 0.0001. In wastewater influent samples, the concentration of crAssphage ranged from 1.45E + 04 to 1.02E + 08 gc/l in WWTP-A and from 3.51E + 05 to 2.39E + 08 gc/l in WWTP-B (Fig. ). Regarding human viruses, HBoVs were detected in effluents wastewater samples at concentration orders of magnitude lower than HAdV and HPyV (Fig. ). Similarly, in effluent samples of WWTPs, the concentration of crAssphage was significantly higher than HBoVs and HPyV, ranging from 1.25E + 04 to 6.29E + 06 gc/l in WWTP-A and 7.49E + 04 to 7.88E + 06 gc/l for WWTP-B (Kruskal–Wallis test, P = 0.0001) (Fig. ). However, no significant difference between crAssphage concentration and HAdV in effluent wastewater samples was identified. Additionally, in Fig. , the overall annual viral concentration is compared between influent and effluent samples of WWTPs. CrAssphage and human enteric viruses concentrations were relatively stable during the study course. Viral reduction during the treatment process The annual mean reduction of all tested viruses was relatively similar, varying between ~ 1 ± 0.64 log 10 for HBoVs, 0.84 ± 0.5 log 10 for HAdV, 1.1 ± 0.8log 10 for HPyV, and 1.32 ± 0.7log 10 for crAssphage. No significant difference between viral reduction was observed. Figure shows ratios of human pathogenic viruses concentrations at influent and effluent samples, normalized over crAssphage concentrations. The ratio showed a slight increase from influents to effluents. These ratios were used to assess the differences in the fate of crAssphage and other human viruses during the wastewater treatment process. Only samples with both targets within the quantifiable range were considered in pair comparison. CrAssphage correlation with human viruses Spearman’s rank correlation coefficients ( r ) were determined between human viruses and crAssphage concentrations in influent and effluent wastewater samples. As noted in Table , a strong positive correlation ( P = 0.001) was found between crAssphage and HPyV in influent samples. Also, a significant correlation ( P < 0.05) was detected between crAssphage and both HAdV and HPyV in the treated samples. 

The concentration of crAssphage in influent samples was significantly higher than those of HAdV, HPyV, and HBoVs (ANOVA), p < 0.0001. In wastewater influent samples, the concentration of crAssphage ranged from 1.45E + 04 to 1.02E + 08 gc/l in WWTP-A and from 3.51E + 05 to 2.39E + 08 gc/l in WWTP-B (Fig. ). Regarding human viruses, HBoVs were detected in effluents wastewater samples at concentration orders of magnitude lower than HAdV and HPyV (Fig. ). Similarly, in effluent samples of WWTPs, the concentration of crAssphage was significantly higher than HBoVs and HPyV, ranging from 1.25E + 04 to 6.29E + 06 gc/l in WWTP-A and 7.49E + 04 to 7.88E + 06 gc/l for WWTP-B (Kruskal–Wallis test, P = 0.0001) (Fig. ). However, no significant difference between crAssphage concentration and HAdV in effluent wastewater samples was identified. Additionally, in Fig. , the overall annual viral concentration is compared between influent and effluent samples of WWTPs. CrAssphage and human enteric viruses concentrations were relatively stable during the study course. 

The annual mean reduction of all tested viruses was relatively similar, varying between ~ 1 ± 0.64 log 10 for HBoVs, 0.84 ± 0.5 log 10 for HAdV, 1.1 ± 0.8log 10 for HPyV, and 1.32 ± 0.7log 10 for crAssphage. No significant difference between viral reduction was observed. Figure shows ratios of human pathogenic viruses concentrations at influent and effluent samples, normalized over crAssphage concentrations. The ratio showed a slight increase from influents to effluents. These ratios were used to assess the differences in the fate of crAssphage and other human viruses during the wastewater treatment process. Only samples with both targets within the quantifiable range were considered in pair comparison. 

Spearman’s rank correlation coefficients ( r ) were determined between human viruses and crAssphage concentrations in influent and effluent wastewater samples. As noted in Table , a strong positive correlation ( P = 0.001) was found between crAssphage and HPyV in influent samples. Also, a significant correlation ( P < 0.05) was detected between crAssphage and both HAdV and HPyV in the treated samples. 

No data are available on crAssphage in the Egyptian environment. The primary objectives of the current work were to assess crAssphage reduction during wastewater treatment and its usefulness as a viral process indicator of the treatment process. Thus, targeted research typically selects pathogens that are more relevant to humans or that are more abundant in wastewater. Some of these viruses have been involved in the present study. Samples were taken from influents and effluents of two WWTPs, and the results of crAssphage genome levels were compared with that of different human enteric viruses. All wastewater samples tested positive for crAssphage (Table ), with no identifiable seasonal variations. Also, both WWTPs showed relatively the same range of viral concentrations (Fig. ) due to using the same treatment technology regardless of their treatment capacity. In raw sewage, the annual crAssphage concentrations varied between 1.45E + 04 and 2.39E + 08 gc/l (Fig. ). The log 10 concentrations of crAssphage in our study are lower than the previously detected values in Florida, USA (9–10 log 10 gc/l) (Ahmed et al. ), Spain (8.4–9.9log 10 gc/l) (García‐Aljaro et al. ), Japan (10.98–12.03log 10 gc/l) (Malla et al. ), Indiana, USA (8.23 ± 0.36 log 10 gc/l) (Wu et al. ), and UK (5.3–9 log 10 gc/l) (Farkas et al. ). The detected crAssphage concentration in this study is relatively the same as the previously reported in Thailand (5.23–7.19 log 10 gc/l) (Kongprajug et al. ). On the other hand, crAssphage concentrations in effluent samples ranged from 1.25E + 04 to 7.88E + 06 gc/l which is relatively the same range as determined in effluent samples examined by Kongprajug et al. ( ). Whereas others from different geographical areas have reported higher concentrations of crAssphage in effluent samples (Ballesté et al. ; Malla et al. ; Tandukar et al. ). The difference in crAssphage between different studies could be attributed to the different geographic distribution of viruses, the capacity of WWTPs, and the difference in industrialized lifestyle (Honap et al. ; Stachler and Bibby ). Moreover, using different concentration techniques, processed water samples, and the quantification method can contribute to the discrepancies in the viral concentrations from different investigations. Additionally, the diversity of crAssphage in the human gut has been recently described (Edwards et al. ). It is likely that such natural diversity in crAssphage was not detected by the CPQ56 assay which was designed based on the prototype crAssphage. A comparison between the level of crAssphage and human viruses showed that in influent and effluent samples, the mean concentration of crAssphage has one order of magnitude higher than HAdV and HPyV and three orders of magnitude higher than HBoVs (Fig. ). Similar trends have been observed in recent reports. Farkas et al. ( ) estimated that all viruses-positive wastewater samples contained approximately 2 log 10 higher crAssphage than other enteric viruses such as NoV, AdV, and HPyV. Also, crAssphage was up to 5 orders of magnitude higher than HPyV in wastewater (Stachler et al. ). The present data showed no seasonal pattern for human viruses and crAssphage. This finding is consistent with other year-long monitoring investigations that have also revealed the constant existence of crAssphage in treated wastewater without seasonal variations (Crank et al. ; Farkas et al. ; Wu et al. ). Meanwhile, the levels of human enteric viruses may have more variations according to the clinical situation of the population. In general, the annual mean reduction of all tested viruses between ~ 1 ± 0.64 log 10 for HBoVs, 0.84 ± 0.5 log 10 for HAdV, 1.1 ± 0.8 log 10 for HPyV, and 1.32 ± 0.7 log 10 for crAssphage. Our results agree with Farkas et al. who found up to 2 log 10 reduction in crAssphage using activated sludge treatment and lower reduction levels (1 log 10 ) by biofilter treatment (Farkas et al. ). Tandukar et al. ( ) observed that crAssphage had the greatest removal ratio (3.3 ± 1.0 log 10 ) among studied enteric viruses such as HPyV, NoVGII, EV, and AiV. Accordingly, Tandukar et al. ( ) argued that crAssphage cannot be used as an indication of viral reduction throughout wastewater treatment. Another study by Wu et al. ( ) reported that the log 10 reduction of crAssphage (2.88 ± 0.68) during wastewater treatment was relatively higher than HAdV (2.24 ± 0.53) or HPyV (1.51 ± 0.37). Although crAssphage had a greater initial concentration in the main influent, the variation in removal is likely limited to crAssphage since it was eliminated in a higher fraction than HAdV or HPyV after secondary treatment (Wu et al. ). Ultimately, the log 10 removal rate of HBoV, HAdV, and HPyV during activated sludge treatment was reported as 0.35–1 log 10 , 0.8–3.7 log 10 , and 1.0–3.7 log 10 , respectively (Hamza et al. , ; Kitajima et al. ; Sangkaew et al. ; Schmitz et al. ). While crAssphage log reductions are less variable than that of other viruses, the results suggest that crAssphage has a high potential as a process indicator for pathogenic viral reduction during wastewater treatment. The ratio of viruses over crAssphage (Fig. ) has been slightly increased from inlet to outlet samples indicating slightly lower removal of human viruses than crAssphage (Wilcox test, P > 0.05). Notably, crAssphage was detected in all samples, and lower detected rates have been identified for other human viruses. Data normalization over crAssphage has been proposed before to assess the performance of the wastewater treatment process. For instance, Wu et al. ( ) reported that ratios of HAdV/CPQ56 and HPyV/CPQ56 increased during secondary treatment, indicating that both viruses were removed relatively smaller than crAssphage. However, both viruses had the same removal mechanism owing to the correlation between crAssphage and HAdV and HPyV. A correlation between viral human fecal indicators and viral pathogens in wastewater is required to obtain an accurate picture of the viral risk posed by human feces. The present study compared the concentration of crAssphage with HBoVs, HPyV, and HAdV. In influent samples, the co-occurrence analysis between crAssphage and human viruses revealed a strong positive correlation between crAssphage and HPyV. However, HAdV and HPyV correlated with crAssphage in effluent samples (Table ). This finding is consistent with a report of crAssphage concentration correlating with HPyV and HAdV through a WWTP (Wu et al. ). Similarly, Crank et al. ( ) observed a positive correlation between crAssphage and DNA viruses (HPyV, HBoV) in raw sewage samples and no correlation was found between crAssphage and HEV. Although the virus enrichment approach may affect this association, the correlation between crAssphage and HPyV was stable regardless of the concentration method (Crank et al. ). Additionally, concentrations of crAssphage in raw wastewater correlated positively with the concentrations of HAdV, HPyV, and NoVGII ( p < 0.05), suggesting the applicability of crAssphage as a suitable indicator to estimate human enteric virus concentrations in raw wastewater. Likewise, Farkas et al. ( ) found a positive correlation between HPyV and crAssphage in both influent and effluent samples. It should be noted that locality and crAssphage marker selection in qPCR assays were likely to contribute to the observed correlations (Sabar et al. ). Future studies should investigate which crAssphage markers correlate well with each water-related pathogen in different locations. The ideal viral indicator to assess the performance in wastewater should be prevalent at a high concentration in raw sewage and has similar or more persistence in wastewater treatment than the pathogenic viruses of the reduction target. CrAssphage possesses several properties that would make it a potential viral process indicator during wastewater treatment. In raw sewage, it was the most abundant of the fecal markers utilized in the current investigation, making it easier to determine. The virus was more persistent during the treatment process than human viruses enabling the performance assessment. Also, high values of crAssphage could be found in the effluent samples, promoting the evaluation of the treatment process in terms of log reduction. Alternatively, crAssphage meets Bonde’s criteria for an ideal indicator of waterborne pathogens, which include (i) being present when the pathogens are present, (ii) occurring in greater numbers than pathogens, and (iii) being more resistant to disinfectants and to aqueous environments than the pathogens (NASEM ). 

The current study aimed to assess crAssphage reduction in WWTPs and to evaluate its usefulness as a viral process indicator during the treatment process. When crAssphage was compared to human viruses, crAssphage was highly abundant in both raw and treated wastewater samples without a significant difference in the removal rate. Importantly, crAssphage is associated with different human viruses in raw and treated wastewater samples. Also, the high co-occurrence and comparable destiny of crAssphage to human viruses such as HAdV and HPyV during the treatment process shows that crAssphage and human viral pathogens have similar removal mechanisms. These findings provide additional evidence of the usefulness crAssphage as a process indicator for wastewater treatment. Additionally, the constant high prevalence, abundance, and association with human pathogenic viruses including HAdV and HPyV in wastewater support its use as a conservative viral indicator of human fecal pollution. Since this study compared the fate of crAssphage and human DNA viruses in WWTPs, further evaluation including RNA viruses should be performed. 

