
==== Front
Environ Sci Technol
Environ Sci Technol
es
esthag
Environmental Science & Technology
0013-936X
1520-5851
American Chemical Society

39229966
10.1021/acs.est.4c03726
Article
Metagenomic Profiling of Internationally Sourced Sewage Influents and Effluents Yields Insight into Selecting Targets for Antibiotic Resistance Monitoring
https://orcid.org/0000-0002-1579-155X
Garner Emily *†
Maile-Moskowitz Ayella ‡
https://orcid.org/0000-0002-5326-5710
Angeles Luisa F. §
https://orcid.org/0000-0002-4101-5095
Flach Carl-Fredrik ∥
https://orcid.org/0000-0001-6512-7713
Aga Diana S. §
Nambi Indumathi ⊥
https://orcid.org/0000-0002-5496-0328
Larsson D. G. Joakim ∥
https://orcid.org/0000-0002-5651-5906
Bürgmann Helmut #
https://orcid.org/0000-0003-1148-4322
Zhang Tong ¶
https://orcid.org/0000-0003-2654-5132
Vikesland Peter J. *‡
https://orcid.org/0000-0002-3191-6244
Pruden Amy *‡
† Wadsworth Department of Civil and Environmental Engineering, West Virginia University, Morgantown, West Virginia 26505, United States
‡ Department of Civil and Environmental Engineering, Virginia Tech, Blacksburg, Virginia 24061, United States
§ Department of Chemistry, University at Buffalo, Buffalo, New York 14260, United States
∥ Institute of Biomedicine, Department of Infectious Diseases, Centre for Antibiotic Resistance Research in Gothenburg (CARe), University of Gothenburg, Västra Götaland, SE-405 30 Gothenburg, Sweden
⊥ Department of Civil Engineering, Indian Institute of Technology, Madras, Chennai 600036, India
# Eawag: Swiss Federal Institute of Aquatic Science and Technology, Kastanienbaum CH-6047, Switzerland
¶ Department of Civil Engineering, The University of Hong Kong, Pokfulam 999077, Hong Kong
* Email: emily.garner@mail.wvu.edu.
* Email: pvikes@vt.edu.
* Email: apruden@vt.edu.
04 09 2024
17 09 2024
58 37 1654716559
22 04 2024
27 08 2024
26 08 2024
© 2024 The Authors. Published by American Chemical Society
2024
The Authors
https://creativecommons.org/licenses/by/4.0/ Permits the broadest form of re-use including for commercial purposes, provided that author attribution and integrity are maintained (https://creativecommons.org/licenses/by/4.0/).

It has been debated whether wastewater treatment plants (WWTPs) primarily act to attenuate or amplify antibiotic resistance genes (ARGs). However, ARGs are highly diverse with respect to their resistance mechanisms, mobilities, and taxonomic hosts and therefore their behavior in WWTPs should not be expected to be universally conserved. We applied metagenomic sequencing to wastewater influent and effluent samples from 12 international WWTPs to classify the behavior of specific ARGs entering and exiting WWTPs. In total, 1079 different ARGs originating from a variety of bacteria were detected. This included ARGs that could be mapped to assembled scaffolds corresponding to nine human pathogens. While the relative abundance (per 16S rRNA gene) of ARGs decreased during treatment at 11 of the 12 WWTPs sampled and absolute abundance (per mL) decreased at all 12 WWTPs, increases in relative abundance were observed for 40% of the ARGs detected at the 12th WWTP. Also, the relative abundance of mobile genetic elements (MGE) increased during treatment, but the fraction of ARGs known to be transmissible between species decreased, thus demonstrating that increased MGE prevalence may not be generally indicative of an increase in ARGs. A distinct conserved resistome was documented in both influent and effluent across samples, suggesting that well-functioning WWTPs generally attenuate influent antibiotic resistance loads. This work helps inform strategies for wastewater surveillance of antibiotic resistance, highlighting the utility of tracking ARGs as indicators of treatment performance and relative risk reduction.

Metagenomic sequencing across a transect of wastewater treatment plant influents and effluents reveals which antibiotic resistance genes (ARGs) tend to decrease or amplify during treatment, informing surveillance and risk reduction strategies.

antibiotic resistance
resistome
wastewater treatment
metagenomics
monitoring
Office of International Science and Engineering 10.13039/100000089 1545756 VetenskapsrÃ¥det 10.13039/501100004359 2022-00945 VetenskapsrÃ¥det 10.13039/501100004359 2018-05771 Svenska ForskningsrÃ¥det Formas 10.13039/501100001862 2021-00922 Svenska ForskningsrÃ¥det Formas 10.13039/501100001862 2018-00833 Schweizerischer Nationalfonds zur FÃ¶rderung der Wissenschaftlichen Forschung 10.13039/501100001711 186531 Schweizerischer Nationalfonds zur FÃ¶rderung der Wissenschaftlichen Forschung 10.13039/501100001711 167116 Water Research Foundation 10.13039/100005612 4961 U.S. Environmental Protection Agency 10.13039/100000139 84082801 Office of Advanced Cyberinfrastructure 10.13039/100000105 2004751 document-id-old-9es4c03726
document-id-new-14es4c03726
ccc-price
==== Body
pmcIntroduction

Bacterial resistance to antibiotics is a serious global health threat and the role of water, sanitation, and hygiene in controlling its spread is increasingly recognized.1 Expansion and improvement of sewage collection and treatment via a wastewater treatment plant (WWTP) is one means of curbing the spread of resistant microorganisms via environmental pathways.2 WWTPs often employ biological treatment, such as activated sludge, wherein aeration is applied to stimulate organic matter biodegradation by flocs of heterotrophic bacteria, which are subsequently removed via settling with the aqueous effluent returned to the environment. WWTPs are highly effective in reducing the concentrations (cells/L) of bacteria leaving the plant, including pathogens. Removal of antibiotic resistant bacteria (ARB), such as Escherichia coli and Enterococcus spp., by activated sludge treatment have likewise been reported in the range of 2–5-log,3,4 but ARB of clinical concern can persist in WWTP effluents and pose health risks.5 There remains a lack of consensus about whether WWTPs act to increase versus attenuate the relative abundance, i.e., abundance normalized to a marker of the total bacterial population such as the 16S rRNA gene, of ARB or antibiotic resistance genes (ARGs).5−7 Shifts in relative abundances of ARB and ARGs could reflect selective pressures occurring within the WWTP, but also will be strongly shaped by microbial ecological shifts that reflect aeration of the secondary bioreactor.

Previous research has demonstrated that the relative abundance of some ARGs can be reduced during wastewater treatment, such as ermB,8blaOXA-58,9tetC,9 and tetM.9 Heterotrophic bacteria resistant to vancomycin, gentamicin, erythromycin, cephalexin, tetracycline, and sulfadiazine have been reported to decline.10 Yet, other studies have found minimal impact or even enrichment of ARGs, such as vanA,8blaVIM,8 and blaSHV-34.9 Ju et al. observed increased relative abundance of 12 resistance classes, but also biocide and heavy metal resistance classes.11 Another study reported increases in drug-resistant Acinetobacter spp.12 Several studies have documented differences in removal or enrichment of ARB and ARGs among WWTPs.9,11 Further, critical differences in social, clinical, regulatory, and environmental factors around the world are expected to influence the types and magnitudes of ARB and ARGs in influent sewage and thus the efficacy of the WWTP as a barrier for ARB/ARG removal.13−15

There are several drivers that may enrich ARB and ARGs within a WWTP. For example, horizontal gene transfer (HGT) facilitates transfer of ARGs to new cell hosts that grow at higher rates in the activated sludge environment. HGT is particularly relevant in the parts of the WWTP with high bacterial cell density (e.g., activated sludge)16 and might be exacerbated by sewage coconstituents; such as organic compounds, residual antibiotics, and heavy metals, which have been shown to increase rates of conjugation or natural transformation.17−20 Biological wastewater processes have been characterized as housing a high co-occurrence of bacterial species capable of donating their ARGs via HGT, along with corresponding mobilizing agents, such as insertion sequences.21 Selection, HGT, or induction of mutations that increase resistant microbiota have been associated with constituents commonly present in wastewater; such as antibiotics, metals, herbicides, and antidepressants.22−25 Shifts in microbial community structure have been reported to be a key driver of shifts in ARGs in soil,26 and this may also be an important driving force in WWTPs.11

High throughput metagenomic DNA sequencing can be used to efficiently profile and compare known ARGs across WWTPs.13,14,27 However, with over 3500 different ARGs documented to date,28 efforts are needed to establish a baseline understanding of which ARGs are present in the influent versus effluent, which tend to attenuate versus increase during treatment, and whether influent ARG composition significantly affects treatment outcomes. Such knowledge will help improve our understanding of the extent to which WWTPs effectively act as a barrier to the spread of antibiotic resistance as well as to identify anomalies that could indicate problems with treatment or flag emerging public health concerns.

The objective of this study was to broadly categorize ARGs according to their reduction or enrichment during wastewater treatment. The following key hypotheses were tested: (1) relative (per 16S rRNA gene) and absolute (per mL) abundance of total ARGs decrease as a result of wastewater treatment; (2) the final effluent resistome composition mirrors regional patterns previously documented in the influent;14 (3) removal efficiencies for individual ARGs vary widely across and within resistance classes and WWTPs; (4) wastewater treatment reduces the abundance of ARGs carried by bacteria causing human disease; (5) biological wastewater treatment increases ARG diversity in the final effluent relative to the influent; (6) wastewater treatment leads to a reproducible shift in resistome composition across geographically diverse WWTPs. By leveraging metagenomic sequencing across 12 diverse WWTPs from six countries in three continents generated using a consistent field collection and analysis protocol, this study provides insight into the extent to which treatment impacts on the resistome are conserved globally across facilities and proposes specific metrics for comparing treatment efficacy and identifying potential human health hazards. While previous studies have used metagenomic sequencing to comprehensively profile the antibiotic resistome in sewage13,29 or to monitor changes across individual treatment trains,4,27 this study is the first to examine the effects of biological wastewater treatment on the resistome across a broad range of international facilities.

Materials and Methods

Site Description and Sample Collection

Sampling was conducted at two WWTPs from each of six countries (India (IND), Hong Kong, China (HKG), Philippines (PHL), United States of America (USA), Switzerland (CHE), and Sweden (SWE)), between March 2016 and January 2017. All WWTPs employed secondary treatment, with capacities ranging from 2.6 to 66 million gallons per day. Both Hong Kong plants (HKG-1, HKG-2) and one USA plant (USA-1) employed UV-disinfection of effluent; one Philippines WWTP (PHL-1), one USA WWTP (USA-2), and both India plants (IND-1 and IND-2) disinfected using free chlorine; and one Swiss WWTP (CHE-2) disinfected via ozonation followed by sand filtration. The remaining four WWTPs did not disinfect final effluent. Additional characteristics of the sampled WWTPs have been published previously,14 along with a detailed examination of their influent resistomes, and are described in additional detail in Table S1. Sample collection and processing were conducted using standardized protocols validated for sample preservation and stability during international shipment.30,31 Influent and final effluent grab samples for molecular analysis were collected at each WWTP in sterile polypropylene containers and samples for antibiotic analysis were collected in acid-washed, baked amber glass bottles.

DNA Extraction and Quantitative Polymerase Chain Reaction

Samples were collected, processed, preserved and shipped according to Li et al.31 Within 12 h of collection, triplicate aliquots of each sample were concentrated onto 0.22 μm mixed cellulose ester membranes (Millipore, Billerica, MA) until clogging. In cases where clogging occurred prior to the maximum amount of water collected passing through the filter, less than the full volume was filtered. The actual volume filtered was recorded and used for normalization. Filters were preserved in 50% ethanol and shipped to Virginia Tech on ice packs. Upon arrival, filters were frozen at −20 °C until DNA extraction. Filters were aseptically torn into 1 cm2 pieces using sterile forceps and transferred to extraction tubes. DNA was extracted using the FastDNA SPIN Kit for Soil (MP Biomedicals, Solon, Ohio). Quantitative polymerase chain reaction (qPCR) was used to quantify 16S rRNA genes and the sul1 sulfonamide ARG in triplicate reactions using previously published primers (Table S2).32,33 Triplicate standard curves of 10-fold serial diluted standards of each target gene ranging from 101 to 107 gene copies/μL were included on each 96-well plate, along with a triplicate negative control.

Metagenomic Sequencing and Analysis

Composite samples were prepared by pooling triplicates by equal DNA mass. Composites were prepared for sequencing using TrueSeq library preparation (Illumina, San Diego, CA) and sequenced on an Illumina HiSeq 2500 using 2 × 100 paired-end reads at the Virginia Tech Biocomplexity Institute Genomic Sequencing Center (Blacksburg, VA). Reads were uploaded to MetaStorm,34 where quality filtering was performed using Trimmomatic according to default parameters.35 Reads were assembled in MetaStorm using IDBA-UD.36 Reads and scaffolds were annotated using default parameters in MetaStorm against the Comprehensive Antibiotic Resistance Database (CARD; version 2.0.1) to identify ARGs,37 the Silva rRNA database for 16S rRNA genes,38 and the ACLAME database (version 0.4) for mobile genetic elements (MGEs).39 Taxonomy was assigned to reads using MetaPhlan2.40 To calculate relative ARG abundance, ARG read counts were normalized to 16S rRNA gene counts.41 To calculate absolute abundance, 16S rRNA gene normalized abundances were multiplied by 16S rRNA gene abundances per mL, as measured by qPCR. Mobile ARGs were defined based on inclusion in the ResFinder database (version 2.1),42 and were manually curated from CARD annotations. Assembled scaffolds were uploaded to NanoARG43 for identification of ARGs, MGEs, and metal resistance genes using default parameters. Metagenomes are publicly available via National Center for Biotechnology Information’s BioProject PRJNA527877.

Antibiotic Analysis

Sample analysis for antibiotics was conducted as previously described and antibiotic concentrations in these samples have been previously published.30 Briefly, wastewater samples (0.5 L) were acidified and filtered using 0.45 μm glass microfiber filters to remove microorganisms and particulate matter. Na2EDTA (2 mL, 5% v/v in water) and surrogate standards (50 μL of 1000 μg/L surrogate mix solution) were added to each sample. Solid phase extraction was performed by conditioning Oasis HLB cartridges with acetonitrile and deionized water before the water samples were loaded at a rate of 3–5 mL/min. Cartridges were dried under vacuum and shipped to the University at Buffalo for elution and liquid chromatography with tandem mass spectrometry (LC–MS/MS) analysis using an Agilent 1200 LC system (Palo Alto, CA).

Statistical Analysis

Differences in ARG abundance between groups were tested using a Wilcoxon rank sum test in R.44 Percent mobility was defined as described by Ju et al., determined by calculating the percentage of scaffolds containing both ARGs and MGEs out of the total number of ARG-containing scaffolds.11 Correlations between antibiotics and ARGs were tested using a Spearman rank sum test in R.44 To compare resistome profiles across samples, ARG abundances were standardized to a percentage of the total ARGs from each sample and a Bray–Curtis resemblance matrix was generated in Primer-E (version 6.1.13). Nonmetric multidimensional scaling (NMDS), analysis of similarities (ANOSIM), and similarity clustering were conducted in Primer-E. Shannon diversity of microbial genera and ARGs, as well as symmetric Procrustes and PROTEST analysis were conducted in R using vegan.45 Linear discriminant analysis Effect Size (LEfSe) was applied to identify ARGs that best explain differences between sample groups.46 Figures were produced using ggplot247 and Microsoft Excel.

Results

Detection of ARGs and Quantitative Assessment of Metagenomic Data

Metagenomic sequencing yielded over 326 million paired-end reads across all samples, with an average of 1.36 × 106 reads per sample (range: 0.92–1.78 × 106). Across the data set, 1079 different ARGs were identified, with up to 554 different ARGs detected in a single sample (SWE-2 influent; Table S3). The most abundant ARG classes in the influent, on average, were multidrug, macrolide-lincosamide-streptogramin (MLS), beta-lactam, tetracycline, and aminoglycoside (Figure 1A). Similarly, many of the same classes of ARGs dominated the effluent resistome, of which multidrug, MLS, beta-lactam, aminoglycoside, and peptide were most abundant. Relative abundance of sul1 normalized to 16S rRNA genes correlated with absolute abundances of sul1 genes per unit volume measured using qPCR (R2 = 0.48, p = 0.031) (Figure S1). The strong alignment of these results supports the use of subsequent quantitative comparisons made via metagenomic data. However, results are presented in terms of both relative abundance (per 16S rRNA gene copies) and absolute abundance (per mL; Figure 1) to provide insight into the effect of changing overall microbial concentrations on the observed resistome.

Figure 1 ARGs in WWTP influent and effluent based on (A) relative abundance annotated against the CARD database via MetaStorm and normalized to 16S rRNA genes, (B) absolute abundance of ARGs annotated against the CARD database via MetaStorm and normalized to sample volume, (C) potentially mobile ARGs, identified as those included in the ResFinder database of acquired ARGs, and (D) ubiquitous ARGs (i.e., detected in all samples).

Impact of Treatment on Total ARG Abundance

While many ARG classes remained dominant from the influent to the final effluent, the overall relative abundance of the summed total ARGs (total ARGs) decreased at 11 of 12 WWTPs (p = 0.01 across the data set). Only one WWTP, PHL-1, exhibited an increase in relative total ARG during treatment. Absolute abundance (per mL) of total ARGs exhibited a similar trend to that observed for relative abundance, decreasing at all 12 WWTPs by ∼0.5–3.5 log (Figure 1B). Hereafter, we focus on relative abundance as an indicator of the tendency of ARGs to be enriched or attenuated across WWTP microbial communities.

Patterns in the Fate of Individual ARGs

The relative abundance of each ARG is provided in Table S3. While both the relative and absolute abundance of total ARGs decreased consistently during treatment, many individual ARGs increased in relative abundance from the influent to final effluent (Figure 2A). Notably, the greatest percentage of ARGs (40.0%) increased at the PHL-1 WWTP, which was also the only WWTP that did not yield a decrease in relative total ARG abundance. Increases in ARG abundance at the remaining WWTPs ranged from 8.2 to 32.8%: 32.8 and 22.7% of detected ARGs increased in Swiss WWTPs, 32.5 and 26.3% of ARGs in Hong Kong WWTPs, 21.1 and 25.6% in Indian WWTPs, 28.2% in the remaining Philippines WWTP, 24.3 and 31.1% in Swedish WWTPs, and 24.4 and 8.2% in USA WWTPs. Increases in abundance of these enriched ARGs were striking in many instances, sometimes reaching >3-log enrichment (Figure 2B). Across all ARGs and all WWTPs, there was a notable bimodal distribution in the log change in relative abundance, as evidenced by clustering of datapoints in Figure 2B,C: 45.4% of data points indicated between −4.0 and −3.0 log change in relative abundance of ARGs, while 40.3% were between −1.0 and +1.0 log. This bimodal distribution was driven by two common scenarios: (A) ARGs abundant in influent wastewater, but undetectable in the final effluent, and (B) ARGs that were undetectable in either the influent or final effluent and only detected at low levels in the other.

Figure 2 (A) Relative abundance of all detected ARGs in each WWTP’s influent vs final effluent. (B) Log change in relative abundance of individual ARGs during treatment at each WWTP. (C) Log change in relative abundance of individual ARGs during treatment by antibiotic resistance class.

Some notable trends in reduction or enrichment of ARGs were observed as a function of resistance class (i.e., the class of antibiotics to which ARGs encode resistance; Figure 3A). For example, the relative abundance of total tetracycline and MLS ARGs decreased with treatment across all WWTPs. Aminoglycoside, beta-lactam, peptide, phenicol, and quinolone ARGs were reduced in relative abundance at nearly all WWTPs. On the contrary, aminocoumarin, multidrug, rifamycin, sulfonamide, and trimethoprim ARGs increased in relative abundance at the majority of WWTPs. While these general trends are notable, it is important to point out that substantial variations in changes were observed among individual ARGs within each class (Figure 2C).

Figure 3 Change in abundance following treatment of (A) ARGs according to class of resistance and (B) select clinically relevant ARGs.

Across all WWTPs, there was a significant increase in ARG diversity from the influent to the final effluent (Wilcox, p = 0.0005; Figure S2A), consistent with the biological reactors harboring a distinct microbial community that carry a distinct array of ARGs that are additionally detected in final effluents. It was noted that some ARGs that were highly abundant in the influent were among those that experienced the greatest overall reduction during treatment. Twenty-two ARGs fell within both the 15th percentile of the highest influent concentrations and the 15th percentile of greatest reductions in relative abundance during treatment (Figure S3): two aminoglycoside (APH(3″)-lb, APH(6)-ld), two beta-lactam (CfxA4, CfxA6), six MLS (ErmB, ErmF, mefA, mel, mphD, msrE), four multidrug (adeJ, adeK, adeM, CRP), one peptide (pmrE), one phenicol (cat), two quinolone (qacH, QnrS2), one sulfonamide (sul1), and three tetracycline (tet(39), tetQ, tetW) ARGs. While decreases in sul1 concentrations were notable, decreasing by 0.00981 ± 0.02380 gene copies per 16S rRNA gene across all plants, the noted increase in sulfonamide ARGs at the class level (Figure 1A,B) was largely driven by sul2 and sul4, with increases of 0.01435 ± 0.03717 and 0.00461 ± 0.00810 gene copies per 16S rRNA gene, respectively.

Conserved and Differential Aspects of Influent and Effluent Resistomes

The relative similarities in ARG profiles across globally distributed WWTPs is noteworthy (Figures 1A,B). The core resistome was evaluated to identify the portion of the resistome that is “ubiquitous” across all influent and effluent sites (Figure 1D). The core influent resistome included 233 subtype ARGs shared across all influent samples, while the core effluent included 51 subtype ARGs. Overall, 49 ARGs at the subtype level were ubiquitous across all influent and effluent samples. This group was overwhelmingly dominated by multidrug (n = 31) ARGs. There were no ARGs at the subtype level that were detected exclusively in the influent or final effluent.

LEfSe analysis was applied to determine which ARGs were the key drivers of differences in the influent and effluent resistomes (Figure S4). Two tetracycline (tetQ, tetA(46)), three beta-lactam (CfxA6, CfxA4, CTX-M-107), two MLS (mphB, ErmG), one quinolone (patB), one fosfomycin (FosA5), four multidrug (lsaB, lsaC, mdtO, ramA), and one glycopeptide (vanSG) ARG were key features of the influent resistomes. Key features of the final effluent resistome included two multidrug (mexN, oprA), one rifamycin (rphA), and one trimethoprim ARG (dfrA2d).

NMDS was used to visualize differences in the resistome composition of each sample (Figure 4). Influent and effluent samples formed distinct clusters, sharing 55% similarity (ANOSIM, R = 0.926, p < 0.001). Effluent resistomes formed a tighter cluster than influent resistomes, with the PHL-1 effluent an apparent outlier. At a high level, the shared ecology of the biological treatment environment shifted the resistome composition from influent to effluent in a similar manner, regardless of geographical location and initial influent resistome composition.

Figure 4 NMDS plot depicting the resistome of influent versus final effluent samples from twelve global WWTPs. Similarity was calculated according to the Bray–Curtis metric. Samples clustered at 55% similarity.

Fate of Clinically-Relevant ARGs in the WWTPs

Key clinically relevant ARGs were selected to assess their removal by the 12 WWTPs (Figure 3B). While OXA-type carbapenemases were the most common clinically relevant ARG detected, they were effectively removed at all WWTPs sampled, except the PHL-1 WWTP. GES, SHV, and TEM-type beta-lactam ARGs were also typically removed during treatment, but removal patterns of CARB and CTX-M beta-lactam ARGs varied among WWTPs. New Delhi metallo (NDM) beta-lactamase ARGs were not detected in any influent or effluent sample. MCR-1, which confers resistance to colistin, was detected in both influent and effluent, though at very low abundance (≤0.00144 copies/16S rRNA gene). qnrS, a plasmid-mediated quinolone ARG, was effectively removed at 11 of 12 WWTPs, with USA-1 the exception.

A total of 5,740,060 scaffolds were generated, averaging 735 bp in length and 179,377 scaffolds per sample (range: 64,576–256,387). Of these, 1214 scaffolds were associated with the WHO list of priority antibiotic resistant pathogens or the ESKAPE pathogens (i.e., Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa, and Enterobacter species), which were further examined for carriage of ARGs using NanoARG (Figure 5).

Figure 5 Abundance of ARGs, by class, identified on assembled scaffolds aligning with pathogen reference genomes using default parameters in NanoARG.43

Multidrug ARGs were found on a substantial percentage of the scaffolds associated with pathogens in the influent and effluent (respectively) for A. baumannii (25, 29%), Haemophilus influenza (33, 40%), K. pneumoniae (35, 25%), Neisseria gonorrheae (25, 100%), P. aeruginosa (32, 41%), and S. aureus (30, 100%). Among the effluent samples, only two scaffolds were annotated as N. gonorrheae and three as S. aureus, so the 100% carriage of multidrug resistance is likely an artifact of low detection rate. Among the influent scaffolds, 28% annotated as E. faecium contained glycopeptide ARGs, a pathogen-resistance combination of high clinical concern. In the effluent, zero of the four scaffolds annotated as E. faecium contained glycopeptide ARGs. In contrast, beta-lactam resistance increased from influent to effluent on A. baumannii (from 18 to 31%) and P. aeruginosa (from 7 to 9%) scaffolds.

All scaffold annotations are provided in Table S4. Scaffolds were further analyzed if they originated from a putative pathogen and contained three or more ARGs, MGEs, metal resistance genes, and, in some cases, repeated detection of conserved regions across multiple samples (Figure 5). Four such assembled regions were derived from A. baumannii, two of which conferred resistance to multiple antibiotics and were found in the final effluent of SWE-1 (Figure S5A,D). One scaffold also included multiple transposase genes and genes conferring resistance to silver and zinc. The remaining A. baumannii scaffolds (Figure S5B,C) were found in influent samples collected from SWE-1 and PHL-1 and carried four different ARGs each. One assembled region derived from E. faecium was particularly noteworthy as the same region was detected in influent samples from CHE-2, SWE-1, SWE-2, and USA-1 (Figure S5E). The region confers resistance to bacitracin, a peptide antibiotic, as well as glycopeptide antibiotics and zinc. This region also contained a conjugative transposon gene. The recurrent detection of this genetic region in different countries suggests that it may be a strain or genetic region globally carried by human populations. The same genetic region was not detected in any final effluent sample. A genetic region associated with H. influenzae was derived from the influent of USA-1 (Figure S5F) and carried a transposase as well as ARGs conferring resistance to beta-lactam and aminoglycoside antibiotics. Finally, a genetic region found on a scaffold originating from P. aeruginosa in a HKG-1 influent sample (Figure S5G) carried ARGs associated with MLS, peptide, aminoglycoside, and glycopeptide resistance, as well as arsenic resistance, recombinase, and transposase genes.

Evaluation of Potential Factors Driving Resistome Shifts from Influent to Effluent

Potential for HGT

The portion of potentially mobile ARGs ranged from 16.8 to 57.1% of the overall resistome of each sample (Figure 1C). This mobile resistome was dominated by beta-lactam (23.2% of samples on average), MLS (23.1%), and tetracycline (14.6%) ARGs and significantly decreased in relative abundance as a result of treatment (p ≤ 0.0001), from 1.10 to 0.49 ARGs per 16S rRNA gene, on average. In influent samples, the most abundant mobile ARGs were msrE (MLS), tetQ (tetracycline), mphD (MLS), cat (phenicol), sul1 (sulfonamide), ErmF (MLS), CfxA4 (MLS), CfxA6 (MLS), and qnrS (quinolone), in descending order. In final effluent, the most abundant mobile ARGs were sul1 (sulfonamide), dfrB6 (trimethoprim), sul2 (sulfonamide), aadA11 (aminoglycoside), cat (phenicol), msrE (MLS), tet(C) (tetracycline), and ErmF (MLS). MGEs were significantly enriched in final effluent samples compared to the influent (Wilcox, p = 0.002316) (Figure S6). This is in contrast to the above observation that the mobile fraction of the resistome was reduced during treatment.

The frequencies of detection of all ARGs on MGE-associated scaffolds are presented in Table S5. The ARGs most frequently found on MGE-associated scaffolds were macB (MLS, 13,014 scaffolds), TaeA (pleuromutilin, 5565 scaffolds), evgS (multidrug, 5372 scaffolds), and tetA(48) (tetracycline, 4887 scaffolds). Clinically relevant ARGs found on MGE-associated scaffolds included MCR-1 (peptide, 73 scaffolds), qnrS (quinolone, 18 scaffolds), and a variety of OXA ARGs (beta-lactam, 479 scaffolds). sul1, a sulfonamide ARG commonly targeted as an indicator of mobile resistance, was found on 62 MGE-associated scaffolds. Percent mobility (i.e., percentage of scaffolds containing both ARGs and MGEs out of the total number of ARG-containing scaffolds) among samples ranged from 9.9 to 13.6%, with no significant difference between influent and final effluent or between geographic regions (Table S6). Excluding ARGs with only a single scaffold occurrence, the ARGs with the highest percent mobility included oleB (MLS, 76.1%), srmB (MLS, 74.0%), carA (MLS, 69.1%), CfxA2 (beta lactam, 68.8%), vgaE (multidrug, 67.8%), sul1 (sulfonamide, 66.7%), and msrC (MLS, 66.7%). Other ARGs of note with high percent mobility include tetW (tetracycline, 39.2%), sul2 (sulfonamide, 34.3%), TEM (beta-lactam, 16.7%), QnrS6 (quinolone, 14.3%), OXA (beta lactam, 12.3%), MCR-1 (peptide, 11.8%), and VIM (beta lactam, 11.1%) (Table S7).

Potential for Selection by Antibiotics During Treatment

Concentrations of sulfonamide, macrolide, quinolone, and tetracycline antibiotics were determined from influent and effluent samples collected at each WWTP (Figures S7 and S8).30 Total antibiotic concentrations in influent varied widely from 171 (SWE-2) to 47,978 ng/L (IND-1). Macrolide antibiotics were particularly prevalent in India and Hong Kong WWTPs, tetracyclines were most abundant in Hong Kong WWTPs, and the USA and Swiss WWTPs were dominated by quinolones and sulfonamides. Treatment reduced total antibiotics at 11 of 12 WWTPs, with total removal ranging from 30 to 98%. Total antibiotic concentrations increased at one WWTP, PHL-2, from 1918 to 2023 ng/L, though these concentrations were much lower than those measured in the influent at other Asian WWTPs.

To examine the potential for coselection, co-occurrence of ARGs on scaffolds was examined. The complete list of ARGs co-occurring on assembled scaffolds is presented in Table S8. vanR and vanS were the most common ARGs to occur on the same scaffold (194 co-occurrences), which is consistent with these ARGs being components of a vancomycin-resistance operon. A variety of co-occurrences of ARGs that are not part of the same operon were observed, which can potentially be coselected by selective agents. Many of the most abundant of these co-occurred with vanR and vanS. cpxR (multidrug), PvrR (aminoglycoside), and mtrA (multidrug) were associated with vanS (42, 34, 26 co-occurrences respectively) and baeS (multidrug), smeS (multidrug), PvrR (aminoglycoside), and cpxA (multidrug) were associated with vanR (31, 30, 28, 25 co-occurrences respectively). This suggested that resistance to glycopeptides may potentially be coselected for by a variety of other antibiotic compounds.

Influence of Microbial Community on Resistome

The microbial community experienced marked shifts as a result of treatment that were largely conserved across the sampled WWTPs (Figure S9A). Influent samples were dominated by Gammaproteobacteria, Clostridia, and Bacteriodia, while final effluent samples were primarily dominated by Betaproteobacteria and Gammaproteobacteria. There was a moderate concordance between the structure of the microbial community and that of the resistome (PROTEST: Procrustes Sum of Squares (m122) = 0.7326, correlation in symmetric Procrustes rotation = 0.5171, significance = 0.004.; Figure S10). In addition, numerous specific taxonomic classes were significantly and positively correlated with several ARG classes (Figure S9B). While ARG diversity consistently increased from the influent to the effluent of each WWTP, this trend did not appear to be driven by a change in microbial diversity, as there was no consistent increase in microbial diversity (Figure S2B).

Discussion

Metagenomic sequencing provided insight into the occurrence and fate of a wide array of ARGs across 12 globally distributed WWTPs in six countries. Detailed analysis aided in delineating typical responses of specific ARGs and ARG classes to biological treatment. Consistent with our hypothesis (1), total ARG relative abundance decreased in 11 of 12 WWTPs sampled, and corresponding absolute abundances decreased at all 12 WWTPs. However, in alignment with our hypothesis (3), the removal efficiencies for individual ARGs varied widely across and within resistance classes and WWTPs. Consistent with our hypothesis (4), clinically relevant ARGs generally decreased during treatment at the majority of WWTPs, consistent with the findings of previous studies.4,27 Others have recently reported striking trends in resistome compositions of raw sewage collected globally,13,14 noting abundance of total and clinically relevant ARGs are higher in Asian than US/European sewages. The present study expands this finding by noting the same pattern is retained in the corresponding WWTP final effluents, even while globally sourced effluent resistomes as a group were distinct from influents, which confirms our hypothesis (2). While this study did not address temporal variability in the resistome composition at each site, several previous studies have demonstrated that sewage resistomes are relatively stable over time, with minimal variation in the broader resistome composition during the short-term.27,29

WWTP PHL-1 was consistently noted as an outlier herein. For example, relative abundance of total ARGs as well as several classes of resistance that decreased across the majority of WWTPs increased during treatment at PHL-1, including aminoglycoside, beta-lactam, fosfomycin, glycopeptide, and phenicol resistance. Similarly, beta-lactam ARGs of GES, OXA, SHV, and TEM types all increased during treatment at PHL-1, despite decreases noted at the majority of studied plants. This finding highlights the value in examining changes in relative abundance at multiple levels, including the overall resistome, among antibiotic resistance classes, and at the individual gene level, as PHL-1 was not noted as an outlier when examining the overall resistome (Figure 1). One possible explanation for the atypical effect of PHL-1 on ARGs was that this facility employed an attached growth biofilm reactor for secondary treatment. All other facilities relied primarily on suspended growth processes, with the only other attached growth process examined being a trickling filter treating approximately 30% of the total flow at SWE-1. The ecological processes governing the microbial community composition and ARG profile in biofilm-based systems are unique from activated sludge processes due to factors such as sorption of antibiotics and micropollutants onto the biofilm matrix48 and stratification of ARB and ARGs in biofilms with filter depth.49

While there was wide variability in the relative and absolute abundance of individual ARGs across the 12 WWTPs, there were also notable similarities. Consistent with our hypothesis (6), the overall resistome shifted in a reproducible manner from the influent to effluent that was conserved across all global sites. Of the 1079 ARGs detected overall, 49 were detected consistently across all samples. These 49 ARGs have limited clinical relevance, as none of them are categorized as Rank 1 ‘current threats’ or Rank 2 'future threats'.50 There are likely two drivers for the ubiquity of ARGs that persist in the effluent: (1) they are consistently found in the human gut and are poorly removed via typical secondary wastewater treatment processes or (2) they are widespread in water and other natural environments. Additional study is needed to comprehensively compare WWTP influent and effluent to other relevant environments, including those local to each WWTP, to differentiate which of the ubiquitous ARGs fall into each group. ARGs in the first group could serve as conserved indicators for monitoring the effectiveness of wastewater treatment processes or for assessing antibiotic resistance pollution associated with WWTP effluent. ARGs in the second group would be poorly suited for such a purpose due to their ubiquity even in environments not impacted by WWTPs or human fecal pollution.

Multidrug efflux ARGs dominated the total ARG profiles (Figure 1A) and effluent trends more strongly reflected influent trends when the analysis focused on the Resfinder database of mobile ARGs, which largely excluded this class (Figure 1C). Multidrug efflux pump ARGs are highly conserved within species, due in part to their location on the chromosome, and their ability to serve the cell not only through efflux of antibiotics, but also through efflux of heavy metals and other toxic compounds, secretion of virulence factors, and transport of compounds such as bile salts and fatty acids across the cell membrane.51 Given the wide utility of these ARGs, they are found in many environments and their detection is not likely concerning for human health, thus making them poor monitoring targets.

Despite the tendency of taxonomic diversity to decrease (Figure S2), the diversity of ARGs from influent to final effluent increased at all 12 WWTPs, confirming our hypothesis (5). This observation is consistent with the idea that biological processes, such as activated sludge, which is known to harbor a highly diverse range of ARGs,52,53 can seed additional ARGs to the final effluent. However, diversity indices represent only high-level trends, requiring closer examination of the behavior of individual ARGs and ARG classes. In particular, 22 ARGs were noted to be highly abundant in influent yet also experienced the greatest overall removal during treatment (Figure S3). Thus, in developing strategies for WWTP surveillance, such ARGs should be included to confirm expected treatment performance, while including ARGs that can be prone to persist or increase during treatment.

While the limited enrichment of only a handful of ARGs observed herein is inconsistent with the notion that WWTPs are a “hotspot” for the uncontrolled spread of ARGs, it nevertheless highlights the potential for a substantial number of clinically important ARGs to escape treatment. As new antibiotic resistant strains emerge, careful monitoring of ARGs in both WWTP influent and effluent can identify potentially clinically relevant ARG targets that may become enriched during biological treatment, or simply are released into the environment in significant numbers, which may prompt a closer examination of potential exposure pathways. Attention to pathogens carrying ARGs will be most informative of transmission risks.54,55 At the same time, it is important to recognize that metagenomic detection limits are relatively high56 and likely will not capture rare, but important, evolutionary events contributing to clinically important resistance.21,57 Further, metagenomic monitoring will detect DNA from inactivated organisms, which can still be useful if the goal is informing wastewater-based surveillance (WBS), but not as much so if the goal is assessment of transmission risks. For example, we suspect the detection of resistant N. gonorrheae would likely not reflect a viable pathogen.

The extensive information derived from metagenomic sequencing of WWTP samples yields insights into some of the key mechanisms driving observed resistome shifts through wastewater treatment. By far the dominant factor will be shifts in taxonomic hosts due to changing ecological conditions, making it difficult to discern contributions of HGT, selection, and coselection by antibiotics and other toxic compounds. While ARGs encoded chromosomally and on plasmids can both cause resistance in clinical infections, acquired or mobile ARGs (i.e., those documented to be transferred by HGT) are of particular concern given their ability to spread between species.58 Interestingly, while the mobile fraction of the resistome was reduced in relative abundance as a result of treatment across the 12 WWTPs, the abundance of MGE-associated ARGs increased. While wastewater treatment may drive increases in plasmids due to HGT or vertical transfer, this will not necessarily increase ARGs unless there is subsequent growth of the recipient.59 Several previous studies have demonstrated that microbial community composition strongly shapes environmental resistomes.11,26,60 This study further validated this result, finding that there was a moderate, but significant, correlation between taxonomic and resistome profiles from across the data set. Given that this study compiled data from 12 WWTPs from across the globe, the observed concordance between taxonomy and ARG profile is notable.

Antibiotics persisting throughout the wastewater treatment process have been identified as a potential selective force encouraging the preferential survival of ARB over nonresistant counterparts.25,61 Some of the measured antibiotic quantities were in excess of predicted no effect concentrations (PNEC), suggesting that there is potential for selection within the wastewater environment.62 For example, 92% of influent samples and 58% of effluent samples were in excess of 64 ng/L of the quinolone antibiotic ciprofloxacin, 33% of influent samples and 17% of effluent samples were in excess of 250 ng/L of the macrolide antibiotic clarithromycin, and 17% of influent samples and 8% of effluent samples were in excess of 500 ng/L of the quinolone antibiotic norfloxacin. The percentage of samples in excess of the PNEC for erythromycin (8%) and azithromycin (17%) remained the same from influent to effluent.

This study demonstrates that it is important to examine multiple targets of antibiotic resistance in monitoring campaigns and that it may be misleading to summarize resistome data via a single summary statistic, such as overall removal of total ARGs. While metagenomic sequencing for wastewater resistome characterization is the most comprehensive approach for ARG monitoring, metagenomic sequencing could also be valuable for informing selection of potential ARG monitoring targets.63 It is important to optimize sampling and monitoring strategies based on specific objectives and to ensure that monitoring targets are selected accordingly.64 Potential end goals include:65 (1) characterizing the prevalence of ARB in a given human population via WWTP influent or sewage collection system monitoring (i.e., WBS),13,14,66 (2) monitoring for the risks of evolution of new pathogenic strains of resistant organisms, also involving selection and HGT during biological treatment processes,67,68 (3) assessing the effectiveness of treatment processes for ARG removal via sampling of WWTP effluent, often alongside influent,27,69 and (4) assessing potential human health risks associated with exposures, particularly to resistant pathogens, downstream of WWTP effluent discharge via sampling receiving water bodies.70,71 In WBS and effluent monitoring for human health risk applications, it is important to select clinically relevant ARGs that are associated with pathogens linked to human infection. Linking ARGs to specific hosts is thus necessary. To identify selection pressures that could drive resistance evolution within WWTPs, simple analyses of relative abundances of ARGs is insufficient since increases in ARGs can be caused by taxonomic shifts unrelated to the antibiotic selection pressure.

Monitoring for the emergence of new strains presents numerous challenges. The use of scaffolds assembled from short reads helped identify potential ARG host bacteria. However, the taxonomic origin of ARGs may be impacted by the ability of conjugative transposons to facilitate movement of chromosomal genes to plasmids. When monitoring effluent to assess treatment process effectiveness, ARGs that are prevalent in influent as well as effluent, such as sul1, are strong candidates, though targets that are also prevalent in natural environments and have high background concentrations should be excluded. For all of these goals, metagenomic sequencing is valuable for assessing broad changes in the resistome that may be overlooked when narrowly selecting monitoring targets. When or where a metagenomics approach is not yet feasible or practical, the results of this and other metagenomic surveys can be useful for identifying potential monitoring targets.

Supporting Information Available

The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.est.4c03726.Characteristics of sampled WWTPs; primers and cycling conditions used for qPCR; abundances of the sul1 resistance gene via qPCR versus metagenomics; Shannon diversity in the influent versus final effluent at each WWTP calculated for ARGs and taxonomic classification at the genus level; subset of ARGs that fell within both the 15th percentile of highest influent concentrations and the 15th percentile of reductions in relative abundance during treatment; discriminatory ARGs in influent vs final effluent; key assembled scaffolds of interest; abundance of MGEs; concentrations of antibiotics measured in collected samples; relative abundance of ARGs and corresponding classes of antibiotics in influent versus final effluent; microbial taxonomy in influent and effluent; symmetric Procrustes analysis between microbial community and resistome; percentage mobility (% M) among samples (PDF)

Relative abundance of all ARGs across all samples as determined by metagenomic sequencing; all annotations of ARGs, MGEs, and metal resistance genes on assembled scaffolds; frequency of detection of ARGs on assembled scaffolds; percentage mobility (% M) among ARGs; count of co-occurring scaffold-associated ARGs (XLSX)

Supplementary Material

es4c03726_si_001.pdf

es4c03726_si_002.xlsx

The authors declare no competing financial interest.

Acknowledgments

We wish to thank all waterworks that allowed us to sample at their WWTPs and aided in sample collection. We also thank all who contributed to the collection of samples used in this study: Dr. Caroline Jaraula, Matthew Blair, Matthew Chan, Alexandria Cook, Kris Mapili, Jacob Metch, and Maria Virginia Riquelme. This work was funded by the U.S. National Science Foundation (NSF) Partnerships for International Research and Education Award #1545756, NSF CSSI Award # 2004751, NSF NRT Award #2122798, Water Research Foundation Project 4961, and U.S. Environmental Protection Agency Grant #84082801. H.B. acknowledges additional funding form Swiss National Science Foundation Grants # 167116 and 186531. D.G.J.L. acknowledges additional funding from the Swedish Research Council VR (2022-00945, 2018-05771). C.F.F. acknowledges additional funding from the Swedish Research Council Formas (2018-00833 and 2021-00922).
==== Refs
References

World Health Organization. WHO Strategic Priorities on Antimicrobial Resistance, 2022. https://www.who.int/publications-detail-redirect/9789240041387 (accessed 2023–01–24).
Graham D. W. ; Giesen M. J. ; Bunce J. T. Strategic Approach for Prioritising Local and Regional Sanitation Interventions for Reducing Global Antibiotic Resistance. Water 2018, 11 (1 ), 27 10.3390/w11010027.
Hazra M. ; Durso L. M. Performance Efficiency of Conventional Treatment Plants and Constructed Wetlands towards Reduction of Antibiotic Resistance. Antibiotics 2022, 11 (1 ), 114 10.3390/antibiotics11010114.35052991
Keenum I. ; Calarco J. ; Majeed H. ; Hager-Soto E. E. ; Bott C. ; Garner E. ; Harwood V. J. ; Pruden A. To What Extent Do Water Reuse Treatments Reduce Antibiotic Resistance Indicators? A Comparison of Two Full-Scale Systems. Water Res. 2024, 254 , 121425 10.1016/j.watres.2024.121425.38492480
Marutescu L. G. ; Popa M. ; Gheorghe-Barbu I. ; Barbu I. C. ; Rodríguez-Molina D. ; Berglund F. ; Blaak H. ; Flach C.-F. ; Kemper M. A. ; Spießberger B. ; Wengenroth L. ; Larsson D. G. J. ; Nowak D. ; Radon K. ; de Roda Husman A. M. ; Wieser A. ; Schmitt H. ; Pircalabioru Gradisteanu G. ; Vrancianu C. O. ; Chifiriuc M. C. Wastewater Treatment Plants, an “Escape Gate” for ESCAPE Pathogens. Front. Microbiol. 2023, 14 , 1193907 10.3389/fmicb.2023.1193907.37293232
Bürgmann H. ; Frigon D. ; Gaze W. ; Manaia C. ; Pruden A. ; Singer A. C. ; Smets B. ; Zhang T. Water & Sanitation: An Essential Battlefront in the War on Antimicrobial Resistance. FEMS Microbiol. Ecol. 2018, 94 , fiy101 10.1093/femsec/fiy101.
Manaia C. M. ; Rocha J. ; Scaccia N. ; Marano R. ; Radu E. ; Biancullo F. ; Cerqueira F. ; Fortunato G. ; Iakovides I. C. ; Zammit I. ; Kampouris I. ; Vaz-Moreira I. ; Nunes O. C. Antibiotic Resistance in Wastewater Treatment Plants: Tackling the Black Box. Environ. Int. 2018, 115 (April ), 312–324. 10.1016/j.envint.2018.03.044.29626693
Alexander J. ; Knopp G. ; Dötsch A. ; Wieland A. ; Schwartz T. Ozone Treatment of Conditioned Wastewater Selects Antibiotic Resistance Genes, Opportunistic Bacteria, and Induce Strong Population Shifts. Sci. Total Environ. 2016, 559 , 103–112. 10.1016/j.scitotenv.2016.03.154.27058129
Laht M. ; Karkman A. ; Voolaid V. ; Ritz C. ; Tenson T. ; Virta M. ; Kisand V. Abundances of Tetracycline, Sulphonamide and Beta-Lactam Antibiotic Resistance Genes in Conventional Wastewater Treatment Plants (WWTPs) with Different Waste Load. PLoS One 2014, 9 (8 ), e103705 10.1371/journal.pone.0103705.25084517
Yuan Q. ; Guo M. ; Wei W. ; Yang J. Reductions of Bacterial Antibiotic Resistance through Five Biological Treatment Processes Treated Municipal Wastewater. Environ. Sci. Pollut. Res. 2016, 23 , 19495–19503. 10.1007/s11356-016-7048-8.
Ju F. ; Beck K. ; Yin X. ; Maccagnan A. ; McArdell C. S. ; Singer H. P. ; Johnson D. R. ; Zhang T. ; Bürgmann H. Wastewater Treatment Plant Resistomes Are Shaped by Bacterial Composition, Genetic Exchange, and Upregulated Expression in the Effluent Microbiomes. ISME J. 2019, 13 (2 ), 346–360. 10.1038/s41396-018-0277-8.30250051
Zhang Y. L. ; Marrs C. F. ; Simon C. ; Xi C. W. Wastewater Treatment Contributes to Selective Increase of Antibiotic Resistance among Acinetobacter Spp. Sci. Total Environ. 2009, 407 (12 ), 3702–3706. 10.1016/j.scitotenv.2009.02.013.19321192
Hendriksen R. S. ; Munk P. ; Njage P. ; van Bunnik B. ; Mcnally L. ; Lukjancenko O. ; Röder T. ; Nieuwenhuijse D. ; Pedersen S. K. ; Kjeldgaard J. ; Kaas R. S. ; Clausen P. T. L. C. ; Vogt J. K. ; Leekitcharoenphon P. ; van de Schans M. G. M. ; Zuidema T. ; de Roda Husman A. M. ; Rasmussen S. ; Petersen B. ; Bego A. ; Rees C. ; Cassar S. ; Coventry K. ; Collignon P. ; Allerberger F. ; Rahube T. O. ; Oliveira G. ; Ivanov I. ; Vuthy Y. ; Sopheak T. ; Yost C. K. ; Ke C. ; Zheng H. ; Baisheng L. ; Jiao X. ; Donado-Godoy P. ; Coulibaly K. J. ; Jergović M. ; et al. Global Monitoring of Antimicrobial Resistance Based on Metagenomics Analyses of Urban Sewage. Nat. Commun. 2019, 10 (1 ), 1124 10.1038/s41467-019-08853-3.30850636
Prieto Riquelme M. V. ; Garner E. ; Gupta S. ; Metch J. ; Zhu N. ; Blair M. F. ; Arango-Argoty G. ; Maile-Moskowitz A. ; Li A. ; Flach C.-F. ; Aga D. S. ; Nambi I. M. ; Larsson D. G. J. ; Bürgmann H. ; Zhang T. ; Pruden A. ; Vikesland P. J. Demonstrating a Comprehensive Wastewater-Based Surveillance Approach That Differentiates Globally Sourced Resistomes. Environ. Sci. Technol. 2022, 56 (21 ), 14982–14993. 10.1021/acs.est.1c08673.35759608
Vikesland P. ; Garner E. ; Gupta S. ; Kang S. ; Maile-Moskowitz A. ; Zhu N. Differential Drivers of Antimicrobial Resistance across the World. Acc. Chem. Res. 2019, 52 (4 ), 916–924. 10.1021/acs.accounts.8b00643.30848890
Kim S. ; Yun Z. ; Ha U.-H. ; Lee S. ; Park H. ; Kwon E. E. ; Cho Y. ; Choung S. ; Oh J. ; Medriano C. A. ; Chandran K. Transfer of Antibiotic Resistance Plasmids in Pure and Activated Sludge Cultures in the Presence of Environmentally Representative Micro-Contaminant Concentrations. Sci. Total Environ. 2014, 468–469 , 813–820. 10.1016/j.scitotenv.2013.08.100.
Prudhomme M. ; Attaiech L. ; Sanchez G. ; Martin B. ; Claverys J. P. Antibiotic Stress Induces Genetic Transformability in the Human Pathogen Streptococcus Pneumoniae. Science 2006, 313 (5783 ), 89–92. 10.1126/science.1127912.16825569
Jiao Y. N. ; Chen H. ; Gao R. X. ; Zhu Y. G. ; Rensing C. Organic compounds stimulate horizontal transfer of antibiotic resistance genes in mixed wastewater treatment systems. Chemosphere 2017, 184 (June ), 53–61. 10.1016/j.chemosphere.2017.05.149.28578196
Zhang P. Y. ; Xu P. P. ; Xia Z. J. ; Wang J. ; Xiong J. ; Li Y. Z. Combined Treatment with the Antibiotics Kanamycin and Streptomycin Promotes the Conjugation of Escherichia Coli. FEMS Microbiol. Lett. 2013, 348 (2 ), 149–156. 10.1111/1574-6968.12282.24111668
Zhang Y. ; Gu A. Z. ; Cen T. ; Li X. ; He M. ; Li D. ; Chen J. Sub-Inhibitory Concentrations of Heavy Metals Facilitate the Horizontal Transfer of Plasmid-Mediated Antibiotic Resistance Genes in Water Environment. Environ. Pollut. 2018, 237 , 74–82. 10.1016/j.envpol.2018.01.032.29477117
Berglund F. ; Ebmeyer S. ; Kristiansson E. ; Larsson D. G. J. Evidence for Wastewaters as Environments Where Mobile Antibiotic Resistance Genes Emerge. Commun. Biol. 2023, 6 (1 ), 321–411. 10.1038/s42003-023-04676-7.36966231
Jin M. ; Lu J. ; Chen Z. ; Nguyen S. H. ; Mao L. ; Li J. ; Yuan Z. ; Guo J. Antidepressant Fluoxetine Induces Multiple Antibiotics Resistance in Escherichia Coli via ROS-Mediated Mutagenesis. Environ. Int. 2018, 120 (July ), 421–430. 10.1016/j.envint.2018.07.046.30125859
Kurenbach B. ; Marjoshi D. ; Amábile-Cuevas C. F. ; Ferguson G. C. ; Godsoe W. ; Gibson P. ; Heinemann J. a. Sublethal Exposure to Commercial Formulations of the Herbicides Dicamba, 2,4-Dichlorophenoxyacetic Acid, and Glyphosate Cause Changes in Antibiotic Susceptibility in Escherichia coli and Salmonella enterica serovar Typhimurium. mBio 2015, 6 (2 ), e0009–e0015. 10.1128/mBio.00009-15.
Baker-Austin C. ; Wright M. S. ; Stepanauskas R. ; McArthur J. V. Co-Selection of Antibiotic and Metal Resistance. Trends Microbiol. 2006, 14 (4 ), 176–182. 10.1016/j.tim.2006.02.006.16537105
Rizzo L. ; Manaia C. ; Merlin C. ; Schwartz T. ; Dagot C. ; Ploy M. C. ; Michael I. ; Fatta-Kassinos D. Urban Wastewater Treatment Plants as Hotspots for Antibiotic Resistant Bacteria and Genes Spread into the Environment: A Review. Sci. Total Environ. 2013, 447 , 345–360. 10.1016/j.scitotenv.2013.01.032.23396083
Forsberg K. J. ; Patel S. ; Gibson M. K. ; Lauber C. L. ; Knight R. ; Fierer N. ; Dantas G. Bacterial Phylogeny Structures Soil Resistomes across Habitats. Nature 2014, 509 , 612–616. 10.1038/nature13377.24847883
Majeed H. J. ; Riquelme M. V. ; Davis B. C. ; Gupta S. ; Angeles L. ; Aga D. S. ; Garner E. ; Pruden A. ; Vikesland P. J. Evaluation of Metagenomic-Enabled Antibiotic Resistance Surveillance at a Conventional Wastewater Treatment Plant. Front. Microbiol. 2021, 12 , 657954 10.3389/fmicb.2021.657954.34054755
Alcock B. P. ; Huynh W. ; Chalil R. ; Smith K. W. ; Raphenya A. R. ; Wlodarski M. A. ; Edalatmand A. ; Petkau A. ; Syed S. A. ; Tsang K. K. ; Baker S. J. C. ; Dave M. ; McCarthy M. C. ; Mukiri K. M. ; Nasir J. A. ; Golbon B. ; Imtiaz H. ; Jiang X. ; Kaur K. ; Kwong M. ; Liang Z. C. ; Niu K. C. ; Shan P. ; Yang J. Y. J. ; Gray K. L. ; Hoad G. R. ; Jia B. ; Bhando T. ; Carfrae L. A. ; Farha M. A. ; French S. ; Gordzevich R. ; Rachwalski K. ; Tu M. M. ; Bordeleau E. ; Dooley D. ; Griffiths E. ; Zubyk H. L. ; Brown E. D. ; Maguire F. ; Beiko R. G. ; Hsiao W. W. L. ; Brinkman F. S. L. ; Van Domselaar G. ; McArthur A. G. CARD 2023: Expanded Curation, Support for Machine Learning, and Resistome Prediction at the Comprehensive Antibiotic Resistance Database. Nucleic Acids Res. 2023, 51 (D1 ), D690–D699. 10.1093/nar/gkac920.36263822
Brinch C. ; Leekitcharoenphon P. ; Duarte A. S. R. ; Svendsen C. A. ; Jensen J. D. ; Aarestrup F. M. Long-Term Temporal Stability of the Resistome in Sewage from Copenhagen. mSystems 2020, 5 (5 ), 10 10.1128/mSystems.00841-20.
Singh R. R. ; Angeles L. F. ; Butryn D. M. ; Metch J. W. ; Garner E. ; Vikesland P. J. ; Aga D. S. Towards a Harmonized Method for the Global Reconnaissance of Multi-Class Antimicrobials and Other Pharmaceuticals in Wastewater and Receiving Surface Waters. Environ. Int. 2019, 124 (November 2018 ), 361–369. 10.1016/j.envint.2019.01.025.30660849
Li A.-D. ; Metch J. W. ; Wang Y. ; Garner E. ; Zhang A. N. ; Riquelme M. V. ; Vikesland P. J. ; Pruden A. ; Zhang T. Effects of Sample Preservation and DNA Extraction on Enumeration of Antibiotic Resistance Genes in Wastewater. FEMS Microbiol. Ecol. 2018, 94 , 1–11. 10.1093/femsec/fix189.
Pei R. T. ; Kim S. C. ; Carlson K. H. ; Pruden A. Effect of River Landscape on the Sediment Concentrations of Antibiotics and Corresponding Antibiotic Resistance Genes (ARG). Water Res. 2006, 40 (12 ), 2427–2435. 10.1016/j.watres.2006.04.017.16753197
Suzuki M. T. ; Taylor L. T. ; DeLong E. F. Quantitative Analysis of Small-Subunit rRNA Genes in Mixed Microbial Populations via 5 ’-Nuclease Assays. Appl. Environ. Microbiol. 2000, 66 (11 ), 4605–4614. 10.1128/AEM.66.11.4605-4614.2000.11055900
Arango-Argoty G. ; Singh G. ; Heath L. S. ; Pruden A. ; Xiao W. ; Zhang L. MetaStorm: A Public Resource for Customizable Metagenomics Annotation. PLoS One 2016, 11 (9 ), e01624422 10.1371/journal.pone.0162442.
Bolger A. M. ; Lohse M. ; Usadel B. Trimmomatic: A Flexible Trimmer for Illumina Sequence Data. Bioinformatics 2014, 30 (15 ), 2114–2120. 10.1093/bioinformatics/btu170.24695404
Peng Y. ; Leung H. C. M. ; Yiu S. M. ; Chin F. Y. L. IDBA-UD: A de Novo Assembler for Single-Cell and Metagenomic Sequencing Data with Highly Uneven Depth. Bioinformatics 2012, 28 (11 ), 1420–1428. 10.1093/bioinformatics/bts174.22495754
Jia B. ; Raphenya A. R. ; Alcock B. ; Waglechner N. ; Guo P. ; Tsang K. K. ; Lago B. A. ; Dave B. M. ; Pereira S. ; Sharma A. N. ; Doshi S. ; Courtot M. ; Lo R. ; Williams L. E. ; Frye J. G. ; Elsayegh T. ; Sardar D. ; Westman E. L. ; Pawlowski A. C. ; Johnson T. A. ; Brinkman F. S. L. ; Wright G. D. ; McArthur A. G. CARD 2017: Expansion and Model-Centric Curation of the Comprehensive Antibiotic Resistance Database. Nucleic Acids Res. 2016, 45 , D566–D573. 10.1093/nar/gkw1004.27789705
Quast C. ; Pruesse E. ; Yilmaz P. ; Gerken J. ; Schweer T. ; Yarza P. ; Peplies J. ; Glockner F. O. The SILVA Ribosomal RNA Gene Database Project: Improved Data Processing and Web-Based Tools. Nucleic Acids Res. 2012, 41 , D590–D596. 10.1093/nar/gks1219.23193283
Leplae R. ; Hebrant A. ; Wodak S. J. ; Toussaint A. ACLAME: A CLAssification of Mobile Genetic Elements. Nucleic Acids Res. 2004, 32 (90001 ), D45–D49. 10.1093/nar/gkh084.14681355
Segata N. ; Waldron L. ; Ballarini A. ; Narasimhan V. ; Jousson O. ; Huttenhower C. Metagenomic Microbial Community Profiling Using Unique Clade-Specific Marker Genes. Nat. Methods 2012, 9 (8 ), 811–814. 10.1038/nmeth.2066.22688413
Li B. ; Yang Y. ; Ma L. ; Ju F. ; Guo F. ; Tiedje J. M. ; Zhang T. Metagenomic and Network Analysis Reveal Wide Distribution and Co-Occurrence of Environmental Antibiotic Resistance Genes. ISME J. 2015, 9 , 2490–2502. 10.1038/ismej.2015.59.25918831
Zankari E. ; Hasman H. ; Cosentino S. ; Vestergaard M. ; Rasmussen S. ; Lund O. ; Aarestrup F. M. ; Larsen M. V. Identification of Acquired Antimicrobial Resistance Genes. J. Antimicrob. Chemother. 2012, 67 (11 ), 2640–2644. 10.1093/jac/dks261.22782487
Arango-Argoty G. A. ; Dai D. ; Pruden A. ; Vikesland P. ; Heath L. S. ; Zhang L. NanoARG: A Web Service for Detecting and Contextualizing Antimicrobial Resistance Genes from Nanopore-Derived Metagenomes. Microbiome 2019, 7 (1 ), 88 10.1186/s40168-019-0703-9.31174603
R Core Team R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2017.
Kindt R. ; Legendre P. ; Hara B. O. ; Simpson G. L. ; Stevens M. H. H. ; Wagner H. The Vegan Package, 2008; Vol. 10 ( (01 ), ).
Segata N. ; Izard J. ; Waldron L. ; Gevers D. ; Miropolsky L. ; Garrett W. S. ; Huttenhower C. Metagenomic Biomarker Discovery and Explanation. Genome Biol. 2011, 12 (6 ), R60 10.1186/gb-2011-12-6-r60.21702898
Wickham H. Ggplot2: Elegant Graphics for Data Analysis; Springer-Verlag: New York, 2016.
Arya V. ; Philip L. ; Murty Bhallamudi S. Performance of Suspended and Attached Growth Bioreactors for the Removal of Cationic and Anionic Pharmaceuticals. Chem. Eng. J. 2016, 284 , 1295–1307. 10.1016/j.cej.2015.09.070.
Petrovich M. ; Rosenthal A. ; Griffin J. ; Wells G. Spatially Resolved Abundances of Antibiotic Resistance Genes and intI1 in Wastewater Treatment Biofilms. Biotechnol. Bioeng. 2019, 116 (3 ), 543–554. 10.1002/bit.26887.30512194
Zhang A.-N. ; Gaston J. M. ; Dai C. L. ; Zhao S. ; Poyet M. ; Groussin M. ; Yin X. ; Li L.-G. ; van Loosdrecht M. C. M. ; Topp E. ; Gillings M. R. ; Hanage W. P. ; Tiedje J. M. ; Moniz K. ; Alm E. J. ; Zhang T. An Omics-Based Framework for Assessing the Health Risk of Antimicrobial Resistance Genes. Nat. Commun. 2021, 12 (1 ), 4765 10.1038/s41467-021-25096-3.34362925
Martinez J. L. ; Sánchez M. B. ; Martínez-Solano L. ; Hernandez A. ; Garmendia L. ; Fajardo A. ; Alvarez-Ortega C. Functional Role of Bacterial Multidrug Efflux Pumps in Microbial Natural Ecosystems. FEMS Microbiol. Rev. 2009, 33 (2 ), 430–449. 10.1111/j.1574-6976.2008.00157.x.19207745
Parsley L. C. ; Consuegra E. J. ; Kakirde K. S. ; Land A. M. ; Harper W. F. ; Liles M. R. Identification of Diverse Antimicrobial Resistance Determinants Carried on Bacterial, Plasmid, or Viral Metagenomes from an Activated Sludge Microbial Assemblage. Appl. Environ. Microbiol. 2010, 76 (11 ), 3753–3757. 10.1128/AEM.03080-09.20382816
Yang Y. ; Li B. ; Ju F. ; Zhang T. Exploring Variation of Antibiotic Resistance Genes in Activated Sludge over a Four-Year Period through a Metagenomic Approach. Environ. Sci. Technol. 2013, 47 (18 ), 10197–10205. 10.1021/es4017365.23919449
Ashbolt N. J. ; Amézquita A. ; Backhaus T. ; Borriello P. ; Brandt K. K. ; Collignon P. ; Coors A. ; Finley R. ; Gaze W. H. ; Heberer T. ; Lawrence J. R. ; Larsson D. J. ; McEwen S. A. ; Ryan J. J. ; Schönfeld J. ; Silley P. ; Snape J. R. ; Van den Eede C. ; Topp E. Human Health Risk Assessment (HHRA) for Environmental Development and Transfer of Antibiotic Resistance. Environ. Health Perspect. 2013, 121 (9 ), 993–1001. 10.1289/ehp.1206316.23838256
Garner E. ; Organiscak M. ; Dieter L. ; Shingleton C. ; Haddix M. ; Joshi S. ; Pruden A. ; Ashbolt N. ; Medema G. ; Hamilton K. Towards Risk Assessment for Antibiotic Resistant Pathogens in Recycled Water: A Systematic Review and Summary of Research Needs. Environ. Microbiol. 2021, 23 (12 ), 7355–7372. 10.1111/1462-2920.15804.34632683
Davis B. C. ; Brown C. ; Gupta S. ; Calarco J. ; Liguori K. ; Milligan E. ; Harwood V. J. ; Pruden A. ; Keenum I. Recommendations for the Use of Metagenomics for Routine Monitoring of Antibiotic Resistance in Wastewater and Impacted Aquatic Environments. Crit. Rev. Environ. Sci. Technol. 2023, 53 (19 ), 1731–1756. 10.1080/10643389.2023.2181620.
Larsson D. G. J. ; Flach C.-F. Antibiotic Resistance in the Environment. Nat. Rev. Microbiol. 2022, 20 (5 ), 257–269. 10.1038/s41579-021-00649-x.34737424
von Wintersdorff C. J. H. ; Penders J. ; van Niekerk J. M. ; Mills N. D. ; Majumder S. ; van Alphen L. B. ; Savelkoul P. H. M. ; Wolffs P. F. G. Dissemination of Antimicrobial Resistance in Microbial Ecosystems through Horizontal Gene Transfer. Front. Microbiol. 2016, 7 , 173 10.3389/fmicb.2016.00173.26925045
Perry J. ; Wright G. The Antibiotic Resistance “Mobilome”: Searching for the Link between Environment and Clinic. Front. Microbiol. 2013, 4 , 138 10.3389/fmicb.2013.00138.23755047
Wang C. ; Dong D. ; Strong P. J. ; Zhu W. ; Ma Z. ; Qin Y. ; Wu W. Microbial Phylogeny Determines Transcriptional Response of Resistome to Dynamic Composting Processes. Microbiome 2017, 5 (1 ), 103 10.1186/s40168-017-0324-0.28814344
Bengtsson-Palme J. ; Hammarén R. ; Pal C. ; Östman M. ; Björlenius B. ; Flach C. F. ; Fick J. ; Kristiansson E. ; Tysklind M. ; Larsson D. G. J. Elucidating Selection Processes for Antibiotic Resistance in Sewage Treatment Plants Using Metagenomics. Sci. Total Environ. 2016, 572 , 697–712. 10.1016/j.scitotenv.2016.06.228.27542633
Bengtsson-Palme J. ; Larsson D. J. Concentrations of Antibiotics Predicted to Select for Resistant Bacteria: Proposed Limits for Environmental Regulation. Environ. Int. 2016, 86 , 140–149. 10.1016/j.envint.2015.10.015.26590482
Tarek M. ; Garner E. A Proposed Framework for the Identification of Indicator Genes for Monitoring Antibiotic Resistance in Wastewater: Insights from Metagenomic Sequencing. Sci. Total Environ. 2023, 854 , 158698 10.1016/j.scitotenv.2022.158698.36108825
Larsson D. G. J. ; Flach C.-F. ; Laxminarayan R. Sewage Surveillance of Antibiotic Resistance Holds Both Opportunities and Challenges. Nat. Rev. Microbiol. 2023, 21 (4 ), 213–214. 10.1038/s41579-022-00835-5.36470999
Liguori K. ; Keenum I. ; Davis B. C. ; Calarco J. ; Milligan E. ; Harwood V. J. ; Pruden A. Antimicrobial Resistance Monitoring of Water Environments: A Framework for Standardized Methods and Quality Control. Environ. Sci. Technol. 2022, 56 (13 ), 9149–9160. 10.1021/acs.est.1c08918.35732277
Karkman A. ; Berglund F. ; Flach C.-F. ; Kristiansson E. ; Larsson D. G. J. Predicting Clinical Resistance Prevalence Using Sewage Metagenomic Data. Commun. Biol. 2020, 3 (1 ), 711–810. 10.1038/s42003-020-01439-6.33244050
Hutinel M. ; Fick J. ; Larsson D. J. ; Flach C.-F. Investigating the Effects of Municipal and Hospital Wastewaters on Horizontal Gene Transfer. Environ. Pollut. 2021, 276 , 116733 10.1016/j.envpol.2021.116733.33631686
Wei Z. ; Feng K. ; Wang Z. ; Zhang Y. ; Yang M. ; Zhu Y.-G. ; Virta M. P. J. ; Deng Y. High-Throughput Single-Cell Technology Reveals the Contribution of Horizontal Gene Transfer to Typical Antibiotic Resistance Gene Dissemination in Wastewater Treatment Plants. Environ. Sci. Technol. 2021, 55 (17 ), 11824–11834. 10.1021/acs.est.1c01250.34415164
Mao D. ; Yu S. ; Rysz M. ; Luo Y. ; Yang F. ; Li F. ; Hou J. ; Mu Q. ; Alvarez P. J. J. Prevalence and Proliferation of Antibiotic Resistance Genes in Two Municipal Wastewater Treatment Plants. Water Res. 2015, 85 , 458–466. 10.1016/j.watres.2015.09.010.26372743
Rodriguez-Mozaz S. ; Chamorro S. ; Marti E. ; Huerta B. ; Gros M. ; Sànchez-Melsió A. ; Borrego C. M. ; Barceló D. ; Balcázar J. L. Occurrence of Antibiotics and Antibiotic Resistance Genes in Hospital and Urban Wastewaters and Their Impact on the Receiving River. Water Res. 2015, 69 , 234–242. 10.1016/j.watres.2014.11.021.25482914
Xu J. ; Xu Y. ; Wang H. ; Guo C. ; Qiu H. ; He Y. ; Zhang Y. ; Li X. ; Meng W. Occurrence of Antibiotics and Antibiotic Resistance Genes in a Sewage Treatment Plant and Its Effluent-Receiving River. Chemosphere 2015, 119 , 1379–1385. 10.1016/j.chemosphere.2014.02.040.24630248
