Integrating Metagenomic and Bayesian Analyses to Evaluate the Performance and Confidence of CrAssphage as an Indicator for Tracking Human Sewage Contamination in China

Sources and Prevalence of Cyclospora cayetanensis in Southeastern U.S. Growing Environments

Recent cyclosporiasis outbreaks associated with fresh produce grown in the United States highlight the need to better understand Cyclospora cayetanensis prevalence in U.S. agricultural environments. In this study, C. cayetanensis occurrence was assessed in municipal wastewater sludge, on-farm portable toilets, irrigation pond water, and spent packing house dump tank water in a Southeastern Georgia growing region over two years. Detection of the C. cayetanensis 18S rRNA qPCR gene target in pond samples was 0%, 28%, and 42% (N = 217) depending on the detection definition used, and ≤1% in dump tank samples (N = 46). However, no qPCR detections were confirmed by sequencing, suggesting false detection occurred due to cross-reactions. C. cayetanensis qPCR detections were confirmed in 9% of wastewater sludge samples (N = 76). The human-specific fecal markers HF183 and crAssphage were detected in 33% and 6% of pond samples, respectively, and 4% and 0% of dump tank samples, respectively. Despite community Cyclospora shedding and evidence of human fecal contamination in irrigation water, there was no correlation between C. cayetanensis and HF183 qPCR detections, further supporting that 18S gene target qPCR amplifications were due to cross-reactions. When evaluating C. cayetanensis qPCR environmental detection data, the impact of assay specificity and detection criteria should be considered. Moreover, additional sequence-based testing may be needed to appropriately interpret Cyclospora qPCR environmental data.

Assessment of crAssphage as a human fecal source tracking marker in the lower Great Lakes

CrAssphage or crAss-like phage ranks as the most abundant phage in the human gut and is present in human feces-contaminated environments. Due to its high human specificity and sensitivity, crAssphage is a potentially robust source tracking indicator that can distinguish human fecal contamination from agricultural or wildlife sources. Its suitability in the Great Lakes area, one of the world's most important water systems, has not been well tested. In this study, we tested a qPCR-based quantification method using two crAssphage marker genes (ORF18-mod and CPQ_064) at Toronto recreational beaches along with their adjacent river mouths. Our results showed a 71.4 % (CPQ_064) and 100 % (ORF18-mod) human sensitivity for CPQ_064 and ORF18-mod, and a 100 % human specificity for both marker genes. CrAssphage was present in 57.7 % or 71.2 % of environmental water samples, with concentrations ranging from 1.45 to 5.14 log10 gene copies per 100 mL water. Though concentrations of the two marker genes were strongly correlated, ORF18-mod features a higher human sensitivity and higher positive detection rates in environmental samples. Quantifiable crAssphage was mostly present in samples collected in June and July 2021 associated with higher rainfall. In addition, rivers had more frequent crAssphage presence and higher concentrations than their associated beaches, indicating more frequent and greater human fecal contamination in the rivers. However, crAssphage was more correlated with E. coli and Enterococcus at the beaches than in the rivers, suggesting human fecal sources may be more predominant in driving the increases in E. coli and Enterococcus at the beaches when impacted by river plumes.

Integrating global microbiome data into antibiotic resistance assessment in large rivers

Rivers are important in spreading antimicrobial resistance (AMR). Assessing AMR risk in large rivers is challenged by large spatial scale and numerous contamination sources. Integrating river resistome data into a global framework may help addressing this difficulty. Here, we conducted an omics-based assessment of AMR in a large river (i.e. the Pearl River in China) with global microbiome data. Results showed that antibiotic resistome in river water and sediment was more diversified than that in other rivers, with contamination levels in some river reaches higher than global baselines. Discharge of WWTP effluent and landfill waste drove AMR prevalence in the river, and the resistome level was highly associated with human and animal sources. Detection of 54 risk rank I ARGs and emerging mobilizable mcr and tet(X) highlighted AMR risk in the river reaches with high human population density and livestock pollution. Florfenicol-resistant floR therein deserved priority concerns due to its high detection frequency, dissimilar phylogenetic distance, mobilizable potential, and presence in multiple pathogens. Co-sharing of ARGs across taxonomic ranks implied their transfer potentials in the community. By comparing with global genomic data, we found that Burkholderiaceae, Enterobacteriaceae, Moraxellaceae and Pseudomonadaceae were important potential ARG-carrying bacteria in the river, and WHO priority carbapenem-resistant Enterobacteriaceae, A. baumannii and P. aeruginosa should be included in future surveillance. Collectively, the findings from this study provide an omics-benchmarked assessment strategy for public risk associated with AMR in large rivers.

Metagenomics-assembled analysis revealed the characteristics of antibiotic resistome and community coalescence in the soils irrigated with different irrigation materials

Antibiotic resistance has received widespread attention in recent years. Soil irrigation and fertilization are routine agricultural practices, but also lead to the spread of antibiotic resistance genes (ARGs) in soil-crop system such as via resistome coalescence. Despite community coalescence being ubiquitous and important in natural ecosystems, little research has been done to investigate resistome coalescence during soil irrigation activities. In this study, the characteristics of antibiotic resistome and community coalescence in the soils irrigated with different irrigation materials (wastewater, wastewater-river water, and wastewater-manure) have been revealed by utilizing microcosm experiments and high-throughput sequencing-based metagenomic assembly approaches. Results showed irrigation and coalescence changed soil quality and resistome. Totally, 789 unique ARGs were identified in the irrigation system, including some emerging ARGs. The abundance and diversity of ARGs increased in the coalesced soils, mainly due to the newly imported ARGs from irrigation materials. Relatively, the soils irrigated with wastewater and manure showed higher level of ARGs. Irrigation with the mixtures containing river water caused greater loss of indigenous taxa, while the community structure of mixing treatment with manure changed more dramatically. Interestingly, the succession of community in coalesced soils was influenced by transient competition for resources and ecological niche width, and the highest abundance and diversity of microorganisms and ARGs were found in the initial phase of coalescence, followed by a gradual succession towards the original community. With increasement of wastewater in the irrigation materials, the soil community showed a stepwise change rather than linear change. Notably, natural deposit of irrigation materials reduced their impacts on the ARGs in the coalesced soils. Findings provide new insights into the resistome coalescence during agricultural practices for reducing the spread risks of ARGs.

Fate, risk and sources of antibiotic resistome and its attenuation dynamics in the river water-sediment system: Field and microcosm study

Antibiotic resistance genes (ARGs) in rivers have received widespread attentions. Deciphering the fate and spread mechanisms of ARGs in river system can contribute to the design of effective strategies for reducing resistome risk in the environment. Although some studies have reported the prevalence and distribution of ARGs in rivers worldwide, few have systematically explored their fates, sources and risks in river water-sediment system. Also, the role of natural sunlight on the attenuation and fate of ARGs in river remains to be demonstrated. To fill the gaps, field investigation and microcosm experiment have been conducted in this study to reveal the fate, risk, source-sink relationship and attenuation dynamics of ARGs in an urban river water-sediment system, by utilizing high-throughput sequencing-based metagenomic assembly analysis and machine-learning-based source tracking tool. In all, 527 unique ARGs belonging to 29 antimicrobial classes were identified in the river. Relatively, the level of ARGs in the sediments were significantly higher (p < 0.05) than that in the waters. Variance partitioning analysis indicated the biotic and abiotic factors co-governed the riverine resistome, totally explaining 76% and 83% variations of ARGs in the waters and sediments, respectively. Microcosm experiment under natural light and dark condition showed that light induced the decay of ARGs in the waters and might promote their transfers from waters to sediments, which were also confirmed by the attenuation dynamics of bacteria in the experimental water-sediment system. Notably, the co-occurrences of ARGs with MGEs and VFs on the same contigs suggested resistome risk in the river, and relatively, the risk scores in the sediments were significantly higher (p < 0.05) than those in the waters. Source apportionment with metagenomic resistome signatures showed the Wenyu River was the most dominant contributor of ARGs in the downstream, with average contributions of 44.5% and 40.8% in the waters and sediments, respectively.

An integrated multimedia fate modeling framework for identifying mitigation strategy of antibiotic ecological risks: A case study in a peri-urban river

Antibiotics have been heavily used over the past decades, resulting in their frequent detections in rivers and increasing ecological risks. Recognizing characteristics of antibiotic ecological risks (AERs) and making effective strategies to mitigate the AERs are essential to ensure the safety of aquatic ecosystem and public health. In this study, an integrated technological framework has been proposed toward identifying management options for reducing AERs by jointly utilizing multimedia fugacity modelling and ecotoxicological risk assessment, and applied to characterize the AERs in a peri-urban river in Beijing. Specifically, a level III fugacity model has been successfully established to simulate the fate of antibiotics in the environment, and the manageable parameters have been screened out via sensitivity analysis of the model. Then the validated fugacity model has been used for scenario modellings to optimize mitigation strategies of AERs. Results show most of the antibiotics considered are frequently detected in the river, and pose medium or high risks to aquatic organisms. Relatively, the macrolides and fluoroquinolones present higher ecotoxicological risks than sulfonamides and tetracyclines. Furthermore, the mixture risk quotient and predictive equation of concentration addition suggest joint and synergistic/antagonistic effects of AERs for multiple or binary antibiotics in the environment. Largely, the concentrations of antibiotics in the river are determined by the source emissions into water and soil. Scenario modellings show the improvement of antibiotic removal rates would be considered preferentially to mitigate the AERs. Also, controlling human consumption is conducive to reducing the risks posed by tetracyclines, macrolides and trimethoprim, while controlling animal consumption would benefit the reduction for sulfonamides. Overall, the joint strategy presents the greatest reduction of AERs by reducing antibiotic consumption and together improving sewage treatment rate and antibiotic removal rate. The study provides us a useful guideline to make ecological risk-based mitigation strategy for reducing AERs in environment.

Risk characteristics of resistome coalescence in irrigated soils and effect of natural storage of irrigation materials on risk mitigation

Irrigation and fertilization are the routinely agricultural practices but also cause resistome coalescence, by which the entire microbiomes from irrigation materials invade soil microbial community, to transfer antibiotic resistance genes (ARGs) in the coalesced soils. Although studies have reported the effect of irrigation or fertilization on the prevalence and spread of ARGs in soils, risk characteristics of resistome coalescence in irrigation system remain to be demonstrated and few has shown whether natural storage of irrigation materials will reduce resistance risks. To fill the gaps, two microscopic experiments were conducted for deeply exploring resistance risks in the soils irrigated with wastewater and manure fertilizer from a perspective of community coalescence by metagenomic analysis, and to reveal the effect of natural storage of irrigation materials on the reduction of resistance risks in the coalesced soils. Results showed irrigation and coalescence significantly increased the abundance and diversity of ARGs in the soils, and introduced some emerging resistance genes into the coalesced community, including mcr-type, tetX, qacB, and an array of genes conferring resistance to carbapenem. Procrustes analysis demonstrated microbial community was significantly correlated with the ARGs in coalesced soils, and variance partitioning analysis quantified its dominant role on shaping resistome profile in the environment. Besides ARGs, abundant and diverse mobile genetic elements (MGEs) were also identified in the coalesced soils and co-existed on the ARG-carrying contigs, implying potential transfer risk of ARGs in the irrigation system. Further, the analysis of metagenome-assembled genomes (MAGs) confirmed the risk by recovering 358 ARGs-carrying MAGs and identifying the resistant bacteria that co-carried multiple ARGs and MGEs. As expected, the natural storage of irrigation water and manure fertilizer reduced about 27%-54% of ARGs, MGEs and virulence factors in the coalesced soils, thus caused the soils to move towards lower resistance risks to a certain extent.

Microcosm experiments revealed resistome coalescence of sewage treatment plant effluents in river environment

Sewage treatment plant (STP) effluents are important contributors of antibiotic resistance (AR) pollution in rivers. Effluent discharging into rivers causes resistome coalescence. However, their mechanisms and dynamic processes are poorly understood, especially for the effects of dilution, diffusion, and sunlight-induced attenuation on coalescence. In this study, we have constructed microcosmic experiments based on in-situ investigation to explore these issues. The first batch experiment revealed the effects of dilution and diffusion. The coverage of water coalesced resistomes ranged 66.26∼152.18 × /Gb and was positively correlated with effluent volume (Mann-Kendall test, p < 0.01). Principal coordinate analysis (PCoA) and source tracking analysis demonstrated that dilution and diffusion stepwise reduced AR pollution. The second batch experiment explored the temporal dynamics and sunlight attenuation on coalesced resistomes. Under natural light, the coverage and diversity of water resistomes posed decreasing trends, primarily attributed to drastic erasure of effluent traces. The proportion of effluent-specific ARGs in coalesced resistomes significantly declined over time (Spearman's r = -0.83 and -0.94 in coverage and richness). While under dark condition, the coverage and diversity increased. Sunlight radiation intensified the interactions between water and sediment resistomes, as evidenced by more shared ARGs and less dissimilarities across niches. Network analysis, metagenome-assembled genome (MAG) analysis and variation partitioning analysis (VPA) showed that microbiome controlled resistome coalescence, explaining 56.5% and 58.4% of resistomes in water and sediment, respectively. Biotic and abiotic factors synergistically explained 40% of water resistomes. This study offers a comprehensive understanding of AR transmission and provides theoretical bases for grasping AR pollution and developing effective suppression strategies.

Improvement of crAssphage detection/quantification method and its extensive application for food safety

Water-borne diseases are usually caused by the fecal-oral transmission of human fecal pathogens. Traditionally, coliforms and enterococci are widely used as indicator bacteria, but they do not allow to differentiate between human and animal fecal contamination. Owing to its presence only in the human gut environment, crAssphage has been suggested as an alternative indicator of human fecal contamination to overcome the above challenges. In this study, 139 human and 89 animal fecal samples (e.g., chicken, cow, dog, pig, pigeon, and mouse) were collected. For the rapid detection of human crAssphage in fecal samples, quantitative real-time PCR (qPCR) was performed using five different oligonucleotide primer/probe combinations. These included three previously reported oligonucleotide primer/probe combinations (RQ, CPQ056, and CrAssBP) and two newly developed combinations (ORF00018-targeting CrAssPFL1 and ORF00044-targeting CrAssPFL2). The detection rate (crAssphage-positive rate) in human fecal samples were 23.0, 30.2, 28.8, 20.1, and 30.9%, respectively, suggesting CrAssPFL2 showed the highest detection rate. Furthermore, the lowest copy numbers (436.16 copy numbers) could be detected using the CrAssPFL2 combination. Interestingly, no difference in crAssphage detection rates was found between healthy people and intestinal inflammatory patients. As expected, no crAssphage was detected in any animal fecal samples, indicating its human specificity. Furthermore, qPCR analysis of sewage samples collected from five different sewage treatment plants revealed that they were all contaminated with 10^5.71 copy numbers/mL of crAssphage on average. The simulation test of crAssphage-contaminated food samples also confirmed that the detection limit was from 10^7.55 copy numbers of crAssphage in foods. Therefore, the newly developed and optimized qPCR would be useful for the sensitive detection of crAssphage while identifying the source of human fecal contamination.

CrAssphage May Be Viable Markers of Contamination in Pristine and Contaminated River Water

Viruses are the most biologically abundant entities and may be ideal indicators of fecal pollutants in water. Anthropogenic activities have triggered drastic ecosystem changes in rivers, leading to substantial shifts in chemical and biological attributes. Here, we evaluate the viability of using the presence of crAssphage as indicators of fecal contamination in South African rivers. Shotgun analysis revealed diverse crAssphage viruses in these rivers, which are impacted by chemical and biological pollution. Overall, the diversity and relative abundances of these viruses was higher in contaminated sites compared to pristine locations. In contrast to fecal coliform counts, crAssphage sequences were detected in pristine rivers, supporting the assertion that the afore mentioned marker may be a more accurate indicator of fecal contamination. Our data demonstrate the presence of diverse putative hosts which includes members of the phyla Bacteroidota, Pseudomonadota, Verrucomicrobiota, and Bacillota. Phylogenetic analysis revealed novel subfamilies, suggesting that rivers potentially harbor distinct and uncharacterized clades of crAssphage. These data provide the first insights regarding the diversity, distribution, and functional roles of crAssphage in rivers. Taken together, the results support the potential application of crAssphage as viable markers for water quality monitoring. IMPORTANCE Rivers support substantial populations and provide important ecosystem services. Despite the application of fecal coliform tests and other markers, we lack rapid and reproducible approaches for determining fecal contamination in rivers. Waterborne viral outbreaks have been reported even after fecal indicator bacteria (FIB) were suggested to be absent or below regulated levels of coliforms. This indicates a need to develop and apply improved indicators of pollutants in aquatic ecosystems. Here, we evaluate the viability of crAssphage as indicators of fecal contamination in two South African rivers. We assess the abundance, distribution, and diversity of these viruses in sites that had been predicted pristine or contaminated by FIB analysis. We show that crAssphage are ideal and sensitive markers for fecal contamination and describe novel clades of crAss-like phages. Known crAss-like subfamilies were unrepresented in our data, suggesting that the diversity of these viruses may reflect geographic locality and dependence.

CrAssphage as an indicator of human-fecal contamination in water environment and virus reduction in wastewater treatment

Viral indicators of human-fecal contamination in wastewaters and environmental waters have been getting much attention in the past decade. Cross-assembly phage (crAssphage) is the most abundant DNA virus in human feces. Recently, the usefulness of crAssphage as a microbial source tracking and water quality monitoring tool for human-fecal contamination has been highlighted. Here, we conducted a comprehensive review on crAssphage in water, focusing on detection methodology, concentration range in various waters and wastewaters, specificity to human-fecal contamination, and reduction in wastewater treatment systems. This review highlights that crAssphage is globally distributed in wastewaters and various fecal-contaminated water bodies at high concentrations without seasonal fluctuations. CrAssphage is highly specific to human-fecal contamination and is rarely found in animal feces. It also has a good potential as a performance indicator to ensure virus reduction in wastewater treatment systems. Accordingly, crAssphage could be an effective tool for monitoring of human-fecal contamination and potential presence of fecal pathogenic microbes in environmental waters. Bridging the research gaps highlighted in this review would make crAssphage a powerful tool to support the control of water-related health risks.

Persisting antibiotic resistance gene pollution and its association with human sewage sources in tropical marine beach waters

Antibiotic-resistant bacteria (ARB) and antibiotic resistance genes (ARGs) are pollutants of worldwide concern that threaten human health and ecosystems. Anthropogenic activities and wastewater could be ARB and ARG pollution sources; however, research on ARG abundance and microbial source tracking (MST) of contamination in tropical marine waters is limited. This study examined spatiotemporal variations of six ARGs (bla_NDM, bla_TEM, bla_VIM, mcr-1, sul1, and tetQ) against the widely used antibiotic groups and a class 1 integron-integrase gene (intI1) at two Thai tropical recreational beaches (n = 41). Correlations between ARGs and sewage-specific MST markers (i.e., crAssphage and human polyomaviruses [HPyVs]) and fecal indicator bacteria (i.e., total coliforms, fecal coliforms, and enterococci) were also investigated. Bla_TEM, intI1, sul1, and tetQ were ubiquitous at both beaches (85.4-100% detection rate); intI1 was the most abundant (3-6 orders in log₁₀ copies/100 mL), followed by bla_TEM (2-4 orders), sul1 (2-3 orders), and tetQ (2-4 orders). Bla_NDM was found in 7.3% (up to 4 orders), and no mcr-1 was detected. Interestingly, bla_VIM was prevalent at one beach (2-5 orders; n = 17), but found in only one sample at the other (4 orders). Temporal, but not spatial, differences were noticed; bla_TEM was at higher levels in the wet season. IntI1 correlated with sul1 and tetQ (Spearman's rho = 0.47-0.97), suggesting potential horizontal gene transfer. CrAssphage, but not HPyVs, correlated with intI1, sul1, and tetQ (Spearman's rho = 0.50-0.74). Higher numbers of ARGs tended to co-occur in samples with higher crAssphage concentrations, implying sewage contribution to the marine water, with a persisting ARG background. This study provides insight into the ARG pollution status of tropical coastal waters and suggests crAssphage as a proxy for ARG pollution, which could facilitate effective management policies to minimize ARG dissemination in marine environments.